1. 引言
温度是影响岩石物理力学性能的重要因素,诸多岩石工程都涉及高温岩石问题 [1]。例如矿山、化石能源和地下热能深部开采、核废料深部地质处置、城市地下空间开发 [2] 等。因此,研究高温作用后岩石的力学性质在工程应用上具有重大意义。
国内外很多学者对高温后和高温下的岩石力学性质开展了试验研究 [3]。张志峰 [4] 利用压汞法测试了25℃~1200℃高温热处理后花岗岩样品的孔隙特征,并研究了不同温度影响下岩石孔隙的分形结构和孔隙率演化模型;黄真萍 [5] 和韩观胜 [6] 分别以石灰岩和砂岩为研究对象,研究了不同高温遇水冷却岩样的力学性质;Yuliang Zhang [7] 研究了高温下砂岩孔隙特性和力学性能的变化,发现高温会增加砂岩累计孔隙体积和孔隙率,降低力学性质;梁永庆 [8] 研究了5种温度作用后砂岩的纵波速度以及力学特性,结果表明高温作用使砂岩纵波速度减小;杨礼宁 [9] 结合扫描电镜技术对高温作用后砂岩试样的物理力学性质和微观结构的变化规律进行观察,结果表明了砂岩试样的变形位移总体上随处理温度升高而增大,加热温度越高时砂岩微观孔隙越多;Weiqiang Zhang [10] 研究了高温作用后砂岩和花岗岩的矿物成分、内部结构和含水量的变化,以及其导致的物理力学性质的变化;杨圣奇 [11] 通过巴西劈裂及常规三轴压缩试验对高温后两种晶粒花岗岩的破坏力学特征进行了研究,研究表明花岗岩高度随着温度升高不断增高,而密度和拉伸强度不断降低;Chengdong Su [12] 制备了不同加热温度的煤样,采用超声波、扫描电镜、核磁共振、单轴压缩试验等多种检测方法,深入研究了高温对煤结构和力学性质的影响;高春杨 [13] 测量了高温下不含水的地幔岩电导率,结果显示岩石的电导率随温度升高而大幅度增大;Yan-jun Shen [14] 通过系统实验探讨分析了不同温度梯度的冷却冲击对高温花岗岩开裂的影响;田振兴 [15] 对花岗岩进行了高温–淬火循环处理后再进行单轴压缩试验和巴西劈裂试验,分析其抗压强度、弹性模量和抗拉强度变化规律;Shuting Miao [16] 采用声发射(AE)、数字图像相关(DIC)和光学显微镜观察相结合的方法,对25℃~800℃条件北山花岗岩物理力学特性和压裂性能进行了实验研究,实验结果表明低温下的花岗岩破坏以拉伸微裂为主,而高温下的花岗岩逐渐以剪切破裂为主;熊健 [17] 对高温影响下页岩岩石的声学特性进行了试验研究,研究了高温处理后页岩岩样物理性质和动弹性参数的变化规律,并讨论了温度对页岩岩石声波特性的影响;Yan [18] 对高温处理后具有不同裂纹倾斜度的3个有缺陷砂岩试样进行了单轴压缩试验,并利用摄影、声发射和数字图像相关技术研究了裂纹的萌生、扩展和合并行为;Wang [19] 通过高温冷却后,将花岗岩浸泡在不同养护温度的水中以及施加冲击荷载的方法模拟井筒高温围岩条件,并对圆形花岗岩试件在径向冲击载荷作用下、热处理和水养护条件下的力学特性进行了试验研究;Tu-bing YIN [20] 开展了花岗岩在高温和轴压作用下的冲击加载试验,并用扫描电镜观察了不同温度下花岗岩的内部结构特征,结果表明轴向压力下的热损伤大于纵波速度所描述的热损伤,随着温度升高,岩石内部结构出现原生裂隙和新生裂隙的扩展和缓慢增加;Xiaohu Zhang [21] 在扫描几何数据的基础上建立三维模型,评估岩石温度、水流速度、粗糙度和裂缝面孔径大小对换热系数的影响,模拟试验表明流速对岩石粗糙度影响最大,其次是孔径大小和岩石粗糙度。
上述研究对在高温作用下岩石的静力学特征研究较多,而应力和温度耦合作用下岩石的动力学机理相对比较复杂,研究较少。因此,为了揭示对岩石在力–热耦合作用下的动力学特征,开展相关研究是十分有必要的。本文利用f 50 mm分离式霍普金森压杆(Split Hopkinson pressure bar, SHPB)实验装置和配套高温环境箱,对常温(25℃),100℃,300℃,500℃,700℃共5种温度环境中黄砂岩试样,通过选用0.2 MPa,0.3 MPa,0.4 MPa,0.5 MPa,0.6 MPa等五种不同冲击载荷进行冲击动力学试验,进一步开展了对含交叉节理黄砂岩动力学分析,对黄砂岩应力峰值、破坏形态等关键力学参数进行系统地研究,探讨了温度和冲击载荷耦合作用下变化的规律。
2. 试验方法
2.1. 试样制作
所选岩石试样为采自江苏省徐州市矿区的含交叉节理黄砂岩,呈暗黄色,主要成分有石英、长石、云母、方解石等。根据国际岩石动力学委员会的建议 [22],所选用于冲击压缩试验的圆柱形试样直径为50 mm左右,厚度为25 mm。对试样两端和侧面进行打磨(如图1),使得两端面的平行度不超过0.05 mm [23]。为了试验的准确性,避免因岩石试样不同而引起不必要的误差,岩石试样均取自同地同一岩块。
2.2. 实验装置
实验采用实时SHPB高温实验系统。该系统由直径50 mm SHPB实验装置和配套高温环境箱组成,如下图2所示。
SHPB实验装置由气室、入射杆、透射杆、吸收杆、发射管道、子弹、阻尼器等组成。其中,入射杆和透射杆均为MnCr合金高强度钢制成,其弹性模量为206 GPa,屈服强度为1200 MPa。入射杆与透射杆侧面对称地贴有两个相同的应变片,以实时记录两杆上的变形状况。
2.3. 试验过程
将岩石试样进行编号处理,分为5组,每个组5个试样。SHPB试验时,检查调整入射杆、透射杆和吸收杆是否处于同一水平线上,调试LK2400N数据采集装置和应变片;并调整阻尼器,使其达到最大缓冲效果。
(a) (b)
Figure 1. Rock sample. (a) Both ends of the joint; (b) Lateral joint
图1. 岩石试样。(a) 两端节理;(b) 侧面节理
Figure 2. SHPB experimental setup schematic diagram
图2. SHPB实验装置原理图
试验操作时,将岩石试样放入配套高温环境箱中,利用温控器设置指定温度后,加温使其达到指定温度,再恒温保存30 min,使试样维持该温度下的状态。在夹岩石的入射杆与透射杆端均匀擦抹凡士林,以增大摩擦效应。在子弹撞击入射杆侧贴上圆柱形波形整形器,从而使试样更早地达到均匀应力状态。当试样加热到预定温度时,启动高温同步组装系统,推动透射杆向试样移动,同时开启空气炮发射子弹,实施高温黄砂岩SHPB冲击试验。子弹撞击入射杆后,在入射杆中产生应力脉冲;通过有效控制,使得应力波到达入射杆与试样接触面时,入射杆、试样和透射杆刚好紧密接触,从而避免对杆传热引起的测量误差。
2.4. 二波法处理
二波法 [24] 是基于一维假定及均匀假定的测试数据处理方法,被广泛应用于各种材料测试的结果处理。根据图2的SHPB测试装置,按照一维应力波理论,测试过程中试样的平均应变率
和平均应变
,试件两端的面应力的平均值
可由下式得出。
(1)
若将试件引入均匀性假定,则上式可变为经典的二波法测试数据处理公式,如下。
(2)
上式中,C、E、A分别为压杆的弹性波波速、弹性模量及截面面积,
为试件的长度和截面面积,
分别为测得的入射波、反射波和透射波。
3. 结果及分析
3.1. 动态应力平衡关系
图3为加载下试样的典型动态应力平衡图。由图3可以看出,在临近入射杆的一端,试样端面应力为入射应力和反射应力叠加,而临近透射杆的一端是由透射应力引起的。两端动应力几乎相等,表明试样在动态加载过程中,基本上处于应力平衡 [25] 状态,满足均质性假说,保证了实验可靠性。
Figure 3. Test chart for stress state uniformity of sample
图3. 试样应力状态均匀性检验图
3.2. 不同温度下黄砂岩的应力应变曲线
常温(25℃)~700℃黄砂岩试件采用不同冲击载荷进行SHPB冲击压缩试验,得到动态应力–应变曲线,如图4所示。
如下图所示,常温(25℃)和高温(100℃~700℃)条件下,黄砂岩在不同冲击载荷作用下,动态应力应变曲线可分为4个阶段。第一阶段:初始压密阶段。由于外力作用,黄砂岩原有微裂缝产生闭合。在此过程中,曲线斜率较小,很小的应力可以引起较多应变;试样被压实,将原有空隙去除;曲线斜率逐渐增大。第二阶段:弹性变形阶段。此阶段中,曲线可看做一条直线,其斜率保持不变,即黄砂岩的动态弹性模量。第三阶段:塑性变形阶段。此阶段中,曲线的斜率发生变化,呈减小走势,使得曲线走向呈现上凹状态。不同温度,不同冲击载荷下,曲线斜率变化速度不一样。当曲线斜率到达0时,到达峰值应力。第四阶段:破坏阶段。此阶段中,曲线快速下降,曲线斜率变为负值。
在25℃、100℃、300℃、500℃、700℃温度下,不同冲击载荷下,试样的应力应变曲线大致相同。25℃时,曲线有一定的右倾趋势,且随着冲击载荷的增大,试样的变形能力逐渐增强。100~500℃,曲线表现出一定的左移倾向。500℃下,当冲击荷载为0.6 MPa时,因试样间的不均质性出现的偶然误差而导致试样应力应变曲线明显出现偏差;曲线塑性变形较为明显,整个试样表现出显著的弹塑性状态。700℃下压密阶段和塑性变形阶段都很小,弹性变形占比远大于塑形变形和初始压密阶段,故而曲线峰值区域较为尖锐。从图中观察,25℃、100℃、300℃、500℃、700℃条件下,试样弹性变形阶段曲线长度均大于塑性变形阶段长度。
3.3. 动态峰值应力的变化规律
25℃~700℃黄砂岩试样动态峰值应力随冲击载荷的变化情况,如图5所示。
Figure 5. Diagram of dynamic peak stress and impact load
图5. 动态峰值应力和冲击载荷的关系图
由图5可知,不同温度条件下,黄砂岩试样的动态峰值应力均随着冲击载荷的增大而增大,即冲击载荷对黄砂岩的单轴抗载荷有加强作用。
动态峰值应力与冲击载荷呈二次函数关系,拟合关系式为
(3)
式中a、b、c为拟合相关参数。不同温度黄砂岩试样拟合参数取值如表1所示。
25℃时,冲击载荷从0.2 MPa增加到0.6 MPa,试样的真实应力从95.18 MPa增加到243.10 MPa,增加幅度达155.41%,说明冲击载荷对常温条件下的黄砂岩的单轴抗压强度具有明显的加强作用。
Table 1. Fitting parameters for the peak stress and impact load
表1. 峰值应力和冲击载荷拟合参数取值
此外,不同温度对试样产生的热力学效应不同,300℃、500℃、700℃时的真实应力关系曲线均位于25℃曲线上方,此温度范围内,真实应力均高于同种冲击载荷下的常温真实应力。而300℃对试样强度影响最为显著,在冲击载荷0.2 MPa~0.6 MPa中,真实应力分别达到134.44 MPa、241.49 MPa、251.99 MPa、292.94 MPa,比同冲击载荷下常温的试样真实应力分别增加42.32%、76.86%、40.00%、23.54%、20.50%。
100℃真实应力关系曲线位于最下方,即曲线值最小。在冲击载荷0.2 MPa~0.6 MPa中,真实应力分别达到92.96 MPa、125.34 MPa、142.05 MPa、191.14 MPa、209.94 MPa,比同冲击载荷下常温的试样真实应力分别减少2.32%、8.21%、21.08%、−8.69%、13.64%。说明100℃对黄砂岩试样单轴抗压强度具有弱化作用。
不同冲击载荷作用下的黄砂岩试样动态峰值应力随着温度变化情况如图6所示。从图中可以看出,不同冲击载荷作用下,黄砂岩试样动态峰值应力随温度的变化曲线总体呈现W型。且随着冲击载荷的增大,曲线整体呈上移趋势。25℃~100℃时,除0.5 MPa曲线,曲线均表现为下降走势。100℃~500℃,曲线表现为先增大后减小,且变化幅度较大,100℃~300℃曲线整体上升,在300℃时均到达曲线的最高点;300℃~500℃,曲线下降;除0.6 MPa曲线,在100℃~500℃曲线上升速率大于下降速率。500℃~700℃时,除0.6 MPa外,其余曲线较为平缓。说明在500℃~700℃时,温度对黄砂岩试样的单轴抗压强度影响较小。
Figure 6. Diagram of dynamic peak stress and temperature
图6. 动态峰值应力和温度的关系图
同种冲击载荷下,25℃ ≤ T ≤ 100℃和500℃ ≤ T ≤ 700℃时,黄砂岩试样动态峰值应力受温度的影响明显小于100℃ ≤ T ≤ 500℃时的情况。其中,300℃时黄砂岩的动态峰值应力最大,大于其他温度的同种冲击载荷的情况。说明黄砂岩试样动态峰值应力受温度作用的影响明显,并且在300℃时最为显著。
3.4. 试件破碎形态
本文选取25℃不同冲击荷载下的破碎情况和0.6 MPa不同温度下的破碎情况进行分析讨论,破碎情况如图7所示。
相同温度下,随着冲击载荷的增大,试样的破裂面逐渐增多,破坏程度加剧,碎块尺寸减小。25℃时,随着冲击载荷的增大,碎块逐渐变细,当冲击载荷为0.4 MPa时,黄砂岩试样已有大面积粉末状态,当荷载为0.5 MPa时,黄砂岩试样已变成完全粉末状态。
相同冲击载荷下,随着温度的升高,黄砂岩试样破碎程度先减小,到达临界温度后,破碎程度呈增大趋势。0.6 MPa冲击荷载条件下,可以明显看出,25℃~300℃时,黄砂岩试样破碎状态分别为粉末状态、破碎细颗粒状态、破碎粗颗粒状态,破碎程度逐渐下降;300℃~500℃,黄砂岩试样的破碎程度缓缓加剧,由破碎粗颗粒状态变为破碎细颗粒状态;但整体破碎程度不如25℃~300℃;300℃为临界温度。700℃条件下,黄砂岩试样保持高度完整状态。同时,随着温度的升高,黄砂岩试样逐渐由黄色变成棕红色;且温度越高,颜色变化效应越明显。
对于以上现象我们进行了以下分析。相同温度下,冲击载荷越大,初始动能越高,吸收的能量越多,裂隙扩展越大,碎块数目越多,破碎程度越大。这是因为岩石内部大小不同的裂隙逐渐发育、扩展、贯通最终导致的。温度对黄砂岩试样的力学性质有两方面的影响,一方面,温度升高使得黄砂岩试样热熔变形,试样内部原生裂隙逐渐愈合,裂隙数量变少,密实性得到提升,从而提升了试样的承载能力和抗变形能力;另一方面,不同热膨胀率引起跨颗粒边界的热膨胀不协调引起结构热应力,试样内部产生微裂缝,同时胶结物的形成也会影响试样的变形,使得试样的承载能力和抗变形能力下降 [26]。
4. 结论
(1) 黄砂岩在不同冲击载荷作用下,动态应力应变曲线可分为4个阶段,分别为压密阶段、弹性变形阶段、塑性变形阶段、损伤破坏阶段。从整体上看,黄砂岩应力应变曲线经过先上升,再平缓,达到峰值后,再下降的发展趋势。在25℃、100℃、300℃、500℃、700℃温度下,不同冲击载荷下,试样的应力应变曲线大致相同。各温度作用下,试样弹性变形阶段曲线长度均大于塑性变形阶段长度。
(2) 不同温度条件下,黄砂岩试样的动态峰值应力均随着冲击载荷的增大而增大。300℃、500℃、700℃曲线均位于25℃曲线上方,此温度范围内,真实应力均高于同种冲击载荷下的常温真实应力。根据100℃真实应力曲线关系,说明100℃对黄砂岩试样单轴抗压强度具有弱化作用。
(3) 不同冲击载荷作用下,黄砂岩试样动态峰值应力随温度的变化曲线总体呈现W型;且随着冲击载荷的增大,
曲线整体呈上移趋势。同种冲击载荷下,25℃ ≤ T ≤ 100℃和500℃ ≤ T ≤ 700℃时,黄砂岩试样动态峰值应力受温度的影响明显小于100℃ ≤ T ≤ 500℃时的情况。黄砂岩试样动态峰值应力受温度的影响明显,在300℃时最显著。
(4) 相同温度下,随着冲击载荷的增大,试样的破裂面逐渐增多,破坏程度加剧,碎块尺寸减小。相同冲击载荷下,随着温度的升高,黄砂岩试样破碎程度先减小,到达临界温度300℃后,破碎程度呈缓慢增大趋势,整体上25℃~300℃破碎程度比300℃~700℃破碎程度高。
项目基金
国家自然科学基金面上项目编号(12072363)。
NOTES
*通讯作者。