*HEADING
: PORE PRESSURE COUPLED CONTACT (C3D8P) (NO SLIDING)
*RESTART,WRITE
*NODE,nset=nall
       1,          1.,        12.5
       2,          2.,        12.5
       3,          2.,         15.
       4,          1.,         15.
       5,          1.,        12.5,          1.
       6,          2.,        12.5,          1.
       7,          2.,         15.,          1.
       8,          1.,         15.,          1.
**
    1001,          1.,        10.0000001,
    1002,          2.,        10.0000001
    1003,          2.,        12.5
    1004,          1.,        12.5
    1005,          1.,        10.0000001,          1.,
    1006,          2.,        10.0000001,          1.,
    1007,          2.,        12.5,          1.
    1008,          1.,        12.5,          1.
**
*NSET,NSET=BOTTOM
1001,1002,1005,1006
*NSET,NSET=TOP
3,4,7,8
*ELEMENT, TYPE=C3D8P, ELSET=SOLID
1, 1, 2, 3, 4, 5, 6, 7, 8
2, 1001, 1002, 1003, 1004, 1005, 1006, 1007, 1008
**2, 9, 10, 20, 19,  13,  14, 24, 23
*SOLID SECTION, ELSET=SOLID, MATERIAL=SOLMAT
*MATERIAL, NAME=SOLMAT
*ELASTIC
30.E6,0.0
*PERMEABILITY,SPECIFIC=1.0
 1.E-1, 
*SURFACE INTERACTION, NAME=GAP
1.0, 
** surfaces
**
*surface, NAME=ASURF
1, S3
*surface, NAME=BSURF
2, S5
*CONTACT PAIR, INTERACTION=GAP
BSURF, ASURF
**
*INITIAL CONDITIONS,TYPE=RATIO
NALL,1.
**
*STEP
*SOILS, CONSOLIDATION
 0.1 , 1.
*BOUNDARY
BOTTOM,1,3
TOP,1,1
TOP,3,3
BOTTOM,8,8,1
TOP,8,8,2
NALL,1,1
NALL,3,3
*DLOAD
1,P5,10
*CONTACT PRINT,F=10
PFL,PFLA,PTL,PTLA
TPFL,TPTL
**
*CONTACT FILE,F=10
CAREA,
PFL, PFLA, PTL, PTLA
TPFL,TPTL
*NODE PRINT,F=10
POR,
*NODE FILE,F=10
POR,
*EL PRINT, F=0
*END STEP