A study on explosive load history of rock blasting considering rock failure zones
-
摘要: 针对岩石爆破爆炸荷载历程中未联合考虑岩石爆破动态过程和炮孔周围岩体破坏分区的不足,开展了考虑岩体破坏分区的岩石爆破爆炸荷载历程及其适用性研究。联合岩石爆破动态过程和岩体破坏分区的理论解,推导了考虑岩体破坏分区的岩石爆破爆炸荷载理论公式,比较了考虑岩体破坏分区的岩石爆破爆炸荷载历程与实测炮孔爆炸压力曲线,开展了单孔爆破现场试验和相应条件下3种爆炸荷载工况的数值模拟,并对爆破振动现场实测和数值模拟结果进行了对比。研究结果表明:考虑岩体破坏分区的爆炸荷载历程包括上升段和衰减段Ⅰ、Ⅱ、Ⅲ,上升段持续时间极短,衰减段持续时间较长且主要由填塞情况控制;考虑岩体破坏分区的爆炸荷载历程理论计算结果与实测炮孔爆炸压力曲线的变化趋势一致,验证了考虑岩体破坏分区的岩石爆破爆炸荷载理论公式的可靠性;考虑岩体破坏分区的爆炸荷载工况下,单孔爆破振动波形的数值模拟结果与现场实测结果的主要特征一致,该荷载工况下质点峰值振速计算结果与现场实测值偏差率最小,绝大部分在7%以内,表明了其应用于数值模拟的优越性;考虑岩体破坏分区的爆炸荷载可随岩石爆破系统条件的变化而动态调整,其可靠性好、适应性强、应用效果佳。Abstract: Due to the deficiency that dynamic processes of rock blasting and rock failure zones around a blasthole are not simultaneously considered, the explosion load history of rock blasting considering rock failure zones and its reliability were investigated. Combining theoretical solutions of the dynamic processes of rock blasting and the rock failure zones around a blasthole, a theoretical formula of the explosive load history considering rock failure zones was derived, and a comparison was made between the derived explosive load history and a measured explosion pressure curve inside a blasthole. Both the field test on an ideal site and the numerical simulation including three explosion load conditions of single hole blasting were carried out, and the field and numerical results of blasting vibration were compared. The results show that the explosive load history considering rock failure zones consists of an ascending stage and three attenuation stages Ⅰ, Ⅱ, and Ⅲ, among which the ascending stage lasts for an extremely short time, while the attenuation stages last for a long time and are controlled by the stemming conditions. The change tendency of the calculated explosive load history considering rock failure zones is consistent with that of the measured explosion pressure curve, indicating the reliability of the explosive load history considering rock failure zones. The numerical results of single hole blasting vibration waveforms under the theoretical explosive load condition are consistent with the filed results, and the deviation ratios between the calculated peak particle velocity (PPV) results under the theoretical explosive load condition and the field PPV results are the smallest, most of which are within 7%, indicting the explosive load history considering rock failure zones has strong reliability. The explosive load history considering rock failure zones can be adjusted as the rock blasting system changes, and it has wide adaptability and good application potentials. The research results may help provide a theoretical basis for realizing efficient and accurate calculation about rock blasting.
-
Key words:
- rock blasting /
- explosive load /
- theoretical formula /
- failure zone /
- rock breaking process
-
近年来空间碎片环境日益恶化,截止到2019年,尺寸为1~10 cm的空间碎片数量约为75万个,尺寸为1~10 mm的碎片数量超过1.5亿个,未来50年空间碎片数量将以每年10%的速度增长[1]。航天器和空间碎片碰撞时的平均相对速度为10 km/s,一旦撞击将导致航天器表面被损伤或破坏,引起部件、分系统甚至整个航天器失效。
为提高航天器在恶劣空间碎片环境中的生存能力,自20世纪80年代以来,NASA、ESA、JAXA等基于传统Whipple防护结构进行了大量的超高速撞击实验,研制出多种高性能防护结构,选择的材料主要有高强度铝合金板、铝网、泡沫铝、蜂窝板、Nextel陶瓷布、Kevlar纤维布等[2-5]。我国航天器空间碎片防护研究工作起步较晚,缺少高性能防护材料,面对空间碎片防护工程需求日趋强烈的现状,一些学者开展了高性能防护材料的研制工作,并得到了初步的工程应用。贾古寨等[6]研究制备了高性能玄武岩纤维织物,并进行了超高速撞击实验研究,获得了其超高速撞击损伤与防护机理;王应德等[7]研究制备了高性能SiC纤维织物,开发了国产高性能SiC纤维织物填充防护结构,并建立了对应撞击极限;侯明强等[8]提出将阻抗梯度材料应用于空间碎片防护,并系统研究了阻抗梯度材料及其防护结构的超高速撞击特性,获得了材料和结构的最优设计方案。尽管航天器空间碎片防护结构材料多样,但均为惰性材料,由于防护机理单一,很大程度上制约了防护结构对大尺寸空间碎片的防护能力。
PTFE/Al是一种冲击引发的活性材料,可在冲击加载下被引发反应释放化学潜能,并对目标造成穿甲、燃烧、内爆等多种形式的综合毁伤效应。学者们围绕活性材料动态力学特性、冲击反应释能特性、高效毁伤战斗部应用等进行了大量的理论分析与实验研究,但均以增强战斗部的毁伤效果为研究背景[9-13]。因此,开展PTFE/Al活性材料超高速撞击实验研究,获取防护结构的损伤与防护特性,对航天器空间碎片防护能力的提高具有重要意义。本文中利用二级轻气炮开展活性材料防护结构的超高速撞击实验,结合防护结构在不同撞击条件下的碎片云与后板损伤特性,探讨活性材料防护性能,拟合得到撞击极限曲线。
1. 超高速撞击实验
1.1 PTFE/Al活性材料薄板的制备
实验中作为防护屏的PTFE/Al是一种零氧配比的活性材料,采用冷压烧结制备工艺,经混合、干燥、压制、烧结而成。其中PTFE的质量分数为73.5%,Al的质量分数为26.5%,表1中给出了主要原材料的相关参数。由于PTFE/Al需要作为空间碎片防护结构的一部分,因此要求材料的力学强度越高越好。通过前期对制备工艺的研究可知,当成型压力为80 MPa时,结合具有熔融、结晶平台的烧结工艺(见图1),获得的材料具有最理想的力学性能。
表 1 主要原料参数Table 1. Parameters of raw materials原料 规格 生产厂家 聚四氟乙烯 粒径26 μm 美国Dupont公司 铝粉 粒径10 μm,纯度99.9% 河南远洋铝业 高纯氩气 纯度≥99.99% 北京亚男伟业 无水乙醇 纯度≥99.7% 北京化工厂 实验中将防护屏制备成直径110 mm的圆板,为了准确控制防护屏面密度,材料参数需要满足公式:
ρAlhAl=ρ0h0 (1) 式中:
ρAl 、hAl 分别为铝合金的密度及厚度,ρ0 、h0 为对应PTFE/Al活性材料的密度、厚度。PTFE/Al活性材料模压成型后,还需要经历烧结工艺,由于烧结过程中基体PTFE经历晶体熔化、再结晶,并伴有孔洞闭合、内应力释放等机制,烧结前后材料体积会根据成型压力的大小有所变化,所以很难单独确定
ρ0 或者h0 。制备活性材料板所需质量的计算公式为:M=πR2ρ0h0 (2) 式中:R为模具半径。
将式(1)代入式(2),可得:
M=πR2ρAlhAl (3) 为了保证面密度相同,只需要确定所需活性材料的质量即可,采用这种方法制备的薄板面密度误差小于1%。为了尽可能地接近工程应用背景,同时考虑活性材料制备工艺的可行性,实验中选用面密度为0.84 g/cm2的PTFE/Al活性材料防护屏,对应粉末质量为81.2 g。
1.2 实验系统与方案
采用二级轻气炮进行实验,发射管口径为14.5 mm,最高发射速度可达7.1 km/s。超高速撞击实验靶板结构选用经典的Whipple防护结构(见图2),撞击角度为0°,防护间距为100 mm,球形弹丸材料为LY-12铝。
采用磁测速系统测量弹丸速度,并采用激光阴影摄像机记录弹丸的超高速撞击过程。为了保证弹丸完整性,采用气动分离系统实现弹丸与弹托的分离,气动弹托及弹托拦截靶拦截前后对比如图3所示。
1.3 实验结果
共进行9发超高速撞击实验,除实验1的防护屏为铝合金外,其余实验的防护屏均为PTFE/Al活性材料,面密度均为0.84 g/cm2。球形弹丸直径分别为5.0、6.0、6.4 mm,弹丸撞击速度范围为2.3~6.1 km/s。实验结果见表2,主要包括弹丸参数、Whipple防护结构参数以及后板损伤情况。
表 2 实验参数及损伤情况Table 2. Hypervelocity impact test configurations and damage results实验 弹丸参数 Whipple防护结构参数 后板损伤 直径/mm 质量/g 速度/(km·s−1) 前板材料 面密度/(g·cm−2) 后板材料 1 6.4 0.38 5.06 LY-12铝 0.84 LY-12铝 穿孔撕裂 2 6.4 0.38 5.03 PTFE/Al 0.84 LY-12铝 轻微层裂 3 5.0 0.18 3.79 PTFE/Al 0.84 LY-12铝 鼓包 4 5.0 0.18 3.88 PTFE/Al 0.84 LY-12铝 鼓包 5 5.0 0.18 4.00 PTFE/Al 0.84 LY-12铝 鼓包 6 6.0 0.31 3.71 PTFE/Al 0.84 LY-12铝 穿孔 7 6.4 0.38 6.08 PTFE/Al 0.84 LY-12铝 轻微鼓包 8 5.0 0.18 2.65 PTFE/Al 0.84 LY-12铝 鼓包 9 6.0 0.31 2.31 PTFE/Al 0.84 LY-12铝 鼓包 2. 实验结果分析
2.1 碎片云特性
文献[14]采用高速摄像机记录弹丸的超高速撞击过程,证明活性材料在撞击瞬间发生了可靠的冲击起爆反应。由于冲击起爆反应过程伴有剧烈的发光现象,强光会导致图像的过渡曝光,使碎片云形貌湮没,所以高速摄像机无法有效记录含能防护结构碎片云演变过程。激光阴影摄影系统可以很好地解决这一问题[15]。
图4~5给出了弹丸分别撞击铝合金防护屏(实验1)、含能活性防护屏(实验2)后碎片云的激光阴影照片,摄影系统的时间序列间隔为5 μs。由于阴影照片的视场直径只有80 mm,所以图片中未显示防护结构的后板。
由图4可知,碎片云的初始时刻由于没有充分膨胀,碎片密度大,激光光源无法穿透,所以照片只能显示碎片云的外部轮廓,不能分辨出块状碎片。随着碎片云飞行过程中沿轴向、径向的不断膨胀,碎片云密度减小,逐渐有光源可以透过,可以看到碎片云中包含有非常多的大尺寸碎片,这些碎片由铝合金防护屏和弹丸材料共同组成,可导致后板成坑甚至穿孔失效。
不同于图4中碎片云的椭球形膨胀过程,图5中弹丸撞击活性材料防护屏后形成的碎片云膨胀过程更接近于球形,且在整个拍摄过程中碎片云都不透光,没有观测到大尺寸碎片。这是因为活性材料超高速撞击后发生爆炸反应,生成高温、高压、高密度的气体产物与弹丸破碎产生的碎片共同组成碎片云,高密度的气体产物导致激光光源无法穿过,同时高压气体产物在真空中会迅速膨胀,使碎片云的轮廓逐渐接近于球形。这一现象也间接证明PTFE/Al活性材料在超高速撞击下发生了爆炸反应。此时,由于活性材料防护屏超高速撞击破碎后大部分发生爆轰或爆燃反应,少量没有反应的PTFE/Al碎片动能非常低,使得碎片云中具有侵彻能力的碎片主要来自于弹丸材料,从而大幅减少有害碎片数量,导致后板损伤情况减弱。
2.2 后板损伤特性
空间碎片防护结构后板损伤特性的分析对于获取碎片云的形成信息具有重要意义[16],通过观察分析后板损伤特点,可以在一定程度上反推碎片云的结构与形态。实验中回收的典型后板损伤情况如图6~9所示。
实验1~2的防护屏材料分别为铝合金、活性材料,防护屏面密度等其他实验参数保持一致,后板损伤情况如图6~7所示。实验1中后板正面中间有一个直径25 mm左右的严重毁伤区,背面明显发生层裂,甚至出现长10 mm左右的撕裂区,直径3 mm以上的弹坑数目多达46个,而实验2后板则未出现大面积毁伤区域,后板损伤程度大幅降低,直径3 mm以上的弹坑数目只有13个,这是因为活性材料超高速撞击后的反应产物基本不具备侵彻能力,直观反映就是后墙弹坑数目的大幅降低。
为了更准确地分析活性材料防护结构中后板损伤情况,得到撞击速度对弹丸破碎情况的影响规律,对3种尺寸的弹坑数量进行了统计,分别为dc>4 mm(大弹坑),3 mm≤dc≤4 mm(中弹坑),2 mm≤dc<3 mm(小弹坑),详细结果见表3。
表 3 后板损伤情况统计Table 3. Damage statistics of rear wall实验
编号撞击速度/
(km·s−1)弹坑数目 弹坑总数
(dc≥2 mm)(dc>4 mm) (3 mm≤dc≤4 mm) (2 mm≤dc<3 mm) 2 5.03 6 5 12 23 3 3.79 2 2 0 4 4 3.88 5 0 3 8 5 4.00 3 3 1 7 7 6.08 0 6 9 15 8 2.65 1 0 0 1 9 2.31 1 0 0 1 实验8~9中后板只有1个4 mm以上弹坑,表明弹丸与活性材料防护屏作用后并没有发生整体破碎,仍以比较完整的弹体撞击后板,说明撞击速度为2.65、2.31 km/s时,碰撞与爆炸两种作用下产生的冲击力不足以使弹丸发生整体破碎,对应的撞击速度处于弹道段。
实验3~5中弹丸速度均为约4 km/s,此时弹坑总数明显增多,弹坑分布近似于中心圆域,说明此时在碰撞与爆炸两种作用下弹丸已经发生破碎,撞击速度处于破碎段,对应的撞击速度应该高于弹丸的临界破碎速度。在空间碎片防护结构撞击极限方程中,临界破碎速度定义为弹丸刚好发生破碎时的速度,也是弹道区域破碎区的分界点[17]。
实验2中撞击速度增大为5.03 km/s,此时无论是小弹坑数目还是弹坑总数都出现大幅增加,说明相较于实验3~5中4 km/s左右的撞击速度,此时弹丸材料发生更充分的破碎。弹坑分布也开始发生变化,不再集中在后板中心区域,而是出现不太规则的环状分布,中心区域发生轻微的整体凹陷。
实验7中,撞击速度进一步增大到6.08 km/s,此时,直径大于4 mm的大弹坑数目为零,说明弹丸破碎更加充分;但相较于实验2,弹坑总数反而大幅降低,弹坑分布呈圆环状,直径约52 mm,中心圆形区域发生明显的整体凹陷,如图10所示。由碎片云形成理论可知,如果发生了熔化、气化,这时碎片云是由固、液、气三相组成的气泡。这时碎片云对后板的破坏主要表现为“气泡”对主结构层的整体破坏,起决定作用的是后板单位面积上的冲量。冲量足够大时,可能会使后板出现整体的撕裂破坏,冲量较小时,破坏模式表现为局部凹陷。从后板的损伤情况可以推测,6.08 km/s撞击条件下,弹丸中部与活性材料防护屏充分接触的部分在碰撞与爆炸的共同作用下发生液化、气化,而在弹丸边缘,稀疏波作用下活性材料不足以充分爆轰,释放能量明显小于中心区域,导致弹丸不能完全液化、气化,而是形成中间为液、气两相气泡,周围为固相颗粒的混合碎片云。分析表明,对于实验中的防护结构,6.08 km/s的撞击速度已经处于熔化/气化区。
图11中给出了弹坑总数随速度的变化趋势,从图7可以清晰地观察到弹坑数目的变化主要分为3个阶段,当弹丸速度位于弹道区时,弹丸保持完整,弹坑数目为1;当弹丸进入破碎区后,开始产生破碎,随着速度的增大,破碎程度迅速增加,表现为弹坑数目迅速增多;当速度进一步增大时,由于弹丸材料开始发生熔化、气化,对后板的损伤转变为整体冲量破坏,弹坑数目开始减少。
2.3 撞击极限分析
撞击极限是航天器空间碎片高速撞击风险评估、防护结构和防护材料优化设计的重要依据,是指结构在撞击后发生失效与不发生失效的临界状态,通常的定义形式包括弹丸的直径及撞击速度、板的厚度、防护间距,舱壁强度等。防护结构后板穿孔或层裂则定义为结构失效,失效区位于撞击极限曲线上方,未失效区位于曲线下方。
目前应用最广泛的是Christiansen方程[17],可用于描述铝合金Whipple双层结构撞击极限。根据弹丸不同速度下破坏机理的不同,将速度分为3个阶段:弹道区、破碎区、熔化/气化区,并直接给出弹丸撞击极限尺寸。一般情况下,对于铝合金Whipple双层结构分段速度为别为3、7 km/s。方程具体形式如下:
当v0 cos θ≤3 km/s时:
dp=[tw(σ40)0.5+tb0.6(cosθ)5/3ρ1/2pv2/30]18/19 (4) 当3 km/s<v0 cos θ<7 km/s时:
dp=[tw(σ40)0.5+tb1.248ρ1/2pcosθ]18/19(1.75−v0cosθ4)+[1.071t2/3wρ−1/3pρ−1/9bS1/3(σ70)1/3](v0cosθ4−0.75) (5) 当v0cosθ≥7 km/s时:
dp=3.918t2/3wρ−1/3pρ−1/9b(v0cosθ)−2/3s1/3(σ70)1/3 (6) 式中:
dp 为临界弹丸直径,cm;ρp 为弹丸密度,g/cm3;ρb 为防护屏密度,g/cm3;tw 为后板厚度,cm;S为防护间距,cm;v0 为弹丸撞击速度,km/s;θ 为撞击速度与靶板法线方向夹角,(°);σ 为后板屈服强度,MPa。实验2弹丸直径为6.4 mm,撞击速度为5.03 km/s,此时虽然后板没有发生穿孔,但后表面产生了小尺寸层裂区域,判断为发生轻微层裂,防护结构恰好处于临界防护状态,则6.4 mm为防护结构的临界弹丸直径。由方程(5)可知,此速度下铝合金防护结构对应的临界弹丸直径为5.04 mm,计算可知活性材料防护结构防护能力提高约27%,活性材料防护结构撞击极限曲线刚好通过该实验点。实验7保持弹丸直径不变,将撞击速度提高到6.08 km/s,此时防护结构后板只是轻微的鼓包,没有发生失效,说明活性材料防护结构的防护能力在该速度段内与速度成正比,实验点应处于活性材料防护结构撞击极限曲线的下方。
实验3~5表明,弹丸直径为5 mm时,在3.79~4.0 km/s速度范围内,防护结构均处于有效防护状态,为了能够在这个速度段内找到含能防护结构的失效点,实验6中将弹丸直径增大到6 mm,速度降为3.71 km/s,此时防护结构发生失效。对比铝合金防护结构临界弹丸直径可知,4.0 km/s左右时活性材料防护结构防护能力提升约25%~50%,活性材料防护结构撞击极限曲线应该处于实验点3~5和6对应的数据点之间。
为了进一步探索活性材料防护结构在弹道区(≤3 km/s)的防护能力,实验8的撞击速度降为2.65 km/s,对应的弹丸直径为5 mm,此时防护结构未发生失效。由式(4)可知,活性材料防护结构防护能力提高了31%。实验9撞击速度为2.31 km/s,弹丸直径增大为6 mm,但此时防护结构依然没有失效,防护能力提高大于45%。由此可知,在弹道段,撞击速度的降低有利于防护能力的提高,这是因为活性材料超高速撞击下完全反应时释放的能量是一定的,当弹丸速度降低时,释放能量与弹丸动能之比增大,防护能力提高。但这是以活性材料完全反应为前提的,由于PTFE/Al活性材料的冲击起爆特性与撞击速度密切相关,当弹丸速度过低以致其穿靶过程中活性材料没有来得及反应时,活性材料防护屏也就失去了其防护能力的立足之本,加上其无论是波阻抗还是屈服强度都远低于铝合金,所以在弹道段总有一个速度点使PTFE/Al活性材料与铝合金防护能力相当,当低于这个速度时铝合金防护能力更强,高于这个速度时PTFE/Al活性材料更强。利用最小二乘法,结合实验8~9对弹道段撞击极限曲线进行拟合,过防护状态临界点(实验4)并结合实验5~9对破碎段撞击极限进行拟合,获得活性材料面密度为0.84 g/cm2时对应的弹道段与破碎段的撞击极限曲线,如图12所示,两者的交点为弹道段与破碎段的分界点,黑色对比曲线为同面密度铝合金防护结构对应的撞击极限曲线。通过对比可知,活性材料防护结构防护能力获得大幅提升。
2.4 防护机理分析
弹丸与防护屏超高速撞击作用下的破碎、熔化甚至气化现象,是Whipple结构防护机理的核心,弹靶冲击压力则是弹丸与防护屏破碎、相变的能量源,活性材料的应用直接改变了这种能量源。首先,弹丸撞击活性材料产生的冲击压力不仅包括冲击压缩引起的压力,还包括因化学能释放而贡献的压力,两种压力的共同作用增加了弹丸的破碎与相变程度。同时,活性材料的冲击起爆过程伴随有高温、高压产物的剧烈膨胀,会对弹丸产生反向冲量而降低弹丸的轴向动能。其次,后板弹坑分析表明,相较于铝合金结构,活性材料结构后墙弹坑数量大幅减少,这是因为活性材料冲击起爆过程伴随有化学能的迅速释放,导致自身产物多为气态,对后板没有侵彻能力,后板的侵彻破坏只来自于碎片云中弹丸破碎产生的碎片。
3. 结 论
利用二级轻气炮完成PTFE/Al活性材料的超高速撞击实验,速度范围为2.3~6.1 km/s,覆盖弹道区、破碎区。碎片云特性及后板损伤特性分析表明,由于活性材料超高速撞击发生冲击起爆反应,爆炸产物尺寸非常细小且多为气态,此时碎片云中具有侵彻能力的碎片仅由弹丸破碎产生,大大减弱了对后板的破坏。撞击极限分析显示,相较于同面密度铝合金,活性材料具有更加优异的防护性能,提升比例与撞击速度密切相关。活性材料良好的防护特性使其可能成为未来航天器防护材料的新选择。
-
表 1 常见的指数函数类爆炸荷载
Table 1. Typical explosive loads in exponential function forms
药卷直径dc/mm ρe/(kg·m−3) D/(m·s−1) 爆热Q/(MJ·kg−1) 70 1200 4000 3.991 ρ/(kg∙m−3) E/GPa μ σc/MPa σt/MPa ϕ/(°) ψ 2670 59.5 0.23 129.1 10.3 45 2 ρs/(kg·m−3) 弹性模量Es/GPa 泊松比μs φs/(°) fd 1800 0.2 0.30 28 0.055 密度/(kg·m−3) 爆速/(m·s−1) 爆热/(kJ·g−1) 装药直径/mm 821 3800 4.0 95 表 6 单孔爆破试验岩体参数
Table 6. Rock mass parameters for single hole blasting test
密度/(kg·m−3) 弹性模量/GPa 泊松比 单轴抗压强度/MPa 单轴抗拉强度/MPa 2400 60 0.26 116 10.9 表 7 单孔爆破试验钻孔装药参数
Table 7. Blasting parameters for single hole blasting test
炮孔直径/mm 炮孔长度/m 炮孔倾角/(°) 药卷直径/mm 装药长度/m 单孔药量/kg 填塞长度/m 115 9.5 90 90 4.5 30 5.0 -
[1] 汪旭光. 中国工程爆破与爆破器材的现状及展望 [J]. 工程爆破, 2007, 13(4): 1–8. DOI: 10.3969/j.issn.1006-7051.2007.04.001.WANG X G. Current status and future prospect of engineering blasting and explosive materials in China [J]. Engineering Blasting, 2007, 13(4): 1–8. DOI: 10.3969/j.issn.1006-7051.2007.04.001. [2] WANG M L, LI X, LI Q H, et al. Study on blasting technology for open-pit layering of complex mine adjacent to high and steep slope [J]. Frontiers in Earth Science, 2021, 9: 773872. DOI: 10.3389/feart.2021.773872. [3] 王志亮, 毕程程, 李鸿儒. 混凝土爆破损伤的SPH-FEM耦合法数值模拟 [J]. 爆炸与冲击, 2018, 38(6): 1419–1428. DOI: 10.11883/bzycj-2017-0209.WANG Z L, BI C C, LI H R. Numerical simulation of blasting damage in concrete using a coupled SPH-FEM algorithm [J]. Explosion and Shock Waves, 2018, 38(6): 1419–1428. DOI: 10.11883/bzycj-2017-0209. [4] 杨建华, 代金豪, 姚池, 等. 爆破开挖扰动下锚固节理岩质边坡位移突变特征与能量机理 [J]. 爆炸与冲击, 2022, 42(3): 035201. DOI: 10.11883/bzycj-2021-0126.YANG J H, DAI J H, YAO C, et al. Displacement mutation characteristics and energy mechanisms of anchored jointed rock slopes under blasting excavation disturbance [J]. Explosion and Shock Waves, 2022, 42(3): 035201. DOI: 10.11883/bzycj-2021-0126. [5] HUSTRULID W. Blasting principles for open pit mining: theoretical foundations [M]. Rotterdam: Balkema, 1999: 992. [6] 孙鹏昌. 基于模态分析的岩石高边坡爆破振动影响评价 [D]. 武汉: 武汉大学, 2021: 47–60.SUN P C. Impact evaluation of blasting vibration on high rock slope based on modal analysis [D]. Wuhan: Wuhan University, 2021: 47–60. [7] CASTEDO R, NATALE M, LÓPEZ L M, et al. Estimation of Jones-Wilkins-Lee parameters of emulsion explosives using cylinder tests and their numerical validation [J]. International Journal of Rock Mechanics and Mining Sciences, 2018, 112: 290–301. DOI: 10.1016/j.ijrmms.2018.10.027. [8] KUZMENKO A A, VOROBEV V D, DENISYUK I I, et al. Seismic effects of blasting in rock [M]. Rotterdam: A. A. Balkema, 1993: 35–37. [9] 王文龙. 钻眼爆破 [M]. 北京: 煤炭工业出版社, 1984: 178–179.WANG W L. Drilling and blasting [M]. Beijing: China Coal Industry Publishing House, 1984: 178–179. [10] SHARPE J A. The production of elastic waves by explosion pressures: I. theory and empirical field observations [J]. Geophysics, 1942, 7(2): 144–154. DOI: 10.1190/1.1445002. [11] DUVALL W I. Strain-wave shapes in rock near explosions [J]. Geophysics, 1953, 18(2): 310–323. DOI: 10.1190/1.1437875. [12] CHO S H, MIYAKE H, KIMURA T, et al. Effect of the waveform of applied pressure on rock fracture process in one free-face [J]. Science and Technology of Energetic Materials, 2003, 64(3): 116–125. [13] CHO S H, KANEKO K. Influence of the applied pressure waveform on the dynamic fracture processes in rock [J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(5): 771–784. DOI: 10.1016/j.ijrmms.2004.02.006. [14] STARFIELD A M, PUGLIESE J M. Compression waves generated in rock by cylindrical explosive charges: a comparison between a computer model and field measurements [J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1968, 5(1): 65–77. DOI: 10.1016/0148-9062(68)90023-5. [15] JONG Y, LEE C, JEON S, et al. Numerical modeling of the circular-cut using particle flaw code [C]//Proceedings of the 31st Annual Conference on Explosives and Blasting Technique. Orlando, 2005: 289–298. [16] BLAIR D, MINCHINTON A. On the damage zone surrounding a single blasthole [J]. Fragblast, 1997, 1(1): 59–72. DOI: 10.1080/13855149709408390. [17] BLAIR D P. A comparison of Heelan and exact solutions for seismic radiation from a short cylindrical charge [J]. Geophysics, 2007, 72(2): E33–E41. DOI: 10.1190/1.2424543. [18] CHO S H, KANEKO K. Rock fragmentation control in blasting [J]. Materials Transactions, 2004, 45(5): 1722–1730. DOI: 10.2320/matertrans.45.1722. [19] TRIVIÑO L F, MOHANTY B, MUNJIZA A. Seismic radiation patterns from cylindrical explosive charges by analytical and combined finite-discrete element methods [C]//Proceedings of the 9th International Symposium on Rock Fragmentation by Blasting. Boca Raton, FL, USA: CRC Press, 2009: 415–426. [20] 陶振宇, 董振华, 卢文波. 敞口炮孔压力变化历程的理论分析与计算 [J]. 武汉水利电力大学学报, 1994, 27(4): 394–399.TAO Z Y, DONG Z H, LU W B. A theoretical analysis and calculation of the pressure course of unstemmed borehole [J]. Journal of Wuhan University of Hydraulic and Electric Engineering, 1994, 27(4): 394–399. [21] 卢文波, 陶振宇. 爆生气体驱动的裂纹扩展速度研究 [J]. 爆炸与冲击, 1994, 14(3): 264–268.LU W B, TAO Z Y. A study of fracture propagation velocity driven by gases of explosion products [J]. Explosion and Shock Waves, 1994, 14(3): 264–268. [22] 卢文波. 岩石爆破中应力波的传播及其效应研究 [D]. 武汉: 武汉水利电力大学, 1994: 15–24.LU W B. Propagation and effect of stress wave in rock blasting [D] Wuhan: Wuhan University of Hydraulic and Electric Engineering, 1994: 15–24. [23] LU W B, YANG J H, CHEN M, et al. An equivalent method for blasting vibration simulation [J]. Simulation Modelling Practice and Theory, 2011, 19(9): 2050–2062. DOI: 10.1016/j.simpat.2011.05.012. [24] 杨建华, 卢文波, 陈明, 等. 岩石爆破开挖诱发振动的等效模拟方法 [J]. 爆炸与冲击, 2012, 32(2): 157–163. DOI: 10.11883/1001-1455(2012)02-0157-07.YANG J H, LU W B, CHEN M, et al. An equivalent simulation method for blasting vibration of surrounding rock [J]. Explosion and Shock Waves, 2012, 32(2): 157–163. DOI: 10.11883/1001-1455(2012)02-0157-07. [25] 卢文波, 杨建华, 陈明, 等. 深埋隧洞岩体开挖瞬态卸荷机制及等效数值模拟 [J]. 岩石力学与工程学报, 2011, 30(6): 1089–1096.LU W B, YANG J H, CHEN M, et al. Mechanism and equivalent numerical simulation of transient release of excavation load for deep tunnel [J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(6): 1089–1096. [26] 李庆文, 乔兰, 陈璐. 基于精确爆破载荷分析的安全距离判据 [J]. 工程力学, 2015, 32(10): 123–129. DOI: 10.6052/j.issn.1000-4750.2014.03.0248.LI Q W, QIAO L, CHEN L. Based on the accurate blasting loading to estimate the safety criterion [J]. Engineering Mechanics, 2015, 32(10): 123–129. DOI: 10.6052/j.issn.1000-4750.2014.03.0248. [27] 钱七虎. 岩石爆炸动力学的若干进展 [J]. 岩石力学与工程学报, 2009, 28(10): 1945–1968. DOI: 10.3321/j.issn:1000-6915.2009.10.001.QIAN Q H. Some advances in rock blasting dynamics [J]. Chinese Journal of Rock Mechanics and Engineering, 2009, 28(10): 1945–1968. DOI: 10.3321/j.issn:1000-6915.2009.10.001. [28] 陈士海, 王明洋, 赵跃堂, 等. 岩石爆破破坏界面上的应力时程研究 [J]. 岩石力学与工程学报, 2003, 22(11): 1784–1788. DOI: 10.3321/j.issn:1000-6915.2003.11.006.CHEN S H, WANG M Y, ZHAO Y T, et al. Time-stress history on interface between cracked and uncracked zones under rock blasting [J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(11): 1784–1788. DOI: 10.3321/j.issn:1000-6915.2003.11.006. [29] CHAPMAN D L, OXON A B. VI. On the rate of explosion in gases [J]. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 1899, 47(284): 90–104. DOI: 10.1080/14786449908621243. [30] JOUGUET E. On the propagation of chemical reactions in gases [J]. Journal de Mathématiques Pures et Appliquées, 1905, 1: 347–425. [31] HENRYCH J. The dynamics of explosion and its use: developments in civil engineering [M]. Amsterdam: Elsevier, 1979: 59–64. [32] Л. Д. 朗道, E. M. 栗弗席兹. 连续介质力学(第二册) [M]. 彭旭麟, 译. 北京: 人民教育出版社, 1960: 500–509. [33] WANG L L. Foundations of stress waves [M]. Oxford: Elsevier, 2007: 337–347. [34] 李世愚, 和泰名, 尹祥础. 岩石断裂力学 [M]. 北京: 科学出版社, 2015: 23–25.LI S Y, HE T M, YIN X C. Rock fracture mechanics [M]. Beijing: Science Press, 2015: 23–25. [35] 冷振东, 卢文波, 陈明, 等. 岩石钻孔爆破粉碎区计算模型的改进 [J]. 爆炸与冲击, 2015, 35(1): 101–107. DOI: 10.11883/1001-1455(2015)01-0101-07.LENG Z D, LU W B, CHEN M, et al. Improved calculation model for the size of crushed zone around blasthole [J]. Explosion and Shock Waves, 2015, 35(1): 101–107. DOI: 10.11883/1001-1455(2015)01-0101-07. [36] LU W B, LENG Z D, CHEN M, et al. A modified model to calculate the size of the crushed zone around a blast-hole [J]. Journal of the Southern African Institute of Mining and Metallurgy, 2016, 116(5): 413–422. DOI: 10.17159/2411-9717/2016/v116n5a7. [37] RAKISHEV B, RAKISHEVA Z B. Basic characteristics of the stages of rock massif destruction by explosive crushing [C]//Proceedings of the 3rd Asia-Pacific Symposium on Blasting Techniques. Kunming, 2011: 65–69. [38] MENCACCI S, CHAVEZ R. The measurement and analysis of detonation pressure during blasting [C]//Proceedings of 2005 European Federation of Explosives Engineers. Brighton, 2005: 231–237. [39] SUN P C, LU W B, ZHOU J R, et al. Comparison of dominant frequency attenuation of blasting vibration for different charge structures [J]. Journal of Rock Mechanics and Geotechnical Engineering, 2022, 14(2): 448–459. DOI: 10.1016/j.jrmge.2021.07.002. [40] SINGH P K, ROY M P. Damage to surface structures due to blast vibration [J]. International Journal of Rock Mechanics and Mining Sciences, 2010, 47(6): 949–961. DOI: 10.1016/j.ijrmms.2010.06.010. [41] LANGEFORS U, KIHLSTRÖM B. The modern technique of rock blasting [M]. 3rd ed. New York: Wiley, 1979: 28–56. 期刊类型引用(1)
1. 安德隆,肖建光,马俊杨,徐瑞泽. 聚四氟乙烯基活性材料化学反应分子动力学模拟. 中北大学学报(自然科学版). 2024(02): 222-228 . 百度学术
其他类型引用(3)
-