我出来学习只办1件事：1.学习PFC_pfc中平行粘结法向强度-CSDN博客

来源公众号：材料科研SCI

PFC2D模拟

层状岩体PFC模拟的细观参数分析

【目的】研究层状岩体的非均质性，通过数值模拟的方法建立模型。

【方法】基于复合岩体思想，在PFC2D中加入等间距的层理弱面以模拟层岩的各向异性。

【结论】数值模型可以模拟层状岩体的非均质性。

细观参数的敏感性分析

本研究选用岩石模拟中常用的平行黏结模型，前人研究表明，细观参数的取值与数值试件的宏观力学参数具有一定的关系。数值试件弹性模量主要受颗粒弹性模量和平行黏结弹性模量的影响，且成正比例关系；峰值强度主要受平行黏结法向强度和切向强度的影响，随着黏结强度的提高，峰值强度也随之升高；泊松比主要受颗粒的黏结法向刚度和切向刚度比值的影响，比值越大，泊松比越大，但是鲜有人研究细观参数对三轴力学特性的影响。

因此，本研究进行数值试件三轴压缩模拟，研究细观参数对黏聚力和摩擦角的影响。采用PFC2D生成宽（W）为50mm，高（H）为100mm的数值试件，并进行三轴压缩模拟，颗粒的基本细观参数如表1所示，如无特殊说明，颗粒间的黏结距离值均取0。

细观参数对力学性能的影响

1、颗粒半径的影响

颗粒半径是PFC2D数值试件的基本参数，选择试件的宽度（W）与平均颗粒半径（R）的比值（S=W/R）为10、20、30、50、100、150，其余细观参数取值如表1所示，分别进行三轴压缩模拟。

模拟结果如图1所示。

由图1可以看出，当S从10增加到20时，黏聚力和摩擦角均急剧降低；当S从20增加到50时，黏聚力和摩擦角均出现上升现象；当S大于50时，黏聚力和摩擦角的变化幅度不大，基本稳定为一个恒定值。根据REV理论，采用PFC2D进行数值模拟时，试件的最小尺寸不得小于平均颗粒半径的50倍。

2、摩擦系数的影响

为研究颗粒间摩擦系数对数值试件三轴力学特性的影响，分别取颗粒间的摩擦系数f为0.1、0.2、0.3、0.5、1.0、2.0、3.0、4.0，其余细观参数取值如表1所示，分别进行三轴压缩模拟。

模拟结果如图2所示。

由图2可以看出，随着摩擦系数的增加，试件的摩擦角逐渐升高，上升速率逐渐降低，当摩擦系数大于2.0时，摩擦角趋于稳定；与此相反，随着摩擦系数的增加，试件的黏聚力逐渐降低，且下降速率逐渐降低，当摩擦系数大于2.0时，黏聚力亦趋于稳定。通过回归分析，黏聚力c和摩擦角φ与摩擦系数f的关系为式（1）、式（2）。