基于加载速率效应的煤样损伤演化研究

    李庆文, 高森林, 黄筱

    李庆文,高森林,黄筱. 基于加载速率效应的煤样损伤演化研究[J]. 煤矿安全,2023,54(11):105−115. DOI: 10.13347/j.cnki.mkaq.2023.11.001
    引用本文: 李庆文,高森林,黄筱. 基于加载速率效应的煤样损伤演化研究[J]. 煤矿安全,2023,54(11):105−115. DOI: 10.13347/j.cnki.mkaq.2023.11.001
    LI Qingwen, GAO Senlin, HUANG Xiao. Study on damage evolution of coal samples based on loading rate effect[J]. Safety in Coal Mines, 2023, 54(11): 105−115. DOI: 10.13347/j.cnki.mkaq.2023.11.001
    Citation: LI Qingwen, GAO Senlin, HUANG Xiao. Study on damage evolution of coal samples based on loading rate effect[J]. Safety in Coal Mines, 2023, 54(11): 105−115. DOI: 10.13347/j.cnki.mkaq.2023.11.001

    基于加载速率效应的煤样损伤演化研究

    基金项目: 国家自然科学基金资助项目(51704144);辽宁省博士科研启动基金资助项目(2019-BS-120);辽宁省自然科学基金指导资助项目(20180550297)
    详细信息
      作者简介:

      李庆文(1987—),男,辽宁绥中人,讲师,硕士研究生导师,博士,主要从事岩石力学、新材料与新型组合结构及离散元-有限元耦合细观模拟等方面的教学及科研工作。E-mail: lgjzlqw@163.com

    • 中图分类号: TD313

    Study on damage evolution of coal samples based on loading rate effect

    • 摘要:

      为了明确加载速率对煤样损伤行为的影响,利用SAM-2000岩石力学试验系统进行不同加载速率下单轴压缩试验,分析了破坏形态、应力—应变、弹性模量、峰值应力和峰值应变等随加载速率的变化规律;推导了能量损伤本构模型,并验证本构模型的合理性。研究结果表明:煤样在压缩破坏过程中可分为压密、弹性、屈服和破坏4个阶段,煤样呈现剪切破坏;随着加载速率的变化,峰值强度和加载速率呈正相关;弹性模量初始呈现近70°的线性增长,后缓慢增加;峰值应变随着加载速率的增加先呈现出减小而后增大趋势;煤样的弹性能、耗散能随着加载速率的增大表现出线性增长的关系,比吸能随加载速率增大而增大,峰值应力前主要以弹性能为主,耗散能为辅,破坏阶段反之;加载速率不会影响煤样能量转化规律趋势和损伤破坏过程。

      Abstract:

      In order to clarify the influence of loading rate on the damage behavior of coal samples, the uniaxial compression tests of heterogeneous coal samples under different loading rates were carried out based on SAM-2000 rock mechanics test system, and the variation of the failure morphology, stress-strain, elastic modulus, peak stress and peak strain with loading rates were analyzed; the energy damage constitutive model was derived and the rationality of the constitutive model was verified. The study result shows that coal samples can be divided into four stages during the compressive failure process: compaction stage, elastic stage, plastic stage and failure stage; with the change of loading rate, peak strength is positively correlated with loading rate. Elastic modulus firstly increases approximately 70° linearly, and then gradually decreases slowly. With the increase of loading rate, the peak strain of coal pillars first decreases and then increases; the elastic energy and dissipation energy of coal samples show a linear increase with the increase of loading rate. And the specific energy absorption increases with the increase of loading rate. Before the peak stress, the elastic energy is mainly dominated, and the dissipation energy is supplemented, and the failure stage is the opposite. The trend of energy conversion and the damage process of coal samples will not be influenced by loading rates.

    • 煤本质上是非均质有机岩石,属于1种矿物质集合体,天然状态下其内部结构呈现出较多孔洞、裂隙等[1]。在深地资源开采过程中,由于地质环境复杂多变等,大量煤柱会作为安全承载结构留于矿区。受到采场工作面推进速度的影响,如爆破、掘进等,留设的煤柱会受到不同加载速率作用影响[2-4]。矿震、岩爆等动力灾害已成为深地煤矿开采领域亟须解决的关键问题[5];因此,国内外众多学者研究不同加载速率对煤岩力学行为、能量存储和耗散规律方面的影响。夏冬等[6]进行了煤样在不同应变率下单轴压缩破坏试验,结果表明煤样破坏机制与应变率相关;HUANG等[7]进行了不同加载速率下的煤岩组成的单轴压缩试验,得出了轴向应变率对煤和岩石力学行为的影响规律;王晓东等[8]基于室内单轴压缩试验,得出不同加载速率下花岗岩力学行为和能量变化规律;王业平等[9]研究了不同加载速率下含瓦斯突出煤体损伤特征,得出了不同速率下煤体损伤演化规律;LU等[10]基于单轴压缩条件和声发射,发现随着加载速率的增加,试件力学行为对非均匀性的敏感性逐渐降低;ZHU等[11]制备长方体煤样,利用超声波探测仪等手段,研究了不同加载速率对试件渐进损伤破坏特性的影响;李海涛等[12-13]开展了不同加载速率下煤样单轴压缩试验,深入分析了煤样的力学行为;马振乾等[14]探究了5种不同加载速率对煤样侧向变形的影响;LI等[15]、GONG[16]开展不同应变率和围压下的单轴压缩试验,探究了煤样变形和力学性能;ZHANG等[17]利用岩石三轴仪对煤样进行了不同应变率和不同承压卸载速率下的三轴载荷破坏实验,基于特征参数构建了损伤模型;王凯等[18]开展了不同含水率下煤样单轴压缩试验,推导出煤样含水率分段式损伤本构模型。综上,关于不同加载速率下煤样单轴压缩试验研究,主要集中在力学性能、破坏形态和能量等方面,而不同加载速率条件下煤样的损伤演化本构模型规律却鲜有研究。为此,通过开展不同加载速率下单轴压缩试验,比较不同加载速率下的破坏形态和力学性能,深入分析能量转化过程,推导出不同加载速率下分段式损伤本构模型,并利用数据进行验证。

      煤样试件使用SAM-2000型微机控制电液伺服岩石三轴试验机进行单轴压缩试验,试验装置主要是由电子自动液压伺服系统、外荷载加载系统、全数字监控系统和应力-应变采集系统组成。单轴压缩试验机轴向荷载最高可达2000 kN,其中轴向外荷载由试验机自动加载并记录数据,应变则是由粘贴在煤样上的应变片采集数据而得,荷载和应变应同时同步存储于计算机上。

      试验煤样选自内蒙古自治区鄂尔多斯市塔然高勒矿区红庆梁煤矿某工作面煤柱,根据GB/T 50266—2013《工程岩体试验方法标准》[19]现场原位取样与标准煤样制备方法,通过剪、切、磨等方法制成50 mm×100 mm(直径×高)的标准圆柱体煤样。试件平均含水率约为8%,煤样加工精准度应满足标准单轴压缩试验要求。为降低煤样试件离散性对试验结果的影响,对所有煤样试件进行超声波试验检测,从中择取离散性小的试件,以减小试验误差。

      试验择选出12个尺寸为ϕ50 mm×100 mm的煤样圆柱试件,考虑不同加载速率(0.01、0.1、1.0、10.0 mm/min)分成4组,每组3个煤样试件。试验基于SAM-2000单轴压缩试验,加载方式采用位移控制,加载至煤样完全破坏。煤样主要试验结果见表1

      表  1  煤样主要试验结果
      Table  1.  Main test results of coal samples
      组别试件
      直径/
      mm
      高/
      mm
      密度/
      (g·cm−3
      波速/
      (km·s−1
      弹性模量/MPa抗压强度/MPa
      弹性
      模量
      平均值峰值
      强度
      平均值
      1 A-1 49.60 100.19 1.29 1.282 4 632.57 4 776.89 15.79 19.04
      A-2 49.70 100.32 1.30 1.292 3 567.16 16.14
      A-3 49.70 100.48 1.32 1.773 6 130.94 25.20
      2 B-1 49.59 100.37 1.32 1.377 3 599.95 4 653.94 19.03 24.28
      B-2 49.84 100.44 1.31 1.401 5 105.76 27.52
      B-3 49.68 100.28 1.32 1.403 5 256.10 26.29
      3 C-1 49.55 100.54 1.30 1.512 4 849.18 5 712.25 20.95 25.63
      C-2 49.63 100.39 1.30 1.527 5 932.97 23.99
      C-3 49.61 100.24 1.33 1.538 6 354.59 31.94
      4 D-1 49.42 100.36 1.34 1.596 6 961.22 6 536.69 30.81 32.56
      D-2 49.86 100.38 1.31 1.612 6 585.36 34.81
      D-3 49.55 100.39 1.33 1.614 6 063.48 32.07
      下载: 导出CSV 
      | 显示表格

      为便于后续讨论,准确区分每个试件,根据试件参数不同对试件进行命名,其命名方式则是由字母与数字构成。第1个字母表示加载速率,即A、B、C、D分别代表0.01、0.1、1.0、10.0 mm/min的加载速率;其后,连字符后的数字则用来区分同组的不同试件。以“C-3”为例,C-3表示加载速率为1.0 mm/min、煤样圆柱3号试件。

      煤样破坏形态随不同加载速率变化,而其力学特性也随不同加载速率而改变。在试验过程中,对所有试件进行拍照和录像记录,自动测量并记录轴应力和轴应变,并对试验结果进行分析。煤样强度和变形行为同岩石变形性质相差无几,与煤样内部微裂纹的闭合、萌生、扩展和贯通有关[20]。随着加载速率的增大,峰值应力与峰值处应变也呈现增大趋势。

      不同加载速率下煤样单轴压缩应力-应变曲线如图1。通过图1可将煤样压缩变形过程大致分为4个阶段,即压密、弹性、屈服和破坏阶段。

      图  1  不同加载速率下煤样单轴压缩应力-应变曲线
      Figure  1.  Stress-strain curves of coal samples under uniaxial compression at different loading rates

      图1可以看出:当煤样处于压密阶段时,煤样中的原生裂纹和微裂隙逐渐被压密,应力-应变曲线向下凹,即应力-应变曲线斜率随应力增大而增大,表明原生裂纹和微裂隙的闭合开始较快,而后逐渐减慢;当煤样处于弹性阶段时,应力-应变呈现出线性关系,在较大程度上煤样表现为可恢复的变形,在持续的外荷载下,煤样逐渐开始产生新的次生裂纹与微裂隙,且与原生裂隙逐渐贯通融合,在主裂隙带上扩展,继续产生新微裂纹,形成裂纹局部化;当煤样处于屈服阶段时,煤样破坏将发生本质变化,随着应力集中效应,首先煤样内部较为薄弱的位置开始破坏,直到煤样完全破坏,煤样从体积压缩转化为扩容,其峰值点称为抗压强度;当煤样处于破坏阶段时,煤样达到抗压强度后内部结构已完全破坏,从表观看煤样可大致保持原态,应力-应变曲线开始快速下降,但轴向应力未到0,说明煤样仍具有一定强度。

      单轴压缩代表性试件破坏形态如图2

      图  2  代表性试件破坏模式
      Figure  2.  Representative failure modes of specimens

      图2可看出:煤样破坏形态可分为2种:①第1种:煤样在受到外荷载作用且加载速率较小时,煤样表面出现微裂缝,随着荷载的增大,微裂缝逐步扩展,呈现剪切破坏,裂缝分为斜裂缝和直裂缝,产生少量碎屑沿裂缝迸出,导致煤样破坏严重,不能维持其原有形态;②第2种:煤样在受到外荷载作用且加载速率较大时,微裂缝逐步扩展贯穿整个煤样,裂缝形态呈“树叶”状,产生小块碎屑崩出,以破碎状为主,表明加载速率严重影响煤样的损伤过程,当加载速率为0.01 mm/min时,煤样内部原生裂纹和裂隙有足够时间去扩展和贯通,反之则不然。

      结果表明,煤样长期处于地应力环境中,使煤样破坏变形变得缓慢,逐渐变形直至破坏。相反,煤样在高压冲击荷载作用下,煤样径向扩展变形大,直至破坏。

      随着加载速率的增大,煤样轴向应力逐渐增大。从应力-应变关系中可以观察出,加载速率越大,弹性模量也越大,煤样破坏时间越短[21]

      弹性模量F、峰值应力σρ和峰值应变εP与加载速率v的关系如图3

      图  3  弹性模量、峰值应力和峰值应变与加载速率的关系
      Figure  3.  The relationship between elastic modulus, peak stress and peak strain and loading rate

      当加载速率为0.01 mm/min,煤样的平均峰值强度为19.04 MPa;当加载速率为0.1 mm/min,煤样的平均峰值强度为24.28 MPa;当加载速率为1.0 mm/min,煤样的平均峰值强度为25.63 MPa;当加载速率为10.0 mm/min,煤样的平均峰值强度为32.56 MPa。

      图3(a)可知:随着加载速率的增加,峰值强度和加载速率呈正相关。

      图3(b)可知:加载速率为0.01 mm/min时的弹性模量微大于加载速率为0.1 mm/min时的弹性模量,随着加载速率的增大,弹性模量初始呈现近70°的线性增长,后缓慢增加。

      图3(c)可知:当加载速率较小时,峰值应变随加载速率的增大而陡然增大,随着加载速率的增大,峰值应变随着加载速率的增加先呈现出减小而后增大趋势。

      图3可知:弹性模量、峰值应力和峰值应变的试验值与拟合曲线基本吻合,R2=0.99。

      在外荷载下岩石加载变形破坏过程中,本质上是能量输入、积聚、耗散和释放的过程。假设不考虑温度的影响,外荷载对岩石做功一部分以弹性能的形式产生,一部分则以塑性变形能和损伤能等耗散能的形式产生;当存储弹性能达到极限平衡状态时,岩石发生变形破坏,存储的一部分能量以破坏损伤形式表现出来,其余的能量以热能、动能以及各种辐射能等形式表现出来。其中研究煤样在不同加载速率下变形破坏过程的弹性能和耗散能的演化规律,从能量角度揭示煤柱的损伤机制,对于认识煤样的损伤演化具有重要意义[22]

      根据能量守恒定律(热力学第一定律)[23],可得到如下关系:

      $$ U=U^{\rm{d}}+U^{\rm{e}}+U^0 $$ (1)

      式中:U为总能量;Ud为煤样耗散能,主要用于塑性变形和内部损失变形;Ue为弹性应变能,主要是存储在煤样内部的能量;U0为以热辐射、热交换等方式释放的能量。

      热辐射、热交换产生的能量很小,可以忽略不计,因此可得[24]

      $$ U=U^{\rm{d}}+U^{\rm{e}} $$ (2)

      煤样单元能量[25]为:

      $$ U=\int_{0}^{\varepsilon_{1}} \sigma_{1} {\rm{d}} \varepsilon_{1}+\int_{0}^{\varepsilon_{2}} \sigma_{2} {\rm{d}} \varepsilon_{2}+\int_{0}^{\varepsilon_{3}} \sigma_{3} {\rm{d}} \varepsilon_{3} $$ (3)
      $$ U^{\rm{e}}=0.5 \sigma_1 \varepsilon_1^{\rm{e}}+0.5 \sigma_2 \varepsilon_2^{\rm{e}}+0.5 \sigma_3 \varepsilon_3^{\rm{e}} $$ (4)

      式中:σ1σ2σ3分别为第一主应力、第二主应力、第三主应力;ε1ε2ε3分别为σ1σ2σ3对应的轴向应变。

      由胡克定律,将式(4)改写为:

      $$ U^{{\rm{e}}}=\frac{1}{2 E_{0}}\left[\sigma_{1}^{2}+\sigma_{2}^{2}+\sigma_{3}^{2}-2 \mu\left(\sigma_{1} \sigma_{2}+\sigma_{2} \sigma_{3}+\sigma_{1} \sigma_{2}\right)\right] $$ (5)

      式中:E0为初始弹性模量;μ为泊松比。

      对于单轴压缩而言,σ2=0,σ3=0,弹性能计算公式可简化为:

      $$ U^{\rm{e}}=\frac{1}{2 E_0} \sigma_1^2 $$ (6)

      根据定积分计算面积的概念,总能量U计算公式为:

      $$ U=\sum_{i=1}^n \frac{1}{2}\left(\sigma_1^i+\sigma_1^{i+1}\right)\left(\varepsilon_1^{i+1}-\varepsilon_1^i\right) $$ (7)

      式中:$ \sigma_{1}^{i} $为应力−应变曲线中每一点应力值;$ \varepsilon_{1}^{i} $为应力-应变曲线中每一点应变值。

      因此,煤样试验过程中的耗散能为:

      $$ U^{{\rm{d}}}=U-U^{{\rm{e}}} $$ (8)

      比吸能是试件单位质量所吸收的能量,即在外荷载压缩过程中所吸收的总能量U与试件质量m的比值,其式为[26]

      $$ E_{{\rm{s}}}=U / m $$ (9)

      式中:Es为比吸能;U为试件所吸收的总能量;m为试件质量。

      不同加载速率下煤样的比吸能见表2,不同加载速率与比吸能关系如图4

      表  2  不同加载速率下煤样的比吸能
      Table  2.  Specific energy absorption of coal samples at different loading rates
      试件号质量/g总能量/(J·cm−3比吸能/(10−3 J·cm−3·g−1
      A-2251.60.078 80.310
      B-2256.10.126 60.490
      C-3256.30.106 70.416
      D-2255.40.097190.381
      下载: 导出CSV 
      | 显示表格
      图  4  不同加载速率与比吸能关系
      Figure  4.  Relationship between different loading rates and specific energy absorption

      当煤样处于0.01 ~10 mm/min 范围时,将0.03 mm/min视为界限加载速率,即在同质量的煤样中速度达到界限加载速率时,煤样吸收的能量达到峰值。基于煤本身的非均质性,天然裂纹性,界限加载速率可用于反映在煤矿施工过程中巷道围岩加卸载速率,例如巷道开挖速率、采煤工作面的推进速度等。由表2可知:随着加载速率的增大,煤样吸收的总能量越多,导致其比吸能逐渐增大,而试件C-3、D-2因煤样本身离散性较大,呈现变小的结果。

      根据能量计算,可得到煤样在不同加载速率下应力-应变过程中的能量演化特征曲线,不同加载速下率煤样能量演化如图5

      图  5  不同加载速下率煤样能量演化
      Figure  5.  Energy evolution of coal samples under different loading rates

      随着加载速率增大煤样总能量先激增再缓降;其中当速率处于0.01~0.1 mm/min范围时,煤样中的总能量激增,约提高1.4倍;当处于0.1~1.0 mm/min范围时,此过程是1个过渡区,煤样最大总能量位于此过程,约是前1个过程的1.2倍;当处于1.0~10.0 mm/min范围时,煤样峰值总能量缓慢下降。

      随着加载速率增大煤样的弹性能水平缓慢增加后猛增,接着下降;当煤样处于0.01~0.1 mm/min范围时,初始弹性能猛增,约提高了2倍;处于0.1~1.0 mm/min范围时,此过渡区约提高3.7%;处于1.0~10.0 mm/min范围时,弹性能缓慢增加,约提高4.2%。

      总能量U和耗散能Ud随轴应变的增大而增大。弹性应变能Ue曲线与应力-应变曲线吻合较好,在峰值应力之前,弹性应变能Ue逐渐增大;在峰值应力后,弹性应变能Ue逐渐变小。

      不同加载速率下的能量演化特征曲线中,可将煤样试验过程分为初始压密、弹性、屈服和破坏阶段。如图5:在压密阶段,总能量U、弹性应变能Ue和耗散能Ud均增大,此时弹性应变能Ue小于耗散能Ud,总能量U转化为耗散能Ud的形式耗散和弹性应变能Ue的形式存储起来,天然煤样自身具有较多裂隙,在外界能量作用下,使内部部分裂隙闭合,颗粒间摩擦滑动消耗部分能量;在弹性阶段,煤样内部原生裂纹已压密,缓慢出现新的次生裂纹,煤样从外界吸收的能量表现在弹性应变能Ue形式存储于煤样中,其内部存储的弹性应变能Ue随弹性变形而增加,煤样吸收弹性应变能Ue的速率恒定,因此弹性应变能呈斜直线;在屈服阶段,煤样内部的次生裂纹以较快速度萌生、扩展和贯通形成新生裂纹;接近峰值应力时,耗散能Ud逐渐开始增加,弹性应变能Ue增长速率开始放缓并趋于平缓,表明煤样内部结构发生了大的变化;在峰值应力后,耗散能Ud迅速增大至最大值,弹性应变能Ue快速下降至最小值。

      总体上看,煤样的损伤破坏是不断吸收和释放能量的过程。在不同加载速率下,煤样能量演化特征曲线表现出相同的规律,表明加载速率不会影响煤样能量转化规律趋势和损伤破坏过程。

      煤样损伤演化实质上就是能量的吸收与释放,将煤样单元能量损伤值D定义为[27]任意时刻耗散能与耗散能最大值的比值,且微元体的强度损伤特征符合统计学规律和Weibull分布规律。

      $$ D=\frac{U^{{\rm{d}}}}{U_{\max }^{{\rm{d}}}} $$ (10)

      式中:$ U^{{\rm{d}}} $为煤样每一时刻所对应产生的耗散能; $ U_{\text {max }}^{{\rm{d}}} $为最大耗散能。

      采用不同加载速率下代表性试件试验数据绘制出损伤值D与应力-应变曲线关系图,能量损伤拟合曲线如图6

      图  6  能量损伤拟合曲线
      Figure  6.  Energy damage fitting curves

      图6可知:损伤值D图样形式走向大概一致,在不同加载速率下,煤样都有1个较小的损伤初始值,这是由于煤样本身性质及外界因素(开采扰动等)所造成;初始阶段损伤值D都较小,主要是弹性应变能吸收存储能量的1个过程,试件还未开始破坏,耗散能Ud还未发挥作用;随着时间推移,达到损伤突变点,煤样试件开始破坏,能量由弹性应变能Ue逐渐转化成耗散能Ud,且耗散能Ud越来越大,直到停止,煤样完全破坏;损伤突变点位于应力-应变峰值点处左右。

      图6可知:当加载速率小于0.01 mm/min时,加载速率对煤样力学特性影响甚微,以初始损伤为主;当加载速率大于0.01 mm/min时,损伤值D开始逐渐增大,因为初始损伤加快微裂纹的发展、扩展和贯通。表明不同加载速率对于损伤煤样强度影响显著[28]

      根据煤样能量损伤曲线可拟合出轴应变ε1与损伤值D的关系曲线,在不同加载速率下损伤值D相关性很好,煤样损伤值D拟合均满足Weibull分布,表达式为:

      $$ D=1-{\rm{e}}^{\left[-\left(\frac{\varepsilon_1}{D_0}\right)^n\right]} $$ (11)

      式中:ε1为轴应变;D0为初始损伤值;n为不同加载速率下指数。

      不同加载速率下参数见表3。参数与加载速率拟合曲线如图7

      表  3  不同加载速率下参数
      Table  3.  Parameters at different loading rates
      加载速率/(mm·min−1)Don
      0.010.007 706.823 46
      0.100.008 539.515 30
      1.000.006 3233.347 92
      10.000.005 8728.195 12
      下载: 导出CSV 
      | 显示表格
      图  7  参数与加载速率拟合曲线
      Figure  7.  Fitting curves of parameters and loading rate

      根据图7,得到参数Don随加载速率的变化关系拟合关系式,其中Do初始损伤值是由煤样本身性质决定,n随着加载速率的增大而逐渐缓慢增加;随着加载速率的增大,前期损伤值D变化较小,后期阶段损伤值D激增,煤样能量耗散极具增大,达到最大耗散能,煤样破坏,能量转化过程结束。

      对于损伤演化计算方法一般根据岩石某一点强度出发,进一步推出岩石损伤本构模型方程。假设煤样在压密阶段不发生损伤破坏,以压密阶段最后1个点为拐点(εAσA),采用分段式损伤本构模型计算出的理论模型曲线与试验值进行对比[29]

      根据鲁祖德[30]对岩石裂隙理论分析的结果表明,煤样在压密阶段应力-应变关系式经修正为:

      $$ \sigma=\sigma_{{\rm{A}}}\left(\varepsilon / \varepsilon_{{\rm{A}}}\right)^{1.5} $$ (12)

      式中:σA为压密阶段最大应力;εA为压密阶段最大应变;σ为某一时刻应力;ε为某一时刻应变。

      假设煤样在弹性、屈服和破坏阶段发生不间断损伤变形。根据Kachanov[31]对损伤演化值的定义研究,煤样损伤演化曲线和方程,及相关研究结果[32-33]煤样在单轴压缩下后阶段损伤演化本构模型为:

      $$ \sigma=\sigma_{{\rm{A}}}+E\left(\varepsilon-\varepsilon_{{\rm{A}}}\right)(1-D) $$ (13)

      联合式(12)和式(13)可得到不同加载速率下基于可释放的弹性应变能和耗散能的损伤演化本构模型为:

      $$ \sigma= \begin{cases}\sigma_{\rm{A}}\left(\varepsilon / \varepsilon_{\rm{A}}\right)^{1.5} & \left(\varepsilon \leqslant \varepsilon_{\rm{A}}\right) \\ \varepsilon_{\rm{A}}+E\left(\varepsilon-\varepsilon_{\rm{A}}\right)(1-D) & \left(\varepsilon>\varepsilon_A\right)\end{cases} $$ (14)

      根据单轴压缩试验数据对损伤演化本构模型进行验证。将参数代入式(13),可绘制出不同加载速率下煤样分段式损伤演化本构模型的理论模型曲线(图略)。从试验曲线拟合度看,损伤演化本构模型方程可以很好地拟合煤样在不同加载速率下的应力-应变曲线,并且还可以从拟合的应力-应变曲线中清晰分析出压密、弹性、屈服和破坏阶段。总之,损伤演化本构模型理论曲线与试验值曲线基本吻合。

      采用分段式损伤演化本构方程将优于连续式损伤演化本构方程,因为连续式损伤演化本构方程在压密阶段会产生很大误差,导致在不同加载速率下试验曲线与理论模型曲线拟合度变小。分段式损伤演化本构模型临界点(εAσA)是压密阶段与弹性阶段的分界点,从弹性阶段开始考虑煤样开始损伤,增加了峰值点前的拟合度,导致分段式损伤演化本构模型更加适合用于考虑分析不同加载速率下煤样的破坏过程和能量转化问题。

      1)当加载速率较小时(0.01、0.1 mm/min),煤样主要处于压缩,当达到峰值应力时,煤样体积开始在环向膨胀,少量碎屑飞出,煤样呈剪切破坏,裂缝为直裂缝和斜裂缝;当加载速率较大时(1.0、10.0 mm/min),煤样压缩时间较短,少量小块迸出,裂缝呈“树叶”状;随着加载速率增大,弹性模量逐渐增大,表明加载速率对煤样破坏过程中的应力-应变曲线有很大影响。

      2)在不同加载速率下的煤样能量损伤演化特征曲线规律在压密、弹性、屈服和破坏阶段大致相同,表明加载速率不会影响煤样能量转化规律趋势和损伤破坏过程;但随着轴应变增加其煤样吸收的能量也逐渐增加。

      3)根据耗散能定义损伤值D,对能量损伤曲线进行拟合,得出在不同加载速率下的损伤值D拟合结果均满足Weibull分布,且拟合度较高。建立参数Don与加载速率的关系,通过函数进行拟合。Do初始损伤值是由煤样本身性质决定,n随着加载速率的增大而逐渐缓慢增加。

      4)推导出在不同加载速率下煤样破坏的分段式损伤本构模型,通过分段式本构模型进一步提高峰值点前的拟合精度;基于弹性模量与不同加载速率的拟合函数,得到不同加载速率对本构模型的影响。

    • 图  1   不同加载速率下煤样单轴压缩应力-应变曲线

      Figure  1.   Stress-strain curves of coal samples under uniaxial compression at different loading rates

      图  2   代表性试件破坏模式

      Figure  2.   Representative failure modes of specimens

      图  3   弹性模量、峰值应力和峰值应变与加载速率的关系

      Figure  3.   The relationship between elastic modulus, peak stress and peak strain and loading rate

      图  4   不同加载速率与比吸能关系

      Figure  4.   Relationship between different loading rates and specific energy absorption

      图  5   不同加载速下率煤样能量演化

      Figure  5.   Energy evolution of coal samples under different loading rates

      图  6   能量损伤拟合曲线

      Figure  6.   Energy damage fitting curves

      图  7   参数与加载速率拟合曲线

      Figure  7.   Fitting curves of parameters and loading rate

      表  1   煤样主要试验结果

      Table  1   Main test results of coal samples

      组别试件
      直径/
      mm
      高/
      mm
      密度/
      (g·cm−3
      波速/
      (km·s−1
      弹性模量/MPa抗压强度/MPa
      弹性
      模量
      平均值峰值
      强度
      平均值
      1 A-1 49.60 100.19 1.29 1.282 4 632.57 4 776.89 15.79 19.04
      A-2 49.70 100.32 1.30 1.292 3 567.16 16.14
      A-3 49.70 100.48 1.32 1.773 6 130.94 25.20
      2 B-1 49.59 100.37 1.32 1.377 3 599.95 4 653.94 19.03 24.28
      B-2 49.84 100.44 1.31 1.401 5 105.76 27.52
      B-3 49.68 100.28 1.32 1.403 5 256.10 26.29
      3 C-1 49.55 100.54 1.30 1.512 4 849.18 5 712.25 20.95 25.63
      C-2 49.63 100.39 1.30 1.527 5 932.97 23.99
      C-3 49.61 100.24 1.33 1.538 6 354.59 31.94
      4 D-1 49.42 100.36 1.34 1.596 6 961.22 6 536.69 30.81 32.56
      D-2 49.86 100.38 1.31 1.612 6 585.36 34.81
      D-3 49.55 100.39 1.33 1.614 6 063.48 32.07
      下载: 导出CSV

      表  2   不同加载速率下煤样的比吸能

      Table  2   Specific energy absorption of coal samples at different loading rates

      试件号质量/g总能量/(J·cm−3比吸能/(10−3 J·cm−3·g−1
      A-2251.60.078 80.310
      B-2256.10.126 60.490
      C-3256.30.106 70.416
      D-2255.40.097190.381
      下载: 导出CSV

      表  3   不同加载速率下参数

      Table  3   Parameters at different loading rates

      加载速率/(mm·min−1)Don
      0.010.007 706.823 46
      0.100.008 539.515 30
      1.000.006 3233.347 92
      10.000.005 8728.195 12
      下载: 导出CSV
    • [1] 李庆文,高森林,胡露露,等. 不同加载速率下非均质煤样能量耗散损伤本构关系[J]. 煤炭学报,2022,47(S1):90−102.

      LI Qingwen, GAO Senlin, HU Lulu, et al. Constitutive relation of energy dissipation damage of heterogeneous coal samples under different loading rates[J]. Journal of China Coal Society, 2022, 47(S1): 90−102.

      [2] 杨胜利,王兆会,蒋威,等. 高强度开采工作面煤岩灾变的推进速度效应分析[J]. 煤炭学报,2016,41(3):586−594.

      YANG Shengli, WANG Zhaohui, JIANG Wei, et al. Advancing rate effect on rock and coal failure format in high intensity mining face[J]. Journal of China Coal Society, 2016, 41(3): 586−594.

      [3] 宋录生,赵善坤,刘军,等. “顶板-煤层”结构体冲击倾向性演化规律及力学特性试验研究[J]. 煤炭学报,2014,39(S1):23−30.

      SONG Lusheng, ZHAO Shankun, LIU Jun, et al. Experimental research on rules of rock burst tendency evolution and mechanical properties of “roof-coal” structure body[J]. Journal of China Coal Society, 2014, 39(S1): 23−30.

      [4] 肖晓春,丁鑫,赵鑫,等. 加载速率影响的煤体破裂过程声-电荷试验研究[J]. 岩土力学,2017,38(12):3419−3426.

      XIAO Xiaochun, DING Xin, ZHAO Xin, et al. Experimental study on acoustic emission and charge signals during coal failure process at different loading rates[J]. Rock and Soil Mechanics, 2017, 38(12): 3419−3426.

      [5]

      WANG L, CHENG L B, CHENG Y P, et al. Characteristics and evolutions of gas dynamic disaster under igneous intrusions and its control technologies[J]. Journal of Natural Gas Science & Engineering, 2014, 18: 164−174.

      [6] 夏冬,杨天鸿,王培涛,等. 干燥及饱和岩石循环加卸载过程中声发射特征试验研究[J]. 煤炭学报,2014,39(7):1243−1247.

      XIA Dong, YANG Tianhong, WANG Peitao, et al. Experimental study of acoustic emission characteristics of dry and saturated rocks during cyclic loading and unloading process[J]. Journal of China Coal Society, 2014, 39(7): 1243−1247.

      [7]

      HUANG B X, LIU J W. The effect of loading rate on the behavior of samples composed of coal and rock[J]. International Journal of Rock Mechanics and Mining Sciences, 2013, 61: 23−30. doi: 10.1016/j.ijrmms.2013.02.002

      [8] 王晓东,王坤. 加载速率效应对花岗岩破裂的力学性能及能量转化机制的影响[J]. 煤矿安全,2020,51(4):31−35.

      WANG Xiaodong, WANG Kun. Effect of loading rate on mechanical properties and energy conversion mechanism of granite fracture[J]. Safety in Coal Mines, 2020, 51(4): 31−35.

      [9] 王业平,刘相文,陈德任,等. 不同加载速率下含瓦斯突出煤体受载破坏损伤特征研究[J]. 煤矿安全,2020,51(6):25−30.

      WANG Yeping, LIU Xiangwen, CHEN Deren, et al. Damage characteristics of outburst coal sample containing gas under different loading rates[J]. Safety in Coal Mines, 2020, 51(6): 25−30.

      [10]

      LU Z G, JU W J, GAO F Q, et al. Influence of loading rate on the failure characteristics of composite coal–rock specimens under quasi-static loading conditions[J]. Rock Mechanics and Rock Engineering, 2022, 55(2): 909−921. doi: 10.1007/s00603-021-02699-2

      [11]

      ZHU C Q, LI S B, LUO Y, et al. Progressive damage process and failure characteristics of coal under uniaxial compression with different loading rates[J]. Shock and Vibration, 2021, 2021: 3360738.

      [12] 李海涛,蒋春祥,姜耀东,等. 加载速率对煤样力学行为影响的试验研究[J]. 中国矿业大学学报,2015,44(3):430−436. doi: 10.13247/j.cnki.jcumt.000326

      LI Haitao, JIANG Chunxiang, JIANG Yaodong, et al. Mechanical behavior and mechanism analysis of coal samples based on loading rate effect[J]. Journal of China University of Mining & Technology, 2015, 44(3): 430−436. doi: 10.13247/j.cnki.jcumt.000326

      [13] 李海涛,宋力,周宏伟,等. 率效应影响下煤的冲击特性评价方法及应用[J]. 煤炭学报,2015,40(12):2763−2771. doi: 10.13225/j.cnki.jccs.2015.0221

      LI Haitao, SONG Li, ZHOU Hongwei, et al. Evaluation method and application of coal burst performance under the effect of loading rate[J]. Journal of China Coal Society, 2015, 40(12): 2763−2771. doi: 10.13225/j.cnki.jccs.2015.0221

      [14] 马振乾,梁旭超,池恩安,等. 不同加载速率下煤的侧向变形特征研究[J]. 应用力学学报,2020,37(5):2007−2012.

      MA Zhenqian, LIANG Xuchao, CHI Enan, et al. Experimental study on lateral deformation characteristics of coal under different loading rates[J]. Chinese Journal of Applied Mechanics, 2020, 37(5): 2007−2012.

      [15]

      LI Y, YANG S Q, LIU Z L, et al. Study on mechanical properties and deformation of coal specimens under different confining pressure and strain rate[J]. Theoretical and Applied Fracture Mechanics, 2022, 118: 103287. doi: 10.1016/j.tafmec.2022.103287

      [16]

      GONG S. Investigation of tensile and fracture mechanical properties of bituminous coal at different strain rates[J]. Journal of Materials Research and Technology, 2021, 15: 834−845. doi: 10.1016/j.jmrt.2021.08.076

      [17]

      ZHANG M B, CUI L, HU W J, et al. Acoustic emission experimental research of the damage characteristics of raw coal under different loading and unloading rates[J]. Shock and Vibration, 2020, 2020: 9063929.

      [18] 王凯,蒋一峰,徐超. 不同含水率煤体单轴压缩力学特性及损伤统计模型研究[J]. 岩石力学与工程学报,2018,37(5):1070−1079.

      WANG Kai, JIANG Yifeng, XU Chao. Mechanical properties and statistical damage model of coal with different moisture contents under uniaxial compression[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(5): 1070−1079.

      [19] GB/T 50266—2013工程岩体试验方法标准[S].
      [20]

      XIAO W J, ZHANG D M, CAI Y, et al. Study on loading rate dependence of the coal failure process based on uniaxial compression test[J]. Pure and Applied Geophysics, 2020, 177(10): 4925−4941. doi: 10.1007/s00024-020-02513-0

      [21] 李庆文,胡露露,曹行,等. CFRP布均匀约束煤圆柱轴压性能[J]. 复合材料学报,2022,39(11):5611−5624.

      LI Qingwen, HU Lulu, CAO Hang, et al. Axial compressive behavior of CFRP uniformly wrapped coal in circular columns[J]. Acta Materiae Compositae Sinica, 2022, 39(11): 5611−5624.

      [22] 刘鹏飞,范俊奇,郭佳奇,等. 三轴应力下花岗岩加载破坏的能量演化和损伤特征[J]. 高压物理学报,2021,35(2):44−53.

      LIU Pengfei, FAN Junqi, GUO Jiaqi, et al. Damage and energy evolution characteristics of granite under triaxial stress[J]. Chinese Journal of High Pressure Physics, 2021, 35(2): 44−53.

      [23]

      HAWKINS A B, MCCONNELL B J. Sensitivity of sandstone strength and deformability to changes in moisture content[J]. Quarterly Journal of Engineering Geology and Hydrogeology, 1992, 25(2): 115−130. doi: 10.1144/GSL.QJEG.1992.025.02.05

      [24] 陈子全,何川,吴迪,等. 深埋碳质千枚岩力学特性及其能量损伤演化机制[J]. 岩土力学,2018,39(2):445−456.

      CHEN Ziquan, HE Chuan, WU Di, et al. Mechanical properties and energy damage evolution mechanism of deep-buried carbonaceous phyllite[J]. Rock and Soil Mechanics, 2018, 39(2): 445−456.

      [25]

      SOLECKI R, CONANT R J. Advanced mechanics of materials[M]. London: Oxford University Press, 2003.

      [26] 徐佳佳,方海,韩娟,等. 格构腹板增强泡沫夹芯复合材料准静态压缩吸能特性[J]. 复合材料学报,2022,39(8):3965−3981.

      XU Jiajia, FANG Hai, HAN Juan, et al. Energy absorption behavior of foam-filled sandwich composite materials reinforced by lattice webs under quasi-static compression[J]. Acta Materiae Compositae Sinica, 2022, 39(8): 3965−3981.

      [27] 谢和平,鞠杨,黎立云. 基于能量耗散与释放原理的岩石强度与整体破坏准则[J]. 岩石力学与工程学报,2005(17):3003−3010. doi: 10.3321/j.issn:1000-6915.2005.17.001

      XIE Heping, JU Yang, LI Liyun. Criteria for strength and structural failure of rocks based on energy dissipation and energy release principles[J]. Chinese Journal of Rock Mechanics and Engineering, 2005(17): 3003−3010. doi: 10.3321/j.issn:1000-6915.2005.17.001

      [28]

      HUANG P, ZHANG J X, DAMASCENE N J, et al. Effect of loading rate on mechanical behavior of coal samples with initial damage accumulation[J]. Mechanics of Time-dependent Materials, 2022, 26(2): 309−322. doi: 10.1007/s11043-021-09489-x

      [29] 赵红鹤,高富强,杨小林. 基于不同分布的分段式岩石损伤本构模型[J]. 矿业研究与开发,2015,35(4):64−67.

      ZHAO Honghe, GAO Fuqiang, YANG Xiaolin. Study on segmentary damage constitutive model based on different distribution[J]. Mining Research and Development, 2015, 35(4): 64−67.

      [30] 鲁祖德. 裂隙岩石水-岩作用力学特性试验研究与理论分析[D]. 武汉: 中国科学院武汉岩土力学研究所, 2010.
      [31]

      KACHANOV L M. Rupture time under creep conditions[J]. International Journal of Fracture, 1999, 97(1): 11−18.

      [32] 刘保县,赵宝云,姜永东. 单轴压缩煤岩变形损伤及声发射特性研究[J]. 地下空间与工程学报,2007(4):647−650.

      LIU Baoxian, ZHAO Baoyun, JIANG Yongdong. Study of deformation-damage and acoustic emission character of coal-rock under uniaxial compression[J]. Chinese Journal of Underground Space and Engineering, 2007(4): 647−650.

      [33] 王超,袁旭朋,汪开旺. 煤和砂岩加载弹塑性损伤本构方程[J]. 煤田地质与勘探,2017,45(2):105−111. doi: 10.3969/j.issn.1001-1986.2017.02.019

      WANG Chao, YUAN Xupeng, WANG Kaiwang. Constitutive equation of loading-induced elastic-plastic damage of coal and sandstone[J]. Coal Geology & Exploration, 2017, 45(2): 105−111. doi: 10.3969/j.issn.1001-1986.2017.02.019

    • 期刊类型引用(2)

      1. 王金成. 煤层注水对煤矿粉尘和瓦斯灾害的影响分析. 煤. 2025(02): 54-57 . 百度学术
      2. 郭庆彪,谢扬,汪锋,郑美楠. 废弃采空区封存CO_2地表形变特征模拟研究. 煤矿安全. 2024(05): 1-10 . 本站查看

      其他类型引用(0)

    图(7)  /  表(3)
    计量
    • 文章访问数:  27
    • HTML全文浏览量:  4
    • PDF下载量:  10
    • 被引次数: 2
    出版历程
    • 收稿日期:  2022-09-15
    • 网络出版日期:  2023-11-21
    • 刊出日期:  2023-11-19

    目录

    /

    返回文章
    返回