*HEADING
ADAPTIVE MESHING EXAMPLE 
OBLIQUE IMPACT OF A COPPER ROD(3D MODEL)
Units - N, m ,second
**
*RESTART,WRITE,NUMBER INTERVAL=30
*INCLUDE,INPUT=ale_rodimpac_inclined_nodelem.inp
**
** Representative material properties for copper
**
*MATERIAL,NAME=COPPER
*ELASTIC
 124.E9, 0.34
** Young's modulus unit: Pa
*PLASTIC,HARDENING=JOHNSON COOK
 90.E6, 292.E6, 0.31, 1.09, 1058., 25.
** A & B unit: Pa
*RATE DEPENDENT,TYPE=JOHNSON COOK
 0.025, 1.0
**unit:reference strain rate has a unit of s^(-1)
*INELASTIC HEAT FRACTION
 0.5,
** dimensionless
*SPECIFIC HEAT
 383.,
** unit: J/kgK
*EXPANSION
 5.E-5,
** unit: K^(-1)
*DENSITY
 8.96E3,
** unit: kg/m^3
*SOLID SECTION,ELSET=ROD,MATERIAL=COPPER,
CONTROLS=SECT
*SECTION CONTROLS,NAME=SECT,
KINEMATICS=ORTHOGONAL,HOURGLASS=VISCO
*BOUNDARY
 SPLANE1, 2
 200000,  1, 6
*INITIAL CONDITIONS,TYPE=TEMPERATURE
 ROD0, 25.
*INITIAL CONDITIONS,TYPE=VELOCITY
 BODY, 3, -340.
**
** Run the problem for 120 microseconds
**
*SURFACE,TYPE=CYLINDER,NAME=RSURF
 -20.E-3,0.,-13.3945E-3,20.E-3,0.,-13.3945E-3
 -20.E-3,1.,-13.3945E-3
 START, 0.,     0.
 LINE,  40.E-3, 23.094E-3
*SURFACE,TYPE=ELEMENT,NAME=RODSURF,
REGION TYPE=SLIDING
 ROD,
*RIGID BODY,REFNODE=200000,
ANALYTICAL SURFACE =RSURF
*STEP
*DYNAMIC,EXPLICIT,ADIABATIC 
 ,120.E-6
*CONTACT PAIR,INTERACTION=INTER
 RSURF,RODSURF
*SURFACE INTERACTION,NAME=INTER
*FRICTION
 0.2,
*FILE OUTPUT,TIMEMARKS=YES,NUM=1
*EL FILE, ELSET=QA_TEST
 PEEQ,
*NODE FILE, NSET=QA_TEST
 U,
*ENERGY FILE
**
*****
*NSET, NSET=QA_TEST_REFN
200000,
*ELSET, ELSET=QA_TEST
OUTER,
*NSET, NSET=QA_TEST
OUTER,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTERVAL=4
*ELEMENT OUTPUT, ELSET=QA_TEST

PEEQ,
*NODE OUTPUT, NSET=QA_TEST
U,
*NODE OUTPUT, NSET=QA_TEST_REFN
RF,
*OUTPUT, HIST, FREQ=999
*ENERGY OUTPUT, VAR=PRESELECT
******
*ADAPTIVE MESH,ELSET=ROD,FREQUENCY=5,
MESH SWEEPS=3
**** output for figures 
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=3 
*NODE OUTPUT 
U,
*OUTPUT, HIST, FREQ=1 
*ENERGY OUTPUT, VAR=PRESELECT 
*END STEP