王 永,林中朝,張 玉,梁昌洪

(西安電子科技大學(xué)天線與微波技術(shù)重點(diǎn)實(shí)驗(yàn)室,陜西西安 710071)

并行矩量法分析波導(dǎo)縫隙天線的輻射特性

王 永,林中朝,張 玉,梁昌洪

(西安電子科技大學(xué)天線與微波技術(shù)重點(diǎn)實(shí)驗(yàn)室,陜西西安 710071)

采用并行矩量法對(duì)工程中遇到的Ka波段波導(dǎo)縫隙天線進(jìn)行仿真計(jì)算.矩量法中難以實(shí)現(xiàn)端口匹配,導(dǎo)致其難以有效地處理行波天線問題.為了有效地解決這一難題,采用等效原理和模式匹配法在端口面上構(gòu)造出波端口模型,并進(jìn)一步將波導(dǎo)與波端口分離構(gòu)造出區(qū)域分解模型,將波導(dǎo)與波端口的互作用等效為兩個(gè)區(qū)域的外加激勵(lì),從而使用迭代方法對(duì)此類問題進(jìn)行求解.同時(shí),并行計(jì)算技術(shù)的引入大大提高了矩量法求解電大尺寸問題的能力.數(shù)值仿真結(jié)果與實(shí)測(cè)結(jié)果的對(duì)比表明,波端口模型的使用可以有效地解決波導(dǎo)縫隙天線這一類電磁仿真難題.

矩量法;并行算法;波導(dǎo)縫隙天線;波端口

近年來(lái),機(jī)載平臺(tái)的發(fā)展對(duì)機(jī)載天線的要求越來(lái)越高.大功率、高增益是現(xiàn)代工程應(yīng)用中機(jī)載天線的普遍趨勢(shì).目前,只能采用波導(dǎo)縫隙天線陣構(gòu)建大功率天線.為了減少機(jī)載天線的設(shè)計(jì)、測(cè)試成本,前期的電磁仿真設(shè)計(jì)越來(lái)越顯得重要.然而,機(jī)載波導(dǎo)縫隙天線由于要求高增益,所以其電尺寸一般都比較大.常規(guī)的電磁仿真軟件(如HFSS,FEKO等)都難以在短期內(nèi)計(jì)算出結(jié)果.Ansoft公司的HFSS軟件是基于有限元法的電磁仿真軟件,由于內(nèi)存和計(jì)算時(shí)間的限制,Ka波段的波導(dǎo)縫隙陣通常無(wú)法仿真計(jì)算.FEKO是基于矩量法的電磁仿真軟件,其通過out-of-core[1]技術(shù)解決了內(nèi)存不足的問題,但out-of-core將硬盤當(dāng)內(nèi)存使用的做法大大降低了計(jì)算速度,一個(gè)Ka波段波導(dǎo)縫隙陣的計(jì)算通常需要幾天時(shí)間.這迫使人們必須找到一種新的可實(shí)現(xiàn)的數(shù)值仿真方法.

筆者采用計(jì)算精度較高的矩量法作為數(shù)值仿真方法[2].通常矩量法分析波導(dǎo)縫隙天線存在3個(gè)難點(diǎn):矩量法(Method of Moments,Mo M)是一種嚴(yán)格計(jì)算互耦的方法,與時(shí)域有限差分法(Finite Difference Time Domain,FDTD)和有限元法(Finite Element Method,FEM)等方法不同,矩量法中沒有理想吸收邊界的概念.這就使得矩量法中難以模擬出匹配負(fù)載,無(wú)法實(shí)現(xiàn)行波天線的計(jì)算.波導(dǎo)縫隙天線的縫隙形狀受到很大限制.傳統(tǒng)矩量法分析縫隙天線時(shí),都是在縫隙開口處使用等效原理,即在開口面上放置等效磁流作為等效面內(nèi)外耦合的中介[3].波導(dǎo)縫隙被等效為一個(gè)腔體,這就需要導(dǎo)出腔體格林函數(shù).但是,當(dāng)縫隙形狀不規(guī)則時(shí),其格林函數(shù)很難找到矩量法的計(jì)算量非常大,而Ka波段的波導(dǎo)縫隙天線電尺寸通常也比較大,內(nèi)存、CPU等計(jì)算資源束縛著矩量法的使用.對(duì)于第個(gè)問題(即波端口的建模),德國(guó)不萊梅大學(xué)的Bunger等[4]在2000年時(shí)給出解決方法.對(duì)于第個(gè)問題,筆者采用RWG基函數(shù)[5]將波導(dǎo)縫隙也剖分為電流面片,即在縫隙處不用等效原理,依然使用自由空間的格林函數(shù)進(jìn)行全波仿真,避免了求解腔體格林函數(shù).對(duì)于第個(gè)問題,筆者采用并行矩量法[6],在上海超級(jí)計(jì)算中心的“魔方”超級(jí)計(jì)算機(jī)上使用大量CPU核實(shí)現(xiàn)了核內(nèi)(in-core)求解,大大地加快了計(jì)算速度.同時(shí)作為對(duì)Bunger等方法的改進(jìn),筆者提出的方法可將波導(dǎo)與波端口進(jìn)行區(qū)域分解,從而可以采用基于區(qū)域分解的迭代法求解,使得可計(jì)算的問題規(guī)模進(jìn)一步擴(kuò)大.改進(jìn)后的迭代法可以使波端口建模與波導(dǎo)上的基函數(shù)無(wú)關(guān),易于實(shí)現(xiàn).

1 矩量法中波端口的建模

由于矩量法中不存在理想吸收邊界,所以無(wú)反射的波端口的建模對(duì)波導(dǎo)饋電和行波天線的計(jì)算至關(guān)重要.波端口建模是基于孔徑耦合的等效原理[7]和模式匹配法[8-9].圖1給出了一個(gè)二端口波導(dǎo)的等效模型.

圖1 二端口波導(dǎo)等效模型

為了模擬波導(dǎo)無(wú)反射條件(即接匹配負(fù)載),可假設(shè)波端口外部接了一個(gè)半無(wú)限長(zhǎng)波導(dǎo).因?yàn)橥獠拷恿税霟o(wú)限長(zhǎng)波導(dǎo),則入射到波端口處的場(chǎng)將被導(dǎo)到無(wú)窮遠(yuǎn)處而不會(huì)被反射回來(lái).將原波導(dǎo)內(nèi)區(qū)域記為區(qū)域Ⅰ,而假想的無(wú)限長(zhǎng)波導(dǎo)區(qū)域記為區(qū)域Ⅱ.這樣假設(shè)以后,原問題變?yōu)榘霟o(wú)限長(zhǎng)波導(dǎo)(區(qū)域Ⅱ)與原波導(dǎo)(區(qū)域Ⅰ)在波端口處的孔徑耦合問題.

根據(jù)孔徑耦合等效原理,首先假設(shè)波端口處被理想導(dǎo)體封閉,然后在理想導(dǎo)體面內(nèi)側(cè)放置磁流M.因?yàn)椴ǘ丝诿媸羌傧朊?實(shí)際中并不存在,所以必須保證波端口內(nèi)外兩側(cè)電磁場(chǎng)切向連續(xù).因而需要在其外側(cè)放置磁流-M,記波端口向內(nèi)側(cè)的法向量為.波端口內(nèi)外面上邊界條件為

由于波端口被理想導(dǎo)體封閉,則波端口內(nèi)側(cè)面上還應(yīng)該存在電流Jd.波導(dǎo)內(nèi)壁上也應(yīng)該存在電流,記為Jc.這樣以后,波導(dǎo)內(nèi)問題就可以使用自由空間格林函數(shù),而與波導(dǎo)形狀無(wú)關(guān).對(duì)于區(qū)域Ⅱ問題,因?yàn)槭前霟o(wú)限長(zhǎng)波導(dǎo),所以不能使用自由空間的格林函數(shù),只能在波端口外側(cè)面上放置磁流-M,通過波導(dǎo)內(nèi)的腔體格林函數(shù)求解.

1.1 波導(dǎo)內(nèi)區(qū)域積分方程的建立

由圖1可以看出,波導(dǎo)內(nèi)的電磁場(chǎng)由3個(gè)源產(chǎn)生:波端口上的電流Jd、磁流M及波導(dǎo)內(nèi)壁上的電流Jc.其表達(dá)式為

其中,Z0為波阻抗.L和K算子分別定義為

1.2 模式匹配法建立內(nèi)外耦合積分方程

端口外側(cè)半無(wú)限長(zhǎng)波導(dǎo)內(nèi)的場(chǎng)由激勵(lì)(即入射波Ei,Hi)和波端口外表面上的磁流-M產(chǎn)生的(EⅡ, HⅡ)疊加而成.因?yàn)橥鈧?cè)為半無(wú)限長(zhǎng)波導(dǎo),所以無(wú)法數(shù)值計(jì)算.此時(shí)可采用模式匹配法解析計(jì)算外側(cè)場(chǎng).波導(dǎo)不同的模式即是波導(dǎo)格林函數(shù)的特征向量,模式匹配法也即等價(jià)于采用腔體格林函數(shù)計(jì)算半無(wú)限長(zhǎng)波導(dǎo)內(nèi)的場(chǎng).因?yàn)椴▽?dǎo)的各特征模式已經(jīng)天然地滿足了波導(dǎo)壁上的邊界條件,所以可以只在波端口面上放置激勵(lì)源-M,而不用考慮半無(wú)限長(zhǎng)波導(dǎo)壁上的電流分布.

假設(shè)入射場(chǎng)為

其中,ηj為第j個(gè)模式的波阻抗,ej為波導(dǎo)第j個(gè)模式的電場(chǎng)歸一化切向矢量.其歸一化方式為

式(7)的積分區(qū)間為整個(gè)波端口面.

首先,考慮到波端口被理想導(dǎo)體封閉,半無(wú)限長(zhǎng)波導(dǎo)內(nèi)電磁場(chǎng)(EⅡ,HⅡ)由入射波、全反射波和磁流產(chǎn)生的場(chǎng)疊加而成.磁流產(chǎn)生的場(chǎng)可由波導(dǎo)各模式組合而成,即

其中,bi為未知系數(shù).因?yàn)樵诓ǘ丝诿嫔蟴=0,所以,端口面外側(cè)電磁場(chǎng)可寫為

其次,考慮到波端口外側(cè)面邊界條件EⅡ×(-^n)=-M,所以由式(9)可得

由于式(7)波導(dǎo)各模式相互正交,所以將式(10)兩邊同乘以×ei,然后在波端口面上積分,則有

這樣就將未知系數(shù)bi用端口磁流-M表示出來(lái).將式(11)帶入式(9),得到波端口外側(cè)電磁場(chǎng)的最終表達(dá)式為

其中,ej表示以第j個(gè)模式激勵(lì).

前面已經(jīng)證明波端口內(nèi)外兩側(cè)電場(chǎng)切向連續(xù),所以還須保證內(nèi)外兩側(cè)磁場(chǎng)切向連續(xù),即HⅠ.因?yàn)閮?nèi)側(cè)面上邊界條件,所以.再由矢量旋轉(zhuǎn)關(guān)系,可得HⅡ互相垂直).代入式(12),得積分方程

最后,整理得波導(dǎo)問題的積分方程組為

2 基于區(qū)域分解的迭代法求解

式(14)是適合直接求解的積分方程組.為了減少對(duì)原有矩量法程序的修改,下面提出一種基于區(qū)域分解思想的迭代解法[10],該方法可在邏輯上將波導(dǎo)壁與波端口看做分離的兩個(gè)區(qū)域.

針對(duì)單根天線,區(qū)域分解思想是將端口磁流對(duì)波導(dǎo)上電流的作用看成是區(qū)域Ⅰ的外加激勵(lì).同時(shí),將波導(dǎo)電流對(duì)端口磁流的作用也看做是區(qū)域Ⅱ的外加激勵(lì).

設(shè)置初始電磁流都為零,首先給區(qū)域Ⅱ一個(gè)激勵(lì)模式,解出端口磁流;然后將這個(gè)磁流作為區(qū)域Ⅰ的激勵(lì)算出電流分布;再用算出的電流更新區(qū)域Ⅱ的激勵(lì).如此循環(huán)下去,直到電磁流的解穩(wěn)定.需要指出的是,通過數(shù)值實(shí)驗(yàn)證明,這種迭代對(duì)于波導(dǎo)縫隙天線是收斂的.但對(duì)于單純的波導(dǎo)內(nèi)問題(如濾波器),其收斂性難以保證.

此外,對(duì)于波導(dǎo)縫隙陣天線[11],由于不同天線單元之間通常都是物理上分離的,所以也可以進(jìn)一步將各天線單元分解到不同區(qū)域,從而可以使天線陣問題通過先逐根分析每個(gè)單元,再逐根考慮各單元之間的耦合,最后迭代出結(jié)果.即使在相同內(nèi)存下,可求解的問題規(guī)模也能夠進(jìn)一步增大.因?yàn)樘炀€單元物理上就是分離的,所以迭代收斂性較好.

3 數(shù)值算例

3.1 行波波導(dǎo)縫隙天線陣的計(jì)算

為了驗(yàn)證筆者提出算法的正確性,此處計(jì)算了一個(gè)包含10個(gè)單元的X波段波導(dǎo)縫隙陣天線.其中每個(gè)單元為一個(gè)單側(cè)開有10個(gè)縫隙的行波天線.波導(dǎo)單元的尺寸為22.86 mm×10.16 mm×266.58 mm,天線工作頻率為9.375 GHz.各天線單元的饋電幅度與相位列于表1,天線陣的排布方式如圖2所示.

表1 各天線單元的饋電幅度與相位

由于頻率較高,天線的電尺寸較大,所以矩量法計(jì)算時(shí)剖分的網(wǎng)格較多.以1/10波長(zhǎng)尺寸的三角形網(wǎng)格剖分時(shí),共產(chǎn)生95 388個(gè)三角形面片.對(duì)應(yīng)的矩陣方程有143 482個(gè)未知量,需要內(nèi)存306.8 GB.如此規(guī)模的問題在一般PC機(jī)上串行計(jì)算顯然是不現(xiàn)實(shí)的.筆者在上海超級(jí)計(jì)算中心使用1 024個(gè)CPU核進(jìn)行計(jì)算,填充矩陣耗時(shí)39 min,求解矩陣方程耗時(shí)30 min.該天線陣的歸一化遠(yuǎn)場(chǎng)方向圖如圖3所示.

圖2 10單元行波陣天線模型

圖3 10單元波導(dǎo)縫隙天線陣遠(yuǎn)場(chǎng)輻射方向圖

圖3中實(shí)線為筆者提出方法的數(shù)值計(jì)算結(jié)果,虛線為實(shí)際測(cè)量結(jié)果.可以看出,筆者提出方法的數(shù)值計(jì)算結(jié)果與實(shí)測(cè)結(jié)果吻合良好.尤其在主瓣區(qū)域,幾乎完全重合,從而證明了筆者構(gòu)造的波端口模型在模擬計(jì)算行波天線時(shí)的有效性和準(zhǔn)確性,為計(jì)算更大規(guī)模的波導(dǎo)縫隙天線陣提供了理論支持.3.2 Ka波段波導(dǎo)縫隙天線的分析

此處選取的計(jì)算模型為一根一側(cè)開有104個(gè)縫的二端口波導(dǎo)縫隙天線,如圖4所示.天線為一端口饋電,另一端口匹配的行波天線.在構(gòu)造出波端口模型前,由于缺少理想匹配負(fù)載模型,矩量法難以模擬行波天線.以前的解決方法是在匹配端口上接高損耗的吸收媒質(zhì),將傳輸?shù)狡ヅ涠丝谏系牟〒p耗掉來(lái)模擬匹配.然而吸收媒質(zhì)的吸收效果與頻率、波導(dǎo)尺寸、傳輸波的模式等因素都有關(guān)系,每一個(gè)天線要獨(dú)立設(shè)計(jì)吸收媒質(zhì)的材料參數(shù)和所需要的媒質(zhì)層數(shù),該方法無(wú)法通用.即使設(shè)計(jì)出吸收媒質(zhì),它的吸收效果也難以達(dá)到100%.實(shí)踐證明,即使損耗媒質(zhì)的吸收達(dá)到98%,依然會(huì)有少量的反射波回到波導(dǎo)內(nèi),對(duì)計(jì)算產(chǎn)生影響.

圖4 Ka波段波導(dǎo)縫隙天線模型

選取天線尺寸為7.112 mm×3.556 mm×591.6 mm,天線工作頻率為35 GHz.經(jīng)矩量法剖分后產(chǎn)生132 474個(gè)未知量,需要261 GB內(nèi)存.采用并行矩量法在“魔方”超級(jí)計(jì)算機(jī)中對(duì)該波導(dǎo)縫隙天線的輻射特性進(jìn)行仿真計(jì)算,使用256 CPU核在2 h內(nèi)即可得到計(jì)算結(jié)果.該天線的遠(yuǎn)場(chǎng)方向圖如圖5所示.

圖5中實(shí)線為筆者采用波端口模型的計(jì)算結(jié)果,虛線為采用吸收媒質(zhì)模擬匹配的計(jì)算結(jié)果.虛線的主瓣較窄,是因?yàn)榉较驁D采樣點(diǎn)的緣故.虛線為每1°采一個(gè)點(diǎn),實(shí)線為每0.1°采一個(gè)點(diǎn),所以虛線有跳過的點(diǎn),導(dǎo)致顯示的主瓣較窄.在-45°、45°和105°附近出現(xiàn)的3個(gè)較大副瓣是交叉極化產(chǎn)生的.需要注意的是,在主瓣右側(cè)吸收媒質(zhì)計(jì)算出的結(jié)果比波端口方法多出了一個(gè)小瓣.通過分析,筆者認(rèn)為這是吸收媒質(zhì)無(wú)法做到100%吸收而產(chǎn)生的反射波輻射.這也證明了波端口方法更加適合模擬匹配負(fù)載.

圖5 Ka波段波導(dǎo)縫隙天線輻射方向圖

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

在構(gòu)造的波端口模型基礎(chǔ)上,筆者提出了適合于迭代求解的區(qū)域分解模型計(jì)算波導(dǎo)縫隙天線問題.由于波端口實(shí)現(xiàn)了匹配負(fù)載,使得行波天線的計(jì)算成為可能.并行計(jì)算和迭代方法的引入解決了大規(guī)模波導(dǎo)縫隙天線的計(jì)算難題,配合“魔方”平臺(tái)可以解決機(jī)載天線布局等工程中迫切需要解決的問題.筆者分析Ka波段波導(dǎo)縫隙天線只是解決工程問題的開始,下一步嘗試將天線加載到飛機(jī)平臺(tái)上去,解決更大規(guī)模的電磁仿真問題.

[1]Zhang Y,Sarkar T K.Parallel Solution of Integral Equation Based EM Problems in the Frequency Domain[M]. Hoboken:Wiley-IEEE Press,2009.

[2]Harrington R F.Field Computation by Moment Methods[M].New York:IEEE Press Series on Electromagnetic Wave Theory,1993.

[3]Lyon R W,Sangster A J.Efficient Moment Method Analysis of Radiating Slots in a Thick-walled Rectangular Waveguide[C]//IEE Proceedings H:Microwaves,Optics and Antennas,1981,128(4):197-205.

[4]Bunger R,Arndt F.Moment-method Analysis of Arbitrary 3-D Metallic N-port Waveguide Structures[J].IEEE Transactions on Microwave Theory and Techniques,2000,48(4):531-537.

[5]Rao S M,Wilton D,Glisson A W.Electromagnetic Scattering by Surfaces of Arbitrary Shape[J].IEEE Transactions on Antennas and Propagation,1982,30(3):409-418.

[6]Yan Y,Zhao W X,Zhao X W,et al.Parallel Computation of RCS of Electrically Large Platform with Coatings Modeled with NURBS Surfaces[J].International Journal of Antennas and Propagation,2012,2012:348052.

[7]Harrington R F,Mautz J R.A Generalized Network Formulation for Aperture Problems[J].IEEE Transactions on Antennas and Propagation,1976,24(6):870-873.

[8]Auda H,Harrington R F.A Moment Solution for Waveguide Junction Problems[J].IEEE Transactions on Microwave Theory and Techniques,1983,31(7):515-520.

[9]趙偉,趙永久,路宏敏.矩形波導(dǎo)不連續(xù)性的通用模式匹配法分析[J].西安電子科技大學(xué)學(xué)報(bào),2008,35(5):894-898. Zhao Wei,Zhao Yongjiu,Lu Hongmin.General Analysis Method of Rectangular Waveguide Discontinuities by Mode Matching[J].Journal of Xidian University,2008,35(5):894-898.

[10]Peng Z,Wang X C,Lee J F.Integral Equation Based Domain Decomposition Method for Solving Electromagnetic Wave Scattering from Non-penetrable Objects[J].IEEE Transactions on Antennas and Propagation,2011,59(9):3328-3338.

[11]李建瀛,梁昌洪.波導(dǎo)窄邊縱縫的理論分析及頻掃縫陣[J].西安電子科技大學(xué)學(xué)報(bào),1998,25(5):689-692. Li Jianying,Liang Changhong.Theoretical Analysis of the Longitudial Slot in the Narrow Wall of a Rectangular Waveguide and Frequence Scan Array[J].Journal of Xidian University,1998,25(5):689-692.

(編輯:郭 華)

Analysis of the radiation of slotted waveguide antennas using the parallel Mo M

WANG Yong,LIN Zhongchao,ZHANG Yu,LIANG Changhong
(Science and Technology on Antenna and Microwave Lab.,Xidian Univ.,Xi’an 710071,China)

This paper numerically analyzes the Ka band slotted waveguide antennas which are involved in engineering using parallel method of moments(Mo M).For lack of the match terminator,it is difficult for the Mo M to simulate the traveling wave antenna.To deal with this problem,the equivalent principle on the wave ports and the mode match methods are used to construct the model of wave ports.As an improvement,the waveguide and the wave ports are divided into different domains.The affections between the two domains are assumed as the external excitations of each domain.Then,the iteration method can be used to solve the equation.At the same time,the introduction of the parallel calculation technique can improve the ability of the Mo M to analyze the electrically large targets.The agreement of the measurement result and numerical result proves that the model of wave ports proposed in this paper can help to solve the problem of the simulation of slotted waveguide antenna validly.

method of moments(Mo M);parallel algorithm;slotted waveguide antenna;wave port

TN820

A

1001-2400(2015)05-0075-05

2014-05-20< class="emphasis_bold">網(wǎng)絡(luò)出版時(shí)間:

時(shí)間:2014-12-23

國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)資助項(xiàng)目(2012AA01A308);國(guó)家自然科學(xué)基金資助項(xiàng)目(61301069,61072019);教育部新世紀(jì)優(yōu)秀人才支持計(jì)劃資助項(xiàng)目(NCET-13-0949);陜西省青年科技新星資助項(xiàng)目(2013KJXX-67);中央高校基本科研業(yè)務(wù)費(fèi)重點(diǎn)資助項(xiàng)目(JY10000902002)

王 永(1987-),男,西安電子科技大學(xué)博士研究生,E-mail:yongmom@163.com.

http://www.cnki.net/kcms/detail/61.1076.TN.20141223.0946.013.html

10.3969/j.issn.1001-2400.2015.05.013