*HEADING
SIMPLE  TESTS, FINITE SLIDING, [CGAX8, RAX2]             :--
**RESTART,WRITE, FREQ=30
*NODE
*ELEMENT,TYPE=CGAX8,ELSET=LEFTEL
*ELEMENT,TYPE=RAX2,ELSET=RIGID
*NSET, NSET=LT_N
*SOLID SECTION,ELSET=LEFTEL,MATERIAL=SOLMAT
*RIGID BODY, ELSET=RIGID, REFNODE=100
*MATERIAL,NAME=SOLMAT
*ELASTIC
30.E6,0.30
*SURFACE INTERACTION,NAME=SLIDE
 1.,
*FRICTION, ANISOTROPIC
0.1,0.2
*surface,NAME=ASURF
LEFTEL,S2
*surface,NAME=RSURF
RIGID, SPOS
*CONTACT PAIR,INTERACTION=SLIDE,TYPE=SURFACE
ASURF,RSURF
*BOUNDARY
LT_N,1,2
LT_N,5,5
100,1,6
**
**    step 1  :  interference fit, (radial intereference)
**
*STEP, INC=20, NLGEOM, unsymm=yes
           SHRINK FIT
*STATIC
0.1,1.,0.1,0.2
*CONTACT INTERFERENCE,SHRINK
ASURF, RSURF
*NODE PRINT, F=0
*EL PRINT, F=0
*OUTPUT,FIELD,FREQ=0
*ELEMENT OUTPUT
*OUTPUT,FIELD,FREQ=0
*NODE OUTPUT
*OUTPUT,FIELD,FREQ=50
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF,MASTER=RSURF
*CONTACT PRINT, SLAVE=ASURF , MASTER=RSURF, F=50
*ENDSTEP
**
**     step 2 : linear perturbation (about shrink fit); sliding upwards
**
*STEP,PERT,  NLGEOM, unsymm=yes
 PERT :- SLIDING UPWARDS
*STATIC
*BOUNDARY, OP=MOD
LT_N,2,,0.1
*NODE PRINT, F=0
*EL PRINT, F=0
*OUTPUT,FIELD,FREQ=0
*ELEMENT OUTPUT
*OUTPUT,FIELD,FREQ=0
*NODE OUTPUT
*OUTPUT,FIELD,FREQ=50
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF,MASTER=RSURF
*CONTACT PRINT, SLAVE=ASURF , MASTER=RSURF, F=50
*ENDSTEP
**
**     step 3 : linear perturbation (about shrink fit); twist
**
*STEP, PERT, NLGEOM, unsymm=yes
  PERT: TWISTING
*STATIC
*BOUNDARY, OP=MOD
LT_N,5,,0.1
*OUTPUT,FIELD,FREQ=50
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF,MASTER=RSURF
*CONTACT PRINT, SLAVE=ASURF , MASTER=RSURF, F=50
*ENDSTEP
**
**     step 4: UNSYMM, linear perturbation, slide and twist
** 	          with shrink-fit base state.
**
*STEP,  PERT, NLGEOM, UNSYMM=YES
*STATIC
*BOUNDARY
LT_N,5,,0.1
LT_N,2,,0.1
*OUTPUT,FIELD,FREQ=50
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF,MASTER=RSURF
*CONTACT PRINT, SLAVE=ASURF , MASTER=RSURF, F=50
*ENDSTEP
**
**       step 5: NLGEOM, UNSYMM, twist and sliding
**
*STEP, NLGEOM, UNSYMM=YES
*STATIC
0.1,1.0
*CONTROLS,PARAMETERS=FIELD
1.0E-6,
*BOUNDARY
LT_N,2,,0.1
LT_N,5,,0.1
*ENDSTEP