*HEADING
DSLOAD, OP=NEW  FOR COUPLED TEMP-DISPLACEMENT-PORE PRESSURE ELEMENTS [CPE4RPT]
LOAD LABELS USED: P, HP
Results should match with mechanical load steps of ece4pfdl.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
*SURFACE, TYPE=ELEMENT, NAME=SURF1
1, S1
*SURFACE, TYPE=ELEMENT, NAME=SURF2
1, S2
*SURFACE, TYPE=ELEMENT, NAME=SURF3
1, S3
*SURFACE, TYPE=ELEMENT, NAME=SURF4
1, S4
*ELEMENT,TYPE=CPE4RPT, 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 4:  *DSLOAD  ->  Pressure loading on a surface
** STEPS 5 TO 8:  *DLOAD   ->  Hydrostatic pressure loading on a surface
** ####################################################
**############
**   STEP 1
**############
*STEP
Pressure on bottom face
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
SURF1,P,100.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*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
Pressure on right face
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
SURF2,P,100.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 3
**############
*STEP
Pressure on top face
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
SURF3,P,100.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 4
**############
*STEP
Pressure on left face
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
SURF4,P,100.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 5
**############
*STEP
Hydrostatic Pressure on bottom face
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
SURF1,HP,100.,5.,0.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 6
**############
*STEP
Hydrostatic Pressure on right face
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
SURF2,HP,100.,5.,0.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 7
**############
*STEP
Hydrostatic Pressure on top face
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
SURF3,HP,100.,5.,0.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
**############
**   STEP 8
**############
*STEP
Hydrostatic Pressure on left face
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
SURF4,HP,100.,5.,0.
*BOUNDARY
NALL,1,2
NALL,8,8
NALL, 11, 11
*END STEP
***************