*HEADING
 FINITE LENGTH AXIAL FLAW IN PRESSURIZED CYLINDER - SHALLOW CRACK  
*NODE,SYSTEM=C
 1, 10.5,-90.,0.
 97, 10.5,-90.,4.5
 99, 10.5,-90.,6.75
 103, 10.5,-90.,15.75
 9601, 10.5,234.,0.
 9697, 10.5,234.,4.5
 9699, 10.5,234.,6.75
 9703, 10.5,234.,15.75
 11201, 10.5,90.,0.
 11297, 10.5,90.,4.5
 11299, 10.5,90.,6.75
 11303, 10.5,90.,15.75
 99991, 0.,0.,0.
 99992, 0.,0.,4.5
 99993, 0.,0.,6.75
 99994, 0.,0.,15.75
*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.2500000
    3,  0.2491304
    5,  0.2465033
    7,  0.2420615
    9,  0.2357023
   11,  0.2272648
   13,  0.2165064
   15,  0.2030582
   17,  0.1863390
   18,  0.1764695
   19,  0.1653595
   20,  0.1527383
   21,  0.1381927
   22,  0.1210307
   23,  0.0999132
   24,  0.0714131
   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,20700.
N3,3,41400.
N2,3,82800.
*DLOAD
SHELL,P,5985.
*BOUNDARY
 EDGE1,XSYMM
 EDGE2,XSYMM
 BOTTOM,ZSYMM
 11201,2
*NODE PRINT,NSET=ND1
U, 
RF, 
*EL PRINT,ELSET=CRACK
JK, 
*ENERGY PRINT
*NODE FILE
U,RF
*OUTPUT,FIELD
*NODE OUTPUT
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,INC=200
*STATIC
.05,1.,,.05
*CLOAD
N1,3,34500.
N3,3,69000.
N2,3,138000.
*DLOAD
SHELL,P,9975.
*BOUNDARY
 EDGE1,XSYMM
 EDGE2,XSYMM
 BOTTOM,ZSYMM
 11201,2
*END STEP