*HEADING
Mass Diffusion Analysis with quadratic elements using
Quadratic Stress analysis
MODEL BUILT IN mm,N,s
**
*PREPRINT,ECHO=NO,MODEL=NO,HIST=NO
*NODE,NSET=ALL,INPUT=difftocrack_node.inp
**
**
*ELEMENT, TYPE=DC2D8, ELSET=ALL,INPUT=difftocrack_quad_elements.inp
**
**
*NODE,NSET=TIP
*NGEN,NSET=TIP
**
** tipbound
**
*NSET, NSET=TIPBOUND
*NFILL,NSET=TIPMESH,SINGULAR=1,TWO STEP
TIP,TIPBOUND,50,100
*NSET,NSET=RING2,GEN
*NFIL,NSET=TIPMESH,TWO STEP
TIP,RING2,2,100
*ELEMENT,TYPE=DC2D8,ELSET=ALL
*ELGEN,ELSET=ALL
**
** front
**
*NSET,NSET=FRONT,GEN
*NSET, NSET=FRONT
**
** RIGHT
**
*NSET, NSET=RIGHT
**
** left
**
*NSET, NSET=LEFT
**
** crack
**
*NSET,NSET=CRACK,GEN
*NSET, NSET=CRACK
**
** top
**
*NSET, NSET=TOP
**
** top
**
*ELSET, ELSET=TOP
**
** all
**
*SOLID SECTION, ELSET=ALL, MATERIAL=STEEL
          1.,
**
** steel
**
*MATERIAL, NAME=STEEL
**
**
*ELASTIC, TYPE=ISO
      2.E+5,         0.3
*DIFFUSIVITY
*KAPPA,TYPE=PRESS
*SOLUBILITY
*PHYSICAL CONS,ABS=-273
*INITIAL CONDITIONS,TYPE=PRESS
ALL,0.
*INITIAL CONDITIONS,TYPE=CONC
ALL,50.
*INITIAL CONDITIONS,TYPE=TEMP
ALL,52.0
*STEP,UNSYMM=YES,AMPLITUDE=STEP
*MASS DIFFUSION,STEADY STATE
*BOUNDARY
RIGHT,11,11,264.9388
*PRESSURE STRESS,FILE=difftocrack_stress
*EL PRINT,FREQ=0
*NODE PRINT,NSET=FRONT
NNC11,
*NODE FILE,NSET=FRONT,FREQ=2
NNC,
*RESTART,WRITE,OVERLAY,FREQ=1
*EL FILE,FREQ=2
CONC,
SOL,
ESOL,
ISOL,
*END STEP