*HEADING
**
** cylinder indentation test with power-law rate dependence, 
** plane strain, regular mesh 10x30
** material: PUR60 (B), SI units
**
*NODE
  1,   0.,     0.
 31,   0.090,  0.
*NGEN,NSET=BOT
 1, 31, 1
*NODE
 311,  0.,     0.030
 341,  0.090,  0.030
*NGEN,NSET=TOP
 311, 341, 1
*NFILL
 BOT, TOP, 10, 31
*NSET,NSET=LFT,GEN
 1, 311, 31
*ELEMENT,TYPE=CPE4R
 1, 1, 2, 33, 32
*ELGEN,ELSET=FOAM
 1, 30, 1, 1, 10, 31, 30
*ELSET,ELSET=TOP,GENERATE
 271, 300, 1
*ELSET,ELSET=OUT
 1, 271
**
*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=ISOTROPIC
 1.0, 0.0
*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
*RATE DEPENDENT
 4638., 2.285
*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=VEL
 0., 0., 0.01, 5.4
** 
** STEP: Step-1
*STEP
 Step-1
*DYNAMIC,EXPLICIT
 ,0.01
**
*BOUNDARY,TYPE=VELOCITY,AMPLITUDE=VEL
 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=2.E-5
*NODE OUTPUT,NSET=PUNCHREF
 U, RF
*ELEMENT OUTPUT,ELSET=OUT
 MISES, PRESS, PEEQ
**
** Output requests for qa testing
**
*NSET,NSET=QA_TEST
PUNCHREF,
*ELSET,ELSET=QA_TEST
OUT,
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TEST
U, RF
*ELEMENT OUTPUT,ELSET=QA_TEST
PEEQ,
*OUTPUT,HISTORY,VAR=PRESELECT,TIME INTERVAL=0.01
*FILE OUTPUT,NUMBER INTERVAL=1
*EL FILE,ELSET=OUT
 PEEQ,
*NODE FILE,NSET=PUNCHREF
 U, RF
*END STEP