*HEADING
 DYNAMIC LOADING OF AN ELASTOMERIC FOAM WITH A HEMISPHERICAL PUNCH
 (using faceted rigid surface)
******************************************************************* 
** 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)
*******************************************************************
*PREPRINT,MODEL=YES,HISTORY=NO
*RESTART,WRITE
*NODE
   1,   0., 300.
  19,   0.,   0.
 481, 300., 300.
 499, 300.,   0.
 601, 600., 300.
 619, 600.,   0.
*NSET,NSET=N1
 1, 
*NSET,NSET=N19
 19, 
*NSET,NSET=N481
 481, 
*NSET,NSET=N499
 499, 
*NFILL,NSET=TOP,BIAS=.85
 N1,N481,12,40
*NGEN,NSET=TOP
 481,601,40
*NFILL,NSET=BOT,BIAS=.85
 N19,N499,12,40
*NGEN,NSET=BOT
 499,619,40
*NFILL,NSET=ALLN,BIAS=.9
 TOP,BOT,9,2
*NSET,NSET=CENTER,GENERATE
 1,19
*ELEMENT,TYPE=CAX4R
 1,3,43,41,1
*ELGEN,ELSET=FOAM
 1,15,40,10,9,2,1
*ELSET,ELSET=CENT,GENERATE
 1,141,10
 61,69
*ELSET,ELSET=ETOP,GENERATE
 1,101,10
*ELSET,ELSET=UPPER,GEN
 1,141,10
*NODE,NSET=PUNCH
 1000,  0.,510.
 1001,  0.,310.
 1031,200.,510.
*NGEN,NSET=PUNCH,LINE=C
 1001,1031,1,1000
*NSET,NSET=NOUT
 1000 
*ELEMENT,TYPE=RAX2,ELSET=PUNCH
 1001,1001,1002
*ELGEN,ELSET=PUNCH
 1001,30,1,1
*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
 1000,1
 1000,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,TYPE=ELEMENT,NAME=IMPACTOR
 PUNCH,SNEG
*RIGID BODY,REF NODE=1000,ELSET=PUNCH
****************************************************
*elset,elset=QA_TEST
1
*NSET,NSET=QA_TEST,ELSET=QA_TEST
*STEP
 Displace punch 250 mm downwards and return punch to original position.
*DYNAMIC,EXPLICIT
 ,2.0
*FIXED MASS SCALING, FACTOR=10
*BOUNDARY,AMPLITUDE=SMOOTH
 1000,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=NOUT
 U,V,A,RF
*ELEMENT OUTPUT,ELSET=UPPER
 S,LE,LEP,NE,NEP
*ENERGY OUTPUT,VARIABLE=PRESELECT
** Output Request for QA Purposes **
*OUTPUT,FIELD,NUMBER INTERVAL=2
*NODE OUTPUT,NSET=QA_TEST
U
*ELEMENT OUTPUT,ELSET=QA_TEST
LE
*END STEP