各向异性位势问题边界元法中几乎奇异积分的解析算法

导出了一种解析积分算法,精确计算了二维各向异性位势问题边界元法中近边界点的几乎奇异积分。对线性单元,几乎奇异积分可用解析公式直接计算。对二次单元,可将其细分为几个线性单元,采用该解析公式间接近似计算。当内点离积分单元较远时,仍然保持常规高斯数值积分模式;而当内点离其较近时,高斯积分结果失效,采用该解析积分取代高斯数值积分。数值算例证明了该算法的有效性和精确性。二次元比线性元计算结果更精确。     

二维位势边界元法高阶单元几乎奇异积分半解析算法

准确计算几乎奇异积分是边界元法难题之一。目前,对于一般的高阶单元的几乎奇异积分尚缺乏通用高效的计算方法。本文在单元局部坐标系中表征了二维高阶单元的几何特征,提出了源点相对高阶单元的接近度概念。针对二维位势边界元法的3节点二次等参单元,构造出与单元积分核具有相同几乎奇异性的近似奇异核函数。从二维位势几乎奇异积分单元积分核中扣除近似奇异核函数,把几乎奇异积分项转换为规则积分和奇异积分两部分之和,规则积分部分用常规Gauss数值积分计算,奇异积分部分由导出的解析公式计算,从而建立了二维位势问题高阶单元几乎强奇异和超奇异积分的半解析算法。算例结果表明了本文半解析算法的有效性和计算精度。     

Similar theory and the applicable ranges of one-equation and two-equation models

Here, we are discussing the applicable ranges of one-equation and two-equation turbulence models with the viewpoint of similar theory. The criterions of determining the applicable ranges of these models are given.     

一种处理弹性动力边界元法中奇异积分的方法

利用边界元法求解瞬态弹性动力学问题时,时域基本解函数的分段连续性和奇异性为该问题的求解带来很大的困难。为了解决时域基本解中的奇异性问题,本文依据柯西主值的定义,对经过时间解析积分之后的时域基本解进行奇异值分解,将其分成奇异和正则积分两部分;其中正则部分可通过采用常规高斯积分方法来计算,而奇异部分具有简单的形式,可以利用解析积分计算。经过上述操作之后,就可以达到直接消除时域基本解中奇异积分的目的。和传统方法相比,本文方法并不依赖静力学基本解来消除奇异性,是一种直接求解方法。最后给定两个数值算例来验证本文提出方法的正确性和可行性,结果表明使用本文算法可以解决弹性动力学边界积分方程中的奇异性问题。     

基于坐标变换精确计算近奇异积分的双向sinh变换法

提出一种与坐标变换相结合的双向sinh变换方法用于精确计算近奇异积分.首先利用坐标变换法对近奇异积分进行双向分离,再针对分离出的两个方向的积分形式,根据复变函数极点理论,构造双向sinh变换.与传统结合环向变换和极坐标变换的单向sinh变换以及迭代sinh变换方法相比,双向sinh变换法的计算精度更高,并且最近点的位置...     

An improved exponential transformation for nearly singular boundary element integrals in elasticity problems

This paper presents an improved exponential transformation for nearly singular boundary element integrals in elasticity problems. The new transformation is less sensitive to the position of the projection point compared with the original transformation. In our work, the conventional distance function is modified into a new form in the polar coordinate system. Based on the refined distance function, an improved exponential transformation is proposed in the polar coordinate system. Moreover, to perform integrations on irregular elements, an adaptive integration scheme considering both the element shape and the projection point associated with the improved transformation is proposed. Furthermore, when the projection point is located outside the integration element, another nearest point is introduced to subdivide the integration elements into triangular or quadrilateral patches of fine shapes. Numerical examples are presented to verify the proposed method. Results demonstrate the accuracy and efficiency of our method.     

三维变系数热传导问题边界元分析中几乎奇异积分计算

在边界积分的数值计算过程中,当源点离积分单元很近时,边界积分就会具有几乎奇异性,此时不能直接用高斯数值积分公式计算几乎奇异积分。本文以三维非均质热传导问题为例,介绍了一种计算几乎奇异边界积分的新方法。首先,采用Newton-Raphson迭代算法确定积分单元上离源点最近的点;然后,将积分单元上任意一点的坐标在最近点处展开成泰勒级数,并计算源点到积分单元任意点的距离;最后,将距离函数代入几乎奇异边界积分中,并运用指数变换方法导出积分单元上几乎奇异积分的计算公式。文中给出了两个非均质热传导问题的算例来验证所述方法的正确性、有效性和稳定性。     

A higher-order boundary element method for three-dimensional potential problems

A method of eliminating the singularities involved in boundary element methods for three-dimensional potential problems is presented and the non-singular expressions of integrals on an element on which the singular point is situated are given for linear and quadratic interpolation functions. Numerical examples are compared with analytical solutions to show that the higher-order interpolations have better precision.     

求解位势问题的虚边界元法

本文提出了求解位势问题的虚边界元法，建立了位势问题的虚边界元的离散方程式，推导了离散化求系数的积分解析式。该方法与传统边界元法相比具有不存在奇异积分和边界附近精度较高等优点，可用来计算真空静电场、稳定温度场、流体绕流、介质中的渗流等各类位势问题。大量算例均获得了满意的结果。     

边界元法分析狭长体结构

针对边界元法分析狭长结构时遇到的几乎奇异积分难以计算的困难，将几乎奇异积分划分为两种类型，分别通过分部积分交换把引起积分几乎奇异的参量移至积分号之外，从而建立了一个新的正则化算法，解决了边界积分方程中几乎奇异积分的计算难题。文中用边界元法计算了弹性力学平面问题的狭长结构，算例证明了本法的有效性。     

Electro-mechanical coupling analysis of MEMS structures by boundary element method

In this paper, we present the applications of Boundary Element Method (BEM) to simulate the electro-mechanical coupling responses of Micro-Electro-Mechanical systems (MEMS). The algorithm is programmed in our research group based on BEM modeling for electrostatics and elastostatics. Good agreement is shown while the simulation results of the pull-in voltages are compared with the theoretical/experimental ones for some examples. The project supported by the 973 Program (G199033108) and the national Natural Science Foundation of China (10125211)     

高阶边界元奇异积分的一种通用高效计算方法

针对三维边界元法中曲面单元上的(弱、强、超)奇异积分提出了一种通用高效的计算方法。经极坐标变换,将奇异积分转化为常规积分;采用数值方法计算Cauchy主值积分和Hadamard有限项积分系数;引入保角变换和反曲变换消除因单元畸形或因积分点靠近单元边界而引起的周向积分奇异性。该方法可以统一处理(弱、强、超)奇异积分,并且只需要知道核函数的奇异阶数和少数几个点上的被积函数值,不依赖于积分和函数的具体选取;所需的积分点少,精度高,并且受单元畸形程度影响较小,稳定性好。采用该方法计算了声学和弹性力学中的典型奇异积分,并结合二阶Nystrm方法求解了弹性力学的边界积分方程,验证了方法的高精度和高效性。本文数值积分程序可向作者索取。     

NOVEL REGULARIZED BOUNDARY INTEGRAL EQUATIONS FOR POTENTIAL PLANE PROBLEMS

The universal practices have been centralizing on the research of regulariza-tion to the direct boundary integal equations (DBIEs). The character is elimination of singularities by using the simple solutions. However, up to now the research of regular ization to the first kind integral equations for plane potential problems has never been found in previous literatures. The presentation is mainly devoted to the research on the regularization of the singular boundary integral equations with indirect unknowns. A novel view and idea is presented herein, in which the regularized boundary integral equations with indirect unknowns without including the Cauchy principal value (CPV) and Hadamard-finite-part (HFP) integrals are established for the plane potential problems. With some numerical results, it is shown that the better accuracy and higher efficiency, especially on the boundary, can be achieved by the present system.     

High order multiplication perturbation method for singular perturbation problems

This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations （ODEs） are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.     

自由面势流问题边界元法的数值色散误差

以连续及离散Fourier分析研究自由面势流问题边界元法的数值色散误差,并从理论上探讨有关计算中数值色散误差的改善问题.研究表明:对于该问题的数值色散误差而言,重要的在于以问题相应的离散算子考察计及各种数值手段后的总体色散误差,而非仅考虑该数值手段自身的数值色散误差大小.高阶面元、自由面域外奇点或适当的耦合方法是降低有关问题算子总体色散误差的较好选择.     

模拟相变问题的精细积分边界元法

将精细积分边界元法和界面追踪法相结合求解相变问题。因为边界元法只需要将待求解空间域的边界离散,方便连续追踪移动界面位置和重构网格,所以边界元法适合应用于移动边界问题的模拟。首先,利用精细积分边界元法在固相区域和液相区域分别求解相应的瞬态热传导控制方程,从而求得温度场和边界热流密度。然后,根据固-液相变界面上的能量平衡方程,利用热流密度求得相变界面的移动速度,再采用界面追踪法预测移动相变界面的位置变化。最后,给出了几个数值算例,并通过与参考解的对比验证本文方法的准确性。     

三维非规则非均匀边界元网格的简便的高精度算法

对三维直接边界元中一阶奇异积分、一阶近奇异积分以及非奇异积分进行统一处理，给出了一种提高积分计算精度的简便有效的方法，对非规则非均匀边界元网格可获得比一般方法高得多的计算精度，非常适合边界形状比较复杂的三维实际问题的边界元分析．     

ELASTIC ANALYSIS OF ORTHOTROPIC PLANE PROBLEMS BY THE SPLINE FICTITIOUS BOUNDARY ELEMENT METHOD

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｈｅｃｏｎｖｅｎｔｉｏｎａｌｂｏｕｎｄａｒｙｅｌｅｍｅｎｔｍｅｔｈｏｄ (ＢＥＭ) ,ｓｉｎｇｕｌａｒｐｏｉｎｔｓａｒｅｄｉｓｔｒｉｂｕｔｅｄａｌｏｎｇｔｈｅｂｏｕｎｄａｒｙｏｆｔｈｅｄｏｍａｉｎｕｎｄｅｒｓｔｕｄｙ ,ｗｈｉｃｈｌｅａｄｓｔｏｓｉｎｇｕｌａｒｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎｓｄｕｅｔｏｔｈｅｓｉｎｇｕｌａｒｉｔｙｏｆｆｕｎｄａｍｅｎｔａｌｓｏｌｕｔｉｏｎｓ .Ｔｈｅｒｅｆｏｒｅ ,ｓｉｎｇｕｌａｒｉｎｔｅｇｒａｔｉｏｎｍｕｓｔｂｅｈａｎｄｌｅｄｗｈｅｎｓｏｌ…     

Helmholtz边界积分方程中奇异积分间接求解方法

提出了间接求解传统Helmholtz边界积分方程CBIE的强奇异积分和自由项系数,以及Burton-Miller边界积分方程BMBIE中的超强奇异积分的特解法。对于声场的内域问题,给出了满足Helmholtz控制方程的特解,间接求出了CBIE中的强奇异积分和自由项系数。对于声场外域对应的BMBIE中的超强奇异积分,按Guiggiani方法计算其柯西主值积分需要进行泰勒级数展开的高阶近似,公式繁复,实施困难。本文给出了满足Helmholtz控制方程和Sommerfeld散射条件的特解,提出了间接求出超强奇异积分的方法。推导了轴对称结构外场问题的强奇异积分中的柯西主值积分表达式,并通过轴对称问题算例证明了本文方法的高效性。数值结果表明,对于内域问题,采用本文特解法的计算结果优于直接求解强奇异积分和自由项系数的结果,且本文的特解法可避免针对具体几何信息计算自由项系数,因而具有更好的适用性。对于外域问题,两者精度相当,但本文的特解法可避免对核函数进行高阶泰勒级数展开,更易于数值实施。     

Coupling of high order multiplication perturbation method and reduction method for variable coefcient singular perturbation problems

Based on the precise integration method(PIM), a coupling technique of the high order multiplication perturbation method(HOMPM) and the reduction method is proposed to solve variable coefcient singularly perturbed two-point boundary value problems(TPBVPs) with one boundary layer. First, the inhomogeneous ordinary diferential equations(ODEs) are transformed into the homogeneous ODEs by variable coefcient dimensional expansion. Then, the whole interval is divided evenly, and the transfer matrix in each sub-interval is worked out through the HOMPM. Finally, a group of algebraic equations are given based on the relationship between the neighboring sub-intervals, which are solved by the reduction method. Numerical results show that the present method is highly efcient.