《工程科学学报》录用稿，htps:/doi.org/10.13374/i,issn2095-9389.2021.04.11.001©北京科技大学2020 工程科学学报，第37卷，第1期：1-8,2015年1月 Chinese Joumal of Engineering,Vol.37,No.1:1-8,January 2015 D0I:10.13374/.issn2095-9389.2015.01.001;htp://journal..ustb.edu.cn 基于DIC的含3D打印起伏节理试样破裂特性及损伤 本构研究 王本鑫12,3)，金爱兵1,2)，赵怡晴1,2☒，孙浩1,2)，刘加柱12) 1)北京科技大学金属矿山高效开采与安全教有部重点实验室，北京1000832)北京科技大学士木与资源工程学觉Z北京100083 3)Department of Architecture and Civil Engineering.Chalmers University of Technology.Gothenburg4196.Sweden ☒通信作者，E-mail:yiqingz小ao@126.com 摘要受地质构造的影响，岩体工程中经常赋存起伏结构面（如扭转皱褶）， 由于形态复杂，目前起伏节理岩体的破裂及损伤 本构研究仍不充分.采用3D打印技术制作不同倾角的起伏节理模型，通过单轴压缩试验和数字图像相关技术(DIC)对起伏节 理试样的力学及破裂特性进行研究，并基于断裂力学原理，首次提出采用DIC位移场求解节理尖端应力强度因子(SF)进而探 究损伤本构特性的思路。结果表明：通过分析最小强度确定了起伏节理对试样的损伤上限为46.6%，起伏节理试样单轴强度对 倾角的敏感性大于平直节理试样：起裂发生在峰值应力附近，破裂过程可分为破裂路径上微裂隙的产生和同步贯通，破裂模 式表现为多条裂隙张剪组合模式：峰前SF随荷载增加而增加，峰辰同 一荷载下K>K,节理左右两端均匀以剪切形式扩展： 起伏节理对试样的损伤与倾角呈正弦关系，节理和荷载对试样的总损伤与应变均呈“$”型曲线. 关键词起伏节理；DIC,应力强度因子；破裂特性；损伤本构 分类号TU45.1 Study on fracture characteristics and damage constitutive model of 3D printing undulating joint samples based on DIC WANG Ben-xin2)JIN Ai-bihs 2fA0Yh-gimg',2☒，SUN Hao,2）,LIUJia-h,2) 1) Key Laboratory of Ministry ient Mining and Safety of Metal Mines,University of Science and Technology Beijing.Beijing 100083,China School of Civil an ing,University of Science and Technology Beijing,Beijing 100083,China Department of Architecfre and Civil Engineering,Chalmers University of Technology,Gothenburg 41296,Sweden Corresponding author,E-mail:yiqingzhao@126.com ABSTRACT Due to the influence of geological structure,various forms of joint structural planes are often present in rock mass engineering.The undulating structural planes(such as torsional fold surface)are a special kind of them.And they will seriously affect the stablity of the rock mass and bring potential hazards to the rock mass engineering.Because of the complexity of their shape,the research on fracture and damage constitutive law of rock mass with undulating joints is not conducted fully.Undulating joints with different dip angles were fabricated using 3D printing technology.Uniaxial compression test and digital image correlation (DIC) 收稿日期：2021-04-10 基金项目：国家自然科学基金资助项目(52004017)：中国博士后科学基金资助项目(2020M670138):中央高校基本科研业务费 专项资金资助项目(FRF-TP-19-026A1)
工程科学学报，第 37 卷，第 1 期：18，2015 年 1 月 Chinese Journal of Engineering, Vol.37, No. 1: 18, January 2015 DOI: 10.13374/j.issn2095-9389.2015.01.001; http://journal.ustb.edu.cn 基于 DIC 的含 3D 打印起伏节理试样破裂特性及损伤 本构研究1 王本鑫 1,2,3)，金爱兵 1,2)，赵怡晴 1,2)，孙 浩 1,2)，刘加柱 1,2) 1) 北京科技大学金属矿山高效开采与安全教育部重点实验室, 北京 100083 2) 北京科技大学土木与资源工程学院, 北京 100083 3) Department of Architecture and Civil Engineering, Chalmers University of Technology, Gothenburg 41296, Sweden 通信作者, E-mail: yiqingzhao@126.com 摘 要 受地质构造的影响，岩体工程中经常赋存起伏结构面(如扭转皱褶)，由于形态复杂，目前起伏节理岩体的破裂及损伤 本构研究仍不充分. 采用 3D 打印技术制作不同倾角的起伏节理模型，通过单轴压缩试验和数字图像相关技术(DIC)对起伏节 理试样的力学及破裂特性进行研究，并基于断裂力学原理，首次提出采用 DIC 位移场求解节理尖端应力强度因子(SIF)进而探 究损伤本构特性的思路. 结果表明：通过分析最小强度确定了起伏节理对试样的损伤上限为 46.6%，起伏节理试样单轴强度对 倾角的敏感性大于平直节理试样；起裂发生在峰值应力附近，破裂过程可分为破裂路径上微裂隙的产生和同步贯通，破裂模 式表现为多条裂隙张剪组合模式；峰前 SIF 随荷载增加而增加，峰后同一荷载下 KII＞KI，节理左右两端均匀以剪切形式扩展； 起伏节理对试样的损伤与倾角呈正弦关系，节理和荷载对试样的总损伤与应变均呈“S”型曲线. 关键词 起伏节理; DIC; 应力强度因子; 破裂特性; 损伤本构 分类号 TU45.1 Study on fracture characteristics and damage constitutive model of 3D printing undulating joint samples based on DIC WANG Ben-xin1,2) , JIN Ai-bing1,2) , ZHAO Yi-qing1,2), SUN Hao1,2) , LIU Jia-zhu1,2) 1) Key Laboratory of Ministry of Education for Efficient Mining and Safety of Metal Mines, University of Science and Technology Beijing, Beijing 100083, China 2) School of Civil and Resource Engineering, University of Science and Technology Beijing, Beijing 100083, China 3) Department of Architecture and Civil Engineering, Chalmers University of Technology, Gothenburg 41296, Sweden Corresponding author, E-mail: yiqingzhao@126.com ABSTRACT Due to the influence of geological structure, various forms of joint structural planes are often present in rock mass engineering. The undulating structural planes (such as torsional fold surface) are a special kind of them. And they will seriously affect the stablity of the rock mass and bring potential hazards to the rock mass engineering. Because of the complexity of their shape, the research on fracture and damage constitutive law of rock mass with undulating joints is not conducted fully. Undulating joints with different dip angles were fabricated using 3D printing technology. Uniaxial compression test and digital image correlation (DIC) 收稿日期:2021-04-10 基金项目:国家自然科学基金资助项目(52004017)；中国博士后科学基金资助项目(2020M670138)；中央高校基本科研业务费 专项资金资助项目(FRF-TP-19-026A1) 《工程科学学报》录用稿，https://doi.org/10.13374/j.issn2095-9389.2021.04.11.001 ©北京科技大学 2020 录用稿件，非最终出版稿
2 technology were used to study the mechanical and fracture characteristics of undulating joint specimens.And based on the principle of fracture mechanics,an idea using DIC displacement field to solve the stress intensity factor (SIF)at the joint tips and further to study the damage constitutive law is firstly proposed.The results show that:the upper limit of undulating joint damage to specimens is determined with 46.6%through the minimum strength analysis.The sensitivity of undulating joint specimen uniaxial strength to joint dip angle is greater than that of straight joint specimen.The initiation of fracture occurs near the peak stress.The fracture process can be divided into the initiation and synchronous penetration of microcracks on the fracture path.And the fracture mode shows a combination mode of multiple tension and shear fractures.The stress intensify factor(SIF)increases with loading in the pre-peak stage,and the cracks propagate in shear at the joint left and right tips in the post peak stage,because of Ku>Ki under the same stress.The undulating joint damage to specimen with the dip angle is in a sinusoidal curve,and the relationships between the total damage coupled by joint and load with strain are all "S"curves KEY WORDS undulating joint;DIC;stress intensity factor;fracture characteristics;damage constitutive model 岩体工程的失稳破坏受很多因素影响四，除了岩体是否完整以及岩体的连续性，岩体中节理形貌特征 也是影响岩体力学及破裂特性的重要因素)，目前虽然规则节理岩体力学性质已有较多研究，但受限于实际 表征的困难性和复杂性，含起伏等不规则节理岩体力学性质方面的研究仍不够深久.起伏节理的形貌特征与 平直节理有很大差异，室内试验中插缝和切槽等方式不能适用于起伏节理武样的制作.室内条件下如何进行 起伏节理模型的准确表征和精准制作是研究起伏节理岩体力学特性的基础和前提 3D打印技术的兴起4，为解决上述问题提供了有效途径共操作简单灵活，能够制作各种形状的节理 模型，实现同种模型批量制作，保证模型参数的一致性，提高节理角度和位置精度-)目前己有较多学者将 3D打印技术应用于岩体力学及破裂特性的研究2绵涛等基于3D打印技术分别进行了含粗糙节理 网络和粗糙节理试样的单轴压缩和直剪试验，证明了3打技术可用于含粗糙节理试样力学特性研究.金 爱兵等4对含3D打印单一节理和交叉节理试样的破裂机理进行了研究，通过理论分析探究了起裂角与节 理倾角之间的关系，结果与前人研究具有良好的一致催再次证明了3D打印技术可用于探究节理岩体的破 裂特性 荷载条件下岩体破裂和损伤本构特性的研究一直以来都是岩体力学研究的重点.虽然已有较多关于节 理岩体损伤本构方面的研究62)，但是大多以平直节理为前提，起伏节理岩体的本构模型研究涉及较少，需 进一步探究起伏节理岩体的损伤本构模型 应力强度因子(Stress intensit,factor,SIF)2l作为求解节理对岩体损伤的基本参数，可根据断裂力学原理 采用数学方法通过节理尖端的位移场计算求得2].数字图像相关技术(Digital Image Correlation method,DIC) 能够实现对加载过程中试样的金程粒移场进行提取.基于DIC技术将室内试验数据和节理岩体损伤本构模 型相结合，为节理尖端S一起伏节理试样破裂机理和损伤本构分析提供新的途径. 本文采用3D打印技术制作了不同倾角的起伏节理模型，浇筑成试样后进行单轴压缩试验.为了探究起 伏节理试样的破裂机理和损伤本构特性，基于DIC技术对加载过程中试样变形进行监测分析，提取应变及 位移场进行破裂演化规律和破裂模式表征.根据断裂力学方法利用节理尖端位移场计算节理尖端的$F,进 而通过数据拟合对起伏节理试样的损伤本构模型进行求解分析，该方法亦可用于其他节理岩体的破裂和损 伤本构特性的研究 1试验方案 1.1含3D打印起伏节理试样的制作 含3D打印起伏节理试样的制作分为两个部分：起伏节理的制作和试样的浇筑，制作流程如图1所示. 首先采用CAD软件进行起伏节理数值模型的构建，并导出为STL格式，然后将模型导入打印软件并设置打 印参数，最后进行模型的打印.模型采用XYZ Printing3D打印机基于熔融堆积原理制作而成
2 technology were used to study the mechanical and fracture characteristics of undulating joint specimens. And based on the principle of fracture mechanics, an idea using DIC displacement field to solve the stress intensity factor (SIF) at the joint tips and further to study the damage constitutive law is firstly proposed. The results show that: the upper limit of undulating joint damage to specimens is determined with 46.6% through the minimum strength analysis. The sensitivity of undulating joint specimen uniaxial strength to joint dip angle is greater than that of straight joint specimen. The initiation of fracture occurs near the peak stress. The fracture process can be divided into the initiation and synchronous penetration of microcracks on the fracture path. And the fracture mode shows a combination mode of multiple tension and shear fractures. The stress intensify factor (SIF) increases with loading in the pre-peak stage, and the cracks propagate in shear at the joint left and right tips in the post peak stage, because of KII＞KI under the same stress. The undulating joint damage to specimen with the dip angle is in a sinusoidal curve, and the relationships between the total damage coupled by joint and load with strain are all "S" curves. KEY WORDS undulating joint; DIC; stress intensity factor; fracture characteristics; damage constitutive model 岩体工程的失稳破坏受很多因素影响[1]，除了岩体是否完整以及岩体的连续性[2]，岩体中节理形貌特征 也是影响岩体力学及破裂特性的重要因素[3] . 目前虽然规则节理岩体力学性质已有较多研究，但受限于实际 表征的困难性和复杂性，含起伏等不规则节理岩体力学性质方面的研究仍不够深入. 起伏节理的形貌特征与 平直节理有很大差异，室内试验中插缝和切槽等方式不能适用于起伏节理试样的制作. 室内条件下如何进行 起伏节理模型的准确表征和精准制作是研究起伏节理岩体力学特性的基础和前提. 3D 打印技术的兴起[4-6]，为解决上述问题提供了有效途径. 其操作简单灵活，能够制作各种形状的节理 模型，实现同种模型批量制作，保证模型参数的一致性，提高节理角度和位置精度[7-8] . 目前已有较多学者将 3D 打印技术应用于岩体力学及破裂特性的研究[9-12] . 王培涛等[13]基于 3D 打印技术分别进行了含粗糙节理 网络和粗糙节理试样的单轴压缩和直剪试验，证明了 3D 打印技术可用于含粗糙节理试样力学特性研究. 金 爱兵等[14-15]对含 3D 打印单一节理和交叉节理试样的破裂机理进行了研究，通过理论分析探究了起裂角与节 理倾角之间的关系，结果与前人研究具有良好的一致性，再次证明了 3D 打印技术可用于探究节理岩体的破 裂特性. 荷载条件下岩体破裂和损伤本构特性的研究一直以来都是岩体力学研究的重点. 虽然已有较多关于节 理岩体损伤本构方面的研究[16-21]，但是大多以平直节理为前提，起伏节理岩体的本构模型研究涉及较少，需 进一步探究起伏节理岩体的损伤本构模型. 应力强度因子(Stress intensity factor，SIF)[22]作为求解节理对岩体损伤的基本参数，可根据断裂力学原理 采用数学方法通过节理尖端的位移场计算求得[23] . 数字图像相关技术(Digital Image Correlation method，DIC) 能够实现对加载过程中试样的全程位移场进行提取. 基于 DIC 技术将室内试验数据和节理岩体损伤本构模 型相结合，为节理尖端 SIF、起伏节理试样破裂机理和损伤本构分析提供新的途径. 本文采用 3D 打印技术制作了不同倾角的起伏节理模型，浇筑成试样后进行单轴压缩试验. 为了探究起 伏节理试样的破裂机理和损伤本构特性，基于 DIC 技术对加载过程中试样变形进行监测分析，提取应变及 位移场进行破裂演化规律和破裂模式表征. 根据断裂力学方法利用节理尖端位移场计算节理尖端的 SIF，进 而通过数据拟合对起伏节理试样的损伤本构模型进行求解分析. 该方法亦可用于其他节理岩体的破裂和损 伤本构特性的研究. 1 试验方案 1.1 含 3D 打印起伏节理试样的制作 含 3D 打印起伏节理试样的制作分为两个部分：起伏节理的制作和试样的浇筑，制作流程如图 1 所示. 首先采用 CAD 软件进行起伏节理数值模型的构建，并导出为 STL 格式，然后将模型导入打印软件并设置打 印参数，最后进行模型的打印. 模型采用 XYZ Printing 3D 打印机基于熔融堆积原理制作而成. 录用稿件，非最终出版稿
3 Model Model importing and Model printing Undulating specimen Specimen pouring Undulating joint 图1含3D打印起伏节理试样制作流程示意图 Fig.1 Schematic diagram of production process of specimens with 3D printing undulating joints 预制节理为聚乳酸PLA材料刚度远小于试样基质材料，相同应力条件下节理变远大于基质变形，在 剪切作用下能发生断裂，所以3D打印节理变形不会影响试样本有的破裂特性人起伏带理制作时打印了宽度 为50m的底座，底座宽度与试模宽度相等，浇筑时两者可卡牢固定.节理心何参数如图2所示，节理厚 度为lmm,节理长度取两端的连线距离，用1表示，仁15mm:起伏高度为项点到两端连线的距离，用h表 示，文中节理的最大起伏高度hmax=4.5mm,故最大起伏角=arctan(2h/05)50.2°：节理倾角为两端连线 延长线与水平方向的夹角，用a表示，试验设计了0°，30°，45 60°和90°5种有代表性的倾角.试样为长 方体形，长×宽×高=50mm×50mm×100mm,均在相同尺寸的钢制模具中浇筑而成.浇筑材料为水泥砂浆， 按水泥、砂和水以质量比4：2：1混合均匀后配制而成.经测试洗筑材料的单轴强度、弹性模量和泊松比分 别为26.7MPa、7.51GPa和0.25 50mm 起伏节理试件几何参数示意图 diagram of geometric parameters of undulating joint specimen 1.2试验条件及加载方案 为了探究单轴压缩条件下起伏节理试样的力学变形特性、破裂演化规律及损伤本构特性，采用YW600 型微机控制电液伺服岩石压力试验机对试样进行加载，采用Correlated Solutions非接触全场应变监测系统对 试样表面全程变形过程进行监测，基于DIC技术通过Vic-3D图像处理软件对试样表面的全场位移和应变进 行分析，最后根据位移和应变场探究破裂和损伤本构特性，单轴压缩测试系统如图3所示. Load
3 Model construction Model importing and parameter setting Model printing Undulating specimen Specimen pouring Undulating joint 图 1 含 3D 打印起伏节理试样制作流程示意图 Fig.1 Schematic diagram of production process of specimens with 3D printing undulating joints 预制节理为聚乳酸 PLA 材料刚度远小于试样基质材料，相同应力条件下节理变形远大于基质变形，在 剪切作用下能发生断裂，所以 3D 打印节理变形不会影响试样本有的破裂特性. 起伏节理制作时打印了宽度 为 50 mm 的底座，底座宽度与试模宽度相等，浇筑时两者可卡牢固定. 节理几何参数如图 2 所示，节理厚 度为 1mm，节理长度取两端的连线距离，用 l 表示，l=15 mm；起伏高度为顶点到两端连线的距离，用 h 表 示，文中节理的最大起伏高度 hmax=4.5 mm，故最大起伏角 i=arctan(2hmax/0.5l)≈50.2°；节理倾角为两端连线 延长线与水平方向的夹角，用 α 表示，试验设计了 0°，30°，45°，60°和 90° 5 种有代表性的倾角. 试样为长 方体形，长×宽×高=50 mm×50 mm×100 mm，均在相同尺寸的钢制模具中浇筑而成. 浇筑材料为水泥砂浆， 按水泥、砂和水以质量比 4：2：1 混合均匀后配制而成. 经测试浇筑材料的单轴强度、弹性模量和泊松比分 别为 26.7 MPa、7.51 GPa 和 0.25. α 100mm l 50mm h i 图 2 起伏节理试件几何参数示意图 Fig.2 Schematic diagram of geometric parameters of undulating joint specimen 1.2 试验条件及加载方案 为了探究单轴压缩条件下起伏节理试样的力学变形特性、破裂演化规律及损伤本构特性，采用 YAW-600 型微机控制电液伺服岩石压力试验机对试样进行加载，采用 Correlated Solutions 非接触全场应变监测系统对 试样表面全程变形过程进行监测，基于 DIC 技术通过 Vic-3D 图像处理软件对试样表面的全场位移和应变进 行分析，最后根据位移和应变场探究破裂和损伤本构特性，单轴压缩测试系统如图 3 所示. Loading device Speckle specimen Servo control system Camera Light Light Vic-Snap Collect system Tripod Camera beam 录用稿件，非最终出版稿
图3单轴压缩测试系统示意图 Fig.3 Schematic diagram of uniaxial compressive test system YAW-600型试验机最大试验力600kN,加载过程采用位移控制，速率为0.0015mm/s,位移分辨率为3 um,峰后应力降至峰值应力的50%左右时停止加载.Correlated Solutions非接触全场应变监测系统采集的单 张图像大小12MB,全程采用1帧/s的速率进行图像采集，该速率可有效观测到裂隙的演化过程，同时可 以减小采集数据量，便于后期数据处理分析.Vic-3D图像分析软件基于DIC技术可进行最小分辨率50με的 数据处理工作，应变测量范围在0.005%~2000%，比普通处理方法的精度高、测量范围广.为了避免因试样 数量太少造成试验结果误差较大，每种工况条件的试样均进行三次试验，取强度与平均值较为接近的试样 进行研究 2试验结果分析 2.1不同倾角起伏和平直节理试样力学特性对比分析 不同倾角条件下起伏节理试样的应力一应变曲线均有与完整试样相同的压密、弹性、屈服和应变软化四 个变形阶段（图4），但由于起伏节理的存在，应力一应变曲线与完整试样均有所木同，其中0°试样除峰值 强度略有下降外，其压密和弹性变形阶段与完整试样的应力一应变曲线最为接近，由此可知，α=0°时对力学 性质影响最小.由图4中完整试样的应力一应变曲线可知，水泥砂浆林科具有完整的压密、弹性、屈服、应 变软化和残余强度阶段，力学及变形特性与真实岩石十分相似，采用水泥砂浆进行本文研究是合理的。 30r Straight 0 25 20 "Intact 1 10 0.40.60.81.012 Strain/% 图4试样的应力一应变曲线 Fig. Stress strain curves of joint specimens 对比图4前期研究中不同倾角平直节理试样的单轴压缩应力一应变曲线，两种试样在=30°和60°时 的峰值应力较为接近，其他倾角条件下差异均较大.=90时，起伏节理试样的峰值应力接近最小值，但平直 节理试样峰值应力最大(45时，起伏节理试样的峰值应力较大，但平直节理试样的峰值应力最小.这是 由起伏节理复杂性造成的，入节理的起伏使其相对于平直节理具有了更强的各向异性，所以不同倾角条件下 起伏节理试样力学姓质的规律性较差 起伏节理试样强度在14.27~22.07MPa之间波动，相对于完整试样强度减小了17.3%~46.6%.随节理倾 角的变化，起伏节理试样单轴强度呈倒N型，平直节理试样呈“V”型（图5）.起伏节理试样的单轴强度在0° 时最大，而平直节理试样在a=90时最大：a=30°时起伏节理试样的单轴强度最小，且=60°和90时单轴强 度与之相差较小，三者单轴强度与45时对应的平直节理试样最小单轴强度非常接近，差异性小于5%， 相对于完整试样强度最多减小46.6%，故可认为两种节理对试样造成的损伤上限为46.6%.某一强度条件下 (如图5中红色虚线所示)，平直节理试样最多对应两种倾角，起伏节理试样最多对应3种倾角，所以起伏节 理试样强度受倾角影响的敏感性大于平直节理试样，这与起伏节理形态比平直节理复杂，试样的各向异性 更强，强度离散性更大有关，因此会出现相同倾角条件下两种试样的强度差异较大的情况
4 图 3 单轴压缩测试系统示意图 Fig.3 Schematic diagram of uniaxial compressive test system YAW-600 型试验机最大试验力 600 kN，加载过程采用位移控制，速率为 0.0015 mm/s，位移分辨率为 3 μm，峰后应力降至峰值应力的 50%左右时停止加载. Correlated Solutions 非接触全场应变监测系统采集的单 张图像大小 12 MB，全程采用 1 帧/s 的速率进行图像采集，该速率可有效观测到裂隙的演化过程，同时可 以减小采集数据量，便于后期数据处理分析. Vic-3D 图像分析软件基于 DIC 技术可进行最小分辨率 50 με 的 数据处理工作，应变测量范围在 0.005%~2000%，比普通处理方法的精度高、测量范围广. 为了避免因试样 数量太少造成试验结果误差较大，每种工况条件的试样均进行三次试验，取强度与平均值较为接近的试样 进行研究. 2 试验结果分析 2.1 不同倾角起伏和平直节理试样力学特性对比分析 不同倾角条件下起伏节理试样的应力—应变曲线均有与完整试样相同的压密、弹性、屈服和应变软化四 个变形阶段(图 4)，但由于起伏节理的存在，应力—应变曲线与完整试样均有所不同，其中 α=0°试样除峰值 强度略有下降外，其压密和弹性变形阶段与完整试样的应力—应变曲线最为接近，由此可知，α=0°时对力学 性质影响最小. 由图 4 中完整试样的应力—应变曲线可知，水泥砂浆材料具有完整的压密、弹性、屈服、应 变软化和残余强度阶段，力学及变形特性与真实岩石十分相似，采用水泥砂浆进行本文研究是合理的。 0.0 0.2 0.4 0.6 0.8 1.0 1.2 0 5 10 15 20 25 30 Straight 0° Straight 30° Straight 45° Straight 60° Straight 90° Intact Stress/MPa Strain/% Undulating 0° Undulating 30° Undulating 45° Undulating 60° Undulating 90° 图 4 试样的应力—应变曲线 Fig.4 Stress strain curves of joint specimens 对比图 4 前期研究中[14]不同倾角平直节理试样的单轴压缩应力—应变曲线，两种试样在 α=30°和 60°时 的峰值应力较为接近，其他倾角条件下差异均较大. α=90°时，起伏节理试样的峰值应力接近最小值，但平直 节理试样峰值应力最大；α=45°时，起伏节理试样的峰值应力较大，但平直节理试样的峰值应力最小. 这是 由起伏节理复杂性造成的，节理的起伏使其相对于平直节理具有了更强的各向异性，所以不同倾角条件下 起伏节理试样力学性质的规律性较差. 起伏节理试样强度在 14.27~22.07 MPa 之间波动，相对于完整试样强度减小了 17.3%~46.6%. 随节理倾 角的变化，起伏节理试样单轴强度呈倒“N”型，平直节理试样呈“V”型(图5). 起伏节理试样的单轴强度在α=0° 时最大，而平直节理试样在 α=90°时最大；α=30°时起伏节理试样的单轴强度最小，且 α=60°和 90°时单轴强 度与之相差较小，三者单轴强度与 α=45°时对应的平直节理试样最小单轴强度非常接近，差异性小于 5%， 相对于完整试样强度最多减小 46.6%，故可认为两种节理对试样造成的损伤上限为 46.6%. 某一强度条件下 (如图 5 中红色虚线所示)，平直节理试样最多对应两种倾角，起伏节理试样最多对应 3 种倾角，所以起伏节 理试样强度受倾角影响的敏感性大于平直节理试样，这与起伏节理形态比平直节理复杂，试样的各向异性 更强，强度离散性更大有关，因此会出现相同倾角条件下两种试样的强度差异较大的情况. 录用稿件，非最终出版稿
5 Undulating joint specime 0一Strarght joint specimen 12 0 30 60 90 Q/b 图5不同倾角起伏节理和平直节理试样的单轴强度对比 Fig.5 Comparison of uniaxial strength of undulating and straight joints with different dip angles 2.2起伏节理试样破裂模式及破裂演化分析 采用DIC技术对试样的全场应变处理分析，获取压缩过程中起伏节理试样表面的全程位移场，进一步 计算Lagrange应变张量，获取应变场，通过对比位移云图和e2-Lagrange主应变云图，发现e1-Lagrange主 应变云图能有效表征试样的破裂模式，如图6所示，图中的应变集中带与试样被裂后裂隙的分布情况有很 好的对应关系，故采用e1-Lagrange(下文用e1表示)主应变云图对试样的破裂模武进行分析 e%) (9%) % 15.05 326 1.60 006 =30p =60 图6不同倾角起伏理试样的1主应变云图 Fig.6 e principal strain nephogram of undulating joint specimens with different dip angles 由图7()可见，起伏节理试样中的裂隙主要由纵向张裂隙（红线）和倾斜剪裂隙（蓝线）组成，其中，剪裂 隙多分布在节理附近，并向试样两端演桃为沿最大主应力方向的张裂隙.=30°、60°和90°试样，裂隙呈正 “Y形和倒Y”形分布，a0°试样的破裂模式为两条连接节理端部的纵向裂隙，a45试样下部为剪裂隙，上 部为连接节理上端的张裂隙和剪裂隙分支.对比前期研究中不同倾角条件下平直节理试样的破裂模式（图 7(b)4,起伏节理试样的破裂模式更为复杂，平直节理试样的破裂以单条贯通裂隙为主，裂隙类型较为单 一，如=30°和60°试样约切破裂，=90°试样为张拉破裂，而起伏节理试样裂隙数量多大于等于2条，且 多以张剪破裂组合为主 (au=0' =30 5 =60” =90
5 0 30 60 90 12 16 20 24 Uniaxial strength/MPa α/° Undulating joint specimen Straight joint specimen Difference<5% 图 5 不同倾角起伏节理和平直节理试样的单轴强度对比 Fig.5 Comparison of uniaxial strength of undulating and straight joints with different dip angles 2.2 起伏节理试样破裂模式及破裂演化分析 采用 DIC 技术对试样的全场应变处理分析，获取压缩过程中起伏节理试样表面的全程位移场，进一步 计算 Lagrange 应变张量，获取应变场，通过对比位移云图和 e2-Lagrange 主应变云图，发现 e1-Lagrange 主 应变云图能有效表征试样的破裂模式，如图 6 所示，图中的应变集中带与试样破裂后裂隙的分布情况有很 好的对应关系，故采用 e1-Lagrange(下文用 e1表示)主应变云图对试样的破裂模式进行分析. e1(%) 5.05 2.45 -0.15 α=0° e1(%) 3.26 1.60 -0.06 α=30° 3.78 1.80 -0.18 α=45° e1(%) 4.90 2.40 -0.10 α=60° e1(%) 3.18 1.57 -0.04 α=45° e1(%) 图 6 不同倾角起伏节理试样的 e1 主应变云图 Fig.6 e1 principal strain nephogram of undulating joint specimens with different dip angles 由图 7(a)可见，起伏节理试样中的裂隙主要由纵向张裂隙(红线)和倾斜剪裂隙(蓝线)组成，其中，剪裂 隙多分布在节理附近，并向试样两端演化为沿最大主应力方向的张裂隙. α=30°、60°和 90°试样，裂隙呈正 “Y”形和倒“Y”形分布，α=0°试样的破裂模式为两条连接节理端部的纵向裂隙，α=45°试样下部为剪裂隙，上 部为连接节理上端的张裂隙和剪裂隙分支. 对比前期研究中不同倾角条件下平直节理试样的破裂模式(图 7(b)) [14]，起伏节理试样的破裂模式更为复杂，平直节理试样的破裂以单条贯通裂隙为主，裂隙类型较为单 一，如 α=30°和 60°试样为剪切破裂，α=90°试样为张拉破裂，而起伏节理试样裂隙数量多大于等于 2 条，且 多以张剪破裂组合为主. α=0° Tensile crack Shear crack (a) α=30° α=45° α=60° α=90° 录用稿件，非最终出版稿