基于氣動(dòng)耦合的導(dǎo)彈適配器初始分離彈道數(shù)值分析

牛鈺森,姜　毅,李　靜,陳　苗,史少巖,于邵禎

(1.北京理工大學(xué) 宇航學(xué)院,北京 100081;2.海軍裝備研究院,北京 100161)

?

基于氣動(dòng)耦合的導(dǎo)彈適配器初始分離彈道數(shù)值分析

牛鈺森1,姜毅1,李靜1,陳苗1,史少巖1,于邵禎2

(1.北京理工大學(xué) 宇航學(xué)院,北京 100081;2.海軍裝備研究院,北京 100161)

摘要:為了準(zhǔn)確預(yù)測(cè)與導(dǎo)彈分離后導(dǎo)彈適配器的初始分離彈道,使用基于局部重構(gòu)方法的動(dòng)網(wǎng)格技術(shù),將適配器的運(yùn)動(dòng)與流場(chǎng)的變化耦合。通過(guò)CFD方法實(shí)時(shí)求解適配器飛行過(guò)程中的氣動(dòng)載荷,并對(duì)適配器的六自由度運(yùn)動(dòng)微分方程進(jìn)行求解。使用基于尺寸函數(shù)的局部重構(gòu)方法對(duì)適配器周圍網(wǎng)格進(jìn)行更新。通過(guò)數(shù)值計(jì)算得到了適配器在初始分離彈道階段的氣動(dòng)載荷系數(shù)曲線,以及適配器的分離速度、角速度曲線和分離飛行軌跡,并與試驗(yàn)結(jié)果進(jìn)行了對(duì)比,兩者高度一致。結(jié)果表明,在分離初速的作用下，適配器與導(dǎo)彈的距離逐漸增大，在氣動(dòng)力作用下適配器發(fā)生偏轉(zhuǎn)，直到初始分離彈道結(jié)束適配器未與導(dǎo)彈發(fā)生碰撞。

關(guān)鍵詞:適配器;CFD;六自由度;動(dòng)網(wǎng)格;氣動(dòng)耦合

安裝在導(dǎo)彈與發(fā)射箱之間的適配器能夠?yàn)閷?dǎo)彈提供支撐,并且在發(fā)射過(guò)程中起到減震與導(dǎo)向的作用。對(duì)于安裝在導(dǎo)彈外表面的實(shí)心適配器,在導(dǎo)彈發(fā)射過(guò)程中隨著導(dǎo)彈一起運(yùn)動(dòng),出箱后通過(guò)爆炸螺栓或者彈簧機(jī)構(gòu)與導(dǎo)彈解鎖分離。在初始分離彈道階段,適配器的飛行軌跡和姿態(tài)受氣動(dòng)載荷的作用時(shí)刻發(fā)生變化,有可能與導(dǎo)彈發(fā)生碰撞造成發(fā)射事故。因此,在發(fā)射之前有必要對(duì)適配器的初始分離彈道進(jìn)行計(jì)算分析。

傳統(tǒng)的非耦合方法計(jì)算數(shù)個(gè)攻角下適配器的氣動(dòng)系數(shù),得到一系列離散的氣動(dòng)載荷數(shù)據(jù),擬合出氣動(dòng)載荷隨姿態(tài)的變化曲線,然后將曲線作為外力輸入,通過(guò)求解剛體的六自由度運(yùn)動(dòng)微分方程計(jì)算適配器的飛行姿態(tài)與軌跡。以上方法雖具有可操作性,但是為了得到精確的計(jì)算結(jié)果,需要計(jì)算大量的工況以得到更多的數(shù)據(jù)點(diǎn),耗費(fèi)時(shí)間和精力。并且,在真實(shí)的環(huán)境中適配器周圍的流場(chǎng)會(huì)隨著適配器與導(dǎo)彈間的距離發(fā)生變化,這種變化過(guò)程是傳統(tǒng)的非耦合計(jì)算方法不能實(shí)現(xiàn)的。

本文使用基于重構(gòu)方法的六自由度動(dòng)網(wǎng)格技術(shù),將適配器的運(yùn)動(dòng)狀態(tài)與流場(chǎng)的變化相耦合。每一時(shí)刻適配器受到的氣動(dòng)載荷由CFD程序直接求解,然后由剛體的六自由度運(yùn)動(dòng)微分方程計(jì)算出適配器的位移與角度變化。適配器的運(yùn)動(dòng)狀態(tài)更新后成為新的初始條件影響下一時(shí)刻氣動(dòng)載荷的計(jì)算。如此循環(huán)迭代可以真正實(shí)時(shí)求解適配器在初始分離彈道階段的飛行姿態(tài)與軌跡。使用耦合方法計(jì)算剛體在氣動(dòng)力驅(qū)動(dòng)作用下的運(yùn)動(dòng)過(guò)程,已在國(guó)內(nèi)外關(guān)于飛機(jī)外掛物的分離彈道研究中得到了廣泛應(yīng)用[2-4],其正確性得到了驗(yàn)證。

1適配器模型

適配器安裝在導(dǎo)彈中部靠近重心的位置處,總共4個(gè),并以導(dǎo)彈中心線為軸線按90°陣列。發(fā)射筒與地面垂直,導(dǎo)彈被彈射出筒,整個(gè)發(fā)射系統(tǒng)模型如圖1所示。

圖1　發(fā)射系統(tǒng)示意圖

適配器整體呈長(zhǎng)方形,靠近彈頭一端有向內(nèi)的楔形切面。內(nèi)外表面均呈圓弧形,以便與導(dǎo)彈和發(fā)射筒壁面貼合。適配器中心位置有放置彈簧結(jié)構(gòu)的盲孔,其外形如圖2所示。

圖2　導(dǎo)彈適配器幾何外形

發(fā)射前適配器被發(fā)射筒和彈體約束,彈簧處于壓縮狀態(tài)。導(dǎo)彈彈射出筒的過(guò)程中適配器與導(dǎo)彈一起運(yùn)動(dòng),當(dāng)適配器離開(kāi)發(fā)射筒,失去了筒壁的約束作用,彈簧開(kāi)始伸展使適配器產(chǎn)生沿導(dǎo)彈徑向的速度。適配器的彈出距離x與彈簧工作時(shí)間t的關(guān)系為

(1)

式中:l為彈簧的壓縮量,k為彈簧的剛度,m為適配器的質(zhì)量。由此計(jì)算出彈簧完全伸展開(kāi)需要的時(shí)間為10 ms,在此時(shí)間內(nèi)適配器的速度很小,因此氣動(dòng)載荷對(duì)適配器產(chǎn)生的沖量很小,可以忽略不計(jì)。因?yàn)檫m配器彈出的時(shí)間極短,所以在此時(shí)間內(nèi)導(dǎo)彈與適配器沿軸線方向的速度也變化很小,可以認(rèn)為兩者的速度保持相同。以彈簧完全展開(kāi)時(shí)為計(jì)算的初始時(shí)刻,此時(shí)導(dǎo)彈與適配器沿軸線方向的速度為20 m/s,適配器沿徑向的速度為5.8 m/s。將適配器下落到導(dǎo)彈尾部高度時(shí)指定為初始分離彈道的結(jié)束時(shí)刻。

計(jì)算適配器的平移運(yùn)動(dòng)選擇全局坐標(biāo)系較為方便,計(jì)算適配器的旋轉(zhuǎn)運(yùn)動(dòng)選擇隨體坐標(biāo)系較為方便。全局坐標(biāo)系以大地為基準(zhǔn),垂直于地面向上為X軸正方向,垂直于X軸且在流場(chǎng)模型的對(duì)稱平面內(nèi)為Y軸,Z軸由右手螺旋定則得到。適配器的隨體坐標(biāo)系以適配器中心孔為基準(zhǔn);經(jīng)過(guò)重心垂直于中心孔旋轉(zhuǎn)軸且在適配器對(duì)稱平面內(nèi)為X′軸,指向適配器楔形面為正方向;垂直于X′軸經(jīng)過(guò)重心與中心孔旋轉(zhuǎn)軸平行為Y′軸,指向中心孔開(kāi)口側(cè)為正方向;Z′軸由右手螺旋定則確定。

設(shè)隨體坐標(biāo)系由全局坐標(biāo)系先繞Z軸旋轉(zhuǎn)φ角,再繞Y軸旋轉(zhuǎn)θ角,然后繞X軸旋轉(zhuǎn)ψ角得到。則從全局坐標(biāo)系到隨體坐標(biāo)系的轉(zhuǎn)換矩陣R為

(2)

r11=cosθcosψ, r12=cosθsinψ, r13=-sinθ

r21=sinφsinθcosψ-cosφsinψ

r22=sinφsinθcosψ+cosφcosψ, r23=sinφcosθ

r31=cosφsinθcosψ+sinφsinψ

r32=cosψsinθsinψ-sinφcosψ, r33=cosφcosθ

使用Gambit建模工具以四面體非結(jié)構(gòu)網(wǎng)格對(duì)整個(gè)流場(chǎng)計(jì)算域進(jìn)行網(wǎng)格劃分,網(wǎng)格單元數(shù)量為200萬(wàn)。為保證適配器氣動(dòng)載荷計(jì)算的精度,對(duì)適配器周圍區(qū)域進(jìn)行局部加密。并且在計(jì)算過(guò)程加密區(qū)域隨適配器一起運(yùn)動(dòng),以確保適配器壁面附近的網(wǎng)格密度如圖3所示。

圖3　適配器與導(dǎo)彈周圍網(wǎng)格

2理論模型

2.1運(yùn)動(dòng)微分方程

適配器在空中飛行的過(guò)程中,除受到的氣動(dòng)力和重力外沒(méi)有其他約束條件。因此,具有6個(gè)方向上的自由度。將適配器的運(yùn)動(dòng)分解為兩部分,即質(zhì)心的質(zhì)點(diǎn)運(yùn)動(dòng)以及剛體的轉(zhuǎn)動(dòng)。

在全局坐標(biāo)系下建立質(zhì)心的質(zhì)點(diǎn)運(yùn)動(dòng)微分方程:

(3)

式中:FX、FY、FZ分別為X、Y、Z方向上適配器受到的作用力,vX、vY、vZ分別為這3個(gè)方向上的速度分量。在適配器的隨體坐標(biāo)系下建立剛體轉(zhuǎn)動(dòng)的運(yùn)動(dòng)微分方程:

L·ω=∑M-ω×(L·ω)

(4)

式中:L為適配器的慣性張量,ω為角速度向量,M為作用在適配器上的力矩。Admas-Moulton多步積分法非常適合求解運(yùn)動(dòng)微分方程,因此使用Admas-Moulton隱式三步四階格式對(duì)以上微分方程進(jìn)行數(shù)值求解。計(jì)算過(guò)程中式(3)與式(4)中的力與力矩通過(guò)積分適配器表面的壓強(qiáng)與切應(yīng)力得到,而壓強(qiáng)與切應(yīng)力由CFD程序解算流場(chǎng)控制方程得到。

2.2網(wǎng)格重構(gòu)方法

目前常用的動(dòng)網(wǎng)格方法有:光滑方法、分層方法和重構(gòu)方法。光滑方法適用于網(wǎng)格變形較小的情況,分層方法適用于網(wǎng)格運(yùn)動(dòng)方向恒定的情況。而適配器在空中飛行的過(guò)程中,速度方向時(shí)刻發(fā)生變化,且有滾轉(zhuǎn)、俯仰和偏航運(yùn)動(dòng),周圍網(wǎng)格的變形量較大,因此選用重構(gòu)動(dòng)網(wǎng)格方法最為合適。

基于恒定尺寸閾值的局部網(wǎng)格重構(gòu)方法靈活性較差,有可能造成剛體運(yùn)動(dòng)方向前方的網(wǎng)格尺寸過(guò)小,而運(yùn)動(dòng)方向后方的網(wǎng)格尺寸過(guò)大,降低計(jì)算精度甚至引起迭代發(fā)散。并且隨著剛體的運(yùn)動(dòng)可能導(dǎo)致計(jì)算域內(nèi)網(wǎng)格單元不斷增多,使計(jì)算速度逐漸降低。而采用基于尺寸函數(shù)的網(wǎng)格重構(gòu)方法可以有效避免以上問(wèn)題。其基本原理是,根據(jù)網(wǎng)格單元中心點(diǎn)與邊界的距離以及邊界網(wǎng)格單元的尺寸,計(jì)算此單元網(wǎng)格重構(gòu)時(shí)的尺寸基準(zhǔn)。重構(gòu)的網(wǎng)格尺寸為邊界網(wǎng)格單元尺寸的加權(quán)平均,其權(quán)值是重構(gòu)網(wǎng)格單元與邊界網(wǎng)格單元距離的反比,即:

sp=sb×γ

(5)

式中:

(6)

(7)

式中:FI與邊界上的單元格尺寸相關(guān);LI為網(wǎng)格單元中心點(diǎn)距離節(jié)點(diǎn)的距離;db是一個(gè)無(wú)量綱量,表征網(wǎng)格單元與邊界單元之間的距離;α為尺寸函數(shù)變化率,控制網(wǎng)格單元尺寸相對(duì)于邊界網(wǎng)格單元的大小;β為尺寸函數(shù)比例,控制網(wǎng)格單元尺寸隨db的變化程度。當(dāng)運(yùn)動(dòng)邊界靠近時(shí),db減小,sp也減小,此時(shí)重構(gòu)單元網(wǎng)格加密,當(dāng)運(yùn)動(dòng)邊界遠(yuǎn)離時(shí),db增大,sp也增大,此時(shí)重構(gòu)單元網(wǎng)格稀疏。

2.3流場(chǎng)控制方程

在基于歐拉描述的流場(chǎng)中,任意空間位置處的標(biāo)量φ在當(dāng)?shù)氐臅r(shí)間變化率是由:流入流出此控制體的對(duì)流流量,溫度梯度、濃度梯度等帶來(lái)的擴(kuò)散流量以及其他現(xiàn)象產(chǎn)生的源項(xiàng)引起的。因此標(biāo)量φ的輸運(yùn)方程可表示為

(8)

p=ρRT

(9)

h=∫cp(T)dT

(10)

式中:ρ為當(dāng)?shù)亓鲌?chǎng)的密度,u為流場(chǎng)的速度矢量,為拉普拉斯算子,Γ為擴(kuò)散系數(shù),Sφ為標(biāo)量φ的產(chǎn)生源項(xiàng)。當(dāng)φ=1時(shí),可以得到流場(chǎng)的質(zhì)量守恒方程;當(dāng)φ分別取X、Y、Z方向上的速度標(biāo)量u、v、w時(shí),可以得到動(dòng)量守恒方程;當(dāng)φ取比焓h時(shí),可以得到能量守恒方程。以上5個(gè)方程再加上描述氣體狀態(tài)的理想氣體狀態(tài)的式(9)以及描述比焓與溫度之間關(guān)系的式(10),就可以構(gòu)成封閉的控制方程組,對(duì)基本標(biāo)量ρ、u、v、w、p、T與h進(jìn)行求解。

2.4Spalart-Allmaras湍流模型

在適配器的氣動(dòng)流場(chǎng)中,適配器的楔形面與立面的交接過(guò)渡處存在曲率較大的流線以及較強(qiáng)的壓強(qiáng)梯度,當(dāng)適配器的攻角較大時(shí)還會(huì)存在貼體流動(dòng)分離現(xiàn)象。對(duì)于這樣的氣動(dòng)湍流問(wèn)題,Spalart-Allmaras模型經(jīng)過(guò)一定的優(yōu)化,是專門(mén)為氣動(dòng)流場(chǎng)設(shè)計(jì)的湍流模型。且在本文研究的適配器初始分離彈道問(wèn)題中導(dǎo)彈發(fā)動(dòng)機(jī)沒(méi)有點(diǎn)火,因此不存在燃?xì)馍淞鬟@樣的自由剪切流,選用此模型作為湍流粘性模型對(duì)于本文研究的問(wèn)題十分合適。Spalart-Allmaras模型是一方程湍流粘性模型,通過(guò)輸運(yùn)方程求解湍流粘性νT,其表達(dá)式為

(11)

式中:σν為常數(shù),Sν為湍流粘性源項(xiàng)。

3計(jì)算結(jié)果

初始時(shí)刻流場(chǎng)內(nèi)為標(biāo)準(zhǔn)大氣環(huán)境,壓強(qiáng)為101.325kPa,溫度為293.15K。流場(chǎng)中有沿Y軸正方向速度為3m/s的橫風(fēng)。按照初始時(shí)刻全局坐標(biāo)系下適配器質(zhì)心位置給適配器命名,質(zhì)心在Y軸正方向的為“SPQ-1”,質(zhì)心在Y軸負(fù)方向的為“SPQ-2”,質(zhì)心在Z軸正方向的為“SPQ-3”,質(zhì)心在Z軸負(fù)方向的為“SPQ-4”。

3.1適配器氣動(dòng)載荷變化

初始時(shí)刻適配器“SPQ-1”、“SPQ-2”迎風(fēng)面的壓力分布云圖如圖4所示。

圖4　初始時(shí)刻適配器迎風(fēng)面壓力分布

由圖4可以看出，“SPQ-2”迎風(fēng)面平均壓強(qiáng)要高于“SPQ-1”。選取適配器楔形面在隨體坐標(biāo)系O′Y′Z′平面上的投影面積S為參考面積,選取適配器最長(zhǎng)邊的一半為參考長(zhǎng)度l,由公式(12)可得適配器無(wú)量綱參數(shù)氣動(dòng)載荷系數(shù)。

(12)

式中:CF為作用在適配器上的氣動(dòng)力系數(shù),CM為作用在適配器上的氣動(dòng)力矩系數(shù),q∞為無(wú)窮遠(yuǎn)處自由流場(chǎng)的動(dòng)壓。適配器氣動(dòng)力系數(shù)CFX、CFY、CFZ在初始分離彈道階段隨時(shí)間變化曲線如圖5所示。

圖5　氣動(dòng)力系數(shù)曲線

適配器有沿X軸正方向的速度,因此受到沿X軸負(fù)方向的氣動(dòng)力作用。可以看出4個(gè)適配器X方向的氣動(dòng)力變化趨勢(shì)一致。從初始時(shí)刻到0.15s這段時(shí)間內(nèi)適配器氣動(dòng)力系數(shù)CFX幅值不斷增大,因?yàn)槭軋D5(b)、圖5(c)所示的氣動(dòng)力矩作用,在此時(shí)間內(nèi)適配器繞O′Z′轉(zhuǎn)動(dòng),隨著旋轉(zhuǎn)角度的增大,適配器在OYZ平面內(nèi)的投影面積不斷增大,使X方向上的氣動(dòng)力作用面積增大,使FX幅值上升。0.15s左右4個(gè)適配器繞O′Z′旋轉(zhuǎn)角度達(dá)到90°,適配器處于水平狀態(tài)如圖11(c)所示,之后適配器繼續(xù)旋轉(zhuǎn),OYZ平面內(nèi)的投影面積不斷減小,同時(shí)受FX作用vX不斷減小,2個(gè)因素共同作用下FX幅值下降如圖5(a)所示。在0.3s左右時(shí)適配器的旋轉(zhuǎn)角度達(dá)到180°,如圖11(d)所示,此時(shí)OYZ平面內(nèi)的投影面積達(dá)到最小,且vX小于初始時(shí)刻,因此FX幅值最小。在Y、Z方向上氣動(dòng)力的變化過(guò)程與此類似。

不同時(shí)刻適配器表面的壓力分布如圖6所示。初始時(shí)刻,適配器處于豎直狀態(tài),此時(shí)X方向氣動(dòng)力主要作用在適配器前端的楔形面上,如圖6(a)所示。楔形面外法線與O′X′軸有夾角,因此楔形面受到的氣動(dòng)力有沿O′Y′軸正方向分量,在適配器質(zhì)心處產(chǎn)生的力矩沿O′Z′軸正向,這是致使適配器旋轉(zhuǎn)的主要因素。隨體坐標(biāo)系下適配器受到的氣動(dòng)力矩系數(shù)CMZ′隨時(shí)間變化曲線如圖7所示。

圖6　適配器表面壓力云圖

圖7　氣動(dòng)載荷系數(shù)CMZ′曲線

可以看出，在最初的一段時(shí)間內(nèi)適配器受到的氣動(dòng)力矩有所上升,這是因?yàn)檫m配器繞O′Z′軸旋轉(zhuǎn)后增大了X方向的氣動(dòng)力作用面積,適配器迎風(fēng)面位于質(zhì)心之前的區(qū)域提供了更多的沿O′Z′軸正向的氣動(dòng)力矩,但是位于質(zhì)心之后的區(qū)域產(chǎn)生的是沿O′Z′軸負(fù)向的氣動(dòng)力矩,兩者疊加氣動(dòng)力矩幅值有較小的上升。當(dāng)適配器轉(zhuǎn)過(guò)一定的角度時(shí),質(zhì)心后部區(qū)域產(chǎn)生的負(fù)向氣動(dòng)力矩超過(guò)前部正向氣動(dòng)力矩。并且隨著適配器轉(zhuǎn)角的增大，楔形面外法線與來(lái)流之間的夾角不斷增大,由此產(chǎn)生的正向氣動(dòng)力矩逐漸減小,在2個(gè)因素的共同作用下致使適配器沿O′Z′軸正向的氣動(dòng)力矩隨轉(zhuǎn)角的增大而減小。在0.15s作用時(shí)適配器旋轉(zhuǎn)至水平位置,此時(shí)CMZ′達(dá)到最小值,之后適配器繼續(xù)旋轉(zhuǎn)適配器背面和底部成為迎風(fēng)面,從而產(chǎn)生了正向氣動(dòng)力矩使CMZ′有所上升,但是此時(shí)適配器的運(yùn)動(dòng)速度已經(jīng)減小很多,因此上升幅度不大。

3.2適配器分離彈道軌跡

由圖5、圖7可以看出，適配器“SPQ-3”與“SPQ-4”的氣動(dòng)載荷系數(shù)體現(xiàn)出較好的對(duì)稱性與一致性,而適配器“SPQ-1”與“SPQ-2”的氣動(dòng)載荷系數(shù)有較大的差別。這是因?yàn)檫m配器“SPQ-1”的分離速度沿Y軸正方向與風(fēng)的方向相同,而適配器“SPQ-2”的分離速度沿Y軸負(fù)方向與風(fēng)的方向相反,因此兩者相對(duì)于空氣的運(yùn)動(dòng)速度不同,造成兩者的氣動(dòng)載荷也不相同。而氣動(dòng)載荷的非對(duì)稱性又加劇了兩者Y方向上速度的差異。適配器質(zhì)心運(yùn)動(dòng)速度隨時(shí)間變化曲線如圖8所示。

可以看出,4個(gè)適配器在X方向上速度變化趨勢(shì)保持一致。適配器“SPQ-2”直接受到從Y軸負(fù)方向吹來(lái)的風(fēng)影響,vY的減小速度快于適配器“SPQ-1”。同時(shí)受風(fēng)力作用,本來(lái)在Y方向沒(méi)有速度的適配器“SPQ-3”與“SPQ-4”也逐漸產(chǎn)生了速度。在Z方向上沒(méi)有風(fēng)力作用,因此適配器“SPQ-3”與“SPQ-4”的速度vZ變化趨勢(shì)完全一致,且適配器“SPQ-1”與“SPQ-2”在Z方向上的速度始終為0。適配器質(zhì)心運(yùn)動(dòng)軌跡呈拋物線形如圖9所示。

圖8　適配器質(zhì)心速度曲線

圖9　適配器質(zhì)心軌跡

適配器與導(dǎo)彈分離0.575s后導(dǎo)彈底部在X方向上超過(guò)所有適配器,初始分離彈道結(jié)束。設(shè)此時(shí)適配器質(zhì)心與導(dǎo)彈軸線的垂直距離為s,結(jié)束時(shí)所有適配器的速度與位移見(jiàn)表1。

表1　初始分離彈道終點(diǎn)適配器運(yùn)動(dòng)參數(shù)

分離過(guò)程中繞隨體坐標(biāo)系O′Z′軸正方向的旋轉(zhuǎn)是主要的旋轉(zhuǎn)運(yùn)動(dòng),致使適配器前端楔形面向彈體方向翻轉(zhuǎn),旋轉(zhuǎn)角度ψ與角速度ωZ隨時(shí)間的變化曲線如圖10所示。

初始分離彈道結(jié)束時(shí)所有適配器的旋轉(zhuǎn)角速度與角度見(jiàn)表1。數(shù)值計(jì)算得到的適配器運(yùn)動(dòng)姿態(tài)與實(shí)驗(yàn)中高速攝影記錄的不同時(shí)刻適配器姿態(tài)變化對(duì)比如圖11所示。

由數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果均可以看出,初始分離彈道結(jié)束時(shí)適配器未與導(dǎo)彈發(fā)生碰撞,適配器的分離是安全的。此外,數(shù)值計(jì)算預(yù)測(cè)的適配器飛行姿態(tài)和軌跡與高速攝影記錄的實(shí)驗(yàn)結(jié)果保持一致。

圖10　適配器旋轉(zhuǎn)運(yùn)動(dòng)變化曲線

圖11　適配器飛行姿態(tài)

5結(jié)束語(yǔ)

本文建立了適配器與導(dǎo)彈的初始分離彈道模型,劃分了非結(jié)構(gòu)網(wǎng)格,使用基于六自由度動(dòng)網(wǎng)格技術(shù)將適配器的分離運(yùn)動(dòng)與氣動(dòng)流場(chǎng)相耦合,得到了初始分離彈道階段適配器的運(yùn)動(dòng)軌跡與姿態(tài),并與實(shí)驗(yàn)結(jié)果進(jìn)行了對(duì)比,可以得出以下結(jié)論:

①計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果高度一致,表明此種研究方法具有較高的可信度和準(zhǔn)確性;

②在本文研究的工況中, 適配器在初始分離階

段沒(méi)有與導(dǎo)彈發(fā)生碰撞,分離過(guò)程是安全的;

③受氣動(dòng)載荷的作用適配器產(chǎn)生旋轉(zhuǎn)運(yùn)動(dòng),且適配器前端楔形面向彈體方向翻轉(zhuǎn);

④沿著風(fēng)速方向?qū)ΨQ分布的2個(gè)適配器,受到的氣動(dòng)載荷不同,兩者的運(yùn)動(dòng)狀態(tài)也不同。

本文使用的研究方法對(duì)于研究適配器與導(dǎo)彈的分離過(guò)程,預(yù)防適配器與導(dǎo)彈的碰撞,以及地面設(shè)備和人員的防護(hù)有指導(dǎo)作用。

參考文獻(xiàn)

趙華,王敏杰,楊為,等.箱式發(fā)射導(dǎo)彈適配器.戰(zhàn)術(shù)導(dǎo)彈技術(shù),2007(4):42-50.

ZHAO Hua,WANG Min-jie,YANG Wei,et al.Adapters for canister-launched missile.Tactical Missile Technology,2007(4):42-50.(in Chinese)

MAHMOOD T,AIZUD N,ZAHIR S.CFD simulations of the store separation and its effect on the attitude of the parent vehicle//International Bhurban Conference on Applied Sciences.Piscataway,NJ,USA:IEEE,2012:198-202.

PANAGIOTOPOULOS E E,KYPARISSIS S D.CFD transonic store separation trajectory predictions with comparison to wind tunnel investigations.International Journal of Engineering,2010(3):538-553.

張培紅,王明,鄧有奇,等.武器分離及艙門(mén)開(kāi)啟過(guò)程數(shù)值模擬研究.空氣動(dòng)力學(xué)學(xué)報(bào),2013,31(3):277-293.

ZHANG Pei-hong,WANG Ming,DENG You-qi,et al.Numerical simulation of store separation and door operation.Acta Aerodynamic Sincia,2013,31(3):277-293.(in Chinese)

BURDEN R L,FAIRES J D.Numerical analysis.Youngstown:Youngstown State University,2010.

POPE S B.Turbulent flows.UK:Cambridge University Press,2000.

ANDERSON J D.Fundamentals of aerodynamics.5th ed.New York:McGraw-Hill,2009.

Numerical Analysis of Initial Separation Trajectory of

Missile Adapter Based on Aerodynamic Coupling

NIU Yu-sen1,JIANG Yi1,LI Jing1,CHEN Miao1,SHI Shao-yan1,YU Shao-zhen2

(1.School of Aerospace Engineering,BIT,Beijing 100081,China;

2.Naval Acadme of Armament,Beijing 100161,China)

Abstract:In order to accurately predict the initial separation trajectory of the adapter after releasing from the missile,the motion of the adapter was coupled with the variation of the flow field by dynamic mesh technology based on local reconstruction.The aerodynamic loads of the adapter were computed in real time by CFD routine.The 6DOF motion differential equations of the adapter were solved.The models of missile and four adapters were built,and the computational domain was meshed by using tetrahedron cells.The grid cells around the adapter were updated by using local remeshing method based on size function.The initial separation trajectory of the adapter was calculated.The results accord with the experimental result.The curves of motion velocities and aerodynamic loads coefficients were drawn.Results show that the distance between adapter and missile increases because of the initial separation velocity,and the adapters rotate because of the aerodynamic force.The adapter doesn’t collide with missile untile the end of initial separation trajectory.

Key words:adapter;CFD;six degree of freedom;dynamic mesh;aerodynamic coupling.

中圖分類號(hào):TJ303.4

文獻(xiàn)標(biāo)識(shí)碼:A

文章編號(hào):1004-499X(2015)04-0052-07

作者簡(jiǎn)介:牛鈺森(1987- ),男,博士研究生,研究方向?yàn)楸靼l(fā)射理論與技術(shù)。E-mail:shenzhou1987@aliyun.com。

收稿日期:2015-06-11