*HEADING
**
** Reference for test case: vuinteraction_rbody_normal
**
*NODE,NSET=REFNOD
*RIGID BODY,REF NODE=1,ANALYTICAL SURFACE=SB1,ISOTHERMAL=YES
*HEATCAP,ELSET=PTCAPAC
*ELEMENT,TYPE=HEATCAP,ELSET=PTCAPAC
*BOUNDARY
1,1,6
**
*NODE,NSET=NB2
*ELEMENT,ELSET=EB2,TYPE=C3D8R
*ELGEN,ELSET=EB2
*NSET,NSET=LOAD_N
**
*ELSET,ELSET=EALL
EB2,
**
*SOLID SECTION,ELSET=EALL,MATERIAL=MAT1
*MATERIAL,NAME=MAT1
*ELASTIC
1.E6,0.2
*DENSITY
1.E1,
**
*SURFACE,TYPE=CYLINDER,NAME=SB1
START,0.,0.
LINE,6.,0.
*SURFACE,TYPE=ELEMENT,NAME=SB2
EB2, S3
**
*AMPLITUDE,NAME=SMOOTH,DEFINITION=SMOOTH STEP
0.,0.,1.E-3,-0.1,5.E-3,1.
**
*NSET, NSET=QA_TEST_ALLNODES
REFNOD, NB2
**
*SURFACE INTERACTION,NAME=Default_PP
*CONTACT DAMPING,DEFINITION=CRITICAL DAMPING FRACTION
**
** Pressure-overclosure slope = E(interface) / thick(interface)
** Shear traction slope       = G(interface) / thick(interface)
**
*SURFACE BEHAVIOR,NO SEPARATION,PRESSURE-OVERCLOSURE=LINEAR
 1.E6,
*FRICTION,ROUGH,SHEAR TRACTION SLOPE=416666.7
**
*STEP
*DYNAMIC,EXPLICIT,DIRECT USER CONTROL
5.E-5,5.E-3
**
*BOUNDARY,AMPLITUDE=SMOOTH
LOAD_N,2,2,0.001
**
*CONTACT
*CONTACT INCLUSIONS
SB1, SB2
*CONTACT PROPERTY ASSIGNMENT
SB1, SB2, Default_PP
**
*OUTPUT, FIELD, NUMBER INTERVAL=1
*NODE OUTPUT,NSET=QA_TEST_ALLNODES
U,
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.005
**
*END STEP