*HEADING
: SIMPLE CONTACT TEST, RIGID
*RESTART,WRITE
*NODE,NSET=NALL
1,0.
2,5.
3,5.,5.
4,0.,5.
1000,0.
*ELEMENT,TYPE=CPE4,ELSET=SOLID
1,1,2,3,4
*rigid body,analytical surface=BSURF,REFNODE=1000
*surface,TYPE=SEGMENT,NAME=BSURF
START,-100.
LINE,  100.
*surface,type=node,NAME=ASURF
1,2.5
2,2.5
*CONTACT PAIR,INTERACTION=SMOOTH
ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
*SOLID SECTION,ELSET=SOLID,MATERIAL=SOLMAT
*MATERIAL,NAME=SOLMAT
*ELASTIC
30.E6,.3
*DENSITY
10., 
*BOUNDARY
1,1
4,1
1000,1,2
1000,6
*STEP
*STATIC
1.,1.
*CONTACT PRINT,SECONDARY=ASURF
*EL PRINT
 S,
*NODE PRINT,FREQ=0
*DLOAD
1,P3,5000.
*END STEP
*STEP,NLGEOM
*STATIC
 .5,1.
*DLOAD
1,P3,5000.
*END STEP
**
** step 5: Steady state dynamics, direct analysis
**
*STEP
*STEADY STATE DYNAMICS,DIRECT
20,80,5
*BOUNDARY
1,1
4,1
1000,1,2
1000,6
*DLOAD
1,P3,5000.
*OUTPUT,FIELD
*NODE OUTPUT
U,
RF,
*ELEMENT OUTPUT
E,
S,
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASURF
*END STEP
**
** step 6: frequency extraction
**
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS, SIM=NO
4,,,,
*BOUNDARY
1,1
4,1
1000,1,2
1000,6
*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
 1,P3,5000.
*OUTPUT,FIELD
*NODE OUTPUT
U,
RF,
*ELEMENT OUTPUT
S,
E,
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASURF
*END STEP