有限差分法、边界元法和离散元法(共3页).doc

《有限差分法、边界元法和离散元法(共3页).doc》由会员分享，可在线阅读，更多相关《有限差分法、边界元法和离散元法(共3页).doc（3页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上有限差分法已经发展的一些近似数值分析方法中，最初常用的是有限差分法，它可以处理一些相当困难的问题。但对于几何形状复杂的边界条件，其解的精度受到限制，甚至发生困难。作为60年代最重要的科技成就之一的有单元法。在理论和工程应用上都_得到迅速发展，几乎所有用经典力学解析方法难以解决的工程力学问题郁可以用有限元方法求解。它将连续的求解域离散为一组有限个单元的组合体，解析地模拟或逼近求解区域。由于单元能按各种不同的联结方式组合在一起，且单元本身又可有不同的几何形状，因此可以适应几何形状复杂的求解域。相限元的另一特点是利用每一单元内假设的近似函数来表示全求解区域上待求的未知场函数

2、。单元内的近似函数由未知场函数在各个单元结点上数值以及插值函数表达，这就使未知场函数的结点值成为新的未知量，把一个连续的无限自由度问题变成离散的有限自由度问题，只要结点来知量解出，便可以确定单元组合体上的场函数。随着单元数目的增加，近似解收敛于精确解。但是有限元方法常常需要很大的存贮容量，甚至大得无法计算；由于相邻界面上只能位移协调，对于奇异性问题(应力出现间断)的处理比较麻烦。这是有限单元法的不足之处。边界元法边界元法是在有限元法之后发展起来的一种较精确有效的工程数值分析方法。与有限元法在连续体域内划分单元的基本思想不同，边界元法是在定义域的边界上划分单元，用满足控制方程的函数去逼近边界条件

3、，通过对边界分元插值离散，化为代数方程组求解。降低了问题的维数，可用较简单的单元准确地模拟边界形状，利用微分算子的解析的基本解作为边界积分方程的核函数，而具有解析与数值相结合的特点，通常具有较高的精度。边界元法的主要缺点是它的应用范围以存在相应微分算子的基本解为前提，对于非均匀介质等问题难以应用，故其适用范围远不如有限元法广泛，而且通常由它建立的求解代数方程组的系数阵是非对称满阵，对解题规模产生较大限制。上述两种数值方法的主要区别在于，边界元法是“边界”方法，而有限元法是“区域”方法，但都是针对连续介质而言，只能获得某一荷载或边界条件下的稳定解。对于节理裂隙发育的岩体或颗粒散体的处理则要麻烦得

4、多，更无法进行大变形、分离、回转及塌落过程的模拟。这就使得人们去探索和寻求适合模拟节理岩体和颗粒散体运动变形特性的有效数值方法。离散元法离散元法是由Cundall P A (1971) 首先提出并应用于岩土体稳定性分析的一种数值分析方法。它是一种动态的数值分析方法,可以用来模拟边坡岩体的非均质、不连续和大变形等特点,因而,也就成为目前较为流行的一种岩土体稳定性分析数值方法。该方法在进行计算时,首先将边坡岩体划分为若干刚性块体(目前已可以考虑块体的弹性变形) ,以牛顿第二运动定律为基础,结合不同本构关系,考虑块体受力后的运动及由此导致的受力状态和块体运动随时间的变化。它允许块体间发生平动、转动,

5、甚至脱离母体下落,结合CAD 技术可以在计算机上形象地反应出边坡岩体中的应力场、位移及速度等力学参量的全程变化。该方法对块状结构、层状破裂或一般碎裂结构岩体比较适合。caoyb19822009-05-24 09:09的数值方法有很多，可以分为两大类：一类是连续介质力学数值方法，另一类是非连续介质力学数值方法。其中连续介质力学数值方法将岩体简化成数学意义上的连续体来进行分析，主要有：有限差分法、有限元法、边界元法、无单元法等等。 !Hq$7j_ 有限差分法（Finite Difference Method，FDM）是求解偏微分方程的最重要的数值方法之一，其主要思想是将微分方程近似地用相应的差分方

6、程来替代，从而将求解偏微分方程的问题转化为求解代数方程组的问题。 ,Ta k, 有限元分析（finite element analysis ，fea）是使用有限元方法来分析静态或动态的物理物体或物理系统，是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的)近似解，然后推导求解这个域总的满足条件(如结构的平衡条件)，从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。然而，当岩体中

7、存在着大量结构面时，将为上述接触单元和接触模型的使用带来困难，处理连续介质力学的有限单元法对于需要考虑大量结构面的非连续介质力学问题如岩石边坡稳定等，就显现出其局限性。 3,INo 边界积分方程 边界元法（ Boundary Integral Equation-BoundaryElement Method）简称边界元法（BEM）是继有限元之后发展起来的一种有效的数值分析方法。其基本思想是以边界积分方程为数学基础，同时采用与有限元法类似的离散技术，通过将边界离散为边界元，将边界 +xY 积分方程离散为线性代数方程组，再由数值方法求解线性代数方程组，从而得到原问题的边界积分方程解。 U9g;t+/

8、 无单元法（Element-free/Mesh-less Method）是Nayroles80等于1992年针对有限元法的一些缺点，如网格依赖性、奇异边界等提出的一种新的数值方法。其基本方法是采用节点信息及其局部支撑域上的权函数实行局部精确逼近，通过配点法或伽辽金法得到微分方程弱形式，再选用合适的积分方案聚合整体平衡方程，从而实现对问题的求解。 j.=&qYc0 扩展有限元法87（eXtended Finite Element Method，X-FEM）是近年来新发展起来的数值方法之一，该方法基于单位分解（PU）理论对传统有限元方法进行了扩展，引入完全独立于网格划分的非连续位移模式来表征裂纹尖

9、端不连续界面的演化，因此计算过程中不需要预设开裂路径和调整计算网格。 c8z6-6i0 连续介质力学数值方法均采用连续体假定，必须满足应力平衡和位移协调条件，因此，在模拟岩石工程计算中，如滑坡体的变形破坏、坝基、隧道工程的失稳等方面均有一定的局限。随着计算技术的发展，出现了很多新的处理非连续介质力学的方法 bu0i # 极限平衡法是岩土工程稳定分析中最为广泛使用的一种方法。该法以摩尔库仑的抗剪强度理论为基础，对滑动岩体进行力平衡分析，结合结构面的强度参数得到抗滑稳定安全系数刚体弹簧元法（Rigid Body-Spring Model, RBSM）首先将结构体离散化为一系列块体，每个块体包含六个自由度，块体与块体之间用弹簧连接，每个块体本身均是一个刚性体。该方法本质上是由有限元法演绎而来的， #aI(fQZe 与有限元不同的是，刚体弹簧元在块体形心处插值，用块体形心的位移作为基本未知量，用分片的刚体位移去逼近实际整体位移场 F*uLTH 离散单元法（Discrete/Distinct Element Method，DEM）是1971 年由Cundall115提出的一种分析节理岩体的数值计算方法，最初是为了模拟岩质边坡的破坏过程。 专心-专注-专业