Numerical simulation on the deflection behavior of large caliber conical nose projectile at oblique high-speed water entry
-
摘要: 结合某大口径锥头弹体高速倾斜入水试验,采用任意拉格朗日-欧拉(arbitrary Lagrange-Euler,ALE)流固耦合方法对弹体倾斜入水偏转行为进行数值模拟,研究了弹体以500 m/s高速倾斜入水过程中不同受力模式、载荷变化特征以及弹体发生偏转的力学机理,分析了入水角度对弹体偏转规律的影响。结果表明:在俯仰力矩作用下,弹体均发生抬头方向偏转,且偏转速度呈现先增大后减小的趋势,偏转程度在不同入水角度范围内呈现不同的变化趋势。当入水角度小于15°时,弹体会发生“跳弹”现象;当入水角度为30°~60°时,弹体偏转趋势基本一致,均由初始倾斜状态逐渐转动至水平状态、竖直状态并最终以弹头入水反方向的“出水”姿态向水下运动;当入水角度为75°时,弹体转动至水平状态后,并未继续偏转至竖直状态,弹头以朝斜上方的姿态向水下运动;弹体的入水侵深随入水角度的增大而增大,且增大趋势近似满足指数函数关系。Abstract: Integrated with a high-velocity oblique water entry test of a large caliber conical nose projectile, the deflection behavior of the projectile obliquely entering water was studied based on the arbitrary Lagrange-Euler (ALE) fluid-structure coupling method. Firstly, based on the experiment of a projectile impacting an inclined water tank at a high velocity, a finite element model was established to simulate the corresponding response characteristics, and the rationality of the numerical model and related method was verified. Secondly, the variation of the contact force mode and load characteristics of the projectile when the water entry velocity was 500 m/s was analyzed, and the corresponding mechanism was discussed. In addition, the influence of water entry angle on the deflection behavior of the projectile was investigated. The related analysis demonstrates that the projectile will deflect upward due to the effect of pitch moment, and the deflection velocity increases first and then decreases gradually during the entry process. The variation trends of deflection degrees are different within different entry angle ranges. When the entry angle is less than 15º, the projectile jumps usually out of the water. When the entry angle is in the range from 30º to 60º, the deflection trend of the projectile is almost the same, i.e., the projectile rotates from the initial tilted state to a horizontal state, then to a vertical state, and moves finally downwards with its nose in the opposite direction to the initial water entry direction. When the entry angle increases to 75º, the projectile cannot continue to rotate to a vertical state after it rotates to a horizontal state, but instead moves downwards in a tilted state with its nose facing upwards. Under different water entry angles, the axial force exerted on the projectile is negative, leading to a continuous decrease in the projectile velocity. Comparatively, the transverse force is positive, and the peak value decreases with increasing water entry angle. Moreover, the penetration depth of the projectile increases with the increase of entry angle, and it almost shows an exponential relationship.
-
作为一种活性金属, 铝粉应用十分广泛, 由于高效的燃烧率, 常被用做推进剂和烟火剂。由于金属氧化而引起的爆炸占了总粉尘爆炸的19%, 其中绝大多数是铝粉爆炸[1]。因此如何安全使用铝粉对铝粉工业中的有效防爆具有重要的意义。
为了在封闭容器中形成实验所必需的粉尘悬浮状态, 有各种扬尘方式[2-6], 包括气流扬尘和散落式扬尘, 无论哪种方式都会在粉尘云中诱导产生湍流(称为扬尘湍流)。具有适当强度的扬尘湍流既是粉尘稳定悬浮必不可少的条件, 又加剧了粉尘等容燃烧过程。扬尘湍流具有瞬态性, 它的初始强度及随时间的衰减特性取决于扬尘装置的设计、高压贮气室的容积和初始压力、燃烧室的容积和几何形状等。湍流速度对粉尘爆炸强度以及粉尘爆炸特性参数的测试有很大影响。Y.K.Pu[7]将甲烷-空气、甲烷-空气-惰性粉尘预混气通过相同的扬尘机构扬至等容燃烧器中燃烧, 获得了气相和两相的湍流燃烧结果, 实验结果表明, 粉尘颗粒的存在对甲烷空气混合物中的湍流燃烧结果影响小于10%。C.S.Tai等[8]进一步研究, 表明有无粉尘对初始扬尘湍流强度影响很小。采用纯空气的扬尘湍流参数评估粉尘空气混合物的扬尘湍流特性, 已成为研究粉尘空气混合物初始扬尘湍流特性对粉尘燃烧爆炸特性影响的一种有效手段, 并广为应用[9-11]。
对于给定的扬尘装置, 粉尘点火时刻的扬尘湍流残存强度与点火延迟时间有关。点火延迟时间(即扬尘电磁阀开启时刻和点火时刻间的延迟时间)常被用来定性表征点火时刻所对应的扬尘湍流残存强度[12-13]。
在点火延迟时间的设置上, 标准20L实验装置一般点火延迟时间统一为60ms; 而对非标准装置, 胡俊等[14]建议在粉尘燃烧的点火延迟时间设置上, 应避免扬尘湍流强度由零达到最大值的这段时间, 以减小由湍流强度在急剧上升过程中的大幅度脉动给点火过程带来的不稳定性; 而Q.M.Liu等[15]直接选取了粉尘爆炸点火延迟时间范围内的中间值。
无论通过哪种方式确定点火延迟时间, 都是在单一的粉尘浓度下进行的, 而且也没考虑粒径的影响。尉存娟等[16]、袁旌杰等[17]已做了关于粉尘粒径的实验研究, 但目前还没有就粉尘浓度对点火延迟时间的影响的研究。另外, 扬尘湍流作为影响粉尘爆炸发展过程的重要因素之一, 在进行粉尘爆炸的安全防护和设计时必须考虑。本文中, 以铝粉为研究对象, 利用实验室自建的5L爆炸罐, 采用高压电极放电研究点火延迟时间, 对不同浓度铝粉的爆炸压力和压力上升速率的影响规律进行研究。
1. 实验
1.1 装置
粉尘爆炸实验装置主要由点火系统、控制系统、采集系统、喷粉系统和5L的圆柱形爆炸容器组成, 如图 1所示。点火系统采用常规电容储能放电, 使用中心点火方式。点火电极为钨电极, 曲率半径0.5mm, 电极间距3mm。控制系统是基于PLC的控制电路, 主要控制电磁阀的开启时间、喷粉时间以及点火延迟时间。喷粉系统由空压机、储气室、电磁阀、单向阀、储粉室和蘑菇形喷头组成, 用来扩散粉尘。5L的爆炸容器高160mm, 内径199mm。
采集系统与文献[18]中的类似, 由安装在爆炸容器上的Kistler压力传感器、电荷放大器、信号处理器、采集卡和瞬态爆炸参数测试软件组成。
1.2 原料
样品为片状铝粉, 平均粒径10~16μm, 厚度1.5~2.0μm。图 2为实验所用片状铝粉的扫描电子显微镜图像。
2. 结果与讨论
在铝粉的爆炸浓度范围内, 选取5组分别为200、262、310、500和1 000g/m3, 通过改变点火延迟时间测定铝粉爆炸的最大爆炸压力。实验用点火电容14μF, 点火电压1 800V, 电容器储能22.68J。吹粉压力0.6MPa, 电磁阀开启时间10ms。实验结果列于图 3中, 为保证数据的准确性, 每次实验至少重复3次, 图中的数据是3次实验的平均。由图可以很清楚地看到, 除262g/m3浓度, 其他浓度的最大爆炸压力和最大压力上升速率都在同一延迟时间达到最大。262g/m3浓度的最大爆炸压力在30ms时最大, 而最大压力上升速率在20ms时达到最大。这种不一致性, 在其他的金属粉尘爆炸中也可以看到, 对于表面氧化过程比中心燃烧过程, 这种现象似乎更普遍[19]。
点火延迟时间在变化过程中存在一个最佳值, 在此点火延迟时间下, 容器内粉尘云散布达到最佳状态, 在电极周围悬浮的粉尘浓度也达到最大值, 在此点火延迟时间点火, 会使爆炸的最大爆炸压力取得最大值。由图 3可知, 铝粉浓度为200、262、310、500和1 000g/m3时, 最佳点火延迟时间分别为20、30、60、60和60ms。
从图 3可以看到, 随着铝粉浓度的增加, 最佳点火延迟时间逐渐增大, 当铝粉浓度大于310g/m3时, 最佳点火延迟时间保持在60ms不变。悬浮在空间的粉尘云是一个不断运动的集合体, 粉尘受重力的影响会发生沉降, 沉降的速度与粒度有一定的关系, 但是粒子之间相互碰撞的布朗运动又阻止它们向下沉降, 即会抵消粒子的沉降。因此对于同一种粉尘, 其沉降速度一定, 当粉尘浓度较低时, 粒子间的布朗运动较弱, 随着点火延迟时间的增加, 粉尘沉降明显。如图 3(a)所示, 铝粉浓度为200g/m3时, 在点火延迟时间20~60ms之间, 最大爆炸压力就从0.546MPa下降到0.355MPa。此后, 随着粉尘浓度的增加, 爆炸容器内的粉尘量增加, 粒子间的布朗运动增强, 粉尘混合均匀需要的时间增加, 因此最佳延迟时间增大。当粉尘量增加到一定程度, 粒子间的布朗运动达到最大, 粉尘悬浮在空间的时间不再变化, 因此粉尘混合均匀的时间也不变, 即最佳点火延迟时间保持不变。这同样可以由实验数据进行验证, 当铝粉浓度为1 000g/m3时, 在点火延迟时间分别为40、50、60和70ms时, 平均最大爆炸压力分别为0.777、0.794、0.820和0.771MPa, 最大爆炸压力随点火延迟时间变化不大。
图 4是不同浓度铝粉在最佳点火延迟时间和固定点火延迟时间为60ms时最大爆炸压力和最大压力上升速率的对比。从图中可以看到, 在500g/m3浓度以下, 铝粉在最佳点火延迟时间下得到的最大爆炸压力和最大压力上升速率明显高于固定点火延迟时间为60ms时得到的。这种差异直观地反映了装置点火延迟时间对粉尘最大爆炸压力和最大压力上升速率的显著影响。因此装置内粉尘云的分布状态, 不完全由脉冲气流气相湍流度决定, 还与粉尘的浓度有关, 浓度决定着粉尘的分散与沉降, 是一个气固两相共同作用的过程, 不能仅凭气相湍流度或单一的粉尘浓度设置固定的装置点火延迟时间。如果忽略了这种影响, 则必然导致测得的粉尘最大爆炸压力偏小, 这不利于粉尘爆炸的风险分析, 也不利于防护装置的设计。因此从安全角度考虑, 对于不同浓度的粉尘不应采用统一的点火延迟时间, 而应找出各自的最佳点火延迟时间。
最大爆炸压力pmax和最大压力上升速率(dp/dt)max随铝粉浓度变化的典型趋势是:pmax起初随着粉尘浓度的增加有较明显的上升。当粉尘浓度高于500g/m3以后, pmax逐渐趋于最大值, 并且在一定范围内相对稳定; (dp/dt)max则随着粉尘浓度的增加而增加。这与O.Dufaud等[20]在20L爆炸罐点火延迟时间为60ms时得到的实验曲线, 在1 000g/m3浓度以前是一致的。本文中, 在最佳点火延迟时间测得的最大爆炸压力0.820MPa、最大压力上升速率71.3MPa/s, 都远低于O.Dufaud等[20]测得的1.1MPa和200MPa/s, 这可能由爆炸容器的尺度效应和点火能的不同引起, O.Dufaud用了2个5kJ点火器, 高强度的点火能会引起pmax和(dp/dt)max的增加。
3. 结论
装置点火延迟时间对铝粉爆炸压力和压力上升速率的影响十分显著, 最佳点火延迟时间下测得的最大爆炸压力和最大压力上升速率明显高于点火延迟时间固定为60ms时的。最佳点火延迟时间不仅与粉尘粒度有关, 而且与浓度有关。粒度相同但浓度不同的铝粉的最佳点火延迟时间也不同, 浓度越低, 最佳点火延迟时间越小, 浓度较大时, 最佳点火延迟时间保持在60ms不变。因此, 仅以气相湍流度和单一浓度的实验值为依据确定固定点火延迟时间, 测得的粉尘爆炸压力和压力上升速率可能明显偏低, 不能真正反应不同浓度粉尘的爆炸威力, 所以可靠测定粉尘对应爆炸参数应先确定其最佳点火延迟时间。
-
材料 ρ/(kg·m–3) E/GPa ν cp/(J·kg–1·K–1) Tr/K Tm/K ˙ε/s–1 A/MPa B/MPa n C G50 7 620 205 0.28 469.0 300 1765 1 1 445 1 326 0.356 0.005 7A04 2 850 69.35 0.33 921.0 293 878 1 602.5 732.1 0.753 0.014 材料 m D1 D2 D3 D4 D5 c0/(m·s–1) S1 γ0 a0 G50 1.120 0.100 0.760 1.57 0 0 4280 1.99 2.00 0.46 7A04 1.015 0.059 0.246 –2.41 –0.1 –0.1 5240 1.40 1.97 0.48 -
[1] VON KARMAN T H. The impact on seaplane floats during landing: NACA technical note No. 321 [R]. Washington: NACA, 1929. [2] LOGVINOVICH G V. Hydrodynamics of flows with free boundaries [M]. Kiev: Naukova Dumka, 1969. [3] MAY A, WOODHULL J C. Drag coefficients of steel spheres entering water vertically [J]. Journal of Applied Physics, 1948, 19(12): 1109–1121. DOI: 10.1063/1.1715027. [4] 陈先富. 弹丸入水空穴的试验研究 [J]. 爆炸与冲击, 1985, 5(4): 70–73.CHEN X F. Experimental studies on the cavitation phenomena as a pellet entering water [J]. Explosion and Shock Waves, 1985, 5(4): 70–73. [5] 张伟, 郭子涛, 肖新科, 等. 弹体高速入水特性实验研究 [J]. 爆炸与冲击, 2011, 31(6): 579–584. DOI: 10.11883/1001-1455(2011)06-0579-06.ZHANG W, GUO Z T, XIAO X K, et al. Experimental investigations on behaviors of projectile high-speed water entry [J]. Explosion and Shock Waves, 2011, 31(6): 579–584. DOI: 10.11883/1001-1455(2011)06-0579-06. [6] 郭子涛, 张伟, 郭钊, 等. 截卵形弹水平入水的速度衰减及空泡扩展特性 [J]. 爆炸与冲击, 2017, 37(4): 727–733. DOI: 10.11883/1001-1455(2017)04-0727-07.GUO Z T, ZHANG W, GUO Z, et al. Characteristics of velocity attenuation and cavity expansion induced by horizontal water-entry of truncated-ogive nosed projectiles [J]. Explosion and Shock Waves, 2017, 37(4): 727–733. DOI: 10.11883/1001-1455(2017)04-0727-07. [7] 刘思华, 王占莹, 李利剑, 等. 头型对射弹高速入水稳定性的影响 [J]. 航空学报, 2023, 44(21): 528437. DOI: 10.7527/S1000-6893.2023.28437.LIU S H, WANG Z Y, LI L J, et al. Influence of nose shapes on high-speed water entry stability of projectile [J]. Acta Aeronautica et Astronautica Sinica, 2023, 44(21): 528437. DOI: 10.7527/S1000-6893.2023.28437. [8] 王云, 袁绪龙, 吕策. 弹体高速入水弯曲弹道实验研究 [J]. 兵工学报, 2014, 35(12): 1998–2002. DOI: 10.3969/j.issn.1000-1093.2014.12.010.WANG Y, YUAN X L, LV C. Experimental research on curved trajectory of high-speed water-entry missile [J]. Acta Armamentarii, 2014, 35(12): 1998–2002. DOI: 10.3969/j.issn.1000-1093.2014.12.010. [9] SHI Y, HUA Y, PAN G. Experimental study on the trajectory of projectile water entry with asymmetric nose shape [J]. Physics of Fluids, 2020, 32(12): 122119. DOI: 10.1063/5.0033906. [10] 马庆鹏, 魏英杰, 王聪, 等. 不同头型运动体高速入水空泡数值模拟 [J]. 哈尔滨工业大学学报, 2014, 46(11): 24–29. DOI: 10.11918/j.issn.0367-6234.2014.11.004.MA Q P, WEI Y J, WANG C, et al. Numerical simulation of high-speed water entry cavity of cylinders [J]. Journal of Harbin Institute of Technology, 2014, 46(11): 24–29. DOI: 10.11918/j.issn.0367-6234.2014.11.004. [11] GAO J G, CHEN Z H, HUANG Z G, et al. Numerical investigations on the oblique water entry of high-speed projectiles [J]. Applied Mathematics and Computation, 2019, 362: 124547. DOI: 10.1016/j.amc.2019.06.061. [12] CHANG Y N, TONG A Y. A numerical study on water entry of cylindrical projectiles [J]. Physics of Fluids, 2021, 33(9): 093304. DOI: 10.1063/5.0059892. [13] 肖海燕, 罗松, 朱珠, 等. 高速射弹小角度入水弹道特性研究 [J]. 北京理工大学学报, 2019, 39(8): 784–791. DOI: 10.15918/j.tbit1001-0645.2019.08.003.XIAO H Y, LUO S, ZHU Z, et al. Trajectory and cavitation characteristics of high-speed projectiles at small angle of water entry [J]. Transactions of Beijing Institute of Technology, 2019, 39(8): 784–791. DOI: 10.15918/j.tbit1001-0645.2019.08.003. [14] 黄振贵, 王瑞琦, 陈志华, 等. 90°锥头弹丸不同速度下垂直入水冲击引起的空泡特性 [J]. 爆炸与冲击, 2018, 38(6): 1189–1199. DOI: 10.11883/bzycj-2018-0115.HUANG Z G, WANG R Q, CHEN Z H, et al. Experimental study of cavity characteristic induced by vertical water entry impact of a projectile with a 90° cone-shaped head at different velocities [J]. Explosion and Shock Waves, 2018, 38(6): 1189–1199. DOI: 10.11883/bzycj-2018-0115. [15] 胡明勇, 张志宏, 刘巨斌, 等. 低亚声速射弹垂直入水的流体与固体耦合数值计算研究 [J]. 兵工学报, 2018, 39(3): 560–568. DOI: 10.3969/j.issn.1000-1093.2018.03.018.HU M Y, ZHANG Z H, LIU J B, et al. Fluid-solid coupling numerical simulation on vertical water entry of projectile at low subsonic speed [J]. Acta Armamentarii, 2018, 39(3): 560–568. DOI: 10.3969/j.issn.1000-1093.2018.03.018. [16] GUO Z T, ZHANG W, XIAO X K, et al. An investigation into horizontal water entry behaviors of projectiles with different nose shapes [J]. International Journal of Impact Engineering, 2012, 49: 43–60. DOI: 10.1016/j.ijimpeng.2012.04.004. [17] SONG Z J, DUAN W Y, XU G D, et al. Experimental and numerical study of the water entry of projectiles at high oblique entry speed [J]. Ocean Engineering, 2020, 211: 107574. DOI: 10.1016/j.oceaneng.2020.107574. [18] 李佳川, 魏英杰, 王聪, 等. 不同扰动角速度高速射弹入水弹道特性 [J]. 哈尔滨工业大学学报, 2017, 49(4): 131–136. DOI: 10.11918/j.issn.0367-6234.201512058.LI J C, WEI Y J, WANG C, et al. Water entry trajectory characteristics of high-speed projectiles with various turbulent angular velocity [J]. Journal of Harbin Institute of Technology, 2017, 49(4): 131–136. DOI: 10.11918/j.issn.0367-6234.201512058. [19] 汪振, 吴茂林, 戴文留. 大口径弹体高速入水载荷特性研究 [J]. 弹道学报, 2020, 32(1): 15–22. DOI: 10.12115/j.issn.1004-499X(2020)01-003.WANG Z, WU M L, DAI W L. Study on load characteristics of high-speed water-entry of large caliber projectile [J]. Journal of Ballistics, 2020, 32(1): 15–22. DOI: 10.12115/j.issn.1004-499X(2020)01-003. [20] 孙玉松, 周穗华, 张晓兵, 等. 基于多介质ALE方法的大型弹体入水载荷特性研究 [J]. 海军工程大学学报, 2019, 31(6): 101–106. DOI: 10.7495/j.issn.1009-3486.2019.06.019.SUN Y S, ZHOU S H, ZHANG X B, et al. On water-impact load of heavy projectiles base on multi-material ALE method [J]. Journal of Naval University of Engineering, 2019, 31(6): 101–106. DOI: 10.7495/j.issn.1009-3486.2019.06.019. [21] 张斌, 李继承, 陈建良, 等. 构型弹体跌落冲击载荷及结构响应特性 [J]. 爆炸与冲击, 2023, 43(3): 033201. DOI: 10.11883/bzycj-2022-0098.ZHANG B, LI J C, CHEN J L, et al. Loading characteristics and structural response of a warhead during drop impact [J]. Explosion and Shock Waves, 2023, 43(3): 033201. DOI: 10.11883/bzycj-2022-0098. [22] 张伟, 肖新科, 魏刚. 7A04铝合金的本构关系和失效模型 [J]. 爆炸与冲击, 2011, 31(1): 81–87. DOI: 10.11883/1001-1455(2011)01-0081-07.ZHANG W, XIAO X K, WEI G. Constitutive relation and fracture model of 7A04 aluminum alloy [J]. Explosion and Shock Waves, 2011, 31(1): 81–87. DOI: 10.11883/1001-1455(2011)01-0081-07. 期刊类型引用(1)
1. 韩向东,王西泉,李超,陈俊彪,李军智,封慧勇,魏雪梅. 多因素下弹体跨介质入水试验系统构建. 兵工学报. 2024(S2): 199-207 . 百度学术
其他类型引用(2)
-