肇 格,張 軍,胡杰民

(國(guó)防科技大學(xué)ATR國(guó)防科技重點(diǎn)實(shí)驗(yàn)室,長(zhǎng)沙410073)

1 引 言

低頻方法在雷達(dá)散射截面(RCS)的計(jì)算上有著很大的局限性,特別是對(duì)于電大尺寸具有復(fù)雜結(jié)構(gòu)的目標(biāo),該方法往往是行不通的。而對(duì)于普通的高頻方法,例如物理光學(xué)方法,由于其是建立在很多種近似上的,所以在計(jì)算很多問題時(shí)很難達(dá)到準(zhǔn)確性的要求。射線追蹤與物理光學(xué)(PO)相結(jié)合的方法包含了散射場(chǎng)中多次散射的高階項(xiàng),是計(jì)算具有多次散射效應(yīng)目標(biāo)RCS 的一種有效方法,其有效性已經(jīng)得到了大量的驗(yàn)證[1-3]。但這種方法也只被用在了角反射器的計(jì)算上。

根據(jù)衛(wèi)星目標(biāo)的三維面元模型建立起來的Kd-tree 是一種根據(jù)k 維空間中的點(diǎn)集對(duì)空間進(jìn)行分割的數(shù)據(jù)結(jié)構(gòu),常用于范圍查找和最近鄰查找等,是一種特殊的二叉空間分割樹[4-5]。文獻(xiàn)[5]中作者將其用于射線追蹤方法使得對(duì)于射線的追蹤效率大大提高。本文將Kd-tree 方法應(yīng)用于射線追蹤,仿真結(jié)果表明該方法在計(jì)算效率方面達(dá)到了良好的效果,結(jié)合PO 方法就可以得到給定模型的RCS 預(yù)估值。計(jì)算結(jié)果表明,結(jié)合了Kd-tree 和射線追蹤的方法可以有效地計(jì)算復(fù)雜電大尺寸空間目標(biāo)的RCS預(yù)估值。

2 衛(wèi)星目標(biāo)模型的剖分原則

衛(wèi)星目標(biāo)模型建好之后,將其剖分成小的三角面元,具體的剖分應(yīng)當(dāng)遵循以下原則:假定每個(gè)小面元都只有被照射和未被照射兩種情況。因此,為了比較準(zhǔn)確地反映多次散射的情況,剖分時(shí)每個(gè)面元都不能太大,以使最小結(jié)構(gòu)都能被面元模型很好地描述。在滿足上述條件下,為了提高計(jì)算精度,這里采用九點(diǎn)插值方法對(duì)每一面元的物理光學(xué)積分進(jìn)行計(jì)算。

3 Kd-tree 結(jié)構(gòu)的建立

Kd-tree 結(jié)構(gòu)是自上而下以一種遞歸的方式構(gòu)造的,其中根節(jié)點(diǎn)對(duì)應(yīng)于包含了目標(biāo)所有面元的長(zhǎng)方體包圍盒。然后長(zhǎng)方體包圍盒按一個(gè)軸方向被分割平面分成兩個(gè)子長(zhǎng)方體,組織成一個(gè)二叉樹,即一個(gè)根節(jié)點(diǎn)被分成了兩個(gè)內(nèi)節(jié)點(diǎn)。剖分平面是選擇使射線與面元相交概率最大的平面,且該平面可以是垂直于x 、y 、z 軸的任何一個(gè)平面。若某一子長(zhǎng)方體包圍盒中所含的面元數(shù)大于給定的值,且Kd -tree 的深度小于給定的深度,則需對(duì)該子包圍盒進(jìn)一步剖分。上述剖分過程遞歸地進(jìn)行,直至Kd -tree 的每一個(gè)葉子節(jié)點(diǎn)所對(duì)應(yīng)的子空間所包含的面元數(shù)小于給定的值或該Kd-tree 的深度大于給定的深度為止。另外,如果分割過程中面元與分割平面相交,則使面元分別屬于兩個(gè)相應(yīng)的節(jié)點(diǎn)。Kd-tree建立之后,還需使每一葉子節(jié)點(diǎn)對(duì)應(yīng)的包圍盒的6個(gè)面與其他內(nèi)節(jié)點(diǎn),葉子節(jié)點(diǎn)或空節(jié)點(diǎn)建立指向關(guān)系,具體細(xì)節(jié)可以參見文獻(xiàn)[4] 。Kd-tree 的建立可以大大減小反射線追蹤時(shí)從Kd-tree 自上而下的追蹤過程。二維的Kd-tree 如圖1 所示。

圖1 二維Kd-tree 模型Fig.1 The two-dimensional model of Kd-tree

4 PO 積分與路徑的射線追蹤

4.1 應(yīng)用物理光學(xué)積分計(jì)算每一面元一階等效電流對(duì)于遠(yuǎn)區(qū)散射場(chǎng)的貢獻(xiàn)

本文所介紹的方法是以面元為單位進(jìn)行的,即計(jì)算出每一面元各階等效電流對(duì)于遠(yuǎn)區(qū)散射場(chǎng)的貢獻(xiàn),然后將所有面元的貢獻(xiàn)相加即可得總的遠(yuǎn)區(qū)散射場(chǎng)。其中,每一面元對(duì)于遠(yuǎn)區(qū)散射場(chǎng)的貢獻(xiàn)可用物理光學(xué)積分表示為

其中,Hsm為面元m 對(duì)于遠(yuǎn)區(qū)散射場(chǎng)的貢獻(xiàn),Jm1、Jm2等為面元m 上的一階、二階和更高階電流。

應(yīng)用切平面近似,面元m 上的一階等效電流可表示為

其中, m 為面元m 的法向矢量, Hi為面元m 處的磁場(chǎng)強(qiáng)度。

4.2 根據(jù)入射波判斷初次入射明暗面

該步驟需對(duì)所有面元一一進(jìn)行判斷,判斷其是否被入射波直接照射到。判斷某一面元m 是否被入射波照射到的方法如下。

設(shè)入射波矢量為 i,面元m 的法向矢量為 n ,則m 面元被入射波照亮應(yīng)滿足以下兩個(gè)條件:

(1)入射波矢量為 i 與面元m 的法向矢量為 m(mx,my,mz)滿足: i × m ≤0;

(2)入射波照射到面元m 的過程未被其他面元遮擋。

其中,條件2 的判斷方法如下:

設(shè)面元中心點(diǎn)坐標(biāo)rm(xm, ym,zm),入射波矢量 i(ix,iy,iz),該面元上的入射線方程為

將該直線方程與其他面元逐一進(jìn)行判斷,如果該直線方程與其他所有面元無交點(diǎn),則面元m 未被其他面元遮擋,即面元m 被入射波直接照射到,否則面元m 在入射波照射過程被其他面元遮擋,未被入射波直接照射到。

以面元n 為例,判斷式(3)直線與面元n 是否有交點(diǎn)的方法如下。

設(shè)面元n 的中心點(diǎn)坐標(biāo)為rn(xn, yn,zn),面元法向矢量為 n(nx, ny,nz),則該面元所在的平面可表示為

聯(lián)立式(3)與式(4)可求得式(3)直線與面元n所在平面的交點(diǎn)P(x0,y0,z0)。這里:

其中:

t0=[ nx(xn-xm)+ny(yn-ym)+nz(zn-zm)]/[ nx·ix +ny·i y+nz·iz] 。

因?yàn)槲覀兣袛嗟氖巧渚€到達(dá)面元m 之前的情況,所以t 0應(yīng)該小于0。此時(shí),如果該交點(diǎn)P(x0,y0,z0)在面元n 內(nèi),則式(3)與面元n 有交點(diǎn),否則式(3)與面元n 無交點(diǎn)。判斷一點(diǎn)P(x0, y0,z0)是否在三角形面元n 內(nèi)的方法如下:

3 個(gè)頂點(diǎn)N1、N2、N3與面元n 的法向矢量n (nx,ny,nz)成右手螺旋法則的順序排列,如圖2 所示,則P 點(diǎn)位于面元n 內(nèi)部的充要條件是式(6)中的3 個(gè)式同時(shí)滿足。

圖2 射線與面元相交判斷Fig.2 Intersect judge of ray and surface

4.3 面元處二階等效電流對(duì)遠(yuǎn)區(qū)散射場(chǎng)的貢獻(xiàn)的計(jì)算

仍以經(jīng)面元n 反射的反射線照亮面元m 為例。面元n 處的入射場(chǎng)可表示為

則經(jīng)面元n 反射的反射場(chǎng)可表示為

其中,Hi// = e//·Hi,Hi⊥= e ⊥·Hi, e ⊥= i × n, e// = i× e⊥, er//= k × e⊥, Γ//和Γ⊥分別為平行極化和垂直極化的反射系數(shù)。

5 仿真結(jié)果

算例1:三面角反射器是驗(yàn)證高頻多次彈跳散射的標(biāo)準(zhǔn)[6]。這里所取的三面角反射器中直角三角形的直角邊長(zhǎng)為1 m,如圖3 所示。θ在φ=45°的平面從0°到90°變化,入射波頻率為6 GHz,VV 極化。

圖3 三面角反射器Fig.3 Three sides of corner reflectors

圖4 給出了三面角反射器的后向散射系數(shù),并將其與MoM 方法的結(jié)果進(jìn)行了比較。從圖中可以看出,本文所提出的方法與MoM 方法所得的結(jié)果符合得很好。在計(jì)算效率方面,MoM 方法耗時(shí)超過24 h,而本文所提出的方法計(jì)算所有90 個(gè)角度用時(shí)5 min。

圖4 本文方法與MoM 方法結(jié)果比較Fig.4 Comparison between the proposed method and the MoM

算例2:將該方法應(yīng)用于衛(wèi)星目標(biāo)的計(jì)算。衛(wèi)星模型如圖5 所示,在模型上面的平臺(tái)上為一天線面結(jié)構(gòu),該天線面結(jié)構(gòu)中有圓柱和方體結(jié)構(gòu),模型被剖分為29 569個(gè)面元。入射波頻率為3 GHz,VV 極化。

圖5 衛(wèi)星模型Fig.5 Satellite model

圖6 給出了本文方法所得結(jié)果與實(shí)驗(yàn)結(jié)果和單純PO 方法所得結(jié)果的比較。從圖中我們可以看出,單純的PO 方法結(jié)果與實(shí)驗(yàn)結(jié)果相差很大,這表明對(duì)于衛(wèi)星天線面這種復(fù)雜結(jié)構(gòu),各部分之間多次散射的效應(yīng)十分明顯,而單純的PO 方法未能包含這種效應(yīng)。本文的改進(jìn)方法則能有效涵蓋這種多次散射效應(yīng),所以與實(shí)驗(yàn)結(jié)果符合得比較好。

圖6 衛(wèi)星RCS 計(jì)算結(jié)果比較Fig.6 Comparison of satellite RCS calculation

6 結(jié) 論

本文提供了一種將PO 方法與Kd-tree 空間分割方法及射線追蹤相結(jié)合的改進(jìn)方法。該方法因?yàn)橛行У赜?jì)入了目標(biāo)各部分之間的多次散射效應(yīng),且為一種高頻方法,所以對(duì)于計(jì)算電大尺寸具有較強(qiáng)耦合效應(yīng)的復(fù)雜目標(biāo)是非常行之有效的。為了方便起見,本文沒有考慮邊緣繞射效應(yīng)所產(chǎn)生的影響,如果用物理繞射理論(PTD)對(duì)之加以修正,可能會(huì)得到更令人滿意的結(jié)果。值得一提的是,其他大多數(shù)方法的計(jì)算量會(huì)隨著入射波頻率的增大而顯著增加,而本文方法的計(jì)算量則不會(huì)受到入射波頻率的影響,這使得該方法在高頻散射計(jì)算時(shí)更能顯示出更多的優(yōu)勢(shì)。

[1] Knott E F.RCS reduction of dihedral corners[ J] .IEEE Transactions on Antennas and Propagation,1977, 25(3):406-409.

[2] Ross R A.Application of geometrical diffraction theory to reflex scattering centers[C]//Proceedings of 1968 IEEE International Symposium on Antennas and Propagation.Monterey,CA:IEEE,1968:94-99.

[3] Anderson W C.Consequence of nonorthogonality on the scattering properties of dihedral reflectors[ J] .IEEE Transactions on Antennas and Propagation, 1987,35(10):1154-1159.

[4] Havran V.Heuristic Ray Shooting Algorithm[D] .Prague:University Czech Technical, 2000.

[5] Tao Y B, Lin H, Bao H J.Kd-tree based fast ray tracing for RCS prediction[ J] .Progress in Electromagnetics Research, 2008, 81:329-341.

[6] Baldauf J, Lee S W, Lin L, et al.High frequency scattering from trihedral corner reflectors and other benchmark targets:SBR vs.experiments[ J] .IEEE Transactions on Antennas and Propagation, 1991,39(9):1345-1351.

[7] 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.