*HEADING
: PORE PRESSURE COUPLED CONTACT (C3D8P) (FINITE SLIDING)
*RESTART,WRITE,F=2
*NODE,NSET=ALL
  1, 1., 0.,-0.01
  7, 7., 0.,-0.01
 11, 1., 1.,-0.01
 17, 7., 1.,-0.01
 21, 1., 0., 1.01
 27, 7., 0., 1.01
 31, 1., 1., 1.01
 37, 7., 1., 1.01
101, 3., 1.
105, 5., 1.
111, 3., 2.
115, 5., 2.
121, 3., 1., 1.0
125, 5., 1., 1.0
131, 3., 2., 1.0
135, 5., 2., 1.0
*NGEN,NSET=ALL
  1,  7
 11, 17
 21, 27
 31, 37
101,105
111,115
121,125
131,135
*NSET,NSET=BOTTOM_NODES,GEN
1,7
21,27
*NSET,NSET=TOP_NODES,GEN
111,115
131,135
**
*ELEMENT,TYPE=C3D8PT,ELSET=DEAD
1,1,3,13,11,21,23,33,31
2,3,5,15,13,23,25,35,33
3,5,7,17,15,25,27,37,35
*SOLID SECTION,ELSET=DEAD,MATERIAL=DEAD
*MATERIAL,NAME=DEAD
*ELASTIC
30.E6,
*PERMEABILITY,SPECIFIC=1.0
 1.E-4, 
*CONDUCTIVITY
11.9
*SPECIFIC HEAT
40.
*DENSITY
1000.
**
*ELEMENT,TYPE=C3D8PT,ELSET=ACTIVE
101,101,103,113,111,121,123,133,131
102,103,105,115,113,123,125,135,133
*SOLID SECTION,ELSET=ACTIVE,MATERIAL=ACTIVE
*MATERIAL,NAME=ACTIVE
*ELASTIC
30.E6,
*PERMEABILITY,SPECIFIC=1.0
 1.E-4, 
*CONDUCTIVITY
11.9
*SPECIFIC HEAT
40.
*DENSITY
1000.
*surface,NAME=MAS
1,S5
2,S5
3,S5
*surface,NAME=SLA
101,S3
102,S3
*CONTACT PAIR,INTERACTION=ISLS
SLA, MAS
*SURFACE INTERACTION,NAME=ISLS
1.0,
**
*INITIAL CONDITIONS,TYPE=RATIO
ALL, 1.0
*INITIAL CONDITIONS,TYPE=PORE PRESSURE
ALL, 100.0
**INITIAL CONDITIONS,TYPE=TEMPERATURE
**ALL, 0.0
**
*STEP,INC=10,NLGEOM,AMP=STEP
Step 1
*SOIL, CONSOLIDATION
1,1.
*BOUNDARY
BOTTOM_NODES,2,3
TOP_NODES,3,3
1,1,1
7,1,1
21,1,1
27,1,1
113,1,1
133,1,1
*DLOAD
101,P5,100.
102,P5,100.
*DFLOW
101,S5,3.e-4
102,S5,3.e-4
1,S3,-1.e-4
2,S3,-1.e-4
3,S3,-1.e-4
**
*EL PRINT,F=100
S, 
E, 
*PRINT,CONTACT=YES
*CONTACT PRINT
CAREA, TPFL
*CONTACT FILE
CAREA, TPFL
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP
Step 2
*SOIL, CONSOLIDATION
1,1.
*BOUNDARY,FIXED
TOP_NODES,2,2
*BOUNDARY
113,1,1, 1.5
133,1,1, 1.5
*CONTACT PRINT
CAREA, TPFL
*CONTACT FILE
CAREA, TPFL
*END STEP