*heading
clay rig
*restart,write
**node definitions
*node
1,0,0,0.
2,0,110.6,0.
3,55,0.0,0.
4,55,110.6,0.
99,18,110.6,0.
10,0,0,0.
30,55,0.0,0.
999,18,110.6,0
*nset,nset=epj
1,3,10,30
*nset,nset=base
10,30
*element,type=b23,elset=beams
1,1,2
3,3,4
*element,elset=rbeams,type=rb2d2
901,2,99
902,4,99
*element,elset=joint2d,type=joint2d
101,10,1
103,30,3
**
**
*beam general section,elset=beams,section=general
1.0,1.6104e-06,0.,1.6104e-06,3.2208e-06

1.54e15,3.61e15
**
**
*rigid body,elset=rbeams,ref node=999
**
**orientations and transformations
**
*orientation,name=ori1
0,1,0,-1,0,0
*transform,nset=epj
0,1,0,-1,0,0
**
**
**property definitions
**
*epjoint,elset=joint2d,orien=ori1,section=spud
20.0,0.0
*joint elasti,ndim=2,moduli=spud can
1.56e4,2.34e4,6.38e3,0.5
*joint plasti,type=member
100,100,1.e15,15.
**
**
**
*step,inc=1000
vertical loading
*static
0.1,5.0
*boundary
10,1,6,0
30,1,6,0
*cload
99,2,-100
*el file,elset=joint2d
e, 
s, 
pe, 
*OUTPUT,FIELD
*ELEMENT OUTPUT,elset=joint2d
e,
s,
pe,
*node file
u, 
rf, 
*OUTPUT,FIELD
*NODE OUTPUT
u,
rf,
*el print,elset=joint2d
s, 
pe, 
*end step
**
**
*step,inc=1000
monotonic horizontal loading
*static
0.1,1.0
*cload
99,1,14.
*end step
**
**
**
**