2022年离散元精心整理.pptxVIP

- 1 二 基本原理-球形颗粒元离散元法 岩土工程争论所 刘 军 CHYT — HZOM — Z3738 这是第1页 - 2 依据离散化模型中所接受的单元种类分别介绍离散元法的基本原理： 颗粒元 二维圆盘单元 三维圆球单元 块体元 多边形单元 多面体单元 二 基本原理 CHYT — HZOM — Z3738 这是第2页 - 3 基本假设 假定速度和加速度在每个时间步长内为常量 ； 选取的时间步长应当足够小以至于在单个时间步长内扰动的传播不会超过当前与之相邻的粒子 ； 二 基本原理-球形颗粒元离散元法 CHYT — HZOM — Z3738 这是第3页 - 4 运动描述 处于一个抱负散体中的任意一个颗粒,具有6个自由度,3个平动自由度与三个转动自由度,可通过Newton其次定律分别描述； 二 基本原理-球形颗粒元离散元法 CHYT — HZOM — Z3738 这是第4页 - 5 运动描述 平动方程： 基本原理-球形颗粒元离散元法 式中， 与 分别为颗粒 的质量和速度； 为时间， 为颗粒的重力， 与 分别为颗粒 与 的接触力与粘性接触阻尼力， 为全部与颗粒接触的颗粒总数； CHYT — HZOM — Z3738 这是第5页 - 6 运动描述 接触力的分解： 颗粒 与 间的接触力可分解为法向与切向接触力，即 基本原理-球形颗粒元离散元法 同理，粘性接触阻尼力也可分解为法向与切向重量形式，即 CHYT — HZOM — Z3738 这是第6页 - 7 运动描述 接触力产生的力矩： 基本原理-球形颗粒元离散元法 颗粒间的接触力作用在两个颗粒的接触点上，而不是作用在颗粒的中心，所以这些接触力（除法向接触力 外）将会对颗粒产生力矩 ， 式中， 为从颗粒 的质心指向接触点的矢量，其幅值为 （颗粒的半径）； CHYT — HZOM — Z3738 这是第7页 - 8 运动描述 转动方程： 基本原理-球形颗粒元离散元法 转动方程可以表示为 式中， 与 分别为颗粒 的转动惯量与角速度，对于球形颗粒 为 CHYT — HZOM — Z3738 这是第8页 - 9 接触模型 综述： 二 基本原理-球形颗粒元离散元法 关于接触力的运算模型已有大量的争论成果，目前仍然是一个活跃的争论领域，特殊是对于切向力的运算方法； 对于抱负散体颗粒（无粘连）：接受Hertz理论描述法向作用，而接受Mindlin与Deresiewicz理论描述切向作用； 对于存在粘连的散体颗粒：法向接触力依据在Hertz理论基础上考虑粘连力的JKR（Johnson-Kendall-Roberts）理论确定，切向接触力增量就依据把Savkoor和Briggs理论与Mindlin和Deresiewicz理论相结合形成的理论确定； CHYT — HZOM — Z3738 这是第9页 - 10 接触模型 两个处于接触颗粒单位法向和切向向量： 基本原理-球形颗粒元离散元法 单位法向向量 单位切向向量 单位切向量之所以通过两个颗粒的相对速度来运算，是由于接触力与粘性阻尼力的方向与相对速度的方向相同； CHYT — HZOM — Z3738 这是第10页 - 11 接触模型 两个处于接触颗粒接触点的相对速度： 基本原理-球形颗粒元离散元法 法向相对速度为 切向相对速度为 或者写为 CHYT — HZOM — Z3738 这是第11页 - 12 接触模型 法向接触力运算模型 —Hertz模型： 基本原理-球形颗粒元离散元法 为颗粒i与j接触时的侵入深度 式中 CHYT — HZOM — Z3738 这是第12页 - 13 接触模型 法向接触力运算模型— Cundall模型： 基本原理-球形颗粒元离散元法 式中， 为法向弹簧刚度； CHYT — HZOM — Z3738 这是第13页 - 14 接触模型 法向接触力运算模型 —法向粘性接触阻尼力： 二 基本原理-球形颗粒元离散元法 式中， 为法向粘性接触阻尼系数； CHYT — HZOM — Z3738 这是第14页 - 15 接触模型 切向接触力运算模型 —综述： 二 基本原理-球形颗粒元离散元法 处于接触中的两个颗粒的切向作用，从本质上讲，是一种摩擦行为，依据摩擦机理，摩擦力包括：滑动摩擦，滚动摩擦与静摩擦，其中滑动摩擦与静摩擦属于切向摩擦力；滚动摩擦是由于法向接触应力的不均匀分布产生的；介