无奇异积分的径向积分边界单元法.doc

《无奇异积分的径向积分边界单元法.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、无奇异积分的径向积分边界单元法高效伟东南大学工程力学系，江苏南京210096Email:xwgao@seu.edu.cn摘要本文提出了一种无奇异积分的边界单元法。该方法用基于紧支径向基函数的形函数为权函数，用加权余量法建立伽辽金弱形式的边界－区域积分方程，并由径向积分法将区域积分转换成等效的边界积分，形成只需要边界离散的纯边界元算法。本方法兼有边界单元法和无网格法的主要特点，导出的系统方程的系数矩阵为带状稀疏矩阵，对线性问题、非线性问题以及变系数问题具有统一的表达形式，是一种很有发展潜力的方法。关键

2、词:边界单元法，无网格法，径向积分法，奇异积分1引言数值计算方法已经变成工程技术人员解决大型工程问题的重要工具，现最流行的数值方法是有限单元法[1],以其强大的解题能力和通用灵活性而普遍被采用。然而随着其应用范围的不断扩大，固有的一些缺陷也日益突出，如在金属成形、优化设计、隐藏区域的识别、渗流问题自由面的确定等移动边界问题中，基于拉格朗日法的有限元网格可能产生畸变和重叠，严重地影响解的精度。边界单元法是继有限单元法后又一种被广泛应用的数值方法[2]，其主要特点是只需要将问题的边界离散成单元，因而准备

3、数据简单、便于复杂边界问题的建模[3,4]。此外，边界单元法在解决上述移动边界问题中具有突出的优势[5]，计算中移动边界的位移与原边界节点的坐标相加就自然形成了新的边界单元信息[6]，不需要专门重构单元，也不会有网格重叠的问题。然而，传统的边界单元法存在着固有的弱点：(1)所用基本解导致有奇异积分的出现，数值计算时需要消除积分奇异性；(2)所形成的系数矩阵是非对称满阵，限制了解题规模和速度；(3)在解决非均质和非线性问题时由于很难求得问题的基本解，因而不得不用对应线性问题的基本解来建立非均质[7]和

4、非线性[8]问题的积分方程，结果有区域积分出现在积分方程中。为了计算区域积分，传统的方法是将问题的区域离散成内部网格,用网格积分法来计算区域积分[9,10]，这种做法消除了边界单元法只需将边界离散成单元的特点。正是这些弱点严重地影响了边界单元法的发展。鉴于有限元和传统非线性边界元需要内部网格的缺陷，近年来许多学者致力于无网格算法的研究[11-13]。纯无网格法不需要借助于任何单元，只需要在计算域中布置一系列的离散点即可，因而具有很强的解题灵活性。此外，无网格法通常采用基于紧支基函数的形函数[12]，

5、形成的系统方程的系数矩阵为带状稀疏矩阵，占用较小的内存空间。然而无网格法还发展的很不成熟，现提出的十几种无网格法中存在的缺点和不足表现在:1)缺少坚实的理论基础和严格的数学证明,因此计算精度、守恒性等一直没有明确的答案;2)需要布置较多的内部点,布点方案不仅限制计算规模而且严重影响计算精度;3)由于没有边界单元（或内部网格）可以利用，对于具有复杂几何边界条件的移动边界问题，在迭代计算中重新产生以及自动判断内部、外部还是边界点时会有困难。为了利用边界元和无网格法所具有的优点,国际上许多学者致力于无内部

6、网格边界单元法的研究，以便只通过对问题的边界离散来求解非线性和非均质问题。其中最常用的方法是用某种特殊技术将出现在传统积分方程中的区域积分转换成等效的边界积分，从而避免对内部区域的离散。其中应用最广泛的方法是Brebbia等人提出的双重互易法（DRM）[14]和多重互易法（MRM）[15]。这类方法的基本思想是利用微分算子的特解，将区域积分转换成边界积分。虽然DRM和MRM的应用范围很广泛，但对于一些复杂的基函数和微分算子，求特解是很困难的。此外，如果在同一积分方程中存在有不同核函数的区域积分，则很

7、难求得统一的特解，以致该方法无法实施。基于这种缺陷，高效伟在2002年提出了径向积分法（RIM）[16,17]。该方法通过首先计算一个径向积分的途径，将一个区域积分转换成边界积分。众所周知，著名的高斯散度公式和格林公式只可将含微分算子的区域积分转换成边界积分，但RIM是一种建立在纯数学处理上的方法，不借助任何特解和微分算子，也不需要使用基本解则可将任何一个区域积分转换成边界积分。不仅能对二维和三维区域积分给出统一的转换公式，还可以显式地消除积分中存在的奇异性[16]。由于基于RIM的径向积分边界元算

8、法能够不需要内部网格而有效地求解非线性和非均质问题，因此该方法提出后已得到了广泛的应用[18,19]和好评[20,21]。其中，Albuquerque等人[20]用算例对RIM与DRM进行了数值比较，证明了RIM比DRM更精确、更稳定；而Hematiyan[21]评论到，双重互易法(DRM）是用得最多的方法，径向积分法(RIM)则是最有力的方法。这在一定程度上说明了径向积分法的发展潜力。虽然径向积分边界单元法求解非线性和非均质问题时不需要内部网格，但仍然包含有奇异积分