*HEADING
 FINITE LENGTH AXIAL FLAW IN PRESSURIZED CYLINDER - DEEP CRACK   
*NODE,SYSTEM=C
 1, 10.5,-90.,0.
 97, 10.5,-90.,9.6
 99, 10.5,-90.,14.4
 103, 10.5,-90.,33.6
 9601, 10.5,234.,0.
 9697, 10.5,234.,9.6
 9699, 10.5,234.,14.4
 9703, 10.5,234.,33.6
 11201, 10.5,90.,0.
 11297, 10.5,90.,9.6
 11299, 10.5,90.,14.4
 11303, 10.5,90.,33.6
 99991, 0.,0.,0.
 99992, 0.,0.,9.6
 99993, 0.,0.,14.4
 99994, 0.,0.,33.6
*NGEN,NSET=BOTTOM,LINE=C
 1,9601,200,99991,,,,,,-1.
 9601,11201,200,99991,,,,,,-1.
*NGEN,NSET=MID1,LINE=C
 97,9697,200,99992,,,,,,-1.
 9697,11297,200,99992,,,,,,-1.
*NGEN,NSET=MID2,LINE=C
 99,9699,200,99993,,,,,,-1.
 9699,11299,200,99993,,,,,,-1.
*NGEN,NSET=TOP,LINE=C
 103,9703,200,99994,,,,,,-1.
 9703,11303,200,99994,,,,,,-1.
*NGEN
 1,97
 97,99
 99,103
*NSET,NSET=EDGE1,GENERATE
 25,103
*NGEN,NSET=EDGE2
 11201,11297
 11297,11299
 11299,11303
*NFILL
 BOTTOM,MID1,96
 MID1,MID2,2
 MID2,TOP,4
*NSET,NSET=N1
 103,11303
*NSET,NSET=N2
 4903,9903,10303,10703,11103
*NSET,NSET=N3
 9703,10103,10503,10903
*NSET,NSET=ND1
 1, 
*ELEMENT,TYPE=S8R,ELSET=SHELL
 9, 1,5,405,401,3,205,403,201
 13, 17,19,419,417,18,219,418,217
 17, 25,29,429,425,27,229,427,225
 18, 29,37,437,429,33,237,433,229
 22, 61,77,477,461,69,277,469,261
 23, 77,97,497,477,87,297,487,277
 24, 401,409,1209,1201,405,809,1205,801
 26, 1201,1209,2809,2801,1205,2009,2805,2001
 28, 417,421,1221,1217,419,821,1219,817
 31, 1217,1221,2821,2817,1219,2021,2819,2017
 34, 429,437,1237,1229,433,837,1233,829
 38, 1229,1237,2837,2829,1233,2037,2833,2029
 42, 461,477,1277,1261,469,877,1269,861
 43, 1261,1277,2877,2861,1269,2077,2869,2061
 44, 477,497,1297,1277,487,897,1287,877
 45, 1277,1297,2897,2877,1287,2097,2887,2077
 46, 2801,2817,6017,6001,2809,4417,6009,4401
 48, 2817,2821,6021,6017,2819,4421,6019,4417
 50, 2821,2829,6029,6021,2825,4429,6025,4421
 58, 2877,2897,6097,6077,2887,4497,6087,4477
 47, 6001,6017,9617,9601,6009,7817,9609,7801
 49, 6017,6021,9621,9617,6019,7821,9619,7817
 51, 6021,6029,9629,9621,6025,7829,9625,7821
 59, 6077,6097,9697,9677,6087,7897,9687,7877
 1002, 97,99,9699,9697,98,4899,9698,4897
 1101, 9601,9697,10097,10001,9649,9897,10049,9801
 1102, 9697,9699,10099,10097,9698,9899,10098,9897
*ELGEN,ELSET=SHELL
 9,4,4,1
 13,4,2,1
 18,4,8,1
 24,2,8,1
 26,2,8,1
 28,3,4,1
 31,3,4,1
 34,4,8,1
 38,4,8,1
 46,2,28,6
 52,3,16,2
 47,2,28,6
 53,3,16,2
 1002,3,2,1
 1102,3,2,1
 1101,4,400,100
 1102,4,400,100
 1103,4,400,100
 1104,4,400,100
*ELSET,ELSET=SHELL1,GENERATE
 1,59
*MATERIAL,NAME=MSHELL
*ELASTIC
 30.E6,.3
*PLASTIC
7.E4,0.
1.E5,.0976
*SHELL SECTION,ELSET=SHELL,MATERIAL=MSHELL
 1., 
*ELEMENT,TYPE=LS3S
 1, 25,24,23
 5, 17,15,13
*ELGEN,ELSET=CRACK
 1, 4,-2,1
 5, 4,-4,1
*SURFACE FLAW,SIDE=NEGATIVE
    1,  0.8000000
    3,  0.7972174
    5,  0.7888106
    7,  0.7745967
    9,  0.7542472
   11,  0.7272475
   13,  0.6928203
   15,  0.6497863
   17,  0.5962848
   18,  0.5647025
   19,  0.5291503
   20,  0.4887626
   21,  0.4422166
   22,  0.3872983
   23,  0.3197221
   24,  0.2285218
   25,  0.0000000
*SHELL SECTION,ELSET=CRACK,MATERIAL=MSHELL
 1., 
*MPC
QUADRATIC, 403,401,405,409
QUADRATIC, 407,401,405,409
QUADRATIC, 411,409,413,417
QUADRATIC, 415,409,413,417
QUADRATIC, 418,417,419,421
QUADRATIC, 420,417,419,421
QUADRATIC, 422,421,423,425
QUADRATIC, 424,421,423,425
QUADRATIC, 2805,2801,2809,2817
QUADRATIC, 2813,2801,2809,2817
QUADRATIC, 2823,2821,2825,2829
QUADRATIC, 2827,2821,2825,2829
QUADRATIC, 2833,2829,2837,2845
QUADRATIC, 2841,2829,2837,2845
QUADRATIC, 2849,2845,2853,2861
QUADRATIC, 2857,2845,2853,2861
QUADRATIC, 9609,9601,9649,9697
QUADRATIC, 9617,9601,9649,9697
QUADRATIC, 9619,9601,9649,9697
QUADRATIC, 9621,9601,9649,9697
QUADRATIC, 9625,9601,9649,9697
QUADRATIC, 9629,9601,9649,9697
QUADRATIC, 9637,9601,9649,9697
QUADRATIC, 9645,9601,9649,9697
QUADRATIC, 9653,9601,9649,9697
QUADRATIC, 9661,9601,9649,9697
QUADRATIC, 9669,9601,9649,9697
QUADRATIC, 9677,9601,9649,9697
QUADRATIC, 9687,9601,9649,9697
QUADRATIC, 7897,9697,4897,97
QUADRATIC, 6097,9697,4897,97
QUADRATIC, 4497,9697,4897,97
QUADRATIC, 2897,9697,4897,97
QUADRATIC, 2097,9697,4897,97
QUADRATIC, 1297,9697,4897,97
QUADRATIC, 897,9697,4897,97
QUADRATIC, 497,9697,4897,97
QUADRATIC, 297,9697,4897,97
*RESTART,WRITE,FREQUENCY=5
*STEP
*STATIC
1.,1.
*CLOAD
N1,3,1.15E4
N3,3,2.30E4
N2,3,4.60E4
*DLOAD
SHELL,P,3325.
*BOUNDARY
 EDGE1,XSYMM
 EDGE2,XSYMM
 BOTTOM,ZSYMM
 11201,2
*NODE PRINT,NSET=ND1,SUMMARY=NO,FREQUENCY=10
U, 
RF, 
*EL PRINT,ELSET=CRACK
JK, 
*ENERGY PRINT
*NODE FILE,NSET=ND1
U,RF
*OUTPUT,FIELD
*NODE OUTPUT,NSET=ND1
U,RF
*EL FILE,ELSET=CRACK
JK,
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=CRACK
JK,
*EL FILE,ELSET=SHELL
LOADS,
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=SHELL
LOADS,
*ENERGY FILE
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*END STEP
*STEP
*STATIC
.1,1.,,.1
*CLOAD
N1,3,.2875E5
N3,3,.575E5
N2,3,1.15E5
*DLOAD
SHELL,P,8312.5
*END STEP