*HEADING
 STATIC INDENTATION OF ELASTOMERIC FOAM 
 WITH A HEMISPHERICAL PUNCH AND CAX4 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=CAX4H
 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=HYPER
*MATERIAL, NAME=HYPER
*HYPERELASTIC
80.0, 20.0
*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
****************************************************
*STEP,NLGEOM,INC=200,AMPLITUDE=RAMP
 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=9999
*NODE OUTPUT
 U,V,A,RF
*ELEMENT OUTPUT
 S,LE,NE
*OUTPUT,HISTORY,FREQUENCY=0
*Adaptive Mesh, Elset=Foam, Frequency=1, Mesh Sweeps=20,control=xx
*adaptive mesh control,name=xx,geometric=yes
1.,1.
*END STEP