Stress Intensity Factors for Partially Autofrettaged Pressurized Thick-Walled Cylinders Containing Closely and Densely Packed Cracks

2010 ◽  
Vol 132 (5) ◽  
Author(s):  
Q. Ma ◽  
C. Levy ◽  
M. Perl

Due to acute temperature gradients and repetitive high-pressure impulses, extremely dense internal surface cracks can be practically developed in highly pressurized thick-walled vessels, typically in gun barrels. In the authors’ previous studies, networks of typical radial and longitudinal-coplanar, semi-elliptical, internal surface cracks have been investigated assuming both ideal and realistic full autofrettage residual stress fields (ε=100%). The aim of the present work is to extend the analysis twofold: to include various levels of partially autofrettaged cylinders and to consider configurations of closely and densely packed radial crack arrays. To accurately assess the stress intensity factors (SIFs), significant computational efforts and strategies are necessary, especially for networks with closely and densely packed cracks. This study focuses on the determination of the distributions along the crack fronts of KIP, the stress intensity factor due to internal pressure KIA, the negative stress intensity factor resulting from the residual stress field due to ideal or realistic autofrettage, and KIN, the combined SIF KIN=KIP−|KIA|. The analysis is performed for over 1000 configurations of closely and densely packed semicircular and semi-elliptical networked cracks affected by pressure and partial-to-full autofrettage levels of ε=30–100%, which is of practical benefit in autofrettaged thick-walled pressure vessels. The 3-D analysis is performed via the finite element method and the submodeling technique employing singular elements along the crack front and the various symmetries of the problem. The network cracks will include up to 128 equally spaced cracks in the radial direction: with relative longitudinal crack spacing, 2c/d, from 0.1 to 0.99; autofrettage level of 30–100%; crack depth to wall thickness ratios, a/t, from 0.01 to 0.4; and, cracks with various ellipticities of crack depth to semicrack length, a/c, from 0.2 to 2. The results clearly indicate that the combined SIFs are considerably influenced by the three dimensionality of the problem and the Bauschinger effect (BE). The Bauschinger effect is found to have a dramatic effect on the prevailing combined stress intensity factors, resulting in a considerable reduction of the fatigue life of the pressure vessel. While the fatigue life can be finite for ideal autofrettage, it is normally finite for realistic autofrettage for the same crack network. Furthermore, it has been found that there are differences in the character of the SIFs between closely packed and densely packed crack networks, namely, more dramatic drop-offs in KIA and KIN at the crack-inner bore interface for densely packed cracks further influenced by crack depth.

Author(s):  
Do-Jun Shim ◽  
Matthew Kerr ◽  
Steven Xu

Recent studies have shown that the crack growth of PWSCC is mainly driven by the weld residual stress (WRS) within the dissimilar metal weld. The existing stress intensity factor (K) solutions for surface cracks in pipe typically require a 4th order polynomial stress distribution through the pipe wall thickness. However, it is not always possible to accurately represent the through thickness WRS with a 4th order polynomial fit and it is necessary to investigate the effect of the WRS fitting on the calculated stress intensity factors. In this paper, two different methods were used to calculate the stress intensity factor for a semi-elliptical circumferential surface crack in a pipe under a given set of simulated WRS. The first method is the Universal Weight Function Method (UWFM) where the through thickness WRS distribution can be represented as a piece-wise cubic fit. In the second method, the through thickness WRS profiles are represented as a 4th order polynomial curve fit (both using the entire wall thickness data and only using data up to the crack-tip). In addition, three-dimensional finite element (FE) analyses (using the simulated weld residual stress) were conducted to serve as a reference solution. The results of this study demonstrate the potential sensitivity of stress intensity factors to 4th order polynomial fitting artifacts. The piece-wise WRS representations used in the UWFM was not sensitive to these fitting artifacts and the UWFM solutions were in good agreement with the FE results.


1980 ◽  
Vol 102 (4) ◽  
pp. 342-346 ◽  
Author(s):  
J. C. Newman ◽  
I. S. Raju

The purpose of this paper is to present stress-intensity factors for a wide range of semi-elliptical surface cracks on the inside of pressurized cylinders. The ratio of crack depth to crack length ranged from 0.2 to 1; the ratio of crack depth to wall thickness ranged from 0.2 to 0.8; and the ratio of wall thickness to vessel radius was 0.1 to 0.25. The stress-intensity factors were calculated by a three-dimensional finite-element method. The finite-element models employ singularity elements along the crack front and linear-strain elements elsewhere. The models had about 6500 degrees of freedom. The stress-intensity factors were evaluated from a nodal-force method. An equation for the stress-intensity factors was obtained from the results of the present analysis. The equation applies over a wide range of configuration parameters and was within about 5 percent of the present results. A comparison was also made between the present results and other analyses of internal surface cracks in cylinders. The results from a boundary-integral equation method were in good agreement (± 2 percent) and those from another finite-element method were in fair agreement (± 8 percent) with the present results.


1997 ◽  
Vol 32 (5) ◽  
pp. 351-363 ◽  
Author(s):  
X B Lin ◽  
R A Smith

Stress intensity factors for internal semi-elliptical surface cracks in autofrettaged cylinders with and without internal pressures applied are presented. The three-dimensional finite element based displacement method with the crack tip square-root singularity of stresses and strains simulated is used to evaluate the stress intensity factors along the crack front. Both allowing and disallowing the overlapping of crack faces are considered in this investigation, the latter being simulated by considering crack surface contact through a kind of interface element introduced into the cracked area. The residual stress distribution assumed to act on the crack face is obtained according to Tresca's yield criterion with the material assumed to be elastic-perfectly plastic. Three different overstrains of autofrettage are chosen. The results show that the stress intensity factor is generally underestimated if the crack contact that has actually occurred is ignored, which may lead to a danger in the assessment of either fracture strength or fatigue life. Implications of the stress intensity factor results are also briefly discussed, particularly for the prediction of fatigue lives, and it is shown that the full autofrettage treatment might be the most beneficial for increasing the fatigue life of cracks initiated from the inner core.


1982 ◽  
Vol 104 (4) ◽  
pp. 293-298 ◽  
Author(s):  
I. S. Raju ◽  
J. C. Newman

The purpose of this paper is to present stress-intensity factor influence coefficients for a wide range of semi-elliptical surface cracks on the inside or outside of a cylinder. The crack surfaces were subjected to four stress distributions: uniform, linear, quadratic, and cubic. These four solutions can be superimposed to obtain stress-intensity factor solutions for other stress distributions, such as those caused by internal pressure and by thermal shock. The results for internal pressure are given herein. The ratio of crack depth to crack length from 0.2 to 1; the ratio of crack depth to wall thickness ranged from 0.2 to 0.8; and the ratio of wall thickness to vessel radius was 0.1 or 0.25. The stress-intensity factors were calculated by a three-dimensional finite-element method. The finite-element models employ singularity elements along the crack front and linear-strain elements elsewhere. The models had about 6500 degrees of freedom. The stress-intensity factors were evaluated from a nodal-force method. The present results were also compared to other analyses of surface cracks in cylinders. The results from a boundary-integral equation method agreed well (±2 percent), and those from other finite-element methods agreed fairly well (±10 percent) with the present results.


Author(s):  
Hsin Jen Hoh ◽  
John H. L. Pang ◽  
Kin Shun Tsang

Offshore pipelines and risers transfer oil and gas across long distances, from seabed to production facility to the surface. The long pipelines are formed by welding together pipe segments. The welded joints formed are a source of stress concentration and defects from which fatigue cracks can grow. Hence, it is imperative that the effect of the weld geometry on the stress concentration be understood so that appropriate measures can be taken to assess the potential remaining service life of the welded structure. The effects can be understood by the linear elastic fracture mechanics approach, where the stress intensity factors quantify the stress concentration. While the classical equations of Newman and Raju have been long available for semi-elliptical surface cracks in plates, no similarly elegant stress intensity factor solutions are available for pipes. There have been solutions in tabular form which can be cumbersome in practice. Moreover, solutions of welded pipe geometries have not been developed. The objectives of the current work are to develop closed-form solutions for stress intensity factors for external semi-elliptical surface cracks in plates. The welded pipe geometry will also be studied to develop solutions for the weld toe magnification factors of welded pipe geometries. The stress intensity factors can be used to determine the propagation rate of cracks in pipe or welded pipe geometries. The stress intensity factors are obtained by the J-integral output of the three-dimensional finite element method. First, a plate with a circular crack is modelled. The initial step transforms the model to a plate with a semi-elliptical crack with the appropriate crack aspect ratio and width. A second transformation follows to transform the geometry to pipe form. The main parameters studied are the relative crack depth to thickness, crack aspect ratio, radius and thickness. The developed stress intensity factor solutions can be reduced to the classical equations. The new solutions show good agreement compared to previous work. A similar approach is developed to study the welded pipe geometry to develop weld toe magnification factor solutions. The weld toe magnification factor solutions for certain geometries are presented as a function of the relative crack depth. The stress intensity factor solutions are then applied to predict the crack growth rates of cracks in pipe geometries. The prediction was conducted by a program written to assess the fatigue life of single and multiple cracks in pipes and welded pipes. The fatigue life assessment of welded pipes using the weld toe magnification factor solutions shows how significantly the weld geometry affects fatigue life.


Author(s):  
D. J. Shim ◽  
S. Tang ◽  
T. J. Kim ◽  
N. S. Huh

Stress intensity factor solutions are readily available for flaws found in pipe to pipe welds or shell to shell welds (i.e., circumferential/axial crack in cylinder). In some situations, flaws can be detected in locations where an appropriate crack model is not readily available. For instance, there are no practical stress intensity factor solutions for circular-arc cracks which can form in circular welds (e.g., nozzle to vessel shell welds and storage cask closure welds). In this paper, stress intensity factors for circular-arc cracks in finite plates were calculated using finite element analysis. As a first step, stress intensity factors for circular-arc through-wall crack under uniform tension and crack face pressure were calculated. These results were compared with the analytical solutions which showed reasonable agreement. Then, stress intensity factors were calculated for circular-arc semi-elliptical surface cracks under the lateral and crack face pressure loading conditions. Lastly, to investigate the applicability of straight crack solutions for circular-arc cracks, stress intensity factors for circular-arc and straight cracks (both through-wall and surface cracks) were compared.


Author(s):  
Mayumi Ochi ◽  
Kiminobu Hojo ◽  
Itaru Muroya ◽  
Kazuo Ogawa

Alloy 600 weld joints have potential for primary water stress corrosion cracks (PWSCC). At the present time it has been understood that PWSCC generates and propagates in the Alloy 600 base metal and the Alloy 600 weld metal and there has been no observation of cracking the stainless and the low alloy steel. For the life time evaluation of the pipes or components the crack extension analysis is required. To perform the axial crack extension analysis the stress intensity database or estimation equation corresponding to the extension crack shape is needed. From the PWSCC extension nature mentioned above, stress intensity factors of the conventional handbooks are not suitable because most of them assume a semi-elliptical crack and the maximum aspect ratio crack depth/crack half length is one (The evaluation in this paper had been performed before API 579-1/ASME FFS was published). Normally, with the advance of crack extension in the thickness direction at the weld joint, the crack aspect ratio exceeds one and the K-value of the conventional handbook can not be applied. Even if those equations are applied, the result would be overestimated. In this paper, considering characteristics of PWSCC’s extension behavior in the welding material, the axial crack was modeled in the FE model as a rectangular shape and the stress intensity factors at the deepest point were calculated with change of crack depth. From the database of the stress intensity factors, the simplified equation of stress intensity factor with parameter of radius/thickness and thickness/weld width was proposed.


Author(s):  
Christian Malekian ◽  
Eric Wyart ◽  
Michael Savelsberg ◽  
Anne Teughels ◽  
Pierre-Eric Fouquet ◽  
...  

Most of the literature about fracture mechanics considers cracks having an elliptical shape with a flaw aspect ratio a/l lower or equal to 0.5 where ‘a’ is the crack depth and ‘l’ the total length of the crack. This is also case in the ASME XI Appendix A where Stress Intensity Factors KI formulations are given for a large range of crack depths and for a flaw aspect ratio a/l between 0 and 0.5. The limitation to 0.5 corresponds to a semi-circular shape for surface cracks and to a circular shape for subsurface cracks. This limitation does not seem to be inspired by a theoretical limitation nor by a computational limit. Moreover, it appears that limiting the ratio a/l to 0.5 may generate in some cases some unnecessary conservatism in flaw analysis. The present article specifically deals with the more unusual narrow cracks having a/l >0.5, in the case of surface cracks in infinite flat plates. Several Finite-Elements calculations are performed to compute KI for a large range of crack depths and for 4 typical load cases (uniform, linear, quadratic and cubic). The results can be presented with the same formalism as in the ASME XI Appendix A, such that the work can provide an extension of the ASME coefficients in table A-3320-1&2. By doing the study, one had the opportunity to compare the results obtained by two different Finite-Elements softwares (Systus and Ansys), each one with a different cracked mesh. In addition, a comparison has been made for some cases with results obtained by a XFEM approach (eXtended Finite-Element Method), where the crack does not need to be meshed in the same way as in classical Finite-Elements. The results indicate how the KI can be reduced when considering the real flaw aspect ratio instead of the conventional semi-circular flaw shape. They also show that, for specific theoretical stress distributions, it is not always possible to reduce the analysis of KI to only 2 points, namely the crack surface point and the crack deepest point. The crack growth evaluation of such unusual crack shape should still be investigated to verify whether simple rules can be established to estimate the evolution of the crack front.


2011 ◽  
Vol 214 ◽  
pp. 192-196 ◽  
Author(s):  
Al Emran Ismail ◽  
Ahmad Kamal Ariffin ◽  
Shahrum Abdullah ◽  
Mariyam Jameelah Ghazali ◽  
Ruslizam Daud

This study presents a numerical investigation on the stress intensity factors (SIF), K of surface cracks in round bars that were obtained under pure torsion loadings or mode III. ANSYS finite element analysis (FEA) was used to determine the SIFs along the crack front of surface cracks embedded in the solid circular bars. 20-node isoparametric singular elements were used around the crack tip by shifting the mid-side node ¼-position close to a crack tip. Different crack aspect ratio, a/b were used ranging between 0.0 to 1.2 and relative crack depth, a/D were ranged between 0.1 to 0.6. Mode I SIF, KI obtained under bending moment was used to validate the proposed model and it was assumed this proposed model validated for analyzing mode III problems. It was found that, the mode II SIF, FII and mode III SIF, FIII were dependent on the crack geometries and the sites of crack growth were also dependent on a/b and a/D.


Sign in / Sign up

Export Citation Format

Share Document