*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, TIMEMARKS=NO
*PREPRINT,ECHO=NO,MODEL=NO,HISTORY=NO
**
**  Generate a (10 x 2) x (36 x 2) mesh for the bar
**
*NODE, NSET=CORNERS
    1, .0000, .0000
   21, .0032, .0000
 1513, .0000, .0324
 1533, .0032, .0324
*NGEN, NSET=HEAD
 1, 21, 1
*NGEN, NSET=TAIL
 1513, 1533, 1
*NFILL, NSET=ROD
 HEAD, TAIL, 72, 21
*NSET, NSET=AXIS, GENERATE
 1, 1513, 21
*NSET, NSET=BODY, GENERATE
 22, 1533, 1
*ELEMENT, TYPE=CAX4R
 1, 1, 2, 23, 22
*ELGEN, ELSET=ROD
 1, 20, 1, 1, 72, 21, 20
**
**  Representative material properties for copper
**
*SOLID SECTION, ELSET=ROD, MATERIAL=COPPER
*MATERIAL, NAME=COPPER
*ELASTIC
 110.E9, .3
*PLASTIC
 314.E6,
*DENSITY
 8.97E3,
**
*ELSET, ELSET=ELOUT
 1,10,20
*NSET, NSET=NOUT
 1,11,21
**
**  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.
**
**  Run the problem for 80 microseconds
**
*ELSET,ELSET=QA_TEST_ELOUT
ELOUT,
*NSET,NSET=QA_TEST_NOUT
NOUT,
**
*STEP
*DYNAMIC, EXPLICIT
 ,80.E-6
*FILE OUTPUT,TIMEMARKS=YES,NUM=1
*EL FILE, ELSET=ELOUT
 PEEQ,
*NODE FILE, NSET=NOUT
 U,
*ENERGY FILE
**
***OUTPUT,HISTORY,TIME INTERVAL=80.E-9
***NODE OUTPUT,NSET=NOUT
**U,V
***ELEMENT OUTPUT,ELSET=ELOUT
**PEEQ,MISES
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLSE,ALLAE,ALLVD,ALLCD,ALLPD,ALLFD,ALLWK,ETOTAL,DT
*OUTPUT, FIELD, TIMEMARKS=YES, NUMBER INTERVAL=1
*ELEMENT OUTPUT, ELSET=QA_TEST_ELOUT
 PEEQ,
*NODE OUTPUT, NSET=QA_TEST_NOUT
 U,
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=8E-05
*END STEP