*HEADING
: SURFACE-BASED PRESSURE PENETRATION TEST
WITH RIGID SURFACE, CAX4
*RESTART,WRITE, FREQ=1
*NODE
*nset,nset=print
*NSET,NSET=NALL
*ELEMENT,TYPE=CAX4,ELSET=BOT
*ELEMENT,TYPE=MAX1,ELSET=TOP
*NSET,NSET=TOP
*RIGID BODY,ELSET=TOP,REF NODE=100
*MEMBRANE SECTION,ELSET=TOP,MATERIAL=MAT
1
*SOLID SECTION,ELSET=BOT,MATERIAL=MAT
*MATERIAL,NAME=MAT
*ELASTIC
1.E5,0.3
*DENSITY
1.0E-3,
*surface,NAME=TOPS
TOP,SPOS
*surface,NAME=BOTS
1,S1
2,S2
3,S4
4,S2
5,S3
6,S3
7,S1
8,S1
*SURFACE INTERACTION,NAME=ISLS
1.0,
*FRICTION
0.,
*CONTACT PAIR,INTERACTION=ISLS, SMALL SLIDING
BOTS,TOPS
*BOUNDARY
100,1,2
100,6,6
NALL,2,2
*STEP,INC=20
*STATIC
0.1,1
*BOUNDARY
10,1,1,-3.75E-3
11,1,1,-5.25E-3
12,1,1,-7.35E-3
13,1,1,-10.29E-3
14,1,1,-14.406E-3
15,1,1,-10.29E-3
16,1,1,-7.35E-3
17,1,1,-5.25E-3
*NODE PRINT,NSET=TOP
U,
*node print,nset=print
rf
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=BOTS,MASTER=TOPS
CSTRESS,
*CONTACT PRINT,SLAVE=BOTS,MASTER=TOPS
CSTRESS,
*ENDSTEP
*STEP, INC=20
*STATIC
0.01,1
*PRESSURE PENETRATION,SLAVE=BOTS,MASTER=TOPS
*CONTACT PRINT,SLAVE=BOTS,MASTER=TOPS,TOTAL=YES,FREQ=1
CSTRESS,PPRESS,CDISP
*OUTPUT,FIELD,FREQ=1
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=BOTS,MASTER=TOPS
CSTRESS,PPRESS,CDISP
*PRINT,CONTACT=YES
*ENDSTEP
*STEP, INC=20
*STATIC
0.01,1
*PRESSURE PENETRATION,SLAVE=BOTS,MASTER=TOPS
*CONTACT PRINT,SLAVE=BOTS,MASTER=TOPS,TOTAL=YES,FREQ=1
CSTRESS,PPRESS,CDISP
*OUTPUT,FIELD,FREQ=1
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=BOTS,MASTER=TOPS
CSTRESS,PPRESS,CDISP
*PRINT,CONTACT=YES
*ENDSTEP