*HEADING
: SIMPLE CONTACT TEST, FINITE SLIDING CONTACT WITH ZERO FRICTION, DEFORMABLE
*AMPLITUDE,NAME=CONST
0,1,1,1
*RESTART,WRITE
*NODE
 1,0.0,0.0
 7,7.5,0.0
 11,,2.
 17,7.5,2.0
101,0.5,2.1
 103,2.25,2.1
 111,0.5,4.1
 113,2.25,4.1
 1001,,1.
 1007,7.5,1.
 1101,0.5,3.1
 1103,2.25,3.1
1102,1.375,3.1
*NGEN,NSET=BASE
 1,7
*NGEN,NSET=TOP
 111,113
*NGEN
 11,17
 101,103
 1001,1007
 1101,1103
*ELEMENT,TYPE=CPS6M,ELSET=STEEL
 1,1,3,11,2,1002,1001
 2,3,13,11,1003,12,1002
 3,3,5,13,4,1004,1003
 4,5,15,13,1005,14,1004
 5,5,7,15,6,1006,1005
 6,7,17,15,1007,16,1006
 11,101,103,111,102,1102,1101
 12,103,113,111,1103,112,1102
*ELSET,ELSET=ETOP
11,12
*ELSET,ELSET=EBOT
1,2,3,4,5,6
*surface, NAME=ASURF
11,S1
*surface, NAME=BSURF
2,S2
4,S2
6,S2
*CONTACT PAIR,INTERACTION=ROUGH,SMOOTH=0.1
ASURF,BSURF
*SURFACE INTERACTION, NAME=ROUGH
*SOLID SECTION,ELSET=STEEL,MATERIAL=STEEL
*MATERIAL,NAME=STEEL
*ELASTIC
  30.E6,0.3
*DENSITY
 10.0, 
******************************
*STEP,NLGEOM
CONTACT
*STATIC
0.67, 1.0
*CONTROLS,PARAMETERS=FIELD
1.0E-6, 
*BOUNDARY
 BASE,2
 1,1
 11,1
 1001,1
 TOP,2,,-.15
 TOP,1
*PRINT,CONTACT=YES
*CONTACT PRINT,SLAVE=ASURF
*EL PRINT
S,E
*NODE PRINT
U,RF
*END STEP
*STEP,NLGEOM
 MOVE TOP ACROSS
*STATIC
  1.,10.,1.
*BOUNDARY,OP=NEW
 BASE,2
 1,1
 11,1
 1001,1
 TOP,2,,-.15
 TOP,1,,4.
*END STEP
**
** step 3: Steady state dynamics, direct analysis
**
*STEP
*STEADY STATE DYNAMICS,DIRECT
20,100,5
*BOUNDARY
 BASE,2
 1,1,1
 11,1,1
 1001,1,1
 TOP,1,1
 TOP,2,2,0.1
*END STEP
**
** step 4: frequency extraction 
**
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS
54,,,-1.,,15
*BOUNDARY
 BASE,2
 1,1,1
 11,1,1
 1001,1,1
 TOP,1,1
*BOUNDARY,BASE NAME=MOVE
 TOP,2,2
*NODE PRINT,F=0
*EL PRINT,F=0
*END STEP
**
** step 5: Steady state dynamics, subspace projection
**
*STEP
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
20,100,5
*BASE MOTION,BASE NAME=MOVE,AMPLITUDE=CONST,DOF=2,TYPE=DISPLACEMENT,SCALE=0.1
*OUTPUT,FIELD,OP=NEW
*ELEMENT OUTPUT
*NODE OUTPUT
U,
*END STEP