申寧 魏兵

(1. 西安電子科技大學(xué)物理與光電工程學(xué)院,西安 710071; 2. 西安電子科技大學(xué)信息感知技術(shù)協(xié)同創(chuàng)新中心, 西安 710071; 3. 西安郵電大學(xué)電子工程學(xué)院,西安 710121)

引 言

隨著短脈沖和超寬帶雷達在工程中的廣泛應(yīng)用,目標時域電磁散射的計算分析越來越受到重視.常用的時域數(shù)值方法[1]有時域矩量法(Time Domain-Method of Moment,TD-MOM),時域有限元法(Time Domain-Finite Element Method,TD-FEM),時域有限差分(Finite Difference Time Domain,FDTD)方法等.數(shù)值方法計算精度高,但是要消耗較多的內(nèi)存,適用于求解電尺寸相對較小目標的電磁問題.電大尺寸目標電磁特性的研究常常采用高頻方法.常見的時域高頻方法有時域物理光學(xué)(Time Domain Physical Optics,TDPO)法[2],時域等效邊緣電流(Time Domain Equivalent Electric Current,TDEEC)方法,時域幾何繞射(Time Domain Geometric Theory of Diffraction,TDGTD)方法,時域彈跳射線(Time Domain Shooting and Bouncing Rays,TDSBR)方法等等.與時域數(shù)值方法相比,時域高頻近似方法具有計算內(nèi)存需求小、計算速度快的優(yōu)點.

TDSBR是彈跳射線(Shooting and Bouncing Rays,SBR)法和TDPO方法相結(jié)合的時域算法.SBR法將入射的平面波等效為大量離散的射線管來模擬入射電磁波在目標中傳播的情況,并采用TDPO方法積分求解散射場,特別適合于計算目標幾何結(jié)構(gòu)之間的多次反射場,具有物理概念清晰、精度高、容易實現(xiàn)等優(yōu)點.傳統(tǒng)的SBR方法為滿足準確度要求,射線管的尺寸至少為入射電磁波波長的十分之一.對于電大復(fù)雜目標,射線數(shù)量龐大,計算效率低.為提高SBR方法的計算效率,許多學(xué)者提出了改進方案.Suk[3]等于2001年提出了多分辨率網(wǎng)格法,多分辨率加速技術(shù)提高了計算效率.然而該算法中射線管太大可能導(dǎo)致SBR方法無法捕捉到目標的一些細小幾何結(jié)構(gòu),有損計算精度.2013年殷宏成[4]等提出了一種自適應(yīng)射線管分裂算法,將射線管依據(jù)其在實際追蹤過程中的情況分裂成若干個子射線管,這種方式不僅能保證計算精度,還能顯著降低待追蹤射線管數(shù)量.另外為了加速射線和目標求交測試的效率,Gao等人利用KD-tree進行加速求交[5],Tao等人用GPU對射線跟蹤進行了加速[6].

本文將基于解析積分的TDPO算法和SBR方法相結(jié)合,給出了一種基于解析積分的TDSBR算法.與傳統(tǒng)的TDSBR算法相比,本文方法可以明顯減少射線管數(shù)量,在保證計算精確度的前提下提高電大尺寸目標電磁特性的計算效率.

1 時域彈跳射線(TDSBR)法

TDSBR方法是一種求解復(fù)雜目標多次散射的經(jīng)典時域高頻方法.TDSBR方法基本原理是:首先是射線管的生成,在垂直于入射方向上定義虛擬孔徑面(等相位面).考慮到計算的精確性,虛擬孔徑面需以波長的十分之一來劃分射線管,本文射線管形狀取為三角形.其次,對每一射線管的角點射線進行追蹤,直到射出目標為止,記錄角點射線的最終出射位置.然后,追蹤射線管中心點發(fā)出的射線,與目標相交時利用幾何光學(xué)法計算反射場.同時考慮射線管發(fā)散、極化和介質(zhì)反射系數(shù)等因素,并記錄射線管中心點的出射場;最后,采用解析TDPO積分計算該射線管在接收方向上的散射場,將所有射線管的遠區(qū)散射場疊加,即得到目標總的遠區(qū)散射場.彈跳射線方法主要包括三個主要部分:射線管的生成、射線管的追蹤、口徑面相關(guān)的積分計算,如圖1所示.

圖1 時域彈跳射線方法

1.1 射線場強的追蹤

在角點射線跟蹤結(jié)束后,確定有效射線管,選取其中心射線的場強進行跟蹤計算,直到中心射線不再與目標相交或達到終止條件才停止追蹤.射線管在傳播過程中,第i+1次的交點ri+1處的場值可以通過第i次的交點ri處的場值遞推:

E(ri+1)=(DF)i(R)iE(ri)e-j,φ.

(1)

式中:(R)i為ri處的反射系數(shù)矩陣;由傳播路徑引起的相位變化

φ=kr·(ri+1-ri)+ki·ri;

(2)

(DF)i為ri處的散度因子,用于描述射線管束中場強的變化,由幾何光學(xué)可得到.在射線管發(fā)出時,面元被賦予了一定的能量,傳播過程中總能量不變,隨著橫截面積發(fā)生變化,隨之改變的是單位面積的能量,散度因子(DF)i就是用來描述這種變化[7],

(3)

1.2 遠場解析積分計算

在入射場已知的情況下,完成射線追蹤后,可由TDPO積分求得每一個射線管在遠區(qū)產(chǎn)生的散射場.

TDPO算法求解時,積分區(qū)域就是射線管最后一次投影在目標的面積.傳統(tǒng)的TDPO算法為了計算精度,需要將射線管離散劃分取入射電磁波波長的十分之一,射線管投影面元中每一點的入射場大小看作是一致的,或者采用高斯積分的辦法,用三角形射線管面元上多個特征點的場來擬合整個射線管的遠場散射.本文采用解析的TDPO算法[8]計算,可明顯減少射線管數(shù)量,提高計算效率.

考慮多次反射時的情況,如圖2所示,TDPO方法觀察點的場是目標上射線管最后一次反射時照亮面積積分的累加結(jié)果,遠區(qū)散射電場可以由式(4)計算

(4)

圖2 TDPO算法各參量示意圖

(5)

式中:ki為射線管初始入射方向;kri為射線管第i次反射時的方向,也是第i+1次反射時的入射方向;rre為源點所在的空間位置.

對于單個射線管,式(4)中的入射磁場方向、射線管面元法向以及散射方向都是已知的,即散射場公式可以改寫成

(6)

式(6)中ds在積分計算時需向各個坐標面投影，如果投影到xoy平面，時間延遲用位置坐標表示，即

t-τ1-τ2-τ3=t+ax+by+c.

(7)

式中a、b、c為常數(shù)，可由式(5)計算得到.

入射波為微分高斯脈沖時,經(jīng)過積分運算,最終得到散射場表達式如下[8]:

(8)

其中

(9)

式中:τ為入射波脈沖寬度;L1、L2和a1x+b1、a2x+b2分別為三角形射線管投影到xoy坐標面后沿著x軸的頂點坐標位置和對應(yīng)的三角形邊的直線方程．對單個射線管而言,在式(8)的求解過程中脈沖寬度、時間變量均由外界給定,與射線管本身尺寸大小無關(guān).因此,在保證計算精度的情況下盡可能使射線管劃分尺寸放大,從而減少對計算機內(nèi)存的需求.

2 算例分析

圖3 二面角的后向散射時域波形

在上述情形下,二面角的后向RCS的解析解為[9]

(10)

式中:L、M為矩形二面角反射器的長和寬;λ為入射波波長.圖4為二面角的后向RCS隨頻率的變化情形.圖中實線、三角形和圓點分別表示本文方法、傳統(tǒng)方法和解析結(jié)果.可見,幾種方法的計算結(jié)果吻合.

圖4 二面角的后向RCS

算例2:金屬角反射器的后向RCS.角反射器由三個5 m×5 m的平面組成.入射波為微分高斯脈沖,入射方向θ=45°,φ=45°(見圖5),后向散射.傳統(tǒng)TDSBR算法需要按照λ/10離散大約300 000個射線管,本文方法按照λ/6離散大約108 000個射線管.如圖5所示為采用解析TDSBR計算θθ極化時三面角后向散射的時域波形.

圖5 三面角的后向散射時域波形

圖6為角反射器后向RCS隨頻率的變化曲線與商業(yè)軟件FEKO計算結(jié)果的對比.由圖可見,兩種方法的計算結(jié)果吻合.此時本文方法射線管的數(shù)量僅為傳統(tǒng)TDSBR射線管數(shù)量的1/3.

圖6 三面角的后向RCS

算例3:雙脊喇叭天線的后向散射.喇叭天線的口徑面為12 cm×10 cm,長30 cm.模型離散為3 317個面元,如圖7所示.入射波為微分高斯脈沖,沿著雙脊喇叭天線口徑垂直入射,后向散射.射線跟蹤時反射次數(shù)上限取為3次.考慮到模型的復(fù)雜性,按照λ/10離散約6 000個射線管.圖8為采用解析TDSBR計算θθ極化時后向散射的時域波形.圖9為采用相同射線管劃分時,解析TDSBR、傳統(tǒng)TDSBR方法計算后向RCS隨頻率的變化曲線.從圖中可以看出傳統(tǒng)TDSBR與解析TDSBR計算結(jié)果在入射波頻率為11 GHz以前基本吻合,之后傳統(tǒng)TDSBR計算結(jié)果出現(xiàn)強烈的不穩(wěn)定性.可見本文方法比傳統(tǒng)方法更為穩(wěn)定.

圖7 雙脊喇叭天線模型

圖8 雙脊喇叭天線的后向時域散射波形

圖9 雙脊喇叭天線的后向RCS

圖10為解析TDSBR計算按不同尺度劃分射線管時RCS的對比,從圖中可以看出按照λ/8劃分射線管的計算結(jié)果與λ/10,λ/15劃分是基本一致.圖11 為傳統(tǒng)TDSBR計算按不同尺度劃分射線管時RCS的對比,從圖中可以看出射線管λ/15劃分時,高頻區(qū)域的震蕩才消失.由此可見,對復(fù)雜模型采用解析TDSBR計算結(jié)果比較穩(wěn)定,與傳統(tǒng)TDSBR比較,可用較少的射線管得到較精確的計算結(jié)果.

圖10 解析TDSBR不同射線劃分RCS的對比

圖11 傳統(tǒng)TDSBR不同射線劃分RCS的對比

3 結(jié) 論

時域彈跳射線法的優(yōu)點是在處理電大尺寸目標時計算速度快、所需內(nèi)存少、一次計算能得到寬頻帶RCS.但是龐大的射線管數(shù)目和耗時的射線與目標求交測試,使得計算效率低.文中的時域彈跳射線法采用解析TDPO積分計算目標的遠區(qū)散射場,算例表明在保證計算精度的前提條件下,使得射線管的劃分可不必取入射電磁波波長的十分之一,明顯減少了射線管的數(shù)量,計算穩(wěn)定性較好,提高了計算效率.

[1] 葛德彪, 魏兵. 電磁波時域計算方法[M]. 西安:西安電子科技大學(xué)出版社,2014.

[2] 楊凌霞, 葛德彪, 魏兵. 電大目標電磁散射的TDPO分析 [J]. 電波科學(xué)學(xué)報, 2007, 22(4): 552-556.

YANG L X, GE D B, WEI B. Analysis of large-scale EM scattering problems by using TD-PO method[J]. Chinese journal of radio science, 2007, 22(4): 552-556. (in Chinese)

[3] SUK S H, SEO T I, PARK H S, et al. Multi-resolution grid algorithm in the SBR and its application to the RCS calculation[J]. Microwave and optical technology letters, 2001, 29(6):394-397.

[4] 殷紅成, 朱國慶, 董純柱, 等. 基于自適應(yīng)射線管分裂的多次反射計算方法[J]. 系統(tǒng)工程與電子技術(shù), 2013,35(4):700-705.

YIN H C, ZHU G Q, DONG C Z, et al. Efficient multi-reflection computational method based on adaptive ray tube splitting[J]. Systems engineering and electronics, 2013, 35(4):700-705. (in Chinese)

[5] GAO P C, TAO Y B, LIN H. Fast RCS prediction using multiresolution shooting and bouncing ray method on the GPU[J]. Progress in electromagnetics research, 2010, 107(4): 187-202.

[6] TAO Y B, LIN H, BAO H J. GPU-based shooting and bouncing ray method for fast RCS prediction[J]. IEEE transactions on antennas and propagation, 2010, 58(2): 494-502.

[7] 張磊, 侯兆國, 董純柱, 等.含曲面介質(zhì)結(jié)構(gòu)復(fù)雜目標電磁散射計算的射線追蹤方法[J]. 電波科學(xué)學(xué)報, 2016,31(3): 546-552.

ZHANG L, YAO Z G, DONG C Z, et al. Ray tracing method for electromagnetic scattering computation from complex target with curved surface dielectric structure[J]. Chinese journal of radio science, 2016, 31(3):546-552. (in Chinese)

[8] 覃璐瑤, 魏兵, 邊成, 等. 一種基于解析積分的TDPO算法及其在散射問題中的應(yīng)用[J]. 電波科學(xué)學(xué)報, 2016,31(3): 508-511.

QIN L Y, WEI B, BIAN C, et al. A TDPO algorithm based on analytic integral and its application in scattering problems[J]. Chinese journal of radio science, 2016, 31(3): 508-511. (in Chinese)

[9] 阮穎錚. 雷達截面與隱身技術(shù)[M]. 北京: 國防工業(yè)出版社, 1998.

[10] 丁建軍, 陳磊, 劉志偉, 等. 基于時域彈跳射線法分析電大尺寸目標的散射[J]. 系統(tǒng)工程與電子技術(shù), 2010,32(9): 1846-1849.

DING J J, CHEN L, LIU Z W, et al. Analysis of electromagnetic scattering of electrically large objects with time domain shooting and bouncing rays method[J]. Systems engineering and electronics, 2010, 32(9): 1846-1849. (in Chinese)

[11] 丁建軍, 劉志偉, 徐侃, 等. 基于高頻方法分析電大尺寸目標的散射,系統(tǒng)工程與電子技術(shù), 2010, 32(11): 2309-2312.

DING J J, LIU Z W, XU K, et al. Application of high frequency methods for scattering analysis of electrically large objects[J]. Systems engineering and electronics, 2010, 32(11): 2309-2312. (in Chinese)