向迎春 曲長文 平殿發(fā) 趙維青
(海軍航空工程學(xué)院電子信息工程系1) 煙臺 264001)(中國人民解放軍91599部隊(duì)2) 萊陽 265200)
艦船雷達(dá)波隱身技術(shù)是提高艦船在現(xiàn)代海戰(zhàn)中生存能力的有效手段。發(fā)達(dá)國家海軍目前在艦船建造中普遍采用隱身外形和涂覆隱身材料等新技術(shù)來反雷達(dá)探測。艦船雷達(dá)波隱身技術(shù)已從試驗(yàn)階段步入實(shí)用階段,歐美等國已相繼有多艘采用雷達(dá)隱身設(shè)計(jì)的艦船投入現(xiàn)役??紤]代價(jià)因數(shù)以及進(jìn)一步發(fā)展隱身外形和隱身材料的困難性,研究新的反雷達(dá)隱身技術(shù)就顯得尤其重要。有源對消隱身就是艦船雷達(dá)波隱身技術(shù)發(fā)展的一個重要方向。
艦船對電磁波的反射基于艦船表面感應(yīng)電流的二次輻射。若已知源的輻射場分布,反過來又可以得到散射體的特性和場源分布。而從艦船對雷達(dá)的隱身看,往往只需雷達(dá)所在的一個小的立體角范圍內(nèi)實(shí)現(xiàn)電磁波對消,此時也只需給出這一特定區(qū)域的遠(yuǎn)區(qū)電磁場數(shù)值即可。雷達(dá)散射截面(RCS)就是對艦船在一定方向上反射的功率與入射的功率密度進(jìn)行歸一化處理后得出的度量。它表征艦船表面感應(yīng)電流的遠(yuǎn)區(qū)輻射場,其數(shù)學(xué)表示式為[1]
從理論上講,艦船在不同方向的 RCS數(shù)值可以精確地測到,而瞬態(tài)的雷達(dá)入射場也可以通過平臺上的傳感器精確測到。這樣根據(jù)式(1)的定義,艦船在雷達(dá)方向上的散射場就確定了,即
根據(jù)式(2)在雷達(dá)方向與散射場相干但相位相反的場總是可以產(chǎn)生出來的,因此,有源電磁波對消在理論上是成立的。提高艦船的雷達(dá)隱身性能,就是減小雷達(dá)能接收到的散射電磁波強(qiáng)度,也就是減弱艦船對雷達(dá)波的散射即雷達(dá)散射截面??梢?通過有源或無源的手段,降低投照到艦船的信號強(qiáng)度,則從艦船散射回雷達(dá)的信號強(qiáng)度相應(yīng)地更弱了;或者只降低艦船散射回波,這種技術(shù)途徑也可以達(dá)到雷達(dá)散射截面縮減(RCSR)的目的。
1)對消幅相條件
復(fù)雜艦船可被分解為 N個離散散射體或散射中心的組合,其σ取決于各散射中心的散射截面積的相關(guān)和[2]:
式中,σn是第n個散射中心的復(fù)數(shù) RCS,φn是該散射中心的相對相位。在眾多的散射中心中,對于特定的工作頻率和入射角,一般存在著若干起主要作用的散射中心。集中減縮主要散射中心對散射所做的貢獻(xiàn),可有效地實(shí)現(xiàn)RCSR。人為引入與主要散射中心相對應(yīng)的對消源是有源對消的出發(fā)點(diǎn)。艦船平臺上的傳感器將入射雷達(dá)波的參數(shù)送入控制系統(tǒng),自適應(yīng)地產(chǎn)生并在所需方向發(fā)射一列合適的相干波,改變艦船的散射分布,有效地降低雷達(dá)方向的散射功率密度。
設(shè)艦船原有的雷達(dá)散射截面積為σ0,人為通過有源方式引入一個等效的散射中心,其等效散射截面積為σ1,它們的相位分別為φ0和φ1,疊加后的艦船散射截面積σ為
即
控制σ1和φ1,就可能進(jìn)行尋優(yōu),以實(shí)現(xiàn)最佳對消效果。當(dāng)達(dá)到最佳對消參數(shù):
2)相位幅值差對隱身效果的影響
若相位與幅值同時存在誤差時,則歸一化雷達(dá)截面ˉσ隨相位與幅值誤差的變化關(guān)系為
圖1 相位和幅值偏離最佳值時的歸一化RCSR曲面
從圖中可以看出當(dāng)相位與幅值相對最佳條件都有變化時,要對被保護(hù)艦船達(dá)到一定的隱身效果,則相位與幅值差偏離最佳條件的范圍被進(jìn)一步的縮小,而對相位的要求就更苛刻。在圖中坐標(biāo)原點(diǎn)處幅值和相位均為最佳條件,此時雷達(dá)截面積為0,實(shí)現(xiàn)了完全隱身。
事實(shí)上,對海面艦船的散射研究尤其是掠射條件下的散射特性分析長期以來一直是計(jì)算電磁學(xué)的一個重要研究課題。同時從上面的原理分析也可以看到,艦船有源對消隱身的實(shí)現(xiàn)并不容易。它有很多關(guān)鍵技術(shù)需要突破:艦船散射特性的精確預(yù)估,對消場的生成與控制以及對消數(shù)據(jù)庫的建立等,而首先要做的工作就是艦船散射特性的精確預(yù)估。
為了能夠準(zhǔn)確計(jì)算海面艦船的散射特性必需要盡可能考慮海洋背景對雷達(dá)波散射的影響。艦船的RCS預(yù)估思想為:以粗糙面為海面模型,以三維外形圖為艦船模型,在此基礎(chǔ)上研究海面艦船散射。
1)艦船的三維建模?,F(xiàn)在流行的圖形建模軟件有 Auto CAD、3DS Max、Ansys、GiD 等。構(gòu)造艦船外形,可以根據(jù)艦船設(shè)計(jì)圖紙和設(shè)備圖紙精確構(gòu)造,模型的精度可達(dá)毫米量級。而對外方艦船,則可以通過照片等資料構(gòu)造出近似外形。
2)海面的模擬。海面的模擬,模型有統(tǒng)計(jì)模型、分形模型和混沌模型。統(tǒng)計(jì)模型中海面的運(yùn)動用隨機(jī)函數(shù)來描述,Monte Carlo方法是最為常用的統(tǒng)計(jì)方法;Pierson-Moskowitz海譜和Phillip譜在海面粗糙面的統(tǒng)計(jì)構(gòu)造中也有著廣泛的運(yùn)用。分形模型構(gòu)造的粗糙面結(jié)構(gòu)具有自相似性,同時兼顧了大范圍和小范圍的無序的特點(diǎn),接近實(shí)際的自然粗糙面。此外,Shaw和Dougan[3]提出可以結(jié)合相應(yīng)的海域模型以及邊界條件來研究海面散射,他們提出的分類中將海面散射的研究分為七個層次,最終達(dá)到對海面散射的精確分析。
3)對三維艦船和模擬海面進(jìn)行RCS預(yù)估。為了計(jì)算艦船的RCS,需要求出散射艦船在遠(yuǎn)場的散射場。在海面艦船的散射計(jì)算中,將艦船以及其周圍海面看作整體進(jìn)行研究。
1)RCS的一般預(yù)估方法。RCS預(yù)估方法有有限元法(FEM)、矩量法(MOM)等頻域技術(shù);也有時域有限差分法(FDTD)、時域積分方程法等時域技術(shù)以及屬于高頻技術(shù)的幾何光學(xué)法(GO)、物理光學(xué)(PO)、幾何繞射(GTD)、物理繞射(PTD)、射線追蹤(RT)等。對海面散射的研究一般是基于海面積分方程,采用MOM離散得到矩陣方程求解,求解精度高但往往受限于艦船電尺寸大帶來的龐大運(yùn)算量與存儲量,可行性不高。時域有限差分法具有良好的寬帶散射特性分析能力,但同矩量法一樣受限于運(yùn)算量與存儲量。高頻方法簡單高效但精度有限,而且屬于點(diǎn)頻方法,寬帶散射特性分析能力弱。因此國內(nèi)外都在發(fā)展各種快速方法和混和方法。
2)復(fù)雜目標(biāo)RCS計(jì)算工具。最早是將復(fù)雜目標(biāo)劃分為多個簡單幾何體部件分別計(jì)算的部件分解法;Lee等發(fā)展XPatch軟件半臺,不僅計(jì)入物理光學(xué)和物理繞射,還基于彈跳射線法(SBR)計(jì)算多次散射,現(xiàn)已發(fā)展成較為成熟的由高頻近似方法估算復(fù)雜目標(biāo)RCS和成像特征的工具;Rius等人利用圖形硬件加速器發(fā)展了一種實(shí)時計(jì)算目標(biāo)RCS的工具GRECO,該工具雖然僅能給出一次散射與繞射場,但因其計(jì)算速度快而多為采用[4]。這些軟件工具或方法往往受到發(fā)布機(jī)構(gòu)對用戶范圍的各種限制,而其方法本身也并不能完全適應(yīng)快速計(jì)算目標(biāo)與環(huán)境復(fù)合散射的多種問題。因此,混和方法日益受到人們的重視,很多軟件也開始逐漸集成利用各種算法進(jìn)行優(yōu)化計(jì)算,這也是未來電磁場EDA軟件發(fā)展的趨勢。
3)粗糙海面與艦船的RCS的預(yù)估現(xiàn)狀。美國佐治亞技術(shù)研究所、諾斯羅普公司已經(jīng)研制了廣泛應(yīng)用于工業(yè)中的RCS預(yù)估程序。美國海軍研究實(shí)驗(yàn)室(NRL)也在研制RCS預(yù)估程序,并且包括海面影響產(chǎn)生的多路徑效應(yīng)。李守秀[5]等給出了粗略估計(jì)艦船雷達(dá)散射截面的經(jīng)驗(yàn)公式,要得到高精度的雷達(dá)散射截面必須用計(jì)算電磁學(xué)方法計(jì)算。
雷達(dá)波在遠(yuǎn)區(qū)時其入射波可視為平面波,針對海面艦船在平面波照射下的散射問題,徐樂[6]等采用介質(zhì)半空間格林函數(shù)在海面建立電場積分方程,并采用廣義前向后向法(GFB)算法結(jié)合矩量法(MOM)計(jì)算了海面以及海面艦船的感應(yīng)電流分布以及雙站雷達(dá)散射截面,為平面波入射時海面散射問題提供更為嚴(yán)謹(jǐn)?shù)姆治瞿P汀?/p>
在解決艦船與海面耦合散射問題時,陳勇[7]等提出了一種基于混合面元投影(HPP)和物理光學(xué)法(PO)的計(jì)算艦船與海面耦合散射的快速算法。計(jì)算平臺為XEON,3.4GHz CPU、8GB內(nèi)存的工作站,計(jì)算中每個角度平均耗時約15s說明了此方法的高效性。
在粗糙海面與艦船的RCS的快速計(jì)算方面,徐豐[4]等提出雙向射線追蹤法(BART),實(shí)現(xiàn)目標(biāo)與粗糙面共存問題的復(fù)合散射計(jì)算。用BART方法計(jì)算縮比船目標(biāo)與下墊平表面后向散射180個觀測角度,其結(jié)果與矩量法(MOM)和多層快速多極子(M LFMM)比較,MOM方法耗時在2d以上,MLFMM耗時約3.6h,BART方法耗時約9s。
在建造一艘有隱身要求的艦船,設(shè)計(jì)方案能否達(dá)到隱身要求不可能靠實(shí)際測量來驗(yàn)定;而為檢驗(yàn)一艘艦船添加的隱身技術(shù)作實(shí)際測定時也將耗資巨大,因而必須借助于RCS的數(shù)學(xué)模型,通過計(jì)算機(jī)仿真實(shí)現(xiàn)。同時在艦船有源對消系統(tǒng)中,對消數(shù)據(jù)庫的生成和更新也要求RCS的快速精確預(yù)估。因此,RCS的準(zhǔn)確快速預(yù)估對海面艦船的隱身尤其是在有源對消隱身中起著決定性作用。
1)分析艦船雷達(dá)散射閃爍點(diǎn)。我國目前的中型水面艦船采用的是常規(guī)排水型船體,且作為主戰(zhàn)艦船必須布置大量的武器裝備,而這些武器裝備中,眾多的武器、觀通等裝備并未進(jìn)行隱身設(shè)計(jì),在這些艦船的總體設(shè)計(jì)中并未因要達(dá)到隱身性的目的而犧牲其總體性能,因此我們現(xiàn)在還只是對艦船進(jìn)行局部的、最大可能的隱身改造[9]。大型水面艦船的雷達(dá)波散射源可以按以下方式進(jìn)行結(jié)構(gòu)劃分:(1)舷體(包括海面);(2)上層建筑;(3)桅桿及布置在桅桿上的天線;(4)甲板上的武器裝備(包括艦載機(jī))[8]。一般通過消滅一類幾何體(角形反射器)、改造二類(二面角)和三類(平板)幾何體來達(dá)到隱身目的。通過分析艦船RCS雷達(dá)散射特性,尋找艦船雷達(dá)散射的“閃爍點(diǎn)”及其在各角度區(qū)域中的變化,采用有源對消方法來實(shí)現(xiàn)艦船隱身。
2)設(shè)計(jì)對消場發(fā)生器分布。在研制新型艦船時,通過RCS預(yù)估程序?qū)υO(shè)計(jì)中的艦船的雷達(dá)散射特性分析,布置合理的對消場發(fā)生器。
3)獲取有利的對敵突擊陣位。通過對艦船RCS的預(yù)估,可以獲得己方艦船的雷達(dá)散射特性并結(jié)合敵雷達(dá)性能預(yù)知敵艦對我艦的探知狀態(tài),從而采取相應(yīng)對策。通過分析我艦船在各角度下的雷達(dá)散射特性,在接敵時可盡可能規(guī)避敵方探測雷達(dá),以較小的散射特性隱蔽突擊。
進(jìn)入21世紀(jì)以后,艦船雷達(dá)隱身的發(fā)展體現(xiàn)在兩個方面:擴(kuò)展隱身技術(shù)領(lǐng)域和綜合各種隱身技術(shù)即隱身波段向毫米波、亞毫米波、紅外、激光和長波段擴(kuò)展和綜合運(yùn)用隱身外形設(shè)計(jì)、吸波材料、電子對抗措施并利用反輻射導(dǎo)彈/高能激光器或粒子束武器等。目前,我海軍的現(xiàn)役艦船尤其是驅(qū)護(hù)艦,其雷達(dá)隱身技術(shù)還非常缺乏,在雷達(dá)隱身技術(shù)方面還有很大改進(jìn)和提高的余地。隱身艦船雖還沒有經(jīng)過實(shí)戰(zhàn)的考驗(yàn),但可以預(yù)料,在未來海戰(zhàn)中將發(fā)揮越來越重要的作用。國外在此方面已走在了我們前面,我們必須盡快入手,從而縮短與西方海軍強(qiáng)國的差距。因此,有源對消隱身和雷達(dá)散射截面預(yù)估等技術(shù)值得進(jìn)一步研究。
[1]E.F.Knott.Radar Cross Section[M].Dedham,MA:Artech House,Inc.,1985
[2]E.F.Knott.RCSR Guidelines Handbook[R].Final Technical Report on EES/GIT Project A-1560-001,Engineering Experiment Station,Georgia Institute of Technology,1976,4
[3]W T Shaw,A J Dougan.Green's function refinement as an approach to radar backscatter:general theory and applications to LGA scattering from the ocean[J].IEEE Trans.On Antennas and Propagat.,1998,46(1):57~66
[4]徐豐,金亞秋.計(jì)算粗糙海面與大型艦船復(fù)合散射的雙向解析射線追蹤法[J].自然科學(xué)進(jìn)展,2008,18(7):814~825
[5]李守秀,王林.艦船雷達(dá)散射截面的計(jì)算公式[J].海軍航空工程學(xué)院學(xué)報(bào),2007,22(3):326~328
[6]徐樂,謝擁軍,史小衛(wèi).海上艦船在平面波照射下的半空間散射研究[J].電子學(xué)報(bào),2006,34(12):2322~2325
[7]陳勇,董純柱,王超,等.基于HPP/PO的艦船與海面耦合散射快速算法[J].系統(tǒng)工程與電子技術(shù),2008,30(4):449~592
[8]吳楠,黃松高.大型水面艦船雷達(dá)波散射特征分析[J].裝備環(huán)境工程,2008,5(1):40~43
[9]黃祥兵,陳霖,王志國.艦船雷達(dá)隱身技術(shù)研究[J].船海工程,2002(2):1~4