三维位势问题边界元法中几乎奇异积分的正则化-计算物理.PDFVIP

计 算 物 理 第22卷第6期 v01．22．No．6 PHYSICS 2005年11月 CHINESEJOURNALOFCOMPUTATIONAL N。．j2005 [文章编号】1001—246X(2005)06．0501．06 三维位势问题边界元法中几乎奇异积分的正则化 周焕林1， 牛忠荣1， 王秀喜2 (1．合肥工业大学工程力学系，安徽合肥230009； 2．中国科学技术大学中科院材料力学行为和设计重点实验室，安徽合肥230026) [摘 要】 采用一种半解析正则化算法，计算了三维位势问题边界元法中近边界点的几乎强奇异和几乎超奇异 面积分．该算法适用于三角形线性等参元．对高次单元将其细分为几个三节点三角形单元即可应用该算法．由 于几乎奇异性，与内点邻近的单元上的积分，采用半解析正则化积分算法计算；而远处单元的积分仍保持常规 高斯积分．对三维热传导算例，计算了近边界点的温度和热流．数值结果证明了该算法的有效性和精确性． [关键词】 边界元法；几乎奇异积分；正则化；三维位势问题；热传导 [中图分类号】0241．83；0343．2[文献标识码】 A 0 引言 边界元法作为一种重要的数值方法在许多领域得到了广泛应用n】．奇异积分和几乎奇异积分的处理一 述．余德浩H3较早提出了在自然边界元法中求解一类超奇异积分的级数展开法．董春迎等舾1对超奇异积分的 计算作了简要综述． 当源点为域内近边界点时，其物理参量的计算通常由于几乎奇异积分导致失效，这种现象称为边界层效 应．当研究区域含有狭长细薄几何形状时，常规边界元法同样面临几乎奇异积分处理的困难，此类问题称为 薄体问题．近来研究表明，边界元法能够有效处理薄体结构问题【6“]，而且相对有限元法，边界元法所用单元 数目少，计算精确度高． 许多研究者做了大量工作处理二维或三维问题中的几乎奇异积分，文[8]对此做了全面综述．对三维问 性等参元上强奇异面积分的解析结果．Devey等u¨在研究稳态弹性动力学问题边界积分方程中一类奇异积 分时，采用三角形单元将奇异积分核函数做Taylor级数展开，将级数首项的奇异部分转化为线积分计算，非 奇异面积分做数值积分．Ma等n23将源点到场点的距离函数r采用泰勒展开，通过距离函数变换法，将几乎 奇异积分转化为非奇异积分，采用常规高斯积分计算． 文[13]对三维问题提出了一种计算近边界点参量的半解析正则化算法，该算法将单元面积分核函数在 局部极坐标系矽下表示，采用分部积分，获得对变量』D做单层积分的解析表达式．几乎强奇异和超奇异面积 分被转化为沿单元围道上一系列线积分，而高斯数值积分能够有效计算这些线积分，从而半解析地计算出了 三维弹性力学问题近边界点参量的几乎强奇异和几乎超奇异积分．本文将此算法拓展用于三维位势问题，计 算了近边界点的位势及其导数． 1三维位势问题内点位势及其导数 考虑三维位势问题，区域内点J，的位势及其导数边界积分方程可写成 [收稿日期]2004—07—29“修回日期】2005—01—28 [基金项目]国家自然科学基金资助项目 【作者简介】周焕林(1973一)，男，安徽宿松，博士，副教授，主要从事计算力学边界元法的研究工作，合肥工业大学7号信箱230009 计 算 物 理 第22卷 “(J，)：lu。(x，y)q(x)dP—Ig。(工，j，)u(x)dr， (1) (2) 其中g(工)为场点工的法向位势梯度，-(x)2场A 基本解对边界外法线的梯度．Y。为源点Y的坐标．离散积分计算式(1)，(2)用到的基本核函数有