*heading
	Cylinder - 24-node full integration elements
        tested with def-vs-def finite sliding contact
*preprint, model=yes, history=yes
*node, nset=bot
*node
*ngen, line=c, nset=bot
*ncopy, old set=bot, change number=100, shift, new set=mid
*ncopy, old set=bot, change number=200, shift, new set=top
*nset, nset=n_outer
bot,mid,top
*element, type=ccl24, elset=e_outer
*elgen, elset=e_outer
**
**  definition of inner cylinder
**
*node, nset=inner_bot
*node, nset=rigid
**
*ngen, line=c, nset=inner_bot
*ncopy, old set=inner_bot, change number=200, shift, new set=inner_top
*element, type=mcl6, elset=e_inner
*elgen, elset=e_inner
*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=out, gen
*nset, nset=out_bot, gen
*transform, type=c, nset=n_outer
**
** 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.05, 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, direct
0.1, 1.0
*boundary, fixed
out, 2, 2
*boundary, op=mod
rigid, 6, 6, 3.1416
*end step