*HEADING
 AXISYMMETRIC ANALYSIS OF PENNY SHAPED CRACK  
 IN ROUND BAR  A/R=.5
*SYSTEM
10., 
*NODE
1, 
7, 10.
87,10., 10.
247, -10., 10.
327, -10.
91, 10., 30.
251, -10.,  30.
321, 
*NGEN
1,321, 20
7,87,20
247,327,20
*NGEN,NSET=END1
87,247,20
*NGEN,NSET=END2
91,251,20
*NFILL
END1,END2,4,1
*NGEN,LINE=P
1,7,,,2.5
21,27,,,2.5,.625
41,47,,,2.5,1.25
61,67,,,2.5,1.875
81,87,,,2.5,2.5
101,107,,,1.875,2.5
121,127,,,1.25,2.5
141,147,,,.625,2.5
161,167,,,0.,2.5
181,187,,,-.625,2.5
201,207,,,-1.25,2.5
221,227,,,-1.875,2.5
241,247,,,-2.5,2.5
261,267,,,-2.5,1.875
281,287,,,-2.5,1.25
301,307,,,-2.5,.625
321,327,,,-2.5,0.
*ELEMENT,TYPE=CAX8
1,  1,3,43,41,  2,23,42,21
*ELGEN,ELSET=ALL
1,3,2,1,8,40,10
23,3,2,1,4,40,10
*SOLID SECTION,ELSET=ALL,MATERIAL=ELAS
*NSET,NSET=XAXIS,GENERATE
1,7
*NSET,NSET=YAXIS,GENERATE
247,251
247,327,20
*ELSET,ELSET=PULL
25,35,45,55
*NSET,NSET=CRACK,GENERATE
1,321,20
*ELSET,ELSET=DRAW,GENERATE
1,71,10
2,72,10
3,73,10
*MATERIAL,NAME=ELAS
*ELASTIC
30.E6, .3
*BOUNDARY
XAXIS,2
YAXIS,1
*STEP
*STATIC
*DLOAD
PULL, P2,-100.
*CONTOUR INTEGRAL,FREQ=1,CONTOURS=3,SYMM,OUTPUT=BOTH
CRACK, 1.0, 0.0
*CONTOUR INTEGRAL,FREQ=1,CONTOURS=3,SYMM,OUTPUT=BOTH,TYPE=K FACTORS
CRACK, 1.0, 0.0
*CONTOUR INTEGRAL,FREQ=1,CONTOURS=3,SYMM,OUTPUT=BOTH,TYPE=T-STRESS
CRACK, 1.0, 0.0
*ENERGY PRINT
*EL PRINT
S,E
*NODE PRINT
U,RF
*NODE FILE
U,
*OUTPUT,FIELD
*NODE OUTPUT
U,
*END STEP