*heading
cylindrical second-order membrane
tested with def-vs-rigid finite sliding contact
*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
102,0.707106781,0.707106781,0.5
103,0.,1.,0.5
104,-0.707106781,0.707106781,0.5
105,-1.,0.,0.5
106,-0.707106781,-0.707106781,0.5
107,0.,-1.,0.5
108,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.
*node
1000, 0,0,0
**
*element,type=mcl9,elset=eall
1,1,11,31,3,101,21,103,2,102
2,3,31,51,5,103,41,105,4,104
3,5,51,71,7,105,61,107,6,106
4,7,71,11,1,107,81,101,8,108
*membrane section,elset=eall,material=mone
0.05
*material,name=mone
*elastic
30.e6,0.3
**
** contact definition
**
*surface,name=asurf
eall,spos
*surface,name=rsurf,type=revolution
0.,0.,0.,0.,0.,1.
start, 1.2,  3.0
line,  1.2, -2.0
*rigid body,analytical surface=rsurf,ref node=1000
*contact pair,interaction=rough,type=surface
asurf,rsurf
*surface interaction,name=rough
*surface behavior,penalty
,,1.e3
**
*nset,nset=ncyl
ntop, nmid, nbot
*transform,nset=ncyl,type=c
0,0,0, 0,0,1.0
*equation
2
nmid,1,1.0, nbot,1,-1.0
**
** history definition
**
**
** Step 1: establish contact
**
*step,unsymm=yes,CONVERT SDI=NO
*static
0.1,1.
*boundary
1, 2,2,0.0
nbot,3,3,0.0
1000, 1, 2, 0.0
1000, 3, 3, 0.0
1000, 4, 6, 0.0
*contact interference,shrink
asurf, rsurf
*controls,parameter=field
1.0e-6,
*print,contact=yes
**
*el print
s,
e,
*contact print
*output, field, freq=6
*element output, elset=eall
s
*node output
u
*contact output, slave=asurf, variable=preselect
*end step
**
** step 2: change surface friction
**
*step, nlgeom,CONVERT SDI=NO
*static
1.0,1.0
*change friction, interaction=rough
*friction
0.3,
*end step
**
** step 3: slide rigid surface to initiate sliding
**
*step, nlgeom, unsymm=yes,CONVERT SDI=NO
*static
0.1, 1.0
*boundary
1000, 3, 3, 1.0
*end step