*HEADING
HERTZ CONTACT PROBLEM - DISPLACEMENT CONTROL, 2D
USING PENALTY CONTACT
*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
N1,N3,5,100
N2,N4,5,100
*NFILL,BIAS=.9
N5,N6,9,100
**
*NSET,NSET=CONTACT,GENERATE
**
*ELEMENT,TYPE=CPE4R,ELSET=SOLID
**
**
*ELGEN,ELSET=SOLID
*ELSET,ELSET=BOT,GENERATE
*SOLID SECTION,ELSET=SOLID,MATERIAL=A1
1.0,
*MATERIAL,NAME=A1
*ELASTIC
 30.E6,0.3
*DENSITY
 7.5E-4,
**
*ELEMENT,TYPE=R2D2,ELSET=RIGID
*NSET,NSET=LEFT,GEN
*NSET,NSET=TOP,GEN
**
*AMPLITUDE,NAME=SMOOTH,DEFINITION=SMOOTH STEP
0.0,0.0,0.01,1.0,
**
*SURFACE,TYPE=ELEMENT,NAME=RIGSUR
RIGID,SPOS
*SURFACE,TYPE=ELEMENT,NAME=BOTSUR
BOT,S1
*RIGID BODY,ELSET=RIGID,REF NODE=REFNODE
*STEP,NLGEOM=NO
*DYNAMIC,EXPLICIT
,.011
*BOUNDARY
LEFT,1,1
REFNODE,1,6
*BOUNDARY,TYPE=DISPLACEMENT, AMPLITUDE=SMOOTH
TOP,2,2,-.4
*CONTACT PAIR,INTERACTION=NO_FRI,MECHANICAL CONSTRAINT=PENALTY
BOTSUR,RIGSUR
*SURFACE INTERACTION,NAME=NO_FRI
**
*RESTART,WRITE,NUMBER=1
**
*OUTPUT,FIELD,VARIABLE=PRESELECT,NUMBER INTERVAL=4
*OUTPUT,FIELD,OP=ADD,NUMBER INTERVAL=4
*CONTACT OUTPUT
CSTRESS, CFORCE
*OUTPUT,HISTORY,TIME INTERVAL=1.E-5
*NODE OUTPUT,NSET=TOP
RF,
**
** 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
*END STEP