*HEADING
: SIMPLE CONTACT TEST, FINITE SLIDING CONTACT WITH FRICTION, DEFORMABLE
*RESTART,WRITE,FREQ=999
*AMPLITUDE,NAME=CONST
0,1,1,1
*AMPLITUDE,NAME=ZERO
0,0,1,0
*NODE
*NGEN,NSET=BASE
*NGEN,NSET=TOP
*NGEN
*ELEMENT,TYPE=CPE8,ELSET=STEEL
*ELGEN,ELSET=STEEL
*ELSET,ELSET=ETOP
*ELSET,ELSET=EBOT
*surface, NAME=ASURF
ETOP,S1
*surface, NAME=BSURF
EBOT,S3
*CONTACT PAIR, INTERACTION=ROUGH,SMOOTH=0.1
ASURF,BSURF
*SURFACE INTERACTION, NAME=ROUGH
*FRICTION
0.2,
*SOLID SECTION,ELSET=STEEL,MATERIAL=STEEL
*MATERIAL,NAME=STEEL
*ELASTIC
  30.E6,0.3
*DENSITY
  7800.00,
******************************
*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
*CONTACT FILE, SLAVE=ASURF,F=5
*EL PRINT
S,E
*NODE PRINT
U,RF
*EL FILE,FREQ=5
S,E
*NODE FILE,FREQ=5
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
Frequency Extraction
*FREQUENCY,EIGENSOLVER=LANCZOS, SIM=NO
5,
*CONTACT FILE, SLAVE=ASURF
*EL FILE
S,E
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*ENDSTEP
**
** step 4: Steady state dynamics, direct analysis
**
*STEP
*STEADY STATE DYNAMICS,DIRECT
20,100,5
*BOUNDARY
 BASE,2
 1,1
 11,1
 1001,1
 TOP,2,,0.01
 TOP,1
*NODE FILE
U,
RF,
*EL FILE
E,
S,
*CONTACT FILE
*END STEP
**
** step 5: frequency extraction
**
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS, SIM=NO
30,,,-1.
*BOUNDARY
 1,1
 11,1
 1001,1
 TOP,1
 TOP,2,2
*BOUNDARY,BASE NAME=ZERO
 BASE,2
*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 6: Steady state dynamics, subspace projection
**
*STEP
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
20,100,5
*SELECT EIGENMODES,GENERATE
*BASE MOTION,BASE NAME=ZERO, AMPLITUDE=ZERO ,DOF=2,TYPE=DISPLACEMENT
*BASE MOTION,AMPLITUDE=CONST,DOF=2,TYPE=DISPLACEMENT,SCALE=0.01
*NODE FILE
TU,
RF,
*EL FILE
E,
S,
*CONTACT FILE
*END STEP