A new boundary element method for mixed boundary value problems involving cracks and holes: Interactions between rigid inclusions and cracks

This paper derives a new boundary integral equation (BIE) formulation for plane elastic bodies containing cracks and holes and subjected to mixed displacement/ traction boundary conditions, and proposes a new boundary element method (BEM) based upon this formulation. The basic unknown in the formulation is a complex boundary function H(t), which is a linear combination of the boundary traction and boundary displacement density. The present BIE formulation can be related directly to Muskhelishvili's formalism. Singular interpolation functions of order r ^–1/2 (where r is the distance measured from the crack tip) are introduced such that singular integrand involved at the element level can be integrated analytically. By applying the BEM, the interaction between a rigid circular inclusion and a crack is investigated in details. Our results for the stress intensity factor are comparable with those given by Erdogan and Gupta (1975) and Gharpuray et al. (1990) for a crack emanating from a stiff inclusion, and with those by Erdogan et al. (1974) for a crack in the neighborhood of a stiff inclusion.     

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.     

BEM for crack-hole problems in thermopiezoelectric materials

The boundary element formulation for analysing interaction between a hole and multiple cracks in piezoelectric materials is presented. Using Green's function for hole problems and variational principle, a boundary element model (BEM) for a 2-D thermopiezoelectric solid with cracks and holes has been developed and used to calculate stress intensity factors of the crack-hole problem. In BEM, the boundary condition on the hole circumference is satisfied a priori by Green's function, and is not involved in the boundary element equations. The method is applicable to multiple-crack problems in both finite and infinite solids. Numerical results for stress and electric displacement intensity factors at a particular crack tip in a crack-hole system of piezoelectric materials are presented to illustrate the application of the proposed formulation.     

Stress analyses around holes in composite laminates using boundary element method

A three-dimensional (3D) boundary element method (BEM) is developed for the analysis of composite laminates with holes. Instead of using Kelvin-type Green's functions of anisotropic infinite space, 3D layered Green's functions with the materials of each layer being generally anisotropic, derived recently in the Fourier transform domain, are implemented into a 3D BEM formulation. A novel numerical algorithm is designed to calculate layered Green's functions efficiently. It should be noted that since layered Green's functions satisfy exactly the continuity conditions along the interfaces and top and bottom free surfaces a priori, the model becomes truly 2D and discretization is only needed along the hole surface and prescribed traction and/or displacement boundaries. To test the validity and accuracy of the proposed method, the present layered BEM formulation is applied to the problem of an infinite anisotropic plate with a circular hole where the analytical solution is available. It is found that even with a very coarse mesh, the present BEM can predict the hoop stress very accurately along the hole surface. The BEM formulation is then applied to analyze two composite laminates (90/0)_s and (−45/45)_s, under a remote in-plane strain, that have been studied previously with different approaches. For the (90/0)_s case, the hoop stresses along the hole surface predicted by the present layered BEM formulation are in very close agreement with the previous results. For the (−45/45)_s case, however, it is found that a nearly converged solution (less than 5% convergence by doubling the mesh) by the present method is at significant variance with the previous ones that are lack-of-convergence checks. It can be expected that for designing the bolted joints of composites with many layers, a computational tool developed based on the present techniques would be robust and offer a much better solution with regard to accuracy, versatility and design cycle time.     

Material properties of piezoelectric composites by BEM and homogenization method

Applications of boundary element method (BEM) to piezoelectric composites in conjunction with homogenization approach for determining their effective material properties are discussed in this paper. The composites considered here consist of inclusion and matrix phases. The homogenization model for composites with inhomogeneities is developed and introduced into a BE formulation to provide an effective means for estimating overall material constants of two-phase composites. In this model, a representative volume element (RVE) is used whose volume average stress and strain are calculated by the boundary tractions and displacements of the RVE. Thus BEM is suitable for performing calculations on average stress and strain fields of the composites. Numerical results for a piezoelectric plate with circular inclusions are presented to illustrate the application of the proposed micromechanics––BE formulation.     

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).     

Fatigue crack growth through particulate clusters in polycarbonate material

The interaction of a crack with a perfectly bonded inclusion or a cluster of inclusions in polycarbonate matrix was investigated through both numerical simulations and fatigue tests. Stress intensity factors (K_I) were evaluated by boundary element method for several particle sizes, position and finally for inclusion cluster as a precursor study for the experiments. The numerical simulation has shown the crack tendency to circumvent the inclusions with consequential reduction of the growth rate. Fatigue crack growth tests were carried out on several particle-filled specimens at constant value of the applied stress intensity factor range (ΔK_Iapp) highlighting the crack delay due to the presence of the stiff second phase. The experiments demonstrated that the inclusion effect on the crack growth rate can be explained with a model based on the crack shielding effect in which the particle would act to reduce the effective stress intensity factor at crack tip (K_Ieff). Finally, the crack growth rate was predicted with an analytical model, and then compared to that obtained by the fatigue testing. Possible explanations for differences are discussed.     

Boundary element analysis for an interface crack between dissimilar elastoplastic materials

The boundary element method (BEM) is presented for elastoplastic analysis of cracks between two dissimilar materials. The boundary integral equations and integral representation of stress rates are written in such a form that all integrals can be evaluated by the regular Gaussian quadrature rule. An advanced multidomain BEM formulation is suggested for the solution of analysed problems where the substantial reduction of stiffness matrix is observed. The elastoplastic behaviour is modelled through the use of an approximation for the plastic component of the stresses. The boundary and the yielding zone are discretized by elements with quadratic approximations. In numerical examples the path independence of the J- and L-integrals for a straight interface crack and a circular arc-shaped interface crack are investigated, respectively. The influence of the different values of Young's modulus on the J-integral, shape and size of plastic zones is treated too.     

A numerical Green's function and dual reciprocity BEM method to solve elastodynamic crack problems

A computational model based on the numerical Green's function (NGF) and the dual reciprocity boundary element method (DR-BEM) is presented for the study of elastodynamic fracture mechanics problems. The numerical Green's function, corresponding to an embedded crack within the infinite medium, is introduced into a boundary element formulation, as the fundamental solution, to calculate the unknown external boundary displacements and tractions and in post-processing determine the crack opening displacements (COD). The domain inertial integral present in the elastodynamic equation is transformed into a boundary integral one by the use of the dual reciprocity technique. The dynamic stress intensity factors (SIF), computed through crack opening displacement values, are obtained for several numerical examples, indicating a good agreement with existing solutions.     

On the conventional boundary integral equation formulation for piezoelectric solids with defects or of thin shapes

In this paper, the conventional boundary integral equation (BIE) formulation for piezoelectric solids is revisited and the related issues are examined. The key relations employed in deriving the piezoelectric BIE, such as the generalized Green's identity (reciprocal work theorem) and integral identities for the piezoelectric fundamental solution, are established rigorously. A weakly singular form of the piezoelectric BIE is derived for the first time using the identities for the fundamental solution, which eliminates the calculation of any singular integrals in the piezoelectric boundary element method (BEM). The crucial question of whether or not the piezoelectric BIE will degenerate when applied to crack and thin shell-like problems is addressed. It is shown analytically that the conventional BIE for piezoelectricity does degenerate for crack problems, but does not degenerate for thin piezoelectric shells. The latter has significant implications in applications of the piezoelectric BIE to the analysis of thin piezoelectric films used widely as sensors and actuators. Numerical tests to show the degeneracy of the piezoelectric BIE for crack problems are presented and one remedy to this degeneracy by using the multi-domain BEM is also demonstrated.     

The 3‐D BEM implementation of a numerical Green's function for fracture mechanics applications

The use of Green's functions has been considered a powerful technique in the solution of fracture mechanics problems by the boundary element method (BEM). Closed‐form expressions for Green's function components, however, have only been available for few simple 2‐D crack geometry applications and require complex variable theory. The present authors have recently introduced an alternative numerical procedure to compute the Green's function components that produced BEM results for 2‐D general geometry multiple crack problems, including static and dynamic applications. This technique is not restricted to 2‐D problems and the computational aspects of the 3‐D implementation of the numerical Green's function approach are now discussed, including examples. Copyright © 2000 John Wiley & Sons, Ltd.     

2D Green's functions of defective magnetoelectroelastic solids under thermal loading

Thermomagnetoelectroelastic problems for various defects embedded in an infinite matrix are considered in this paper. Using Stroh's formalism, conformal mapping, and perturbation technique, Green's functions are obtained in closed form for a defect in an infinite magnetoelectroelastic solid induced by the thermal analog of a line temperature discontinuity and a line heat source. The defect may be of an elliptic hole or a Griffith crack, a half-plane boundary, a bimaterial interface, or a rigid inclusion. These Green's functions satisfy the relevant boundary or interface conditions. The proposed Green's functions can be used to establish boundary element formulation and to analyzing fracture behaviour due to the defects mentioned above.     

The boundary contour method for piezoelectric media with linear boundary elements

This paper presents a development of the boundary contour method (BCM) for piezoelectric media. First, the divergence‐free property of the integrand of the piezoelectric boundary element is proved. Secondly, the boundary contour method formulation is derived and potential functions are obtained by introducing linear shape functions and Green's functions (Computer Methods in Applied Mechanics and Engineering 1998; 158 : 65) for piezoelectric media. The BCM is applied to the problem of piezoelectric media. Finally, numerical solutions for illustrative examples are compared with exact ones and those of the conventional boundary element method (BEM). The numerical results of the BCM coincide very well with the exact solution, and the feasibility and efficiency of the method are verified. Copyright © 2003 John Wiley & Sons, Ltd.     

On a regularized method of fundamental solutions coupled with the numerical Green's function procedure to solve embedded crack problems

The method of fundamental solutions (MFS) is applied to solve linear elastic fracture mechanics (LEFM) problems. The approximate solution is obtained by means of a linear combination of fundamental solutions containing the same crack geometry as the actual problem. In this way, the fundamental solution is the very same one applied in the numerical Green's function (NGF) BEM approach, in which the singular behavior of embedded crack problems is incorporated. Due to severe ill-conditioning present in the MFS matrices generated with the numerical Green's function, a regularization procedure (Tikhonov's) was needed to improve accuracy, stabilization of the solution and to reduce sensibility with respect to source point locations. As a result, accurate stress intensity factors can be obtained by a superposition of the generalized fundamental crack openings. This mesh-free technique presents good results when compared with the boundary element method and estimated solutions for the stress intensity factor calculations.     

Wave propagation in a simplified modelled poroelastic continuum: fundamental solutions and a time domain boundary element formulation

In finite element formulations for poroelastic continua a representation of Biot's theory using the unknowns solid displacement and pore pressure is preferred. Such a formulation is possible either for quasi‐static problems or for dynamic problems if the inertia effects of the fluid are neglected. Contrary to these formulations a boundary element method (BEM) for the general case of Biot's theory in time domain has been published (Wave Propagation in Viscoelastic and Poroelastic Continua: A Boundary Element Approach. Lecture Notes in Applied Mechanics. Springer: Berlin, Heidelberg, New York, 2001.). If the advantages of both methods are required it is common practice to couple both methods. However, for such a coupled FE/BE procedure a BEM for the simplified dynamic Biot theory as used in FEM must be developed. Therefore, here, the fundamental solutions as well as a BE time stepping procedure is presented for the simplified dynamic theory where the inertia effects of the fluid are neglected. Further, a semi‐analytical one‐dimensional solution is presented to check the proposed BE formulation. Finally, wave propagation problems are studied using either the complete Biot theory as well as the simplified theory. These examples show that no significant differences occur for the selected material. Copyright © 2005 John Wiley & Sons, Ltd.     

Fracture mechanics analysis of anisotropic plates by the boundary element method

This paper presents an analysis of mixed mode fracture mechanics problems arising in anisotropic composite laminates. The boundary element method (BEM) and the J _k integral are presented as accurate techniques to compute the stress intensity factors K _I and K _II of two dimensional anisotropic bodies. Using function of a complex variable a decoupling procedure is derived to obtain the stress intensity factors. The procedure is based on the computation of the J ₁-integral and of the ratio of relative displacements at the crack faces, near the crack tip. Applications are presented for unidirectional and symmetric laminates of glass, boron and graphite-epoxy materials. Numerical examples of problems of pure mode I and mixed mode deformations are given, in order to demonstrate the accuracy of the method.     

Mixed Mode I,II and III analysis of multiple cracks in plane anisotropic solids by the BEM: a dislocation and point force approach

A direct boundary element method (BEM) for plane anisotropic elasticity is formulated for the generalized plane strain. It deals with the general case when the in-plane and out-of-plane deformations are coupled, including the special case when they are decoupled. The formulation is based on the distributions of point forces and dislocation dipoles following the physical interpretation of Somigliana's identity. We adopt Lekhnitskii-Eshelby-Stroh formalism for anisotropic elasticity and represent the point force and the dislocation, their dipoles, and continuous distributions systematically; the duality relations between the point force and the dislocation solutions are fully exploited. The analytical formulas for the displacement and the traction BEM are applied to the mixed mode crack analysis for multiply cracked anisotropic bodies. We extend the physical interpretation of Somigliana's identity to cracked bodies and represent the crack by the continuous distribution of dislocation dipoles. The mixed mode stress intensity factors (K_I, K_II and K_III) are determined accurately with the help of the conservation integrals of anisotropic elasticity.     

p-version hierarchical three dimensional curved shell element for elastostatics

This paper presents a hierarchical three dimensional curved shell finite element formulation based on the p-approximation concept. The element displacement approximation can be of arbitrary and different polynomial orders in the plane of the shell (ξ, η) and the transverse direction (ξ). The curved shell element approximation functions and the corresponding nodal variables are derived by first constructing the approximation functions of orders p_ξ, p_η and p_ξ and the corresponding nodal variable operators for each of the three directions ξ, η and ξ and then taking their products (sometimes also known as tensor product). This procedure gives the approximation functions and the corresponding nodal variables corresponding to the polynomial orders p_ξ, p_η and p_ξ. Both the element displacement functions and the nodal variables are hierarchical; therefore, the resulting element matrices and the equivalent nodal load vectors are hierarchical also, i.e. the element properties corresponding to the polynomial orders p_ξ, p_η and p_ξ are a subset of those corresponding to the orders (p_ξ + 1), (p_η +1) and (p_ξ +1). The formulation guarantees C° continuity or smoothness of the displacement field across the interelement boundaries. The geometry of the element is described by the co-ordinates of the nodes on its middle surface (ξ = 0) and the nodal vectors describing its bottom (ξ = ?1) and top (ξ = +1) surfaces. The element properties are derived using the principle of virtual work and the hierarchical element approximation. The formulation is equally effective for very thin as well as very thick plates and curved shells. In fact, in many three dimensional applications the element can be used to replace the hierarchical three dimensional solid element without loss of accuracy but significant gain in modelling convenience. Numerical examples are presented to demonstrate the accuracy, efficiency and overall superiority of the present formulation. The results obtained from the present formulation are compared with those available in the literature as well as analytical solutions.     

A boundary element method for calculating the K field for cracks along a bimaterial interface

A multi-region boundary element method (BEM) based on the modified crack closure method (CCM) is developed to obtain the energy release rate G for cracks in homogeneous materials and along a bimaterial interface. The energy release rate obtained using the CCM are compared with that obtained using the crack opening displacement (COD) method. A combination of these methods allows us to determine the phase angle and therefore the complex stress intensity factor K for crack problems. We access the accuracy of our BEM by comparing its results with known analytic solutions and previous FEM results in the literature. Computations are also carried out for the asymmetric double cantilever beam (ADCB) specimen, which has been used to determine fracture toughness of polymer/polymer and polymer/nonpolymer interfaces. An auxiliary K _A field method to evaluate K is also discussed.This work is supported by the Cornell Materials Science Center (MSC) which is funded by the National Science Foundation (DMR-MRL program).     

Dislocation and point-force-based approach to the special Green's Function BEM for elliptic hole and crack problems in two dimensions

In this paper we give the theoretical foundation for a dislocation and point-force-based approach to the special Green's function boundary element method and formulate, as an example, the special Green's function boundary element method for elliptic hole and crack problems. The crack is treated as a particular case of the elliptic hole. We adopt a physical interpretation of Somigliana's identity and formulate the boundary element method in terms of distributions of point forces and dislocation dipoles in the infinite domain with an elliptic hole. There is no need to model the hole by the boundary elements since the traction free boundary condition there for the point force and the dislocation dipole is automatically satisfied. The Green's functions are derived following the Muskhelishvili complex variable formalism and the boundary element method is formulated using complex variables. All the boundary integrals, including the formula for the stress intensity factor for the crack, are evaluated analytically to give a simple yet accurate special Green's function boundary element method. The numerical results obtained for the stress concentration and intensity factors are extremely accurate. © 1997 John Wiley & Sons, Ltd.