离散元课件(一).pptVIP

离散单元法的基本理论及应用 The Theory of DEM (Discrete Element Method) and Its Applications 离散单元法的基本理论及应用 历史由来及研究现状 基本原理 程序设计及商业软件介绍 应用 离散单元法的基本理论及应用 历史由来及研究现状 产生背景 发展现状 存在的问题 离散单元法的基本理论及应用 基本原理 Cundall二维圆盘单元离散元法 三维球体单元(颗粒元)离散元法 多边形单元离散元法 多面体单元离散元法 接触发现算法 Cundall公共面法 其他方法 离散单元法的基本理论及应用 程序设计及商业软件介绍 程序设计方法 UDEC 3DEC PFC2D与PFC3D 离散单元法的基本理论及应用 应用 理想散体岩土介质中的应用 堆积问题 磨矿 铁路道砟 离散单元法的基本理论及应用 应用 在连续岩土介质中的应用 岩石力学性质分析 颗粒破碎模拟 边坡稳定性分析 洞室围岩稳定性分析 桩与土介质相互作用模拟 分层岩土介质中侵彻与爆炸模拟 一 历史由来及研究现状 产生背景 散粒岩土材料在自然界中普遍存在 从本质上讲，岩土材料都是由离散的、尺寸不一、形状各异的颗粒或块体组成的，例如，土就是松散颗粒的堆积物，同样，天然岩体也是由被结构面切割而成的大小不一、形态各异的岩石块体所组成。散粒岩土材料的力学特性有着重要的工程应用，如泥砂的沉淀，土堤、土（岩）坡、铁路道渣等的稳定性研究，散粒岩土材料的力学特性研究是岩土力学中最基本的、也是最重要的问题之一。 一 历史由来及研究现状 产生背景 用连续介质力学研究散粒岩土材料力学特性的不足 连续介质力学把散粒体作为一个整体来考虑，研究的重点放在建立粒子集合的本构关系，从粒子集合整体的角度研究散粒体介质的力学行为。 不足：不能体现颗粒间的复杂相互作用及高度非线性行为；不能真实刻画散体材料的流动变形特征。 一 历史由来及研究现状 用多体动力学描述散粒体的力学行为的困难 ? （1）对于被研究的多粒子系统而言，已经存在的接触不断地分开，而新的接触频繁的形成，在多体动力学中，接触的分开与形成都需要改变控制方程； ? （2）即使接触网络保持相同，在每一个接触中，也可能发生在依附与滑动间的过渡，而这种过渡也会导致系统运动方程的改变。 因而：多体动力学方法只能描述少数散体体系的力学行为，对于大量散体组成的岩土材料则相当困难。 一 历史由来及研究现状 分子动力学方法的引入 ? 分子动力学模拟是一种用来计算一个经典多体体系的平衡和传递性质的方法。所谓的经典意味着颗粒体系的运动遵守经典力学定律。该方法最初是用来描述分子运动的（当处理一些较轻的原子或分子时，才需要考虑量子效应）。? 分子动力学方法模拟分子的运动时，邻近分子间存在吸引或排斥力。该方法可以模拟大量分子的运动。 去除分子间作用力，把分子动力学中的小尺度粒子作为散体岩土材料中的颗粒，并入颗粒间及颗粒与边界间的相互作用描述，即是Cundall离散元法的最初思路。 一 历史由来及研究现状 早期的离散单元法 ? 离散元法的思想源于较早的分子动力学（Molecular Dynamics）。 1971年Cundall提出适于岩石力学的离散元法； 1979年Cundall和Strack又提出适于土力学的离散元法，并推出二维圆盘（Disc）程序BALL和三维圆球程序TRUBAL（后来发展为商业软件PFC-2D/3D），形成较系统的模型与方法，被称为软颗粒模型。 1988年Cundall所在的ITASCA咨询公司推出针对三维块体元的3DEC程序。至此，离散元的理论体系基本形成。 一 历史由来及研究现状 早期的离散单元法 ? Cundall称之为“Distinct Element Method”，随着该方法的推广，有的学者称其为“Discrete Element Method”，缩写形式均为DEM。 最初，离散元的研究对象主要是岩石等非连续介质的力学行为，它的基本思想是把不连续体分离为刚性元素的集合，使各个刚性元素满足运动方程，用时步迭代的方法求解各刚性元素的运动方程，继而求得不连续体的整体运动形态。离散元法允许单元间的相对运动，不一定要满足位移连续和变形谐调条件，计算速度快，所需存储空间小，尤其适合求解大位移和非线性问题。 一 历史由来及研究现状 早期的离散单元法 ? 主要思路：用Newton定律描述颗粒运动，通过颗粒间及颗粒与边界间的相互作用传递载荷，求解方法是解藕的。 理论难点：接触力模型（Contact Force Model)与接触发现算法 (Contact Detection Algorithm)。