*heading
critical load - initial embedment checks for
a. sand model (cylindrical, conical < critical embed, conical > critical embed
b. clay model (same 3 cases)
*****************************************************************************
                     Explanation for node/element set names
The first three letters indicate whether the material is sand(snd) or clay(cly)
The next three letters indicate the section of the spud can: cylindrical(cyl)
or conical(con).
The seventh letter, indicates wheter the embedment is greater than (g) or less
than (l) critical (for conical spud cans only).
If the last letter is `f', it indicates that the material properties for that
element are field variable dependent.
*****************************************************************************
**
**node definitions
**
*node,nset=sndcyl
1,1.0
2,1.0
*node,nset=sndcong
21,2.0
22,2.0
*node,nset=sndconl
31,3.0
32,3.0
*node,nset=clycyl
41,4.0
42,4.0
*node,nset=clycong
51,5.0
52,5.0
*node,nset=clyconl
61,6.0
62,6.0
*node,nset=sndcylf
101,7.0
102,7.0
*node,nset=sndcongf
121,8.0
122,8.0
*node,nset=sndconlf
131,9.0
132,9.0
*node,nset=clycylf
141,10.0
142,10.0
*node,nset=clycongf
151,11.0
152,11.0
*node,nset=clyconlf
161,12.0
162,12.0
**
**element definitions
**
**
*element,type=joint3d, elset=sndcyl
11,1,2
*element,type=joint3d, elset=sndcong
21,21,22
*element,type=joint3d, elset=sndconl
31,31,32
*element,type=joint3d, elset=clycyl
41,41,42
*element,type=joint3d, elset=clycong
51,51,52
*element,type=joint3d, elset=clyconl
61,61,62
*element,type=joint3d, elset=sndcylf
111,101,102
*element,type=joint3d, elset=sndcongf
121,121,122
*element,type=joint3d, elset=sndconlf
131,131,132
*element,type=joint3d, elset=clycylf
141,141,142
*element,type=joint3d, elset=clycongf
151,151,152
*element,type=joint3d, elset=clyconlf
161,161,162
**
**
**nset definitions 
**
*nset,nset=nall
sndcyl,sndcong,sndconl,clycyl,clycong,clyconl,sndcylf,sndcongf
*nset,nset=nall
nall,sndconlf,clycylf,clycongf,clyconlf
**
**elset definitions
**
**orientations
**
*orientation,name=ori1
0,1,0,1,0,0
**transformations
**
*transform,nset=nall
0,1,0,1,0,0
**
**property definitions
**
**sand model ... cylindrical
**
*epjoint,elset=sndcyl,orien=ori1,section=spud
10.9,0.
*joint elasticity,ndim=3,moduli=spud can
5.14e4, 3.87e3, 2.04e4, 0.2, 2.e4
*joint plasticity,type=sand
0.,1,0.5,33.0,10.0
**
*epjoint,elset=sndcylf,orien=ori1,section=spud
10.9,0.
*joint elasticity,ndim=3,moduli=spud can
5.14e4, 3.87e3, 2.04e4, 0.2, 2.e4
*joint plasticity,type=sand,dependencies=1
0.,1,0.4,33.0,10.0,30.,2.
0.,1,0.5,33.0,10.0,30.,3.
**
*initial conditions, type=spudpreload
sndcyl, 50600
sndcylf, 50600
**
**
**sand model... conical > critical
*epjoint,elset=sndcong,orien=ori1,section=spud
10.9,150.0
*joint elasticity,ndim=3,moduli=spud can
5.14e4, 3.87e3, 2.04e4, 0.2, 2.e4
*joint plasticity,type=sand
0.,1,0.5,33.0,10.0
**
*epjoint,elset=sndcongf,orien=ori1,section=spud
10.9,150.0
*joint elasticity,ndim=3,moduli=spud can
5.14e4, 3.87e3, 2.04e4, 0.2, 2.e4
*joint plasticity,type=sand,dependencies=1
0.,1,0.4,33.0,10.0,30.,2.
0.,1,0.5,33.0,10.0,30.,3.
**
**
*initial conditions, type=spudpreload
sndcong, 50000.
sndcongf, 50000.
**
**
**sand model... conical < critical
*epjoint,elset=sndconl,orien=ori1,section=spud
10.9,150.0
*joint elasticity,ndim=3,moduli=spud can
5.14e4, 3.87e3, 2.04e4, 0.2, 2.e4
*joint plasticity,type=sand
0.,1,0.5,33.0,10.0
** 
*epjoint,elset=sndconlf,orien=ori1,section=spud
10.9,150.0
*joint elasticity,ndim=3,moduli=spud can
5.14e4, 3.87e3, 2.04e4, 0.2, 2.e4
*joint plasticity,type=sand,dependencies=1
0.,1,0.4,33.0,10.0,30.,2.
0.,1,0.5,33.0,10.0,30.,3.
**
*initial conditions,type=spudpreload
sndconl, 15000.,
sndconlf, 15000.,
**
**
**clay model : cylindrical
**
*epjoint,elset=clycyl,orien=ori1,section=spud
20.0,0.0
*joint elasticity,ndim=3,moduli=spud can
1.56e4,2.34e4,6.38e3,0.5,2.e4
*joint plasticity,type=clay
150,7.204e4,1.978e3
**
*epjoint,elset=clycylf,orien=ori1,section=spud
20.0,0.0
*joint elasticity,ndim=3,moduli=spud can
1.56e4,2.34e4,6.38e3,0.5,2.e4
*joint plasticity,type=clay,dependencies=1
150,8.204e4,1.978e3,,30.,2.
150,7.204e4,1.978e3,,30.,3.
**
*initial conditions,type=spud preload
clycyl, 1.3e5,
clycylf, 1.3e5,
**
**
**clay model > critical
**
*epjoint,elset=clycong,orien=ori1,section=spud
20.0,150.0,,8.5e5
*joint elasticity,ndim=3,moduli=spud can
1.56e4,2.34e4,6.38e3,0.5,2.e4
*joint plasticity,type=clay
50.0,-2.395e-5,8.77e-6,2.9294
**
*epjoint,elset=clycongf,orien=ori1,section=spud
20.0,150.0,,8.5e5
*joint elasticity,ndim=3,moduli=spud can
1.56e4,2.34e4,6.38e3,0.5,2.e4
*joint plasticity,type=clay,dependencies=1
50.0,-3.395e-5,8.77e-6,2.9294,30.,2.
50.0,-2.395e-5,8.77e-6,2.9294,30.,3.
**
*initial conditions,type=spud preload
clycong, 8.5e5,
clycongf, 8.5e5,
**
**
**clay model < critical
**
*epjoint,elset=clyconl,orien=ori1,section=spud
20.0,150.0,,1.3e5
*joint elasticity,ndim=3,moduli=spud can
1.56e4,2.34e4,6.38e3,0.5,2.e4
*joint plasticity,type=clay
50.0,-2.395e-5,8.77e-6,2.9294
**
*epjoint,elset=clyconlf,orien=ori1,section=spud
20.0,150.0,,1.3e5
*joint elasticity,ndim=3,moduli=spud can
1.56e4,2.34e4,6.38e3,0.5,2.e4
*joint plasticity,type=clay,dependencies=1
50.0,-3.395e-5,8.77e-6,2.9294,30.,2.
50.0,-2.395e-5,8.77e-6,2.9294,30.,3.
**
*initial conditions,type=spud preload
clyconl, 1.3e5,
clyconlf, 1.3e5,
**
*initial conditions, type=field,variable=1
nall,3.
**
**
*step
steps to check whether peeq output correctly, and that embedment
has the correct effect on elastic modulus.
*static
1.0,1.
*boundary
1,1,1,0
21,1,1,0
31,1,6,0
41,1,1,0
51,1,1,0
61,1,6,0
101,1,1,0
121,1,1,0
131,1,6,0
141,1,1,0
151,1,1,0
161,1,6,0
*cload
32,1,-1.2312e4
32,2,-0.0879e4
32,3,-0.0879e4
32,4,-200.0
32,5,-7.48412e3
32,6,-7.48412e2
62,1,-0.6125e3
62,2,-0.6125e2
62,3,-0.6125e3
62,4,-200.0
62,5,-4.01372e2
62,6,-4.01372e2
132,1,-1.2312e4
132,2,-0.0879e4
132,3,-0.0879e4
132,4,-200.0
132,5,-7.48412e3
132,6,-7.48412e2
162,1,-0.6125e3
162,2,-0.6125e2
162,3,-0.6125e3
162,4,-200.0
162,5,-4.01372e2
162,6,-4.01372e2
*elset,elset=check
clyconl,clyconlf,sndconl,sndconlf
*nset,nset=check
32,62,132,162
*elprint,elset=check
s,
e,
*elprint
peeq, 
*elfile,elset=check
s,
e, 
*OUTPUT,FIELD
*ELEMENT OUTPUT,elset=check
s,
e,
*elfile
pe, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
pe,
*nodefile,nset=check
u, 
*OUTPUT,FIELD
*NODE OUTPUT,nset=check
u,
*nodefile
rf, 
*OUTPUT,FIELD
*NODE OUTPUT
rf,
*node print,nset=check
u, 
rf, 
*end step
**static perturbation
*step,perturbation
*static
*cload
32,1,-1.2312e4
32,2,-0.0879e4
32,3,-0.0879e4
32,4,-200.0
32,5,-7.48412e3
32,6,-7.48412e2
62,1,-0.6125e3
62,2,-0.6125e2
62,3,-0.6125e3
62,4,-200.0
62,5,-4.01372e2
62,6,-4.01372e2
132,1,-1.2312e4
132,2,-0.0879e4
132,3,-0.0879e4
132,4,-200.0
132,5,-7.48412e3
132,6,-7.48412e2
162,1,-0.6125e3
162,2,-0.6125e2
162,3,-0.6125e3
162,4,-200.0
162,5,-4.01372e2
162,6,-4.01372e2
*elset,elset=check
clyconl,clyconlf,sndconl,sndconlf
*nset,nset=check
32,62,132,162
*elprint,elset=check
s,
e,
*elprint
peeq, 
*elfile,elset=check
s,
e, 
*OUTPUT,FIELD
*ELEMENT OUTPUT,elset=check
s,
e,
*elfile
pe, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
pe,
*nodefile,nset=check
u, 
rf, 
*OUTPUT,FIELD
*NODE OUTPUT,nset=check
u,
rf,
*node print,nset=check
u, 
rf, 
*end step