6.0数值方法.ppt

6.1 概述 6.2 有限元法 6.3 有限差分法 6.4 离散元法 有限单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。它是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。 由于其灵活处理边界、模拟非线性，广泛应用于岩石工程领域 有限元法分析计算的基本思想： 物体离散化 单元特性分析 选择位移模式 分析单元的力学性质 计算等效节点力 单元组集 求解未知节点位移 物体离散化 将某个工程结构离散为由各种单元组成的计算模型，这一步称作单元剖分。 离散后单元于单元之间利用单元的节点相互连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算进度而定。 用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获 得的结果就与实际情况相符合。 分析单元的力学性质 根据 单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，这是单元分析中的关键一步。此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式，从而导出单元刚 度矩阵，这是有限元法的基本步骤之一。 选择位移模式 位移法：选择节点位移作为基本未 知量称为位移法； 力 法：选择节点力作为基本未 知量时称为力法； 混合法：取一部分节点力和一部分节点位移作为基本未知量时称为混合法。 位移法易于实现计算自动化，所以，在有限单元法中位移法应用范围最广。 计算等效节点力 物体离散化后，假定力是通过节点从一个单元 传递到另一个单元。但是，对于实际的连续体，力是从单元的公共边传递到另一个单元中去的。因而，这种作用在单元边界上的表面力、体积力和集中力都需要等效的移到节点上去，也就是用等效的节点力来代 替所有作用在单元上得力。 单元组集 利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来，形成整体的有限元方程 求解未知节点位移 可以根据方程组的具体特点来选择合适的 计算方法。 有限差分法主要思想是将待解决问题的基本方程组和边界条件（一般为微分方程）近似的改用差分方程（代数方程）来表示，由有一定规则的空间离散点的处场变量（应力、位移）的代数表达式代替。 有限差分法和有限元法都产生一组待解决的方程组，但两者产生方程组的方式不一样 有限元： 插值函数 有限差分：差分代替微分 6.3 有限差分法 设f(x,y)为弹性体内某一个连续函数，在邻近点0处，函数f可以展开为泰勒级数： 在结点3和1处，x分别等于x0-h及x0+h，即x-x0分别等于-h和h。将其代入上式得 6.3 有限差分法 假定h充分小，可以略去三次幂或更次幂的各项，上式可简化为： 联立求解得差分公式： 同理 6.3 有限差分法 利用上式可以推导出混合二阶导数的差分公式： 平衡方程 平衡方程可转化为由各结点的待求量表示的差分方程，结合边界条件即可求解。 平衡方程 （运动方程） 应力-应变关系 （本构方程） 新的速度和 位移 新的应力 或力 商业软件： FLAC2D、FLAC3D 6.3 有限差分法 FLAC—土坡滑动问题 FLAC—土坡滑动问题 FLAC—土坡滑动问题 FLAC—土坡滑动问题 FLAC—土坡滑动问题 FLAC—土坡滑动问题 FLAC—土坡滑动问题 FLAC—土坡滑动问题 * 第六章 岩石力学数值分析方法简介 数值计算是人类认识世界的新手段。自理论分析与科学实验之后，数值模拟已成为人类认识世界最重要的手段。 它主要用来解决以下两类问题：不可能进行实验的问题，进行实验代价太大的问题。同时它又融和了理论分析和科学试验的特点，数值模拟 / 仿真已经不再局限于科学计算，正广泛被用在科学研究，工程与生产领域。  6.1 概述 6.1 概述 解析方法的障碍 岩体是地质结构体，非均质、非连续、非线性、以及复杂的加卸载条件和边界条件， 岩石力学问题通常无法用解析方法简单求解 数值方法的优势 数值方法适用范围广， 能模拟岩体的复杂结构和力学特性， 分析各种边值问题和施工过程 对工程进行预测和预报 6.1 概述 主要用于研究岩土工程活动和自然环境变化过程中岩体及其加固结构的力学行为和工程活动对周围环境的影响。 常用方法:有限元、边界元法、有限差分法、加权余量法、离散元法、刚体元法、不连续变形分析方法、流形方法。 前四种基于连续介质力学方法，后三种基于非连续介质力学方法 6.1 概述 有限元基于最小总势能变分原理，能方便处理各种非线性问题，灵活模拟岩土工程中复杂的施工过程，在岩石力学领域应用最为广泛 边界元法以表述拜特(Betti)互等定理的积分方程为基础，建立了直接法的基本