*HEADING
: SIMPLE CONTACT TEST, FINITE SLIDING CONTACT WITHOUT FRICTION, RIGID 
** Contact Features Tested: 
**	*FRICTION, *CHANGE FRICTION (mu=0 and mu=0.3)
*NODE
1, 0.
2, 5.
3, 0., 5.
4, 5., 5.
1000,0.
*ELEMENT,TYPE=CPE4,ELSET=SOLID
 1, 1,2,4,3
*SOLID SECTION,ELSET=SOLID,MATERIAL=MAT1
 .25, 
*MATERIAL,NAME=MAT1
*ELASTIC
 30.E6,.3
*DENSITY
 10.0, 
*rigid body,analytical surface=BSURF,REFNODE=1000
*surface,TYPE=SEGMENTS,NAME=BSURF
 START,0.,0.
 LINE,5.5,0.
*surface,NAME=ASURF
1,S1
*CONTACT PAIR,INTERACTION=SMOOTH
ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
.25, 
*FRICTION
0.0, 
*STEP
*STATIC
*CONTACT PRINT,SECONDARY=ASURF
*CONTACT FILE,SECONDARY=ASURF
*EL FILE
 S,
*NODE FILE
U, 
*EL PRINT
 S,
*NODE PRINT,FREQ=0
*DLOAD
 1,P3,250.
*BOUNDARY
 1,1
 1000,1,2
 1000,6
*OUTPUT,FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASURF
*END STEP
*STEP,NLGEOM
  LARGE DISPLACEMENT ANALYSIS
*STATIC
 .5,1.
*DLOAD
 1,P3,250.
*BOUNDARY
 1,1
 1000,1,2
 1000,6
*CONTACT FILE,SECONDARY=ASURF
*EL FILE
 S,
*NODE FILE
U, 
*OUTPUT,FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASURF
*END STEP
**
** step 3: change surface friction
**
*STEP,NLGEOM
*STATIC
1.0,1.0
*CHANGE FRICTION,INTERACTION=SMOOTH
*FRICTION
0.3, 
*ENDSTEP
**
** step 4: slide to the left w/ friction
**
*STEP,NLGEOM
*STATIC
*BOUNDARY
3,1,,0.1
4,1,,0.1
*EL FILE
S, 
*NODE FILE
U, 
*CONTACT FILE,SECONDARY=ASURF,MAIN=BSURF
*CONTACT PRINT,SECONDARY=ASURF,MAIN=BSURF
*PRINT,CONTACT=YES
*OUTPUT,FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASURF
*ENDSTEP
**
** step 5: Steady state dynamics, direct analysis
**
*STEP
*STEADY STATE DYNAMICS,DIRECT
20,80,5
*DLOAD
 1,P3,100.
*BOUNDARY
 1,1
 1000,1,2
 1000,6
*NODE FILE
U, 
RF, 
*EL FILE
E, 
S, 
*CONTACT FILE
*OUTPUT,FIELD
*NODE OUTPUT
U,
RF,
*ELEMENT OUTPUT
E,
S,
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASURF
*END STEP
**
** step 6: frequency extraction ->all modes are extracted
**
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS, SIM=NO
3,,,,10
*BOUNDARY
 1,1
 1000,1,2
 1000,6
*CONTACT FILE,FREQ=0
*NODE FILE,F=0
*EL FILE,F=0
*NODE PRINT,F=0
*EL PRINT,F=0
*MODAL FILE
*END STEP
**
** step 7: Steady state dynamics, subspace projection
**
*STEP
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
20,80,5
*DLOAD
 1,P3,100.
*NODE FILE
U, 
RF, 
*EL FILE
E, 
S, 
*CONTACT FILE
*OUTPUT,FIELD
*NODE OUTPUT
U,
RF,
*ELEMENT OUTPUT
S,
E,
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASURF
*END STEP