*HEADING
HERTZ CONTACT PROBLEM - DISPLACEMENT CONTROL, C3D10M
*Preprint, Contact=Yes, Model=Yes, History=Yes
*RESTART,WRITE,NUMBER=1
*NODE,NSET=NALL
*NODE,NSET=RIGNODE
*NODE,NSET=REFNODE
*ELEMENT, TYPE=C3D10M, ELSET=EALL
*NSET, NSET=TOP1,GEN
*NSET, NSET=TOP2,GEN
*NSET, NSET=TOP3,GEN
*NSET, NSET=TOP
TOP1,TOP2,TOP3
*NSET, NSET=LEFT1,GEN
*NSET, NSET=LEFT2,GEN
*NSET, NSET=LEFT3,GEN
*NSET, NSET=LEFT
LEFT1,LEFT2,LEFT3,450,801,1251
*ELSET,ELSET=BOT
**
** Contact definition
**
*ELEMENT,TYPE=R3D4,ELSET=RIGID
*AMPLITUDE,NAME=SMOOTH,DEFINITION=SMOOTH STEP
0.0,0.0,0.01,1.0,
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*MATERIAL,NAME=A1
*ELASTIC
 30.E6,0.3
*DENSITY
7.5E-4,
**
*NSET, NSET=CONTACT
*NSET,NSET=OUTFIL,GENERATE
**
*RIGID BODY,ELSET=RIGID,REF NODE=REFNODE
*************************************************************
**General contact definition
*Contact
*Contact Inclusions, AllEl
*Contact Property Assignment
 , , Hertz
*SURFACE INTERACTION, NAME=Hertz
*************************************************************
*STEP, NLGEOM=NO
*DYNAMIC,EXPLICIT
,.011
*FIXED MASS SCALING,FACTOR=100,ELSET=EALL
*BOUNDARY
LEFT,1,1
REFNODE,1,6
NALL,3,3
*BOUNDARY,TYPE=DISPLACEMENT, AMPLITUDE=SMOOTH
TOP,2,2,-.4
**
*OUTPUT,FIELD,VAR=PRESELECT,NUMBER INTERVAL=5
**
** Output requests for qa testing
**
*NSET,NSET=QA_TEST1
REFNODE,
*NSET,NSET=QA_TEST2
OUTFIL,
*OUTPUT,FIELD,NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TEST1
RF,
*NODE OUTPUT, NSET=QA_TEST2
U,
*OUTPUT,HISTORY,VAR=PRESELECT,TIME INTERVAL=0.00275
*FILE OUTPUT,NUMBER=1
*NODE FILE,NSET=REFNODE
RF,
*NODE FILE,NSET=OUTFIL
U,
***OUTPUT,HISTORY,TIME INTERVAL=1.E-5
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLWK,ALLAE,ALLVD,ETOTAL,
**DT,
***NODE OUTPUT,NSET=TOP
**RF,
***NODE OUTPUT,NSET=CONTACT
**RF,
*END STEP