The first-kind and the second-kind boundary integral equation systems for solution of frictionless contact problems

An original approach to the numerical solution of displacement boundary integral equation (BIE) and traction hypersingular boundary integral equation (HBIE) by the boundary element method (BEM) for contact problems is given. The main point is to show, how the contact conditions are used to formulate the first-kind and the second-kind BIE systems in the case of frictionless two-body elastic contact. The solution of the first-kind BIE is performed by symmetric Galerkin BEM; the second-kind BIE is solved by an appropriate collocation BEM. The contact problem in itself is solved by the method of subsequent approximations of contact region. Both forms of BIE system are compared in several numerical examples. This comparison is made for different kinds of contact problem. The major emphasis is put on the evaluation of contact pressure. The obtained results are compared with referenced numerical and with the analytical ones.     

Traction BIE solutions for flat cracks

The paper deals with the numerical solution techniques for the traction boundary integral equation (BIE), which describes the opening (and sliding) displacements of the surface of the traction loaded crack or arbitrary planform embedded in an elastic infinite body (buried crack problem). The traction BIE is a singular integral equation of the first kind for the displacement gradients. Its solution poses a number of numerical problems, such as the presence of derivatives of the unknown function in the integral equation, the modeling of the crack front displacement gradient singularity, and the regularization of the equation's singular kernels. All of the above problems have been addressed and solved. Details of the algorithm are provided. Numerical results of a number of crack configurations are presented, demonstrating high accuracy of the method.     

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.     

Self‐regular boundary integral equation formulations for Laplace's equation in 2‐D

The purpose of this work is to demonstrate the application of the self‐regular formulation strategy using Green's identity (potential‐BIE) and its gradient form (flux‐BIE) for Laplace's equation. Self‐regular formulations lead to highly effective BEM algorithms that utilize standard conforming boundary elements and low‐order Gaussian integrations. Both formulations are discussed and implemented for two‐dimensional potential problems, and numerical results are presented. Potential results show that the use of quartic interpolations is required for the flux‐BIE to show comparable accuracy to the potential‐BIE using quadratic interpolations. On the other hand, flux error results in the potential‐BIE implementation can be dominated by the numerical integration of the logarithmic kernel of the remaining weakly singular integral. Accuracy of these flux results does not improve beyond a certain level when using standard quadrature together with a special transformation, but when an alternative logarithmic quadrature scheme is used these errors are shown to reduce abruptly, and the flux results converge monotonically to the exact answer. In the flux‐BIE implementation, where all integrals are regularized, flux results accuracy improves systematically, even with some oscillations, when refining the mesh or increasing the order of the interpolating function. The flux‐BIE approach presents a great numerical sensitivity to the mesh generation scheme and refinement. Accurate results for the potential and the flux were obtained for coarse‐graded meshes in which the rate of change of the tangential derivative of the potential was better approximated. This numerical sensitivity and the need for graded meshes were not found in the elasticity problem for which self‐regular formulations have also been developed using a similar approach. Logarithmic quadrature to evaluate the weakly singular integral is implemented in the self‐regular potential‐BIE, showing that the magnitude of the error is dependent only on the standard Gauss integration of the regularized integral, but not on this logarithmic quadrature of the weakly singular integral. The self‐regular potential‐BIE is compared with the standard (CPV) formulation, showing the equivalence between these formulations. The self‐regular BIE formulations and computational algorithms are established as robust alternatives to singular BIE formulations for potential problems. Copyright © 2001 John Wiley & Sons, Ltd.     

A dual BIE approach for large‐scale modelling of 3‐D electrostatic problems with the fast multipole boundary element method

A dual boundary integral equation (BIE) formulation is presented for the analysis of general 3‐D electrostatic problems, especially those involving thin structures. This dual BIE formulation uses a linear combination of the conventional BIE and hypersingular BIE on the entire boundary of a problem domain. Similar to crack problems in elasticity, the conventional BIE degenerates when the field outside a thin body is investigated, such as the electrostatic field around a thin conducting plate. The dual BIE formulation, however, does not degenerate in such cases. Most importantly, the dual BIE is found to have better conditioning for the equations using the boundary element method (BEM) compared with the conventional BIE, even for domains with regular shapes. Thus the dual BIE is well suited for implementation with the fast multipole BEM. The fast multipole BEM for the dual BIE formulation is developed based on an adaptive fast multiple approach for the conventional BIE. Several examples are studied with the fast multipole BEM code, including finite and infinite domain problems, bulky and thin plate structures, and simplified comb‐drive models having more than 440 thin beams with the total number of equations above 1.45 million and solved on a PC. The numerical results clearly demonstrate that the dual BIE is very effective in solving general 3‐D electrostatic problems, as well as special cases involving thin perfect conducting structures, and that the adaptive fast multipole BEM with the dual BIE formulation is very efficient and promising in solving large‐scale electrostatic problems. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

Antiplane crack analysis of a functionally graded material by a BIEM

Elastostatic analysis of an antiplane crack in a functionally graded material (FGM) is performed by using a hypersingular boundary integral equation method (BIEM). An exponential law is applied to describe the spatial variation of the shear modulus of the FGM. A Galerkin method is applied for the numerical solution of the hypersingular traction BIE. Both unidirectional and bidirectional material gradations are investigated. Stress intensity factors for an infinite and linear elastic FGM containing a finite crack subjected to an antiplane crack-face loading are presented and discussed. The influences of the material gradients and the crack orientation on the stress intensity factors are analyzed.     

Boundary element‐free method (BEFM) and its application to two‐dimensional elasticity problems

In this study, we first discuss the moving least‐square approximation (MLS) method. In some cases, the MLS may form an ill‐conditioned system of equations so that the solution cannot be correctly obtained. Hence, in this paper, we propose an improved moving least‐square approximation (IMLS) method. In the IMLS method, the orthogonal function system with a weight function is used as the basis function. The IMLS has higher computational efficiency and precision than the MLS, and will not lead to an ill‐conditioned system of equations. Combining the boundary integral equation (BIE) method and the IMLS approximation method, a direct meshless BIE method, the boundary element‐free method (BEFM), for two‐dimensional elasticity is presented. Compared to other meshless BIE methods, BEFM is a direct numerical method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied easily; hence, it has higher computational precision. For demonstration purpose, selected numerical examples are given. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical solution for degenerate scale problem for exterior multiply connected region

Based on some previous publications, this paper investigates the numerical solution for degenerate scale problem for exterior multiply connected region. In the present study, the first step is to formulate a homogenous boundary integral equation (BIE) in the degenerate scale. The coordinate transform with a magnified factor, or a reduced factor h is performed in the next step. Using the property ln(hx)=ln(x)+lg(h), the new obtained BIE equation can be considered as a non-homogenous one defined in the transformed coordinates. The relevant scale in the transformed coordinates is a normal scale. Therefore, the new obtained BIE equation is solvable. Fundamental solutions are introduced. For evaluating the fundamental solutions, the right-hand terms in the non-homogenous equation, or a BIE, generally take the value of unit or zero. By using the obtained fundamental solutions, an equation for evaluating the magnified factor “h” is obtained. Finally, the degenerate scale is obtainable. Several numerical examples with two ellipses in an infinite plate are presented. Numerical solutions prove that the degenerate scale does not depend on the normal scale used in the process for evaluating the fundamental solutions.     

Boundary integral equation fracture mechanics analysis of plane orthotropic bodies

In this paper, a quick and efficient means of determining stress intensity factors, K _I and K _II, for cracks in generally orthotropic elastic bodies is presented using the numerical boundary integral equation (BIE) method. It is based on the use of quarter-point singular crack-tip elements in the quadratic isoparametric element formulation, similar to those commonly employed in the BIE fracture mechanics studies in isotropic elasticity. Analytical expressions which enable K _Iand K _II to be obtained directly from the BIE computed crack-tip nodal traction, or from the computed nodal displacements, of these elements are derived. Numerical results for a number of test problems are compared with those established in the literature. They are accurate even when only a very modest number of boundary elements are used.     

Boundary integral Equations Method in two- and three-dimensional problems of elastodynamics

In this paper the boundary integral equations method (BIEM) are considered for elastodynamic initial boundary value problems. It's known two approaches are discerned for account time. First of one is a combination of BIEM with Laplace (Fourier) transformation. This approach was suggested and realized by Cruse T.E. and Rizzo F. J. By them BIE in Laplace transformation space were obtained, investigated and some concrete problems were solved. This method was developed also by Manolis G. D., Beskos D. and other scholars for some dynamic problems solving.The second approach using retarding potentials was considered by Brebbia C. A., Fujiki K., Fukui T., Kato S., Kishima T., Kobayashi S., Nishimura N., Niwa Y., Manolis G. D. Mansur W.J. (for 2D elastodynamics), Chutoryansky N.M. (for 3D elastodynamics). Detailed review of abroad scholars elaborating BIEM was made by Beskos D. [7].This paper discusses BIEM for 2 and 3D elastodynamics on the base of the second approach. The fundamental solutions, integral representations and boundary integral equations are constructed by means distributions theory for the general case of anisotropic elastic media. It's suggested some new results concerning special regularization of singularities on the wave fronts of the integral equations kernels. The illustrative numerical examples concern the scattering of elastic waves on cavities embedded in an infinite isotropic medium. So, it's shown the numerical results of waves diffraction on the one and two cavities of arched and rectangular forms in 2 and 3D cases. These results show quite stability of the elaborating algorithm.     

Boundary element‐free method (BEFM) for two‐dimensional elastodynamic analysis using Laplace transform

In this paper, we present a direct meshless method of boundary integral equation (BIE), known as the boundary element‐free method (BEFM), for two‐dimensional (2D) elastodynamic problems that combines the BIE method for 2D elastodynamics in the Laplace‐transformed domain and the improved moving least‐squares (IMLS) approximation. The formulae for the BEFM for 2D elastodynamic problems are given, and the numerical procedures are also shown. The BEFM is a direct numerical method, in which the basic unknown quantities are the real solutions of the nodal variables, and the boundary conditions can be implemented directly and easily that leads to a greater computational precision. For the purpose of demonstration, some selected numerical examples are solved using the BEFM. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

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.     

Meshless boundary integral equations with equilibrium satisfaction

A known feature of any mixed interpolation boundary integral equations (BIE)-based methods is that equilibrium is not generally guaranteed in the numerical solution. Here, a complete meshless technique, based on the boundary element-free method (BEFM) with complete equilibrium satisfaction for 2D elastostatic analysis is proposed. The BEFM adopted is a meshless method based on boundary integral equations such as local boundary integral equation (LBIE) method and boundary node method (BNM), differing from them with respect to the integration domain and the approximation scheme.     

Three‐dimensional BEM for scattering of elastic waves in general anisotropic media

Based on the full‐space Green's functions, a three‐dimensional time‐harmonic boundary element method is presented for the scattering of elastic waves in a triclinic full space. The boundary integral equations for incident, scattered and total wave fields are given. An efficient numerical method is proposed to calculate the free terms for any geometry. The discretization of the boundary integral equation is achieved by using a linear triangular element. Applications are discussed for scattering of elastic waves by a spherical cavity in a 3D triclinic medium. The method has been tested by comparing the numerical results with the existing analytical solutions for an isotropic problem. The results show that, in addition to the frequency of the incident waves, the scattered waves strongly depend on the anisotropy of the media. Copyright © 2003 John Wiley & Sons, Ltd.     

Boundary-integral equation analysis of twisted internally cracked axisymmetric bimaterial elastic solids

 Green's function is obtained for the infinite bimaterial elastic solid, containing an internal circular interface crack, loaded by a unit tangential co-axial circular source. An axisymmetric direct boundary integral equation (BIE) is used for the analysis of a finite bimaterial axisymmetric body containing an internal circular interface crack and a finite homogeneous cracked cylinder, both under torsional loading. Using the proposed technique, no discretization of the crack surface is necessary. Numerical results for both examples as obtained by the proposed method are presented and discussed. Received: 29 October 2001 / Accepted: 29 May 2002     

The computation of the J integral with the traction boundary integral equation

The solutions of the displacement boundary integral equation (BIE) are not uniquely determined in certain types of boundary conditions. Traction boundary integral equations that have unique solutions in these traction and mixed boundary cases are established. For two‐dimensional linear elasticity problems, the divergence‐free property of the traction boundary integral equation is established. By applying Stokes' theorem, unknown tractions or displacements can be reduced to computation of traction integral potential functions at the boundary points. The same is true of the J integral: it is divergence‐free and the evaluation of the J integral can be inverted into the computation of the J integral potential functions at the boundary points of the cracked body. The J integral can be expressed as the linear combination of the tractions and displacements from the traction BIE on the boundary of the cracked body. Numerical integrals are not needed at all. Selected examples are presented to demonstrate the validity of the traction boundary integral and J integral. Copyright © 2004 John Wiley & Sons, Ltd.     

A method of combining rigid body displacements for the determination of rigid body rotations of structural members in the slope‐deflection method

Abstract

In this paper we investigate the stresses and displacements around an advancing face of a tunnel in an elastic half‐space as well as the induced three‐dimensional ground surface subsidence. An indirect boundary integral formulation for the analysis of tunnels in an elastic half‐space is derived and the half‐space fundamental solution (Mindlin's solution) is adopted to satisfy the traction‐free condition on the ground surface. For the evaluation of the non‐singular integral in the numerical implementation, a simple and efficient algorithm to automatically search for integration orders is developed by basing on the upper bound of error in Gaussian quadrature. Numerical results for some tunnel problems are studied to demonstrate the accuracy and versatility of the proposed method. In addition, the stresses, displacements, and the induced ground surface subsidence for an advancing circular tunnel in an elastic half‐space are calculated and presented in detail.     

3-D dynamic interaction between a penny-shaped crack and a thin interlayer joining two elastic half-spaces

Three-dimensional (3-D) elastodynamic interaction between a penny-shaped crack and a thin elastic interlayer joining two elastic half-spaces is investigated by an improved boundary integral equation method or boundary element method. The penny-shaped crack is embedded in one of the half-spaces, perpendicular to the interlayer and subjected to a time-harmonic tensile loading on its surfaces. Effective “spring-like” boundary conditions are applied to approximate the effects of the thin layer in the mathematical model. Integral representations for the displacement and the stress components are derived by using modified Green’s functions, which satisfy the “spring-like” boundary conditions identically. Then, application of the dynamic loading condition on the crack-surfaces results in a boundary integral equation (BIE) for the crack-opening-displacement over the crack-surfaces only. A solution procedure is developed for solving the BIE numerically. Numerical results for the mode-I dynamic stress intensity factor (SIF) are presented and discussed to show the variations of the mode-I dynamic SIF with the angular coordinate of the crack-front points, the dimensionless wave number, the material mismatch and the crack-layer distance.