离散元、有限元和数值流形法的对比.docxVIP

离散元方法也被称为散体单元法最早是 1971 年由 Cundall 提出的一种不连续数值方法 模型离散元理论是由分析离散单元的块间接触入手找出其接触的本构关系建立接触的物理 力学模型并根据牛顿第二定律建立力、加速度、速度 及其位移之间的关系对非连续、离散 的单元进行模拟仿真。 离散元法是专门用来解决不连续介质问题的数值模拟方法。 该方法把节理岩体视为由离 散的岩块和岩块间的节理面所组成 , 允许岩块平移、 转动和变形 , 而节理面可被压缩、 分离或 滑动。因此 , 岩体被看作一种不连续的离散介质。其内部可存在大位移、旋转和滑动乃至块 体的分离 , 从而可以较真实地模拟节理岩体中的非线性大变形特征。离散元法的一般求解过 程为 : 将求解空间离散为离散元单元阵 , 并根据实际问题用合理的连接元件将相邻两单元连 接起来 ; 单元间相对位移是基本变量 , 由力与相对位移的关系可得到两单元间法向和切向的 作用力 ; 对单元在各个方向上与其它单元间的作用力以及其它物理场对单元作用所引起的外 力求合力和合力矩 , 根据牛顿运动第二定律可以求得单元的加速度 ; 对其进行时间积分 , 进而 得到单元的速度和位移。 从而得到所有单元在任意时刻的速度、加速度、角速度、 线位移和 转角等物理量。 离散单元法的特点 岩体或颗粒组合体被模拟成通过角或边的相互接触而产生相互作用。 块体之间边界的相互作用可以体现其不连续性和节理的特性。 使用显式积分迭代算法允许有大的位移、转动和使用。 在岩体计算力学方面， 由于离散单元能更真实地表达节理岩体的几何特点， 便于处理所有非 线性变形和破坏都集中在节理面上的岩体破坏问题， 被广泛应用于模拟边坡、 滑坡和节理岩体地下水渗流等力学过程 离散单元法的求解过程 离散元法具体的求解过程分为显式解法和隐式解法下面分别介绍 其适用范围。 显式解法 显式解法用于动力问题的求解或动态松弛法的静力求解显式算法无 须建立像有限元法那样的大型刚度矩阵只需将单元的运动分别求出计算比较简单数据量较 少并且允许单元发生很大的平移和转动可以用来求解一些含有复杂物理力学模型的非线性 问题时间积分采用中心差分法由于条件收敛的限制使得计算步长不能太大因而增加了计算 时间。 隐式解法 隐式解法用于求解静力问题的静态松弛法隐式解法的动态松弛法式直接找 导块体失去平衡后达到再平衡的力位移关系建立隐式方法解联立方程组并通过迭代求解以 完全消除块体的残余力和力矩。 有限元方法的基本概念 将介质复杂几何区域离散为具有简单几何形状的单元而单元内的材料性质和控制方程通过 单元节点的未知量来进行表达再通 过单元集成、外载和约束条件的处理得到方程组求解该 方程组就可以得到该介质行为的近似表达。 有限元法的优点 可以用有限的、 相互关联的单元模拟无限的复杂体， 无论多么复杂的几何体都能用相应的单 元简化，从而建模分析计算出结果。 使复杂的、感觉无处下手的工程问题简单化， 这是最大 的优点。 有限元法采用矩阵形式表达，编程性高 对于线弹性问题当实际结构位移场函数连续光滑时能够得到收敛解。 对于任意复杂结构理论上总是可以通过细分单元的方法获得足够近似的模拟。 长期大量工程应用积累了丰富的经验。 有限元的插值是基于网格的， 所以需要人为做好单元， 这很耗时间， 但是单元就好像人们修了路一样， 计算的时候可以节省很多时间， 效率比较高。 同时，这也是有限元法的一个缺点，大变形问题中的网格畸变问题，本质在于单元插值造成的 缺点 ：精确度浮动性比较大。基于建模的水平和边界条件、载荷工况的模拟是否真实等等。 数值流形方法简介 由石根华博士提出 , 在最近几年发展起来的数值流形方法是一种一般的数值方法。这种方法 采用有限覆盖体系 ( 一套数学覆盖 , 一套物理覆盖两者相互独立地定义 , 但又有一定的依赖关 系), 特别适合模拟断续介质材料的变形和物体的大位移运动。 而以往的有限元方法及非连续 变形分析 (DDA)可以作为数值流形方法的特殊情况。 数值流形方法具有两套网格 , 即物理网格和数学网络 . 物理网格是由分析域的边界、 节理、块 体及不同材料区域的界面所组成 , 它是不能人为选择的材料条件 . 而数学网格可以任意选择 , 近期研究基本上均采用了有限元网格作为数学网格 . 数学网格可直接转换为有限数学覆盖 ; 由数学覆盖与物理网络形成物理覆盖系统 , 而物理覆盖的交集 ( 公共区域 ) 称为流形意义下的 单元 . 数值流形方法是在每个物理覆盖上建立覆盖函数 ( 覆盖位移函数 ), 在几个覆盖的公共 区域内 ( 单元 ), 将其上所有覆盖位移函数加权求和即可形成适应于该域的总体位移函数 , 以 此根据总势能变分原理建立求解岩体力学问题 ( 包括一般结构体 ) 的数值流