*HEADING
**
** cylinder indentation test, plane strain, fine mesh 15x45
** material: PUR60 (B), SI units
**
*NODE
  1,   0.,     0.
 46,   0.090,  0.
*NGEN,NSET=BOT
 1, 46, 1
*NODE
 691,  0.,     0.030
 736,  0.090,  0.030
*NGEN,NSET=TOP
 691, 736, 1
*NFILL
 BOT, TOP, 15, 46
*NSET,NSET=LFT,GEN
 1, 691, 46
*ELEMENT,TYPE=CPE4R
 1, 1, 2, 48, 47
*ELGEN,ELSET=FOAM
 1, 45, 1, 1, 15, 46, 45
*ELSET,ELSET=TOP,GENERATE
 631, 675, 1
*ELSET,ELSET=OUT
 1, 631
**
*SECTION CONTROLS,NAME=FOAM,HOURGLASS=ENHANCED
*SOLID SECTION,ELSET=FOAM,MATERIAL=FOAM,CONTROLS=FOAM
 0.180,
*MATERIAL,NAME=FOAM
*ELASTIC
 7.5E6, 0.0
*CRUSHABLE FOAM,HARDENING=VOLUMETRIC
 1.0, 0.1
*CRUSHABLE FOAM HARDENING
 0.2000E6, 0.0000
 0.2577E6, 0.0094
 0.2760E6, 0.0258
 0.3053E6, 0.0452
 0.3267E6, 0.0655
 0.3623E6, 0.1084
 0.3891E6, 0.1540
 0.4250E6, 0.2405
 0.4568E6, 0.3812
 0.4738E6, 0.4600
 0.5170E6, 0.6391
 0.5862E6, 0.8570
 0.6503E6, 0.9857
 0.7470E6, 1.1324
 0.9820E6, 1.2965
 1.4702E6, 1.4808
 2.7262E6, 1.6609
 5.3911E6, 1.9000
*DENSITY
 60.,
** 
** SURFACE DEFINITIONS
*NODE,NSET=PUNCHREF
 1000,  0.0,    0.1125
*SURFACE,TYPE=SEGMENT,NAME=PUNCH
 START, 0.0825, 0.1125
 CIRCL, 0.0,    0.0300,   0.0,  0.1125
*RIGID BODY,REF NODE=PUNCHREF,ANALYTICAL SURFACE=PUNCH 
*SURFACE,TYPE=ELEMENT,NAME=FOAM
 TOP, S3
** 
** BOUNDARY CONDITIONS
*BOUNDARY
 LFT, XSYMM
 BOT, ENCASTRE
 PUNCHREF, 1
 PUNCHREF, 6
*AMPLITUDE,NAME=DISP,TIME=TOTAL TIME,DEFINITION=SMOOTH STEP
 0., 0., 0.05, 0.027
** 
** STEP: Step-1
*STEP
 Step-1
*DYNAMIC,EXPLICIT
 ,0.05
**
*BOUNDARY,AMPLITUDE=DISP
 PUNCHREF, 2, 2, -1.
**
** INTERACTION
*CONTACT PAIR,INTERACTION=INTPROP
 PUNCH, FOAM
*SURFACE INTERACTION,NAME=INTPROP
*FRICTION,TAUMAX=1.15E5
 0.2,
** 
** OUTPUT REQUESTS
*OUTPUT,FIELD,NUMBER INTERVAL=5
*NODE OUTPUT
 U,
*ELEMENT OUTPUT
 S, PE, PEEQ, LE
*OUTPUT,HISTORY,TIME INTERVAL=1.0E-4
*NODE OUTPUT,NSET=PUNCHREF
 U, RF
*ELEMENT OUTPUT,ELSET=OUT
 MISES, PRESS, PEEQ
**
** OUTPUT REQUESTS FOR QA PURPOSES
**
*NSET,NSET=QA_TEST
553,599,645,691,
*ELSET,ELSET=QA_TEST
586,631,
*OUTPUT,FIELD,NUMBER INTERVAL=1
*NODE OUTPUT,NSET=QA_TEST
 U
*ELEMENT OUTPUT,ELSET=QA_TEST
 PEEQ,
*OUTPUT,HISTORY,VAR=PRESELECT,TIME INTERVAL=1.25e-2
*FILE OUTPUT,NUMBER INTERVAL=1
*EL FILE,ELSET=QA_TEST
 PEEQ, 
*NODE FILE,NSET=PUNCHREF
 U, RF
**
*END STEP