*HEADING
DIFFUSION ANALYSIS                            
*PREPRINT,ECHO=NO,MODEL=NO,HISTORY=NO
*NODE, NSET=A
1,0.0,-0.005
*NODE, NSET=B
161,0.1,-0.005
*NODE, NSET=C
2737,0.1,0.005
*NODE, NSET=D
2577,0.0,0.005
*NFILL,NSET=LEFT
A,D,16,161
*NFILL,NSET=RIGHT
B,C,16,161
*NFILL,NSET=PLATE
LEFT,RIGHT,160,1
*NSET, NSET=TOP,GENERATE
2577,2737
*NSET,NSET=BOTTOM,GENERATE
1,161
*NSET,NSET=CONL,GENERATE
1,2577,322
*NSET,NSET=CONR,GENERATE
161,2737,322
** GENERATE A NODE SET FOR OBSERVATION AT X=7.5cm
*NSET,NSET=RPOINT,GENERATE
121,2697,161
*NSET,NSET=LPOINT,GENERATE
17,2593,161
*ELEMENT,TYPE=DC2D4,ELSET=SC2D
1,1,3,325,323
*ELGEN, ELSET=SC2D
1,80,2,1,8,322,80,1
*ELSET,ELSET=LPOINT,GENERATE
9,569,80
*ELSET,ELSET=RPOINT,GENERATE
61,621,80
*SOLID SECTION,ELSET=SC2D,MATERIAL=STEEL
1.0D0, 
*PHYSICAL CONST,ABSOLUTEZ=0.
*MATERIAL,NAME=STEEL
*ELASTIC
2e11,0.3
*CONDUCTIVITY
1e3, 
*SOLUBILITY
 1., 
*DIFFUSIVITY
 3.6E-3,
*KAPPA,TYPE=TEMP
0.000000000E+00, 0.000000E+00, 0.000000E+00
0.000000000E+00, 0.150000E+03, 0.000000E+00
0.000000000E+00, 0.000000E+00, 0.123000E+03
0.266500000E+01, 0.150000E+03, 0.123000E+03
0.000000000E+00, 0.000000E+00, 0.273000E+03
0.591500000E+01, 0.150000E+03, 0.273000E+03
0.000000000E+00, 0.000000E+00, 0.373000E+03
0.808166667E+01, 0.150000E+03, 0.373000E+03
0.000000000E+00, 0.000000E+00, 0.423000E+03
0.916500000E+01, 0.150000E+03, 0.423000E+03
0.000000000E+00, 0.000000E+00, 0.473000E+03
0.102483333E+02, 0.150000E+03, 0.473000E+03
*KAPPA,TYPE=PRESS
0.0,0.0
1e-6,150
**
*INITIAL CONDITIONS, TYPE=TEMPERATURE,FILE=thermomechdiffusion_tempdisp.odb,STEP=1,INC=0
*INITIAL CONDITIONS,TYPE=CONCENTRATION
PLATE,50.0D-0
*INITIAL CONDITIONS,TYPE=PRESSURE,FILE=thermomechdiffusion_tempdisp,STEP=1,INC=0
**
*STEP,INC=100,AMPLITUDE=RAMP
STEP 1, CONCENTRATION REACHES STEADY STATE
*MASS DIFFUSION,STEADY STATE
0.001,.001
*BOUNDARY
LEFT,11,11,50.
RIGHT,11,11,100.
*RESTART,WRITE,FREQ=20
*NODE PRINT,FREQ=0
*EL PRINT,FREQ=0
*EL FILE
ESOL, 
SOL, 
CONC, 
NFLUX, 
SOL, 
ISOL, 
MFL, 
*NODE FILE
NNC, 
CFL, 
RFL, 
*EL FILE, POSITION=AVERAGED AT NODES
TEMP, 
*END STEP
*STEP,INC=100,AMPLITUDE=RAMP
STEP 2, PRESSURE APPLIED
*MASS DIFFUSION,DCMAX=0.5,END=SS
100.,1.0E6, , ,1e-6
*PRESSURE STRESS,FILE=thermomechdiffusion_tempdisp,BINC=1,BSTEP=1,ESTEP=2
*RESTART,WRITE,FREQ=100
*NODE PRINT,FREQ=0
*NODE PRINT,FREQ=0
*EL FILE,FREQ=5
ESOL, 
SOL, 
CONC, 
NFLUX, 
SOL, 
ISOL, 
MFL, 
*NODE FILE,FREQ=5
NNC, 
CFL, 
RFL, 
*EL FILE, POSITION=AVERAGED AT NODES,FREQ=5
TEMP, 
*END STEP
**
*STEP,INC=1000,AMPLITUDE=RAMP
STEP 3, TEMPERATURE APPLIED
*MASS DIFFUSION,DCMAX=0.5,END=SS
100.,1.0E6, , ,1e-6
*TEMPERATURE,FILE=thermomechdiffusion_tempdisp.odb,BINC=1,BSTEP=1,ESTEP=2
*RESTART,WRITE,FREQ=20
*EL FILE,FREQ=5
ESOL, 
SOL, 
CONC, 
NFLUX, 
SOL, 
ISOL, 
MFL, 
*NODE PRINT,FREQ=0
*EL PRINT,FREQ=0
*NODE FILE,FREQ=5
NNC, 
CFL, 
RFL, 
*EL FILE, POSITION=AVERAGED AT NODES,FREQ=5
TEMP, 
*END STEP