*HEADING
HYDROGEN DIFFUSION OF A SHELL SECTION       
*RESTART,WRITE,FREQUENCY=1000
*NODE
1,
2001,,10.
41,5.,
2041,5.,10.
121,205.
2121,205.,10.
*NGEN,NSET=LHEND
1,2001,1000
*NGEN,NSET=MID
41,2041,1000
*NGEN,NSET=RHEND
121,2121,1000
*NFILL,NSET=NALL,TWO STEP,BIAS=1.125 
LHEND,MID,40,1
*NFILL,NSET=NALL,TWO STEP,BIAS=.875
MID,RHEND,80,1
*ELEMENT,TYPE=DC2D8
1,1,3,2003,2001,2,1003,2002,1001
*ELGEN,ELSET=WELD
1,20,2,2
*ELEMENT,TYPE=DC2D8
41,41,43,2043,2041,42,1043,2042,1041
*ELGEN,ELSET=BASE
41,40,2,2 
*ELSET,ELSET=ELS1
39,41
*SOLID SECTION,MATERIAL=WELD,ELSET=WELD
1.,
*MATERIAL,NAME=WELD
*SOLUBILITY
2.48858D+01, 2.73150D+02
3.46472D+01, 2.98150D+02
4.58298D+01, 3.23150D+02
5.82346D+01, 3.48150D+02
7.16596D+01, 3.73150D+02
8.59120D+01, 3.98150D+02
1.00815D+02, 4.23150D+02
1.16210D+02, 4.48150D+02
1.31960D+02, 4.73150D+02
1.47945D+02, 4.98150D+02
1.64063D+02, 5.23150D+02
1.80229D+02, 5.48150D+02
1.96372D+02, 5.73150D+02
2.12432D+02, 5.98150D+02
2.28360D+02, 6.23150D+02
2.44118D+02, 6.48150D+02
2.59672D+02, 6.73150D+02
2.74998D+02, 6.98150D+02
2.90077D+02, 7.23150D+02
2.92704D+02, 7.27550D+02
*DIFFUSIVITY
 1.62094D-07,, 2.73150D+02
 1.29398D-06,, 2.98150D+02
 7.49035D-06,, 3.23150D+02
 3.36944D-05,, 3.48150D+02
 1.23910D-04,, 3.73150D+02
 3.86932D-04,, 3.98150D+02
 1.05615D-03,, 4.23150D+02
 2.57731D-03,, 4.48150D+02
 5.72350D-03,, 4.73150D+02
 1.17322D-02,, 4.98150D+02
 2.24546D-02,, 5.23150D+02
 4.05055D-02,, 5.48150D+02
 6.94021D-02,, 5.73150D+02
 1.13680D-01,, 5.98150D+02
 1.78977D-01,, 6.23150D+02
 2.72085D-01,, 6.48150D+02
 4.00960D-01,, 6.73150D+02
 5.74695D-01,, 6.98150D+02
 8.03461D-01,, 7.23150D+02
 8.50242D-01,, 7.27550D+02
*SOLID SECTION,MATERIAL=BASE,ELSET=BASE
1.,
*MATERIAL,NAME=BASE
*SOLUBILITY
 1.78188D-01, 3.23150D+02
 3.67776D-01, 3.48150D+02
 6.88842D-01, 3.73150D+02
 1.19243D-00, 3.98150D+02
 1.93457D-00, 4.23150D+02
 2.97365D-00, 4.48150D+02
 4.36782D-00, 4.73150D+02
 6.17278D-00, 4.98150D+02
 8.43995D-00, 5.23150D+02
 1.12152D+01, 5.48150D+02
 1.45379D+01, 5.73150D+02
 1.84407D+01, 5.98150D+02
 2.29491D+01, 6.23150D+02
 2.80820D+01, 6.48150D+02
 3.38515D+01, 6.73150D+02
 4.02639D+01, 6.98150D+02
 4.73200D+01, 7.23150D+02
 4.86282D+01, 7.27550D+02
*DIFFUSIVITY
 7.85674D-03,,  2.73150D+02
 3.34289D-02,,  2.98150D+02
 1.13043D-01,,  3.23150D+02
 3.17681D-01,,  3.48150D+02
 7.65060D-01,,  3.73150D+02
 1.61326D-00,,  3.98150D+02
 3.02697D-00,,  4.23150D+02
 5.12273D-00,,  4.48150D+02
 7.92165D-00,,  4.73150D+02
 1.13403D+01,,  4.98150D+02
 1.52239D+01,,  5.23150D+02
 1.93975D+01,,  5.48150D+02
 2.37067D+01,,  5.73150D+02
 2.80357D+01,,  5.98150D+02
 3.23085D+01,,  6.23150D+02
 3.64808D+01,,  6.48150D+02
 4.05306D+01,,  6.73150D+02
 4.44494D+01,,  6.98150D+02
 4.82372D+01,,  7.23150D+02
 4.88906D+01,,  7.27550D+02
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,727.55
*INITIAL CONDITIONS,TYPE=CONCENTRATION
NALL,.1
*AMPLITUDE,NAME=AMP1,TIME=TOTAL TIME,
VALUE=ABSOLUTE
0.,727.55,1.E-5,727.55,21.5,298.15,1000.,298.15
**
**
**
**
*STEP
STEADY-STATE SOLUTION UNDER OPERATING CONDITIONS
*MASS DIFFUSION,STEADY STATE
1.E-5,1.E-5 
*BOUNDARY                          
LHEND,11,11,.12247
RHEND,11,11,0.
*EL PRINT,FREQUENCY=1000
CONC,MFL,TEMP
*EL PRINT,FREQUENCY=1000
MFL1,MFL2,MFLM
*EL PRINT,FREQUENCY=1000,
POSITION=AVERAGED AT NODES,ELSET=WELD
CONC, 
*EL PRINT,FREQUENCY=1000,
POSITION=AVERAGED AT NODES,ELSET=BASE
CONC, 
*EL PRINT,FREQUENCY=1000
NFLUX, 
*NODE PRINT,FREQUENCY=1000 
NNC,NT,RFL
*OUTPUT,FIELD,Frequency=10
*ELEMENT OUTPUT
CONC,MFL,NFLUX 
*NODE OUTPUT
NNC,RFL
*END STEP              
*STEP
TRANSIENT ANALYSIS DUE TO SHUTDOWN -- PHASE 1
*MASS DIFFUSION,DCMAX=.01, END=PERIOD
.1,2.7,,,.0001
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*OUTPUT,FIELD,Frequency=10
*ELEMENT OUTPUT
CONC,MFL,NFLUX 
*NODE OUTPUT
NNC,RFL
*END STEP
*STEP
TRANSIENT ANALYSIS DUE TO SHUTDOWN -- PHASE 2
*MASS DIFFUSION,DCMAX=.01
.1,2.5,,,.0001
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*OUTPUT,FIELD,Frequency=10
*ELEMENT OUTPUT
CONC,MFL,NFLUX 
*NODE OUTPUT
NNC,RFL
*END STEP
*STEP
TRANSIENT ANALYSIS DUE TO SHUTDOWN -- PHASE 3
*MASS DIFFUSION,DCMAX=.01
.1,5.0,,,.0001
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*OUTPUT,FIELD,Frequency=10
*ELEMENT OUTPUT
CONC,MFL,NFLUX 
*NODE OUTPUT
NNC,RFL
*END STEP
*STEP,INC=1000
TRANSIENT ANALYSIS DUE TO SHUTDOWN -- PHASE 4
*MASS DIFFUSION,DCMAX=.05
.1,11.3,,,.0001
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*OUTPUT,FIELD,Frequency=100
*ELEMENT OUTPUT
CONC,MFL
*NODE OUTPUT
NNC
*END STEP