*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=NALL
*ELEMENT,TYPE=CPE4P,ELSET=SOLID
*SOLID SECTION,ELSET=SOLID,MATERIAL=MAT1
 .25,
*MATERIAL,NAME=MAT1
*ELASTIC
 30.E6,.3
*DENSITY
 10.0,
*PERMEABILITY
1.e-4,
*INITIAL CONDITIONS,TYPE=RATIO
NALL,1.0
*rigid body,analytical surface=bsurf,REFNODE=1000
*surface,TYPE=SEGMENTS,NAME=BSURF
 START,0.,0.
 LINE,5.5,0.
*surface,NAME=ASURF
1,S1
*CONTACT PAIR,INTERACTION=SMOOTH
ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
.25,
*FRICTION
0.0,
*STEP,INC=10
*SOIL,CONSOLIDATION
*CONTACT PRINT,SLAVE=ASURF
*CONTACT FILE,SLAVE=ASURF
*EL FILE
 S,
*NODE FILE,F=100
U,
POR,
*EL PRINT
 S,
*NODE PRINT,FREQ=100
POR,
*DLOAD
 1,P3,250.
*BOUNDARY
 1,1
 1000,1,2
 1000,6
 3,8,8,1
 4,8,8,1
*OUTPUT,FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF
*END STEP
*STEP,NLGEOM
  LARGE DISPLACEMENT ANALYSIS
*SOIL, CONSOLIDATION
*DLOAD
 1,P3,250.
*BOUNDARY
 1,1
 1000,1,2
 1000,6
 3,8,8,1
 4,8,8,1
*CONTACT FILE,SLAVE=ASURF
*EL FILE
 S,
*NODE PRINT,F=100
U,
POR,
*NODE FILE, FREQ=100
U,
POR,
*END STEP
**
** step 3: Steady state dynamics, direct analysis
**
*STEP
*STEADY STATE DYNAMICS,DIRECT
20,80,5
*DLOAD
 1,P3,100.
*BOUNDARY
 1,1
 1000,1,2
 1000,6
*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 4: frequency extraction ->all modes are extracted
**
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS
5,
*BOUNDARY
 1,1
 1000,1,2
 1000,6
*MODAL FILE
*NODE FILE,F=0
*EL FILE,F=0
*NODE PRINT,F=0
*EL PRINT,F=0
*CONTACT FILE,FREQ=0
*END STEP
**
** step 5: Steady state dynamics, subspace projection
**
*STEP
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
20,80,5
*DLOAD
 1,P3,100.
*NODE FILE
U,
RF,
*EL FILE
E,
S,
*CONTACT FILE
*END STEP