*HEADING
HERTZ CONTACT PROBLEM - DISPLACEMENT CONTROL, 3D
USING PENALTY CONTACT
*Preprint, Contact=Yes, Model=Yes, History=Yes
*NODE,NSET=REFNODE
**
*NODE
**
**
**
*NGEN,LINE=C
*NSET,NSET=N1,GENERATE
*NSET,NSET=N2,GENERATE
*NGEN,NSET=N3
*NSET,NSET=N510
*NSET,NSET=N519
*NFILL,NSET=N4,BIAS=.9
N510,N519,9,1
*NSET,NSET=N5,GENERATE
*NGEN,NSET=N6
**
*NFILL,BIAS=.9
N1,N3,5,100
N2,N4,5,100
*NFILL
N5,N6,9,100
**
*NSET,NSET=FRONT,GENERATE
*NCOPY, CHANGE NUMBER=2000,SHIFT,OLDSET=FRONT,NEWSET=BACK
**
*NSET,NSET=CONTACT,GENERATE
**
*ELEMENT,TYPE=C3D8R,ELSET=SOLID
**
**
*ELGEN,ELSET=SOLID
*ELSET,ELSET=BOTTOM,GENERATE
*SOLID SECTION,ELSET=SOLID,MATERIAL=A1
1.0,
*MATERIAL,NAME=A1
*ELASTIC
 30.E6,0.3
*DENSITY
 7.5E-4,
**
*ELEMENT,TYPE=R3D4,ELSET=RIGID
*NSET,NSET=LEFT,GEN
*NSET,NSET=TOP,GEN
**
*AMPLITUDE,NAME=SMOOTH,DEFINITION=SMOOTH STEP
0.0,0.0,0.01,1.0,
**
*RIGID BODY,ELSET=RIGID,REF NODE=REFNODE
*************************************************************
**General contact definition
*Contact
*Contact Inclusions, All Element Based
*Contact Property Assignment
 , ,NO_FRIC
*SURFACE INTERACTION,NAME=NO_FRIC
*************************************************************
*STEP,NLGEOM=NO
*DYNAMIC,EXPLICIT
,.011
*BOUNDARY
LEFT,1,1
REFNODE,1,6
FRONT,3,3
BACK,3,3
*BOUNDARY,TYPE=DISPLACEMENT, AMPLITUDE=SMOOTH
TOP,2,2,-.4
**
*RESTART,WRITE,NUMBER=1
**
*OUTPUT,FIELD,VAR=PRESELECT,NUMBER INTERVAL=5
**
** Output requests for qa testing
**
*NSET,NSET=QA_TEST1
REFNODE,
*NSET,NSET=QA_TEST2
N3,
*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=N3
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