*HEADING
 DYNAMIC LOADING OF AN ELASTOMERIC FOAM WITH A HEMISPHERICAL PUNCH
******************************************************************* 
** Measure the punch reaction force during the following step:
**     Step 1: Displace punch 250 mm downwards and 
**             return punch to original position.
**             (Units: N, mm, sec)
*******************************************************************
 TEST 9-DIGIT LABELS
*PREPRINT,MODEL=YES,HISTORY=YES,CONTACT=YES
*RESTART,WRITE
*NODE
900000001,   0., 300.
900000019,   0.,   0.
900000481, 300., 300.
900000499, 300.,   0.
900000601, 600., 300.
900000619, 600.,   0.
*NSET,NSET=N1
900000001, 
*NSET,NSET=N19
900000019, 
*NSET,NSET=N481
900000481, 
*NSET,NSET=N499
900000499, 
*NFILL,NSET=TOP,BIAS=.85
 N1,N481,12,40
*NGEN,NSET=TOP
900000481,900000601,40
*NFILL,NSET=BOT,BIAS=.85
 N19,N499,12,40
*NGEN,NSET=BOT
900000499,900000619,40
*NFILL,NSET=ALLN,BIAS=.9
 TOP,BOT,9,2
*NSET,NSET=CENTER,GENERATE
900000001,900000019
*NSET,NSET=QA_TEST
ALLN, 
*ELEMENT,TYPE=CAX4R
900000001,900000003,900000043,900000041,900000001
*ELGEN,ELSET=FOAM
900000001,15,40,10,9,2,1
*ELSET,ELSET=CENT,GENERATE
900000001,900000141,10
900000061,900000069
*ELSET,ELSET=ETOP,GENERATE
900000001,900000101,10
*ELSET,ELSET=UPPER,GEN
900000001,900000141,10
*NODE,NSET=PUNCH
900001000,0.,510.
*SOLID SECTION,ELSET=FOAM,MATERIAL=FOAM,ORIENTATION=RECT
*ORIENTATION,NAME=RECT,SYSTEM=RECTANGULAR
 1.,0.,0., 0.,1.,0.
 3,0.
****************************************************
*MATERIAL,NAME=FOAM
*HYPERFOAM,N=2,TEST DATA INPUT,POISSON=0.0,MODULI=INSTANTANEOUS
*UNIAXIAL TEST DATA
 -.0217,   -.05
 -.0317,   -.10
 -.0367,   -.15
 -.0402,   -.20                              
 -.0433,   -.25
 -.0467,   -.30
 -.0504,   -.35
 -.0542,   -.40
 -.0604,   -.45
 -.0668,   -.50
 -.0759,   -.55
 -.0909,   -.60
 -.1083,   -.65
 -.1410,   -.70
 -.1933,   -.75
 -.2896,   -.80
*SIMPLE SHEAR TEST DATA      
  .0140, .08, .0046
  .0334, .16, .0166
  .0533, .24, .0366
  .0853, .32, .0573
  .1280, .40, .0817
  .1653, .48, .1098
  .2080, .56, .1394
  .2560, .64, .1666
  .2987, .72, .1904
*DENSITY
 1.E-11, 
****************************************************
*BOUNDARY
 BOT,1,2
 CENTER,1
900001000,1
900001000,6
****************************************************
*AMPLITUDE,NAME=SMOOTH,DEFINITION=TABULAR,TIME=TOTAL TIME
 0.0,0.0,1.0,-250.,2.0,0.
****************************************************
*SURFACE,TYPE=ELEMENT,NAME=TARGET
 UPPER,S3
*SURFACE, NAME=IMPACTOR,TYPE=SEGMENTS
 START,141.42,641.42
 CIRCL,0.,300., 0.,500.
*RIGID BODY,ANALYTICALSURFACE=IMPACTOR,REFNODE=900001000
****************************************************
*STEP
 Displace punch 250 mm downwards and return punch to original position.
*DYNAMIC,EXPLICIT
 ,2.0
*FIXED MASS SCALING, FACTOR=10
*BOUNDARY,AMPLITUDE=SMOOTH
900001000,2,2,1.
*SURFACE INTERACTION,NAME=IMP_TARG
*FRICTION
 0.8,
*CONTACT PAIR,INTERACTION=IMP_TARG
 TARGET,IMPACTOR
*OUTPUT,FIELD,TIME MARKS=YES,NUMBER INTER=10
*NODE OUTPUT
 U,V,A,RF
*ELEMENT OUTPUT
 S,LE,LEP,NE,NEP
*OUTPUT,HISTORY,TIME INTERVAL=8.E-4
*NODE OUTPUT,NSET=PUNCH
 U,V,A,RF
*ELEMENT OUTPUT,ELSET=UPPER
 S,LE,LEP,NE,NEP
*********************
*OUTPUT,HISTORY,TIME INTERVAL=2.0
*ENERGY OUTPUT,VARIABLE=PRESELECT
*************************
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TEST
U,
***********************
*END STEP