*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=C3D4P,ELSET=ELEMS
 1, 7, 9,25, 1
 2, 1, 9,21, 3
 3, 9,27,25,21
 4, 1,21,25,19
 5, 1,21, 9,25
*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
*SOILS,CONSOLIDATION
 1. , 1.
*EL PRINT
 LOADS,
*NODE PRINT,TOTALS=YES
 RF, 
*NODE PRINT
 POR, 
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS, 8
*DLOAD,OP=NEW
 ELEMS,CENTRIF,100.,0.,2.5,-.5,0.,0.,1.
*output, field
*element output
 s
*node output
 u,rf
*END STEP
*STEP
*SOILS, CONSOLIDATION
 1.,1.
*DLOAD, OP=NEW
 ELEMS, BX, 100.
 ELEMS, BY, 100.
 ELEMS, BZ, 100.
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*DLOAD, OP=NEW
 ELEMS, GRAV, 100.,1,1,1
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*DLOAD, OP=NEW
 1,P2,100.
 2,P2,100.
 4,P3,100.
 3,P3,100.
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*DLOAD, OP=NEW
 4,P2,100.
 2,P4,100.
 3,P1,100.
 1,P1,100.
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*DLOAD, OP=NEW
 2,P3,100.
 3,P2,100.
 1,P4,100.
 4,P4,100.
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*DLOAD, OP=NEW
 1,HP2, 100., 5., 0.
 2,HP2, 100., 5., 0.
 4,HP3, 100., 5., 0.
 3,HP3, 100., 5., 0.
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*DLOAD, OP=NEW
 4,HP2, 100., 5., 0.
 2,HP4, 100., 5., 0.
 3,HP1, 100., 5., 0.
 1,HP1, 100., 5., 0.
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*DLOAD, OP=NEW
 2,HP3, 100., 5., 0.
 3,HP2, 100., 5., 0.
 1,HP4, 100., 5., 0.
 4,HP4, 100., 5., 0.
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS2, 8
*DLOAD, OP=NEW
*FLOW, OP=NEW
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY,OP=NEW
 NALL, 1,3
 CORNERS1, 8
*FLOW, OP=NEW, AMP=STEP-FUNCTION
 4,Q3, 14.7, 2.E-5
 3,Q3, 14.7, 2.E-5
*output, field
*element output
 por,s,sat
*node output
 u,por,rf
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS5, 8
*FLOW, OP=MOD
 4,Q3, 14.7, 0.
 3,Q3, 14.7, 0.
 4,Q2, 14.7, 2.E-5
 2,Q4, 14.7, 2.E-5
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS6, 8
*FLOW, OP=NEW
 2,Q3, 14.7, 2.E-5
 3,Q2, 14.7, 2.E-5
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS3, 8
*FLOW, OP=NEW
 3,Q1, 14.7, 2.E-5
 1,Q1, 14.7, 2.E-5
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS4, 8
*FLOW, OP=NEW
 1,Q4, 14.7, 2.E-5
 4,Q4, 14.7, 2.E-5
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS2, 8
*FLOW, OP=NEW
*DFLOW, OP=NEW
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS1, 8
*CFLOW, OP=NEW, AMP=STEP-FUNCTION
 21,,2.250e-5
 27,,2.250e-5
 19,,2.250e-5
 25,,2.250e-5
*output, field
*element output
 s,por,sat
*node output
 u,rf,cff,por
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS5, 8
*CFLOW, OP=NEW
*DFLOW, OP=NEW
 4,S2, 3.E-5
 2,S4, 3.E-5
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS6, 8
*DFLOW, OP=NEW
 2,S3, 3.E-5
 3,S2, 3.E-5
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS3, 8
*DFLOW, OP=NEW
 3,S1, 3.E-5
 1,S1, 3.E-5
*END STEP
*STEP
*SOILS,CONSOLIDATION
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS4, 8
*DFLOW, OP=NEW
 1,S4, 3.E-5
 4,S4, 3.E-5
*END STEP