Boundary element‐free method for fracture analysis of 2‐D anisotropic piezoelectric solids

This paper considers a 2‐D fracture analysis of anisotropic piezoelectric solids by a boundary element‐free method. A traction boundary integral equation (BIE) that only involves the singular terms of order 1/r is first derived using integration by parts. New variables, namely, the tangential derivative of the extended displacement (the extended displacement density) for the general boundary and the tangential derivative of the extended crack opening displacement (the extended displacement dislocation density), are introduced to the equation so that solution to curved crack problems is possible. This resulted equation can be directly applied to general boundary and crack surface, and no separate treatments are necessary for the upper and lower surfaces of the crack. The extended displacement dislocation densities on the crack surface are expressed as the product of the characteristic terms and unknown weight functions, and the unknown weight functions are modelled using the moving least‐squares (MLS) approximation. The numerical scheme of the boundary element‐free method is established, and an effective numerical procedure is adopted to evaluate the singular integrals. The extended ‘stress intensity factors’ (SIFs) are computed for some selected example problems that contain straight or curved cracks, and good numerical results are obtained. Copyright © 2006 John Wiley & Sons, Ltd.     

A general Boundary Element Analysis of 2-D Linear Elastic Fracture Mechanics

This paper presents a boundary element method (BEM) analysis of linear elastic fracture mechanics in two-dimensional solids. The most outstanding feature of this new analysis is that it is a single-domain method, and yet it is very accurate, efficient and versatile: Material properties in the medium can be anisotropic as well as isotropic. Problem domain can be finite, infinite or semi-infinite. Cracks can be of multiple, branched, internal or edged type with a straight or curved shape. Loading can be of in-plane or anti-plane, and can be applied along the no-crack boundary or crack surface. Furthermore, the body-force case can also be analyzed. The present BEM analysis is an extension of the work by Pan and Amadei (1996a) and is such that the displacement and traction integral equations are collocated, respectively, on the no-crack boundary and on one side of the crack surface. Since in this formulation the displacement and/or traction are used as unknowns on the no-crack boundary and the relative crack displacement (i.e. displacement discontinuity) as unknown on the crack surface, it possesses the advantages of both the traditional displacement BEM and the displacement discontinuity method (DDM) and yet gets rid of the disadvantages associated with these methods when modeling fracture mechanics problems. Numerical examples of calculation of stress intensity factors (SIFs) for various benchmark problems were conducted and excellent agreement with previously published results was obtained. This revised version was published online in August 2006 with corrections to the Cover Date.     

Boundary element analysis of fracture mechanics in anisotropic bimaterials

This paper presents a boundary element formulation for the analysis of linear elastic fracture mechanics problems involving anisotropic bimaterials. The most important feature associated with the present formulation is that it is a single domain method, and yet it is accurate, efficient and versatile. In this formulation, the displacement integral equation is collocated on the uncracked boundary only, and the traction integral equation is collocated on one side of the crack surface only. The complete Green's functions for anisotropic bimaterials are also derived and implemented into the boundary integral formulation so that discretization along the interface can be avoided except for the interfacial crack part. A special crack-tip element is introduced to capture exactly the crack-tip behavior.Numerical examples are presented for the calculations of stress intensity factors for a straight crack with various locations in infinite bimaterials. It is found that very accurate results can be obtained by the proposed method even with relatively coarse discretization. Numerical results also show that material anisotropy can greatly affect the stress intensity factor.     

Fracture mechanics analysis of cracked 2-D anisotropic media with a new formulation of the boundary element method

A new formulation of the boundary element method (BEM) is proposed in this paper to calculate stress intensity factors for cracked 2-D anisotropic materials. The most outstanding feature of this new approach is that the displacement and traction integral equations are collocated on the outside boundary of the problem (no-crack boundary) only and on one side of the crack surfaces only, respectively. Since the new BEM formulation uses displacements or tractions as unknowns on the outside boundary and displacement differences as unknowns on the crack surfaces, the formulation combines the best attributes of the traditional displacement BEM as well as the displacement discontinuity method (DDM). Compared with the recently proposed dual BEM, the present approach doesn't require dua elements and nodes on the crack surfaces, and further, it can be used for anisotropic media with cracks of any geometric shapes. Numerical examples of calculation of stress intensity factors were conducted, and excellent agreement with previously published results was obtained. The authors believe that the new BEM formulation presented in this paper will provide an alternative and yet efficient numerical technique for the study of cracked 2-D anisotropic media, and for the simulation of quasi-static crack propagation.     

Cracks emanating from a square hole in rectangular plate in tension

A numerical analysis of cracks emanating from a square hole in a rectangular plate in tension is performed using a hybrid displacement discontinuity method (a boundary element method). Detailed solutions of the stress intensity factors (SIFs) of the plane elastic crack problem are given, which can reveal the effect of geometric parameters of the cracked body on the SIFs. By comparing the calculated SIFs of the plane elastic crack problem with those of the centre crack in a rectangular plate in tension, in addition, an amplifying effect of the square hole on the SIFs is found. The numerical results reported here also prove that the boundary element method is simple, yet accurate, for calculating the SIFs of complex crack problems in finite plate.     

Analysis of elliptical cracks perpendicular to the interface of two joined transversely isotropic solids

This paper presents a boundary element analysis of elliptical cracks in two joined transversely isotropic solids. The boundary element method is developed by incorporating the fundamental singular solution for a concentrated point load in a transversely isotropic bi-material solid of infinite space into the conventional displacement boundary integral equations. The multi-region method is used to analyze the crack problems. The traction-singular elements are employed to capture the singularity around the crack front. The values of stress intensity factors (SIFs) are obtained by using crack opening displacements. The results of the proposed method compare well with the existing exact solutions for an elliptical crack parallel to the isotropic plane of a transversely isotropic solid of infinite extent. Elliptical cracks perpendicular to the interface of transversely isotropic bi-material solids of either infinite extent or occupying a cubic region are further examined in detail. The crack surfaces are subject to the uniform normal tractions. The stress intensity factor values of the elliptical cracks of the two types are analyzed and compared. Numerical results have shown that the stress intensity factors are strongly affected by the anisotropy and the combination of the two joined solids.     

A unified and integrated numerical‐analytical approach for evaluation of Jk‐integrals of an interfacial crack in orthotropic and isotropic bimaterials

A new unified and integrated method for the numerical‐analytical calculation of J_k‐integrals of an in‐plane traction free interfacial crack in homogeneous orthotropic and isotropic bimaterials is presented. The numerical algorithm, based on the boundary element crack shape sensitivities, is generic and flexible. It applies to both straight and curved interfacial cracks in anisotropic and/or isotropic bimaterials. The shape functions of semidiscontinuous quadratic quarter point crack tip elements are correctly scaled to adapt the singular oscillatory near tip field of tractions. The length of crack is designated as the design variable to compute the strain energy release rate precisely. Although an analytical equation relating J₁ and stress intensity factors (SIFs) exists, a similar relation for J₂ in debonded anisotropic solids for decoupling SIFs is not available. An analytical expression was recently derived by this author for J₂ in aligned orthotropic/orthotropic bimaterials with a straight interface crack. Using this new relation and the present computed J_k values, the SIFs can be decoupled without the need for an auxiliary equation. Here, the aforementioned analytical relation is reconstructed for cubic symmetry/isotropic bimaterials and used with the present numerical algorithm. An example with known analytical SIFs is presented. The numerical and analytical magnitudes of J_k for an interface crack in orthotropic/orthotropic and cubic symmetry/isotropic bimaterials show an excellent agreement.     

Cracks emanating from a rhombus hole in infinite and finite plate subjected to internal pressure

This note deals with the stress intensity factors (SIFs) of cracks emanating from a rhombus hole in a rectangular plate subjected to internal pressure by means of the displacement discontinuity method with crack-tip elements (a boundary element method) proposed recently by the author. Moreover, an empirical formula of the SIFs of the crack problem is presented and examined. It is found that the empirical formula is very accurate for evaluating the SIFs of the crack problem.     

Cracks emanating from circular hole or square hole in rectangular plate in tension

A numerical analysis of cracks emanating from a circular hole (Fig. 1) or a square hole (Fig. 2) in rectangular plate in tension is performed by means of the displacement discontinuity method with crack-tip elements (a boundary element method) presented recently by the author. Detail solutions of the stress intensity factors (SIFs) of the two plane elastic crack problems are given, which can reveal the effect of geometric parameters of the cracked bodies on the SIFs. By comparing the SIFs of the two crack problems with those of the center crack in rectangular plate in tension (Fig. 3), in addition, an effect of the circular hole or the square hole on the SIFs of the center crack is discussed in detail. The numerical results reported here also illustrate that the boundary element method is simple, yet accurate for calculating the SIFs of complex crack problems in finite plate.     

Boundary element formulation for 3D transversely isotropic cracked bodies

The boundary traction integral representation is obtained in elasticity when the classical displacement representation is differentiated and combined according to Hooke's law. The use of both traction and displacement integral representations leads to a mixed (or dual) formulation of the BEM where the discretization effort for crack problems is much smaller than in the classical formulation. A boundary element analysis of three‐dimensional fracture mechanics problems of transversely isotropic solids based on the mixed formulation is presented in this paper. The hypersingular and strongly singular kernels appearing in the formulation are regularized by using two terms of the displacement series expansion and one term of the traction expansion, at the collocation point. All the remaining integrals are analytically evaluated or transformed by means of Stokes' theorem into regular or weakly singular integrals, which are numerically computed. The method is general and can be used for elements of any shape including quarter‐point crack front elements. No change of co‐ordinates is required for the integration. The formulation as presented in this paper is something as clear, general and easy to handle as the classical BE formulation. It is used in combination with three‐dimensional quadratic and quarter‐point elements to obtain accurate results for several different crack problems. Cracks in boundless and finite transversely isotropic domains are studied. The meshes are simple and include only discretization of the crack and the external boundary. The obtained results are in good agreement with those existing in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

A BEM analysis of fracture mechanics in 2D anisotropic piezoelectric solids

This paper presents a single-domain boundary element method (BEM) analysis of fracture mechanics in 2D anisotropic piezoelectric solids. In this analysis, the extended displacement (elastic displacement and electrical potential) and extended traction (elastic traction and electrical displacement) integral equations are collocated on the outside boundary (no-crack boundary) of the problem and on one side of the crack surface, respectively. The Green's functions for the anisotropic piezoelectric solids in an infinite plane, a half plane, and two joined dissimilar half-planes are also derived using the complex variable function method. The extrapolation of the extended relative crack displacement is employed to calculate the extended `stress intensity factors' (SIFs), i.e., K_I, K_II, K_III and K_IV. For a finite crack in an infinite anisotropic piezoelectric solid, the extended SIFs obtained with the current numerical formulation were found to be very close to the exact solutions. For a central and inclined crack in a finite and anisotropic piezoelectric solid, we found that both the coupled and uncoupled (i.e., the piezoelectric coefficient e_ijk=0) cases predict very similar stress intensity factors K_I and K_II when a uniform tension σ_yy is applied, and very similar electric displacement intensity factor K_IV when a uniform electrical displacement D_y is applied. However, the relative crack displacement and electrical potential along the crack surface are quite different for the coupled and uncoupled cases. Furthermore, for a inclined crack within a finite domain, we found that while a uniform σ_yy (=1 N m⁻²) induces only a very small electrical displacement intensity factor (in the unit of Cm^−3/2), a uniform D_y (=1 C m⁻²) can produce very large stress intensity factors (in the unit of Nm^−3/2).     

An effective method of stress intensity factor calculation for cracks emanating from a triangular or square hole under biaxial loads

This paper is concerned with stress intensity factors for cracks emanating from a triangular or square hole under biaxial loads by means of a new boundary element method. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfied and the 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 the corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. The method is called a Hybrid Displacement Discontinuity Method (HDDM). Numerical examples are included to show that the method is very efficient and accurate for calculating stress intensity factors for plane elastic crack problems. In addition, the present numerical results can reveal the effect of the biaxial loads on stress intensity factors.     

A hypersingular time‐domain BEM for 2D dynamic crack analysis in anisotropic solids

A hypersingular time‐domain boundary element method (BEM) for transient elastodynamic crack analysis in two‐dimensional (2D), homogeneous, anisotropic, and linear elastic solids is presented in this paper. Stationary cracks in both infinite and finite anisotropic solids under impact loading are investigated. On the external boundary of the cracked solid the classical displacement boundary integral equations (BIEs) are used, while the hypersingular traction BIEs are applied to the crack‐faces. The temporal discretization is performed by a collocation method, while a Galerkin method is implemented for the spatial discretization. Both temporal and spatial integrations are carried out analytically. Special analytical techniques are developed to directly compute strongly singular and hypersingular integrals. Only the line integrals over an unit circle arising in the elastodynamic fundamental solutions need to be computed numerically by standard Gaussian quadrature. An explicit time‐stepping scheme is obtained to compute the unknown boundary data including the crack‐opening‐displacements (CODs). Special crack‐tip elements are adopted to ensure a direct and an accurate computation of the elastodynamic stress intensity factors from the CODs. Several numerical examples are given to show the accuracy and the efficiency of the present hypersingular time‐domain BEM. Copyright © 2008 John Wiley & Sons, Ltd.     

An efficient dual boundary element technique for a two-dimensional fracture problem with multiple cracks

An efficient dual boundary element technique for the analysis of a two-dimensional finite body with multiple cracks is established. In addition to the displacement integral equation derived for the outer boundary, since the relative displacement of the crack surfaces is adopted in the formulation, only the traction integral equation is established on one of the crack surfaces. For each crack, a virtual boundary is devised and connected to one of the crack surfaces to construct a closed integral path. The rigid body translation for the domain enclosed by the closed integral path is then employed for evaluating the hypersingular integral. To solve the dual displacement/traction integral equations simultaneously, the constant and quadratic isoparametric elements are taken to discretize the closed integral paths/crack surfaces and the outer boundary, respectively. The present method has distinct computational advantages in solving a fracture problem which has arbitrary numbers, distributions, orientations and shapes of cracks by a few boundary elements. Several examples are analysed and the computed results are in excellent agreement with other analytical or numerical solutions.     

Dynamic crack interactions in magnetoelectroelastic composite materials

In this paper, the dynamic interactions among cracks embedded in a two-dimensional (2-D) piezoelectric-piezomagnetic composite material are analyzed by means of a hypersingular formulation of the boundary element method. In the numerical solution procedure, extended crack opening displacements and extended traction jumps across the crack are considered as basic unknowns, so that only the traction boundary integral equations are needed on the crack surfaces. Quadratic discontinuous boundary elements are implemented together with discontinuous quarter-point elements placed next to the crack tips to ensure a proper representation of the square root asymptotic behavior. Several impermeable cracks configurations subjected to time-harmonic incident L-waves are analyzed in order to characterize the effects of the magnetoelectromechanical coupling on the dynamic crack interactions and to illustrate the dependence on such coupling of the fracture parameters: stress intensity factors, electric displacement intensity factor and magnetic induction intensity factor.     

Rectangular tensile sheet with single edge crack or edge half-circular-hole crack

By using the displacement discontinuity method with crack-tip elements (a boundary element method) proposed recently by the author, this note presents the stress intensity factors (SIFs) of a rectangular tensile plate with single edge crack. Further this note studies the SIFs of crack emanating from an edge half-circular hole. By comparing the calculated SIFs of the single edge half-circular-hole crack with those of the single edge crack, a shielding effect of the half-circular hole on the SIFs of the single edge crack is discussed. It is found that the boundary element method is simple, yet accurate for calculating the SIFs of complex crack problems in finite plate.     

A new boundary integral equation method for cracked 2-D anisotropic bodies

By using integration by parts to the traditional boundary integral formulation, a traction boundary integral equation for cracked 2-D anisotropic bodies is derived. The new traction integral equation involves only singularity of order 1/r and no hypersingular term appears. The dislocation densities on the crack surface are introduced and the relations between stress intensity factors and dislocation densities near the crack tip are induced to calculate the stress intensity factors. The boundary element method based on the new equation is established and the singular interpolation functions are introduced to model the singularity of the dislocation density (in the order of ) for crack tip elements. The proposed method can be directly used for the 2-D anisotropic body containing cracks of arbitrary geometric shapes. Several numerical examples demonstrate the validity and accuracy of BEM based on the new boundary integral equation.     

Rectangular Tensile Sheet with Double Edge Notch Cracks

This paper deals with the rectangular tensile sheet with symmetric double edge notch cracks. Such a crack problem is called an edge notch crack problem for short. By using a hybrid displacement discontinuity method (a boundary element method), two edge notch models are analyzed in detail. By changing the geometrical forms and parameters of the edge notch, and by comparing the stress intensity factors (SIFs) of the edge notch crack problem with those of the double edge cracked plate tension specimen (DECT), which is a model frequently used in fracture mechanics, the effect of the geometrical forms and parameters of the edge notch on the SIFs of the DECT specimen, is revealed. Some geometric characterestic parameters are introduced here, which are used to formulate the notch length and the branch crack length, which are to be determined in mechanical machining of the DECT specimen So we can say that the geometric characterestic parameters and the formulae used to determine the notch length and the branch crack length presented in this paper perhaps have some guidance role for mechanical machining of the DECT specimen.     

Analyzing interaction between coplanar square cracks using an efficient boundary element‐free method

This paper examines the interaction between coplanar square cracks by combining the moving least‐squares (MLS) approximation and the derived boundary integral equation (BIE). A new traction BIE involving only the Cauchy singular kernels is derived by applying integration by parts to the traditional boundary integral formulation. The new traction BIE can be directly applied to a crack surface and no displacement BIE is necessary because all crack boundary conditions (both upper and lower ones) are incorporated. A boundary element‐free method is then developed by combining the derived BIE and MLS approximation, in which the crack opening displacement is first expressed as the product of weight functions and the characteristic terms, and the unknown weight is approximated with the MLS approximation. The efficiency of the developed method is tested for isotropic and transversely isotropic media. The interaction between two and three coplanar square cracks in isotropic elastic body is numerically studied and the case of any number of coplanar square cracks is deduced and discussed. Copyright © 2012 John Wiley & Sons, Ltd.     

Three-dimensional analysis of thermally loaded cracks

The stress intensity factors for cracks in three-dimensional, thermally stressed structures are computed by using the boundary element method. While many boundary and volume-integral-based formulations are available for the treatment of thermoelastic problems in solids, the present analysis is based on a recently developed boundary-only formulation. The accuracy of the solutions in the present work is improved by using special elements at the crack front that accurately model the variation of displacement, temperature fields and singularity of traction, flux fields.