*HEADING
 RIGID SPHERE IMPACTING ON THICK PLATE WITH PENETRATION 
 AND ELEMENT EROSION
*PREPRINT,ECHO=NO,MODEL=NO,HISTORY=NO
*NODE
  1,-.0040,  0.0,0.0
  9, .0040,  0.0,0.0
 37,-.0031,.0031,0.0
 41,   0.0,.0040,0.0
 45, .0031,.0031,0.0
 46, .0040,  0.0,0.0
 50, .0031,.0031,0.0
 54,   0.0,.0040,0.0
 58,-.0031,.0031,0.0
 62,-.0040,  0.0,0.0
318, .02,0.0,0.0
334,-.02,0.0,0.0
488, .04,0.0,0.0
504,-.04,0.0,0.0
505, .05,0.0,0.0
521,-.05,0.0,0.0
*NGEN,NSET=A
1,9,1
*NGEN,NSET=B
37,41,1
41,45,1
*NGEN,NSET=C
46,50,1
50,54,1
54,58,1
58,62,1
*NGEN,NSET=D,LINE=C
318,334,1,,0.,0.,0.,0.,0.,1.
*NGEN,NSET=E,LINE=C
488,504,1,,0.,0.,0.,0.,0.,1.
*NGEN,NSET=TOP,LINE=C
505,521,1,,0.,0.,0.,0.,0.,1.
*NFILL,NSET=TOP
A,B,4,9
C,D,16,17
*NFILL,BIAS=.9,NSET=TOP
D,E,10,17
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=1000
0.,0.,-.001
0.,0.,0.,0.,0.,1.,0.0
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=2000
0.,0.,-.002
0.,0.,0.,0.,0.,1.,0.0
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=3000
0.,0.,-.003
0.,0.,0.,0.,0.,1.,0.0
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=4000
0.,0.,-.004
0.,0.,0.,0.,0.,1.,0.0
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=5000
0.,0.,-.005
0.,0.,0.,0.,0.,1.,0.0
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=6000
0.,0.,-.006
0.,0.,0.,0.,0.,1.,0.0
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=7000
0.,0.,-.007
0.,0.,0.,0.,0.,1.,0.0
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=8000
0.,0.,-.008
0.,0.,0.,0.,0.,1.,0.0
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=9000
0.,0.,-.009
0.,0.,0.,0.,0.,1.,0.0
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=10000
0.,0.,-.010
0.,0.,0.,0.,0.,1.,0.0
**
*NSET,NSET=YSYM,GEN
   1,   9, 1
  46, 488,17
  62, 504,17
1001,1009,1
1046,1488,17
1062,1504,17
2001,2009,1
2046,2488,17
2062,2504,17
3001,3009,1
3046,3488,17
3062,3504,17
4001,4009,1
4046,4488,17
4062,4504,17
5001,5009,1
5046,5488,17
5062,5504,17
6001,6009,1
6046,6488,17
6062,6504,17
7001,7009,1
7046,7488,17
7062,7504,17
8001,8009,1
8046,8488,17
8062,8504,17
9001,9009,1
9046,9488,17
9062,9504,17
10001,10009,1
10046,10488,17
10062,10504,17
*NSET,NSET=FRONT,GEN
   1, 385,1
**
*ELEMENT,TYPE=C3D8R,ELSET=TOP
 1, 1002,1003,1012,1011,  2, 3,12,11
19, 1008,1046,1047,1017,  8,46,47,17
20, 1017,1047,1048,1026, 17,47,48,26
21, 1026,1048,1049,1035, 26,48,49,35
22, 1035,1049,1050,1051, 35,49,50,51
23, 1034,1035,1051,1052, 34,35,51,52
24, 1033,1034,1052,1053, 33,34,52,53
25, 1032,1033,1053,1054, 32,33,53,54
26, 1031,1032,1054,1055, 31,32,54,55
27, 1030,1031,1055,1056, 30,31,55,56
28, 1029,1030,1056,1057, 29,30,56,57
29, 1059,1029,1057,1058, 59,29,57,58
30, 1060,1020,1029,1059, 60,20,29,59
31, 1061,1011,1020,1060, 61,11,20,60
32, 1062,1002,1011,1061, 62, 2,11,61
33, 1046,1063,1064,1047, 46,63,64,47
*ELGEN,ELSET=TOP
 1,  6,1,1,  3, 9, 6
33, 16,1,1, 26,17,16
*ELEMENT,TYPE=CIN3D8,ELSET=TINF
449, 489,488,1488,1489, 506,505,1505,1506
*ELGEN,ELSET=TINF
449, 16,1,1
**
*ELCOPY,OLD SET=TOP,NEW SET=PLATE,SHIFT NODES=1000,ELEMENT SHIFT=1000
*ELCOPY,OLD SET=TOP,NEW SET=PLATE,SHIFT NODES=2000,ELEMENT SHIFT=2000
*ELCOPY,OLD SET=TOP,NEW SET=PLATE,SHIFT NODES=3000,ELEMENT SHIFT=3000
*ELCOPY,OLD SET=TOP,NEW SET=PLATE,SHIFT NODES=4000,ELEMENT SHIFT=4000
*ELCOPY,OLD SET=TOP,NEW SET=PLATE,SHIFT NODES=5000,ELEMENT SHIFT=5000
*ELCOPY,OLD SET=TOP,NEW SET=PLATE,SHIFT NODES=6000,ELEMENT SHIFT=6000
*ELCOPY,OLD SET=TOP,NEW SET=PLATE,SHIFT NODES=7000,ELEMENT SHIFT=7000
*ELCOPY,OLD SET=TOP,NEW SET=PLATE,SHIFT NODES=8000,ELEMENT SHIFT=8000
*ELCOPY,OLD SET=TOP,NEW SET=PLATE,SHIFT NODES=9000,ELEMENT SHIFT=9000
*ELSET,ELSET=PLATE
TOP, 
**
*ELCOPY,OLD SET=TINF,NEW SET=INF,SHIFT NODES=1000,ELEMENT SHIFT=1000
*ELCOPY,OLD SET=TINF,NEW SET=INF,SHIFT NODES=2000,ELEMENT SHIFT=2000
*ELCOPY,OLD SET=TINF,NEW SET=INF,SHIFT NODES=3000,ELEMENT SHIFT=3000
*ELCOPY,OLD SET=TINF,NEW SET=INF,SHIFT NODES=4000,ELEMENT SHIFT=4000
*ELCOPY,OLD SET=TINF,NEW SET=INF,SHIFT NODES=5000,ELEMENT SHIFT=5000
*ELCOPY,OLD SET=TINF,NEW SET=INF,SHIFT NODES=6000,ELEMENT SHIFT=6000
*ELCOPY,OLD SET=TINF,NEW SET=INF,SHIFT NODES=7000,ELEMENT SHIFT=7000
*ELCOPY,OLD SET=TINF,NEW SET=INF,SHIFT NODES=8000,ELEMENT SHIFT=8000
*ELCOPY,OLD SET=TINF,NEW SET=INF,SHIFT NODES=9000,ELEMENT SHIFT=9000
*ELSET,ELSET=INF
TINF, 
**
*SOLID SECTION,ELSET=PLATE,MATERIAL=RHA
*SOLID SECTION,ELSET=INF,  MATERIAL=INF
**
**  Sphere with radius of .01 m.
**
*NODE,NSET=SPHERE,INPUT=sphere_n.inp
*ELEMENT,TYPE=R3D4,ELSET=SPHERE,INPUT=sphere_e.inp
*ELEMENT,TYPE=MASS,ELSET=MASS
50000,599991
*MASS,ELSET=MASS
**  Mass based on a density of 37240 Kg/m^3 and radius of .01 m.
.078, 
**
*MATERIAL,NAME=RHA
*DENSITY
7800., 
*ELASTIC
207.8E9,.3
*PLASTIC
1220.E6,0.
2440.E6,1.
*DAMAGE INITIATION, CRITERION=DUCTILE
1.0, 
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT 
0.0
***SHEAR FAILURE
**1.0,   
*MATERIAL,NAME=INF
*DENSITY
7800., 
*ELASTIC
206.8E9,.3
*INITIAL CONDITIONS,TYPE=VELOCITY
599991,1, 500.
599991,3,-866.6
*BOUNDARY
YSYM,2,2
599991, 2
599991, 4,6
*RESTART,WRITE,NUM=8
*SURFACE,TYPE=ELEMENT,NAME=SPHERE
SPHERE,SPOS
*NSET,NSET=ERODE_NSET,GEN
   1, 385,1
1001,1385,1
2001,2385,1
3001,3385,1
4001,4385,1
5001,5385,1
6001,6385,1
7001,7385,1
8001,8385,1
9001,9385,1
10001,10385,1
*SURFACE,TYPE=NODE,NAME=ERODE
ERODE_NSET,
*RIGID BODY,ELSET=SPHERE,REF NODE=599991
**
**
*STEP
*DYNAMIC,EXPLICIT
,40.E-6
*CONTACT PAIR
SPHERE,ERODE
***OUTPUT,HISTORY,TIME INTERVAL=0.0
*NSET, NSET=N1
599991,
***NODE OUTPUT,NSET=N1
**V,
*ELSET, ELSET=TOUT
 TOP, TINF
******
*NSET, NSET=QA_TEST
N1, YSYM
*ELSET, ELSET=QA_TEST
TOP,
*FILE OUTPUT, NUM=2, TIMEMARKS=YES
*NODE FILE, NSET=QA_TEST
 U,
*EL FILE, ELSET=QA_TEST
 PEEQ,
*OUTPUT,FIELD,OP=NEW,NUMBER INTERVAL=5,TIMEMARKS=NO
*ELEMENT OUTPUT, ELSET=QA_TEST
 PEEQ,
*NODE OUTPUT, NSET=QA_TEST
 U, 
*OUTPUT,FIELD,OP=NEW,NUMBER INTERVAL=5,TIMEMARKS=NO
*ELEMENT OUTPUT
 S, PEEQ, STATUS 
*NODE OUTPUT
 U, 
*OUTPUT,HISTORY,OP=NEW,TIME INTERVAL=8.E-8
*ENERGY OUTPUT
 ALLAE,ALLIE,ALLKE,ALLPD,ALLSE,ALLVD,ALLWK,ETOTAL
*END STEP