• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

冰孔约束条件下的弹丸倾斜入水实验研究

张东晓 鹿麟 闫雪璞 高词松 胡彦晓 陈凯敏

胡宏伟, 宋浦, 王建灵, 郭炜, 徐洪涛, 金鹏刚, 任松涛. 炸药水中爆炸冲击因子的新型计算方法[J]. 爆炸与冲击, 2014, 34(1): 11-16.
引用本文: 张东晓, 鹿麟, 闫雪璞, 高词松, 胡彦晓, 陈凯敏. 冰孔约束条件下的弹丸倾斜入水实验研究[J]. 爆炸与冲击, 2023, 43(10): 103304. doi: 10.11883/bzycj-2023-0143
Hu Hong-wei, Song Pu, Wang Jian-ling, Guo Wei, Xu Hong-tao, Jin Peng-gang, Ren Song-tao. A new calculation method for shock factor of underwater explosion[J]. Explosion And Shock Waves, 2014, 34(1): 11-16.
Citation: ZHANG Dongxiao, LU Lin, YAN Xuepu, GAO Cisong, HU Yanxiao, CHEN Kaimin. Experimental study on oblique water-entry of projectile constrained by ice hole[J]. Explosion And Shock Waves, 2023, 43(10): 103304. doi: 10.11883/bzycj-2023-0143

冰孔约束条件下的弹丸倾斜入水实验研究

doi: 10.11883/bzycj-2023-0143
基金项目: 国家自然科学基金(52201385);山西省自然科学基金(20210302123023) ;山西省回国留学人员科研资助项目(2020-106)
详细信息
    作者简介:

    张东晓(2000- ),男,硕士研究生,zhangdongxiao_nuc@163.com

    通讯作者:

    鹿 麟(1988- ),男,博士,副教授, lulin2016@nuc.edu.com

  • 中图分类号: O385;TJ012.3

Experimental study on oblique water-entry of projectile constrained by ice hole

  • 摘要: 基于高速摄影技术,开展了冰孔约束条件下的弹丸倾斜入水实验;通过对比分析无冰环境与冰孔约束条件下的弹丸入水运动过程,并将入水运动过程分为空泡扩张、空泡闭合以及空泡溃灭三个阶段进行了研究,得到了冰孔约束条件下弹丸入水的空泡演化特性;通过对比同一直径冰孔约束条件不同入水初速下弹丸的空泡演化过程及速度变化规律,总结得出了入水初速对于冰孔约束条件下弹丸空泡演化特性以及入水运动特性的影响规律。研究结果表明:空泡扩张阶段,冰孔约束条件下产生的喷溅较为分散,弹丸背水面产生部分隆起;此外,冰孔约束条件下的空泡扩张受到阻碍,空泡最大直径减小。空泡闭合阶段,冰孔约束条件下的空泡闭合时间提前,并且撞击冰板的反射流冲击空泡侧壁使空泡发生局部冲击溃灭。空泡溃灭阶段,冰孔约束条件下的溃灭尾迹由局部冲击溃灭、脱落溃灭和正常溃灭组成;空泡溃灭产生的尾迹旋涡较小。随着入水初速的提高,空泡的长度和最大直径明显增大,局部冲击溃灭的宽度增加;冰孔约束条件会使得弹丸在空泡扩张阶段的速度衰减幅度增大,空泡的闭合时间提前,开始溃灭的时刻延后。
  • 冲击波作为爆破战斗部的主要毁伤方式之一,是评估武器毁伤威力的一项重要指标[1]。毁伤威力评估时,需要对战斗部炸点周围的冲击波场进行分布式测量,同时要求多测点实现同步触发。实战环境下战斗部落点及起爆时刻具有一定的随机性,且战斗部爆炸前具有一定的速度,战斗部的运动速度会改变爆炸冲击波的压力场分布。传统的触发方法如断线触发[2]、光触发[3]、无线触发[1]等均难以实现战斗部实战环境下冲击波超压的可靠触发,因此对战斗部动爆压力场的特性分析主要是通过仿真计算,并结合少量试验数据结果和爆炸相似律获得经验公式[4-6],缺少实战环境下的试验研究。

    本文中,提出了一种基于地震波可靠触发的战斗部空中爆炸冲击波超压测试方法,对着靶速度为0、535和980 m/s的战斗部空中爆炸冲击波进行了测试分析。结果表明,基于地震波信号触发测试方法能可靠获取战斗部动爆冲击波超压峰值。试验成果可为实战复杂环境下基于实测数据研究动爆冲击波特性提供依据。

    为了验证基于地震波实现冲击波超压测试触发的可行性,设计了基于加速度信号触发的冲击波超压测试系统并进行了试验。测试系统主要包括传感器、信号调理电路、电源管理、无线通信和信号采集存储等五部分,测试系统组成如图1所示。进行冲击波超压测试时,信号调理电路对传感器获得的冲击波超压信号进行滤波、降噪,在进入FPGA (field programmable gate array)之前进行信号有无和是否达到触发阈值的判断,A/D控制模块将模拟信号转为数字信号后存储在外部同步动态存储器(synchronous dynamic random-access memory, SDRAM)中,最后通过USB (universal serial bus)或无线通信模块将数据上传到上位机上,在上位机上完成数据最终的显示、分析和处理。

    图  1  测试系统
    Figure  1.  Test system

    基于地震波信号的冲击波超压测试触发方法借助配置在各冲击波超压测试节点上的加速度计与信号调理电路,利用地震波传播速度比冲击波传播速度快的特点,在冲击波到达各测试节点之前,提前感应到的加速度信号,经专用调理电路处理后,触发该测试节点的冲击波超压信号存储。其触发原理如图2所示,冲击波超压信号采集缓存区划分为循环采集和时序采集两部分,通过加速度信号启动冲击波超压信号的循环采集完成第1步触发,通过预设超压阈值的比较进行第2步触发,满足触发条件时,立即固化循环采集区,并开始冲击波超压的时序采集。

    图  2  地震波触发冲击波超压测试原理
    Figure  2.  The principle for measuring shock wave overpressure triggered by seismic wave

    爆炸时,形成以炸点为中心向四周传播的地震波,地震波传播速度最快的是纵波,其使地面发生上下振动,在地壳中的传播速度:u0=5.5~7 km/s。

    (1)根据金尼-格雷厄姆公式,计算冲击波超压峰值(pp,MPa):

    pppair=808[1+(fdR4.5)2]1+(fdR0.048)21+(fdR0.032)21+(fdR1.35)2,fd=3pairp0T0Tair
    (1)

    式中:R为比例距离,是观测点到距爆炸中心的距离r(m)与炸药TNT当量W(kg)的立方根之比,即R=r/3W,本文中0.053m/kg1/3R500m/kg1/3pair为试验现场大气压;p0为标准大气压,p0=101.325kPaTair为试验现场大气温度;T0为标准大气温度,T0=288.16K

    (2)根据Rankine-Hugoiot方程,爆炸点空气冲击波传播速度u与冲击波超压峰值pp之间的关系为:

    u=(ppp0×γ+12γ)+1c0
    (2)

    式中:p0为大气压,γ为空气的比热比,c0为波阵面前的空气声速。

    (3)冲击波和地震波传播到观测点的平均时间差Δt为:

    Δt=r(1u1u0)
    (3)

    联立式(1)~(3)可得不同观测距离处冲击波传播平均速度及其与地震波传播平均时间差的关系,如图3所示。由图3(a)可以看出,冲击波超压值随爆距的增大而减小,当距离r5m时,冲击波的传播平均速度u995m/s,该值远小于地震波的传播速度;由图3(b)可以看出,当距离r5m时,冲击波与地震波传播到观测点的平均时间差Δt4.1ms,而地震波信号触发该测点的冲击波信号只需要几十微秒,在该时间差内能够完成触发该测试节点的冲击波信号的存储。因此,当等效TNT装药量不大于100 kg、爆心距不小于5 m时,可以通过安装在测试终端的加速度计采集地震波信号,作为冲击波超压测试的可靠触发信号。

    图  3  冲击波超压峰值和传播时间差随距离的变化
    Figure  3.  Variation of shock wave overpressure peak and propagation time difference with distance

    图4是弹丸爆炸时在距爆心5 m处获取的加速度和冲击波超压信号,可以看出加速度信号的触发阈值先于冲击波信号到达测试节点,提前时长为12.42 ms,这进一步验证了基于地震波信号的冲击波测试触发方法原理的可行性。

    图  4  试验测试的加速度曲线和超压峰值曲线
    Figure  4.  Test curves of acceleration and shock wave overpressure peak

    为分析不同速度战斗部爆炸时冲击波场的分布规律,建立分布式动爆试验测试系统,采用球形裸装药,装药量为1.2 kg。比例距离R分别为4.71、9.41和14.12 m/kg1/3,战斗部着靶速度v0分别为0、535和980 m/s,火炮火药发射获得着靶速度,距目标炸点50 m处顺序放置4台天幕靶,每台间隔20 m,通过区截法获取弹丸在天幕靶处的飞行速度,再结合炮口靶获取的出炮口速度以及制式弹的外弹道模型,计算出目标炸点处弹丸速度。采用模块装药,速度为535和980 m/s弹丸的装药分别为3×B模块和6×B模块,每个B模块装药2.35 kg,引信采用瞬触发引信。战斗部与地靶平面的水平夹角β为45°,测点分布如图5所示,以爆炸中心为原点,建立三维坐标系Oxyz,分别在地靶平面距爆心5、10、15 m处安装冲击波超压测试装置。在Oxy平面上(即地靶平面),爆心到测试点的连线与x轴的夹角为θ,受试验条件限制,共24个测点。理论上测点越多,对爆炸冲击波的重建越有利。

    图  5  战斗部爆炸测点分布
    Figure  5.  Measuring points of blasting warhead explosion

    图6为比例距离R=4.71 m/kg1/3时不同方向上测得的冲击波超压时域曲线。由图6可知,以静爆条件下(v0=0 m/s)的冲击波超压时域曲线为参照,当θ=0°时,动爆冲击波超压远高于静爆冲击波超压;当θ增大到45°时,动爆冲击波超压有所下降,但是仍然高于静爆冲击波超压;当θ增大到90°时,动爆冲击波超压继续降低,v0对冲击波超压的影响变小,不同速度战斗部爆炸的冲击波超压趋于一致,与静爆冲击波压力相当;随着θ的进一步增大,动爆冲击波压力进一步降低,开始低于静爆冲击波压力;当θ增大到180°时,测点处于与战斗部速度完全相反的方向,为压力最低点,此时的冲击波超压远低于静爆冲击波超压。此外,在动爆冲击波超压高于静爆冲击波超压的方向上,即θ在0°~90°和270°~360°范围内时,v0越大,压力越高,冲击波到达时间越短;相反地,在动爆冲击波超压低于静爆冲击波超压的方向上,即θ在90°~270°范围内时,v0越大,压力越低,冲击波到达时间越长。

    图  6  不同方向的冲击波超压时域曲线
    Figure  6.  Shock wave overpressure-time curves in different directions

    读取冲击波的超压峰值,得到不同速度战斗部的爆炸冲击波超压峰值对比曲线,如图7所示。由图7可知:(1)比例距离相同时,战斗部爆炸的冲击波超压峰值随θ的增大近似呈余弦衰减,当θ=0°时,超压峰值最大,θ=180°时,超压峰值最小,并且战斗部着靶速度v0越大,超压峰值衰减得越快;(2)以静爆冲击波超压峰值为参照,在与战斗部速度方向相同的区域(0°~90°和270°~360°),动爆冲击波存在较大的压力升,超压峰值大于静爆状态下的超压峰值,使得冲击波场呈现出局部高压区,而在与战斗部速度方向相反的区域(90°~270°),则存在较大的压力降。由此,可将运动战斗部的爆炸冲击波场分为压力升和压力降两个区域,分界点大约在θ=90°处。

    图  7  不同速度战斗部的冲击波超压峰值曲线
    Figure  7.  Shock wave overpressure peak curves of the blasting warheads with different velocities

    根据战斗部动爆冲击波超压峰值pp,d与战斗部相对靶平面的速度v0及静爆冲击波超压峰值pp,s的关系[4]

    pp,dpp,s=(1+0.31+Rv0c0cosθ)2
    (4)

    战斗部在靶平面的投影速度v0为:

    v0=v0cosβ
    (5)

    由式(1)、(4)、(5)联合可得战斗部动爆冲击波超压峰值的经验公式计算值,实测战斗部动爆冲击波超压峰值pp,e与经验公式计算结果pp,d对比如表1所示,其中ε=(pp,epp,d)/pp,e。由表1可知,实测战斗部静爆冲击波超压峰值与经验公式计算结果一致性较好,除个别点外,实测战斗部动爆冲击波超压峰值与理论值较为接近,相对误差小于20%,且相对误差随着战斗部速度的增大而增大。此外,越靠近战斗部运动方向轴线(θ=0°和θ=180°)的实测冲击波超压峰值相对误差越大,垂直于战斗部运动方向轴线(θ=90°和θ=270°)的实测冲击波超压峰值相对误差较小。

    表  1  试验结果与理论值对比
    Table  1.  Comparison between experimental and theoretical results
    R/(m·kg−1/3θ/(°)v0=0 m/sv0=535 m/sv0=980 m/s
    pp,e/kPapp,d/kPaε/%pp,e/kPapp,d/kPaε/%pp,e/kPapp,d/kPaε/%
    4.710321.39317.351.26434.86355.5718.23483.16389.0219.48
    45343.49317.357.61384.34344.1510.46392.65367.286.46
    90342.92317.357.45285.10317.35−11.31288.32317.35−10.07
    135335.52317.355.41270.90291.64−7.65258.87271.08−4.72
    180320.94317.351.12228.45281.30−23.14187.89252.98−34.64
    225294.20317.35−7.87261.93291.64−11.34215.23271.08−25.95
    270296.57317.35−7.01285.98317.35−10.97294.97317.35−7.59
    315319.73317.350.74368.50344.156.61402.27367.288.70
    9.410112.04104.866.41127.36111.7012.30149.08117.5421.15
    45110.61104.865.19110.62109.670.86118.86113.764.30
    90100.10104.86−4.7696.52104.86−8.6596.21104.86−9.00
    135103.81104.86−1.0288.96100.16−12.5986.6996.34−11.13
    180106.88104.861.8990.3598.25−8.7575.8892.91−22.44
    225102.12104.86−2.6994.81100.16−5.6488.0696.34−9.40
    27098.28104.86−6.7099.25104.86−5.6698.10104.86−6.90
    315111.00104.865.53118.45109.677.41123.47113.767.86
    14.12063.8763.041.3174.8765.8512.0578.2368.2412.77
    4563.2363.040.3066.3565.022.0169.5666.694.12
    9057.7263.04−9.2159.1663.04−6.5559.0363.04−6.78
    13559.0863.04−6.7057.7561.08−5.7850.6359.48−17.48
    18060.5163.04−4.1751.7260.28−16.5645.7358.04−26.91
    22557.7563.04−9.1653.6961.08−13.7748.0659.48−23.77
    27057.7863.04−9.1059.5263.04−5.9257.5263.04−9.58
    31564.7363.042.6169.4265.026.3473.9866.699.85
    下载: 导出CSV 
    | 显示表格

    表1中的冲击波超压峰值为插值点,利用MATLAB数据处理软件的薄板样条插值方法(thin-plate-spline interpolation)对实测数据进行插值(该插值方法可以使得三维超压曲面弯曲能量最小),得到战斗部速度分别为0、535、980 m/s的爆炸冲击波超压峰值场分布和等压曲线,如图810所示。由图8可知,战斗部静爆冲击波超压峰值在各个方向基本相同;由图910可知,战斗部动爆冲击波超压峰值在战斗部运动速度方向增强,在战斗部运动相反方向减弱,且战斗部速度越大,增强和减弱的程度越大。

    图  8  战斗部速度为0 m/s的冲击波超压峰值场分布
    Figure  8.  Shock wave overpressure field for the blasting warhead with the velocity of 0 m/s
    图  9  战斗部速度为535 m/s的冲击波超压峰值场分布
    Figure  9.  Shock wave overpressure field for the blasting warhead with the velocity of 535 m/s
    图  10  战斗部速度为980 m/s的冲击波超压峰值场分布
    Figure  10.  Shock wave overpressure field for the blasting warhead with the velocity of 980 m/s

    提出了一种基于地震波可靠触发的战斗部空中爆炸冲击波超压测试方法,并对速度为0、535和980 m/s的战斗部空中爆炸冲击波进行了测试研究,通过测试结果和经验公式计算值的对比分析,以及重建的战斗部动爆冲击波超压三维可视化模型,可以得出以下结论:

    (1)本文中提出的测试方法能可靠获取战斗部动爆冲击波超压峰值;

    (2)战斗部动爆冲击波超压峰值在战斗部运动速度方向增强,在战斗部运动相反方向减弱,且战斗部速度越大,增强和减弱的程度越大。

  • 图  1  实验系统示意图

    Figure  1.  Schematic diagram of experimental system

    图  2  实验现场布置

    Figure  2.  Experimental site

    图  3  实验弹丸模型

    Figure  3.  Experimental projectile model

    图  4  实验工况示意图

    Figure  4.  Schematic diagram of experimental conditions

    图  5  空泡扩张阶段空泡演化照片

    Figure  5.  Photos of cavity evolution in water-entry cavity expansion stage

    图  6  空泡扩张阶段空泡细节图

    Figure  6.  Diagram of cavity detail in cavity expansion stage

    图  7  空泡扩张阶段(3.72 ms)仰视照片

    Figure  7.  Photos of cavity expansion stage (3.72 ms) from bottom view

    图  8  空泡闭合阶段空泡演化照片

    Figure  8.  Photos of cavity evolution in cavity contraction stage

    图  9  空泡闭合阶段仰视照片

    Figure  9.  Photos of cavity closure stage from bottom view

    图  10  空泡溃灭阶段空泡演化照片

    Figure  10.  Photo of cavity evolution in cavity collapse stage

    图  11  不同初速下弹丸的入水过程

    Figure  11.  Projectile water-entry processes at different initial velocities

    图  12  闭合以及溃灭时刻的空泡细节

    Figure  12.  Diagram of cavity detail at cavity closure and collapse moment

    图  13  不同入水初速下的空泡演化轮廓图

    Figure  13.  Cavity evolution contours at different initial velocities

    图  14  不同工况下的弹丸速度变化曲线

    Figure  14.  Projectile velocity attenuation under different working conditions

    图  15  不同工况下弹丸的加速度变化曲线

    Figure  15.  Projectile acceleration curves under different working conditions

    图  16  不同工况下弹丸的姿态角变化曲线

    Figure  16.  Projectile attitude angle curves under different working conditions

  • [1] VON KARMAN T. The impact on seaplane floats during landing: NACA-TN-321 [R]. National Advisory Committee for Aeronatics, 1929: 309-313.
    [2] FOROUZANI H, SARANJAM B, KAMALI R. A study on the motion of high speed supercavitating projectiles [J]. Journal of Applied Fluid Mechanics, 2018, 11(6): 1727–1738. DOI: 10.29252/jafm.11.06.28807.
    [3] ERFANIAN M R, ANBARSOOZ M, RAHIMI N, et al. Numerical and experimental investigation of a three dimensional spherical-nose projectile water entry problem [J]. Ocean Engineering, 2015, 104: 397–404. DOI: 10.1016/j.oceaneng.2015.05.024.
    [4] AKBARI M A, MOHAMMADI J, FEREIDOONI J. Stability of oblique water entry of cylindrical projectiles [J]. Journal of Applied Fluid Mechanics, 2021, 14(1): 301–314. DOI: 10.47176/jafm.14.01.31682.
    [5] AKBARI M A, MOHAMMADI J, FEREIDOONI J. A dynamic study of the high-speed oblique water entry of a stepped cylindrical-cone projectile [J]. Journal of the Brazilian Society of Mechanical Sciences and Engineering, 2021, 43(1): 2. DOI: 10.1007/s40430-020-02727-2.
    [6] SHI H H, ITOH M, TAKAMI T. Optical observation of the supercavitation induced by high-speed water entry [J]. Journal of Fluids Engineering, 2000, 122(4): 806–810. DOI: 10.1115/1.1310575.
    [7] 高英杰, 孙铁志, 张桂勇, 等. 回转体高速倾斜入水的流场特性及结构响应 [J]. 爆炸与冲击, 2020, 40(12): 123301. DOI: 10.11883/bzycj-2020-0014.

    GAO Y J, SUN T Z, ZHANG G Y, et al. Flow characteristics and structure response of high-speed oblique water-entry for a revolution body [J]. Explosion and Shock Waves, 2020, 40(12): 123301. DOI: 10.11883/bzycj-2020-0014.
    [8] LU L, WANG C, LI Q, et al. Numerical investigation of water-entry characteristics of high-speed parallel projectiles [J]. International Journal of Naval Architecture and Ocean Engineering, 2021, 13(5): 450–465. DOI: 10.1016/j.ijnaoe.2021.05.003.
    [9] LU L, GAO C S, QI X B, et al. Numerical study on the water-entry characteristics of asynchronous parallel projectiles at an oblique impact angle [J]. Ocean Engineering, 2023, 271: 113697. DOI: 10.1016/j.oceaneng.2023.113697.
    [10] BERGSMA J M, BOUHUYS C W, SCHAAP T, et al. On the measurement of submersion ice resistance of ships, using artificial ice [C]// Proceedings of the Twenty-fourth International Ocean and Polar Engineering Conference. Busan: ISOPE, 2014.
    [11] ZONG Z, YANG B Y, SUN Z, et al. Experimental study of ship resistance in artificial ice floes [J]. Cold Regions Science and Technology, 2020, 176: 103102. DOI: 10.1016/j.coldregions.2020.103102.
    [12] 张军, 蔡晓伟, 宣建明, 等. 弹体穿越冰水混合物流动过程的数值模拟 [J]. 弹道学报, 2020, 32(3): 35–40. DOI: 10.12115/j.issn.1004-499X(2020)03-008.

    ZHANG J, CAI X W, XUAN J M, et al. Numerical simulation of flow filed of projectile passing through ice water mixture [J]. Journal of Ballistics, 2020, 32(3): 35–40. DOI: 10.12115/j.issn.1004-499X(2020)03-008.
    [13] YOU C, SUN T Z, ZHANG G Y, et al. Numerical study on effect of brash ice on water exit dynamics of ventilated cavitation cylinder [J]. Ocean Engineering, 2022, 245: 110443. DOI: 10.1016/j.oceaneng.2021.110443.
    [14] 张健宇. 航行体冰孔约束出水空泡演化及载荷特性研究 [D]. 大连: 大连理工大学, 2021.

    ZHANG J Y. Study on cavity evolution and load characteristics of water exit of underwater vehicle constrained by ice environment [D]. Dalian: Dalian University of Technology, 2021.
    [15] 蔡晓伟, 宣建明, 王宝寿, 等. 细长体穿越冰-水混合物的出水流场数值模拟 [J]. 兵工学报, 2020, 41(S1): 79–90. DOI: 10.3969/j.issn.1000-1093.2020.S1.012.

    CAI X W, XUAN J M, WANG B S, et al. Numerical simulation of thin body passing through the ice-water mixture flow field [J]. Acta Armamentarii, 2020, 41(S1): 79–90. DOI: 10.3969/j.issn.1000-1093.2020.S1.012.
    [16] WANG H, LUO Y C, CHEN Z H, et al. Influences of ice-water mixture on the vertical water-entry of a cylinder at a low velocity [J]. Ocean Engineering, 2022, 256: 111464. DOI: 10.1016/j.oceaneng.2022.111464.
    [17] LOGVINOVICH G V. Hydrodynamics of free-boundary flows [R]. Jersualem: IPST Press, 1972.
  • 期刊类型引用(0)

    其他类型引用(3)

  • 加载中
图(16)
计量
  • 文章访问数:  227
  • HTML全文浏览量:  54
  • PDF下载量:  68
  • 被引次数: 3
出版历程
  • 收稿日期:  2023-04-19
  • 修回日期:  2023-08-07
  • 刊出日期:  2023-10-27

目录

/

返回文章
返回