Green's functions for the stress intensity factor evolution in finite cracks in orthotropic materials

The elastodynamic response of an infinite orthotropic material with finite crack under concentrated loads is examined. Solution for the stress intensity factor history around the crack tips is found. Laplace and Fourier transforms are employed to solve the equations of motion leading to a Fredholm integral equation on the Laplace transform domain. The dynamic stress intensity factor history can be computed by numerical Laplace transform inversion of the solution of the Fredholm equation. Numerical values of the dynamic stress intensity factor history for some example materials are obtained. This solution can be used as a Green's function to solve dynamic problems involving fini te cracks.     

Dynamic fracture mechanics study of an electrically impermeable mode III crack in a transversely isotropic piezoelectric material under pure electric load

The dynamic response of an electrically impermeable Mode III crack in a transversely isotropic piezoelectric material under pure electric load is investigated by treating the electric loading process as a transient impact load, which may be more appropriate to mimic the real service environment of piezoelectric materials. The stress intensity factor, the mechanical energy release rate, and the total energy release rate are derived and expressed as a function of time for a given applied electric load. The theoretical results indicate that a purely electric load can fracture the piezoelectric material if the stress intensity factor or the mechanical energy release rate is used as a failure criterion.     

Dynamic stress intensity factors around three cracks in an infinite elastic plane subjected to time-harmonic stress waves

The time-harmonic problem for an infinite elastic plane weakened by three parallel cracks has been solved. In this problem, two cracks are situated symmetrically on either side of a central crack and incident stresses impinge perpendicular to the cracks. Using the Fourier transform technique, the boundary conditions are reduced to four simultaneous integral equations. To solve the equations, the differences of displacements inside the cracks are expanded in a series. The unknown coefficients in the series are solved by the Schmidt method. The dynamic stress intensity factors are calculated numerically for several crack configurations. This revised version was published online in August 2006 with corrections to the Cover Date.     

Arbitrarily oriented microcracks near the tip of an interface macrocrack in bonded dissimilar anisotropic materials

It is well known that microcracking in brittle materials results in a reduction of the stress intensity factor (SIF) and energy release rate (ERR). The reduced SIF or ERR represents crack tip shielding which is of significant interest to micromechanics and material science researchers. However, the effect of microcracking on the SIF and ERR is a complicated subject even for isotropic homogeneous materials, and becomes much more formidable in case of interface cracks in bonded dissimilar solids. To unravel the micromechanics of interface crack tip shielding in bonded dissimilar anisotropic solids, an interface crack interacting with arbitrarily oriented subinterface microcracks in bonded dissimilar anisotropic materials is studied. After deducing the fundamental solutions for a subinterface crack under concentrated normal and tangential tractions, the present interaction problem is reduced to a system of integral equations which is then solved numerically. A J‐integral analysis is then performed with special attention focused on the J₂‐integral in a local coordinate system attached to the microcracks. Theoretical and numerical results reassert the conservation law of the J‐integral derived for isotropic materials 1 , 2 also to be valid for bonded dissimilar anisotropic materials. It is further concluded that there is a wastage when the remote J‐integral transmits across the microcracking zone from infinity to the interface macrocrack tip. In order to highlight the influence of microstructure on the interfacial crack tip stress field, the crack tip SIF and ERR in several typical cases are presented. It is interesting to note that the Mode I SIF at the interface crack tip is quite different from the ERR in bonded dissimilar anisotropic materials.     

Stress intensity factor analysis for an interface crack between dissimilar isotropic materials under thermal stress

Thermal stresses, one of the main causes of interfacial failure between dissimilar materials, arise from different coefficients of linear thermal expansion. Two efficient numerical procedures in conjunction with the finite element method (FEM) for the stress intensity factor (SIF) analysis of interface cracks under thermal stresses are presented. The virtual crack extension method and the crack closure integral method are modified using the superposition method. The SIF analyses of some interface crack problems under mechanical and thermal loads are demonstrated. Very accurate mode separated SIFs are obtained using these methods.     

On the dynamic crack propagation in an interphase with spatially varying elastic properties under inplane loading

This article provides a comprehensive theoretical investigation on a finite crack with constant length (Yoffe type crack) propagating in an interfacial layer with spatially varying elastic properties under inplane loading. The analytical formulations are developed using Fourier transforms and solving the resulting singular integral equations in terms of the opening and sliding displacements of the crack. The dynamic stress intensity factors and energy release rate are analyzed to study the dynamic fracture property of this inherent mixed mode crack problem. Numerical examples are provided to show the effects of the material properties, the thickness of the interfacial layer, the crack position and speed upon the dynamic fracture behaviour, and the singularity transition between the current crack and the corresponding interfacial crack for thin interphase.     

Analysis of the driving force for a generally oriented crack in a functionally graded strip sandwiched between two homogeneous half planes

The driving forces for a generally oriented crack problem embedded in a Functionally Graded strip sandwiched between two half plane are analyzed using singular integral equations with Cauchy kernels, and integrated using Lobatto-Chebyshev collocation. Mixed-mode Stress Intensity Factors (SIF) and Strain Energy Release Rates (SERR) are calculated. The Stress Intensity Factors are compared for accuracy with previously published results. Parametric studies are conducted for various non-homogeneity ratios, crack lengths, crack orientation and thickness of the strip. It is shown that the SERR is more complete and should be used for crack propagation analysis.     

Computation of stress intensity factors of interface cracks based on interaction energy release rates and BEM sensitivity analysis

Stress intensity factors of bimaterial interface cracks are evaluated based on the interaction energy release rates. The interaction energy release rate is defined based on the energy release rates of a cracked body, corresponding to two independent loading conditions, actual field and an auxiliary field, and is related to the sensitivities of the potential energies for crack extensions. The potential energy of a cracked body is expressed with a domain integral, which is converted to a boundary integral expression by applying the divergence theorem. By differentiating this expression with the crack length, a boundary integral expression for the interaction energy release rate is obtained. The boundary integral representation for the interaction energy release rate involves the displacement, the traction, and their sensitivity coefficients with respect to the crack length. The boundary element sensitivity analyses are used to calculate these quantities accurately. A regularized boundary integral equation relating the boundary displacement and traction is differentiated with respect to an arbitrary shape parameter to derive the regularized boundary integral equation for the sensitivity coefficients of the boundary displacement and traction. The proposed approach is applied to several cracks in dissimilar media and the results are compared with those obtained by the conventional approach based on the extrapolation method. The analytical displacement and stress solutions for an interface crack between two infinite dissimilar media subjected to uniform stresses at infinity are used to give the auxiliary field, in which the values of the stress intensity factors are known. It is demonstrated that the present method can give accurate results for the stress intensity factors of various bimaterial interface cracks under coarse mesh discretizations.     

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.     

Weight functions for interface cracks in dissimilar anisotropic piezoelectric materials

Weight functions proposed for interface cracks in dissimilar isotropic materials (Gao, 1991; Chen and Hasebe, 1994) are extended to treat those in piezoelectric materials. The difficulties in separating the eight distinct complex arguments are overcome. The pseudo-orthogonal properties of the eigenfunction expansion form found in isotropic dissimilar cases(Chen and Hasebe, 1994) are proved to be valid in the present cases although the mathematical manipulations performed here seem much more complicated than those in isotropic dissimilar materials. Several path-independent integrals are obtained and all the coefficients in the eigenfunction expansion form, including the K _I, K _II, K _III and K _e, could be calculated by the weight functions introduced in this paper. It is concluded that the weight functions presented here provide a powerful tool to calculate the dominant parameters at the interface crack tip without any special treatment to the singular stress field of the near-tip region.     

A sandwich three-point bend specimen for testing mode-I interlaminar fracture toughness for fiber-reinforced composite materials

A sandwich three-point bend specimen has recently been proposed to test mode-I interlaminar fracture toughness for fiber-reinforced composite materials. The test composite consist of a thin layer bonded by two lateral reusable steel bars (Sohn et al. 1995). Some time earlier this specimen configuration was used to test fracture toughness of adhesives (Zdaniewsk et al. 1987). However, formulae for analysing its fracture mechanics parameters such as stress intensity factor and energy release rate can not be found in the literature. The lack of adequate formulae may explain why suitable quantitative analysis using this specimen configuration has not been achieved. In this paper, a simple and effective homogenisation method is used to change the bi-material system, which represents the specimen, into single uniform test material. This physical homogenisation is carried out by geometric change of the cross section of lateral steel parts based on equal deflection rigidity. For the transformed specimen configuration of single uniform material, the corresponding stress intensity factor solution from handbooks is available. Two formulae of stress intensity factor for the sandwich three-point bend specimen are given as upper limit and lower limit respectively, they are plotted with varying elastic tensile modulus mismatch. Then the relation between stress intensity factor and energy release rate, with special consideration of orthotropy of the tested composite material, is used to derive its energy release rate. The specimen and its formulae can also be applied to test other materials such as wood, welded joints (Burstow and Ainsworth, 1995), as well as to test dynamic fracture toughness. This revised version was published online in July 2006 with corrections to the Cover Date.     

Dynamic stress intensity factor due to concentrated loads on a propagating semi-infinite crack in orthotropic materials

The elastodynamic response of an infinite orthotropic material with a semi-infinite crack propagating at constant speed under the action of concentrated loads on the crack faces is examined. Solution for the stress intensity factor history around the crack tip is found for the loading modes I and II. Laplace and Fourier transforms along with the Wiener-Hopf technique are employed to solve the equations of motion. The asymptotic expression for the stress near the crack tip is analyzed which lead to a closed-form solution of the dynamic stress intensity factor. It is found that the stress intensity factor for the propagating crack is proportional to the stress intensity factor for a stationary crack by a factor similar to the universal function k(v) from the isotropic case. Results are presented for orthotropic materials as well as for the isotropic case.     

A continuum shape sensitivity method for fracture analysis of isotropic functionally graded materials

This paper presents a new continuum shape sensitivity method for calculating mixed-mode stress-intensity factors for a stationary crack in two-dimensional, linear- elastic, isotropic FGMs with arbitrary geometry. The method involves the material derivative concept taken from continuum mechanics, the mutual potential energy release rate, and direct differentiation. Since the governing variational equation is differentiated prior to discretization, resulting sensitivity equations are independent of approximate numerical techniques, such as the finite element method, boundary element method, mesh-free method, or others. The discrete form of the mutual potential energy release rate is simple and easy to calculate, as it only requires multiplication of displacement vectors and stiffness sensitivity matrices. By judiciously selecting the velocity field, the method only requires displacement response in a subdomain close to the crack tip, thus making the method computationally efficient. Seven finite-element based numerical examples, which comprise mode-I and mixed-mode deformations and/or single or multiple interacting cracks, are presented to evaluate the accuracy of the fracture parameters calculated by the proposed method. Comparisons have been made between stress-intensity factors predicted by the proposed method and available reference solutions in the literature, generated either analytically or numerically using various other fracture integrals or analyses. Excellent agreement is obtained between the results of the proposed method and previously obtained solutions. Therefore, shape sensitivity analysis provides an attractive alternative to fracture analysis of cracks in homogeneous and non-homogeneous materials.     

An interaction integral method for computing mixed mode stress intensity factors for curved bimaterial interface cracks in non-uniform temperature fields

In finite element analysis the interaction integral has been a useful tool for computing the stress intensity factors for fracture analysis. This work extends the interaction integral to account for non-uniform temperatures in the calculation of stress intensity factors for three dimensional curvilinear cracks either in a homogeneous body or on a bimaterial interface. First, the derivation of the computational algorithm, which includes the additional terms developed by the non-zero gradient of the temperature field, is presented in detail. The algorithm is then implemented in conjunction with commercial finite element software to calculate the stress intensity factors of a crack undergoing non-uniform temperatures on both a homogeneous and a bimaterial interface. The numerical results displayed path independence and showed excellent agreement with available analytical solutions.     

Mode I stress intensity factors for curved cracks in gears by a weight functions method

The most recent trend in power transmission design considers the damage-tolerant approach as one of the methods to obtain safe, reliable and light systems. This means that components containing cracks must be considered and analysed to understand the conditions that cause critical cracks and defects and their dimensions.
In this paper a cracked tooth of an automotive gearwheel is considered. A numerical procedure (based on the slice synthesis weight function method) to calculate the stress intensity factors of curved cracks due to bending loads is illustrated. The results are compared with those obtained by expensive finite element calculations. The agreement is satisfactory and the proposed technique proves to be very attractive from the point of view of time saving.
One example of an application to fatigue design practice is provided, namely the analysis of fatigue crack propagation in surface-treated gears. The results show the role played by residual stresses induced by carburizing and shot peening.     

Application of finite element method for determination of stress intensity factor and plastic zone geometry in an aluminium alloy sheet under uniaxial tension

High strength materials have gained prominence in the fields of aero-structures, space missiles, ship-building, pressure vessels etc. However, high strength materials are often characterised by low values of crack resistance or fracture toughness. Knowledge of stress intensity factor (SIF) is essential to predict their fracture toughness. SIF values can be obtained both theoretically and experimentally. Theoretical methods include analytical techniques as well as the finite element method (FEM). The former is used for simpler geometries and the latter for complicated geometries of engineering structures. The SIF as a function of crack size in an aluminium alloy 2024-T₃ (Al-4·5% Cu, 1·5% Mg, 0·6% Mn) sheet was determined by a computer method. These values were obtained directly from the stresses as well as indirectly from strain energy release rateG andJ integral. The results agree well with the normalised values obtained from an ASTM formula. The size and shape of the plastic zone at the crack tip have been determined as a function of nominal stress for a fixed crack length. The plastic zone has the form of two ellipsoids with their maximum spreads oriented around 69° to the crack axis.     

Determination of stress intensity factors for bimaterial interface stationary rigid line inclusions by boundary element method

The stress intensity factors for a rigid line inclusion lying along a bimaterial interface are calculated by the boundary element method with the multiregion and the discontinuous traction singular elements. The relationships between the stress intensity factors and the inclusion surface stresses are derived. The numerically computed stress intensity factors for the bimaterial interface rigid line inclusion in the infinite body are proved to be in good agreement within 3% when compared with the previous exact solutions. In the finite bimaterial models, the stress intensity factors for the center and edge rigid line inclusions at the interface are computed with the variation of the rigid line inclusion length and the shear modulus ratio under the uniaxial and biaxial loading conditions.     

Experimental and finite element analyses on stress intensity factors of an elliptical surface crack in a circular shaft under tension and bending

Experimental backtracking technique and finite element analysis have been employed to evaluate the stress intensities along the front of an elliptical surface crack in a cylindrical rod. The finite element solution covers a wide range of crack shapes loaded under end-free and end-constrained axial tension and pure bending. Convenient closed form stress intensity expressions along the whole crack front for each of the loading cases have been given in terms of the crack aspect ratio, crack depth ratio and place ratio.The closed form solutions have been compared against a number of representative solutions collected from the literature. It has been found that different finite element results for the interior points are generally in good mutual agreement, while solutions derived from other methods may sometimes indicate different trends. At the surface interception point agreement is less good because of a complication in the interpretation of stress intensity there.Experimental backtracking results on the end-constrained axial tension case corroborate well with the closed form solution presented. It suggests that the current closed form solution is adequate in describing the stress intensities along the whole crack front of real surface cracks in cylindrical rods.     

A direct method for calculating the stress intensity factor in BEM

The proposed algorithm employs singular crack tip elements in which the stress intensity factor appears as a degree of freedom. The additional degrees of freedom are compensated by constraint conditions which originate from imposing continuity across elements and a contour integration formula. The two benchmark problems indicate the proposed algorithm can accurately predict the stress intensity factor and the distribution of the primary and secondary variables in fracture problems.