初始扰动对射弹尾拍运动及弹道特性影响分析

赵成功，王 聪，孙铁志，张孝石

(哈尔滨工业大学 航天学院，哈尔滨150001)



初始扰动对射弹尾拍运动及弹道特性影响分析

赵成功，王 聪，孙铁志，张孝石

(哈尔滨工业大学 航天学院，哈尔滨150001)

为研究射弹的初始扰动对射弹的尾拍运动、超空泡形态及弹道特性的影响，基于计算流体动力学程序CFX，通过二次开发使用CEL语言将刚体动力学方程嵌入计算模块并实现刚体动力学方程与流体(unsteady reynolds averaged navier-stokes,URANS)方程的耦合求解，建立了超空泡射弹在不同初始扰动下自由减速运动的计算流体动力学模型.在相同的初始速度下，计算并对比分析射弹在不同初始扰动下的空泡形态特征、弹尾压力分布特征及弹道特性.计算结果表明：射弹发生尾拍运动时弹体尾部与空泡壁面发生碰撞，这种碰撞破坏了空泡形态的对称性；初始扰动对射弹的水平速度与位移影响较小，但对垂直位移影响较大；初始扰动越大，射弹的俯仰角、角速度和角加速度的振幅越大且频率越高.

多相流；超空泡射弹；尾拍运动；初始扰动；CFD

水下高速运动航行体表面附近的水域会因压力低于水的饱和蒸汽压而发生空化现象，进而在航行体表面形成空泡，当速度较高时会形成完全包裹航行体表面的超空泡[1].超空泡的形成能够明显降低航行体水下运动时的阻力，可大幅提高航行体的水下运动速度[2]，但是由于空泡流动本身的复杂性使得超空泡技术在实际应用中面临诸多挑战，以超空泡航行体的稳定性和控制问题尤为突出[3-4].

随着计算机技术的飞速发展，基于N-S方程的计算流体动力学(CFD)方法已成为揭示一些复杂流动现象的有效途径.本文通过对流体动力计算软件CFX进行二次开发，运用程序自带的CEL语言将射弹在垂直平面内运动的刚体运动方程离散化并添加到程序的计算求解过程中，结合动网格技术，实现了多相流动非定常雷诺平均N-S方程(URANS)与刚体运动学方程的耦合求解，建立了水下超空泡射弹纵向平面的尾拍运动动力学计算模型并分析了初始扰动对空泡形态及射弹尾拍运动和弹道特性的影响.本文建立的计算模型考虑了射弹与空泡表面复杂的作用力、射弹的自由运动、空泡的时变特征及湍流对射弹的流体动力影响.

1 数值计算模型

1.1 基本控制方程

本文研究水下运动的超空泡射弹，涉及的是自然空化流动问题，流体域的两相介质为水蒸汽和水.连续性方程为

动量方程为

1.2 湍流模型

湍流是一种强非线性流体状态，是空化流动的主要特征之一.空泡发展及闭合等均伴随着湍流现象.本文采用的湍流模型是基于Baselinek-ω模型的SST(shear stress transport)湍流模型.该湍流模型能够较好地解决湍流剪切应力的传输问题，对于逆压梯度下的流动分离问题有更精确的预测，基本方程为[17]

SST湍流模型的涡黏度限制方程为

式中：μ为流体黏性系数；μt为湍流黏度；νt=μt/ρ；Pk为湍动能中由黏性力产生的部分；Pkb为湍动能中由浮力产生的部分；F1、F2分别为混合函数；常数β′=0.09、α1=5/9；S为应变率的不变测度.

1.3 空化模型

用来描述流场空化程度强弱的特征参数为空化数，为

不计热传输和非平衡相变效应，空化流动中水蒸汽的体积分数输运方程为

式中：pv、p分别为气泡内压力和环境压力；αv为蒸汽的体积分数；rnuc=5×10-4为成核点的体积分数；RB为气泡半径；Fe=50，Fc=0.01均为经验常数.

2 模型尺寸及网格划分

本文计算采用的射弹模型为具有大长细比的锥柱结合体，如图 1所示.

射弹结构包括3段，分别为头部柱段、肩部锥段和后体柱段.头部柱段直径4 mm，长度3 mm；肩部锥段的半锥角为4°；后体柱段直径D=10 mm；弹体总长L= 170 mm；质量为0.09 kg；质心距离弹体头部中心点95.8 mm；射弹绕通过质心且垂直于弹体轴线的坐标轴的转动惯量为1.72×10-4kg·m2.

计算流域为圆柱体，入口距弹体头部长度取2L=340 mm，出口距弹体尾部取20L=3 400 mm，直径取60D=600 mm，x轴位于弹体轴线上，计算中考虑重力影响.整个计算域均采用六面体结构化网格离散，单元数量约为163万，弹体附近网格如图 2所示.计算中采用CFX中CEL语言将刚体运动学方程离散后并嵌入程序包中与U-RNAS方程耦合求解，瞬态计算开始前，首先通过调整流场参数得到覆盖整个弹体的超空泡流动稳态数值解作为后续计算的初始值，以利于计算收敛.射弹初始位置位于水面下5 m深度，初始速度v0= 150 m/s为水平方向，初始自然空化数σv=0.013 1.

图2 弹体模型及计算网格

针对图1所示的射弹模型，建立了粗糙(101万)、中等精细(163万)、精细(210万)3种不同精度的网格，3种网格数量下射弹头部附近压力系数对比如图3所示.

图3 不同网格数量下射弹头部附近压力系数对比

从图3中可以看出，3种网格精度下射弹头部附近压力系数变化差异较小，其中精细网格(210万)和中等精细网格(163万)的计算结果基本一致，而粗糙网格(101万)在射弹头部与精细网格和中等精细网格的计算结果有一定差异，综合考虑，为了节省计算资源，并且达到所需精度，本文选用中等精细网格(163万)进行计算.

为验证本文所建立数值计算模型的有效性，进行了v=150 m/s，水深5 m时的仿真计算，空泡形态与Lognovinch独立膨胀原理所得的空泡形态对比如图4所示.

图4 空泡形态计算结果对比

由图4可以看出，本文数值模型所得空泡形态与独立膨胀原理所得空泡形态基本一致，在空泡尾部由于二者的处理方法不同而略有差异，这说明了本文所建立数值计算模型的有效性.

3 结果与分析

3.1 空泡形态分析

图5 射弹尾拍运动时空泡形态变化(T为尾拍振动周期，T=0.025 s)

Fig.5 Variations of cavity profile during the tail-slapping of projectile(Tas tail-slapping vibration period,T=0.025 s)

由图 5可以看出，射弹的尾拍运动引起空泡形态的改变，破坏空泡原有的光滑对称壁面.弹体尾部与空泡壁面碰撞的位置，空泡形态发生显著改变，碰撞点处发生动量交换，空泡壁面凸出，造成空泡形态不对称.

以射弹尾部截面为0点， 图 5 (c)中发生撞击时弹体尾部沿轴向的压力分布如图 6所示.

图6 弹体尾部压力分布

由图 6可以看出，当射弹尾部撞击空泡壁面时，射弹尾部表面的压力由原来空泡内较低且均匀的空化压力改变为图6所示的在弹体尾部0

3.2 尾拍运动特性分析

图7 不同初始扰动下射弹的俯仰角变化

从图 7中可以看出，不同初始扰动下射弹俯仰角的变化不同，具体表现为俯仰角的幅值和变化周期均不相同.初始扰动角速度越大，俯仰角变化的最大幅值越大，变化频率较快.在t=0.055 s后，俯仰角不再呈现周期性变化，因为此时空泡开始溃灭于弹体尾部并逐渐减小至消失，如图 8所示.

不同初始扰动角速度下射弹尾拍运动的俯仰角速度和俯仰角加速度变化历程分别如图 9, 10所示.

图8 溃灭于弹体尾部的空泡形态

Fig.8 Cavity profile of projectile tail

图9 不同初始扰动下射弹的俯仰角速度变化

Fig.9 Angular velocity variations of pitch angle under different initial disturbances

图10 不同初始扰动下射弹的俯仰角加速度变化

Fig.10 Angular acceleration variations of pitch angle under different initial disturbances

尾拍振动的周期越来越短，频率越来越高，随着后期射弹速度的衰减，空化数变化，超空泡长度变短并逐渐溃灭于弹体上，射弹尾拍振动时刺透空泡壁面浸入流体的深度增加，造成尾拍反力的增加，进而影响尾拍运动俯仰角加速度的增加；当空泡开始溃灭于弹体上并逐渐消失时，射弹由尾拍运动变为双空泡运动(如图 8所示)，尾拍反力急剧减小并消失，故射弹的俯仰角加速度的变化也呈现先增加后急剧减小的趋势；且初始扰动角速度越大，射弹尾拍运动的俯仰角速度和角加速度的幅值也越大，变化频率越高.

3.3 弹道特性分析

图 11给出了射弹在不同初始扰动下以相同初速度自由运动的水平速度变化情况.初始扰动越大，射弹的速度衰减越快，但整体而言射弹水平速度变化的差别并不明显.

图11 不同初始扰动下射弹的水平速度变化

Fig.11 Variations of horizontal velocity under different initial disturbances

不同初始扰动下射弹的阻力系数变化如图 12所示.阻力系数的变化主要包括两个比较明显的阶段，当t<0.06 s，阻力系数主要呈现为周期振动特性，这是由于射弹尾拍运动时弹体尾部刺透空泡壁面并浸入水中增加了流体阻力；当t>0.06 s，阻力系数开始急剧增加，主要是由于空泡长度变短接触到弹体尾部并开始逐渐溃灭于弹体上，因此阻力系数急剧增加.不同初始扰动下尾拍运动时弹体尾部浸入水中的深度不同，所以引起的阻力变化幅值也不同，初始扰动越大，阻力系数变化的振幅越大，频率越高.

图12 不同初始扰动下射弹的阻力系数变化

Fig.12 Variations of drag coefficient under different initial disturbances

图 13, 14分别给出了不同初始扰动下射弹自由衰减运动的水平位移和垂直位移的变化历程.从图中可以看出，不同初始扰动对射弹的水平位移影响较小，而对射弹的垂直位移影响较大.

初始扰动分别为2 rad/s和8 rad/s的两种情况下，在相同运动时间内射弹的垂直位移相差约1倍；随着初始扰动的增加，射弹的垂直位移增大；射弹垂直位移的增加是重力和扰动的共同作用.

图13 不同初始扰动下射弹的水平位移变化

Fig.13 Variations of horizontal displacement under different initial disturbances

图14 不同初始扰动下射弹的垂直位移变化

Fig.14 Variations of vertical displacement under different initial disturbances

4 结 论

1)不同初始扰动下射弹尾拍运动的俯仰角、俯仰角速度、俯仰角加速度、阻力系数都表现出周期性变化的特点，且初始扰动越大，变化的幅值越大，频率越高，尾拍振动的频率随着时间发展越来越高.

2)不同初始扰动对射弹速度衰减和水平位移变化的影响较小，但对射弹垂直位移的影响较大.

3)超空泡射弹尾拍运动时弹体尾部与空泡壁面发生碰撞并浸入流体中，破坏空泡形态的对称性，射弹尾部碰撞点附近产生较大的压力峰值.

[1] 孟庆昌，张志宏，刘巨斌，等. 水下高速航行体超空泡流动研究进展[J]. 船海工程, 2006，35(6): 26-29.DOI: 10.3963/j.issn.1671-7953.2006.06.008.

MENG Qingchang, ZHANG Zhihong, LIU Jubin, et al. Research development of high-speed super-cavitating flow of the high-speed underwater vehicle[J]. Ship & Ocean Engineering, 2006，35(6): 26-29.DOI: 10.3963/j.issn.1671-7953.2006.06.008.

[2] SAVCHENKO Y N. Control of supercavitation flow and statbility of supercavitating motion of bodies [C] // VKI Special Course on Supercavitating Flows 2001.Brussels: Rto2avt and vki, 2001: 313-341.

[3] STINEBRING D R, COOK R B, DZIELSKI J E, et al. High-speed supercavitating vehicles [C] // Proceedings of the AIAA Guidance, Navigation and Control Conference and Exhibit 2006. Keystone, CO: AIAA, 2006: 3092-3099.

[4] KIRSCHNER I N, KRING D C, STOKES A W, et al. Control strategies for supercavitating vehicles [J]. Journal of Vibration and Control, 2002, 8(2): 219-242. DOI: 10.1177/107754602023818.

[5] 孟庆昌，张志宏，顾建农，等. 超空泡射弹尾拍分析与计算[J].爆炸与冲击，2009，29(1)：56-60. DOI: 10.3321/j.issn:1001-1455.2009.01.011.

MENG Qingchang, ZHANG Zhihong, GU Jiannong, et al. Analysis and calculation for tail-slaps of supercavitating projectiles[J]. Explosion and Shock Waves，2009，29(1)：56-60. DOI: 10.3321/j.issn:1001-1455.2009.01.011.

[6] RUZZENE M, SORANNA F. Impact dynamics of elastic stiffened supercavitating underwater vehicles[J]. Journal of Vibration and Control, 2004, 10(2): 243-267. DOI: 10.1177/1077546304035607.

[7] RUZZENE M, KAMADA R, BOTTASSO C L, et al. Trajectory optimization strategies for supercavitating underwater vehicles [J]. Journal of Vibration and Control, 2008, 14(5): 611-644.DOI: 10.1177/1077546307076899.

[8] SAVCHENKO Y N. Experimental investigation of supercavitating motion of bodies [R]. Brussels: VKI, 2001.

[9] VASIN A D. Immersion of a cylinder in a fluid through a cylindrical free surface[J]. Fluid Dynamics, 2001, 36(2): 169-177.DOI: 10.1023/A:1019299930896.

[10]KUBENKO V D. Impact of a long thin body on a cylindrical cavity in liquid: a plane problem[J]. International Applied Mechanics, 2006, 42(6): 32-53. DOI: 10.1007/s10778-006-0131-y.

[11]KUBENKO V D, GAVRILENKO O V. Impact of a spherical rigid body on the surface of a cavity in a compressible liquid an axisymmetric problem [J]. International Applied Mechanics, 2008, 44(1): 8-15. DOI: 10.1007/s10778-008-0020-7.

[12]KUBENKO V D, GAVRILENKO O V. Impact interaction of cylindrical body with a surface of cavity during supercavitation motion in compressible fluid [J]. Journal of Fluids and Structures, 2009, 25(5): 794-814. DOI: 10.1016/j.jfluidstructs.2008.07.007.

[13]金大桥，王聪，魏英杰，等. 通气超空泡水下射弹实验研究[J]. 工程力学，2011，28(9)：214-217.

JIN Daqiao, WANG Cong , WEI Yingjie, et al. Experimental study of ventilated supercavity by underwater projectile[J]. Engineering Mechanics, 2011, 28(9): 214-217.

[14]曹伟，王聪，魏英杰，等. 自然超空泡形态特性的射弹试验研究[J]. 工程力学，2006，23(12)：175-187.DOI: 10.3969/j.issn.1000-4750.2006.12.031.

CAO Wei, WANG Cong, WEI Yingjie, et al. High-speed projectile experimental investigations on the characteristics of natural supercaviation[J]. Engineering Mechanics, 2006，23(12)：175-187.DOI: 10.3969/j.issn.1000-4750.2006.12.031.

[15]何乾坤，魏英杰，尤天庆，等. 空泡摆动对超空泡航行体尾拍影响分析[J]. 北京航空航天大学学报，2012, 38(4):509-512.

HE Qiankun, WEI Yingjie, YOU Tianqing, et al. Analysis of tail-slaps of supercavitating vehicle influenced by distortion of cavity shape[J]. Journal of Beijing University of Aeronautics and Astronautics，2012, 38(4):509-512.

[16]冯光，颜开. 超空泡航行体水下弹道的数值计算[J]. 船舶力学, 2005，9(2)：1-8. DOI: 10.3969/j.issn.1007-7294.2005.02.001.

FENG Guang, YAN Kai. Numerical calculation of underwater trajectory of supercavitating bodies[J]. Journal of Ship Mechanics, 2005，9(2)：1-8. DOI: 10.3969/j.issn.1007-7294.2005.02.001.

[17]MENTER F R. Two-equation eddy-viscosity turbulence models for engineering applications [J]. AIAA-Journal, 1994, 32(8): 1598-1605.

(编辑 张 红)

Analysis of tail-slapping and ballistic characteristics of supercavitating projectiles under different initial disturbances

ZHAO Chenggong, WANG Cong, SUN Tiezhi, ZHANG Xiaoshi

(School of Astronautics, Harbin Institute of Technology, Harbin 150001, China)

In order to investigate the influence of initial disturbances on the cavity profile, tail-slapping and ballistic characteristics, the computational fluid dynamics(CFD) model that can be used to solve the coupling equations of rigid body dynamics and unsteady reynolds averaged navier-stokes (URANS) of fluid under different initial disturbances is established based on the CFD program CFX and its expression language CEL. The cavity profile pressure distribution of the projectile tail and ballistic characteristic when the tail-slapping occurs are studied using this model. The result shows that the tail-slapping of the projectile destroys the symmetry of cavity profile; the initial disturbances have slight influences on the velocity and displacements inx-direction but have obvious influences on the vertical displacements; the amplitudes and frequencies of the pitch angle, angular velocity and angular acceleration increase with the increase of the initial disturbances.

multiphase flow; supercavitating projectile; tail-slapping; initial disturbance; CFD

10.11918/j.issn.0367-6234.2016.10.010

2015-11-16

中央高校基本科研业务费专项基金(HIT.NSRIF.201159)

赵成功(1986—)，男，博士研究生；

王 聪(1966—)，男，教授，博士生导师

王 聪，alanwang@hit.edu.cn

O351

A

0367-6234(2016)10-0071-06