*heading
clay rig
**node definitions
*node
1,0.0,0.0,31.75
331,0,0,-31.75
2,0,110.6,31.75
332,0,110.6,-31.75
3,55,0.0,0.
4,55,110.6,0.
99,18,110.6,0.
100,18,55.6,0.
10,0,0,31.75
20,0,0,-31.75
30,55,0.0,0.
999,18,82.5,0.
*nset,nset=epj
1,3,10,30,20
*nset,nset=base
10,30,20
*element,type=b31,elset=beams
1,1,2
2,331,332
3,3,4
*element,elset=rbeams,type=rb3d2
901,2,99
904,332,99
902,4,99
903,99,100
*element,elset=joint2d,type=joint3d
101,10,1
102,20,331
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
**
**
**property definitions
**
*epjoint,elset=joint2d,orien=ori1,section=spud
20.0,0.0
*joint elasti,ndim=3,moduli=spud can
1.56e4,2.34e4,6.38e3,0.5,40000
*joint plasti,type=clay
150,7.204e4,1.978e3
**
*initial conditions, type=spud preload
joint2d,1.3e5
*step
vertical loading
*static
0.1,1.0
*boundary
10,1,6,0
20,1,6,0
30,1,6,0
*cload
99,2,-198000
*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.001,1.,,0.05
*cload
100,1,100000
*end step