1、二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法 岩土工程研究所岩土工程研究所刘刘 军军 根据离散化模型中所采用的单元种类分别介绍离根据离散化模型中所采用的单元种类分别介绍离散元法的基本原理:散元法的基本原理:颗粒元二维圆盘单元三维圆球单元块体元多边形单元多面体单元二二 基本原理基本原理基本假设假定速度和加速度在每个时间步长内为常量;选取的时间步长应该足够小以至于在单个时间步长内扰动的传播不会超过当前与之相邻的粒子。二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法运动描述运动描述 处于一个理想散体中的任意一个颗粒,具有6个自由度,3个平动自由度与三个转动自由度,可通过N
2、ewtonNewton第第二定律二定律分别描述。二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法运动描述运动描述平动方程平动方程:二二基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法式中,与 分别为颗粒 的质量和速度。为时间,为颗粒的重力,与 分别为颗粒 与 的接触力与粘性接触阻尼力,为所有与颗粒接触的颗粒总数。运动描述运动描述接触力的分解接触力的分解:颗粒 与 间的接触力可分解为法向与切向接触力,即二二基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法同理,粘性接触阻尼力也可分解为法向与切向分量形式,即 运动描述运动描述接触力产生的力矩接触力产生的力矩:二二基本原理基本
3、原理-球形颗粒元离散元法球形颗粒元离散元法颗粒间的接触力作用在两个颗粒的接触点上,而不是作用在颗粒的中心,所以这些接触力(除法向接触力 外)将会对颗粒产生力矩 ,式中,为从颗粒 的质心指向接触点的矢量,其幅值为 (颗粒的半径)。运动描述运动描述转动方程转动方程:二二基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法转动方程可以表示为式中,与 分别为颗粒 的转动惯量与角速度,对于球形颗粒 为 接触模型接触模型综述综述:二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法关于接触力的计算模型已有大量的研究成果,目前仍旧是一个活跃的研究领域,特别是对于切向力的计算方法。对对于于理理想想散
4、散体体颗颗粒粒(无无粘粘连连):采用Hertz理理论论描述法向作用,而采用Mindlin与Deresiewicz理论描述切向作用;对对于于存存在在粘粘连连的的散散体体颗颗粒粒:法法向向接接触触力力根据在Hertz理理论论基础上考虑粘连力的JKR(Johnson-Kendall-Roberts)理论确定,切切向向接接触触力力增量则根据把Savkoor和Briggs理论与Mindlin和Deresiewicz理论相结合形成的理论确定。接触模型接触模型两个处于接触颗粒单位法向和切向向量两个处于接触颗粒单位法向和切向向量:二二基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法单位法向向量 单位切向
5、向量 单位切向量之所以通过两个颗粒的相对速度来计算,是因为接触力与粘性阻尼力的方向与相对速度的方向相同。接触模型接触模型两个处于接触颗粒接触点的相对速度两个处于接触颗粒接触点的相对速度:二二基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法法向相对速度为 切向相对速度为 或者写为 接触模型接触模型法向接触力计算模型法向接触力计算模型 HertzHertz模型模型:二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法 为颗粒i与j接触时的侵入深度式中接触模型接触模型法向接触力计算模型法向接触力计算模型 CundallCundall模型模型:二二基本原理基本原理-球形颗粒元离散元法球形
6、颗粒元离散元法式中,为法向弹簧刚度。接触模型接触模型法向接触力计算模型法向接触力计算模型 法向粘性接触阻尼力法向粘性接触阻尼力:二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法式中,为法向粘性接触阻尼系数。接触模型接触模型切向接触力计算模型切向接触力计算模型 综述综述:二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法处于接触中的两个颗粒的切向作用,从本质上讲,是一种摩擦行为,按照摩擦机理,摩擦力包括:滑动摩擦、滚动摩擦与静摩擦,其中滑动摩擦与静摩擦属于切向摩擦力;滚动摩擦是由于法向接触应力的不均匀分布产生的。介绍两个切向接触力模型:Coulomb准则准则Mindlin与
7、与Deresiewicz切向接触力模型切向接触力模型 接触模型接触模型切向接触力计算模型切向接触力计算模型 Coulomb准则准则:二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法 在离散元模拟中,一般用Coulomb准则这种简单的形式描述,静摩擦的详细刻画需要涉及切向位移甚至可能要考虑时间依赖效应。式中,为静摩擦系数,切向摩擦力的方向为与相对滑动的趋势相反。接触模型接触模型切向接触力模型切向接触力模型 Mindlin与与Deresiewicz模型模型:二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法式中 为颗粒 与 间的累积切向位移矢量接触模型接触模型切向接触力模型切
8、向接触力模型 阻尼力阻尼力:二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法式中,为切向粘性接触阻尼系数。二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法计算模型总结计算模型总结运动方程运动方程接触力的计算接触力的计算法向接触力法向接触力 切向接触力切向接触力 二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法运动方程的求解运动方程的求解一般采用两种方法求解运动方程:一般采用两种方法求解运动方程:中心差分法中心差分法Verlet积分法积分法二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法运动方程的求解运动方程的求解中心差分法中心差分法 运动方程可由
9、Verlet显式积分求解。通过积分可获得粒子的新位置,积分时需要粒子的当前及上一步长的位置数据,而不需要粒子的速度数据。二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法运动方程的求解运动方程的求解Verlet积分法积分法二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法运动方程的求解运动方程的求解Verlet积分法积分法定义定义:从而导出Verlet方程为接触发现算法接触发现算法 在一个由众多颗粒组成的体系中,直接判别颗粒是否接触需要耗费大量的计算时间,因而,为了节约计算时间,提高计算效率,一般不直接判别任意两个颗粒间是否存在接触,而是分两个步骤判别颗粒间的接触是否存在:
10、首先,对一个颗粒,判别其潜在的邻居个数,然后,准确确定该颗粒与每个邻居是否接触。虽然在确定邻居数目时也要耗费一定的计算时间,但是仍旧比逐个准确判别颗粒间接触是否存在要节约时间。因而,接触发现算法的效率在多颗粒体系力学行为模拟中至关重要。二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法接触发现算法接触发现算法介绍三种针对球形颗粒的接触发现算法:Verlet邻居目录邻居目录法法连接单元法连接单元法 边界盒法边界盒法 二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法接触发现算法接触发现算法Verlet邻居目录邻居目录法法二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散
11、元法 当需要判别体系中某个颗粒的邻居数量时,在该粒子周围构建一个球(称之为参考球,称该颗粒为核心颗粒),参考球半径为体系中最大粒子半径的若干倍,那么参考球所包围的所有粒子为该球中心粒子的邻居。参考球半径的选取取决于粒子的运动速度及体系中粒子的密度。对于每个粒子,都可生成一个邻居粒子的目录。为了得到邻居目录,对每一个粒子而言,所有标号大于该粒子的粒子都必须被检验,判断是否位于该粒子的参考球中,而对于标号小于该粒子标号的粒子则没有必要被检验,因为邻居是互相的,没有必要对一个邻居对检验两次。Verlet邻居目录邻居目录法法及粒子存储目录及粒子存储目录二二 基本原理基本原理-球形颗粒元离散元法球形颗粒
12、元离散元法接触发现算法接触发现算法Verlet邻居目录邻居目录法法 二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法接触发现算法接触发现算法Verlet邻居目录邻居目录法法对于n个颗粒组成的体系,用Verlet邻居目录法需要 次计算,也就是说计算次数仍旧为 量级。然而,并不需要在每个时间步长上都对邻居目录进行更新。更新的频率取决于体系中粒子的密度、粒子的运动速度以及参考球的尺寸。参考球的半径也可以根据颗粒体系的稠密程度及运动速度进行动态调整,并且参考球半径与邻居目录的更新频率呈反比关系,参考球半径越小,邻居目录的更新频率越高;但是参考球半径越大,则有更多的粒子位于球体内,所以判别是
13、否为邻居就需要较长的时间。接触发现算法接触发现算法连接单元法连接单元法二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法将颗粒体系所占据的空间划分成规则的网格,对于三维问题,可以划分为 个立方体单元,对于二维体系,则划分为 个正方形单元,对于颗粒体系所占据空间形状不规则时,也可采用其他形状单元划分。但是所有单元的尺寸必须大于粒子的尺寸。与Verlet邻居目录法的主要区别在于:相邻单元法中的单元不依附于粒子,单元不随粒子的运动而运动。如果粒子就当前的位置被分配到某个单元,显然,只有在同一个单元或直接相邻单元内的粒子间才可能发生相互作用,也就是说,只有相邻单元内的粒子才能成为邻居。接触发
14、现算法接触发现算法连接单元法连接单元法二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法例如,对一个二维体系而言,只有在9个不同的单元内可能包含邻居粒子,对于三维体系,则只有在27个不同单元内可能包含邻居粒子。与前面介绍的Verlet邻居目录法相同,每个粒子对只需要检验一次,这样,没必要对9个单元都进行检验,只需检验中心单元及邻居单元的一半即可,即在二维情况下,只需检验5个单元,在三维情况下,只需检验14个单元。接触发现算法接触发现算法连接单元法连接单元法二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法2D体系中的连接单元法体系中的连接单元法接触发现算法接触发现算法边界盒
15、法边界盒法二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法这个方法与前述两种方法不同。首先,在每一个粒子周围构建一个边界盒,边界盒的尺寸按这样的方式选取:使每个粒子刚好放进它使每个粒子刚好放进它的边界盒内的边界盒内。边界盒的边为直线,并且与体系的坐标轴平行。在判别颗粒的邻居时,把边界盒投影到体系的坐标轴上。通过边界盒在坐标轴上投影的起点和终点来判别是否为邻居。接触发现算法接触发现算法边界盒法边界盒法二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法 a a 围绕围绕在每在每个粒个粒子周子周围的围的边界边界盒盒 b b两个不两个不同时刻同时刻粒子的粒子的边界盒边界盒在轴上在
16、轴上的投影的投影 在判别颗粒的邻居时,把边界盒投影到体系的坐标轴上。例如,图a表明了粒子及边界盒的位置,图b为图a中的边界盒在体系x轴的投影。通过边界盒在坐标轴上投影的起点和终点来判别是否为邻居,出于这个原因,投影的起点和终点序列被存储在目录中。接触发现算法接触发现算法边界盒法边界盒法二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法对于3D体系,必须把边界盒子在三个坐标轴上投影,所以生成三个目录。如果一个粒子的边界盒在某个坐标轴上投影的起点和终点间包含另外一个粒子边界盒投影的起点、终点或起点和终点,就说明这两个粒子的边界盒在该坐标轴上的投影发生了重叠。如果两个边界盒的投影在每一个坐标轴上都发生重叠,那么就说明这两个边界盒发生了接触。接触发现算法接触发现算法边界盒法边界盒法二二 基本原理基本原理-球形颗粒元离散元法球形颗粒元离散元法存储目录的更新存储目录的更新检查一个边界盒在每一个坐标轴的投影是否在另外一个边界盒投影的起点和终点之间仍旧需要花费很多计算时间。但是,尽管这些目录在每个计算时间步长上都必须被更新,所用的计算时间仍旧可以被减少到与体系中粒子数量成比例的量级,因为在每
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1