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

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

摘要:GPS 接收天線,反射面饋源等天線的設(shè)計(jì)需要對于天線的相位特征加以考慮。仿真軟件在相位方向圖的計(jì)算上可能存在一定誤差,繼而影響到對于相位特征有較高要求的天線設(shè)計(jì)。本篇文章通過比較對于同一模型在不同仿真仿真軟件下與不同求解設(shè)臵下的計(jì)算結(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ù)量級。然而仿真軟件的相位方向圖計(jì)算能否達(dá)到這一精度則少有文章提及。

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

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

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

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

可見在有限大的計(jì)算邊界下,二元陣的相位方向圖會因?yàn)橛?jì)算邊界大小而變化。當(dāng)計(jì)算邊界越大時(shí),相位方向圖的波動(dòng)越小,其函數(shù)表達(dá)式約接近于無限大邊界下的階躍函數(shù)。

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

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

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

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

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

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

度為0.48mm,0.16mm 與0.43mm。

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

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

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

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

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

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

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

可見,在CST 微波工作室仿真環(huán)境下,計(jì)算邊界大小對于相位方向圖計(jì)算結(jié)果并沒有直接影響。增加PML 層數(shù)可以在一定范圍內(nèi)減小相位方向圖的計(jì)算誤差。同時(shí)相位方向圖的計(jì)算結(jié)果與求解頻率關(guān)系密切,擴(kuò)展求解頻率范圍可以減小相位方向圖的均方根誤差,即令計(jì)算邊界更接近于理想的輻射邊界。同一天線模型在不同的求解頻率下的均方根相差對應(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è)置下對于相位方向圖的均方根誤差計(jì)算結(jié)果,相差可以達(dá)到0.75mm?？梢姴煌抡孳浖c同一軟件的不同設(shè)置對于相位方向圖的計(jì)算存在一定影響。在CST 微波工作室下擴(kuò)展求解頻率可以在一定程度上減小計(jì)算相位方向圖的誤差。而對于工作在L 波段的衛(wèi)星接收天線而言,相位方向圖均方根誤差0.01mm 的精度可能已經(jīng)超過了仿真軟件對于相位方向圖計(jì)算的誤差限。

參考文獻(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] 王秉中, 計(jì)算電磁學(xué). 科學(xué)出版社, 2002,第四章

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