A numerical green's function approach for boundary elements applied to fracture mechanics

The most accurate boundary element formulation to deal with fracture mechanics problems is obtained with the implementation of the associated Green's function acting as the fundamental solution. Consequently, the range of applications of this formulation is dependent on the availability of the appropriate Green's function for actual crack geometry. Analytical Green's functions have been presented for a few single crack configurations in 2-D applications and require complex variable theory. This work extends the applicability of the formulation through the introduction of efficient numerical means of computing the Green's function components for single or multiple crack problems, of general geometry, including the implementation to 3-D problems as a future development. Also, the approach uses real variables only and well-established boundary integral equations.     

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.     

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.     

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.     

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.     

BEM for crack‐inclusion problems of plane thermopiezoelectric solids

The problem of interactions between an inclusion and multiple cracks in a thermopiezoelectric solid is considered by boundary element method (BEM) in this paper. First of all, a BEM for the crack–inclusion problem is developed by way of potential variational principle, the concept of dislocation, and Green's function. In the BE model, the continuity condition of the interface between inclusion and matrix is satisfied, a priori, by the Green's function, and not involved in the boundary element equations. This is then followed by expressing the stress and electric displacement (SED) and elastic displacements and electric potential (EDEP) in terms of polynomials of complex variables ξ_t and ξ_k in the transformed ξ‐plane in order to simulate SED intensity factors by the BEM. The least‐squares method incorporating the BE formulation can, then, be used to calculate SED intensity factors directly. Numerical results for a piezoelectric plate with one inclusion and a crack are presented to illustrate the application of the proposed formulation. Copyright © 2000 John Wiley & Sons, Ltd.     

A numerical Green's function BEM formulation for crack growth simulation

This paper presents a crack growth prediction analysis based on the numerical Green's function (NGF) procedure and on the minimum strain energy density criterion for crack extension, also known as S-criterion. In the NGF procedure, the hypersingular boundary integral equation is used to numerically obtain the Green's function which automatically includes the crack into the fundamental infinite medium. When solving a linear elastic fracture mechanisms (LEFM) problem, once the NGF is obtained, the classical boundary element method can be used to determine the external boundary unknowns and, consequently, the stress intensity factors needed to predict the direction and increment of crack growth. With the change in crack geometry, another numerical analysis is carried out without need to rebuilding the entire element discretization, since only the crack built in the NGF needs update. Numerical examples, contemplating crack extensions for two-dimensional LEFM problems, are presented to illustrate the procedure.     

Boundary element modeling of cracks in piezoelectric solids

This study focuses on the application of boundary element methods for linear fracture mechanics of two-dimensional piezoelectric solids. A complete set of piezoelectric Green's functions, based on the extended Lekhnitskii's formalism and distributed dislocation modeling, are presented. Special Green's functions are obtained for an infinite medium containing a conducting crack or an impermeable crack. The numerical solution of the boundary integral equation and the computation of fracture parameters are discussed. The concept of crack closure integral is utilized to calculate energy release rates. Accuracy of the boundary element solutions is confirmed by comparing with analytical solutions reported in the literature. The present scheme can be applied to study complex cracks such as branched cracks, forked cracks and microcrack clusters.     

A boundary integral approach for plane analysis of electrically semi-permeable planar cracks in a piezoelectric solid

A plane electro-elastostatic problem involving arbitrarily located planar stress free cracks which are electrically semi-permeable is considered. Through the use of the numerical Green's function for impermeable cracks, the problem is formulated in terms of boundary integral equations which are solved numerically by a boundary element procedure together with a predictor–corrector method. The crack tip stress and electric displacement intensity factors can be easily computed once the boundary integral equations are properly solved.     

The singular function boundary integral method for a two-dimensional fracture problem

The singular function boundary integral method (SFBIM) originally developed for Laplacian problems with boundary singularities is extended for solving two-dimensional fracture problems formulated in terms of the Airy stress function. Our goal is the accurate, direct computation of the associated stress intensity factors, which appear as coefficients in the asymptotic expansion of the solution near the crack tip. In the SFBIM, the leading terms of the asymptotic solution are used to approximate the solution and to weight the governing biharmonic equation in the Galerkin sense. The discretized equations are reduced to boundary integrals by means of Green's theorem and the Dirichlet boundary conditions are weakly enforced by means of Lagrange multipliers. The numerical results on a model problem show that the method converges extremely fast and yields accurate estimates of the leading stress intensity factors.     

A Green's function time‐domain boundary element method for the elastodynamic half‐plane

The transient Green's function of the 2‐D Lamb's problem for the general case where point source and receiver are situated beneath the traction‐free surface is derived. The derivations are based on Laplace‐transform methods, utilizing the Cagniard–de Hoop inversion. The Green's function is purely algebraic without any integrals and is presented in a numerically applicable form for the first time. It is used to develop a Green's function BEM in which surface discretizations on the traction‐free boundary can be saved. The time convolution is performed numerically in an abstract complex plane. Hence, the respective integrals are regularized and only a few evaluations of the Green's function are required. This fast procedure has been applied for the first time. The Green's function BEM developed proved to be very accurate and efficient in comparison with analogue BEMs that employ the fundamental solution. Copyright © 1999 John Wiley & Sons, Ltd.     

Local discretization error bounds using interval boundary element method

In this paper, a method to account for the point‐wise discretization error in the solution for boundary element method is developed. Interval methods are used to enclose the boundary integral equation and a sharp parametric solver for the interval linear system of equations is presented. The developed method does not assume any special properties besides the Laplace equation being a linear elliptic partial differential equation whose Green's function for an isotropic media is known. Numerical results are presented showing the guarantee of the bounds on the solution as well as the convergence of the discretization error. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

Elastoplastic BIE analysis of cracked plates and related problems. Part 1: Formulation

The boundary-integral equation formulation for two-dimensional, planar fracture mechanics based on the use of a special Green's function forms the basis of this analytical paper. The Green's function method is extended to problems of anelastic strain distributions (e.g. elastoplasticity, thermal gradients, residual strains) through a volume (area) integral. The role of the elastic Green's function for the crack problem on the distribution of elastoplastic strains is reviewed. Further, new analytical results for elastic stress intensity factor models for the residual strain and thermal gradient problems are presented. Part 2 of this paper outlines the numerical solution strategy and results for several test problems.     

Numerical Green's functions for some electroelastic crack problems

A plane electroelastic problem involving planar cracks in a piezoelectric body is considered. The deformation of the body is assumed to be independent of time and one of the Cartesian coordinates. The cracks are traction free and are electrically either permeable or impermeable. Numerical Green's functions which satisfy the boundary conditions on the cracks are derived using the hypersingular integral approach and applied to obtain a boundary integral solution for the electroelastic crack problem considered here. As the conditions on the cracks are built into the Green's functions, the boundary integral solution does not contain integrals over the cracks. It is used to derive a boundary element procedure for computing the crack tip stress and electrical displacement intensity factors.     

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.     

On the Theory of Direct and Inverse Problems in the Elastic Rectangle: Antiplane Case

In this article, we study the antiplane deformation of the boundary surface of a rectangular domain in the presence of a void and a shear force on the outer boundary surface. For a formulated inverse problem, we develop some analytical results and use them to solve the problem numerically for various elliptic geometrical configurations. The analytical method allows us to give an efficient representation for Green's function in the rectangular domain. Then we derive the same Green's function by an alternative method based on Fourier series expansions. Finally, for a number of configurations, we demonstrate the comparison between real and reconstructed defects.     

Three-dimensional extended displacement discontinuity method for vertical cracks in transversely isotropic piezoelectric media

The conventional displacement discontinuity method is extended to study a vertical crack under electrically impermeable condition, running parallel to the poling direction and normal to the plane of isotropy in three-dimensional transversely isotropic piezoelectric media. The extended Green's functions specifically for extended point displacement discontinuities are derived based on the Green's functions of extended point forces and the Somigliana identity. The hyper-singular displacement discontinuity boundary integral equations are also derived. The asymptotical behavior near the crack tips along the crack front is studied and the ordinary 1/2 singularity is obtained at the tips. The extended field intensity factors are expressed in terms of the extended displacement discontinuity on crack faces. Numerical results on the extended field intensity factors for a vertical square crack are presented using the proposed extended displacement discontinuity method.