Universal Weight Function Consistent Method to Fit Polynomial Stress Distribution for Calculation of Stress Intensity Factor

Author(s):  
Douglas A. Scarth ◽  
Steven X. Xu

Procedures for analytical evaluation of flaws in nuclear pressure boundary components are provided in Section XI of the ASME B&PV Code. The flaw evaluation procedure requires calculation of the stress intensity factor. Engineering procedures to calculate the stress intensity factor are typically based on a polynomial equation to represent the stress distribution through the wall thickness, where the polynomial equation is fitted using the least squares method to discrete data point of stress through the wall thickness. However, the resultant polynomial equation is not always an optimum fit to stress distributions with large gradients or discontinuities. Application of the weight function method enables a more accurate representation of the stress distribution for the calculation of the stress intensity factor. Since engineering procedures and engineering software for flaw evaluation are typically based on the polynomial equation to represent the stress distribution, it would be desirable to incorporate the advantages of the weight function method while still retaining the framework of the polynomial equation to represent the stress distribution when calculating the stress intensity factor. A method to calculate the stress intensity factor using a polynomial equation to represent the stress distribution through the wall thickness, but which provides the same value of the stress intensity factor as is obtained using the Universal Weight Function Method, is provided in this paper.

2012 ◽  
Vol 134 (6) ◽  
Author(s):  
Douglas A. Scarth ◽  
Steven X. Xu

Procedures for the analytical evaluation of flaws in nuclear pressure boundary components are provided in Section XI of the ASME B&PV Code. The flaw evaluation procedure requires calculation of the stress intensity factor. Engineering procedures to calculate the stress intensity factor are typically based on a polynomial equation to represent the stress distribution through the wall thickness, where the polynomial equation is fitted using the least squares method to discrete data point of stress through the wall thickness. However, the resultant polynomial equation is not always an optimum fit to stress distributions with large gradients or discontinuities. Application of the weight function method enables a more accurate representation of the stress distribution for the calculation of the stress intensity factor. Since engineering procedures and engineering software for flaw evaluation are typically based on the polynomial equation to represent the stress distribution, it would be desirable to incorporate the advantages of the weight function method while still retaining the framework of the polynomial equation to represent the stress distribution when calculating the stress intensity factor. A method to calculate the stress intensity factor using a polynomial equation to represent the stress distribution through the wall thickness, but which provides the same value of the stress intensity factor as is obtained using the Universal Weight Function Method, is provided in this paper.


Author(s):  
Yinsheng Li ◽  
Hiroto Itoh ◽  
Kunio Hasegawa ◽  
Steven X. Xu ◽  
Douglas A. Scarth

Many solutions of the stress intensity factor have been proposed in recent years. However, most of them take only third or fourth-order polynomial stress distributions into account. For complicated stress distributions which are difficult to be represented as third or fourth-order polynomial equations over the stress distribution area such as residual stress distributions or thermal stress distributions in dissimilar materials, it is important to further improve the accuracy of the stress intensity factor. In this study, a weight function method with segment-wise polynomial interpolation is proposed to calculate solutions of the stress intensity factor for complicated stress distributions. By using this method, solutions of the stress intensity factor can be obtained without employing finite element analysis or difficult calculations. It is therefore easy to use in engineering applications. In this method, the stress distribution area is firstly divided into several segments and the stress distribution in each segment is curve fitted to segment-wise polynomial equation. The stress intensity factor is then calculated based on the weight function method and the fitted stress distribution in each segment. Some example solutions for both infinite length cracks and semi-elliptical cracks are compared with the results from finite element analysis. In conclusion, it is confirmed that this method is applicable with high accuracy to the calculation of the stress intensity factor taking actual complicated stress distributions into consideration.


2014 ◽  
Vol 136 (2) ◽  
Author(s):  
Yinsheng Li ◽  
Kunio Hasegawa ◽  
Steven X. Xu ◽  
Douglas A. Scarth

Many solutions of the stress intensity factor have been proposed in recent years. However, most of them take only third or fourth-order polynomial stress distributions into account. For complicated stress distributions which are difficult to be represented as third or fourth-order polynomial equations over the stress distribution area such as residual stress distributions or thermal stress distributions in dissimilar materials, it is important to further improve the accuracy of the stress intensity factor. In this study, a weight function method with segment-wise polynomial interpolation is proposed to calculate solutions of the stress intensity factor for complicated stress distributions. By using this method, solutions of the stress intensity factor can be obtained without employing finite element analysis or difficult calculations. It is therefore easy to use in engineering applications. In this method, the stress distribution area is firstly divided into several segments and the stress distribution in each segment is curve fitted to segment-wise polynomial equation. The stress intensity factor is then calculated based on the weight function method and the fitted stress distribution in each segment. Some example solutions for both infinite length cracks and semi-elliptical cracks are compared with the results from finite element analysis. In conclusion, it is confirmed that this method is applicable with high accuracy to the calculation of the stress intensity factor taking actual complicated stress distributions into consideration.


Author(s):  
Steven X. Xu ◽  
Douglas A. Scarth ◽  
Russell C. Cipolla

Linear elastic fracture mechanics based flaw evaluation procedures in ASME Section XI require calculation of the stress intensity factor. The method to calculate the stress intensity factor that is provided in the 2010 Edition of Appendix A of ASME Section XI is to fit the stress distribution ahead of the crack tip to a polynomial equation, and then use standardized influence coefficients. In PVP2011-57911, an alternate method for calculation of the stress intensity factor for a surface flaw that makes explicit use of the Universal Weight Function Method and does not require a polynomial fit to the actual stress distribution was proposed for implementation into Appendix A of Section XI. The alternate method provides closed-form solutions for the stress intensity factor. A numerical approximation is the assumed piece-wise linear variation of stress between discrete locations where stresses are known. For highly nonlinear stress distributions, piece-wise cubic interpolation of stress over intervals between discrete locations where stresses are known is an improvement over piece-wise linear interpolation. Investigation of a cubic interpolation of stress between discrete locations where stresses are known has been conducted. Closed-form equations for calculation of the stress intensity factor for a surface flaw were developed using the Universal Weight Function Method and generic piece-wise cubic interpolation of stress over intervals. Example calculations are provided to compare the results of stress intensity factors using piece-wise cubic interpolation with piece-wise linear interpolation.


Author(s):  
Rui Sun ◽  
Zongwen An ◽  
Hong-Zhong Huang ◽  
Qiming Ma

Propagation of a critical unstable crack under the action of static or varying stresses is determined by the intensity of strain field at tips of the crack. Stress intensity factor (SIF) is an important parameter in fracture mechanics, which is used as a criterion to judge the unstable propagation of a crack and plays an important role in calculating crack propagation life. SIF is related to both geometrical form and loading condition of a structure. In the paper, a weight function method is introduced to study crack propagation of center through cracks and edge cracks in a finite-size plate. In addition, finite element method, linear regression, and polynomial interpolating technique are used to simulate and verify the proposed method. Comparison studies among the proposed and current methods are performed as well. The results show that the weight function method can be used to calculate SIF easily.


Author(s):  
Douglas A. Scarth ◽  
Steven X. Xu

A method for calculating the stress intensity factor for linear elastic fracture mechanics based flaw evaluation is provided in Appendix A-3000 of ASME Section XI. In the 2010 Edition of ASME Section XI, the calculation of stress intensity factor for a surface crack is based on characterization of stress field with a cubic equation and use of influence coefficients. The influence coefficients are currently only provided for flat plate geometry in tabular format. The ASME Section XI Working Group on Flaw Evaluation is in the process of rewriting Appendix A-3000. Proposed major updates include the implementation of explicit use of Universal Weight Function Method for calculation of the stress intensity factor for a surface flaw and the inclusion of closed-form influence coefficients for cylinder geometry. The explicit use of weight function method eliminates the need for fitting polynomial equations to the actual through-thickness stress distributions at crack location. In this paper, the proposed Appendix A procedure is applied to calculate the stress intensity factors in closed-form for an axial ID surface flaw in a cylinder subjected to a set of nonlinear hoop weld residual stress profiles. The calculated stress intensity factor results are compared with the results calculated based on the current method in Appendix A using cubic equations to represent the stress distribution. Three-dimensional finite element analyses were performed to verify the accuracy of the stress intensity factor results calculated based on the current and proposed Appendix A procedures. The results in this paper support the implementation of the proposed stress intensity factor calculation procedure into ASME Code.


Author(s):  
Steven X. Xu ◽  
Douglas A. Scarth ◽  
Russell C. Cipolla

Analytical evaluation procedures for determining the acceptability of flaws detected during in-service inspection of nuclear power plant components are provided in Section XI of the ASME Boiler and Pressure Vessel Code. Linear elastic fracture mechanics based evaluation procedures in ASME Section XI require calculation of the stress intensity factor. A method for calculating the stress intensity factor is provided in Appendix A of ASME Section XI. This method consists of a two-step process. In the first step, the stress distribution, as calculated in the absence of the flaw, is obtained at the flaw location. For a surface flaw, the stress distribution at the flaw location is then fitted to a third-order polynomial equation. In the second step, the fitted polynomial representation of the stress distribution is used with standardized influence coefficients to calculate the stress intensity factor. An alternate method for calculation of the stress intensity factor for a surface flaw that makes explicit use of the Universal Weight Function Method and does not require a polynomial fit to the actual stress distribution is proposed in this paper for implementation into Appendix A of ASME Section XI. Universal Weight Function coefficients are determined from standardized influence coefficients through closed-form equations. Closed-form equations for calculation of the stress intensity factor are provided. The technical basis and verification for this alternate method for calculation of the stress intensity factor are described in this paper.


Sign in / Sign up

Export Citation Format

Share Document