第5单元： 边界元法课件.pptVIP

边界元法Boundary Element Method 第一章 绪 论 工程中边界元法 边界积分方程＋离散法 当代三大数值解法之一 边界元法特征 奇异降维数值解法 2阶微分方程－转换－边界积分方程 降维 降阶 奇异 矩阵非对称 满阵 边界元法发展史 始于1978年 Brebbia [1] C. A.等（英国南安普敦大学） 《边界元法的基础和应用》 跻身于 当代计算力学数值解析理论与工具 （有限差分法、有限元法、边界元法）之林 应用领域 1.简史 边界元法有了迅速的发展 （高性能、大容量微机出现） 大型和非线性科学工程 结构静力学 裂纹的生成及扩展 流体运动 骨骼生长 接触问题等研究领域 国内简史 在国内，1978年起步 杜庆华院士 率先推动工程中边界元法 冯康、胡海昌、何广乾院士等 加入到边界元法的研究者行列 我国边界元法研究得到了迅速的发展 研究起点和热点 我国大部分工程中边界元法 固体力学方面开始 后迅速转入非线性问题领域 出版自然边界元法1993 弹塑性BEM 弹塑性问题的边界元法 Swedlow Cruse 和Mendelson上世纪70年代初分别给出二维弹塑性问题的边界积分方程 Mukherje和Kumar提出平面问题公式中的错误并作了修正 Tells和Brebbia给出了边界应力的计算方法和奇异积分的数值处理方案。 岑章志博士采用初应力形式给出三维弹塑性问题边界元法，取得奇异积分的计算法及加速收敛方法 弹塑性BEM 陈政清博士给出弹塑性大变形边界元法 完成了拉伸试件颈缩定量数值模拟 肖宏博士建立三维弹塑性有限形变边界元法和轧制过程模拟边界元法源程序 给出板带轧制过程变形－面力－应力场 很好的处理奇异问题 规模局限性 典型边界元法计算结构（边界积分方程－影响系数数值积分－矩阵方程及消去法求解）局限性 系数积分计算和方程组的求解时间长，占用大量的计算机内存和主机CPU的时间 抵消边界元法降维特性 对单一三维弹性体静力学问题，当离散网格节点数超过2000时，因计算时间过长而失去可行性 工程接触问题的矛盾更加突出，如果再叠加弹塑性和摩擦非线性，在微机上将无法得到收敛解 研究热点与方向 边界元法发展至关重要问题 扩大运算规模提高计算效率 普遍注意到边界元法求解过程 运算量主要集中在两个方面 影响系数计算和方程组求解法 1）分域法 2）网络并行法 3）多极展开法 2. 多极边界元法 多极展开法 （Fast Multipole Methods，简称FMM方法） 求和项的近似计算方法 计算的近似程度和展开的项数密切相关 直接计算的计算量和n2成正比 上世纪80年代后期 Greengard等 提出快速多极展开法，基于球谐函数在空间中的多级展开，采用递归算法结构使计算量级降为一阶 1999年底 清华大学的姚振汉教授在国内介绍了多极展开法在边界元法的应用 多极展开法 FMM法来源于计算电荷集中各个电荷间相互作用力的方法，如二维静电模等离子体粒子模拟 Rokhlin将多极展开法用于Laplace’s边界积分方程的迭代求解 在二维问题中运算等价点是 三维椭圆边值问题中等价点是 最近，Greengard and Rokhlin提出了三维多极展开法，并且将其运算等价点降为 多极边界元法简史 G. etal给出二维双调和方程的多极边界元法 Peirce and Napier给出二维线弹性Navier方程多极边界元法，计算效率较差 Sangani and Mo对三维低雷诺数方程（Stokes’s方程）的积分方程适应多极展开法，或者使多极展开法适应积分方程，上述适应变换是繁琐而往往使问题特殊化，因为刚性质点，不能扩展到三维线弹性问题 Greenbaum etal的多极边界元法方法不能扩展到三维问题 国内多极边界元法 近几年在国内 银燕等将多极展开法用于无网格等离子体静电模粒子的模拟 王浩刚对多极展开法中不变项的计算进行了优化 清华大学姚振汉教授课题组率先将泰勒级数展开式适用于2维弹性问题边界元法中，模拟了复合材料不同粒子变形－应力场 燕山大学申光宪教授课题组将域多极展开法扩展至边界展开法，给出3维弹塑性摩擦接触多极边界元法，完成四辊轧机板带轧制模拟计算 3. 轧制工程中边界元法 轧制理论 板带和轧辊耦合轧制理论 重要的研究方向－轧制过程数值模拟 有限差分法 主要的方法有：有限元法 边界元法 板带和轧辊轧制理论回顾 起始于20世纪40年代 Hitchcock给出轧辊弹性压扁公式 20世纪50年代 St