*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, TYPE=SURFACE TO SURFACE
SSLIDE, SBASE
*SURFACE INTERACTION, NAME=FRIC
***
*** PARAMETERS: mu_i  = 0.1
***             mu_f  = 0.2
***             alpha =-0.2
***             beta  = 0.0
***
***        SDV(1) = Current friction coefficient for visualization purposes only
***        SDV(2) = Solution dependent variable that depends on equivalent contact slip
***
*FRICTION,USER=COEFFICIENT,PROPERTIES=4,DEPVAR=2
0.1,0.2,-0.2,0.0
*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
***
***	ESTABLISH CONTACT
***
*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
*OUTPUT,FIELD,FREQ=1
*ELEMENT OUTPUT, VARIABLE=PRESELECT
*NODE OUTPUT, VARIABLE=PRESELECT
*CONTACT OUTPUT
 CSTRESS,CDISP,SDV1
*END STEP
***
***     ACCUMULATE SLIP AND MODIFY FRICTION COEFFICIENT
***
*STEP,INC=30,UNSYMM=YES
SLIDE TO THE RIGHT
*STATIC
.1,1.,,0.2
*BOUNDARY
TOP,2,2,-.15
TOP,1,1,4.
*END STEP
***
***   DUMMY STEP TO TRANSFER STATE VARIABLES
***
*STEP,INC=30,UNSYMM=YES
SLIDE TO THE RIGHT
*STATIC
1.,1.
*END STEP
**
**
** STEP 4: 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
*OUTPUT,FIELD, FREQ=1
*NODE OUTPUT
U
*ELEMENT OUTPUT
S
*CONTACT OUTPUT
CSTRESS,CDISP,SDV1
*END STEP
**
** STEP 5: frequency extraction 
**
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS, SIM=NO
15,
*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
*OUTPUT,FIELD, FREQ=1
*NODE OUTPUT
U
*ELEMENT OUTPUT
S
*CONTACT OUTPUT
CSTRESS,CDISP,SDV1
*END STEP
**
** SSD SUBSPACE 6: 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
*OUTPUT,FIELD, FREQ=1
*NODE OUTPUT
U
*ELEMENT OUTPUT
S
*CONTACT OUTPUT
CSTRESS,CDISP,SDV1
*END STEP