摘 要：航空發(fā)動(dòng)機(jī)排氣系統(tǒng)是作戰(zhàn)飛機(jī)后向電磁散射的主要貢獻(xiàn)源，在排氣系統(tǒng)涂敷雷達(dá)吸波材料能夠有效降低其后向雷達(dá)散射截面積。為了研究雷達(dá)吸波材料涂敷方案對軸對稱塞式噴管電磁散射特性的影響，設(shè)計(jì)了7種不同的涂敷方案，獲得了噴管后向雷達(dá)散射截面積分布。數(shù)值研究結(jié)果表明：在噴管壁面和塞錐壁面涂敷雷達(dá)吸波材料均能夠降低噴管后向的RCS；在噴管出口區(qū)域涂敷能夠在使用較少的雷達(dá)吸波材料的前提下，獲得較好的RCS縮減效果；在塞錐前部區(qū)域涂敷雷達(dá)吸波材料能夠提供最大的單位面積縮減效率。

關(guān)鍵詞：雷達(dá)吸波材料；軸對稱；塞式噴管；雷達(dá)散射截面積；涂敷方案

中圖分類號：V231" " " " " " "文獻(xiàn)標(biāo)識碼：A" " " " " " " 文章編號：1007 - 9734 （2024） 06 - 0031 - 07

0 引 言

隨著戰(zhàn)爭形態(tài)的改變，下一代作戰(zhàn)飛行器需要進(jìn)一步提升自身的機(jī)動(dòng)能力和隱身能力，從而提升作戰(zhàn)飛行器總體的生存率和隱蔽性［1］。航空發(fā)動(dòng)機(jī)的尾噴管是典型的電大尺寸腔體結(jié)構(gòu)，是飛行器后向電磁散射的主要貢獻(xiàn)源［2］。目前，常用于航空發(fā)動(dòng)機(jī)噴管雷達(dá)散射截面積（RCS， radar cross section）縮減的措施主要是外形修形和涂敷雷達(dá)吸波材料（RAM， radar absorbing material）。雷達(dá)吸波材料通過將入射到吸波材料上的電磁波所攜帶的電磁能轉(zhuǎn)化成熱能耗散或者使電磁波因干涉而消失［3］。Mosallaei［4］等人采用遺傳算法得到吸波材料在矩形、球形等規(guī)則腔體的最優(yōu)涂敷組合形式，獲得到最佳RCS縮減效果。何小祥［5］等人建立了非完全理想導(dǎo)體邊界的電磁散射問題理論模型，并將其開發(fā)的方法應(yīng)用于內(nèi)壁涂敷吸波材料的電大尺寸腔體的電磁散射分析中。王龍［6］等研究了不同電參數(shù)的雷達(dá)吸波材料涂敷對S彎進(jìn)氣道的電磁特性影響，獲得了電參數(shù)對進(jìn)氣道RCS的影響規(guī)律。西北工業(yè)大學(xué)楊青真［7］［8］團(tuán)隊(duì)研究了雷達(dá)吸波材料涂敷位置對S彎噴管、球面收斂二元矢量噴管電磁散射特性的影響，以及不同涂敷位置對雷達(dá)吸波材料RCS縮減作用的影響。

塞式噴管作為一種機(jī)械調(diào)節(jié)式矢量噴管，具有較好的氣動(dòng)性能和紅外抑制特性［9］，同時(shí)在電磁散射抑制方面也具有一定的優(yōu)勢；與傳統(tǒng)鐘形噴管相比，塞式噴管可以作為一種高度補(bǔ)償類噴管應(yīng)用于固體火箭發(fā)動(dòng)機(jī)上。國內(nèi)外學(xué)者針對塞式噴管的氣動(dòng)、電磁散射和紅外輻射性能開展了一系列的研究。在氣動(dòng)特性研究方面，國內(nèi)學(xué)者采用實(shí)驗(yàn)和計(jì)算模擬的方式研究了塞式噴管設(shè)計(jì)參數(shù)、冷卻措施對流場流動(dòng)性能的影響。王旭［10］等人采用計(jì)算流體力學(xué)的方法分析了塞式噴管幾何設(shè)計(jì)參數(shù)—塞錐角度和幾何矢量角對噴管氣動(dòng)性能的影響，數(shù)值分析結(jié)果表明高空狀態(tài)下塞式噴管的推力系數(shù)與地面狀態(tài)下相比有更大的損失；增加塞錐后體長度能夠縮短塞錐尾部低壓區(qū)。周兵［11］等人通過試驗(yàn)測試的方式獲得了不同吹風(fēng)比條件下，二元塞式噴管塞錐壁面多斜孔氣膜冷卻對塞錐壁面溫度、噴管出口截面處溫度以及對噴管紅外輻射特性的影響。試驗(yàn)結(jié)果表明，增大吹風(fēng)比可以有效提高氣膜的冷卻效率，降低噴管出口中心區(qū)域溫度，進(jìn)而降低整個(gè)噴管后向的紅外輻射強(qiáng)度。征建生［12］針對塞錐所處工作環(huán)境的特點(diǎn)及其結(jié)構(gòu)特點(diǎn)，提出了塞錐前緣弧面段沖擊冷卻與直段氣膜冷卻相結(jié)合的冷卻結(jié)構(gòu)、塞錐后緣出口截面以內(nèi)采用沖擊冷卻與出口截面以外采用氣膜冷卻相結(jié)合的冷卻方式，并對其影響塞錐表面冷卻特性和噴管氣動(dòng)性能進(jìn)行了數(shù)值研究，研究結(jié)果表明：噴管總壓恢復(fù)系數(shù)和推力系數(shù)隨著吹風(fēng)比、開孔率和矢量偏轉(zhuǎn)的增加而下降。盛超［13］等人開展了熱態(tài)下塞式矢量噴管內(nèi)流動(dòng)特性的試驗(yàn)，獲得了不同落壓比下塞式噴管幾何矢量角對塞式噴管總壓恢復(fù)系數(shù)、推力系數(shù)和氣動(dòng)矢量角的影響規(guī)律。李泓瑾［14］等人采用數(shù)值仿真得到隨海拔高度變化時(shí)不同擴(kuò)張比的塞式噴管推力性能的基本變化規(guī)律，數(shù)值結(jié)果表明塞式噴管的推力優(yōu)勢會(huì)隨著海拔高度的提升有明顯的下降；在低海拔下，推力優(yōu)勢明顯。李修明［15］等人針對導(dǎo)彈推進(jìn)系統(tǒng)性能需求，設(shè)計(jì)了一套環(huán)喉形固體塞式噴管，通過移動(dòng)膨脹唇部位置實(shí)現(xiàn)推力可調(diào)，對其推力性能進(jìn)行了數(shù)值模擬預(yù)估和冷態(tài)試驗(yàn)測試；結(jié)果表明測試結(jié)果與數(shù)值模擬結(jié)果得到的推力性能吻合，塞式噴管推力可實(shí)現(xiàn)4∶1的調(diào)節(jié)變化。在電磁散射方面，陳玲玲［16］等人研究了塞錐錐頂角對塞式噴管后向電磁散射特性的影響，數(shù)值計(jì)算結(jié)果表明塞錐的存在會(huì)增大0°探測角下噴管的雷達(dá)散射截面積；但是噴管后向整體的RCS均值得到了10%以上的降低。

綜合已有文獻(xiàn)發(fā)現(xiàn)，現(xiàn)有的研究主要集中在研究塞式噴管設(shè)計(jì)參數(shù)對噴管氣動(dòng)性能特性和目標(biāo)特性的研究，對RCS縮減措施尤其是雷達(dá)吸波材料涂敷對電磁散射特性影響的研究較少。為此，本文采用彈跳射線法對設(shè)計(jì)的7種雷達(dá)吸波材料涂敷方案的噴管后向RCS進(jìn)行了計(jì)算分析。

1 計(jì)算模型與方法

1.1" 計(jì)算模型

本文中塞錐的前錐半錐頂角β固定為30°，塞錐后錐半錐頂角α固定為15°，塞錐半徑比固定為0.25（塞錐最大半徑r與內(nèi)涵進(jìn)口半徑R之比）。本文計(jì)算所用的軸對稱塞式噴管如圖1所示。

1.2" 雷達(dá)吸波材料涂敷方案設(shè)計(jì)

為研究介質(zhì)在不同涂敷區(qū)域?qū)S對稱塞式噴管RCS的影響，且在仿真模擬中便于在不同區(qū)域進(jìn)行介質(zhì)涂敷，在設(shè)計(jì)涂敷方案時(shí)，需要考慮中心錐外表面和噴管殼體內(nèi)外壁面，因此本文對噴管殼體及塞錐均進(jìn)行了區(qū)域劃分，分別劃分出編碼為1、2、3、4、5、6的六個(gè)區(qū)域。對噴管殼體進(jìn)行區(qū)域劃分，分別在距進(jìn)氣口截面300mm處、600mm處、1000mm處、1400mm處、1700mm處建立平面，并對殼體進(jìn)行拆分體操作，殼體區(qū)域劃分如圖2所示。

對塞錐進(jìn)行區(qū)域劃分，根據(jù)塞錐自身幾何特征，分別在距塞錐后錐頂點(diǎn)200mm處、400mm處、460mm處、535mm處、635mm處建立平面，并對塞錐進(jìn)行拆分體操作，塞錐區(qū)域劃分如圖3所示。

本文根據(jù)所劃分區(qū)域共設(shè)計(jì)8個(gè)涂敷方案，分別為涂敷外殼1、2區(qū)域，涂敷外殼3、4區(qū)域，涂敷外殼5、6區(qū)域，涂敷塞錐1、2區(qū)域，涂敷塞錐3、4區(qū)域，涂敷塞錐5、6區(qū)域，以下分別稱為Model 2、Model 3、Model 4、Model 5、Model 6、Model 7。一組對比方案不予涂敷，稱為Model 0；一組對比方案每個(gè)區(qū)域進(jìn)行全涂敷，稱為Model 1。以全涂敷方案的涂敷面積為基準(zhǔn)面積１，其他涂敷方案的面積與全涂敷相比確定其涂敷面積比例。噴管外壁面涂敷方案和塞錐涂敷方案分別如表1和表2所示。

2 計(jì)算方法

2.1" 彈跳射線法

本文采用彈跳射線算法模擬涂敷吸波材料的電大尺寸腔體電磁散射特性計(jì)算，研究不同涂敷位置對RCS的影響。首先用幾何光學(xué)射線來模擬入射電磁波在物體表面的多次反射，追蹤射線在傳播過程中的場強(qiáng)和相位變化，最后當(dāng)射線離開物體時(shí)利用物理光學(xué)求解散射場，解決了物理光學(xué)法沒有考慮電磁波在目標(biāo)表面多次反射的問題，該算法具有簡練容易實(shí)現(xiàn)、物理推導(dǎo)清晰、誤差小等優(yōu)點(diǎn)[17]。

根據(jù)彈跳射線法，目標(biāo)的總散射場為[18][19]：

[Esr，ω=jωZ0e-jrωc4πrcs×ss×n×Hr，ω+1Z0Er，ω×nejs?riωcds]" " （1）

式中E（[r]，ω）和H（[r]，ω）分別為目標(biāo)表面處的總電場和總磁場。

2.2" 雷達(dá)吸波材料RCS縮減能力表征參數(shù)

在本文中，為準(zhǔn)確描述吸波材料不同涂敷方案對于軸對稱塞式噴管RCS均值縮減能力的影響，定義如下兩個(gè)特征參數(shù)：縮減效果A[5]、單位面積提供的縮減效果B。

A= [σ0-σnσ0-σall] （2）

B= [Ak] （3）

式中：[k]為不同涂敷方案下涂敷面積所占比例；[σn]為不同涂敷方案下模型的RCS均值；[σ0]為未涂敷模型的RCS均值；[σall]為全涂敷方案下模型的RCS均值。

3 計(jì)算結(jié)果與分析

本文中，入射雷達(dá)波頻率為10GHz，入射角范圍為0°～30°，入射角設(shè)置如圖4所示。

3.1噴管壁面涂敷方案對RCS的影響分析

圖5所示為水平極化方式下不同噴管壁面涂敷設(shè)計(jì)方案RCS角向分布曲線。由圖可知，在大部分入射角下，在噴管外壁面涂敷雷達(dá)吸波材料能夠降低塞式噴管后向的RCS幅值，Model 1具有最優(yōu)的RCS縮減能力。在0°～5°入射角范圍內(nèi)，三種壁面涂敷方案與Model 0的RCS幅值差異較小。在5°～10°入射角范圍內(nèi)，Model 2的RCS角向分布規(guī)律與Model 0接近且兩者的RCS幅值差較小，Model 3和Model 4與Model 1的角向分布規(guī)律存在比較大的差異。在該入射角范圍內(nèi)，這兩種涂敷方案的RCS縮減能力較為明顯，這主要是因?yàn)樵谛∪肷浣窍?，入射電磁波直接照射的區(qū)域主要是喉道及噴管出口區(qū)域，Model 2中在該區(qū)域無涂敷。當(dāng)入射角大于10°之后，Model 2、Model 3和Model 4的后向RCS角向分布與Model 0接近，且RCS幅值差異較小；Model 2在更多的入射角下的RCS大于Model 3和Model 4，即Model 2在大入射角范圍內(nèi)的RCS縮減能力較差，這主要是因?yàn)槿肷涞組odel 2涂敷區(qū)域的電磁波角度與壁面法向方向偏離較大，降低了吸波材料的RCS縮減能力。

圖6是垂直極化方式下不同噴管壁面涂敷設(shè)計(jì)方案RCS角向分布曲線。由圖可知，在靠近噴管入口區(qū)域涂敷吸波材料對噴管RCS的角向分布規(guī)律和幅值的影響均較小。在5°～15°入射角范圍內(nèi)，Model 2和Model 3展現(xiàn)了較好的RCS縮減能力。當(dāng)入射角大于15°之后，在壁面區(qū)域涂敷吸波材料的RCS縮減效果較差。

表3為兩種極化方式下不同壁面區(qū)域雷達(dá)吸波材料涂敷方案的RCS均值及特征參數(shù)。由表可知，在兩種極化方式下全涂敷方案都可以達(dá)到較好的RCS縮減效果。三種壁面涂敷方案，在水平極化方式下，Model 3的RCS均值最小，Model 4的單位面積提供的RCS縮減效率最高；在垂直極化方式下，Model 4的RCS均值最小，單位面積提供的RCS縮減效率也是最高的。綜合來看，三種涂敷方案中，Model 4的RCS縮減效果較好。

3.2" 塞錐表面涂敷方案對RCS的影響

圖7為水平極化方式下不同塞錐區(qū)域涂敷設(shè)計(jì)方案的RCS角向分布曲線。由圖可知，在大部分入射角范圍內(nèi)，Model 6、Model 7的RCS角向分布規(guī)律與Model 0接近，與Model 0的RCS幅值差較??；在5°～10°入射角范圍內(nèi)，Model 5與其他涂敷方案的RCS幅值差較大。在部分入射角下，塞錐區(qū)域涂敷設(shè)計(jì)方案的RCS幅值會(huì)大于Model 0，這可能是因?yàn)樵谶@些入射角下，入射雷達(dá)波在經(jīng)過塞錐與壁面的反射后會(huì)在塞錐前部區(qū)域匯聚，這部分電磁波所攜帶的能量被吸波材料吸收，降低了返回的電磁波的電場強(qiáng)度，進(jìn)而降低了后向的RCS幅值。在入射角大于15°之后，在塞錐任何區(qū)域涂敷雷達(dá)吸波材料對軸對稱噴管后向RCS的分布規(guī)律和幅值的影響都較小。

圖8為垂直極化方式下不同塞錐區(qū)域涂敷設(shè)計(jì)方案的噴管RCS角向分布曲線。

由圖可知，在5°入射角附近， RCS角向分布曲線在此入射角下出現(xiàn)一個(gè)波谷；Model 6、Model 7在此入射角下則出現(xiàn)一個(gè)波峰；在5°～15°入射角范圍內(nèi)，Model 5的RCS幅值低于Model 0、Model 6和Model 7，這主要是因?yàn)樵谠撊肷浣欠秶鷥?nèi)，Model 5所涂敷的塞錐區(qū)域處于入射電磁波經(jīng)過壁面反射之后的區(qū)域。當(dāng)入射角大于15°時(shí)，在塞錐區(qū)域涂敷雷達(dá)吸波材料還能保持一定的RCS縮減能力，Model 6和Model 7的RCS分布規(guī)律接近，與Model 5的分布規(guī)律存在差異。

表4為兩種極化方式下三種不同塞錐涂敷方案的RCS均值和縮減效果特征參數(shù)。由表可知，在兩種極化方式下，Model 5的RCS均值最小，Model 6的RCS均值最大；在塞錐區(qū)域涂敷雷達(dá)吸波材料，垂直極化方式下涂敷方案的RCS縮減能力優(yōu)于水平極化方式下。所有涂敷方案中，Model 5單位面積的RCS縮減效率最大。

4 結(jié)論與分析

本文采用彈跳射線法對七種不同涂敷方案的軸對稱塞式噴管的后向RCS特性進(jìn)行了分析，結(jié)果表明：

（1）在噴管壁面和塞錐區(qū)域涂敷雷達(dá)吸波材料能夠達(dá)到最好的RCS縮減效果。

（2）在兩種極化方式下，在噴管喉道附近區(qū)域涂敷雷達(dá)吸波材料能夠達(dá)到的RCS縮減效果較好，RCS縮減效果可以達(dá)到全涂敷方案的49%，但是涂敷面積僅占37%。

（3）在兩種極化方式下，在塞錐前部區(qū)域涂敷雷達(dá)吸波材料設(shè)計(jì)方案的RCS縮減效果可以達(dá)到全涂敷方案的31%，涂敷面積僅占1.4%；Model 5方案的單位面積縮減效率最高。

參考文獻(xiàn)：

［1］桑建華.飛行器隱身技術(shù)［M］.北京：航空工業(yè)出版社，2013.

［2］艾俊強(qiáng)，周莉，楊青真.S彎隱身噴管［M］.北京：國防工業(yè)出版社，2017.

［3］姬金祖，黃沛霖，馬云鵬，等.隱身原理［M］.北京：北京航空航天大學(xué)出版社，2018.

［4］MOSALLAEI H，RAHMAT-SAMII Y. RCS reduction of canonical targets using genetic algorithm synthesized RAM ［J］.IEEE Transactions on Antennas amp; Propagation，2000，48（10）：1594-1606.

［5］何小祥，徐金平，顧長青.介質(zhì)涂敷電大腔體電磁散射IPO研究［J］.電子與信息學(xué)報(bào)，2005，27（1）：136-138.

［6］王龍，鐘易成，張堃元.金屬/介質(zhì)涂覆的S形擴(kuò)壓器電磁散射特性［J］.物理學(xué)報(bào)，2012，61（23）：193-201.

［7］高翔，施永強(qiáng)，楊青真，等.介質(zhì)涂覆位置對雙S彎排氣系統(tǒng)電磁散射特性影響研究［J］.物理學(xué)報(bào)，2015，64（2）：106-115.

［8］郭霄，楊青真，施永強(qiáng)，等.介質(zhì)涂覆位置對球面收斂噴管電磁散射特性影響［J］.航空學(xué)報(bào)，2017，38（4）：245-253.

［9］陳靜，單勇，張靖周，等.引氣結(jié)構(gòu)對二元塞式噴管冷卻和紅外輻射特性的影響［J/OL］.推進(jìn)技術(shù)，2023，44（11）：1-14.

［10］王旭，張靖周，單勇.塞錐形狀和偏轉(zhuǎn)角對軸對稱塞式噴管氣動(dòng)性能的影響［J］.航空發(fā)動(dòng)機(jī)，2016，42（1）：11-15.

［11］周兵，吉洪湖.塞錐氣膜冷卻對二元塞式噴管紅外特征的影響［J］.航空動(dòng)力學(xué)報(bào)，2016，31（12）：2895-2903.

［12］征建生.二元矢量塞式噴管塞錐表面冷卻特性研究［D］.南京：南京航空航天大學(xué)，2017.

［13］盛超，滕狀，李慶林，等.塞式矢量噴管熱態(tài)內(nèi)流特性試驗(yàn)［J］.航空發(fā)動(dòng)機(jī)，2020，46（6）：29-33.

［14］李泓瑾，謝侃，王寧飛.塞式噴管推力性能研究［A］.中國科協(xié)航空發(fā)動(dòng)機(jī)產(chǎn)學(xué)聯(lián)合體.第六屆空天動(dòng)力聯(lián)合會(huì)議暨中國航天第三專業(yè)信息網(wǎng)第四十二屆技術(shù)交流會(huì)暨2021航空發(fā)動(dòng)機(jī)技術(shù)發(fā)展高層論壇論文集（第二冊）［C］.2022.

［15］李修明，王一白，童悅，等.固體火箭發(fā)動(dòng)機(jī)推力可調(diào)塞式噴管數(shù)值模擬與冷流試驗(yàn)［J］.上海航天，2019，36（S1）：29-34.

［16］陳玲玲，楊青真，陳立海，等.軸對稱塞式噴管電磁散射特性數(shù)值模擬［J］.推進(jìn)技術(shù)，2018，39（1）：86-91.

［17］冉靜.典型復(fù)合目標(biāo)電磁散射特性分析［D］.西安：西安電子科技大學(xué)，2022.

［18］馮恬恬，郭立新.基于彈跳射線法的電大尺寸部分涂覆目標(biāo)的RCS計(jì)算［A］//中國電子學(xué)會(huì).2020年全國微波毫米波會(huì)議論文集（上冊）［C］.2020.

［19］蘇遠(yuǎn).金屬/介質(zhì)目標(biāo)電磁散射分析的彈跳射線法及在GPU計(jì)算平臺上的實(shí)現(xiàn)［D］.南京：南京理工大學(xué)，2018.

責(zé)任編校：裴媛慧，陳 強(qiáng)

Influence of Radar Absorbing Material on the RCS of Axisymmetric Plug Nozzle

GUO Xiao1，XIA Xinru1，CHEN Haobo2，LI Heng3，ZHAO Zhujun1

（1.School of Aeronautical Engineering，Zhengzhou University of Aeronautics，Zhengzhou 450046，China；

2.AVIC Jonhon Optronic Technology Co.Ltd.，Luoyang 471000，China；

3.Sichuan Gas Turbine Research Establishment，Aero Engine Corporation of China，Chengdu 610500，China）

Abstract：The aeroengine exhaust system" is the main contributor to the backward electromagnetic scattering of combat aircraft.It can effectively reduce the cross-sectional area of backward radar scattering by coating the radar absorbing material.Seven kinds of radar absorbing material （RAM） coating schemes were designed to get out the influence on radar cross section（RCS） of axisymmetric plug nozzle.The backward radar cross section distribution was studied by shooting and bouncing ray method.The numerical simulation results show that it can reduce the RCS by coating RAM on the nozzle wall surface and the plug wall surface.The better backward RCS reduction effect can be reached that coating RAM on the outlet aera under the premise by using less RAM.It can provide the largest reduction efficiency per unit area by coating RAM on the front plug area.

Key words：radar absorbing material;axisymmetric;plug nozzle;radar cross section；coating scheme