运动炸药近地爆炸冲击波场特性研究(2).pdf
第35卷 第3期 2018年9月 爆 破 BLASTING Vol. 35 No. 3 Sep. 2018 doi10. 3963/ j. issn. 1001 -487X. 2018. 03. 008 运动炸药近地爆炸冲击波场特性研究* 畅 博, 谷鸿平, 牛晨伟, 吕永柱, 敬怡东 ( 西安近代化学研究所, 西安710065) 摘 要 为研究运动炸药近地爆炸冲击波场分布特性, 运用AUTODYN软件对不同运动速度炸药近地爆炸 过程进行数值仿真, 地面超压峰值仿真结果与等效经验公式计算值符合程度较好。通过对扰动压力场、 爆炸 初始压力场以及地面特征位置峰值超压的对比分析, 得到炸药运动速度对近地爆炸冲击波场分布的影响规 律。研究结果表明 炸药运动速度对地面冲击波峰值超压有增益作用, 随着炸药运动速度的增加该增益作用 增强; 炸药运动速度对地面峰值超压的增益作用随比例距离的增加逐渐减弱,3倍比例距离外可认为无明显 差异。 关键词 近地爆炸;运动炸药;冲击波场;爆炸力学;数值计算 中图分类号 TJ55 文献标识码 A 文章编号 1001 -487X(2018)03 -0049 -06 Characteristics Investigation on Shock Wave Field of Near Ground Blasting of Moving Explosive Charge CHANG Bo,GU Hong-ping,NIU Chen-wei,LV Yong-zhu,JING Yi-dong (Xi′an Modern Chemistry Research Institute,Xi′an 710065,China) Abstract To study the characteristics of near ground blasting shock wave field of explosive charge moving at dif- ferent velocities,the simulations were carried out by using AUTODYN software. The simulation results of peak over- pressure near ground was in good agreement with the computing results obtained by equivalent empirical equation. Then the disturbance pressure field,initial explosive pressure field,peak overpressure at characteristic ground loca- tion were compared and analyzed,and the influencing rules of moving velocity on distribution of shock wave near ground were obtained. The simulation results showed that the charge moving velocity rted positive influence on the peak near ground overpressure,meanwhile the positive influence increased along with the increasing of charge moving velocity. This positive influence decreased along with the increasing of characteristic ground position,and no differ- ence can be found at the position for more than triple scale distance. Key words blast near ground;moving explosive charge;shock wave field;explosive mechanics;numerical simulation 收稿日期2018 -04 -10 作者简介畅 博(1990 -) , 男, 工程师、 硕士, 主要从事战斗部设计 与仿真科研工作, (E-mail)changbo204s@ sina. com。 基金项目装备发展部预先研究项目 炸药爆炸过程瞬时(10 -6s量级)释放出巨大的 能量, 转变为高温(103K量级)和高压(1010Pa)的 爆轰产物[ 1], 产生如冲击波超压、 火球热效应等各 种目标破坏效应, 其中冲击波超压是破坏时间最长, 破坏区域最广, 针对人员、 设备和防护设施最具杀伤 力的攻击方式。自从二战初期人们发现由于地面介 质的反射作用, 炸药近地爆炸有着比空中爆炸更强 的破坏效应, 近地爆炸冲击波特性研究成为了常规 武器弹药有效评估手段, 能够为武器弹药的设计开 发提供技术数据支撑, 因此研究炸药近地爆炸典型 条件的冲击波能量输出结构至关重要。 杜红棉等人理论分析了近地爆炸冲击波场传播 万方数据 特性, 并通过仿真和试验验证了冲击波传播特性理 论分析的正确性[ 2]; 刘伟等人进行了近地面 TNT球 形和圆柱形炸药爆炸冲击波场的试验和数值仿真研 究, 两者符合程度较好[ 3]; 汪维等人对 1000 kg TNT 炸药在一定地面高度爆炸冲击波场分布进行二维数 值仿真, 得到了近场与远场峰值超压与炸高的影响 规律,并预测了冲击波超压峰值最大值出现位 置[ 4]。综上可看出近年来国内研究重点多集中于 炸药近地爆炸冲击波场特性分析, 但关于运动炸药 近地爆炸冲击波场特性研究内容较少。随着高速武 器系统的不断发展, 武器设计和研发人员越来越关 注战斗部在高速状态下爆炸毁伤效果。运动炸药近 地爆炸试验研究由于费用和测试手段的限制, 困难 程度较高, 往往也存在试验控制和关键信息捕捉等 问题, 难以对整体过程有全面深入认识; 同时通过纯 理论方法得出问题的完整解析解非常困难, 因此采 用数值仿真方法进行研究具有很大的优越性。 利用AUTODYN显式非线性动力学有限差分软 件, 重点对炸药不同运动速度近地爆炸地面冲击波 分布特性进行研究, 得到炸药运动速度对近地爆炸 冲击波场分布的影响规律, 可为爆破战斗部设计提 供一定的数据参考。 1 近地爆炸冲击波传播特性 当炸药TNT当量和爆炸炸高满足条件H/ 3 ω≤ 0. 35时, 就可视为近地爆炸[ 5]。当爆炸冲击波传播 达到刚性地面时, 地面位置处受压缩空气质点不断 聚集, 密度和压力急剧上升形成反射冲击波。由于 前向入射波的压缩作用, 入射波后的空气介质压力 和温度升高, 因此经地面反射的波阵面速度大于入 射波阵面速度, 不断追赶前向初始冲击波( 未经地面 反射的入射波) 阵面。当入射冲击波与地面夹角达到 某一临界值( 马赫反射临界角) 时, 地面反射波追赶上 前向初始冲击波形成马赫反射。随着地面位置到爆 心距离的不断增加, 初始冲击波与地面反射波两者交 点逐渐远离地面, 马赫波阵面高度不断升高[ 6,7]。近 地爆炸冲击波传播过程示意图如图1所示。 根据初始冲击波与地面夹角的不同范围, 可将 近地爆炸冲击波反射过程分为正反射、 规则反射和 马赫反射三个阶段[ 2,8]。其中正反射和规则反射区 域距离爆心很近, 该区域地面冲击波超压理论计算 过程非常复杂, 同时试验过程中传感器不易于采集 冲击波压力, 因此工程试验一般将近地爆炸马赫反 射区域作为重点关注区域。 图1 近地爆炸冲击波传播过程示意图 Fig. 1 Schematic diagram of shock wave dispersion for blast near ground 2 数值仿真 2. 1 仿真模型 炸药近地爆炸过程仿真计算模型采用二维轴对 称模型, 计算域由炸药和空气域两部分组成, 炸药为 长径比1. 5、 圆柱形的30 kg TNT炸药, 空气域尺寸 为10 m 6. 25 m, 为了简化计算将空气域下边界即 地面处理为零位移刚性地面, 除刚性地面及对称边 界外其它边界设置为Flow-Out自由流出边界。计 算模型采用Euler 2D Multi-Material二维多物质欧拉 法进行计算, 网格尺寸为1. 25 cm 1. 25 cm, 使用 mm-mg-ms单位制进行建模, 仿真计算模型如图2 所示。为研究炸药近地爆炸冲击波场的分布特性, 在地面位置距炸药中心Z = 1至Z = 3区间内均匀 设置1 ~5共五个高斯点, 均位于马赫反射区内; 炸 药起爆方式为两端点起爆;炸药中心距地约为 1. 2 m, 满足近地爆炸的炸高和当量条件。 结合所建立的运动炸药近地爆炸仿真模型, 为 方便表述对近地爆炸冲击波场分布的影响, 规定炸 药轴向Y轴负向(0方向)为炸药运动前向, 反之炸 药轴向Y轴正向(180方向)为炸药运动后向, 炸药 侧向为90方向。定义比例距离Z = R/ 3 W, 式中R 为地面测点位置到炸药中心之间的距离,W为炸药 的等效TNT药量[ 9]。 2. 2 材料模型参数 TNT炸药和空气分别采用JWL状态方程和I- deal Gas状态方程进行描述, 两者状态方程参数如 表1和表2所示[ 10,11]。 p = A 1 - ω R1 V e-R1V+ B 1 - ω R2 V e-R2V+ ωE 0 V 式中A、B、R1、R2、ω为状态方程独立参数;P、 E0、V分别为爆轰产物压力、 炸药比内能和比容。 05爆 破 2018年9月 万方数据 P = A(γ - 1) ρe 式中γ为理想气体常数;ρ为空气初始密度;e 为空气比内能。 图2 数值仿真模型示意图 Fig. 2 Schematic diagram of the numerical model 表1 TNT炸药状态方程参数 Table 1 Material parameters of TNT explosive A/ GPaB/ GPaR1R2E0/ GPaωV 3. 73773. 7474. 150. 90. 060. 31 表2 空气状态方程参数 Table 2 Air material parameters ρ/ (kgm -3) γe/(Jkg -1) 1. 2251. 42. 068E5 3 结果分析 3. 1 不同速度炸药爆炸前扰动压力场 飞行状态下的战斗部炸药在前向产生压缩波, 当亚音速状态时, 产生的是弱的压缩波, 当超音速飞 行状态, 前向产生的是一系列的压缩波, 近地爆炸冲 击波形成的本质是压缩波的聚集、 反射与叠加。在 仿真研究炸药运动速度对近地爆炸冲击波场影响 时, 设置1 ms炸药起爆延时时间, 使得空气中产生 较为稳定的流场, 与真实试验工况相接近, 提高仿真 计算结果的可信度。 通过仿真计算得到的340 m/ s(1 Ma) 、680 m/ s (2 Ma) 、1020 ms(3 Ma)运动速度下炸药运动形成 的初始压力场分布云图, 如图3所示。由图3可以 看出 (1) 在炸药未爆轰前跨音速和超音速运动时, 炸药运动前向产生由无数较强压缩波叠加形成的边 界波 激波。激波波阵面前后空气介质的物理性 质发生突变, 波后空气介质压强急剧升高, 密度和温 度也随之升高。而炸药运动的后向和侧向由于炸药 高速运动产生空化, 周围空气随炸药的运动流动匀 化出现稀疏波, 总的来说由于炸药运动产生的初速 压力场分布较为复杂; (2)随着炸药运动速度的增 加, 炸药运动前向产生的激波和侧向及后向产生的 稀疏波扰动影响区域愈发扁平贴近炸药表面, 这是 因为跨音速和超音速运动时, 炸药运动形成的扰动 影响仅能在马赫锥面范围内, 而被扰动的马赫锥面 角度随扰动速度的增加而减小[ 12]。 图3 炸药爆炸前压力云图 Fig. 3 Pressure nephogram before explosion 3. 2 不同速度炸药近地爆炸初始压力场 运动速度为0 m/ s、340 m/ s(1 Ma) 、680 m/ s (2 Ma) 、1020 ms(3 Ma) 炸药起爆0.2 ms后( 冲击波前 沿尚未到达地面) 冲击波压力分布云图如图4所示。 图4 不同运动速度炸药爆炸压力分布云图 Fig. 4 Pressure nephogram of explosive charge moving at different velocities (1) 研究对象为两端面点起爆的柱形TNT炸 药, 由于两点起爆导致在炸药中心位置发生前后端 面冲击波的叠加作用, 炸药侧面90附近出现局部 高压区, 同时在炸药的两端面附近由于稀疏波作用 出现局部低压区, 整体冲击波场呈现近球形的菱形 波形式扩展。 (2)炸药静爆条件下冲击波场为规则对称结 构, 而炸药动爆时冲击波场在炸药运动前向0附近 出现静爆时未曾出现的局部高压区, 炸药后向180 附近稀疏波局部低压区范围扩大,同时炸药侧面 90附近的局部高压区范围向着炸药运动前向移动。 (3) 与静爆冲击波场比较, 炸药动爆时炸药运 动前向(0 ~ 90范围)冲击波压力增加, 其中0方 向冲击波压力增加幅度最大,而运动速度的反向 (90 ~180范围) 冲击波压力小于静爆冲击波压力, 180方向为压力最小点, 小于静爆冲击波压力。 根据动爆冲击波场特性建立动爆冲击波场简化 模型, 联立Hugoniot方程和线性动量守恒方程, 假设 冲击波速度可由静爆冲击波速度与炸药运动速度的 矢量和求得, 分析后得出影响动爆初始冲击波场分 布的主要因素有炸药运动速度V、 特征位置与炸药 15第35卷 第3期 畅 博, 谷鸿平, 牛晨伟, 等 运动炸药近地爆炸冲击波场特性研究 万方数据 中心夹角θ、 特征位置的比例距离Z等[ 13]。通过分 析不同角度方向超压场分布认为, 当特征位置与运 动前向夹角在0 ~ 90范围内时, 炸药运动速度对 冲击波峰值超压起增强作用, 冲击波到达时间相应 地缩短, 同时在该范围内冲击波超压增量随与运动 前向夹角角度的增加而逐渐减小; 而当特征位置与 运动前向夹角在90 ~ 180范围内时, 运动炸药冲 击波峰值超压小于静爆时冲击波峰值超压。 图5和图6为不同运动速度炸药爆炸前向和后 向冲击波前沿运动趋势图, 结合图4不同运动速度 炸药近地爆炸初始压力分布云图, 可以看出 静爆时 初始冲击波场为规则对称分布, 而动爆初始冲击波 场相较于静爆时前向冲击波前沿外凸, 且随着炸药 运动速度的增加, 前向冲击波前沿相对外凸程度增 大; 相反地, 动爆时后向冲击波场相较于静爆时冲击 波前沿内凹, 且随着炸药运动速度的增加, 后向冲击 波前沿相对内凹程度增大; 动爆初始压力场后向冲 击波前沿相对内凹程度小于前向冲击波前沿相对外 凸程度, 总体随着炸药运动速度的增加, 炸药近地爆 炸初始冲击波场向着地面方向运动。 图5 炸药爆炸前向冲击波运动对比直方图 Fig. 5 Comparison histogram of blast forward-facing shock wave motion 图6 炸药爆炸后向冲击波运动对比直方图 Fig. 6 Comparison histogram of blast backward-facing shock wave motion 3. 3 不同运动速度炸药近地爆炸地面超压分布 由于运动炸药近地爆炸试验的实施难度大, 目 前较完整的公开试验数据较少。选用国防工程设 计规范 中使用的萨道夫斯基公式计算空中静爆冲 击波超压峰值, 其中冲击波峰值超压单位为MPA。 Δp = 0. 085 Z + 0. 30 Z2 + 0. 8 Z3 ,1 ≤ Z ≤ 10 由终点效应理论可知 运动炸药爆炸能量要比 静爆大得多, 爆炸冲击波作用场产生明显变化。当 装药运动的方向与爆炸产物飞散的方向一致时爆炸 效应最大, 并且随两者速度矢量的夹角的增加而减 小[ 14]。根据能量相似原理, 可将炸药运动动能引起 的能量增加视为炸药装药量增加。炸药近地爆炸过 程既不能视作无限空气中爆炸处理, 同时地面在近 地爆炸过程中受到高温高压爆轰产物作用产生破 坏, 消耗一部分能量, 也不能视作完全刚性地面, 因 此当取ωe=(1. 7 ~1. 8)ω时,TNT空中静爆经验预 估公式可用于TNT炸药近地爆炸地面峰值超压近 似计算。综上运动炸药近地爆炸地面超压等效经验 公式为 ωDe= 1. 8 ωse= 1. 8 Qm+ 1 2 V2 Qm ω ZDe= R ωDe Δp De = 0. 085 ZDe + 0. 30 Z2 De + 0. 8 Z3 De 式中Qm为TNT炸药单位质量能量密度;ω为 TNT炸药质量;V为TNT炸药运动速度;R为地面位 置距离炸药爆心的水平距离; Δp De为运动炸药近地 爆炸地面峰值超压等效药量经验公式计算值。 由表3可知, 除个别点以外, 运动炸药近地爆炸 地面峰值超压仿真结果与等效药量经验公式符合程 度较好, 相对误差小于20%。考虑到经验公式的本 身误差以及等效药量的近似简化等因素, 对比结果 表明所使用的仿真计算模型方法在一定范围内具有 较好的适用性。 由图7可以看出, 在炸药药量和起爆方式都相 同的前提下, 不同运动速度炸药近地爆炸地面峰值 超压随着比例距离的增加不断减小, 但随着炸药运 动速度的增加, 地面冲击波峰值超压衰减趋势变陡。 比例距离Z = 1 ~ 3范围内,当炸药运动速度由 0 m/ s上升到1020 m/ s(3 Ma) , 地面冲击波峰值超 压最大下降量分别为1. 54 MPa、1. 94 MPa、 2. 44 MPa和2. 98MPa; 比例距离Z =1 ~3范围内地 面冲击波超压由于远场冲击波场的匀化作用, 不同 运动速度(0 ~3 Ma) 炸药近地爆炸峰值超压增益幅 度由Z =1. 5时的46. 4%下降至Z =3时的9. 8%, 25爆 破 2018年9月 万方数据 可以预见随着比例距离的增加, 炸药运动速度对近地 爆炸冲击波场分布影响将逐渐减小直至无明显差异。 表3 仿真结果与经验公式对比 Table 3 Comparison of results for simulation and empirical equation VelocityZSimulation Empirical equation Compare with empirical equation/ % 1. 01. 6811. 987-18. 22 1. 50. 7320. 6935. 34 0 m/ s2. 00. 3690. 3437. 13 2. 50. 2350. 20512. 96 3. 00. 1450. 1375. 43 1. 02. 0842. 0153. 30 1. 50. 8260. 70215. 01 340 m/ s2.00.3890.34710.81 2.50.2220.2076.78 3.00.1450.1440.84 1.02.5962.09819.17 1.50.9380.72922.26 680 m/ s2.00.4250.36015.38 2.50.2390.21410.41 3.00.1520.1435.77 1.03.1432.23628.85 1.51.0730.77427.84 1020 m/ s2.00.4640.38117.98 2.50.2610.22613.44 3.00.1650.1518.65 图7 炸药近地爆炸地面峰值超压曲线 Fig. 7 Peak overpressure curves of explosion charge blast near ground 分析认为, 随着炸药运动速度的增加, 爆炸初始 压力场近地一侧( 与炸药运动前向夹角0 ~90)压 力值升高, 在相同地面测点处入射冲击波入射角相 同的前提下, 相同地面测点处马赫反射冲击波峰值 超压值随炸药运动速度增加而升高; 同时运动炸药 的近地爆炸初始压力场呈现不规则分布特性, 在与 炸药运动前向夹角0 ~ 90范围内, 随着角度增加 爆炸初始压力值相对于静爆初始压力值增量逐渐减 小; 而地面马赫反射区较远测点位置的入射冲击波, 对应于爆炸初始压力场与炸药运动前向夹角较大的 区域, 且当地面测点比例距离增加时入射冲击波入 射角逐渐减小。因此随着地面测点比例距离的增 加, 不同炸药运动速度引起的地面测点峰值超压增 益值逐渐减小。 运动炸药近地爆炸冲击波场分布特性由炸药运 动速度和炸药炸高等因素决定, 本研究由于仿真计 算量的限制, 重点就运动速度对近地爆炸冲击波场 分布影响进行研究, 未涉及炸高等因素, 运动炸药近 地爆炸冲击波场分布特性相似律还需进一步深入 研究。 4 小结 (1) 运动炸药近地爆炸与静止炸药近地爆炸冲 击波场分布不同, 炸药运动速度对地面冲击波峰值 超压有增益作用, 随着炸药运动速度的增加该增益 作用增强; (2) 由于运动炸药近地爆炸初始压力场的分布 特性以及远场冲击波匀化作用, 炸药运动速度对地 面峰值超压的增益作用随着比例距离的增加逐渐减 弱,3倍比例距离外可认为无明显差异。 参考文献(References) [1] 张俊秀, 刘光烈.爆炸及其应用作用[M].北京 兵器 工业出版社,1998. [2] 杜红棉, 曹学友, 何志文, 等.近地爆炸空中和地面冲 击波特性分析和验证[J].弹箭与制导学报,2014, 34(4) 65-68. [2] DU Hong-mian,CAO Xue-you,HE Zhi-wen,et al. Analy- sis and validation for characteristics of air and ground shock wave near field explosion[J]. Journal of Projec- tiles,Rockets,Missiles and Guidance,2014,34(4) 65- 68.(in Chinese) [3] 刘 伟, 郑毅, 秦 飞.近地面TNT爆炸的试验研究和 数值模拟[J].爆破,2012,29(1) 5-9. [3] LIU Wei,ZHENG Yi,QIN Fei. Experimental and numeri- cal simulation of TNT explosion on the ground[J]. Blas- ting,2012,29(1) 5-9.(in Chinese) [4] 汪 维, 张 舵, 卢芳云, 等.大当量TNT装药爆轰的 远场数值模拟及超压预测[J].弹箭与制导学报, 2010,30(1) 127-130. [4] WANG Wei,ZHANG Duo,LU Fang-y un,et al. Simula- tion of thousands kilograms of tnt explosion shock wave in air at distance far away and forecast of the peak overpres- 35第35卷 第3期 畅 博, 谷鸿平, 牛晨伟, 等 运动炸药近地爆炸冲击波场特性研究 万方数据 sure[J]. Journal of Projectiles,Rockets,Missiles and Guidance,2010,30(1) 127-130.(in Chinese) [5] 赵蓓蕾, 崔村燕, 陈景鹏, 等.近地爆炸地面冲击波传 播规律的数值研究[J].四川兵工学报,2015,36(9) 45-48. [5] ZHAO Bei-lei,CUI Cun-yan,CHEN Jing-peng,et al. Nu- merical study of propagation law of ground shock wave in near surface explosion[J]. Journal of Sichuan Ordnance, 2015,36(9) 45-48.(in Chinese) [6] 王伟力.曲面冲击波马赫反射特性理论研究[J].稀有 金属材料与工程,1992,21(5) 27-30. [6] WANG Wei-li. The theoretical study on the characteristic of mach reflection of curve shock wave[J]. Rare Metal Materials and Engineering,1992,21(5) 27-30.(in Chi- nese) [7] 辛春亮.高能炸药爆炸能量输出结构的数值分析[D]. 北京 北京理工大学,2008. [7] XIN Chun-liang. Numerical analysis of explosion energy output structure for high energy explosives[D]. Beijing Beijing Institute of Technology,2008.(in Chinese) [8] 宁建国, 王 成, 马天宝.爆炸与冲击动力学[M].北 京 国防工业出版社,2010. [9] 谷鸿平, 柳 雯, 李广嘉.起爆方式对柱形炸药空爆冲 击波场的影响[J].爆破,2016,33(4) 34-38. [9] GU Hong-ping,LIU Wen,LI Guang-jia. Influences of ini- tiation modes on air blast wave fields produced by cylin- drical explosives[J]. Blasting,2016,33(4) 34-38.(in Chinese) [10] 杨亚东, 李向东, 王晓鸣.爆炸冲击波空中传播特征 参量的优化拟合[J].爆破器材,2014,43(1) 13-18. [10] YANG Ya-dong,LI Xiang-dong,WANG Xiao-ming. Op- timum fitting for characteristic parameters of blast shock- waves traveling in air[J]. Explosive Materials,2014, 43(1) 13-18.(in Chinese) [11] 时党勇, 李裕春, 张胜民.基于ANSYS/ LS-DYNA 8. 1 进行显示动力分析[M].北京清华大学出版社, 2005. [12] 陈新华, 聂万胜.液体推进剂爆炸危害性评估方法及 应用[M].北京 国防工业出版社,2005. [13] 蒋海燕, 李芝绒, 张玉磊, 等.运动装药空中爆炸冲击 波特性研究研究[J].高压物理学报,2017,31(3) 286-294. [13] JIANG Hai-yan,LI Zhi-rong,ZHANG Yu-lei,et al. Char- acteristics of air blast wave field for explosive charge moving at different velocities[J]. Chinese Journal of High Pressure Physics,2017,31(3) 286-294.(in Chi- nese) [14] 张光莹, 周 旭, 黄咏政, 等.动爆冲击波特性分析方 法研究[C]∥第四届全国计算爆炸力学会议论文集, 2008282-287. [14] ZHANG Guang-ying,ZHOU Xu,HUANG Yong-zheng,et al. The study of analysis of the shock-wave char- acteristic of moving explosive[C]∥The 4th National Symposium onComputationalExplosionMechanics, 2008282-287.(in Chinese) (上接第20页) [14] 闫长斌.爆破作用下岩体累积损伤效应及其稳定性 研究[D].长沙 中南大学,2006. [14] YAN Chang-bin. Study on cumulative damage effects and stability of rock mass under blasting loading[D]. Chang- shaCentral South University,2006.(in Chinese) [15] 肖建清, 冯夏庭, 林大能.爆破循环对围岩松动圈的 影响[J].岩石力学与工程学报,2010,29(11) 2248- 2255. [15] XIAO Jian-qing,FENG Xia-ting,LIN Da-neng. Influ- ence of blasting round on excavation damaged zone of surrounding rock[J]. Chinese Journal of Rock Mechan- ics and Engineering,2010,29(11) 2248-2255.(in Chi- nese) [16] HOEK E,BROWN E T. Underground excavations in rock [M]. Hertford,EnglandAustin &Sons Ltd,1980. [17] HOEK E,CARRANZA-TORRES C,CORKUM B. Hoek- Brown failure criterion[C]∥Proceedings of NARMS- Tac[S. l.] [s. n.] ,2002267-273. [18] 朱传云, 喻胜春.爆破引起岩体损伤的判别方法研究 [J].工程爆破,2001,7(1) 12-16. [18] ZHU Chuan-yun,YU Sheng-chun. Study on the criterion of rock mass damage caused by blasting[J]. Engineering Blasting,2001,7(1) 12-16.(in Chinese) [19] BARTON N. Some new Q-value correlations to assist in site characterization and tunnel design[J]. International Journal of Rock Mechanics and Mining Sciences,2002, 39(2) 185-216. [20] 向明智.金沙江白鹤滩电站工程岩体质量分类[D]. 成都 成都理工大学,2009. [20] XIANG Ming-zhi. Engineering rock mass classification of Baihetan hydroelectric station at Jin Sha River[D]. ChengduChengdu University of Technology,2009.(in Chinese) 45爆 破 2018年9月 万方数据