*HEADING
EC3FMFDC: DFLUX FOR MASS DIFFUSION ELEMENTS  [DC3D15]
*PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES
*NODE
  1, 0.0, 0.0, 0.0
  2, 7.0, 0.0, 0.0
  3, 7.0, 7.0, 0.0
  4, 0.0, 7.0, 0.0
  5, 0.0, 0.0, 7.0
  6, 7.0, 0.0, 7.0
  7, 7.0, 7.0, 7.0
  8, 0.0, 7.0, 7.0
  9, 3.5, 0.0, 0.0
 10, 7.0, 3.5, 0.0
 11, 3.5, 7.0, 0.0
 12, 0.0, 3.5, 0.0
 13, 3.5, 0.0, 7.0
 14, 7.0, 3.5, 7.0
 15, 3.5, 7.0, 7.0
 16, 0.0, 3.5, 7.0
 17, 0.0, 0.0, 3.5
 18, 7.0, 0.0, 3.5
 19, 7.0, 7.0, 3.5
 20, 0.0, 7.0, 3.5
 21, 3.5, 3.5, 0.0
 22, 3.5, 3.5, 7.0
 23, 3.5, 0.0, 3.5
 24, 7.0, 3.5, 3.5
 25, 3.5, 7.0, 3.5
 26, 0.0, 3.5, 3.5
 27, 3.5, 3.5, 3.5
*ELEMENT,TYPE=DC3D15,ELSET=ALLELS
 1, 1,2,4,5,6,8,9,21,12,13,22,16,17,18,20
 2, 2,3,4,6,7,8,10,11,21,14,15,22,18,19,20
*SOLID SECTION,MATERIAL=A1, ELSET=ALLELS
*MATERIAL,NAME=A1
*DIFFUSIVITY
 3.77E-5, 
*SOLUBILITY
 1., 
*NSET,NSET=FIX1
1,2,3,4,9,10,11,12,21
*NSET,NSET=FIX2
5,6,7,8,13,14,15,16,22
*NSET,NSET=FIX3
1,2,6,5,9,18,13,17
*NSET,NSET=FIX4
2,3,7,6,10,19,14,18
*NSET,NSET=FIX5
3,4,8,7,11,20,15,19
*NSET,NSET=FIX6
4,1,5,8,12,17,16,20
*BOUNDARY
 FIX1,11
*STEP
*MASS DIFFUSION, STEADY STATE
*NODE PRINT
 NNC,RFL
*OUTPUT,FIELD
*NODE OUTPUT
 NNC,RFL
*EL PRINT,FREQ=0
*OUTPUT,FIELD
*ELEMENT OUTPUT
*DFLUX
 ALLELS, BF, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX2, 11
*DFLUX, OP=NEW
 1, S1, .3
 2, S1, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX1, 11
*DFLUX,OP=NEW
 1, S2, .3
 2, S2, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX5, 11
*DFLUX,OP=NEW
 1, S3, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX6, 11
*DFLUX,OP=NEW
 2, S3, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX3, 11
*DFLUX,OP=NEW
 2, S4, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX4, 11
*DFLUX,OP=NEW
 1, S5, .3
*END STEP