高寶春 彭 戈 胡一峰

中國(guó)艦船研究設(shè)計(jì)中心，上海 201108

基于矩量法的線天線電磁散射和電磁輻射分析

高寶春 彭 戈 胡一峰

中國(guó)艦船研究設(shè)計(jì)中心，上海 201108

針對(duì)線電流模型分析半徑／波長(zhǎng)較大的線天線的缺點(diǎn)，首先利用解析法描述線天線的結(jié)構(gòu)，將線天線模擬為細(xì)帶模型，線面連接問題轉(zhuǎn)化為面面連接問題，選取RWG函數(shù)對(duì)天線模型進(jìn)行剖分，采用矩量法分析不同半徑下偶極子天線的電流分布，輻射場(chǎng)，方向圖，對(duì)散射模式和輻射模式的情況進(jìn)行了比較。隨后計(jì)算了有限平板上10m單極子天線的天線特性。計(jì)算結(jié)果表明，當(dāng)半徑／波長(zhǎng)較小時(shí)，電流基本沿天線軸向流動(dòng)，線電流模型適用于分析線天線的問題，半徑／波長(zhǎng)較大時(shí)，橫向電流不可忽略，線電流模型不適用于分析線天線的問題。

矩量法；RWG基函數(shù)；線天線；電磁散射；電磁輻射

1 引言

線天線是最常用的天線形式之一，大到飛機(jī)、艦船，小到手機(jī)等移動(dòng)通信設(shè)備都能見到線天線的使用，因此準(zhǔn)確的分析線天線的電磁特性具有重要的實(shí)際意義。基于電場(chǎng)積分方程的矩量法能夠準(zhǔn)確的分析線天線的電磁特性［1］。

常見的線天線形式主要有偶極子天線和單級(jí)子天線兩種，本文將分別對(duì)這2種天線形式進(jìn)行研究。在分析線天線的問題時(shí)，通常將線天線作為一維分段模型加以研究，假設(shè)電流只沿天線的軸向方向流動(dòng)，在橫向方向上沒有電流分布。在實(shí)際情況中，當(dāng)天線的半徑較小時(shí)，橫向電流可以忽略不計(jì)，但當(dāng)半徑較大時(shí)，橫向電流將不可忽略不計(jì)，此時(shí)，電流分布已不能近似為沿軸向方向的流動(dòng)［2］。

在用矩量法分析線面連接的問題時(shí)，通常需要構(gòu)造三類基函數(shù)，分別是線元之間，面元之間，線元與面元之間［3］。三類基函數(shù)的構(gòu)造給分析和計(jì)算帶來了較大的麻煩。部分文獻(xiàn)［4］給出了一種用解析法描述天線結(jié)構(gòu)的例子并受到了較大的關(guān)注［5-8］。本文采用這種方法，建立金屬細(xì)帶模型模擬線天線的結(jié)構(gòu)，導(dǎo)體面和金屬細(xì)帶模型均通過RWG邊元對(duì)其剖分，將線面連接問題轉(zhuǎn)化為面面連接問題，給分析和計(jì)算帶來了較大的方便。在此基礎(chǔ)之上，分析了偶極子天線的散射特性與輻射特性，并對(duì)較少受關(guān)注的近場(chǎng)場(chǎng)強(qiáng)進(jìn)行了計(jì)算。隨后，對(duì)有限平面上的單根大功率單極子天線進(jìn)行了仿真，給出了單極子天線的輸入阻抗，并引入了組合單極子模型進(jìn)行比較，分析兩種模型下的近場(chǎng)場(chǎng)強(qiáng)與方向性圖。

2 理論分析

2.1 電磁散射理論分析

對(duì)于電磁散射問題，假設(shè)有一入射場(chǎng)Ei照射在某理想導(dǎo)體表面產(chǎn)生感應(yīng)電流J。Es為J在自由空間的輻射場(chǎng)，n為理想導(dǎo)體表面外法線單位矢量。Es可由下式表示：

根據(jù)理想導(dǎo)體的邊界條件 n×（Ei＋Es）＝0，得到：

采用矩量法計(jì)算電場(chǎng)積分方程時(shí)，導(dǎo)體表面采用平面三角形面元進(jìn)行剖分，基函數(shù)選擇如圖1所示的RWG基函數(shù)［9］，導(dǎo)體表面電流J可表示為：

其中，N為未知量總數(shù)；In為待求電流密度系數(shù)。用伽略金法檢驗(yàn)后，按照矩量法的標(biāo)準(zhǔn)形式得到矩陣方程：

當(dāng)確定［Zmn］和［Vm］后，可通過求解矩陣方程得到［Im］。結(jié)合電場(chǎng)積分方程，可得出Zmn和Vm的表達(dá)式如下［9］：

式（4）～式（8）中，r為場(chǎng)點(diǎn)坐標(biāo)；r′為源點(diǎn)坐標(biāo)；m、n對(duì)應(yīng)于2個(gè)邊元；r為邊元m的2個(gè)三角T的中心點(diǎn)，分別是邊元m的2個(gè)三角T的自由頂點(diǎn)到中心點(diǎn)的矢量。

2.2 電磁輻射理論分析

當(dāng)我們考慮線天線的輻射模式時(shí)，激勵(lì)源將由入射波變?yōu)殡妷涸?，需要將饋電模型引入天線結(jié)構(gòu)以便考慮電壓源的影響。在實(shí)際工程中，有多種方式設(shè)置激勵(lì)源，其中較為常用的是由傳輸線通過兩個(gè)靠近的端子饋電。對(duì)于這種激勵(lì)模式，可以很好的用δ函數(shù)縫隙電壓源模型進(jìn)行模擬［10］。當(dāng)用縫隙電壓源激勵(lì)時(shí)，假設(shè)縫隙寬度可以忽略不計(jì)，則縫隙內(nèi)的電場(chǎng)可以用δ函數(shù)近似表示為：

式中，V0為外加加壓；n為電場(chǎng)方向。

將間隙與邊元結(jié)構(gòu)的內(nèi)部邊m聯(lián)系起來，對(duì)應(yīng)此內(nèi)部邊只有一個(gè)RWG邊元m。除此邊元以外，其他地方，入射場(chǎng)均為0。對(duì)于邊元m上的激勵(lì)電壓，可通過下式給出：

而后，可通過式（6）計(jì)算得出［Im］。

2.3 細(xì)帶模型建模

根據(jù)相關(guān)文獻(xiàn)論述［11］，具有電小尺寸的圓柱線天線，當(dāng)天線半徑r遠(yuǎn)小于波長(zhǎng)λ時(shí)，其電磁特性可用一條金屬細(xì)帶等效。天線半徑與金屬細(xì)帶寬度之間的關(guān)系為r＝0.25 s。其中，r為天線的半徑；s為金屬細(xì)帶的寬度。進(jìn)行等效后，可通過建立一個(gè)金屬細(xì)帶模型代替圓柱線天線，并進(jìn)行三角形網(wǎng)格剖分，用RWG函數(shù)描述金屬細(xì)帶上的電流分布?？梢宰C明，每個(gè)邊元的基函數(shù)可近似成長(zhǎng)度為的有限小電偶極子。式中，rc±為邊元的兩個(gè)三角T±的中心點(diǎn)，RWG邊元可以保證沿金屬細(xì)帶有均勻的軸向電流，滿足細(xì)線理論的假設(shè)［10］。

2.4 電場(chǎng)磁場(chǎng)計(jì)算

建立金屬細(xì)帶模型后，要計(jì)算空間某點(diǎn)處的電場(chǎng)、磁場(chǎng)，可以通過計(jì)算每個(gè)邊元對(duì)應(yīng)的有效電偶極子在空間某點(diǎn)處的電場(chǎng)、磁場(chǎng)并進(jìn)行求和，得出該點(diǎn)處的電場(chǎng)、磁場(chǎng)大小。位于原點(diǎn)的長(zhǎng)度為d的有限小電偶極子在空間某點(diǎn)處R的輻射磁場(chǎng)和電場(chǎng)以矢量形式可以表示為［12］：

式中，r、θ、φ分別為球坐標(biāo)系中的單位坐標(biāo)。上述方程是精確表達(dá)式，沒有進(jìn)行遠(yuǎn)場(chǎng)近似，可以計(jì)算離開偶極子外任意點(diǎn)處的電場(chǎng)、磁場(chǎng)。但考慮到模型中等效的有效小電偶極子具有一定的長(zhǎng)度，實(shí)際能計(jì)算的范圍應(yīng)該大于電偶極子長(zhǎng)度的量級(jí)。如果要計(jì)算非常近的距離的電場(chǎng)、磁場(chǎng)強(qiáng)度，則需要在模型網(wǎng)格剖分時(shí)劃分更細(xì)的網(wǎng)格，以滿足計(jì)算精度的需求。

考慮各個(gè)邊元的合成效果，則空間某點(diǎn)R處的電場(chǎng)和磁場(chǎng)為：

3 數(shù)值計(jì)算與分析

為了驗(yàn)證上述算法，分別進(jìn)行了不同情況下的偶極子天線的散射仿真，輻射仿真和單極子天線的輻射仿真。首先在MATLAB環(huán)境下生成均勻網(wǎng)格化的天線模型?？梢酝ㄟ^編寫程序，定義天線和平板的長(zhǎng)、寬、剖分密度，生成網(wǎng)格時(shí)，使用德洛內(nèi)三角化［10］。實(shí)際工作中可以很方便的生成任意長(zhǎng)、任意寬的三角化天線網(wǎng)格，通過控制剖分密度達(dá)到控制未知量的數(shù)目，滿足計(jì)算精度和計(jì)算速度的要求，具有廣泛的適用性。

3.1 偶極子天線散射模型仿真

仿真頻率設(shè)置為50 MHz，偶極子半徑分別設(shè)置為波長(zhǎng)的 10－5、10－4、10－3、10－2， 等效細(xì)帶寬度分別為波長(zhǎng)的4 × 10－5、4 × 10－4、4 × 10－3、4 × 10－2，細(xì)帶長(zhǎng)度為3 m。入射信號(hào)是垂直于細(xì)帶模型平面的平面波，該平面波只有沿細(xì)帶軸向方向的電場(chǎng)分量，大小為1 V／m。分別得到細(xì)帶模型的表面電流強(qiáng)度分布，輻射電場(chǎng)強(qiáng)度和方向性圖如圖2所示。

在圖2中可以看出，隨著天線半徑的增加，表面電流強(qiáng)度分布逐漸減小，這是因?yàn)樵谄涞刃Ы饘偌?xì)帶模型中，隨著天線半徑的增加，金屬細(xì)帶的寬度逐漸增加，表面電流在橫向方向上的分量逐步增大。從圖3和圖4中可以看出天線半徑的變化對(duì)于電場(chǎng)和方向性的變化影響較小。

3.2 偶極子天線輻射仿真

偶極子天線模型采用3.1節(jié)中所采用的模型，激勵(lì)信號(hào)由入射平面波改為饋電電壓為1V的電壓信號(hào)，饋電點(diǎn)位于天線中點(diǎn)，分別得到細(xì)帶模型的表面電流強(qiáng)度分布，輻射電場(chǎng)強(qiáng)度和方向性圖如圖5所示。

從圖5中可以看出，當(dāng)天線半徑較小時(shí)，細(xì)帶模型上的電流分布基本呈現(xiàn)正弦分布，這是因?yàn)樘炀€半徑較小時(shí)，天線上的電流基本沿軸向流動(dòng)。隨著天線半徑的增大，當(dāng)半徑與波長(zhǎng)的比值為0.01時(shí)，細(xì)帶模型上的電流分布已不在呈現(xiàn)正弦分布，究其原因是此時(shí)的電流分布已不能近似為沿軸向分布，橫向方向上的電流分布增大所致。這種情況與前面表述的偶極子天線電磁散射的情況有所不同，在電磁散射情況下，當(dāng)半徑與波長(zhǎng)的比值等于0.01時(shí)，雖然電流的幅度有較明顯的變化，但仍然呈現(xiàn)正弦分布。這兩者不同的原因在于，在散射環(huán)境下，整個(gè)金屬細(xì)帶上的電流分布更加平均，對(duì)細(xì)帶寬度的要求比輻射情況下更低。所以在進(jìn)行天線輻射仿真時(shí)，應(yīng)注意金屬帶寬度與長(zhǎng)度的比值不能過小。

在圖6、圖7中，我們舍棄半徑與波長(zhǎng)的比值等于0.01的情況，分別計(jì)算了輻射電場(chǎng)從近場(chǎng)到遠(yuǎn)場(chǎng)的變化以及方向性圖。在圖7中，引入在FEKO中仿真進(jìn)行的偶極子天線的方向性圖作為比較依據(jù)?？梢钥闯?，在三種情況下，電場(chǎng)分布數(shù)據(jù)的吻合度都很高，只有細(xì)微的差別。方向性圖的吻合度也很高，并且與在FEKO中的仿真結(jié)果極為相似，驗(yàn)證了這種算法的準(zhǔn)確性。

3.3 單極子天線輻射仿真

10 m鞭天線是艦船上常用的通信天線，屬于單級(jí)子天線，通常有2根天線分別布置于上層甲板的左右舷。分析10 m單級(jí)子天線的本身特性和組合布置情況下的特性具有重要的實(shí)際意義。本文分別仿真了兩種情況下的單極子天線，第一種情況是置于20×20 m的平板中心的一根長(zhǎng)為10 m的單極子天線，底部饋電。第二種情況是在20×20 m的平板上分別放置兩個(gè)單極子天線，位置分別位于平面中心的左側(cè)7.5 m處和右側(cè)7.5 m處，對(duì)左側(cè)的單極子天線底部饋電，右側(cè)的單極子天線不饋電，饋電電壓設(shè)置為300 V，模型如圖8所示。

常用的10 m單極子天線的工作頻率為5～28 MHz，計(jì)算10 m單極子天線的輸入阻抗如圖9所示。

表1比較了兩種模型在7.5 MHz時(shí)的輸入阻抗和增益值，可以看出，組合單極子模型的增益更大，輸入阻抗的實(shí)部增大。圖10所示為兩種模型在7.5 MHz時(shí)的輻射電場(chǎng)從近場(chǎng)到遠(yuǎn)場(chǎng)的變化關(guān)系，可以看出兩者的輻射電場(chǎng)有細(xì)微的差距，隨著距離的增加，兩條曲線逐漸重合。圖11給出了兩種模型的方向性圖，從圖中可以明顯的看出組合單極子的增益更大，方向性更強(qiáng)。

表1 7.5 MHz時(shí)2個(gè)模型的比較Tab.1 Comparison of the two models at 7.5 MHz

4 結(jié)論

基于矩量法，采用RWG基函數(shù)對(duì)3 m長(zhǎng)偶極子天線和有限平板上的10 m長(zhǎng)單極子天線的特性進(jìn)行了分析。分析偶極子天線的半徑發(fā)生變化時(shí)，其散射和輻射模式下的表面電流強(qiáng)度分布，場(chǎng)強(qiáng)分布和方向性圖，并與FEKO中的方向性圖試驗(yàn)結(jié)果進(jìn)行比較，證明了這種方法的正確性。隨后進(jìn)行了單極子天線和組合單極子天線的特性比較，為實(shí)際工程中的仿真計(jì)算提供有效的參考依據(jù)。

［1］HARRINGTON R F.Field Computation by Moment Methods［M］.New York： Macmillan，1968.

［2］闕肖峰，聶在平，宗顯政.復(fù)雜金屬載體上線天線的MoM分析［J］.微波學(xué)報(bào)，2006，22（5）：16－20，38

［3］梁榮江.基于Matlab分析矩量法在線天線中的應(yīng)用［J］.西安郵電學(xué)院學(xué)報(bào)，2009.14（1）：49－52，56

［4］MAKAROV S.MoM antenna simulations with MATLAB：RWG basis functions ［J］.IEEE Antennas and Propagation Magazine，2001，43（5）：100－107.

［5］張?jiān)迫A.用基于rao－wilton－glisson基函數(shù)的矩量法分析線天線 ［J］.系統(tǒng)工程與電子技術(shù)，2005，27 （6）：1105－1108.

［6］董健，柴舜連，毛鈞杰.任意形狀線、面、體組成導(dǎo)體目標(biāo)的電磁建模［J］.電子學(xué)報(bào)，2005，33（9）：1656－1659.

［7］芮翔，徐忠偉，程嗣怡.基于RWG基函數(shù)的平面螺旋天線建模及仿真研究［J］.現(xiàn)代防御技術(shù)，2009，37（4）：109－113.

［8］廖艷蘋，楊莘元，馬惠珠，侯艷麗.有限偶極子天線陣的電磁散射理論研究［J］.兵工學(xué)報(bào)，2007，28（2）：182－186.

［9］RAO S M，WILTON D R，Glisson A W.Electromagnetic Scattering by Surfaces of Arbitrary Shape ［J］.IEEE Trans Antennas and Propagation，1982，30（3）：409－418.

［10］MAKAROV S N.通信天線建模與MATLAB仿真分析.許獻(xiàn)國(guó)，譯［M］.北京：北京郵電大學(xué)出版社，2006.

［11］JOHN L V.Antenna engineering handbook，fourth edition［M］.New York：McGraw－Hill，2007.

［12］BALANIS C A.Antenna theory：analysis and design，third edition［M］.New York：Wiley，2005.

Analysis on the Electromagnetic Scattering and Radiation of Wire Antenna by Moment Method

Gao Bao-chun Peng Ge Hu Yi-feng
China Ship Development and Design Center， Shanghai Division，Shanghai 201108，China

Analytical method was chosen to describe the structure of wire antenna.The wire antenna was modeled as a thin strip， wire－surface junction problem was transformed to surface－surface junction problem.RWG function was chosen to discretize the surface of antenna.Current distribution， radiation field and directivity of different radius dipole antennas were analyzed by using method of moment，and then the results of radiation pattern and scattering pattern were compared.The characteristics of 10m monopole antenna on the finite plate were calculated.The results show that： when the radius and or wavelength are small， the current flows along the antenna axial， the line current model is applicable to the problems of wire antennas， but when they are large， horizontal currents cannot be ignored， the line current model is unsuitable for the analysis of wire antenna problems.

method of moment； RWG basis function； wire antenna； electromagnetic scattering； electromagnetic radiation

U674.7

A

1673－3185（2011）03－49－06

10.3969/j.issn.1673-3185.2011.03.011

2010-09-03

“十一·五”海軍裝備預(yù)先研究項(xiàng)目（40******101）；“十一·五”海軍武器裝備研制項(xiàng)目（09***／***3）

高寶春（1987－），男，碩士研究生。研究方向：艦船電磁兼容，計(jì)算電磁學(xué)。E-mail：gaobaochun＠yahoo.cn