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.     

Weight functions for bimaterial interface cracks

An efficient approach using the analytically decoupled near-tip displacement solution for bimaterial interface cracks presented in this paper involves: (1) the calculation of the decoupled strain energy release rates G _I and G _II associated respectively with the decoupled stress intensity factors K _I and K _II and (2) the extension of Rice's displacement derivative representation of Bueckner's weight function vectors beyond the homogeneous media. It is shown that the stress intensity factors for a bimaterial interface crack predicted by the present approach agree very well with those solutions available in the literature. The computational efficiency is enhanced through the use of singular elements in the crack-tip neighborhood.As reported in the homogeneous case, the calculated weight function for a bimaterial interface crack is load-independent but depends strongly on geometry and constraint conditions. Due to the coupling nature of the stress intensity factors of a bimaterial interface crack, the invariant characteristics of the dimensionless weight function vectors are different from those of a crack in homogeneous material. In addition, the elastic constants of two constituents can significantly alter the weight function behavior for a cracked bimaterial medium.Due to the load-independent characteristic of the weight functions, the stress intensity factors for a bimaterial interface crack can be obtained accurately and inexpensively by performing the sum of worklike products between the applied loads and the weight functions for the cracked bimaterial body under any loading conditions once the weight functions are explicitly predetermined. The same calculation can also be applied for the identical cracked bimaterial medium with different constraint conditions by including the self-equilibrium forces that contain both the external loads and the reaction forces induced at the constraint locations. Moreover, the physical interpretation of the weight functions can provide a guidance for damage tolerant design application.     

Stress intensity factors and T-stresses for offset double edge-cracked plates under mixed-mode loadings

A particular weight function method is used in this study to determine the stress intensity factors (SIFs) and T-stresses for offset double edge-cracked plates (ODECPs). By using reference loading conditions prescribed on the crack flanks for finite element analyses, the coefficients of weight functions are derived and compiled in the form of tables. With the weight functions, the SIFs and T-stresses for several loading cases are calculated. The results compare well with those obtained using the displacement method. Applying the derived weight functions, the SIFs and T-stresses for ODECPs under several loading conditions are determined. The results can be used as references for related applications.     

Weight functions for the determination of stress intensity factor and T‐stress for semi‐elliptical cracks in finite thickness plate

This paper presents the application of weight function method for the calculation of stress intensity factors (K) and T‐stress for surface semi‐elliptical crack in finite thickness plates subjected to arbitrary two‐dimensional stress fields. New general mathematical forms of point load weight functions for K and T have been formulated by taking advantage of the knowledge of a few specific weight functions for two‐dimensional planar cracks available in the literature and certain properties of weight function in general. The existence of the generalised forms of the weight functions simplifies the determination of specific weight functions for specific crack configurations. The determination of a specific weight function is reduced to the determination of the parameters of the generalised weight function expression. These unknown parameters can be determined from reference stress intensity factor and T‐stress solutions. This method is used to derive the weight functions for both K and T for semi‐elliptical surface cracks in finite thickness plates, covering a wide range of crack aspect ratio (a/c) and relative depth (a/t) at any point along the crack front. The derived weight functions are then validated against stress intensity factor and T‐stress solutions for several linear and nonlinear two‐dimensional stress distributions. These derived weight functions are particularly useful for the development of two‐parameter fracture and fatigue models for surface cracks subjected to fluctuating nonlinear stress fields, such as these resulting from surface treatment (shot peening), stress concentration or welding (residual stress).     

A variational technique for boundary element analysis of 3D fracture mechanics weight functions: static

The dual boundary element method coupled with the weight function technique is developed for the analysis of three-dimensional elastostatic fracture mechanics mixed-mode problems. The weight functions used to calculate the stress intensity factors are defined by the derivatives of traction and displacement for a reference problem. A knowledge of the weight functions allows the stress intensity factors for any loading on the boundary to be calculated by means of a simple boundary integration without singularities. Values of mixed-mode stress intensity factors are presented for an edge crack in a rectangular bar and a slant circular crack embedded in a cylindrical bar, for both uniform tensile and pure bending loads applied to the ends of the bars. © 1998 John Wiley & Sons, Ltd.     

Weight functions and stress intensity factors for quarter-elliptical corner cracks in fastener holes

Approximate weight functions for a quarter‐elliptical crack in a fastener hole were derived from a general weight function form and two reference stress intensity factors. Closed‐form expressions were obtained for the coefficients of the weight functions. The derived weight functions were validated against numerical data by comparison of stress intensity factors calculated for several nonlinear stress fields. Good agreements were achieved. These derived weight functions are valid for the geometric range of 0.5 ≤a/c≤ 1.5 and 0 ≤a/t≤ 0.8 and R/t= 0.5; and are given in forms suitable for computer numerical integration. The weight functions appear to be particularly suitable for fatigue crack growth prediction of corner cracks in fastener holes and fracture analysis of such cracks in complex stress fields.     

A direct determination of weight functions for arbitrary three-dimensional cracks – a combined analytical and numerical approach

A combined analytical and numerical method is proposed for the determination of the weight functions of stress intensity factors of cracks in an arbitrary three-dimensional elastic body. Having defined the weight functions for a given geometry of a structure, the stress intensity factors for arbitrary loading conditions can be obtained by a simple inner product of the weight function and a traction vector. Traditionally weight functions are defined in the two ways; the one is defined by the hyper-singular term of the eigen-function expansion of the displacement field of a cracked body, and the other is defined by the variation of displacement field with respect to a virtual extension of a crack. In the present paper, the weight functions for stress intensity factors are defined by applying the Maxwell-Betti's reciprocal theorem to an original problem and the auxiliary problems subjected to three kinds of force-couples acting on the crack surfaces near the limiting periphery of an arbitrary three-dimensional crack. In the present formulation, weight functions can be calculated by using a general-purpose finite element code combined with analytical expressions near the condensation point, where hyper-singularities exist. The validity of the method is confirmed by two- and three-dimensional illustrative examples.     

Stress intensity factors for surface fatigue crack in a round bar under cyclic axial loading

In this paper, the surface fatigue crack growth shape for an initial straight-fronted edge crack in an elastic bar of circular cross-section is determined through experiments under pure fatigue axial loading. Three different initial notch depths are discussed. The relations of the aspect ratio (b/c) and relative crack depth (b/D) are obtained, and it is shown that there is a great difference in the growth of cracks with different initial front shapes and crack depths. Further, using the three-dimensional finite element method, the stress intensity factors (SIFs) are determined under remote uniform tension loading. Since the relationship of b/c and b/D changes during the fatigue crack growth, the SIFs are determined for different surface crack configurations.     

Description of loadings and screenings of cracks with the aid of universal weight functions

The stress intensity vector K _i is defined as the limiting behaviour of the stress near the tip of a crack, the stress components being proportional to r ^?1/2 for any external loading. Internal stresses caused by dislocations show the same power dependence at the crack tip; the stress intensity associated with a loading can thus be screened (or amplified) by a plastic zone. Since for any particular specimen and crack geometry the stress intensity vector must be a functional of the loading and screening which are of vectorial character (lines of forces f _i or dislocations with a Burgers vector b _i) one can define two tensorial weight functions, one for screenings, D _si(x, a), and one for forces, F _si(x, a), so that the stress intensity K _s can be found by integration over the product of weight functions and dislocation or force density. In order to find the weight functions the displacement field and the Airy stress vector must be known for some, completely arbitrary, loading or screening.     

Effect of curvature at the crack tip on the stress intensity factor for curved cracks

In this paper, the numerical solution of the hypersingular integral equation using the body force method in curved crack problems is presented. In the body force method, the stress fields induced by two kinds of standard set of force doublets are used as fundamental solutions. Then, the problem is formulated as a system of integral equations with the singularity of the form r ^–2. In the numerical calculation, two kinds of unknown functions are approximated by the products of the fundamental density functions and power series. The calculation shows that the present method gives rapidly converging numerical results for curved cracks under various geometrical conditions. In addition, a method of evaluation of the stress intensity factors for arbitrary shaped curved cracks is proposed using the approximate replacement to a simple straight crack.     

Stiffness derivative finite element technique to determine nodal weight functions with singularity elements

The use of the stiffness derivative technique coupled with “quarter-point” singular crack-tip elements permits very efficient finite element determination of both stress intensity factors and nodal weight functions. Two-dimensional results are presented in this paper to demonstrate that accurate stress intensity factors and nodal weight functions can be obtained from relatively coarse mesh models by coupling the stiffness derivative technique with singular elements.The principle of linear superposition implies that the calculation of stress intensity factors and nodal weight functions with crack-face loading, σ($\text{r}{}_{\mathrm{s}}$), is equivalent to loading the cracked body with remote loads, which produces σ($\text{r}{}_{\mathrm{s}}$) on the prospective crack face in the absence of crack. The verification of this equivalency is made numerically, using the virtual crack extension technique. Load independent nodal weight functions for two-dimensional crack geometry is demonstrated on various remote and crack-face loading conditions. The efficienct calculation of stress intensity factors with the use of the “uncracked” stress field and the crack-face nodal weight functions is also illustrated.In order to facilitate the utilization of the discretized crack-face nodal weight functions, an approach was developed for two-dimensional crack problems. Approximations of the crack-face nodal weight functions as a function of distance, ($\text{r}{}_{\mathrm{s}}$), from crack-tip has been sucessfully demonstrated by the following equation: $\text{h(a, r}{}_{\mathrm{s}}\text{,) =}\text{A(a)}\text{√r}{}_{\mathrm{s}}\text{+ B(a) + C(a)√r}{}_{\mathrm{s}}\text{+ D(a)r}{}_{\mathrm{s}}$.Coefficients $\text{A(a)}$, $\text{B(a)}$, $\text{C(a)}$ and $\text{D(a)}$, which are functions of crack length ($\text{a}$), can be obtained by least-squares fitting procedures. The crack-face nodal weight functions for a new crack geometry can be approximated using cubic spline interpolation of the coefficients $\text{A, B, C}$ and $\text{D}$ of varying crack lengths. This approach, demonstrated on the calculation of stress intensity factors for single edge crack geometry resulted in total loss of accuracy of less than 1%.     

Micro-crack initiation at tip of a rigid line inhomogeneity in piezoelectric materials

For the square-root singularity shear stress found at the tip of a rigid line inhomogeneity (an anti-crack) in piezoelectric media, one possible way of releasing high strain energy is to initiate a micro-crack at the inhomogeneity tip. In our current study, a dislocation pileup model for micro-crack initiation at the inhomogeneity tip is proposed based on Zener-Stroh crack initiation mechanism. An interesting and important physical result that emerges from the analysis is that the critical stress intensity factor for the anti-crack (the line inhomogeneity) can be related to the fracture toughness of a conventional Griffith crack in the same material. Analytical results further show that under mechanical loading, the critical stress and electric displacement intensity factors of an anti-crack are only related to the corresponding intensity factors of stress and electric displacement of the crack, respectively. While if the anti-crack is under displacement loading (with net dislocation pile-up at the inhomogeneity tip), the critical stress and electric displacement intensity factors of an anti-crack depend on both of the total mechanical dislocations b_T and electricity dislocations b_D.     

Crack propagation from a pre-existing flaw at a notch root. I. Introduction and general form of the stress intensity factors at the initial crack tip

This paper and its companion are devoted to the study of crack kinking from some small pre-existing crack originating from a notch root (the notch root radius being zero). Both the notch boundaries and the initial crack are allowed to be curved; also, the geometry of the body and the loading are totally arbitrary. The ingredients required are knowledge of the stress intensity factors at the initial crack tip and use of a suitable mixed mode propagation criterion. This paper is devoted to the first point, and more specifically to establishing the general (that is, not yet fully explicit) form of the formulae giving these stress intensity factors. The method used is based on changes of scale (homogeneity properties of the equations of elasticity) on the one hand, and on continuity of the displacement and stresses at a given, fixed point with respect to the crack length on the other hand. The formulae derived for the stress intensity factors at the tip of the small crack are of universal value: they apply to any situation, whatever the geometry of the body, the notch and the crack and whatever the loading, the stress intensity factors depending always only upon the `stress intensity factor of the notch' (the multiplicative coefficient of the singular stress field near the notch root in the absence of the crack), the length of the crack, the aperture angle of the notch and the angle between its bisecting line and the direction of the crack.     

Engineerng method of constructing a weight function for plane normal-separation cracks in three-dimensional bodies

On the basis of an analysis of the known solutions for a circular and for a semiinfinite normal-separation crack with straight front in an infinite body we constructed, in analogy to the method of Burns and Oore, a weight function for elliptical cracks. Its use for finding the stress intensity factors K_I at the point of minimal curvature under conditions of uniform loading leads to a maximal error of 10%; at the point of minimal curvature of the crack front the error increases with decreasing ratio of the semiaxes of the ellipse. With the aid of this solution the weight function is found in bounded bodies with a quarter-elliptical, semielliptical, and elliptical crack. A comparison of the data obtained by this method with the values of K_I calculated by the finite-element method under nonuniform loading showed that the suggested method is very accurate.Translated from Problemy Prochnosti, No. 10, pp. 14–22, October, 1992.     

Determination of approximate point load weight functions for embedded elliptical cracks

This paper presents the application of weight function method for the calculation of stress intensity factors in embedded elliptical cracks under complex two-dimensional loading conditions. A new general mathematical form of point load weight function is proposed based on the properties of weight functions and the available weight functions for two-dimensional cracks. The existence of this general weight function form has simplified the determination of point load weight functions significantly. For an embedded elliptical crack of any aspect ratio, the unknown parameters in the general form can be determined from one reference stress intensity factor solution. This method was used to derive the weight functions for embedded elliptical cracks in an infinite body and in a semi-infinite body. The derived weight functions are then validated against available stress intensity factor solutions for several linear and non-linear stress distributions. The derived weight functions are particularly useful for the fatigue crack growth analysis of planer embedded cracks subjected to fluctuating non-linear stress fields resulting from surface treatment (shot peening), stress concentration or welding (residual stress).     

Variations of stress intensity factor of a semi-elliptical surface crack subjected to mixed mode loading

Maximum stress intensity factors of a surface crack usually appear at the deepest point of the crack, or a certain point along crack front near the free surface depending on the aspect ratio of the crack. However, generally it has been difficult to obtain smooth distributions of stress intensity factors along the crack front accurately due to the effect of corner point singularity. It is known that the stress singularity at a corner point where the front of 3 D cracks intersect free surface is depend on Poisson's ratio and different from the one of ordinary crack. In this paper, a singular integral equation method is applied to calculate the stress intensity factor along crack front of a 3-D semi-elliptical surface crack in a semi-infinite body under mixed mode loading. The body force method is used to formulate the problem as a system of singular integral equations with singularities of the form r ⁻³ using the stress field induced by a force doublet in a semi-infinite body as fundamental solution. 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 mixed modes stress intensity factors along the crack front accurately. Distributions of stress intensity factors are indicated in tables and figures with varying the elliptical shape and Poisson's ratio.     

Weight function for a single edge cracked geometry with clamped ends

A single edge cracked geometry with clamped ends is well suited for fracture toughness and fatigue crack growth testing of composites and thin materials. Analysis of fiber bridging phenomenon in the composites and determination of stress intensity factors due to non-uniform stress distributions such as residual and thermal stresses generally require the use of a weight function. This paper describes the development and verification of a weight function for the single edge cracked geometry with clamped ends. Finite element analyses were conducted to determine the stress intensity factors (K) and crack opening displacements (COD) due to different types of stress distributions. The weight function was developed using the K and COD solution for a constant stress distribution. K and COD predicted using this weight function correlated well with the finite element results for non-uniform crack surface stress distributions.     

Edge distance effects on residual stress distribution around a cold expanded hole in Al 2024 alloy

Cold expansion process is a well-known technique for improving the fatigue life of aerospace structures by introducing a compressive residual stress around the fastener holes. However, there are concerns about the residual stress distribution around those holes which are located near the free edges of structure. The purpose of this study is to investigate the influence of edge distance ratio (e/D) on the residual stress distribution around a cold expanded hole in Al 2024 alloy. A two-dimensional finite element simulation was carried out with various degrees of cold expansion and various values of e/D. It was found that for edge distance ratios less than e/D = 3, there are considerable effects on the residual stress profile. Also, the dependency of residual stress distribution on the degree of expansion was reduced significantly for small e/Ds. The results revealed that the bulging of the plate free edge increases with reduction of e/D but in worse cases the maximum bulging at the plate free edge was lower than 3% of the hole radius. The weight function method was then used for determining stress intensity factors for a crack emanating from a cold expanded hole.     

On the fatigue and fracture behavior of necked double shear lugs for aircraft applications

The fatigue and fracture behavior of double shear lugs subjected to axial loading is investigated. The focus is on specific shapes, so-called waisted or necked lugs. These structural components used in aircraft interior are prone to fatigue loads. Three different sizes of necked double shear lugs made of high strength aluminum 2024-T351 and steel 17–4 PH are tested using constant amplitude cyclic loadings with a load ratio R = 0.01. Measurement data is used to identify the number of cycles to crack initiation and final fracture. Fatigue tests show that cracks initiate either at the inside or outside surface of necked lugs. However, no clear dependency on the load amplitude, lug size and material could be found. Numerical simulations using both conventional finite element method (FEM) and extended finite element method (XFEM) are performed to calculate the stress intensity factors (SIFs) for multiple crack lengths of straight and necked double shear lugs. Calculated stress intensity factors for straight lugs fit well to stress intensity factors reported in literature. Stress intensity factor curves of inside and outside cracks of necked lugs plotted with respect to crack length, cross each other, which could have an influence on the fracture behavior observed in fatigue tests.     

Variation of elastic T-stresses along slightly wavy 3D crack fronts

The non-singular terms in the series expansion of the elastic crack-tip stress field, commonly referred to as the elastic T-stresses, play an important role in fracture mechanics in areas such as the stability of a crack path and the two-parameter characterization of elastic-plastic crack-tip deformation. In this paper, a first order perturbation analysis is performed to study some basic properties of the T-stress variation along a slightly wavy 3D crack front. The analysis employs important properties of Bueckner-Rice 3D weight function fields. Using the Boussinesq-Papkovitch potential representation for the mode I weight function field, it is shown that, for coplanar cracks in an infinite isotropic and homogeneous linear elastic body, the mean T-stress along an arbitrary crack front is independent of the shape and size of the crack. Further, a universal relation is discovered between the mean T-stress and the stress field at the same crack front location under the same loading but in the absence of a crack. First-order-accurate solutions are given for the T-stress variation along a slightly wavy crack front with nearly circular or straight confifurations. Specifically, cosine wave functions are adopted to describe smooth polygonal and slightly undulating planar crack shapes. The results indicate that T ₁₁, the 2D T-stress component acting normal to the crack front, increases with the curvature of the crack front as it bows out but T ₃₃, acting parallel to the crack front, decreases with the crack front curvature.