*HEADING
: SURFACE-BASED PRESSURE PENETRATION TEST 
WITH RIGID SURFACE, CAX4
*RESTART,WRITE, FREQ=1
*NODE
1,1.0,12.0
2,1.0,10.5
3,1.0,9.0
4,1.0,7.5
5,1.0,6.0
6,1.0,4.5
7,1.0,3.0
8,1.0,1.5
9,1.0,0
10,2.0,12.0
11,2.0,10.5
12,2.0,9.0
13,2.0,7.5
14,2.0,6.0
15,2.0,4.5
16,2.0,3.0
17,2.0,1.5
18,2.0,0
19,1.0,14.0
20,1.0,12.0
21,1.0,10.0
22,1.0,8.0
23,1.0,6.0
24,1.0,4.0
25,1.0,2.0
26,1.0,0
27,1.0,-2.0
100,1.0,6.0
*nset,nset=print
100
*NSET,NSET=NALL
1,2,3,4,5,6,7,8,9,10,
11,12,13,14,15,16,17,18
*ELEMENT,TYPE=CAX4,ELSET=BOT
1,1,2,11,10
2,11,2,3,12
3,4,13,12,3
4,13,4,5,14
5,15,14,5,6
6,16,15,6,7
7,7,8,17,16
8,8,9,18,17
*ELEMENT,TYPE=MAX1,ELSET=TOP
10,19,20
11,20,21
12,21,22
13,22,23
14,23,24
15,24,25
16,25,26
17,26,27
*NSET,NSET=TOP
10,11,12,13,14,15,16,17,18
*RIGID BODY,ELSET=TOP,REF NODE=100
1.0,
*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, NOTHICK
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
18,1,1,-3.75E-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
1,19,800
9,27,800
*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
1,19,900
9,27,900
*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