*HEADING
**
** spherical indentation test, regular mesh 10x30
** material: PUR60 (B), SI units
**
*RESTART,WRITE,FREQ=999
*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=CAX4R
 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
**
*SOLID SECTION,ELSET=FOAM,MATERIAL=FOAM
*MATERIAL,NAME=FOAM
*ELASTIC
 7.5E6, 0.0
*CRUSHABLE FOAM,HARDENING=ISOTROPIC
 1.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
*CONTACT PAIR,INTERACTION=INTPROP
 FOAM, PUNCH
*SURFACE INTERACTION,NAME=INTPROP
*FRICTION,TAUMAX=1.15E5
 0.2,
** 
** BOUNDARY CONDITIONS
*BOUNDARY
 LFT, 1
 BOT, ENCASTRE
 PUNCHREF, 1
 PUNCHREF, 6
** 
*AMPLITUDE,NAME=VEL
 0., 0., 0.01, 5.4
** STEP: Step-1
*STEP,INC=500,NLGEOM,UNSYMM=YES
 Step-1
*STATIC
 0.00001, .01
**
** APPLIED LOAD
*BOUNDARY,TYPE=VELOCITY,AMPLITUDE=VEL
 PUNCHREF, 2, 2, -1.
**BOUNDARY
** PUNCHREF, 2, 2, -0.0270
**
** Relax tolerance from a default of 0.005 to 0.2 because the time and
** space averaging produce values that are very small, thus requiring
** convergence with unnecessarily small residuals. Relax displacement
** tolerance from 0.01 to 0.1
*CONTROLS,PARAMETERS=FIELD
 0.2, 0.1
**
*MONITOR,NODE=PUNCHREF,DOF=2
*ENERGY PRINT,FREQUENCY=999
*NODE FILE,NSET=PUNCHREF
 U, RF
*OUTPUT,FIELD,FREQ=100
*ELEMENT OUTPUT
 S, PE, PEEQ, LE
*NODE OUTPUT
 U,
*OUTPUT,HISTORY,FRE=1
*ELEMENT OUTPUT,ELSET=OUT
 MISES, PRESS, PEEQ
*NODE OUTPUT,NSET=PUNCHREF
 U, RF
*END STEP
*******************************************
** STEP: Step-2
*STEP,INC=500,NLGEOM,UNSYMM=YES
 Step-2
*STATIC
 0.0001, .01
**
** APPLIED LOAD
*BOUNDARY,TYPE=VELOCITY,AMPLITUDE=VEL
 PUNCHREF, 2, 2, 1.
**BOUNDARY
** PUNCHREF, 2, 2, 0
*END STEP