*HEADING
Mass Diffusion Analysis with quadratic elements using
Quadratic Singular 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
 1,2.0,0.
81,2.0,0.
*NGEN,NSET=TIP
1,81,1
**
** tipbound
**
*NSET, NSET=TIPBOUND
    5001,    5002,    5003,    5004,    5005,    5006,    5007,    5008,
    5009,    5010,    5011,    5012,    5013,    5014,    5015,    5016,
    5017,    5018,    5019,    5020,    5021,    5022,    5023,    5024,
    5025,    5026,    5027,    5028,    5029,    5030,    5031,    5032,
    5033,    5034,    5035,    5036,    5037,    5038,    5039,    5040,
    5041,    5042,    5043,    5044,    5045,    5046,    5047,    5048,
    5049,    5050,    5051,    5052,    5053,    5054,    5055,    5056,
    5057,    5058,    5059,    5060,    5061,    5062,    5063,    5064,
    5065,    5066,    5067,    5068,    5069,    5070,    5071,    5072,
    5073,    5074,    5075,    5076,    5077,    5078,    5079,    5080,
    5081,
*NFILL,NSET=TIPMESH,SINGULAR=1,TWO STEP
TIP,TIPBOUND,50,100
*ELEMENT,TYPE=DC2D8,ELSET=ALL
2001,1,201,203,3,101,202,103,2
*ELGEN,ELSET=ALL
2001,25,200,50,40,2,1
**
** front
**
*NSET,NSET=FRONT,GEN
1,5001,100
*NSET, NSET=FRONT1
    5001,    5522,    5533,    5554,    5565,    5586,    5597,    5618,
    5629,    5650,    5661,    5682,    5693,    5714,    5725,    5746,
    5757,    5778,    5789,    5810,    5821,    6802,    6813,    6834,
    6845,    6866,    6877,    6898,    6909,    6930,    6941,    6962,
    6973,    6994,    7005,    7026,    7037,    7058,    7069,    7090,
    7101,    7122,    7133,    7154,    7165,    7186,    7197,    7218,
    7229,    7250,    7261,    7282,    7293,    7314,    7325,    7346,
    7357,    7378,    7389,    7410,     7421
*NSET,NSET=FRONT
FRONT1,
**
** RIGHT
**
*NSET, NSET=RIGHT
    7421,    7422,    7423,    7424,    7425,    7426,    7427,    7428,
    7429,    7430,    7431,    7432,    7433,    7434,    7435,    7436,
    7437,    7438,    7439,    7440,    7441,    8642,    8653,    8674,
    8685,    8706,    8717,    8738,    8749,    8770,    8781,    8802,
    8813,    8834,    8845,    8866,    8877,    8898,    8909,    8930,
    8941,    8962,    8973,    8994,    9005,    9026,    9037,    9058,
    9069,    9090,    9101,    9122,    9133,    9154,    9165,    9186,
    9197,    9218,    9229,    9250,    9261,    9282,    9293,    9314,
    9325,    9346,     9357
**
** left
**
*NSET, NSET=LEFT
    6481,    6782,    6783,    6784,    6785,    6786,    6787,    6788,
    6789,    6790,    6791,    6792,    6793,    6794,    6795,    6796,
    6797,    6798,    6799,    6800,    6801,    8051,    8071,    8081,
    8101,    8111,    8131,    8141,    8161,    8171,    8191,    8201,
    8221,    8231,    8251,    8261,    8281,    8291,    8311,    8321,
    8341,    8351,    8371,    8381,    8401,    8411,    8431,    8441,
    8461,    8471,    8491,    8501,    8521,    8531,    8551,    8561,
    8581,    8591,    8611,    8621,    8641,    9387,    9392,    9402,
    9407,    9417,    9422,    9432,    9437,    9447,    9452,    9462,
    9467,    9477,    9482,    9492,    9497,    9507,    9512,    9522,
    9527,    9537,    9542,    9552,    9557,    9567,    9572,    9582,
    9587,    9597,    9602,    9612,    9617,    9627,    9632,    9642,
    9647,    9657,    9662,    9672,    9677,    9687,    9688,    9689,
    9690,    9691,    9692,    9745,     9746
**
** crack
**
*NSET,NSET=CRACK,GEN
81,5081,100
*NSET, NSET=CRACK
    5081,    6511,    6531,    6541,    6561,    6571,    6591,    6601,
    6621,    6631,    6651,    6661,    6681,    6691,    6711,    6721,
    6741,    6751,    6771,    6781,     6801
**
** top
**
*NSET, NSET=TOP
    9357,    9358,    9359,    9360,    9361,    9362,    9363,    9364,
    9365,    9366,    9367,    9368,    9369,    9370,    9371,    9372,
    9373,    9374,    9375,    9376,    9377,    9678,    9679,    9680,
    9681,    9682,    9683,    9684,    9685,    9686,     9687
**
** top
**
*ELSET, ELSET=TOP
    1221,    1222,    1223,    1224,    1225,    1226,    1227,    1228,
    1229,    1230,    1326,    1327,    1328,    1329,     1330
**
** all
**
*SOLID SECTION, ELSET=ALL, MATERIAL=STEEL
          1.,
**
** steel
**
*MATERIAL, NAME=STEEL
**
**
*ELASTIC, TYPE=ISO
      2.E+5,         0.3
*DIFFUSIVITY
34.096361E-6,
*KAPPA,TYPE=PRESS
0.0,0.0
3.921890,1E3
*SOLUBILITY
0.18872283, 
*PHYSICAL CONS,ABS=-273
*INITIAL CONDITIONS,TYPE=PRESS
ALL,0.
*INITIAL CONDITIONS,TYPE=CONC
ALL,264.9388
*INITIAL CONDITIONS,TYPE=TEMP
ALL,52.0
*STEP,UNSYMM=YES,AMPLITUDE=STEP
*MASS DIFFUSION,STEADY STATE
0.1,1.0
*BOUNDARY
RIGHT,11,11,264.9388
*PRESSURE STRESS,FILE=difftocrack_quarterpstress
*EL PRINT,FREQ=0
*NODE PRINT,NSET=FRONT
NNC11, 
*NODE FILE,NSET=FRONT1,FREQ=2
NNC, 
*RESTART,WRITE,OVERLAY,FREQ=1
*ELSET,ELSET=OUTPUT,GENERATE
500,1000
*EL FILE,ELSET=OUTPUT,FREQ=2
CONC, 
SOL, 
ESOL, 
ISOL, 
*END STEP