基于微平面模型和正则化的脆性岩石损伤破裂模拟方法

doi:10.12068/j.issn.1005-3026.2022.08.011

东北大学学报（自然科学版） ›› 2022, Vol. 43 ›› Issue (8): 1141-1148.DOI: 10.12068/j.issn.1005-3026.2022.08.011

基于微平面模型和正则化的脆性岩石损伤破裂模拟方法

袁阳¹，徐涛¹，周广磊²，勒治华¹

(1. 东北大学资源与土木工程学院，辽宁沈阳110819； 2. 山东科技大学能源与矿业工程学院，山东青岛266590)

修回日期:2021-05-01 接受日期:2021-05-01 发布日期:2022-08-11
通讯作者: 袁阳
作者简介:袁阳(1993-)，男，江西新余人，东北大学博士研究生; 徐涛(1975-)，男，湖北随州人，东北大学教授，博士生导师.
基金资助:
国家自然科学基金资助项目(51974062，42172312); “十三五”国家重点研发计划项目(2017YFC1503101).

Simulation Method of Damage and Fracture for Brittle Rock Based on Microplane Model and Regularization

YUAN Yang¹， XU Tao¹， ZHOU Guang-lei²， LE Zhi-hua¹

1. School of Resources & Civil Engineering， Northeastern University， Shenyang 110819， China; 2. College of Energy and Mining Engineering， Shandong University of Science and Technology， Qingdao 266590， China.

Revised:2021-05-01 Accepted:2021-05-01 Published:2022-08-11
Contact: XU Tao
About author:-
Supported by:
-

摘要/Abstract

摘要： 为解决有限元计算在模拟脆性岩石损伤过程方面的网格依赖性问题，本文基于微平面模型和正则化方法建立了改进的微平面损伤模型，将单轴压缩条件下的岩石实验结果与数值模型进行了对比.引入不同非均质系数对三点弯曲岩石试样的损伤发展进行模拟，结果表明非均质系数越高损伤区域越集中.此外，通过引入高阶梯度正则化方法对不同特征长度的三点弯曲实验进行数值模拟.特征长度为1mm时，在不同网格尺寸下的正则化微平面模型的损伤演化结果基本相同，而没有采用正则化方法的微平面损伤模型在不同网格尺寸下出现明显的网格依赖性.最后，通过对比普通模型和正则化模型的载荷-位移曲线，发现网格密度对采用高阶梯度正则化后的微平面损伤模型损伤演化影响较低.

关键词: 微平面模型；非均质；损伤演化；正则化；数值模拟

Abstract: To solve the mesh dependence of finite element method in simulating the failure of brittle rock， an improved microplane damage model was built based on the microplane model and regularization method. And the numerical simulations were compared with the experimental results under uniaxial compression test. Heterogeneous indexes were introduced for the damage process simulation of three points bending tests. The results showed that the damage area was more concentrated at the high heterogeneous index. In addition， the three points bending tests with different characteristic length were simulated by regularization based on higher-order deformation gradients. The results showed that the damage evolution of microplane damage models with regularization are analogous when the characteristic length is 1mm， while the microplane damage model without regularization has a significant mesh-dependence at different mesh size. Finally， the load-displacement relationship of normal model and regularization model was contrasted. Mesh density has little effect on the microplane damage model based on regularization with higher-order deformation gradients.

Key words: microplane model; heterogeneous; damage evolution; regularization; numerical simulation

中图分类号:

TD853

袁阳，徐涛，周广磊，勒治华. 基于微平面模型和正则化的脆性岩石损伤破裂模拟方法[J]. 东北大学学报（自然科学版）, 2022, 43(8): 1141-1148.

YUAN Yang， XU Tao， ZHOU Guang-lei， LE Zhi-hua. Simulation Method of Damage and Fracture for Brittle Rock Based on Microplane Model and Regularization[J]. Journal of Northeastern University(Natural Science), 2022, 43(8): 1141-1148.

参考文献

[1]冯西桥.脆性材料的细观损伤理论和损伤结构的安定分析［D］.北京:清华大学，1995.(Feng Xi-qiao.Micro-failure theory for brittle materials and shakedown analysis of structures with damage［D］.Beijing:Tsinghua University，1995.)
[2]Batdorf S B，Budiansky B.A mathematical theory of plasticity based on the concept of slip［J］.Soviet Applied Mechanics，1949，12(12):1102-1112.
[3]Bazant Z，Caner F，Carol I，et al.Microplane model M4 for concrete I:formulation with work-conjugate deviatoric stress［J］.Journal of Engineering Mechanics，2000，126(9):944-953.
[4]Caner F，Bazant Z.Microplane model M7 for plain concrete I:formulation［J］.Journal of Engineering Mechanics，2013，139(12):1714-1723.
[5]Carol I，Bazant Z.Damage and plasticity in microplane theory［J］.International Journal of Solids and Structures，1997，34(29):3807-3835.
[6]陈新，杨强，李德建.岩体裂隙网络各向异性损伤力学效应研究［M］.北京:科学出版社，2016.(Chen Xin，Yang Qiang，Li De-jian.Study on anisotropic damage mechanics effect of fracture network in rock mass［M］.Beijing:Science Press，2016.)
[7]Murakami S，Liu Y.Mesh-dependence in local approach to creep fracture［J］.International Journal of Damage Mechanics，1995，4:230-250.
[8]Liu Y，Murakami S，Yamada T，et al.Mesh dependence in local approach to creep fracture analysis of perforated plate［J］.Transactions of the Japan Society of Mechanical Engineers，1995，61:2030-2036.
[9]Borst R，Sluys L J，Muhlhaus H，et al.Fundamental issues in finite element analysis of localization of deformation［J］.Engineering Computations，1993，10(2):99-121.
[10]Mondal S，Olsen-Kettle L，Gross L.Regularization of continuum damage mechanics models for 3-D brittle materials using implicit gradient enhancement［J］.Computers and Geotechnics，2020，122:103505.
[11]崔焕平，崔燕平，王宗敏.混凝土非线性有限元分析中的网格尺寸效应［J］.混凝土，2007(6):27-29.(Cui Huan-ping，Cui Yan-ping，Wang Zong-min.Mesh size effect in nonlinear finite element analysis of concrete［J］.Concrete，2007(6):27-29.)
[12]Bazant Z，Pijaudier-Cabot G.Nonlocal continuum damage，localization instability and convergence［J］.Journal of Applied Mechanics，1988，55(2):287-293.
[13]Bazant Z，Jirásek M.Nonlocal integral formulations of plasticity and damage:survey of progress［J］.Journal of Engineering Mechanics，2002，128(11):1129-1239.
[14]Bazant Z，Obolt J.Nonlocal microplane model for fracture，damage，and size effect in structures［J］.Journal of Engineering Mechanics，1990，116(11):2485-2505.
[15]崔溦，杨娜娜，宋慧芳.基于非局部微平面模型M7的混凝土非线性有限元分析［J］.建筑结构学报，2017，38(2):126-133.(Cui Wei，Yang Na-na，Song Hui-fang.Nonlinear finite element analysis of concrete structure based on nonlocal microplane model M7［J］.Journal of Building Structures，2017，38(2):126-133.)
[16]杨娜娜.基于非局部理论的混凝土M7微平面本构模型研究及应用［D］.天津:天津大学，2016.(Yang Na-na.Research and application of microplane constitutive model M7 based on nonlocal theory for plain concrete［D］.Tianjin:Tianjin University，2016.)
[17]Bazant Z，Di L G.Nonlocal microplane model with strain-softening yield limits［J］.International Journal of Solids and Structures，2004，41(24/25):7209-7240.
[18]Caner F，Bazant Z.Microplane model M7 for plain concrete II:calibration and verification［J］.Journal of Engineering Mechanics，2013，139(12):1724-1735.
[19]Xu T，Tang C A，Zhao J，et al.Modelling the time-dependent rheological behaviour of heterogeneous brittle rocks［J］.Geophysical Journal International，2012，189(3):1781-1796.
[20]Chen X，Bazant Z.Microplane damage model for jointed rock masses［J］.International Journal for Numerical and Analytical Methods in Geomechanics，2014，38(14):1431-1452.
[21]Bazant Z，Oh B.Microplane model for progressive fracture of concrete and rock［J］.Journal of Engineering Mechanics，1985，111(4):559-582.
[22]Bazant Z，Xiang Y，Prat P.Microplane model for concrete I:stress-strain boundaries and finite strain［J］.Journal of Engineering Mechanics，1996，122(3):245-254.
[23]Kachanov L M.Rupture time under creep conditions［J］.International Journal of Fracture，1999，97(1/2/3/4):11-18.
[24]朱万成，唐春安，杨天鸿，等.岩石破裂过程分析用(RFPA~(2D))系统的细观单元本构关系及验证［J］.岩石力学与工程学报，2003，22(1):24-29.(Zhu Wan-cheng，Tang Chun-an，Yang Tian-hong，et al.Constitutive relationship of mesoscopic elements used in RFPA~2D and its validations［J］.Chinese Journal of Rock Mechanics and Engineering，2003，22(1):24-29.)
[25]梁正召.三维条件下的岩石破裂过程分析及其数值试验方法研究［D］.沈阳:东北大学，2005.(Liang Zheng-zhao.Three-dimensional failure process analysis of rock and associated numerical tests［D］.Shenyang:Norteastern University，2005.)
[26]Xu T，Yuan Y，Heap M J，et al.Microwave-assisted damage and fracturing of hard rocks and its implications for effective mineral resources recovery［J］.Minerals Engineering，2021，160:106663.
[27]Yuan Y，Xu T，Heap M J，et al.A three-dimensional mesoscale model for progressive time-dependent deformation and fracturing of brittle rock with application to slope stability［J］.Computers and Geotechnics，2021，135:104160.

基于微平面模型和正则化的脆性岩石损伤破裂模拟方法

Simulation Method of Damage and Fracture for Brittle Rock Based on Microplane Model and Regularization

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	勒治华，于庆磊，郑浩田，曹永胜. 侧限条件下充填散体与岩柱相互作用机理[J]. 东北大学学报（自然科学版）, 2021, 42(1): 124-130.
[2]	姜立春，谢波. 沉积型铝土矿采空区矿柱-顶板支撑系统滑移突变失稳分析[J]. 东北大学学报:自然科学版, 2020, 41(12): 1767-1774.
[3]	任凤玉，刘洋，何荣兴，张晶. 放矿条件下塌陷区内尾砂穿流特性[J]. 东北大学学报:自然科学版, 2020, 41(6): 858-863.
[4]	温震江，高谦，王永定，何建元. 不同浓度料浆流变特性与混合骨料级配相关性试验[J]. 东北大学学报:自然科学版, 2020, 41(5): 642-648.
[5]	闫保旭，朱万成，侯晨，贾瀚文. 充填体应力分布理论分析及数值模拟对比研究[J]. 东北大学学报:自然科学版, 2019, 40(12): 1773-1778.
[6]	谭宝会，张志贵，何荣兴，朱强. 放矿椭球体排列理论的合理性探讨及实验研究[J]. 东北大学学报:自然科学版, 2019, 40(7): 1014-1019.
[7]	任凤玉，谭宝会，付煜，朱强. 诱导冒落法回采挂帮矿引发的边坡失稳[J]. 东北大学学报:自然科学版, 2019, 40(2): 273-278.
[8]	张飞，杨天鸿，胡高建. 复杂应力扰动下围岩稳定性评价与采场参数优化[J]. 东北大学学报:自然科学版, 2018, 39(5): 699-704.
[9]	任思潼，顾晓薇，付煜，李传迎. 露天边帮下紧邻尾矿库采空区地表充填治理方法[J]. 东北大学学报:自然科学版, 2018, 39(5): 722-725.
[10]	刘洋，任凤玉，何荣兴，曹建立. 基于放矿下临界散体柱理论的地表塌陷范围预测[J]. 东北大学学报:自然科学版, 2018, 39(3): 416-420.
[11]	谭宝会，任凤玉，何荣兴，张志贵. 垂直分区与组合放矿在倾斜中厚矿体中的应用[J]. 东北大学学报:自然科学版, 2018, 39(1): 133-137.
[12]	雷刚，李元辉，张志贵，李睿. 无底柱分段崩落法矿石覆盖层回收实验[J]. 东北大学学报:自然科学版, 2017, 38(3): 400-405.
[13]	柯愈贤，王新民，张钦礼，刘恩彦. 基于全尾砂充填体非线性本构模型的深井充填强度指标[J]. 东北大学学报:自然科学版, 2017, 38(2): 280-283.
[14]	由希，张国建. 覆岩散体粒径空间分布对岩石混入的影响[J]. 东北大学学报:自然科学版, 2016, 37(11): 1604-1608.
[15]	姜立春，肖康，吴爱祥，李令鑫. 基于结构模态的纵向采空区群爆破动力响应[J]. 东北大学学报:自然科学版, 2016, 37(11): 1615-1619.