*HEADING
 RIGID SPHERE IMPACTING ON THICK PLATE WITH
 PENETRATION AND ELEMENT EROSION
 SECTION CONTROLS USED (KINEMA=CENTROID,
 HOURGLASS=ENHANCED)
*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,CONTROL=B
*SOLID SECTION,ELSET=INF,  MATERIAL=INF,CONTROL=B
*SECTION CONTROLS, KINEMA=CENTROID,
HOURGLASS=ENHANCED, NAME=B
**
**  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
206.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
*NSET, NSET=N1
599991,
*OUTPUT,FIELD,OP=NEW,NUMBER INTERVAL=5,TIMEMARKS=NO
*ELEMENT OUTPUT
 STATUS 
*NODE OUTPUT
 U
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=N1
U,V
*ELSET, ELSET=TOUT
 TOP, TINF
*FILE OUTPUT, NUM=2, TIMEMARKS=YES
*NODE FILE, NSET=YSYM
 U,
*NODE FILE, NSET=N1
 U, V
*EL FILE, ELSET=TOUT
 STATUS,
*END STEP