*HEADING
: SIMPLE CONTACT TEST, FINITE SLIDING CONTACT WITH ZERO FRICTION, DEFORMABLE
*AMPLITUDE,NAME=CONST
0,1,1,1
*AMPLITUDE,NAME=ZERO
0,0,1,0
*RESTART,WRITE
*NODE
1,0.0,0.0
7,0.0,7.5
11,2.,0.0
17,2.0,7.5
101,2.1,7.0
103,2.1,5.25
111,4.1,7.0
113,4.1,5.25
1001,1.0,0.0
1007,1.0,7.5
1101,3.1,7.0
1103,3.1,5.25
*NGEN,NSET=INNER
1,7
*NGEN,NSET=OUTER
111,113
*NGEN
11,17
101,103
1001,1007,2
1101,1103
*ELEMENT,TYPE=CAX8,ELSET=STEEL
1,1,11,13,3,1001,12,1003,2
11,101,103,113,111,102,1103,112,1101
*ELGEN,ELSET=STEEL
1,3,2
*ELSET,ELSET=ETOP
11, 
*ELSET,ELSET=EBOT
3,2,1
*surface,NAME=ASURF
ETOP,S1
*surface,NAME=BSURF
EBOT,S2
*CONTACT PAIR,INTERACTION=SMOOTH,TYPE=SURFACE
ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
*SURFACE BEHAVIOR, DIRECT
*SOLID SECTION,ELSET=STEEL,MATERIAL=STEEL
*MATERIAL,NAME=STEEL
*ELASTIC
 30.E6,0.3
*DENSITY 
10.0, 
**********
*STEP,NLGEOM,UNSYMM=YES
  CONTACT
*STATIC
0.7, 1.0
*CONTROLS,PARAMETERS=FIELD
1.0E-6, 
*BOUNDARY
 1,2
 11,2
 1001,2
 INNER,1
 OUTER,1,1,-.15
 OUTER,2
*PRINT,CONTACT=YES
*CONTACT PRINT,SLAVE=ASURF
*EL PRINT
S,E
*NODE PRINT
U,RF
*END STEP
*STEP,NLGEOM
  MOVE OUTER ACROSS
*STATIC
 1.,10.,1.
*BOUNDARY,OP=NEW
 INNER,1
1,2
11,2
1001,2
OUTER,1,1,-.15
OUTER,2,2,-4.
*END STEP
**
** step 3: Steady state dynamics, direct analysis
**
*STEP
*STEADY STATE DYNAMICS,DIRECT
20,100,5
*BOUNDARY,OP=NEW
 INNER,1
 1,2
 7,2
 11,2
 17,2
 OUTER,1,1
 OUTER,2,2,0.01
*END STEP
**
** step 4: frequency extraction 
**
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS
32,
*BOUNDARY,OP=NEW
 INNER,1
 OUTER,1,1
 OUTER,2,2
*BOUNDARY,BASE NAME=ZERO
 1,2
 7,2
 11,2
 17,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=ZERO,AMPLITUDE=ZERO,DOF=2,TYPE=DISPLACEMENT
*BASE MOTION,AMPLITUDE=CONST,DOF=2,TYPE=DISPLACEMENT,SCALE=0.01
*END STEP