嵌地線天線的輻射特性分析

廖飛龍 張東升 范啟蒙 王劉軍

(1.陸軍軍官學院,合肥 230031;2.電子工程學院,合肥 230031)

?

嵌地線天線的輻射特性分析

廖飛龍1張東升1范啟蒙2王劉軍1

(1.陸軍軍官學院,合肥 230031;2.電子工程學院,合肥 230031)

采用混合位積分方程結合矩量法對嵌地線天線進行了分析,討論了Sommerfeld積分、阻抗矩陣元素及輻射場的計算等關鍵技術,仿真計算了嵌地線天線的輻射特性,并和FEKO軟件的計算結果進行比對,驗證了分析方法的正確性.與FEKO等電磁軟件相比,文中方法在分析跨嵌電磁目標時,其計算效率具有明顯的優(yōu)勢.

平面分層媒質;混合位積分方程;嵌地線天線;Sommerfeld積分;離散復鏡像法

DOI 10.13443/j.cjors.2016061501

引 言

平面分層介質結構在單片微波與毫米波集成電路、微帶天線、埋地及嵌地目標的電磁特性分析等領域中,具有廣泛的應用．采用混合位積分方程[1]結合矩量法分析這類問題具有高效、準確等優(yōu)點,因而一直為國內外學者所青睞．當前的研究中,對于電磁目標處于同一層媒質的情況討論較多[2-6],而對于跨界面目標的電磁特性分析卻較為少見[7-9],最主要的困難在于電磁目標跨界面時空域格林函數中所包含Sommerfeld積分的計算．平面分層媒質中Sommerfeld積分的計算,最主流的方法之一是離散復鏡像法(Discrete Complex Image Method,DCIM)[10-11],它通常只對電磁目標處于同一層媒質時有效,當電磁目標處于不同層時,即跨層情況,Sommerfeld積分的計算往往需要進行繁瑣的插值運算[9]．目前相當多的商用電磁軟件僅有2.5維計算能力[12-14],只有IE3D、FEKO等少數軟件能夠進行全三維仿真,而IE3D、FEKO在計算空域格林函數時都采用數值積分和查表內插的方法,計算比較耗時．

本文首先根據混合位積分方程(Mixed-Potential Integral Equation,MPIE),結合矩量法,得到阻抗矩陣元素及其包含的譜域格林函數的表達式．然后分析阻抗矩陣元素計算的方法,并利用Formulation C[1]中標量位格林函數在分界面上對于場點縱坐標連續(xù)的性質,通過算子運算將對標量位格林函數的微分運算轉移至表達式較為簡單的電流基函數中,降低了矩陣元素計算的復雜性和計算量．在阻抗元素的計算中,對不同場源點位置情況下格林函數的計算采用了不同的方法,較好地解決了跨界面電磁目標分析中Sommerfeld積分計算困難的問題,然后推導了輻射場的計算,最后利用Matlab編程仿真了嵌地天線的相關電磁特性．

1 嵌地線天線的矩量法分析

1.1 阻抗矩陣的計算

平面分層媒質中(如圖1所示),根據MPIE[1],第n層中任意點r′處電流源J(r′)在第m層中任意點r處產生的電場可以表示為

圖1 平面分層媒質結構示意圖

(1)

(2)

在矩量法中,阻抗矩陣元素的計算可以歸結為源電流元產生的場和檢驗電流元的卷積運算,電流元Jn(r′)與Jm(r)之間的互阻抗元素按照MPIE可以表示為

Zmn=〈En(r′,r),Jm(r)〉

Jn(r′)·Jm(r)dl′dl+

Jn(r′)·Jm(r)dl′dl.

(3)

′·Jn(r′)dl′.

(4)

討論式(4)中∫l[·Jm]dl積分部分,對于線天線,電流沿導線軸線方向一維分布,故梯度與散度的運算實際上就是一維求導運算,有

∫l.

(5)

對于第m模,如果上下半模處于同一層媒質,則

∫

(6)

如果第m模上下半模處于分界面的兩側,則

∫

(7)

這里利用了Formulation C中標量位格林函數在分界面處關于場點縱坐標連續(xù)的性質,即

(8)

從式(6)和式(7)可以看出,無論第m模的上下半模處于同一層媒質,還是分別處于平面媒質分界面的兩側,都有

∫

(9)

所以式(3)的第二項積分可以寫為

Jn(r′)dl·dl′.

(10)

式(10)無需對標量位格林函數進行微分運算,而將微分算子轉移至表達式相對簡單的電流基函數上,這將使計算量大大減小,并且此式對于電流基函數和檢驗函數具有較好的互易性,便于阻抗矩陣的優(yōu)化計算．

綜合前面的分析,阻抗元素的表達式可寫為

Jm(r)′·Jn(r′)dl′dl.

(11)

1.2 Sommerfeld積分的計算

(12)

求出積分的近似解析表達式．

對于嵌地線天線,源電流和檢驗電流所處的相對位置不同,其格林函數具有不同形式的表達式．在求取阻抗矩陣時,將Sommerfeld積分的計算分下面四種情況分別采用不同的方法．

1) 源電流和檢驗電流都在上半空間．阻抗元素計算過程中,場、源點均在上半空間,采用傳統(tǒng)的二級DCIM[11]計算格林函數中所包含的Sommerfeld積分．

2) 源電流在下半空間,檢驗電流在上半空間．阻抗元素計算過程中,源點在下半空間,場點在上半空間,可采用二維離散復鏡像法(2-Dimensional Discrete Complex Image Method,2D-DCIM)[15]計算格林函數中所包含的Sommerfeld積分,以避免繁瑣的插值運算．

3) 源電流在上半空間,檢驗電流在下半空間．阻抗元素計算過程中,源點在上半空間,場點在下半空間,需要采用文獻[16]提出的改進的二維DCIM計算格林函數中所包含的Sommerfeld積分,以避免介質的有耗特性給計算帶來的錯誤．

4) 源電流和檢驗電流都在下半空間．場源點均在下半空間,如果下層媒質是無耗的,可采用傳統(tǒng)的二級DCIM方法計算格林函數中所包含的Sommerfeld積分;如果下層媒質是有耗的,則需要采用文獻[16]提出的改進的2D-DCIM方法計算位函數中所包含的Sommerfeld積分,以保證計算的準確性．

計算過程中,無論是采用傳統(tǒng)的二級DCIM、2D-DCIM還是改進的2D-DCIM,如果電磁目標的電尺寸較大(lg(k0ρ)>1),都需要采用增強DCIM[17]的采樣路徑,以保證計算精度．

2 輻射場的計算

對于嵌地線天線的輻射場,根據實際工程情況,主要討論天線在上半空間產生的輻射場,包括地上部分在上半空間產生的場和埋地部分在上半空間產生的場的總和．關于地上部分在上半空間產生的場,可通過文獻[18]的方法直接得到,這里主要討論埋地部分在上半空間產生的輻射場．

先計算遠區(qū)磁矢位表達式．假設任意電偶極子的方向為(cos α,cos β,cos γ),其產生的磁矢位為

(13)

這里源點在下層,場點在上層,即m=1,n=2．下面計算矢量位各分量的近似表達式,由于z′?z,可將關于z′的指數部分歸為函數慢變部分,將其對積分的貢獻歸為駐相點的值從積分中提出[18],剩余部分用Sommerfeld恒等式得到即可．于是

(14)

分別計算A12(r)的各個分量,根據式(13)有

(15)

(16)

將地上部分和埋地部分產生的場相加可得總的輻射場．

3 數值計算

編程分析了如圖2所示的λ天線,此天線架設在地面上空,兩邊用地釘固定,左端接地處用電壓源饋電．在目前的實際工程中,都是將大地當成無限大理想金屬導體板處理,也沒有考慮地釘的作用,這樣處理跟實際問題并不相符．本文將地面作為平面分層媒質,對地釘和天線一起進行剖分建模,展開矩量法分析計算．天線尺寸如圖3所示,地釘長度為0.5 m,半徑為0.005 m,天線的半徑為0.015 m,地質相對介常數為εr2=4,地質電導率為σ2=0.001 S/m．

圖2 λ天線架設示意圖

(a) E面方向圖

(b) H面方向圖圖3 天線的方向圖

首先仿真了f=10 MHz時天線的方向圖,并與FEKO軟件的仿真結果進行了比對．從圖3可以看出,本文方法計算結果與軟件的結果吻合良好,證實了本文方法的正確性．

為了解地質參數對天線性能的影響,計算了天線在不同地質參數的方向圖,如圖4所示．天線工作頻率f=10 MHz,地質相對介電常數為εr2=4,圖4(a)表示φ=0°時,垂直平面內的方向圖,圖4(b)表示θ=90°時,水平面內的方向圖．從圖4可以看出,在保持其他參數不變時,地質參數的變化對天線方向圖具有較大影響,并且隨著地質電導率的增大,地質金屬性增強,向上半空間反射的電場增大,從而使上半空間總的輻射場變大．

比較FEKO及本文方法的計算效率．天線工作頻率為f=20 MHz,在Inter I5處理器、2G內存的計算機上,分別采用兩種方法計算不同地質參數情況下天線的方向圖,所用時間如表1所示．可以看出,隨著地質電導率及相對介電常數的增加,FEKO計算時間急劇增加;而本文方法的計算時間基本穩(wěn)定,并且天線電尺寸越大,結構越復雜,本文方法優(yōu)勢將越明顯．這是因為FEKO在計算Sommerfeld積分時采用了數值積分的方法,每改變一次場源點位置都需要重新計算Sommerfeld積分的值,并且在地質電導率及相對介電常數變大的情況下,Sommerfeld積分的被積函數劇烈振蕩,數值方法對其進行計算的時間就會急劇增加．而本文方法采用解析方法計算Sommerfeld積分,在一定地質條件下獲得格林函數解析表達式以后,只要代入不同位置點的坐標即可得到不同位置坐標的格林函數,故計算穩(wěn)定．

(a) E面方向圖 (b) H面方向圖圖4 改變地質電導率時近地斜架λ天線的方向圖

所用方法地質參數σ2=0.1S/mεr2=4σ2=0.01S/mεr2=4σ2=0.001S/mεr2=4σ2=0.001S/mεr2=20σ2=0.001S/mεr2=80FEKO339s61s32s496s1551s本文計算方法79s72s66s82s93s

4 結 論

跨嵌導體目標的電磁特性分析由來是分層媒質中電磁問題的難點,本文通過混合位積分方程結合矩量法分析嵌地線天線,仿真計算了天線的電磁輻射特性．從計算結果可以看出:地質參數的改變對于天線的電磁輻射特性具有較大影響,天線在架設使用時應充分考慮不同地面的影響;由于FEKO等電磁軟件在處理跨層電磁目標問題時,對于Sommerfeld積分的計算多采用數值積分法,而本文采用解析的方法,在計算效率上具有明顯的優(yōu)勢．為進一步提高嵌地線天線分析的效率,阻抗矩陣的優(yōu)化計算將是后續(xù)研究的方向．

[1]MICHASKI K A, ZHENG D L.Electromagnetic scattering and radiation by surfaces of arbitrary shape in layered media, part I and part II [J].IEEE transactions on antennas and propagation, 1990, 38(3)：335-352.

[2]LEE J, ZOU J, J WAN, et al.Fast evaluating the electromagnetic fields generated by a vertical electric dipole over the lossy ground using Sommerfeld integral without truncation[J].IEEE transactions on electromagnetic compatibility, 2011, 53(4)：977-986.

[3]WU B Y, SHENG X Q.A complex image reduction technique using genetic algorithm for the MoM solution of half-space MPIE [J].IEEE transactions on antennas and propagation, 2015, 63(8)：3727-3731.

[4]SHOORY A, MOINI R, SADEGHI S H H.Direct use of discrete complex image method for evaluating electric field expressions in a lossy half space[J].IET microwave antennas propagation, 2010, 4(2)：258-268.

[5]于繼軍, 盛新慶.地下三維電磁散射的矩量法計算[J].電子與信息學報, 2006, 28(5)：950-954.

YU J J, SHENG X Q.Scattering from 3-D targets in the subsurface using MOM[J].Journal of electronics &information technology, 2006, 28(5)：950-954.(in Chinese)

[6]張云華, 肖伯勛, 朱國強.一種快速計算分層媒質中三維目標電磁散射的方法[J].武漢大學學報, 2007, 53(1)：119-122.

ZHANG Y H, XIAO B X, ZHU G Q.Fast evaluation of electromagnetic scattering from three dimensional objects buried in layered media[J].Journal of Wuhan University(natural science edition), 2007, 53(1)：119-122.(in Chinese)

[7]YE H X, JIN Y Q.Dual GPOF/DCIM for fast computation of Sommerfeld integrals and EM scattering from an object partially embedded in dielectric half-space[J].IEEE transactions on antennas and propagation, 2010, 58(5)：1801-1807.

[8]ARNAUTOVSKI-TOSEVA V, GRCEV L.Image and exact models of a vertical wire penetrating a two-layered earth[J].IEEE transactions on electromagnetic compatibility, 2011, 53(4)：968-976.

[9]LIU Y, LI L W, LEONG T S.Appliction of DCIM to MPIE-MOM analysis of 3-D PEC objects in multilayered media [J].IEEE transactions on antennas and propagation, 2002, 50(11)：157-162.

[10]FANG D G, YANG J J, DELISLE G Y.Discrete complex images for horizontal electric dipoles in a multilayered medium[J].IEE proceedings H-microwaves, antennas and propagation,1988, 135(5)：297-303.

[11]AKSUN M I.A robust approach for the derivation of closed-form Green’s function [J].IEEE transactions on microwave theory and techniques, 1996, 44(5)：651-658.

[12]TANG W H, GEDNEY S D..An efficient application of the discrete complex image method for quasi-3-D microwave circuits in layered media[J].IEEE transactions on microwave theory and techniques, 2007, 50(8)：1723-1729.

[13]尹成友.平面分層媒質空域Green函數快速計算(I) [J].微波學報, 2012, 28(2)：1-6.

YIN C Y.Fast computation of spatial Green’s functions for planar-layered media(I)[J].Journal of microwaves, 2012, 28(2)：1-6.(in Chinese)

[14]尹成友.平面分層媒質空域Green函數快速計算(Ⅱ)[J].微波學報, 2012, 28(3):1-8.

YIN C Y.Fast computation of spatial Green’s functions for planar-layered media(Ⅱ)[J].Journal of microwaves, 2012, 28(3)：1-8.(in Chinese)

[15]YUAN M T, ZHANG Y, DE A, et al.Two-dimensional discrete complex image method(DCIM) for closed-form Green’s function of arbitrary 3D structures in general multilayered media[J].IEEE transactions on antennas and propagation, 2008, 56(5)：1350-1356.

[16]廖飛龍, 尹成友, 杜紅兵, 等.一種有效計算有耗分層媒質中的格林函數的方法[J].微波學報, 2010, 26(1)：27-32.

LIAO F L, YIN C Y, DU H B, et al.A robust approach for derivation of Green’s function for lossy layered media[J].Journal of microwaves, 2010, 26(1)：27-32.(in Chinese)

[17]YUAN M T, SARKAR T K, SALAZAR-PALAM M.A direct discrete complex image method from the closed-form Green’s functions in multilayered media [J].IEEE transactions on microwave theory and techniques, 2006, 54(3)：1025-1032.

[18]CHEW A W.A quick way to approximate a Sommerfeld-Weyl-type integral [J].IEEE transactions on antennas and propagation, 1988, 36(11)：1654-1657.

Radiation characteristics analysis of wire antenna partially embedded in the ground

LIAO Feilong1ZHANG Dongsheng1FAN Qimeng2WANG Liujun1

(1.ArmyOfficerAcademyofPLA,Hefei230031,China;2.ElectronicEngineeringInstituteofPLA,Hefei230031,China)

Mixed-potential integral equation (MPIE) combined with method of moment is used to analyze the wire antenna partially embedded in the ground.Some key techniques such as the calculation of Sommerfeld integrals, impedance matrix element and radiation field are discussed in this paper.Radiation characteristics of wire antennas partially embedded in the ground are calculated, and the results obtained are compared with that of the FEKO software, which confirms the correctness of this method.Compared with FEKO and other electromagnetic software, the approach presented in this paper has obvious advantages in computational efficiency for electromagnetic characteristics analysis of the object penetrating planar layered media.

plannar layered media;mixed-potential electric field integral equation;wire antenna partially embedded in the ground;Sommerfeld Integral;discrete complex image method

廖飛龍, 張東升, 范啟蒙, 等.嵌地線天線的輻射特性分析[J].電波科學學報,2016,31(5):858-863.

10.13443/j.cjors.2016061501

LIAO F L, ZHANG D S, FAN Q M, et al.Radiation characteristics analysis of wire antenna partially embedded in the ground [J].China journal of radio science,2016,31(5):858-863.DOI：10.13443/j.cjors.2016061501

2016-06-15

安徽省自然科學基金項目(1408085QF121)

TN92

A

1005-0388(2016)05-0858-06

廖飛龍 (1981-),男,四川人,博士,講師,主要研究方向為電磁數值計算、天線分析設計等．

張東升 (1968-),男,安徽人,副教授,碩士研究生導師,主要研究方向為無人機飛控技術等．

范啟蒙 (1993-),男,河南人,碩士研究生,主要研究方向為電磁數值計算．

聯(lián)系人：廖飛龍 E-mail：liaofeilong0514@163.com