不同炮孔间距对岩石爆炸裂纹扩展影响的数值分析

李洪伟 雷战 江向阳 刘伟 何志伟 章彬彬

引用本文:
Citation:

不同炮孔间距对岩石爆炸裂纹扩展影响的数值分析

    作者简介: 李洪伟(1979-),男,硕士研究生,副教授,主要从事控制爆破技术研究. E-mail: 1227002529@qq.com;
    通讯作者: 雷战, 15055403215@163.com
  • 中图分类号: 0346.1

Numerical Analysis of Impact of Shot Hole Spacing on Crack Growth in Rock

    Corresponding author: LEI Zhan, 15055403215@163.com ;
  • CLC number: 0346.1

  • 摘要: 为实现对岩石的充分破碎,有效利用炸药能量,基于有限元分析软件ANSYS/LS-DYNA及流固耦合(ALE)算法,研究了不同炮孔间距对岩石(角岩)爆炸裂纹扩展的影响,同时将模拟结果应用在工程实践上加以验证。结果表明:随着两炮孔间距离的增大,单炮孔周围裂纹扩展更加充分,炮孔周围粉碎区增大,萌生的分支小裂纹逐渐减少,主裂纹增多。在两炮孔之间受到相邻炮孔爆炸应力波影响的区域,主裂纹发育扩展较为明显,且随着炮孔间距的增大,主裂纹相互贯通的位置越靠近两炮孔中心连线方向。工程实践表明:数值模拟结果与爆破工程效果具有较好的吻合性,将数值模拟结果用来指导爆破方案设计是可行的,能够为爆破工程提供重要的参考。
  • 图 1  计算模型

    Figure 1.  Calculation model

    图 2  L=40 cm时岩石裂纹扩展过程

    Figure 2.  Growth process of rock crack at L=40 cm

    图 3  L=60 cm时岩石裂纹扩展过程

    Figure 3.  Growth process of rock crack at L=60 cm

    图 4  L=80 cm时岩石裂纹扩展过程

    Figure 4.  Growth process of rock crack at L=80 cm

    图 5  不同孔间距下炮孔中心连线中点处x方向的压力曲线

    Figure 5.  Pressure curve in the x direction at the center of the connecting line of gun hole with different hole spacing

    图 6  起爆网路连接示意图

    Figure 6.  Schematic of detonating network

    图 7  爆破效果图

    Figure 7.  Blasting effect

    表 1  岩石参数[6]

    Table 1.  Material parameters of rock[6]

    ${\rho}$/(kg·m–3)E0/GPa${\mu}$${{\sigma _0}}$/MPaEtan/GPa${{\sigma _{\rm{c}}}}$/MPa${{\sigma _{{\rm{st}}}}}$/MPaC/s–1P
    270068.690.22875401505.62.633.96
    下载: 导出CSV

    表 2  炸药参数[10-12]

    Table 2.  Material parameters of explosive[10-12]

    ${\rho}$/(kg·m–3)D/(m·s–1)pCJ/GPaA/GPaB/GPaR1R2${\omega}$
    120040004.802140.0934.150.950.3
    下载: 导出CSV

    表 3  岩石裂隙扩展平均速度

    Table 3.  The average speed of rock crack extension

    Speed of rock crack extension/(m·s−1V/Cp
    L=40 cmL=60 cmL=80 cm
    2051(t=40 ${{\text{μ}}{\rm{s}}}$)2157(t=60 ${{\text{μ}}{\rm{s}}}$)2286(t=70 ${{\text{μ}}{\rm{s}}}$)0.39
    2197(t=70 ${{\text{μ}}{\rm{s}}}$)2083(t=120 ${{\text{μ}}{\rm{s}}}$)2285(t=140 ${{\text{μ}}{\rm{s}}}$)0.41
    1656(t=230 ${{\text{μ}}{\rm{s}}}$)1923(t=210 ${{\text{μ}}{\rm{s}}}$)2140(t=190 ${{\text{μ}}{\rm{s}}}$)0.44
    1968(Average)2054(Average)2237(Average)0.41
    下载: 导出CSV

    表 4  爆破参数表

    Table 4.  Blasting parameters

    H/mh/ma/mb/mL/ml/mQ/kg
    1.702.01.20.201.500.9
    2.002.01.20.221.781.0
    2.502.51.50.452.052.0
    3.003.01.50.672.333.0
    下载: 导出CSV

    表 5  振动监测表

    Table 5.  Blasting vibration monitoring

    Measuring pointDistance/mvx/(cm·s–1)fx/Hzvy/(cm·s–1)fy/Hzvz/(cm·s–1)fz/Hz
    1163–0.219.16–0.289.160.2718.31
    21540.048.850.0435.710.039.48
    33390.0412.160.057.830.0613.29
    42340.1918.31–0.149.16–0.1118.31
    下载: 导出CSV
  • [1] MA G W, HAO H, WANG F. Simulations of explosion-induced damage to underground rock chambers [J]. Journal of Rock Mechanics and Geotechnical Engineering, 2011, 3(1): 19–29. doi: 10.3724/SP.J.1235.2011.00019
    [2] YUE Z W, QIU P, YANG R S, et al. Stress analysis of the interaction of a running crack and blasting waves by caustics method [J]. Engineering Fracture Mechanics, 2017, 184: 339–351. doi: 10.1016/j.engfracmech.2017.08.037
    [3] 何成龙, 杨军. 主动围压和爆炸加载作用下岩石动态响应研究 [J]. 兵工学报, 2017, 38(12): 2395–2405. doi: 10.3969/j.issn.1000-1093.2017.12.013
    HE C L, YANG J. Research on dynamic response of rock under blast loading and active confining pressure [J]. Acta Armamentarii, 2017, 38(12): 2395–2405. doi: 10.3969/j.issn.1000-1093.2017.12.013
    [4] ZHAO J J, ZHANG Y, RANJITH P G. Numerical simulation of blasting-induced fracture expansion in coal masses [J]. International Journal of Rock Mechanics and Mining Sciences, 2017, 100: 28–39. doi: 10.1016/j.ijrmms.2017.10.015
    [5] ZHANG X L, JIAO Y Y, MA J F. Simulation of rock dynamic failure using discontinuous numerical approach [J]. Computers and Geotechnics, 2018, 96: 160–166. doi: 10.1016/j.compgeo.2017.10.001
    [6] 夏祥, 李海波, 李俊如, 等. 岩体爆生裂纹的数值模拟 [J]. 岩土力学, 2006, 27(11): 1987–1991. doi: 10.3969/j.issn.1000-7598.2006.11.025
    XIA X, LI H B, LI J R, et al. Numerical simulation of blast-induced cracks in rock [J]. Rock and Soil Mechanics, 2006, 27(11): 1987–1991. doi: 10.3969/j.issn.1000-7598.2006.11.025
    [7] 周艳, 叶海旺, 易长平, 等. 台阶爆破相邻炮孔间损伤范围的数值模拟 [J]. 工程爆破, 2014, 20(2): 17–20, 28. doi: 10.3969/j.issn.1006-7051.2014.02.005
    ZHOU Y, YE H W, YI C P, et al. Numerical simulation of damage zone between the adjacent boreholes in bench blasting [J]. Engineering Blasting, 2014, 20(2): 17–20, 28. doi: 10.3969/j.issn.1006-7051.2014.02.005
    [8] 郭东明, 闫鹏洋, 杨仁树, 等. 动静荷载下巷道围岩倾斜裂纹的动焦散试验 [J]. 采矿与安全工程学报, 2016, 33(4): 668–675.
    GUO D M, YAN P Y, YANG R S, et al. Dynamic caustic test of inclined crack in roadwaysurrounding rock under dynamic and static load [J]. Journal of Mining and Safety Engineering, 2016, 33(4): 668–675.
    [9] 钟波波, 李宏, 张永彬. 爆炸荷载作用下岩石动态裂纹扩展的数值模拟 [J]. 爆炸与冲击, 2016, 36(6): 825–831. doi: 10.11883/1001-1455(2016)06-0825-07
    ZHONG B B, LI H, ZHANG Y B. Numerical simulation of dynamic cracks propagation of rock under blasting loading [J]. Explosion and Shock Waves, 2016, 36(6): 825–831. doi: 10.11883/1001-1455(2016)06-0825-07
    [10] 门建兵, 蒋建伟, 王树有. 爆炸冲击数值模拟技术基础 [M]. 北京: 北京理工大学出版社, 2015.
    MEN J B, JIANG J W, WANG S Y. Numerical simulation of explosion impact [M]. Beijing: Beijing Institute of Technology Press, 2015.
    [11] 钮强. 岩石爆破机理[M]. 沈阳: 东北工学院出版社, 1990.
    NIU Q. Rock blasting mechanism [M]. Shenyang: Northeast Institute of Technology Press, 1990.
    [12] LI X F, LI H B, ZHAO J. 3D polycrystalline discrete element method (3PDEM) for simulation of crack initiation and propagation in granular rock [J]. Computers and Geotechnics, 2017, 90: 96–112. doi: 10.1016/j.compgeo.2017.05.023
    [13] YILMAZ O, UNLU T. Three dimensional numerical rock damage analysis under blasting load [J]. Tunnelling and Underground Space Technology, 2013, 38: 266–278. doi: 10.1016/j.tust.2013.07.007
    [14] WANG Z L, KONIETZKY H. Modelling of blast-induced fractures in jointed rock masses [J]. Engineering Fracture Mechanics, 2009, 76(12): 1945–1955. doi: 10.1016/j.engfracmech.2009.05.004
    [15] ZHU W C, BAI Y, LI X B. Numerical simulation on rock failure under combined static and dynamic loading during SHPB tests [J]. International Journal of Impact Engineering, 2012, 49: 142–157. doi: 10.1016/j.ijimpeng.2012.04.002
    [16] 王礼立, 朱兆祥. 应力波基础 [M]. 2版. 北京: 国防工业出版社, 2005.
    WANG L L, ZHU Z X. Stress wave foundation [M]. 2nd ed. Beijing: National Defense Industry Press, 2005.
    [17] 陈宝心, 杨勤荣. 爆破动力学基础 [M]. 武汉: 湖北科学技术出版社, 2005.
    CHEN B X, YANG Q R. Basis of blasting dynamics [M]. Wuhan: Hubei Science and Technology Press, 2005.
  • [1] 傅华刘仓理王文强李涛 . 冲击动力学中离散元与有限元相结合的计算方法研究. 高压物理学报, 2006, 20(4): 379-385 . doi: 10.11858/gywlxb.2006.04.007
    [2] 毛勇建李玉龙陈颖黄含军张青平苗应刚 . 炸药条加载圆柱壳的数值模拟(Ⅰ):流固耦合模拟. 高压物理学报, 2012, 26(2): 155-162. doi: 10.11858/gywlxb.2012.02.006
    [3] 姚蓬飞韩阳姚芬李志强 . 基于固有型内聚力模型模拟双层夹胶玻璃冲击断裂行为. 高压物理学报, 2019, 33(6): 064105-1-064105-11. doi: 10.11858/gywlxb.20190718
    [4] 江松青刘文韬 . 钽的层裂实验数值模拟. 高压物理学报, 2007, 21(2): 113-120 . doi: 10.11858/gywlxb.2007.02.001
    [5] 卞梁王肖钧章杰 . SPH/FEM耦合算法在陶瓷复合靶抗侵彻数值模拟中的应用. 高压物理学报, 2010, 24(3): 161-167 . doi: 10.11858/gywlxb.2010.03.001
    [6] 张祎杨春霞栗保明 . 电磁轨道炮C型固体电枢坡膛段的装填方式研究. 高压物理学报, 2012, 26(1): 48-54. doi: 10.11858/gywlxb.2012.01.007
    [7] 尹晓文田泽韩阳李志强 . 工业管道系统流致疲劳寿命的数值模拟. 高压物理学报, 2018, 32(6): 064102-1-064102-9. doi: 10.11858/gywlxb.20180559
    [8] 种涛唐志平谭福利王桂吉赵剑衡 . 纯铁相变和层裂损伤的数值模拟. 高压物理学报, 2018, 32(1): 014102-1-014102-5. doi: 10.11858/gywlxb.20170528
    [9] 肖桂凤唐志平周昌国胡晓军蔡建 . 氮化硅陶瓷层裂强度的研究. 高压物理学报, 2005, 19(3): 219-224 . doi: 10.11858/gywlxb.2005.03.005
    [10] 张忠珍孙承纬丰树平高顺受 . 短脉冲载荷下金属材料的层裂破坏. 高压物理学报, 1991, 5(2): 137-144 . doi: 10.11858/gywlxb.1991.02.009
    [11] 董玉斌苏林祥陈大年经福谦韩钧万封加波 . 滑移爆轰作用下内爆柱形钢壳层裂的数值模拟. 高压物理学报, 1989, 3(1): 1-10 . doi: 10.11858/gywlxb.1989.01.001
    [12] 祁美兰罗超范端贺红亮 . 高纯铝层裂的自由面速度剖面特征分析. 高压物理学报, 2011, 25(6): 487-492. doi: 10.11858/gywlxb.2011.06.002
    [13] 陈涛蒋招绣辛铭之申海艇王永刚 . 预加温条件下高强铝合金动态屈服及层裂行为的研究. 高压物理学报, 2015, 29(5): 321-328. doi: 10.11858/gywlxb.2015.05.001
    [14] 李庆忠陈永涛胡海波徐永波 . 冲击加载下FeMnNi合金相变和层裂特性实验研究. 高压物理学报, 2010, 24(2): 107-112 . doi: 10.11858/gywlxb.2010.02.005
    [15] 谭福利李永池郭扬曹结东 . 脉冲激光产生层裂的一维和二维数值计算. 高压物理学报, 2005, 19(1): 24-28 . doi: 10.11858/gywlxb.2005.01.005
    [16] 孙宇新张进李永池胡时胜董杰 . 内爆加载下热塑性管壳的应力波演化与层裂效应研究. 高压物理学报, 2005, 19(4): 319-324 . doi: 10.11858/gywlxb.2005.04.006
    [17] 周洪强秦承森张树道张凤国刘文韬王裴 . 一个新的损伤演化方程及其在层裂中的应用. 高压物理学报, 2005, 19(3): 201-205 . doi: 10.11858/gywlxb.2005.03.002
    [18] 陈永涛李庆忠唐小军卢秋虹 . 给定压力下纯铁材料的冲击相变与层裂特性研究. 高压物理学报, 2010, 24(3): 213-218 . doi: 10.11858/gywlxb.2010.03.009
    [19] 李庆忠陈永涛汤铁钢 . FeMnNi合金高压加卸载历程测量和相变层裂分析. 高压物理学报, 2009, 23(2): 117-122 . doi: 10.11858/gywlxb.2009.02.007
    [20] 王志华李志强任会兰赵隆茂 . 有机玻璃的几何形状对爆炸切割影响的研究. 高压物理学报, 2010, 24(1): 43-48 . doi: 10.11858/gywlxb.2010.01.008
  • 加载中
图(7)表(5)
计量
  • 文章访问数:  592
  • 阅读全文浏览量:  538
  • PDF下载量:  14
出版历程
  • 收稿日期:  2018-11-09
  • 录用日期:  2018-11-19
  • 网络出版日期:  2019-07-24
  • 刊出日期:  2019-08-01

不同炮孔间距对岩石爆炸裂纹扩展影响的数值分析

    作者简介:李洪伟(1979-),男,硕士研究生,副教授,主要从事控制爆破技术研究. E-mail: 1227002529@qq.com
    通讯作者: 雷战, 15055403215@163.com
  • 安徽理工大学化工学院,安徽 淮南 232001

摘要: 为实现对岩石的充分破碎,有效利用炸药能量,基于有限元分析软件ANSYS/LS-DYNA及流固耦合(ALE)算法,研究了不同炮孔间距对岩石(角岩)爆炸裂纹扩展的影响,同时将模拟结果应用在工程实践上加以验证。结果表明:随着两炮孔间距离的增大,单炮孔周围裂纹扩展更加充分,炮孔周围粉碎区增大,萌生的分支小裂纹逐渐减少,主裂纹增多。在两炮孔之间受到相邻炮孔爆炸应力波影响的区域,主裂纹发育扩展较为明显,且随着炮孔间距的增大,主裂纹相互贯通的位置越靠近两炮孔中心连线方向。工程实践表明:数值模拟结果与爆破工程效果具有较好的吻合性,将数值模拟结果用来指导爆破方案设计是可行的,能够为爆破工程提供重要的参考。

English Abstract

  • 爆炸作用过程是一个高温高压的瞬态过程,装药爆炸以后爆轰产物的压力高达10 GPa级,强大的冲击波压应力使炮孔周围岩石受压破碎,瞬间完成压缩破碎和裂隙扩展[1]。在工程实践中,例如,露天矿山台阶爆破、巷道掘进爆破、拆除爆破等一般有成百上千个炮孔,炮孔间距偏小会造成爆炸能量的极度浪费,炮孔间距过大会造成岩石破碎不充分,增加工程成本和爆破安全隐患等一系列问题,因此合理选择炮孔间距对爆破起着至关重要的作用[23]。目前炮孔间距的选择主要是以经验公式为参考依据,缺乏一定的理论支撑,研究双炮孔甚至多炮孔间裂纹扩展问题及孔间距的选择成为了研究的热点和难点问题。

    针对不同炮孔间距间岩体破坏效应问题,Zhao等[4]研究了距离炮孔不同位置、不同方向的节理面对动态裂纹扩展方向的影响,并讨论了不同微差起爆时间对煤体裂纹扩展的影响;Zhang等[5]基于不连续数值模拟的方法对动态加载下的岩石破坏过程进行了研究;夏祥等[6]通过ANSYS/LS-DYNA模拟软件对单孔柱状装药下岩体粉碎区和裂隙区的形成过程及范围进行了研究,并通过理论计算对模拟结果进行了验证;周艳等[7]通过建立三合一新站石灰石矿台阶松动爆破模型,对相邻炮孔间的损伤范围进行研究,发现两炮孔间的损伤百分比随孔距的增大而减小,最后逐渐趋于定值;郭东明等[8]采用透射式动态焦散线实验系统对邻近巷道裂纹缺陷受爆炸荷载作用的扩展机理进行了研究;钟波波等[9]运用RFPA-Dynamic分析软件就爆炸载荷加载速率、炮孔到自由边界距离及炮空间空孔大小3方面对岩石动态裂纹的扩展进行了研究。

    以上工作为岩体在爆炸冲击载荷作用下的破坏分析提供了重要的参考。鉴于岩石在爆炸载荷作用下破坏的复杂性,在前人工作的基础上,本研究基于ANSYS/LS-DYNA有限元软件,采用Mises屈服准则,利用流固耦合(Arbitrary Lagrange-Euler,ALE)算法模拟不同炮孔间距离对岩石在爆炸载荷下裂纹扩展的影响过程,得到炮孔间应力波的叠加曲线,以期为充分理解爆炸应力波与裂纹间的相互影响过程提供理论依据,模拟结果对于岩石爆破相关工程具有一定的参考价值。

    • 炸药爆炸时产生的爆轰波和高温高压气体作用在炮孔壁上,在爆炸近区,岩石中激起的应力波强度远大于岩石的动态抗压强度,使炮孔周围岩石呈塑性状态,并且爆破工程中岩石的加载应变率$\dot \varepsilon $在1~105 s–1之间,应变率效应明显,塑性硬化模型本构简单,岩石参数易于获取,故采用含应变率效应的塑性硬化模型比较合适。

      岩石选用LS-DYNA自带的*MAT_PLASTIC_KINEMATIC材料[10]模型,该模型中应变率用Cowper-Symonds模型表示,采用与应变率有关的函数表示屈服应力

      ${\sigma _{\rm{y}}} = \left[ {1 + {{\left( {\frac{{\dot\varepsilon}}{C}} \right)}^{\frac{1}{P}}}} \right]\left( {{\sigma _0} + \beta {E\rm_p}\varepsilon _{\rm{p}}^{{\rm{eff}}}} \right)$

      式中:${\sigma _0}$为岩石初始屈服应力,$\dot\varepsilon$为应变率,CP为应变率参数,$\varepsilon _{\rm{p}}^{{\rm{eff}}}$为岩石有效塑性应变,$\beta $为硬化参数,0 ≤ $\beta $ ≤ 1,Ep为岩石塑性硬化模量。

      Ep的表达式为

      $E_{\rm{p}} = \frac{{{E_{{\rm{tan}}}} \cdot {E_0}}}{{{E_0} - {E_{{\rm{tan}}}}}}$

      式中:Etan为切线模量,E0为杨氏模量。

      岩石基本力学参数选用文献[6]所测角岩动态力学参数,具体参数如表1所示,其中:$\rho $为岩石密度,$\mu $为泊松比,${\sigma _{\rm{c}}}$为岩石静态抗压强度,${\sigma _{{\rm{st}}}}$为岩石静态抗拉强度。

      ${\rho}$/(kg·m–3)E0/GPa${\mu}$${{\sigma _0}}$/MPaEtan/GPa${{\sigma _{\rm{c}}}}$/MPa${{\sigma _{{\rm{st}}}}}$/MPaC/s–1P
      270068.690.22875401505.62.633.96

      表 1  岩石参数[6]

      Table 1.  Material parameters of rock[6]

    • 在ANSYS/LS-DYNA中选用*MAT_HIGH_EXPLOSIVE_BURN材料定义炸药爆轰模型,同时选用JWL状态方程,具体表述为

      $p = A\left( {1 - \frac{\omega }{{{R_1}V}}} \right){{\rm{e}}^{ - {R_1}V}} + B\left( {1 - \frac{\omega }{{{R_2}V}}} \right){{\rm{e}}^{ - {R_2}V}} + \frac{{\omega {E_0}}}{V}$

      式中:p为爆轰产物内部压力,V为爆轰产物相对体积,E0为初始比内能,A、B、R1R2$\omega $为与炸药有关的常数。炸药具体参数[1012]表2所示,其中D为爆速,pCJ为炸药爆轰的CJ压力。

      ${\rho}$/(kg·m–3)D/(m·s–1)pCJ/GPaA/GPaB/GPaR1R2${\omega}$
      120040004.802140.0934.150.950.3

      表 2  炸药参数[10-12]

      Table 2.  Material parameters of explosive[10-12]

    • 炸药在岩石中爆炸时会导致岩体发生大变形甚至断裂破坏,为此本研究选用多物质Euler材料和Lagraner结构相耦合的算法,炸药为欧拉算法,岩石为拉格朗日算法,同时鉴于模型的对称性,建立1/2平面计算模型,在对称面节点施加对称约束,模型尺寸为300 cm×300 cm,计算时间为2000 ${\text{μ}}{\rm{s}}$,炮孔直径D=70 mm,炮孔间距为L

      自由面数目对爆破作用效果的影响较大,当爆炸应力波遇到自由面发生反射时,压缩应力波变为拉伸波,引起岩石的片落和径向裂隙的延伸,同时改变岩石的应力状态和强度极限[13]。为了提高炸药对岩石的破坏作用以及观察从边界反射的拉伸波对岩石的二次甚至多次破坏作用,岩石边界均设置为自由边界。物理计算模型如图1所示。

      图  1  计算模型

      Figure 1.  Calculation model

    • 图2(a)图2(i)分别给出了在双炮孔间距为40 cm时岩石裂纹扩展过程。

      图  2  L=40 cm时岩石裂纹扩展过程

      Figure 2.  Growth process of rock crack at L=40 cm

      图2可以看出:炮孔起爆后,随着应力波向孔外传播,在炮孔周围逐渐形成破碎区,岩石裂纹在孔壁边缘逐渐形成,并呈放射状向四周扩展;在t = 40 ${\text{μ}}{\rm{s}}$ 时,两应力波波阵面相遇,在应力波波阵面切线方向产生合拉应力,在应力波传播过程中,当合成应力波的大小超过岩石抗拉强度时,首先在两个炮孔之间形成径向裂隙,直至裂纹贯穿两个炮孔。

      图2(e)图2(g)所示,当爆炸应力波从边界反射回来时,岩石主要受到拉伸应力波的影响,岩石裂隙进一步发育扩展形成二次破坏裂隙。当两炮孔叠加应力波传播至上下边界时,由于应力波的反射拉伸作用,使距离上下边界一定位置处出现层裂(剥落)现象,如图2(e)所示,层裂主裂纹方向与边界方向基本平行,与文献[14]的模拟结果一致。

    • 图3(a)图3(i)分别给出了在双炮孔间距为60 cm时岩石裂纹扩展过程。通过与图2对比发现,在双炮孔起爆所产生的应力波阵面未相遇前,岩石裂纹的形成和扩展过程与图2基本一致。随着炮孔间距离L的增大,单个炮孔周围裂纹扩展得更加充分(如图3(c)所示),在两个炮孔之间受到爆炸应力波初次影响的区域,主裂纹逐渐扩展并相互贯通形成明显的裂隙,但微裂纹数目减少。

      图  3  L=60 cm时岩石裂纹扩展过程

      Figure 3.  Growth process of rock crack at L=60 cm

      图2对比可知,当L=60 cm时,从边界反射回来的爆炸应力波随着炮孔间距离的增大,岩石裂隙受到拉伸应力波产生的二次破坏作用较小,且在上下边界处引起层裂的主裂纹长度减小,主裂纹距离边界的距离缩短,即层裂厚度降低。

    • 图4(a)图4(i)分别给出了在炮孔间距为80 cm时岩石裂纹的扩展过程。与图2、图3相比较可知,在两炮孔间距离为80 cm时,单炮孔岩石裂纹扩展得更加充分,两炮孔之间的主裂纹逐渐扩展并相互贯通形成明显的裂隙,另外炮孔间因应力波叠加,在叠加处应力波幅值增大所产生的裂隙数目增多,岩石破碎效果较好。

      图  4  L=80 cm时岩石裂纹扩展过程

      Figure 4.  Growth process of rock crack at L=80 cm

      图2图3对比可知,当L=80 cm时,在上下边界处引起层裂的主裂纹长度最长,层裂厚度最大,如图4(i)所示。另外发现从边界反射回来的爆炸应力波随着炮孔间距离的增大,岩石裂隙受到拉伸应力波产生的二次破坏作用最小,说明随着孔间距的增大,岩石受到的爆炸应力波初次破坏作用所占比例增大,受到二次破坏作用所占比例降低。

      考虑到应力波对岩石的初次破坏作用较大,选择应力波从产生至到达岩石边界此段时间,计算岩石裂隙的扩展速度,岩石纵波速度的计算公式[15]

      ${C_{\rm{p}}} = \sqrt {\frac{{{E_0}}}{\rho }} $

      式中:Cp为岩石纵波速度,m/s;E0为杨氏模量,GPa;$\rho $为岩石密度,kg/m3

      通过(4)式计算得到岩石的纵波速度为5044 m/s,岩石裂隙扩展速度V表1所示。在孔距L 为40、60和80 cm时岩石裂隙扩展平均速度V分别为岩石纵波速度Cp的0.39、0.41和0.44倍,与杜良奈[11](Dulaneg)研究认为的V = 0.38Cp基本一致。从表3中可以看出,在炮孔间距一定时,岩石裂隙扩展速度随时间的增加逐渐减小,与岩石中爆炸应力波随距离的增加而衰减相符合。另外,随着炮孔间距的增加,岩石裂隙扩展平均速度呈现逐渐增大的趋势,L = 80 cm时岩石裂隙扩展平均速度VL为40、60 cm时相比,分别提高了13.7%和8.9%。

      Speed of rock crack extension/(m·s−1V/Cp
      L=40 cmL=60 cmL=80 cm
      2051(t=40 ${{\text{μ}}{\rm{s}}}$2157(t=60 ${{\text{μ}}{\rm{s}}}$2286(t=70 ${{\text{μ}}{\rm{s}}}$0.39
      2197(t=70 ${{\text{μ}}{\rm{s}}}$2083(t=120 ${{\text{μ}}{\rm{s}}}$2285(t=140 ${{\text{μ}}{\rm{s}}}$0.41
      1656(t=230 ${{\text{μ}}{\rm{s}}}$1923(t=210 ${{\text{μ}}{\rm{s}}}$2140(t=190 ${{\text{μ}}{\rm{s}}}$0.44
      1968(Average)2054(Average)2237(Average)0.41

      表 3  岩石裂隙扩展平均速度

      Table 3.  The average speed of rock crack extension

    • 岩石爆破中的压碎区是岩石受压所致,故在压碎区选用Mises屈服准则[16],而裂隙区主要是受应力波的拉伸作用和爆轰产物的尖劈作用形成的,假设岩石中任意一点的等效应力强度为

      ${\sigma _{\rm{i}}}{\rm{ = }}\frac{1}{{\sqrt 2 }}\left[ {{{\left( {{\sigma _1}{\rm{ - }}{\sigma _2}} \right)}^2}{\rm{ + }}{{\left( {{\sigma _2}{\rm{ - }}{\sigma _3}} \right)}^2}{\rm{ + }}{{\left( {{\sigma _1}{\rm{ - }}{\sigma _3}} \right)}^2}} \right]$

      当满足(6)式时,岩石被破坏。

      $\left\{\begin{aligned} & {\sigma _{\rm{i}}} \geqslant {\sigma _{{\rm{cd}}}}\quad {\text{压碎区}}\\ & {\sigma _{\rm{t}}} \geqslant {\sigma _{{\rm{td}}}} \quad {\text{裂隙区}} \end{aligned} \right.$

      岩石爆破中岩石动态抗压强度近似取

      ${\sigma _{{\rm{cd}}}} = {\sigma _{\rm{c}}}{\dot \varepsilon ^{\frac{1}{3}}}$

      岩石爆破中岩石动态抗拉强度近似取

      ${\sigma _{{\rm{td}}}} = {\sigma _{{\rm{st}}}}$

      式中:${\sigma _{{\rm{cd}}}}$为岩石单轴动态抗压强度,${\sigma _{{\rm{td}}}}$为岩石单轴动态抗拉强度,${\sigma _{\rm{t}}}$为岩体中任意一点在爆炸载荷作用下所受到的拉应力,${\sigma _{\rm{c}}}$为岩石的单轴抗压强度,${\sigma _{{\rm{st}}}}$为岩石的单轴抗拉强度。

    • 选取距离两炮孔中心L/2处一个单元,其在x方向的应力曲线如图5所示。当L=40 cm时,在10~40 ${\text{μ}}{\rm{s}}$内,所选单元在x方向受到的拉应力逐渐增大(压应力为正,拉应力为负),拉应力峰值为39.9 MPa(A点),到达A点以后拉应力急速减小,在70 ${\text{μ}}{\rm{s}}$D点)时受到的压缩应力为3.27 MPa,此段时间内应力平均加载速率约为1.33 MPa/${\text{μ}}{\rm{s}}$,平均卸载速率约为1.44 MPa/${\text{μ}}{\rm{s}}$;当L=60 cm时,所选单元在x方向受到的拉应力峰值为25.9 MPa(B点),相对于A点峰值应力衰减了35.1%;当L=80 cm时,所选单元在x方向受到的拉应力峰值为12.9 MPa,相对于A、B点峰值应力分别衰减67.7%和50.1%,且在E点以后,单元所受的拉应力和压缩应力动态转换,与岩石爆破动力学的研究结果[17]相符合。

      图  5  不同孔间距下炮孔中心连线中点处x方向的压力曲线

      Figure 5.  Pressure curve in the x direction at the center of the connecting line of gun hole with different hole spacing

      综合以上分析表明,岩石裂隙区主要是爆炸应力波拉伸破坏作用的结果,应力波峰值大小与炸药产生的爆轰压力成正比,与炮孔间距离成反比。虽然在孔间距较小时,两炮孔中心岩石受到的拉伸应力较大,但是易造成岩石的过度粉碎,岩石破坏范围较小,因此工程上应根据岩石性质、炸药性能等因素选取适当的炮孔间距,既能保证岩石被充分破碎,又能保证一定破坏的范围。

      岩石塑性区以外的部分尽管没有像塑性区内那样存在明显的裂隙网,但实际上岩石内本身就含有很多细微的裂隙(原生裂隙),这些细微的裂隙会加速应力波的衰减,使模拟结果往往与实际情况存在一定的差别,因此在数值模拟和理论计算的基础上,需要在实际应用中根据具体情况逐步进行修正,这也是有待于进一步深入研究的方向。

    • 本次爆破工程为文一·名门森林语土石方爆破工程,爆破位置位于安徽省合肥市肥西县紫蓬路与森林大道交口。区域面积约为104 m2,本工程土石方开挖量约2.45×104 m3。经过现场勘察,爆破区域北侧63 m为居民区,南侧124 m为居民区,东侧55 m为森林大道,西侧为施工工地,整体爆破环境良好。

    • 根据数值模拟得到的不同炮孔间距对岩石爆炸裂纹的影响规律,并综合考虑现场机械设备、岩石类型、要求开挖爆破高程、周边环境等诸多因素,选取的爆破参数及设计的起爆网路如下所示。

      (1)爆破参数

      此次爆破为岩石松动爆破,钻孔孔径为90 mm,爆破作用指数n=0.6,最小抵抗线W=1.2 m,单耗k取0.22 kg/m3,炮孔参数见表4,其中a为炮孔间距,b为炮孔排距,H 为炮孔深度,h为炮孔超挖深度,L为装药长度,l为堵塞长度,Q为单孔装药量。

      H/mh/ma/mb/mL/ml/mQ/kg
      1.702.01.20.201.500.9
      2.002.01.20.221.781.0
      2.502.51.50.452.052.0
      3.003.01.50.672.333.0

      表 4  爆破参数表

      Table 4.  Blasting parameters

      (2)爆破网路

      此次爆破采用可操作性强、安全可靠的非电导爆管雷管复式起爆网路,采用孔内9段(MS9)孔外4段(MS4)延期微差逐排起爆,起爆网路示意图见图6

      图  6  起爆网路连接示意图

      Figure 6.  Schematic of detonating network

    • (1)爆破振动

      工程中共设4个爆破振动监测点,采用2台四川拓普测控科技有限公司生产的NUBOX-8016型爆破振动智能监测仪、2台成都中科测控有限公司生产的TC-4850型爆破测振监测仪进行检测。

      爆破振动智能监测仪监测的数据见表5,其中vxvyvz分别为xyz方向的振动速度,fxfyfz分别为xyz方向的振动频率。可以看出,土石方爆破作业中,爆破振动速度远小于《爆破安全规程》(GB 6722–2014)中一般民用建筑物的安全允许质点振动速度标准(f≤10 Hz时振动速度v的范围是1.5~2.0 cm/s;10 Hz<f≤50 Hz时振动速度v的范围是2.0~2.5 cm/s)。

      Measuring pointDistance/mvx/(cm·s–1)fx/Hzvy/(cm·s–1)fy/Hzvz/(cm·s–1)fz/Hz
      1163–0.219.16–0.289.160.2718.31
      21540.048.850.0435.710.039.48
      33390.0412.160.057.830.0613.29
      42340.1918.31–0.149.16–0.1118.31

      表 5  振动监测表

      Table 5.  Blasting vibration monitoring

      (2)现场爆破效果

      起爆后,所有炮孔均被起爆(无盲炮),从中选择4个炮孔查看裂纹扩展情况,具体如图7所示。从图7可以看出, 每个炮孔周围有许多分支裂纹,表明岩石裂纹在孔壁边缘逐渐形成并呈放射状向四周扩展。4个炮孔之间共有8条主裂纹,且裂纹之间相互贯通,在每两个炮孔之间因应力波叠加效应而导致岩石发生部分破碎,这与模拟结果和理论分析一致,从而证明了数值模拟结果的可靠性。

      图  7  爆破效果图

      Figure 7.  Blasting effect

    • (1)爆炸应力波在未相遇前,随炮孔间距离的增加,单炮孔周围裂纹扩展得更加充分,炮孔周围粉碎区增大,萌生的分支小裂纹逐渐减少,主裂纹增多。在两炮孔之间受到相邻炮孔爆炸应力波影响的区域,主裂纹发育扩展较为明显,且随着炮孔间距的增大,主裂纹相互贯通的位置越靠近两炮孔中心连线方向。

      (2)在3种不同炮孔间距条件下,所选取单元x方向爆炸应力波峰值压力衰减规律基本一致,炮孔间距L=80 cm相对于L=60 cm和L=40 cm,应力波峰值分别衰减67.7%和50.1%。

      (3)由于岩石爆破裂隙区的形成主要是爆炸应力波拉伸应力作用的结果,因此工程上应根据岩石性质和炸药性能选择爆炸应力波衰减到岩石动态抗拉强度附近时对应的距离,即为较合理的炮孔间距。

      (4)工程实践结果表明,数值模拟结果用来指导爆破方案设计是可行的,能够为爆破工程提供重要的参考。

参考文献 (17)

目录

    /

    返回文章
    返回