离散元方式与有限元方式的比较.docx

离散元方式与有限元方式的比较 摘要 离散元方式是由分析离散单元的块间接触入手，找出其接触的 本构关系，成立接触的物理力学模型，并依照牛顿第二定律，对非 持续、离散的单元进行模拟仿真。而有限元方式是将介质复杂几何 区域离散为具有简单几何形状的单元，通过单元集成、外载和约束 条件的处置，取得方程组，再求解该方程组就能够够取得该介质行 为的近似表达。 本文中并介绍刚体一弹簧元法及极限平稳法，还有离散元法有 限元法结合之应用，和工程中的离散元方式的应用实例。本文中介 绍的实例有：丽江地震区应力场研究及离散变量结构拓扑优化设计 研究及基于混合离散复合形法的工程优化设计及离散元与壳体有 限元结合的多尺度方式及其应用和昌马水库枢纽工程右岸岩石边 坡稳固性的离散元法分析。 关键词：离散元方式、有限元方式、刚体一弹簧元法、极限平稳法 离散元方式 1.1离散元方式的大体概念⑴ 离散元方式也被称为散体单元法，最先是1971年由Cundall 提出的一种不持续数值方式模型，离散元理论是由分析离散 单元的块间接触入手，找出其接触的本构关系，成立接触的 物理力学模型，并依照牛顿第二定律成立力、加速度、速度 及其位移之间的关系，对非持续、离散的单元进行模拟仿真。 1.2离散元方式的历史背景⑵ 离散元法又称DEM(Discrete Element Method)法，它的思想 源于较早的分子动力学(Molecular Dynamics)o 1971年由Cundall 最先提出，其研究对象是岩石等非持续介质的力学行为。1979 年，Cundall和Strack又提出适于土力学的离散元法。国内显现 了用于土木工程设计的块体离散元分析系统2D-Block和三维离 散单元法软件TRUDEC；在冲击波研究方面，唐志平等成立了 二维和三维细观离散元理论和DM2程序。 1.3离散单元法的特点⑶ ?岩体或颗粒组合体被模拟成通过角或边的彼此接触而产生彼此 作用。 ?块体之间边界的彼此作用能够表现其不持续性和节理的特性。 利用显式积分迭代算法，许诺有大的位移、转动和利用。 1.4离散单元法的求解进程 离散元法具体的求解进程分为显式解法和隐式解法，下面别离介 绍其适用范围。 显式解法【4】： 显式解法用于动力问题的求解或动态松弛法的静力求解，显式算 法不必成立像有限元法那样的大型刚度矩阵，只需将单元的运动 别离求出，计算比较简单，数据量较少，而且许诺单元发生专门 大的平移和转动，能够用来求解一些含有复杂物理力学模型的非 线性问题，时刻积分采纳中心差分法，由于条件收敛的限制，使 得计算步长不能太大，因此增加了计算时刻。 隐式解法【4】： 而隐式解法用于求解静力问题的静态松弛法，隐式解法的动态松 弛法式直接找导块体失去平稳后达到再平稳的力一位移关系，成 立隐式方式解联立方程组，并通过迭代求解以完全排除块体的残 余力和力矩。 有限元方式 2.1有限元方式的大体概念【5】 将介质复杂几何区域离散为具有简单几何形状的单元，而单元 内的材料性质和操纵方程通过单元节点的未知量来进行表达，再通 过单元集成、外载和约束条件的处置，取得方程组，求解该方程组 就能够够取得该介质行为的近似表达。 2.2有限元方式的历史背景⑸ Hrenikoff于1941年采纳框架形变功法计算了弹性问题， Courant于1943年发表了采纳三角形区域内的分片多项式来处置 扭转问题的论文；Turner等人于1956年推导了杆、梁等单元的 刚度矩阵；而“有限单元”这一名称是Clough于1960年提出。 第一本关于有限元方式的书是Zienkiewicz和Cheung于1967年 完成的；1972年Oden完成了有关非线性介质方面的专着，现在， 随着运算机的进展和普及，使得学生和工程师能够充分的利用有 限元方式这一有力的工具。 2.3有限元法的优势【3】 ?关于线弹性问题，当实际结构位移场函数持续滑腻时，能够取得 收敛解。 ?关于任意复杂结构，理论上老是能够通过细分单元的方式取得足 够近似的模拟。 ?刚度矩阵系数带状，在结构不显现软化的时候仍是对称正定的， 求解方便。 ?长期大量工程应用，积存了丰硕的体会。 其它数值方式 3.1刚体一弹簧元法 刚体弹簧元法【3】（Rigid Body Spring Method ） 又称RBSM, 最先由Kawai于1976年提出，当初提出的用意是以较少的自由 度来求解结构问题。它把体系分解为一些由均布在接触面上的弹 簧系统联系起来的刚性元，刚性元本身不发生弹性变形，因此结 构的变形能仅能贮存在接触面的弹簧系统中，由于刚体弹簧元单 元间的作使劲通过单元界面上弹簧传递，故可直接取得界面的作 使劲，因此在极限分析等领域也有着较好的应用。 3.2极限平稳法 极限平稳法【1】是一种近似的计算方式，采纳