*HEADING
 INDENTATION WITH A HEMISPHERICAL PUNCH
 ADAPTIVE MESHING VERIFICATION FOR CRUSHABLE FOAM PLASTICITY
 WITH ISOTROPIC HARDENING
 Units: N, mm
*NODE,NSET=ALLN
   1,   0., 300.
  15,   0.,   0.
 286, 600., 300.
 300, 600.,   0.
*NGEN,NSET=BOT
 15, 300, 15
*NGEN,NSET=CENTER
 1, 15, 1
*NGEN,NSET=LEFT
 1, 15, 1
*NGEN,NSET=RIGHT
 286, 300, 1
*NFILL,NSET=NALL
 LEFT, RIGHT, 19, 15
*ELEMENT,TYPE=CAX4R,ELSET=PLFOAM
 1, 2, 17, 16, 1
*ELGEN,ELSET=PLFOAM
 1, 14, 1, 1, 19, 15, 14
*ELSET,ELSET=UPPER,GEN
 1, 253, 14
*NODE,NSET=PUNCH
 1000, 0., 500.
*NSET,NSET=NOUT
 1000, 
*ELEMENT,TYPE=MASS,ELSET=PMASS
 2001, 1000
*MASS,ELSET=PMASS
 0.2, 
**
*SOLID SECTION,ELSET=PLFOAM,MATERIAL=PLFOAM,CONTROL=SECT
*SECTION CONTROLS,NAME=SECT,HOURGLASS=STIFFNESS
*MATERIAL,NAME=PLFOAM
*ELASTIC
 3000.0, 0.0
*CRUSHABLE FOAM,HARDENING=ISOTROPIC
 1.1, 0.2983
** plastic Poisson's ratio is chosen in such a way 
** that the plastic flow is associative.
*CRUSHABLE FOAM HARDENING
 220.0000,  0.0
 246.5121,  0.1
 272.9380,  0.2
 299.0192,  0.3
 324.5491,  0.4
 349.3453,  0.5
 373.2598,  0.6
 396.1709,  0.7
 418.0111,  0.8
 438.7210,  0.9
 458.2682,  1.0
 493.8430,  1.2
 524.8363,  1.4
 551.5284,  1.6
 574.3073,  1.8
 593.5948,  2.0
 629.3596,  2.5
 651.9947,  3.0
 683.3387,  5.0
 688.3322, 10.0
*DENSITY
 500.E-9,
** 
*BOUNDARY
 BOT,    1, 2
 CENTER, 1, 1
 RIGHT,  1, 1
 NOUT,   1
 NOUT,   6
*AMPLITUDE,NAME=RAMPP,TIME=TOTAL TIME,DEFINITION=SMOOTH STEP 
 0.0, 0.0, 0.6, -250.
*SURFACE,TYPE=SEGMENT,NAME=PUNCH,FILLET RADIUS=0.01
 START, 200., 700.
 LINE,  200., 500.
 CIRCL,   0., 300., 0., 500.
*SURFACE,TYPE=ELEMENT,NAME=TARGET,REGION TYPE=SLIDING
 UPPER,S3
*RIGID BODY,REF NODE=1000,ANALYTICAL SURFACE=PUNCH
*ELSET,ELSET=QA_TEST_UPPER
 UPPER,
*NSET,NSET=QA_TEST_NOUT
 NOUT,
**
*STEP
*DYNAMIC,EXPLICIT
 ,0.6
*BOUNDARY,AMPLITUDE=RAMPP
 1000, 2, 2, 1.
*SURFACE INTERACTION,NAME=IMP_TARG
*FRICTION
 0.4, 
*CONTACT PAIR,INTERACTION=IMP_TARG
 PUNCH,TARGET
*ADAPTIVE MESH,ELSET=PLFOAM
**
*FILE OUTPUT,NUMBER INTERVAL=4,TIMEMARKS=YES
*EL FILE,ELSET=UPPER
 PEEQ,
*NODE FILE,NSET=NOUT
 U
**
*OUTPUT,FIELD,NUMBER INTERVAL=4,TIMEMARKS=YES
*ELEMENT OUTPUT,ELSET=QA_TEST_UPPER
 PEEQ,
*NODE OUTPUT
 U,
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIME INTERVAL=0.1
**
*END STEP