边界元群面多极展开法的探索.pdfVIP

维普资讯 第28卷 第6期 太 原 科 技 大 学 学 报 vo1．28 No．6 2007年 12月 JOURNALOFTAIYUANUNIVERSITYOFSCIENCEANDTECHNOLOGY Dec．2007 文章编号：1673—2057(2007)06—0466—03 边界元群面多极展开法的探索 张志平，冀俊杰，黄庆学，王建梅 (太原科技大学，太原030024) 摘 要：在现有边界元快速多极展开法(FMM—BEM)的基础上，将群面多极展开法和 广义极小残值法应用于三维弹性问题的边界元法中，变革计算结构，以适应大规模数值 计算，提高运算精度。 关键词：边界元法；多极展开法；广义极小残值法 中图分类号：O241．6 文献标识码 ：A 在微机上使用传统边界元法 ¨’2求解三维弹性 1．1 群面多极边界元法的数学理论 问题要求运算规模较小，其制约因素主要是必须建 假设有 Ⅳ个强度为g，q，…，g 的电荷分别位 立方程组的系数矩阵且是满秩的，没有稀疏性和对 于点 ， ，… ，如图1所示。其球坐标分别为(p， 称性。因此，需要 占用大容量的内存 ，而方程组的 求解一般采用高斯消去法或由它演变得到的其他 X(r，8 ) 方法。快速多极展开法 (FMM)与传统边界元法结 合，始于静电场的大规模快速近似计算，将三维弹 性基本解核函数分解成多极展开法格式，引入广义 极小残值法 (TheGeneralizedMinimalResidualAI— gorithm 简称GMRES)，建立快速多极边界元法可 以更新传统边界元法的计算结构，适应大规模运算 需要。而传统的边界元法计算结构，要求解N个未 知数问题所需要的内存量为 D(』、72)和 O ( )次 图 1 多极展开 Fig．1M ultipoleexpansion 运算，当N2000时难 以在微机上完成计算。多 1，卢1)，(p2，2，卢2)，…，(pⅣ，Ⅳ，Ⅳ)，并且点 l， 极展开法是求解静电场问题的有效方法，它体现了 ， … 位于圆心在原点，半径为 a的球 内。由此， 力学圣维南原理将电荷 (或力源)的远近不同影响 就有在任意点X ：(r，，)∈R。，在 ra时，由电 区分开来加以分别处理的科学思想，应用于传统的 边界元法中，更新其计算结构，给出求解三维弹性 荷g，g，…，g产生的位势 ()为： 问题所需要的内存量为D(Ⅳ)方法 。 ()：壹塞M7．U(o，) (1) 1 群面多极边界元法的数学理论与算法 其中：M7=∑g·p · (，) (2) 收稿 日期：2007-01．10 基金项 目：国家 自然科学基金项 目(5057155) 作者简介：张志平 (1969一)，男，在读硕士，研究方向为多极边界元法。 维普资讯 第28卷第6期 张志平，等：边界元群面多极展开法的探索