离散元(二forcopy)PPT课件.ppt

二 基本原理-球形颗粒元离散元法 岩土工程研究所 刘 军 根据离散化模型中所采用的单元种类分别介绍离散元法的基本原理： 颗粒元 二维圆盘单元 三维圆球单元 块体元 多边形单元 多面体单元 二 基本原理 基本假设 假定速度和加速度在每个时间步长内为常量 ； 选取的时间步长应该足够小以至于在单个时间步长内扰动的传播不会超过当前与之相邻的粒子 。 二 基本原理-球形颗粒元离散元法 运动描述 处于一个理想散体中的任意一个颗粒,具有6个自由度,3个平动自由度与三个转动自由度,可通过Newton第二定律分别描述。 二 基本原理-球形颗粒元离散元法 运动描述 平动方程： 基本原理-球形颗粒元离散元法 式中， 与 分别为颗粒 的质量和速度。 为时间， 为颗粒的重力， 与 分别为颗粒 与 的接触力与粘性接触阻尼力， 为所有与颗粒接触的颗粒总数。 运动描述 接触力的分解： 颗粒 与 间的接触力可分解为法向与切向接触力，即 基本原理-球形颗粒元离散元法 同理，粘性接触阻尼力也可分解为法向与切向分量形式，即 运动描述 接触力产生的力矩： 基本原理-球形颗粒元离散元法 颗粒间的接触力作用在两个颗粒的接触点上，而不是作用在颗粒的中心，所以这些接触力（除法向接触力 外）将会对颗粒产生力矩 ， 式中， 为从颗粒 的质心指向接触点的矢量，其幅值为 （颗粒的半径）。 运动描述 转动方程： 基本原理-球形颗粒元离散元法 转动方程可以表示为 式中， 与 分别为颗粒 的转动惯量与角速度，对于球形颗粒 为 接触模型 综述： 二 基本原理-球形颗粒元离散元法 关于接触力的计算模型已有大量的研究成果，目前仍旧是一个活跃的研究领域，特别是对于切向力的计算方法。 对于理想散体颗粒（无粘连）：采用Hertz理论描述法向作用，而采用Mindlin与Deresiewicz理论描述切向作用； 对于存在粘连的散体颗粒：法向接触力根据在Hertz理论基础上考虑粘连力的JKR（Johnson-Kendall-Roberts）理论确定，切向接触力增量则根据把Savkoor和Briggs理论与Mindlin和Deresiewicz理论相结合形成的理论确定。 接触模型 两个处于接触颗粒单位法向和切向向量： 基本原理-球形颗粒元离散元法 单位法向向量 单位切向向量 单位切向量之所以通过两个颗粒的相对速度来计算，是因为接触力与粘性阻尼力的方向与相对速度的方向相同。 接触模型 两个处于接触颗粒接触点的相对速度： 基本原理-球形颗粒元离散元法 法向相对速度为 切向相对速度为 或者写为 接触模型 法向接触力计算模型 —Hertz模型： 基本原理-球形颗粒元离散元法 为颗粒i与j接触时的侵入深度 式中 接触模型 法向接触力计算模型— Cundall模型： 基本原理-球形颗粒元离散元法 式中， 为法向弹簧刚度。 接触模型 法向接触力计算模型 —法向粘性接触阻尼力： 二 基本原理-球形颗粒元离散元法 式中， 为法向粘性接触阻尼系数。 接触模型 切向接触力计算模型 —综述： 二 基本原理-球形颗粒元离散元法 处于接触中的两个颗粒的切向作用，从本质上讲，是一种摩擦行为，按照摩擦机理，摩擦力包括：滑动摩擦、滚动摩擦与静摩擦，其中滑动摩擦与静摩擦属于切向摩擦力；滚动摩擦是由于法向接触应力的不均匀分布产生的。介绍两个切向接触力模型： Coulomb准则 Mindlin与Deresiewicz切向接触力模型 接触模型 切向接触力计算模型 —Coulomb准则： 二 基本原理-球形颗粒元离散元法 在离散元模拟中，一般用Coulomb准则这种简单的形式描述，静摩擦的详细刻画需要涉及切向位移甚至可能要考虑时间依赖效应。 式中， 为静摩擦系数，切向摩擦力的方向为与相对滑动的趋势相反。 接触模型 切向接触力模型 —Mindlin与Deresiewicz模型 ： 二 基本原理-球形颗粒元离散元法 式中 为颗粒 与 间的累积切向位移矢量 接触模型 切向接触力模型 —阻尼力 ： 二 基本原理-球形颗粒元离散元法 式中， 为切向粘性接触阻尼系数。 二 基本原理-球形颗粒元离散元法 计算模型总结 运动方程 接触力的计算 法向接触力 切向接触力 二 基本原理-球形颗粒元离散元法 运动方程的求解 一般采用两种方法求解运动方程： 中心差分法