*heading
cylindrical membrane second-order elements
tested with def-vs-rigid small sliding contact
*preprint,model=yes,history=yes
*node,nset=nbot
**
*node,nset=nmid
**
*node,nset=ntop
*node
**
*element,type=mcl9,elset=eall
*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=bsurf,type=revolution
start, 1.2,  3.0
line,  1.2, -2.0
*rigid body,analytical surface=bsurf,ref node=1000
*contact pair,interaction=rough, small sliding, TYPE=SURFACE TO SURFACE
asurf,bsurf
*surface interaction,name=rough
**
** history definition
**
** Step 1: establish contact
**
*step,nlgeom
*static
0.1,1.
*boundary
1, 2,2,0.0
nbot,3,3,0.0
1000,1,6,0.0
*contact interference,shrink
asurf, bsurf
*el print
s,
e,
*contact print
*controls,parameter=field
1.0e-6,
*print,contact=yes
*output, field, freq=6
*element output, elset=eall
s
*node output
u
*contact output, slave=asurf, variable=preselect
*end step
**
** step 2: slide w/o friction, linear perturbation
**
*step, nlgeom, pert
*static
*boundary, op=mod
1000, 3, 3, 0.01
*contact print
*output, field
*element output, elset=eall
s
*node output
u
*end step
**
** step 3: change surface friction
**
*step, nlgeom
*static
1.0,1.0
*change friction, interaction=rough
*friction
0.3,
*output, field, freq=6
*element output, elset=eall
s
*node output
u
*end step
**
** step 4: slide with friction
**
*step, nlgeom
*static
1.0,1.0
*boundary, op=mod
1000, 3, 3, 0.01
*output, field, freq=6
*element output, elset=eall
s
*node output
u
*contact output, slave=asurf, variable=preselect
*end step