A new boundary integral equation formulation for plane orthotropic elastic media

A new boundary integral equation formulation for solving plane elasticity problems involving orthotropic media is presented in this paper. Based on the real variable fundamental solutions of the considered problems, a limit theorem for the transformation from domain integral equations into boundary integral equations (BIEs) and a novel decomposition technique to the fundamental solutions, the regularized BIEs with indirect unknowns, which do not involve the direct calculation of CPV and HFP integrals, are established. The limiting process is done in global coordinates and no separate numerical treatment for strong and weak singular integrals was necessary. The current method does not need to transform the considered problems into isotropic ones as is normally done in the existing literature, so no inverse transform is required. The numerical implementation is carried out using both discontinuous quadratic elements and exact elements, which is developed to model its boundary with negligible error. The validity of the proposed scheme is demonstrated by three numerical examples. Excellent agreement between the numerical results and exact solutions was obtained even with using small amounts of element.     

The h-p version of the coupling of finite element and boundary element methods for transmission problems in polyhedral domains

Summary. This paper analyzes the rate of convergence of the h-p version of the coupling of the finite element and boundary element method for transmission problems with a linear differential operator with variable coefficients in a bounded polyhedral domain and with constant coefficients in the exterior domain . This procedure uses the variational formulation of the differential equation in and involves integral operators on the interface between and . The finite elements are used to obtain approximate solutions of the differential equation in and the boundary elements are used to obtain approximate solutions of the integral equations. For given piecewise analytic data we show that the Galerkin solution of this coupling procedure converges exponentially fast in the energy norm if the h-p version is used both for finite elements and boundary elements. Received February 10, 1996 / Revised version received April 4, 1997     

自适应快速多极正则化无网格法求解大规模三维位势问题

正则化无网格法（regularized meshless method, RMM）是一种新的边界型无网格数值离散方法．该方法克服了近年来引起广泛关注的基本解方法（method of fundamental solutions, MFS）的虚假边界缺陷,继承了其无网格、无数值积分、易实施等优点．另一方面,RMM方法同MFS方法的插值方程都涉及非对称稠密系数矩阵,运用常规代数方程的迭代法求解时都要求O(N2)量级的乘法计算量和存储量．随着问题自由度的增加,该方法的计算量增加极快,效率较低,一般难以计算大规模问题．为了克服这个缺点,利用对角形式的快速多级算法（fast multipole method, FMM）来加速RMM方法,发展了快速多级正则化无网格法（fast multipole regularized mesheless method, FM-RMM）．该方法无需数值积分并且具有O(N)量级的计算量和存储量,可有效地求解大规模工程问题．数值算例表明,FM-RMM算法可成功在内存为4GB的Core(TM)Ⅱ台式机上求解高达百万级自由度的三维位势问题．     

On boundary-volume element discretization of inhomogenous elastodynamic problems

This paper presents an integral equation formulation and its discretization scheme for the elastodynamic problem in which the material properties are prescribed as arbitrary, continuous and differentiable functions of the spatial coordinates. The formulation is made by using the Green's function for the corresponding problem in homogenous elasticity. From a weighted residual statement of the problem, the governing differential equation is transformed into a set of the integral equations in the inner domain as well as on the boundary. These integral equations are discretized by introducing a finite number of the boundary-volume-time elements, and the solution for the system of linear equations thus obtained is discussed.     

Local convergence of an adaptive scalar method and its application in a nonoverlapping domain decomposition scheme

In this paper, we demonstrate a local convergence of an adaptive scalar solver which is practical for strongly diagonal dominant Jacobian problems such as in some systems of nonlinear equations arising from the application of a nonoverlapping domain decomposition method. The method is tested to a nonlinear interface problem of a multichip heat conduction problem. The numerical results show that the method performs slightly better than a Newton-Krylov method.     

A subregion DRBEM formulation for the dynamic analysis of two-dimensional cracks

The dual reciprocity boundary element method employing the step by step time integration technique is developed to analyse two-dimensional dynamic crack problems. In this method the equation of motion is expressed in boundary integral form using elastostatic fundamental solutions. In order to transform the domain integral into an equivalent boundary integral, a general radial basis function is used for the derivation of the particular solutions. The dual reciprocity boundary element method is combined with an efficient subregion boundary element method to overcome the difficulty of a singular system of algebraic equations in crack problems. Dynamic stress intensity factors are calculated using the discontinuous quarter-point elements. Several examples are presented to show the formulation details and to demonstrate the computational efficiency of the method.     

Mixed boundary element solution for three-dimensional convection-diffusion problem with a velocity profile

A mixed boundary element formulation is presented for convection-diffusion problems with a velocity profile. In this formulation the convection-diffusion equation is considered as a nonlinear diffusion equation with inhomogeneous terms in which the convective term is involved additionally, because the spatial distribution of the drift velocity cannot be straightforwardly expressed in boundary integral form. Accordingly, a corresponding boundary integral equation may be described usually in the form of a so-called hybrid-type boundary integral equation.

In the present paper, mixed boundary elements are employed in a discrete model of the original convection-diffusion system. In the mixed element, potentials are approximated linearly, and their normal derivatives to boundaries are assumed constant. A simple iterative scheme is adopted in order to solve hybrid-type mixed boundary element equations. Simple three-dimensional models are dealt with in numerical experiments. The proposed approach gives more accurate and stable solutions compared with constant boundary elements which have been reported.     

Infinite domain potential problems by a new formulation of singular boundary method

This study proposes a new formulation of singular boundary method (SBM) and documents the first attempt to apply this new method to infinite domain potential problems. The essential issue in the SBM-based methods is to evaluate the origin intensity factor. This paper derives a new regularization technique to evaluate the origin intensity factor on the Neumann boundary condition without the need of sample solution and nodes as in the traditional SBM. We also modify the inverse interpolation technique in the traditional SBM to get rid of the perplexing sample nodes in the calculation of the origin intensity factor on the Dirichlet boundary condition. It is noted that this new SBM retains all merits of the traditional SBM being truly meshless, free of integration, mathematically simple, and easy-to-program without the requirement of a fictitious boundary as in the method of fundamental solutions (MFS). We examine the new SBM by the four benchmark infinite domain problems to verify its applicability, stability, and accuracy.     

Hybrid boundary element formulation for acoustic wave propagation

The so‐called hybrid stress boundary element method (HSBEM) is introduced in a frequency domain formulation for the computation of acoustic radiation and scattering in closed and in finite domains. Different from other boundary element formulations, the HSBEM is based on an extended Hellinger‐Reissner variational principle and leads to a Hermitian, frequency‐dependent stiffness equation. Due to this, the method is very well suited for treating fluid structure interaction problems since the effort for the coupling the structure, discretized by a finite elements, and the fluid, discretized by the HSBEM is strongly reduced. To arrive at a boundary integral formulation, the field variables are separated into boundary variables, which are approximated by piecewise polynomial functions, and domain variables, which are approximated by a superposition of singular fundamental solutions weighed by source strength. This approximation cancels the domain integral over the equation of motion in the hybrid principle and leads to a boundary integral formulation, incorporating singular integrals. Comparing to previous results published by the authors, new considerations concerning the interpretation of singular contributions in the stiffness matrix for exterior domain problems are communicated here.     

Oblique derivative problems for second order nonlinear equations of mixed type with degenerate hyperbolic curve

The present paper deals with oblique derivative problems for second order nonlinear equations of mixed type with degenerate hyperbolic curve, which include the Tricomi problem as a special case. Firstly the formulation of the problems for the equations is given, next the representation and estimates of solutions for the above problems are obtained, finally the existence of solutions for the problems is proved by the successive iteration of solutions of the equations and the fixed-point principle. In this paper, we use the complex analytic method, namely the new partial derivative notations, elliptic complex functions in the elliptic domain and hyperbolic complex functions in the hyperbolic domain are introduced, such that the second order equations of mixed type with degenerate curve are reduced to the first order mixed complex equations with singular coefficients, and then the advantage of complex analytic method can be applied.     

Invariant factors of integral quaternion matrices II

In the author's previous paper a diagonal form was found under unimodular equivalence for matrices over the Hurwitz ring of integral quaternions, and uniqueness was established for the norms of certain constituents of the diagonal elements. In the present paper it is shown that the odd, primitive, parts of the all but one of the diagonal elements may be freely chosen provided that the norm constraint is met.     

Factorization in Integral Domains,IV

Let D be an integral domain such that every nonzero nonunit of D is a finite product of irreducible elements. In this article, we introduce and study several unifying concepts for the theory of nonunique factorization in D. They give a new way to measure, in some sense, how far an half-factorial domain (resp., bounded factorization domain, atomic domain) D is from being a UFD (resp., finite factorization domain, Cohen–Kaplansky domain) based on equivalence relations on the set of irreducible elements of D.     

Numerical solution of thermal convection problems using the multidomain boundary element method

The multidomain dual reciprocity method (MD‐DRM) has been effectively applied to the solution of two‐dimensional thermal convection problems where the momentum and energy equations govern the motion of a viscous fluid. In the proposed boundary integral method the domain integrals are transformed into equivalent boundary integrals by the dual reciprocity approach applied in a subdomain basis. On each subregion or domain element the integral representation formulas for the velocity and temperature are applied and discretised using linear continuous boundary elements, and the equations from adjacent subregions are matched by additional continuity conditions. Some examples showing the accuracy, the efficiency and flexibility of the proposed method are presented. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 469–489, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10016     

发展方程初边值问题的高阶差分-边界积分方程法及误差分析

本文结合差分方法与边界积分方程方法,提出并研究了一类新的求解发展型方程初边值问题的高阶差分与边界积分方程耦合数值方法.对于有界区域问题与无界区域问题给出了数值计算格式及其误差的先验估计.     

Boundary element bending analysis of heated elastic plates

Generally, the boundary integral equations for nonlinear and inhomogeneous problems need to be formulated in terms of boundary and domain integrals. The present paper declares the possibility of the transformation of the domain integral term into the boundary integral for the thermal bending analysis of a thin elastic plate,provided that the temperature can be expressed by a linear variation over the thickness.     

Analytical evaluation of the origin intensity factor of time-dependent diffusion fundamental solution for a matrix-free singular boundary method formulation

The singular boundary method (SBM) with the empirical formulas of the origin intensity factors (OIFs) can be effectively used to simulate one- and two-dimensional time-dependent diffusion problems. However, there is no such empirical formula available for determining the OIFs in three-dimensional problems so that the traditional inverse interpolation technique (IIT) has to be employed in three-dimensional case. This paper presents the analytical evaluation formulas to derive the OIFs and thereby overcome the above shortcomings. The proposed new formulation not only has clear theoretical foundations, but also ensures good stability compared with the IIT. Moreover, the present method can effectively simulate three-dimensional diffusion problems. Consequently, our new formulation, most importantly, is matrix-free and fully explicit due to completely avoiding the IIT. As a result, the proposed SBM formulation is mathematically simple, computationally fast and stable, and requiring very low memory since it does not need to solve any algebraic equations. In stark contrast to the boundary element method, the present SBM only requires integration and background grid to calculate the OIFs, while remaining free of integration and mesh for the rest of the calculation. Five benchmark problems are tested to verify the feasibility and accuracy of the new formulation. Numerical results clearly demonstrate the applicability and accuracy of the proposed SBM for solving three-dimensional transient diffusion problems.     

Boundary integral equation method in thermoelasticity part III: uncoupled thermoelasticity

This paper presents conversion of the volume integral of temperature gradients to a surface integral in the integral formulation of boundary value problems of uncoupled thermoelasticity. Particular classes of problems such as thermal stresses, quasi-static problems of uncoupled thermoelasticity, and stationary problems of thermoelasticity are considered. The improved formulation makes numerical computation more accurate and less formidable. The integral formulations in uncoupled thermoelasticity are given in the Laplace transform domain as well as in the time domain.     

椭圆外区域上的自然边界元法

１．引言 二十年来,自然边界元法已在椭圆问题求解方面取得了许多研究成果。它可以直接用来解决圆内（外）区域、扇形区域、球内（外）区域及半平面区域等特殊区域上的椭圆边值问题［１,２,５］,也可以结合有限元法求解一般区域上的椭圆边值问题,例如基于自然边界归化的耦合算法及区域分解算法就是处理断裂区域问题及外问题的一种有效手段［２－４,６］。 人们在设计求解外问题的耦合算法或者区域分解算法时,通常选取圆周或球面作人工边界。但对具有长条型内边界的外问题,以圆周或球面作人工边界显然并非最佳选择,它将会导致大量的…     

A Galerkin boundary element method based on interpolatory Hermite trigonometric wavelets

A Galerkin boundary element method based on interpolatory Hermite trigonometric wavelets is presented for solving 2-D potential problems defined inside or outside of a circular boundary in this paper. In this approach, an equivalent variational form of the corresponding boundary integral equation for the potential problem is used; the trigonometric wavelets are employed as trial and test functions of the variational formulation. The analytical formulae of the matrix entries indicate that most of the matrix entries are naturally zero without any truncation technique and the system matrix is a block diagonal matrix. Each block consists of four circular submatrices. Hence the memory spaces and computational complexity of the system matrix are linear scale. This approach could be easily coupled into domain decomposition method based on variational formulation. Finally, the error estimates of the approximation solutions are given and some test examples are presented.     

Boundary integral equation formulation for generalized thermoelastic diffusion – Analytical aspects

The present work deals with the formulation of the boundary integral equations for the solution of equations under linear theory of generalized thermoelastic diffusion in a three-dimensional Euclidean space. A mixed initial-boundary value problem is considered in the present context and the fundamental solutions of the corresponding coupled differential equations are obtained in the Laplace transform domain by employing the treatment of scalar and vector potential theory. A reciprocal relation of Betti type is established. Then we formulate the boundary integral equations for generalized thermoelastic diffusion on the basis of these fundamental solutions and the reciprocal relation.