*HEADING
	2-D ELEMENT IN CONTACT WITH ANOTHER 2-D ELEMENT.
	NORMALS ARE NOT MODIFIED. 
** Contact Features Tested:
**	*CLEARANCE
** Dynamic Features Tested:
**      *STEADY STATE DYNAMICS,DIRECT
**      *FREQUENCY,EIGENSOLVER=LANCZOS, SIM=NO
**      *STEADY STATE DYNAMICS,SUBSPACE PROJECTION
*AMPLITUDE,NAME=CONST
0,1,1,1
*RESTART,WRITE, FREQ=4
*NODE, Nset=ODB_NSET
1,
2,1.,
3,1.,1.1,
4,,1.1
11,-1.,1.0
12,2.,1.0
13,2.,3.
14,-1.,3.
*ELEMENT,TYPE=CPE4,ELSET=BOT
1,1,2,3,4
*ELEMENT,TYPE=CPE4,ELSET=TOP
11,11,12,13,14
*SOLID SECTION,ELSET=TOP,MATERIAL=MAT
*SOLID SECTION,ELSET=BOT,MATERIAL=MAT
*MATERIAL,NAME=MAT
*ELASTIC
3.E7,0.0
*DENSITY
10.0, 
*surface,NAME=TOPS
TOP,S1
*surface,NAME=SECONDARY
BOT,S3
*CLEARANCE,SECONDARY=SECONDARY,MAIN=TOPS,TABULAR
3,0.0
4,0.0
*SURFACE INTERACTION,NAME=ISLS
*CONTACT PAIR,INTERACTION=ISLS, SMALL SLIDING
SECONDARY,TOPS
*BOUNDARY
1,1,2
2,2,2
13,2,2
14,1,2
**
** Sets added by addodb
**
*Elset,Elset=ODB_ESET
TOP, BOT
**
**
*Nset, Nset=ODB_NSET
ODB_NSET, 
**
** step 1: Clearance
**
*STEP,INC=50,NLGEOM
*STATIC
0.1,1.0
*BOUNDARY
11,2,2,-0.2
12,2,2,-0.2
*PRINT,CONTACT=YES
*NODE FILE,FREQ=100
U, 
RF, 
*EL FILE,FREQ=100
S, 
E, 
*OUTPUT,FIELD, OP=ADD,Frequency=100
*NODE OUTPUT
U, 
RF, 
*OUTPUT,HISTORY,OP=ADD,Frequency=100
*NODE OUTPUT, Nset=ODB_NSET
U,RF
*OUTPUT,FIELD, OP=ADD,Frequency=100
*ELEMENT OUTPUT
S, 
E, 
*OUTPUT,HISTORY,OP=ADD,Frequency=100
*ELEMENT OUTPUT, Elset=ODB_ESET
S, 
E, 
*END STEP
**
** step 2: Steady state dynamics, direct analysis
**
*STEP
*STEADY STATE DYNAMICS,DIRECT
4,20,5
*BOUNDARY
1,2,,0.01
2,2,,0.01
*NODE FILE
U, 
RF, 
*EL FILE
E, 
S, 
*CONTACT FILE
*OUTPUT,FIELD
*NODE OUTPUT
U, 
RF, 
*OUTPUT,HISTORY
*NODE OUTPUT, Nset=ODB_NSET
U, 
RF, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
E, 
S, 
*OUTPUT,HISTORY
*ELEMENT OUTPUT, Elset=ODB_ESET
E, 
S, 
*OUTPUT,FIELD
*CONTACT OUTPUT, VARIABLE=ALL
*OUTPUT,HISTORY
*CONTACT OUTPUT, VARIABLE=ALL
*END STEP
**
** step 3: frequency extraction ->all modes are extracted
**
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS, SIM=NO
4,0,,,
*BOUNDARY,BASE NAME=SHAKE
1,2,2
2,2,2
*CONTACT FILE,FREQ=0
*MODAL FILE
*NODE FILE,F=0
*EL FILE,F=0
*NODE PRINT,F=0
*EL PRINT,F=0
*OUTPUT,FIELD,Frequency=0
*END STEP
**
** step 4: Steady state dynamics, direct analysis
**
*STEP
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
4,20,5
*BASE MOTION,BASE NAME=SHAKE,AMPLITUDE=CONST,DOF=2,TYPE=DISPLACEMENT,SCALE=0.01
*NODE FILE
U, 
RF, 
*EL FILE
E, 
S, 
*MODAL FILE
BM, 
*CONTACT FILE
*OUTPUT,FIELD
*NODE OUTPUT
U, 
RF, 
*OUTPUT,HISTORY
*NODE OUTPUT, Nset=ODB_NSET
U, 
RF, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
E, 
S, 
*OUTPUT,HISTORY
*ELEMENT OUTPUT, Elset=ODB_ESET
E, 
S, 
*OUTPUT,FIELD
*MODAL OUTPUT,VARIABLE=ALL
BM, 
*OUTPUT,FIELD
*CONTACT OUTPUT, VARIABLE=ALL
*OUTPUT,HISTORY
*CONTACT OUTPUT, VARIABLE=ALL
*END STEP