冲击荷载下煤系砂岩应变率效应及能量耗散特征

余永强1,2,张文龙1,范利丹1,2,龚 健1,杨 蒙1,孙 亮1

(1.河南理工大学 土木工程学院,河南 焦作 454003; 2.河南省地下空间开发及诱发灾变防治国际联合实验室,河南 焦作 454003)

摘 要:在煤矿巷道掘进过程中,巷道围岩在动载作用下变形将会增大,为研究煤系砂岩在冲击荷载作用下的力学特性及能量耗散,以河南陈四楼煤矿巷道围岩中的砂岩为研究对象,利用直径为50 mm的分离式霍普金森压杆试验装置对煤矿砂岩开展单轴单次冲击压缩试验和循环冲击压缩试验,对冲击荷载作用下煤矿砂岩的应变率效应、能量耗散特征和破坏模式等进行分析。研究结果表明:在单轴单次冲击荷载作用下,随着平均应变率的增加,砂岩试样的峰值应力和峰值应变均增大,割线模量逐渐降低,砂岩试样的塑性增加,强度提高;且峰值应变与平均应变率呈线性递增关系,峰值应力近似与平均应变率的1/3次幂呈递增关系;随着平均应变率的增加,砂岩试样的单位体积吸收能呈线性增加趋势,且试样破碎程度不断增大,在压应力持续作用下砂岩试样内部裂纹不断交叉扩展,沿轴向发生劈裂破坏。在循环冲击荷载作用下,随冲击荷载作用次数的增加,砂岩试样的平均应变率和峰值应变均逐渐增大,峰值应力、割线模量和第2类割线模量均随着冲击次数的增加而逐渐降低;在固定冲击气压下进行循环冲击时,随着冲击次数的增加,入射能基本保持不变,反射能和吸收能均逐渐增大,透射能逐渐减小,砂岩试样的单位体积吸收能随冲击次数增加呈先缓慢增加后突然增加的趋势;在循环冲击荷载作用下,砂岩试样吸收能量的能力增强,随着冲击次数的不断增加,试样内部微裂纹分布增多且扩展速度加快,主要发生受拉破坏,沿轴向加载方向产生贯穿于砂岩试样的环状破坏面。

关键词:煤系砂岩;分离式霍普金森压杆;能量耗散;应变率;破坏模式

在煤炭资源开采过程中,煤矿巷道围岩不可避免的遭受到较为频繁的爆破以及采矿冲击地压等强烈的动荷载作用,严重威胁着煤矿开采安全[1-4]。另一方面,在实际工程中多数岩体的破坏并不是在单次冲击荷载作用下发生的,如炮弹对目标的多次打击、工程岩体的多次爆破以及隧道的多次爆破开挖、冲击破岩等,都是在多种形式荷载作用下使岩体逐步产生损伤断裂,最终发生宏观破坏[5-9]。因此,研究动荷载作用下煤矿岩石的应变率效应、能量耗散特征和破坏模式等,对于更好的理解和控制煤矿巷道围岩稳定性的动力学问题以及保证煤矿安全经济生产具有重要意义。

国内外相关学者利用分离式霍普金森压杆系统(Splitting Hopkinson Pressure Bar,简称SHPB)对岩石在动荷载作用下的力学特性和破坏模式进行了一些研究。LI等[10]利用摆锤驱动的分离式SHPB装置对绿色砂岩进行反复冲击压缩试验,研究了材料的声学特性、耗散能、变形特性和微观结构演变。BRAUNAGEL M J等[11]利用SHPB试验装置对花岗岩在循环加载过程中的动态抗压强度及破坏机理进行研究,认为在循环加载条件下试样的抗压强度会发生降低,且岩石的强度随加载速度的增加而增加,在极快的加载速度下岩石的破坏由简单断裂变为破碎状态。朱晶晶等[12]对花岗岩试样开展单轴循环冲击压缩试验,分析花岗岩的力学特性及能量吸收规律,认为随冲击荷载循环作用次数的增加,岩石的变形模量变小、峰值应力降低,试样的屈服应变逐渐增大,并根据试验结果建立了岩石动态统计损伤本构模型。李地元等[13]结合岩石声发射监测研究了花岗岩在循环冲击压缩试验下的累计损伤特性。张智宇等[14]对深部白云岩开展单轴冲击试验,从应力-应变、强度和能量分布等方面对白云岩的动态力学性质进行分析,认为岩石在冲击荷载下发生拉伸破坏。王志亮等[15-16]分析了花岗岩在等幅循环冲击下的峰值应力、应变、平均应变率、损伤值与冲击次数之间的关系,以及轴压对试样总循环冲击次数的影响;金解放[17-19]、于洋[20]等利用改进的SHPB试验装置对砂岩试样进行循环冲击试验,研究了轴压和循环冲击次数对岩石动态应力-应变曲线、动态强度和变形特性的影响,认为在无围压作用下砂岩发生脆性破坏,有围压时表现出明显的弹塑性特征。唐礼忠等[21]利用改进的SHPB试验装置对矽卡岩岩样进行一维静载及循环冲击共同作用下的加载试验,研究了岩石的动态应力-应变曲线、能量传递以及破坏模式等。许金余等[22]利用带围压装置的SHPB对斜长角闪岩、绢云母石英片岩和砂岩开展循环冲击试验,认为在循环冲击荷载下岩石的杨氏模量明显减小,并以杨氏模量表征岩石的损伤程度。李晓锋等[23]对灰岩、白云岩和砂岩进行动态冲击试验,得到了岩石的动强度因子、耗散能密度及破碎尺寸与应变率间的变化关系。目前,对于煤矿中的砂岩在动载作用下的力学性质和能量耗散虽有一定的研究,但对岩石试样在循环冲击荷载作用下的能量耗散和试样破坏模式研究较少。

在巷道掘进过程中,巷道围岩在冲击动载和循环动载作用下,变形将会逐渐增大,影响巷道的稳定性。针对上述研究中存在的不足,以陈四楼煤矿巷道围岩中的砂岩为研究对象,利用改进的SHPB试验装置对煤矿巷道砂岩开展冲击试验,分析砂岩在冲击荷载作用下的动力学特征,总结其力学特性和能量耗散规律、破坏形态变化特征等,对正确认识煤矿巷道围岩在冲击力学作用下的动态力学特性、准确评价巷道围岩稳定性具有重要意义,为煤矿的合理高效开采及围岩支护稳定性等相关实际工程问题提供理论参考。

1 砂岩SHPB试验

1.1 试样制备

试样取自河南陈四楼煤矿巷道中的岩体,选取完整性和均质性较好的砂岩作为研究对象。根据国际岩石力学学会(ISRM)推荐的试样尺寸[24],将砂岩加工为φ50 mm×25 mm的圆柱体试样,试样两端面的不平行度和不垂直度均小于0.02 mm,试样如图1所示。砂岩试样的静态物理力学参数见表1。

图1 冲击试验砂岩试样
Fig.1 Sandstone samples for impact tests

表1 砂岩静态物理力学参数
Table 1 Static physical mechanics parameters of sandstone

密度/(kg·m-3)纵波波速/(m·s-1)抗压强度/MPa抗拉强度/MPa弹性模量/GPa泊松比2 7034 07896.494.809.640.17

1.2 试验装置及原理

利用如图2所示的河南理工大学冲击实验室改进的SHPB试验装置对砂岩试样进行单轴冲击试验,采用江苏东华生产的DH8302超动态应变仪进行数据采集。试验所采用的子弹为直径37 mm、长度400 mm的钢杆,入射杆为直径50 mm、长度2 400 mm的钢杆,透射杆为直径50 mm、长度1 200 mm的钢杆。入射杆和透射杆均为变截面锥形杆,钢杆中纵波波速为5 190 m/s,弹性模量为210 GPa,吸收杆材质为铝杆。

图2 SHPB试验装置系统和砂岩试样放置位置
Fig.2 SHPB test facility system and placement of sandstone
samples

试验中典型试样两端的应力结果如图3所示,可以看出,入射波加反射波后形成的曲线与透射波曲线基本重合,表明在试验过程中试样左右两端受力相等,满足应力平衡条件[25]

图3 砂岩动态应力平衡检验
Fig.3 Calibration test of dynamic stress balance of sandstone

基于一维应力波和应力-应变均匀性假定,根据入射杆上粘贴的应变片所测得的入射波应变εI(t)及反射波应变εR(t)和透射杆上应变片测得的透射波应变εT(t)等应变信号,根据“三波法”计算公式可分别计算出试样两端面上的应力σ、应变ε以及应变率其中“三波法”计算公式[26]

(1)

(2)

(3)

式中,A0为压杆的横截面积,mm2;AS为试样的横截面积,mm2;E0为压杆的弹性模量,GPa;C0为压杆的纵波波速,m/s;LS为试样的长度,mm。

根据应力均匀性假设,试样长度远小于输入应力波的波长,可认为试样中的应力均匀化、无衰减,由一维应力波理论可得εI(t)+εR(t)=εT(t),选取入射波εI(t)和透射波εT(t)应变信号进行计算,则式(1)~(3)可简化为“二波法”[26]计算公式

(4)

(5)

(6)

1.3 冲击气压确定

对试验采集得到的原始波形进行分析计算,得到砂岩试样在冲击荷载作用下的峰值应力、峰值应变等参数。单次冲击时不同冲击气压P与砂岩的峰值应力σf和平均应变率间的关系如图4所示,可以发现,随冲击气压的增大,平均应变率和峰值应力呈增大的趋势,结合砂岩试样的静态物理力学参数和动载作用下试样的强度,确定循环冲击气压的大小为0.40,0.45和0.50 MPa。

图4 峰值应力和平均应变率与冲击气压关系
Fig.4 Relationship between peak stress,average strain rate
and impact pressure

1.4 试验方法

试验时在子弹与入射杆之间采用宽度为15 mm、厚度0.3 mm的正方形H62黄铜片作为波形整形器,试样放置于入射杆和透射杆之间并在试样与压杆接触面涂抹适量黄油,减小压杆和试样接触面间的端面摩擦效应。试验首先选取不同的冲击气压对砂岩试样开展单轴单次冲击压缩试验。循环冲击试验时,在一次冲击后,若试样未破碎,则在相同的冲击气压下,进行下一次冲击,经多次冲击作用,直至砂岩试样产生宏观力学破坏,停止试验。砂岩试样在单次和循环冲击作用下主要试验结果分别见表2,3。

表2 砂岩单次冲击试验结果
Table 2 Experimental results of sandstone under single impact loads

编号L/mmD/mmρ/(kg·m-3)P/MPaε·/s-1εfσf/MPaDIFE50/GPaE250/GPa124.7949.772 6770.4018.990.004 2100.001.04136.4229.64224.8249.692 6900.5021.810.005 0118.601.23111.8135.84325.1549.772 6740.7034.740.007 0143.901.49106.835.32425.1049.672 7380.8047.990.009 8149.101.5584.1124.34524.9849.672 7211.2075.700.016 4165.101.7142.7529.22

注:L为试样长度;D为试样直径;ρ为试样密度;εf为峰值应变;DIF为应力增长系数;E50为割线模量;为第2类割线模量。

2 SHPB冲击试验结果分析

将砂岩试样共分为2组,分别开展不同冲击气压下的单次冲击试验和固定冲击气压下的循环冲击试验,相同加载类型下重复冲击多次,选取有代表性试验数据进行分析。以砂岩试样在0.45 MPa固定冲击气压下循环冲击时应变片测得的典型冲击波形图为例,如图5所示,可以看出,每次冲击时入射波波形基本重合,能够保证每次冲击时试样所受到的冲击荷载基本相同。

图5 0.45 MPa固定冲击气压循环冲击时砂岩应变波形
Fig.5 Strain wave time history curves of sandstone under cyclic
impact at 0.45 MPa constant impact pressure

2.1 冲击荷载作用下砂岩试样的破坏形态

不同冲击气压下单次冲击时砂岩试样的破坏形态如图6所示,在冲击气压较低时试样的变形较小破碎为几大块,随冲击气压的提高平均应变率逐渐增大,试样的损伤范围由试样边缘分布逐渐扩展到试样中心部位[27],试样的破碎程度增大,由大块破碎过渡为小块粉碎,碎块数量增多且碎块体积逐渐减小。随着冲击气压进一步增大,砂岩试样内部裂纹分布增多且裂纹宽度逐渐增大并贯穿于整个砂岩试样,在应力持续作用下砂岩试样裂纹出现交叉发展,发生受拉破坏和压剪破坏[28]

表3 砂岩循环冲击试验结果
Table 3 Experimental results of sandstone under cyclic impact loads

编号L/mmD/mmρ/(kg·m-3)P/MPanε·/s-1εfσf/MPaDIFE50/GPaE250/GPa111.710.002 9103.901.08376.8248.71214.730.003 3105.201.09144.8932.38313.860.003 4104.101.0895.7834.78125.2249.672 7080.40415.760.003 7104.401.08128.3938.09521.360.004 5101.201.0573.2933.02625.010.005 598.601.02142.3327.79735.890.008 391.730.9548.2117.38132.250.007 2109.201.13138.4639.51224.2549.712 7310.45230.170.006 3105.501.09148.2731.46350.850.010 4103.701.0741.7123.434120.250.024 064.740.6726.0113.24119.080.004 2116.701.21180.8027.55324.7149.742 6800.50228.560.006 4114.901.1975.0425.28359.290.012 883.730.87123.623.96

注:n为冲击次数。

图6 不同冲击气压下单次冲击时砂岩的破碎形态
Fig.6 Fracture morphology of sandstone under single impact
with different impact pressures

固定冲击气压下循环冲击时砂岩试样的破坏形态如图7所示,可以看出,随冲击气压增大,砂岩出现宏观破坏时的冲击总次数减小。在循环冲击荷载作用下砂岩试样主要是在拉应力作用下发生受拉破坏,在0.40 MPa冲击气压下,第1次冲击后砂岩试样基本保持完整状态,随着冲击次数增加,砂岩试样内部损伤不断累积,微裂纹不断产生与扩展,试样沿轴向加载方向分裂为几大块。在0.50 MPa冲击气压下,砂岩试样在第1次冲击后表面出现裂纹,在第3次冲击时砂岩内部损伤严重,微裂纹分布增多且扩展速度加快,裂纹宽度逐渐增大并贯穿整个砂岩试样,试样的承载力迅速下降,破碎程度增大,发生受拉破坏,沿轴向加载方向产生贯穿于砂岩试样的环状破坏面。

图7 循环冲击荷载作用下砂岩的破坏形态
Fig.7 Failure morphology of sandstone under cyclic impact loads

2.2 砂岩试样应力、应变和应变率分析

2.2.1 单次冲击时砂岩试样强度及变形分析

砂岩试样在不同冲击荷载作用下单次冲击时的典型应力-应变曲线如图8所示,可以看出,砂岩试样的峰值应力和峰值应变均随平均应变率的增加而增加。当平均应变率为75.70 s-1时,砂岩试样的峰值应变为0.016 4,峰值应力为165.10 MPa,约为1.71 fc(fc为砂岩试样的静载平均抗压强度,为96.40 MPa),砂岩试样的塑性增加,强度提高。这是因为在更大冲击荷载作用时试样内部产生了更多的微裂纹,破坏形态也由较为完整的大块破碎逐渐过渡到小块甚至粉碎状态,试样的变形增大,吸收能量增加,因此砂岩试样所能承受的荷载增加,表现为强度增大。砂岩试样的动态应力-应变曲线可分为弹性阶段、裂缝扩展阶段、塑性变形阶段和破坏阶段4个阶段,在初始弹性阶段随应变增加应力呈快速线性增加趋势,且在初始弹性阶段曲线基本重合。在冲击荷载作用下试样内部微裂缝不断萌生发育,达到弹性极限应力后进入塑性变形阶段,当达到屈服应力后,试样内部裂纹扩展贯通,应变增大且应力迅速降低,砂岩试样发生破坏[29]

图8 不同冲击气压下砂岩试样的应力-应变曲线
Fig.8 Stress-strain curves of sandstone under different pressures

不同冲击气压下单次冲击时砂岩试样的峰值应力和峰值应变与平均应变率关系如图9所示,可以看出,砂岩的峰值应变随应变率增加基本呈线性增加趋势,峰值应力近似与平均应变率的1/3次幂呈递增关系。

图9 不同冲击气压下峰值应变和应力与平均应变率关系
Fig.9 Relationship between peak strain,peak stress and average
strain rate under different impact pressures

2.2.2 循环冲击时砂岩试样强度及变形分析

循环冲击时砂岩试样的应力-应变曲线如图10所示,可以看出,随冲击次数增加,砂岩试样的峰值应力逐渐降低,应变逐渐增大,前几次冲击时应力应变曲线基本重合。试样内部存在的原始缺陷如微裂纹、空隙等会引起应力集中,在冲击荷载作用下裂纹扩展发育,导致岩石组构传递荷载的能力降低,随着冲击次数的增加,应力-应变曲线初始弹性上升阶段的斜率逐渐减小,即弹性模量逐渐降低,岩石试样的力学性能不断下降[12,22]。随着冲击气压的增大,岩石破坏时所需冲击总次数逐渐减小,在0.40 MPa冲击气压下,岩石试样在多次冲击作用下内部裂纹不断闭合、张开,试样累计损伤增加,最终试样发生破坏。在0.45 MPa与0.50 MPa冲击气压下,试样破坏时所需的冲击总次数相差不大,这可能是由于此时的冲击气压较大,在第1次冲击时对试样造成的损伤程度过大,超过了岩石的临界破坏值,试样受损严重,因此在这2种冲击气压下冲击总次数相差不大。

图10 循环冲击荷载作用下砂岩试样的应力-应变曲线
Fig.10 Stress-strain curves of sandstone under cyclic
impact loads

图11(a)为循环冲击荷载作用下砂岩的平均应变率与冲击次数间的关系。可以看出,循环冲击时冲击气压越高,砂岩破坏时的冲击总次数越少。3种冲击气压下砂岩试样的平均应变率均随着冲击荷载作用次数的增加先缓慢增加,达到一定冲击次数时由于试样内部损伤严重,变形增大,平均应变率突然增大。图11(b)为3种冲击气压循环冲击时砂岩的峰值应变与冲击次数之间的关系,可以看出,砂岩试样的峰值应变随冲击次数的增加逐渐增大。在0.40 MPa冲击气压时,前4次冲击时砂岩试样的峰值应变增加缓慢,在第7次冲击时峰值应变较第6次增加较多,0.50 MPa冲击气压下,砂岩试样在第3次冲击时的峰值应变突然增大。图11(c)为3种冲击气压循环冲击时砂岩的峰值应力与冲击次数之间的关系,砂岩试样的峰值应力随冲击次数的增加总体呈下降趋势,在0.40 MPa冲击气压下,砂岩的峰值应力随着冲击次数的增加呈先增大后减小的趋势,这可能是因为冲击气压较低,第1次冲击后试样内部并未发生严重损伤,试样内部微裂缝被压密从而提高了砂岩的强度,在第2次冲击时砂岩的峰值应力略大于第1次冲击时的峰值应力,随着冲击次数继续增加峰值应力缓慢下降,在第7次冲击时峰值应力突然下降。

图11 峰值应变、峰值应力和平均应变率随冲击次数
增加变化规律
Fig.11 Variation of peak strain,peak stress and average strain
rate along with the increase of impact times

0.45 MPa和0.50 MPa冲击气压下砂岩试样的峰值应力随冲击次数的增加呈先缓慢下降后突然下降的趋势。每次冲击荷载作用后均有一部分能量用于加剧试样损伤,随着冲击次数不断增加砂岩试样内部会产生更多的微裂纹,内部损伤不断累积导致岩石的承载力不断降低,变形和平均应变率增大,达到一定冲击次数时,砂岩峰值应力突然下降。

砂岩试样在3种循环冲击荷载作用下的平均应变率与其对应的峰值应变和峰值应力间的关系如图12所示,其中,分别为0.40,0.45,0.50 MPa冲击气压下对应的动态峰值应变;分别为0.40,0.45,0.50 MPa冲击气压下对应的动态峰值应力。在3种循环冲击气压作用下,砂岩试样的峰值应变随平均应变率的增大呈线性增加趋势,峰值应力随平均应变率的增大呈线性降低趋势。

图12 峰值应变和峰值应力与平均应变率关系
Fig.12 Relationship between peak strain,peak stress and average strain rate

2.2.3 单次冲击和循环冲击时割线模量分析

以割线模量作为岩石动态应力-应变曲线的衡量指标进行对比[30],即

E50=σd50d50

(7)

式中,σd50为动态峰值应力50%时的应力值,MPa;εd50σd50对应的峰值前的轴向应变值。

为进一步观察岩石模量的变化,将岩石的应力-应变曲线简化为2段,如图13所示,并定义第2类割线模量,图13中,εdσd分别为试样在动态冲击作用下试样的峰值应变和峰值应力。即取应力-应变曲线上峰值应力和峰值应力的50%两点之间的斜率,用表示[30]。其中,割线模量对岩石弹性变形阶段的力学特征进行评价,第2类割线模量对岩石塑性变形阶段和峰值阶段的力学特性进行评价,计算结果见表2,3。

图13 岩石割线模量示意
Fig.13 Schematic of rock secant modulus

不同冲击气压下单次冲击时,砂岩试样割线模量和第2类割线模量与平均应变率间的关系如图14所示。可以看出,在单次冲击荷载作用下,随平均应变率的增大,割线模量逐渐降低,第2类割线模量受平均应变率的变化影响较小。在更大冲击气压作用下砂岩内部损伤加剧,微裂纹扩展速度和宽度增加,从而弱化微元体的承载能力,砂岩试样变形增大,平均应变率也逐渐增大。随平均应变率增大,割线模量逐渐降低,式(8)为单次冲击时砂岩试样的割线模量及第2类割线模量与平均应变率间拟合关系。

图14 单次冲击割线模量与平均应变率关系
Fig.14 Relationship between secant modulus and average
strain rate under single impact loads

(8)

循环冲击荷载作用下砂岩试样的割线模量及第2类割线模量随冲击次数增加时的变化趋势如图15所示,可以看出,割线模量及第2类割线模量随冲击次数的增加整体上均呈下降的趋势,割线模量随冲击次数的增加表现出一定的离散性,第2类割线模量在数值上比割线模量小,随冲击次数的增加整体呈下降趋势。随冲击次数的增加,割线模量的降低表明岩石的强度在衰减,内部微裂纹和损伤程度在逐渐增加,弹性能也在发生变化。割线模量降幅明显较大,而第2类割线模量降幅很小。但割线模量对应的是弹性势能,说明弹性势能变化更大,而塑性势能变化较小。

图15 循环冲击时割线模量与冲击次数关系
Fig.15 Relationship between secant modulus and impact
times under cyclic impact loads

2.3 冲击荷载作用下能量耗散特征分析

在SHPB冲击压缩试验过程中,在试样两端面和杆与试样的接触面均涂抹了适量的黄油作为润滑剂,因此在计算时可忽略试样与压杆之间因端部摩擦所消耗的能量。根据SHPB试验加载原理和能量守恒定律,利用入射波εI(t)、反射波εR(t)和透射波εT(t)的应变曲线,通过以下公式可分别求出各项能量[3]:

(9)

(10)

(11)

WS=WI-WR-WT

(12)

(13)

式中,WIWRWTWS分别为冲击过程中的入射能、反射能、透射能和吸收能,J;σI为入射应力,MPa;σR为反射应力,MPa;σT为透射应力,MPa;VS为试样体积,cm3;EV为单位体积吸收能,J/cm3

对入射杆和透射杆的应变数据进行分析和计算,得到砂岩试样在单次冲击和循环冲击过程中的能量变化计算结果,分别见表4,5。

表4 砂岩单次冲击能量计算结果
Table 4 Calculation results of single impact energy of sandstone

试样编号入射能WI/J反射能WR/J透射能WT/J吸收能WS/J单位体积吸收能/(J·cm-3)177.785.6660.3111.810.2452104.896.3381.6516.910.3513149.7311.66110.0328.040.5734175.2011.62116.1647.420.9755253.9136.47142.0175.431.558

表5 砂岩循环冲击能量计算结果
Table 5 Calculation results of cyclic impact energy of sandstone

试样编号冲击次数n入射能WI/J反射能WR/J透射能WT/J吸收能WS/J单位体积吸收能/(J·cm-3)178.985.1567.296.540.134281.662.5868.0811.000.225379.362.5265.8510.990.2251479.063.1364.0311.900.244577.326.3059.7811.240.230680.327.3157.0715.940.326779.336.7747.7524.810.508191.489.8460.6720.970.4462291.136.8563.9720.310.432391.349.7149.532.130.683484.5160.8211.2212.470.2651103.234.9483.2415.050.31332104.256.7474.2723.240.4843105.4717.4746.3541.650.867

2.3.1 不同冲击气压下单次冲击时能量变化规律

在不同冲击气压作用下进行单次冲击时,砂岩试样的平均应变率与入射能、反射能、透射能和吸收能之间的关系如图16所示。可以看出,随平均应变率的增大,入射能、反射能、透射能和吸收能均是增大的,各能量与平均应变率间的拟合关系为

图16 单次冲击时能量与平均应变率关系
Fig.16 Relationship between energy and average strain
rate at single impact

(14)

砂岩试样的平均应变率与单位体积吸收能间的关系如图17所示,可以看出,随平均应变率的增大砂岩单位体积吸收能基本呈线性增加的趋势,结合砂岩试样冲击后的破坏形态,可以认为,当单位体积吸收能大于0.573 J/cm3时砂岩试样发生较为明显的宏观破碎。随冲击气压增大入射能增加,砂岩试样由较完整状态逐渐向破碎及粉碎状态过渡,在这个过程中需要吸收更多的能量用于试样内部微裂纹的萌生、扩展和贯通,因此单位体积吸收能随平均应变率增加逐渐增大。单位体积吸收能与平均应变率间拟合关系为

图17 单位体积吸收能与平均应变率关系
Fig.17 Relationship between absorbed energy per unit volume
and average strain rate

(15)

2.3.2 固定冲击气压下循环冲击时能量变化规律

固定冲击气压下进行循环冲击时的入射能、反射能和透射能与冲击次数间的关系如图18所示。可以看出,冲击气压为0.40,0.45和0.50 MPa时对应的入射能分别约为80,90和100 J。在固定冲击气压下进行循环冲击时,入射能随着冲击次数的增加基本不变,说明每次冲击时入射波是基本相同的,重复性较好。随着冲击次数的增加,反射能逐渐增大,透射能逐渐降低,吸收能逐渐增大。当试样内部损伤严重时吸收能也会出现降低的趋势,如在0.45 MPa冲击气压下第3次冲击时,试样四周边缘产生裂缝,内部损伤严重,在此基础上继续冲击时,试样发生严重破碎且所吸收的能量相对降低,大部分能量以反射能的形式释放。每次冲击时砂岩试样内部会受到反射应力与透射应力的共同作用,随着冲击次数的增加,试样入射端反射系数则不断增大而透射系数不断减小,导致反射波形不断增大而透射波形不断减小,即反射能随冲击次数增加而不断增大,透射能随冲击次数增加而不断减小,吸收能则不断增加从而用于试样内部微裂缝的萌生、扩展和发育,达到一定损伤程度时微裂缝将贯穿整个砂岩试样,试样发生受拉破坏。

图18 循环冲击荷载作用下入、反、透、吸收能随冲击
次数变化规律
Fig.18 Variation of different energy along with cycle
impact times

固定冲击气压下循环冲击时砂岩试样的单位体积吸收能与冲击次数间的变化趋势如图19所示。可以看出,随冲击次数增加,单位体积吸收能呈先缓慢增加后突然增加的趋势,当试样内部损伤严重时,随冲击次数继续增加,由于大部分能量以反射能的形式释放且试样吸收的能量较低,单位体积吸收能也可能出现降低趋势,这是因为试样内部损伤严重导致吸能能力降低,如图19中0.45 MPa冲击气压加载情况所示。在0.40 MPa冲击气压下,第2次冲击时的单位体积吸收能0.225 J/cm3约为第1次0.134 J/cm3的1.68倍,第2~5次冲击时单位体积吸收能变化不大,在第6次冲击时单位体积吸收能突然增加,为0.326 J/cm3,约为第5次0.230 J/cm3的1.42倍,是第1次冲击时的2.43倍。当冲击气压为0.45 MPa和0.50 MPa时试样吸收更多的能量,第1次冲击时砂岩的单位体积吸收能分别为0.446 J/cm3和0.313 J/cm3,均大于0.40 MPa冲击气压下第1次冲击时的单位体积吸收能0.134 J/cm,在更大冲击气压作用下砂岩试样内部损伤加剧。

图19 循环冲击荷载作用下砂岩单位体积吸收能随冲击
次数变化规律
Fig.19 Variation of sandstone absorbed energy per unit volume
along with impact cycle times

2.3.3 固定冲击气压下循环冲击时累计比能量吸收值

定义累计比能量吸收能ϑ(J/cm3)为

ϑ

(16)

因此可得砂岩的累计比能量吸收值与冲击荷载作用次数的关系如图20所示。随着冲击次数增加,砂岩试样的累计比能量吸收值逐渐增加,试样内部损伤不断累积,吸收能增大,达到一定程度时吸收能快速增加,砂岩试样的累计比能量吸收值也逐渐增加。在0.40 MPa冲击气压时,第6次和第7次冲击时的累计比能量吸收值急剧增大,第7次冲击时为1.891 J/cm3;在0.50 MPa冲击气压时,第3次冲击时的累计比能量吸收值为1.665 J/cm3,是第3次冲击时砂岩单位体积吸收能0.867 J/cm3的1.92倍。

图20 砂岩累计比能量吸收值随冲击次数的变化规律
Fig.20 Variation of cumulative specific energy absorption
value of sandstone along with cyclic impact times

在不同冲击气压下,单次冲击时砂岩试样的单位体积吸收能大于0.573 J/cm3时发生破碎;在循环冲击时,0.40,0.45和0.50 MPa冲击气压下最后1次冲击时砂岩的单位体积吸收能分别在0.508,0.265和0.867 J/cm3时发生破坏,此时砂岩的累计比能量吸收值分别为1.891,1.825和1.665 J/cm3,而0.45 MPa冲击气压下前3次冲击累计比能量吸收值为1.560 J/cm3,故总体可认为在循环冲击下,当砂岩试样的累计比能量吸收值大于1.5 J/cm3时砂岩试样发生破坏,在循环冲击荷载下砂岩发生破坏时吸收的总能量增加。

3 结 论

(1)单次冲击荷载下,砂岩的峰值应力和峰值应变均与平均应变率呈正相关关系,砂岩的动态抗压强度最大为165.10 MPa,约为试样静载抗压强度的1.71倍;在循环冲击荷载下,随冲击作用次数的增加,砂岩试样的峰值应变和平均应变率增大,峰值应力减小,且峰值应力和峰值应变皆与平均应变率间存在较好的线性相关性。

(2)单次冲击荷载下,随着平均应变率的增大,砂岩试样的割线模量逐渐降低,循环冲击时随冲击作用次数的增加,砂岩试样的割线模量及第2类割线模量降低。

(3)单次冲击荷载下,砂岩单位体积吸收能随平均应变率增加呈线性增加趋势,当单位体积吸收能大于0.573 J/cm3时砂岩试样发生破碎;循环冲击时,砂岩单位体积吸收能随冲击次数增加先缓慢增加后突然增加,当砂岩试样的累计比能量吸收值大于1.5 J/cm3时砂岩试样发生破坏,并且在循环冲击荷载下砂岩发生破坏时吸收的总能量增加。

(4)单次冲击荷载下,砂岩试样随冲击气压的提高,由大块破碎过渡为小块粉碎,发生受拉破坏和压剪破坏;固定气压下循环冲击时,随冲击作用次数增加,砂岩试样由基本完整状态向大块破碎状态过渡,主要发生受拉破坏,在冲击气压较大时沿轴向加载方向出现贯穿于砂岩试样的环状破坏面。

参考文献(References):

[1] 何满潮,谢和平,彭苏萍,等.深部开采岩体力学研究[J].岩石力学与工程学报,2005,24(16):2803-2813.

HE Manchao,XIE Heping,PENG Suping,et al.Study on rock mechanics in deep mining engineering[J].Chinese Journal of Rock Mechanics and Engineering,2005,24(16):2803-2813.

[2] 谢和平.深部岩体力学与开采理论研究进展[J].煤炭学报,2019,44(5):1283-1305.

XIE Heping.Research review of the state key research development program of China:Deep rock mechanics and mining theory[J].Journal of China Coal Society,2019,44(5):1283-1305.

[3] 李夕兵.岩石动力学基础与应用[M].北京:科学出版社,2014.

[4] 黄理兴.岩石动力学研究成就与趋势[J].岩土力学,2011,32(10):2889-2900.

HUANG Lixing.Development and new achievements of rock dynamics in China[J].Rock and Soil Mechanics,2011,32(10):2889-2900.

[5] LI Haibo,XIA Xiang,LI Jianchun,et al.Rock damage control in bedrock blasting excavation for a nuclear power plant[J].International Journal of Rock Mechanics and Mining Sciences,2011,48(2):210-218.

[6] RAMULU M,CHAKRABORTY A K,SITHARAM T G.Damage assessment of basaltic rock mass due to repeated blasting in a railway tunnelling project-A case study[J].Tunnelling and Underground Space Technology,2009,24(2):208-221.

[7] LUO Xuedong,JIANG Nan,WANG Mingyang,et al.Response of leptynite subjected to repeated impact loading[J].Rock Mechanics and Rock Engineering,2016,49(10):4137-4141.

[8] QI Chengzhi,WANG Mingyang,BAI Jiping,et al.Investigation into size and strain rate effects on the strength of rock-like materials[J].International Journal of Rock Mechanics and Mining Sciences,2016,86(1):132-140.

[9] BAGDE M N,PETROS V.Fatigue properties of intact sandstone samples subjected to dynamic uniaxial cyclical loading[J].International Journal of Rock Mechanics and Mining Sciences,2005,42(2):237-250.

[10] LI S H,ZHU W C,NIU L L,et al.Dynamic characteristics of green sandstone subjected to repetitive impact loading:Phenomena and mechanisms[J].Rock Mechanics and Rock Engineering,2018,51(6):1921-1936.

[11] BRAUNAGEL M J,GRIFFITH W A.The effect of dynamic stress cycling on the compressive strength of rocks[J].Geophysical Researeh Letters,2019,46(12):6479-6486.

[12] 朱晶晶,李夕兵,宫凤强,等.单轴循环冲击下岩石的动力学特性及其损伤模型研究[J].岩土工程学报,2013,35(3):531-539.

ZHU Jingjing,LI Xibing,GONG Fengqiang,et al.Dynamic characteristics and damage model for rock underuniaxial cyclic impact compressive loads[J].Chinese Journal of Geotechnical Engineering,2013,35(3):531-539.

[13] 李地元,孙小磊,周子龙,等.多次冲击荷载作用下花岗岩动态累计损伤特性[J].实验力学,2016,31(6):827-835.

LI Diyuan,SUN Xiaolei,ZHOU Zilong,et al.On the dynamic accumulated damage characteristics of granite subjected to repeated impact load action[J].Journal of Experimental Mechanics,2016,31(6):827-835.

[14] 张智宇,王建国,黄永辉,等.白云岩动力学特性及破坏模式试验研究[J].地下空间与工程学报,2018,14(5):1195-1201.

ZHANG Zhiyu,WANG Jianguo,HUANG Yonghui,et al.Experimental analysis on dynamic characteristics and failure mode of deep dolomite based on high-speed photography[J].Chinese Journal of Underground Space and Engineering,2018,14(5):1195-1201.

[15] 王志亮,杨辉,田诺成.单轴循环冲击下花岗岩力学特性与损伤演化机理[J].哈尔滨工业大学学报,2020,52(2):59-66.

WANG Zhiliang,YANG Hui,TIAN Nuocheng.Mechanical property and damage evolution mechanism of granite under uniaxial cyclic impact[J].Journal of Harbin Institute of Technology,2020,52(2):59-66.

[16] 田诺成,王志亮,熊峰,等.循环冲击荷载下轴压对花岗岩动力学特性的影响[J].哈尔滨工业大学学报,2021,53(5):156-164.

TIAN Nuocheng,WANG Zhiliang,XIONG Feng,et al.Influence of axial pressure on dynamic mechanical properties of granite under cyclic impact loading[J].Journal of Harbin Institute of Technology,2021,53(5):156-164.

[17] 金解放,李夕兵,王观石,等.循环冲击载荷作用下砂岩破坏模式及其机理[J].中南大学学报(自然科学版),2012,43(4):1453-1461.

JIN Jiefang,LI Xibing,WANG Guanshi,et al.Failure modes and mechanisms of sandstone under cyclic impact loadings[J].Journal of Central South University (Science and Technology),2012,43(4):1453-1461.

[18] 金解放,李夕兵,常军然,等.循环冲击作用下岩石应力应变曲线及应力波特性[J].爆炸与冲击,2013,33(6):613-619.

JIN Jiefang,LI Xibing,CHANG Junran,et al.Stress-strain curve and stress wave characteristics of rock subjected to cyclic impact loadings[J].Explosion and Shock Waves,2013,33(6):613-619.

[19] 金解放,李夕兵,殷志强,等.轴压和循环冲击次数对砂岩动态力学特性的影响[J].煤炭学报,2012,37(6):923-930.

JIN Jiefang,LI Xibing,YIN Zhiqiang,et al.Effects of axial pressure and number of cyclic impacts on dynamic mechanical characteristics of sandstone[J].Journal of China Coal Society,2012,37(6):923-930.

[20] 于洋,徐倩,刁心宏,等.循环冲击对围压作用下砂岩特征的影响[J].华中科技大学学报(自然科学版),2019,47(6):127-132.

YU Yang,XU Qian,DIAO Xinhong,et al.Effect of cyclic impact on sandstone characteristics under confining pressures[J].Journal of Huazhong University of Science and Technology (Natural Science Edition),2019,47(6):127-132.

[21] 唐礼忠,王春,程露萍,等.一维静载及循环冲击共同作用下矽卡岩力学特性试验研究[J].中南大学学报(自然科学版),2015,46(10):3898-3906.

TANG Lizhong,WANG Chun,CHENG Luping,et al.Experimental study of mechanical characteristics of skarn under one-dimensional coupled static and cyclic impact loads[J].Journal of Central South University (Science and Technology),2015,46(10):3898-3906.

[22] 许金余,吕晓聪,张军,等.围压条件下岩石循环冲击损伤的能量特性研究[J].岩石力学与工程学报,2010,29(S2):4159-4165.

XU Jinyu,LÜ Xiaocong,ZHANG Jun,et al.Research on energy properties of rock cyclical impact damage under confining pressure[J].Chinese Journal of Rock Mechanics and Engineering,2010,29(S2):4159-4165.

[23] 李晓锋,李海波,刘凯,等.冲击荷载作用下岩石动态力学特性及破裂特征研究[J].岩石力学与工程学报,2017,36(10):2393-2405.

LI Xiaofeng,LI Haibo,LIU Kai,et al.Dynamic properties and fracture characteristics of rocks subject to impact loading[J].Chinese Journal of Rock Mechanics and Engineering,2017,36(10):2393-2405.

[24] ZHOU Y X,XIA K,LI X B,et al.Suggested methods for determining the dynamic strength parameters and mode-I fracture toughness of rock materials[J].International Journal of Rock Mechanics and Mining Sciences,2012,49:105-112.

[25] WANG P,YIN T,LI X,et al.Dynamic properties of thermally treated granite subjected to cyclic impact loading[J].Rock Mechanics and Rock Engineering,2019,52(4):991-1010.

[26] 宋力,胡时胜.SHPB数据处理中的二波法与三波法[J].爆炸与冲击,2005,25(4):368-373.

SONG Li,HU Shisheng.Two-wave and three-wave method in SHPB data processing[J].Explosion and Shock Waves,2005,25(4):368-373.

[27] SHAN Renliang,SONG Yongwei,SONG Liwei,et al.Dynamic property tests of frozen red sandstone using a split hopkinson pressure bar[J].Earthquake Engineering and Engineering Vibration,2019,18(3):511-519.

[28] 王梦想,汪海波,宗琦.冲击荷载作用下煤矿泥岩能量耗散试验研究[J].煤炭学报,2019,44(6):1716-1725.

WANG Mengxiang,WANG Haibo,ZONG Qi.Experimental study on energy dissipation of mudstone in coal mine under impact loading[J].Journal of China Coal Society,2019,44(6):1716-1725.

[29] 王梦想,汪海波,宗琦.煤矿泥岩冲击动态力学特性与破裂破碎特征分析[J].振动与冲击,2019,38(4):137-143.

WANG Mengxiang,WANG Haibo,ZONG Qi.Analysis dynamic mechanical characteristics and fracture breaking characteristics of coal mine mudstone[J].Journal of Vibration and Shock,2019,38(4):137-143.

[30] 宫凤强,李夕兵,刘希灵.三维动静组合加载下岩石力学特性试验初探[J].岩石力学与工程学报,2011,30(6):1179-1190.

GONG Fengqiang,LI Xibing,LIU Xiling.Preliminary experimental study of characteristics of rock subjected to 3D coupled static and dynamic loads[J].Chinese Journal of Rock Mechanics and Engineering,2011,30(6):1179-1190.

Strain rate effect and energy dissipation characteristics of sandstone in coal measures under impact loading

YU Yongqiang1,2,ZHANG Wenlong1,FAN Lidan1,2,GONG Jian1,YANG Meng1,SUN Liang1

(1.School of Civil Engineering,Henan Polytechnic University,Jiaozuo 454003,China; 2.International Joint Research Laboratory of Henan Province for Underground Space Development and Disaster Prevention,Jiaozuo 454003,China)

Abstract:The deformation of surrounding rock will increase under dynamic load during the excavation of coal mine roadway.In order to study the mechanical properties and energy dissipation of sandstone in coal measures under impact loading,the sandstone in the surrounding rock of Chensilou coal mine roadway was collected as the research sample,and a split Hopkinson pressure bar test device,with a diameter of 50 mm,was used in uniaxial single impact compression test and cyclic impact compression test.The strain rate effect,energy dissipation characteristics and failure modes of the sandstone were analyzed.Under the uniaxial single impact loading,with the increase of average strain rate,peak stress and peak strain,the plasticity and strength of the sandstone samples increase,while the secant modulus decreases gradually.Moreover,the peak strain has a linear increase relationship with the average strain rate,and the peak stress approximately increases with the 1/3 power of the average strain rate.With the increase of the average strain rate,the absorption energy per unit volume of the sandstone samples shows a linear increase tendency,and the degree of fragmentation of the samples continues to increase.Under the impact of constant compressive stress,the internal cracks of the sandstone samples keep crossing and expanding,and the splitting failure occurs along the axial direction.Under the cyclic impact loading,the average strain rate and peak strain of sandstone samples gradually develop with the increase of the cyclic loading times,while the peak stress,secant modulus and second type secant modulus gradually decreases.During the cyclic impact under the fixed pressure,with the increase of the impact times,the incident energy basically remains unchanged,and the reflection energy and absorption energy gradually increase,whereas the transmission energy gradually decreases.The absorption energy per unit volume of the sandstone samples increases slowly firstly and then suddenly in-creases with the impact times.Under the cyclic impact loading,the energy absorption capacity of the sandstone samples is enhanced.With the increase of impact times,the distribution of microcracks inside the samples enlarges and the expansion speed accelerates.As a result,the tensile failure mainly occurs,and the annular failure plane running through the sandstone samples generates along the axial loading direction.

Key words:coal measures sandstone;SHPB;energy dissipation;strain rate;failure mode

中图分类号:TD315

文献标志码:A

文章编号:0253-9993(2021)07-2281-13

移动阅读

收稿日期:2020-03-19

修回日期:2020-05-24

责任编辑:郭晓炜

DOI:10.13225/j.cnki.jccs.2020.0440

基金项目:国家自然科学基金资助项目(51874119)

作者简介:余永强(1972—),男,河南商城人,教授。Tel:0391-3987635,E-mail:yyq@hpu.edu.cn

通讯作者:范利丹(1982—),女,河南濮阳人,副教授。Tel:0391-3987639,E-mail:lidanfan@hpu.edu.cn

引用格式:余永强,张文龙,范利丹,等.冲击荷载下煤系砂岩应变率效应及能量耗散特征[J].煤炭学报,2021,46(7):2281-2293.

YU Yongqiang,ZHANG Wenlong,FAN Lidan,et al.Strain rate effect and energy dissipation characteristics of sandstone in coal measures under impact loading[J].Journal of China Coal Society,2021,46(7):2281-2293.