*HEADING
HYDROGEN DIFFUSION OF A SHELL SECTION (FICK'S LAW)
*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
*ELEMENT,TYPE=DC2D8
*ELGEN,ELSET=WELD
*ELEMENT,TYPE=DC2D8
*ELGEN,ELSET=BASE
*ELSET,ELSET=ELS1
*SOLID SECTION,MATERIAL=WELD,ELSET=WELD
1.,
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0.0
*MATERIAL,NAME=WELD
*SOLUBILITY
*DIFFUSIVITY,LAW=FICK
*SOLID SECTION,MATERIAL=BASE,ELSET=BASE
1.,
*MATERIAL,NAME=BASE
*SOLUBILITY
*DIFFUSIVITY,LAW=FICK
*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
*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,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,MFL
*EL FILE,ELSET=ELS1,POSITION=AVERAGED AT NODES,FREQUENCY=1000
CONC,
*NODE FILE,FREQUENCY=1000
NNC,RFL
*END STEP
*STEP
*MASS DIFFUSION,DCMAX=.01
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*END STEP
*STEP
*MASS DIFFUSION,DCMAX=.01
*BOUNDARY
LHEND,11,11,0.
*TEMPERATURE,AMPLITUDE=AMP1
NALL,
*END STEP
*STEP
*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