Elastic T-stress solutions for semi-elliptical surface cracks in finite thickness plates

Three-dimensional finite element analyses have been conducted to calculate the elastic T-stress for semi-elliptical surface cracks in finite thickness plates. Far-field tension and bending loads were considered. The analysis procedures and results were verified using both exact solutions and approximate solutions. The T-stress solutions are presented along the crack front for cracks with a/t values of 0.2, 0.4, 0.6 or 0.8 and a/c values of 0.2, 0.4, 0.6 or 1.0. Based on the present finite element calculations for T-stress, empirical equations for the T-stress at three locations: the deepest, the surface and the middle points of the crack front under tension or bending are presented. The numerical results are approximated by empirical formulae fitted with an accuracy of 1% or better. They are valid for 0.2?a/c?1 and 0?a/t?0.8. These T-stress results together with the corresponding K or J values for surface cracks are suitable for the analysis of constraint effects for surface cracked components.     

The T-stress solutions for through-wall circumferential cracks in cylinders subjected to general loading conditions

The elastic T-stress is a parameter used to define the level of constraint at a crack tip. It is important to provide T-stress solutions for practical geometries to apply the constraint-based fracture mechanics methodology. In the present work, T-stress solutions are provided for circumferential through-wall cracks in thin-walled cylinders. First, cylinders with a circumferential through-wall crack were analyzed using the finite element method. Three cylinder geometries were considered; defined by the mean radius of the cylinder (R) to wall thickness (t) ratios: R/t = 5, 10, and 20. The T-stress was obtained at eight crack lengths (θ/π = 0.0625, 0.1250, 0.1875, 0.2500, 0.3125, 0.3750, 0.4375, and 0.5000, θ is the crack half angle). Both crack face loading and remote loading conditions were considered including constant, linear, parabolic and cubic crack face pressures and remote tension and bending. The results for constant and linear crack face pressure were used to derive weight functions for T-stress for the corresponding cracked geometries. The weight functions were validated against several linear and non-linear stress distributions. The derived weight functions are suitable for T-stress calculations for circumferential cracks in cylinders under complex stress fields.     

Estimations of the T-stress for small cracks at notches

The T-stress is increasingly being recognized as an important additional stress field characterizing parameter in the analyses of cracked bodies. Using T-stress as the constraint parameter, the framework of failure assessments including the constraint effect has been established; and the effect of T-stress on fatigue crack propagation rate has been investigated by several researchers. In this paper, a simple method for determining the T-stress for small notch-emanating cracks is presented. First, the background on the T-stress calculation using the superposition principle and the similarities between the elastic notch-tip stress fields described by two parameters: the stress concentration factor (K_t) and the notch-tip radius (ρ), are summarized. Then, the method of estimating T-stress for both short and long cracks at the notches is presented. The method is used to predict T-stress solutions for cracks emanating from an internal hole in a wide plate, and cracks emanating from an U-shaped edge notch in a finite thickness plate. The results are compared to the T-stress results in the literature, and the T-stresses solutions obtained from finite element analysis. Excellent agreements have been achieved for small cracks. The method presented here can be used for a variety of notch crack geometries and loading conditions.     

Stress intensity factors for a wide range of long-deep semi-elliptical surface cracks, partly through-wall cracks and fully through-wall cracks in tubular members

To calculate the rate of fatigue crack growth in tubular members, one approach is to make use of the fracture mechanics based Paris law. Stress intensity factors (SIF) of the cracked tubular members are prerequisite for such calculations. In this paper, stress intensity factors for circumferential deep semi-elliptical surface crack (a/t > 0.8), semi-elliptical partly through-wall crack and fully through-wall crack cracks in tubular members subjected to axial tension are presented. The work has produced a comprehensive set of equations for stress intensity factors as a function of a/T, c/πR and R/T for deep surface cracks. For the partly through-wall cracks and fully through-wall cracks, two sets of bounding stress intensity factor equations were produced based on which all stress intensity factors within the range of parameters can be obtained by interpolation.     

Surface flaws behavior under tension,bending and biaxial cyclic loading

Fatigue surface crack growth and in-plane and out-of-plane constraint effects are studied through experiments and computations for the aluminum alloy D16T. A tension/bending central notched plate and cruciform specimens under different biaxial loadings with external semi-elliptical surface cracks are studied. The variation of the fatigue crack growth rate and surface crack paths is studied under cyclic tension, bending and biaxial tension–compression loading. For the experimental surface crack paths in the tested specimens, the T-stress, out-of-plane T_z factor, local triaxiality parameter h and the governing parameter for the 3D-fields of the stresses and strains at the crack tip in the form of the I_n-integral are calculated as a function of the aspect ratio by finite element analysis to characterize the constraint effects along the semi-elliptical crack front. The plastic stress intensity factor approach is applied to the fatigue crack growth on the free surface, as well as at the deepest point of the semi-elliptical surface crack front, of the tested tension/bending plate and cruciform specimens. From the results, characteristics of the fatigue surface crack growth rate as a function of the loading conditions are established.     

Three finite element analysis of semi-elliptical crack in high density poly-ethylene pipe subjected to internal pressure

The principal objective of this work is to analyze the severity of semi-elliptical crack defects and to study the degree of damage in the equipment under internal pressure during the crack propagation. The semi-elliptical cracks are considered in this work located in different position in the wall of poly-ethylene pipe. The tree finite element method based on the computation of the J integral was used to analyze the fracture behaviour of these structures. The effect of the position shape and size of the crack on the J integral was highlighted. The effects of strain rate and the temperature on the J integral values were also examined. The obtained results show that, whatever the material (strain rate) for a semi-elliptical crack, the J integral value has not an important variation with respect to the crack size. However, the energy for axial crack is more important compared to circumferential crack. The effect of the depth of the crack becomes important when the ratio (a/t) reaches a critical value of 0.6 (a/t = 0.6), especially when the ratio a/c is weak (semi-elliptical crack). We recall finally, that the temperature effect on circumferential cracks behaviour is more important compared to the axial cracks. It is also shown that in the wall of pipe, the internal cracks are more dangerous than the external cracks.     

Elastic T-stress for cracks in test specimens subjected to non-uniform stress distributions

Finite element analyses have been conducted to calculate elastic T-stress solutions for cracked test specimens. The T-stress solutions are presented for single edge cracked plates, double edge crack plates and centre cracked plates. Uniform, linear, parabolic or cubic stress distributions were applied to the crack face. The results for uniform and linear stress distributions were used to derive weight functions for T-stress for the corresponding specimens. The weight functions for T-stress are then verified against several linear and non-linear stress distributions. The derived weight functions are suitable for the T-stress calculation for cracked specimens under any given stress field.     

3-D stress intensity factors for arrays of inner radial lunular or crescentic cracks in a typical spherical pressure vessel

Many spherical pressure vessels are manufactured by methods such as the integrated hydro-bulge forming (IHBF) method, where the sphere is composed of a series of double curved petals welded along their meridional lines. Such vessels are susceptible to multiple radial cracking along the welds. For fatigue life assessment and fracture endurance of such vessels one needs to evaluate the stress intensity factors (SIFs) distribution along the fronts of these cracks. However, to date, only two 3-D solutions for the SIF for one inner semi-elliptical crack in thin or thick spheres are available, as well as 2-D SIFs for one through-the-thickness crack in thin spherical shells. In the present paper, mode I SIF distributions for a wide range of lunular and crescentic cracks are evaluated. The 3-D analysis is performed, via the FE method employing singular elements along the crack front, for a typical spherical pressure vessel with outer to inner radius ratios of η = R_o/R_i = 1.1. SIFs are evaluated for arrays containing n = 1-20 cracks; for a wide range of crack depth to wall thickness ratio, a/t, from 0.025 to 0.95; and for various ellipticities of the crack, i.e., the ratio of crack depth to semi crack length, a/c, from 0.2 to 1.5. The obtained results clearly indicate that the SIFs are considerably affected by the three-dimensionality of the problem, and the following parameters: the number of cracks in the array-n, the relative crack depth a/t, and the crack ellipticity a/c.     

Elastic T-stress solutions for flat elliptical cracks under tension and bending

Exact solutions for elastic T-stress of a flat elliptical crack in an infinite body under tension and bending are obtained in this paper. Many papers have been devoted to the problems for elliptical cracks in an elastic medium, but all their attention has been concentrated on the determination of stress intensity factors. In the current paper, elastic T-stress solutions are derived by means of the potential method and a specific collection of harmonic functions. The formulas of the elastic T-stress for a penny-shaped crack [Wang X. Elastic T-stress solutions for penny-shaped cracks under tension and bending. Engng Fract Mech 2004;71:2283-98] follow from the present results as a special case. It is obtained that under tension loading, the elastic T-stress is always compressive along the elliptical crack front. In both tension and bending cases, T-stress essentially depends on the Poisson’s ratio of the material, a parametric angle and semi-axes of the ellipse.     

Two-parameter fracture assessment of surface cracked cylindrical shells during collapse

The present study addresses the use of CTOD and T-stress in fracture assessments of surface cracked shell structures. A new software is developed for this purpose, denoted LINK_pipe. It is based on a combination of a quadrilateral assumed natural deviatoric strain thin shell finite element and an improved linespring finite element. Plasticity is accounted for using stress resultants. A power law hardening model is used for shell and linespring materials. A co-rotational formulation is employed to represent nonlinear geometry effects. With this, one can carry out nonlinear fracture mechanics assessments in structures that show instabilities due buckling (local/global), ovalisation and large rigid body motion. Many constraint-measuring parameters have been proposed, with the Q-parameter or the T-stress being the most popular ones. Solid finite element meshing for complex structures such as pipes containing semi-elliptical surface cracks in order to compute Q is at present not a feasible approach. However, shell structures are most conveniently meshed with shell finite elements, and the linespring finite element is a natural way of accounting for surface cracks. The T-stress is readily obtained from the linespring membrane force and bending moment along the surface crack. In this study we present a new approach to analyse cracked shell structures subjected to large geometric changes. By numerical examples it is shown how geometric instabilities and fracture compete as governing failure mode.     

Elastic T stress estimates for circumferential surface‐cracked cylinders*

On the basis of detailed three‐dimensional (3D) elastic finite element (FE) analyses, this paper provides tractable approximations for elastic T stress solutions for circumferential inner‐surface cracks in cylinders. Internal pressure and global bending moment were considered. The FE model and analysis procedure employed in the analysis were verified using existing solutions for both elastic stress intensity factor and T stress. To cover a practical range, three different values of the ratio of the mean radius of cylinder to the thickness, R_m/t, were selected; furthermore, four different values of the ratio of the crack depth to the thickness, a/t, ranging from 0.1 to 0.75 and three different values of θ/π ranging from 0.1 to 0.4 were selected. On the basis of FE analyses results, polynomial approximations were proposed at three different locations: surface point, middle point and deepest point. On the basis of the detailed 3D elastic FE analysis, the solutions presented are believed to be the most accurate, and thus provide valuable information for structural integrity assessment considering a crack‐tip constraint.     

Two-parameter characterization of elastic-plastic crack front fields: Surface cracked plates under tensile loading

In this paper the J-Q two-parameter characterization of elastic-plastic crack front fields is examined for surface cracked plates under uniaxial and biaxial tensile loadings. Extensive three-dimensional elastic-plastic finite element analyses were performed for semi-elliptical surface cracks in a finite thickness plate, under remote uniaxial and biaxial tension loading conditions. Surface cracks with aspect ratios a/c = 0.2, 1.0 and relative depths a/t = 0.2, 0.6 were investigated. The loading levels cover from small-scale to large-scale yielding. In topological planes perpendicular to the crack fronts, the crack stress fields were obtained. In order to facilitate the determination of Q-factors, modified boundary layer analyses were also conducted. The J-Q two-parameter approach was then used in characterizing the elastic-plastic crack front stress fields along these 3D crack fronts. Complete distributions of the J-integral and Q-factors for a wide range of loading conditions were obtained. It is found that the J-Q characterization provides good estimate for the constraint loss for crack front stress fields. It is also shown that for medium load levels, reasonable agreements are achieved between the T-stress based Q-factors and the Q-factors obtained from finite element analysis. These results are suitable for elastic-plastic fracture mechanics analysis of surface cracked plates.     

The in-plane and out-of-plane stress constraint factors and K-T-Tz description of stress fields near the border of a quarter-elliptical corner crack

The elastic T-stress and stress intensity factor K for quarter-elliptical corner cracks have been investigated in elastic plates by detailed three-dimensional finite-element calculations. The distributions of normalized K and T-stress have been obtained along the crack front with aspect ratios (a/c) of 0.2, 0.3, 0.4, 0.5, 0.6, 0.8 and 1.0, and far-field tension and the effect of Poisson's ratio have also been considered. The normalized K increases and the normalized T-stress decreases with the increase of Poisson's ratio v. For v= 0.3, the normalized K gradually increases in the range of crack-face angle φ≥ 22.5° and decreases in the range of φ≤ 22.5° with the increase of a/c. The normalized T-stress increases in the beginning and then decreases with increasing φ except for a/c= 0.2 and a/c= 0.3. By fitting the numerical results with the least squares method, empirical formulae have been given for the convenience of engineering applications. Combining with the corresponding out-of-plane constraint factor T_z, the three-parameter K-T-T_z approach has been provided, which can accurately describe the stress field around the crack front.     

T-stress of cracks loaded by near-tip tractions

T-stress solutions were derived for tractions acting on the crack-faces near a crack tip. Such solutions are of interest for the determination of the leading term of a weight function representation of T-stresses and the computation of an “intrinsic” T-stress for cracks growing in a material with a rising crack growth resistance. First, the type of the Green’s function for T-stresses is theoretically established. Then, results of finite element computations are reported for edge-cracked bars, DCB and CT specimens, which are suited for the determination of the first series term. As an application of the Green’s functions, the T-stresses caused by bridging interactions very close to the crack tip are computed.     

J integral solution for semi-elliptical surface crack in high density poly-ethylene pipe under bending

The goal of this work is to analyse the severity of semi-elliptical crack defects and to study the degree of damage in the poly-ethylene pipe in bending during the crack propagation. The semi-elliptical cracks are considered in this work located in different position in the wall of the pipe. The three finite element method based on the computation of the J integral was used to analyse the fracture behaviour of these structures. The effect of the position, shape and size of the crack on the J integral values was highlighted. The effects of strain rate and the temperature on the J integral values were also examined. The obtained results show that the strain rates have a strong influence on the J integral values especially for circumferential crack at higher bending moment. However, the energy for circumferential crack is more important compared to axial crack. The effect of the depth of the crack becomes important when the ratio (a/t) reaches a critical value of 0.6 (a/t = 0.6), especially when the ratio a/c is weak (semi-elliptical crack, a/c = 0.2) where the J integral values becomes independently of the crack depth, this conclusion is opposite to the above for the poly-ethylene pipe subjected to internal pressure. We recall finally, that the temperature effect on circumferential cracks behaviour is more important compared to the axial cracks at critical crack size (a/c = 0.2 and a/t = 0.6). It is also shown that in the wall of pipe, the internal cracks are more dangerous than the external cracks.     

The effective T-stress estimation and crack paths emanating from U-notches

The concept of the T-stress as a local constraint factor has been extended to U-notch tip stress distribution as the effective T-stress. The effective T-stress has been estimated as the average value of the T-stress in the region corresponding to the effective (characteristic) distance ahead of the notch tip. The T-stress is evaluated by finite element method using the experimental load for crack initiation and computing the difference between principal stresses along ligament. A large range of critical effective T-stress values is investigated for different specimen configurations and notch aspect ratios. Crack stabilisation and crack bifurcation for fracture emanating from notches according to the critical effective T-stress is discussed. A model involving the influence of the critical effective T-stress on void growth for ductile failure in the vicinity of the notch tip has been proposed.     

On the full set of elastic T-stress terms of internal elliptical cracks under mixed-mode loading condition

Analytical expressions for the elastic constant stress terms of the asymptotic field, the so called T-stresses, for internal mixed-mode elliptical cracks in infinite homogeneous and isotropic elastic solids are addressed. To solve the problem the mixed-mode crack problem is divided into sub-problems using the superposition method, and each sub-problem is then solved for the asymptotic stress field. Considering the expansion of the local stress field at the crack front, the elastic T-stress terms are derived for each sub-problem. The results are superimposed to give the analytical expressions of the so far missing elastic T-stresses for mixed-mode elliptical cracks.The effect of the T-stresses on the size and shape of the plastic zone at the crack tip is discussed, and analytical results are compared to the ones from finite element analyses, both for the T-stress components and the size of the plastic zone. For an accurate prediction of the plastic zone all singular and constant terms (T-stresses) in the stress expansion formulae should be considered. It is observed that negative T-stresses increase the size of the plastic zone, while positive ones reduce it.     

On the full set of elastic T-stress terms of internal circular cracks under mixed-mode loading conditions

Analytical expressions for all non-singular stress terms of the asymptotic crack tip field, the so-called T-stresses of internal mixed mode circular (penny shaped) cracks in infinite homogeneous and isotropic elastic solids are addressed. To solve the problem the mixed mode crack problem is divided into sub-problems using the superposition method, and each sub-problem is then solved for the asymptotic stress field. Considering the expansion of the local stress field at the crack front, the elastic T-stress terms are derived for each sub-problem. The results are superimposed to give the analytical expressions of the so far missing elastic T-stresses of internal mixed mode penny shaped cracks.The effect of the T-stresses on the size and shape of the plastic zone at the crack tip is discussed, and analytical results are compared to the ones from finite element analyses, both for the T-stress components and the size of the plastic zone. For an accurate prediction of the plastic zone all the singular terms and the constant terms (T-stresses) in the stress expansion formulae should be considered. It is observed that negative T-stresses increase the size of the plastic zone, while positive ones reduce it.     

T-stress near the tips of a cruciform crack with unequal arms

The T-stress near the tips of a crack of cross shape embedded in an isotropic elastic solid is analyzed. The integral transform technique is employed to convert the associated boundary value problem to a system of singular integral equations. According to the stress difference method, T-stresses can be expressed as a sum of an integral involving crack opening displacement (COD) and applied loading at infinity. Obtained results indicate that, in addition to applied loading, T-stresses at the horizontal (vertical) crack tips depend on the COD of the vertical (horizontal) crack surface. COD plays a leading role in determining T-stresses, in particular for a cruciform crack of two crack-arm lengths of the same order. Moreover, T-stresses for a single-crack limiting case can be recovered from the present results as the length of one arm approaches zero. For a biaxial tension of the same magnitude, T-stresses are present for a cruciform crack, but absent for a single crack. Finally, for several cases of interest, T-stresses for a cruciform crack are evaluated and compared with those for a single crack, and the influence of the ratio of two crack-arm lengths b/a and the COD on the T-stress of a cruciform crack is examined.     

Three-dimensional stress state around quarter-elliptical corner cracks in elastic plates subjected to uniform tension loading

Detailed full-field three-dimensional (3D) finite element analyses have been conducted to study the out-of-plane stress constraint factor T_z around a quarter-elliptical corner crack embedded in an isotropic elastic plate subjected to uniform tension loading. The distributions of T_z are studied in the forward section (0° ? θ ? 90°) of the corner cracks with aspect ratios a/c of 0.2, 0.4, 0.5, 0.6, 0.8 and 1.0. In the normal plane of the crack front line, T_z drops radially from Poisson’s ratio at the crack tip to zero beyond certain radial distances. Strong 3D zones (T_z > 0) exist within a radial distance r/a of about 4.6-0.7 for a/c = 0.2-1.0 along the crack front, despite the stress-free boundary conditions far away. At the same radial distance along the crack front in the 3D zones, T_z increases from zero on one free surface to a peak value in the interior, and then decreases to zero on another free surface. The distributions of T_z near the corner points are also discussed. Empirical formulae describing the 3D distributions of T_z are obtained by fitting the numerical results, which prevail with a sufficient accuracy in the valid range of 0.2 ? a/c ? 1.0 and 0° ? θ ? 90° except very near the free surfaces where T_z is extremely low. Combined with the K-T solution, the transition of approximate plane-stress state near the surfaces to plane-strain state in the interior can be characterized more accurately.