*HEADING
 FRICTIONAL CONTACT BETWEEN A DEFORMABLE BLOCK OF C3D8R ELEMENTS AND 
 AN ANALYTICAL RIGID SURFACE WITH COULOMB FRICTION MODEL
**
** DEFORMABLE ELEMENT
**
*NODE, NSET=NELEM
101, -5.0, 0.0, 0.0
102,  0.0, 0.0, 0.0
103,  0.0, 1.0, 0.0
104, -5.0, 1.0, 0.0
111, -5.0, 0.0, 1.0
112,  0.0, 0.0, 1.0
113,  0.0, 1.0, 1.0
114, -5.0, 1.0, 1.0
**
*NSET, NSET=NTOP
114, 113, 103, 104
*NSET, NSET=NBOT
111, 112, 102, 101
**
*ELEMENT, TYPE=C3D8R, ELSET=ELEM
1001, 101, 111, 112, 102, 104, 114, 113, 103
*SOLID SECTION, ELSET=ELEM, MATERIAL=METAL
 1., 
**
*MATERIAL, NAME=METAL
*ELASTIC
 3.E7,0.33
*DENSITY
 7.3E-4, 
**
*NODE, NSET=REFNODE
201, 5.0,-1.0
**
*ELEMENT, TYPE=MASS, ELSET=RMASS
200, 201
*MASS, ELSET=RMASS
1.0, 
**
*BOUNDARY
201,1,6
**
*INITIAL CONDITIONS, TYPE=VELOCITY
NELEM, 1, 200.0
*INITIAL CONDITIONS, TYPE=STRESS
ELEM, 0.0, -2000.0, 0.0
**
*SURFACE, TYPE=CYLINDER, NAME=RSURF
-10., -10., 1., 0., -10., 1.
-10., -10., 0.
START, 0., 10.
LINE, 20., 10.
*SURFACE, TYPE=ELEMENT, NAME=SELEM
ELEM, S1
*SURFACE, TYPE=ELEMENT, NAME=SLOAD
ELEM, S2
**
*RIGID BODY, REF NODE=201, ANALYTICAL SURFACE =RSURF
**
*SURFACE INTERACTION, NAME=RSUR_PM
*FRICTION, USER=FRICTION, PROPERTIES=3
0.15, 0.05, 0.01
**
*ELSET,ELSET=QA_TEST_All_ELSETS_AUTO
ELEM, RMASS
*NSET,NSET=QA_TEST_All_NSETS_AUTO
NELEM, REFNODE
**
*STEP
*DYNAMIC,EXPLICIT
,0.001
**
*CONTACT
*CONTACT INCLUSIONS, ALL EXTERIOR
*CONTACT PROPERTY ASSIGNMENT
SELEM, RSURF, RSUR_PM
**
*DSLOAD
SLOAD, P, 2000.0
**
**  Some typical output requests.
**
*RESTART,TIMEMARKS=NO, WRITE, NUMBER=4
*FILE OUTPUT, NUMBER=4, TIMEMARKS=YES
*ENERGY FILE
*NODE FILE
 U,
**
*OUTPUT,FIELD,VARIABLE=PRESELECT,NUMBER INTERVAL=4,OP=NEW
*OUTPUT,HISTORY,TIME INTERVAL=5.e-5,OP=NEW
*NODE OUTPUT,NSET=QA_TEST_All_NSETS_AUTO
 U, V
*OUTPUT, FIELD, NUMBER=4, TIMEMARKS=YES
*NODE OUTPUT,NSET=QA_TEST_All_NSETS_AUTO
 U, V
*ELEMENT OUTPUT,ELSET=ELEM
 S, 
*Contact Output
CDISP
*END STEP