*HEADING
: SIMPLE CONTACT TEST, FINITE SLIDING CONTACT WITHOUT FRICTION, RIGID
*RESTART,WRITE
** Contact Features Tested:
**	*FRICTION, *CHANGE FRICTION (mu=0 and mu=0.3)
*NODE
*NGEN
*ELEMENT,TYPE=CPE6MH,ELSET=SOLID
*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, 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 1: linear geometry analysis
**
*STEP
*STATIC
*CONTACT PRINT,SLAVE=ASURF
*CONTACT FILE,SLAVE=ASURF
*EL FILE
 S,
*EL PRINT
 S,
*NODE FILE
U,
*NODE PRINT,FREQ=0
*DLOAD
2,P2,250.
*BOUNDARY
1,1
1000,1,2
1000,6
*OUTPUT,FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
*CONTACT OUTPUT,SLAVE=ASURF,VARIABLE=PRESELECT
*END STEP
**
** Step 2: Nonlinear geometry
**
*STEP,NLGEOM
  LARGE DISPLACEMENT ANALYSIS
*STATIC
 .5,1.
*DLOAD
2,P2,250.
*BOUNDARY
1,1
1000,1,2
1000,6
*EL FILE
 S,
*EL PRINT
 S,
*NODE FILE
U,
*OUTPUT,FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
*END STEP
**
** step 3: change surface friction
**
*STEP,NLGEOM
 CHANGE SURFACE FRICTION
*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
7,1,,0.1
8,1,,0.1
9,1,,0.1
*EL FILE
S,
*NODE FILE
U,
*CONTACT FILE,SLAVE=ASURF,MASTER=BSURF
*CONTACT PRINT,SLAVE=ASURF,MASTER=BSURF
*PRINT,CONTACT=YES
*OUTPUT,FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
*CONTACT OUTPUT,SLAVE=ASURF,MASTER=BSURF,VARIABLE=PRESELECT
*ENDSTEP
**
** step 5: Steady state dynamics, direct analysis
**
*STEP
*STEADY STATE DYNAMICS,DIRECT
20,80,5
*DLOAD
2,P2,250.
*BOUNDARY
1,1
1000,1,2
1000,6
*EL PRINT
 S,
*NODE FILE
U,
RF,
*EL FILE
E,
S,
*CONTACT FILE
*OUTPUT,FIELD
*NODE OUTPUT
U,
RF,
*ELEMENT OUTPUT
S,
E,
*CONTACT OUTPUT,SLAVE=ASURF,MASTER=BSURF,VARIABLE=PRESELECT
*END STEP
**
** step 6: frequency extraction ->all modes are extracted
**
*STEP
*FREQUENCY, eigensolver=lanczos
16, , , , ,14
*BOUNDARY
 1,1
 1000,1,2
 1000,6
*MODAL FILE
*CONTACT FILE,FREQ=0
*NODE FILE,F=0
*EL FILE,F=0
*NODE PRINT,F=0
*EL PRINT,F=0
*END STEP
**
** step 7: Steady state dynamics, subspace projection
**
*STEP
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
20,80,5
*DLOAD
2,P2,250.
*NODE FILE
U,
RF,
*EL FILE
E,
S,
*CONTACT FILE
*OUTPUT,FIELD
*NODE OUTPUT
U,
RF,
*ELEMENT OUTPUT
S,
E,
*CONTACT OUTPUT,SLAVE=ASURF,MASTER=BSURF,VARIABLE=PRESELECT
*END STEP