*HEADING
 FINITE LENGTH AXIAL FLAW IN PRESSURIZED CYLINDER     
*NODE,SYSTEM=C
 1, 10.5,-90.,0.
 97, 10.5,-90.,6.
 99, 10.5,-90.,9.
 103, 10.5,-90.,21.
 9601, 10.5,234.,0.
 9697, 10.5,234.,6.
 9699, 10.5,234.,9.
 9703, 10.5,234.,21.
 11201, 10.5,90.,0.
 11297, 10.5,90.,6.
 11299, 10.5,90.,9.
 11303, 10.5,90.,21.
 99991, 0.,0.,0.
 99992, 0.,0.,6.
 99993, 0.,0.,9.
 99994, 0.,0.,21.
*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.5000000
    3,  0.4982609
    5,  0.4930066
    7,  0.4841229
    9,  0.4714045
   11,  0.4545297
   13,  0.4330127
   15,  0.4061164
   17,  0.3726780
   18,  0.3529390
   19,  0.3307189
   20,  0.3054766
   21,  0.2763854
   22,  0.2420615
   23,  0.1998263
   24,  0.1428261
   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
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