*HEADING
: SURFACE-BASED PRESSURE PENETRATION TEST 
DEFORMABLE vs. DEFORMABLE SURFACES, 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
61,1.0,12.0
62,1.0,10.5
63,1.0,9.0
64,1.0,7.5
65,1.0,6.0
66,1.0,4.5
67,1.0,3.0
68,1.0,1.5
69,1.0,0
70,0.0,12.0
71,0.0,10.5
72,0.0,9.0
73,0.0,7.5
74,0.0,6.0
75,0.0,4.5
76,0.0,3.0
77,0.0,1.5
78,0.0,0
*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=CAX4,ELSET=TOP
9,62,61,70,71
10,63,62,71,72
11,73,64,63,72
12,74,65,64,73
13,74,75,66,65
14,75,76,67,66
15,67,76,77,68
16,68,77,78,69
*NSET,NSET=TOP
70,71,72,73,74,75,76,77,78
*SOLID SECTION,ELSET=BOT,MATERIAL=MAT
*SOLID SECTION,ELSET=TOP,MATERIAL=MAT
*MATERIAL,NAME=MAT
*ELASTIC
1.E5,0.3
*DENSITY
1.0E-3, 
*surface,NAME=TOPS
9,S1
10,S1
11,S2
12,S2
13,S3
14,S3
15,S4
16,S4
*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
NALL,2,2
TOP,1,2
*STEP,INC=20
*STATIC
0.1,1
*BOUNDARY
10,1,1,-6.468E-3
11,1,1,-9.056E-3
12,1,1,-12.68E-3
13,1,1,-17.75E-3
14,1,1,-24.85E-3
15,1,1,-17.75E-3
16,1,1,-12.68E-3
17,1,1,-9.056E-3
18,1,1,-6.468E-3
*NODE PRINT,NSET=TOP
U, 
*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,61,800
9,69,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,61,900
9,69,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