*HEADING
********************************************************
** 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
*NSET,NSET=N1
*NSET,NSET=N19
*NSET,NSET=N481
*NSET,NSET=N499
*NFILL,NSET=TOP,BIAS=.85
 N1,N481,12,40
*NGEN,NSET=TOP
*NFILL,NSET=BOT,BIAS=.85
 N19,N499,12,40
*NGEN,NSET=BOT
*NFILL,NSET=ALLN,BIAS=.9
 TOP,BOT,9,2
*NSET,NSET=CENTER,GENERATE
*ELEMENT,TYPE=CPE4H
*ELGEN,ELSET=FOAM
*ELSET,ELSET=CENT,GENERATE
*ELSET,ELSET=ETOP,GENERATE
*ELSET,ELSET=UPPER,GEN
*NODE,NSET=PUNCH
*SOLID SECTION,ELSET=FOAM,MATERIAL=HYPER,CONTROL=CNTRL-1
*MATERIAL, NAME=HYPER
*HYPERELASTIC
80.0, 20.0
*SECTION CONTROL,NAME=CNTRL-1,DISTORTION CONTROL=NO
*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=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=10,control=xx
*adaptive mesh control,name=xx,geometric=yes
*END STEP