轴对称Poisson方程的Trefftz有限元解法

将径向基函数应用到一类轴对称Poisson方程的数值求解中,提出了一种Trefftz有限元计算格式.非0右端项将问题的特解引入Trefftz单元域内场,致使单元刚度方程涉及区域积分.利用径向基函数对特解近似处理,可消除区域积分,从而保持Trefftz有限元法只含边界积分的优势.为获得特解,选取求解域内所有单元的节点和形心作为基本插值点,而在求解域之外构造一个虚拟边界,在其上布置一定数目的虚拟点作为额外插值点.数值算例验证了该方法的有效性和可行性.     

用于解析函数复分析的共轭边界元法

由2个共轭的实调和函数构建1个复解析函数,其复分析在应用数学和力学领域具有重要的作用.提出了一个加权残数方程组,证明了该方程组为2个共轭函数的域内控制方程、边界条件和边界上Cauchy Riemann(柯西-黎曼)条件的近似解,等效为复解析函数的逼近方程.在离散空间中,由该加权残数方程分别推导出2个位势问题的直接边界积分方程和1个表示Cauchy-Riemann条件的有限差分方程,随后解决了弱奇异线性方程组的求解难题,并提出用Cauchy积分公式求内点值的方法,从而建立了一种用于复分析的常单元共轭边界元法.最后,用3个算例证明了所提出方法适用于域内或域外的幂函数、指数函数或对数函数形式的解析函数,而且其误差与2维位势问题是同等量级的.     

双参数地基上板弯曲问题的边界积分方程

本文应用广义函数的Fourrier积分变换导出了双参数地基上板弯曲问题的基本解,并将基本解展成一致收敛的级数形式.在此基础上,应用广义Rayleigh-Green公式建立了适用于任意形状、任意边界条件情形的两个边界积分方程,为边界元法在这一问题中的应用提供了理论基础.     

含孔薄板弯曲波动的双互易边界元法

采用双互易边界元法对开孔无限大薄板弹性波的散射与动应力集中问题进行理论分析和数值计算．基于功的互等定理,采用静力基本解建立了薄板弯曲波动问题的双互易边界积分方程．作为数值算例,计算了圆孔附近的动应力集中系数,通过与已有结果进行比较,表明该方法简单有效并能够保证计算精度．     

基于径向基函数插值的积分方程求解

将径向基函数(radial basis function,RBF)插值引入积分方程的求解中,具体将待求函数表示为RBF的线性组合,再通过配点法将积分方程离散为线性或非线性方程组,求得权系数后给出待求函数的近似表示.论文选用的RBF是插值性能优异的多重二次曲面(multiquadric,MQ)函数,能在较少节点下取得较高的近似精度;而且RBF定义为距离的函数,在三维或高维插值时仅需改变距离公式,因而便于推广到高维积分方程求解中.在RBF插值矩阵的构造中,元素的积分计算分别通过高斯积分或基于区域剖分的数值求积完成,实现了一维、二维下Fredholm和Volterra方程的求解.算例结果表明:论文方法具有实施方便和精度较高的优点,是一种适合积分方程求解的新方法.     

正交各向异性弹性力学平面问题的样条虚边界元法

采用域外奇点技术并根据问题的边界条件,建立了正交各向异性弹性力学平面问题的非奇异虚边界积分方程,然后采用性态优越的B样条函数去逼近未知虚荷载函数,并采用性能稳定的最小二乘边界子段法去消除边界余量,据此获得积分方程的数值解.数值算例表明:该方法具有相当高的精度和良好的数值稳定性,且计算工作量少.文中引言部分还对域外奇点法的发展作了系统的评述.     

求解平片裂纹问题的有限部积分与边界元法

本文利用位移的Somigliana公式和有限部积分的概念,导出了求解三维弹性力学中的任意形状平片裂纹问题的超奇异积分方程组,进而联合使用有限部积分法与边界元法对所得方程建立了数值法.为验证本文的方法,计算了若干数值例子的裂纹面的位移间断及裂纹前沿的应力强度因子,它们与理论值相比符合很好.     

用径向基函数插值解自共轭椭圆型方程

本文讨论用MQ作为插值的径向基函数,对自共轭椭圆型方程进行插值,证明了插值系数的唯一性,并用投影法证明了用径向基函数解自共轭椭圆型方程的收敛性.     

改进的无奇异局部边界积分方程方法

在局部边界积分方程方法中,当源节点位于分析域的整体边界上时,局部边界积分将出现奇异积分问题,这些奇异积分需要做特别的处理．为此,提出了对域内节点采用局部积分方程,而对边界节点直接采用移动最小二乘近似函数引入边界条件来解决奇异积分问题,这同时也解决了对积分边界进行插值引入近似误差的问题．作为应用和数值实验,对Laplace方程和Helmholtz方程问题进行了分析,取得了很好的数值结果．进而,在Helmholtz方程求解中,采用了含波解信息的修正基函数来代替单项式基函数进行近似．数值结果显示,这样处理是简单高效的,在高波数声传播问题的求解中非常具有前景．     

简支梯形底扁球壳自由振动问题的准Green函数方法

以简支梯形底扁球壳的自由振动问题为例,详细阐明了准Green函数方法的思想.即利用问题的基本解和边界方程构造一个准Green函数,此函数满足了问题的齐次边界条件,采用Green公式,将简支梯形底扁球壳自由振动问题的振形控制微分方程化为两个耦合的第二类Fredholm积分方程.边界方程有多种选择,在选定一种边界方程的基础上,可以通过建立一个新的边界方程来表示问题的边界,以克服积分核的奇异性.最后由积分方程的离散化方程组有非平凡解的条件,求得固有频率.数值结果表明,该方法具有较高的精度.     

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.     

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     

Multi-domain mass conservative dual reciprocity method for the solution of the non-Newtonian Stokes equations

The dual reciprocity method (DRM) is a technique to transform the domain integrals that appear in the boundary element method into equivalent boundary integrals. In this approach the non-linear terms are approximated by an interpolation applied to the non-Newtonian stress tensor for an inelastic fluid. In the present paper we introduce a radial basis function interpolation scheme for the velocity field that satisfies the continuity equation (mass conservative interpolation). The proposed method performs better than the classical interpolation used in the DRM approach to represent such field. The new scheme together with a sub-domain variation of the DRM yields a more accurate solution for inelastic non-Newtonian problems.     

奇异杂交边界点方法求解泊松方程

As a boundary-type meshless method,the singular hybrid boundary node method(SHBNM)is based on the modified variational principle and the moving least square(MLS)approximation,so it has the advantages of both boundary element method(BEM)and meshless method.In this paper,the dual reciprocity method(DRM)is combined with SHBNM to solve Poisson equation in which the solution is divided into particular solution and general solution.The general solution is achieved by means of SHBNM,and the particular solution is approximated by using the radial basis function(RBF).Only randomly distributed nodes on the bounding surface of the domain are required and it doesn't need extra equations to compute internal parameters in the domain.The postprocess is very simple.Numerical examples for the solution of Poisson equation show that high convergence rates and high accuracy with a small node number are achievable.     

A numerical solution for advection‐diffusion equation using dual reciprocity method

This article describes a numerical method based on the boundary integral equation and dual reciprocity method(DRM) for solving the one‐dimensional advection‐diffusion equations. The concept of DRM is used to convert the domain integral to the boundary that leads to an integration free method. The time derivative is approximated by the time‐stepping method. Numerical results are presented for some problems to demonstrate the usefulness and accuracy of the new approach. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Boundary element analysis of uncoupled transient thermo-elastic problems with time- and space-dependent heat sources

A boundary element method (BEM) for the analysis of two- and three-dimensional uncoupled transient thermo-elastic problems involving time- and space-dependent heat sources is presented. The domain integrals are efficiently treated using the Cartesian transformation and the radial integration methods without considering any internal cells. Similar to the dual reciprocity method (DRM), some internal points without any connectivity are considered; however, in contrast to the DRM, any arbitrary mesh-free interpolation method can be used in the present formulation. There is no need to find any particular solutions and the shape functions in the mesh-free interpolation method can be arbitrary and sufficiently complicated. Unlike the DRM, the generated system of equations contains the unknowns only on the boundary. After finding the primary unknowns on the boundary, the temperature, displacement, and stress components at all internal points can directly be found without solving any system of equations. Three examples with different forms of heat sources are presented to demonstrate the efficiency and accuracy of the proposed method. Although the proposed BEM is mathematically more complicated than domain methods, such as the finite element method (FEM), it is more efficient from a modelling viewpoint since only the surface mesh has to be generated in the presented method.     

Multigrid technique for biharmonic interpolation with application to dual and multiple reciprocity method

A biharmonic-type interpolation method is presented to solve 2D and 3D scattered data interpolation problems. Unlike the methods based on radial basis functions, which produce a large linear system of equations with fully populated and often non-selfadjoint and ill-conditioned matrix, the presented method converts the interpolation problem to the solution of the biharmonic equation supplied with some non-usual boundary conditions at the interpolation points. To solve the biharmonic equation, fast multigrid techniques can be applied which are based on a non-uniform, non-equidistant but Cartesian grid generated by the quadtree/octtree algorithm. The biharmonic interpolation technique is applied to the multiple and dual reciprocity method of the BEM to convert domain integrals to the boundary. This makes it possible to significantly reduce the computational cost of the evaluation of the appearing domain integrals as well as the memory requirement of the procedure. The resulting method can be considered as a special grid-free technique, since it requires no domain discretisation. This revised version was published online in June 2006 with corrections to the Cover Date.     

The comparison between the DRBEM and DQM solution of nonlinear reaction–diffusion equation

In this study, both the dual reciprocity boundary element method and the differential quadrature method are used to discretize spatially, initial and boundary value problems defined by single and system of nonlinear reaction–diffusion equations. The aim is to compare boundary only and a domain discretization method in terms of accuracy of solutions and computational cost. As the time integration scheme, the finite element method is used achieving solution in terms of time block with considerably large time steps. The comparison between the dual reciprocity boundary element method and the differential quadrature method solutions are made on some test problems. The results show that both methods achieve almost the same accuracy when they are combined with finite element method time discretization. However, as a method providing very good accuracy with considerably small number of grid points differential quadrature method is preferrable.