蔣明明，袁慶航，趙長(zhǎng)見(jiàn)，卜奎晨

（中國(guó)運(yùn)載火箭技術(shù)研究院，北京，100076）

助推-滑翔式導(dǎo)彈總體參數(shù)設(shè)計(jì)方法初探

蔣明明，袁慶航，趙長(zhǎng)見(jiàn)，卜奎晨

（中國(guó)運(yùn)載火箭技術(shù)研究院，北京，100076）

以兩級(jí)固體推進(jìn)劑火箭發(fā)動(dòng)機(jī)助推、彈頭直接入軌，而后全程在大氣層內(nèi)滑翔飛行的助推-滑翔式導(dǎo)彈為研究對(duì)象，對(duì)導(dǎo)彈總體參數(shù)設(shè)計(jì)進(jìn)行研究，給出一種適用于助推-滑翔式導(dǎo)彈的總體參數(shù)設(shè)計(jì)方法。根據(jù)助推-滑翔式導(dǎo)彈的彈道特點(diǎn)，通過(guò)分段分析彈道特性，推導(dǎo)出導(dǎo)彈總體參數(shù)與關(guān)機(jī)點(diǎn)理想速度間的關(guān)系式。通過(guò)仿真分析，建立滑翔起點(diǎn)參數(shù)與關(guān)機(jī)點(diǎn)參數(shù)間的關(guān)系模型；考慮平衡滑翔條件，得到滑翔射程公式?；谝陨瞎胶湍Ｐ?，給出助推-滑翔式導(dǎo)彈射程與關(guān)機(jī)點(diǎn)參數(shù)之間的解析關(guān)系，初步建立了助推-滑翔式導(dǎo)彈總體參數(shù)的設(shè)計(jì)方法。

助推-滑翔式導(dǎo)彈；武器裝備；總體參數(shù)設(shè)計(jì)

0 引 言

助推-滑翔式導(dǎo)彈是一種利用火箭助推，彈頭在大氣層內(nèi)靠氣動(dòng)升力進(jìn)行滑翔飛行的導(dǎo)彈武器。助推-滑翔式導(dǎo)彈結(jié)合了彈道導(dǎo)彈與巡航導(dǎo)彈的優(yōu)點(diǎn)，可以實(shí)現(xiàn)大氣層內(nèi)高超聲速飛行，實(shí)現(xiàn)遠(yuǎn)距離精確打擊，具有強(qiáng)機(jī)動(dòng)能力和突防能力，是目前世界主要軍事大國(guó)導(dǎo)彈武器研究的主要方向。助推-滑翔式彈道的概念早在1933年就已經(jīng)被提出，但由于科學(xué)技術(shù)發(fā)展水平的限制，基于這種彈道的高超聲速滑翔飛行器一直處于驗(yàn)證階段。近年來(lái)，以美國(guó)和俄羅斯為代表的主要軍事大國(guó)正加緊驗(yàn)證助推-滑翔式飛行器的關(guān)鍵技術(shù)，并列入武器裝備研究計(jì)劃。

導(dǎo)彈總體參數(shù)的設(shè)計(jì)是導(dǎo)彈設(shè)計(jì)的基礎(chǔ)，在方案論證階段快速獲得總體參數(shù)將有助于提高總體設(shè)計(jì)效率、縮短設(shè)計(jì)周期。文獻(xiàn)[1]推導(dǎo)出了彈道式飛行器的航程公式；文獻(xiàn)[2]推導(dǎo)出了單級(jí)彈道式飛行器、高超聲速助推-滑翔式飛行器及高超聲速助推-巡航飛行器的起飛質(zhì)量與射程、推進(jìn)劑比沖和載荷之間的解析關(guān)系；文獻(xiàn)[3]通過(guò)引入射程經(jīng)驗(yàn)公式并進(jìn)行近似質(zhì)量分析，將求解固體彈道導(dǎo)彈最優(yōu)級(jí)間比轉(zhuǎn)化為求解各子級(jí)推進(jìn)劑質(zhì)量比之間的關(guān)系，以最小起飛質(zhì)量為優(yōu)化目標(biāo)，求得了n級(jí)導(dǎo)彈各級(jí)最佳推進(jìn)劑質(zhì)量比；文獻(xiàn)[4]以導(dǎo)彈起飛質(zhì)量最小為目標(biāo)，對(duì)導(dǎo)彈質(zhì)量進(jìn)行近似分析，并忽略速度的阻力損失以及噴管出口壓力損失，將求解n級(jí)液體彈道導(dǎo)彈最優(yōu)總體設(shè)計(jì)參數(shù)轉(zhuǎn)化為求解各級(jí)最優(yōu)推進(jìn)劑質(zhì)量比和推重比，獲得了可行的理論解法。

目前，助推-滑翔式導(dǎo)彈的總體參數(shù)設(shè)計(jì)尚未有系統(tǒng)的理論方法。本文根據(jù)助推-滑翔式導(dǎo)彈的彈道特點(diǎn)，對(duì)其彈道進(jìn)行分段研究，得到導(dǎo)彈射程與關(guān)機(jī)點(diǎn)參數(shù)之間的解析關(guān)系，初步建立了助推-滑翔式導(dǎo)彈總體參數(shù)的設(shè)計(jì)方法，并通過(guò)設(shè)計(jì)實(shí)例驗(yàn)證此方法的可行性，為此類型導(dǎo)彈的總體參數(shù)設(shè)計(jì)提供了參考。

1 助推-滑翔式導(dǎo)彈總體參數(shù)設(shè)計(jì)方法

根據(jù)助推-滑翔式導(dǎo)彈的彈道特點(diǎn)，將其彈道分為助推段、無(wú)動(dòng)力下壓段、拉起段、滑翔段及下壓攻擊段，如圖1所示。

在助推-滑翔式導(dǎo)彈的彈道中，滑翔段射程占導(dǎo)彈射程的絕大部分，助推段及下壓攻擊段射程占導(dǎo)彈射程比例較小，解析分析時(shí)可以忽略這2段對(duì)導(dǎo)彈全射程的影響。

1.1 導(dǎo)彈總體參數(shù)與關(guān)機(jī)點(diǎn)理想速度的解析關(guān)系

雖然助推段射程對(duì)于導(dǎo)彈全射程的影響很小，但其彈道特性直接決定了關(guān)機(jī)點(diǎn)參數(shù)，對(duì)導(dǎo)彈的射程及飛行控制等均有較大影響。因此，本節(jié)以兩級(jí)固體導(dǎo)彈為例，研究導(dǎo)彈起飛質(zhì)量等總體設(shè)計(jì)參數(shù)與關(guān)機(jī)點(diǎn)理想速度間的關(guān)系，為建立助推-滑翔式導(dǎo)彈總體參數(shù)設(shè)計(jì)方法打下基礎(chǔ)。

引入發(fā)動(dòng)機(jī)質(zhì)量比σ1和σ2，導(dǎo)彈各級(jí)起飛質(zhì)量可表示為

式中 Myx為有效載荷即彈頭質(zhì)量；M01和M02分別為導(dǎo)彈一、二級(jí)起飛質(zhì)量；Mp1和Mp2分別為導(dǎo)彈一、二級(jí)推進(jìn)劑質(zhì)量；Mj1和Mj2分別為導(dǎo)彈一、二子級(jí)結(jié)構(gòu)質(zhì)量。

根據(jù)齊奧爾科夫斯基公式，兩級(jí)導(dǎo)彈二級(jí)發(fā)動(dòng)機(jī)關(guān)機(jī)時(shí)刻的理想速度Vk1為

式中 Ispv1和Ispv2分別為導(dǎo)彈一、二級(jí)真空比沖。

引入有效推進(jìn)劑質(zhì)量比μki= Mpi/M0i，則有：

這是一個(gè)求函數(shù)條件極值的問(wèn)題：Vk1給定的情況下，求M01為最小時(shí)的參數(shù)μk1和μk2。即為求下述拉格朗日函數(shù)無(wú)條件極值的問(wèn)題：

式中 Λ為待定拉格朗日乘子；φ 為約束條件：

令Myx+Mj2= a，Mj1= b，Ispv2/Ispv1= c，經(jīng)計(jì)算得：

在得到Vk1之后，可由式（7）求得二級(jí)有效推進(jìn)劑質(zhì)量比μk2，再將μk2帶入式（8），即可求得一級(jí)有效推進(jìn)劑質(zhì)量比μk1。之后，即可計(jì)算Mp1和Mp2等參數(shù)，從而確定出導(dǎo)彈最優(yōu)起飛質(zhì)量。

通過(guò)計(jì)算驗(yàn)證公式推導(dǎo)的正確性計(jì)算條件如下：

a）Vk1= 4000 m/s，有效載荷（即彈頭）質(zhì)量為1600 kg；

b）Ispv1= 2600 (N·s)/kg，Ispv2= 2800 (N·s)/kg；

c）Mj1= 500 kg，Mj2= 200 kg。

導(dǎo)彈起飛質(zhì)量與二級(jí)有效推進(jìn)劑質(zhì)量比、發(fā)動(dòng)機(jī)質(zhì)量比的關(guān)系曲線如圖2所示。

圖2 導(dǎo)彈起飛質(zhì)量隨有效推進(jìn)劑質(zhì)量比變化曲線

從圖2中可以看出，對(duì)應(yīng)不同發(fā)動(dòng)機(jī)質(zhì)量比，均有一個(gè)最優(yōu)二級(jí)有效推進(jìn)劑質(zhì)量比使得導(dǎo)彈起飛質(zhì)量最小，且其值與理論值相同，證明了導(dǎo)彈總體參數(shù)設(shè)計(jì)方法的正確性。

1.2 滑翔起點(diǎn)參數(shù)與關(guān)機(jī)點(diǎn)參數(shù)的關(guān)系

滑翔起點(diǎn)參數(shù)包括滑翔起點(diǎn)的速度和射程（即無(wú)動(dòng)力下壓段及拉起段射程之和），下面分別求解這2個(gè)參數(shù)與關(guān)機(jī)點(diǎn)速度Vk和關(guān)機(jī)點(diǎn)當(dāng)?shù)貜椀纼A角θk之間的關(guān)系。

通過(guò)仿真分析，滑翔起點(diǎn)速度值V0與關(guān)機(jī)點(diǎn)參數(shù)近似滿足以下關(guān)系：

式中 a2和b2為擬合為kθ的多項(xiàng)式函數(shù)。

通過(guò)對(duì)仿真結(jié)果進(jìn)行擬合，無(wú)動(dòng)力下壓段及拉起段的射程之和L1與關(guān)機(jī)點(diǎn)參數(shù)間近似滿足以下關(guān)系：

式中 a1，b1，c可以擬合為Vk的多項(xiàng)式函數(shù)；Lke為以關(guān)機(jī)點(diǎn)參數(shù)為起點(diǎn)計(jì)算出的真空狀態(tài)下的拋物線射程，計(jì)算公式[5]如下：

式中 βe為射程角；rk為關(guān)機(jī)點(diǎn)地心距；νk為關(guān)機(jī)點(diǎn)能量參數(shù)，；μ 為地球引力系數(shù)。

1.3 滑翔段射程公式

在進(jìn)行滑翔段射程分析時(shí)，為了簡(jiǎn)化問(wèn)題，對(duì)導(dǎo)彈滑翔運(yùn)動(dòng)模型作如下假設(shè)：

a）不考慮地球旋轉(zhuǎn)，即ωe= 0；

b）地球?yàn)橐痪|(zhì)圓球，即引力加速度與地心距平方成正比，g = fM/r2；

c）導(dǎo)彈僅在平面內(nèi)運(yùn)動(dòng)，無(wú)滾轉(zhuǎn)和側(cè)滑。導(dǎo)彈滑翔段的運(yùn)動(dòng)方程為

式中 V，θ，h，φ 分別為導(dǎo)彈速度、當(dāng)?shù)貜椀纼A角、飛行高度及射程角；R0，ρ，S，m分別為球形地球半徑、大氣密度、導(dǎo)彈參考面積和質(zhì)量；CL和CD分別為導(dǎo)彈的升力系數(shù)和阻力系數(shù)。

本文中，滑翔段采用平衡滑翔飛行，平衡滑翔指飛行器在大氣層內(nèi)無(wú)動(dòng)力飛行時(shí)，沿航跡每一點(diǎn)處飛行器所受的升力、重力和離心力平衡，即0θ=˙的狀態(tài)。對(duì)式（13）進(jìn)行整理可得平衡滑翔狀態(tài)受力關(guān)系為

高超聲速飛行器滑翔距離一般在上千公里，而滑翔段高程變化一般小于50 km，因此當(dāng)?shù)貜椀纼A角和徑向速度都很小，可以認(rèn)為θ = 0，代入式（16）中可得：

由式（15）除以式（12），并結(jié)合平衡滑翔條件，可得平衡滑翔條件下射程角對(duì)滑翔速度的偏導(dǎo)數(shù)：

將式（12）代入式（13），并令K = CL/CD，整理可得：

假設(shè)平衡滑翔過(guò)程中，升阻比K和高度h近似不變，設(shè)V0和h0分別為滑翔起始點(diǎn)的速度和高度，對(duì)式（19）進(jìn)行積分，可得平衡滑翔射程角φ 與速度V的關(guān)系式：

相應(yīng)的平衡滑翔射程為

由式（21）建立滑翔距離與升阻比、滑翔高度及滑翔速度的關(guān)系。其中，滑翔高度及升阻比均近似為常值，滑翔結(jié)束點(diǎn)速度可根據(jù)落地點(diǎn)速度約束選定，即平衡滑翔射程僅與滑翔起點(diǎn)速度值有關(guān)。

1.4 導(dǎo)彈射程與關(guān)機(jī)點(diǎn)參數(shù)之間的關(guān)系

將式（10）和式（21）相加，即可得到助推-滑翔式導(dǎo)彈全射程公式：

將式（9）及式（10）代入式（22），即建立了助推-滑翔式導(dǎo)彈射程與關(guān)機(jī)點(diǎn)速度Vk和關(guān)機(jī)點(diǎn)當(dāng)?shù)貜椀纼A角θk間的關(guān)系。

1.5 導(dǎo)彈總體參數(shù)與射程指標(biāo)之間的關(guān)系

給定射程指標(biāo)L及式（22）中的各常量后，即可求得助推-滑翔式導(dǎo)彈的關(guān)機(jī)點(diǎn)參數(shù)。之后，即可根據(jù)導(dǎo)彈總體參數(shù)與關(guān)機(jī)點(diǎn)理想速度間的解析關(guān)系，通過(guò)關(guān)機(jī)點(diǎn)參數(shù)求得助推-滑翔式導(dǎo)彈的起飛規(guī)模等總體設(shè)計(jì)參數(shù)。

2 仿真算例

技術(shù)指標(biāo)要求：射程1300 km，彈頭質(zhì)量1600 kg。

一般，助推-滑翔式導(dǎo)彈關(guān)機(jī)點(diǎn)處當(dāng)?shù)貜椀纼A角θk較小，可取θk= 7°，滑翔升阻比K = 2.5，滑翔高度h = 40 km，滑翔結(jié)束點(diǎn)速度V = 1000 m/s；導(dǎo)彈一、二子級(jí)結(jié)構(gòu)質(zhì)量分別為500 kg和200 kg；一、二級(jí)發(fā)動(dòng)機(jī)質(zhì)量比σ1= σ2= 0.9；一、二級(jí)發(fā)動(dòng)機(jī)真空比沖分別為Ispv1= 2600 (N·s)/kg，Ispv2= 2800 (N·s)/kg。

根據(jù)工程經(jīng)驗(yàn)，經(jīng)大量仿真計(jì)算，將a2和b2擬合為θk的五次多項(xiàng)式，a1，b1，c擬合為Vk的三次多項(xiàng)式函數(shù)，所得擬合結(jié)果較為理想。

由式（10），得到導(dǎo)彈關(guān)機(jī)點(diǎn)速度：Vk= 3017 m/s。

考慮到實(shí)際飛行中存在的速度損失，同時(shí)為簡(jiǎn)化問(wèn)題，本文中取速度損失ΔVk= 0.25Vk，則理想速度：Vk1= Vk+ΔVk= 1.25Vk= 3771 m/s。

根據(jù)助推-滑翔式導(dǎo)彈總體參數(shù)設(shè)計(jì)方法，計(jì)算得到助推-滑翔式導(dǎo)彈總體參數(shù)為μk1= 0.32；μk2= 0.63；M01= 10 017 kg；M02= 5 964 kg。

為驗(yàn)證助推-滑翔式導(dǎo)彈總體參數(shù)設(shè)計(jì)方法的可行性，基于上述總體參數(shù)進(jìn)行助推-滑翔式彈道的數(shù)學(xué)仿真驗(yàn)證。取一級(jí)真空推力P1= 280 kN，二級(jí)真空推力P2= 200 kN，并進(jìn)行彈道仿真，仿真結(jié)果如圖3所示。

仿真所得導(dǎo)彈落速為560 m/s，當(dāng)?shù)貜椀纼A角為-80°，導(dǎo)彈射程為1250 km，與設(shè)計(jì)指標(biāo)間的誤差為3.8%，由此可見(jiàn)，本文得到的助推-滑翔式導(dǎo)彈總體參數(shù)設(shè)計(jì)方法可以用于方案論證初期的導(dǎo)彈總體參數(shù)選擇。

Preliminary Investigation of Parameter Design Method for Boost-Glide Missile

Jiang Ming-ming, Yuan Qing-hang, Zhao Chang-jian, Bu Kui-chen
(China Academy of Launch Vehicle Technology, Beijing, 100076)

Based on two class solid propellant missile whose warhead orbits directly and then glides in the aerosphere all along, missile parameter design is studied, and a parameter design method for boost-glide missile is proposed. According to the characteristic of boost-glide trajectory, the trajectory is divided into several parts and each part is studied: the relationship between missile parameters and the theoretical velocity is deducted; through simulation, the model used to describe the relationship between parameters of glide starting point and that of engine shutdown point is established; considering the balanced glide condition, the formula to estimate the glide range is obtained. Based on the models, the analytical expression between the missile range and the parameters of the engine shutdown point is given, and the parameter design method for glide missile is established preliminarily.

Boost-glide missile; Weaponry and equipment; Parameter design

V421.1

A

1004-7182(2016)04-0005-04

10.7654/j.issn.1004-7182.20160402

2015-01-06；

2015-02-13

蔣明明（1990-），男，助理工程師，主要研究方向?yàn)轱w行器設(shè)計(jì)