*HEADING
 THREE-DIMENSIONAL CONICAL CRACK MODEL USING REFINED MESHES 
 --SUBMODEL WITH 13 RINGS
 File name: conicalcrack_3dsubmodel_rms.inp
 uses the results file from conicalcrack_3dglobal.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,1000,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
*NGEN,NSET=TIPI
1501,65501,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,SINGULAR=1,NSET=J
TIPI,TIP,10,-50
*NFILL,NSET=J,BIAS=0.75
TIP,OUTER,16,2
*NSET,NSET=ALL0,GENERATE
1033,65033,4000
*NSET,NSET=EVEN0,GENERATE
3033,63033,4000
**
*NSET,NSET=N0
J, 
*NCOPY,CHANGE NUMBER=100000,OLD SET=N0,MULTIPLE=19,SHIFT,NEW SET=N119
0.,0.,0.
0.,0.,0.,0.,1.,0.,-4.5
*NCOPY,CHANGE NUMBER=2000000,OLD SET=N0,MULTIPLE=1,SHIFT,NEW SET=N20
0.,0.,0.
0.,0.,0.,0.,1.,0.,-90.
*NCOPY,CHANGE NUMBER=100000,OLD SET=ALL0,MULTIPLE=19,SHIFT,NEW SET=ALL119
0.,0.,0.
0.,0.,0.,0.,1.,0.,-4.5
*NCOPY,CHANGE NUMBER=2000000,OLD SET=ALL0,MULTIPLE=1,SHIFT,NEW SET=ALL20
0.,0.,0.
0.,0.,0.,0.,1.,0.,-90.
*NCOPY,CHANGE NUMBER=200000,OLD SET=EVEN0,MULTIPLE=9,SHIFT,NEW SET=EVEN218
0.,0.,0.
0.,0.,0.,0.,1.,0.,-9.
*NCOPY,CHANGE NUMBER=2000000,OLD SET=EVEN0,MULTIPLE=1,SHIFT,NEW SET=EVEN20
0.,0.,0.
0.,0.,0.,0.,1.,0.,-90.
**
*ELEMENT,TYPE=C3D20R
1001,1001,1005,5005,5001,201001,201005,205005,205001,1003,3005,5003,3001,
201003,203005,205003,203001,101001,101005,105005,105001
1,1001,5001,5101,1101,201001,205001,205101,201101,3001,5051,3101,1051,
203001,205051,203101,201051,101001,105001,105101,101101
*ELGEN,ELSET=RINGS
1001,16,4000,4000,8,4,4
1,16,4000,1,5,100,16
**
*ELSET,ELSET=E1
RINGS, 
*ELCOPY,ELEMENT SHIFT=200000,OLD SET=E1,SHIFT NODES=200000,NEW SET=E2
*ELCOPY,ELEMENT SHIFT=200000,OLD SET=E2,SHIFT NODES=200000,NEW SET=E3
*ELCOPY,ELEMENT SHIFT=200000,OLD SET=E3,SHIFT NODES=200000,NEW SET=E4
*ELCOPY,ELEMENT SHIFT=200000,OLD SET=E4,SHIFT NODES=200000,NEW SET=E5
*ELCOPY,ELEMENT SHIFT=200000,OLD SET=E5,SHIFT NODES=200000,NEW SET=E6
*ELCOPY,ELEMENT SHIFT=200000,OLD SET=E6,SHIFT NODES=200000,NEW SET=E7
*ELCOPY,ELEMENT SHIFT=200000,OLD SET=E7,SHIFT NODES=200000,NEW SET=E8
*ELCOPY,ELEMENT SHIFT=200000,OLD SET=E8,SHIFT NODES=200000,NEW SET=E9
*ELCOPY,ELEMENT SHIFT=200000,OLD SET=E9,SHIFT NODES=200000,NEW SET=E10
*ELSET,ELSET=EALL
E1,E2,E3,E4,E5,E6,E7,E8,E9,E10
*MATERIAL,NAME=STEEL
*ELASTIC
30.E6,0.3
*SOLID SECTION, MATERIAL=STEEL, ELSET=EALL
*ELSET,ELSET=TIPCORE,GEN
1,80
200001,200080
400001,400080
600001,600080
800001,800080
1000001,1000080
1200001,1200080
1400001,1400080
1600001,1600080
1800001,1800080
*NSET,NSET=T0,GENERATE
1501,65501,1000
*NSET,NSET=T1,GENERATE
101501,165501,1000
*NSET,NSET=T2,GENERATE
201501,265501,1000
*NSET,NSET=T3,GENERATE
301501,365501,1000
*NSET,NSET=T4,GENERATE
401501,465501,1000
*NSET,NSET=T5,GENERATE
501501,565501,1000
*NSET,NSET=T6,GENERATE
601501,665501,1000
*NSET,NSET=T7,GENERATE
701501,765501,1000
*NSET,NSET=T8,GENERATE
801501,865501,1000
*NSET,NSET=T9,GENERATE
901501,965501,1000
*NSET,NSET=T10,GENERATE
1001501,1065501,1000
*NSET,NSET=T11,GENERATE
1101501,1165501,1000
*NSET,NSET=T12,GENERATE
1201501,1265501,1000
*NSET,NSET=T13,GENERATE
1301501,1365501,1000
*NSET,NSET=T14,GENERATE
1401501,1465501,1000
*NSET,NSET=T15,GENERATE
1501501,1565501,1000
*NSET,NSET=T16,GENERATE
1601501,1665501,1000
*NSET,NSET=T17,GENERATE
1701501,1765501,1000
*NSET,NSET=T18,GENERATE
1801501,1865501,1000
*NSET,NSET=T19,GENERATE
1901501,1965501,1000
*NSET,NSET=T20,GENERATE
2001501,2065501,1000
**
*NSET,NSET=TIP01,GENERATE
1501,64501,1000
*NSET,NSET=TIP11,GENERATE
101501,164501,1000
*NSET,NSET=TIP21,GENERATE
201501,264501,1000
*NSET,NSET=TIP31,GENERATE
301501,364501,1000
*NSET,NSET=TIP41,GENERATE
401501,464501,1000
*NSET,NSET=TIP51,GENERATE
501501,564501,1000
*NSET,NSET=TIP61,GENERATE
601501,664501,1000
*NSET,NSET=TIP71,GENERATE
701501,764501,1000
*NSET,NSET=TIP81,GENERATE
801501,864501,1000
*NSET,NSET=TIP91,GENERATE
901501,964501,1000
*NSET,NSET=TIP101,GENERATE
1001501,1064501,1000
*NSET,NSET=TIP111,GENERATE
1101501,1164501,1000
*NSET,NSET=TIP121,GENERATE
1201501,1264501,1000
*NSET,NSET=TIP131,GENERATE
1301501,1364501,1000
*NSET,NSET=TIP141,GENERATE
1401501,1464501,1000
*NSET,NSET=TIP151,GENERATE
1501501,1564501,1000
*NSET,NSET=TIP161,GENERATE
1601501,1664501,1000
*NSET,NSET=TIP171,GENERATE
1701501,1764501,1000
*NSET,NSET=TIP181,GENERATE
1801501,1864501,1000
*NSET,NSET=TIP191,GENERATE
1901501,1964501,1000
*NSET,NSET=TIP201,GENERATE
2001501,2064501,1000
**
*NSET,NSET=TIP02,GENERATE
2501,65501,1000
*NSET,NSET=TIP12,GENERATE
102501,165501,1000
*NSET,NSET=TIP22,GENERATE
202501,265501,1000
*NSET,NSET=TIP32,GENERATE
302501,365501,1000
*NSET,NSET=TIP42,GENERATE
402501,465501,1000
*NSET,NSET=TIP52,GENERATE
502501,565501,1000
*NSET,NSET=TIP62,GENERATE
602501,665501,1000
*NSET,NSET=TIP72,GENERATE
702501,765501,1000
*NSET,NSET=TIP82,GENERATE
802501,865501,1000
*NSET,NSET=TIP92,GENERATE
902501,965501,1000
*NSET,NSET=TIP102,GENERATE
1002501,1065501,1000
*NSET,NSET=TIP112,GENERATE
1102501,1165501,1000
*NSET,NSET=TIP122,GENERATE
1202501,1265501,1000
*NSET,NSET=TIP132,GENERATE
1302501,1365501,1000
*NSET,NSET=TIP142,GENERATE
1402501,1465501,1000
*NSET,NSET=TIP152,GENERATE
1502501,1565501,1000
*NSET,NSET=TIP162,GENERATE
1602501,1665501,1000
*NSET,NSET=TIP172,GENERATE
1702501,1765501,1000
*NSET,NSET=TIP182,GENERATE
1802501,1865501,1000
*NSET,NSET=TIP192,GENERATE
1902501,1965501,1000
*NSET,NSET=TIP202,GENERATE
2002501,2065501,1000
*MPC
TIE,TIP01,TIP02
TIE,TIP11,TIP12
TIE,TIP21,TIP22
TIE,TIP31,TIP32
TIE,TIP41,TIP42
TIE,TIP51,TIP52
TIE,TIP61,TIP62
TIE,TIP71,TIP72
TIE,TIP81,TIP82
TIE,TIP91,TIP92
TIE,TIP101,TIP102
TIE,TIP111,TIP112
TIE,TIP121,TIP122
TIE,TIP131,TIP132
TIE,TIP141,TIP142
TIE,TIP151,TIP152
TIE,TIP161,TIP162
TIE,TIP171,TIP172
TIE,TIP181,TIP182
TIE,TIP191,TIP192
TIE,TIP201,TIP202
****
*SUBMODEL
ALL0,EVEN0,ALL119,EVEN218,ALL20,EVEN20
*BOUNDARY
N0,3,3
N20,1,1
*STEP
*STATIC
1.0,1.0
*BOUNDARY,SUBMODEL,STEP=1
ALL0,1,2
EVEN0,1,2
ALL20,2,3
EVEN20,2,3
ALL119,1,3
EVEN218,1,3
*CONTOUR INTEGRAL,CONTOURS=13,OUTPUT=BOTH
T0,0.70710678,-0.70710678,0.
T1,0.70492701,-0.70710678,0.055478959
T2,0.69840112,-0.70710678,0.11061587
T3,0.68756936,-0.70710678,0.1650708
T4,0.67249851,-0.70710678,0.21850801
T5,0.65328148,-0.70710678,0.27059805
T6,0.63003676,-0.70710678,0.32101976
T7,0.60290764,-0.70710678,0.36946228
T8,0.5720614,-0.70710678,0.41562694
T9,0.53768821,-0.70710678,0.45922912
T10,0.5,-0.70710678,0.5
T11,0.45922912,-0.70710678,0.53768821
T12,0.41562694,-0.70710678,0.5720614
T13,0.36946228,-0.70710678,0.60290764
T14,0.32101976,-0.70710678,0.63003676
T15,0.27059805,-0.70710678,0.65328148
T16,0.21850801,-0.70710678,0.67249851
T17,0.1650708,-0.70710678,0.68756936
T18,0.11061587,-0.70710678,0.69840112
T19,0.055478959,-0.70710678,0.70492701
T20,0.,-0.70710678,0.70710678
*CONTOUR INTEGRAL,CONTOURS=13,OUTPUT=BOTH,TYPE=K FACTORS
T0,0.70710678,-0.70710678,0.
T1,0.70492701,-0.70710678,0.055478959
T2,0.69840112,-0.70710678,0.11061587
T3,0.68756936,-0.70710678,0.1650708
T4,0.67249851,-0.70710678,0.21850801
T5,0.65328148,-0.70710678,0.27059805
T6,0.63003676,-0.70710678,0.32101976
T7,0.60290764,-0.70710678,0.36946228
T8,0.5720614,-0.70710678,0.41562694
T9,0.53768821,-0.70710678,0.45922912
T10,0.5,-0.70710678,0.5
T11,0.45922912,-0.70710678,0.53768821
T12,0.41562694,-0.70710678,0.5720614
T13,0.36946228,-0.70710678,0.60290764
T14,0.32101976,-0.70710678,0.63003676
T15,0.27059805,-0.70710678,0.65328148
T16,0.21850801,-0.70710678,0.67249851
T17,0.1650708,-0.70710678,0.68756936
T18,0.11061587,-0.70710678,0.69840112
T19,0.055478959,-0.70710678,0.70492701
T20,0.,-0.70710678,0.70710678
*CONTOUR INTEGRAL,CONTOURS=13,OUTPUT=BOTH,TYPE=T-STRESS
T0,0.70710678,-0.70710678,0.
T1,0.70492701,-0.70710678,0.055478959
T2,0.69840112,-0.70710678,0.11061587
T3,0.68756936,-0.70710678,0.1650708
T4,0.67249851,-0.70710678,0.21850801
T5,0.65328148,-0.70710678,0.27059805
T6,0.63003676,-0.70710678,0.32101976
T7,0.60290764,-0.70710678,0.36946228
T8,0.5720614,-0.70710678,0.41562694
T9,0.53768821,-0.70710678,0.45922912
T10,0.5,-0.70710678,0.5
T11,0.45922912,-0.70710678,0.53768821
T12,0.41562694,-0.70710678,0.5720614
T13,0.36946228,-0.70710678,0.60290764
T14,0.32101976,-0.70710678,0.63003676
T15,0.27059805,-0.70710678,0.65328148
T16,0.21850801,-0.70710678,0.67249851
T17,0.1650708,-0.70710678,0.68756936
T18,0.11061587,-0.70710678,0.69840112
T19,0.055478959,-0.70710678,0.70492701
T20,0.,-0.70710678,0.70710678
*EL PRINT,freq=0
S, 
*NODE PRINT,freq=0
U,RF
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=TIPCORE
E,S 
*END STEP