An iterative boundary element method for Cauchy inverse problems

 We are interested in this paper in recovering an harmonic function from the knowledge of Cauchy data on some part of the boundary. A new inversion method is introduced. It reduces the Cauchy problem resolution to the determination of the resolution of a sequence of well-posed problems. The sequence of these solutions is proved to converge to the Cauchy problem solution. The algorithm is implemented in the framework of boundary elements. Displayed numerical results highlight its accuracy, as well as its robustness to noisy data. Received 6 November 2000     

Identification of material properties and cavities in two-dimensional linear elasticity

 In this paper, the simultaneous identification of the Poisson ratio, the shear modulus and a circular cavity embedded in an isotropic linear elastic material from boundary measurements is investigated. The numerical method proposed is based on the least-squares minimisation of the errors between the measured and the calculated tractions on the outer boundary using the boundary element method (BEM). Received: 26 July 2002 / Accepted: 4 February 2003 L. Marin would like to acknowledge the financial support received from the EPSRC. The very constructive comments made by the referees on the first version of this paper are gratefully acknowledged.     

Conjugate gradient-boundary element solution to the Cauchy problem for Helmholtz-type equations

 In this paper, an iterative algorithm based on the conjugate gradient method (CGM) in combination with the boundary element method (BEM) for obtaining stable approximate solutions to the Cauchy problem for Helmholtz-type equations is analysed. An efficient regularising stopping criterion for CGM proposed by Nemirovskii [25] is employed. The numerical results obtained confirm that the CGM + BEM produces a convergent and stable numerical solution with respect to increasing the number of boundary elements and decreasing the amount of noise added into the input data. Received: 5 November 2002 / Accepted: 5 March 2003 L. Marin would like to acknowledge the financial support received from the EPSRC. The authors would like to thank Professor Dinh Nho Hào and Dr. Thomas Johansson for some useful discussions and suggestions.     

A general algorithm for the numerical evaluation of nearly singular boundary integrals of various orders for two- and three-dimensional elasticity

 A general algorithm of the distance transformation type is presented in this paper for the accurate numerical evaluation of nearly singular boundary integrals encountered in elasticity, which, next to the singular ones, has long been an issue of major concern in computational mechanics with boundary element methods. The distance transformation is realized by making use of the distance functions, defined in the local intrinsic coordinate systems, which plays the role of damping-out the near singularity of integrands resulting from the very small distance between the source and the integration points. By taking advantage of the divergence-free property of the integrals with the nearly hypersingular kernels in the 3D case, a technique of geometric conversion over the auxiliary cone surfaces of the boundary element is designed, which is suitable also for the numerical evaluation of the hypersingular boundary integrals. The effects of the distance transformations are studied and compared numerically for different orders in the 2D case and in the different local systems in the 3D case using quadratic boundary elements. It is shown that the proposed algorithm works very well, by using standard Gaussian quadrature formulae, for both the 2D and 3D elastic problems. Received: 20 November 2001 / Accepted: 4 June 2002 The work was supported by the Science Foundation of Shanghai Municipal Commission of Education.     

Domain decomposition in the symmetric boundary element analysis

 Recent developments in the symmetric boundary element method (SBEM) have shown a clear superiority of this formulation over the collocation method. Its competitiveness has been tested in comparison to the finite element method (FEM) and is manifested in several engineering problems in which internal boundaries are present, i.e. those in which the body shows a jump in the physical characteristics of the material and in which an appropriate study of the response must be used. When we work in the ambit of the SBE formulation, the body is subdivided into macroelements characterized by some relations which link the interface boundary unknowns to the external actions. These relations, valid for each macroelement and characterized by symmetric matricial operators, are similar in type to those obtainable for the FEM. The assembly of the macroelements based on the equilibrium conditions, or on the compatibility conditions, or on both of these conditions leads to three analysis methods: displacement, force, and mixed-value methods, respectively. The use of the fundamental solutions involves the punctual satisfaction of the compatibility and of the equilibrium inside each macroelement and it causes a stringent elastic response close to the actual solution. Some examples make it possible to perform numerical checks in comparison with solutions obtainable in closed form. These checks show that the numerical solutions are floating ones when the macroelement geometry obtained by subdividing the body changes. Received 26 January 2001     

A first order method for the Cauchy problem for the Laplace equation using BEM

The purpose is to propose an improved method for inverse boundary value problems. This method is presented on a model problem. It introduces a higher order problem. BEM numerical simulations highlight the efficiency, the improved accuracy, the robustness to noisy data of this new approach, as well as its ability to deblur noisy data.     

Evolutionary optimization in thermoelastic problems using the boundary element method

 The paper is devoted to application of evolutionary algorithms and the boundary element method to shape optimization of structures for various thermomechanical criteria, inverse problems of finding an optimal distribution of temperature on the boundary and identification of unknown boundary. Design variables are specified by Bezier curves. Several numerical examples of evolutionary computation are presented. Received 6 November 2000     

Modeling of 3D transversely piezoelectric and elastic bimaterials using the boundary element method

 This paper presents a numerical model for three-dimensional transversely isotropic bimaterials based on the boundary element formulation. The point force solutions expressed in a united-form for distinct eigenvalues are studied for transversely isotropic piezoelectricity and pure elasticity. A boundary integral formulation is implemented for the modeling of two-phase materials. In this study, the stress distributions are computed for a near interface flaw. The influences of the shape and location of the flaw on the the stress concentration are examined. The accuracy of the numerical procedures is validated through selected example problems and comparison studies. Received 3 October 2001 / Accepted 9 April 2002     

A Trefftz function approximation in local boundary integral equations

 In the present paper the Trefftz function as a test function is used to derive the local boundary integral equations (LBIE) for linear elasticity. Since Trefftz functions are regular, much less requirements are put on numerical integration than in the conventional boundary integral method. The moving least square (MLS) approximation is applied to the displacement field. Then, the traction vectors on the local boundaries are obtained from the gradients of the approximated displacements by using Hooke's law. Nodal points are randomly spread on the domain of the analysed body. The present method is a truly meshless method, as it does not need a finite element mesh, either for purposes of interpolation of the solution variables, or for the integration of the energy. Two ways are presented to formulate the solution of boundary value problems. In the first one the local boundary integral equations are written in all nodes (interior and boundary nodes). In the second way the LBIE are written only at the interior nodes and at the nodes on the global boundary the prescribed values of displacements and/or tractions are identified with their MLS approximations. Numerical examples for a square patch test and a cantilever beam are presented to illustrate the implementation and performance of the present method. Received 6 November 2000     

Regularization of hypersingular boundary integral equations: a new approach for axisymmetric elasticity

A hypersingular boundary integral equation (HBIE) formulation, for axisymmetric linear elasticity, has been recently presented by de Lacerda and Wrobel [Int. J. Numer. Meth. Engng 52 (2001) 1337]. The strongly singular and hypersingular equations in this formulation are regularized by de Lacerda and Wrobel by employing the singularity subtraction technique. The present paper revisits the same problem. The axisymmetric HBIE formulation for linear elasticity is interpreted here in a ‘finite part’ sense and is then regularized by employing a ‘complete exclusion zone’. The resulting regularized equations are shown to be simpler than those by de Lacerda and Wrobel.     

Direct evaluation of Cauchy-principal-value integrals in boundary elements for infinite and semi-infinite three-dimensional domains

Algorithms for the direct numerical evaluation of Cauchy-principal-value integrals in three-dimensional problems have recently been implemented. An application of one of these algorithms to infinite boundary elements and to semi-infinite special fundamental solutions is discussed.     

The boundary contour method for two-dimensional Stokes flow and incompressible elastic materials

 A variant of the boundary element method, called the boundary contour method (BCM), offers a further reduction in dimensionality. Consequently, boundary contour analysis of two-dimensional (2-D) problems does not require any numerical integration at all. While the method has enjoyed many successful applications in linear elasticity, the above advantage has not been exploited for Stokes flow problems and incompressible media. In order to extend the BCM to these materials, this paper presents a development of the method based on the equations of Stokes flow and its 2-D kernel tensors. Potential functions are derived for quadratic boundary elements. Numerical solutions for some well-known examples are compared with the analytical ones to validate the development. Received 28 August 2001 / Accepted 15 January 2002     

An alternating iterative algorithm for the Cauchy problem in anisotropic elasticity

The alternating iterative algorithm proposed by Kozlov et al. [An iterative method for solving the Cauchy problem for elliptic equations. USSR Comput Math Math Phys 1991;31:45–52] for obtaining approximate solutions to the Cauchy problem in two-dimensional anisotropic elasticity is analysed and numerically implemented using the boundary element method (BEM). The ill-posedness of this inverse boundary value problem is overcome by employing an efficient regularising stopping. The numerical results confirm that the iterative BEM produces a convergent and stable numerical solution with respect to increasing the number of boundary elements and decreasing the amount of noise added into the input data.     

A mixed least squares method for solving problems in linear elasticity: theoretical study

 In a previous paper we proposed a mixed least squares method for solving problems in linear elasticity. The solution to the equations of linear elasticity was obtained via minimization of a least squares functional depending on displacements and stresses. The performance of the method was tested numerically for low order elements for classical examples with well known analytical solutions. In this paper we derive a condition for the existence and uniqueness of the solution of the discrete problem for both compressible and incompressible cases, and verify the uniqueness of the solution analytically for two low order piece-wise polynomial FEM spaces. Received: 20 January 2001 / Accepted: 14 June 2002 The authors gratefully acknowledge the financial support provided by NASA George C. Marshall Space Flight Centre under contract number NAS8-38779.     

Boundary element model of microsegregation during volumetric solidification of binary alloy

 In the paper the mathematical model of heat and mass transfer processes proceeding in the domain of casting is discussed. In particular, the volumetric solidification of a binary alloy under the assumption that the temperature is only time-dependent is analyzed. From the practical point of view such situation takes place when the casting is produced in typical moulding sand. The differential equation describing the course of solidification and cooling processes is presented in Sect. 1. In this equation the capacity of an internal heat source results from the Mehl–Johnson–Avrami–Kolmogorov theory [1, 2], at the same time the constant number of nuclei is accepted, while the rate of the solid phase growth is proportional to the second power of undercooling below the liquidus temperature. The macroscopic model is coupled with a microsegregation one (Sect. 2). This process is analyzed at the level of a single grain. The distribution of the alloy component in the control volume corresponding to the final grain radius is found on a basis of the boundary element method using discretization in time [3, 4]. The examples of numerical computations are also presented. Received 6 November 2000     

Finite parts of singular and hypersingular integrals with irregular boundary source points

The subject of this paper is the evaluation of finite parts (FPs) of certain singular and hypersingular integrals, that appear in boundary integral equations (BIEs), when the source point is an irregular boundary point (situated at a corner on a one-dimensional plane curve or at a corner or edge on a two-dimensional surface). Two issues addressed in this paper are: an unified, consistent and practical definition of a FP with an irregular boundary source point, and numerical evaluation of such integrals together with solution strategies for hypersingular BIEs (HBIEs). The proposed formulation is compared with others that are available in the literature and interesting connections are made between this formulation and those of other researchers.     

The least-squares meshfree method for solving linear elastic problems

 A meshfree method based on the first-order least-squares formulation for linear elasticity is presented. In the authors' previous work, the least-squares meshfree method has been shown to be highly robust to integration errors with the numerical examples of Poisson equation. In the present work, conventional formulation and compatibility-imposed formulation for linear elastic problems are studied on the convergence behavior of the solution and the robustness to the inaccurate integration using simply constructed background cells. In the least-squares formulation, both primal and dual variables can be approximated by the same function space. This can lead to higher rate of convergence for dual variables than Galerkin formulation. In general, the incompressible locking can be alleviated using mixed formulations. However, in meshfree framework these approaches involve an additional use of background grids to implement lower approximation space for dual variables. This difficulty is avoided in the present method and numerical examples show the uniform convergence performance in the incompressible limit. Therefore the present method has little burden of the requirement of background cells for the purposes of integration and alleviating the incompressible locking. Received: 16 December 2001 / Accepted: 4 November 2002     

A semi‐analytical curved element for linear elasticity based on the scaled boundary finite element method

This work introduces a semi‐analytical formulation for the simulation and modeling of curved structures based on the scaled boundary finite element method (SBFEM). This approach adapts the fundamental idea of the SBFEM concept to scale a boundary to describe a geometry. Until now, scaling in SBFEM has exclusively been performed along a straight coordinate that enlarges, shrinks, or shifts a given boundary. In this novel approach, scaling is based on a polar or cylindrical coordinate system such that a boundary is shifted along a curved scaling direction. The derived formulations are used to compute the static and dynamic stiffness matrices of homogeneous curved structures. The resulting elements can be coupled to general SBFEM or FEM domains. For elastodynamic problems, computations are performed in the frequency domain. Results of this work are validated using the global matrix method and standard finite element analysis. Copyright © 2016 John Wiley & Sons, Ltd.     

The traction boundary contour method for linear elasticity

This paper presents a further development of the boundary contour method. The boundary contour method is extended to cover the traction boundary integral equation. A traction boundary contour method is proposed for linear elastostatics. The formulation of traction boundary contour method is regular for points except the ends of the boundary element and corners. The present approach only requires line integrals for three‐dimensional problems and function evaluations at the ends of boundary elements for two‐dimensional cases. The implementation of the traction boundary contour method with quadratic boundary elements is presented for two‐dimensional problems. Numerical results are given for some two‐dimensional examples, and these are compared with analytical solutions. This method is shown to give excellent results for illustrative examples. Copyright © 1999 John Wiley & Sons, Ltd.     

A fast multipole boundary element method for solving the thin plate bending problem

A fast multipole boundary element method (BEM) for solving large-scale thin plate bending problems is presented in this paper. The method is based on the Kirchhoff thin plate bending theory and the biharmonic equation governing the deflection of the plate. First, the direct boundary integral equations and the conventional BEM for thin plate bending problems are reviewed. Second, the complex notation of the kernel functions, expansions and translations in the fast multipole BEM are presented. Finally, a few numerical examples are presented to show the accuracy and efficiency of the fast multipole BEM in solving thin plate bending problems. The bending rigidity of a perforated plate is evaluated using the developed code. It is shown that the fast multipole BEM can be applied to solve plate bending problems with good accuracy. Possible improvements in the efficiency of the method are discussed.