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

旋转射弹高速倾斜入水多相流场与弹道数值模拟

朱珠 罗松 卢丙举 于勇

引用本文:
Citation:

旋转射弹高速倾斜入水多相流场与弹道数值模拟

    作者简介: 朱 珠(1981- )女,博士,高级工程师,zzqy801@163.com;
    通讯作者: 于勇, yuyong@bit.edu.cn
  • 中图分类号: O359

Numerical simulation of multiphase flow field and trajectory of high-speed oblique water entry of rotating projectile

    Corresponding author: YU Yong, yuyong@bit.edu.cn
  • CLC number: O359

  • 摘要: 基于VOF多相流模型和有限体积法求解水、汽、气多相流动的RANS方程,结合重叠网格技术和six DOF算法对某一型号舰载射弹倾斜入水过程进行数值模拟研究。首先基于该方法研究了射弹旋转效应对射弹运动特性及流体动力特性的影响,然后对不同入水角下倾斜入水过程进行分析,得到不同倾角下旋转射弹入水空泡形态发展规律、弹体运动特征及流体动力特性变化规律。研究结果表明:射弹的旋转有利于弹体在初始对称面内的弹道稳定性,但会降低弹体侧向稳定性,使射弹受到的阻力系数、俯仰力矩系数变小;入水角越小,形成的空泡越不对称,由射弹运动状态的改变引起的空泡形态变化越明显,在超空泡航行阶段,弹体运动较稳定,不同角度下流体动力系数差别很小,当弹体下表面刺破空泡壁沾湿时,弹体运动状态发生较大变化,流体动力系数迅速增大,此时入水角度过小,弹体容易失稳;弹体的沾湿对空泡形态、弹体运动稳定性和流体动力特性有着重要的影响。
  • 图 1  弹体对称面结构示意图

    Figure 1.  Sketch of projectile’s symmetry plane

    图 2  弹体倾斜入水流域对称面示意图

    Figure 2.  Sketch of the fluid domain’s symmetry plane for oblique water entry of projectile

    图 3  30°入水角z=0平面网格示意图

    Figure 3.  Schematic of the grid in z=0 plane of 30 degrees inclined angle

    图 4  平头圆柱入水的位移与速度变化

    Figure 4.  Displacement and velocity of the flat head cylinder water-entry

    图 5  不同旋转速度下弹体质心轨迹图

    Figure 5.  Centroid trajectory of projectile at different rotation velocities

    图 6  不同旋转速度下弹体质心轨迹xOy平面投影曲线

    Figure 6.  Projection of centroid trajectory of projectile on xOy plane at different rotation velocities

    图 7  不同旋转速度下弹体质心速度变化曲线

    Figure 7.  Velocity curves of projectile at different rotation velocities

    图 8  不同旋转速度下弹体姿态角变化曲线

    Figure 8.  Attitudes change of projectile at different rotation velocities

    图 9  不同旋转速度下弹体阻力系数(Cd)变化曲线

    Figure 9.  Drag coefficients (Cd) evolution of projectile at different rotation velocities

    图 10  不同旋转速度下弹体升力(Cl)系数变化曲线

    Figure 10.  Lift coefficients (Cl) evolution of projectile at different rotation velocities

    图 11  不同旋转速度下弹体偏航力矩系数变化曲线

    Figure 11.  Yawing moment coefficients evolution of projectile at different rotation velocities

    图 12  不同旋转速度下弹体俯仰力矩系数变化曲线

    Figure 12.  Pitching moment coefficients evolution of projectile at different rotation velocities

    图 13  弹体不同角度入水空泡形态图

    Figure 13.  Cavitation shapes of the projectile at different inclined water entry angles

    图 14  不同入水角度下弹体质心轨迹图

    Figure 14.  Centroid trajectory of projectile at different inclined water entry angles

    图 15  不同入水角度下弹体质心速度变化曲线

    Figure 15.  The velocity curves of projectile at different inclined water entry angles

    图 16  不同入水角度下弹体滚转角(ϕ)变化曲线

    Figure 16.  Roll angle (ϕ) curves of projectile at different inclined water entry angles

    图 17  不同入水角度下弹体偏航角(θ)变化曲线

    Figure 17.  Yaw angle (θ) curves of projectile at different inclined water entry angles

    图 18  不同入水角度下弹体俯仰角(ψ)变化曲线

    Figure 18.  Pitch angle (ψ) curves of projectile at different inclined water entry angles

    图 19  不同入水角度下弹体阻力系数变化曲线

    Figure 19.  Drag force coefficients evolution of projectile at different inclined water entry angles

    图 20  不同入水角度下弹体升力系数变化曲线

    Figure 20.  Lift force coefficients evolution of projectile at different inclined water entry angles

    图 21  不同入水角度下弹体滚转力矩系数变化曲线

    Figure 21.  Rolling moment coefficients evolution of projectile at different inclined water entry angles

    图 22  不同入水角度下弹体偏航力矩系数变化曲线

    Figure 22.  Yawing moment coefficients evolution of projectile at different inclined water entry angles

    图 23  不同入水角度下弹体俯仰力矩系数变化曲线

    Figure 23.  Pitching moment coefficients evolution of projectile at different inclined water entry angles

  • [1] WORTHINGTON A M, COLE R S. Impact with a liquid surface studied by the aid of instantaneous photography [J]. Philosophical Transactions of the Royal Society of London, 1900, 189: 175–199.
    [2] WORTHINGTON A M, COLE R S. A study of splashes [M]. New York: Longmans Green and Company, 1908: 25−76.
    [3] SAVCHENKO Y. Supercavitation-problems and perspectives [C] // 4th International Sysmposium on Cavitation. Pasadena, USA: California Institute of Technology, 2001: 1−8.
    [4] SAVCHENKO Y N. Control of supercavitation flow and stability of supercavitating motion of bodies [C] // Vki Lecture Series Supercavitating Flows, 2001.
    [5] SEMENENKO V N. Artificial Supercavitation. Physics and calculation [J]. Materials Research, 1(3): 1−5.
    [6] LOGVINOVICH G V. Hydrodynamics of flows with free boundaries [M]. Kiev: Naukova Dumka, 1969.
    [7] 施红辉, 周浩磊, 吴岩, 等. 伴随超空泡产生的高速细长体入水实验研究 [J]. 力学学报, 2012, 44(1): 49–55. DOI: 10.6052/0459-1879-2012-1-lxxb2011-062.
    SHI Honghui, ZHOU Haolei, WU Yan, 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.
    [8] 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.
    [9] 宋武超, 王聪, 魏英杰, 许昊. 回转体倾斜入水空泡及弹道特性实验 [J]. 北京航空航天大学学报, 2016, 42(11): 2386–2394.
    SONG Wuchao, WANG Cong, WEI Yingjie, XU Hao. Experiment of cavity trajectory characteristics of oblique water entry of revolution bodies [J]. Journal of Beijing University of Aeronautics and Astronautics, 2016, 42(11): 2386–2394.
    [10] 马庆鹏. 高速射弹入水过程多相流场特性研究[D]. 哈尔滨: 哈尔滨工业大学, 2014.
    MA Qingpeng. Investigation of multiphase flow characteristics induced by water entry of high-speed projectiles [D]. Harbin: Harbin Institute of Technology, 2014.
    [11] 何春涛. 超空泡射弹结构参数设计与数值模拟研究[D]. 哈尔滨: 哈尔滨工业大学, 2009.
    HE Chuntao. Structure parameter design and numerical simulation study on supercavitation projectile [D]. Harbin: Harbin Institute of Technology, 2009.
    [12] 孙健. 小尺度回转体入水过程的三维数值模拟[D]. 哈尔滨: 哈尔滨工业大学, 2014.
    SUN Jian. Three dimensional simulation in water entry of small scale rotary body [D]. Harbin: Harbin Institute of Technology, 2014.
    [13] 宋武超, 王聪, 魏英杰, 等. 不同头型回转体低速倾斜入水过程流场特性数值模拟 [J]. 北京理工大学学报, 2017, 37(7): 661–666, 671.
    SONG Wuchao, WANG Cong, WEI Yingjie, et al. Numerical simulation of flow field characteristics of low speed oblique water entry of revolution body [J]. Transaction of Beijing Institute of Technology, 2017, 37(7): 661–666, 671.
    [14] 齐亚飞. 弹体高速入水弹道稳定及空泡特性研究[D]. 哈尔滨: 哈尔滨工业大学, 2016.
    QI Yafei. Research on the trajectory stability and cavity characteristics of high-speed water entry projectiles [D]. Harbin: Harbin Institute of Technology, 2016.
    [15] 赵成功, 王聪, 魏英杰, 等. 细长体水下运动空化流场及弹道特性实验 [J]. 爆炸与冲击, 2017, 37(3): 439–446. DOI: 10.11883/1001-1455(2017)03-0439-08.
    ZHAO Chenggong, WANG Cong, WEI Yingjie, 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.
    [16] 赵成功, 王聪, 孙铁志, 等. 初始扰动对射弹尾拍运动及弹道特性影响分析 [J]. 哈尔滨工业大学学报, 2016, 48(10): 71–76. DOI: 10.11918/j.issn.0367-6234.2016.10.010.
    ZHAO Chenggong, WANG Cong, SUN Tiezhi, et al. Analysis of tail-slapping and ballistic characteristics of supercavitating projectiles under different initial disturbances [J]. Journal of Harbin Institute of Technology, 2016, 48(10): 71–76. DOI: 10.11918/j.issn.0367-6234.2016.10.010.
    [17] 李佳川. 高速射弹入水过程流体动力与弹道特性研究[D]. 哈尔滨: 哈尔滨工业大学, 2015.
    LI Jiachuan. Research on water entry hydrodynamic and trajectory characteristics of high-speed projectiles [D]. Harbin: Harbin Institute of Technology, 2016.
    [18] Fluent Inc. Fluent theory guide [M/DK]. 2017.
    [19] MENTER F R. Two-equation eddy-viscosity turbulence models for engineering applications [J]. AIAA Journal, 1994, 32(8): 1598–1605. DOI: 10.2514/3.12149.
    [20] RAYLEIGH L. On the pressure developed in a liquid during the collapse of a spherical cavity [J]. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science: Series 6, 1917, 34: 94–98. DOI: 10.1080/14786440808635681.
    [21] LEE M, LONGORIA R G, WILSON D E. Cavity dynamics in high-speed water entry [J]. Physics of Fluids, 1997, 9(3): 540–550. DOI: 10.1063/1.869472.
    [22] 郭子涛. 弹体入水特性及不同介质中金属靶的抗侵彻性能研究[D]. 哈尔滨: 哈尔滨工业大学, 2012.
    GUO Zitao. Research on characteristics of projectile water entry and ballistic resistance of targets under different mediums [D]. Harbin: Harbin Institute of Technology, 2012.
  • 加载中
图(23)
计量
  • 文章访问数:  34
  • HTML全文浏览量:  95
  • PDF下载量:  12
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-08-27
  • 录用日期:  2018-12-03
  • 网络出版日期:  2019-09-25
  • 刊出日期:  2019-11-01

旋转射弹高速倾斜入水多相流场与弹道数值模拟

    作者简介:朱 珠(1981- )女,博士,高级工程师,zzqy801@163.com
    通讯作者: 于勇, yuyong@bit.edu.cn
  • 1. 中国船舶重工集团公司第七一三研究所河南省水下智能装备重点实验室,河南 郑州 450015
  • 2. 北京理工大学宇航学院,北京 100081

摘要: 基于VOF多相流模型和有限体积法求解水、汽、气多相流动的RANS方程,结合重叠网格技术和six DOF算法对某一型号舰载射弹倾斜入水过程进行数值模拟研究。首先基于该方法研究了射弹旋转效应对射弹运动特性及流体动力特性的影响,然后对不同入水角下倾斜入水过程进行分析,得到不同倾角下旋转射弹入水空泡形态发展规律、弹体运动特征及流体动力特性变化规律。研究结果表明:射弹的旋转有利于弹体在初始对称面内的弹道稳定性,但会降低弹体侧向稳定性,使射弹受到的阻力系数、俯仰力矩系数变小;入水角越小,形成的空泡越不对称,由射弹运动状态的改变引起的空泡形态变化越明显,在超空泡航行阶段,弹体运动较稳定,不同角度下流体动力系数差别很小,当弹体下表面刺破空泡壁沾湿时,弹体运动状态发生较大变化,流体动力系数迅速增大,此时入水角度过小,弹体容易失稳;弹体的沾湿对空泡形态、弹体运动稳定性和流体动力特性有着重要的影响。

English Abstract

  • 高速运动体入水过程中,运动体头部会发生相变产生空化,空化蒸汽与空气排开水质点形成入水空泡。空泡的存在使得弹体与水的接触面积大大减小,可以达到减阻的效果,如何形成包裹全弹身的超空泡是超空泡减阻技术的关键。同时,超空泡形态和变化规律是弹体以正确的姿态和较高的末动能达到目标地点的关键因素。对高速弹体入水空化问题进行研究,具有重要意义。

    入水问题源于Worthington等[1-2]在20世纪初对于球体入水的大量实验,他们在实验中捕捉到了入水溅射和空泡发展。随后,由于水下鱼雷失稳、弹跳等实际问题的出现使得入水问题的研究得到重视,入水问题的研究对象从球体转向了不同头型的圆柱回转体,船舶和实际的一些弹型。俄罗斯和乌克兰学者对于高速射弹的空泡特性研究开展了大量的基础性工作[3-5]。Logvinovich[6]基于独立膨胀原理对入水空泡形态的求解成为之后研究空泡形态和空泡壁发展的理论基础。施红辉[7-8]、宋武超[9]等实验研究了细长体低速入水时入水角度和入水速度等因素对入水空泡形态的影响。马庆鹏[10]对圆柱体垂直入水问题进行了数值模拟研究,详细分析了入水空化过程的流动机理。何春涛[11]详细研究了超空泡射弹的几何参数对空泡形态和阻力特性的影响,为超空泡射弹的设计提供了参考。随后,孙健[12]、宋武超等[13]对小型回转体倾斜和垂直入水问题进行了三维数值计算,利用六自由度算法和动网格对回转体入水弹道进行了求解。齐亚飞[14]利用实验和数值计算对不同头型的回转圆柱射弹高速入水时的弹道和空泡特性进行了研究,得到了不同头型不同初速度下的弹体衰减规律以及不同头型的入水稳定性规律,同时在实验中确定了临界入水角的存在。赵成功等[15-16]研究了细长体倾斜入水的弹道特性和空泡闭合脱落特性,同时研究了弹体的尾拍运动对弹体稳定性的影响。李佳川[17]利用数值计算模拟了高速射弹入水过程中的初始扰动角速度、弹体密度和入水角度对弹道稳定性的影响。

    目前对于入水问题的研究,低速入水问题研究的比较多,高速入水问题由于条件的限制,实验研究和数值计算都很少,特别是弹体高速旋转倾斜入水空化和弹道稳定性问题的研究,有待进一步深入。本文使用基于某一型号舰载弹体改进后的弹型作为研究对象,利用Fluent平台进行计算,考虑自然空化,利用six DOF算法计算弹体姿态轨迹,研究不同入水角度下高速旋转射弹入水过程中的空化特性、弹道特性和流体动力特性的变化规律。

    • 计算采用ANSYS Fluent流体力学计算软件,采用有限体积法求解N-S方程组,多相流场的描述采用均质多相流模型中的VOF模型[18],湍流流场的描述采用SST k-ω 湍流模型[19]。同时,本文选取基于Rayleigh-Plesset气泡方程[20]建立的Schnerr and Sauer空化模型来描述高速射弹入水过程中的自然空化现象。计算过程为瞬态计算,考虑重力的作用。弹体弹道的求解采用刚体运动的动力学方程和运动学方程[17]。在每个时间步,Fluent会计算出弹体受到的流体动力及力矩,根据动力学方程求得弹体的平移加速度和旋转角加速度,之后积分得到弹体的速度和角速度,进一步计算得到弹体新的质心位置和欧拉角,更新网格之后进行下一时间步的求解,直至求解结束。

    • 本文计算采用的射弹模型是如图1所示的实心回转体,弹体模型由头部锥台段、曲线弹身段、圆柱段及后体斜台段组成,弹体总长Lt,头部空化器直径为Dn,弹体最大直径为Dt,弹体长细比为Lt/Dt=8.746。图2展示了计算域对称面内示意图以及边界条件。计算域长20 m,高6.5 m,宽4 m,水深4 m,空气域2.5 m,重力沿y轴负方向,弹体初始质心位置位于(8 m, 0 m, 0 m)处。初始时刻弹体对称轴与水面夹角为α,弹体初速度408 m/s,对应自然空化数为0.001 177,弹体绕着其对称轴高速旋转,旋转速度为7 500 r/min。计算域x+, x-面为滑移壁面,z+, z-面为symmetry边界条件,计算域上方为压力入口,压强与标准大气压一致,底部为压力出口,压强为对应4 m水深的静水压强。

      图  1  弹体对称面结构示意图

      Figure 1.  Sketch of projectile’s symmetry plane

      图  2  弹体倾斜入水流域对称面示意图

      Figure 2.  Sketch of the fluid domain’s symmetry plane for oblique water entry of projectile

    • 计算中采用了Overset网格技术,划分网格时分为部件网格和背景网格。同时导入背景网格和部件网格进行计算,计算时Fluent会对背景网格进行挖洞操作,然后嵌套部件网格,去除Dead Cell后进行计算。流场参数会在部件网格和背景网格重叠边界处进行插值计算。计算网格采用结构化网格划分方式,部件网格区域为一包裹弹体的圆柱,背景网格区域是2.1节中的长方体计算域。图3显示了30°入水角时初始时刻计算采用的进行挖补操作之后的网格在z=0平面上的分布。

      图  3  30°入水角z=0平面网格示意图

      Figure 3.  Schematic of the grid in z=0 plane of 30 degrees inclined angle

      采用Overset网格技术,设置部件网格计算域与运动体同步运动,结合six DOF算法可以实现运动体运动的六自由度运动解算。重叠网格技术计算中,只涉及部件网格和不同区域的背景网格的挖补插值操作,不会产生网格的变形,因此可以保持整个计算过程中都有质量很高的计算网格。

    • 计算过程中,体积分数的求解采用隐式格式(Implicit),相间交界面的构建采用Compressive格式,计算中设置水为主相,水蒸气和空气为第二相,相间相互作用选取Schnerr and Sauer空化模型,入水湍流流场的模拟采用选取SST k-ω湍流模型;时间离散采用一阶精度,空间离散采用二阶精度;瞬时速度场和压力场的耦合计算采用Coupled算法,压力差值采用PRESTO格式,其余项采用二阶迎风格式,计算采用的时间步长为10−6 s。

    • Lee[21]利用能量守恒原理得到了入水空泡的计算方法。忽略弹体入水产生的粘性热,考虑流体不可压,对弹体运用牛顿第二定律可以得到弹体的入水速度方程

      $\frac{{{\rm d}{v_{\rm{p}}}}}{{{\rm d}t}} = g - \frac{1}{{2m}}{\rho _{\rm l}}{A_0}{C_{{\rm{d}}x}}v_{\rm{p}}^2$

      式中:${v_{\rm{p}}}$是弹体入水过程中的瞬时速度,$m$是弹体的质量,${\rho _{\rm l}}$是水的密度,${A_0}$是弹体的横截面积,$g$是重力加速度,${C_{{\rm{d}}x}}$是阻力系数,阻力系数的值采用经验公式${C_{{\rm{d}}x}} = {C_{{\rm{d}}0}} + \sigma $,其中$\sigma $是弹体运动过程中的自然空化数,其数值跟弹体的速度和入水深度h有关,其表达式为$\sigma = \dfrac{{{p_{_0}} + {\rho _{\rm l}}gh - {p_{\rm{v}}}}}{{{{{\rho _{\rm{l}}}v_{\rm{p}}^2}/2}}}$p0为大气压,pv为饱和蒸气压。最终可以得到弹体入水的速度方程的最终形式为

      ${\dot v_{\rm{p}}} = g - \frac{{{\rho _{\rm l}}{A_0}}}{{2m}}\left({C_{{\rm{d}}0}} + \frac{{{p_{_0}} + {\rho _{\rm{l}}}g\int_0^t {{v_{\rm{p}}}{\rm d}t} - {p_{\rm{v}}}}}{{{{{\rho _{\rm{l}}}v_{\rm{p}}^{\rm{2}}}/2}}}\right)v_{\rm{p}}^2$

      为了验证数值计算的可靠性,首先对郭子涛实验研究[22]的平头圆柱垂直入水问题进行数值模拟计算,实验采用的平头圆柱长25.4 mm,直径为12.56 mm,重25.1 g,初速度为603 m/s。计算采用六自由度算法和重叠网格技术。图4显示了圆柱体入水过程中的速度衰减曲线和入水位移曲线,将仿真结果与实验结果及采用以上弹道理论解获得的结果进行了对比,可以看出,数值计算结果与实验结果和理论计算结果均具有较好的一致性,说明数值计算方法是可信的。

      图  4  平头圆柱入水的位移与速度变化

      Figure 4.  Displacement and velocity of the flat head cylinder water-entry

    • 超空泡射弹发射后一般都具有沿自身轴线的角速度,射弹在高速旋转入水及水下运动时,其弹道稳定性和流体动力特性都会有一定的特殊性。为了得到射弹的旋转对超空泡射弹入水过程的影响规律,这里选取图1中的弹体在初速度408 m/s,入水角度α=20°,初始旋转速度ω0分别为0和7 500 r/min的工况进行数值模拟,将得到的结果记性了对比。弹体的空化器直径为Dn=35 mm,质心到头部的距离是弹体总长的60.21%。

      图5显示了射弹在无初始旋转速度和高速旋转状态下的质心在计算流域内的三维轨迹变化曲线。在相同的入水条件下,射弹旋转使得弹体质心的侧向偏移增大,为了进一步的得到弹体在初始对称平面内的弹道偏转情况,图6显示了质心轨迹在xOy平面上的投影,可以看出在高速旋转状态下弹体质心在xOy平面内偏离初始运动方向的程度大幅度减小,这说明弹体的旋转使得弹体在初始对称平面内的弹道稳定性得到了提高。

      图  5  不同旋转速度下弹体质心轨迹图

      Figure 5.  Centroid trajectory of projectile at different rotation velocities

      图  6  不同旋转速度下弹体质心轨迹xOy平面投影曲线

      Figure 6.  Projection of centroid trajectory of projectile on xOy plane at different rotation velocities

      图7显示了不同旋转速度下的弹体质心速度的变化曲线,可以看出,高速旋转状态下的弹体在入水过程中的质心速度衰减明显小于无初始旋转的弹体。弹体高速旋转使得弹体的阻力减小,末动能可以得到很大的提高。图8显示了两种初始旋转状态下的弹体入水过程中的偏航角和俯仰角变化曲线,从图中可以看出,弹体旋转状态下入水俯仰角变化大幅变小,但偏航角变化明显大于无初始旋转速度入水工况。这说明弹体的旋转对弹体水下运动的姿态有着重要的影响,特别是在弹体表面有部分沾湿时,弹体的旋转使得弹体的俯仰姿态的稳定性得到提高,但是偏航姿态的稳定性减弱。这是由于弹体旋转过程中,由于弹体表面的不对称沾湿产生不对称的压力分布,从而产生侧向的作用力,这作用力力会产生较大的偏航力矩。而旋转状态使得弹体的沾湿区域相对于无旋转或近似无旋转状态有所变小,这使得弹体的俯仰力矩减小。这将在流体动力特性曲线中得到证实。

      图  7  不同旋转速度下弹体质心速度变化曲线

      Figure 7.  Velocity curves of projectile at different rotation velocities

      图  8  不同旋转速度下弹体姿态角变化曲线

      Figure 8.  Attitudes change of projectile at different rotation velocities

      图9图10分别显示了两种初始旋转状态下的弹体阻力系数和升力系数的变化规律,从图中可以看出在入水运动初期,两种状态的弹体阻力系数和升力系数基本没有差别,此时弹体处于超空泡航行阶段,弹体仅头部沾湿,其余部分被空化蒸汽包裹;在入水运动后期,高速旋转的弹体阻力系数和升力系数均小于无旋转弹体,此时弹体前弹身下表面部分沾湿,弹体高速旋转改变了弹体沾湿状态时的受力。

      图  9  不同旋转速度下弹体阻力系数(Cd)变化曲线

      Figure 9.  Drag coefficients (Cd) evolution of projectile at different rotation velocities

      图  10  不同旋转速度下弹体升力(Cl)系数变化曲线

      Figure 10.  Lift coefficients (Cl) evolution of projectile at different rotation velocities

      图11图12显示了两种初始旋转状态下的弹体偏航力矩系数和俯仰力矩系数的变化曲线,从图中可以看出,在入水运动前期,弹体处于超空泡航行阶段时,旋转对弹体的偏航力矩系数和俯仰力矩系数影响很小;在入水运动后期,弹体表面沾湿时,高速旋转状态的弹体受到的偏航力矩系数迅速增大,无旋转弹体的偏航力矩系数一直近似为0,高速旋转状态的弹体俯仰力矩系数的量级远小于无旋转状态的弹体。

      图  11  不同旋转速度下弹体偏航力矩系数变化曲线

      Figure 11.  Yawing moment coefficients evolution of projectile at different rotation velocities

      图  12  不同旋转速度下弹体俯仰力矩系数变化曲线

      Figure 12.  Pitching moment coefficients evolution of projectile at different rotation velocities

    • 为了进一步得到入水角度对射弹入水过程影响的规律,对射弹在初速度408 m/s,旋转速度7 500 r/min,入水角度α分别为20°,25°,30°情形下入水过程中的超空泡形态、弹道稳定性、流体动力特性进行了分析。

    • 图13显示了不同入水角度不同时刻空泡形态的变化,空泡壁面取水蒸气相体积分数为0.5的等值面,图中时间从1~19 ms,时间间隔为3 ms。从图中可以看出,弹体迎水面空泡发展快于背水侧,形成的空泡不对称,随着入水深度的增加,这种不对称性逐渐变小,但当弹体运动状态变化使弹身刺破空泡壁时,空泡的对称性被打破,空泡的形态发生较大变化。弹体自然空化剧烈,空泡内充满蒸汽。随着入水深度的增加,自然空化数增大,空泡尺寸有所减小,弹体下表面有所沾湿,弹体姿态开始变化,其中20°入水工况弹体下表面沾湿和姿态变化较明显。同时可以看出,入水角度越小,撞击水面形成的空泡越不对称,开空泡阶段的空泡形态越不对称。

      图  13  弹体不同角度入水空泡形态图

      Figure 13.  Cavitation shapes of the projectile at different inclined water entry angles

      通过分析不同入水角度弹体产生的空泡形态变化过程,我们可以得知,在本文计算工况下,弹体能够较为稳定的射入水中,但随后弹体姿态的轻微变化导致弹体发生“尾拍”现象,不同于正常的“尾拍”情形,由于空泡直径偏小,空泡壁面与弹体下表面距离较近,同时弹体曲线段曲率较大,弹体发生“尾拍”时并非是尾部段接触空泡壁,而是弹体的曲线弹身段先刺破空泡壁沾水,沾湿区域会受到极大的摩擦阻力和和压力,而这一位置位于弹体质心的前方,使得弹体受到的力和力矩是不稳定的,导致弹体轨迹和姿态发生很大变化。

    • 图14显示了不同入水角度下弹体质心在计算流域内的运动轨迹,可以看出在计算的工况下,入水角度越大,弹体质心轨迹偏离初始弹体对称轴方向的程度越小,不同入水角下弹体都会有一定程度的侧向偏移,入水角越大,侧向偏移量越小,其中20°入水工况下弹体能够稳定入水,但在入水深度达到一定程度后,弹体逐渐失稳,质心轨迹发生很大的偏转。同时,图15显示了弹体入水时的速度衰减规律,可以看出在入水运动前期,各角度下弹体的速度衰减较为一致,此时弹体仅头部沾湿,阻力仅来源于头部,随后弹身下表面沾湿,弹体的阻力变大,弹体速度衰减变快,不同角度下弹体表面沾湿大小不一样,速度衰减产生差异,其中20°工况速度衰减明显快于25°和30°工况,这是由于20°工况下的弹体沾湿面积扩张的更快,导致阻力增加的更快。通过这里可以看出,弹体下表面的沾湿对弹体的轨迹和速度会产生很大的影响,弹体侧表面的不对称沾湿是弹体轨迹偏转和速度急剧减小的重要原因。

      图  14  不同入水角度下弹体质心轨迹图

      Figure 14.  Centroid trajectory of projectile at different inclined water entry angles

      图  15  不同入水角度下弹体质心速度变化曲线

      Figure 15.  The velocity curves of projectile at different inclined water entry angles

      图16图17图18显示了弹体入水过程中的欧拉角变化,对于计算的几种工况,滚转角由于初始转动角速度的存在持续增大,受到的扰动较小,不同入水角度下滚转角的变化几乎一样,只在入水后期产生轻微差别;弹体运动初期稳定入水,弹体处于超空泡航行阶段,弹体的俯仰角和偏航角初始阶段变化很小,随后由于弹体姿态的轻微变化导致弹体前弹身下表面沾湿,使得弹体产生较大的力矩,弹体的俯仰角和偏航角随之发生很大的变化。同时可以看到,入水角度越小,弹体俯仰角和偏航角的变化会越大,这种扰动变化会随着入水深度的增加逐渐累积变大。其中对于20°工况弹体在入水初期运动稳定,但当其运动到一定深度后,弹体俯仰角发生极大变化,由正变负,弹体逐渐发生失稳,弹体运动姿态由入水姿态变为出水姿态,弹体抬头发生翻转,同时弹体产生很大的偏航角,弹体的运动姿态相对于初始状态发生很大变化。弹体的侧表面的不对称沾湿对弹体的运动姿态产生了极大的影响,是影响弹体姿态稳定性的重要因素。

      图  16  不同入水角度下弹体滚转角(ϕ)变化曲线

      Figure 16.  Roll angle (ϕ) curves of projectile at different inclined water entry angles

      图  17  不同入水角度下弹体偏航角(θ)变化曲线

      Figure 17.  Yaw angle (θ) curves of projectile at different inclined water entry angles

      图  18  不同入水角度下弹体俯仰角(ψ)变化曲线

      Figure 18.  Pitch angle (ψ) curves of projectile at different inclined water entry angles

    • 定义弹体的流体动力系数如下:

      ${C_i} = \frac{{{F_i}}}{{0.5{\rho _{\rm{l}}}v_0^2S}}$

      ${m_i} = \frac{{{M_i}}}{{0.5{\rho _{\rm{l}}}v_0^2SL}}$

      式中:${F_i}$代表弹体的阻力${F_{\rm{d}}}$和升力${F_{\rm{l}}}$${M_i}$代表弹体的滚转力矩${M_x}$、偏航力矩${M_y}$和俯仰力矩${M_z}$${C_i}$${m_i}$为对应的力和力矩系数,${\rho _{\rm{l}}}$为水的密度,${v_0}$为弹体的初速度,$S = \dfrac{{{\text{π}} D_{\rm{t}}^{\rm{2}}}}{4}$为弹体的最大横截面积,L为弹体的特征长度,这里取弹体的总长度${L_{\rm{t}}}$。其中,弹体的升力和阻力是弹体受力相对于速度坐标系分解得到,弹体的滚转力矩,偏航力矩和俯仰力矩是弹体力矩相对于弹体坐标系的分量。

      图19图20显示了弹体入水过程中阻力系数和升力系数变化情况,从图中可以看出,弹体撞击水面时,阻力系数和升力系数均在极短的时间内达到峰值,随后随着空泡的形成,弹体的阻力系数和升力系数缓慢下降,随后趋于稳定,稳定运动时的升力系数趋于0,阻力系数很小;但随后弹体沾湿,升阻力系数开始持续增大,从这里可以看出,弹体在姿态发生变化时,表面沾湿形成的升力并未成为恢复升力使弹体恢复到空泡内部运动,弹体沾湿形成的受力使弹体的轨迹变化加剧,沾湿面积增大,从而进一步增大了受力。其中对于20°工况弹体失稳后升阻力变化明显大于25°和30°工况。

      图  19  不同入水角度下弹体阻力系数变化曲线

      Figure 19.  Drag force coefficients evolution of projectile at different inclined water entry angles

      图  20  不同入水角度下弹体升力系数变化曲线

      Figure 20.  Lift force coefficients evolution of projectile at different inclined water entry angles

      各力矩系数的变化规律如图21图22图23所示。在弹体撞击水面时刻,偏航力矩系数和俯仰力矩系数受到轻微的扰动,滚转力矩系数没有受到扰动,在弹体带空泡运动初期,弹体的力矩系数趋于平缓且量级很小,在弹体表面沾湿后,弹体的力矩系数产生较大变化和波动,其中25°,30°入水角下弹体力矩系数变化趋势较为一致,20°入水角情形下弹体的力矩系数量级较大,因此20°工况下弹体的运动姿态变化比另外两种工况快。滚转力矩系数持续增大,而偏航力矩系数和俯仰力矩系数正负交替,这是由于旋转导致的弹体坐标系方向的改变造成的。

      图  21  不同入水角度下弹体滚转力矩系数变化曲线

      Figure 21.  Rolling moment coefficients evolution of projectile at different inclined water entry angles

      图  22  不同入水角度下弹体偏航力矩系数变化曲线

      Figure 22.  Yawing moment coefficients evolution of projectile at different inclined water entry angles

      图  23  不同入水角度下弹体俯仰力矩系数变化曲线

      Figure 23.  Pitching moment coefficients evolution of projectile at different inclined water entry angles

    • (1)当弹体处于超空泡航行阶段时,弹体的旋转对弹体的弹道姿态,流体动力特性无明显影响;但入水深度达到一定值时,此时形成的超空泡不能完全包裹弹体,弹体的旋转使得弹体弹道的侧向稳定性下降,但初始对称平面内的弹道稳定性得到大幅度的提高,弹体的旋转使得弹体的升阻力系数减小,弹体速度衰减变慢,弹体的旋转使得弹体的俯仰力矩系数大幅变小,姿态角变化相对于无旋转弹体大幅减小,但旋转使得弹体的偏航力矩系数增大,偏航角变化相对于无旋转弹体加剧。(2)在本文计算的倾斜入水角度下,弹体入水初期产生的空泡不对称性较小,弹体产生覆盖弹身的超空泡,入水角度越小,空泡不对称性越强;入水中后期弹体轨迹和姿态发生变化,弹体下表面刺破空泡壁面,空泡形态产生变化,其中20°工况的弹体沾湿和空泡形态变化最明显。(3)在本文的入水角度下,弹体均能稳定入水,但弹体的运动轨迹会产生不同程度的偏移,同时弹体的姿态也会发生不同程度的变化,弹体入水角度越大,弹体质心的偏移量越小;相较于另外两种入水角度,20°入水时弹体入水运动一段时间后发生严重失稳,弹体质心轨迹和弹体姿态产生的变化极大,速度衰减急剧变快。(4)三种角度入水角度下,弹体在下表面未沾湿前受到的力和力矩系数差别不大,但是当弹体下表面沾湿后,入水角度越小,弹体受到的力和力矩会越大,其中20°入水角下弹体的力和力矩系数的量级远大于另外两种入水角度下的力和力矩中时间后弹体质心轨迹和弹体姿态产生的变化极大,弹体发生反转。

参考文献 (22)

目录

    /

    返回文章
    返回