*HEADING
HYDROGEN DIFFUSION OF A SHELL SECTION  ***  DC3D20
*RESTART,WRITE,FREQUENCY=1000
*NODE
*NGEN,NSET=LHEND
*NGEN,NSET=MID
*NGEN,NSET=RHEND
*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
*NSET,NSET=NRFL
*ELEMENT,TYPE=DC3D20
*ELGEN,ELSET=WELD
*ELEMENT,TYPE=DC3D20
*ELGEN,ELSET=BASE
*ELSET,ELSET=ELS1
*SOLID SECTION,MATERIAL=WELD,ELSET=WELD
*PHYSICAL CONSTANTS,ABSOLUTE ZERO=-273.15
*MATERIAL,NAME=WELD
*SOLUBILITY
*DIFFUSIVITY
*SOLID SECTION,MATERIAL=BASE,ELSET=BASE
*MATERIAL,NAME=BASE
*SOLUBILITY
*DIFFUSIVITY
*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
*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
*END STEP
*STEP,INC=100
*MASS DIFFUSION,DCMAX=.01
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*END STEP
*STEP,INC=100
*MASS DIFFUSION,DCMAX=.01
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*END STEP
*STEP,INC=100
*MASS DIFFUSION,DCMAX=.01
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*END STEP
*STEP,INC=1000
*MASS DIFFUSION,DCMAX=.05
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*END STEP