*HEADING
 FINITE LENGTH AXIAL FLAW IN PRESSURIZED CYLINDER - SHALLOW CRACK
 C3D20R Submodel
**
*NODE
**
**
*ELEMENT, TYPE=C3D20R, ELSET=EALL
**
** DRIVE
**
*NSET, NSET=DRIVE
******************************************
******************************************
** NODE SET FOR J-INTEGRAL CONTOURS
** --------------------------------
*NSET, NSET=J0, GENERATE
*NSET, NSET=J1, GENERATE
*NSET, NSET=J2, GENERATE
*NSET, NSET=J3, GENERATE
*NSET, NSET=J4, GENERATE
*NSET, NSET=J5, GENERATE
*NSET, NSET=J6, GENERATE
*NSET, NSET=J7, GENERATE
*NSET, NSET=J8, GENERATE
*NSET, NSET=J9, GENERATE
*NSET, NSET=J10, GENERATE
*NSET, NSET=J11, GENERATE
*NSET, NSET=J12, GENERATE
*NSET, NSET=J13, GENERATE
*NSET, NSET=J14, GENERATE
*NSET, NSET=J15, GENERATE
*NSET, NSET=J16, GENERATE
************************************
*SOLID SECTION, ELSET=EALL,MATERIAL=BASE
*MATERIAL, NAME=BASE
*ELASTIC
30.E6,0.3
*PLASTIC
7.E4,0.
1.E5,0.0976
************************************
** ELEMENT SETS FOR PLOTTING
** -------------------------
*ELSET, ELSET=SEG1, GENERATE
*ELSET, ELSET=SEG1A, GENERATE
*ELSET, ELSET=CORNER, GENERATE
*ELSET, ELSET=CRACK, GENERATE
*ELSET, ELSET=SEG2, GENERATE
*ELSET, ELSET=CUTOUT
SEG1A,SEG2
*NSET,NSET=Z_CON,GENERATE
*BOUNDARY
Z_CON,     3
40000,     1
 2001,     1
 2002,     1
 2003,     1
 1001,     1
 1003,     1
    1,     1
    2,     1
****************************
*NORMAL
*MPC
QUADR,   221,   220,   222,   224
QUADR,  2221,  2220,  2222,  2224
C BIQ,  1222,   220,   224,  2224,  2220,   222,  1224,  2222,  1220
QUADR,  4221,  4220,  4222,  4224
C BIQ,  3222,  2220,  2224,  4224,  4220,  2222,  3224,  4222,  3220
QUADR,  6221,  6220,  6222,  6224
C BIQ,  5222,  4220,  4224,  6224,  6220,  4222,  5224,  6222,  5220
QUADR,  8221,  8220,  8222,  8224
C BIQ,  7222,  6220,  6224,  8224,  8220,  6222,  7224,  8222,  7220
QUADR, 10221, 10220, 10222, 10224
C BIQ,  9222,  8220,  8224, 10224, 10220,  8222,  9224, 10222,  9220
QUADR, 12221, 12220, 12222, 12224
C BIQ, 11222, 10220, 10224, 12224, 12220, 10222, 11224, 12222, 11220
QUADR, 14221, 14220, 14222, 14224
C BIQ, 13222, 12220, 12224, 14224, 14220, 12222, 13224, 14222, 13220
QUADR, 16221, 16220, 16222, 16224
C BIQ, 15222, 14220, 14224, 16224, 16220, 14222, 15224, 16222, 15220
QUADR,   223,   220,   222,   224
QUADR,  2223,  2220,  2222,  2224
QUADR,  4223,  4220,  4222,  4224
QUADR,  6223,  6220,  6222,  6224
QUADR,  8223,  8220,  8222,  8224
QUADR, 10223, 10220, 10222, 10224
QUADR, 12223, 12220, 12222, 12224
QUADR, 14223, 14220, 14222, 14224
QUADR, 16223, 16220, 16222, 16224
QUADR,   225,   224,   226,   228
QUADR,  2225,  2224,  2226,  2228
C BIQ,  1226,   224,   228,  2228,  2224,   226,  1228,  2226,  1224
QUADR,  4225,  4224,  4226,  4228
C BIQ,  3226,  2224,  2228,  4228,  4224,  2226,  3228,  4226,  3224
QUADR,  6225,  6224,  6226,  6228
C BIQ,  5226,  4224,  4228,  6228,  6224,  4226,  5228,  6226,  5224
QUADR,  8225,  8224,  8226,  8228
C BIQ,  7226,  6224,  6228,  8228,  8224,  6226,  7228,  8226,  7224
QUADR, 10225, 10224, 10226, 10228
C BIQ,  9226,  8224,  8228, 10228, 10224,  8226,  9228, 10226,  9224
QUADR, 12225, 12224, 12226, 12228
C BIQ, 11226, 10224, 10228, 12228, 12224, 10226, 11228, 12226, 11224
QUADR, 14225, 14224, 14226, 14228
C BIQ, 13226, 12224, 12228, 14228, 14224, 12226, 13228, 14226, 13224
QUADR, 16225, 16224, 16226, 16228
C BIQ, 15226, 14224, 14228, 16228, 16224, 14226, 15228, 16226, 15224
QUADR,   227,   224,   226,   228
QUADR,  2227,  2224,  2226,  2228
QUADR,  4227,  4224,  4226,  4228
QUADR,  6227,  6224,  6226,  6228
QUADR,  8227,  8224,  8226,  8228
QUADR, 10227, 10224, 10226, 10228
QUADR, 12227, 12224, 12226, 12228
QUADR, 14227, 14224, 14226, 14228
QUADR, 16227, 16224, 16226, 16228
QUADR,   229,   228,   230,   232
QUADR,  2229,  2228,  2230,  2232
C BIQ,  1230,   228,   232,  2232,  2228,   230,  1232,  2230,  1228
QUADR,  4229,  4228,  4230,  4232
C BIQ,  3230,  2228,  2232,  4232,  4228,  2230,  3232,  4230,  3228
QUADR,  6229,  6228,  6230,  6232
C BIQ,  5230,  4228,  4232,  6232,  6228,  4230,  5232,  6230,  5228
QUADR,  8229,  8228,  8230,  8232
C BIQ,  7230,  6228,  6232,  8232,  8228,  6230,  7232,  8230,  7228
QUADR, 10229, 10228, 10230, 10232
C BIQ,  9230,  8228,  8232, 10232, 10228,  8230,  9232, 10230,  9228
QUADR, 12229, 12228, 12230, 12232
C BIQ, 11230, 10228, 10232, 12232, 12228, 10230, 11232, 12230, 11228
QUADR, 14229, 14228, 14230, 14232
C BIQ, 13230, 12228, 12232, 14232, 14228, 12230, 13232, 14230, 13228
QUADR, 16229, 16228, 16230, 16232
C BIQ, 15230, 14228, 14232, 16232, 16228, 14230, 15232, 16230, 15228
QUADR,   231,   228,   230,   232
QUADR,  2231,  2228,  2230,  2232
QUADR,  4231,  4228,  4230,  4232
QUADR,  6231,  6228,  6230,  6232
QUADR,  8231,  8228,  8230,  8232
QUADR, 10231, 10228, 10230, 10232
QUADR, 12231, 12228, 12230, 12232
QUADR, 14231, 14228, 14230, 14232
QUADR, 16231, 16228, 16230, 16232
************************************
*AMPLITUDE, NAME=PRESSURE, TIME=TOTAL TIME
0.,0.,1.,5985.,2.,9975.
************************************
*RESTART, WRITE,FREQ=1000
************************************
*SUBMODEL,SHELL TO SOLID,SHELL THICKNESS=1.0
DRIVE,
************************************
*STEP
*STATIC
1.0,1.0
*BOUNDARY,SUBMODEL,STEP=1
DRIVE,
*DLOAD, AMPLITUDE=PRESSURE
 8001,    P1,        1.0
 8002,    P1,        1.0
 8003,    P1,        1.0
 8004,    P1,        1.0
 8005,    P1,        1.0
 8006,    P1,        1.0
 8007,    P1,        1.0
 8008,    P1,        1.0
 8009,    P1,        1.0
 8010,    P1,        1.0
 8011,    P1,        1.0
 8012,    P1,        1.0
 8013,    P1,        1.0
 8014,    P1,        1.0
 8015,    P1,        1.0
 8016,    P1,        1.0
 8017,    P1,        1.0
 8018,    P1,        1.0
 8019,    P1,        1.0
 8020,    P1,        1.0
 8021,    P1,        1.0
 8022,    P1,        1.0
 8023,    P1,        1.0
 8024,    P1,        1.0
 8025,    P1,        1.0
 8026,    P1,        1.0
 8027,    P1,        1.0
 8028,    P1,        1.0
 8029,    P1,        1.0
 8030,    P1,        1.0
 8031,    P1,        1.0
 8032,    P1,        1.0
 1033,    P6,        1.0
 2033,    P6,        1.0
 3033,    P6,        1.0
 4033,    P6,        1.0
 5033,    P6,        1.0
 6033,    P6,        1.0
 7033,    P6,        1.0
 1034,    P6,        1.0
 2034,    P6,        1.0
 3034,    P6,        1.0
 4034,    P6,        1.0
 5034,    P6,        1.0
 6034,    P6,        1.0
 7034,    P6,        1.0
 8035,    P1,        1.0
 8036,    P1,        1.0
 8037,    P1,        1.0
 8038,    P1,        1.0
 8039,    P1,        1.0
 8040,    P1,        1.0
 8041,    P1,        1.0
14041,    P6,        1.0
 8042,    P1,        1.0
14042,    P6,        1.0
 1043,    P6,        1.0
 2043,    P6,        1.0
 3043,    P6,        1.0
 4043,    P6,        1.0
 5043,    P6,        1.0
 6043,    P6,        1.0
 7043,    P6,        1.0
 8044,    P1,        1.0
 8045,    P1,        1.0
 8046,    P1,        1.0
 8047,    P1,        1.0
 8048,    P1,        1.0
 8049,    P1,        1.0
 8050,    P1,        1.0
 8051,    P1,        1.0
14051,    P6,        1.0
 1052,    P6,        1.0
 2052,    P6,        1.0
 3052,    P6,        1.0
 4052,    P6,        1.0
 5052,    P6,        1.0
 6052,    P6,        1.0
 7052,    P6,        1.0
 8053,    P1,        1.0
 8054,    P1,        1.0
 8055,    P1,        1.0
 8056,    P1,        1.0
 8057,    P1,        1.0
14057,    P6,        1.0
 1058,    P6,        1.0
 2058,    P6,        1.0
 3058,    P6,        1.0
 4058,    P6,        1.0
 5058,    P6,        1.0
 6058,    P6,        1.0
 7058,    P6,        1.0
 8059,    P1,        1.0
 8060,    P1,        1.0
 8061,    P1,        1.0
 8062,    P1,        1.0
 8063,    P1,        1.0
14063,    P6,        1.0
*EL PRINT, POSITION=CENTROID,FREQ=40
S,MISES
*NODE PRINT,FREQ=40
U,RF
*NODE FILE,FREQ=1,NSET=J0
U,
*CONTOUR INTEGRAL,NORMAL,OUTPUT=BOTH,CONTOURS=4, SYMM
J0,J1,J2,J3,J4,J5,J6,J7,J8,J9,J10,J11,J12,J13,J14,J15
J16,
*END STEP
********************
*STEP
*STATIC
0.05,1.0,,0.05
*BOUNDARY,SUBMODEL,STEP=2
DRIVE,
*END STEP