*heading
: contact test, finite sliding of ring on cylinder, deformable
*restart,write,freq=99
*node,nset=nfile
 1, 1.3,7.
 36, 1.3,0.
 301, 1.6,7.
 336, 1.6,0.
 401, 1.6,8.
 416, 1.6,6.
 601, 1.8,8.
 616, 1.8,6.
 701, 1.8,7.
 704, 1.8,6.4
 1001, 2.1,7.
 1004, 2.1,6.4
*ngen,nset=base
 36,336,100
*ngen,nset=top
 1,301,100
*ngen,nset=b1
 401,601,100
*ngen,nset=b2
 416,616,100
*ngen,nset=c1
 701,1001,100
*ngen,nset=c2
 704,1004,100
*nfill
 top,base,35
 b1,b2,15
 c1,c2,3
*element,type=cgax4
 1, 2,102,101,1
 501, 402,502,501,401
 801, 702,802,801,701
*elgen,elset=steel
 1, 3,100,1, 35,1,10
 501, 2,100,1, 15,1,10
 801, 3,100,1, 3,1,10
*elset,elset=outer
 803,813,823
*elset,elset=innerall,generate
 1,341,10
 2,342,10
 3,343,10
*elset,elset=middleall,generate
 501,641,10
 502,642,10
*elset,elset=outerall,generate
 801,803,1
 811,813,1
 821,823,1
*nset,nset=innerall,elset=innerall
*nset,nset=middleall,elset=middleall
*nset,nset=outerall,elset=outerall
*elset,elset=e1,generate
 3,293,10
*elset,elset=e2,generate
 501,641,10
*elset,elset=e3,generate
 502,642,10
*elset,elset=e4,generate
 801,821,10
*surface,name=asurf
 e1,s2
*surface,name=bsurf
 e2,s4
*surface,name=csurf
 e3,s2
*surface,name=dsurf
 e4,s4
*contact pair,interaction=smooth,smooth=0.1
 bsurf,asurf
*contact pair,interaction=rough,smooth=0.1
 dsurf,csurf
*surface interaction,name=smooth
*friction
 .2, 
*surface interaction,name=rough
*friction
 .6, 
*solid section,elset=steel,material=steel
*material,name=steel
*elastic
 30.e6,0.3
** ----------------------------------------------------------------
** 
** STEP: Step-1
** 
*step,nlgeom
  pressurize collar
*static
 1.0,1.0
*boundary
 base,2
 top,2
 336,1
 1001,2
 601,2
 innerall,5,5
 middleall,5,5
 outerall,5,5
*dload
 outer,p2,30000.
*contact print
 ctrq,
*el print,freq=0
 s,e
 sinv, 
*output, field
*element output
 s,
*node output
 u,
*output, history
*contact output
 ctrq,cms
*end step
** ----------------------------------------------------------------
** 
** STEP: Step-2
** 
*step,nlgeom
  move slider down
*static
 .025,1.0
*boundary
 base,2
 top,2
 336,1
 1001,2,,-4.5
 601,2,,-5.
 innerall,5,5
 middleall,5,5
 outerall,5,5
*end step
** ----------------------------------------------------------------
** 
** STEP: Step-3
** 
*step,nlgeom
  twist
*static
 .025,1.0
*boundary,op=new
 base,2
 top,2
 336,1
 1001,2,,-4.5
 601,2,,-5.
 innerall,5,5
 outerall,5,5,0.1
*end step