*HEADING
HYDROGEN DIFFUSION OF A SHELL SECTION  ***  DC3D20   
*RESTART,WRITE,FREQUENCY=1000
*NODE
1,
2001,,10.
41,5.,
2041,5.,10.
121,205.
2121,205.,10.
10001,,,.5
12001,,10.,.5
10041,5.,,.5
12041,5.,10.,.5
10121,205.,,.5
12121,205.,10.,.5
20001,,,1.
22001,,10.,1.
20041,5.,,1.
22041,5.,10.,1.
20121,205.,,1.
22121,205.,10.,1.
*NGEN,NSET=LHEND
1,2001,1000
10001,12001,1000
20001,22001,1000
*NGEN,NSET=MID
41,2041,1000
10041,12041,1000
20041,22041,1000
*NGEN,NSET=RHEND
121,2121,1000
10121,12121,1000
20121,22121,1000
*NFILL,NSET=NALL,TWO STEP,BIAS=1.125 
LHEND,MID,40,1
*NFILL,NSET=NALL,TWO STEP,BIAS=.875
MID,RHEND,80,1
*NSET,NSET=NOUT,GENERATE
1,121,1
2001,2121,1
20001,20121,1
22001,22121,1
*NSET,NSET=NRFL
20121,21121,22121
*ELEMENT,TYPE=DC3D20
1,1,3,2003,2001,20001,20003,22003,22001,2,1003,2002,1001,20002,21003,22002,
21001,10001,10003,12003,12001
*ELGEN,ELSET=WELD
1,20,2,2
*ELEMENT,TYPE=DC3D20
41,41,43,2043,2041,20041,20043,22043,22041,42,1043,2042,1041,20042,21043,22042,
21041,10041,10043,12043,12041
*ELGEN,ELSET=BASE
41,40,2,2 
*ELSET,ELSET=ELS1
39,41
*SOLID SECTION,MATERIAL=WELD,ELSET=WELD
*PHYSICAL CONSTANTS,ABSOLUTE ZERO=-273.15
*MATERIAL,NAME=WELD
*SOLUBILITY
2.48858D+01, 0.000D+00
3.46472D+01, 2.500D+01
4.58298D+01, 5.000D+01
5.82346D+01, 7.500D+01
7.16596D+01, 1.000D+02
8.59120D+01, 1.250D+02
1.00815D+02, 1.500D+02
1.16210D+02, 1.750D+02
1.31960D+02, 2.000D+02
1.47945D+02, 2.250D+02
1.64063D+02, 2.500D+02
1.80229D+02, 2.750D+02
1.96372D+02, 3.000D+02
2.12432D+02, 3.250D+02
2.28360D+02, 3.500D+02
2.44118D+02, 3.750D+02
2.59672D+02, 4.000D+02
2.74998D+02, 4.250D+02
2.90077D+02, 4.500D+02
2.92704D+02, 4.544D+02
*DIFFUSIVITY
 1.62094D-07,, 0.000D+00                         
 1.29398D-06,, 2.500D+01
 7.49035D-06,, 5.000D+01
 3.36944D-05,, 7.500D+01
 1.23910D-04,, 1.000D+02
 3.86932D-04,, 1.250D+02
 1.05615D-03,, 1.500D+02
 2.57731D-03,, 1.750D+02
 5.72350D-03,, 2.000D+02
 1.17322D-02,, 2.250D+02
 2.24546D-02,, 2.500D+02
 4.05055D-02,, 2.750D+02
 6.94021D-02,, 3.000D+02
 1.13680D-01,, 3.250D+02
 1.78977D-01,, 3.500D+02
 2.72085D-01,, 3.750D+02
 4.00960D-01,, 4.000D+02
 5.74695D-01,, 4.250D+02
 8.03461D-01,, 4.500D+02
 8.50242D-01,, 4.544D+02
*SOLID SECTION,MATERIAL=BASE,ELSET=BASE
*MATERIAL,NAME=BASE
*SOLUBILITY
 1.78188D-01, 5.000D+01
 3.67776D-01, 7.500D+01
 6.88842D-01, 1.000D+02
 1.19243D-00, 1.250D+02
 1.93457D-00, 1.500D+02
 2.97365D-00, 1.750D+02
 4.36782D-00, 2.000D+02
 6.17278D-00, 2.250D+02
 8.43995D-00, 2.500D+02
 1.12152D+01, 2.750D+02
 1.45379D+01, 3.000D+02
 1.84407D+01, 3.250D+02
 2.29491D+01, 3.500D+02
 2.80820D+01, 3.750D+02
 3.38515D+01, 4.000D+02
 4.02639D+01, 4.250D+02
 4.73200D+01, 4.500D+02
 4.86282D+01, 4.544D+02
*DIFFUSIVITY
 7.85674D-03,, 0.000D+00                          
 3.34289D-02,, 2.500D+01 
 1.13043D-01,, 5.000D+01 
 3.17681D-01,, 7.500D+01 
 7.65060D-01,, 1.000D+02 
 1.61326D-00,, 1.250D+02 
 3.02697D-00,, 1.500D+02 
 5.12273D-00,, 1.750D+02 
 7.92165D-00,, 2.000D+02 
 1.13403D+01,, 2.250D+02 
 1.52239D+01,, 2.500D+02 
 1.93975D+01,, 2.750D+02 
 2.37067D+01,, 3.000D+02 
 2.80357D+01,, 3.250D+02 
 3.23085D+01,, 3.500D+02 
 3.64808D+01,, 3.750D+02 
 4.05306D+01,, 4.000D+02 
 4.44494D+01,, 4.250D+02 
 4.82372D+01,, 4.500D+02 
 4.88906D+01,, 4.544D+02 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,454.4
*INITIAL CONDITIONS,TYPE=CONCENTRATION
NALL,.1
*AMPLITUDE,NAME=AMP1,TIME=TOTAL TIME,VALUE=ABSOLUTE
0.,454.4,1.E-5,454.4,21.5,25.0,1000.,25.0
*STEP
*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,MFL3,MFLM
*EL PRINT,FREQUENCY=1000,POSITION=AVERAGED AT NODES,ELSET=WELD
CONC, 
*EL PRINT,FREQUENCY=1000,POSITION=AVERAGED AT NODES,ELSET=BASE
CONC, 
*NODE PRINT,FREQUENCY=1000
NNC,NT,RFL
*EL FILE,ELSET=ELS1,FREQUENCY=1000
CONC,
*EL FILE,ELSET=ELS1,POSITION=AVERAGED AT NODES,FREQUENCY=1000
CONC, 
*NODE FILE,FREQUENCY=1000
NNC,
*NODE FILE,NSET=NRFL,FREQUENCY=1000
RFL
*OUTPUT, FIELD
*ELEMENT OUTPUT
 CONC,MFL
*NODE OUTPUT
 NNC
*END STEP              
*STEP,INC=100
*MASS DIFFUSION,DCMAX=.01
.1,2.7,,,.0001
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*END STEP
*STEP,INC=100
*MASS DIFFUSION,DCMAX=.01
.1,2.5,,,.0001
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*END STEP
*STEP,INC=100
*MASS DIFFUSION,DCMAX=.01
.1,5.0,,,.0001
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*END STEP
*STEP,INC=1000
*MASS DIFFUSION,DCMAX=.05
.1,11.3,,,.0001
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*END STEP