Criterion of Plate Structure Damage Caused by Underwater Explosion Shock Wave Based on Effective Impulse
-
摘要: 为了准确评估水下爆炸冲击波对平板结构造成的毁伤效果,提出以有效冲量作为毁伤准则的威力参量类别,并给出了考虑球面波斜入射效应的平板结构有效冲量计算方法。新准则通过动量守恒方程计算得到的平板结构实际获得冲量对比毁伤效果,表现为冲击波峰压、时间常数以及板结构特征参数的联合形式。借助数值模拟与文献数据,对比分析了准则的准确度和适用性。结果表明:相对于冲击波峰压、比冲量、能流密度等单一威力参量,新毁伤准则在评估平板结构的毁伤程度时误差更小;在对比不同炸药毁伤威力以及预估未知工况毁伤效果两种应用场景中,新准则的相对误差均在10%以内。新提出的毁伤准则用于对比和评估水下爆炸冲击波对板结构的毁伤效果时具有良好的通用性。Abstract: In order to accurately evaluate damage effect of underwater explosion shock wave on flat plate structure, the power parameter category with the effective impulse as the damage criterion is proposed, and a correction method considering the oblique incidence effect of the spherical wave is given. The new criterion compares the damage effect by the actual impulse of the plate structure, which is calculated by the momentum conservation equation. The new criterion is expressed in the joint form of peak pressure of shock wave, time constant and characteristic parameters of plate structure. With the help of simulation and literature data, the accuracy and applicability of the criteria are compared and analyzed. The results show that the new damage criterion has less error in evaluating the damage degree of flat plate structure, compared with single power parameters such as peak pressure of shock wave, specific impulse and energy flux density. In comparing the damage power of different explosives and predicting the damage effect under unknown conditions, the relative error of the new criterion is within 10%. The proposed damage criterion has good versatility when used to compare and evaluate the damage effect of underwater explosion shock wave on plate structure.
-
Key words:
- underwater explosion /
- shock wave /
- damage criterion /
- flat plate structure /
- effective impulse
-
乙炔是一种无色无味的可燃气体,广泛应用于金属焊接、切割等,同时也是化工生产中合成橡胶、芳烃和高聚物等的重要原材料[1-3]。乙炔在空气中的爆炸极限为2.5%~80%[4],具有很高的爆炸风险和破坏性。在乙炔的生产、运输和储存等环节,由于气体泄漏、通风和回火等原因,极易引发爆炸[5]。近年来,由乙炔引发的爆炸事故频发,例如:2018年,张家口盛华化工有限公司附近发生燃爆事故[6],装载乙炔的货车被爆燃的氯乙烯引爆,造成24人死亡、22人受伤;2018年,印度科钦造船厂因工人焊接导致乙炔泄漏引发爆炸事故,造成5人死亡、13人受伤。乙炔爆炸事故发展迅速、爆炸危险性大,给人们的生命和财产安全带来了巨大的威胁。因此,研究乙炔气体的抑爆技术对于保障安全生产具有重要的意义。
国内外学者在可燃气体抑爆方面做了大量的研究工作,常用的抑爆剂包括气体抑爆剂、液体抑爆剂和固体抑爆剂。Wang等[7]运用氮气和二氧化碳稀释烃类燃料,发现惰性气体的加入能够抑制燃料的爆炸强度,且二氧化碳的抑制效果优于氮气。Luo等[8]研究了气-固抑爆剂(二氧化碳-ABC粉末)对瓦斯气体爆炸的影响,实验结果表明,ABC粉末与二氧化碳存在协同增效,可以降低瓦斯的爆炸威力。Wang等[9]运用实验与数值模拟相结合的方法,研究了超细水雾对气体爆炸的抑制机理,结果表明,超细水雾对火焰温度、火焰传播速度和最大爆炸压力均有显著的衰减作用。Wen等[10]研究了惰性纳米粉末颗粒对甲烷-空气爆炸的抑制作用,发现氢氧化铝(Al(OH)3)、碳酸氢钠(NaHCO3)和二氧化硅(SiO2)都可以用于抑制火焰的传播和超压,并且金属氢氧化物比二氧化硅的抑制效果更好。丁超等[11]基于100 L的爆炸容器探究了超细ABC粉末对瓦斯爆炸的影响,结果表明,抑制爆炸所需的超细ABC粉体存在确定的临界量,且抑爆效果与抑爆装置的触发时间有关。国内外有关乙炔爆炸特性的研究较少,而涉及抑爆效果的研究则更少。在乙炔的爆炸特性方面:宋诗祥[12]利用20 L球形爆炸装置研究了水分子在乙炔-碳化钙气-固两相爆炸中的作用,揭示了其爆炸反应机制;郭璐等[3]探究了杂质对乙炔气体分解爆炸的影响,发现氮气对乙炔爆炸分解有抑制作用,而氧气、硫化氢、铁锈对其具有促进作用。Mizutani等[13]利用1 L密闭容器研究了乙炔的分解爆炸特性,结果表明,最大分解爆燃压力和压力上升速率随着容器温度的降低而升高。在乙炔的抑爆特性方面:王犇等[14]采用20 L球形爆炸测试系统研究了氮气对乙炔分解爆炸特性的影响,结果显示,乙炔的爆炸抑制作用随着氮气体积分数的增加而逐渐增大;Kopylov等[15]研究了丙烷、甲烷和18%С3Н6-40%С3Н8-42%С4Н10混合物对乙炔-空气混合气体燃烧的影响,发现当乙炔的体积分数在2%~8%范围内时,丙烷和18%С3Н6-40%С3Н8-42%С4Н10混合物对乙炔的燃烧具有较强的抑制作用,并且随着混合物中乙炔含量的增加,抑制作用减弱。目前,有关乙炔抑爆研究主要集中在气体和液体抑爆剂方面,而固体抑爆剂对乙炔的抑爆特性研究很少。
SiO2、Al(OH)3和NaHCO3[16-19]是典型的固体抑爆剂。本研究采用20 L球形爆炸测试系统,研究3种固体抑爆剂对乙炔-空气(当量比
φ =1)预混气体的抑爆效果,通过爆炸特性参数(最大爆炸压力pex、最大爆炸压力上升速率(dp/dt)ex)分析固体抑爆剂对乙炔的抑爆效果,结合乙炔气体爆炸机理,阐明固体抑爆剂的抑爆机理,为乙炔工业生产、储运和使用过程中乙炔意外爆炸防治提供理论指导。1. 实验装置与材料
1.1 实验装置
本实验采用的20 L球形爆炸测试系统如图1所示。该系统主要由爆炸容器、点火系统、喷粉系统、配气系统、数据采集系统和控制系统组成。爆炸容器为球形,体积为20 L;点火系统由导线和点火电极构成;喷粉系统由气固两相阀、储粉罐(600 mL)、压力表和压缩气体构成;数据采集系统由压力传感器和计算机构成;控制系统由导线和控制器构成。首先,在储粉罐中装入固体抑爆剂粉体,用真空泵将爆炸容器抽至真空,采用分压法向其中通入乙炔气体;然后,将高压空气通入储粉罐,当储粉罐前的电接点式压力表达到2.0 MPa后,开启气固两相阀,将固体抑爆剂粉体喷入爆炸容器中,此时爆炸容器内部达到常压,固体抑爆剂粉体弥散形成粉尘云;最后,采用控制系统点火,由压力传感器(美国Dytran)采集数据,并传输至计算机,经计算机软件初步处理得到压力时程曲线,存储图像和数据。实验中,乙炔与空气的当量比
φ =1,每种工况下重复3次以上实验。1.2 实验材料
本实验采用的乙炔气体由安徽合肥恒隆电气有限公司提供,纯度为99.9%;空气由安徽合肥恒隆电气有限公司提供,氧气和氮气的体积分数分别为20.98%和79.02%。考虑到粒子的表面积会影响吸热效果[20],为此本研究统一选用微米级抑爆剂粉体。微米级SiO2、Al(OH)3和NaHCO3均由国药集团化学试剂有限公司提供,其粒径分布如图2所示。可见,SiO2、Al(OH)3和NaHCO3粉体的粒径分布接近,平均粒径D50分别为13.0、10.4和11.5 μm。
2. 实验和模拟结果
乙炔-空气爆炸的主要反应式为[21]
C2H2+2.5(O2+3.76N2)→2CO2+H2O+9.4N2 (1) 采用乙炔与空气的当量比
φ =1的乙炔-空气预混气体进行爆炸实验φ = n/nst (2) 式中:n为实验中乙炔与空气的物质的量之比;nst为理论上乙炔完全反应时乙炔与空气的物质的量之比[5]。图3为20 L球形爆炸容器内
φ=1 时乙炔-空气预混气体的爆炸压力时程曲线(未加入固体抑爆剂的实验结果,即对照组),其最大压力pex为0.68 MPa,最大压力上升速率(dp/dt)ex为203.09 MPa/s。2.1 SiO2对乙炔-空气预混气体爆炸特性的影响
SiO2的化学性质稳定,熔点高达(1723±5) ℃,沸点高达2230 ℃,是一种常用的耐火材料和固体抑爆剂。微米级SiO2因表面效应特性,在催化、扩散、光吸收、热阻和自由基捕捉能力等方面具有广泛的应用[22]。图4给出了在不同粉体浓度的SiO2作用下乙炔-空气爆炸压力时程曲线,除100和300 g/m3粉体浓度外,高浓度的SiO2对乙炔-空气爆炸压力起到了抑制作用。图5为不同SiO2粉体浓度对乙炔-空气预混气体的爆炸特征参数的影响。当SiO2粉体浓度小于300 g/m3时,加入SiO2粉体能够快速吸收燃烧产生的水分,使燃烧反应加速,此时SiO2粉体表现出对乙炔爆炸参数的促进作用[23];当SiO2粉体浓度超过300 g/m3时,SiO2粒子的吸热作用和吸收活性自由基作用增强,远高于SiO2粒子吸收水分的影响,从而导致pex和(dp/dt)ex均随着SiO2粉体浓度的增大而不断下降。
2.2 Al(OH)3对乙炔-空气预混气体爆炸特性的影响
Al(OH)3是一种两性氧化物,熔点为300 ℃,受热易分解,在378.7 ℃出现热分解吸热峰[24],具有无毒、阻燃和不产生二次污染的优点,被广泛用于阻燃剂[25-26]。图6显示了在不同Al(OH)3粉体浓度下乙炔-空气的爆炸压力时程曲线。可见,Al(OH)3粉体对乙炔-空气的爆炸压力起到抑制作用。图7显示了不同Al(OH)3粉体浓度下乙炔-空气预混气体的爆炸特征参数。随着Al(OH)3粉体浓度的增加,pex和(dp/dt)ex均呈下降趋势。当Al(OH)3的粉体浓度小于300 g/m3时,pex和(dp/dt)ex下降缓慢,说明Al(OH)3粉体量较少时,对乙炔爆炸特性的影响较小;当Al(OH)3粉体浓度大于300 g/m3时,随着Al(OH)3粉体浓度的增加,pex和(dp/dt)ex的下降速率加快。
2.3 NaHCO3对乙炔-空气预混气体爆炸特性的影响
NaHCO3的化学性质不稳定,分解温度在110~140 ℃[27]之间,易分解生成惰性物质,因此具有良好的灭火作用,是BC类干粉灭火剂的主要填充材料[28]。图8显示了不同NaHCO3粉体浓度下乙炔-空气爆炸压力时程曲线。可以看出,NaHCO3对乙炔-空气的爆炸压力起到抑制作用。图9显示了不同NaHCO3粉体浓度下乙炔-空气预混气体的爆炸特征参数变化情况。随着NaHCO3粉体浓度的增加,pex和(dp/dt)ex均呈不断下降的趋势。这是因为在乙炔-空气爆炸火焰传播过程中,NaHCO3粉体能够使前锋面的火焰出现局部熄灭,由于微米级粒子的比表面积高,因此分解反应速率快,与无NaHCO3粉体的对照组相比,pex和(dp/dt)ex下降明显,抑制火焰发展的效果好。随着粉体浓度的增加,爆炸容器内参与抑制作用的粉体颗粒数目增多,对乙炔-空气爆炸的抑制作用增强,pex和(dp/dt)ex逐渐下降。
2.4 爆炸压力敏感性分析
为了进一步研究固体抑爆剂的影响机制,选取与实验相同当量的乙炔-空气预混气体,对影响爆炸压力的主要反应进行敏感性分析,进而探究固体抑爆剂对乙炔爆炸的抑爆机理。设初始温度为1400 K,初始压力为0.1 MPa,时间为5 ms,影响爆炸压力的主要反应见表1。图10给出了当量比为1的乙炔-空气预混气体爆炸压力的敏感性分析结果。
表 1 影响爆炸压力的主要反应Table 1. Main reactions affecting explosion pressureNo. Reaction No. Reaction R1 H+O2↔O+OH R158 C2H2+O↔HCCO+H R39 HCO+M↔CO+H+M R161 C2H2+OH↔C2H+H2O R41 HCO+O2↔CO+HO2 R166 C2H2+C2H↔C4H2+H R141 C2H+O2↔HCO+CO R194 C2H3+O2↔CH2CHO+O R155 C2H3(+M) ↔C2H2+H(+M) R195 C2H3+O2↔HCO+CH2O 从图10可以看出:R1、R39、R158、R161、R166和R194反应促进爆炸压力,R41、R141、R155和R195反应抑制爆炸压力,其中R1反应占主导地位。加入固体抑爆剂时,固体抑爆剂会与其中的关键自由基结合,从而影响爆炸特征参数。
3. 分析与讨论
图11为3种固体抑爆剂对乙炔-空气预混气体的抑爆特性对比。3种固体抑爆剂中,从pex和(dp/dt)ex的减小程度上看,SiO2粉体对乙炔-空气预混气体爆炸的抑制效果弱于Al(OH)3粉体,NaHCO3粉体对乙炔-空气爆炸的抑制效果最好。乙炔-空气预混气体爆炸是一种剧烈的自由基氧化链式反应,在反应过程中会生成O、H、OH和HCO等高活性自由基[22]。虽然微米级SiO2粒子(固相)不与乙炔-空气爆炸火焰反应,但是由于其具有很高的比表面积,粒子很容易接触和吸收燃烧反应热,从而直接降低燃烧反应强度,影响反应的传热,进而降低爆炸火焰速度和压力波[17]。当SiO2粉体浓度大于300 g/m3时,随着粉体浓度的增加,能够吸收更多的自由基,对乙炔爆炸的抑制效果更好。
Al(OH)3在乙炔-空气爆炸过程中既充当吸热介质,又发生分解反应[18],即
2Al(OH)3→Al2O3+3H2O (3) Al(OH)3分解生成Al2O3(固相)和H2O(气-液二相),产生吸热冷却效应。H2O汽化吸热,导致爆炸容器内的温度下降,生成的水蒸气降低了乙炔分子与氧分子的碰撞概率。同时,Al(OH)3中的自由基OH能够消耗乙炔-空气链式反应中的关键自由基O,抑制反应R158(C2H2+O↔HCCO+H)的进行,使得爆炸压力降低。Al(OH)3的吸热冷却效应和分解生成的H2O汽化吸热使得Al(OH)3的抑爆效果略优于SiO2,表现出更低的pex和(dp/dt)ex。
NaHCO3的性质不稳定,NaHCO3晶体在熔化过程中可吸收热量,并且在乙炔-空气预混气体爆炸过程中发生分解反应[18],即
2NaHCO3→Na2CO3+CO2+H2O (4) 分解反应的产物为惰性CO2气体(气相)和H2O(气-液二相),是气体爆炸中的惰性物质。其中:分解生成的Na2CO3(固相)使燃烧反应的传热受到阻碍;H2O汽化能够吸收大量的热量,影响乙炔粒子之间的传热。同时,生成的CO2气体作为一种惰性物质能够稀释乙炔和氧气的体积分数,阻碍火焰的发展,产生气体惰化。此外,NaHCO3中的HCO消耗乙炔-空气爆炸链式反应中的关键自由基O和OH[18],由敏感性分析可知,R158和R161反应对爆炸压力起促进作用,加入NaHCO3使关键自由基被结合,参与反应的自由基减少,爆炸压力的促进作用减弱。同时,NaHCO3中的HCO参与R41反应,R41反应对爆炸压力起抑制作用,爆炸压力的抑制作用增强,发生化学抑制。晶体熔化吸热、反应吸热、惰性物质和化学抑制的耦合作用使得爆炸容器内的温度和乙炔气体的体积分数降低,从而抑制乙炔-空气的爆炸反应,pex和(dp/dt)ex的降幅最大。因此,NaHCO3粉体的抑制爆炸效果优于SiO2和Al(OH)3粉体。
4. 结 论
(1) 当SiO2粉体浓度低于300 g/m3时,SiO2粉体对乙炔-空气预混气体的爆炸参数具有促进作用,当SiO2粉体浓度超过300 g/m3时,乙炔-空气预混气体的爆炸参数随着SiO2粉体浓度的增加而显著降低。
(2) SiO2粉体对乙炔-空气预混气体爆炸的抑制效果弱于Al(OH)3粉体,NaHCO3粉体对乙炔-空气爆炸的抑制效果最好。
(3) NaHCO3在乙炔-空气预混气体的爆炸过程中兼具气、固、液三相抑爆特点,通过晶体熔化吸热、反应吸热、惰性物质和化学抑制的多种耦合作用,使其表现出最佳的抑爆效果。
-
表 1 数值模拟计算材料模型
Table 1. Material model in numerical simulation
Material EOS Strength model Failure model TNT-2 JWL Air Ideal gas Water Shock Hydro Q235 steel Linear Johnson-Cook Johnson-Cook 表 2 冲击波冲量与时间常数数值模拟结果验证
Table 2. Verification of simulation results of shock wave impulse and time constant
Standoff
distance/mW/g Time decay constants Impulse Numerical
simulation/msEmpirical
formula/msRelative
error/%Numerical
simulation/(kPa·s)Empirical
formula/(kPa·s)Relative
error/%0.4 50 0.0300 0.0315 −4.87 1.976 1.969 0.35 0.5 50 0.0319 0.0332 −4.06 1.607 1.614 −0.45 0.6 50 0.0349 0.0346 0.91 1.369 1.372 −0.26 表 3 冲击波作用下圆板的变形挠度
Table 3. Deflection of circular plate under shock wave
D ω/mm W=50 g W=100 g W=200 g W=400 g W=800 g 10 59.52 71.98 85.32 99.22 113.90 16 40.85 47.53 54.85 62.84 71.83 20 32.46 37.73 43.40 49.92 57.24 26 24.24 28.13 32.77 37.82 43.60 30 20.29 23.75 27.86 32.39 37.54 表 4 炸药的相关系数 α,
k Table 4. Coefficients α,
k for explosivesExplosive type αp kp αθ kθ TNT 1.13 52.4 −0.23 0.084 H-6 1.19 59.2 −0.28 0.088 Pentolite 1.14 56.5 −0.23 0.084 表 5 不同工况下结构变形挠度计算误差对比
Table 5. Comparison of calculation errors of deflection under different working conditions
CI,plate Explosive W/g Standoff distance/m ω/mm Relative error of the same explosive/% Relative error with TNT/% 1.0 TNT 200 0.749 3.53 3.68 TNT 400 1.039 3.66 3.68 H-6 200 0.849 3.61 3.05 −1.95 H-6 400 1.166 3.72 3.05 −1.95 Pentolite 200 0.806 3.34 3.29 5.56 Pentolite 400 1.114 3.45 3.29 5.56 1.6 TNT 200 0.452 6.00 0.67 TNT 400 0.647 6.04 0.67 H-6 200 0.530 6.16 −1.79 −1.14 H-6 400 0.745 6.05 −1.79 −1.14 Pentolite 200 0.493 5.91 −2.71 3.16 Pentolite 400 0.699 5.75 −2.71 3.16 表 6 3种工况下的结构变形挠度
Table 6. Deflections under three working conditions
Standoff distance/m ω/mm W/kg 2.131 49.68 1 1.705 90.86 1 1.421 120.30 1 表 7 平板结构变形挠度预测值与实验结果的对比
Table 7. Comparison between predicted and experimental values of deflection
W/kg Standoff distance/m ω/mm Predicted deflection/mm Prediction error/% 1 1.136 149.94 160.26 6.88 1 0.710 243.81 222.24 −8.85 -
[1] KEIL A H. The response of ships to underwater explosions [R]. New York: Society of Naval Architects and Marine Engineers, 1961: 43. [2] REID W D. The response of surface ships to underwater explosions [R]. Canberra, Australia: Defence Science and Technology Organization, 1996. [3] 任鹏, 田阿利, 张伟, 等. 水下冲击波载荷作用下气背固支圆板动态毁伤实验 [J]. 爆炸与冲击, 2016, 36(5): 617–624. doi: 10.11883/1001-1455(2016)05-0617-08REN P, TIAN A L, ZHANG W, et al. Failure mode of clamped air-back circular panel subjected to underwater shock loading [J]. Explosion and Shock Waves, 2016, 36(5): 617–624. doi: 10.11883/1001-1455(2016)05-0617-08 [4] 冯晓伟, 卢永刚, 李永泽. 飞机目标在爆炸冲击波作用下的毁伤效应评估方法 [J]. 高压物理学报, 2019, 33(4): 045101. doi: 10.11858/gywlxb.20180687FENG X W, LU Y G, LI Y Z. Damage assessment method of aircraft targets under blast wave [J]. Chinese Journal of High Pressure Physics, 2019, 33(4): 045101. doi: 10.11858/gywlxb.20180687 [5] 王树山, 张静骁, 王传昊, 等. 水中爆炸冲击波对靶体结构的毁伤准则研究 [J]. 火炸药学报, 2020, 43(3): 262–270. doi: 10.14077/j.issn.1007-7812.201909015WANG S S, ZHANG J X, WANG C H, et al. Damage criterion of underwater explosion shock wave on target [J]. Chinese Journal of Explosives & Propellants, 2020, 43(3): 262–270. doi: 10.14077/j.issn.1007-7812.201909015 [6] 胡宏伟, 宋浦, 王建灵, 等. 炸药水中爆炸冲击因子的新型计算方法 [J]. 爆炸与冲击, 2014, 34(1): 11–16. doi: 10.11883/1001-1455(2014)01-0011-06HU H W, SONG P, WANG J L, et al. A new calculation method for shock factor of underwater explosion [J]. Explosion and Shock Waves, 2014, 34(1): 11–16. doi: 10.11883/1001-1455(2014)01-0011-06 [7] 卢熹, 王树山, 王新颖. 水中爆炸对鱼雷壳体的毁伤准则和判据研究 [J]. 兵工学报, 2016, 37(8): 1469–1475. doi: 10.3969/j.issn.1000-1093.2016.08.019LU X, WANG S S, WANG X Y. Research on damage criterion of torpedo shell subjected to underwater explosive shock waves [J]. Acta Armamentarii, 2016, 37(8): 1469–1475. doi: 10.3969/j.issn.1000-1093.2016.08.019 [8] 薛再清, 辛春亮, 梁海东, 等. 水下爆炸冲击波及其作用判据研究 [J]. 防护工程, 2016, 38(3): 48–52.XUE Z Q, XIN C L, LIANG H D, et al. Study on shock-wave of underwater explosion and effect criterion [J]. Protective Engineering, 2016, 38(3): 48–52. [9] ZHANG W Z, JIANG W K. An improved shock factor to evaluate the shock environment of small-sized structures subjected to underwater explosion [J]. Shock and Vibration, 2015, 2015: 451583. doi: 10.1155/2015/451583 [10] 姚熊亮, 曹宇, 郭君, 等. 一种用于水面舰船的水下爆炸冲击因子 [J]. 哈尔滨工程大学学报, 2007, 28(5): 501–509. doi: 10.3969/j.issn.1006-7043.2007.05.004YAO X L, CAO Y, GUO J, et al. Research on the response of warships to impulsive factor of underwater explosions [J]. Journal of Harbin Engineering University, 2007, 28(5): 501–509. doi: 10.3969/j.issn.1006-7043.2007.05.004 [11] 罗泽立, 周章涛, 毛海斌, 等. 水下爆炸强冲击波与平板结构相互作用的理论分析方法 [J]. 高压物理学报, 2017, 31(4): 443–452. doi: 10.11858/gywlxb.2017.04.013LUO Z L, ZHOU Z T, MAO H B, et al. Theoretical analysis of the interaction between the plate structure and strong shock wave in underwater explosion [J]. Chinese Journal of High Pressure Physics, 2017, 31(4): 443–452. doi: 10.11858/gywlxb.2017.04.013 [12] 王树山. 终点效应学 [M]. 2版. 北京: 科学出版社, 2019: 33−34.WANG S S. Terminal effects [M]. 2nd ed. Beijing: Science Press, 2019: 33−34. [13] TAYLOR G I. The pressure and impulse of submarine explosion waves on plates [C]//The Scientific Papers Sir Geoffrey Ingram Taylor. Cambridge: Cambridge University Press, 1963: 287−303. [14] 李旭东, 尹建平, 杜志鹏, 等. 多次水下爆炸钢制圆板应变与挠度增长规律分析 [J]. 振动与冲击, 2020, 39(5): 131–136. doi: 10.13465/j.cnki.jvs.2020.05.017LI X D, YIN J P, DU Z P, et al. Growth law analysis for strain and deflection of steel circular plates subjected to multiple underwater explosions [J]. Journal of Vibration and Shock, 2020, 39(5): 131–136. doi: 10.13465/j.cnki.jvs.2020.05.017 [15] 张斐, 张春辉, 张磊, 等. 多次水下爆炸作用下钢板动态响应数值模拟 [J]. 中国舰船研究, 2019, 14(6): 122–129. doi: 10.19693/j.issn.1673-3185.01431ZHANG F, ZHANG C H, ZHANG L, et al. Numerical simulation of dynamic response of steel plate subjected to multiple underwater explosions [J]. Chinese Journal of Ship Research, 2019, 14(6): 122–129. doi: 10.19693/j.issn.1673-3185.01431 [16] COLE R H. Underwater explosions [M]. New Jersey: Princeton University Press, 1948. [17] 张显丕, 刘建湖, 潘建强, 等. 基于效应靶的装药水下近场爆炸威力评估方法 [J]. 兵工学报, 2016, 37(8): 1430–1435. doi: 10.3969/j.issn.1000-1093.2016.08.013ZHANG X P, LIU J H, PAN J Q, et al. An evaluation method for near-field underwater explosion power based on effect target [J]. Acta Armamentarii, 2016, 37(8): 1430–1435. doi: 10.3969/j.issn.1000-1093.2016.08.013 -