*heading
	cylindrical membrane
        tested with def-vs-def finite sliding contact
*preprint,model=yes,history=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=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=n_outer
ntop, nbot
*transform, type=c, nset=n_outer
0.0,0.0,0.0, 0.0,0.0,1.0
**
*node, nset=inner_bot
501, 1.2, 0.0, -1.0
517, 1.2, 0.0, -1.0
*node, nset=rigid
1000, 0.0, 0.0, -1.0
**
*ngen, line=c, nset=inner_bot
501, 517, 1, 1000, 0.0, 0.0, -1.0, 0.0, 0.0, 1.0
*ncopy, old set=inner_bot, change number=100, shift, new set=inner_mid
0.0, 0.0, 1.5

*ncopy, old set=inner_bot, change number=200, shift, new set=inner_top
0.0, 0.0, 3.0

*element,type=mcl6,elset=e_outer
1,1,11,31,3,21,2
2,3,31,51,5,41,4
3,5,51,71,7,61,6
4,7,71,11,1,81,8
*element, type=mcl9, elset=e_inner
11, 501, 701, 703, 503, 601, 702, 603, 502, 602
18, 515, 715, 701, 501, 615, 716, 601, 516, 616
*elgen, elset=e_inner
11, 7, 2, 1
*elset, elset=eall
e_outer, e_inner
*membrane section,elset=eall,material=mone,membrane thickness=distThick
 1.
*distribution table, name=distTableThick
 length
*distribution, name=distThick, table=distTableThick, location=element 
  ,1.
 eall, 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
*orientation,name=contactOrient,system=cylindrical
0,0,0,0,0,1
1,90
*contact pair,interaction=rough,TYPE=SURFACE TO SURFACE
asurf,bsurf,contactOrient
*surface interaction,name=rough
*surface behavior, direct
**
**
** history definition
**
**
** Step 1: establish contact
**
*step,nlgeom,unsymm=yes
*static,stabilize
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
*contact file
*controls,parameter=field
1.0e-6,
*print,contact=yes
*output, field, freq=6
*element output, elset=e_outer
s
*CONTACT OUTPUT
 CSTATUS
*node output
u
*end step
**
** step 2: change surface friction
**
*step, nlgeom
*static
1.0,1.0
*change friction, interaction=rough
*friction, anisotropic
0.2, 0.0
*end step
**
** step 3: slide rigid surface to initiate sliding
**
*step, nlgeom, unsymm=yes
*static
0.1, 1.0
*boundary, op=new
1, 2,2,0.0
nbot,3,3,0.0
1000, 1, 2, 0.0
1000, 4, 6, 0.0
1000, 3, 3, 1.0
*output, field, freq=6,VARIABLE=PRESELECT
*element output, elset=e_outer
s
*CONTACT OUTPUT
 CSTATUS
*node output
u
*end step
**
** step 4: twist rigid surface
**
*step, nlgeom, extrapolation=no
*static
0.1, 1.0
*boundary, op=new
1, 2,2,0.0
nbot,3,3,0.0
1000, 1, 2, 0.0
1000, 4, 5, 0.0
1000, 6, 6, 3.1416
1000, 3, 3, 1.0
*output, field, freq=6, VARIABLE=PRESELECT
*element output, elset=e_outer
s
*CONTACT OUTPUT
 CSTATUS
*node output
u
*end step