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.     

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.     

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.     

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.     

An indirect elastodynamic boundary element method with analytic bases

An indirect time‐domain boundary element method (BEM) is presented here for the treatment of 2D elastodynamic problems. The approximated solution in this method is formulated as a linear combination of a set of particular solutions, which are called bases. The displacement and stress fields of a basis are analytically derived by means of solving Lame's displacement potentials. A semi‐collocation method is proposed to be the time‐stepping algorithm. This method is equivalent to a displacement discontinuity method with piecewise linear discontinuities in both space and time. The resulting time‐stepping scheme is explicit. The BEM is implemented to solve three numerical examples, Lamb's problem, half‐plane with a buried crack and Selberg's problem. Though Lamb's problem is considered a difficult problem for numerical methods, the current numerical results for the surface displacements show accurately the characteristics of the Rayleigh wave. This method is efficient and accurate. Copyright © 2003 John Wiley & Sons, Ltd.     

A hyper-singular numerical Green's function generation for BEM applied to dynamic SIF problems

This paper introduces the extension of the numerical Green's function approach for elastodynamic fracture mechanics problems. The formulation uses the hyper-singular boundary integral equation to obtain the fundamental solution for the cracked unbounded medium. The procedure is general and can be applied to multiple crack problems of general geometry. Applications to time harmonic and transient (through inverse numerical Fourier and Laplace transforms) stress intensity factor (SIF) computations are presented and compared with other numerical and analytical results, showing the good accuracy of the present strategy for these kinds of problems.     

3D non-linear time domain FEM–BEM approach to soil–structure interaction problems

Dynamic soil–structure interaction is concerned with the study of structures supported on flexible soils and subjected to dynamic actions. Methods combining the finite element method (FEM) and the boundary element method (BEM) are well suited to address dynamic soil–structure interaction problems. Hence, FEM–BEM models have been widely used. However, non-linear contact conditions and non-linear behavior of the structures have not usually been considered in the analyses. This paper presents a 3D non-linear time domain FEM–BEM numerical model designed to address soil–structure interaction problems. The BEM formulation, based on element subdivision and the constant velocity approach, was improved by using interpolation matrices. The FEM approach was based on implicit Green's functions and non-linear contact was considered at the FEM–BEM interface. Two engineering problems were studied with the proposed methodology: the propagation of waves in an elastic foundation and the dynamic response of a structure to an incident wave field.     

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.     

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.     

Structural integrity assessment of FBR components using a distributed computing environment

In this work, the structural integrity of the critical components of the fast breeder reactor (FBR) that are subjected thermal striping is assessed using a fracture mechanics approach based on linear elastic fracture mechanics (LEFM). The structural integrity is assessed in terms of the actual life of the component for a particular difference between the hot and cold liquid temperatures at the critical mixing velocities. A generalized procedure is attempted for the computation of fatigue life. It is demonstrated in this work that the analysis procedure adopted is computationally very efficient. Green's function method is used for transient mode I crack propagation analysis. An inherent parallelism in the method is exploited for computational efficiency. A distributed computing environment is, therefore, used to demonstrate the effectiveness of Green's function method for crack propagation analysis for the kind of problem solved in this work. A simple idealization in the form of flat plate geometry is used in a numerical example to show the computational efficiency. The method shows a good scale‐up justifying the benefit of using a distributed computing environment given a large amount of input data for the thermal striping problem.     

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.     

Numerical instabilities of the integral approach to the interior boundary-value problem for the two-dimensional Helmholtz equation

The integral equations arising from the Green's formula, applied to the two-dimensional Helmholtz equation defined in a limited domain, are considered and the presence of instabilities in their numerical solution, when a real Green's function is adopted, is pointed out. A complete study has been carried out for a circular domain and the conditions under which such instabilities can occur in a domain of arbitrary geometry have been investigated. In particular, it has been shown that in every case the use of a complex Green's function is able to avoid their presence.     

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.     

Elastoplastic BIE analysis of cracked plates and related problems. Part 2: Numerical results

Part 1 of this paper reports on the formulation of an advanced boundary—integral equation model for fracture mechanics analysis of cracked plates, subject to elastoplastic behaviour or other, related body force problems. The basis of this formulation contrasts with other BIE elastoplastic formulations in the use of the Green's function for an infinite plane containing a stress free crack. This Green's function formulation assures that the total elastic strain field for the crack problem is accurately imbedded in the numerical model. The second part of this paper reports on the numerical implementation of this algorithm, as currently developed. The anelastic strain field (residual strains, thermal strains, plastic strains, etc.) is approximated as piecewise constant, while the boundary data is modelled with linear interpolations. An iteration solution scheme is adopted which eliminates the need for recalculation of the BIE matrices. The stability and accuracy of the algorithm are demonstrated for an uncracked, notch geometry, and comparison to finite element results is made for the centre-cracked panel. The data shows that even the crude plastic strain model applied is capable of excellent resolution of crack tip plastic behaviour.     

Heritage and early history of the boundary element method

This article explores the rich heritage of the boundary element method (BEM) by examining its mathematical foundation from the potential theory, boundary value problems, Green's functions, Green's identities, to Fredholm integral equations. The 18th to 20th century mathematicians, whose contributions were key to the theoretical development, are honored with short biographies. The origin of the numerical implementation of boundary integral equations can be traced to the 1960s, when the electronic computers had become available. The full emergence of the numerical technique known as the boundary element method occurred in the late 1970s. This article reviews the early history of the boundary element method up to the late 1970s.     

Efficient numerical models for the prediction of acoustic wave propagation in the vicinity of a wedge coastal region

In this paper, numerical frequency domain formulations are developed to simulate the 2D acoustic wave propagation in the vicinity of an underwater configuration which combines two sub-regions: the first one consists of a wedge with rigid seabed and free surface, and the second one is assumed to have a rigid flat bottom and a free flat surface.The problem is solved using two different numerical methods: the Boundary Element Method (BEM) and the Method of Fundamental Solutions (MFS). Two models are developed by using a sub-region technique, where only the vertical interface between sub-regions of different geometries has to be discretized. These formulations incorporate Green's functions that take into account the presence of flat rigid and free surfaces and of a wedge. Green's functions are defined using two approaches: the image source method is used to model the rigid flat bottom and free flat interface, whereas the response provided by the wedge sub-region is based on a normal mode solution. Additionally, a MFS and a BEM model are also implemented which require the discretization of the sloping rigid seabed of the wedge, therefore making use of Green's functions for a rigid flat bottom and a free surface (using the image source method).A detailed discussion on the performance of these formulations is performed, with the aim of finding an efficient formulation to solve the problem. It is found that the model based on the MFS and on the sub-region technique has a significantly lower computational cost and is stable, therefore being the most suitable for the analysis of acoustic wave propagation in the studied configurations.     

BEM calculation of the complex thermal impedance of microelectronic devices

This paper presents a numerical method for modelling the dynamic thermal behaviour of microelectronic structures in the frequency domain. A boundary element method (BEM) based on a Green's function solution is proposed for solving the 3D heat equation in phasor notation. The method is capable of calculating the AC temperature and heat flux distributions and complex thermal impedance for packages composed of an arbitrary number of bar-shaped components. Various types of boundary conditions, including thermal contact resistance and convective cooling, can be taken into account. A simple benchmark case is investigated and a good convergence towards the analytical solution is obtained. Simulation results for a thin plate under convective cooling are compared with a theoretical model and an excellent agreement is observed. In a second example a more complicated three-layer structure is investigated. The BEM is used to analyse the thermal behaviour if delamination of the package occurs, and a physical explanation for the results is given.     

Efficient computation of the Green's function and its derivatives for three-dimensional anisotropic elasticity in BEM analysis

An alternative scheme to compute the Green's function and its derivatives for three dimensional generally anisotropic elastic solids is presented in this paper. These items are essential in the formulation of the boundary element method (BEM); their evaluation has remained a subject of interest because of the mathematical complexity. The Green's function considered here is the one introduced by Ting and Lee [Q. J. Mech. Appl. Math. 1997; 50: 407–26] which is of real-variable, explicit form expressed in terms of Stroh's eigenvalues. It has received attention in BEM only quite recently. By taking advantage of the periodic nature of the spherical angles when it is expressed in the spherical coordinate system, it is proposed that this Green's function be represented by a double Fourier series. The Fourier coefficients are determined numerically only once for a given anisotropic material; this is independent of the number of field points in the BEM analysis. Derivatives of the Green's function can be performed by direct spatial differentiation of the Fourier series. The resulting formulations are more concise and simpler than those derived analytically in closed form in previous studies. Numerical examples are presented to demonstrate the veracity and superior efficiency of the scheme, particularly when the number of field points is very large, as is typically the case when analyzing practical three dimensional engineering problems.