*HEADING
    STATIC INDENTATION OF AN ELASTOMERIC
    FOAM WITH A HEMISPHERICAL PUNCH
    DSA
    Step: Displace punch 250 mm downwards.
*PARAMETER
L=0.0
MU1=0.162449
*DESIGN PARAMETER
L,MU1
*NODE,NSET=ALLN
*NODE,NSET=N9999
*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
*PARAMETER SHAPE VARIATION, PARAMETER=L
*ELEMENT,TYPE=CAX4,ELSET=FOAM
*ELGEN,ELSET=FOAM
*ELSET,ELSET=CENT,GENERATE
*ELSET,ELSET=ETOP,GENERATE
*RIGID BODY,ANALYTICAL SURFACE=BSURF,REF NODE=9999
*SURFACE,TYPE=SEGMENTS,NAME=BSURF
START,141.42,641.42
CIRCL,-1.,300.,0.,500.
*SURFACE,NAME=ASURF
ETOP,S3
*CONTACT PAIR,INTERACTION=ROUGH
ASURF,BSURF
*SURFACE INTERACTION,NAME=ROUGH
*FRICTION
0.8,
*SOLID SECTION,ELSET=FOAM,MATERIAL=FOAM,
ORIENTATION=RECT
*ORIENTATION,NAME=RECT,SYSTEM=RECTANGULAR
*MATERIAL,NAME=FOAM
*HYPERFOAM,N=2
<MU1>, 8.89239, 3.597339e-05, -4.52156, 0.0, 0.0
*BOUNDARY
BOT,1,2
CENTER,1
N9999,1
N9999,6
*ELSET, ELSET=ER, GENERATE
*NSET, NSET=NU, GENERATE
*NSET, NSET=NR, GENERATE
*STEP,NLGEOM,INC=200,AMPLITUDE=RAMP, DSA, UNSYMM=YES
** Step 1: Move down punch.
*STATIC
.0015,1.,,.05
*CONTROLS, PARAMETERS=FIELD
1.E-5,
*BOUNDARY
9999,2,,-250.
**
** ODB OUTPUT REQUESTS
**
*OUTPUT, FIELD, FREQUENCY=999
*ELEMENT OUTPUT
S,
*NODE OUTPUT
U,
*NODE OUTPUT
RF,
*CONTACT OUTPUT
CSTRESS,CDISP
*DESIGN RESPONSE, FREQUENCY=999
*ELEMENT RESPONSE
S,
*NODE RESPONSE
U,
*NODE RESPONSE
RF,
*CONTACT RESPONSE
CSTRESS,CDISP
*END STEP