顾国华;杨溢;王积忠;陈浩
【摘 要】Shock wave of supersonic projectile contains plenty of information about its motion state. Based on the analysis of body of revolution shock wave formation mechanism and establish the relationship model between shock wave formation, the paper establishes the three-dimensional Euler model and fusion algorithm of acoustic signal based on shock capture, and then put into spatial spectrum identification experiment. Compared to background field test, it proposed dynamic g-h-k filter tracking methods to verify the effect of the tracking of the shock wave field inversion. Above all , it proposes a new method of tracking supersonic projectile based on shock wave field.%超声速弹丸产生的激波蕴含丰富的运动状态信息.在分析弹丸激波形成机理和建立激波形成前后参数关系的基础上,建立空间张量的激波捕捉三维Euler模型和声信号融合算法,并进行空间谱辨识实验;为实现反演跟踪,结合外场试验背景,提出动态g-h-k滤波跟踪方法,验证激波声场反演跟踪的有效性,为探索弹丸激波场进行弹丸跟踪提供了新方法.
【期刊名称】《火力与指挥控制》
【年(卷),期】2017(042)012
【总页数】4页(P31-34)
【关键词】超声速;弹丸;激波场;张量;跟踪
【作 者】顾国华;杨溢;王积忠;陈浩ADAU1452
【作者单位】陆军军官学院,合肥 230031;陆军军官学院,合肥 230031;陆军军官学院,合肥 230031;陆军军官学院,合肥 230031
【正文语种】中 文
【中图分类】TP391;TJ011+2
0 引言
空气动力学研究表明，当旋成体（如弹丸）以超声速在空中飞行时将产生激波，不同物体
、不同飞行条件产生的激波不同，激波状况体现了飞行物的气动特征。为研究旋成体气动特性，国内外学者和团队从不同角度和方向在理论和实验上进行了较为深入的研究［1-4］。
从公开报道看，Munk［5］较早提出了旋成体气动特性，提出基于细长体位势流理论并进行数值模拟解算，可求解 0～5°迎角的气动问题。Jorgense［6］将横流理论气动分解为势流项和横流粘性项，并对相关阻力进行修正和分解。后来逐渐发展了激波捕捉方法，如Harten1983年提出的TVD差分格式，该方法是一种具有二阶精度的总变差不增算法；之后发展为对称型TVD。近代高分辨率的激波捕捉算法得以长足发展，大体可归纳为以下4种：
大嘴泉①TVD差分格式的升级，ENO格式和WENO格式。②Godunov间断分解方法。③FCT差分格式。④Roe差分格式。
为探索基于激波场张量的弹丸跟踪，本文提出建立三维Euler方程求解来确立激波捕捉张量，通过激波捕捉阵列多声信号融合处理，实现对高精度反演激波声场位置解算和跟踪。
1 弹丸激波场
由空气动力学原理可知，当弹丸以超声速在空中飞行时，势必造成对弹丸头部周边空气的挤压，而使得弹尾部空气稀疏。因而在头部形成压缩波（高压波），在尾部形成稀疏波。通常将其压缩波称为弹丸激波。是一种非线性大气声学现象，如图1所示。由激波形成的某一区域称为激波场。
拟合值图1 超声速弹丸飞行产生的激波
在二维空间分析，激波的高压区传播速度大于声速，而低压区则低于声速。因而，通常形成起伏状波形，形似英文字母“N”，所以俗称为“N”波。
弹丸“N”波的形成过程大体上可分为4个阶段［7］如图 2所示。
图2 “N”波
2 弹丸激波捕捉
2.1 激波张量分析吉林市苏宁天润城
二维Euler方程是分析空气动力学的基本模型。为了抑制振荡和过冲现象，将通量项分解为
对流通量项和压力项，在文献［8-9］的基础上拓展，采用理想气体模型建立三维Euler方程（1）。采用张量捕捉方法的优势在于信号更全面、更精确，利于分析处理。
其中，
其中，ρ为压力，u，v和 w 分别为 x，y，z方向上的速度分量，e，E为单位体积内的内能和总能。
方程（1）可转化为守恒型形式：
由式（1）、式（2）可得
令可解得：
解得非线性雅可比系数特征值矩阵ΛA如下，同理可解得 ΛB、ΛC。
解得A特征矢量矩阵RA如下，同理可解得B、C 的特征矢量矩阵 RB、RC。
2.2 激波捕捉累积多声信号融合算法
根据文献［10］建立的二维守恒型Euler方程，在任一点（x，y）处，将体积采用有限体积法可得：
其中，Rx，y表示流出控制体的净通量，在某一定点（i，j）处，净通量以通量形式可表示为
在α角度上，令信号空间U1与V1张成，噪声子空间由U2与V2张成，将其进行进行奇异值分解［11］
摊丁入亩其中，为非零奇异值矩阵。
以i方向的通量为例，将通量计算公式改为：
对于由M*N矩阵组成的声矢量传感器阵列
在i-1/2界面处的压强p可表示为
对于三维矢量信号，令第页的矩阵为：
进行奇异值分解可得
式中：∑'1为非零奇异值矩阵。U'1和V'1张成信号子空间，U'2和V'2张成与信号子空间正交的噪声子空间。
oec管理模式考虑到实际接收的数据矩阵长度有限，因此，对循环相关矩阵取最大似然估计
其中，L为快拍数。
3 弹丸运动跟踪
超声速弹丸空中飞行时，运动轨迹影响因素多、不易捕捉、具有一定随机性。试验与反演以某弹丸对空射击背景，航路参数为：等速直线水平飞行1 000 m，高度500 m，捷径400 m，目标航速200 m/s，航向由北向南。通过测试与反演弹目偏差解算跟踪弹丸，检验弹丸运动跟踪滤波效果。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议