考虑岩石微缺陷影响的损伤本构模型.pdf

第31卷第3期2023 年 6 月Vol.31 No.3Jun.，2023Gold Science and Technology507考虑岩石微缺陷影响的损伤本构模型刘志祥，晏孟洋*，张双侠，熊帅，王凯中南大学资源与安全工程学院，湖南 长沙 410083摘 要：为准确描述岩石应力应变曲线非线性过程，提出了考虑岩石微缺陷影响的损伤本构模型。首先对含微缺陷岩石进行分析，把含微缺陷岩石抽象成不含微缺陷的岩石骨架部分和微缺陷部分。岩石微缺陷包括初始空隙和岩石受载新增的微裂纹，微缺陷只能产生应变，不能承受应力。初始空隙产生的应变反映压密阶段的非线性特征，新增的微裂纹产生的应变反映岩石峰后阶段应变软化，用修正系数b来表示新增微缺陷的影响，假设岩石骨架部分的损伤符合Weibull概率分布，最后推导出基于微裂纹的损伤本构模型，并给出参数Vm、n、F0、m和b的确定方法，对模型参数进行讨论。用砂岩和苏长岩试验数据进行验证，结果显示试验数据与理论结果相吻合。关键词：含微缺陷岩石；本构模型；微缺陷；压密阶段；应变软化；修正系数中图分类号：TD853 文献标志码：A 文章编号：1005-2518（2023）03-0507-09 DOI：10.11872/j.issn.1005-2518.2023.03.164引用格式：LIU Zhixiang，YAN Mengyang，ZHANG Shuangxia，et al.Damage Constitutive Model Considering the Effect of Rock Microdefects J.Gold Science and Technology，2023，31（3）：507-515.刘志祥，晏孟洋，张双侠，等.考虑岩石微缺陷影响的损伤本构模型 J.黄金科学技术，2023，31（3）：507-515.岩石是一种天然的地质材料，在成岩过程中，其内部形成了不同形状、尺寸、方向、数目以及不同种类的空隙，称为微缺陷，由于微缺陷的影响，岩石的变形呈现出非线性。随着采矿工程的不断发展，进入深部工程后，岩石受“三高”（高地应力、高地温、高渗透压和扰动）的影响（何满潮等，2005），非线性工程问题不断增加，研究岩石损伤本构模型具有重要的实践意义。国内外学者从多个角度对岩石本构模型进行了研究。如：Chen et al.（2019）对传统统计损伤模型进行修正；Chen et al.（2018）比较了幂函数与Weibull分布函数的统计损伤模型；Gao et al.（2018）提出了考虑热力耦合的本构模型；Miao et al.（2018）建立了水化学环境下的本构模型；Hu et al.（2018）构建了硬岩的非线性本构模型；Huang et al.（2018）建立了弹塑性损伤本构模型；曹文贵等（2003）建立了正态分布损伤本构模型；李海潮等（2017）对应变等价性假设的不足进行了修正，得到新的岩石损伤本构模型；Liu et al.（2018）使用Logistic方程来拟合岩石的损伤过程，并以此建立了岩石损伤本构模型；贾宝新等（2022）提出了基于 Hoek-Brown 准则的高温岩石损伤本构。上述岩石损伤本构模型与特定环境下的岩石应力应变曲线相吻合，但是对应力应变的压密阶段描述略有不足，而含有初始微缺陷的脆性岩石一般有压密阶段。张超等（2020）对两类初始空隙进行分析，建立了考虑初始空隙的岩石损伤本构模型。曹文贵等（2016）考虑空隙压密效应，提出了相应的本构模型。张超等（2021）基于应变等效假设的弹性模量法，构建了岩石损伤本构模型。Cao et al.（2010）通过确定孔隙 收稿日期：2022-11-01；修订日期：2023-03-30基金项目：国家重点研发计划项目“深部厚大矿体大参数高效连续开采技术”（编号：2022YFC290410）和国家自然科学基金项目“海底金属矿开采充填体约束矿柱群力学模型构建与混沌破坏机制”（编号：51974359）联合资助作者简介：刘志祥（1967-），男，湖南宁乡人，博士生导师，从事岩石力学与采矿工程研究工作。*通信作者：晏孟洋（1999-），男，湖南邵阳人，硕士研究生，从事延时损伤力学研究工作。Vol.31 No.3 Jun.，2023采选技术与矿山管理508率来分析空隙和体积变化对岩石行为的影响，然后建立相关的损伤模型。这些学者建立的本构模型考虑到微缺陷会导致岩石产生压密阶段，但忽略了新增微缺陷对脆性岩石应变软化阶段的影响。由此可见，上述学者的研究对于岩石损伤本构模型的发展和进步作出了重要的贡献，提出的岩石损伤本构模型已考虑到微缺陷会导致岩石产生压密阶段，具有很强的针对性，但对于脆性岩石全应力应变曲线的描述仍有一定局限性。因此，本文针对含有初始微缺陷脆性岩石，重点考虑微缺陷对压密阶段和应变软化阶段的影响，进一步完善脆性岩石全应力应变过程损伤本构模型。1 损伤本构模型的建立1.1 含微缺陷脆性岩石变形分析为建立含微缺陷岩石的应力应变全过程的损伤本构模型，选取一体积单元进行分析，如图1所示，该体积单元包括岩石骨架和微缺陷。由于微缺陷的多样性，微缺陷的大小、形状和数目各不相同。对单个微缺陷进行分析，理论十分复杂，而且不能代表所有的微缺陷，所以把岩石中的微缺陷集中起来，将岩石划分为岩石骨架和由微缺陷组成的空隙部分。岩石在外荷载ij的作用下开始发生变形，第一阶段为岩石中的初始空隙闭合阶段。此时岩石中的微缺陷被压密，岩石变形为非线性变形，当初始微缺陷完全关闭后，岩石变形转变为弹性变形。第二阶段，当应力超过萌生应力后，岩石骨架开始损伤，新的微裂纹开始形成，组成新的空隙部分。第三阶段，岩石骨架完全损伤。如图2所示。1.2 微裂纹萌生前岩石本构模型微缺陷会对岩石的力学性质和损伤行为产生重要影响，偏应力作用下的含微缺陷脆性岩石压密阶段将十分明显，含初始微缺陷材料的压缩性明显大于不含初始微缺陷的相同固体材料。初始微缺陷是岩石不可缺少的一部分，假设岩石骨架部分的变形只与其本身性质有关，空隙部分的应变与外荷载ij有关。选取体积单元进行分析。（1）岩石骨架部分变形分析在不考虑蠕变的弹塑性变形中，岩石骨架部分的本构关系可采用经典的弹塑性本构关系式：Rij=eij+pij（1）eij=1+Eij-ERRij（2）pij=f（3）式中：Rij（i、j=1，2，3）为岩石骨架部分的应变；pij为岩石骨架部分的塑性变形；eij为岩石骨架部分的弹性应变；为岩石材料泊松比；E为岩石弹性模量；f为屈服面；为材料常数。由于该本构模型主要针对脆性岩石，不考虑塑性变形，只考虑弹性变形的影响，即=0，可得：pij=0（4）图1岩石微缺陷集中示意图Fig.1Schematic diagram of the concentration of micro-defects in rocks图2岩石损伤与微缺陷变化示意图Fig.2Schematic diagram of rock damage and micro-defect variation2023 年 6 月 第 31 卷 第 3 期509刘志祥等：考虑岩石微缺陷影响的损伤本构模型Rij=eij（5）（2）微缺陷部分变形分析选取岩石的体积单元进行分析，该单元包括足够多的微缺陷和细观结构，所有微缺陷部分所引起的应变之和为cij。cij=k=1Ni()ij（6）式中：N为单元中微缺陷总数；i()ij为第个微缺陷所引起的应变。对于单个的微缺陷的闭合，已知微缺陷的形状尺寸和受力状态的情况下，能得到特定情况下微缺陷闭合时所引起的应变的理论解。许多学者通过对含有点缺陷、线缺陷和体缺陷的弹性介质体进行研究，将各种细观参数作为随机变量，通过平均化的方法得到岩石材料的总体有效参数。但这种理论提供的方法十分复杂。为避免分析材料中微缺陷形状、尺寸和方位的复杂性，从全部微缺陷的角度出发，将体积单元中的所有微缺陷集中考虑，考虑微缺陷整体在外荷载作用下对材料整体的应力应变关系的影响，假设微缺陷部分所引起的应变cij与ij存在函数关系：cij=f(ij)（7）当i=j=1时，cij就变为c11，表示微缺陷部分所引起的轴向应变，所以 c11与偏应力存在函数关系。由于岩石骨架部分在微裂纹萌生前只发生弹性变形，所以岩石骨架部分的应变与应力呈线性关系。研究表明，当轴向应力较小时，微缺陷部分产生的轴向应变c11是随轴向应力的增大而增加的；微缺陷完全闭合时，c11不再随轴向应力的增大而增加，而是趋于一个稳定值。微缺陷产生的轴向应变与轴向应力呈负指数关系（Peng et al.，2015）：c11=Vm 1-exp(-11-33n)（8）式中：c11为初始微缺陷引起的轴向应变；Vm和n为常量参数；11和33为主应力。（3）微裂纹萌生前岩石本构模型确定岩石的轴向总应变11等于岩石骨架部分的轴向应变R11与岩石微缺陷部分应变c11之和。结合式（1）、式（2）、式（4）、式（8）可得：11=11-33E+Vm 1-exp(-11-33n)（9）为验证上述微裂纹萌生前岩石本构模型的正确性，引用文献（Yumlu et al.，1995）中砂岩和苏长岩的常规三轴压缩试验数据。由图 3 可知，当初始微缺陷闭合之前（B 点之前），由于微缺陷影响，岩石应力应变曲线为非线性，岩石轴向应变11等于微缺陷部分应变c11与岩石骨架部分应变R11之和。B点之后微缺陷产生的应变几乎为0，称cc为微缺陷闭合应变，B点为初始微缺陷闭合点。BC段为弹性阶段，曲线呈线性变化，C点之后初始微裂纹开始萌生，岩石出现损伤，称aa为微裂纹萌生应变，C点为微裂纹萌生点。由图 4曲线拟合可得模型参数 Vm和 n，从图 4可以看出，试验数据与拟合效果较好，在相同应力条件下，砂岩的微缺陷闭合应变大于苏长岩的微缺陷闭合应变，砂岩与苏长岩微缺陷模型参数结果见表1。将模型参数代入式（9）中，可得微裂纹萌生前岩石的本构表达式，如图5所示。由图5可知，理论曲线与试验曲线基本吻合，验证微裂纹萌生前岩石本构模型的合理性。1.3 微裂纹萌生后岩石本构模型当岩石所受的应力超过萌生应力时，岩石应力应变曲线进入屈服阶段，新的微裂纹开始出现，岩石出现不可逆损伤。由于压密阶段后，岩石中的初始微缺陷已闭合，岩石骨架屈服后又有新的微裂纹生成，新增微裂纹形成空隙部分，岩石骨架开始损伤。图3砂岩微裂纹萌生前应力应变曲线Fig.3Stress-strain curve of sandstone before microcrack initiationVol.31 No.3 Jun.，2023采选技术与矿山管理510基于等效应变假设（Lemaitre et al.，1985），材料的一维损伤本构可表示为=*E=(1-D)E（10）式中：和分别为应变和应力；*为有效应力；D为损伤变量。对式（12）进行等量变换可得：=*(1-D)（11）=(1-D)E（12）Lemaitre et al.（1985）提出的等效应变假设认为岩石未损伤时，D=0，D=1时岩石完全损伤，实际上岩石完全损伤仍有承载能力。等效应变假设改进为11=*11(1-D)+D(r11+c11)（13）式中：*11为未损伤部分所受应力；r11为已损伤部分所受应力；c11为空隙部分所受应力。未损伤部分只发生弹性变形，表示为*11=E*11+(*22+*33)（14）含初始微缺陷脆性岩石完全损伤后，还有一定的残余强度。已损伤岩石部分所受应力r11等于岩石的残余应力，而岩石的残余应力可通过岩石应力应变曲线确定，即岩石完全破坏后，应变不断增加，应力趋于稳定时的应力值确定为残余强度。空隙部分所受应力为c11，由于空隙部分没有实际的承载能力，所以空隙部分所受应力c11为0。虽然空隙部分所承受的应力为0，但该部分会产生应变，对未损伤部分岩石所受的应变会产生影响。如上所述，岩石的本构模型可改进为11=*11(1-D)+Dr11（15）假设岩石只在轴向发生损伤且满足变形协调条件（曹文贵等，2012），即11=*11=r11（16）22=*22=r22（17）33=*33=r33（18）式中：*11为未损伤部分所受的轴向应变；*22和*33为未损伤部分所受纵向应力；r11为损伤部分所受的轴向应变；r22和r33为损伤部分所受纵向应力。上述公式未考虑到初始微缺陷的闭合与新增微裂纹对应变的影响。岩石的名义轴向应变11并非未损伤岩石的实际应变，未损伤部分轴向应变*11与损伤部分的轴向应变r11应等于岩石名义轴向应变 11减去岩石初始微缺陷闭合所导致的应变c11，可表示为*11=r11=11-c11（19）*11=11-Vm 1-exp()-11-33n（20）在图3中，由于岩石在初始微缺陷闭合点（B点）图5岩石微裂纹萌生前试验与理论应变应力曲线Fig.5Experimental and theoretical strain-stress curves before initiation of rock microcracks 图4岩石微缺陷闭合应变应力试验与拟合曲线Fig.4Experimental and fitted curves of closed strain-stress for rock microdefects表1微缺陷闭合模型参数Table 1Micro-defect closure model parameters岩石种类砂岩苏长岩E/GPa2885cc/1.370.67Vm/1.430.68n4.6414.71R2（拟合系数）0.94680.99242023 年 6 月 第 31 卷 第 3 期511刘志祥等：考虑岩石微缺陷影响的损伤本构模型之后，初始微缺陷完全闭合，初始微缺陷产生的应变c11不再增加，趋于一个常数，其值等于微缺陷闭合应变cc。因此，当新的微裂纹开始萌生之后，初始微缺陷对岩石应变的影响可用微缺陷闭合应变cc来体现，即：r11=11-cc（21）除了初始微缺陷会对岩石的应力应变关系产生影响之外，在C点之后新的微裂纹开始产生、扩张、汇合，也会对应变产生影响。岩石萌生的新的微裂纹所产生的空隙部分没有承载能力，但会产生应变，需对未损伤岩石的名义轴向应变进行进一步修正。假设新增空隙部分的应变为nc，修正后为损伤岩石的名义轴向应力，表示为*11=r11=11-cc-nc（22）新增空隙部分的应变nc暂时无法通过试验手段获得，假设nc与ij呈函数关系，nc代表的是新增空隙部分所产生的应变，它会导致未损伤的岩石骨架部分实际所承受的应变减小，主要对岩石骨架部分实际所承受的应变起弱化作用。所以，引入一个修正参数b来代表新增微裂纹对未损伤骨架部分的弱化作用，即：*11=r11=b(11-cc)（23）式中：b为修正参数，取值范围为01。在应力应变曲线峰值点之前，岩石的微裂纹生成和扩展速率较慢，微裂纹形成的空隙部分对岩石未损伤部分应变影响较小，故在峰值点之前，不考虑新增空隙部分的影响。峰值点之后，进入应变软化阶段，岩石应力跌落较快。新增微裂纹的扩展、增加和汇合会加剧岩石的应变软化现象。新增空隙部分所占比例也急剧增加，新增微裂纹形成的空隙部分对未损伤岩石部分应变具有不可忽视的影响。所以用修正参数来修正岩石应变软化阶段的应力应变曲线。根据上述讨论，对三轴压缩下岩石的本构模型进行修正，可得：11=1-D(*11)E*11+(22+33)+D(*11)r11（24）2 模型参数的确定2.1 损伤变量的确定当11aa时，新的微裂纹生成，岩石主体骨架部分开始出现连续损伤，假设岩石的损伤过程符合Weibull概率分布函数，即：f(x)=mF0(xF0)m-1exp -(xF0)m（25）式中：m、F0为参数；x为自变量。则岩石的损伤变量函数为D(*11)=0*11f(*11)d(*11)=1-exp -(*11F0)m（26）D(*11)=0，0 11 aa 1-exp -b(11-cc)F0m，aa 11（27）得到考虑微缺陷脆性岩石应力应变全过程的损伤本构模型为 11=11-33E+Vm 1-exp()-11-33n，0 11 aa 11=1-D(*11)E*11+(22+33)+D(*11)r11，aa aa（35）式中：M11为模型计算的轴向应力；E11为试验数据的轴向应力；m为数据点数。所以b值确定为b=b(MSEmin)（36）3 损伤本构模型的验证3.1 实例验证为验证上述本构模型的合理性，选取文献中的砂岩和苏长岩常规三轴压缩试验进行验证，先根据已有数据确定砂岩和苏长岩的基本参数，已知砂岩的弹性模量E=28 GPa、泊松比=0.25。苏长岩的弹性模量E=85 GPa、泊松比=0.29。拟合曲线得模型参数Vm、n，修正参数b，通过计算得到F0和m。结果见表2、图6和图7。由图6和图7可知，理论曲线与试验曲线基本吻合，可以描述岩石压密阶段、弹性阶段、应变软化阶段和残余应变阶段的特征，其中，砂岩和苏长岩在不同围压的应力状态下取得不同修正b值，不同b值取值能够反映岩石损伤后产生的空隙部分对岩石峰后阶段应变软化的影响程度，从而验证了该损伤本构模型。3.2 参数分析如图8所示，当b=1时，没有考虑新增微缺陷对图6砂岩试验数据与理论曲线比较Fig.6Comparison between experimental data and theoretical curves of sandstone表2砂岩和苏长岩损伤演化方程参数Table 2Parameters of the damage evolution equation for sandstone and saprolite岩石类型砂岩苏长岩33/MPa358358r11/MPa25.147437.200062.945856.873269.0954101.7840p11/%0.740.850.980.440.500.58p11/MPa107.1221132.1709156.6792249.1216274.0948306.4413b=b（MSEmin）0.700.740.850.900.850.80F0/%0.52710.66600.92740.43070.47070.5192m8.76557.50653.42726.87806.85796.69922023 年 6 月 第 31 卷 第 3 期513刘志祥等：考虑岩石微缺陷影响的损伤本构模型脆性岩石应变软化阶段的影响，岩石未损伤部分所受的应变会偏大，导致应变软化阶段的应力偏大。当应变相同时，b值越小，其损伤本构理论应力越小，应力差表示空隙部分所产生的应变对岩石应力的弱化作用。当b=0.70时，砂岩试验曲线与理论曲线吻合良好，从而验证了用修正系数b来表示空隙部分对岩石未损伤部分岩石弱化作用的合理性。考虑b值对裂纹萌生点之后所有阶段的影响，可得不同b值条件下，围压为3 MPa时砂岩的损伤演化图。如图9所示，b值越小，损伤变量曲线斜率越大，说明损伤速度越快，即砂岩实际的损伤比不考虑空隙部分的影响时要大。4 结论首先对含初始微缺陷的脆性岩石进行变形分析，在加载过程中，岩石中的空隙逐渐被压密，直至完全压缩，岩石变形进入弹性阶段，该过程中岩石的应变由空隙部分产生的应变与岩石骨架部分应变组成，岩石不发生损伤。当应力大于萌生应力时，岩石出现新的微裂纹，并开始出现损伤，持续的加载导致岩石骨架部分完全损伤。以岩石的变形为基础，建立了考虑微缺陷影响的脆性岩石的损伤本构模型，对比试验数据与理论模型，得出以下结论：（1）初始微缺陷的闭合是导致岩石压密阶段应力应变曲线呈非线性变化的主要原因。初始微缺陷的闭合产生的应变与岩石的偏应力呈负指数分布，用砂岩和苏长岩的初始微缺陷应变与偏应力进行拟合，拟合系数大于0.94，得到的考虑初始微缺陷的本构模型在压密阶段与试验数据吻合较好。（2）岩石新增微裂纹形成的空隙部分会影响岩图7苏长岩试验数据与理论曲线比较Fig.7Comparison between experimental data and theoretical curves of norite图93=3 MPa时不同b值下砂岩损伤演化Fig.9Sandstone damage evolution for 3=3 MPa and different b values图83=3 MPa时不同b值砂岩理论与试验曲线Fig.8Theoretical and experimental curves of sandstone for 3=3 MPa and different b valuesVol.31 No.3 Jun.，2023采选技术与矿山管理514石的变形。新增微裂纹主要影响岩石的应变软化阶段，在峰值点之后，新增微缺陷会导致岩石的应力跌落加快，使其应变软化更加显著。通过对比发现当应变相同时，考虑新增微裂纹影响的损伤本构模型理论曲线与试验曲线更符合岩石的应变软化阶段变化规律，未考虑新增微裂纹影响的岩石应力会偏大，未考虑新增微裂纹影响的岩石损伤变量会偏小。（3）本文得到的岩石损伤本构模型能够描述岩石应力应变曲线的压密阶段、弹性阶段、屈服阶段、应变软化阶段和残余强度阶段，考虑到了脆性岩石应力应变曲线的所有过程，且理论曲线与试验数据吻合度较高。参考文献（References）：Cao W，Zhao H，Li X，2010.Statistical damage model with strain softening and hardening for rocks under the influence of voids and volume changesJ.Canadian Geotechnical Journal，47（8）：857-871.Cao Wengui，Zhang Chao，He Min，2016.A statistical damage simulation method for rock strain softening considering the characteristics of void compacting stageJ.Journal of Geotechnical Engineering，38（10）：1754-1761.Cao Wengui，Zhao Heng，Li Xiang，2012.A statistical damage simulation method for the whole process of rock deformation based on residual strength deformation phase characteristics J.Journal of Civil Engineering，45（6）：139-145.Cao Wengui，Zhao Minghua，Tang Xuejun，2003.Study on statistical damage simulation of rock fracture process J.Journal of Geotechnical Engineering，（2）：184-187.Chen K，Tang M，Guo Z，2019.Comparative study on three-dimensional conventional and modified statistical damage constitutive modelsJ.Multiscale and Multidisciplinary Modeling，Experiments and Design，2（4）：259-267.Chen S，Qiao C S，Ye Q，2018.Comparative study on three-dimensional statistical damage constitutive modified model of rock based on power function and Weibull distributionJ.Environmental Earth Sciences，77（3）：108.Gao M B，Li T B，Wei T，et al，2018.A statistical constitutive model considering deterioration for brittle rocks under a coupled thermal-mechanical condition J .Geofluids，2018：1-10.DOI：10.1155/2018/3269423.He Manchao，Xie Heping，Peng Suping，2005.Study on rock mechanics in deep mining engineering J.Chinese Journal of Rock Mechanics and Engineering，24（16）：2803-2813.Hu B，Yang S Q，Xu P，et al，2018.Nonlinear rheological damage model of hard rock J.Journal of Central South University，25（7）：1665-1677.Huang J，Zhao M，Du X，2018.An elasto-plastic damage model for rocks based on a new nonlinear strength criterionJ.Rock Mechanics and Rock Engineering，51（5）：1413-1429.Jia Baoxin，Chen Guodong，Liu Fengpu，2022.Ontogenetic model of rock damage at high temperature and its validationJ.Geotechnics，43（Supp.2）：63-73.Lemaitre J A，1985.Continuous damage mechanics model for ductile materials J.Journal of Engineering Materials and Technology，107（1）：83-89.Li Haichao，Zhang Sheng，2017.Rock damage model based on modified Lemaitre strain equivalence assumption J.Geotechnics，38（5）：1321-1326.Liang Mingchun，Miao Shengjun，Cai Meifeng，2021.Ontogenetic model of rock damage considering shear expansion properties and post-peak morphology J.Journal of Rock Mechanics and Engineering，40（12）：2392-2401.Liu D，He M，Cai M，2018.A damage model for modeling the complete stress-strain relations of brittle rocks under uniaxial compressionJ.International Journal of Damage Mechanics，27（7）：1000-1019.Miao S，Wang H，Cai M，2018.Damage constitutive model and variables of cracked rock in a hydro-chemical environmentJ.Arabian Journal of Geosciences，11：1-14.Peng J，Rong G，Cai M，2015.A model for characterizing crack closure effect of rocks J .Engineering Geology，189：48-57.Yumlu M，Ozbay M U，1995.A study of the behaviour of brittle rocks under plane strain and triaxial loading conditions J.International Journal of Rock Mechanics and Mining Sciences&Geomechanics Abstracts，32（7）：725-733.Zhang Chao，Chen Qiunan，Yang Qijun，2020.Simulation of the whole process of deformation damage of brittle rocks considering initial void closure and its influence J.Journal of Coal，45（3）：1044-1052.Zhang Chao，Yang Chuqing，Bai Yun，2021.Analysis of damage evolution of rock-like brittle materials and its modeling method J.Geotechnics，42（9）：2344-2354.曹文贵，张超，贺敏，2016.考虑空隙压密阶段特征的岩石应变软化统计损伤模拟方法 J.岩土工程学报，38（10）：1754-1761.曹文贵，赵衡，李翔，2012.基于残余强度变形阶段特征的岩2023 年 6 月 第 31 卷 第 3 期515刘志祥等：考虑岩石微缺陷影响的损伤本构模型石变形全过程统计损伤模拟方法 J.土木工程学报，45（6）：139-145.曹文贵，赵明华，唐学军，2003.岩石破裂过程的统计损伤模拟研究 J.岩土工程学报，（2）：184-187.何满潮，谢和平，彭苏萍，等，2005.深部开采岩体力学研究J.岩石力学与工程学报，24（16）：2803-2813.贾宝新，陈国栋，刘丰溥，2022.高温下岩石损伤本构模型及其验证 J.岩土力学，43（增2）：63-73.李海潮，张升，2017.基于修正Lemaitre应变等价性假设的岩石损伤模型 J.岩土力学，38（5）：1321-1326.梁明纯，苗胜军，蔡美峰，2021.考虑剪胀特性和峰后形态的岩石损伤本构模型 J.岩石力学与工程学报，40（12）：2392-2401.张超，陈秋南，杨期君，2020.考虑初始空隙闭合及其影响的脆性岩石变形破坏全过程模拟 J.煤炭学报，45（3）：1044-1052.张超，杨楚卿，白允，2021.岩石类脆性材料损伤演化分析及其模型方法研究 J.岩土力学，42（9）：2344-2354.Damage Constitutive Model Considering the Effect of Rock MicrodefectsLIU Zhixiang，YAN Mengyang，ZHANG Shuangxia，XIONG Shuai，WANG KaiSchool of Resources and Safety Engineering，Central South University，Changsha 410083，Hunan，ChinaAbstract：In order to accurately describe the whole process of rock stress-strain curve，a statistical damage contitutive model considering the effect of rock microdefects was proposed.Firstly，rocks containing microdefects are analysed and abstracted into a rock skeleton part and a defective part without microdefects.The rock microdefects include initial microdefects and new microdefects added to the rock by loading，which can only generate strain but not stress.The strain generated by the initial microdefects is negatively exponential to the stress during the compression-density stage，and the model parameters n and Vm are obtained by fitting the experimental curve of closed strain-stress for the microdefects.After the compression-density stage，the initial micro-defects are completely closed and the resulting strain is a constant.As additional micro-defects are mainly generated in large numbers in the post-peak stage，the effect of additional micro-defects on the post-peak stage is mainly considered.A variable b in the range of 0 to 1 is used to represent the weakening effect of the strain formed by the additional micro-defects on the strain generated in the rock.Different values of b can reflect the degree of strain softening in the rock，and a method for determining the value of b is given to obtain the value of b under different stress states for different rocks.The deformation of the rock skeletal part without micro-defects and whose damage conforms to the Weibull probability distribution，the deformation of the rock skeletal part and the micro-defective part make up the deformation of the rock，which leads to the derivation of the damage constitutive model of the rock，and the specific determination method of the parameters m and F0 of this damage constitutive model is given.The parameters of the model are discussed.The smaller the value of b，the