*HEADING
 FINITE LENGTH AXIAL FLAW IN PRESSURIZED CYLINDER
*NODE,SYSTEM=C
***
***
*NGEN,NSET=BOTTOM,LINE=C
**
*NGEN,NSET=MID1,LINE=C
**
*NGEN,NSET=MID2,LINE=C
**
*NGEN,NSET=TOP,LINE=C
**
*NGEN
**
*NGEN,NSET=EDGE2
*NFILL
 BOTTOM,MID1,96
 MID1,MID2,2
 MID2,TOP,4
*NSET,NSET=N1
*NSET,NSET=N2
*NSET,NSET=N3
*ELEMENT,TYPE=S8R,ELSET=SHELL
**
*ELEMENT,TYPE=S8R,ELSET=SHELL
*ELGEN,ELSET=SHELL
**
*ELGEN,ELSET=SHELL
**
*ELEMENT,TYPE=S8R,ELSET=SHELL
**
*ELSET,ELSET=SHELL1,GENERATE
**
*MATERIAL,NAME=MSHELL
*ELASTIC
 30.E6,.3
*PLASTIC
7.E4,0.
1.E5,.0976
*SHELL SECTION,ELSET=SHELL,MATERIAL=MSHELL
 1.,
*ELEMENT,TYPE=LS6
*ELGEN,ELSET=CRACK
*SURFACE FLAW,SIDE=NEGATIVE
*SHELL SECTION,ELSET=CRACK,MATERIAL=MSHELL
 1.,
*MPC
BEAM, 20025,25
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
**
QUADRATIC, 20403,20401,20405,20409
QUADRATIC, 20407,20401,20405,20409
QUADRATIC, 20411,20409,20413,20417
QUADRATIC, 20415,20409,20413,20417
QUADRATIC, 20418,20417,20419,20421
QUADRATIC, 20420,20417,20419,20421
QUADRATIC, 20422,20421,20423,20425
QUADRATIC, 20424,20421,20423,20425
QUADRATIC, 22805,22801,22809,22817
QUADRATIC, 22813,22801,22809,22817
QUADRATIC, 22823,22821,22825,22829
QUADRATIC, 22827,22821,22825,22829
QUADRATIC, 22833,22829,22837,22845
QUADRATIC, 22841,22829,22837,22845
QUADRATIC, 22849,22845,22853,22861
QUADRATIC, 22857,22845,22853,22861
QUADRATIC, 29609,29601,29649,29697
QUADRATIC, 29617,29601,29649,29697
QUADRATIC, 29619,29601,29649,29697
QUADRATIC, 29621,29601,29649,29697
QUADRATIC, 29625,29601,29649,29697
QUADRATIC, 29629,29601,29649,29697
QUADRATIC, 29637,29601,29649,29697
QUADRATIC, 29645,29601,29649,29697
QUADRATIC, 29653,29601,29649,29697
QUADRATIC, 29661,29601,29649,29697
QUADRATIC, 29669,29601,29649,29697
QUADRATIC, 29677,29601,29649,29697
QUADRATIC, 29687,29601,29649,29697
QUADRATIC, 27897,29697,24897,97
QUADRATIC, 26097,29697,24897,97
QUADRATIC, 24497,29697,24897,97
QUADRATIC, 22897,29697,24897,97
QUADRATIC, 22097,29697,24897,97
QUADRATIC, 21297,29697,24897,97
QUADRATIC, 20897,29697,24897,97
QUADRATIC, 20497,29697,24897,97
QUADRATIC, 20297,29697,24897,97
**
*NSET,NSET=ND1
*RESTART,WRITE,FREQUENCY=5
*STEP
*STATIC
1.,1.,
*CLOAD
N3,3,2.30E4
N2,3,4.60E4
*DLOAD
SHELL,P,3325.
*BOUNDARY
 BOTTOM,ZSYMM
 11201,1,2
 11201,6
*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
N3,3,.575E5
N2,3,1.15E5
*DLOAD
SHELL,P,8312.5
*END STEP