Calculation of stress intensity factors for a longitudinal semi-elliptical crack in a finite-length thick-walled cylinder

The purpose of this paper is to present the effect of finite boundary on the stress intensity factor of an internal semi-elliptical crack in a pressurized finite-length thick-walled cylinder  ( R _i/ t = 4)  . The three-dimensional finite element method, in conjunction with the weight function method, is used for computing the stress intensity factor at the deepest and surface points of an axial semi-elliptical crack in a cylinder. The transition aspect ratios, the aspect ratios in which the maximum stress intensity factor translates from the deepest to the surface points of the crack, are calculated for different relative depths and cylinder lengths. The results show that the stress intensity factor increases as the cylinder length decreases, especially at the corner point of the crack compared with the deepest point. The major advantage of this paper is that a closed-form expression is extracted for the stress intensity factor at the surface point of a semi-elliptical crack, which experiences higher changes due to the effect of the finite boundary of the cylinder.     

Variation of mixed modes stress intensity factors of an inclined semi-elliptical surface crack

In this paper, a singular integral equation method is applied to calculate the stress intensity factor along crack front of a 3D inclined semi-elliptical surface crack in a semi-infinite body under tension. The stress field induced by displacement discontinuities in a semi-infinite body is used as the fundamental solution. Then, the problem is formulated as a system of integral equations with singularities of the form r ^–3. In the numerical calculation, the unknown body force doublets are approximated by the product of fundamental density functions and polynomials. The results show that the present method yields smooth variations of mixed modes stress intensity factors along the crack front accurately for various geometrical conditions. The effects of inclination angle, elliptical shape, and Poisson's ratio are considered in the analysis. Crack mouth opening displacements are shown in figures to predict the crack depth and inclination angle. When the inclination angle is 60 degree, the mode I stress intensity factor F _I has negative value in the limited region near free surface. Therefore, the actual crack surface seems to contact each other near the surface.     

Closed‐form thermal stress intensity factors for an internal circumferential crack in a thick‐walled cylinder

In this paper the method of weight functions is employed to calculate the stress intensity factors for an internal circumferential crack in a thick‐walled cylinder. The pressurized cylinder is also subjected to convection cooling on the inner surface. Finite element method is used to determine an accurate weight function for the crack and a closed‐form thermal stress intensity factor with the aid of the weight function method is extracted. The influence of crack parameter and the heat transfer coefficient on the stress intensity factors are determined. Comparison of the results in the special cases with those cited in the literature and the finite element data shows that the results are in very good agreement.     

Stress intensity factor calculation of steel-lined hoop-wrapped cylinders with internal semi-elliptical circumferential crack

The purpose of this paper is to present mode I stress intensity factor for a circumferential semi-elliptical crack on the inner surface of a hoop-wrapped steel-lined CNG cylinder. The stress intensity factors along the crack front are directly computed by 3D finite element method for a wide range of variations of the crack geometry. Also influence of many parameters such as cylinder internal pressure, composite layer thickness, composite material properties and undertaking Auto-Frettage pressure are studied on the stress intensity factor of the crack and some conclusive results are drawn. For the sake of validation of the results and because of lack of the results for a circumferential semi-elliptical crack in the literature, a semi-elliptical axial crack in a composite hoop-wrapped cylinder has been modeled and the results have been compared with those in the literature showing a good agreement.     

Stress intensity factors of semi-elliptical surface cracks in pressure vessels by global-local finite element methodology

In the present study the problem of calculation of the stress intensity factors (SIF) of semi-elliptical cracks located in the stress concentration areas of a pressure vessel is numerically solved by advanced global-local finite element (FE) analysis. The common characteristic of the cases solved is that the stress field at the crack area varies along the axial, the circumferential, as well as, the through-the-thickness directions. SIF solutions for such problems are not available, neither analytically, nor numerically, as the currently existing solutions in the literature (numerical results, Newman-Raju empirical equations, weight function solutions, etc.) are only valid for uniform stress distribution along the axial and circumferential directions of the pressure vessel and allow variation only through-the-thickness. The crack locations considered are the intersection of the cylinder to a nozzle and the connection of the cylinder with its hemi-spherical ends. The stress intensity factors are presented in a suitable table format for various geometrical configurations of both the pressure vessel and the semi-elliptical crack, thus providing a useful tool for the fracture mechanics design of cracked pressure vessels. The modeling details of the sub-structuring methodology, employed in the analysis, are extensively discussed and the numerical approach is proven to be very efficient for the SIF calculation of pressure vessel semi-elliptical cracks.     

Analysis of thermal stress intensity factors for cracked cylinders using weight function method

In this paper a general weight function was derived to evaluate the thermal stress intensity factors of a circumferential crack in cylinders. The weight function derived is valid for a wide range of thin- to thick-walled cylinders and relative crack depth. Closed-form stress intensity factor based on the weight function method was derived as a function of the Biot number and relative depth and various inner-to-outer radius ratios of cylinders. The accuracy of the analysis has been examined using the finite element method results and were compared to existing solutions for uniform loading in the literature for special geometries, indicating an excellent agreement.     

Stress intensity factor solution for a semi-elliptical crack in an arbitrarily distributed stress field

An equation for the stress intensity factor (SIF) for semi-elliptical crack has been developed. It is based on the Newman-Raju's solution for the crack in a plate under bending or tension. The equation can be applied when a stress distribution is described by a power function. Using the approach outlined, the SIF for a surface crack in a T-butt welded connection has been estimated. The results obtained can be used in a fracture-mechanics-based fatigue analysis.     

Stress intensity factors of two diametrically opposed edge cracks in a thick-walled functionally graded material cylinder

A method is developed to evaluate stress intensity factors for two diametrically-opposed edge cracks emanating from the inner surface of a thick-walled functionally graded material (FGM) cylinder. The crack and the cylinder inner surfaces are subjected to an internal pressure. The thermal eigenstrain induced in the cylinder material due to nonuniform coefficient of thermal expansion after cooling from the sintering temperature is taken into account. First, the FGM cylinder is homogenized by simulating its nonhomogeneous material properties by an equivalent eigenstrain, whereby the problem is reduced to the solution of a cracked homogenized cylinder with an induced thermal and an equivalent eigenstrains and under an internal pressure. Then, representing the cracks by a continuous distribution of edge dislocations and using their complex potential functions, generalized formulations are developed to calculate stress intensity factors for the cracks in the homogenized cylinder. The stress intensity factors calculated for the cracks in homogenized cylinder represents the stress intensity factors for the same cracks in the FGM cylinder. The application of the formulations are demonstrated for a thick-walled TiC/Al₂O₃ FGM cylinder and some numerical results of stress intensity factors are presented for different profiles of material distribution in the FGM cylinder.     

Axisymmetric crack problem of thick-walled cylinder with loadings on crack surfaces

This study is concerned with the fracture of an infinite thick-walled cylinder. The inner surface of the cylinder is stress free and the outer is rigidly fixed. The cylinder having a ring-shaped crack located at the symmetry plane is subjected to distributed compressive load on its surfaces. The Hankel and Fourier transform techniques are used for the solution of the field equations. By applying the boundary conditions, the singular integral equation in terms of crack surface displacement derivative is derived. By using an appropriate quadrature formula, the integral equation is reduced to a system of linear algebraic equations. Numerical results are obtained for the stress intensity factors at the edges of the crack, surfaces of which are subjected to uniform, linear and parabolic load distributions.     

Estimations of stress intensity factor for a short crack in an elastic–plastic notch

This note presents a simple method for estimating the stress intensity factor (SIF) for a short crack emanating from an elastic–plastic notch.     

Method for efficient computation of stress intensity factors from weight functions by singular point elimination

Mode I stress intensity factor K_I can be computed by integration of a function representing a stress profile (e.g., variation of stress with depth), modified by an appropriate weight function. Usually, numerical integration is required. However, widely used weight functions cause the end (s) of integration intervals to be singular points, complicating numerical integration. Approaches for computing K_I that deal with singularities by approximating stress profiles by a linear function near a singular point, or transforming a weight function to a form that enables Gauss–Chebyshev integration, are reviewed. As an alternative to those approaches, this study presents a different method for numerical integration involving weight functions. First, a general, variable transformation method to eliminate singularities is introduced. Elimination of singular point enables elementary integration approaches such as Simpson’s rule, as well more involved methods, such as adaptive-Lobatto integration, to be applied. Benchmark tests using a variety of numerical integration formulas show the singular point elimination method to provide accurate, robust and computationally efficient integrations.     

A general weight function approach to compute mode-II stress intensity factors and crack paths for slightly curved or kinked cracks in finite bodies

A simple procedure is proposed that allows computing the stress intensity factors for slightly curved and kinked cracks in finite bodies. Basis of the method is the computation of the stress field around a straight crack under externally applied tractions. Then, this auxiliary crack is replaced by the crack of interest. The stress intensity factors are computed from the stresses caused by the auxiliary crack using the weight function technique. In a practical application of the method, mode-II stress intensity factors are computed for the edge-cracked half-space. From the usual crack path condition, K_II = 0, the paths of propagating cracks under biaxial loading and the critical biaxiality ratio for global directional instability are computed. The results are in very good agreement with finite element computations.     

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.     

Thermal stress intensity factors of an infinite orthotropic layer with a crack

Stress intensity factors are determined for a crack in an infinite orthotropic layer. The crack is situated parallel to the plane surfaces of the layer. Stresses are solved for two kinds of the boundary conditions with respect to temperature field. In the first problem, the upper surface of the layer is heated to maintain a constant temperature T ₀, while the lower surface is cooled to maintain a constant temperature –T ₀. In the other problem, uniform heat flows perpendicular to the crack. The surfaces of the crack are assumed to be insulated. The boundary conditions are reduced to dual integral equations using the Fourier transform technique. To satisfy the boundary conditions outside the crack, the difference in temperature at the crack surfaces and differences in displacements are expanded in a series of functions that vanish outside the crack. The unknown coefficients in each series are evaluated using the Schmidt method. Stress intensity factors are then calculated numerically for a steel layer that behaves as an isotropic material and for a tyrannohex layer that behaves as an orthotropic material.     

Stress intensity factors and weight functions for deep semi-elliptical surface cracks in finite-thickness plates

ABSTRACT Three-dimensional finite element analyses have been conducted to calculate the stress intensity factors for deep semi-elliptical cracks in flat plates. The stress intensity factors are presented for the deepest and surface points on semi-elliptic cracks with a/t -values of 0.9 and 0.95 and aspect ratios ( a/c ) from 0.05 to 2. Uniform, linear, parabolic or cubic stress distributions were applied to the crack face. The results for uniform and linear stress distributions were combined with corresponding results for surface cracks with a/t = 0.6 and 0.8 to derive weight functions over the range 0.05 ≤  a/c  ≤ 2.0 and 0.6 ≤  a/t  ≤ 0.95. The weight functions were then verified against finite element data for parabolic or cubic stress distributions. Excellent agreements are achieved for both the deepest and surface points. The present results complement stress intensity factors and weight functions for surface cracks in finite thickness plate developed previously.     

Stress intensity factors of an inclined elliptical crack near a bimaterial interface

In this paper the stress intensity factors are discussed for an inclined elliptical crack near a bimaterial interface. The solution utilizes the body force method and requires Green’s functions for perfectly bonded semi-infinite bodies. The formulation leads to a system of hypersingular integral equation whose unknowns are three modes of crack opening displacements. In the numerical calculation, unknown body force densities are approximated by using fundamental density functions and polynomials. The results show that the present method yields smooth variations of stress intensity factors along the crack front accurately. Distributions of stress intensity factors are presented in tables and figures with varying the shape of crack, distance from the interface, and elastic modulus ratio. It is found that the inclined crack can be evaluated by the models of vertical and parallel cracks within the error of 24% even for the cracks very close to the interface.     

Crack tip stress intensity factors for a crack emanating from a sharp notch

Accurate calibrations are provided for the crack tip stress intensity factor for a crack of finite length emanating from the symmetric tip of a sharp notch, of arbitrary angle, in terms of the generalised stress intensity quantifying remote loading of the notch. The solution is applied to example problems and shown to be accurate for cases where the crack is much shorter then the notch depth.     

Elastodynamic stress intensity factors for a semi-infinite crack due to three-dimensional concentrated shear loadings on the crack faces

The dynamic stress intensity factors for a semi-infinite crack in an otherwise unbounded elastic body is investigated. The crack is subjected to a pair of suddenly-applied shear point loads on its faces at a distance l away from the crack tip. This problem is treated as the superposition of two problems. The first problem considers the disturbance by a concentrated shear force acting on the surface of an elastic half space, while the second problem discusses a half space with its surface subjected to the negative of the tangential surface displacements induced by the first problem in the front of the crack edge. A fundamental problem is proposed and solved by means of integral transforms together with the application of the Wiener–Hopf technique and Cagniard–de Hoop method. Exact expressions are then derived for the mode II and III dynamic stress intensity factors by taking integration over the fundamental solution. Some features of the solutions are discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Stress intensity factors for internal circular cracks in fibers under tensile loading

Stress intensity factors for inner circular cracks placed eccentrically in a fiber with round cross section were computed and are presented in this paper in both analytical and graphical form. The crack plane was perpendicular to the fiber axis and remote tensile loading was assumed. The stress intensity factors were numerically computed using the finite element method. Mesh objectivity and some other aspects of computational precision are considered. The asymptotic behaviour when the crack size and the ligament depth vanish were considered in order to formulate accurate interpolation expressions.     

The application of asymptotic solutions to a semi-elliptical crack at the root of a notch

This paper presents an approximate method based on asymptotic solutions for estimating the stress intensity factor K for semi-elliptic surface cracks at stress concentrations. The proposed equation for estimating K makes use of the near-notch and remote-notch solution to interpolate over the entire range from shallow to deep cracks. The near-notch solution is obtained by means of the stress concentration factor. For cracks located in the remote stress field, K is obtained by considering the crack to be located in a smooth plate with a crack depth equal to the sum of the notch depth and the actual crack depth. The accuracy of the predictions is assessed using numerical calculations and solutions found in the literature.