边界元法chapter01—概述

1、边界元法边界元法武汉大学水利水电学院武汉大学水利水电学院王桥王桥教学内容教学内容参考教材参考教材1、姚振汉等姚振汉等. 边界元法边界元法.高等教育出版社，高等教育出版社，20102、高效伟等高效伟等. 高等边界元法高等边界元法-理论与程序理论与程序. 科学出版社科学出版社 ，2015边界元法基本理论边界元法基本理论第一章第一章 概述概述1.1 边界元法的发展边界元法的发展1.2 边界元法的特点边界元法的特点本章主要内容：本章主要内容：1.3 边界元法的应用边界元法的应用1.1 边界元方法的发展边界元方法的发展工程的数值方法工程的数值方法边值问题的数值解边值问题的数值解有限差分法有限差分法FDM

2、有限元法有限元法FEM边界元法边界元法BEM有限体积法有限体积法finite volume method多重网格法多重网格法multigrid method无网格法无网格法meshless method1.1 边界元方法的发展边界元方法的发展有限差分法：FLAC3D FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。通过调整三维网格中的多面体单元来拟合实际的结构。单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发变形和移动（大变形模式）。 FLAC3D采用了显式拉格朗日算法和混合-离散分区技

3、术，能够非常准确地模拟材料的塑性破坏和流动。由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。FLAC3D是采用ANSI C+语言编写的。1.1 边界元方法的发展边界元方法的发展有限元方法： 有限元法（FEA，Finite Element Analysis）的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元

4、不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。 在岩土工程中，主要特点FEM应反映岩土材料的特点，如含结构面、流变、开挖 等。软件有：ANSYS,ABACUS1.1 边界元方法的发展边界元方法的发展可分为两个阶段：第一阶段第一阶段: 19世纪各种积分方程的出现世纪各种积分方程的出现 Helmholtz (1859), Rayleigh(1889), Kirchhoff(1882)提出各种边界积分方程. 第二阶段第二阶段: 20世纪边界元方法的出现世纪边界元方法的出现 对于复杂的问题, 很难得到边界积分方程的解析解. 随着计算机技术的发展, 并借鉴有限元方法, 边界元方

5、法发展成熟. 边界元方法就是边界积分方程的离散形式. 边界元方法的发展边界元方法的发展1.1 边界元方法的发展边界元方法的发展边界元方法的发展边界元方法的发展1976年, Lachat JC 和 Waston JO 提出了有效的3D弹性力学的边界元数值方法(Int. J. Numer. Meth. Engng 1976; 10:991-1005). 弹性力学边界元方法具有里程碑 意义的两篇文章: 1969年, Rizzo FJ.发表了第一篇边界元文章(Quarterly Applied Mathematics 25(1):83-95), 形成有效的分元数值计算方法.研究对象为2D弹性力学问题.

6、1.1 边界元方法的发展边界元方法的发展边界元方法的发展边界元方法的发展一些专集也相继出现, 促进了边界元方法的发展. 比较有名的有: Cruse A. 和 Rizzo FJ. (1975) 出版了边界元专集. Banerjee PK (1979-1986) 出版了边界元专集I-IV. 目前, 有许多国际期刊刊登边界元方面的文章, 如: Engineering Analysis with Boundary Element International Journal of Solids and Structures Journal of Applied Mechanics, ASME. Mech

7、anics of Materials International Journal of Fracture 国内刊物有: 力学学报, 固体力学学报, 计算力学学报1.1 边界元方法的发展边界元方法的发展边界元方法的发展边界元方法的发展边界元发展史1.1 边界元方法的发展边界元方法的发展我国边界元方法的发展我国边界元方法的发展 1978年开始,我们开始学习国外的文献。 清华大学工程力学系的杜庆华、姚振汉教授组织研讨班，推动边界元方法在我国的发展。 冯康、胡海昌、何广乾院士等加入到边界元法的研究者行列1.1 边界元方法的发展边界元方法的发展我国边界元方法的发展我国边界元方法的发展不同学科的发展情况：

8、水利工程： 张楚汉院士（清华大学）发展了动力边界元方法，研究大坝与地基的相互作用。 周维垣教授（清华大学）发展了裂隙岩体渗流的边界元方法. 张有天教授（水科院）发展了岩体地下水的渗流岩土工程： 王泳嘉教授（东北大学）发展了粘弹性边界元方法 朱合华教授（同济大学）地下工程的开挖 谭云亮的FEM和BEM耦合;刘承伦教授间接边界元方法.1.2 边界元方法的特点边界元方法的特点1. 边界元方法的原理：偏微分方程边值问题边界积分方程(采用加权余量或其它方法)-边界元方法(离散方法).2. 特点: 降低求解问题的维数, 3D问题变为2D问题, 2D变为1D问题. 具有较高的精度, 原因: 仅仅对边界进行离

9、散, 域内点的值采用边界上的 已知量计算得到. 采用辅助解, 而该解为解析解. 适合于无限域问题1.2 边界元方法的特点边界元方法的特点3. 局限性: 适应性差, 没有通用程序。 处理弹塑性、动力学、含体积力的问题， 需要在域内划分单元，效率降低。 奇异性。 矩阵非对称，满阵。1.3 边界元方法的应用边界元方法的应用土木建筑工程机械工程海洋工程航天工程环境工程强、弱电工程生物工程推动各工程理论解析和创新技术的发展发挥了巨大作用1.3 边界元方法的应用边界元方法的应用边界元法有了迅速的发展（高性能、大容量微机出现）大型和非线性科学工程结构静力学裂纹的生成及扩展流体运动骨骼生长接触问题等研究领域1.3 边界元方法的应用边界元方法的应用1. 一些岩土工程应用 Map3D程序，见网站：，澳大利亚Mine Modelling Pty Ltd 单一断层分析单一断层分析：断层采用DD单元， 开挖岩体采用块体单元1. 一些岩土工程应用单一断层分析单一断层分析：断层采用DD单元， 开采岩体采用块体单元1.3 边界元方法的应用边界元方法的应用/research/geomech/Software/Software.htmPoly3D程序：美国Stanford大学编写模拟 1999 Turkish earthquakes (after