三维边界元法中拟奇异积分的正则化及其在薄体与涂层结构中的应用-数学专业论文.docxVIP

山东理工大学硕士学位论文摘要 山东理工大学硕士学位论文 摘要 摘 要 边界元法（BEM）作为一种重要的数值方法，因其具有计算精度高、降维等优点， 近年得到了很大发展，已被广泛应用于科学及工程问题的数值分析中。然而，几乎奇 异积分的存在大大地制约了它的应用领域，如三维薄体与涂层结构问题等。虽然，已 有许多计算几乎奇异积分的方法，但是大部分都局限于平面几何单元的采用。近年来， 尽管在曲面单元上的几乎奇异积分的计算方面取得很大的进展，但是这些方法仅在单 个曲面单元上的简单几乎奇异积分计算中得到验证，从未在薄体及涂层结构等实际边 值问题中得到检验。薄体与涂层结构一般厚度很薄，有的在微米级甚至纳米级，它的 数值分析一直是科学计算工作者面临的挑战。近年来，研究者发现，边界元法求解此 类问题的瓶颈是高阶曲面单元上的几乎奇异积分的计算。 本文提出了一种计算三维边界元法中高阶几何单元上的 2D 拟奇异积分的有效方 法。该方法的新颖之处在于：（1）构造了‘精确’的距离公式。此距离公式不仅能够 精确地描述场点（源点）与积分单元上的一般点之间的距离 r ，同时它便于与各种非 线性变换结合；（2）拓展 Sinh 变换到 3D 边界元法中，并与前述的距离公式有机地结 合，有效地降低了距离函数在积分区间上的剧烈震荡，从而充分地将几乎奇异的被积 函数规则化。 所提出的方法已在三维薄体与涂层结构中的边界元分析中得到了验证，取得了很 理的效果。它表明，本文解决了最一般形式的曲面单元，即 8 节点二次曲面单元上的 几乎奇异积分的计算问题，使得准确、高效地计算三维薄体及涂层结构问题及 3D 边 界元法中的边界层效应问题成为可能。因此，本文不仅成功地解决了薄体及涂层结构 计算的困难问题，同时拓展了边界元法的应用领域。 关键词：3D BEM；拟奇异积分；薄体结构；涂层结构 I 山东理工大学硕士学位论文Abs 山东理工大学硕士学位论文 Abstract II II Abstract Due to the merits of high accuracy and reducing the number of unknowns by one order, the boundary element method (BEM) has become an appealing area of research and is widely used in engineering problems. However, the existence of nearly singular integral limits its application in a lot of aspects, such as the 3D thin-body and coating structure problems. Although many methods of evaluating two-dimensional (2D) nearly singular integrals have been developed , most of the techniques are used to the lower-order geometry elements, e.g. usual linear or planar elements. In recent years, the computation of the nearly singular integrals over surface elements has achieved great progress, but only the single element integrals are tested. Generally speaking, most of the engineering processes occur in complex geometrical domains, such as thin body and coating structures whose thickness are in the nano-or micro-meters. However, the numerical analysis of these problems by BEM is always a great challenge. In recent years, scholars found that the bottleneck of the problems is the evaluation of nearly singular integral on higher-order elements. In this paper,