** : push-over analysis for sand model, paqajsandp
*heading
4-leg jack up on sand cylindrical
**
*restart,write,f=20
**node definitions
**
*node,nset=beams
1,0,0
2,0,59.0
3,29.33,0
4,29.33,59.0
5,58.66,0.0
6,58.66,59.0
999,29.33,59.0
33,5,0.
34,6,0.
*node,nset=base
10, 0.  
30, 29.33
50, 58.66
*nset,nset=epj
base,1,3,5
**
**element definitons
**
*element,elset=beaml,type=b21
1,1,2
2,3,4
3,3,4
4,5,6
*element,elset=rbeams,type=rb2d2
11,2,4
12,4,6
*element,type=joint2d, elset=joint2d
101,10,1
102,30,3
103,30,3
104,50,5
**
**
*element,type=t2d2,elset=fake
333,33,34
*solid section,elset=fake,material=m1
1., 
*material,name=m1
*elastic
1.,0.
**
**
**orientations and transformations
**
*orientation,name=ori1
0,1,0,1,0,0
*transform,nset=epj
0,1,0,1,0,0
**
**property definitions
**
*beam general section,elset=beaml,section=general
1.0,0.33333,0.,0.33333,0.6666

3.0e15,2.0e15
*rigid body,elset=rbeams,ref node=999
**
**
*epjoint,elset=joint2d,orien=ori1,section=spud
14.0,
*joint elasti,ndim=2,moduli=spud can
1.63e5,2.92e4,2.10e4,0.2
*joint plasti,type=sand
40000,0.3,0.3,35.0,10.0
**
**
*initial conditions,type=spud preload
joint2d, 775000
**
**
*boundary
33,1,2,0
34,1,2,0.
**
**
**analysis 4 leg jack up
**
**
**
*step,inc=1000
 apply operational vertical load
*static
1.,1.
*print,plasticity=yes
*boundary
10,1,6,0
50,1,6,0
30,1,6,0
*output,hist,freq=0
*output,field,freq=20
*element output,elset=joint2d
pe,
e,
*node output,nset=epj
u
*cload
4,2,-209000
*end step
**
**
**
*step,inc=10000
*static,direct
.01,1
**
**
**** -- load beyond maximum
***cload
**999,1,200000
***
** -- load to 95% of max
*cload
999,1,142500
*boundary
34,1,1,2.
*end step