*HEADING
FLOW, DFLOW, OP=NEW  FOR COUPLED TEMP-DISPLACEMENT-PORE PRESSURE ELEMENTS [CPE4RPT]
LOAD LABELS USED: Q1, Q2, Q3, Q4, S1, S2, S3, S4
Results should match with flow load steps of ece4pfde.inp (plane strain elements with pore-pressure)
**###################################################################
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE
1, 1., 0.
2, 4., 0.
3, 4., 5.
4, 1., 5.
*NSET, NSET=SIDE1
1, 2
*NSET, NSET=SIDE2
2, 3
*NSET, NSET=SIDE3
3, 4
*NSET, NSET=SIDE4
4, 1
*NSET, NSET=NALL
1,2,3,4
*ELEMENT,TYPE=CPE4RPHT, ELSET=ALL_ELE
1, 1,2,3,4
*SOLID SECTION,MATERIAL=ELMAT, ELSET=ALL_ELE
*MATERIAL,NAME=ELMAT
**#################
** Properties for the solid
**#################
*ELASTIC
1.0E+08,0.0
*DENSITY
1.4142, 
*PERMEABILITY, SPECIFIC=1.0
1.0E-05
*CONDUCTIVITY
0.1
*SPECIFIC HEAT
0.39,
**####################
** Initial Conditions
**####################
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL, 0.
*INITIAL CONDITIONS,TYPE=RATIO
NALL,1.
*INITIAL CONDITIONS,TYPE=PORE
NALL,0.
**###################################################
** STEPS 1 TO 4: *FLOW	 ->  Specify seepage coefficient and sink pore pressure on element face to model pore fluid flow as a function of current pore pressure
** STEPS 5 TO 8: *DFLOW	 ->  Specify seepage velocity on element face
**###################################################
**############
**   STEP 1
**############
*STEP
FLOW through face 1
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY,OP=NEW
NALL, 1,2
SIDE3, 8
NALL, 11, 11
*FLOW, OP=NEW
ALL_ELE, Q1, 14.7, 2.E-5
*EL PRINT
LOADS,
*NODE PRINT,TOTALS=YES
RF, 
*NODE PRINT
POR,
*OUTPUT,FIELD
*NODE OUTPUT
RF,
POR
*OUTPUT,HISTORY
*ENERGY OUTPUT,VAR=PRESELECT
*END STEP
**############
**   STEP 2
**############
*STEP
FLOW through face 2
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY,OP=NEW
NALL, 1,2
SIDE4, 8
NALL, 11, 11
*FLOW, OP=NEW
ALL_ELE, Q2, 14.7, 2.E-5
*END STEP
**############
**   STEP 3
**############
*STEP
FLOW through face 3
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY,OP=NEW
NALL, 1,2
SIDE1, 8
NALL, 11, 11
*FLOW, OP=NEW
ALL_ELE, Q3, 14.7, 2.E-5
*END STEP
**############
**   STEP 4
**############
*STEP
FLOW through face 4
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY,OP=NEW
NALL, 1,2
SIDE2, 8
NALL, 11, 11
*FLOW, OP=NEW
ALL_ELE, Q4, 14.7, 2.E-5
*END STEP
**#######################################
**#######################################
**############
**   STEP 5
**############
*STEP
DFLOW through face 1
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE3, 8
NALL, 11, 11
*FLOW, OP=NEW
*DFLOW, OP=NEW
ALL_ELE, S1, 3.E-5
*END STEP
**############
**   STEP 6
**############
*STEP
DFLOW through face 2
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE4, 8
NALL, 11, 11
*DFLOW, OP=NEW
ALL_ELE, S2, 3.E-5
*END STEP
**############
**   STEP 7
**############
*STEP
DFLOW through face 3
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE1, 8
NALL, 11, 11
*DFLOW, OP=NEW
ALL_ELE, S3, 3.E-5
*END STEP
**############
**   STEP 8
**############
*STEP
DFLOW through face 4
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE2, 8
NALL, 11, 11
*DFLOW, OP=NEW
ALL_ELE, S4, 3.E-5
*END STEP