*HEADING
: SIMPLE CONTACT TEST, 3D FINITE SLIDING, FRICTION
*AMPLITUDE,NAME=CONSTA
0,0.048,1,0.048
*AMPLITUDE,NAME=CONSTB
0,0.060,1,0.060
*AMPLITUDE,NAME=CONSTC
0,-0.064,1,-0.064
*ELEMENT,TYPE=C3D8,ELSET=BASE
1, 1,3,13,11,101,103,113,111
2, 3,5,15,13,103,105,115,113
3, 5,7,17,15,105,107,117,115
*ELEMENT,TYPE=C3D8,ELSET=SLIDER
4, 21,23,33,31,121,123,133,131
*ELSET,ELSET=SOLID
BASE,SLIDER
*NODE,NSET=NALL
1,0.0,0.0
7,7.5,0.0
11,0.0,2.0
17,7.5,2.0
21,0.5,2.1,.25
23,2.25,2.1,.25
31,0.5,4.1,.25
33,2.25,4.1,.25
101,0.0,0.0,2.5
107,7.5,0.0,2.5
111,0.0,2.0,2.5
117,7.5,2.0,2.5
121,0.5,2.1,2.25
123,2.25,2.1,2.25
131,0.5,4.1,2.25
133,2.25,4.1,2.25
*NGEN,NSET=NALL
1,7,2
11,17,2
21,23,2
31,33,2
101,107,2
111,117,2
121,123,2
131,133,2
*NSET,NSET=BOTTOM,G
1,7,2
101,107,2
*NSET,NSET=TOP,G
31,33,2
131,133,2
*NMAP,NSET=NALL,TYPE=RECTANGULAR
0,0,0,.8,.0,.6
.48,.6,-.64
*TRANSFORM,NSET=NALL
.8,.0,.6,.48,.6,-.64
*SOLID SECTION, ELSET=SOLID,MATERIAL=STEEL
*surface, NAME=SBASE
BASE, S5
*surface, NAME=SSLIDE
SLIDER, S3
*CONTACT PAIR, INTERACTION=FRIC
SSLIDE, SBASE
*SURFACE INTERACTION, NAME=FRIC
*FRICTION
0.1, 
*MATERIAL,NAME=STEEL
*ELASTIC
 30.E6,0.3
*DENSITY
 10.0, 
*BOUNDARY
1,1
11,1
101,1
111,1
BOTTOM,2,3
TOP,3
*STEP,NLGEOM,INC=20,UNSYMM=YES
 FIRST COMPRESS
*STATIC,DIRECT
0.67,1.
*CONTROLS,PARAMETERS=FIELD
1.0E-6, 
*BOUNDARY
TOP,1
TOP,2,2,-.15
*PRINT,CONTACT=YES
*CONTACT PRINT
*CONTACT FILE
*EL PRINT
S, 
E, 
*NODE PRINT
U,
RF,
*EL FILE,FREQ=5
S, 
E, 
*NODE FILE,FREQ=5
U,RF
*END STEP
*STEP,INC=30,UNSYMM=YES
SLIDE TO THE RIGHT
*STATIC
.2,1.
*CHANGE FRICTION,INTERACTION=FRIC
*FRICTION
0.2, 
*BOUNDARY
TOP,2,2,-.15
TOP,1,1,4.
*END STEP
**
** step 3: Steady state dynamics, direct analysis
**
*STEP
*STEADY STATE DYNAMICS,DIRECT
20,100,5
*BOUNDARY
TOP,1,1 
TOP,2,2,0.1
TOP,3,3
BOTTOM,1,3
*NODE FILE
U, 
RF, 
*EL FILE
E, 
S, 
*OUTPUT,FIELD
*node output
u,
*CONTACT FILE
*END STEP
**
** step 4: frequency extraction 
**
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS
15,,,-1.0,
*BOUNDARY
 BOTTOM,1,3
*BOUNDARY,BASE NAME=MOVA
 TOP,1,1
*BOUNDARY,BASE NAME=MOVB
 TOP,2,2
*BOUNDARY,BASE NAME=MOVC
 TOP,3,3
*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 5: Steady state dynamics, subspace projection
**
*STEP
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
20,100,5
*select eigenmodes,generate
1,15,1
*BASE MOTION,BASE NAME=MOVA,AMPLITUDE=CONSTA,DOF=1,TYPE=DISPLACEMENT
*BASE MOTION,BASE NAME=MOVB,AMPLITUDE=CONSTB,DOF=2,TYPE=DISPLACEMENT
*BASE MOTION,BASE NAME=MOVC,AMPLITUDE=CONSTC,DOF=3,TYPE=DISPLACEMENT
*NODE FILE
TU, 
RF, 
*EL FILE
E, 
S, 
*CONTACT FILE
*END STEP