Numerical simulation study on the mechanism and characteristics of high-speed water entry of hollow projectiles
-
摘要: 为分析空心弹高速入水的机理及其特性,基于雷诺时均Navier-Stokes方程、VOF(volume of fluid)多相流模型、Realizable k-ε湍流模型,引入Schnerr-Sauer空化模型和重叠网格技术对空心弹高速入水进行数值模拟研究,获得了通孔孔径和头部形状对空心弹的空化特性、空泡形态和入水运动特性的影响规律。研究显示数值计算的空泡形态和入水速度、位移曲线与实验结果吻合较好,验证了数值模拟方法的可行性。结果表明:当通孔孔径不同时,通孔孔径越大,空化现象越明显,通孔射流越长,但对空泡半径的影响不大;通孔孔径越小,空泡闭合时间越早,与水面碰撞产生的阻力系数峰值越高,空心弹入水稳定后其阻力系数也越大;无量纲直径在0.575~0.600之间时,空心弹的运动最为稳定。当头部锥角不同时,头部锥角越大,空泡直径越大,空化现象出现得越晚,但空化生成的速度更快;随着头部锥角的增大,阻力系数变大,空心弹的速度衰减变快,相同时间运动的距离较短;头部锥角越大,俯仰角的变化越小,空心弹的运动越稳定。Abstract: To analyze the mechanism and characteristics of high-speed water entry of hollow projectiles, a numerical simulation study of the high-speed water entry of the hollow projectiles was carried out based on the Reynolds average Navier-Stokes equation (RANS), the volume of fluid (VOF) multi-phase flow model, the realizable k-ε flux model, the Schnerr and Sauer aeration model, the six-degrees-freedom (6-DOF) motion simulation method, and the overlapping grid technology. The effects of through-hole aperture and head shape on the cavitation characteristics, cavity morphology, and water entry kinematic properties of the hollow projectile were obtained. The effectiveness of the calculation method was verified by comparing it with the water entry experiment. The results show that the numerically calculated cavity morphology and water entry velocity and displacement curves are in good agreement with the experimental results, which verify the effectiveness of the numerical simulation method. When the through-hole aperture is different, the larger the through-hole aperture, the more pronounced the cavitation phenomenon and the longer the through-hole jet, but the effect on the radius of the cavity is not significant. The smaller the through-hole aperture, the earlier the closure time, the higher the peak drag coefficient resulting from impact with the water surface, and the greater the drag coefficient after the hollow projectile has stabilized in the water. The movement of the hollow projectile is most stable when the dimensionless diameter is between 0.575 and 0.600. When the head cone angle is varied, the larger the head cone angle, the larger the diameter of the cavity, and the later the cavitation phenomenon begins, but the faster the cavitation is generated. As the head cone angle increases, the drag coefficient becomes larger and the velocity of the hollow projectile decays faster, moving a shorter distance at the same time. However, the larger the head cone angle, the smaller the change in pitch angle and the more stable the motion of the hollow projectile.
-
Key words:
- hollow projectile /
- cavity evolution /
- high-speed water entry /
- motion characteristics
-
具有一定初速度的弹丸依靠自身的动能撞击并侵入混凝土靶体过程的研究,在民用和军事上都有着广泛的应用,对于工程防护的设计和军事武器的研发具有重要意义[1-3]。关于混凝土侵彻问题的研究已经较为成熟,主要集中于侵彻阻力、侵彻深度、开坑半径、穿透速度等方面的理论和实验研究[4-7]。文鹤鸣[8]进一步总结了相关计算公式和理论的优缺点,指出进行更深入的理论和实验研究的必要性。这些结果在一定程度上很好地解决了相应的工程问题。但是,这些研究中靶体放置于金属筒体内,侵彻过程受到筒体的约束作用,无法考虑实际应力状态对侵彻过程的影响。
弹丸侵彻混凝土靶体,在靶体表面的开坑过程是侵彻作用的起点,所涉及的力学过程很复杂,具有重要的研究意义,引起了极大的关注[9-10]。由于侵彻初期的开坑过程是与惯性作用(密度ρ)、靶体材料的强度特性(屈服强度σy),以及压缩特性(vp/ct,其中vp为弹速,ct为靶材声速)密切相关的,且以靶体材料的强度特性为主,因此,揭示靶体开坑过程的应力路径依赖性显得尤为重要。基于此,本文中将基于先期研制的真三轴静载混凝土Hopkinson实验设备[11-12]进行改进,研制新型的三轴应力状态作用下混凝土侵彻实验装置,对立方体混凝土试件进行不同应力状态下的低速侵彻实验。通过记录不同应力状态下立方体试件侵彻过程中各面上的动态信号,研究应力状态对低速侵彻性能的影响。
1. 三轴应力状态下混凝土侵彻实验系统
图1(a)~(b)为实验原理图。将50 mm× 50 mm× 50 mm的立方体置于6根方杆中间。在侵彻实验之前,采用液压伺服控制系统在立方体试件三个方向分别施加预先设计的应力状态:
σx 、σy 和σz 。y方向右边杆中空,设有6 mm的弹道。高压气体驱动子弹通过该弹道射出,对立方体试件进行侵彻。立方体试件6个面上的动态响应可通过6根方杆上的应变信号进行分析,由此分析混凝土试件的侵彻特性。图1(c)为实验系统照片。实验系统包括侵彻方向(y方向)、垂直于侵彻方向的水平方向x、垂直于侵彻方向的竖直方向z以及伺服液压系统泵站和实验系统控制台。侵彻方向(y方向)包含高压气炮、子弹入射的中空方杆、支撑方杆、y方向液压缸;垂直于侵彻方向的水平方向x包含x方向左支撑方杆、x方向右支撑方杆、x方向液压缸;垂直于侵彻方向的竖直方向z包含z方向下支撑方杆、z方向上支撑方杆、z方向液压缸。实验系统中6根方杆的截面尺寸均为50 mm×50 mm,x方向的2根方杆长度均为2 000 mm;y方向的2根方杆的长度均为1 500 mm;z方向的2根方杆的长度均为1 000 mm。可确保杆上可记录到完整的波形。3个液压缸可独立地给试验系统的三个方向提供最大180 MPa的静载。利用气炮驱动子弹,弹速最大可达300 m/s。实验子弹为平头实心弹,其直径为5 mm、长度为20 mm。弹体材料采用合金结构钢30CrMnSi。y方向右边杆靠近试件部位预设有5个穿透孔,其中2个孔穿透激光进行子弹速度测量,其余孔释放子弹前方的压缩空气。实验结束后,对靶体的破坏情况进行拍摄和测量,记录开坑深度和表面直径等,以评估开坑效果。
子弹开坑过程的端部作用主要包括:端部阻力、侧向扩孔阻力以及侧面摩擦阻力。因为采用平头弹进行侵彻实验,因此端部阻力可根据y方向左边杆上的应力波形来分析,侧向扩孔阻力可根据x方向和z方向的四根杆上的应力信号进行分析,而侧面摩擦力的测试相对较复杂。利用x方向和z方向方杆进行立方体试件侧面阻力测试的计算原理如图2所示。以z方向为例,试件与杆接触面上的剪力(Qzy)的表达式为[13]:
Qzy=dMdz+ρ0Idωdt (1) 式中:M为截面弯矩,I为绕z轴的转动惯量,ω为截面转动的角速度。
如图2(b)所示,在杆前后表面对称贴有应变计zy1和zy2,它们的应变值基本对称,即:|εzy1|≈|εzy2|。则截面的弯矩为:
M=E0L36εzy1 。其中L为杆的边长,E0为杆的弹性模量。因此,可以得到:dMdz=dMdtdtdz=E0L3s6cfdεzy1dt (2) 式中:cf为弯曲波速度,cf ≡ dz/dt。
又由于
ω=dθdt ,以及dθ=ME0Ids=ME0I(1+εzy1)dz (其中ds为杆弯曲后的弧长),可以得到:ω=dθdt=ME0I(1+εzy1)cf (3) 由式(1)~(3),可以得到侧向剪力(Qzy)和剪切应力(
τzy−up )的表达式:Qzy=E0L33cf(1+εzy1)dεdt (4) τzy−up=QzxA=E0L3cf(1+εzy1)dεdt (5) 式中:A为杆的截面积。由此可以基于z杆前后表面的应变信号计算杆与试件接触界面的阻力。
2. 混凝土材料C30的侵彻实验
2.1 实验方案
实验中采用的C30混凝土,由水、水泥、粗骨料、标准砂以一定配合比制得。实验过程主要由三轴应力状态的施加和y方向侵彻两个阶段组成。首先采用液压伺服控制对C30混凝土试件分别同步施加无静载、双轴静载以及三轴静载等三种应力状态的静载值,再利用气炮驱动子弹,在y方向进行真三轴应力状态下的侵彻实验。通过六根杆上的应变片,适时记录子弹开坑过程中的应变波形。
2.2 实验波形分析
真三轴静载下侵彻实验实测到的原始波形如图3所示。其中:[0 MPa,0 MPa,0 MPa]表示三向静载均为0 MPa。子弹侵彻方向即y方向有二个波:子弹端部开坑阻力信号
εpenet 和带弹道的杆上信号εback ,前者用于计算端部阻力,后者没有明确的规律,随应力状态变化,可正可负。x方向的二个波(εx−left 、εx−right )和z方向的二个波(εz−up 、εz−down )在加载历程上具有较好的一致性,用于计算侧向扩孔阻力,其幅值的差异表明开坑过程混凝土的非均匀性。图3(d)所示为按照图2(a)的形式布置的应变片测试得到的信号。在应变片与试件之间预设一段距离,AB段为z方向方杆前后截面上的压缩信号,对其进行平均,与εz−up 信号基本一致;B点后面信号为弯曲波信号,可用于接触界面摩擦阻力的计算,对应的弯曲波速cf≈2 400 m/s。2.3 侵彻阻力
图4为冲击速度50 m/s下试件在静载[0 MPa,0 MPa,0 MPa]状态下,且x、z轴方杆不与试件接触、x轴方杆与试件接触、或x、z轴方杆均与试件接触等三种情况下各杆上应力波形。着弹点在粗骨料上,除了x、z轴方杆不与试件接触情况下的试件有一定坑深之外,其余试件表面未见明显开坑。如图4(a)所示,即便侵彻方向之外的其他两个方向只是与试件接触,都会对侵彻阻力有较大影响。这种接触会增大侧向杆中的压应力信号(图4(b)~(c)),降低试件与杆接触面的摩擦力(图4(d))。随着冲击速度和侧限条件的改变,试件表面开坑,这些应力会产生一些规律性的变化,即相同冲击速度下,侧限应力的增加会降低端部侵彻阻力的幅值及其脉宽,如图5所示。实验得到x、z方向的侧向阻力(图5(b)~(c))和摩擦力(图5(d)~(e))具有较好的规律,即在不同侧限条件下,其侧向阻力及脉宽会发生很大变化,并且侧面摩擦力会随着侧限约束的增加而减小:当x和z轴的静载力不相等时,其对应的阻力也不一致,但有序发展。因此,可以测试材料在侵彻下的各向异性特性,以及应力状态对侵彻性能的影响。
图6为试件在静载[12 MPa, 12 MPa, 0 MPa]下,不同冲击速度对其侧面摩擦阻力的影响。从图中可以看出,随着冲击速度的增大,侧面摩擦阻力是逐渐减小的,并且阻力的时间脉宽也越宽,可见,子弹冲击速度对侧面摩擦阻力影响很大。
如图7所示,试件在一维侵彻下(无三轴静载,[0 MPa, 0 MPa, 0 MPa]),端部阻力幅值较高,且脉宽较宽,表明具有较长的开坑过程。如图7(a)所示,随着静载的应力状态由无静载到y轴静载(12 MPa),侵彻阻力脉宽较窄,幅值有所增强,曲线更加光滑;如图7(b)所示,随着静载的应力状态由无静载到三轴静载(12 MPa),侵彻阻力脉宽较窄、幅值增强、曲线光滑的特点更加明显。这也表明:真三轴静载状态对混凝土侵彻性能的影响机制发生较大的改变,这需要进一步描述。三轴静载下的样品在单轴施加冲击载荷,由于动态泊松效应,必然引起试件在y、z方向的变形。
在无静载、y轴静载、xy轴静载,以及xyz轴静载等四种应力状态下,C30混凝土试样侵彻破坏的形貌如图8所示,其开坑直径d与深度h如图所示。当无静载作用时,试件周围处于自由状态,侵彻作用下混凝土发生较严重的碎裂(见图8(a)),y轴静载情况与之相似(见图8(b)),并且子弹开坑直径随y轴加静载而增大,但开坑深度反而被限制。当试件在x和y方向施加双向静载时,此时试件z方向的两个面处于自由状态,当在x方向冲击加载试件时,试件的破坏呈现明显的层状劈裂现象,层状劈裂的层面与z轴垂直(见图8(c))。这种层状破坏面表现出试件破坏模式对应力状态的依赖性。当在试件的三个方向都施加一定的真三轴静载时,试件在侵彻作用下不会产生较大规模的显著破坏(见图8(d)),其开坑直径和开坑深度都是逐渐减小的。由于试样各个面都受到限制,其失效发生在材料的内部,如图8(d)中的白色条带。这些侵彻破坏模式表明混凝土侵彻性能的应力状态的依赖性,需要进一步深入工作进行解释。
3. 结 论
介绍了一种三维静载下混凝土材料的侵彻实验系统。通过此系统可以研究混凝土、岩石类材料在不同真三轴静载条件下的侵彻性能。实验分析了不同静载条件下C30混凝土三个方向上的端部阻力、侧向阻力和侧向摩擦力。此实验系统为研究侵彻作用下材料的各向异性特性提供了一种有效的实验技术。
感谢研究生陈丽娜、周李姜参与实验工作。
-
表 1 空心弹模型参数
Table 1. Model parameters of hollow projectiles
弹丸模型 d1/mm d2/mm θ/(°) m/g rC/mm S/mm2 M1 1.6 4 60 13.215 14.535 36.474 M2 2 4 60 12.837 14.612 35.343 M3 2.4 4 60 12.412 14.773 33.961 M4 2.8 4 60 11.938 14.952 32.327 M5 2.4 4 120 12.528 14.151 33.961 M6 2.4 4 180 12.576 14.014 33.961 -
[1] ZHAO Q, CHEN Z H, HUANG Z G, et al. Optimization of the aerodynamic configuration of a tubular projectile based on blind kriging [J]. Scientia Iranica, 2019, 26(1): 311–322. DOI: 10.24200/SCI.2017.20015. [2] 黄振贵, 李艳玲, 陈志华, 等. 空心弹的阻力特性与气动外形数值分析 [J]. 兵工学报, 2013, 34(5): 535–540. DOI: 10.3969/j.issn.1000-1093.2013.05.004.HUANG Z G, LI Y L, CHEN Z H, et al. Numerical investigations on the drag and aerodynamic characteristics of a hollow projectile [J]. Acta Armamentarii, 2013, 34(5): 535–540. DOI: 10.3969/j.issn.1000-1093.2013.05.004. [3] WORTHINGTON A M, COLE R S IV. Impact with a liquid surface studied by the aid of instantaneous photography: Paper II [J]. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 1900, 194(252): 175–199. DOI: 10.1098/rsta.1900.0016. [4] WORTHINGTON A M. A study of splashes [M]. New York: Longmans, Green, and Company, 1908: 25−76. [5] DE BACKER G, VANTORRE M, BEELS C, et al. Experimental investigation of water impact on axisymmetric bodies [J]. Applied Ocean Research, 2009, 31(3): 143–156. DOI: 10.1016/j.apor.2009.07.003. [6] TRUSCOTT T T, TECHET A H, BEAL D N. Shallow angle water entry of ballistic projectiles [C]//Proceedings of the 7th International Symposium on Cavitation. Michigan, USA: Ann Arbor, 2009: 1–14. [7] 陈先富. 弹丸入水空穴的试验研究 [J]. 爆炸与冲击, 1985, 5(4): 70–73.CHEN X F. A method of observing both the motion of the free surface and the front of the mass-ejection of shocked lead [J]. Explosion and Shock Waves, 1985, 5(4): 70–73. [8] 施红辉, 周浩磊, 吴岩, 等. 伴随超空泡产生的高速细长体入水实验研究 [J]. 力学学报, 2012, 44(1): 49–55. DOI: 10.6052/0459-1879-2012-1-lxxb2011-062.SHI H H, ZHOU H L, WU Y, et al. Experiments on water entry of high-speed slender body and the resulting supercavitation [J]. Chinese Journal of Theoretical and Applied Mechanics, 2012, 44(1): 49–55. DOI: 10.6052/0459-1879-2012-1-lxxb2011-062. [9] 宋武超, 王聪, 魏英杰, 等. 回转体倾斜入水空泡及弹道特性实验 [J]. 北京航空航天大学学报, 2016, 42(11): 2386–2394. DOI: 10.13700/j.bh.1001-5965.2015.0690.SONG W C, WANG C, WEI Y J, et al. Experiment of cavity and trajectory characteristics of oblique water entry of revolution bodies [J]. Journal of Beijing University of Aeronautics and Astronasutics, 2016, 42(11): 2386–2394. DOI: 10.13700/j.bh.1001-5965.2015.0690. [10] 赵成功, 王聪, 魏英杰, 等. 细长体水下运动空化流场及弹道特性实验 [J]. 爆炸与冲击, 2017, 37(3): 439–446. DOI: 10.11883/1001-1455(2017)03-0439-08.ZHAO C G, WANG C, WEI Y J, et al. Experiment of cavitation and ballistic characteristics of slender body underwater movement [J]. Explosion and Shock Waves, 2017, 37(3): 439–446. DOI: 10.11883/1001-1455(2017)03-0439-08. [11] 罗驭川, 黄振贵, 高建国, 等. 截锥体头型弹丸低速斜入水实验研究 [J]. 爆炸与冲击, 2019, 39(11): 113902. DOI: 10.11883/bzycj-2018-0498.LUO Y C, HUANG Z G, GAO J G, et al. Experiment research of low-speed oblique water-entry of truncated cone-shaped projectile [J]. Explosion and Shock Waves, 2019, 39(11): 113902. DOI: 10.11883/bzycj-2018-0498. [12] 朱珠, 罗松, 卢丙举, 等. 旋转射弹高速倾斜入水多相流场与弹道数值模拟 [J]. 爆炸与冲击, 2019, 39(11): 113901. DOI: 10.11883/bzycj-2018-0315.ZHU Z, LUO S, LU B J, et al. Numerical simulation of multiphase flow field and trajectory of high-speed oblique water entry of rotating projectile [J]. Explosion and Shock Waves, 2019, 39(11): 113901. DOI: 10.11883/bzycj-2018-0315. [13] 邵志宇, 伍思宇, 曹苗苗, 等. 斜截头弹体入水的弹道特性 [J]. 兵工学报, 2022, 43(6): 1255–1265. DOI: 10.12382/bgxb.2021.0301.SHAO Z Y, WU S Y, CAO M M, et al. Water-entry trajectory of truncated cone-shaped projectile [J]. Acta Armamentarii, 2022, 43(6): 1255–1265. DOI: 10.12382/bgxb.2021.0301. [14] 高旭东, 钱建平, 王晓鸣, 等. 空心弹丸流场数值模拟与阻力特性 [J]. 南京理工大学学报(自然科学版), 2005, 29(2): 158–161.GAO X D, QIAN J P, WANG X M, et al. Flowfield calculation and drag characteristic of hollow projectile [J]. Journal of Nanjing University of Science and Technology, 2005, 29(2): 158–161. [15] 任登凤, 谭俊杰, 张军. 非结构隐式方法在空心弹丸流场模拟中的应用 [J]. 力学与实践, 2006, 28(5): 24–27. DOI: 10.3969/j.issn.1000-0879.2006.05.005.REN D F, TAN J J, ZHANG J. Flowfield calculation of hollow projectile using implicit method based on unstructured meshes [J]. Mechanics in Engineering, 2006, 28(5): 24–27. DOI: 10.3969/j.issn.1000-0879.2006.05.005. [16] 杜宏宝, 蒋锋, 黄振贵, 等. 内锥型空心弹阻塞临界入口锥角仿真研究 [J]. 南京理工大学学报(自然科学版), 2018, 42(6): 642–646, 670. DOI: 10.14177/j.cnki.32-1397n.2018.42.06.002.DU H B, JIANG F HUANG Z G, et al. Simulation of critical inlet angle of inner conical hollow projectile under choke flow [J]. Journal of Nanjing University of Science and Technology, 2018, 42(6): 642–646, 670. DOI: 10.14177/j.cnki.32-1397n.2018.42.06.002. [17] 全鑫, 张红艳, 马铁华, 等. 内锥型空心弹阻塞喉径面积比的仿真研究 [J]. 兵器装备工程学报, 2021, 42(4): 97–101. DOI: 10.11809/bqzbgcxb2021.04.018.QUAN X, ZHANG H Y, MA T H, et al. Simulation research on critical area ratio of throat to entrance of inner conical hollow projectile under choke flow [J]. Journal of Ordnance Equipment Engineering, 2021, 42(4): 97–101. DOI: 10.11809/bqzbgcxb2021.04.018. [18] WESSAM M E, HUANG Z G, CHEN Z H. Aerodynamic characteristics and flow field investigations of an optimal hollow projectile [C]// Proceedings of the 5th International Conference on Mechanical Engineering and Mechanics. Yangzhou: ICMEM, 2014: 181–186. DOI: 10.13140/2.1.2037.9529. [19] SAVCHENKO Y N. Hydrodynamic characteristics of a disc with central duct in a supercavitation flow [M]//NESTERUK I. Supercavitation. Berlin: Springer, 2012: 107–113. DOI: 10.1007/978-3-642-23656-3_6. [20] HOU Y, HUANG Z G, CHEN Z H, et al. Different closure patterns of the hollow cylinder cavities with various water-entry velocities [J]. Ocean Engineering, 2021, 221: 108526. DOI: 10.1016/j.oceaneng.2020.108526. [21] HOU Y, HUANG Z G, CHEN Z H, et al. Experimental investigations on the oblique water entry of hollow cylinders [J]. Ocean Engineering, 2022, 266: 112800. DOI: 10.1016/j.oceaneng.2022.112800. [22] LIU H, ZHOU B, YU J W, et al. Experimental investigation on the multiphase flow characteristics of oblique water entry of the hollow cylinders [J]. Ocean Engineering, 2023, 272: 113902. DOI: 10.1016/j.oceaneng.2023.113902. [23] SCHNERR G H, SAUER J. Physical and numerical modeling of unsteady cavitation dynamics [C]// 4th International Conference on Multiphase Flow. New Orleans, LO, USA: ICMF, 2001. [24] CHEN T, HUANG W, ZHANG W, et al. Experimental investigation on trajectory stability of high-speed water entry projectiles [J]. Ocean Engineering, 2019, 175: 16–24. DOI: 10.1016/j.oceaneng.2019.02.021. 期刊类型引用(2)
1. 王聪,许海雨,马贵辉,孙龙泉. 跨介质航行器动力学前沿技术研究进展. 水下无人系统学报. 2024(03): 384-395 . 百度学术
2. 黄振贵,王浩,蔡晓伟,刘想炎,陈志华,秦健,郝戌龙. 速度对射弹垂直破冰入水空泡流动特性的影响. 兵工学报. 2024(10): 3371-3384 . 百度学术
其他类型引用(1)
-