李寧波，邵 雷，王華吉，雷虎民，龍振國(guó)

(1.空軍工程大學(xué) 防空反導(dǎo)學(xué)院，西安 710051；2.中國(guó)人民解放軍95100部隊(duì)，廣州 510405)

0 引言

近年來(lái)，臨近空間高超聲速武器迅猛發(fā)展使其成為國(guó)家安全的現(xiàn)實(shí)威脅[1-2]，研究反臨近空間的先進(jìn)攔截技術(shù)和防御手段已迫在眉睫。從目前的研究現(xiàn)狀來(lái)看，國(guó)外關(guān)于反臨的研究文獻(xiàn)很少，其主要以進(jìn)攻型高超聲速飛行器為背景，但也可以為反臨研究提供借鑒。國(guó)內(nèi)方面，許多學(xué)者開(kāi)展了積極的探索，提出了有益的思路，對(duì)臨近空間高超聲速目標(biāo)攔截難點(diǎn)的認(rèn)識(shí)逐漸統(tǒng)一，并逐步形成空基攔截和地基攔截兩類主流思想[3-4]。但對(duì)于攔截彈的彈道設(shè)計(jì)目前還沒(méi)有定論，有待研究和完善。

面對(duì)臨近空間高超聲速武器這類新型、高動(dòng)態(tài)目標(biāo)，攔截彈應(yīng)采用復(fù)合制導(dǎo)體制，遠(yuǎn)程、梯次進(jìn)行攔截，從而保證攔截的成功率。由于攔截彈大部分時(shí)間飛行在中制導(dǎo)段，因此中制導(dǎo)段的飛行性能決定了整體性能。傳統(tǒng)防空導(dǎo)彈制導(dǎo)律研究主要基于彈目的運(yùn)動(dòng)關(guān)系設(shè)計(jì)導(dǎo)引律，彈道由導(dǎo)引律決定?？紤]到臨近空間攔截問(wèn)題中各種嚴(yán)格的限制條件，在中制導(dǎo)律的設(shè)計(jì)過(guò)程中必須要考慮到彈道的優(yōu)化設(shè)計(jì)，以提高攔截效率。

彈道設(shè)計(jì)問(wèn)題一直是飛行器研究領(lǐng)域的熱點(diǎn)。文獻(xiàn)[5-6]分別研究了初始及終點(diǎn)位置固定、飛行時(shí)間受約束和考慮地球偏率影響的大氣層外彈道規(guī)劃問(wèn)題。文獻(xiàn)[7]提出了使相對(duì)位置向量與速度向量共線的中制導(dǎo)飛行策略，研究了大氣層外中制導(dǎo)段的彈道。文獻(xiàn)[8]對(duì)防空導(dǎo)彈中制導(dǎo)段彈道優(yōu)化問(wèn)題進(jìn)行了較為全面的論述，同時(shí)也給出了準(zhǔn)最佳彈道的相關(guān)理論和實(shí)現(xiàn)方法。文獻(xiàn)[9]研究了虛擬域動(dòng)態(tài)逆的方法，并將其應(yīng)用于空地反輻射導(dǎo)彈的彈道優(yōu)化設(shè)計(jì)。

雖然許多學(xué)者對(duì)彈道優(yōu)化設(shè)計(jì)問(wèn)題進(jìn)行了深入的探索，但是對(duì)于以反臨攔截彈為背景的研究卻很少。本文以地基攔截臨近空間高超聲速目標(biāo)為背景，根據(jù)臨近空間目標(biāo)特點(diǎn)及攔截作戰(zhàn)要求，設(shè)計(jì)了一種高拋再入的彈道方案，并分析了中末制導(dǎo)交接班問(wèn)題從而為建立中制導(dǎo)終端約束模型的建立提供參考。

1 彈道方案設(shè)計(jì)

1.1 目標(biāo)彈道特性分析

目前，典型的高超聲速助推滑翔飛行器主要以HTV-2和AHW為代表，本文以HTV-2為攔截目標(biāo)進(jìn)行彈道特性分析。

Lu Ping等[10]研究了HTV-2的兩種彈道模式，一種是采用彈道式發(fā)射，助推段完成后的終端狀態(tài)與彈道導(dǎo)彈類似，試驗(yàn)飛行器先飛出大氣層，之后再入至臨近空間高度進(jìn)行平衡滑翔飛行；另一種在助推級(jí)利用發(fā)動(dòng)機(jī)推力直接將飛行軌跡調(diào)整為平衡滑翔模式。兩種模式的飛行彈道相比于彈道導(dǎo)彈，其飛行高度都大大的降低，與彈道導(dǎo)彈長(zhǎng)時(shí)間在大氣層外飛行不同，HTV-2將長(zhǎng)時(shí)間在大氣層內(nèi)的臨近空間高度飛行，且其通常采用空氣動(dòng)力控制[11]。不管采用何種模式，HTV-2都將長(zhǎng)時(shí)間在30～60 km高度的范圍內(nèi)滑翔飛行，且其飛行彈道較為平穩(wěn)。此外，通過(guò)分析速度剖面可發(fā)現(xiàn)，兩者的速度最高都可達(dá)到6000～7000 m/s，而在滑翔段的末端其速度一般只有2000 m/s左右。因此，從速度的角度看，滑翔段的末端部分是一個(gè)適合攔截的窗口。

另一方面，在滑翔段目標(biāo)航跡較為平穩(wěn)，易于探測(cè)發(fā)現(xiàn)并進(jìn)行穩(wěn)定跟蹤，攔截窗口較大，將該段作為目標(biāo)與攔截彈的交會(huì)時(shí)間段，有利于對(duì)目標(biāo)軌跡預(yù)測(cè)計(jì)算，目標(biāo)可用過(guò)載較小。在一定程度上可降低對(duì)攔截彈的過(guò)載要求。另外，從攔截彈的控制角度考慮，臨近空間下層高度范圍內(nèi)存在有稀薄大氣，有利于對(duì)攔截彈實(shí)施氣動(dòng)控制，減少直接力控制的使用。

1.2 彈道與導(dǎo)引律分離設(shè)計(jì)

在制導(dǎo)律設(shè)計(jì)中，防空導(dǎo)彈通常都以最速接近為目標(biāo)。通過(guò)數(shù)學(xué)方法可證明，前置法、半前置法、三點(diǎn)法等遙控指令制導(dǎo)方法，和追蹤法、平行接近法、比例導(dǎo)引法等尋的制導(dǎo)方法，其本質(zhì)上都是根據(jù)不同初始條件，以最速接近為性能指標(biāo)推導(dǎo)得到的。由于臨近空間攔截彈的作戰(zhàn)距離相對(duì)較遠(yuǎn)，如果采用上述的基于彈目相對(duì)運(yùn)動(dòng)關(guān)系的攔截導(dǎo)引策略，將會(huì)導(dǎo)致攔截武器制導(dǎo)指令的頻繁調(diào)整，造成更多能量損耗。

對(duì)高超聲速目標(biāo)攔截問(wèn)題，當(dāng)采用彈道與導(dǎo)引律分離設(shè)計(jì)策略時(shí)，中制導(dǎo)段攔截彈彈道無(wú)需隨目標(biāo)運(yùn)動(dòng)進(jìn)行實(shí)時(shí)調(diào)整，只有當(dāng)終端約束條件發(fā)生大的變化時(shí)才需相應(yīng)的彈道調(diào)整，以保證攔截彈以較好的相對(duì)陣位進(jìn)入末制導(dǎo)，因此中制導(dǎo)段合理設(shè)計(jì)優(yōu)化彈道有利于提高飛行品質(zhì)，從而達(dá)到有效命中目標(biāo)的目的[12]。

另一方面，攔截彈飛行過(guò)程中，其結(jié)構(gòu)強(qiáng)度、熱防護(hù)、發(fā)動(dòng)機(jī)正常工作以及控制穩(wěn)定性等一系列問(wèn)題對(duì)熱流密度、動(dòng)壓、過(guò)載、控制量等過(guò)程特征提出了嚴(yán)格的要求，上述導(dǎo)引律難以將各種約束考慮其中。綜上，為了克服以上的缺點(diǎn)，提升攔截彈的性能，應(yīng)將彈道和導(dǎo)引方法分離設(shè)計(jì)。

1.3 高拋再入彈道

由于地基攔截彈從地面發(fā)射，如果按照傳統(tǒng)彈道設(shè)計(jì)，將會(huì)導(dǎo)致攔截彈的飛行彈道全程處于目標(biāo)的下方，由于下方的空氣阻力相對(duì)較大，因此與目標(biāo)相比，攔截彈處于不利地位，也限制了攔截彈的射程。如果按照拋物線彈道飛行，先向上進(jìn)行爬升，然后再沿拋物線下降，末制導(dǎo)段的彈道按照尋的導(dǎo)引律飛行，此時(shí)大部位的彈道將會(huì)處于目標(biāo)上方，空氣阻力降低，攔截彈射程也將會(huì)大大增加。

此外，由于臨近空間及其上部空域大氣相對(duì)稀薄，其空氣中的紅外吸收物質(zhì)較少，所以導(dǎo)引頭的探測(cè)視角應(yīng)采用自上而下的方式[13]。另一方面，攔截彈長(zhǎng)期在目標(biāo)上方的稀薄大氣中飛行，可降低導(dǎo)引頭和彈體的氣動(dòng)加熱，為導(dǎo)引頭提供較好的初始工作條件。綜上，本文提出攔截彈采用高拋再入的攔截彈道方案，此處“再入”不是指?jìng)鹘y(tǒng)意義上的從大氣層外再入，而是指攔截彈向上爬升后在下降飛行的過(guò)程。

根據(jù)上述設(shè)計(jì)思想，攔截彈的彈道將按照不同的飛行過(guò)程劃分為如下四個(gè)階段：

(1)無(wú)控段：假設(shè)攔截彈采用垂直發(fā)射模式，冷彈升起以后，在空中進(jìn)行點(diǎn)火，垂直飛行一定高度后，制導(dǎo)雷達(dá)捕獲攔截彈，這段時(shí)間僅進(jìn)行攔截彈的穩(wěn)定控制。

(2)轉(zhuǎn)彎段：從啟控點(diǎn)開(kāi)始，到助推發(fā)動(dòng)機(jī)熄火關(guān)機(jī)，與彈體分離，本章的設(shè)計(jì)為兩級(jí)的火箭助推發(fā)動(dòng)機(jī)。該階段攔截彈的速度和彈道傾角的變化范圍較大。

(3)滑翔段：從助推發(fā)動(dòng)機(jī)關(guān)機(jī)分離開(kāi)始，到導(dǎo)引頭捕獲目標(biāo)，順利實(shí)現(xiàn)中末制導(dǎo)交接班為止。該階段攔截彈進(jìn)行無(wú)動(dòng)力滑行，不斷接近目標(biāo)，飛行狀態(tài)較為平穩(wěn)。

(4)末制導(dǎo)段：從順利實(shí)現(xiàn)中末制導(dǎo)交接班開(kāi)始，到精確命中目標(biāo)為止。該階段攔截彈在尋的導(dǎo)引律的作用下，直接飛向目標(biāo)，進(jìn)行碰撞殺傷。

上述過(guò)程中，無(wú)控段垂直上升，沒(méi)有必要對(duì)彈道進(jìn)行優(yōu)化。本文所研究的中制導(dǎo)段為轉(zhuǎn)彎段和滑翔段。圖1為攔截過(guò)程的彈道示意圖。

圖1 攔截彈全程彈道示意圖Fig.1 Schematic diagram of the whole trajectory

2 攔截彈彈道優(yōu)化模型

2.1 攔截彈運(yùn)動(dòng)方程

中制導(dǎo)段攔截彈應(yīng)朝向終端位置約束飛行，理想軌跡位于攔截彈與目標(biāo)所在的縱向平面內(nèi)。另一方面，在實(shí)際飛行的過(guò)程中，由攔截彈偏離預(yù)定彈道平面造成的側(cè)向運(yùn)動(dòng)參數(shù)的變化量相比于縱向參數(shù)的變化都是小量，其實(shí)際飛行彈道與縱向平面內(nèi)的彈道差別不大，所以可將縱向運(yùn)動(dòng)和側(cè)向運(yùn)動(dòng)分開(kāi)來(lái)研究[14]。因此，考慮如下縱向平面內(nèi)的攔截彈質(zhì)點(diǎn)運(yùn)動(dòng)模型：

(1)

(2)

(3)

(4)

式中v為攔截彈速度；θ為彈道傾角；x、y分別為攔截彈在大地慣性坐標(biāo)系下的位置；m為攔截彈質(zhì)量；q為動(dòng)壓；S為參考面積；g為重力加速度；P為作用在攔截彈上的推力。

為了使攔截彈達(dá)到較高的飛行速度，需要大推力助推加速，初步設(shè)計(jì)了兩級(jí)助推方案，推力計(jì)算公式如式(5)所示[15]：

(5)

式中Ve為燃?xì)馑俣龋籌sp為比沖，對(duì)固體火箭發(fā)動(dòng)機(jī)，其值一般為200～300 s。

發(fā)動(dòng)機(jī)的相關(guān)參數(shù)設(shè)計(jì)如表1所示。

表1 攔截彈助推發(fā)動(dòng)機(jī)相關(guān)參數(shù)

阻力系數(shù)及升力系數(shù)Cx、Cy可以分別表示為馬赫數(shù)Ma和攻角α的函數(shù)：

(6)

對(duì)大氣密度的假設(shè)：

ρ=ρ0exp(-y/y0)

(7)

式中ρ0為地面空氣密度，ρ0=1.225 0 kg/m3；ρ為飛行器高度處的空氣密度；y0=7254.3 m。

2.2 目標(biāo)函數(shù)與過(guò)程約束

設(shè)t0表示初始時(shí)刻，tf表示終端時(shí)刻，為保證攔截彈的殺傷效果，一般將終端時(shí)刻的速度最大作為優(yōu)化指標(biāo)J，即：

J=minφ(v(tf)，tf)=min{-vf}

(8)

為保證飛行的穩(wěn)定性，對(duì)攔截彈攻角大小應(yīng)加以限制?？刂屏考s束設(shè)置如下：

‖α‖≤αmax

(9)

動(dòng)壓極限值主要取決于熱防護(hù)材料強(qiáng)度與氣動(dòng)控制鉸鏈力矩。其約束設(shè)置為

q=0.5ρv2≤qmax

(10)

在臨近空間范圍內(nèi)，攔截彈所能提供的可用過(guò)載大小有限，過(guò)載約束設(shè)置為

(11)

考慮到攔截彈的熱流密度約束，其簡(jiǎn)化計(jì)算公式為[16]

(12)

3 反臨攔截彈中末制導(dǎo)交接班分析

高精度的復(fù)合制導(dǎo)技術(shù)是攔截臨近空間高超聲速目標(biāo)的關(guān)鍵技術(shù)之一，復(fù)合制導(dǎo)包含了初、中、末制導(dǎo)以及各段間的交接班。其中，實(shí)現(xiàn)中末制導(dǎo)交接班較為困難，成為防御此類目標(biāo)的瓶頸問(wèn)題。

對(duì)于低速目標(biāo)來(lái)說(shuō)，交接班約束較為寬松，末制導(dǎo)的捕獲狀態(tài)空間較大，且攔截常規(guī)目標(biāo)時(shí)對(duì)脫靶量要求并不嚴(yán)苛，可認(rèn)為中末交接班瞬間完成。對(duì)高超聲速目標(biāo)來(lái)說(shuō)，隨著目標(biāo)速度增大，末制導(dǎo)捕獲狀態(tài)空間會(huì)急劇減小，攔截彈要想以很小的脫靶量擊中目標(biāo)，就必須快速滿足中末制導(dǎo)交接班的各類復(fù)雜約束，主要包括導(dǎo)引頭截獲條件約束和末制導(dǎo)捕獲條件約束。導(dǎo)引頭截獲[17]包括距離截獲、角度截獲和速度截獲，保證攔截彈能“看得到目標(biāo)”，此約束條件形式簡(jiǎn)單，容易描述。末制導(dǎo)捕獲條件約束主要由彈目相對(duì)距離矢量、攔截彈及目標(biāo)速度矢量組成，保證攔截彈能“打得上目標(biāo)”，對(duì)此約束的描述較為困難。本節(jié)側(cè)重研究末制導(dǎo)捕獲條件以及它對(duì)中末交接班的影響，最終目的是為中制導(dǎo)設(shè)計(jì)彈道終端約束條件提供參考。

3.1 彈目相對(duì)運(yùn)動(dòng)模型

中制導(dǎo)的終端時(shí)刻對(duì)應(yīng)的是末制導(dǎo)的初始時(shí)刻，此時(shí)的運(yùn)動(dòng)狀態(tài)主要包括攔截彈所在位置和其速度方向，這些狀態(tài)將決定在末制導(dǎo)的初始時(shí)刻攔截彈是否處在有利的攔截陣位。下面對(duì)彈目相對(duì)運(yùn)動(dòng)模型進(jìn)行討論。

如圖2所示，彈目相對(duì)位置矢量定義：

r=rT-rM=ρer

(13)

式中rM、rT分別為攔截彈和目標(biāo)的位置矢量；er為沿彈目視線方向的單位矢量；ρ為彈目的相對(duì)距離。

彈目的相對(duì)運(yùn)動(dòng)速度與加速度可表示如下：

(14)

式中vM、vT、aM、aT分別為攔截彈與目標(biāo)和速度與加速度。

彈目相對(duì)運(yùn)動(dòng)在三個(gè)方向分量上應(yīng)滿足：

(15)

彈目相對(duì)運(yùn)動(dòng)關(guān)系應(yīng)滿足：

(16)

式中aMt、aMr、aTt、aTr分別為攔截彈與目標(biāo)的加速度沿er與et方向上的分量。

圖2 彈目相對(duì)運(yùn)動(dòng)關(guān)系Fig.2 Relative movement of missile and target

3.2 基于PPN的捕獲條件推導(dǎo)

純比例導(dǎo)引律(PPN)是較為常用的一種比例導(dǎo)引改進(jìn)制導(dǎo)律，其被廣泛地應(yīng)用在各種高精度的末制導(dǎo)律設(shè)計(jì)中?；赑PN推導(dǎo)中末制導(dǎo)交接班的約束條件，對(duì)于研究中制導(dǎo)律的終端彈道約束具有十分重要的意義。下面結(jié)合文獻(xiàn)[18]，給出適用于臨近空間攔截彈的捕獲高速非機(jī)動(dòng)目標(biāo)的約束條件。

為便于描述攔截彈和目標(biāo)在末制導(dǎo)段的初始運(yùn)動(dòng)狀態(tài)對(duì)攔截的影響效果，做如下的定義：

(17)

(18)

式中 下標(biāo)“0”為在末制導(dǎo)捕獲時(shí)刻的狀態(tài)；θvM0、θvT0分別為攔截彈速度指向與彈目視線、目標(biāo)的速度指向與彈目視線的夾角。

通過(guò)分析彈目的幾何關(guān)系可知，確定了θvM0、θvT0，結(jié)合目標(biāo)的位置信息，就可確定攔截彈位置和速度指向(即彈道傾角)的初始狀態(tài)。因此，研究捕獲條件就是確定這兩個(gè)變量的所要滿足的約束條件。

定義攔截彈和目標(biāo)的速度大小的比值：

p=VM/VT

相應(yīng)地定義有效速度大小的比值：

基于PPN導(dǎo)引飛行，攔截彈捕獲高速非機(jī)動(dòng)目標(biāo)(p<1)的條件：

(20)

其中

(21)

式中β為導(dǎo)航比；θc為在能夠?qū)崿F(xiàn)攔截的條件下，目標(biāo)的初始速度方向與彈目視線方向的最大夾角。

具體計(jì)算過(guò)程：通過(guò)目標(biāo)軌跡預(yù)測(cè)獲得目標(biāo)速度指向與目標(biāo)位置信息，以目標(biāo)為參考，即可獲得θvT0的變化范圍，根據(jù)導(dǎo)引頭探測(cè)距離確定可行的攔截器位置變化范圍，再由θvM0和彈目視線角范圍可獲得攔截器的速度指向變化范圍。上述捕獲區(qū)描述的是攔截彈捕獲目標(biāo)時(shí)彈目交會(huì)角所要滿足的充分必要性條件。

(21)

4 仿真算例

4.1 PPN的交接班捕獲窗口定量描述

假定通過(guò)目標(biāo)軌跡預(yù)測(cè)獲得的目標(biāo)速度為2250 m/s；攔截彈的平均速度為1500 m/s；通過(guò)計(jì)算可得目標(biāo)與攔截彈速度比近似為p=0.67，選擇β=4。采用第2.4節(jié)中推導(dǎo)的公式計(jì)算捕獲區(qū)，可獲得PPN的交接班捕獲窗口，數(shù)據(jù)如表2所示。

上述PPN的交接班捕獲窗口數(shù)據(jù)可轉(zhuǎn)化為如圖3、圖4所示。

表2 中末制導(dǎo)交接班約束條件的特征點(diǎn)

圖3 捕獲窗口示意圖Fig.3 Schematic diagram of capture window

圖4 可攔截目標(biāo)所在區(qū)域Fig.4 Interception area of target

由以上可知，在采用PPN攔截高速目標(biāo)過(guò)程中，當(dāng)目標(biāo)速度指向與視線夾角θvT0為138.9°～221.8°時(shí),可有效攔截目標(biāo)，在端點(diǎn)處只有一個(gè)θvM0與之對(duì)應(yīng)；當(dāng)θvT0=176.8°時(shí)，攔截彈可選擇的初始速度指向范圍最大，攔截彈的初始速度的指向范圍越大；當(dāng)θvM0=-80.7°時(shí)，可攔截目標(biāo)的初始速度指向范圍最大。圖3以目標(biāo)為參考點(diǎn)，描述了末制導(dǎo)初始時(shí)刻攔截彈位置與速度指向的約束集合。圖4是以攔截彈為參考點(diǎn)，描述了末制導(dǎo)初始時(shí)刻目標(biāo)位置與攔截彈速度指向的約束集合。

由上述對(duì)交接班問(wèn)題的分析可知，給定目標(biāo)和攔截彈的速度比以及導(dǎo)航比，從幾何關(guān)系上可對(duì)末制導(dǎo)捕獲條件進(jìn)行定量描述，得到實(shí)現(xiàn)交接班時(shí)彈目交會(huì)角所要滿足的范圍，為設(shè)計(jì)中制導(dǎo)終端約束條件提供依據(jù)。由于本節(jié)得到的捕獲區(qū)窗口是中制導(dǎo)終端所要滿足的一個(gè)必要條件，即角度滿足約束才能順利實(shí)現(xiàn)中末制導(dǎo)交接班。但具體的終端約束條件的設(shè)置還與攔截場(chǎng)作戰(zhàn)景有關(guān)，根據(jù)高超聲速目標(biāo)的軌跡預(yù)測(cè)理論研究[19]，預(yù)測(cè)得到的目標(biāo)軌跡實(shí)際上是一條管道，目標(biāo)在某一時(shí)刻的預(yù)測(cè)位置實(shí)際上是一個(gè)范圍，即目標(biāo)在此時(shí)刻可能出現(xiàn)的區(qū)域。因此為了確保對(duì)目標(biāo)的有效攔截，應(yīng)采用多彈協(xié)同攔截的方式，借助多枚攔截彈對(duì)預(yù)測(cè)區(qū)域范圍進(jìn)行覆蓋，并根據(jù)合理的任務(wù)分配計(jì)算單枚攔截彈預(yù)測(cè)命中點(diǎn)，進(jìn)而計(jì)算預(yù)測(cè)交班點(diǎn)信息，確定各個(gè)攔截彈的終端約束條件。由于本文不涉及多彈協(xié)同攔截的任務(wù)分配、預(yù)測(cè)命中點(diǎn)及預(yù)測(cè)交班點(diǎn)信息的計(jì)算等內(nèi)容，所以4.2節(jié)的仿真中直接給出預(yù)測(cè)交班條件。

4.2 彈道與導(dǎo)引律分離設(shè)計(jì)驗(yàn)證

垂直上升段的原則是保證攔截彈達(dá)到一定速度，以便開(kāi)始控制。假定垂直上升段末端速度為500 m/s，高度達(dá)5000 m，轉(zhuǎn)彎段以此狀態(tài)為起點(diǎn)。由于本節(jié)驗(yàn)證分離設(shè)計(jì)的優(yōu)化彈道的優(yōu)勢(shì)，所以這里直接給定中制導(dǎo)的終端位置約束為(250 000，30 000)，終端彈道傾角約束為-30°，采用高斯偽譜法優(yōu)化[20]，仿真結(jié)果如圖5所示。

圖5 高斯偽譜法優(yōu)化彈道Fig.5 Optimal trajectory by Gauss pseudospectral method

為說(shuō)明彈道與導(dǎo)引律分離設(shè)計(jì)的優(yōu)勢(shì)，將基于高斯偽譜法的彈道與采用比例導(dǎo)引的彈道進(jìn)行了比較，三種彈道的初始條件一致，仿真結(jié)果如圖6、圖7和表3所示。其中Gauss表示根據(jù)路徑約束與預(yù)測(cè)交班信息，采用高斯偽譜法優(yōu)化的彈道。PN-Target指基于目標(biāo)運(yùn)動(dòng)參數(shù)，按照比例導(dǎo)引飛行，目標(biāo)初始位置距攔截彈發(fā)射點(diǎn)615 km，以2250 m/s的速度在30 km高度上平飛，當(dāng)彈目相距100 km時(shí)認(rèn)為進(jìn)入交接班過(guò)程，計(jì)算停止；PN-HP指將預(yù)測(cè)交班點(diǎn)作為目標(biāo)位置，認(rèn)為目標(biāo)固定，攔截彈按照比例導(dǎo)引飛行。

圖6 彈道對(duì)比Fig.6 Trajectory comparison

圖7 彈道過(guò)程量對(duì)比Fig.7 Comparison of trajectory process variable

彈道t/sθf(wàn)/(°)vf/(m/s)Gauss128．5-302179．5PN?Target159．7-41．762114．8PN?HP169．3-0．261887．8

由圖6、圖7和表3可知，按照本文助推方案設(shè)計(jì)，三種彈道都為先爬升而后下降的過(guò)程。PN-Target彈道在中末制導(dǎo)交接班時(shí)相對(duì)目標(biāo)的位置和速度指向，都是由制導(dǎo)律本身決定的，無(wú)法按照前述終端約束條件的要求，在交接班時(shí)刻實(shí)現(xiàn)對(duì)攔截彈位置和彈道傾角的控制，也就無(wú)法保證進(jìn)入末制導(dǎo)時(shí)，攔截彈處于較為有利的陣位，且其飛行高度較高，導(dǎo)致飛行時(shí)間較長(zhǎng)。若對(duì)中制導(dǎo)終端位置進(jìn)行控制，即PN-HP彈道，其飛行時(shí)間增加，末速降低，不利于攔截作戰(zhàn)。由彈道過(guò)程量對(duì)比可知，三種彈道動(dòng)壓、過(guò)載區(qū)別不大，基于高斯偽譜法的彈道的熱流密度相對(duì)于比例導(dǎo)引彈道較大，主要由于其在臨近空間下層以相對(duì)高速飛行造成。從控制量的角度看，基于高斯偽譜法的彈道全程攻角較小，而比例導(dǎo)引彈道全程攻角接近最大值，導(dǎo)致較大能量損耗，且存在大范圍跳變，控制上難以實(shí)現(xiàn)。

5 結(jié)論

(1)本文結(jié)合目標(biāo)特性分析，基于彈道與導(dǎo)引律分離原則設(shè)計(jì)了高拋再入彈道方案，建立了彈道優(yōu)化模型，研究了中末制導(dǎo)交接班的條件，從而為完善攔截彈終端約束模型提供了參考。仿真結(jié)果表明，本文設(shè)計(jì)的彈道能夠?qū)K端位置和彈道傾角進(jìn)行控制，且飛行時(shí)間短、末速較大，是較為理想的攔截彈道。

(2)本文只考慮了單枚彈的彈道設(shè)計(jì)問(wèn)題，為確保對(duì)目標(biāo)的有效攔截，應(yīng)采用多彈協(xié)同攔截的方式，后續(xù)研究應(yīng)考慮多枚攔截彈對(duì)預(yù)測(cè)區(qū)域范圍進(jìn)行覆蓋，并根據(jù)合理的任務(wù)分配計(jì)算單枚攔截彈預(yù)測(cè)命中點(diǎn)，進(jìn)而計(jì)算預(yù)測(cè)交班點(diǎn)信息，確定各個(gè)攔截彈的終端約束條件。

參考文獻(xiàn)：

[1] Bertin J J,Cummings R M.Fifty years of hypersonics:where we’ve been,where we’re going[J].Progress in Aerospace Sciences,2003,39:511-536.

[2] 馮志高,關(guān)成啟,張紅文.高超聲速飛行器概論[M].北京:北京理工大學(xué)出版社,2016.

FENG Zhigao,GUAN Chengqi,ZHANG Hongwen.An introduction to hypersonic aircraft[M].Beijing:Beijing Institute of Technology Press,2016.

[3] 孫磊,黃可西,常曉飛,等.臨近空間高超聲速巡航導(dǎo)彈攔截問(wèn)題研究[J].西北工業(yè)大學(xué)學(xué)報(bào),2015(4):615-620.

SUN Lei,HUANG Kexi,CHANG Xiaofei,et al.Researching Interception of Near Space Hypersonic Cruise Missile[J].Journal of Northwestern Polytechnical University,2015(4):615-620.

[4] 金欣,梁維泰,王俊,等.反臨近空間目標(biāo)作戰(zhàn)的若干問(wèn)題思考[J].現(xiàn)代防御技術(shù),2013,41(6):1-7.

JIN Xin,LIANG Weitai,WANG Jun,et al.Discussion on Issues of Anti-Near-Space-Target Operation[J].Modern defense technology,2013,41(6):1-7.

[5] 魏倩,蔡遠(yuǎn)利.J2項(xiàng)攝動(dòng)影響下的大氣層外彈道規(guī)劃改進(jìn)算法[J].控制理論與應(yīng)用,2016,33(9):1245-1251.

WEI Qian,CAI Yuanli.Trajectory planning algorithm of exo-atmosphere aircraft under the influence of the J2 perturbation[J].Control Theory & Applications,2016,33(9):1245-1251.

[6] 魏倩,蔡遠(yuǎn)利.一種基于神經(jīng)網(wǎng)絡(luò)的中制導(dǎo)改進(jìn)算法[J].西安交通大學(xué)學(xué)報(bào),2016,50(7):125-130.

WEI Qian,CAI Yuanli.A Modified Algorithm on the Midcourse Guidance Based on BP Neural Network[J].Journal of Xi’an Jiaotong University,2016,50(7):125-130.

[7] Shaul G,Sergey R.Exo-atmospheric mid-course guidance[C]//AIAA Guidance,Navigation and Control Conference.New York,USA,2015:2015-0088.

[8] BT 斯維特洛夫.防空導(dǎo)彈設(shè)計(jì)[M].北京:宇航出版社,2004.

Светлов В Г.Air defense missile design[M].Beijing:Chinese Astronautics Press,2004.

[9] 李寧波,雷虎民,譚詩(shī)利,等.基于虛擬域動(dòng)態(tài)逆的空地反輻射導(dǎo)彈軌跡優(yōu)化[J].彈道學(xué)報(bào),2016,28(4):7-11.

LI Ningbo,LEI Humin,TAN Shili,et al.Trajectory optimization for air-to-ground ARM based on inverse dynamics in the virtual domain[J].Journal of Ballistics,2016,28(4):7-11.

[10] Lu P,Forbes S,Morrgan B.Gliding guidance of high L/D hypersonic vehicles[R].AIAA 2013-4648.

[11] Steven H W,Col J S,Dale S.The DARPA/AF falcon program:The hypersonic technology vehicle #2(HTV-2)flight demonstration Phase[R].AIAA 2008-2539.

[12] LI Ningbo,LEI Humin,ZHOU Jin,et al.Variable-time-domain online neighboring optimal trajectory modification for hypersonic interceptors[J/OL].International Journal of Aerospace Engineering,2017,doi:10.1155/2017/9456179.

[13] 趙杰,王君,張大元,等.反臨近空間高超聲速飛行器中末交接視角研究[J].飛行力學(xué),2015,33(3):253-256.

ZHAO Jie,WANG Jun,ZHANG Dayuan,et al.Visual angle research for midcourse and terminal guidance hand-over of near space hypersonic vehicles[J].Flight Dynamics,2015,33(3):253-256.

[14] 李新國(guó),方群.有翼導(dǎo)彈飛行動(dòng)力學(xué)[M].西安:西北工業(yè)大學(xué)出版社,2005:153-155.

LI Xinguo,FANG Qun.Flight dynamics of winged missile[M].Xi’an:Northwestern Polytechnical University Press,2005:49-50.

[15] Peter J Mantel.The missile defense equation:factors for decision making[M].AIAA,2004.

[16] Afshin R,Krishna D K,Hekmat A.Particle swarm optimization applied to spacecraft reentry trajectory[J].Journal of Guidance,Control,and Dynamics,2013 36(1):307-310.

[17] 樊會(huì)濤.復(fù)合制導(dǎo)空空導(dǎo)彈截獲目標(biāo)概率研究[J].航空學(xué)報(bào),2010,31(6):1225-1229.

FAN Huitao.Study on Target Acquisition Probability of Air-to-air Missiles with Combined Guidance[J].Acta Aeronautica et Astronautica Sinica,2010,31(6):1225-1229.

[18] Feng Tyan.Analysis of 3D PPN guidance laws for nonmaneuvering target[J].IEEE Transactions on Aerospace and Electronic Systems,2015,51(4):2923-2934.

[19] 翟岱亮,雷虎民,李炯,等.基于自適應(yīng)IMM的高超聲速飛行器軌跡預(yù)測(cè)[J].航空學(xué)報(bào),2016,37(11):3466-3475.

ZHAI Dailiang,LEI Humin,LI Jiong,et al.Trajectory prediction of hypersonic vehicle based on adaptive IMM[J].Acta Aeronautica et Astronautica Sinica,2016,37(11):3466-3475.

[20] LI Ningbo,LEI Humin,SHAO Lei,et al.Trajectory optimization based on multi-interval mesh refinement method[J/OL].Mathematical Problems in Engineering,2017.doi:10.1155/2017/8521368.