*heading
	Cylinder - 12-node full integration elements 
        tested with def-vs-def finite sliding contact
*preprint, model=yes, history=yes
*node, nset=bot
1,  1.0, 0.0
17, 1.0, 0.0
21, 1.5, 0.0
37, 1.5, 0.0
41, 2.0, 0.0
57, 2.0, 0.0
*node
999, 0.0, 0.0
*ngen, line=c, nset=bot
1, 17, 1, 999, , , , 0.0, 0.0, 1.0
21, 37, 1, 999, , , , 0.0, 0.0, 1.0
41, 57, 1, 999, , , , 0.0, 0.0, 1.0
*ncopy, old set=bot, change number=200, shift, new set=top
0.0, 0.0, 1.0

*nset, nset=n_outer
bot,top
*element, type=ccl12, elset=e_outer
1,  3,43,243,203,  1,41,241,201,  2,42,242,202
8,  1,41,241,201, 15,55,255,215, 16,56,256,216
*elgen, elset=e_outer
1, 7, 2, 1
**
**  definition of inner cylinder
**
*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=200, shift, new set=inner_top
0.0, 0.0, 3.0

*element, type=mcl6, elset=e_inner
11, 501, 701, 703, 503, 702, 502
18, 515, 715, 701, 501, 716, 516
*elgen, elset=e_inner
11, 7, 2, 1
*elset, elset=eall
e_outer, e_inner
**
*solid section, elset=e_outer, material=st1
*membrane section, elset=e_inner, material=st1
0.05
*material, name=st1
*elastic
30E6, 0.3
**
** contact definition
**
*surface,name=asurf
e_outer, s6
*surface,name=bsurf
e_inner, sneg
*rigid body, elset=e_inner, ref node=rigid
*contact pair,interaction=rough,TYPE=SURFACE TO SURFACE
asurf,bsurf
*surface interaction,name=rough
**
** node set definition
**
*nset, nset=in, gen
1,16
201,216
*nset, nset=out, gen
41,56
241,256
*nset, nset=out_bot, gen
41,56
*transform, type=c, nset=n_outer
0.0,0.0,0.0, 0.0,0.0,1.0
**
** HISTORY DEFINITION
**
** step 1: Establish contact
**
*step, nlgeom, unsymm=yes
 establish contact
*static
0.1, 1.0
*contact interference,shrink
asurf, bsurf
**
*** Boundary conditions
**
*boundary
out_bot, 3, 3
41, 2, 2
rigid, 1, 2
rigid, 3, 3
rigid, 4, 5
rigid, 6, 6
*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
*node output
u
*end step
**
** step 2: change surface friction
**
*step, nlgeom, unsymm=yes
*static
1.0,1.0
*change friction, interaction=rough
*friction, anisotropic
0.2, 0.0
*output, field, freq=6
*element output, elset=e_outer
s
*node output
u
*end step
**
** step 3: slide with friction
**
*step, nlgeom, unsymm=yes
*static
0.1, 1.0
*boundary, op=mod
rigid, 3, 3, 1.0
*output, field, freq=6
*element output, elset=e_outer
s
*node output
u
*end step
**
** step 4: rotate the rigid surface to initiate sliding
**
*step, nlgeom, extrapolation=no
*static
0.1, 1.0
*boundary, fixed
out, 2, 2
*boundary, op=mod
rigid, 6, 6, 3.1416
*end step