*HEADING
 STATIC INDENTATION OF ELASTOMERIC FOAM 
 WITH A HEMISPHERICAL PUNCH AND CPE4 elements 
********************************************************
** Measure the punch reaction force during the following 2 steps:
**     Step 1: Displace punch 200 mm downwards.
**     Step 2: ALE
**     Step 2: Displace punch 50mm further down.
**             (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=CPE4
 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.
*SOLID SECTION,ELSET=FOAM,MATERIAL=FOAM,CONTROLS=HGC
*SECTION CONTROLS,NAME=HGC,HOURGLASS=STIFFNESS
*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
*VISCOELASTIC,TIME=RELAXATION TEST DATA
*SHEAR TEST DATA,SHRINF=0.5000
 1.0000, 0.0001
 0.9695, 0.001
 0.9417, 0.002
 0.8722, 0.005
 0.7913, 0.010
 0.7043, 0.020
 0.6233, 0.050
 0.5736, 0.100
 0.5271, 0.200
 0.5013, 0.500
 0.5000, 1.000
*DENSITY
 1.E-11,
****************************************************
*BOUNDARY
 BOT,1,2     
 CENTER,1
 1000,1
 1000,6
****************************************************
*SURFACE,TYPE=SEGMENTS,NAME=IMPACTOR
 START,141.42,641.42
 CIRCL,-1.,300.,0.,500.
*SURFACE,TYPE=ELEMENT,NAME=TARGET
 ETOP,S3
*RIGID BODY,REF NODE=1000,ANALYTICAL SURFACE=IMPACTOR
*CONTACT PAIR,INTERACTION=ROUGH
 TARGET,IMPACTOR
*SURFACE INTERACTION,NAME=ROUGH
*FRICTION
 0.8, 
****************************************************
*STEP,NLGEOM,INC=200,AMPLITUDE=RAMP,UNSYMM=YES
 Displace punch 200 mm downwards.
*STATIC
 .0015,1.,,.05
*BOUNDARY
 1000,2,2,-200.
*PRINT,CONTACT=YES
*CONTACT CONTROLS,FRICTION ONSET=DELAY
*CONTACT PRINT,SLAVE=TARGET
**CONTACT FILE,SLAVE=TARGET,FREQUENCY=10
*ENERGY PRINT,FREQUENCY=5
*ENERGY PRINT,ELSET=FOAM,FREQ=5
*ENERGY PRINT,ELSET=CENT,FREQ=5
*EL PRINT,FREQUENCY=50,ELSET=CENT   
 S, 
 E,  
*NODE PRINT,FREQUENCY=25
 U,
 RF,
*OUTPUT,FIELD,FREQUENCY=4
*NODE OUTPUT
 U,V,A,RF
*ELEMENT OUTPUT
 S,LE,NE
*CONTACT OUTPUT,VARIABLE=PRESELECT,NSET=TOP
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=PUNCH
 U,V,A,RF
*ELEMENT OUTPUT,ELSET=CENT
 S,LE,NE
*ENERGY OUTPUT,VARIABLE=PRESELECT
*END STEP
****************************************************
*STEP,NLGEOM,INC=200,AMPLITUDE=RAMP,UNSYMM=YES
ADAPTIVE MESHING
*STATIC
 1.0, 1.0
*PRINT,CONTACT=YES
*CONTACT CONTROLS,FRICTION ONSET=DELAY
*CONTACT PRINT,SLAVE=TARGET
**CONTACT FILE,SLAVE=TARGET,FREQUENCY=10
*ENERGY PRINT,FREQUENCY=5
*ENERGY PRINT,ELSET=FOAM,FREQ=5
*ENERGY PRINT,ELSET=CENT,FREQ=5
*EL PRINT,FREQUENCY=50,ELSET=CENT   
 S, 
 E,  
*NODE PRINT,FREQUENCY=25
 U,
 RF,
*OUTPUT,FIELD,FREQUENCY=4
*NODE OUTPUT
 U,V,A,RF
*ELEMENT OUTPUT
 S,LE,NE
*CONTACT OUTPUT,VARIABLE=PRESELECT,NSET=TOP
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=PUNCH
 U,V,A,RF
*ELEMENT OUTPUT,ELSET=CENT
 S,LE,NE
*ENERGY OUTPUT,VARIABLE=PRESELECT
*Adaptive Mesh, Elset=Foam, Frequency=1, Mesh Sweeps=6,control=xx
*adaptive mesh control,name=xx,geometric=yes
1.,0.1
*END STEP
*******************************************************
*STEP,NLGEOM,INC=200,AMPLITUDE=RAMP,UNSYMM=YES
 Displace punch 50 mm further downwards.
*STATIC
 0.1, 1.0
*BOUNDARY
 1000,2,2,-250.
*PRINT,CONTACT=YES
*CONTACT CONTROLS,FRICTION ONSET=DELAY
*CONTACT PRINT,SLAVE=TARGET
**CONTACT FILE,SLAVE=TARGET,FREQUENCY=10
*ENERGY PRINT,FREQUENCY=5
*ENERGY PRINT,ELSET=FOAM,FREQ=5
*ENERGY PRINT,ELSET=CENT,FREQ=5
*EL PRINT,FREQUENCY=50,ELSET=CENT   
 S, 
 E,  
*NODE PRINT,FREQUENCY=25
 U,
 RF,
*OUTPUT,FIELD,FREQUENCY=4
*NODE OUTPUT
 U,V,A,RF
*ELEMENT OUTPUT
 S,LE,NE
*CONTACT OUTPUT,VARIABLE=PRESELECT,NSET=TOP
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=PUNCH
 U,V,A,RF
*ELEMENT OUTPUT,ELSET=CENT
 S,LE,NE
*ENERGY OUTPUT,VARIABLE=PRESELECT
*END STEP