基于压电陶瓷和DIC的软硬互层岩石内外损伤试验及模型研究.pdf

119DConpiezoelectricceramicsandDICLJJ54(9):119-128.ateinCoalMines移动扫码阅读QI Xianyin,KE Ting,GENG Diandong.Damage test and model study of soft and hard interbedded rock basedJ.煤矿安全,2 0 2 3,5 4(9)：119-12 8.元宪寅，柯婷，耿殿栋.基于压电陶瓷和DIC的软硬互层岩石内外损伤试验及模型研究SafetyinCoalMines2023年9 月Sep.2023煤砺发全第9 期第5 4 卷No.9Vol.54DOI:10.13347/ki.mkaq.2023.09.017基于压电陶瓷和DIC的软硬互层岩石内外损伤试验及模型研究元宪寅12，柯婷1，耿殿栋1（1.长江大学城市建设学院，湖北荆州4 3 4 0 2 3;2.中国科学院武汉岩土力学研究所岩土力学与工程国家重点实验室，湖北武汉4 3 0 0 7 1）摘要：为探究软硬互层岩石在破坏后的损伤演化规律，利用相似模型理论制作3 组类岩石试块(软岩、互层岩石、硬岩)，在对类岩石试块进行单轴压缩的同时辅以压电陶瓷试验和数字图像相关技术(digital imagecorrelation,以下简称DIC)来实时检测岩石的内外损伤，分析其破坏状态和损伤演化规律。试验结果表明：软硬互层岩石的强度、变形特性介于类软岩、类硬岩之间，在达到峰值点前其力学规律类似于软岩有一定的膨胀变形，而在破坏后其力学特性为接近硬岩的脆性破坏；将压电陶瓷时域信号和DIC应变图像与岩石从加载至破坏的过程相结合，发现各类岩石呈现初始平缓、中期加速、末期减缓的“s”形损伤演化规律；压电陶瓷时域信号的衰减幅度与岩石损伤中前期过程正相关，而DIC图像应变变化特征与岩石损伤中后期过程关联密切；软硬互层岩石损伤中前期过程主要是软岩部分的变形及硬岩部分少量的裂纹发育导致的损伤累积变化，而后期则主要是硬岩部分微裂纹相互连接贯通导致的损伤累积，最终表现为软岩部分的粉碎、硬岩部分呈现明显的贯通裂纹。结合压电陶瓷与DIC技术的损伤监测特点，基于尖点突变理论量化岩石损伤曲线中的突变区间和突变交点，以突变交点为岩石内外损伤的分段点，提出内外损伤结合的岩石本构模型。关键词：软硬互层岩石；内外损伤；尖点突变理论；损伤模型；数字图像相关技术中图分类号：TD315文献标志码：A文章编号：10 0 3-4 9 6 X（2 0 2 3）0 9-0 119-10amage test and model study of soft and hard interbedded rock based on piezoelectric ceramics and DQI Xianyin2,KE Ting,GENG Diandong(1.College of Urban Construction,Yangtze University,Jingzhou 434023,China;2.State Key Laboratory of Geomechanics andGeotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences,Wuhan 430071,China)Abstract:In order to explore the damage evolution law of soft and hard interbedded rock after failure,three groups of rock-like testblocks(soft rock,interbedded rock and hard rock)were made by using the similar model theory.In addition to uniaxial compressionof rock-like test blocks,piezoelectric ceramic test and digital image correlation(DIC)technology are used to detect rock internal andexternal damages in real time,and analyze its failure state and damage evolution law.The test results show that the strength and de-formation characteristics of soft and hard interbedded rock are between those of soft and hard rock.Before reaching the peak point,收稿日期：2 0 2 3-0 2-0 3责任编辑：朱蕾基金项目：湖北省自然科学基金资助项目（2 0 2 0 CFB367）作者简介：元宪寅（19 8 6 一），男，山东泰安人，副教授，硕士研究生导师，博士，主要从事油气岩石力学以及数值仿真方面的教学和研究工作。E-mail:120SafetyinCoalMinesSep.20232023年9 月No.9Vol.54煤矿全第9 期第5 4 卷the mechanical law is similar to that of soft rock with certain expansion deformation,but after failure,the mechanical property is sim-ilar to that of hard rock with brittle failure.By combining the time domain signal of piezoelectric ceramics and DIC strain imageswith the process of rock loading to failure,it is found that all kinds of rocks present an“s shaped damage evolution law,which isgentle at the beginning,accelerated at the middle stage and slowed down at the end.The time domain signal attenuation amplitude ofpiezoelectric ceramics is positively correlated with the early stage of rock damage,while the strain variation characteristics of DICimages are closely correlated with the middle and late stage of rock damage.In the early stage,the damage accumulation is mainlycaused by the deformation of the soft rock and the development of a small number of cracks in the hard rock,while in the late stage,the damage accumulation is mainly caused by the interconnecting and transfixating of microcracks in the hard rock.The damage ac-cumulation is finally manifested as the crushing of the soft rock and the obvious transfixating cracks in the hard rock.Combined withthe damage monitoring characteristics of piezoelectric ceramics and DIC technology,the abrupt change interval and abrupt changeintersection in rock damage curves were quantified based on the abrupt change theory,and the abrupt change intersection was takenas the subsection point of rock internal and external damage.A rock constitutive model combining internal and external damage wasproposed.Key words:soft and hard interbedded rock;interior and exterior damage;cusp catastrophe theory;damage model;digital image cor-relation(DIC)不同均质度的岩石在成形过程中，经过挤压、脱水、重结晶等过程被沉积压实后形成软硬互层岩石。与单一岩石不同的是在外力荷载作用下，软硬互层岩石力学特性受内部不同均质度岩石的力学性质和层理接触面的影响，变形和强度呈各向异性，破坏机理和损伤演化模式也较为复杂，在实际生活中由软硬互层岩石引起的工程问题较多，如倾倒变形和滑坡等，因此，如何对软硬互层岩石在破坏过程中的损伤进行监测是十分有必要的。近年来，国内外学者基于室内实验和数值模拟的方法对软硬互层岩石从层厚比、倾角等不同角度的宏观破坏模式到微观损伤模式进行了大量的研究2-5 。LUO 等2 采用DIC 技术发现水平软硬层夹层岩石的破坏是软层材料与硬层材料在压应力作用下耦合破坏的过程；黄锋等3 发现层状岩体的破坏最先从硬层开始再逐渐扩展到软岩的规律；高振恒4 利用压电陶瓷智能骨料传感器接收方法实现了对混凝土材料的内部结构损伤和裂缝识别；杨霞5 基于压电陶瓷的健康检测技术，对GFRP管混凝土新型组合结构各阶段的健康和损伤状态进行了识别，发现识别的结果可为实际工程的实时损伤评估提供参考依据。这些研究中均采用单一损伤检测方法，以此为基础得到的损伤数据与宏观破坏过程并不契合，很难对软硬互层岩石的损伤过程做出准确的量化分析。在对岩石材料进行损伤检测的过程中可以发现，不同监测方法的精度不同，适用的层面和范围也不一致，压电陶瓷方法对于材料内部损伤更为敏感，而DIC方法则对材料外部损伤的表征更为精确。为此，通过开展室内单轴压缩试验并辅以压电陶瓷和DIC同时监测材料内外损伤，基于尖点突变理论6-7 定量分析材料损伤特性，通过结合内外损伤信息来表征岩石完整破坏过程的损伤演化规律，并代人试验数据来验证模型的合理性。1基于压电陶瓷和DIC的单轴压缩损伤试验1.1相似材料试样制备由于从自然界中钻取合适层状复合岩石难度较大，因此众多学者8-10 利用相似理论制作岩石试块用于研究岩石的物理力学性质。为此，采用相似理论制作类软岩、软硬交互类岩石和类硬岩3组类岩石试件，试件流程如图1。I试件制作小网配合比实验试块制作标准试样压电陶瓷系统DIC系统岩样室内试验光源、镜头横菜云台、三脚架岩石内部损伤监测岩石外部损伤监测图1试验流程Fig.1Experimentalflow在岩石试块制作过程中，通过降低硅粉、石英砂等用量，使其具有相对较低的弹性模量和抗121SafetyinCoal Mines2023年9 月Sep.2023煤矿发全No.9Vol.54第9 期第5 4 卷压强度，类软岩；采用水泥、石英砂、硅粉作为相似材料，使试件具有较大的弹性模量和抗压强度，得到类硬岩；再通过采用类软岩和类硬岩的相似材料配比，二者1：1叠加，使其具备软硬交互岩石的物理力学性质以及软硬胶结面，即得到软硬交互岩石。但需要注意的是由于胶结面的存在，软硬交互岩石的力学特性并不是单一岩石的简单叠加，在相似岩样制备过程中不仅要满足几何尺寸的相似，还要综合考虑其边界下应力条件、物理力学特性的相似，因此依据相似理论原理1 可知，相似岩样与原岩间的关系为：C。Cf=C,=C=Co=1Cc=C。=Co t =Cs c =CE(1)式中:CL、Co、Cp、Cr、Ce、Cu、Co、C。、Co t、Cc、Ce 分别为几何、应力、密度、内摩擦角、应变、泊松比、摩擦系数、黏聚力、抗拉强度、抗压强度、弹性模量相似比尺，其中：CL=3.2,C,=1.39,C=4.448。由于软硬互层岩石的物理力学特性会受到胶结面以及2 类岩石的耦合影响，根据LEMAITRE等12 的研究，在考虑胶结面的情况下软硬交互岩石的主要物理参数弹性模量E可表示为：E,E2(LI+L2)Ep=(2)(E,L2+E2L)式中：LI、L 为软硬交互岩石中类硬、类软岩的高度，mm；Ei、Ez 为软硬交互岩石中的类硬、类软岩的弹性模量，MPa。因此，为使类岩石试件拥有与原岩相似的结构和破坏特征，其主要力学参数如弹性模量应符合上述规律。在多次配比试验后，得到的类软岩配比为：水泥：石英砂：硅粉：水：消泡剂：减水剂=1：0.8：0.1：0.3：0.0 0 3：0.0 0 3；类硬岩配比为：水泥：石英砂：水：消泡剂=1：0.7：0.4：0.0 0 3，在此配合比下相应力学参数符合式（1)、式(2)，相似材料配比方案及力学性能参数见表1。1.2试验过程本次单轴压缩试验采用压电陶瓷与DIC系统结合HYAS-1000C岩石三轴仪进行损伤全过程监测。试验过程中为确保试验过程中信息采集的统一性，在正式开始时保证压力机、DIC端和压电表1相似材料配比方案及力学性能参数Table 1 Proportioning scheme and mechanical propertyparameters of similar materials岩石类型单轴抗压强度/MPa弹性模量/GPa砂岩178.8626.21类硬岩46.586.82泥岩143.4618.86类软岩37.364.91软硬交互类岩石42.325.84陶瓷端的同步记录。其中压电陶瓷无损检测方法对软硬互层相似岩石在受压过程中的内部损伤进行实时检测，压电陶瓷检测系统包括电脑端、压电传感器、信号发射器和数据采集器(NIUSB-6361)。将1对智能骨料放在试块左右对称两侧，分别作为激励端和接收端，在骨料和试块的表面均匀涂抹硅脂，起到黏结和填充二者接触面保持能量稳定接收的作用，来实现单轴压缩过程中电信号的正常传送和接收。而DIC的检测系统对软硬互层相似岩石在受压过程中的外部损伤进行实时检测，包括2 台高速相机，照明设备、标定板和XT-DIC三维应变分析系统等，其中相机的最大采集频率为12 5 0 0 0 fps，拍摄间隔为1s1张，压缩过程中对岩样表面进行图像实时采集，最后再用DIC分析系统对图像进行处理，得到表面位移、应变云图信息和数据。1.3压电信号图将试验测得的压电信号经MATLAB软件做去噪处理后得到时域信号图，3 种不同岩石时域信号图如图2。每条曲线的周期均为1s，试块的损伤破坏过程分为4 个阶段，分别为初始压密阶段、弹性变形阶段、塑性屈服阶段、峰后破坏阶段，通过分析不同加载阶段压电信号值的变化幅度和衰减趋势来判断试块内部的损伤特征。由图2 可看出：3 组类岩石类软岩、软硬互层类岩石、类硬岩的压电陶瓷峰值电压分别为0.0 14 8、0.0171、0.0 2 5 3 V，出现在岩石初始压密阶段，此电压反映了岩石的强度特征，即峰值电压越大，岩石强度越高13 ；随着轴向荷载增大，岩石进入弹塑性破坏阶段，与之对应的，压电陶瓷电压信号衰减，这是由于轴向荷载增大，岩石内部产生新生裂隙，裂隙阻断应力波传输，造成压电陶瓷电压信号衰减，因此电压信号的衰减幅度可反映岩石内部的损伤情况。122Safety in Coal MinesSep.20232023年9 月No.9煤防发全Vol.54第9 期第5 4 卷800.020.02软岩-压密阶段广软岩-弹性阶段审田审00700.020.020.20.40.60.8 1.000.20.40.6 0.81.0时间/s60时间/s0.020.02软岩-塑性阶段【软岩-蜂后阶段00500.020.0200.20.40.60.81.000.20.40.60.81.0时间/s时间/s4030塑性阶峰后阶20104压密阶段弹性阶段段段00.2 0.40.60.81.01.21.4 1.6应变/%(a)D-0-11000.02软硬交互岩石0.02软硬交互岩石压密阶段弹性阶段09000.020.0200.20.40.60.81.00.20.40.60.88001.0时间/s时间/s0.02厂软硬交互岩石0.02软硬交互岩石70塑性阶段峰后阶段00600.020.0200.20.40.60.81.000.20.40.60.81.0时间/s时间/s504030塑性阶破坏阶段2010压密阶段弹性阶段段00.20.40.60.81.01.21.41.6应变/%(b)D-50-11000.030.03硬岩-压密阶段硬岩-弹性阶段90000.030.20.03800.40.60.81.00.20.40.60.8001.0时间/s时间/s0.03硬岩-塑性阶段0.03硬岩-峰后阶段70宝审00600.030.20.40.80.030.61.000.20.40.60.81.0时间/s时间/s504030塑性阶破坏阶段2010压密阶段弹性阶段段00.20.40.60.81.01.2应变/%(c)D-100-1图2 3 种不同岩石时域信号图Fig.2Time domain signal diagrams ofthreedifferentrocks从压密阶段至塑性阶段的衰减幅值来看，类软岩、软硬互层类岩石、类硬岩的衰减幅值分别为0.0 113、0.0 13 6、0.2 18 V，结果表明：类硬岩信号衰减幅度最大，其压缩过程中裂隙表征最明显，其破坏为典型的脆性破坏特征；而类软岩则表现出一定的膨胀破坏特征，信号衰减幅值最小；而软硬交互类岩石则介于两者之间但更接近于软岩，该现象表明软岩交互岩石在峰前阶段的力学特性更接近于类软岩。从塑性阶段至峰后破坏阶段的衰减幅值来看软岩、软硬互层类岩石、类硬岩的衰减幅值为0.00195、0.0 0 2 17、0.0 0 2 9 9 V，很明显：峰后破坏阶段压电陶瓷的信号变化情况几乎可以忽略不计，即压电陶瓷对于岩石峰后破坏的损伤预估效果并不理想，单纯依靠压电陶瓷的信号衰减评估岩石损伤过程存在相当的难度14 1.4DIC表面应变云图DIC全场应变云图可直观的分析试样局部裂纹的萌生、扩展及相互贯通的过程，3 种不同岩石主应变场演化如图3。由图3(a)可以看到：在初始压密阶段，类软岩内部受轴向压力压缩后，发生明显压缩变形，局部应变分布极不均匀，这是由于软岩内部颗粒孔隙大，颗粒分布不均匀导致的；在弹性阶段，局部应变在上一阶段基础上，局部变形进一步增大，不均匀性增加，但是并没有明显的可表征明显裂隙产生的大变形产生；在塑性阶段，标志性的局部大应变产生，此时最大局部应变达到1.19%，岩石表面出现了明显拉伸裂隙；在峰后破坏阶段，类软岩的表面拉伸裂隙贯通，此时岩石发生明显拉伸破坏。软硬交互类岩石与类硬岩的破坏过程如图3(b)和图3(c)，与类软岩试样不同的是：软硬交互类岩石和类硬岩在初始压缩阶段，局部应变分布相对均匀，这是二者内部孔隙空间相对较小，颗粒硬度大导致的15 ；在峰后破坏阶段，软硬交互类岩石与类硬岩在局部应变为1.4 9%和1.05%时，二者均产生了明显的单裂隙剪切破坏，由此可以看出软硬交互类岩石在破坏特点上更接近于类硬岩，但与之不同的是软硬互层岩石破坏过程中所产生的裂隙虽也表现出脆性破坏特征，但其裂纹的数量及开度都较类硬岩更少，且岩石的破坏是以软岩部分的膨胀开始，逐渐延伸至硬岩部分随后出现一定的脆性破坏。DIC应变云图得到的应变数据属于局部应变特征，而单轴实验中试验机得到的应变曲线为宏观应变特征，因此DIC应变云图无法直接与试验机应变数据比对验证，只能通过间接方法转化成宏观应变数据才能进一步的分析岩石的损伤演化过程。在试验进行的中前期即岩石的初始压缩、弹性阶段，DIC得到的局部应变数据离散性较大，局部应变分布不均匀，很难通过定量方法得到与123SatetyinCoalMines2023年9 月Sep.2023煤码发全第9 期第5 4 卷No.9Vol.541000.049%0.119%1.484%2.058%9080700.0189%0.034%0.702%0.980%6050-0.0125-0.099%403020塑性破坏10压密阶段弹性阶段阶段 阶段00.20.40.60.81.01.21.41.6应变/%(a)D-0-11201.627%1:357%2.343%2.815%1000.773%0.626%1.083%1.191%80600.0800.1050.178-0.432%40塑性阶段破坏阶段20压密阶段弹性阶段00.20.40.60.81.01.21.41.6应变/%(b)D-50-11200.071%0.210%1.237%5.221%1000.033%0.088%0.585%2.536%80600.005.034-0.0679-0.148%40塑性阶段破坏阶段20压密阶段弹性阶段00.20.40.60.81.01.21.4应变/%(c)D-100-1图3 3 种不同岩石主应变场演化Fig.3Evolution of principal strain field inthree different rocks宏观应变相匹配的变化趋势；而在中后期，在岩石表面出现明显裂隙后，局部应变分布极具规律性，相应的定量分析得到的宏观应变与试验机得到的宏观应变数据一致性较高16 。因此，DIC试验方法适合观察岩石的压密与弹性阶段的局部变形与塑性破坏阶段的宏观变形特征。2基于内外损伤结合的软硬互层岩石本构模型在结合压电陶瓷和DIC方法的单轴试验过程中，压电陶瓷在试验的中前期(初始压密、弹性阶段)表征合理，而DIC方法在实验的中后期（塑性阶段、峰后破坏阶段)表征更为有效，因此有必要结合二者的优点，针对岩石的内外损伤，提出1种合理的软硬互层岩石损伤本构演化模型。2.1基于压电陶瓷试验的内损伤表征方法目前基于压电陶瓷的损伤状态分析主要采用小波包能量法17 。高振恒等4 提出了基于小波包能量的损伤指标D，将健康状态下接收的能量Eh作为基准信号，E,为某一加载时刻对应的损伤状态下接收的能量，此时材料的损伤指数可表示为：Eh-E;D;=(3)Eh式中：E为无损状态下接收到的最大幅值能量；E,为在单轴压缩损伤状态下以15 s为单位接收到的累计信号值。D为O时，材料为健康无损伤状态，当D,达到1时，表明材料功能失效处于完全损伤状态。由前文分析可知，随着单轴压缩过程中应力的增长，损伤变大，E,会对应的发生衰减，同时损耗能量E会增大，其中，E,=E-E，那么损伤指数D,可表示为：ED:=(4)Eh上述材料的损伤过程是由于材料内部单元微裂隙的萌生、扩展和贯通导致的，因此可采用统计损伤力学理论18 ，假设微观岩石单元的强度均符合Weibull分布，其概率密度函数为：p(-(%）)m1m-1P(a)=8exp(5)式中：P()为岩石单元的概率密度函数；为岩石应变；m、F为表征材料物理力学性质的参数。若将单轴压缩下岩石的内部损伤变量Dint定义为由压电陶瓷测得的单位时间里损耗能量E与无损状态下健康能量E的比值，可得到：EDint=E(6)E压电陶瓷所测得单位时间内损耗能量即单位面积内岩石单元微裂隙发育所对应的应变能E,为：CFEi=EnP(e)de(7)J0124Safety in Coal MinesSep.20232023年9 月煤防发全No.9Vo1.54第9 期第5 4 卷将式（5）、式(7)代人式(6），可得岩石的内部损伤变量Dint为：intD int=1-exp(8)int式中：m_int、F_i n t 值为基于压电陶瓷试验用以表征材料物理力学性质的Weibull分布参数，反映的是岩石材料在外界荷载下其内损伤的响应特征。通过式(6)计算得到3 组试块的损伤实测值，代人式(8)拟合分别得到3 组不同的m_int、F_i n t值及对应的损伤拟合曲线，3 组岩石的内损伤-应变曲线如图4。损伤变量值可以对岩石的损伤严重程度做定量分析，3 组类岩石类软岩、软硬互层类岩石、类硬岩在应变为0.0 0 4 之前损伤曲线增长较缓，此时对应单轴压缩过程中岩石的压密阶段，试块内部有细小微裂纹萌生，内部损伤不明显，无应力集中区；应变在0.0 0 4 之后损伤变量呈现明显增长趋势，此时出现应变集中区，裂纹稳定扩展；其中软岩和硬岩在应变为0.0 0 6 和0.0 0 8 时损伤加速发展，可看作试件发生失稳破坏的临界点，随着单轴压力的增大，岩样内部损伤加速，裂纹加速扩展形成宏观裂纹；软硬交互岩石在损伤变量达到0.812后，曲线增长趋势减缓，对应峰后残余应力发展阶段，试件逐渐失去承载力，当损伤接近1时，表明试块受到严重破坏，失去基本功效。由此可见：3 组试块在破坏的过程中，应变的增加、试块表面裂纹的扩展都是损伤在积累的体现19 ,随着损伤程度的增加，损伤变量也随之增长；在软硬互层岩石压缩中前期，损伤曲线拟合较好，而在中后期拟合效果较差，预测损伤曲线与压电陶瓷实测值偏差较大。2.2基于DIC的岩石外损伤表征方法DIC方法所测得的应变为局部应变数据，无法直接用于材料整体损伤评估。在张皓2 0 的研究基础上可发现，DIC方法测得的观测应变场内所有点的前5%10%的较大应变点的平均值可反映微裂隙的发育情况，而观测区所有应变的平均值接近岩石试验机直接测得宏观应变值，二者之差定义为应变偏差。M一之N18=(x);(9)M=1=11.0D-0-3实际损伤变量-D-0-3拟合损伤变量0.8D=1exp(-(e/1.024 48).62 5)R2=0.996 910.60.40.200.20.40.60.81.01.21.41.61.8应变/%(a)D-0-11.0D-50-3实际损伤变量D-50-3拟合损伤变量0.8D=1-exp(-(e/0.988 4)2.685 4)R2=0.982 60.60.40.200.20.40.60.8 1.01.21.41.6应变/%(b)D-50-11.0D-100-1实际损伤变量-D-100-1拟合损伤变量0.8D=1exp(-(e/0.953 1)2.57 8)R2=0.996770.60.40.200.20.40.60.81.01.21.41.6应变/%(c)D-100-1图4 3 组岩石的内损伤-应变曲线Fig.4Internal damage-strain curves of three groups ofrocksN式中：之(x),为同一水平应变场内所有点的N=1M1应变平均值；(sx);为所有点中前10%的较大Mi=1应变点的平均值；M=10%N。相比于直接采用较大应变平均值，应力偏差可更好地反映岩石内部的损伤程度，定义外部损伤变量Dext为：125Safety in Coal MinesSep.20232023年9 月No.9煤防发全Vo1.54第9 期第5 4 卷Dext(10)=式中：max为的最大值即岩石临界破坏状态的值。同样采取Weibull统计分布理论，DIC所测得应变偏差即单位面积内岩石所有可观测微裂隙发育变形程度，即：F=EmaxP(e)ds(11)将式(11)代入式(10)可得外损伤变量为：m_extDext-exp(12)8max_ext式中：m_ext、F_e x t 值为基于DIC试验用以表征材料物理力学性质的Weibull分布参数，反映的是岩石材料在外界荷载下其外损伤的响应特征。3组类岩石类软岩、软硬互层类岩石、类硬岩岩石试块的外损伤变量随应变的变化趋势如图5。3组类岩石试块的外损伤变量的趋势基本一致，在初始压密阶段裂纹基本不发育，损伤D值较小趋近于0；当应变分别达到0.7 9%、0.9 1%、0.6%时出现损伤加速点，损伤D呈直线式剧烈增长，此时在外部荷载的作用下，试块表面裂纹加速扩展。其中软硬互层类岩石损伤加速点出现的时间相比软岩和硬岩较晚，是因为软硬互层类岩石前期先在强度较小的软岩处开始发生破坏，而强度较大的硬岩对软岩的竖向破坏变形具有约束作用，随着外部压力的不断增大，软岩处的裂纹逐渐延伸扩展至硬岩处，实现软硬岩裂纹的贯通2 1。可以看到，以DIC方法为基础的损伤预测曲线对于软硬互层岩石破坏的中前期过程拟合效果较差。2.3基于尖点突变理论的内外损伤结合模型验证结合压电陶瓷与DIC测得的内外损伤演化趋势，可发现内损伤在前期荷载作用下损伤值是处于缓慢增长的趋势；而由于本研究用的是圆柱体试块，内部损伤只有积累到一定程度才会有试块表面微裂纹的迅速扩展，所以外损伤在前期没有明显的变化，损伤值接近零，在某一应变点之后才剧烈增长，因此，采用压电陶瓷方法的内损伤评估方法在前期更符合材料的损伤损伤趋势，而后期损伤主要体现在岩石破坏的外表面，采用DIC损伤评估方法来判断会更准确，所以若选择在某一应变点之前用内损伤模型表征，之后用外损伤表征可更直接、全面的判断岩石破坏的整个损伤演化过程。1.0D-0-2实际损伤变量-D-0-2拟合损伤变量0.8D-1-exp(-(c/0.931 8)11.615)0.6R2=0.99250.40.200.20.40.60.81.01.21.41.61.8应变/%(a)D-0-11.0D-50-2实际损伤变量0.8-D-50-2拟合损伤变量D=1-exp(-(e/1.094 5)15.50 8)R2=0.971520.60.40.200.20.40.60.81.01.21.4 1.6应变/%(b)D-50-11.0D-100-1实际损伤变量-D-100-1拟合损伤变量0.8D=1exp(-(e/0.912 0)5.952)R2=0.996.90.60.40.200.20.40.60.81.01.21.41.6应变/%(c)D-100-1图5 3 组岩石的外损伤-应变曲线Fig.5External damage-strain curves ofthree groups of rocks岩石的损伤过程是1个从稳定破坏到突变破坏的过程，当轴向荷载加载到一定水平时，损伤变量会发生明显的突变，破坏的过程可以用尖点突变模型来分析论证，即内外损伤的分界点可采用尖点突变理论来确定。由前文可知，损伤与应变之间的变化趋势为1条非线性曲线，在数学上可采用一元高次方程形式来表示：126Safetyin Coal MinesSep.20232023年9 月No.9煤砺发全Vo1.54第9 期第5 4 卷D(s)=ao+aiel+a2?+ag+a4e4(13)式中：,为与应变相关的待定参数；i可取值1、2、3、4。根据Tschirnhaus变换将式(13)转换成尖点突变模型的标准形式，令=x-A，A=a s/4 a 4，则式(13)可转化为：D(x)=x4+ux?+vx+c(14)3a3a22a1a2a3a33U=(15)8a42Va4a42a428a43对式(14)求导满足D(x)=0，D (x)=0，联立这2 个方程，可得控制阈值：=8 u 3 +2 7 v 2(16)其中当 0 时，表示岩石损伤较小或未发生明显损伤，=0 时，岩石处于损伤变化临界值状态，8时，损伤变量以DIC方法拟合的损伤参数进行预测3组岩石损伤突变曲线如图6。从图6 可看出：随着应变的增大试块的损伤也逐渐增长，结合应力应变曲线，3 类岩石的突变区间都发生在应力应变曲线中的塑性屈服阶段内，突变初始点分别对应峰值强度的5 7.2 6%、5 8.10%和6 5.2 0%，可发现3 种试块的突变区间逐渐向峰值应力靠近；其中硬岩发生突变空间的位置最接近于峰值强度，这是因为硬岩强度最大，加载过程属于脆性破坏，从初始压密阶段到损伤突变持续时间较长，而软岩最早发生出现突变区间也是由于其抗压强度较小；软硬交互类岩石的突变区间出现在应变为1.0 1%1.15%之间，较于其他2软岩压电陶瓷损伤变量1.030软岩DIC损伤变量一软岩应力应变曲线0.8突变交点 0200.60.4100.2000.20.40.60.81.0 1.21.41.6应变/%（a）软岩损伤突变曲线软硬交互岩石压电陶瓷损伤变量一软硬交互岩石DIC损伤变量501.0二软硬交互岩石应力应变曲线400.8突变交点 0300.6突变区0.4200.210间000.20.4 0.60.81.01.2 1.41.6应变/%（b）软硬交互岩石损伤突变曲线硬岩压电陶瓷损伤变量501.0硬岩DIC损伤变量一硬岩应力应变曲线400.8突变交点40300.620变区0.210区间000.20.4 0.60.81.01.21.41.6应变/%（c）硬岩损伤突变曲线图6 3 组岩石损伤突变曲线Fig.6Calculation results of sudden change characteristicvalues of hard rock damage variables组发生突变的位置稍微靠后，是因为软硬交互岩石的破坏属于软岩与硬岩之间的相互作用，硬岩对软岩的破坏和变形具有抑制作用2 2 ，所以软硬交互类岩石从裂纹的萌生、扩展、贯通到发生损伤突变的时间较长。本研究基于尖点突变理论找到内外损伤变量的突变交点，以该点为分段点建立新的内外损伤结合的分段模型，该点之前用内损伤模型表征，该点之后用外损伤模型表征来构建新的分段损伤曲线，3 组岩石内外损伤结合模型如图7。127SafetyinCoal Mines2023年9 月Sep.2023煤防发全Vol.54第5 4 卷No.9第9 期35软岩分段损伤模型1.0软岩应力应变曲线300.8250.620150.4突变交点100.25000.20.40.60.81.01.21.41.61.8应变/%（a）软岩损伤突变曲线软硬交互岩石分段损伤模型501.0厂软硬交互岩石应力应变曲线突变交点400.8300.60.420档0.210000.20.40.60.81.01.21.41.6应变/%（b）软硬交互岩石损伤突变曲线50硬岩损伤分段模型1.0硬岩应力应变曲线400.8300.60.4突变交点0.210000.20.4 0.60.81.01.21.41.6应变/%（c）硬岩损伤突变曲线图7 3 组岩石内外损伤结合模型Fig.7Calculation results of sudden change characteristicvalues of hard rock damage variables由新的损伤曲线可清晰看出：软硬互层岩石损伤发育过程中，前中期的线型与类软岩近似，这是因为在该阶段岩石损伤的产生主要体现在软岩部分的变形及硬岩部分少量的裂纹发育；而后期曲线的变化趋势近似于类硬岩，这是因为软岩部分在较低应力时就已基本破坏完全，当达到峰值点后，岩石破坏已主要为硬岩部分微裂纹相互连接贯通导致的损伤累积。因此软岩互层岩石在最终破坏后表现为软岩部分的粉碎、硬岩部分呈现明显的贯通裂纹。采用本研究的改进损伤模型后，相比于单一的压电陶瓷或者DIC损伤评估方法，在单向受压的岩石破坏过程中，在压密与弹性阶段，岩石内部初始孔隙或裂隙所造成的损伤若采用单一DIC评估方法，其中前期损伤值过低，无法反映岩石内部既有缺陷