*heading
	Cylinder - 18-node full integration elements
        tested with def-vs-def small 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=nall
bot,mid,top
*element, type=ccl18, elset=e_outer
**
*elgen, elset=e_outer
*elset, elset=asurf, gen
**
**  definition of inner cylinder
**
*node, nset=inner_bot
*node, nset=rigid
**
*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=mcl9, 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
asurf, s5
*surface,name=bsurf
e_inner, sneg
*rigid body, elset=e_inner, ref node=rigid
*contact pair,interaction=rough,small sliding
asurf,bsurf
*surface interaction,name=rough
**
** node set definition
**
*nset, nset=in, gen
*nset, nset=out, gen
*nset, nset=out_bot, gen
*transform, type=c, nset=nall
**
** HISTORY DEFINITION
**
** step 1: Shrink fit with *CONTACT INTERFERENCE.
** 	   Establish base state.
**
*step, nlgeom
 establish contact
*static
0.1, 1.0
*contact interference,shrink
asurf, bsurf
**
*** Boundary conditions
**
*boundary
out_bot, 3, 3
41, 2, 2
rigid, 1, 6
*el print
s,
e,
*contact print
*controls,parameter=field
1.0e-6,
*print,contact=yes
*output, field, freq=6
*element output, elset=e_outer
s
*node output
u
*contact output, slave=asurf, variable=preselect
*end step
**
** step 2: slide w/o friction, linear perturbation
**
*step, pert, nlgeom
*static
*boundary, op=mod
rigid, 3, 3, 0.01
*contact print
*output, field
*element output, elset=e_outer
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=e_outer
s
*node output
u
*end step
**
** step 4: slide with friction
**
*step, nlgeom
*static
1.0, 1.0
*boundary, op=mod
rigid, 3, 3, 0.01
*output, field, freq=6
*element output, elset=e_outer
s
*node output
u
*contact output, slave=asurf, variable=preselect
*end step