顧國華，楊 溢，王積忠，陳 浩

（陸軍軍官學(xué)院，合肥 230031）

0 引言

空氣動力學(xué)研究表明，當(dāng)旋成體（如彈丸）以超聲速在空中飛行時將產(chǎn)生激波，不同物體、不同飛行條件產(chǎn)生的激波不同，激波狀況體現(xiàn)了飛行物的氣動特征。為研究旋成體氣動特性，國內(nèi)外學(xué)者和團(tuán)隊從不同角度和方向在理論和實(shí)驗上進(jìn)行了較為深入的研究［1-4］。

從公開報道看，Munk［5］較早提出了旋成體氣動特性，提出基于細(xì)長體位勢流理論并進(jìn)行數(shù)值模擬解算，可求解 0～5°迎角的氣動問題。Jorgense［6］將橫流理論氣動分解為勢流項和橫流粘性項，并對相關(guān)阻力進(jìn)行修正和分解。后來逐漸發(fā)展了激波捕捉方法，如Harten1983年提出的TVD差分格式，該方法是一種具有二階精度的總變差不增算法；之后發(fā)展為對稱型TVD。近代高分辨率的激波捕捉算法得以長足發(fā)展，大體可歸納為以下4種：

①TVD差分格式的升級，ENO格式和WENO格式。②Godunov間斷分解方法。③FCT差分格式。④Roe差分格式。

為探索基于激波場張量的彈丸跟蹤，本文提出建立三維Euler方程求解來確立激波捕捉張量，通過激波捕捉陣列多聲信號融合處理，實(shí)現(xiàn)對高精度反演激波聲場位置解算和跟蹤。

1 彈丸激波場

由空氣動力學(xué)原理可知，當(dāng)彈丸以超聲速在空中飛行時，勢必造成對彈丸頭部周邊空氣的擠壓，而使得彈尾部空氣稀疏。因而在頭部形成壓縮波（高壓波），在尾部形成稀疏波。通常將其壓縮波稱為彈丸激波。是一種非線性大氣聲學(xué)現(xiàn)象，如圖1所示。由激波形成的某一區(qū)域稱為激波場。

在二維空間分析，激波的高壓區(qū)傳播速度大于聲速，而低壓區(qū)則低于聲速。因而，通常形成起伏狀波形，形似英文字母“N”，所以俗稱為“N”波。

彈丸“N”波的形成過程大體上可分為4個階段［7］如圖 2所示。

2 彈丸激波捕捉

2.1 激波張量分析

二維Euler方程是分析空氣動力學(xué)的基本模型。為了抑制振蕩和過沖現(xiàn)象，將通量項分解為對流通量項和壓力項，在文獻(xiàn)［8-9］的基礎(chǔ)上拓展，采用理想氣體模型建立三維Euler方程（1）。采用張量捕捉方法的優(yōu)勢在于信號更全面、更精確，利于分析處理。

其中，

其中，ρ為壓力，u，v和 w 分別為 x，y，z方向上的速度分量，e，E為單位體積內(nèi)的內(nèi)能和總能。

方程（1）可轉(zhuǎn)化為守恒型形式：

由式（1）、式（2）可得

解得非線性雅可比系數(shù)特征值矩陣ΛA如下，同理可解得 ΛB、ΛC。

解得A特征矢量矩陣RA如下，同理可解得B、C 的特征矢量矩陣 RB、RC。

2.2 激波捕捉累積多聲信號融合算法

根據(jù)文獻(xiàn)［10］建立的二維守恒型Euler方程，在任一點(diǎn)（x，y）處，將體積采用有限體積法可得：

其中，Rx，y表示流出控制體的凈通量，在某一定點(diǎn)（i，j）處，凈通量以通量形式可表示為

在α角度上，令信號空間U1與V1張成，噪聲子空間由U2與V2張成，將其進(jìn)行進(jìn)行奇異值分解［11］

以i方向的通量為例，將通量計算公式改為：

對于由M*N矩陣組成的聲矢量傳感器陣列

在i-1/2界面處的壓強(qiáng)p可表示為

進(jìn)行奇異值分解可得

式中：∑'1為非零奇異值矩陣。U'1和V'1張成信號子空間，U'2和V'2張成與信號子空間正交的噪聲子空間。

考慮到實(shí)際接收的數(shù)據(jù)矩陣長度有限，因此，對循環(huán)相關(guān)矩陣取最大似然估計

其中，L為快拍數(shù)。

3 彈丸運(yùn)動跟蹤

超聲速彈丸空中飛行時，運(yùn)動軌跡影響因素多、不易捕捉、具有一定隨機(jī)性。試驗與反演以某彈丸對空射擊背景，航路參數(shù)為：等速直線水平飛行1 000 m，高度500 m，捷徑400 m，目標(biāo)航速200 m/s，航向由北向南。通過測試與反演彈目偏差解算跟蹤彈丸，檢驗彈丸運(yùn)動跟蹤濾波效果。

經(jīng)激波聲場3 500個采樣點(diǎn)后提取彈丸運(yùn)動軌跡如圖3所示。采用ENO格式處理后的的軌跡如圖4所示，采用Godunov法處理后的軌跡如圖5所示，采用本文提出的激波張量捕捉法處理后的波形如下頁圖6所示。

仿真實(shí)驗中，通過比較分析，在3 500點(diǎn)采樣數(shù)據(jù)中，相對于初始數(shù)據(jù)，采用3種方法處理后的殘差曲線統(tǒng)計，采用ENO格式法的跟蹤誤差為18.5%，采用Godunov法的跟蹤誤差為15.6%，采用本文提出的激波張量捕捉法的跟蹤誤差減小到8.9%。

4 結(jié)論

當(dāng)前，采用彈丸激波場對彈丸進(jìn)行跟蹤的研究還不夠深入。本文探討了基于激波場張量的彈丸跟蹤仿真，實(shí)驗和反演表明，提出的激波張量捕捉法收斂效果好，精度較高，有效解決了復(fù)雜激波聲場中彈丸跟蹤精度和收斂速度的矛盾，為超聲速運(yùn)動目標(biāo)定位和跟蹤提供了一種新方法。

［1］URSZULA L，KRYSTIAN S.Wavelet based shock wave and muzzle blast classification for different supersonic projectiles［J］.Expert Systems with Applications，2014，41（11）：5097-5104.

［2］KNIGHT D，LONGO J，DRIKAKIS D，et al.Assessment of CFD capability for prediction of hypersonic shock interactions［J］.Progress in Aerospace sciences，2012（48）：8-26.

［3］GNANI F，LO K H，ZARE B H，et al.Shock wave diffraction in the presence of a supersonic co-flow jet［J］.Shock Waves，2016（3）：1-10.

［4］趙占龍.彈箭線性與非線性氣動特性計算［D］.南京:南京理工大學(xué)，2012.

［5］MUNK M M.The aerodynamic force on airship hulls［R］.NACA-TR-184，1923.

［6］JORGENSE L H.Prediction of static aerodynamic characteristics for slender bodies alone and with lifting surfaces to very high angles of attack［R］.NASA TR-R474，1977.

［7］張飛猛.對空射擊聲學(xué)靶脫靶量測試系統(tǒng)的精度分析［J］.兵工學(xué)報，2000，21（1）:23-26.

［8］張德良.計算流體力學(xué)教程［M］.北京:高等教育出版社，2010：259-260，289.

［9］WAGNER B，SCHMIDT W.Theoretical investigations on shock wave-boundary layer interaction in cryogenic nitrogen［J］.Springer International Publishing，2015，365（3）:341-350.

［10］賈勇，孫剛，劉蘇.一種改進(jìn)混合迎風(fēng)格式在翼型流場激波捕捉中的應(yīng)用［J］.力學(xué)季刊，2007，28（2）:223-227.

［11］鄭彤.復(fù)雜噪聲背景下基于聲矢量傳感器陣列的汽車噪聲源方位估計方法研究［D］.長春：吉林大學(xué)，2012.

［12］王勇，王學(xué)德，許國龍.超聲速彈丸高空流場數(shù)值模擬與稀薄氣體效應(yīng)研究［J］.兵器裝備工程學(xué)報，2016，37（10）:13-18.

［13］李樹軍，李開端.基于自適應(yīng)α-β濾波的目標(biāo)跟蹤算法研究［J］.科學(xué)技術(shù)與工程，2012（12）：9367-9369.

［14］顧國華，陳棟，王彩，等.基于聲學(xué)靶傳感器的彈著點(diǎn)測試研究與實(shí)現(xiàn)［J］.電子測量技術(shù)，2007，30（2）:154-156.