*HEADING
EC38MFDC: DFLUX FOR MASS DIFFUSION ELEMENTS  [DC3D8]
*restart,write
*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.
*NSET,NSET=FIX1
 1,2,3,4
*NSET,NSET=FIX2
 5,6,7,8
*NSET,NSET=FIX3
 1,2,6,5
*NSET,NSET=FIX4
 2,3,7,6
*NSET,NSET=FIX5
 3,4,7,8
*NSET,NSET=FIX6
 1,4,5,8
*ELEMENT,TYPE=DC3D8, ELSET=EALL
 1, 1,2,3,4,5,6,7,8
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*MATERIAL,NAME=A1
*DIFFUSIVITY
 3.77E-5, 
*SOLUBILITY
 1., 
*BOUNDARY
 FIX1,11
*surface,NAME=BOT
1,S1
*STEP
*MASS DIFFUSION, STEADY STATE
*NODE PRINT
 NNC,RFL
*EL PRINT,FREQ=0
*DFLUX
 1, BF, .3
*SECTION PRINT,NAME=BOT,SURFACE=BOT
*SECTION FILE, NAME=BOT,SURFACE=BOT
*OUTPUT,FIELD
*NODE OUTPUT,NSET=ALL
 NNC,RFL
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=ALL
 NNC,RFL
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX2, 11
*DFLUX, OP=NEW
 1, S1, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX1, 11
*DFLUX,OP=NEW
 1, 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
 1, S4, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX3, 11
*DFLUX,OP=NEW
 1, S5, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX4, 11
*DFLUX,OP=NEW
 1, S6, .3
*END STEP