作　　者： (孙锐）; (胡宗军）; (牛忠荣）; (周焕林）;

机构地区： 合肥工业大学土木与水利工程学院,合肥230009

出　　处： 《中国科学：物理学、力学、天文学》 2017年第9期18-27,共10页

摘　　要： 针对边界元法中几乎奇异积分计算难题,本文提出一种基于6节点三角形等参数单元的三维高阶单元半解析算法.通过对三维声场基本解中的三角函数进行T a y l o r级数展开,分离出基本解中的奇异积分项.根据单元的几何特性,构造出与奇异积分核函数具有相同奇异性的近似奇异核函数,对奇异积分项应用扣除法,将奇异积分核函数分为规则核函数和近似奇异核函数两项.规则核函数积分无奇异性,应用常规G a u s s数值积分就能够准确计算;近似奇异核函数积分由导出的半解析公式计算,即在局部极坐标系ρθ下分离积分变量,导出对变量ρ积分的解析计算列式,应用常规G a u s s数值积分计算变量θ积分,从而建立一种三维声场边界元法几乎奇异积分的半解析算法.算例结果表明,本文高阶单元半解析算法比双线性元算法更加有效且算法稳定,能够有效、准确地计算距离单元非常近的近边界点处的声压. For the nearly singular integral of three-dimensional acoustic boundary element method （BEM）, based on the 6-noded triangular isoparametric element, a new semi-analytical algorithm of three-dimensional high order element is proposed in this paper. Using Taylor expansion of trigonometric functions in the three dimensional acoustic fundamental solutions, the singular part of the fundamental solutions is separated. Based on the geometric characteristics of the 6-noded triangular element, an approximate singular kernel function is constructed which has the same singularity as singular integral kernel function. Subtracting the approximate kernel function from the kernel function of the singular integral, the latter is decomposed into a regular kernel function and an approximate singular kernel function. The integral of the regular kernel function can be calculated accurately by using the conventional Gauss numerical quadrature. The integral of the new singular part is calculated by the semi-analytic formula derived in this paper. In the surface of the integral element, the local coordinate systempθ is established and the approximate singular integral is transformed into the integrals of variables p and θ which are already separated in pO system. The integral with respect to polar variable p is expressed by the analytic formulations first. Then the new singular integral which is a surface integral is transformed into the line integral with respect to variable θ, which can be evaluated by the Gaussian quadrature. Consequently, the new semi-analytic algorithm is established to calculate the nearly singular surface integrals in 3D acoustic BEM. Some examples are given in the last part of this paper to show the accuracy and the effectiveness of the present algorithm. The computed results demonstrate that the semi-analytic algorithm with high order element presented in this paper is more effective than linear regularization BEM to solve nearly singular integrals for 3D acoustic BEM.

关 键 词： 边界元法 声辐射 几乎奇异积分 半解析算法 节点三角形单元