*HEADING
 STATIC INDENTATION OF ELASTOMERIC FOAM 
********************************************************
** Measure the punch reaction force during the following:
**     Step 1: Displace punch 200 mm downwards.
********************************************************
*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=CPE3H
  1, 1, 3,43
150, 1,43,41
*ELGEN,ELSET=FOAM
  1,15,40,10,9,2,1
150,15,40,10,9,2,1
*ELSET,ELSET=CENT,GENERATE
 1,141,10
 61,69
*ELSET,ELSET=ETOP,GENERATE
 150,290,10
*ELSET,ELSET=UPPER,GEN
 1,141,10
*NODE,NSET=PUNCH
 1000,0.,510.
*SOLID SECTION,ELSET=FOAM,MATERIAL=FOAM
*MATERIAL,NAME=FOAM
*HYPERELASTIC
80.0, 20.0, 0.01
*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.
*Adaptive Mesh, Elset=Foam, Frequency=1, Mesh Sweeps=20
*END STEP