*HEADING
2-D ELEMENT IN CONTACT WITH ANOTHER 2-D ELEMENT.
NORMALS ARE MODIFIED.
** Contact Features Tested:
**	*CLEARANCE
** Dynamic Features Tested:
**      *STEADY STATE DYNAMICS,DIRECT
**      *FREQUENCY,EIGENSOLVER=LANCZOS
**      *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.05,
4,,1.05
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=SLAVE
BOT,S3
*CLEARANCE,SLAVE=SLAVE,MASTER=TOPS,TABULAR
3,0.0,1.0,-1.0
4,0.0,1.0,-1.0
*SURFACE INTERACTION,NAME=ISLS
*CONTACT PAIR,INTERACTION=ISLS, SMALL SLIDING, TYPE=SURFACE TO SURFACE
SLAVE,TOPS
*BOUNDARY
1,1,2
2,1,2
11,1,1
13,2,2
14,1,2
12,1,1
**
** 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, OP=ADD
*NODE OUTPUT
U,
RF,
*OUTPUT,HISTORY,OP=ADD
*NODE OUTPUT, Nset=ODB_NSET
U,
RF,
*OUTPUT,FIELD, OP=ADD
*ELEMENT OUTPUT
E,
S,
*OUTPUT,HISTORY,OP=ADD
*ELEMENT OUTPUT, Elset=ODB_ESET
E,
S,
*OUTPUT,FIELD, OP=ADD
*CONTACT OUTPUT, VARIABLE=ALL
*OUTPUT,HISTORY,OP=ADD
*CONTACT OUTPUT, VARIABLE=ALL
*END STEP
**
** step 3: frequency extraction ->all modes are extracted
**
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS
3,0,,,
*BOUNDARY,BASE NAME=SHAKE
1,2,2
2,2,2
*CONTACT FILE,FREQ=0
*NODE FILE,F=0
*EL FILE,F=0
*MODAL FILE
*NODE PRINT,F=0
*EL PRINT,F=0
*OUTPUT,FIELD, OP=NEW,Frequency=0
*NODE OUTPUT
*ELEMENT OUTPUT
*OUTPUT,HISTORY,OP=NEW,Frequency=0
*NODE OUTPUT, Nset=ODB_NSET
*ELEMENT OUTPUT, Elset=ODB_ESET
*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, OP=ADD
*NODE OUTPUT
U,
RF,
*OUTPUT,HISTORY,OP=ADD
*NODE OUTPUT, Nset=ODB_NSET
U,
RF,
*OUTPUT,FIELD, OP=ADD
*ELEMENT OUTPUT
E,
S,
*OUTPUT,HISTORY,OP=ADD
*ELEMENT OUTPUT, Elset=ODB_ESET
E,
S,
*OUTPUT,FIELD, OP=ADD
*MODAL OUTPUT,VARIABLE=ALL
BM,
*OUTPUT,FIELD, OP=ADD
*CONTACT OUTPUT, VARIABLE=ALL
*OUTPUT,HISTORY,OP=ADD
*CONTACT OUTPUT, VARIABLE=ALL
*END STEP