Weight function,stress intensity factor and crack opening displacement solutions to periodic collinear edge hole cracks

Periodic collinear edge hole cracks and arbitrary small cracks emanating from collinear holes, which are two typical multiple site damages occurred in the aircraft structures, are studied by using the weigh function method. An explicit closed form weight function for periodic edge hole cracks in an infinite sheet is obtained and further used to calculate the stress intensity factor and crack opening displacement for various loading cases. Compared to finite element method, the present weight function is accurate and highly efficient. The interactions of the holes and cracks on the stress intensity factor and crack opening displacement are quantitatively determined by using the present weight function. An approximate weight function method is also proposed for arbitrary small cracks emanating from multiple collinear holes. This method is very useful for calculating the stress intensity factor for arbitrary small cracks.     

Wide‐range and accurate mixed‐mode weight functions and stress intensity factors for cracked discs

The Wu‐Carlsson displacement‐based weight function method is extended to obtain the mode I and mode II weight functions for the edge‐ and centre‐cracked discs. Compared with Fett's direct adjustment weight functions for the edge‐cracked discs, the present weight functions are more accurate and are applicable for a wider range of crack lengths. Using the present weight functions, extensive and highly accurate mixed‐mode stress intensity factors are obtained for the cracked discs subjected to diametrically compressive forces. Assuming perfect contact and using Coulomb friction law and the present weight functions, the mode II stress intensity factors for the cracked discs with consideration of friction are obtained and widely compared with the corresponding results from finite element analyses.     

Weight functions and strip yield solution for two equal-length collinear cracks in an infinite sheet

The problem of two equal-length collinear cracks in an infinite sheet is treated using the weight function method. Exact weight functions for the inner and outer crack tips are derived based on the crack opening displacement solution for a reference load case. These weight functions are used to calculate stress intensity factors for different load cases, plastic zone sizes and crack tip opening displacements of the strip yield model. The approach is validated by the perfect agreement between the present strip yield model solutions and Collins and Cartwright’s analytical results based on the direct complex stress function formulation.     

Stress intensity factors for cracks in thin plates

This paper presents stress intensity factor solutions for several crack configurations in plates. The loadings considered include internal pressure, and also combined bending and tension. The dual boundary element method is used to model the plate and mixed mode stress intensity factors are evaluated by a crack surface displacement extrapolation technique and the J-integral technique. Several cases including centre crack, edge crack and cracks emanating from a hole in finite width plates are presented.     

A new weight function approach using indirect boundary integral method

A new weight function approach to determine SIF (stress intensity factor) using the indirect boundary integral method has been presented. The crack opening displacement field was represented by one boundary integral in the form of a single-layer potential whose kernel was modified from the fundamental solution. The proposed method enables the calculation of SIF using only one SIF formula without any modification of the crack geometries symmetric in a two-dimensional plane, e.g. a center crack in a plate with or without an internal hole, double edge cracks, circumferential crack or radial cracks in a pipe. The application procedure for this variety of crack geometries is very simple and straghtforward with only one SIF formula. The necessary information in the analysis is two reference SIFs. The analysis results, using several examples, verified that the present closed-form solution was in good agreement with those of the literature and applicable to various crack geometries.     

Prediction of fatigue crack growth path in mechanical joints using a weight function approach

Cracks often initiate from the mechanical joints which are widely used in structural components. It has been reported that cracks in mechanical joints are under mixed‐mode condition and there is a critical angle at which mode I stress intensity factor becomes maximum. The crack propagates in an arbitrary direction and the prediction of fatigue crack growth path is needed to provide against crack propagation and examine safety. In this study, mixed‐mode fatigue crack growth tests are performed for horizontal and critical inclined cracks in mechanical joints. Fatigue crack growth paths are predicted using a weight function approach and maximum tangential stress criterion.     

A new weight function approach using indirect boundary integral method

A new weight function approach to determine SIFs (stress intensity factors) using the indirect boundary integral method has been presented. The crack opening displacement field was represented by one boundary integral term in the form of a single-layer potential whose kernel was modified from the fundamental solution. The proposed method enables the calculation of SIFs using only one SIF solution, without any modification for the crack geometries symmetric in the two-dimensional plane, e.g. a center crack in a plate with or without an internal hole, double edge cracks, circumferential cracks or radial cracks in a pipe. The application procedure for this variety of crack geometries is very simple and straightforward with only one SIF solution. The necessary information in the analysis is two reference SIFs. The analysis results using several examples verified that the present closed-form solution was in good agreement with those of the literature and applicable to various crack geometries.     

Calculation of mixed mode stress intensity factors for an elliptical subsurface crack under arbitrary normal loading

Normal loading causes mixed fracture modes in an elliptical subsurface crack because of the nonsymmetrical geometry with respect to the crack face. In this paper, mixed mode weight functions (MMWFs) for elliptical subsurface cracks in an elastic semi‐infinite space under normal loading are derived. Reference mixed mode stress intensity factors (MMSIFs), calculated by finite element analysis, under uniform normal loading are used to derive MMWFs. The cracks have aspect ratios and crack depth to crack length ratios of 0.2–1.0 and 0.05 to infinity, respectively. MMWFs are used to calculate MMSIFs for any point of the crack front under linear and nonlinear two‐dimensional (2D) loadings. So, in order to evaluate the fatigue crack growth phenomenon under complicated 2D stress distributions, MMWFs can be easily used. The comparison between the MMSIFs obtained from the MMWFs and finite element analysis indicates high accuracy.     

Universal features of weight functions for cracks in mode I

An analysis of known analytical and numerical weight functions for cracks in mode I has revealed that they all have a similar singular term and that it is possible to approximate them with one universal expression with three unknown parameters. The unknown parameters can be determined directly from reference stress intensity factor expressions without using the crack opening displacement function. The universal weight function expression, with suitable reference stress intensity factors, was used to derive the weight functions for internal and external radial cracks in a thick cylinder. These weight functions were then further used to calculate the stress intensity factors for radial cracks in a cylinder subjected to several nonlinear stress fields and were compared to available numerical data.     

Stress intensity factors for corner cracks in single‐edge notch bend specimen by a three‐dimensional weight function method

A three‐dimensional (3D) weight function method is employed to calculate stress intensity factors of quarter‐elliptical corner cracks at a semi‐circular notch in the newly developed single‐edge notch bend specimen. Corner cracks covering a wide range of geometrical parameters under pin‐loading and remote tension conditions are analysed. Stress intensity factors from the 3D weight function analysis agree well with ABAQUS‐Franc3D finite element results. An engineering similitude approach previously developed for the half‐elliptical surface crack in single‐edge notch bend specimen is also applied to the present corner crack configuration. The results compare well with those from the present weight function analysis.     

Dynamic stress intensity factors due to scattering of harmonic waves by a crack in an infinite medium

An analytical method for calculating dynamic stress intensity factors in the mixed mode (combination of opening and sliding modes) using complex functions theory is presented. The crack is in infinite medium and subjected to the plane harmonic waves. The basis of the method is grounded on solving the two‐dimensional wave equations in the frequency domain and complex plane using mapping technique. In this domain, solution of the resulting partial differential equations is found in the series of the Hankel functions with unknown coefficients. Applying the boundary conditions of the crack, these coefficients are calculated. After solving the wave equations, the stress and displacement fields, also the J‐integrals are obtained. Finally using the J‐integrals, dynamic stress intensity factors are calculated. Numerical results including the values of dynamic stress intensity factors for a crack in an infinite medium subjected to the dilatation and shear harmonic waves are presented.     

Stress intensity factors of cracks emanating from a surface square hole in infinite body in tension

This paper deals with such a kind of surface crack problem with a same depth (called a liked‐plane crack problem for short). Based on the previous investigations on an internal rectangular crack and a surface rectangular crack in an infinite solid in tension and the hybrid displacement discontinuity method, a numerical approach for the liked‐plane crack problem is presented. Numerical examples are given to illustrate the numerical approach is simple, yet accurate for calculating the stress intensity factors (SIFs) of the liked‐plane crack problem. Specifically, SIFs of a pair of cracks emanating from a surface square hole in an infinite body in tension are investigated in detail.     

On the computation of two-dimensional stress intensity factors using the boundary element method

General two-dimensional linear elastic fracture problems are investigated using the boundary element method. The √r displacement and 1/√r traction behaviour near a crack tip are incorporated in special crack elements. Stress intensity factors of both modes I and II are obtained directly from crack-tip nodal values for a variety of crack problems, including straight and curved cracks in finite and infinite bodies. A multidomain approach is adopted to treat cracks in an infinite body. The body is subdivided into two regions: an infinite part with a finite hole and a finite inclusion. Numerical results, compared with exact solution whenever possible, are accurate even with a coarse discretization.     

Determination of crack surface displacements for cracks emanating from a circular hole using weight function method

Cracks emanating from a circular hole are of significant engineering importance, especially in aerospace industry. Accurate determination of key fracture mechanics parameters is essential for damage tolerance design and fatigue life predictions. The purpose of this paper is to provide an efficient and accurate closed‐form weight function approach to the calculation of crack surface displacements for radial crack(s) emanating from a circular hole in an infinite and finite‐width plate. Results were presented for two loading conditions: remote applied stress and uniform stress segment applied to crack surfaces, and extensively compared to recent studies using other methods in the literature. Both single and double radial cracks were considered, and also the effect of finite plate width on crack surface displacements has been investigated. A brief assessment was made on an engineering estimation of displacements based on a correction of stress intensity factor ratio. It has been demonstrated that the Wu‐Carlsson closed‐form weight functions are very efficient, accurate and easy‐to‐use for calculating crack surface displacements for arbitrary load conditions. The method will facilitate fatigue crack closure and other fracture mechanics analyses where accurate crack surface displacements are required.     

Some experimental observations on crack closure and crack-tip plasticity

This study presents a methodology for evaluating crack closure and the effect of crack-tip plasticity on stress intensity. Full-field displacement maps obtained by digital image correlation are used to obtain the mixed-mode, crack-driving force. The methodology allows the quantification of the effect of a range of contact phenomena: effects arising from interlocking, plastic deformation of crack face asperities and wedging generated as a consequence of sliding displacements of fatigue cracks have been identified. By evaluating the effective crack-tip stress intensity factor, crack opening levels can be quantified for both mode I and mode II. Moreover, the approach can take into account plasticity effects local to the crack in determining the stress intensity factor. All the information can be extracted in a non-contacting fashion with equipment that can be easily incorporated into industrial environments.     

An elastic strip with periodic surface cracks under thermal shock

The problem of two periodic edge cracks in an elastic infinite strip located symmetrically along the free boundaries under thermal shock is investigated. It is assumed that the infinite strip is initially at constant temperature. Suddenly the surfaces containing the edge cracks are quenched by a ramp function temperature change. Very high tensile transient thermal stresses arise near the cooled surface resulting in severe damage. The degree of the severity for a subcritical crack growth mode is measured by determining the stresses intensity factors. The thermoelastic problem is treated as uncoupled quasi-static. The superposition technique is used to solve the problem. The thermal stresses obtained from the uncracked strip with opposite sign are utilized as the only external loads to formulate the perturbation problem. By expressing the displacement components in terms of finite and infinite Fourier transforms, a hypersingular integral equation is derived with the crack surface displacement as the unknown function. Numerical results for stress intensity factors are carried out and presented as a function of time, cooling rate, crack length, and periodic crack spacing.     

Numerical and experimental evaluation of SIF for threaded connectors

This paper presents a methodology for fatigue crack growth analysis in tubular threaded connectors. A solution for stress intensity factor for semi-elliptical surface cracks emanating from a thread root in a screw connector is also discussed in the paper. The solution is based on a mixed approach incorporating weight function and finite element methods. The weight functions used are the universal functions for cracks in mode I and these are linked with a thread through-thickness stress distribution obtained from finite element analysis to produce a stress intensity factor for a crack at the critical tooth of a thread. The resulting crack growth data are then validated experimentally.     

The residual stress intensity factors for cold‐worked cracked holes: a technical note

Cold‐working of riveted holes reduces the stress intensity factor associated with cracks that may develop at the hole boundary, by creating a compressive residual stress field. The residual stress field is determined using the finite‐element method and the reduction of the stress intensity factor for different values of the interference is evaluated with the weight function method, in the case of an infinite plate made from an elastic–perfectly plastic material, and having a hole with two symmetrical cracks. Once the weight function of the structure is known, further calculation of the stress intensity factors for different loadings such as a remote uniform stress, or a point load that simulates the action of the rivet can be performed without difficulty.     

Stress intensity factors for cracks emanating from a triangular or square hole in an infinite plate by boundary elements

This paper concerns stress intensity factors of cracks emanating from a triangular or square hole in an infinite plate subjected to internal pressure calculated by means of a boundary element method, which consists of constant displacement discontinuity element presented by Crouch and Starfied and crack tip displacement discontinuity elements proposed by the author. In the boundary element implementation the left or the right crack tip displacement discontinuity element is placed locally at corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Numerical examples are included to show that the method is very efficient and accurate for calculating stress intensity factors of plane elasticity crack problems. Specifically, the numerical results of stress intensity factors of cracks emanating from a triangular or square hole in an infinite plate subjected to internal pressure are given.     

Thermal striping fatigue damage in multiple edge-cracked geometries

Thermal fatigue striping damage may be caused when incompletely mixed hot and cold fluid streams pass over the surface of a component or structure containing a defect. Stress intensity factor (SIF) fluctuations are developed in response to the surface temperature fluctuations. An existing methodology for the analysis of striping damage in geometries containing a single edge‐crack geometry is extended to such an analysis of multiple edge cracks. SIFs are calculated as functions of crack depth, when an edge‐cracked plate and semi‐infinite solid, each containing multiple cracks, are subjected to thermal striping. The effect of various restraint conditions and striping frequencies on the SIF values for a stainless steel plate is examined. The degree of conservatism is shown when an assessment of thermal fatigue striping damage is based on a single, rather than multiple, crack analysis. Accurate curve fits are developed resulting in practical weight functions for an edge‐cracked plate and semi‐infinite solid.