*HEADING
EIA2SRS1_PPEN: SURFACE-BASED PRESSURE PENETRATION TEST 
DEFORMABLE SURFACE (SAX1) vs. DEFORMABLE 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
21,1.0,12.0
22,1.0,10.5
23,1.0,9.0
24,1.0,7.5
25,1.0,6.0
26,1.0,4.5
27,1.0,3.0
28,1.0,1.5
29,1.0,0
*NSET,NSET=BOT
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=SAX1,ELSET=TOP
9,21,22
10,22,23
11,24,23
12,25,24
13,25,26
14,26,27
15,28,27
16,29,28
*NSET,NSET=TOP
21,22,23,24,25,26,27,28,29
*SHELL SECTION,ELSET=TOP,MATERIAL=MAT
1.0,3
*SOLID SECTION,ELSET=BOT,MATERIAL=MAT
*MATERIAL,NAME=MAT
*ELASTIC
1.E5,0.3
*DENSITY
1.0E-3, 
*surface,NAME=TOPS
9,SPOS
10,SPOS
11,SNEG
12,SNEG
13,SPOS
14,SPOS
15,SNEG
16,SNEG
*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
TOPS,BOTS
*BOUNDARY
BOT,1,2
TOP,2,2
*STEP,INC=20
*STATIC
0.1,1
*DLOAD
9,P,800
10,P,1066.7
11,P,-1333.3
12,P,-1600
13,P,1600
14,P,1333.3
15,P,-1066.7
16,P,-800
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=TOPS,MASTER=BOTS
CSTRESS,
*CONTACT PRINT,SLAVE=TOPS,MASTER=BOTS
CSTRESS,
*ENDSTEP
*STEP, INC=20
*STATIC
0.01,1
*PRESSURE PENETRATION,SLAVE=TOPS,MASTER=BOTS
21,1,800
29,9,800
*CONTACT PRINT,SLAVE=TOPS,MASTER=BOTS,TOTAL=YES,FREQ=1
CSTRESS,PPRESS,CDISP
*OUTPUT,FIELD,FREQ=1
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=TOPS,MASTER=BOTS
CSTRESS,PPRESS,CDISP
*PRINT,CONTACT=YES
*ENDSTEP
*STEP, INC=20
*STATIC
0.01,1
*PRESSURE PENETRATION,SLAVE=TOPS,MASTER=BOTS
21,1,900
29,9,900
*CONTACT PRINT,SLAVE=TOPS,MASTER=BOTS,TOTAL=YES,FREQ=1
CSTRESS,PPRESS,CDISP
*OUTPUT,FIELD,FREQ=1
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=TOPS,MASTER=BOTS
CSTRESS,PPRESS,CDISP
*PRINT,CONTACT=YES
*ENDSTEP