*HEADING
HERTZ CONTACT PROBLEM - DISPLACEMENT CONTROL, 3D
*Preprint, Contact=Yes, Model=Yes, History=Yes
*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
**
5001,-1.0,-10.0,-.1
5002,-1.0,-10.0,1.1
5003,10.0,-10.0,1.1
5004,10.0,-10.0,-.1
9999, 0.0,-11.0,0.5
*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,BIAS=.9
N1,N3,5,100
N2,N4,5,100
*NFILL
N5,N6,9,100
**
*NSET,NSET=FRONT,GENERATE
1,1500,1
*NCOPY, CHANGE NUMBER=2000,SHIFT,OLDSET=FRONT,NEWSET=BACK
0.0,0.0,1.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0
**
*NSET,NSET=CONTACT,GENERATE
1,7
2001,2007
**
*ELEMENT,TYPE=C3D8R,ELSET=SOLID
 1, 1,2,102,101,2001,2002,2102,2101,
501,501,502,602,601,2501,2502,2602,2601
**
509, 509, 510,511, 609,2509,2510,2511,2609
609, 609, 511,512, 709,2609,2511,2512,2709
709, 709, 512,513, 809,2709,2512,2513,2809
809, 809, 513,514, 909,2809,2513,2514,2909
909, 909, 514,515,1009,2909,2514,2515,3009
1009,1009,515,516,1109,3009,2515,2516,3109
1109,1109,516,517,1209,3109,2516,2517,3209
1209,1209,517,518,1309,3209,2517,2518,3309
1309,1309,518,519,1409,3309,2518,2519,3409
**
*ELGEN,ELSET=SOLID
 1,18,1,1,5,100,100
501,8,1,1,9,100,100
*ELSET,ELSET=BOTTOM,GENERATE
1,9,1
*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
5001,5001,5002,5003,5004
*NSET,NSET=LEFT,GEN
1,1401,100
2001,3401,100
*NSET,NSET=TOP,GEN
19,519,100
2019,2519,100
1401,1409,1
3401,3409,1
**
*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 Exterior
*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