二维弱奇异积分高精度数值求积公式的构造.pdfVIP

第37卷 第4期 东 华 理 工 大 学 学 报(自然科学版) VoL37 No．4 2014年 12月 JOURNAL OF EAST CHINA INSTITUTE OF TECHNOLOGY DeC．2014 doi：10．3969／j．issn．1674-3504．2014．04．016 二维弱奇异积分高精度数值求积公式的构造 曾 光 ，黄 晋2， 雷 莉 ， 宁德圣 (1．东华理工大学理学院，江西 南昌 330013；2．电子科技大学数学科学学院，四川成都 611731) 摘 要：在欧拉一麦克劳林展开式和一维弱奇异积分的求积公式的基础上，推导出了二维弱奇异积分的求积公式及其误 差的渐进展开式。此类求积公式只需赋值，不需计算二重积分，故计算量小。利用这类积分公式进行计算可以得到十分 精确的结果，使得收敛阶大为提高，为讨论更为复杂地多维弱奇异积分方程奠定了基础。 关键词：弱奇异积分；求积公式；高精度；欧拉一麦克劳林展开式 中图分类号：O186 文献标识码 ：A 文章编号：1674-3504(2014)04-0447-04 曾光。黄晋，雷莉，等．2014．二维弱奇异积分高精度数值求积公式的构造[J]．东华理工大学学报：自然科学版，37(4)：447— 450． ZengGuang，HuangJin，LeiLi，eta1．2014．ANewConstructionofHighAccurateNumericalQuadratureFormulasfor2DWeak SingularIntegral[J]．JournalofEastChinaInstituteofTechnology(NaturalScience)，37(4)：447-450． 考虑二维弱奇异积分 展开式和后验误差估计。此外还有Duffy变换消去 法和奇异减方法：这是由Atkinson(1989，1997)提 I(s,t)：．J dxdy0 出，Duffy变换消去法是利用该变换把区间的奇异 1，a 6，ctd． (1) 性转化到端点部分进行计算，不过同样存在端点的 这类积分来源于大量工程问题诸如弹性力学、 超奇异性计算；奇异减方法的不足：虽然奇异消去 断裂力学等问题中。其中寻找有效，精度高的数值 了，但却产生了一个高阶函数项和一个已知的奇异 方法求解弱奇异积分和弱奇异积分方程成为计算 部分，仍然存在着一个奇异部分的计算。 数学领域的研究热点之一 (阮周生等，2010；Lifanov 为了克服以上方法的不足，本文主要工作是基 eta1．，2004；Sidieta1．，1988；Lynesseta1．，1967)。 于欧拉一麦克劳林展开式和 Sidi推导出的一维弱 对于各类一维奇异积分与奇异积分方程的计算已 奇异积分的求积公式基础上，推导出了二维弱奇异 有很多学者做了大量的工作，主要包括配置法和有 积分的求积公式及其误差的渐进展开式，此类求积 限元 Galerkin方法 ，这是由Sloan等(1988)提出的， 公式只需赋值，不需计算二重积分。因此，计算量 配置法是对非奇异部分用一组正交基的线性组合 小，同时收敛阶大为提高。利用这类积分公式进行 表示，然后用配点求积来处理，每个离散元素需要 计算可以得到十分精确的结果。这为后面讨论多 计算一重奇异积分；有限元 Galerkin方