椭球颗粒体系宏、细观特性的3维离散元分析

颗粒材料由大量不同形状的离散颗粒组成,其力学特性与构成颗粒体系的颗粒形状密切相关。目前,关于颗粒形状对颗粒材料剪切强度的影响已有大量的物理试验或数值模拟试验研究,而非球形颗粒材料在复杂应力路径下的宏、细观力学特性认识还不够充分。本文采用离散单元法(DEM)对按一定级配生成的不同伸长率的椭球颗粒体系进行等σ_(3)(第3主应力)及等b(中主应力比)的真三轴加载数值模拟试验,研究椭球颗粒体系的3维宏、细观力学特性,并验证常用3维强度准则对椭球颗粒体系的适用性。结果表明:在3维加载路径下,不同于圆球颗粒体系宏、细观演化规律的一致性,椭球颗粒体系的宏观应力面与细观组构面演化规律不一致,宏观应力比与强接触网络组构偏值/均值比线性不为1,这是因为椭球颗粒体系弱接触网络对宏观应力的贡献;椭球颗粒体系宏观体应变、细观法向接触力分布、接触数目及配位数分布等宏、细观指标与加载路径无关;随着椭球伸长率LAR增大,Mohr–Coulomb准则、Lade–Duncan准则、Mastsuoka–Nakai准则及施维成3维强度准则对φ_(max)(峰值内摩擦角)–b关系预测准确度逐渐增大,其中,Lade–Duncan准则能更好地拟合不同加载路径下椭球颗粒体系的数据点。研究为真实应力状态下椭球颗粒材料的宏、细观特性的认识提供了一定的理论基础。

粒间摩擦对岩土颗粒材料三维力学行为的影响机制

颗粒间摩擦是影响岩土颗粒材料力学特性的重要因素之一,现有研究主要集中在粒间摩擦对颗粒集合体堆积、宏观剪切强度等的影响,而其在复杂应力路径下对岩土颗粒材料的影响机制尚不明确。采用离散单元法(DEM)进行了颗粒材料的真三轴加载路径数值试验,通过改变颗粒间摩擦系数μ来研究粒间摩擦对颗粒材料宏观力学特性的影响。对比了4种三维强度准则,发现当0.2<μ≤0.5时,Lade-Duncan和Matsuoka-Nakai准则能较好地拟合不同加载路径下的数据点。采用平均接触力将颗粒体系的接触网络划分为强、弱两个子网络,分析了强、弱接触网络的配位数分布、接触力和组构各向异性等。结果表明在峰值剪应力状态,随着颗粒间滑动摩擦系数μ的增加,强接触网络中参与形成"力链"的颗粒数目基本不变,但强接触网络中颗粒间法向接触力明显增加以及法向接触力各向异性加强,这是导致宏观抗剪强度增大的主要原因。而弱接触网络的配位数分布随μ值明显变化,其对颗粒体系剪胀程度的增加贡献显著。