边界元素法之介绍.pptVIP

内 容 提 要 引言 由偏微分方程建立边界积分方程 求解边界积分方程的边界元法 位势问题边界元法的应用例题 边界元法的各种应用 快速多极边界元法的大规模计算 边界型无网格法简介 进一步研究方向 快速多极算法发展概况 传统的求解 N particle 问题算法的计算量： O(N 2) 静电势 静电力 Barnes, Hut (1986) ： 树结构算法 (tree code) 计算量：O (NlogN) J. Barnes, P. Hut. A hierarchical O (NlogN) force calculation algorithm. Nature, 324, (1986), 446-449. 快速多极算法发展概况 Rokhlin, Greengard (1987)： 快速多极算法 (Fast multipole method, FMM)： O(N) 采用了： 多极展开 (multipole expansion) 局部展开 (local expansion) Carrier等 (1988) 自适应快速多极算法 (Adaptive FMM)： O(N) 快速多极算法发展概况 Greengard, Rokhlin (1997): 3D新型快速多极算法 (New Version FMM) O(N) 引入 指数展开 (Exponential expansion) Hrycak, Rokhlin (1998)： 2D新型快速多极算法 O(N) 二维含大量随机分布圆形夹杂的复合材料大规模计算 三维含大量随机分布球形颗粒的复合材料大规模计算 沿X轴单向拉伸时颗粒/基体界面法向应力分布图 三维含大量随机分布和随机方向短纤维的复合材料大规模计算 并行计算的计算规模 界面法向面力分布（隐去基体显示） 二维无限弹性体中3000裂纹 (900,000 DOF, p = 20, 6.15 h) 域内变量用基本解来插值，从而可以将弱形式中的域内积分化为边界积分。 边界变量采用移动最小二乘插值。 象边界元法那样降低维数， 又是真正的无网格法， 只需要建立在边界上分布的点的信息。 边界型无网格法简介 用FM-HBNM求解大规模三维位势问题 620,708 节点，53,840 秒 含复杂形状空腔 边界型无网格法简介 碳纳米管复合材料的一种代表性胞元 边界型无网格法简介 得到的等效导热系数和碳纳米管(CNT)长度的关系 边界型无网格法简介 碳纳米管复合材料的另一种代表性胞元 边界型无网格法简介 得到的等效导热系数和碳纳米管波纹度的关系 边界型无网格法简介 边界元法和各种无网格法等在许多工程领域都属于有限元法的一个补充。 它的发展既受到有限元法的启发，又受到有限元法的制约。只有在其优势领域，作为补充才有实际意义。如果只是有限元能计算的它也能算，而且结果吻合，那就证明这种补充没有必要。 进一步研究方向 快速多极边界元法已经证明可以有效地进行有限元法都难以完成的大规模计算。当务之急是为它寻找工程上有重要意义的实际应用。 由于线性问题的有限元方法已经非常成熟，而边界型数值方法在非线性分析方面的研究还很少，因此快速多极边界元法在非线性分析中的应用值得进行研究。我们在弹塑性分析快速多极边界元法方面已经初步取得成功。 进一步研究方向 谢谢诸位！ 边界元法的各种应用 二维弹性力学问题模拟光弹性实验 (本组, 1998) 边界元法的各种应用 三维弹性应力分析 (本组, 2000) 滚筒轴边界元模型：5个子域，8184结点，1721单元 边界元法的各种应用 三维弹性应力分析 (本组, 2000) 滚筒轴有限元模型：13000结点，64000单元 边界元法的各种应用 三维弹性应力分析 (本组, 2000) 边界元解：左端局部Mises应力 (单位：MPa) 有限元解：左端局部Mises应力 (MPa) 便宜了2000倍 关键是改进了算法 规模增大400倍 提高速度9倍 PC机 1个Intel Pentium III 快速多极边界元法 20,000,000个常值面元 每小时1,200,000常值元 大规模并行机 1900个Intel i860 传统边界元法 50,000个常值面元 每小时140,000常值元 计算机 CPU 算法 计算能力 计算速度 评 论 Coyote AutoMEMS 软 件 Sandia BEM软件(1994)