*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
 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
*NSET,NSET=RING2,GEN
201,281
*NFIL,NSET=TIPMESH,TWO STEP
TIP,RING2,2,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=FRONT
    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
**
** 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_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