[发明专利]一种离散元模拟中的细观参数标定方法

说明书

本发明公开了一种离散元模拟中的细观参数标定方法，法向刚度、颗粒摩擦系数和黏性阻尼率三种细观参数，具体步骤如下：步骤1，标定法向刚度Kn：步骤1.1：先假定一个法向刚度；步骤1.1.1：再做数值侧限压缩实验并得到竖向压力‑位移曲线；步骤1.2：判断步骤1.1.1得到的竖向压力‑位移曲线是否和参考实验(室内实验)吻合；步骤1.3：如吻合，即假定法向刚度就是颗粒材料的真实法向刚度，确定法向刚度Kn，如不吻合，改变法向刚度的值并重复步骤1.1.1与步骤1.2；步骤2，标定颗粒摩擦系数：步骤2.1：采用步骤1.3得出的法向刚度Kn。优点在于：将主要标定对象有：法向刚度，颗粒之间摩擦系数和黏性阻尼率能较好地进行标定，方便进行离散元技术的研究。

技术领域

本发明涉及离散元数值仿真模拟技术领域，尤其涉及一种离散元模拟中的细观参数标定方法。

背景技术

离散元法是专门用来解决不连续介质问题的数值模拟方法。该方法把节理岩体视为由离散的岩块和岩块间的节理面所组成,允许岩块平移、转动和变形,而节理面可被压缩、分离或滑动。因此,岩体被看作一种不连续的离散介质。其内部可存在大位移、旋转和滑动乃至块体的分离,从而可以较真实地模拟节理岩体中的非线性大变形特征。离散元法的一般求解过程为:将求解空间离散为离散元单元阵,并根据实际问题用合理的连接元件将相邻两单元连接起来；单元间相对位移是基本变量,由力与相对位移的关系可得到两单元间法向和切向的作用力；对单元在各个方向上与其它单元间的作用力以及其它物理场对单元作用所引起的外力求合力和合力矩,根据牛顿运动第二定律可以求得单元的加速度；对其进行时间积分,进而得到单元的速度和位移。从而得到所有单元在任意时刻的速度、加速度、角速度、线位移和转角等物理量。

离散元技术在岩土、矿冶、农业、食品、化工、制药和环境等领域有广泛地应用，可分为分选、凝聚、混合、装填和压制、推铲、储运、粉碎、爆破、流态化等过程。颗粒离散元法在上述领域均有不少应用：料仓卸料过程的模拟；堆积、装填和压制；颗粒混合过程的模拟。离散元法在岩土工程、地质工程和能源开采领域也具有广泛的应用价值。

在离散元数值仿真模拟中，细观参数的标定是最关键的问题之一，主要细观参数为：法向刚度、颗粒之间摩擦系数和黏性阻尼率。这三个量一旦确定，离散元数值模拟的结果就可以很好的代表实际试验材料的力学试验。而现有技术中没有好的方法来标定。

为解决上述问题，我们提出了一种离散元模拟中的细观参数标定方法。

发明内容

本发明的目的是为了解决背景技术中的问题，而提出的一种离散元模拟中的细观参数标定方法。

为了实现上述目的，本发明采用了如下技术方案：一种离散元模拟中的细观参数标定方法，法向刚度、颗粒摩擦系数和黏性阻尼率三种细观参数，具体步骤如下：

步骤1，标定法向刚度Kn：

步骤1.1：先假定一个法向刚度；

步骤1.1.1：再做数值侧限压缩实验并得到竖向压力-位移曲线；

步骤1.2：判断步骤1.1.1得到的竖向压力-位移曲线是否和参考实验(室内实验)吻合；

步骤1.3：如吻合，即假定法向刚度就是颗粒材料的真实法向刚度，确定法向刚度Kn，如不吻合，改变法向刚度的值并重复步骤1.1.1与步骤1.2；

步骤2，标定颗粒摩擦系数：

步骤2.1：采用步骤1.3得出的法向刚度Kn，先假定一个摩擦系数；

步骤2.1.1：再做数值休止角试验并得到数值休止角；

步骤2.2：判断数值休止角是否和实验测定的休止角一致；

步骤2.3：如一致，即假定摩擦系数就是材颗粒真实的摩擦系数，如果不一致，改变颗粒摩擦系数的值并重复步骤2.1.1与步骤2.2；