Effective Methods to Determine Stress Intensity Factors for 2D and 3D Cracks

Author(s):  
Xian-Kui Zhu ◽  
Brian N. Leis

Increasing concern for crack assessment in the pipeline industry motivates analysis to quantify the crack driving force, with the linear-elastic fracture mechanics stress intensity factor, denoted K, viable for many vintage pipeline applications. This paper presents a brief review of numerical methods developed for calculating K via the finite element analysis (FEA) as a background to identify the “best” approaches for such purposes. The existing methods can be categorized into three groups: the displacement-based methods, the stress-based methods, and the energy-based methods. The first group involves the displacement extrapolation method, the quarter-point displacement method, and the displacement correction method. The second group involves the stress extrapolation method and the force method. The third group includes the J-integral method, the stiffness derivative method, the virtual crack extension method, the virtual crack closure technique (VCCT) and ABAQUS direct K output method. Based on the review, four methods were selected and evaluated for a central-cracked plate (CCP) specimen based on the FEA calculations via ABAQUS. The “best” methods are then applied in an analysis of K for through-wall cracks in a line pipe — important reference geometry for leak-versus-rupture analysis.

2020 ◽  
Vol 36 (2) ◽  
pp. 235-243
Author(s):  
N. Zhu ◽  
E. Oterkus

ABSTRACTThis paper introduces a new approach to calculate stress intensity factors based on a combination of Displacement Extrapolation Method and Peridynamic Theory. After obtaining the displacement field from Peridynamic Theory, by appropriately selecting nodes at the crack tip region and their displacements yield stress intensity factors at the crack tips. To demonstrate the capability of the proposed approach, three different benchmark problems are considered including plate with a central crack, plate with an edge crack and plate with a slanted crack. Results evaluated from the current approach are compared against analytical and finite element analysis results, and good agreement is obtained between three different approaches. This shows that coupled Displacement Extrapolation Method and Peridynamic Theory approach can be an alternative method to calculate stress intensity factors.


Author(s):  
Xian-Kui Zhu

Crack assessments for pressure vessels often need to quantify the crack driving force — stress intensity factor K with the linear-elastic fracture mechanics methods. Different numerical methods have been developed to calculate the stress intensity factors for complex cracks. Of which, four typical methods, i.e., the displacement extrapolation method, the virtual crack closure technique (VCCT), the J-integral conversion method, and the direct K output method are selected and evaluated in this paper using the finite element analysis (FEA) and ABAQUS software. The evaluations are performed based on the benchmark FEA calculations in the linear-elastic conditions for the central-cracked panel (CCP) specimen in the two-dimensional (2D) plane strain conditions. The “best method” is then determined and used to calculate the stress intensity factor for the CCP specimen with a through-thickness crack in the three-dimensional (3D) conditions. The results show that ABAQUS can simply determine very accurate K values for both 2D and 3D cracks.


Author(s):  
Kisaburo Azuma ◽  
Yinsheng Li ◽  
Kunio Hasegawa

The interaction of multiple flaws in close proximity to one another may increase the stress intensity factor of the flaw in structures and components. This interaction effect is not distributed uniformly along the crack front. For instance, the strongest interaction is generally observed at the point closest to a neighboring flaw. For this reason, the closest point could show a higher value of the stress intensity factor than all other points in some cases, even if the original value at the point of the single flaw is relatively low. To clarify the condition when the closest point shows the maximum stress intensity factor, we investigated the interaction of two similar elliptical flaws in an infinite model subjected to remote tension loading. The stress intensity factor of the elliptical flaws was obtained by performing finite element analysis of a linear elastic solid. The results indicated that the interaction factors along the crack front can be expressed by a simple empirical formula. Finally, we show the relationship between geometrical features of the flaw and the stress intensity factor at the closest point to a neighboring flaw.


2010 ◽  
Vol 303-304 ◽  
pp. 63-83
Author(s):  
Ehsan Mahdavi ◽  
Mahmoud Mosavi Mashhadi ◽  
M. Amidpour

It is well known that the crack growth rate fatigue and stress corrosion cracking can be approximated by a power function of the stress intensity factor. In this study, stress intensity factor for elliptical crack under the uniform tension in linear elastic fracture mechanics (LEFM) is investigated therefore for this purpose, a pressure vessel modeled by finite element. A crack modeled on the pressure vessel and then the stress intensity factor for crack propagation in different methods is evaluated. Finite element analysis calculates stress intensity factor in the values of the J-integral are based on the stress intensity factors, JK, and by evaluating the contour integral directly, JA. The stability of crack growth is considered so the ductile crack extension is determined by pursuing the equilibrium between loading and crack resistance. Using especial method of meshing caused to have accurate results. This method causes to decrease run time and considerable accuracy. Then stress intensity factor is calculated for different position of the crack such as crack front and then compared to each other.


2013 ◽  
Vol 3 (5) ◽  
pp. 506-510
Author(s):  
A. Boulenouar ◽  
N. Benseddiq ◽  
M. Mazari

When the loading or the geometry of a structure is not symmetrical about the crack axis, rupture occurs in mixed mode loading and the crack does not propagate in a straight line. It is then necessary to use kinking criteria to determine the new direction of crack propagation. The aim of this work is to present a numerical modeling of crack propagation under mixed mode loading conditions. This work is based on the implementation of the displacement extrapolation method in a FE code and the strain energy density theory in a finite element code. At each crack increment length, the kinking angle is evaluated as a function of stress intensity factors. In this paper, we analyzed the mechanical behavior of inclined cracks by evaluating the stress intensity factors. Then, we presented the examples of crack propagation in structures containing inclusions and cavities.


Sign in / Sign up

Export Citation Format

Share Document