摘要

考虑理想弹塑性弯扭-Mohr-Coulomb-最大拉应力破坏准则,以及拉伸失效、压剪失效、转动屈服和扭转屈服四种失效模式,提出一种新的离散元弹塑性本构模型,其植入二次开发的三维颗粒流程序,数值试验研究宏细观参数关系,模拟Lac du Bonnet(LDB)花岗岩试样压剪和拉伸强度试验。研究表明:颗粒转动弹塑性弯扭矩对宏观内摩擦角和粘聚力具有显著影响;宏观粘聚力随细观粘聚力和屈服弯扭矩增加而趋于恒定;宏观内摩擦角随屈服弯扭矩和细观摩擦系数增加而增加;岩石宏观单轴抗拉强度随细观极限拉应力增加而趋于恒定;宏观弹性模量和泊松比主要受细观弹性模量和刚度比影响。弹塑性颗粒流本构模型模拟的LDB花岗岩的强度参数和变形参数以及压拉强度比,裂缝分布与室内试验结果吻合较好。LDB花岗岩单轴压剪模拟试验在达到峰值应力前,试样内部裂缝以细观拉伸裂缝为主,从试样端部向中部发展成近似X状;在接近或超过峰值应力后,试样内部开始出现压剪裂缝、转动屈服裂缝和扭转屈服裂缝,从试样一端向另一端发展成单斜状,其倾角与宏观破坏倾角接近,发展成宏观压剪破坏面,从而揭示了宏观裂缝形成的细观机理。