*HEADING
EC38VFDF: DFLUX, DECURRENT FOR COUPLED THERMAL-ELECTRICAL
ELEMENTS  [DC3D8E]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE,NSET=ALL
 1, 0., 0., 0.
 2, 7., 0., 0.
 3, 7., 0., -7.
 4, 0., 0., -7.
 5, 0., 7., 0.
 6, 7., 7., 0.
 7, 7., 7., -7.
 8, 0., 7., -7.
*NCOPY, CHANGE NUMBER=10, OLD SET=ALL, NEW SET=ALL2, SHIFT
8.0, 0.0, 0.0

*NCOPY, CHANGE NUMBER=10, OLD SET=ALL2, NEW SET=ALL3, SHIFT
8.0, 0.0, 0.0

*NCOPY, CHANGE NUMBER=10, OLD SET=ALL3, NEW SET=ALL4, SHIFT
8.0, 0.0, 0.0

*NODE,NSET=ALL5
41,37.      , 0.      , 0.
42,41.949747, 4.949747, 0.
43,41.949747, 4.949747,-7.
44,37.      , 0.      ,-7.
45,32.05026 , 4.949747, 0.
46,37.      , 9.899494, 0.
47,37.      , 9.899494,-7.
48,32.05026 , 4.949747,-7
*NSET,NSET=FIX1
 1,2,3,4
 11,12,13,14
 21,22,23,24
 31,32,33,34
 41,42,43,44
*NSET,NSET=FIX2
 5,6,7,8
 15,16,17,18
 25,26,27,28
 35,36,37,38
 45,46,47,48
*NSET,NSET=FIX3
 1,2,6,5
 11,12,16,15
 21,22,26,25
 31,32,36,35
 41,42,46,45
*NSET,NSET=FIX4
 2,3,7,6
 12,13,17,16
 22,23,27,26
 32,33,37,36
 42,43,47,46
*NSET,NSET=FIX5
 3,4,7,8
 13,14,17,18
 23,24,27,28
 33,34,37,38
 43,44,47,48
*NSET,NSET=FIX6
 1,4,5,8
 11,14,15,18
 21,24,25,28
 31,34,35,38
 41,44,45,48
*ELEMENT,TYPE=DC3D8E, ELSET=EALL
 1, 1,2,3,4,5,6,7,8
*ELEMENT, TYPE=DC3D8E, ELSET=ALL2
 11, 11,12,13,14,15,16,17,18
*ELEMENT, TYPE=DC3D8E, ELSET=EALL3
 21, 21,22,23,24,25,26,27,28
*ELEMENT, TYPE=DC3D8E, ELSET=EALL4
 31, 31,32,33,34,35,36,37,38
*ELEMENT, TYPE=DC3D8E, ELSET=EALL5
 41, 41,42,43,44,45,46,47,48
**local orientation matching global system 
*ORIENTATION,NAME=RECT
 1.0, 0.0, 0.0, 0.0, 1.0, 0.0
 1, 0.0
*ORIENTATION, NAME=ORIW
0.0,1.0, 0.0, -1.0, 0.0, 0.0
3, 
***ORIENTATION,SYSTEM=RECTANGULAR,NAME=RECT,DEFINITION=COORDINATES
**1.,1.,0.,-1.,1.,0.
3,
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*SOLID SECTION,MATERIAL=A1, ELSET=ALL2, ORIENTATION=ORIW
*SOLID SECTION,MATERIAL=ISO, ELSET=EALL3
*SOLID SECTION,MATERIAL=ORTH, ELSET=EALL4,ORIENT=RECT
***SOLID SECTION,MATERIAL=ANIS, ELSET=EALL5, ORIENTATION=RECT
*SOLID SECTION,MATERIAL=ANIS, ELSET=EALL5,ORIENT=RECT
*MATERIAL,NAME=A1
* ELECTRICAL CONDUCTIVITY
 3.77E-5, 
*CONDUCTIVITY
 3.77E-5, 
*JOULE HEAT FRACTION
 0.0, 
*DENSITY
 82.9, 
*SPECIFIC HEAT
 .39, 
*MATERIAL,NAME=ISO
* ELECTRICAL CONDUCTIVITY,TYPE=ISO
 3.77E-5, 
*CONDUCTIVITY
 3.77E-5, 
*JOULE HEAT FRACTION
 0.0, 
*DENSITY
 82.9, 
*SPECIFIC HEAT
 .39, 
*MATERIAL,NAME=ANIS
*****
**USED k^0*[M] to define anisotropic values
**     where k^0 = 3.77E-5
**           [M] = |3, -1, 0|
**                 |-1, 3, 0|
**                 |0,  0, 1|
***** 
* ELECTRICAL CONDUCTIVITY,TYPE=ANISO
 11.31E-5,-3.77E-5,11.31E-5,0.,0.,3.77E-5
*CONDUCTIVITY,TYPE=ANISO
 11.31E-5,-3.77E-5,11.31E-5,0.,0.,3.77E-5
*JOULE HEAT FRACTION
 0.0, 
*DENSITY
 82.9, 
*SPECIFIC HEAT
 .39, 
*MATERIAL,NAME=ORTH
*****
**USED k^0*[M] to define anisotropic values
**     where k^0 = 3.77E-5
**           [M] = |2, 0, 0|
**                 |0, 4, 0|
**                 |0, 0, 1|
** Thus, the results obtained for ANISO must equal results obtained for ORTHO
***** 
* ELECTRICAL CONDUCTIVITY,TYPE=ORTHO
 7.54E-5,15.08E-5,3.77E-5
*CONDUCTIVITY,TYPE=ORTHO
 7.54E-5,15.08E-5,3.77E-5
*JOULE HEAT FRACTION
 0.0, 
*DENSITY
 82.9, 
*SPECIFIC HEAT
 .39, 
*BOUNDARY
 FIX1,11
 FIX1, 9
*RESTART, WRITE
*STEP,UNSYMM=YES
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*SOLUTION TECHNIQUE,TYPE=SEPARATED
*NODE PRINT
 NT,RFL
EPOT, RECUR
**EL PRINT,FREQ=0
*EL PRINT
EPG,
ECD, 
*DFLUX
 1, BF, .3
 11, BF, .3
 21, BF, .3
 31, BF, .3
 41, BF, .3
*DECURRENT,OP=NEW
 1, CBF, .3
 11, CBF, .3
 21, CBF, .3
 31, CBF, .3
 41, CBF, .3
*END STEP
*STEP,UNSYMM=YES
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY, OP=NEW
 FIX2, 11
 FIX2,  9
*DFLUX, OP=NEW
 1, S1, .3
 11, S1, .3
 21, S1, .3
 31, S1, .3
 41, S1, .3
*DECURRENT,OP=NEW
 1, CS1, .3
 11, CS1, .3
 21, CS1, .3
 31, CS1, .3
 41, CS1, .3
*END STEP
*STEP,UNSYMM=YES
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY, OP=NEW
 FIX1, 11
 FIX1,  9
*DFLUX,OP=NEW
 1, S2, .3
 11, S2, .3
 21, S2, .3
 31, S2, .3
 41, S2, .3
*DECURRENT,OP=NEW
 1, CS2, .3
 11, CS2, .3
 21, CS2, .3
 31, CS2, .3
 41, CS2, .3
*END STEP
*STEP,UNSYMM=YES
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY, OP=NEW
 FIX5, 11
 FIX5,  9
*DFLUX,OP=NEW
 1, S3, .3
 11, S3, .3
 21, S3, .3
 31, S3, .3
 41, S3, .3
*DECURRENT,OP=NEW
 1, CS3, .3
 11, CS3, .3
 21, CS3, .3
 31, CS3, .3
 41, CS3, .3
*END STEP
*STEP,UNSYMM=YES
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY, OP=NEW
 FIX6, 11
 FIX6,  9
*DFLUX,OP=NEW
 1, S4, .3
 11, S4, .3
 21, S4, .3
 31, S4, .3
 41, S4, .3
*DECURRENT, OP=NEW
 1, CS4, .3
 11, CS4, .3
 21, CS4, .3
 31, CS4, .3
 41, CS4, .3
*END STEP
*STEP,UNSYMM=YES
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY, OP=NEW
 FIX3, 11
 FIX3,  9
*DFLUX,OP=NEW
 1, S5, .3
 11, S5, .3
 21, S5, .3
 31, S5, .3
 41, S5, .3
*DECURRENT, OP=NEW
 1, CS5, .3
 11, CS5, .3
 21, CS5, .3
 31, CS5, .3
 41, CS5, .3
*END STEP
*STEP,UNSYMM=YES
*COUPLED THERMAL-ELECTRICAL, STEADY STATE
*BOUNDARY, OP=NEW
 FIX4, 11
 FIX4, 9
*DFLUX,OP=NEW
 1, S6, .3
 11, S6, .3
 21, S6, .3
 31, S6, .3
 41, S6, .3
*DECURRENT, OP=NEW
 1, CS6,.3
 11, CS6,.3
 21, CS6,.3
 31, CS6,.3
 41, CS6,.3
*END STEP