*HEADING
: SURFACE-BASED PRESSURE PENETRATION TEST 
WITH RIGID SURFACE, SAX1
*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
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=NALL
1,2,3,4,5,6,7,8,9
*ELEMENT,TYPE=SAX1,ELSET=BOT
1,1,2
2,2,3
3,4,3
4,5,4
5,5,6
6,6,7
7,8,7
8,9,8
*ELEMENT,TYPE=RAX2,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
1,2,3,4,5,6,7,8,9
*RIGID BODY,ELSET=TOP,REF NODE=100
1.0, 
*SHELL GENERAL SECTION,ELSET=BOT,MATERIAL=MAT
1.0, 
*MATERIAL,NAME=MAT
*ELASTIC
1.E5,0.3
*DENSITY
1.0E-3, 
*surface,NAME=TOPS
TOP,SPOS
*surface,NAME=BOTS
1,SNEG
2,SNEG
3,SPOS
4,SPOS
5,SNEG
6,SNEG
7,SPOS
8,SPOS
*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
*AMPLITUDE,NAME=AMPP
0,0,1,1
*STEP,INC=20
*STATIC
0.1,1
*DLOAD
1,P,-600
2,P,-800
3,P,1000
4,P,1200
5,P,-1200
6,P,-1000
7,P,800
8,P,600
*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,AMP=AMPP
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,OP=MOD
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