*HEADING
AXISYMMETRIC CONICAL CRACK MODEL--GLOBAL WITH 1 RING      
*SYSTEM
10.,0.,0.
*NODE,SYSTEM=C
1833,0.,0.,0.
9833,0.,0.,0.
57833,0.,0.,0.
57865,0.,0.,0.
65833,0.,0.,0.
1001,15.,-45.,0.
65001,15.,-45.,0.
1033,8.,-45.,0.
65033,8.,-45.,0.
33033,25.,-45.,0.
41033,25.,-18.,0.
57033,10.,-18.,0.
9033,10.,-72.,0.
25033,25.,-72.,0.
41065,25.,0.,0.
57065,10.,0.,0.
*SYSTEM
0.,0.,0.
*NODE,SYSTEM=C
9865,0.,0.,0.
9065,15.,-90.,0.
25065,30.,-90.,0.
25865,170.,-90.,0.
25965,340.,-90.,0.
25833,170.,-60.,0.
25933,340.,-60.,0.
33833,170.,-45.,0.
33933,340.,-45.,0.
41833,170.,-30.,0.
41933,340.,-30.,0.
41865,170.,0.,0.
41965,340.,0.,0.
**
**CRACK TIP REGION
*NGEN,NSET=TIP
1001,65001,1000
*NGEN,NSET=OUTER1
1033,9033,1000
*NGEN,NSET=OUTER2
9033,25033,1000
*NGEN,NSET=OUTER3
25033,33033,1000
*NGEN,NSET=OUTER4
33033,41033,1000
*NGEN,NSET=OUTER5
41033,57033,1000
*NGEN,NSET=OUTER6
57033,65033,1000
*NSET,NSET=OUTER
OUTER1,OUTER2,OUTER3,OUTER4,OUTER5,OUTER6
*NFILL,NSET=JREGION,SINGULAR=1
TIP,OUTER,2,16
**
**SECTION 9
*NGEN,NSET=BOT9
9033,9065,1
*NGEN,NSET=TOP9
9833,9865,1
*NFILL,NSET=ALL9
BOT9,TOP9,16,50
**
**SECTION 25
*NGEN,NSET=BOT25
25033,25065,1
*NGEN,NSET=TOP25,LINE=C
25833,25865,1,9865
*NFILL,NSET=ALL25,BIAS=0.8
BOT25,TOP25,16,50
**
**SECTION 41
*NGEN,NSET=BOT41
41033,41065,1
*NGEN,NSET=TOP41,LINE=C
41833,41865,1,9865
*NFILL,NSET=ALL41,BIAS=0.8
BOT41,TOP41,16,50
**
**SECTION 925
*NFILL,NSET=ALL925
BOT9,BOT25,16,1000
**
**SECTION 2541
*NSET,NSET=BOT2533,GENERATE
25033,33033,1000
*NGEN,NSET=TOP2533,LINE=C
25833,33833,1000,9865
*NSET,NSET=BOT3341,GENERATE
33033,41033,1000
*NGEN,NSET=TOP3341,LINE=C
33833,41833,1000,9865
*NSET,NSET=BOT2541
BOT2533,BOT3341
*NSET,NSET=TOP2541
TOP2533,TOP3341
*NFILL,NSET=ALL2541,BIAS=0.8
BOT2541,TOP2541,16,50
**
**SECTION 4157
*NGEN,NSET=BOT57
57033,57065,1
*NFILL,NSET=ALL4157
BOT41,BOT57,16,1000
**
**SECTION 57
*NGEN,NSET=TOP57
57833,57865,1
*NFILL,NSET=ALL57,BIAS=1.0
BOT57,TOP57,16,50
**
**SECTION 5765
*NGEN,NSET=TOP5765
57833,65833,1000
*NFILL,NSET=ALL5765,BIAS=1.0
OUTER6,TOP5765,16,50
**
**SECTION 19
*NGEN,NSET=TOP19
1833,9833,1000
*NFILL,NSET=ALL19
OUTER1,TOP19,16,50
**
**INFINITE ELEMENT REGION
*NGEN,NSET=INF25,LINE=C
25933,25965,8,9865
*NGEN,NSET=INF2533,LINE=C
25933,33933,4000,9865
*NGEN,NSET=INF3341,LINE=C
33933,41933,4000,9865
*NGEN,NSET=INF41,LINE=C
41933,41965,8,9865
*NSET,NSET=INF
INF25,INF2533,INF3341,INF41
**
**CRACK TIP REGION ELEMENTS
*ELEMENT,TYPE=CAX8R
1001,1001,1033,5033,5001,1017,3033,5017,3001
*ELGEN,ELSET=RING
1001,16,4000,4000
**
**ELEMENTS SECTION 9
*ELEMENT,TYPE=CAX8R
9033,9033,9233,9241,9041,9133,9237,9141,9037
*ELGEN,ELSET=SECT9
9033,4,8,8,4,200,200
*ELSET,ELSET=TOP9,GENERATE
9633,9657,8
**
**ELEMENTS SECTION 25
*ELEMENT,TYPE=CAX8R
25233,25233,25033,25041,25241,25133,25037,
25141,25237
*ELGEN,ELSET=SECT25
25233,4,8,8,4,200,200
**
**ELEMENTS SECTION 19
*ELEMENT,TYPE=CAX8R
1033,1033,1233,5233,5033,1133,3233,5133,3033
*ELGEN,ELSET=SECT19
1033,2,4000,4000,4,200,200
**
**ELEMENT SECTION 925
*ELEMENT,TYPE=CAX8R
13033,13033,9033,9041,13041,11033,9037,
11041,13037
*ELGEN,ELSET=SECT925
13033,4,8,8,4,4000,4000
**
**ELEMENT SECTION 2541
*ELEMENT,TYPE=CAX8R
29233,29233,29033,25033,25233,29133,27033,
25133,27233
*ELGEN,ELSET=SECT2541
29233,4,4000,4000,4,200,200
**
**ELEMENT SECTION 41
*ELEMENT, TYPE=CAX8R
41241,41241,41041,41033,41233,41141,41037,
41133,41237
*ELGEN,ELSET=SECT41
41241,4,8,8,4,200,200
**
**ELEMENT SECTION 4157
*ELEMENT,TYPE=CAX8R
41041,41041,45041,45033,41033,43041,45037,
43033,41037
*ELGEN,ELSET=SECT4157
41041,4,8,8,4,4000,4000
**
**ELEMENT SECTION 57
*ELEMENT, TYPE=CAX8R
57041,57041,57241,57233,57033,57141,57237,
57133,57037
*ELGEN,ELSET=SECT57
57041,4,8,8,4,200,200
**
**ELEMENT SECTION 5765
*ELEMENT,TYPE=CAX8R
57033,57033,57233,61233,61033,57133,59233,
61133,59033
*ELGEN,ELSET=SECT5765
57033,2,4000,4000,4,200,200
**
**INFINITE ELEMENTS
*ELEMENT,TYPE=CINAX5R
41941,41841,41833,41933,41941,41837
*ELGEN,ELSET=INF41
41941,4,8,8
*ELEMENT,TYPE=CINAX5R
29933,29833,25833,25933,29933,27833
*ELGEN,ELSET=INF2541
29933,4,4000,4000
*ELEMENT,TYPE=CINAX5R
25933,25833,25841,25941,25933,25837
*ELGEN,ELSET=INF25
25933,4,8,8
*ELSET,ELSET=INF
INF41,INF2541,INF25
**
*ELSET,ELSET=E1
RING,SECT9,SECT25,SECT41,SECT57,SECT19,SECT925,
SECT2541,SECT4157,SECT5765,INF
*MATERIAL,NAME=STEEL
*ELASTIC
30.E6,0.3
*SOLID SECTION, MATERIAL=STEEL, ELSET=E1
**
**ADD BOUNDARY CONDITIONS
*NSET,NSET=R9,GENERATE
9065,9865,50
*NSET,NSET=R925,GENERATE
9065,25065,1000
*NSET,NSET=R25,GENERATE
25065,25865,50
*NSET,NSET=YAXIS
R9,R925,R25
*NSET,NSET=L57,GENERATE
57065,57865,50
*NSET,NSET=L4157,GENERATE
41065,57065,1000
*NSET,NSET=L41,GENERATE
41065,41865,50
*NSET,NSET=XAXIS
TOP9,L57,L4157,L41
*BOUNDARY
YAXIS,XSYMM
**
**MPC'S TO TIE REDUNDANT NODES
*NSET,NSET=TIP1,GENERATE
1001,64001,1000
*NSET,NSET=TIP2,GENERATE
2001,65001,1000
*MPC
TIE,TIP1,TIP2
*MPC
TIE,57865,57833
*MPC
TIE,57833,65833
*MPC
TIE,1833,9833
**
*STEP
APPLY PRESSURE LOAD
*STATIC
1.0,1.0
*DLOAD
TOP9,P2,10.
*CONTOUR INTEGRAL,CONTOURS=1,OUTPUT=BOTH
TIP,0.707107,-0.707107
*CONTOUR INTEGRAL,CONTOURS=1,OUTPUT=BOTH,
TYPE=K FACTORS
TIP,0.707107,-0.707107
*CONTOUR INTEGRAL,CONTOURS=1,OUTPUT=BOTH,
TYPE=T-STRESS
TIP,0.707107,-0.707107
*EL PRINT
S, 
*NODE PRINT
U,RF
*NODE FILE
U, 
*OUTPUT,FIELD
*NODE OUTPUT
U,
*ENDSTEP