*HEADING
ECA3VFDM: *DSFLUX,*DSECURRENT FOR COUPLES THERMAL-ELECTRICAL ELEMENTS
 [DCAX3E]
*NODE,NSET=ALL
 1, 1.
 2, 8.
 3, 8.,7.
 4, 1.,7.
*NSET,NSET=FIX1
 1,2
*NSET,NSET=FIX2
 2,3
*NSET,NSET=FIX3
 3,4
*NSET,NSET=FIX4
 1,4
*ELEMENT,TYPE=DCAX3E,ELSET=ALLELS
 1, 1,2,3
 2, 1,3,4
*surface, NAME=SIDE_1
 1, S1
*surface, NAME=SIDE_2
 1, S2
*surface, NAME=SIDE_3
 2, S2
*surface, NAME=SIDE_4
 2, S3
*SOLID SECTION,MATERIAL=A1, ELSET=ALLELS
*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
*STEP
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*DSFLUX,OP=NEW
*DSECURRENT, OP=NEW
*NODE PRINT
 NT,RFL
 EPOT,RECUR
*OUTPUT,FIELD
*NODE OUTPUT
 NT,RFL
 EPOT,RECUR
*EL PRINT,FREQ=0
*OUTPUT,FIELD
*ELEMENT OUTPUT
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL,STEADY STATE
*BOUNDARY,OP=NEW
 FIX3, 11
 FIX3,  9
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL,STEADY STATE
*BOUNDARY,OP=NEW
 FIX4, 11
 FIX4,  9
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL,STEADY STATE
*BOUNDARY,OP=NEW
 FIX2, 11
 FIX2,  9
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY,OP=NEW
 FIX1, 11
 FIX1,  9
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY,OP=NEW
 FIX3, 11
 FIX3,  9
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY,OP=NEW
 FIX4, 11
 FIX4,  9
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY,OP=NEW
 FIX2, 11
 FIX2,  9
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY,OP=NEW
 FIX3,11
 FIX3, 9
*DSFLUX,OP=NEW
 SIDE_1, S, .3
*DSECURRENT, OP=NEW
 SIDE_1, CS, .3
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY,OP=NEW
 FIX4,11
 FIX4, 9
*DSFLUX,OP=NEW
 SIDE_2, S, .3
*DSECURRENT,OP=NEW
 SIDE_2, CS, .3
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY,OP=NEW
 FIX1,11
 FIX1, 9
*DSFLUX,OP=NEW
 SIDE_3, S, .3
*DSECURRENT, OP=NEW
 SIDE_3, CS, .3
*END STEP
*STEP
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY,OP=NEW
 FIX2,11
 FIX2, 9
*DSFLUX,OP=NEW
 SIDE_4, S, .3
*DSECURRENT, OP=NEW
 SIDE_4, CS,.3
*END STEP