*HEADING
HERTZ CONTACT PROBLEM - DISPLACEMENT CONTROL, 2D
USING PENALTY CONTACT w/ 10 TIMES DEFAULT PENALTY STIFFNESS
*NODE,NSET=REFNODE
9999,0.0,0.0
**
*NODE
1,   0.0,-10.0
10,3.420214,-9.3969262
19, 10.0,  0.0
**
501, 0.0, -7.8
510, 3.420214,-6.5
519, 7.8,  0.0
**
1401,0.0,0.0
1409,6.93333,0.0
**
2001,-1.0,-10.0
2002,10.0,-10.0
9999, 0.0,-11.0
*NGEN,LINE=C
1,10,1,,0.0,0.0
10,19,1,,0.0,0.0
*NSET,NSET=N1,GENERATE
1,10
*NSET,NSET=N2,GENERATE
10,19
*NGEN,NSET=N3
501,510,1
*NSET,NSET=N510
510,
*NSET,NSET=N519
519,
*NFILL,NSET=N4,BIAS=.9
N510,N519,9,1
*NSET,NSET=N5,GENERATE
501,509,1
*NGEN,NSET=N6
1401,1409
**
*NFILL
N1,N3,5,100
N2,N4,5,100
*NFILL,BIAS=.9
N5,N6,9,100
**
*NSET,NSET=CONTACT,GENERATE
1,7
**
*ELEMENT,TYPE=CPE4R,ELSET=SOLID
 1, 1,2,102,101
501,501,502,602,601
**
509, 509, 510,511,609
609, 609, 511,512,709
709, 709, 512,513,809
809, 809, 513,514,909
909, 909, 514,515,1009
1009,1009,515,516,1109
1109,1109,516,517,1209
1209,1209,517,518,1309
1309,1309,518,519,1409
**
*ELGEN,ELSET=SOLID
 1,18,1,1,5,100,100
501,8,1,1,9,100,100
*ELSET,ELSET=BOT,GENERATE
1,18,1
*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
2001,2001,2002
*NSET,NSET=LEFT,GEN
1,1401,100
*NSET,NSET=TOP,GEN
19,519,100
1401,1409,1
**
*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,CPSET=CC1
BOTSUR,RIGSUR
*CONTACT CONTROLS,CPSET=CC1,SCALE PENALTY=10.
*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