*HEADING
AXISYMMETRIC CONICAL CRACK MODEL--SUBMODEL WITH 13 RINGS 
 conicalcrack_axisubmodel_rms.inp
*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,14.9,-45.,0.
65001,14.9,-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.
1501,15.,-45.,0.
65501,15.,-45.,0.
*NGEN,LINE=C,NSET=TIP
1001,65001,2000,1501,15.,-45,0.0,0.0,0.0,1.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.
*NODE,SYSTEM=R
1033,15.656,-5.6579
65033,15.658,-5.6559
**
**CRACK TIP REGION
*NGEN,NSET=TIPI
1501,65501,2000
*NGEN,NSET=OUTER1
1033,9033,2000
*NGEN,NSET=OUTER2
9033,25033,2000
*NGEN,NSET=OUTER3
25033,33033,2000
*NGEN,NSET=OUTER4
33033,41033,2000
*NGEN,NSET=OUTER5
41033,57033,2000
*NGEN,NSET=OUTER6
57033,65033,2000
*NSET,NSET=OUTER
OUTER1,OUTER2,OUTER3,OUTER4,OUTER5,OUTER6
*NFILL,SINGULAR=1,NSET=JREGION
TIPI,TIP,10,-50
*NFILL,NSET=JREGION,BIAS=0.75
TIP,OUTER,16,2
**
**CRACK TIP REGION ELEMENTS
*ELEMENT,TYPE=CAX8R
1001,1001,1005,5005,5001,1003,3005,5003,3001
1,1001,5001,5101,1101,3001,5051,3101,1051
*ELGEN,ELSET=RINGS
1001,16,4000,4000,8,4,4
1,16,4000,1,5,100,16
**
*MATERIAL,NAME=STEEL
*ELASTIC
30.E6,0.3
*SOLID SECTION, MATERIAL=STEEL, ELSET=RINGS
*ELSET,ELSET=TIPCORE,GEN
1,80
**
**MPC'S TO TIE REDUNDANT NODES
*NSET,NSET=TIP1,GENERATE
1501,63501,2000
*NSET,NSET=TIP2,GENERATE
3501,65501,2000
*MPC
TIE,TIP1,TIP2
*SUBMODEL
OUTER, 
**
*STEP
*STATIC
1.,1.
*BOUNDARY,SUBMODEL,STEP=1
OUTER,1,2
*CONTOUR INTEGRAL,CONTOURS=5,OUTPUT=BOTH
TIPI,0.707107,-0.707107
*CONTOUR INTEGRAL,CONTOURS=5,OUTPUT=BOTH,TYPE=K FACTORS
TIPI,0.707107,-0.707107
*CONTOUR INTEGRAL,CONTOURS=5,OUTPUT=BOTH,TYPE=T-STRESS
TIPI,0.707107,-0.707107
*EL PRINT
S, 
*NODE PRINT
U,RF
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=TIPCORE
E,S 
*ENDSTEP