*HEADING
: SIMPLE CONTACT TEST WITHOUT FRICTION, SMALL SLIDING, RIGID SURFACE
*RESTART,WRITE
*NODE, NSET=BLOCK
*NODE
*NSET,NSET=TOP
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*NSET,NSET=MID
*ELEMENT, TYPE=C3D10MP, ELSET=BLOCK
*surface,NAME=ASURF
1,S2
3,S2
*CONTACT PAIR,INTERACTION=SMOOTH
ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
*FRICTION
0.0,
*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.
**
** Step 1:  Bring block into contact and push against rigid surface
**
*STEP,AMP=RAMP,NLGEOM
*SOILS,CONSOLIDATION
*CONTROLS,PARAMETERS=FIELD
 0.000001,
*BOUNDARY
 1, 1
 2, 1
 9, 1
 9, 3
 9999, 1,6
 TOP,2,, -1.00001
*CONTACT PRINT,SLAVE=ASURF,MASTER=BSURF
*OUTPUT,FIELD,FREQ=9999
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF,MASTER=BSURF
*OUTPUT,FIELD,FREQ=9999
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD,FREQ=9999
*NODE OUTPUT
U,
*END STEP
**
** Step 2: Apply normal load on top surface of block
**
*STEP,AMP=RAMP,NLGEOM
*SOILS,CONSOLIDATION
*PRINT,CONTACT=YES
*BOUNDARY,OP=NEW
 1, 1
 2, 1
 9, 3
 9, 1
 9999, 1,6
*DLOAD
 2, P3, 100.E3
 4, P4, 100.E3
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF,MASTER=BSURF
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD
*NODE OUTPUT
U,
*END STEP
**
** Step 3: Slide block by fixing one end and pulling on the other
**
*STEP,AMP=RAMP,NLGEOM,INC=20
*SOILS,CONSOLIDATION
*BOUNDARY,OP=NEW
 LEFT,1,,-1.0
 RIGHT,1
 MID,3
 9999, 1,6
TOP,8,8,1.
*OUTPUT,FIELD,FREQ=9999
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF,MASTER=BSURF
*OUTPUT,FIELD,FREQ=9999
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD,FREQ=9999
*NODE OUTPUT
U,
*END STEP
**
** step 4: change surface friction
**
*STEP,NLGEOM
*SOILS,CONSOLIDATION
*CHANGE FRICTION,INTERACTION=SMOOTH
*FRICTION
0.0,
*ENDSTEP
**
** step 5: slide to the left w/ friction
**
*STEP,NLGEOM
*STATIC
*BOUNDARY,OP=NEW
 LEFT,1,,-1.0
 RIGHT,1,,-1.0
 MID,3
 9999, 1,6
TOP,8,8,1.0
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD
*NODE OUTPUT
U,
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF,MASTER=BSURF
*CONTACT PRINT,SLAVE=ASURF,MASTER=BSURF
*PRINT,CONTACT=YES
*ENDSTEP