*HEADING
EC23VFDM: *DSFLUX,*DSECURRENT FOR COUPLED THERMAL-ELECTRICAL
ELEMENTS  [DC2D3E]
*RESTART,WRITE
*NODE,NSET=ALL1
 1, 0.
 2, 7.
 3, 7., 7.
 4, 0., 7.
*NCOPY, CHANGE NUMBER=10, OLD SET=ALL1, NEW SET=ALL2, SHIFT
8.0, 0.0, 0.0

*NSET,NSET=FIX1
 1,2
11,12
*NSET,NSET=FIX2
 2,3
12,13
*NSET,NSET=FIX3
 3,4
13,14
*NSET,NSET=FIX4
 1,4
11,14
*ELEMENT,TYPE=DC2D3E, ELSET=ALL1
 1,1,2,3
 2,1,3,4
*ELEMENT, TYPE=DC2D3E, ELSET=ALL2
11, 11,12,13
12, 11,13,14
*ELSET,ELSET=ELEMENT-1
1,11
*ELSET,ELSET=ELEMENT-2
2,12
*surface, NAME=SIDE_1A
 1, S1
*surface, NAME=SIDE_2A
 1, S2
*surface, NAME=SIDE_3A
 2, S2
*surface, NAME=SIDE_4A
 2, S3
*surface, NAME=SIDE_1B
 11, S1
*surface, NAME=SIDE_2B
 11, S2
*surface, NAME=SIDE_3B
 12, S2
*surface, NAME=SIDE_4B
 12, S3
*ORIENTATION, NAME=ORIW
1.0,0.0, 0.0, 0.0, 1.0, 0.0
3, 
*SOLID SECTION,MATERIAL=A1, ELSET=ALL1
*SOLID SECTION, MATERIAL=A1, ELSET=ALL2, ORIENTATION=ORIW
*surface,NAME=SIDE_SURFACE
2,S2
*MATERIAL,NAME=A1
*CONDUCTIVITY
 3.77E-5, 
*ELECTRICAL CONDUCTIVITY
 3.77E-5, 
*JOULE HEAT FRACTION
 0.0, 
*SPECIFIC HEAT
 .39, 
*DENSITY
 82.9, 
*BOUNDARY
 FIX1, 11
 FIX1, 9
*AMPLITUDE,NAME=RAMP-FUNCTION,TIME=STEP
0.,0.,1.,1.
*STEP
*COUPLED THERMAL-ELECTRICAL,STEADY STATE
*NODE PRINT
 NT,RFL
EPOT, RECUR
*OUTPUT,FIELD
*NODE OUTPUT
 NT,RFL
EPOT, RECUR
*EL PRINT,FREQ=0
*OUTPUT,FIELD
*ELEMENT OUTPUT
*OUTPUT,FIELD
*NODE OUTPUT
 NT,RFL
EPOT, RECUR
*DSFLUX,OP=NEW
 SIDE_1A, S, .3
 SIDE_1B, S, .3
*DSECURRENT,OP=NEW
 SIDE_1A, CS, .3
 SIDE_1B, CS, .3
*SECTION PRINT, NAME=SECTION-ID34,SURFACE=SIDE_SURFACE
SOAREA,SOH,SOE
*SECTION FILE,  NAME=SECTION-ID34,SURFACE=SIDE_SURFACE
SOAREA,SOH,SOE
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY,OP=NEW
 FIX4,11
 FIX4, 9
*DSFLUX,OP=NEW
 SIDE_2A, S, .3
 SIDE_2B, S, .3
*DSECURRENT,OP=NEW
 SIDE_2A, CS, .3
 SIDE_2B, CS, .3
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY,OP=NEW
 FIX1,11
 FIX1, 9
*DSFLUX,OP=NEW
 SIDE_3A, S, .3
 SIDE_3B, S, .3
*DSECURRENT,OP=NEW
 SIDE_3A, CS, .3
 SIDE_3B, CS, .3
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY,OP=NEW
 FIX2,11
 FIX2, 9
*DSFLUX,OP=NEW,AMP=RAMP-FUNCTION
 SIDE_4A, S, .3
 SIDE_4B, S, .3
*DSECURRENT, OP=NEW,AMP=RAMP-FUNCTION
 SIDE_4A, CS, .3
 SIDE_4B, CS, .3
*END STEP