三维边界元法高阶元几乎奇异积分半解析法

胡宗军; 牛忠荣; 程长征

doi:10.6052/0459-1879-13-353

三维边界元法高阶元几乎奇异积分半解析法

合肥工业大学土木与水利工程学院, 合肥 230009

基金项目: 国家自然科学基金资助项目（11272111，11372094）.

详细信息

作者简介:
胡宗军，副教授，主要研究方向：计算力学.E-mail：huzongjun1975@163.com

中图分类号: O343.2
计量
- 文章访问数: 1243
- HTML全文浏览量: 98
- PDF下载量: 994
出版历程
- 收稿日期: 2013-10-21
- 修回日期: 2013-12-16
- 刊出日期: 2014-05-17

A NEW SEMI-ANALYTIC ALGORITHM OF NEARLY SINGULAR INTEGRALS IN HIGH ORDER BOUNDARY ELEMENT ANALYSIS OF 3D POTENTIAL

School of Civil Engineering, Hefei University of Technology, Hefei 230009, China

Funds: The project was supported by the National Natural Science Foundation of China (11272111, 11372094).

摘要

摘要: 分析了三维边界元法高阶曲面单元几何特征，定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系，构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数，分离出积分核中主导的奇异函数部分，将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算，奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ，对ρ积分建立解析计算列式，对θ积分应用常规Gauss数值积分计算，从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例，采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度，并与线性单元正则化算法计算结果对比，结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效.
- 三维 /
- 位势 /
- 边界元法 /
- 高阶单元 /
- 几乎奇异积分 /
- 半解析法
Abstract: By analyzing the geometric feature of 8-noded quadrilateral element in three dimensional boundary element method (3D BEM), the relative distance is first defined as the approach degree from a source point to the high order surface element. And then a local polar coordinate ρθ is built which origin point is the project point of the source point on the element surface. The approximate singular kernel function is constructed corresponding to the nearly singular integral on high order surface elements in 3D potential BEM by a series of deduction, which has the same singularity as the nearly singular kernel function. The leading singular part is separated by subtracting the approximate kernel function from the original kernel function. Thus the nearly singular surface integrals on high order elements are transformed into the sum of both the non-singular integrals and singular integrals. The former can be efficiently computed by the Gaussian quadrature. The integral variables ρ and θ of the later are separated in the local polar coordinate. The singular surface integrals with respect to polar variable ρ are firstly expressed by the analytic formulations. Then the surface integrals are transformed into the line integrals with respect to variable θ, which can be evaluated by the Gaussian quadrature. Consequently, the new semi-analytic algorithm is established to calculate the nearly strongly and hyper-singular surface integrals on high order element in 3D potential BEM. Some numerical examples about the high order BE analysis for 3D heat conduction problems are given to demonstrate the efficiency and accuracy of the present semi-analytic algorithm. In comparison with the published regularization algorithm which is applied to calculating the nearly singular integrals on 3-noded triangular element, the present semi-analytic algorithm with 8-noded quadrilateral element can evaluate the potentials and potential gradients of inner points more close to the boundary. Moreover, the semi-analytic algorithm can be applied to more efficiently analyze thin structures in 3D potentials.
- 3D /
- potential /
- BEM /
- high order element /
- nearly singular integrals /
- semi-analytic method

HTML全文

参考文献(16)

Sladek V, Sladek J, Tanaka M. Numerical integration of logarithmic and nearly logarithmic singularity in BEMs. Applied Mathematical Modelling, 2001, 25: 901-922

Padhi GS, Shenoi RA, Moy SS, et al. Analytic integration of kernel shape function product integrals in the boundary element method. Computers and Structures, 2001, 79: 1325-1333

Shhiah YC, Shi YX. Heat conduction across thermal barrier coatings of anisotropic substrates. International Communications in Heat and Mass Transfer, 2006, 33: 827-835

Niu ZR, Cheng CZ, Zhou HL, et al. Analytic formulations for calculating nearly singular integrals in two-dimensional BEM. Engineering Analysis with Boundary Elements, 2007, 31: 949-964

牛忠荣,王秀喜,周焕林. 边界元法计算近边界点参量的一个通用算法. 力学学报,2001,33(2):275-283 (Niu Zhongrong, Wang Xiuxi, Zhou Huanlin. A general algorithm for calculating the quantities at interior points close to the boundary by BEM. Acta Mechanica Sinica, 2001, 33(2): 275-283 (in Chinese))

张耀明,孙翠莲,谷岩. 边界积分方程中近奇异积分计算的一种变量替换法. 力学学报,2008,40(2):207-214 (Zhang Yaoming, Sun Cuilian, Gu Yan. The evaluation of nearly singular integrals in the boundary integral equations with variable transformation. Acta Mechanica Sinica, 2008,40(2): 207-214 (in Chinese))

Dehghan M, Hosseinzadeh H. Calculation of 2D singular and near singular integrals of boundary elements method on complex space C. Applied Mathematical Modelling, 2012, 36: 545-560

Dehghan M, Hosseinzadeh H. Obtaining the upper bound of discretization error and critical boundary integrals of circular arc element method. Mathematical and Computer Modelling, 2012, 55: 517-529

Milroy J, Hinduja S, Davey K. The 3-D thermoelastic boundary element method: Analytical integration for linear isoparametric triangular elements. Appl Math Model, 1997, 21: 763-782

Liu YJ. Analysis of shell-like structures by the boundary element method based on 3-D elasticity: Formulation and verification. International Journal for Numerical Methods in Engineering, 1998, 41(3): 541-558 3.0.CO;2-K" target=_blank>

牛忠荣,王秀喜,周焕林. 三维边界元法中几乎奇异积分的正则化算法. 力学学报,2004,26(1):49-56 (Niu Zhongrong, Wang Xiuxi, Zhou Huanlin. A regularization algorithm for the nearly singular integrals in 3-D BEM. Acta Mechanica Sinica, 2004, 26(1): 49-56 (in Chinese))

Pasternak I, Sulym H. Self-regular integral equation for axisymmetric elasticity. Engineering Analysis with Boundary Elements, 2009, 33: 1001-1004

Nintcheu FS. Semi-analytic treatment of nearly-singular Galerkin surface integral. Applied Numerical Mathematics, 2010, 60: 974-993

Qin XY, Zhang JM, Xie GZ, et al. A general algorithm for the numerical evaluation of nearly singular integrals on 3D boundary element. Journal of Computational and Applied Mathematics, 2011, 235: 4174-4186

Xie GZ, Zhou FL, Zhang JM, et al. New variable transformations for evaluating nearly singular integrals in 3D boundary element method. Engineering Analysis with Boundary Elements, 2013, 37: 1169-1178

Barbara M, Peter R, David E. A new method for the numerical evaluation of nearly singular integrals on triangular elements in the 3D boundary element method. Journal of Computational and Applied Mathematics, 2013, 245: 148-161