基于正交试验法的气泡帷幕削波特性研究

杜明燃; 陈宇航; 陆少锋; 梁进; 李基锐; 王尹军; 王天照; 陈智凡

doi:10.11858/gywlxb.20230684

基于正交试验法的气泡帷幕削波特性研究

doi: 10.11858/gywlxb.20230684

杜明燃^{1, 2,},
陈宇航¹,
陆少锋^{2, 3},
梁进^{2, 3},
李基锐^{2, 3},
王尹军⁴,
王天照¹,
陈智凡¹

1.
安徽理工大学化学工程学院, 安徽淮南　232001
2.
广西壮族自治区水下破岩工程研究中心, 广西南宁　530200
3.
广西新港湾工程有限公司, 广西南宁　530200
4.
矿冶科技集团有限公司, 北京　100160

基金项目: 安徽省高校科学研究项目（KJ2021A0431）；广西重点研发计划(桂科AB22035001)；安徽省自然科学基金（1908085QA33）；安徽理工大学引进人才科研启动基金（11881）

详细信息

作者简介:
杜明燃（1987－），男，博士，副教授，主要从事炸药性能和爆破技术研究. E-mail：dumingranaust@163.com

中图分类号: O389; TD235
计量
- 文章访问数: 262
- HTML全文浏览量: 153
- PDF下载量: 28
出版历程
- 收稿日期: 2023-06-20
- 修回日期: 2023-07-03
- 网络出版日期: 2023-11-26
- 刊出日期: 2023-12-15

Bubble Curtain Clipping Characteristics Based on Orthogonal Test Method

DU Mingran^{1, 2
,},
CHEN Yuhang¹,
LU Shaofeng^{2, 3},
LIANG Jin^{2, 3},
LI Jirui^{2, 3},
WANG Yinjun⁴,
WANG Tianzhao¹,
CHEN Zhifan¹

1.
School of Chemical Engineering, Anhui University of Science and Technology, Huainan 232001, Anhui, China
2.
Research Institution of Underwater Rock-Cutting Engineering at Guangxi Zhuang Autonomous Region, Nanning 530200, Guangxi, China
3.
Guangxi New Harbor Engineering Co., Ltd., Nanning 530200, Guangxi, China
4.
BGRIMM Technology Group, Beijing 100160, China

摘要

摘要: 气泡帷幕削波技术的影响因素众多。为了得到气泡帷幕在实际工程应用中的最佳组合方案，利用AUTODYN软件设计3因素3水平正交试验，以冲击波峰值压力的平均削减率作为评价指标，研究了气泡帷幕层数、气泡帷幕爆心距和药包深度对气泡帷幕削波效果的影响和敏感性。结果表明：气泡帷幕的削波能力存在极大值，此后继续增加气泡帷幕层数反而会出现负效应；气泡帷幕层数对气泡帷幕削波效应的影响最大，药包深度次之，气泡帷幕爆心距的影响最小；气泡帷幕爆心距越小，药包深度越大，气泡帷幕的削波效果越好；气泡帷幕层数为2层，气泡帷幕爆心距为1 m，药包深度为10.5 m时，削波效果最好。
- 水下爆炸 /
- 气泡帷幕 /
- 正交试验 /
- 峰值压力 /
- 敏感性
Abstract: There are many influencing factors of bubble curtain clipping technology. In order to obtain the optimal combination of bubble curtain in engineering applications, AUTODYN software was used to design a three factor and three level orthogonal test. The average reduction rate of shock wave peak pressure was used as an evaluation index to study the influence and sensitivity of the number of bubble curtain layers, the burst distance of the bubble curtain and the charge depth on the bubble curtain clipping effect. The results indicate that the clipping ability of the bubble curtain has a great value, and there will be a negative effect if the number of bubble curtain layers continues to increase. The number of bubble curtain layers has the greatest effect on the clipping ability, followed by the charge depth, and the effect of bubble curtain explosion center distance is the smallest. The smaller the burst center distance and the greater the depth of the packet, the better the bubble curtain chipping effect. When the number of bubble curtain layers is 2, the bubble curtain burst distance is 1 m, and the charge depth is 10.5 m, the clipping effect is the best.
- underwater explosion /
- bubble curtain /
- orthogonal test /
- peak pressure /
- sensitivity

HTML全文

港口码头、水利水电设施和航道疏浚等基础建设工程在我国的发展规划中占据重要地位，然而，其作业环境复杂，难以通过机械拆除手段实现工程目标。水下爆破技术具有适用范围广、价格低廉等优势，能够很好地适应复杂环境，被广泛地应用于工程建设^[1–3]。水下爆炸与其他介质中的爆炸不同，所产生的冲击波具有峰值压力更高、影响范围更广、毁伤效应更强的特点^[4–5]。水下爆破作业过程的主要危害因素有水中冲击波、水面涌浪、爆破飞石等^[6]。其中，水中冲击波具有毁伤作用强、作用时间长的特点^[7]，如何减小水中冲击波的危害是目前水下爆破领域亟待解决的重要问题之一^[8–9]。

气泡帷幕的概念最早由加拿大工程师Adolph提出，用于Oratario水电站的水下爆破，由于其对水中冲击波的削减效果优异，因而得到广泛认可和应用^[10]。为进一步研究气泡帷幕的削波作用，国内学者也开展了相关研究。刘欣等^[11]通过改变气泡帷幕与被保护对象之间的距离，研究了气泡帷幕对水中冲击波衰减效果的影响，发现气泡帷幕与被保护物的间距越小，水中冲击波的削减效果越明显。胡伟才等^[12]利用LS-DYNA软件，分析了气泡帷幕的数量、间距和防护距离对桥梁振动速度削减作用的影响，发现气泡帷幕数量的影响最大，气泡帷幕间距次之，气泡帷幕防护距离的影响最小。李泽华等^[13]研究了气泡与水中冲击波的相互作用过程，发现气泡的强力压缩和破碎是削减水中冲击波的最佳方式。谢达建等^[14]利用LS-DYNA软件建立了水下钻孔爆破模型，结合长江九朝段炸礁工程，分析了气泡帷幕的距离对削波效果的影响，发现在距离被保护对象较近处设置气泡帷幕时的防护效果更好。陆少锋等^[15]研究了不同供风量形成的气泡帷幕对水中冲击波的削减效应，发现气泡帷幕的削波效果随着供风量的增加而增强。

综上所述，气泡帷幕的削波效果与众多因素有关，而气泡帷幕层数、气泡帷幕爆心距和药包深度作为工程实践中的重要因素，对气泡帷幕削波效果的影响以及影响程度尚待深入研究。为此，本研究利用AUTODYN软件建立水下爆炸模型，设计现场试验对模型进行验证，通过3因素3水平正交试验，对气泡帷幕层数、气泡帷幕爆心距和药包深度对削波效果的影响进行敏感性分析，以期获得削波效果最优的组合方案，为气泡帷幕在水下爆炸中的应用提供参考和理论依据。

1. 模型的建立与验证

1.1 材料模型与状态方程

1.1.1 炸药材料模型

炸药采用JWL状态方程^[16]描述，其表达式为

$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{E}{V}$

(1)

式中：p为冲击波压力，A、B、R₁、R₂和ω为JWL状态方程参数，E为炸药内能，V为当前的相对体积。工程中常用2号岩石乳化炸药，具体参数见表1，其中：ρ为密度，D为爆速，p_C-J为C-J爆轰压力。

表 1 炸药的材料参数

Table 1. Material parameters of explosive

ρ/(g·cm⁻³)	D/(km·s⁻¹)	p_C-J/GPa	A/GPa	B/GPa	R₁	R₂	ω	E/(kJ·g⁻¹)
1.30	4.5	9.80	214.4	0.182	4.15	0.95	0.15	4.19

下载: 导出CSV

| 显示表格

1.1.2 水材料模型

水采用多项式状态方程^[16]描述。当水被压缩时（压缩比μ>0），其状态方程为

${p}_{\rm w}={A}_{1}\mu +{A}_{2}\mu +{A}_{3}\mu +({B}_{0}+{B}_{1}\mu ){\rho }_{0}e$

(2)

当水膨胀时（μ<0），其状态方程为

${p_{\rm w}} = {T_1}\mu +{T_2}\mu +{B_0}{\rho _0}e$

(3)

当水既不压缩也不膨胀时（μ=0），其状态方程可简化为

${p}_{\rm w}={B}_{0}{\rho }_{0}e$

(4)

式中：p_w为水的压力；压缩比μ=ρ/ρ₀−1，ρ为水的当前密度，ρ₀为水的初始密度；A₁、A₂、A₃、B₀、B₁、T₁、T₂均为常数，如表2所示；e为水的比内能，e=(p₀ +ρ₀gH)/B₀ρ₀，其中p₀为大气压，g为重力加速度，H为水深（即药包深度）。根据药包深度调节水的比内能e，当H=4.5 m时，水的比内能为519.38 J/kg。

表 2 水的材料参数

Table 2. Material parameters of water

ρ/(g·cm⁻³)	A₁/GPa	A₂/GPa	A₃/GPa	B₀	B₁	T₁/GPa	T₂/GPa	ω
1.000	2.2	9.54	14.57	0.28	0.28	2.2	0	0.150

下载: 导出CSV

| 显示表格

1.1.3 气泡帷幕材料模型

由刘欣等^[11]和谢达建等^[14]的研究可知，数值模拟中可以利用空气层近似替代气泡帷幕。选用AUTODYN材料库中的AIR空气模型和理想气体状态方程^[17]，其表达式为

$p = (\gamma - 1)\frac{{{\rho _{\rm a}}}}{{{\rho }_{\rm a0}}}E_{\rm a}$

(5)

式中：E_a为空气的比内能；γ为绝热指数，取1.4；ρ_a和ρ_a0分别为空气的当前密度和初始密度，ρ_a0取1.225 kg/m³。

1.2 气泡帷幕计算模型

药包在无限水域中爆炸时，球形药包在平面内关于x轴和y轴对称，为了减少运算量，利用AUTODYN建立如图1所示的1/4轴对称计算模型，模型尺寸为20 m×20 m。炸药和水均采用欧拉网格划分，网格大小设为10 mm。炸药选用300 g的2号岩石乳化炸药，起爆点设在药包球心，选用Transmit边界条件，采用mm-mg-ms单位制^[18–19]。通过调节空气层的厚度实现不同气泡帷幕层数的模拟：气泡帷幕层数为1（N=1）时，空气层厚度为10 cm，气泡帷幕层数为2（N=2）时，空气层厚度为20 cm，依此类推。

图 1 数值计算模型示意图

Figure 1. Schematic diagram of the numerical simulation model

下载: 全尺寸图片幻灯片

1.3 模型验证

为验证数值模拟的合理性，在某河道进行了气泡帷幕水下爆炸冲击波测试。试验设置3种工况：工况1为1层气泡帷幕，工况2为2层气泡帷幕，工况3为3层气泡帷幕。每种工况测试2次，共计6次试验。试验区域的水深为9～10 m，平均水流速度为0.1 m/s。传感器和药包均设置在4.5 m水深处，且位于同一水平线上。传感器选用PCB-138A01型压力传感器。冲击波测试仪为4通道Blast-PRO型测试仪，测试过程中，触发电平设置为0.2%，量程设置为10 V，记录时长为0.4 s。炸药选用300 g 2号岩石乳化炸药。气泡帷幕设置在距药包中心6 m处。选取C₁和C₂两个测点，其中，测点C₁位于药包左侧12 m处，测点C₂位于药包右侧12 m处。试验布局如图2所示。现场试验与数值模拟结果列于表3，其中p_max为冲击波峰值压力，δ为冲击波峰值压力削减率。以1层气泡帷幕为例，测点C₁和C₂的压力时程曲线如图3所示。

图 2 气泡帷幕试验布局

Figure 2. Layout of the bubble curtain test

下载: 全尺寸图片幻灯片

表 3 不同气泡帷幕层数下水下爆炸冲击波峰值压力

Table 3. Peak pressure of underwater blast shock wave for different bubble curtain layers

Test No.	N	On-site test			Simulation
		p_max/MPa		δ/%	p_max/MPa		δ/%
		C₁	C₂	δ/%	C₁	C₂	δ/%
1	1	1.883	0.154	91.82	1.750	0.291	83.37
2	1	1.789	0.228	87.26	1.750	0.291	83.37
3	2	1.675	0.157	90.63	1.750	0.116	93.37
4	2	1.486	0.108	92.73	1.750	0.116	93.37
5	3	1.862	0.247	86.73	1.750	0.248	85.83
6	3	1.764	0.196	88.89	1.750	0.248	85.83

下载: 导出CSV

| 显示表格

图 3 设置1层气泡帷幕时测点C₁和C₂处的压力时程曲线

Figure 3. Shock wave pressure-time curves at measurement points C₁ and C₂ for one layer of bubble curtain

下载: 全尺寸图片幻灯片

从表3可以看出：设置1层气泡帷幕时，2次现场试验测得C₁处的峰值压力分别为1.883和1.789 MPa，C₂处的峰值压力分别为0.154和0.228 MPa，削减率分别为91.82%和87.26%；数值模拟得到C₁和C₂处的峰值压力分别为1.750和0.291 MPa，削减率为83.37%。以峰值压力的削减率为误差评判指标，数值模拟与试验结果的相对误差分别为9.20%和4.46%。

设置2层气泡帷幕时，现场试验得到的削减率为90.63%和92.73%，数值模拟得到削减率为93.37%，相对误差为3.02%和0.69%。设置3层气泡帷幕时，现场试验得到的削减率为86.73%和88.89%，数值模拟得到的削减率为85.83%，相对误差为1.04%和3.44%。

由上述分析可知，采用空气层替代本试验工况中的气泡帷幕时，误差较小，说明本数值模型可以很好地模拟实际工况，因此后续正交试验均采用该方法进行模拟运算。

2. 试验设计与分析

为减少试验次数，同时保证试验的可靠性，选用L9(3⁴)正交表，利用AUTODYN软件设计了3因素3水平正交试验方案。选取气泡帷幕层数（N）、气泡帷幕爆心距（D）、药包深度（H）3个因素，分别设为因素A、B、C，其中：气泡帷幕层数分为3种水平，即N取1、2、3；气泡帷幕爆心距有3种水平，分别取1、3、5 m；药包深度有3种水平，分别为2.5、6.5、10.5 m。各因素之间无交互作用。此外，设置未加气泡帷幕的空白对照组。如图1所示，距药包中心6、9和12 m分别设置测点1、测点2和测点3，取冲击波峰值压力削减率δ作为评价指标，正交试验设计因素和因素水平见表4，正交试验方案见表5。

表 4 正交试验设计因素和水平

Table 4. Orthogonal test design factors and levels

Level	Factor
Level	N	D/m	H/m
1	1	1	2.5
2	2	3	6.5
3	3	5	10.5

下载: 导出CSV

| 显示表格

表 5 正交试验方案

Table 5. Orthogonal test protocols

No.	Level			Test programme	No.	Level			Test programme
No.	N	D	H	Test programme	No.	N	D	H	Test programme
1	1	1	1	A₁B₁C₁	6	2	3	1	A₂B₃C₁
2	1	2	2	A₁B₂C₂	7	3	1	3	A₃B₁C₃
3	1	3	3	A₁B₃C₃	8	3	2	1	A₃B₂C₁
4	2	1	2	A₂B₁C₂	9	3	3	2	A₃B₃C₂
5	2	2	3	A₂B₂C₃

下载: 导出CSV

| 显示表格

2.1 空白组结果与分析

利用水的多项式状态方程，可以计算出药包深度为2.5、6.5和10.5 m时，水的比内能分别为449.38、589.38和729.38 J/kg。由此通过修改AUTODYN材料库中水介质的比内能^[20–22]模拟炸药的不同水深环境。在空白对照组中，水下爆炸冲击波峰值压力（p_max）及其到达时间（t_p）如表6所示，压力时程曲线如图4所示。

表 6 空白对照组数据

Table 6. Data of blank control group

No.	H/m	Point	p_max/MPa	t_p/ms	No.	H/m	Point	p_max/MPa	t_p/ms
1	2.5	1	4.108	3.952	3	10.5	1	4.187	3.954
		2	2.456	6.001			2	2.533	6.001
		3	1.531	8.006			3	1.708	8.002
2	6.5	1	4.148	3.954
		2	2.494	6.002
		3	1.688	8.016

下载: 导出CSV

| 显示表格

图 4 无气泡帷幕时冲击波压力时程曲线

Figure 4. Time course curves of shock wave pressure without bubble curtain

下载: 全尺寸图片幻灯片

从表6和图4可以看出：冲击波压力时程曲线的波形一致；在浅水区域，随着水下爆炸药包深度的增加，冲击波到达测点的时间基本保持不变，而冲击波峰值压力增大。以测点1为例，当药包深度为2.5、6.5和10.5 m时，冲击波到达峰值压力的时间分别为3.952、3.954和3.954 ms，峰值压力分别为4.108、4.148和4.187 MPa。

2.2 正交试验结果与分析

按照表5设计正交试验，将所得数据绘制压力时程曲线，如图5所示。可以看出，气泡帷幕层数、气泡帷幕爆心距和药包深度对气泡帷幕的削波效果均有影响。以正交试验1为例：对于空白对照试验1（未设置气泡帷幕，H=2.5 m），测点1、测点2和测点3的冲击波峰值压力分别为4.108、2.456和1.531 MPa；对于正交试验1（N=1，D=1 m，H=2.5 m），测点1、测点2和测点3的冲击波峰值压力分别为0.431、0.314和0.233 MPa，削减率分别为89.51%、87.21%和84.78%，平均削减率为87.17%。同理，求得其余8个试验的冲击波峰值压力削减率，如表7所示。

图 5 有气泡帷幕时的冲击波压力时程曲线

Figure 5. Time course curves of shock wave pressure with bubble curtain

下载: 全尺寸图片幻灯片

表 7 正交试验数据

Table 7. Orthogonal test data

Test No.	Programme	N	D/m	H/m	p_max/MPa			δ/%
Test No.	Programme	N	D/m	H/m	Point 1	Point 2	Point 3	δ/%
1	A₁B₁C₁	1	1	2.5	0.341	0.314	0.295	86.55
2	A₁B₂C₂	1	3	6.5	0.320	0.278	0.261	88.56
3	A₁B₃C₃	1	5	10.5	0.355	0.286	0.264	88.26
4	A₂B₁C₂	2	1	6.5	0.122	0.113	0.106	95.42
5	A₂B₂C₃	2	3	10.5	0.132	0.118	0.107	95.31
6	A₂B₃C₁	2	5	2.5	0.146	0.121	0.109	94.80
7	A₃B₁C₃	3	1	10.5	0.275	0.260	0.250	89.51
8	A₃B₂C₁	3	3	2.5	0.372	0.332	0.308	85.77
9	A₃B₃C₂	3	5	6.5	0.412	0.305	0.274	87.20

下载: 导出CSV

| 显示表格

如表7所示，气泡帷幕防护技术能够很好地削弱水下爆炸冲击波的峰值压力，削减率可达95.42%；在所设置的9组正交试验中，正交试验4（N=2，D=1 m，H=6.5 m）的削减率最大，防护效果最好。设置1、2、3层气泡帷幕时，平均削减率分别为86.55%～88.56%、94.80%～95.42%、85.77%～89.51%。由此可见，气泡帷幕的冲击波峰值压力削减率存在最大值，到达最佳防护效果之后，继续增加气泡帷幕数量反而出现负效应。

另外，以正交试验1、2、3为例，测点1处，D为1、3、5 m时，冲击波到达峰值压力的时间分别为3.507、2.058、0.652 ms，说明气泡帷幕爆心距越小，水下爆炸冲击波到达峰值压力的时间越长。

2.3 极差分析

设K_i表示因素j取第i个水平时的评价指标之和，k_i表示因素j取第i个水平时评价指标的平均值。通过极差R_j分析各个因素对评价指标的影响程度，R_j越大，则表明因素j对评价指标的影响越明显^[23–25]。

设评价指标为冲击波峰值压力的平均削减率δ_ave，根据表7中的正交试验数据，得到极差分析结果，如表8所示，水平与指标的关系如图6所示。

表 8 极差分析结果

Table 8. Results of the variance analysis

Factor	K₁	K₂	K₃	k₁	k₂	k₃	R
N	263.37	285.53	262.48	87.79	95.18	87.49	7.68
D	271.48	269.64	270.26	90.49	89.88	90.09	0.61
H	267.12	271.18	273.08	89.04	90.39	91.03	1.99

下载: 导出CSV

| 显示表格

图 6 水平与指标的关系

Figure 6. Relationship between level and indicators

下载: 全尺寸图片幻灯片

由表8和图6可知，气泡帷幕层数、气泡帷幕爆心距和药包深度的极差分别为7.68、0.61、1.99，说明气泡帷幕层数对气泡帷幕削波效应的影响最大，药包深度次之，气泡帷幕爆心距的影响最小。图6(a)显示，削减率随气泡帷幕层数的增加先增大后减小，进一步验证了当气泡帷幕削波效果达到最大后继续增加气泡帷幕层数反而会降低削波效果的结论。另外，气泡帷幕爆心距越小，药包深度越大，则气泡帷幕削波效果越好。根据表8可以推断，当气泡帷幕层数为2，气泡帷幕爆心距为1 m，药包深度为10.5 m时，削波效果最好。

2.4 方差分析

在方差分析中，取显著性水平P为0.05，查表得F_0.05为19。通过计算因素j的离差平方和S_j、自由度μ_j和均方 ${\bar S_j}$ ，得到因素j的F_j值，与显著性水平F_0.05进行对照，若F_j>F_0.05，则说明因素j对气泡帷幕削波效果的影响显著；反之，则影响不显著。因素j的均方 ${\overline{S}}_{j}$ 越大，其对气泡帷幕削波效果的影响越大。

${S_j} = \frac{1}{n}\sum\limits_{i=1}^n {{K^2_{ij}}} - \frac{1}{h}{\left(\sum\limits_{i = 1}^h {{y_i}} \right)^2}$

(6)

$S = \sum\limits_{i = 1}^h {{y^2_i}} - \frac{1}{h}{\left(\sum\limits_{i = 1}^h {{y_i}} \right)^2}$

(7)

${S_{\rm e}} = S - \sum\limits_{j = 1}^3 {{S_j}}$

(8)

式中：K_ij为第j个因素下第i个水平所对应的试验指标；y_i为第i次试验的试验指标；n为水平数，n=3；h为正交试验次数，h=9；S为总离差平方和；S_e为误差平方和。

${f _{\rm e}} = {f _{\rm t}} - \sum\limits_{j = 1}^3 {{f _j}}$

(9)

式中：f_e为误差项自由度；f_t为总自由度，f_t=h−1；f_j=n−1。

${\bar S_j} = \frac{{{S_j}}}{{{f _j}}}$

(10)

${\bar S_{\rm e}} = \frac{{{S_{\rm e}}}}{{{f _{\rm e}}}}$

(11)

${F_j} = \frac{{{{\bar S}_j}}}{{{{\bar S}_{\rm e}}}}$

(12)

式中： ${\bar S_{\rm e}}$ 为误差项的均方。F统计量服从第一自由度为f_j、第二自由度为f_e的F分布，显著性水平取0.05时，对F统计量进行检验。

利用式(9)～式(17)计算离差平方和S、自由度f、均方 $\bar S$ 和F值^[26–28]，所得结果见表9。

表 9 方差分析结果

Table 9. Results of variance analysis

Sources of variance	S	f	$\bar S$	F	Significance-effect
N	113.68	2	56.84	38.81	Yes
D	3.48	2	1.74	1.19	No
H	6.18	2	3.09	2.11	No
Error	2.93	2	1.46
Sum	126.27	8

下载: 导出CSV

| 显示表格

根据表9中的 $\bar S$ 可以看出，正交试验所选因素对气泡帷幕削波效果的影响由大到小依次为气泡帷幕层数、药包深度、气泡帷幕爆心距。通过方差分析可以看出：气泡帷幕层数的显著性水平F为38.81，大于F_0.05，表明气泡帷幕层数对削波效果有显著影响；药包深度和气泡帷幕爆心距的显著性水平分别为1.19和2.11，均小于F_0.05，表明二者对气泡帷幕削波效果的影响不大，与极差分析所得结论一致。

3. 结　论

基于正交试验极差分析，利用AUTODYN显式有限元分析程序，建立了自由水域内药包爆炸的1/4轴对称计算模型，探究了气泡帷幕层数、气泡帷幕爆心距和药包深度对气泡帷幕削波效应的影响和敏感性，得出以下主要结论。

(1) 气泡帷幕可以有效地削减水中冲击波峰值压力，削减率可达95.42%，从而有效降低水中冲击波对爆破区域附近被保护对象的毁伤作用。

(2) 气泡帷幕的削波效果与气泡帷幕层数、气泡帷幕爆心距和药包深度均有关。气泡帷幕爆心距越小，药包深度越大，则气泡帷幕的削波效果越好；而气泡帷幕层数与气泡帷幕削波效果并不成正比，当气泡帷幕层数为1、2、3时，冲击波峰值压力的平均削减率分别为87.79%、95.17%和87.49%，在实际工程中应合理选择气泡帷幕层数。

(3) 正交试验分析显示：气泡帷幕层数对气泡帷幕削波效应的影响最大，药包深度次之，气泡帷幕爆心距的影响最小；当气泡帷幕层数为2，气泡帷幕爆心距为1 m，药包深度为10.5 m时，削波效果最好。

图 1 数值计算模型示意图

Figure 1. Schematic diagram of the numerical simulation model

下载: 全尺寸图片幻灯片

图 2 气泡帷幕试验布局

Figure 2. Layout of the bubble curtain test

下载: 全尺寸图片幻灯片

图 3 设置1层气泡帷幕时测点C₁和C₂处的压力时程曲线

Figure 3. Shock wave pressure-time curves at measurement points C₁ and C₂ for one layer of bubble curtain

下载: 全尺寸图片幻灯片

图 4 无气泡帷幕时冲击波压力时程曲线

Figure 4. Time course curves of shock wave pressure without bubble curtain

下载: 全尺寸图片幻灯片

图 5 有气泡帷幕时的冲击波压力时程曲线

Figure 5. Time course curves of shock wave pressure with bubble curtain

下载: 全尺寸图片幻灯片

图 6 水平与指标的关系

Figure 6. Relationship between level and indicators

下载: 全尺寸图片幻灯片

表 1 炸药的材料参数

Table 1. Material parameters of explosive

ρ/(g·cm⁻³)	D/(km·s⁻¹)	p_C-J/GPa	A/GPa	B/GPa	R₁	R₂	ω	E/(kJ·g⁻¹)
1.30	4.5	9.80	214.4	0.182	4.15	0.95	0.15	4.19

下载: 导出CSV

表 2 水的材料参数

Table 2. Material parameters of water

ρ/(g·cm⁻³)	A₁/GPa	A₂/GPa	A₃/GPa	B₀	B₁	T₁/GPa	T₂/GPa	ω
1.000	2.2	9.54	14.57	0.28	0.28	2.2	0	0.150

下载: 导出CSV

表 3 不同气泡帷幕层数下水下爆炸冲击波峰值压力

Table 3. Peak pressure of underwater blast shock wave for different bubble curtain layers

Test No.	N	On-site test			Simulation
		p_max/MPa		δ/%	p_max/MPa		δ/%
		C₁	C₂	δ/%	C₁	C₂	δ/%
1	1	1.883	0.154	91.82	1.750	0.291	83.37
2	1	1.789	0.228	87.26	1.750	0.291	83.37
3	2	1.675	0.157	90.63	1.750	0.116	93.37
4	2	1.486	0.108	92.73	1.750	0.116	93.37
5	3	1.862	0.247	86.73	1.750	0.248	85.83
6	3	1.764	0.196	88.89	1.750	0.248	85.83

下载: 导出CSV

表 4 正交试验设计因素和水平

Table 4. Orthogonal test design factors and levels

Level	Factor
Level	N	D/m	H/m
1	1	1	2.5
2	2	3	6.5
3	3	5	10.5

下载: 导出CSV

表 5 正交试验方案

Table 5. Orthogonal test protocols

No.	Level			Test programme	No.	Level			Test programme
No.	N	D	H	Test programme	No.	N	D	H	Test programme
1	1	1	1	A₁B₁C₁	6	2	3	1	A₂B₃C₁
2	1	2	2	A₁B₂C₂	7	3	1	3	A₃B₁C₃
3	1	3	3	A₁B₃C₃	8	3	2	1	A₃B₂C₁
4	2	1	2	A₂B₁C₂	9	3	3	2	A₃B₃C₂
5	2	2	3	A₂B₂C₃

下载: 导出CSV

表 6 空白对照组数据

Table 6. Data of blank control group

No.	H/m	Point	p_max/MPa	t_p/ms	No.	H/m	Point	p_max/MPa	t_p/ms
1	2.5	1	4.108	3.952	3	10.5	1	4.187	3.954
		2	2.456	6.001			2	2.533	6.001
		3	1.531	8.006			3	1.708	8.002
2	6.5	1	4.148	3.954
		2	2.494	6.002
		3	1.688	8.016

下载: 导出CSV

表 7 正交试验数据

Table 7. Orthogonal test data

Test No.	Programme	N	D/m	H/m	p_max/MPa			δ/%
Test No.	Programme	N	D/m	H/m	Point 1	Point 2	Point 3	δ/%
1	A₁B₁C₁	1	1	2.5	0.341	0.314	0.295	86.55
2	A₁B₂C₂	1	3	6.5	0.320	0.278	0.261	88.56
3	A₁B₃C₃	1	5	10.5	0.355	0.286	0.264	88.26
4	A₂B₁C₂	2	1	6.5	0.122	0.113	0.106	95.42
5	A₂B₂C₃	2	3	10.5	0.132	0.118	0.107	95.31
6	A₂B₃C₁	2	5	2.5	0.146	0.121	0.109	94.80
7	A₃B₁C₃	3	1	10.5	0.275	0.260	0.250	89.51
8	A₃B₂C₁	3	3	2.5	0.372	0.332	0.308	85.77
9	A₃B₃C₂	3	5	6.5	0.412	0.305	0.274	87.20

下载: 导出CSV

表 8 极差分析结果

Table 8. Results of the variance analysis

Factor	K₁	K₂	K₃	k₁	k₂	k₃	R
N	263.37	285.53	262.48	87.79	95.18	87.49	7.68
D	271.48	269.64	270.26	90.49	89.88	90.09	0.61
H	267.12	271.18	273.08	89.04	90.39	91.03	1.99

下载: 导出CSV

表 9 方差分析结果

Table 9. Results of variance analysis

Sources of variance	S	f	$\bar S$	F	Significance-effect
N	113.68	2	56.84	38.81	Yes
D	3.48	2	1.74	1.19	No
H	6.18	2	3.09	2.11	No
Error	2.93	2	1.46
Sum	126.27	8

下载: 导出CSV

参考文献(28)

[1]	高明涛, 李昕, 周晶. 水下钻孔爆破水中冲击波的数值模拟研究 [J]. 水电能源科学, 2009, 27(4): 138–141. doi: 10.3969/j.issn.1000-7709.2009.04.043 GAO M T, LI X, ZHOU J. Numerical simulation of shock wave in water of underwater drilling blasting [J]. Water Resources and Power, 2009, 27(4): 138–141. doi: 10.3969/j.issn.1000-7709.2009.04.043
[2]	杨建, 刘静, 张登泰, 等. 气泡帷幕对港口水域爆破波的削弱 [J]. 中国水运, 2020, 20(10): 90–92. YANG J, LIU J, ZHANG D T, et al. Weakening of blast waves in port waters by air bubble curtains [J]. China Water Transport, 2020, 20(10): 90–92.
[3]	司剑峰. 深水钻孔爆破的冲击波衰减规律及防护研究 [D]. 武汉: 武汉科技大学, 2021. SI J F. Research on attenuation law of shock wave and protection in deep-water drilling and blasting [D]. Wuhan: Wuhan University of Science and Technology, 2021.
[4]	俞统昌, 王晓峰, 王建灵. 炸药的水下爆炸冲击波性能 [J]. 含能材料, 2003, 11(4): 182–186. doi: 10.3969/j.issn.1006-9941.2003.04.002 YU T C, WANG X F, WANG J L. Underwater shockwave performance of explosives [J]. Energetic Materials, 2003, 11(4): 182–186. doi: 10.3969/j.issn.1006-9941.2003.04.002
[5]	伍俊, 杨益, 庄铁栓. 水中爆炸作用机理及毁伤效应研究综述 [J]. 火炸药学报, 2016, 39(1): 1–13. doi: 10.14077/j.issn.1007-7812.2016.01.001 WU J, YANG Y, ZHUANG T S. A review of research on action mechanism and damage effect of underwater explosion [J]. Chinese Journal of Explosives & Propellants, 2016, 39(1): 1–13. doi: 10.14077/j.issn.1007-7812.2016.01.001
[6]	陈岩武, 孙远翔, 王成. 水下爆炸载荷下舰船双层底部结构的毁伤特性 [J]. 兵工学报, 2023, 44(3): 670–681. doi: 10.12382/bgxb.2022.0390 CHEN Y W, SUN Y X, WANG C. Damage characteristics of ship’s double bottom structure subjected to underwater explosion [J]. Acta Armamentarii, 2023, 44(3): 670–681. doi: 10.12382/bgxb.2022.0390
[7]	赵为, 梁作民, 郭成喜. 水下近场爆破安全控制 [J]. 水运工程, 2007, 33(9): 159–164. doi: 10.3969/j.issn.1002-4972.2007.09.038 ZHAO W, LIANG Z M, GUO C X. Safety control of underwater near-site explosion [J]. Port & Waterway Engineering, 2007, 33(9): 159–164. doi: 10.3969/j.issn.1002-4972.2007.09.038
[8]	张轶凡, 刘亮涛, 王金相, 等. 水下爆炸冲击波和气泡载荷对典型圆柱壳结构的毁伤特性 [J]. 兵工学报, 2023, 44(2): 345–359. doi: 10.12382/bgxb.2021.0598 ZHANG Y F, LIU L T, WANG J X, et al. Damage characteristics of underwater explosion shock wave and bubble load on typical cylindrical shell structure [J]. Acta Armamentarii, 2023, 44(2): 345–359. doi: 10.12382/bgxb.2021.0598
[9]	王高辉, 高政, 卢文波, 等. 考虑初始应力的混凝土重力坝水下爆炸毁伤特性研究 [J]. 振动与冲击, 2022, 41(11): 133–140. doi: 10.13465/j.cnki.jvs.2022.11.017 WANG G H, GAO Z, LU W B, et al. Damage characteristics of underwater explosion of concrete gravity dam considering initial stress [J]. Journal of Vibration and Shock, 2022, 41(11): 133–140. doi: 10.13465/j.cnki.jvs.2022.11.017
[10]	张志波, 李春军, 李红勇, 等. 气泡帷幕在水下爆破减震工程中的应用 [J]. 爆破, 2003, 20(2): 75–76, 89. doi: 10.3963/j.issn.1001-487X.2003.02.028 ZHANG Z B, LI C J, LI H Y, et al. Application of air bubble purdah in the damping measure in the underwater blasting [J]. Blasting, 2003, 20(2): 75–76, 89. doi: 10.3963/j.issn.1001-487X.2003.02.028
[11]	刘欣, 顾文彬, 陈学平. 气泡帷幕对水中冲击波衰减特性的数值模拟研究 [J]. 爆破, 2015, 32(3): 79–84. doi: 10.3963/j.issn.1001-487X.2015.03.014 LIU X, GU W B, CHEN X P. Numerical simulation study of attenuation characteristics of water shock wave under bubble curtain [J]. Blasting, 2015, 32(3): 79–84. doi: 10.3963/j.issn.1001-487X.2015.03.014
[12]	胡伟才, 吴立, 舒利, 等. 不同设置方式下气泡帷幕对水中冲击波衰减特性的影响 [J]. 科学技术与工程, 2018, 18(17): 33–38. doi: 10.3969/j.issn.1671-1815.2018.17.006 HU W C, WU L, SHU L, et al. Influence of water shock wave on attenuation characteristics under bubble curtain with different settings [J]. Science Technology and Engineering, 2018, 18(17): 33–38. doi: 10.3969/j.issn.1671-1815.2018.17.006
[13]	李泽华, 白春华, 刘庆明, 等. 气泡帷幕减弱水中冲击波强度的研究 [J]. 中国安全科学学报, 1999, 10(5): 72–76. doi: 10.16265/j.cnki.issn1003-3033.1999.05.015 LI Z H, BAI C H, LIU Q M, et al. Study on weakening the shock wave in water by bubble heavy curtain [J]. China Safety Science Journal, 1999, 10(5): 72–76. doi: 10.16265/j.cnki.issn1003-3033.1999.05.015
[14]	谢达建, 吴立, 洪江, 等. 气泡帷幕对水下爆破冲击波的削弱作用研究 [J]. 人民长江, 2018, 49(8): 72–77. doi: 10.16232/j.cnki.1001-4179.2018.08.014 XIE D J, WU L, HONG J, et al. Study on weakening effect of bubble curtain on water shock wave in underwater blasting [J]. Yangtze River, 2018, 49(8): 72–77. doi: 10.16232/j.cnki.1001-4179.2018.08.014
[15]	陆少锋, 梁进, 覃才勇, 等. 供风量对水下爆炸冲击波气泡帷幕消波效应的影响 [J]. 工程爆破, 2022, 28(5): 143–148. doi: 10.19931/j.EB.20210013 LU S F, LIANG J, QIN C Y, et al. Influence of air supply rate on wave attenuation effect of bubble curtain for underwater explosion shock wave [J]. Engineering Blasting, 2022, 28(5): 143–148. doi: 10.19931/j.EB.20210013
[16]	胡亮亮, 黄瑞源, 李世超, 等. 水下爆炸冲击波数值仿真研究 [J]. 高压物理学报, 2020, 34(1): 015102. doi: 10.11858/gywlxb.20190773 HU L L, HUANG R Y, LI S C, et al. Shock wave simulation of underwater explosion [J]. Chinese Journal of High Pressure Physics, 2020, 34(1): 015102. doi: 10.11858/gywlxb.20190773
[17]	刘世聪, 王秋生, 娄浩然. 装药深度及空气域尺寸对水下爆炸的影响分析 [J]. 水下无人系统学报, 2019, 27(6): 664–672. doi: 10.11993/j.issn.2096-3920.2019.06.010 LIU S C, WANG Q S, LOU H R. Effects of charge depth and air domain size on underwater explosion [J]. Journal of Unmanned Undersea Systems, 2019, 27(6): 664–672. doi: 10.11993/j.issn.2096-3920.2019.06.010
[18]	黄洪, 卢熹, 王健. 柱形装药水下爆炸近场冲击波数值仿真 [J]. 水下无人系统学报, 2021, 29(4): 471–476. doi: 10.11993/j.issn.2096-3920.2021.04.015 HUANG H, LU X, WANG J. Near field shock wave numerical simulation of cylindrical charge underwater explosion [J]. Journal of Unmanned Undersea Systems, 2021, 29(4): 471–476. doi: 10.11993/j.issn.2096-3920.2021.04.015
[19]	鲁天龙. 气泡帷幕周围流场运动特性数值模拟研究 [D]. 长沙: 长沙理工大学, 2020. LU T L. Numerical simulation of flow field motion characteristics around bubble curtain [D]. Changsha: Changsha University of Technology, 2020.
[20]	王思, 胡晶, 张雪东, 等. 不同水深水下爆炸数值及离心试验研究 [J]. 哈尔滨工业大学学报, 2020, 52(6): 78–84. doi: 10.11918/202001082 WANG S, HU J, ZHANG X D, et al. Numerical analysis and centrifugal test of underwater explosion effect at different water depths [J]. Journal of Harbin Institute of Technology, 2020, 52(6): 78–84. doi: 10.11918/202001082
[21]	盛振新, 刘荣忠, 郭锐. 壳体厚度和爆炸深度对水下爆炸冲击波的影响 [J]. 火炸药学报, 2011, 34(3): 45–47, 64. doi: 10.3969/j.issn.1007-7812.2011.03.012 SHENG Z X, LIU R Z, GUO R. Effect of shell thickness and explosion depth on underwater explosive shock wave [J]. Chinese Journal of Explosives & Propellants, 2011, 34(3): 45–47, 64. doi: 10.3969/j.issn.1007-7812.2011.03.012
[22]	陈兴, 周兰伟, 李福明, 等. 爆炸深度对装药水下载荷的影响 [J]. 兵器装备工程学报, 2021, 42(8): 79–84. doi: 10.11809/bqzbgcxb2021.08.013 CHEN X, ZHOU L W, LI F M, et al. Research on influence of depth of charge on underwater explosion load [J]. Journal of Ordnance Equipment Engineering, 2021, 42(8): 79–84. doi: 10.11809/bqzbgcxb2021.08.013
[23]	尹岳降, 李瑞泽, 陈明, 等. 基于正交试验法的爆破块度分布影响因素敏感性分析 [J]. 爆破, 2019, 36(4): 37–42. doi: 10.3963/j.issn.1001-487X.2019.04.005 YIN Y J, LI R Z, CHEN M, et al. Sensitivity analysis of influencing factors of blasting fragmentation distribution based on orthogonal experiment method [J]. Blasting, 2019, 36(4): 37–42. doi: 10.3963/j.issn.1001-487X.2019.04.005
[24]	张世豪, 韩晶, 焦国太, 等. 单个装药混凝土爆破毁伤效果敏感性分析 [J]. 工程爆破, 2014, 20(3): 5–9. doi: 10.3969/j.issn.1006-7051.2014.03.002 ZHANG S H, HAN J, JIAO G T, et al. Sensitivity analysis of blasting damage effect in concrete based on single charge [J]. Engineering Blasting, 2014, 20(3): 5–9. doi: 10.3969/j.issn.1006-7051.2014.03.002
[25]	史秀志, 王怀勇, 刘金明, 等. 基于粗糙集的影响爆破振动特征参量因素的敏感性分析 [J]. 爆破器材, 2009, 38(2): 1–4, 7. doi: 10.3969/j.issn.1001-8352.2009.02.001. SHI X Z, WANG H Y, LIU J M, et al. Sensitivity analysis of the factors impact of blasting vibration characteristic parameters based on rough set theory [J]. Explosive Materials, 2009, 38(2): 1–4, 7. doi: 10.3969/j.issn.1001-8352.2009.02.001.
[26]	胡玉林, 姜俊杰, 罗健琛, 等. 基于多因素方差分析探究灌水施肥量对番茄产量的影响 [J]. 浙江农业科学, 2023, 64(5): 1042–1045. doi: 10.16178/j.issn.0528-9017.20230126 HU Y L, JIANG J J, LUO J C, et al. Exploring the effect of irrigation and fertilizer application on tomato yield based on multi-factor analysis of variance [J]. Journal of Zhejiang Agricultural Sciences, 2023, 64(5): 1042–1045. doi: 10.16178/j.issn.0528-9017.20230126
[27]	刘世涛, 郝兵元, 杨冉, 等. 基于方差分析法坚硬顶板下矿压显现分析 [J]. 煤炭技术, 2022, 41(7): 20–23. doi: 10.13301/j.cnki.ct.2022.07.005 LIU S T, HAO B Y, YANG R, et al. Analysis of mineral pressure under rigid top plate based on variance analysis [J]. Coal Technology, 2022, 41(7): 20–23. doi: 10.13301/j.cnki.ct.2022.07.005
[28]	王志建, 龙顺忠, 李颖宏. 基于正交试验的感应控制参数组合优化 [J]. 浙江大学学报(工学版), 2023, 57(6): 1128–1136. doi: 10.3785/j.issn.1008-973X.2023.06.008 WANG Z J, LONG S Z, LI Y H. Combination optimization of induction control parameters based on orthogonal test [J]. Journal of Zhejiang University (Engineering Science), 2023, 57(6): 1128–1136. doi: 10.3785/j.issn.1008-973X.2023.06.008

施引文献

资源附件(0)

访问统计

点击查看大图

图(6) / 表(9)

计量

文章访问数: 262
HTML全文浏览量: 153
PDF下载量: 28

1. 模型的建立与验证
1.1 材料模型与状态方程
1.2 气泡帷幕计算模型
1.3 模型验证
2. 试验设计与分析
2.1 空白组结果与分析
2.2 正交试验结果与分析
2.3 极差分析
2.4 方差分析
3. 结　论

基于正交试验法的气泡帷幕削波特性研究

doi: 10.11858/gywlxb.20230684

作者简介: 杜明燃（1987－），男，博士，副教授，主要从事炸药性能和爆破技术研究. E-mail：dumingranaust@163.com

计量

出版历程

Bubble Curtain Clipping Characteristics Based on Orthogonal Test Method

1. 模型的建立与验证

1.1 材料模型与状态方程

1.1.1 炸药材料模型

1.1.2 水材料模型

1.1.3 气泡帷幕材料模型

1.2 气泡帷幕计算模型

1.3 模型验证

2. 试验设计与分析

2.1 空白组结果与分析

2.2 正交试验结果与分析

2.3 极差分析

2.4 方差分析

3. 结 论

计量

出版历程

目录

1. 模型的建立与验证

1.1 材料模型与状态方程

1.2 气泡帷幕计算模型

1.3 模型验证

2. 试验设计与分析

2.1 空白组结果与分析

2.2 正交试验结果与分析

2.3 极差分析

2.4 方差分析

3. 结 论

作者简介:
杜明燃（1987－），男，博士，副教授，主要从事炸药性能和爆破技术研究. E-mail：dumingranaust@163.com

3. 结　论

3. 结　论