*HEADING
EC26MFDC: DFLUX FOR MASS DIFFUSION ELEMENTS  [DC2D6]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE,NSET=ALL
 1, 0.
 2, 7.
 3, 7., 7.
 4, 0., 7.
 5, 3.5, 0.
 6, 7., 3.5
 7, 3.5, 7.
 8, 0., 3.5
 9, 3.5, 3.5
*NSET,NSET=FIX1
 1,2,5
*NSET,NSET=FIX2
 2,3,6
*NSET,NSET=FIX3
 3,4,7
*NSET,NSET=FIX4
 1,4,8
*ELEMENT,TYPE=DC2D6, ELSET=EALL
 1,1,2,3,5,6,9
 2,1,3,4,9,7,8
**local orientation matching global system 
*ORIENTATION,NAME=RECT
 1.0, 0.0, 0.0, 0.0, 1.0, 0.0
 1, 0.0
*SOLID SECTION,MATERIAL=A1, ELSET=EALL,ORIENT=RECT
*MATERIAL,NAME=A1
*DIFFUSIVITY,TYPE=ANISO
 3.77E-5,3.77E-6,7.54E-5,3.77E-6,3.77E-6,11.31E-5
*SOLUBILITY
 1., 
*BOUNDARY
 FIX1, 11
*STEP
*MASS DIFFUSION, STEADY STATE
*NODE PRINT
 NNC,RFL
*EL PRINT,FREQ=0
*DFLUX,OP=NEW
 1, BF, .3
*OUTPUT,FIELD
*NODE OUTPUT
 NNC,RFL
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=ALL
 NNC,RFL
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY,OP=NEW
 FIX3,11
*DFLUX,OP=NEW
 1, S1, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY,OP=NEW
 FIX4,11
*DFLUX,OP=NEW
 1, S2, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY,OP=NEW
 FIX1,11
*DFLUX,OP=NEW
 2, S2, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY,OP=NEW
 FIX2,11
*DFLUX,OP=NEW
 2, S3, .3
*END STEP