Comparison of The Stress Intensity Factor Influence Coefficients For Axial Id Surface Cracks in Cylinders of Rse-M and Api 579-1
Comparison of The Stress Intensity Factor Influence Coefficients For Axial Id Surface Cracks in Cylinders of Rse-M and Api 579-1
PVP2015
July 19-23, 2015, Boston, Massachusetts, USA
PVP2015-45236
ABSTRACT NOMENCLATURE
Analytical evaluation procedures for determining the a Crack depth
acceptability of flaw detected during in-service inspection of Aj Coefficients of polynomial representation of stress
nuclear power plant components are provided in Appendix 5.4 distribution over crack depth
of the French RSE-M Code. Linear elastic fracture mechanics Bj Coefficients of polynomial representation of stress
based evaluation procedures require calculation of the stress distribution through the wall thickness
intensity factor (SIF). In Appendix 5.4 of the RSE-M Code, c One-half of the crack length
influence coefficients needed to compute the SIF are provided Gj Influence coefficients for the stress intensity factor
for a wide range of surface axial or circumferential flaws in ij Influence coefficients for the stress intensity factor
cylinders, the through-wall stress field being represented by a KI Stress intensity factor
cubic equation. On the other hand, Appendix C of API 579-1 Q Flaw shape parameter
FFS procedure provides also a very complete set of influence Ri Inner radius of a cylinder
coefficients. t Wall thickness
The paper presents the comparison of the influence (x) Actual stress distribution
coefficients from both documents, focused on axial ID semi- (x)A Polynomial representation of stress distribution over crack
elliptical surface flaws in cylinders. The cylinder and crack depth
geometries are represented by three ratios: Ri/t, a/t, and a/c, (x)B Polynomial representation of stress distribution through
where Ri, t, a, and c are respectively the inner radius, the wall wall thickness
thickness, the crack depth and one-half of the crack length. The
solutions for the coefficients G0 and G1 at the deepest point and INTRODUCTION
at the surface point are investigated.
At the deepest point, the agreement between the solutions is The influence coefficients given in Appendix 5.4 of the
good, the relative difference being lower than 2 %, except for RSE-M Code were developed at the end of the 1990’s thanks to
the plate (Ri/t = ) at a/c = 0.125 and 0.0625 and a/t = 0.8 a very large finite element computation grid. Appendix C of API
(around 5 %). At the surface point, the agreement between both 579-1 FFS procedure provides also a very complete set of
solutions is not so good. At this point, the relative differences influence coefficients, developed by Anderson for the Welding
depend strongly on the a/c ratio, being larger for elongated Research Committee.
cracks (with low a/c ratios). However, it must be recalled that The paper presents the comparison of the influence
the absolute values of the coefficients are low at the surface coefficients from both documents, focused on axial ID semi-
point for elongated cracks, and that for these cracks the critical elliptical surface cracks in cylinders (Fig. 1). The solutions at the
point regarding the stress intensity factor is the deepest point. deepest point and at the surface point are investigated. The
comparison is made directly on the influence coefficients, since
In the RSE-M Code, the SIF KI is calculated using a slightly More precisely, the following geometries are covered:
different relation: Ri/t ratios: 1, 3, 5, 10, 20, 60, 100,
3 a/t ratios: 0, 0.2, 0.4, 0.6, 0.8
K I A ji j a (3) a/c ratios: 0, 0.03125, 0.0625, 0.125, 0.25, 0.5, 1, 2
j0
where ij is an influence coefficient relative to the exponent “j” of The case of a/c = 0 corresponds to axially infinite ID
the loading. It is obvious that the two kinds of influence surface crack and the case of Ri/t = to the plate.
coefficients are linked by the relation: Approximate higher order influence coefficients G2 to G4
Gj ij Q (4) can be computed from G0 and G1 [7].
For each Ri/t ratio selected, the coefficients G0 and G1 at the where E and are respectively the Young’s modulus and
deepest point and at the surface point are compared graphically the Poisson’s ratio (200,000 MPa and 0.3).
on Fig. 4 to Fig. 23, as a function of the a/t ratio. One can The evolution of the coefficients G0 and G1 along the crack
observe that at the deepest point, the values are very close front is presented in Fig. 25 and Fig. 26. One can see that the
between API and RSE-M. At the surface point, some differences agreement between API-579 (equation (8)) and the FE solution
are evident between API and RSE-M solutions: is very good all along the crack front, except near the surface
point for the G0 coefficient. The sixth order polynomial fit used
for a/c = 0.0625, the RSE-M values are above the API
by API-579 cannot handle precisely the evolution of G0 near the
ones, whereas for a/c = 1 the RSE-M values are below
surface.
the API ones,
CONCLUSIONS
for the Ri/t ratios of 10, 20 and , and a/c = 1, the API
values show an irregular variation. The influence coefficients G0 and G1 for axial ID semi-
elliptical surface cracks in cylinders from API 579-1 and RSE-M
However, it should be recalled first that the computation of
Code Appendix 5.4 have been compared.
the stress intensity factor (and hence the computation of the
influence coefficient) at the surface point is not easy and At the deepest point, the agreement between the solutions is
accurate and second, that the values of G0 and G1 at this point good, the relative difference being lower than 2 %, except for
are generally lower than at the deepest point, particularly for G 1 the plate (Ri/t = ) at a/c = 0.125 and 0.0625 and a/t = 0.8
(always lower than 0.3). The plane strain condition that prevails (around 5 %).
along the crack front far from the free surface vanishes as the
At the surface point, the agreement between both solutions
point of the crack front comes closer to this surface [10].
is not so good. At this point, the differences depend strongly on
More precisely, the relative difference between the two the a/c ratio. They are generally below 5 % for a/c = 1, 0.5 and
solutions was evaluated by: 0.25, close to 10 % for a/c = 0.125 and close to 20 % for
a/c = 0.0625. The maximum relative difference (-32 %) is
G API observed for a/c = 0.0625 (0.014 against 0.020). However, it
Re lative difference (%) 1 x 100 (9)
G RSEM must be recalled that the absolute values of the coefficients are
low at the surface point for elongated cracks, and that for these
REFERENCES
1. “Rules for In-service Inspection of Nuclear Power Plant
Components”, Section XI, Division 1, ASME Boiler and
Pressure Vessel Code, 2013 Edition
2. “Rules for In-service Inspection of Nuclear Power Plant
Components”, 2010 Edition and 2013 Addendum, AFCEN,
Paris
3. “Proposed Change and Rewrite of K I Calculation Procedure
in Appendix A”, ASME Section XI Working Group on
Flaw Evaluation, February 11, 2014
4. Fitness-for-Service, API 579 Second Edition, API 579-
1/ASME FFS-1, June 5, 1997
5. Xu, S.X., Darrell, D.R., et al, 2014, “Closed-form relations
for stress intensity factor influence coefficients for axial ID
surface flaws in cylinders for Appendix A of ASME Section
XI”, Proceedings of ASME 2014 PV&P Conference, paper
PVP#2014-28222
6. Anderson, T.L. et al, “Development of stress intensity factor
solutions for surface and embedded cracks in API 579”,
WRC Bulletin 471, Welding Research Council, Inc. New
York, NY, May 2002
7. Cipolla, R.C., and Lee, D.R., 2004, “Technical basis for
equations for stress intensity factor coefficients in ASME
Section XI Appendix A”, Proceedings of ASME 2004
PV&P Conference, Vol. 480, paper PVP#2004-2708, pp.
301-312
8. Chapuliot, S., 2000, “Formulaire de KI pour les tubes
comportant un défaut de surface semi-elliptique longitudinal
ou circonférentiel, interne ou externe”, Commissariat à
l’Energie Atomique, Internal Report CEA-R-5900
9. Miura, N. et al, 2008, “Comparison of stress intensity factor
solutions for axial and circumferential cracks”, Nuclear
Engineering and Design, Vol. 238, pp. 423-434
10. Pook, L.P., “Crack profiles and corner point singularities”,
2000, Fatigue and Fracture of Engineering Materials and
Structures, Vol. 23, pp. 141-150
API a/c = 1
RSE-M a/c = 1
1,8
API a/c = 0.5
RSE-M a/c = 0.5
API a/c = 0.25
1,6 RSE-M a/c = 0.25
G0 - deepest point
API a/c = 0
RSE-M a/c = 0
1,4
1,2
1,0
0,8
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
a/t
FIG. 1 AN AXIAL ID SEMI-ELLIPTICAL SURFACE FLAW IN A FIG. 4 CYLINDER WITH Ri/t = 1. COMPARISON OF THE
CYLINDER COEFFICIENT G0 AT THE DEEPEST POINT
G0 - deepest point
0,90 API a/c = 0
API a/c = 0.0625
RSE-M a/c = 0
RSE-M a/c = 0.0625
0,85 2,5
API a/c = 0
2,0
0,75
0,70
1,5
0,65
0,60 1,0
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
a/t a/t
FIG. 5 CYLINDER WITH Ri/t = 1. COMPARISON OF THE FIG. 8 CYLINDER WITH Ri/t = 5. COMPARISON OF THE
COEFFICIENT G1 AT THE DEEPEST POINT COEFFICIENT G0 AT THE DEEPEST POINT
Axial ID surface flaws in cylinders - Ri/t = 1 - G0 - surface point Axial ID surface flaws in cylinders - Ri/t = 5 - G1 - deepest point
1,4 1,8
API a/c = 1
1,2 RSE-M a/c = 1
1,6
API a/c = 0.5
RSE-M a/c = 0.5
1,0 API a/c = 0.25
1,4
RSE-M a/c = 0.25
G1 - deepest point
G0 - surface point
0,0 0,6
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
a/t a/t
FIG. 6 CYLINDER WITH Ri/t = 1. COMPARISON OF THE FIG. 9 CYLINDER WITH Ri/t = 5. COMPARISON OF THE
COEFFICIENT G0 AT THE SURFACE POINT COEFFICIENT G1 AT THE DEEPEST POINT
Axial ID surface flaws in cylinders - Ri/t = 1 - G1 - surface point Axial ID surface flaws in cylinders - Ri/t = 5 - G0 - surface point
0,30 1,6
1,4
0,25
1,2
0,20
1,0
G1 - surface point
G0 - surface point
API a/c = 1
API a/c = 1
RSE-M a/c = 1
0,15 0,8 RSE-M a/c = 1
API a/c = 0.5
API a/c = 0.5
RSE-M a/c = 0.5 0,6
RSE-M a/c = 0.5
0,10
API a/c = 0.25
API a/c = 0.25
RSE-M a/c = 0.25 0,4
RSE-M a/c = 0.25
0,05 API a/c = 0.0625
API a/c = 0.0625
0,2
RSE-M a/c = 0.0625
RSE-M a/c = 0.0625
0,00 0,0
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
a/t a/t
FIG. 7 CYLINDER WITH Ri/t = 1. COMPARISON OF THE FIG. 10 CYLINDER WITH Ri/t = 5. COMPARISON OF THE
COEFFICIENT G1 AT THE SURFACE POINT COEFFICIENT G0 AT THE SURFACE POINT
1,4
0,25
1,2
0,20
1,0 API a/c = 1
G1 - surface point
G0 - surface point
API a/c = 1
RSE-M a/c = 1
RSE-M a/c = 1
0,15 0,8 API a/c = 0.5
API a/c = 0.5
RSE-M a/c = 0.5
RSE-M a/c = 0.5 0,6
API a/c = 0.25
0,10
API a/c = 0.25
RSE-M a/c = 0.25
RSE-M a/c = 0.25 0,4
API a/c = 0.0625
0,05 API a/c = 0.0625
0,2 RSE-M a/c = 0.0625
RSE-M a/c = 0.0625
0,00 0,0
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
a/t a/t
FIG. 11 CYLINDER WITH Ri/t = 5. COMPARISON OF THE FIG. 14 CYLINDER WITH Ri/t = 10. COMPARISON OF THE
COEFFICIENT G1 AT THE SURFACE POINT COEFFICIENT G0 AT THE SURFACE POINT
Axial ID surface flaws in cylinders - Ri/t = 10 - G0 - deepest point Axial ID surface flaws in cylinders - Ri/t = 10 - G1 - surface point
6,0 0,30
G1 - surface point
4,0 API a/c = 1
API a/c = 0.0625
RSE-M a/c = 0.0625 RSE-M a/c = 1
3,5 0,15
API a/c = 0 API a/c = 0.5
3,0 RSE-M a/c = 0 RSE-M a/c = 0.5
0,10 API a/c = 0.25
2,5
RSE-M a/c = 0.25
2,0 API a/c = 0.0625
0,05
RSE-M a/c = 0.0625
1,5
1,0 0,00
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
a/t a/t
FIG. 12 CYLINDER WITH Ri/t = 10. COMPARISON OF THE FIG. 15 CYLINDER WITH Ri/t = 10. COMPARISON OF THE
COEFFICIENT G0 AT THE DEEPEST POINT COEFFICIENT G1 AT THE SURFACE POINT
Axial ID surface flaws in cylinders - Ri/t = 10 - G1 - deepest point Axial ID surface flaws in cylinders - Ri/t = 20 - G0 - deepest point
2,4 8
1,0
2
0,8
0,6 1
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
a/t a/t
G0 - deepest point
RSE-M a/c = 0.25
2,1
API a/c = 0.0625
API a/c = 0.0625
RSE-M a/c = 0.0625
RSE-M a/c = 0.0625 7
API a/c = 0
API a/c = 0
RSE-M a/c = 0
1,6 RSE-M a/c = 0
5
1,1
3
0,6 1
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
a/t a/t
G1 - deepest point
RSE-M a/c = 0.25
1,0
G0 - surface point
0,6
0,0 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
a/t
a/t
1,4
0,25
1,2
0,20 1,0
G0 - surface point
G1 - surface point
API a/c = 1
API a/c = 1
RSE-M a/c = 1
RSE-M a/c = 1 0,8
0,15 API a/c = 0.5
API a/c = 0.5
0,6 RSE-M a/c = 0.5
RSE-M a/c = 0.5
0,10 API a/c = 0.25
API a/c = 0.25
0,4 RSE-M a/c = 0.25
RSE-M a/c = 0.25
API a/c = 0.0625
API a/c = 0.0625
0,05 0,2
RSE-M a/c = 0.0625 RSE-M a/c = 0.0625
0,0
0,00 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
a/t
a/t
FIG. 19 CYLINDER WITH Ri/t = 20. COMPARISON OF THE FIG. 22 CYLINDER WITH Ri/t = (PLATE). COMPARISON OF
COEFFICIENT G1 AT THE SURFACE POINT THE COEFFICIENT G0 AT THE SURFACE POINT
0,30
1,4
0,25 1,2 FE
API-579
RSE-M
1,0
0,20
G0 coefficient
G1 - surface point
0,8
API a/c = 1
0,15
RSE-M a/c = 1 0,6
API a/c = 0.5
0,8
0,7
0,6
G1 Coefficient
0,5
EF
0,4 API-579
RSE-M
0,3
0,2
0,1
0,0
0 10 20 30 40 50 60 70 80 90
Elliptic angle - surface point towards deepest point (°)