*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
1,0.,300.
19,0.,0.
481,300.,300.
499,300.,0.
601,600.,300.
619,600.,0.
*NODE,NSET=N9999
9999,0.,500.
*NSET,NSET=N1
1
*NSET,NSET=N19
19
*NSET,NSET=N481
481
*NSET,NSET=N499
499
*NFILL,NSET=TOP,BIAS=.85
N1,N481,12,40
*NGEN,NSET=TOP
481,601,40
*NFILL,NSET=BOT,BIAS=.85
N19,N499,12,40
*NGEN,NSET=BOT
499,619,40
*NFILL,NSET=ALLN,BIAS=.9
TOP,BOT,9,2
*NSET,NSET=CENTER,GENERATE
3,19
*PARAMETER SHAPE VARIATION, PARAMETER=L
1, 0.000,  0.00000, 0.000
41, 0.000,  0.01463, 0.000
81, 0.000,  0.03185, 0.000
121, 0.000,  0.05210, 0.000
161, 0.000,  0.07592, 0.000
201, 0.000,  0.10396, 0.000
241, 0.000,  0.13693, 0.000
281, 0.000,  0.17573, 0.000
321, 0.000,  0.22137, 0.000
361, 0.000,  0.27507, 0.000
401, 0.000,  0.33823, 0.000
441, 0.000,  0.41257, 0.000
481, 0.000,  0.50000, 0.000
521, 0.000,  0.66667, 0.000
561, 0.000,  0.83333, 0.000
601, 0.000,  1.00000, 0.000
*ELEMENT,TYPE=CAX4,ELSET=FOAM
1,3,43,41,1
*ELGEN,ELSET=FOAM
1,15,40,10,9,2,1
*ELSET,ELSET=CENT,GENERATE
1,141,10
61,69
*ELSET,ELSET=ETOP,GENERATE
1,101,10
*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
1.,0.,0., 0.,1.,0.
3,0.
*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
1, 59, 1
*NSET, NSET=NU, GENERATE
41,601,40
*NSET, NSET=NR, GENERATE
59,619,40
*STEP,NLGEOM,INC=200,AMPLITUDE=RAMP, DSA, UNSYMM=YES,Convert Sdi=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