黃浩然 王宏偉 常 旭

（中國(guó)人民解放軍空軍航空大學(xué)飛行器與動(dòng)力系，吉林 長(zhǎng)春130022）

0 引言

對(duì)飛機(jī)雷達(dá)散射源的分析可知，飛機(jī)是一個(gè)很復(fù)雜的目標(biāo)，由多個(gè)部件組成。在雷達(dá)波照射下，每個(gè)部件都會(huì)產(chǎn)生散射波，有的部件可能同時(shí)產(chǎn)生散射機(jī)理不同的散射源，形成多個(gè)散射源[1]。由于飛機(jī)整機(jī)RCS是各散射源綜合的結(jié)果，不難理解，降低了每個(gè)散射源的RCS值就有可能使整機(jī)的RCS值降低。為此，在飛行器的設(shè)計(jì)中，應(yīng)設(shè)法降低有關(guān)的部件的RCS值。而飛機(jī)的垂尾作為飛機(jī)上的重要散射源之一，降低其RCS對(duì)提高飛機(jī)的隱身性能起著重要作用[2]。因此，需要研究垂尾的雷達(dá)散射特性，并尋找方法來(lái)降低其RCS。

本文首先建立垂尾的三維模型，利用矩量法(MoM)計(jì)算模型RCS，重點(diǎn)研究垂尾前緣后掠角、展長(zhǎng)及傾角對(duì)RCS的影響，并擬合垂尾RCS隨這三個(gè)變量的變化曲線，建立了垂尾RCS與這三個(gè)變量之間的函數(shù)式。并由此提出減縮垂尾雷達(dá)散射截面積RCS的方法。

1 MoM方法原理

積分方程一般采用矩量法(method of moment,MoM)求解。矩量法是由R.F.Harrington于1968年提出的一種嚴(yán)格數(shù)值方法，具有較高的求解精度。其數(shù)學(xué)本質(zhì)是一種求解線性方程組的方法。

矩量法的基本原理是用許多離散的子域來(lái)代表整個(gè)連續(xù)區(qū)域。在每個(gè)子域中，未知函數(shù)用帶有未知系數(shù)的基函數(shù)來(lái)表示，因此無(wú)限個(gè)自由度的問(wèn)題就被轉(zhuǎn)化成了有限個(gè)自由度的問(wèn)題，然后運(yùn)用點(diǎn)匹配法、線匹配法或伽略金法(Galerkin)等方法進(jìn)行檢驗(yàn)，得到一組代數(shù)方程(即矩陣方程)，最后通過(guò)求解這一矩陣方程獲得數(shù)值解。

矩量法求解積分方程主要包括以下幾個(gè)步驟：

（1）區(qū)域的離散化或目標(biāo)的劃分；

（2）選擇合適的基函數(shù)和檢驗(yàn)函數(shù)；

（3）填充阻抗矩陣；

（4）求解矩陣方程。

例如：算子方程為

L{f（x）}=g（x）（1）

式中：L——線性算子；

f（x）——待求的未知函數(shù)；

g（x）——已知函數(shù)。

為了求解該算子方程，需要先將未知函數(shù)f（x）展開(kāi)為一系列已知函數(shù)(基函數(shù))疊加的形式，即

式中：ɑn——第n個(gè)基函數(shù)的待求系數(shù)；

bn——第n個(gè)基函數(shù)。

精確解通常需要無(wú)窮項(xiàng)的求和，對(duì)于近似解，式（2）為有限項(xiàng)求和，項(xiàng)數(shù)N(即未知量數(shù)目)由計(jì)算精度確定?；瘮?shù)之間應(yīng)該是線性無(wú)關(guān)的，所以必須選擇合適的基函數(shù)，使其能精確地模擬電流分布。式（1）的殘差可以表示為

求解該算子方程的目標(biāo)之一就是使得殘差足夠小，以保證求解的精度。因此需要選擇合適的權(quán)函數(shù)tm（x），與基函數(shù)做內(nèi)積。內(nèi)積被定義為

式（4）可以寫為線性方程組的形式，即

式中：Zmn——阻抗元素；

vm——激勵(lì)源。

通過(guò)直接法(如LU分解，即三角分解)或迭代法(如Krylov子空間迭代方法)求解式（5），即可求得所有基函數(shù)的未知系數(shù)ɑn，從而根據(jù)式（1）得到待求的未知函數(shù)。

2 垂尾散射特性分析

影響垂尾雷達(dá)散射截面積的幾何參數(shù)主要有：垂尾前緣后掠角χ、垂尾展長(zhǎng)bc和垂尾傾角θ。在這里，我們利用Feko軟件分別對(duì)以上幾何參數(shù)不同的垂尾模型進(jìn)行電磁仿真，求出其RCS。

2.1 前緣后掠角對(duì)垂尾RCS的影響

首先建立垂尾的三維模型，模型具體參數(shù)為：翼型NACA0009，根部弦長(zhǎng)cc=3.1m，展長(zhǎng)bc=3.7m，后緣后掠角20°，前緣后掠角從30°到47.5°。利用矩量法(MoM)計(jì)算的RCS曲線如圖1所示。

將圖1中各垂尾模型-20°～20°俯仰角范圍內(nèi)的RCS平均值(圖2中的星點(diǎn)）利用最小二乘法進(jìn)行擬合，得到曲線，如圖2所示。

由圖2可以看出，在俯仰角±20°范圍內(nèi)，隨著前緣后掠角的增大，垂尾RCS減小。根據(jù)計(jì)算結(jié)果，可以將垂尾RCS與前緣后掠角的關(guān)系擬合成圖中的曲線，表達(dá)式如下：

S=-0.0023χ2+0.0114χ-26.3780 （6）

通過(guò)式（6）可知，隨著前緣后掠角的增大，垂尾RCS均值減小加快。這表明，當(dāng)進(jìn)一步增大前緣后掠角時(shí)，RCS減縮收益會(huì)更大，需要指出，這可能會(huì)付出較大的氣動(dòng)性能代價(jià)。

2.2 展長(zhǎng)對(duì)垂尾RCS的影響

保持垂尾模型的翼型和根部弦長(zhǎng)不變，前緣后掠角和后緣后掠角分別固定為40°和20°。展長(zhǎng)從3.4m增大到4.1m，計(jì)算垂尾模型的RCS曲線，如圖3所示。并將垂尾模型在俯仰角范圍為-20°～20°內(nèi)的RCS平均值繪制成圖4。

從圖4中可以看出，垂尾RCS的均值（圖中的星點(diǎn)）隨展長(zhǎng)的分布是不穩(wěn)定的。將RCS值與展長(zhǎng)間的關(guān)系近似擬合為圖4中的二次曲線，可得出垂尾RCS值隨展長(zhǎng)的增大而增大，其曲線表達(dá)式如下：

S=0.8757b2c-7.1768bc-15.8524 （7）

通過(guò)式（7）可知，對(duì)所研究翼型，展長(zhǎng)小于3.7時(shí)，RCS均值隨展長(zhǎng)均似線性減?。淮笥?.7時(shí)，兩者呈非線性關(guān)系，RCS均值減小減慢。這表明，當(dāng)展長(zhǎng)大于某個(gè)量值后，進(jìn)一步增大展長(zhǎng)，RCS減縮收益減小，于是對(duì)于展長(zhǎng)的取值，應(yīng)權(quán)衡氣動(dòng)、強(qiáng)度等其他要求。

2.3 傾角對(duì)垂尾RCS的影響

保持垂尾模型的翼型和根部弦長(zhǎng)不變，前緣后掠角和后緣后掠角分別固定為40°和20°，展長(zhǎng)固定為3.7m。傾角從0°變化到45°，計(jì)算垂尾模型的RCS曲線，如圖5所示。并將垂尾模型在俯仰角范圍為-20°～20°內(nèi)的RCS平均值繪制成圖6。

從圖6中可以看出，垂尾RCS均值隨傾角θ的增大而減小，RCS值與傾角間的關(guān)系可近似擬合為圖6中的二次曲線，表達(dá)式如下：

通過(guò)二次曲線及式（8）可知，RCS均值隨垂尾傾角的增大近似線性減小，但需要指出，垂尾傾角的增大會(huì)嚴(yán)重影響飛機(jī)的航向穩(wěn)定性及操縱性。

3 垂尾RCS減縮方法

3.1 優(yōu)化設(shè)計(jì)垂尾參數(shù)

通過(guò)以上對(duì)垂尾電磁散射特性的分析可知，垂尾前緣后掠角χ增大值,RCS減?。淮刮舱归L(zhǎng)bc增大，RCS值增大；垂尾傾角θ增大，RCS減小。于是，可以在滿足垂尾氣動(dòng)特性、強(qiáng)度要求等條件下，使垂尾前緣后掠角θ盡可能大，展長(zhǎng)bc盡可能小，傾角θ盡可能大，從而使垂尾雷達(dá)散射截面積RCS在滿足必要?dú)鈩?dòng)特性、強(qiáng)度要求等條件的情況下達(dá)到最小。

3.2 垂尾采用吸波材料

另一種減縮垂尾RCS的方法是采用吸波材料。采用吸波材料的方式一般有兩種，一種是在垂尾的表面涂覆吸波材料；一種方式是采用結(jié)構(gòu)型吸波材料。

4 結(jié)論語(yǔ)

本文利用矩量法計(jì)算垂尾的RCS，得出了垂尾前緣后掠角χ、展長(zhǎng)bc和傾角θ和垂尾RCS間的關(guān)系：χ增大值,RCS減??；bc增大，RCS值增大；θ增大，RCS減小。從而得出了在滿足垂尾氣動(dòng)特性、強(qiáng)度要求等條件的情況下，可以通過(guò)恰當(dāng)?shù)脑龃螃郑瑴p小bc，增大θ來(lái)優(yōu)化垂尾的幾何參數(shù)來(lái)達(dá)到使其RCS最小。另外，可通過(guò)采用吸波材料來(lái)降低垂尾的RCS。

［1］張考,馬東立.軍用飛機(jī)生存力與隱身設(shè)計(jì)[M].國(guó)防工業(yè)出版社,2002.

［2］桑建華.飛行器隱身技術(shù)[M].航空工業(yè)出版社,2013.

［3］李毅.雷達(dá)隱身目標(biāo)電磁散射計(jì)算與實(shí)驗(yàn)研究[D].國(guó)防科學(xué)技術(shù)大學(xué),2007.

［4］焦子涵,張彬乾,沈冬.翼型幾何參數(shù)對(duì)隱身特性的影響研究[J].機(jī)械科學(xué)與技術(shù),2012.12.

［5］李啟鵬,王和平,孫珍,付偉.鴨翼電磁散 射特性分析與RCS減縮方法研究[J].航空計(jì)算技術(shù),2010.

［6］Weinmann F.Ray Tracing With PO/PTD for RCS Modeling of Large Complex Objects[J].Antennas and Propagation'IEEE'2006(6)：1797-1806.