*HEADING
SFILM, DSFLUX, SRADIATE, OP=NEW  FOR COUPLED TEMP-DISPLACEMENT-PORE PRESSURE ELEMENTS [CPE4PHT]
LOAD LABELS USED: F, S, R
Results should match with heat load steps of ece4tfdk.inp (coupled temperature-displacement plane strain elements)
**###################################################################
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE
1, 1., 0.
2, 8., 0.
3, 8., 7.
4, 1., 7.
*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=CPE4PHT, ELSET=EALL
1, 1,2,3,4
*SOLID SECTION,MATERIAL=ELMAT, ELSET=EALL
*MATERIAL,NAME=ELMAT
**#################
** Properties for the solid
**#################
*ELASTIC
30.E6, .3
*EXPANSION
0., 
*CONDUCTIVITY
3.77E-5, 
*DENSITY
82.9, 
*SPECIFIC HEAT
0.39,
*PERMEABILITY, SPECIFIC=1.0
1.0E-05
**####################
** Properties for the pore fluid
**####################
*EXPANSION, PORE FLUID
0., 
*CONDUCTIVITY, PORE FLUID
3.77E-5, 
*DENSITY, PORE FLUID
82.9, 
*SPECIFIC HEAT, PORE FLUID
.39, 
**####################
** Initial Conditions
**####################
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,0.
*INITIAL CONDITIONS,TYPE=RATIO
NALL,1.0
*INITIAL CONDITIONS,TYPE=PORE
NALL,0.
**
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-460., STEFAN BOLTZMANN = 5.E-8
**
** #####################################################################################
** Loading (OP=NEW)
** STEPS 1 to 4:   *DSFLUX    ->  Surface heat flux; heat surface flux per unit area specified on an element surface
** STEPS 5 to 8:   *SFILM     ->  Surface film condition; film coefficient and sink temperature provided on an element surface
** STEPS 9 to 12:  *SRADIATE  ->  Surface radiation; emissivity and sink temperature specified on an element surface
** #####################################################################################
**
**
**############
**   STEP 1
**############
*STEP
Surface flux on face 1
*SOILS
1., 1.
*DSFLUX, OP=NEW
SURF1, S, 0.3
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE3, 11
NALL, 8, 8, 0.0
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*EL PRINT
LOADS,FLUXS
FILM,RAD
CE, 
*NODE PRINT
RF, 
NT, 
RFL, 
*OUTPUT,FIELD
*NODE OUTPUT
RF,
NT,
RFL,
*END STEP
**############
**   STEP 2
**############
*STEP
Surface flux on face 2
*SOILS
1., 1.
*DSFLUX, OP=NEW
SURF2, S, 0.3
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE4, 11
NALL, 8, 8, 0.0
*END STEP
**############
**   STEP 3
**############
*STEP
Surface flux on face 3
*SOILS
1., 1.
*DSFLUX, OP=NEW
SURF3, S, 0.3
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE1, 11
NALL, 8, 8, 0.0
*END STEP
**############
**   STEP 4
**############
*STEP
Surface flux on face 4
*SOILS
1., 1.
*DSFLUX, OP=NEW
SURF4, S, 0.3
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE2, 11
NALL, 8, 8, 0.0
*END STEP
**############
**   STEP 5
**############
*STEP
SFILM face 1
*SOILS
1., 1.
*DSFLUX, OP=NEW
*SFILM, OP=NEW
SURF1, F, 75., .103
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE3, 11
NALL, 8, 8, 0.0
*END STEP
**############
**   STEP 6
**############
*STEP
SFILM face 2
*SOILS
1., 1.
*SFILM, OP=NEW
SURF2, F, 75., .103
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE4, 11
NALL, 8, 8, 0.0
*END STEP
**############
**   STEP 7
**############
*STEP
SFILM face 3
*SOILS
1., 1.
*SFILM, OP=NEW
SURF3, F, 75., .103
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE1, 11
NALL, 8, 8, 0.0
*END STEP
**############
**   STEP 8
**############
*STEP
SFILM face 4
*SOILS
1., 1.
*SFILM, OP=NEW
SURF4, F, 75., .103
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE2, 11
NALL, 8, 8, 0.0
*END STEP
**############
**   STEP 9
**############
*STEP
Surface radiation on face 1
*SOILS
1., 1.
*SFILM, OP=NEW
*SRADIATE, OP=NEW
SURF1, R, 75., 1.E-6
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE3, 11
NALL, 8, 8, 0.0
*END STEP
**############
**   STEP 10
**############
*STEP
Surface radiation on face 2
*SOILS
1., 1.
*SRADIATE, OP=NEW
SURF2, R, 75., 1.E-6
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE4, 11
NALL, 8, 8, 0.0
*END STEP
**############
**   STEP 11
**############
*STEP
Surface radiation on face 3
*SOILS
1., 1.
*SRADIATE, OP=NEW
SURF3, R, 75., 1.E-6
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE1, 11
NALL, 8, 8, 0.0
*END STEP
**############
**   STEP 12
**############
*STEP
Surface radiation on face 4
*SOILS
1., 1.
*SRADIATE, OP=NEW
SURF4, R, 75., 1.E-6
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE2, 11
NALL, 8, 8, 0.0
*END STEP