*HEADING
: SIMPLE CONTACT TEST, FINITE SLIDING CONTACT WITHOUT FRICTION, RIGID
** Contact Features Tested:
**      permeable-impermeable pore pressure contact for pore pressure coupled contact
*NODE, NSET=BLOCK
*NODE
*NSET,NSET=TOP
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT, TYPE=C3D8P, ELSET=BLOCK
*SOLID SECTION, ELSET=BLOCK, MATERIAL=SOFT
*MATERIAL, NAME=SOFT
*ELASTIC
 30.E6, 0.3
*DENSITY
10.,
*PERMEABILITY
1.e-4,
*INITIAL CONDITIONS,TYPE=RATIO
BLOCK,1.0
*rigid body,analytical surface=BSURF, REFNODE=9999
*surface,TYPE=CYLINDER, NAME=BSURF
 START, 0., 0.
 LINE, 10., 0.
*surface,NAME=ASURF
BLOCK,S1
*CONTACT PAIR,INTERACTION=SMOOTH
ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
*FRICTION
0.0,
**
** Step 1:  Bring block into contact and push against rigid surface
**
*STEP,AMP=RAMP
*SOIL,CONSOLIDATION
*CONTROLS,PARAMETERS=FIELD
 0.000001,
*BOUNDARY
 1, 1
 2, 1
 2, 3
 9999, 1,6
 TOP, 2,, -1.00001
 TOP,8,8,1
*PRINT, CONTACT=YES
*CONTACT PRINT,SLAVE=ASURF
*CONTACT FILE,SLAVE=ASURF
*EL FILE
 S,
*NODE FILE,F=100
U,
POR,
*NODE PRINT
U,POR
*OUTPUT,FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF
*END STEP
**
** Step 2: Apply normal load on top surface of block
**
*STEP,AMP=RAMP
*SOIL, CONSOLIDATION
*BOUNDARY,OP=NEW
 1, 1
 2, 1
 2, 3
 9999, 1,6
 TOP,8,8,1
*DLOAD
 BLOCK, P2, 100.E3
*PRINT, CONTACT=YES
*CONTACT FILE,SLAVE=ASURF
*EL FILE
 S,
*NODE FILE,F=100
U,POR
*NODE PRINT
U,POR
*END STEP
**
** Step 3: Slide block by fixing one end and pulling on the other
**
*STEP,AMP=RAMP,NLGEOM
*SOIL, CONSOLIDATION
*BOUNDARY,OP=NEW
 LEFT,1,,-2.0
 RIGHT,1
 6,3
 9999, 1,6
 TOP,8,8,1
*CONTACT FILE,SLAVE=ASURF
*EL FILE
 S,
*NODE FILE,F=100
U,
*END STEP
**
** step 4: Steady state dynamics, direct analysis
**
*STEP
*STEADY STATE DYNAMICS,DIRECT
20,80,5
*BOUNDARY,OP=NEW
 1, 1
 2, 1
 2, 3
 9999, 1,6
 TOP,8,8,1
*DLOAD
 BLOCK, P2, 100.E3
*NODE FILE
U,
RF,
*EL FILE
E,
S,
*CONTACT FILE
*OUTPUT,FIELD
*NODE OUTPUT
U,
RF,
*ELEMENT OUTPUT
E,
S,
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF
*END STEP
**
** step 5: frequency extraction
**
*STEP
*FREQUENCY, eigensolver=lanczos
17, ,
*BOUNDARY,OP=NEW
 1, 1
 2, 1
 2, 3
 9999, 1,6
 TOP,8,8,1
*MODAL FILE
*NODE FILE,F=0
*EL FILE,F=0
*NODE PRINT,F=0
*EL PRINT,F=0
*OUTPUT,FIELD,OP=NEW
*NODE OUTPUT
U,
*CONTACT FILE,FREQ=0
*END STEP
**
** step 6: Steady state dynamics, subspace projection
**
*STEP
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
20,80,5
*DLOAD
 BLOCK, P2, 100.E3
*NODE FILE
U,
*EL FILE
E,
S,
*OUTPUT,FIELD,OP=NEW
*NODE OUTPUT
U,
*ELEMENT OUTPUT
*CONTACT FILE
*END STEP