*HEADING
 DLOAD, FLOW, DFLOW and CFOW OR PORE PRESSURE ELEMENTS  [C3D4P]
*PREPRINT,model=yes
*NODE,NSET=NALL
  1, 1.0000, 0.0000,  0.0000                                    
  3, 4.0000, 0.0000,  0.0000                                    
  7, 1.0000, 0.0000, -1.0000                                    
  9, 4.0000, 0.0000, -1.0000                                    
 19, 1.0000, 5.0000,  0.0000                                    
 21, 4.0000, 5.0000,  0.0000                                    
 25, 1.0000, 5.0000, -1.0000                                    
 27, 4.0000, 5.0000, -1.0000                                    
*NSET,NSET=NS1
 1 , 7  , 25
*NSET,NSET=NS2
 1 , 3  , 9 , 7
*NSET,NSET=NS3
 1 , 3  , 21 , 19
*NSET,NSET=CORNERS1
 1, 3, 7, 9
*NSET,NSET=CORNERS2
 25, 27, 19, 21
*NSET,NSET=CORNERS3
 19, 21, 1, 3
*NSET,NSET=CORNERS4
 21, 27, 3, 9
*NSET,NSET=CORNERS5
 27, 25, 7, 9
*NSET,NSET=CORNERS6
 25, 19, 1, 7
*NSET,NSET=CORNERS
 CORNERS1,
 CORNERS2,
 CORNERS3,
 CORNERS4,
 CORNERS5,
 CORNERS6,
*ELEMENT,TYPE=C3D6P,ELSET=ELEMS
 1,  1, 3, 9,19,21,27
 2,  9, 7, 1,27,25,19
*surface, NAME=SIDE_1
 1, S1
 2, S1
*surface, NAME=SIDE_2
 1, S2
 2, S2
*surface, NAME=SIDE_3
 1, S3
*surface, NAME=SIDE_4
 1, S4
*surface, NAME=SIDE_5
 2, S3
*surface, NAME=SIDE_6
 2, S4
*SOLID SECTION,MATERIAL=A1, ELSET=ELEMS
*MATERIAL,NAME=A1
*ELASTIC
 1.E8, 
*PERMEABILITY,SPECIFIC=1.0
 1.E-5, 
*DENSITY
 1.7321, 
*AMPLITUDE,NAME=STEP-FUNCTION,TIME=STEP
 0.,1.,1.,1.
*INITIAL CONDITIONS,TYPE=RATIO
 NALL,1.
*STEP
STEP 1
*SOILS,CONSOLIDATION
 1. , 1.
*EL PRINT
 LOADS,
*NODE PRINT,TOTALS=YES
 RF, 
*NODE PRINT
 POR, 
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS, 8
*DSLOAD, OP=NEW
 SIDE_1, P, 100.
 SIDE_2, P, 100.
*output, field
*element output
 loads,s
*node output
 rf,por
*END STEP
*STEP
STEP 2
*SOILS,CONSOLIDATION
 1., 1.
*DSLOAD, OP=NEW
 SIDE_3, P, 100.
 SIDE_5, P, 100.
*END STEP
*STEP
STEP 3
*SOILS,CONSOLIDATION
 1., 1.
*DSLOAD, OP=NEW
 SIDE_4, P, 100.
 SIDE_6, P, 100.
*END STEP
*STEP
STEP 4
*SOILS,CONSOLIDATION
 1., 1.
*DSLOAD, OP=NEW
 SIDE_1, HP, 100., 5., 0.
 SIDE_2, HP, 100., 5., 0.
*END STEP
*STEP
STEP 5
*SOILS,CONSOLIDATION
 1., 1.
*DSLOAD, OP=NEW
 SIDE_3, HP, 300., 5., 0.
 SIDE_5, HP, 100., 5., 0.
*END STEP
*STEP
STEP 6
*SOILS,CONSOLIDATION
 1., 1.
*DSLOAD, OP=NEW
 SIDE_4, HP, 100., 5., 0.
 SIDE_6, HP, 100., 5., 0.
*END STEP
*STEP
STEP 7
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS2, 8
*DSLOAD, OP=NEW
*SFLOW, OP=NEW
*END STEP
*STEP
STEP 8
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY,OP=NEW
 NALL, 1,3
 CORNERS1, 8
*SFLOW, OP=NEW, AMP=STEP-FUNCTION
 SIDE_2, Q, 14.7, 2.E-5
*END STEP
*STEP
STEP 9
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS5, 8
*SFLOW, OP=NEW
 SIDE_3, Q, 14.7, 2.E-5
*END STEP
*STEP
STEP 10
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS6, 8
*SFLOW, OP=NEW
 SIDE_4, Q, 14.7, 2.E-5
*END STEP
*STEP
STEP 11
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS3, 8
*SFLOW, OP=NEW
 SIDE_5, Q, 14.7, 2.E-5
*END STEP
*STEP
STEP 12
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS4, 8
*SFLOW, OP=NEW
 SIDE_6, Q, 14.7, 2.E-5
*END STEP
*STEP
STEP 13
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS2, 8
*SFLOW, OP=NEW
*DSFLOW, OP=NEW
*END STEP
*STEP
STEP 14
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS1, 8
*DSFLOW, OP=NEW, AMP=STEP-FUNCTION
 SIDE_2, S, 3.E-5
*END STEP
*STEP
STEP 15
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS5, 8
*DSFLOW, OP=NEW
 SIDE_3, S, 3.E-5
*END STEP
*STEP
STEP 16
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS6, 8
*DSFLOW, OP=NEW
 SIDE_4, S, 3.E-5
*END STEP
*STEP
STEP 17
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS3, 8
*DSFLOW, OP=NEW
 SIDE_5, S, 3.E-5
*END STEP
*STEP
STEP 18
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS4, 8
*DSFLOW, OP=NEW
 SIDE_6, S, 3.E-5
*END STEP