*HEADING
J-INTEGRAL FOR SEMI-ELLIPTICAL SURFACE CRACK        
IN A SEMI-INFINITE BODY
  A/C=.4      A/T=.1
     A,  MINOR CRACK DEPTH
     T,  MINOR BODY THICKNESS
     C,  MAJOR CRACK WIDTH
*PREPRINT,HISTORY=NO,MODEL=NO
*RESTART,WRITE
*NODE,INPUT=jintegral3d_node.inp
*NSET,NSET=XYPL,GENERATE
10,2090,130
70,90,10
90,110,10
110,130,10
590,610,10
610,630,10
630,650,10
70,590,130
80,600,130
90,610,130
100,620,130
110,630,130
120,640,130
130,650,130
590,1630,130
1630,2150,130
2200,2650,150
2280,2730,150
2200,2280,10
2280,2300,10
2300,2320,10
2320,2340,10
2350,2430,10
2500,2580,10
2650,2730,10
2730,2750,10
2750,2770,10
2770,2790,10
2290,2740,150
2300,2750,150
2310,2760,150
2320,2770,150
2330,2780,150
2340,2790,150
*NSET,NSET=XYPL,GENERATE
10,70,10
140,200,10
270,330,10
400,460,10
530,590,10
660,720,10
790,850,10
920,980,10
1050,1110,10
1180,1240,10
1310,1370,10
1440,1500,10
1570,1630,10
1700,1760,10
1830,1890,10
1960,2020,10
2090,2150,10
*NSET,NSET=YZPL,GENERATE
18,2098,130
78,98,10
98,118,10
118,138,10
598,618,10
618,638,10
638,658,10
78,598,130
88,608,130
98,618,130
108,628,130
118,638,130
128,648,130
138,658,130
598,1638,130
1638,2158,130
2208,2658,150
2288,2738,150
2208,2288,10
2288,2308,10
2308,2328,10
2328,2348,10
2358,2438,10
2508,2588,10
2658,2738,10
2738,2758,10
2758,2778,10
2778,2798,10
2298,2748,150
2308,2758,150
2318,2768,150
2328,2778,150
2338,2788,150
2348,2798,150
*NSET,NSET=YZPL,GENERATE
18,78,10
148,208,10
278,338,10
408,468,10
538,598,10
668,728,10
798,858,10
928,988,10
1058,1118,10
1188,1248,10
1318,1378,10
1448,1508,10
1578,1638,10
1708,1768,10
1838,1898,10
1968,2028,10
2098,2158,10
*NSET,NSET=RAD0,GENERATE
10,18,1
20,28,1
30,38,1
40,48,1
50,58,1
60,68,1
70,78,1
*NSET,NSET=RAD1,GENERATE
140,148,1
160,168,1
180,188,1
200,208,1
*NSET,NSET=RAD2,GENERATE
270,278,1
280,288,1
290,298,1
300,308,1
310,318,1
320,328,1
330,338,1
*NSET,NSET=RAD3,GENERATE
400,408,1
420,428,1
440,448,1
460,468,1
*NSET,NSET=RAD4,GENERATE
530,538,1
540,548,1
550,558,1
560,568,1
570,578,1
580,588,1,
590,598,1
*NSET,NSET=RAD5,GENERATE
660,668,1
680,688,1
700,708,1
720,728,1
*NSET,NSET=RAD6,GENERATE
790,798,1
800,808,1
810,818,1
820,828,1
830,838,1
840,848,1,
850,858,1
*NSET,NSET=RAD7,GENERATE
920,928,1
940,948,1
960,968,1
980,988,1
*NSET,NSET=RAD8,GENERATE
1050,1058,1
1060,1068,1
1070,1078,1
1080,1088,1
1090,1098,1
1100,1108,1,
1110,1118,1
*NSET,NSET=RAD9,GENERATE
1180,1188,1
1200,1208,1
1220,1228,1
1240,1248,1
*NSET,NSET=RADA,GENERATE
1310,1318,1
1320,1328,1
1330,1338,1
1340,1348,1
1350,1358,1
1360,1368,1,
1370,1378,1
*NSET,NSET=RADB,GENERATE
1440,1448,1
1460,1468,1
1480,1488,1
1500,1508,1
*NSET,NSET=RADC,GENERATE
1570,1578,1
1580,1588,1
1590,1598,1
1600,1608,1
1610,1618,1
1620,1628,1,
*NSET,NSET=RADD,GENERATE
1700,1708,1
1720,1728,1
1740,1748,1
*NSET,NSET=RADE,GENERATE
1830,1838,1
1840,1848,1
1850,1858,1
1860,1868,1
1870,1878,1
1880,1888,1
*NSET,NSET=RADF,GENERATE
1960,1968,1
1980,1988,1
2000,2008,1
*NSET,NSET=RADG,GENERATE
2090,2098,1
2100,2108,1
2110,2118,1
2120,2128,1
2130,2138,1
2140,2148,1
*NSET,NSET=LEV0,GENERATE
80,88,1
90,98,1
100,108,1
110,118,1
120,128,1
130,138,1
*NSET,NSET=LEV1,GENERATE
220,228,1
240,248,1
260,268,1
*NSET,NSET=LEV2,GENERATE
340,348,1
350,358,1
360,368,1
370,378,1
380,388,1
390,398,1
*NSET,NSET=LEV3,GENERATE
480,488,1
490,498,1
500,508,1
520,528,1
*NSET,NSET=LEV4,GENERATE
600,608,1
610,618,1
620,628,1
630,638,1
640,648,1
650,658,1
*NSET,NSET=LEV5,GENERATE
2220,2228,1
2240,2248,1
2260,2268,1
2280,2288,1
2300,2308,1
2320,2328,1
2340,2348,1
*NSET,NSET=LEV6,GENERATE
2360,2368,1
2370,2378,1
2380,2388,1
2390,2398,1
2400,2408,1
2410,2418,1
2420,2428,1
2430,2438,1
2440,2448,1
2450,2458,1
2460,2468,1
2470,2478,1
2480,2488,1
2490,2498,1
*NSET,NSET=LEV7,GENERATE
2520,2528,1
2540,2548,1
2560,2568,1
2580,2588,1
2600,2608,1
2620,2628,1
2640,2648,1
*NSET,NSET=LEV8,GENERATE
2660,2668,1
2670,2678,1
2680,2688,1
2690,2698,1
2700,2708,1
2710,2718,1
2720,2728,1
2730,2738,1
2740,2748,1
2750,2758,1
2760,2768,1
2770,2778,1
2780,2788,1
2790,2798,1
*NSET,NSET=J0,GENERATE
10,2090,130
*NSET,NSET=J1,GENERATE
11,2091,260
*NSET,NSET=J2,GENERATE
12,2092,130
*NSET,NSET=J3,GENERATE
13,2093,260
*NSET,NSET=J4,GENERATE
14,2094,130
*NSET,NSET=J5,GENERATE
15,2095,260
*NSET,NSET=J6,GENERATE
16,2096,130
*NSET,NSET=J7,GENERATE
17,2097,260
*NSET,NSET=J8,GENERATE
18,2098,130
*NSET,NSET=M0,GENERATE
140,2090,130
*NSET,NSET=M1,GENERATE
271,2091,260
*NSET,NSET=M2,GENERATE
142,2092,130
*NSET,NSET=M3,GENERATE
273,2093,260
*NSET,NSET=M4,GENERATE
144,2094,130
*NSET,NSET=M5,GENERATE
275,2095,260
*NSET,NSET=M6,GENERATE
146,2096,130
*NSET,NSET=M7,GENERATE
277,2097,260
*NSET,NSET=M8,GENERATE
148,2098,130
*MPC
PIN,M0,10
PIN,M1,11
PIN,M2,12
PIN,M3,13
PIN,M4,14
PIN,M5,15
PIN,M6,16
PIN,M7,17
PIN,M8,18
*NSET,NSET=NALL
XYPL,YZPL,
RAD0,RAD1,RAD2,RAD3,RAD4,RAD5,RAD6,RAD7,RAD8,RAD9,
RADA,RADB,RADC,RADD,RADE,RADF,RADG,
LEV0,LEV1,LEV2,LEV3,LEV4,LEV5,LEV6,LEV7,LEV8,
J0,J1,J2,J3,J4,J5,J6,J7,J8
*ELEMENT,TYPE=C3D20R
100,12,32,292,272,10,30,290,270,22,162,282,142,20,160,280,
140,11,31,291,271
*ELGEN,ELSET=MAST
100,8,260,100
*ELGEN,ELSET=FAN1
100,6,20,10
100,4,2,1
110,4,2,1
120,4,2,1
130,4,2,1
140,4,2,1
150,4,2,1
150,2,260,10
160,4,2,1
*ELEMENT,TYPE=C3D20R,ELSET=FAN1
170,632,652,2492,2472,630,650,2490,2470,642,2342,2482,2322,640,
2340,2480,2320,631,651,2491,2471
180,2472,2492,2792,2772,2470,2490,2790,2770,2482,2642,2782,2622,
2480,2640,2780,2620,2471,2491,2791,2771
*ELGEN,ELSET=FAN1
170,4,2,1
180,4,2,1
*ELGEN
170,2,-20,80
180,2,-20,80
250,2,-20,20
260,2,-20,20
*ELGEN,ELSET=FAN2
200,5,20,10
200,4,2,1
210,4,2,1
220,4,2,1
230,4,2,1
240,4,2,1
250,4,2,1
260,4,2,1
270,4,2,1
280,4,2,1
*ELGEN,ELSET=FAN3
300,3,20,10
300,4,2,1
310,4,2,1
320,4,2,1
*ELEMENT,TYPE=C3D20R,ELSET=FAN3
330,852,592,2432,2412,850,590,2430,2410,722,2282,2422,2262,720,
2280,2420,2260,851,591,2431,2411
340,2412,2432,2732,2712,2410,2430,2730,2710,2422,2582,2722,2562,
2420,2580,2720,2560,2411,2431,2731,2711
*ELGEN,ELSET=FAN3
330,4,2,1
340,4,2,1
*ELGEN,ELSET=FAN4
400,3,20,10
400,4,2,1
410,4,2,1
420,4,2,1
*ELEMENT,TYPE=C3D20R,ELSET=FAN4
430,1112,852,2412,2392,1110,850,2410,2390,982,2262,2402,2242,980,
2260,2400,2240,1111,851,2411,2391
440,2392,2412,2712,2692,2390,2410,2710,2690,2402,2562,2702,2542,
2400,2560,2700,2540,2391,2411,2711,2691
*ELGEN,ELSET=FAN4
430,4,2,1
440,4,2,1
*ELGEN,ELSET=FAN5
500,3,20,10
500,4,2,1
510,4,2,1
520,4,2,1
*ELEMENT,TYPE=C3D20R,ELSET=FAN5
530,1372,1112,2392,2372,1370,1110,2390,2370,1242,2242,2382,2222,
1240,2240,2380,2220,1371,1111,2391,2371
540,2372,2392,2692,2672,2370,2390,2690,2670,2382,2542,2682,2522,
2380,2540,2680,2520,2371,2391,2691,2671
*ELGEN,ELSET=FAN5
530,4,2,1
540,4,2,1
*ELGEN,ELSET=FAN6
600,3,20,10
600,4,2,1
610,4,2,1
620,4,2,1
*ELEMENT,TYPE=C3D20R,ELSET=FAN6
630,1632,1372,2372,2352,1630,1370,2370,2350,1502,2222,2362,2202,
1500,2220,2360,2200,1631,1371,2371,2351
640,2352,2372,2672,2652,2350,2370,2670,2650,2362,2522,2662,2502,
2360,2520,2660,2500,2351,2371,2671,2651
*ELGEN,ELSET=FAN6
630,4,2,1
640,4,2,1
*ELGEN,ELSET=FAN7
700,3,20,10
700,4,2,1
710,4,2,1
720,4,2,1
*ELGEN,ELSET=FAN8
800,3,20,10
800,4,2,1
810,4,2,1
820,4,2,1
*ELSET,ELSET=LOAD,GENERATE
640,643
540,543
440,443
340,343
280,283
260,263
180,183
*ELSET,ELSET=PRINT
120,121,122,123,150,151,152,153
*ELSET,ELSET=PLOTA
LOAD, 
*ELSET,ELSET=PLOTA,GENERATE
100,820,10
*ELSET,ELSET=PLOTB,GENERATE
100,120,10
200,220,10
300,320,10
400,420,10
500,520,10
600,620,10
700,720,10
800,820,10
*ELSET,ELSET=EALL,GENERATE
1,850
*surface, NAME=SURFACE
 LOAD, S5
*MATERIAL,NAME=MAT
*ELASTIC
30.E6,.3
*SOLID SECTION,MATERIAL=MAT,ELSET=EALL
*BOUNDARY
RAD0,2
LEV0,2
138,3
YZPL,1
**
*STEP
*STATIC
*DSLOAD
SURFACE,P,-100.
*CONTOUR INTEGRAL,CONTOURS=3,SYMM,OUTPUT=BOTH,NORMAL
0.0,1.0,0.0
J0,J1,J2,J3,J4,J5,J6,J7,J8
*CONTOUR INTEGRAL,CONTOURS=3,SYMM,OUTPUT=BOTH,NORMAL,TYPE=K FACTORS
0.0,1.0,0.0
J0,J1,J2,J3,J4,J5,J6,J7,J8
*CONTOUR INTEGRAL,CONTOURS=3,SYMM,OUTPUT=BOTH,NORMAL,TYPE=T-STRESS
0.0,1.0,0.0
J0,J1,J2,J3,J4,J5,J6,J7,J8
*EL PRINT,ELSET=PRINT,POSITION=CENTROID
S,MISES
*ENERGY PRINT
*NODE PRINT
U,
*ENERGY FILE
*OUTPUT,FIELD, OP=NEW
*ELEMENT OUTPUT, ELSET=PRINT
S, 
*NODE OUTPUT,NSET=NALL
U, 
*OUTPUT,HISTORY,OP=NEW
*ENERGY OUTPUT,VARIABLE=ALL
*ELEMENT OUTPUT, ELSET=PRINT
S, 
*NODE OUTPUT,NSET=NALL
U, 
*END STEP