*HEADING
: SIMPLE CONTACT TEST, FINITE SLIDING WITHOUT FRICTION, RIGID 
** Contact Features Tested: 
**	*FRICTION, *CHANGE FRICTION (mu=0 and mu=0.3)
*NODE
1,0.
2,1.
1000,0.
*ELEMENT,TYPE=B31,ELSET=BEAM
1,1,2
*BEAM SECTION,SECTION=RECT,ELSET=BEAM,MATERIAL=BEAMAT
.25,.25
*MATERIAL,NAME=BEAMAT
*ELASTIC
30.E6,.3
*rigid body,analytical surface=bsurf,REFNODE=1000
*surface,TYPE=CYLINDER,NAME=BSURF
0.,-1.,0.,    1.,-1.,0.
0.,-1.,-1.
START,0.,1.
LINE,10.,1.
*surface,NAME=ASURF
BEAM,
*CONTACT PAIR,INTERACTION=SMOOTH
ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
*FRICTION
0.0, 
**
** step 1: Apply surface pressure
**
*STEP
*STATIC
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASURF
CSTRESS,CDISP
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD
*NODE OUTPUT
U,
*NODE PRINT,FREQ=0
*EL PRINT
 S,
*CONTACT PRINT,SECONDARY=ASURF
CSTRESS,CDISP
CAREA,
CFN,CFS
CMN,CMS
CFT,CMT
XN,XS,XT
*DLOAD
1,PY,-500.
*BOUNDARY
1,1
1,3
1,4,6
1000,1,6
*END STEP
**
** step 2: Apply surface pressure, large disp.
**
*STEP,NLGEOM
 LARGE DISPLACEMENT ANALYSIS
*STATIC
.5,1.
*DLOAD
1,PY,-500.
*BOUNDARY
1,1
1,3
1,4,6
1000,1,6
*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
1,1,,0.1
2,1,,0.1
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD
*NODE OUTPUT
U,
*ENDSTEP