*HEADING
 COPPER ROD IMPACT (AXISYMMETRIC)
**
**  This example simulates the high velocity impact of a copper rod
**  onto a rigid wall.  Similar tests have been performed to determine
**  the material constants for high-pressure equations of state, and
**  the test is sometimes described as the Taylor bar experiment.
**  Extremely high plastic strains develop at the crushed end of the
**  rod, resulting in severe local mesh distortion.
**
*RESTART, WRITE, NUMBER INTERVAL=4
**
**  Generate a 10 x 36 mesh for the bar
**
*NODE, NSET=CORNERS
   1, .0000, .0000
  11, .0032, .0000
 397, .0000, .0324
 407, .0032, .0324
*NGEN, NSET=HEAD
   1, 11, 1
*NGEN, NSET=TAIL
 397, 407, 1
*NFILL, NSET=ROD
 HEAD, TAIL, 36, 11
*NSET, NSET=AXIS, GENERATE
 1, 397, 11
*NSET, NSET=BODY, GENERATE
 12, 407, 1
*NSET,NSET=NOUT
397, 
*ELEMENT, TYPE=CAX4R
 1, 1, 2, 13, 12
*ELGEN, ELSET=ROD
 1, 10, 1, 1, 36, 11, 10
**
**  Representative material properties for copper
**
*MATERIAL, NAME=COPPER
*ELASTIC
 110.E9, .3
*PLASTIC
 314.E6, 
*DENSITY
 8.97E3, 
*SOLID SECTION, ELSET=ROD, MATERIAL=COPPER,CONTROLS=ENHANCED
*SECTION CONTROLS,NAME=ENHANCED,HOURGLASS=ENHANCED
**
**  Zero radial displacements are imposed along the symmetry axis.
**  To simulate the impact of the rod on a (frictionless) rigid wall,
**  zero axial displacements are prescribed at one end of the rod,
**  while all other nodes are subjected to an initial velocity field.
**
*BOUNDARY
 HEAD, 2
 AXIS, 1
*INITIAL CONDITIONS, TYPE=VELOCITY
 BODY, 2, -227.
**
*NSET,NSET=NOUT
HEAD,407
*ELSET,ELSET=EOUT,GENERATE
1,10,1
**
**  Run the problem for 80 microseconds
**
*ELSET,ELSET=QA_TEST_EOUT
EOUT,
*NSET,NSET=QA_TEST_NOUT
NOUT,
**
*STEP
*DYNAMIC, EXPLICIT
 , 80.E-6
*FILE OUTPUT,TIMEMARKS=YES,NUM=1
*EL FILE,ELSET=EOUT
PEEQ,
*NODE FILE,NSET=NOUT
U, 
*ENERGY FILE
***OUTPUT,HISTORY,TIME INTERVAL=0.
***ENERGY OUTPUT
**DT,
**ALLKE,
***NODE OUTPUT,NSET=NOUT
**U,
*OUTPUT,FIELD,OP=NEW, NUMBER=4
*NODE OUTPUT
U, 
*OUTPUT,FIELD,OP=NEW,TIMEMARKS=YES,NUM=1
*ELEMENT OUTPUT
PEEQ,MISES,PE,LE,NE
*OUTPUT,HISTORY,OP=NEW,TIME INTERVAL=1.6E-7
*ENERGY OUTPUT
ALLKE, 
*NODE OUTPUT,NSET=NOUT
U, 
*INCREMENTATION OUTPUT
DT, 
*OUTPUT, FIELD, TIMEMARKS=YES, NUMBER INTERVAL=1
*ELEMENT OUTPUT, ELSET=QA_TEST_EOUT
PEEQ,
*NODE OUTPUT, NSET=QA_TEST_NOUT
U, 
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=8E-05
*END STEP