*heading
*preprint,model=yes
*node,nset=nbot
1,1.,0.,0.
2,0.707106781,0.707106781,0.
3,0.,1.,0.
4,-0.707106781,0.707106781,0.
5,-1.,0.,0.
6,-0.707106781,-0.707106781,0.
7,0.,-1.,0.
8,0.707106781,-0.707106781,0.
**
*node,nset=nmid
101,1.,0.,0.5
201,0.707106781,0.707106781,0.5
301,0.,1.,0.5
401,-0.707106781,0.707106781,0.5
501,-1.,0.,0.5
601,-0.707106781,-0.707106781,0.5
701,0.,-1.,0.5
801,0.707106781,-0.707106781,0.5
*node,nset=ntop
11,1.,0.,1.
21,0.707106781,0.707106781,1.
31,0.,1.,1.
41,-0.707106781,0.707106781,1.
51,-1.,0.,1.
61,-0.707106781,-0.707106781,1.
71,0.,-1.,1.
81,0.707106781,-0.707106781,1.
**
*nset,nset=nall
nbot,ntop,nmid
*element,type=mcl9,elset=one
1,1,11,31,3,101,21,301,2,201
2,3,31,51,5,301,41,501,4,401
3,5,51,71,7,501,61,701,6,601
4,7,71,11,1,701,81,101,8,801
*membrane section,elset=one,material=mone,ori=or1,
poisson=0.
0.05
*orientation,system=cylindrical,name=or1
0.,0.,0.,0.,0.,1.
1,45.
*material,name=mone
*elastic,type=engineering constants
 1.E2, 1.E8, 1.E2, 0., 0., 0., 1.E2, 1.E2,
 1.E2,
*initial conditions,type=stress
one,0.001,0.001
*transform,type=c,nset=nall
0.,0.,0.,0.,0.,1.
*nset,nset=top
11,21,31,41,51,61,71,81
*node,nset=mstr
1000,0.,0.,1.
*mpc
beam,top,mstr
*step,nlgeom
*static
1.,1.
*boundary
nall,1,3,
1000,1,3,
*output,history,freq=0
*output,field
*element output
s,
e,
sth,
*node output
u,
rf,
*end step
*step,nlgeom
*static
1.,1.
*boundary,op=new
nbot,1,3,0.
*cload
mstr,3,2.
*end step