离散元课件二.pptVIP

第1页，共35页，星期日，2025年，2月5日根据离散化模型中所采用的单元种类分别介绍离散元法的基本原理：颗粒元二维圆盘单元三维圆球单元块体元多边形单元多面体单元二基本原理第2页，共35页，星期日，2025年，2月5日基本假设假定速度和加速度在每个时间步长内为常量；选取的时间步长应该足够小以至于在单个时间步长内扰动的传播不会超过当前与之相邻的粒子。二基本原理-球形颗粒元离散元法

第3页，共35页，星期日，2025年，2月5日运动描述处于一个理想散体中的任意一个颗粒,具有6个自由度,3个平动自由度与三个转动自由度,可通过Newton第二定律分别描述。二基本原理-球形颗粒元离散元法

第4页，共35页，星期日，2025年，2月5日运动描述平动方程：基本原理-

球形颗粒元离散元法

式中，与分别为颗粒的质量和速度。为时间，为颗粒的重力，与分别为颗粒与的接触力与粘性接触阻尼力，为所有与颗粒接触的颗粒总数。第5页，共35页，星期日，2025年，2月5日运动描述接触力的分解：颗粒与间的接触力可分解为法向与切向接触力，即基本原理-

球形颗粒元离散元法

同理，粘性接触阻尼力也可分解为法向与切向分量形式，即第6页，共35页，星期日，2025年，2月5日运动描述接触力产生的力矩：基本原理-

球形颗粒元离散元法

颗粒间的接触力作用在两个颗粒的接触点上，而不是作用在颗粒的中心，所以这些接触力（除法向接触力外）将会对颗粒产生力矩，式中，为从颗粒的质心指向接触点的矢量，其幅值为（颗粒的半径）。第7页，共35页，星期日，2025年，2月5日运动描述转动方程：基本原理-

球形颗粒元离散元法

转动方程可以表示为式中，与分别为颗粒的转动惯量与角速度，对于球形颗粒为第8页，共35页，星期日，2025年，2月5日接触模型综述：二基本原理-球形颗粒元离散元法

关于接触力的计算模型已有大量的研究成果，目前仍旧是一个活跃的研究领域，特别是对于切向力的计算方法。对于理想散体颗粒（无粘连）：采用Hertz理论描述法向作用，而采用Mindlin与Deresiewicz理论描述切向作用；对于存在粘连的散体颗粒：法向接触力根据在Hertz理论基础上考虑粘连力的JKR（Johnson-Kendall-Roberts）理论确定，切向接触力增量则根据把Savkoor和Briggs理论与Mindlin和Deresiewicz理论相结合形成的理论确定。第9页，共35页，星期日，2025年，2月5日接触模型两个处于接触颗粒单位法向和切向向量：基本原理-

球形颗粒元离散元法

单位法向向量单位切向向量单位切向量之所以通过两个颗粒的相对速度来计算，是因为接触力与粘性阻尼力的方向与相对速度的方向相同。第10页，共35页，星期日，2025年，2月5日接触模型两个处于接触颗粒接触点的相对速度：基本原理-

球形颗粒元离散元法

法向相对速度为切向相对速度为或者写为第11页，共35页，星期日，2025年，2月5日接触模型法向接触力计算模型—Hertz模型：基本原理-球形颗粒元离散元法为颗粒i与j接触时的侵入深度式中第12页，共35页，星期日，2025年，2月5日接触模型法向接触力计算模型—Cundall模型：基本原理-

球形颗粒元离散元法

式中，为法向弹簧刚度。第13页，共35页，星期日，2025年，2月5日接触模型法向接触力计算模型—法向粘性接触阻尼力：二基本原理-球形颗粒元离散元法

式中，为法向粘性接触阻尼系数。第14页，共35页，星期日，2025年，2月5日接触模型切向接触力计算模型—综述：二基本原理-球形颗粒元离散元法

处于接触中的两个颗粒的切向作用，从本质上讲，是一种摩擦行为，按照摩擦机理，摩擦力包括：滑动摩擦、滚动摩擦与静摩擦，其中滑动摩擦与静摩擦属于切向摩擦力；滚动摩擦是由于法向接触应力的不均匀分布产生的。介绍两个切向接触力模型：Coulomb准则Mindlin与Deresiewicz切向接触力模型第15页，共35页，星期日，2025年，2月5日接触模型切向接触力计算模型—Coulomb准则：二基本原理-球形颗粒元离散元法

在离散元模拟中，一般用Coulomb准则这种简单