Experiment of explosive consumption by blasting pretreated boulders with overlying stratum of rock-soil
-
摘要: 盾构机在风化花岗岩地层掘进会遇到高强度孤石,为降低施工风险、避免刀盘及刀具严重损坏,常采用地面钻孔爆破对前方孤石进行预处理,但其炸药单耗一般由经验公式确定,导致爆破效果不理想。结合实际工程对爆破块度小于30 cm的特殊要求,采用爆破模型实验对炸药单耗进行研究。实验结果表明:满足特殊块度要求的陆地钻孔爆破炸药单耗是常规钻孔爆破的5.4~6.5倍;上覆地层厚度与炸药单耗呈线性递增关系;实验得出的炸药单耗是瑞典经验公式计算值的3.4~4.9倍。由此提出满足工程需要的炸药单耗修正公式,应用于实际工程并通过验证取得良好效果,为类似工程孤石爆破预处理炸药单耗的计算提供依据和参考。Abstract: During the tunneling construction in hard granite strata utilizing shield machines, pretreatment of boulders by blasting was commonly used to prevent the abrading cutter from serious damages and reduce construction risks, but the explosive consumption was usually determined based on empirical formula, resulting in that the blasting effect was often far from satisfactory. In the present work, aiming to meet the special requirement that rock fragments from the blasting should be less than 30 cm in practice, we investigate the explosive consumption using a model test of blasting. The results show that the explosive consumption of land blasting is 5.4-6.5 times that of conventional blasting when the requirement of rock fragmentation is satisfied; the explosive consumption increases linearly with the depth of overlying stratum; and the explosive consumption of model test was 3.4-4.9 times that of the empirical formula. Thus, the revised empirical formula is derived and proved valid based on the on-the-site practice. The revised formula can serve as calculation basis and reference value for similar projects of tunneling construction.
-
Key words:
- mechanics of explosion /
- explosive consumption /
- blasting /
- boulder /
- shield /
- pre-treatment /
- model experiment
-
应变增长现象是爆炸容器响应研究中的一个重要问题,主要是指容器壳体的最大变形不是出现在载荷最强的第1个周期内,而是在稍后某时刻出现的现象。在应变增长现象中,容器出现超出预期的变形,对容器安全形成威胁。现有研究表明,导致应变增长的原因主要有2个方面。首先是容器内部的周期性载荷与壳体响应频率接近或为壳体振动频率的整数倍时,可能引发共振并导致应变增长[1-2]。对于共振引起的应变增长现象,若周期性载荷峰值的到时位于应变曲线上升沿前段,所引发的应变增长现象更加严重[3]。另外,容器上法兰、开孔等结构会引起壳体的弯曲振动,弯曲振动与壳体呼吸振动的叠加也是导致应变增长的重要原因[4-7]。
为满足装药、测试等需求,容器必须开设开口,开口依靠螺栓联结。螺栓的响应状态直接关系容器开口是否安全以及密封效果。现有关于应变增长现象的研究主要针对容器壳体,很少涉及容器开口的螺栓。Duffey等研究了球形容器开口螺栓预紧力设计方法[8],数值模拟结果显示,螺栓应力曲线出现了应变增长现象[9],但该现象未引起重视,后续也没有进一步的实验以及深入研究。实际上,螺栓的应变增长会导致应力提高、密封面间隙增大,从螺栓强度、密封性两方面对容器使用产生影响。
本文基于球形爆炸容器开口的螺栓响应实验数据探索应变增长现象;通过数值模拟,分析容器开口处爆炸载荷增强现象、周期性载荷与端盖螺栓结构响应共振问题,讨论端盖质量、螺栓预紧力对螺栓应变增长的影响和相应的避免方法。
1. 螺栓应变增长现象实验结果
实验中使用的球形爆炸容器如图1所示,容器内径523 mm,壁厚20 mm,材料为Q345R。容器赤道线上方均布4个压力传感器安装孔。装药开口位于容器上方,开口接管内径120 mm,壁厚30 mm。为保证实验中端盖、螺栓始终处于弹性范围内,端盖厚度设计为30 mm,其内表面距容器中心405 mm;螺栓选用16根M16的高强螺栓。实验采用中心球形装药,主装药成分为60RDX/40TNT,爆炸当量为670 g TNT。实验主装药的中心装有微米级PETN粉压制而成、当量为1 g TNT的起爆药球。
实验中以螺栓应变为对象研究开口结构的响应规律。选取均布的4根螺栓,按照图2的方法,沿螺栓轴向对称粘贴2个应变片,应变片上沿与螺栓六方头下沿距离为10 mm;螺栓拧紧时注意调整2个应变片的位置,使其分别朝向半径内侧和外侧。实验容器法兰和端盖上的螺栓孔直径均为18 mm,螺栓安装时注意使应变片与螺栓孔内侧的间隙均匀,避免实验中应变片与螺栓孔碰撞受损。螺栓孔外侧设有引线槽,便于应变片的导线引出。实验前,为保证螺栓预紧力一致,使用力矩扳手拧紧螺栓,设定最大力矩为170 N∙m,预紧后螺栓轴向拉应力约为300 MPa[10]。
图3为4根螺栓内侧(inside)、外侧(outside)应变曲线的原始数据,其中螺栓内侧应变曲线上存在3个较明显的峰值,螺栓外侧应变曲线变化平缓,说明螺栓在响应中既发生了沿螺栓轴向的拉伸变形,也发生了沿法兰半径方向的弯曲变形。其原因为爆炸载荷作用于端盖的中心区域并导致端盖弯曲变形,使端盖的中心挠度大、边缘挠度小;进一步地,螺栓在端盖的带动下发生相同趋势的弯曲变形。
对每根螺栓内侧、外侧的应变曲线求平均值,可以得到图3中螺栓中轴线(axial)处的应变曲线,对应螺栓在结构响应中的轴向拉伸变形。螺栓中轴线应变的第2个峰值高于第1个峰,且为全部振动过程中的最大值,说明螺栓发生应变增长;同样的现象也出现在螺栓内侧的应变曲线上。
应变增长系数是分析应变增长严重程度的常用量化标准,为应变曲线最大值与第1个峰值的比值。表1中统计了4根螺栓中轴线和内侧应变曲线的最大值(第2个峰值)和第1个峰值,即应变曲线的前两个峰值,并计算了应变增长系数。螺栓内侧和中轴线处的最大应变增长系数分别为1.83和2.10,平均应变增长系数分别为1.47和1.60。
表 1 螺栓应变统计Table 1. Statics of bolt strain编号 第1应变峰值/10−6 第2应变峰值/10−6 应变增长系数 内侧 中轴线 内侧 中轴线 内侧 中轴线 1 762 373 1 391 691 1.83 1.85 2 900 371 1 357 692 1.51 1.87 3 1 057 610 1 148 615 1.09 1.01 4 869 329 1 367 690 1.57 2.10 平均 897 421 1 316 672 1.47 1.60 螺栓发生应变增长现象,既可能导致强度隐患,也可能导致密封面间隙超出预期。因此有必要开展针对性的研究,找出引发螺栓应变增长的原因并在设计中避免。
2. 数值计算模型
为研究螺栓应变增长现象的原因,需以实验数据为基础,通过数值模拟进行深入分析。为简化计算模型、节约计算时间,本文的数值模拟工作分为两步进行:首先通过二维计算模型确定端盖受到的冲击波载荷历程,然后再通过三维模型求解端盖和螺栓的动态响应。
图4为使用AUTODYN软件[11]建立的压力载荷计算模型。模型根据容器尺寸建立,其中空气和炸药采用欧拉网格,用于模拟爆炸流场的演化过程;钢制容器采用拉格朗日网格,通过流固耦合法分析爆炸载荷与容器间的作用。模型中空气采用理想气体状态方程,初始密度ρ0=1.29 kg/m3,初始比内能e0=2.068×105 J/kg,绝热指数γ=1.4。炸药采用JWL状态方程描述:
p=A(1−ω′R1V)exp(−R1V)+B(1−ω′R2V)exp(−R2V)+ω′e (1) 式中:p为压力,V为相对体积,A、B、R1、R2和ω'为状态方程参数。表2给出了炸药参数,其中,D为爆速,pCJ为炸药CJ爆轰压力。
表 2 JWL状态方程参数Table 2. Parameters in JWL Equation of Stateρ0/(kg∙m−3) e0/(J∙kg−1) D/(m∙s−1) pCJ/GPa A/GPa B/GPa R1 R2 ω' 1 630 4.26×106 6 930 21 371.2 3.231 4.15 0.95 0.3 图5为使用LS-DYNA软件[12]建立的法兰响应计算模型,模型尺寸与实验容器完全相同。为减少计算时间,将模型建立为1/16圆周对称模型;并按图添加法兰下方的固定边界、法兰和端盖的旋转对称边界;考虑法兰、端盖和螺栓之间的接触,静、动摩擦系数设置为0.1。由于动态响应中,雷管罩主要通过质量惯性效应影响端盖运动,将雷管罩等结构简化为图示的柱形结构。柱形结构与端盖连接为一个整体,其外径与雷管罩外径一致。模型中的全部结构采用线弹性模型,材料密度为7.83 g/cm3,弹性模量为200 GPa,泊松比为0.3。求解前,首先对螺栓施加300 MPa的轴向预紧力[13];然后将图4模型算得的端盖中心压力载荷施加在端盖内表面,求解动态响应过程。
3. 容器开口处的载荷增强现象
研究螺栓应变增长现象,应首先分析端盖处的载荷特性。对于实验容器,由于端盖中心有雷管罩、雷管、导爆索等结构,不具备安装压力传感器的条件,因此本文采用图4的计算模型,间接获得端盖处载荷历程。将图4模型得到的容器内壁载荷历程与实验中使用压电传感器测得的载荷历程进行比较,如图6所示。实测得到的压力载荷峰值约为88.2 MPa,与数值模拟结果基本一致。载荷峰值过后,实验测得的压力载荷曲线逐渐衰减到负值;而实际上载荷峰值后,受到温度升高和爆炸产物的影响,容器内压力应与准静态压力相当[13],故载荷曲线的变化趋势应与数值模拟曲线相似。因此认为载荷第1个脉冲峰值后的曲线为异常数据,其原因可能是传感器中压电晶体在冲击波作用下强烈振动,导致传感器损坏。综上,通过数值模拟得到的载荷峰值与实验结果基本一致,后续的载荷历程也与实际物理过程相符合,计算模型基本可靠,可用于开展进一步的分析。
使用图4的模型计算端盖(closure)处的载荷历程,并与容器内壁(shell)进行比较,如图7所示。由于端盖的球心距大于容器内壁,因此端盖处载荷第1个峰值到时略晚于容器内壁。同时,受到开口结构的影响,端盖处载荷第1个峰值达到103.8 MPa,比容器内壁高约41%。在后续的载荷曲线上,容器内壁和端盖处均可见冲击波在容器内来回反射形成的周期性脉冲,且脉冲峰值到时能够相互对应。但与容器内壁载荷不同,受到开口结构的影响,端盖处载荷曲线的第奇数个峰的峰值较高,第偶数个峰值较低,因此可认为端盖处载荷的周期约为容器内壁的2倍。按照上述周期计算方法,容器内壁前2个周期的冲量约为10.7 kPa∙s;而端盖处载荷第1个周期的冲量约为13.9 kPa∙s,比容器内壁高约30%。综上,受到开口结构影响,与容器内壁相比,端盖处的载荷峰值、冲量均大幅增强,以下简称为端盖处的载荷增强现象。
为验证通过数值模拟发现的载荷增强现象是否存在,并分析载荷增强现象对螺栓应变增长的影响,将图7中容器内壁压力载荷曲线(pshell)、端盖处的压力载荷曲线(pclosure)分别代入图5的模型中,得到螺栓动态响应曲线如图8所示,图中实验应变曲线为4根螺栓内侧应变曲线(inside strain)、中轴线应变曲线(axial strain)的平均值。由图可见,如果不考虑载荷增强现象,得到的螺栓内侧、中轴线处应变曲线的第1个峰值不足实验数据的1/2,且为动态响应历程中应变的最大值,即螺栓不会发生应变增长现象。当考虑载荷增强现象后,螺栓内侧、中轴线处应变曲线的前两个应变峰值及其变化趋势与实验数据均有较高的一致性;且应变曲线的第2个峰值高于第1个峰值,即在数值计算结果中出现了螺栓变形的应变增长现象。
综上,数值计算结果与实验数据的对比表明:(1)容器开口处的载荷增强现象是真实存在的;(2)载荷增强现象会导致的螺栓变形增大幅度超过100%,也是导致螺栓应变增长的重要原因。
4. 关于螺栓应变增长原因的讨论
如图9所示,从时间历程上看,螺栓应变曲线的前2个峰值分别是在被增强压力载荷的第1周期和第2周期作用下产生的;螺栓应变曲线的第3个峰低于前2个峰,说明螺栓应变增长主要受被增强压力载荷前2个周期的影响。为分析螺栓应变增长与被增强压力载荷第2个周期的关系,在计算中忽略压力载荷第2周期的作用,使端盖、螺栓在压力载荷第1周期结束后自由振动。由图9的计算结果,忽略被增强压力载荷第2个周期后,应变曲线的第2个峰值显著降低且不再发生应变增长现象,说明螺栓应变增长与被增强载荷的第2个周期密切相关。
进一步地,如图10(a)所示,将计算模型中被增强压力载荷的第2个周期到时分别提前100 μs、延后100μs和200 μs,观察端盖和螺栓动态响应的变化。图10(b)为法兰内径处的间隙变化曲线,由图可见,端盖的响应为周期性的“分离”、“合拢”运动。在实验条件下,压力载荷第2个周期到达时,恰处于端盖“合拢”运动结束、“分离”运动开始的时刻;对比载荷和位移曲线的时间历程可知,此时压力载荷第2个周期全部做正功,使得端盖中心位移曲线的第2个峰值高于第1个峰值,并带动螺栓发生应变增长,如图10(c)所示。如果压力载荷第2周期到时提前100 μs、延后100 μs或延后200 μs,载荷第2周期到时将分别位于“合拢”运动后期、“分离”运动中期、“分离”运动结束时刻,作用效果分别是提前“分离”运动开始时间、延长“分离”运动时间、延长“合拢”运动时间;从能量的角度看,上述3种情况的载荷第2周期均有做负功的阶段,因此端盖中心位移有不同程度的降低,图10(c)中的螺栓应变曲线也没有应变增长现象。综上,螺栓发生应变增长还要求被增强压力载荷第2周期的到时与端盖第2周期运动开始的时刻接近,即被增强压力载荷与结构的共振是引发的螺栓应变增长的重要原因。
结构尺寸的变化会改变结构响应频率,从而对共振和螺栓应变增长产生影响。对于端盖和螺栓组成的振动系统,端盖质量是影响振动频率的重要因素。为分析端盖质量对螺栓应变增长的影响,逐渐减小图5的模型中的雷管罩质量进行对比计算,结果如图11所示。图中计算法兰间隙的位置与图10一致,螺栓应变曲线为螺栓轴向的拉伸应变曲线。由图11可见,当雷管罩质量增大,法兰间隙曲线的第1个峰值逐渐降低、第1峰值的到达时间和响应周期逐渐增大。其原因是载荷第2周期的到时逐渐接近位移曲线第2个峰的起始时刻,即逐渐接近共振条件。计算图11中各工况应变曲线第2个峰与第1个峰的比值并记为f,可见随着雷管罩质量增大,f值由0.78增大到1.12,应变增长现象产生并加剧。综上,特定的端盖质量会引发结构响应与压力载荷的共振,并导致螺栓发生应变增长;而改变端盖质量即可避免螺栓应变增长的发生。
预紧力也是影响端盖-螺栓系统响应周期的重要因素。为分析预紧力对螺栓应变增长的影响,改变计算模型中预紧力,得到图12所示的计算结果。预紧力由100 MPa增大到700 MPa,法兰间隙曲线的第1个峰值及其到时、系统响应周期都逐渐减小,同时法兰间隙曲线的第2个峰的起始时间也逐渐提前。当螺栓预紧力为300 MPa、即实验条件下,被增强压力载荷第2周期的到时与法兰间隙曲线第2个峰值的起始时间接近,满足共振形成的条件。由图12应变曲线可知,当螺栓预紧力为100、300、500、700 MPa时,螺栓的应变增长系数分别为0.99、1.12、0.94、0.43,即仅当螺栓预紧力为300 MPa时,螺栓发生了应变增长。因此,在特定的预紧力下,结构响应周期可满足与被增强压力载荷共振的条件,并引发螺栓应变增长;在设计中改变预紧力,即可避免螺栓应变增长现象。
综上,螺栓应变增长的主要原因是被增强的压力载荷与结构响应的共振,实验中发现的螺栓应变增长现象,是在被增强载荷的周期性、端盖与雷管罩总质量、螺栓预紧力大小等因素满足共振的情况下发生的。事实上,共振形成的条件是比较苛刻的,只要在设计中根据压力载荷的周期特性,调整端盖质量、预紧力等结构参数,避免结构响应与压力载荷的共振,就可以避免螺栓应变增长现象。
5. 结 论
(1)通过实验观察到球形爆炸容器螺栓的应变增长现象,实验中螺栓内侧和中轴线处的最大应变增长系数分别为1.83和2.10,平均应变增长系数分别为1.47和1.60。
(2)容器开口处存在载荷增强现象,增强载荷的前两个周期是影响螺栓应变增长的主要因素。
(3)被增强压力载荷与结构响应的共振是引发螺栓应变增长的重要原因,共振发生的条件是被增强压力载荷的第二周期脉冲到时与端盖第二周期运动开始时刻接近;端盖质量、螺栓预紧力等结构参数的变化,可能引起被增强压力载荷与结构响应共振,并导致螺栓的应变增长现象。
-
表 1 模型实验中选取的现场物理量、量纲及数值
Table 1. Dimensions and values of on-the-site physical quantities in model test
编号 符号 名称 单位 量纲 数值 1 Q 炸药量 kg M - 2 v 炸药爆速 m/s LT-1 4 000 3 ρ 炸药密度 kg/m3 ML-3 1 200 4 R 装药半径 m L 0.045 5 h 装药长度 m L 0.90 6 ρr 岩石密度 kg/m3 ML-3 2 800 7 vr 岩石纵波波速 m/s LT-1 4 000 8 σr 岩石强度 N/m2 ML-1T-2 1×106 9 L 岩石尺寸(正方体) m L 1.00 10 D 岩石破碎的平均块度 m L 0.30 11 H 岩石埋深 m L - 12 ρc 上覆层密度 kg/m3 ML-3 2 600 13 vc 上覆层纵波波速 m/s LT-1 2 400 表 2 相似律各参数含义及数值
Table 2. Parameters and values of similarity law
参数 数值 参数含义 π1 0.05 模型和实际的装药尺寸相似比 π2 0.045 模型与实际爆破的尺寸相似比 π3 0.15 模型与实际爆破的爆破块度相似比 π4 0.143 实验模型和几何原型的相似比 π5 0.005 2 炸药的爆轰压力与岩石强度的相似比 π6 0.43 炸药与岩石的波阻抗相似比 π7 1.8 岩石与覆盖层的波阻抗相似比 π8 62.8 装药量与装药半径的相似比 表 3 混凝土试块配合比(kg/m3)
Table 3. Proportions of components in concrete sample (kg/m3)
标号 水泥 砂 碎石(5~10 mm) 碎石(10~20 mm) 粉煤灰 聚羧酸高效减水剂 水 C50 376 659 468 703 94 4.7 145 表 4 试块爆破装药量及爆后块度统计
Table 4. Charge weight of samples and fragmentation statistics
试块 Q/g D/cm q/(kg·m-3) H/cm 实验场地 qb/(kg·m-3) 实验现象 1-1 20 2.00 5.93 0 实验用水池 3.26 在砂层表面形成明显爆坑 1-2 11 4.50 3.26 1-3 8 5.54 2.37 2-1 12 4.63 3.56 30 实验用水池 3.85 覆盖层表面出现了直径约30 cm左右的塌腔 2-2 13 4.54 3.85 2-3 14 4.40 4.15 3-1 14 4.67 4.15 50 实验用水池 4.44 表面出现直径约45 cm的环状裂纹且试块上部出现空腔 3-2 13 5.33 3.85 3-3 15 4.48 4.44 4-1 16 4.46 4.74 70 工地料场 4.74 表面无裂纹,试块周边出现空区 4-2 17 4.28 5.04 4-3 18 3.98 5.33 5-1 15 4.85 4.44 100 工地料场 5.04 表面无裂纹,试块周边出现空区 5-2 18 4.09 5.33 5-3 17 4.45 5.04 6-1 22 3.78 6.52 150 工地料场 5.93 表面无裂纹,试块周边出现空区 6-2 20 4.65 5.93 6-3 19 4.31 5.63 表 5 模型实验和现场爆破的炸药单耗
Table 5. Explosive consumption in model experiment and on-the-tsite blasting
Hm/m Hf/m qb/(kg·m-3) 0 0 3.26 0.3 2.1 3.85 0.5 3.5 4.44 0.7 4.9 4.74 1.0 7.0 5.04 1.5 10.5 5.93 表 6 水下爆破炸药单耗经验公式对比
Table 6. Comparison of explosive consumption empirical formula in underwater explosion
公式来源 表达式 考虑的主要因素 中国水利系统 q水=q陆+0.01H水+0.02H介质+0.03H梯 水深、炸药埋深、梯段高度 日本炸药协会 {Lα=HCα(水压修正)Lβ=H0Cβ(上覆土岩修正) 水深、上覆层厚度 瑞典经验公式 q=q1+q2+q3+q4 水深、梯段高度、覆盖层高度 《工程爆破实用手册》 Q=KWaH[1.45+0.45Exp(-0.33H0/W)] 水深、抵抗线、孔距、梯段高度 表 7 瑞典经验公式和本文拟合公式的炸药单耗对比
Table 7. Comparison of explosive consumption between results of Sweden empirical formula and fitting formula
H/m q/(kg·m-3) η 式(5) 瑞典经验公式 0 1.10 3.39 3.39 5 1.37 4.63 3.38 10 1.52 5.87 3.86 15 1.67 7.11 4.26 20 1.82 8.25 4.53 25 1.97 9.59 4.87 表 8 孤石段与正常段盾构掘进参数对比表
Table 8. Comparison of boulder's driving parameters with those of normal section
掘进段 v/(mm·min-1) t/min F/kN 孤石段 15.5 121 18 739 正常段 13.8 108 18 035 -
[1] 刘建国.深圳地铁盾构隧道技术研究与实践[M].北京:人民交通出版社, 2011. [2] 杜力编译.无自由面爆破破碎岩盘[J].爆破, 1990(4): 66-69. http://www.cqvip.com/qk/91438X/199301/1286481.htmlDu Li Translated. Blasting rock without free surface[J]. Blasting, 1990(4): 66-69. http://www.cqvip.com/qk/91438X/199301/1286481.html [3] 徐朝辉.台山核电站取水隧洞基岩及孤石处理技术[J].陕西水利, 2010(6):81-82. http://d.old.wanfangdata.com.cn/Periodical/shanxsl201006045Xu Zhaohui. Treatment technology of bedrock and boulders in Taishan nuclear intake tunnel[J]. Shaanxi Water Resources, 2010(6):81-82. http://d.old.wanfangdata.com.cn/Periodical/shanxsl201006045 [4] 洪开荣.水下盾构隧道硬岩处理与对接技术[J].隧道建设, 2012, 32(3):361-365. http://d.old.wanfangdata.com.cn/Periodical/sdjs201203021Hong Kairong. Case study on hard rock treatment technology and shield docking technology in boring of underwater tunnels[J]. Tunnel Construction, 2012, 32(3):361-365. http://d.old.wanfangdata.com.cn/Periodical/sdjs201203021 [5] 路耀邦, 刘洪震, 游永锋.海底盾构隧道孤石爆破预处理关键技术[J].现代隧道, 2012, 49(5):117-122. doi: 10.3969/j.issn.1009-6582.2012.05.020Lu Yaobang, Liu Hongzhen, You Yongfeng. Key techniques for the pretreatment of boulder blasting in an under-sea shield-driven tunnel[J]. Modern Tunnelling Technology, 2012.49(5):117-122. doi: 10.3969/j.issn.1009-6582.2012.05.020 [6] 海瑞克股份公司.台山核电引水隧道盾构机说明书[M].海瑞克股份公司, 2010. [7] 刘美山, 吴新霞, 张恒伟.混凝土水下爆破炸药单耗试验分析[J].爆破, 2007, 24(1):10-15. doi: 10.3963/j.issn.1001-487X.2007.01.003Liu Meishan, Wu Xinxia, Zhang Hengwei. Experimental analysis on specific charge of underwater explosion of concrete[J]. Blasting, 2007, 24(1):10-15. doi: 10.3963/j.issn.1001-487X.2007.01.003 [8] 刘慧, 冯叔瑜.炸药单耗对爆破块度分布影响的理论探讨[J].爆炸与冲击, 1997, 17(4):359-362. http://www.bzycj.cn/article/id/10486Liu Hui, Feng Shuyu. Theoretical research of the effect on the blasting fragentation distribution from the explosive specific charge[J]. Explosion and Shock Waves, 1997, 17(4):359-362. http://www.bzycj.cn/article/id/10486 [9] 陈运轩.爆破块度效应对炸药单耗的影响[J].爆破, 1996, 13(3):19-22. http://www.cnki.com.cn/Article/CJFDTOTAL-BOPO603.004.htmChen Yunxuan. Effect on the explosive specific charge considering blasting fragmentation[J]. Blasting, 1996, 13(3):19-22. http://www.cnki.com.cn/Article/CJFDTOTAL-BOPO603.004.htm [10] 马建军, 熊祖钊, 段卫东, 等.工程爆破模拟试验的相似律[J].武汉科技大学学报:自然科学版, 2001, 24(2):170-173. http://d.old.wanfangdata.com.cn/Periodical/whkjdxxb200102020Ma Jianjun, Xiong Zuzhao, Duan Weidong, et al. Similarity law in simulation test of engineering blasting[J]. Jounal of Wuhan University of Science &Technology: Natural Science Edition, 2001, 24(2):170-173. http://d.old.wanfangdata.com.cn/Periodical/whkjdxxb200102020 [11] 杨振声.工程爆破的模型试验与模型律[J].工程爆破, 1995, 1(2):1-15. http://www.cnki.com.cn/Article/CJFDTOTAL-GCBP502.000.htmYang Zhensheng. Modeling experiment of exgineering blasting and model theory[J]. Engineering Blasting, 1995, 1(2):1-15. http://www.cnki.com.cn/Article/CJFDTOTAL-GCBP502.000.htm [12] 秦绍兵, 秦绍红.井下中深孔爆破模型试验的相似性研究[J].工程爆破, 2001, 7(2):9-12. doi: 10.3969/j.issn.1006-7051.2001.02.003Qin Shaobing, Qin Shaohong. Study on the similarity of model experiment of underground medium-length-hole blasting[J]. Engineering Blasting, 2001, 7(2):9-12. doi: 10.3969/j.issn.1006-7051.2001.02.003 [13] 王伟, 李小春.无临空面深层岩体爆破增渗模型试验相似律研究及相似条件的实现[J].岩土力学, 2009, 30(7):1892-1898. doi: 10.3969/j.issn.1000-7598.2009.07.003Wang Wei, Li Xiaochun. Study of similar law of model test of blasting echanced permeability in deep rock mass without free surface and its implementation[J]. Rock and Soil Mechanics, 2009, 30(7):1892-1898. doi: 10.3969/j.issn.1000-7598.2009.07.003 [14] 左启东.模型试验的理论和方法[M].北京.水利水电出版社, 1984. [15] Ma Qinyong, Cai Meifeng. Determination of similarity of explosives for a model experiment[J]. Combustion, Explosion and Shock Waves, 2003, 39(5):606-609. doi: 10.1023/A:1026178222336 [16] 李泉.几种水下钻孔爆破炸药单耗计算公式的分析与比较[J].爆破, 2012, 29(1):94-97. http://d.old.wanfangdata.com.cn/Periodical/bp201201025Li Quan. Analysis and comparision of several calculation formulas of specific charge underwater drilling blasting[J]. Blasting, 2012, 29(1):94-97. http://d.old.wanfangdata.com.cn/Periodical/bp201201025 [17] Arpaz E, Uysal O, Tola Y, et al. Comparison of blasting-induced ground vibration predictors in Seyitomer coal mine[M]. Harmonising Rock Engineering and the Enviroment, 2011, 25(2): 153-175. -