Influence of Plexiglass-Air Interlayer Structure on Foundation Vibration of Small Pool Underwater Explosion
-
摘要: 为了研究有机玻璃-空气层结构对水下爆炸振动的影响, 采用NUBOX-6016爆破振动测试仪监测小水池(直径5.5 m、高3.62 m)中该结构下水下爆炸引起的地基振动信号, 研究不同空气层厚度对最大振动速度的影响; 基于Matlab软件编写相关程序对测得的振动信号进行Hilbert-Huang变换(HHT), 分析不同空气层厚度对振动信号全局频率的影响。结果表明:对于有机玻璃-空气层结构, 随着空气层厚度的增加, 最大振动速度呈先减小后增大的趋势, 当空气层厚度为120 mm时隔振效果最佳; 通过对振动信号的HHT分析得到全局频率所对应的幅值, 5~15 Hz低频区间段的幅值衰减较明显, 且振动作用时间缩短, 能够有效防止水下爆炸与建/构筑物之间产生共振现象。所得试验结果及分析对水下爆破工程防护及军事舰艇防雷仓结构设计等具有一定参考价值。
-
关键词:
- 水下爆炸 /
- 有机玻璃-空气层结构 /
- 爆炸振动 /
- Hilbert-Huang变换
Abstract: To study the influence of the plexiglass-air interlayer structure on underwater explosion vibration, the foundation vibration signals for this structure caused by underwater explosion in the small pool (5.5 m in diameter and 3.62 m in height), the influence of different air interlayer thicknesses on the maximum vibration velocity was studied, HHT (Hilbert-Huang Transform) was used to analyze the vibration test signals by writing relevant programs based on Matlab software, studying the influence of different air interlayer thickness on the global frequency of vibration signals.Experimental results show that:under the condition of plexiglass-air interlayer structure, with the increase of air interlayer thickness, the maximum vibration velocity decreases first and then increases, and the vibration isolation effect is the best when the thickness is 120 mm.The amplitude corresponding to the global frequency can be obtained by means of HHT analysis, the amplitude attenuation of the 5-15 Hz low frequency section is obvious, the action time is shortened, which can effectively prevent the resonance phenomenon between the building and structure.The experimental results and analysis have certain theoretical reference value for underwater blasting engineering protection and structural design of lightning protection bunker for military ship. -
随着海洋资源的开发与水下爆破技术的发展, 水下爆破安全防护问题逐渐受到人们的重视。在陆地岩石中爆破, 地震波能量一般不足爆破总能量的10%;而在水下介质中爆破时, 地震波能量可达20%。水下爆破时结构的振动效应较陆地岩石中更复杂, 一般是低频波与高频波共同作用的结果, 其主要特征是振动频带宽、高频丰富、振幅大[1-3]。在浅水域近距离处开展水下钻孔爆破工程时, 通过采用水下预裂爆破、气泡帷幕防护、空气层防护、多段延时起爆等措施削减地震波和水中冲击波[4]。由于空气的波阻抗比水小很多, 当水下冲击波传播到空气与水的交界面时, 水下冲击波将在介质交界面处发生反射, 其透射系数较小, 少部分能量透射到空气介质中, 可有效缓冲水下冲击波和衰减地震能量[5-6]。樊自建等[7]、贾虎等[8-9]、寇晓枫等[10]开展实验研究了空气隔层对水中冲击波的衰减效果, 得出结论:空气隔层能有效衰减冲击波峰值压力、冲击波冲量及总能量, 对保护对象有较好的防护效果。姚熊亮等[11-12]、刘欣等[13]建立了多介质水下爆炸数值模型进行模拟计算, 其研究结果表明, 空气层结构不仅可以防护接触爆炸产生的冲击波, 而且对非接触爆炸也有很好的防护效果。
本课题组前期针对不同段别雷管的水下爆炸振动信号进行了Hilbert-Huang变换(Hilbert-Huang Transform, HHT)分析[14], 但未进行有机玻璃-空气层结构对水下爆炸振动影响的研究, 且国内外相关文献在此方面报道甚少。基于此考虑, 本工作在小水池中进行有机玻璃-空气层结构约束下的水下爆炸振动测试, 研究空气层厚度对水下爆炸最大振动速度的影响, 通过HHT分析其地基振动信号, 探究空气层厚度对水下爆炸振动全局频率的影响, 为爆炸水池减振、隔振以及减弱爆炸振动对建筑物的破坏作用提供理论支撑。
1. 振动信号HHT分析理论
HHT法适用于处理非线性、非平稳信号, 主要由经验模态分解(Empirical Mode Decomposition, EMD)和Hilbert变换两部分组成[15-17]。相对于小波分析, EMD分解克服了小波基选择问题, 避免了由于基函数不同所带来的影响。另外, 原始信号经EMD分解后得到的固有模态函数(Intrinsic Mode Function, IMF)具有不同特征尺度, 该过程是自适应的、高效的。通过EMD分解得到的各IMF分量经Hilbert变换后, 得到各IMF分量的瞬时频谱, 综合所有IMF分量的瞬时频谱得到Hilbert谱。
爆炸振动原始信号中包含3个方向的信号, 依据《爆破安全规程》(GB 6722-2014)中的13.2.2条规定, 选取振动速度值最大的垂向振动信号(记为X(t))进行分析。基于HHT法对爆炸振动信号进行EMD分解以及Hilbert变换, 过程如下。
原始信号X(t)经n次EMD分解, 可表示为
X(t)=n∑i=1Ci(t)+rn(t) (1) 式中: Ci(t)为经EMD分解得到的第i阶IMF分量, rn(t)为原始信号减去n个IMF分量之和的残差。
对得到的IMF分量做Hilbert变换
H[Ci(t)]=1πPV∫∞−∞Ci(t′)t−t′dt′ (2) 式中:PV表示柯西主值。
构造解析信号Z(t)
Z(t)=Ci(t)+jH[Ci(t)]=ai(t)ejφi(t) (3) 式中:ai(t)为幅值函数, φi(t)为相位函数。
ai(t)=√C2i(t)+H2[Ci(t)] (4) φi(t)=tan−1H[Ci(t)]Ci(t) (5) 原始信号分解后可表示为
X(t)=Ren∑i=1ai(t)ejφi(t)=Ren∑i=1ai(t)e∫ωi(t)dt (6) Hilbert谱可表示为
H(ω,t)=Ren∑i=1ai(t)e∫ωi(t)dt (7) 式中:ω表示角频率。
根据以上理论, 采用Matlab编写程序代码, 加载相应模块即可作出相关Hilbert谱。
2. 有机玻璃-空气层结构条件下小水池水下爆炸振动试验
2.1 试验装置与仪器
设计了4种类型的有机玻璃-空气层结构, 空气层厚度(直径)分别为40、80、120和160 mm。中心装药管为壁厚1 mm的有机玻璃管, 外层为壁厚2 mm的有机玻璃壳体, 装药中心与复合结构中心在同一高度, 如图 1所示。
采用NUBOX-6016型爆破振动测试仪记录振动信号, 其传感器使用TP-3V-10144型号, 采用单次采样方式, 采样频率和时间分别为5 kHz、2 s, 延时点数-256, 触发方式为通道出窗触发, 触发条件为0.05 cm/s, 耦合方式为直流耦合。
2.2 试验系统
在直径为5.5 m、高3.62 m的小型爆炸水池中进行水下爆炸试验。入水深度的选取是基于消除爆炸水池边界效应的考虑。一般装药置于水池中心总深度的2/3处时, 来自水面和池底的边界效应可互相抵消[18], 因此装药中心的入水深度选取为2.4 m。选用3台测振仪, 分别布设在距爆炸水池筒体外壁1、2、3 m处(地基振动速度在3 m处已衰减到较小数值, 而1、2、3 m处测点为等间距的典型观测点, 分别标记为1#、2#、3#)。试验系统如图 2所示。
2.3 试验结果与分析
根据《爆破安全规程》(GB 6722-2014)规定, 选取x、y、z 3个方向中振动速度最大的z方向(与地面垂直方向)振动信号进行分析。测试得到不同距离、不同厚度空气层条件下的水下爆炸振动参数, 如表 1所示。为了探究有机玻璃壳体对水下爆炸振动速度的影响, 进行填充水和空气层的对比试验。表 2为不同距离、不同厚度水隔层情况下的水下爆炸振动参数。
表 1 不同距离、不同厚度空气层对水下爆炸振动的影响Table 1. Influence of distances and thickness of air interlayer on underwater explosion vibrationDistance/m Direction Main vibration frequency/Hz Maximum vibration velocity/(cm·s-1) ∅40 mm ∅80 mm ∅120 mm ∅160 mm ∅40 mm ∅80 mm ∅120 mm ∅160 mm 1 z 56.763 56.763 56.763 56.763 0.086 0.070 0.054 0.070 2 z 56.763 56.763 56.763 44.556 0.118 0.097 0.082 0.101 3 z 56.763 56.763 56.763 56.763 0.080 0.071 0.061 0.086 表 2 不同距离、不同厚度水隔层对水下爆炸振动的影响Table 2. Influence of distances and thickness of water interlayer on underwater explosion vibrationDistance/m Direction Main vibration frequency/Hz Maximum vibration velocity/(cm·s-1) ∅40 mm ∅80 mm ∅120 mm ∅160 mm ∅40 mm ∅80 mm ∅120 mm ∅160 mm 1 z 62.256 62.256 62.256 56.763 0.086 0.169 0.184 0.185 2 z 62.256 62.256 62.256 62.256 0.155 0.145 0.160 0.161 3 z 62.256 62.256 62.256 62.256 0.101 0.097 0.107 0.116 由表 1和表 2可看出:相对于水隔层, 不同厚度空气层条件下各测点处主振频率皆有所降低; 部分情况下出现主振频率值不变的现象, 可能由于爆炸当量较小及测振仪灵敏度较小等原因, 因此需要进一步对信号进行频谱分析。图 3、图 4分别为不同厚度空气层、水隔层结构对各测点振动速度的影响。
从图 3、图 4可看出, 对比相同厚度的空气层和水隔层结构, 空气层结构对水下爆炸振动的衰减更明显。相同厚度水隔层结构条件下:3 m处的振动速度衰减程度最大, 振动速度值最小; 2 m处速度值次之; 1 m处速度值最大。而相同厚度空气层结构条件下:1 m处的振动速度衰减程度最大, 振动速度值最小; 3 m处速度值次之; 2 m处速度值最大。可见空气层结构对近区1 m处的振动速度影响较大。
3. 振动信号的HHT分析
目前国内使用的振动信号采集仪主要采集最大振幅所对应的频率, 即主振频率, 无法对振动信号整个过程的振动频率进行监测; 但爆破工程中需要保护的对象有其自身的固有频率, 为了避免共振造成的破坏, 需要测试振动信号的全局频率, 以观测危险频段的振幅值是否需要进一步防护。本试验爆炸当量较小, 振动信号相对较弱, 同时可能由于测振仪灵敏度较小, 出现部分情况主振频率值不变的现象。因此, 有必要对振动信号进行Hilbert谱分析。
3.1 EMD分解
通过Matlab软件编写相应程序, 将振动测试仪采集到的振动信号读取到矩阵中, 并根据试验设置的采集频率、采集时长等设置运行参数。对事先装载的EMD工具箱, 只需调用命令函数即可使用EMD分解算法。以空气层厚度为40 mm、距离1 m处的垂直方向振动信号为例(见图 5), 利用EMD分解该爆炸振动信号得到各IMF分量, 如图 6所示。
从图 6可看出, 采用EMD分解得到12个IMF分量C1~C12和一个余量r12。其中:C1和C2分量频率最高, 所含能量较低, 代表振动测试过程中周围环境的高频噪声; C3分量代表信号的高频成分; C4~C9分量幅值最大, 所含的能量最多, 是信号的优势频率段, 包含信号的最显著特征; C10~C12为信号的低频成分; r12幅值很小, 可能是信号微弱的趋势或仪器的飘零。
对水下爆炸振动信号经EMD分解后得到的各IMF分量与原始信号进行相关分析, 得到相关分析参数, 如表 3所示。从表 3可看出:各IMF分量的方差值几乎全为零, 且包含信号主要信息的IMF分量段的贡献率达到96.77%。可见, EMD分解过程具有准确性、高效性和自适应性, 保证了信号经分解后仍保留水下爆炸振动信号的非平稳特性, 表明HHT分析方法适用于水下爆炸振动信号处理。
表 3 信号经EMD分解后相关参数Table 3. Related parameters of the signal by EMD decompositionIMF Variance Proportion/% C1 0 1.646 3 C2 0 0.862 9 C3 0 0.626 8 C4 0.000 1 70.307 9 C5 0 20.200 7 C6 0 2.208 9 C7 0 1.254 8 C8 0 1.081 4 C9 0 1.719 6 C10 0 0.015 0 C11 0 0.060 6 C12 0 0.015 0 3.2 振动信号Hilbert谱分析
对垂向爆炸振动信号经EMD分解得到的IMF分量进行Hilbert变换, 得到频率-时间-振幅三维Hilbert谱。以距离1 m、水隔层厚80 mm的情况为例, 对其振动信号的IMF分量进行Hilbert变换得到Hilbert谱, 如图 7所示。
从图 7可看出:幅值较大(即颜色较深)的样点主要分布于7 000时间采样点(即1.4 s)之前, 且800时间采样点(即0.16 s)前振动幅值最大, 相对应的频率为60~65 Hz, 与测试仪采集到的主振频率接近; 频率为5~15 Hz处对应的幅值较大, 且作用时间较长, 该频率段与房屋等建筑物的固有频率较为接近, 因此可能对建筑物造成更大的破坏, 在实际工程爆破中需要加以防护。为此, 采用不同厚度的空气层结构进行防护, 并对采集到的振动信号进行Hilbert谱分析, 如图 8所示。
从图 8可看出, 在有机玻璃-空气层结构下, 最大振幅所对应的主振频率为45~55 Hz, 与测振仪采集到的主振频率接近。相比于水隔层, 有机玻璃-空气层结构下5~15 Hz区间段的幅值减弱较明显, 且作用时间得到缩短, 但随着空气层厚度的增加, 该频带所对应的幅值有所增加。可见, 针对本试验条件, 为了减小水下爆炸引起的地基振动, 并非空气层越厚越好, 当其厚度为120 mm时防护效果较好。
4. 结论
(1) 对于有机玻璃-空气层结构, 随着空气层厚度的增加, 水下爆炸最大振动速度呈先减小后增大趋势, 且当空气层厚度为120 mm时, 振动速度值最小, 对地基振动的衰减效果较好。
(2) HHT法能有效解决振动测试仪对小药量振动频率监测不敏感等问题, 同时可以得到全局频率所对应的幅值, 直观得出所需要防护的振动频率。
(3) 通过HHT法分析得出, 有机玻璃-空气层结构作用下, 能有效衰减5~15 Hz低频区间段的振动速度, 且能够缩短作用时间, 有效防止与建/构筑物之间产生共振。
-
表 1 不同距离、不同厚度空气层对水下爆炸振动的影响
Table 1. Influence of distances and thickness of air interlayer on underwater explosion vibration
Distance/m Direction Main vibration frequency/Hz Maximum vibration velocity/(cm·s-1) ∅40 mm ∅80 mm ∅120 mm ∅160 mm ∅40 mm ∅80 mm ∅120 mm ∅160 mm 1 z 56.763 56.763 56.763 56.763 0.086 0.070 0.054 0.070 2 z 56.763 56.763 56.763 44.556 0.118 0.097 0.082 0.101 3 z 56.763 56.763 56.763 56.763 0.080 0.071 0.061 0.086 表 2 不同距离、不同厚度水隔层对水下爆炸振动的影响
Table 2. Influence of distances and thickness of water interlayer on underwater explosion vibration
Distance/m Direction Main vibration frequency/Hz Maximum vibration velocity/(cm·s-1) ∅40 mm ∅80 mm ∅120 mm ∅160 mm ∅40 mm ∅80 mm ∅120 mm ∅160 mm 1 z 62.256 62.256 62.256 56.763 0.086 0.169 0.184 0.185 2 z 62.256 62.256 62.256 62.256 0.155 0.145 0.160 0.161 3 z 62.256 62.256 62.256 62.256 0.101 0.097 0.107 0.116 表 3 信号经EMD分解后相关参数
Table 3. Related parameters of the signal by EMD decomposition
IMF Variance Proportion/% C1 0 1.646 3 C2 0 0.862 9 C3 0 0.626 8 C4 0.000 1 70.307 9 C5 0 20.200 7 C6 0 2.208 9 C7 0 1.254 8 C8 0 1.081 4 C9 0 1.719 6 C10 0 0.015 0 C11 0 0.060 6 C12 0 0.015 0 -
[1] 库尔·P. 水中爆炸[M]. 罗耀杰, 等译. 北京: 国防工业出版社, 1960.COLE P. Underwater explosion[M]. Translated by LUO Y J, et al. Beijing: National Defense Industry Press, 1960. [2] RAJENDRAN R, LEE J M.Blast loaded plates[J].Marine Structures, 2009, 22(2):99-127. doi: 10.1016/j.marstruc.2008.04.001 [3] 张社荣, 孔源, 王高辉.水下和空中爆炸冲击波传播特性对比分析[J].振动与冲击, 2014, 33(13):148-153. http://www.cnki.com.cn/Article/CJFDTotal-ZDCJ201413027.htmZHANG S R, KONG Y, WANG G H.Comparative analysis on propagation characteristics of shock wave induced by underwater and air explosions[J].Journal of Vibration and Shock, 2014, 33(13):148-153. http://www.cnki.com.cn/Article/CJFDTotal-ZDCJ201413027.htm [4] 高建华, 陆林, 何洋扬.浅水中爆炸及其破坏效应[M].北京:国防工业出版社, 2010.Gao J H, LU L, HE Y Y.Explosion and its damage effect in shallow water[M].Beijing:National Defense Industry Press, 2010. [5] WANG G H, ZHANG S R, YU M, et al.Investigation of the shock wave propagation characteristics and cavitation effects of underwater explosion near boundaries[J].Applied Ocean Research, 2014, 46(2):40-53. https://www.sciencedirect.com/science/article/pii/S014111871400011X [6] PETROV N V, SCHMIDT A A.Multiphase phenomena in underwater explosion[J].Experimental Thermal and Fluid Science, 2015, 60:367-373. doi: 10.1016/j.expthermflusci.2014.05.008 [7] 樊自建, 沈兆武, 马宏昊, 等.空气隔层对水中冲击波衰减效果的实验研究[J].中国科学技术大学学报, 2007, 37(10):1306-1311. doi: 10.3969/j.issn.0253-2778.2007.10.025FAN Z J, SHEN Z W, MA H H, et al.Experimental study on attenuation of underwater shock wave by air interlayer[J].Journal of University of Science and Technology of China, 2007, 37(10):1306-1311. doi: 10.3969/j.issn.0253-2778.2007.10.025 [8] 贾虎, 郑伟花, 罗强, 等.爆炸气泡帷幕对水中冲击波能量的衰减特性[J].含能材料, 2015, 23(10):1015-1019. doi: 10.11943/j.issn.1006-9941.2015.10.018JIA H, ZHENG W H, LUO Q, et al.Attenuation characteristics of underwater explosion bubble curtain on the shock[J].Chinese Journal of Energetic Materials, 2015, 23(10):1015-1019. doi: 10.11943/j.issn.1006-9941.2015.10.018 [9] 贾虎, 沈兆武.空气隔层对水中冲击波的衰减特性[J].爆炸与冲击, 2012, 32(1):61-66. doi: 10.11883/1001-1455(2012)01-0061-06JIA H, SHEN Z W.An investigation into attenuation of underwater shockwave by air interlayer[J].Explosion and Shock Waves, 2012, 32(1):61-66. doi: 10.11883/1001-1455(2012)01-0061-06 [10] 寇晓枫, 王高辉, 卢文波, 等.空气隔层对水下爆炸冲击波的缓冲效应[J].振动与冲击, 2017, 36(3):7-13. http://d.old.wanfangdata.com.cn/Periodical/zdycj201703002KOU X F, WANG G H, LU W B, et al.Mitigation effects of air interlayer on underwater explosion shock wave[J].Journal of Vibration and Shock, 2017, 36(3):7-13. http://d.old.wanfangdata.com.cn/Periodical/zdycj201703002 [11] 姚熊亮, 杨文山, 初文华, 等.水中空气隔层衰减冲击波性能研究[J].高压物理学报, 2011, 25(2):165-172. doi: 10.11858/gywlxb.2011.02.013YAO X L, YANG W S, CHU W H, et al.Research on performance of the underwater air buffer weakening shock wave[J].Chinese Journal of High Pressure Physics, 2011, 25(2):165-172. doi: 10.11858/gywlxb.2011.02.013 [12] 姚熊亮, 刘文韬, 张阿漫, 等.水下爆炸气泡及其对结构毁伤研究综述[J].中国舰船研究, 2016, 11(1):36-45. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgjcyj201601006YAO X L, LIU W T, ZHANG A M, et al.Review of the research on underwater explosion bubbles and the corresponding structural damage[J].Chinese Journal of Ship Research, 2016, 11(1):36-45. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgjcyj201601006 [13] 刘欣, 顾文彬, 陈学平.气泡帷幕对水中冲击波衰减特性的数值模拟研究[J].爆破, 2015, 32(3):79-84. doi: 10.3963/j.issn.1001-487X.2015.03.014LIU 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 [14] 汤有富, 汪泉, 朱恺波, 等.基于HHT变换的小水池水下爆炸振动分析[J].工程爆破, 2017, 23(1):29-33.TANG Y F, WANG Q, ZHU K B, et al.Underwater explosion vibration analysis based on HHT transform in small pond[J].Engineering Blasting, 2017, 23(1):29-33. [15] 张义平. 爆破震动信号的HHT分析与应用研究[D]. 长沙: 中南大学, 2006.ZHANG Y P. HHT analysis and its application of blasting vibration signals[D]. Changsha: Central South University, 2006. [16] HUANG N E, SHEN Z, LONG S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[C]//Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences. The Royal Society, 1998, 454(1971): 903-995. [17] PENG Z K, PETER W T, CHU F L.An improved Hilbert-Huang transform and its application in vibration signal analysis[J].Journal of Sound and Vibration, 2005, 286(1):187-205. https://www.sciencedirect.com/science/article/pii/S0022460X04007928 [18] 张立, 汪大立.水下爆炸炸药能测量消除边界效应的研究[J].爆破器材, 1995(2):1-6.ZHANG L, WANG D L.A study on elimination of boundary effects in underwater explosion testing of explosive[J].Explosive Materials, 1995(2):1-6. -