*heading
	cylindrical membrane
        tested with def-vs-def small sliding contact
*preprint,model=yes,history=yes
*node,nset=nbot
**
*node,nset=ntop
**
*node, nset=inner_bot
*node, nset=rigid
*nset, nset=n21
**
*ngen, line=c, nset=inner_bot
*ncopy, old set=inner_bot, change number=100, shift, new set=inner_mid
*ncopy, old set=inner_bot, change number=200, shift, new set=inner_top
*element,type=mcl6,elset=e_outer
*element, type=mcl9, elset=e_inner
*elgen, elset=e_inner
*elset, elset=eall
e_outer, e_inner
*membrane section,elset=eall,material=mone
0.05
*material,name=mone
*elastic
30.e6,0.3
**
** contact definition
**
*surface,name=asurf
e_outer,spos
*surface,name=bsurf
e_inner, sneg
*rigid body, elset=e_inner, ref node=rigid
*contact pair,interaction=rough,small sliding,smooth=0.0
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, 2, 0.0
1000, 3, 3, 0.0
1000, 4, 6, 0.0
*contact interference,shrink
asurf, bsurf
*el print
s,
e,
*contact print
*controls,parameter=field
1.0e-6,
*print,contact=yes
*output, history
*node output,nset=n21
u1,u2,u3
*contact output, slave=asurf
cpress,
cshear1,
*end step
**
** step 2: slide surface
**
*step, nlgeom, pert
*static
*boundary, op=mod
1000, 3, 3, 0.1
*contact print
*output, history
*node output,nset=n21
u1,u2,u3
*end step
**
** step 3: change surface friction
**
*step, nlgeom
*static
1.0,1.0
*change friction, interaction=rough
*friction
0.3,
*output, history
*node output,nset=n21
u1,u2,u3
*contact output, slave=asurf
cpress,
*end step
**
** step 4: slide rigid surface to initiate sliding
**
*step, nlgeom
*static
1.0, 1.0
*boundary, op=mod
1000, 3, 3, 0.1
*output, history
*node output,nset=n21
u1,u2,u3
*contact output, slave=asurf
cpress,
cshear1,
*end step