*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
*NGEN,NSET=A
*NGEN,NSET=B
*NGEN,NSET=C
*NGEN,NSET=D,LINE=C
*NGEN,NSET=E,LINE=C
*NGEN,NSET=TOP,LINE=C
*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
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=2000
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=3000
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=4000
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=5000
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=6000
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=7000
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=8000
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=9000
*NCOPY,OLD SET=TOP,SHIFT,CHANGE NUMBER=10000
**
*NSET,NSET=YSYM,GEN
*NSET,NSET=FRONT,GEN
**
*ELEMENT,TYPE=C3D8R,ELSET=TOP
*ELGEN,ELSET=TOP
*ELEMENT,TYPE=CIN3D8,ELSET=TINF
*ELGEN,ELSET=TINF
**
*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
*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.
*SHEAR FAILURE
*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
*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
*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