离散元经典课件.ppt

基本思路在解决连续介质力学问题时，除了边界条件外，还有3个方程必须满足：平衡方程、变形协调方程（弹性力学中的几何方程）、本构方程（弹性力学中的物理方程）；对于离散单元法而言，由于介质被假定为离散体（颗粒或块体）的集合，所以离散体单元间不存在变形协调的约束，但是其运动必须遵循经典力学定律；三基本原理-多边形单元离散元法运动描述处于一个理想散体中的任意一个多边形单元,具有3个自由度,2个平动自由度与一个转动自由度,可通过Newton第二定律分别描述。三基本原理-多边形单元离散元法运动描述运动方程（平动和转动）：三基本原理-

多边形单元离散元法

式中，为质量矩阵；为第j个单元给第i个单元的总作用力。接触模型Cundall模型（弹簧阻尼器模型）：三基本原理-

多边形单元离散元法

与颗粒元相同，不再重复介绍。接触模型简化的解析模型：基本原理-

多边形单元离散元法

式中，与分别为法向与切向阻尼系数，为静摩擦系数，为材料的弹性模量，为重叠部分的面积，为特征长度。应用举例基本原理-多边形单元离散元法

应用举例基本原理-多边形单元离散元法

应用举例—UDEC基本原理-多边形单元离散元法

综述问题的复杂性：计算量大；接触的力学模型块体元的计算机生成；模拟结果的可视化；但更为复杂的是块体间的接触发现算法。基本原理-多面体单元离散元法

块体接触的直接判别接触方式：判别块体间是否存在接触的最简单方法是检查接触发生的所有可能性。对于三维块体的接触有很多方式：点-点，点-边，点-面，边-边，边-面，面-面。基本原理-多面体单元离散元法

块体接触的直接判别直接判别的计算效率举例：基本原理-多面体单元离散元法

如果块体有个顶点、条边、个面，第二个块体有个顶点、条边、个面，则用直接法判别接触存在的计算次数为用直接法，两个立方体间存在676种接触的能性。块体接触的直接判别直接判别计算效率的改善：基本原理-多面体单元离散元法

事实上，并不需要如此多次判别，因为有些接触类型可以并入到其他类型中，这样只需要检测点-面、边-边接触两种类型即可，而其他类型可用下面的方法通过点-面、边-边接触两种类型来体现，即（1）当在同样位置存在三个或更多个点-面接触时，说明在该位置存在点-点接触；（2）当两个点-面接触在同一位置同时存在时，表明在该位置存在点-边接触；（3）当两个块体间存在两个边-边接触时，表明存在边-面接触；（4）当三个或更多个边-边接触存在或三个或更多个点-面接触存在时，表明两个块体间存在面-面接触。接触模型两个处于接触颗粒单位法向和切向向量：基本原理-

球形颗粒元离散元法

单位法向向量单位切向向量单位切向量之所以通过两个颗粒的相对速度来计算，是因为接触力与粘性阻尼力的方向与相对速度的方向相同。接触模型两个处于接触颗粒接触点的相对速度：基本原理-

球形颗粒元离散元法

法向相对速度为切向相对速度为或者写为接触模型法向接触力计算模型—Hertz模型：基本原理-球形颗粒元离散元法为颗粒i与j接触时的侵入深度式中接触模型法向接触力计算模型—Cundall模型：基本原理-

球形颗粒元离散元法

式中，为法向弹簧刚度。接触模型法向接触力计算模型—法向粘性接触阻尼力：二基本原理-球形颗粒元离散元法

式中，为法向粘性接触阻尼系数。接触模型切向接触力计算模型—综述：二基本原理-球形颗粒元离散元法

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

在离散元模拟中，一般用Coulomb准则这种简单的形式描述，静摩擦的详细刻画需要涉及切向位移甚至可能要考虑时间依赖效应。式中，为静摩擦系数，切向摩擦力的方向为与相对滑动的趋势相反。接触模型切向接触力模型—Mindlin与Deresiewic