*HEADING
DLOAD, OP=NEW  FOR COUPLED TEMP-DISPLACEMENT-PORE PRESSURE ELEMENTS [CPE4RPT]
LOAD LABELS USED: BX, BY, GRAV, CENTRIF, P1, P2, P3, P4, HP1, HP2, HP3, HP4
Results should match with mechanical 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
**####################
** Initial Conditions
**####################
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL, 0.
*INITIAL CONDITIONS,TYPE=RATIO
NALL,1.
*INITIAL CONDITIONS,TYPE=PORE
NALL,0.
** ####################################################
** Loading (OP=NEW)
** STEPS 1 TO 2  :  *DLOAD  ->  Body force
** STEP 3        :  *DLOAD  ->  Gravity loading
** STEP 4        :  *DLOAD  ->  Centrifugal loading
** STEPS 5 TO 8  :  *DLOAD  ->  Pressure loading on element face n
** STEPS 10 TO 13:  *DLOAD  ->  Hydrostatic pressure loading on element face n
** ####################################################
**############
**   STEP 1
**############
*STEP
Body force in X direction
*SOILS
1.0,1.0
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*DLOAD, OP=NEW
ALL_ELE,BX,100.
*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
Body force in Y direction
*SOILS
1.0,1.0
*DLOAD, OP=NEW
ALL_ELE,BY,100.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 3
**############
*STEP
Gravity loading
*SOILS
1.0,1.0
*DLOAD, OP=NEW
ALL_ELE,GRAV,100.0, 1.0,1.0,0.0
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 4
**############
*STEP
Centrifugal loading
*SOILS
1.0,1.0
*DLOAD, OP=NEW
ALL_ELE,CENTRIF,100., 0.0,2.5,0., 0.,0.,1.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 5
**############
*STEP
Pressure on bottom face
*SOILS
1.0,1.0
*DLOAD, OP=NEW
1,P1,100.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 6
**############
*STEP
Pressure on right face
*SOILS
1.0,1.0
*DLOAD, OP=NEW
1,P2,100.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 7
**############
*STEP
Pressure on top face
*SOILS
1.0,1.0
*DLOAD, OP=NEW
1,P3,100.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 8
**############
*STEP
Pressure on left face
*SOILS
1.0,1.0
*DLOAD, OP=NEW
1,P4,100.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 9
**############
*STEP
Hydrostatic Pressure on bottom face
*SOILS
1.0,1.0
*DLOAD, OP=NEW
1,HP1,100.,5.,0.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 10
**############
*STEP
Hydrostatic Pressure on right face
*SOILS
1.0,1.0
*DLOAD, OP=NEW
1,HP2,100.,5.,0.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 11
**############
*STEP
Hydrostatic Pressure on top face
*SOILS
1.0,1.0
*DLOAD, OP=NEW
1,HP3,100.,5.,0.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 12
**############
*STEP
Hydrostatic Pressure on left face
*SOILS
1.0,1.0
*DLOAD, OP=NEW
1,HP4,100.,5.,0.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
***************