對于仿真軟件相位方向圖計算的討論

2016-12-27 by:CAE仿真在線來源:互聯(lián)網(wǎng)

摘要:GPS 接收天線,反射面饋源等天線的設(shè)計需要對于天線的相位特征加以考慮。仿真軟件在相位方向圖的計算上可能存在一定誤差,繼而影響到對于相位特征有較高要求的天線設(shè)計。本篇文章通過比較對于同一模型在不同仿真仿真軟件下與不同求解設(shè)臵下的計算結(jié)果,說明仿真軟件在天線相位方向圖上的誤差,并主要分析該誤差在CST 微波工作室下產(chǎn)生的原因。

1 引言

由于高精度的全球定位系統(tǒng)需要用戶段的接收天線能夠準(zhǔn)確地測量載波相位,因此指標(biāo)要求天線在近上半空間(俯仰角在±80°范圍內(nèi))的相位誤差小于10°,在L1 頻段(1.575GHz)該指標(biāo)對應(yīng)長度為5mm[1]。在一些文獻(xiàn)中,對于仿真模型的輻射場相位誤差描述則達(dá)到了0.1mm 甚至0.01mm 數(shù)量級。然而仿真軟件的相位方向圖計算能否達(dá)到這一精度則少有文章提及。

本篇文章首先通過二元點(diǎn)源陣與半波陣子天線理論模型,說明計算邊界大小對于相位方向圖的影響,繼而分別在CST 微波工作室,HFSS 與Feko

Suite 電磁仿真軟件中建立半波振子天線模型,并將仿真的相位方向圖結(jié)果與理論計算值進(jìn)行比較。最后則分析CST 微波工作室相位方向圖計算誤差產(chǎn)生的原因。

2 計算邊界對于二元點(diǎn)源陣相位方向圖的影響

IEEE 對于天線相位方向圖與相位中心定義如下[2]:相位方向圖為天線輻射場矢量的相位空間分布,而相位中心則是與天線相關(guān)的一個點(diǎn),在以該點(diǎn)為球心半徑延伸至遠(yuǎn)區(qū)場的球面整體或是球面的一部分上,天線輻射場的相位幾乎是不變的。著作[3]中以如圖1 所示的等幅同相的二元理想各向同性點(diǎn)源陣為例,說明相位方向圖的結(jié)果與計算坐標(biāo)系原點(diǎn)選取的關(guān)系,并指出當(dāng)計算輻射場的坐標(biāo)系原點(diǎn)選取在二元陣的中點(diǎn)時,二元陣的相位方向圖為階躍函數(shù):若陣元間距d 小于二分之波長,那么二元陣的相位方向圖函數(shù)為一常數(shù),。當(dāng)天線陣元間距大于二分之波長時,天線幅度方向圖會有零點(diǎn)產(chǎn)生,在幅度方向圖零點(diǎn)位置相位方向圖會出現(xiàn)180°跳變。

可見在有限大的計算邊界下,二元陣的相位方向圖會因?yàn)橛嬎氵吔绱笮《兓?。?dāng)計算邊界越大時,相位方向圖的波動越小,其函數(shù)表達(dá)式約接近于無限大邊界下的階躍函數(shù)。

3 理論與仿真計算半波陣子E 面相位方向圖比較

繼而分析有限長度線源的相位方向圖。令長度為d 的線源沿z 軸放置,坐標(biāo)原點(diǎn)位于線源幾何中心。那么在球坐標(biāo)系下任一場點(diǎn)(X,Y,Z)距離線源上

任一點(diǎn)(0,0,z)距離r1 可寫作

那么由公式即可繪出半波陣子天線在不同計算半徑r 下的E 面的相位方向圖并與通過CST Microwave Studio, Feko Suite 與HFSS 三種仿真軟件

計算的相位方向圖進(jìn)行比較。仿真軟件中的振子天線工作在L1 頻段(1.575GHz)附近,振子臂半徑為1mm,臂長經(jīng)調(diào)節(jié)為41mm( 0.21入)。比較的結(jié)果如圖4 所示。

可見于第二節(jié)中點(diǎn)源陣相位方向圖隨半徑r 的變化趨勢類似,線源天線理論計算的相位方向圖波動也隨著半徑增大而減小。而由于理想的輻射邊界應(yīng)在形狀上近似球面且在尺寸上盡量大,由圖3 中可見,CST 軟件計算得的相位方向圖與理論計算在函數(shù)曲線上更為相似,這也印證了文章[3]的觀點(diǎn),即CST 的計算邊界更近似于球面。通過CST,Feko與HFSS 軟件計算得的相位方向圖均方根誤差分別為0.92°,0.31°與0.82°,在L1 頻段分別對應(yīng)長

度為0.48mm,0.16mm 與0.43mm。

由于在理想的無窮遠(yuǎn)邊界下,半波陣子天線的E 面相位誤差應(yīng)趨于0。因此對于仿真軟件計算得的E 面均方根誤差可以認(rèn)為是仿真軟件造成的計算方向圖的誤差。因此將在下節(jié)中討論如何減小這種仿真軟件在相位方向圖計算上的誤差。

4 仿真軟件相位方向圖計算誤差產(chǎn)生分析

在HFSS 軟件中的輻射邊界是用來模擬電磁波朝無窮遠(yuǎn)輻射的開放問題,而這類問題在天線設(shè)計中十分常見。一般使用的輻射邊界采用第一階吸收邊界條件(ABC)[4]。文獻(xiàn)[5]中的相應(yīng)章節(jié)有對于該類邊界條件的討論,并給出該類吸收邊界對于不同極化方式,不同入射角的電磁波的反射系數(shù)。

而在CST 中計算輻射場則是使用的開放(openadd space)邊界條件,[6]中說明該邊界條件可以視作PML 邊界加上了一些額外的空間用來計算遠(yuǎn)區(qū)場,

并提到PML 層數(shù)為為取默認(rèn)的四層就已足夠。在[7]中則說明增加PML 層數(shù)可以使得局部和總體誤差都單調(diào)地減小,但是PML 層數(shù)過多又會導(dǎo)致計算

量的劇增。因此折衷考慮吸收邊界效果與計算量,可以選取層數(shù)4 至8 層。

因此通過調(diào)節(jié)CST 微波工作室仿真環(huán)境下半波陣子天線包括求解頻率,計算邊界大小(默認(rèn)為計算頻段中頻波長的八分之一),PML 邊界層數(shù)(默認(rèn)為4 層)的計算設(shè)置,來比較天線E 面相位方向圖在θ為10°~170°范圍內(nèi)的均方根誤差,通過均方根誤差值來比較不同設(shè)置下相位方向圖的計算結(jié)果,并將均方根誤差值對應(yīng)到L1 頻段波長。比較結(jié)果在表1 中給出。仿真采用時域求解器,在不同頻率下計算的網(wǎng)格數(shù)單元數(shù)均為56400 個。

可見,在CST 微波工作室仿真環(huán)境下,計算邊界大小對于相位方向圖計算結(jié)果并沒有直接影響。增加PML 層數(shù)可以在一定范圍內(nèi)減小相位方向圖的計算誤差。同時相位方向圖的計算結(jié)果與求解頻率關(guān)系密切,擴(kuò)展求解頻率范圍可以減小相位方向圖的均方根誤差,即令計算邊界更接近于理想的輻射邊界。同一天線模型在不同的求解頻率下的均方根相差對應(yīng)到距離可以達(dá)到0.75mm(編號1 與編號7 的數(shù)據(jù))。經(jīng)過仿真還發(fā)現(xiàn)增加網(wǎng)格單元數(shù),對于相位方向圖結(jié)果影響不大。

5 結(jié)論

不同仿真軟件對于工作在L1 頻段的半波振子天線相位方向圖均方根誤差的結(jié)果差值可以達(dá)到0.3mm。而在CST 微波工作室下的同一半波陣子天線模型,在不同的求解設(shè)置下對于相位方向圖的均方根誤差計算結(jié)果,相差可以達(dá)到0.75mm。可見不同仿真軟件與同一軟件的不同設(shè)置對于相位方向圖的計算存在一定影響。在CST 微波工作室下擴(kuò)展求解頻率可以在一定程度上減小計算相位方向圖的誤差。而對于工作在L 波段的衛(wèi)星接收天線而言,相位方向圖均方根誤差0.01mm 的精度可能已經(jīng)超過了仿真軟件對于相位方向圖計算的誤差限。

參考文獻(xiàn)

[1] Tranquilla J M, Colpitts B G. ―GPS antenna design characteristics for high-precision applications‖[J]. Journal of Surveying Engineering, 1989, 115(1): 2-14.

[2] IEEE Standard Definitions of Terms for Antennas Sponsored by Antenna Standards Committee of the IEEE Antennas and Propagation Society

[3] John D. Kraus, ―Antennas‖ New York : McGraw-Hill, c1988. 2nd ed..pp 137-140

[4] Investigation of antenna phase pattern by the electromagnetic simulation software, Ze Shen, Zhenxin Cao,2014 ICCP&HSIC

[5] HFSS Help-Assigning boundaries-Assigning Radiation Boundaries;

[6] 金建銘(著) 王建國(譯), 電磁場有限元方法, 西安電子科技大學(xué)出版社,1998,P308~P317。

[7] CST Help Content-Boundary Conditions

[8] 王秉中, 計算電磁學(xué). 科學(xué)出版社, 2002,第四章

開放分享：優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才