*HEADING
: SIMPLE CONTACT TEST, FINITE SLIDING CONTACT WITHOUT FRICTION, RIGID
** Contact Features Tested: 
**	*FRICTION, *CHANGE FRICTION (mu=0 and mu=0.3)
*NODE
1,0.
3,5.
7,0.,5.
9,5.,5.
11,0.,0.,-1.
13,5.,0.,-1.
17,0.,5.,-1.
19,5.,5.,-1.
1000,0.
*ELEMENT,TYPE=C3D8,ELSET=SOLID
1, 1,3,13,11,7,9,19,17
*SOLID SECTION,ELSET=SOLID,MATERIAL=SOLMAT
*MATERIAL,NAME=SOLMAT
*ELASTIC
30.E6, .3
*rigid body,analytical surface=BSURF,REFNODE=1000
*surface,TYPE=CYLINDER,NAME=BSURF
0.,0.,0., 1.,0.,0.
0.,0.,-1.
START,-100.
LINE,  100.
*surface,NAME=ASURF
SOLID,S1
*CONTACT PAIR,INTERACTION=SMOOTH
ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
1.25, 
*FRICTION
0.0, 
*BOUNDARY
LEFT,1
TOP,3
1000,1,6
*NSET,NSET=LEFT
1,11,7,17
*NSET,NSET=TOP
11,13,17,19
**
** step 1: Apply surface pressure
**
*STEP
*STATIC
1.,1.
*CONTACT PRINT,SECONDARY=ASURF
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASURF
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD
*NODE OUTPUT
U,
*EL PRINT,FREQ=999
 S,
*NODE PRINT,FREQ=0
*DLOAD
1,P2,5000.
*END STEP
**
** step 2: Apply surface pressure, large disp.
**
*STEP,NLGEOM
*STATIC
 .5,1.
*DLOAD
1,P2,5000.
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASURF
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD
*NODE OUTPUT
U,
*END STEP
**
** step 3: Change surface friction
**
*STEP,NLGEOM
*STATIC
1.0,1.0
*CHANGE FRICTION,INTERACTION=SMOOTH
*FRICTION
0.3, 
*ENDSTEP
**
** step 4: slide to the left w/ friction
**
*STEP,NLGEOM
*STATIC
*BOUNDARY,OP=NEW
1000,1,6
TOP,1,1,1.0
*DLOAD
1,P2,5000.
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD
*NODE OUTPUT
U,
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASURF,MAIN=BSURF
*CONTACT PRINT,SECONDARY=ASURF,MAIN=BSURF
*PRINT,CONTACT=YES
*ENDSTEP