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.     

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.     

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.     

Dynamic stress intensity factors of two 3D rectangular permeable cracks in a transversely isotropic piezoelectric material under a time‐harmonic elastic P‐wave

The dynamic behavior of two 3D rectangular permeable cracks in a transversely isotropic piezoelectric material is investigated under an incident harmonic stress wave by using the generalized Almansi's theorem and the Schmidt method. The problem is formulated through double Fourier transform into three pairs of dual integral equations with the displacement jumps across the crack surfaces as the unknown variables. To solve the dual integral equations, the displacement jumps across the crack surfaces are directly expanded as a series of Jacobi polynomials. Finally, the relations among the dynamic stress field and the dynamic electric displacement filed near the crack edges are obtained, and the effects of the shape of the rectangular crack, the characteristics of the harmonic wave, and the distance between two rectangular cracks on the stress and the electric intensity factors in a piezoelectric composite material are analyzed. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

Reduced basis approximation and a posteriori error estimation for stress intensity factors

We develop reduced basis approximations, rigorous a posteriori error bounds, and offline–online computational procedures for the accurate, fast, and reliable prediction of stress intensity factors or strain energy release rates for ‘Mode I’ linear elastic crack problems relevant to fracture mechanics. We demonstrate the efficiency and rigor of our numerical method in several examples. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

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.     

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.     

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.     

聚合物梯度材料黏弹性断裂的双控制参数

针对组分材料体积分数任意分布的聚合物功能梯度材料,研究其在蠕变加载条件下Ⅰ型裂纹应力强度因子（SIFs）和应变能释放率的时间相依特征。由Mori-Tanaka方法预测等效松弛模量,在Laplace变换域中采用梯度有限元法和虚拟裂纹闭合方法计算断裂参数,由数值逆变换得到物理空间的对应量。分析边裂纹平行于梯度方向的聚合物功能梯度板条,分别考虑均匀拉伸和三点弯曲蠕变加载。结果表明,聚合物梯度材料应变能释放率随时间增加,其增大的程度与黏弹性组分材料体积分数正相关;材料的非均匀黏弹性性质产生应力重新分布,导致裂纹尖端应力场强度随时间变化,当裂纹位于黏弹性材料含量较低的一边时,应力强度因子随时间增加,反之,随时间减小。而且,材料的应力强度因子与时间相依的变化范围和体积分数分布以及加载方式有关,当体积分数接近线性分布时,变化最明显,三点弯曲比均匀拉伸的变化大。SIFs随时间的延长增加或减小、加剧或减轻裂纹尖端部位的“衰坏”,表明黏弹性功能梯度裂纹体的延迟失稳需要联合采用应力强度因子与应变能释放率作为双控制参数。