*HEADING
[] COUPLED PORE-DISPLACEMENT (C3D20P) (NO SLIDING)
*RESTART,WRITE
*NODE,nset=nall
       1,          1.,         12
       2,          2.,         12
       3,          2.,         14.
       4,          1.,         14.
       5,          1.,         12,          1.
       6,          2.,         12,          1.
       7,          2.,         14.,          1.
       8,          1.,         14.,          1.
**
       9,          1.5,        12
      10,          2.,         13
      11,          1.5 ,       14          
      12,          1.,         13
      13,          1.5,        12,         1          
      14,          2.,         13,         1
      15,          1.5 ,       14,         1     
      16,          1.,         13,         1
      17,          1.,         12,         0.5
      18,          2.,         12,         0.5
      19,          2.,         14.,        0.5
      20,          1.,         14.,        0.5
********
    1001,          1.,        10.0000001,
    1002,          2.,        10.0000001
    1003,          2.,        12
    1004,          1.,        12
    1005,          1.,        10.0000001,          1.,
    1006,          2.,        10.0000001,          1.,
    1007,          2.,        12,          1.
    1008,          1.,        12,          1.
**
      1009,          1.5,        10
      1010,          2.,         11
      1011,          1.5 ,       12          
      1012,          1.,         11
      1013,          1.5,        10,         1          
      1014,          2.,         11,         1
      1015,          1.5 ,       12,         1     
      1016,          1.,         11,         1
      1017,          1.,         10,         0.5
      1018,          2.,         10,         0.5
      1019,          2.,         12.,        0.5
      1020,          1.,         12.,        0.5
**
*NSET,NSET=BOTTOM
1001,1002,1005,1006,1017,1013,1018,1009
*NSET,NSET=TOP
3,4,7,8,11,15,19,20
*ELEMENT, TYPE=C3D20P, ELSET=SOLID
1,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20
2, 1001, 1002, 1003, 1004, 1005, 1006, 1007, 1008,
1009,1010,1011,1012,1013,1014,1015,1016,1017,1018,1019,1020
*SOLID SECTION, ELSET=SOLID, MATERIAL=SOLMAT
*MATERIAL, NAME=SOLMAT
*ELASTIC
30.E6,0.0
*PERMEABILITY,SPECIFIC=1.0
 1.E-4, 
*DENSITY
  1.4142, 
*SURFACE INTERACTION, NAME=GAP
*SURFACE BEHAVIOR, AUGMENTED
** surfaces
**
*surface, NAME=ASURF
1, S3
*surface, NAME=BSURF
2, S5
*CONTACT PAIR, INTERACTION=GAP,SUPPLEMENTARY CONSTRAINTS=SELECTIVE
ASURF, BSURF
**
*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