金濤斌，劉正文，羅 雨，羅毅彪

(1.呼和浩特市城市軌道交通建設(shè)管理有限責(zé)任公司，呼和浩特 010000；2.北京石油化工學(xué)院 機(jī)械工程學(xué)院，北京 102617)

0 引言

20世紀(jì)50年代以來，頻率選擇表面(Frequency Selective Surfaces，F(xiàn)SS)的工作特性和工程應(yīng)用得到了電磁領(lǐng)域的廣泛重視和深入研究。通常是將單層或多層的二維金屬單元結(jié)構(gòu)進(jìn)行周期性排列以構(gòu)成FSS，而生成金屬單元結(jié)構(gòu)的方式通常有兩類，其中一類是將金屬貼片單元蝕刻在介質(zhì)底板上，以構(gòu)成金屬單元；另一類是直接在金屬板上開刻缺陷單元，以構(gòu)成金屬單元。FSS實(shí)質(zhì)上就是一種空間濾波器，廣泛應(yīng)用于微波工程領(lǐng)域。周期性矩形缺陷接地結(jié)構(gòu)(Periodic Rectangle Defected Ground Structure，PRDGS)作為一類FSS，它是在微波電路的金屬板上開刻周期性的矩形縫隙單元，以抑制某些頻段電磁波的傳播，從而改變頻率特性，在頻譜上形成十分顯著的阻帶[1-4]。PRDGS能被用來實(shí)現(xiàn)微波濾波器，其阻帶特性可使得天線在自身性能、抑制諧波、輸出功率及功放效率等方面獲得顯著的提升，因此在天線及微波電路領(lǐng)域具有非常廣大的應(yīng)用和市場(chǎng)前景。

對(duì)PRDGS傳輸特性進(jìn)行的計(jì)算分析，一般有三類方法。第一種為頻域方法，如頻域有限元方法(Frequency Domain Finite Element Method，F(xiàn)DFEM)。這類方法的發(fā)展比較成熟，通常我們?cè)趩我活l率條件下應(yīng)用此種頻域方法來分析計(jì)算PRDGS是完全可行且簡(jiǎn)便的，但以該方法分析處理寬頻帶的電磁信號(hào)時(shí)，需要在很多不同頻率點(diǎn)上耗費(fèi)過多運(yùn)算量，這將直接導(dǎo)致整體計(jì)算效率低下，從而限制該訪求的應(yīng)用范圍[5-6]。由于PRDGS在寬頻帶范圍具有非常顯著的阻帶，如果采用頻域方法去計(jì)算PRDGS的傳輸特性，需要計(jì)算很多的頻率點(diǎn)，從而大大降低計(jì)算效率。因此，計(jì)算PRDGS的寬頻帶特性，一般不宜采用頻域方法。第二種為等效電路模型法。這種方法的優(yōu)點(diǎn)在于電路模型直觀、物理概念清晰，而這種方法的不足在于提取等效電路參數(shù)比較復(fù)雜。對(duì)于具有簡(jiǎn)單結(jié)構(gòu)的PRDGS，采用等效電路模型法進(jìn)行計(jì)算，是直觀可行的。但是對(duì)于結(jié)構(gòu)比較復(fù)雜的PRDGS，采用等效電路模型法進(jìn)行計(jì)算，會(huì)導(dǎo)致電路參數(shù)的提取較為困難，且計(jì)算精度不高[7-8]。第三種為時(shí)域方法，如時(shí)域有限差分法(Finite Difference Time Domain Method，F(xiàn)DTD)或無條件穩(wěn)定的時(shí)域有限差分(Crank-Nicolson FDTDM，CN-FDTD)。這種方法是在獲得足夠的時(shí)域信息后，再通過時(shí)頻變換的簡(jiǎn)單處理，即可獲得寬頻帶范圍內(nèi)的可用頻域信息，因而應(yīng)用該方法的優(yōu)點(diǎn)是時(shí)域下對(duì)具有寬頻帶特性的瞬變電磁場(chǎng)可直接分析處理[9-12]。

有限元方法(Finite Element Method，F(xiàn)EM)雖然也可以完成時(shí)域下的計(jì)算，但其在每一個(gè)時(shí)間步都需要反復(fù)求解大型方程組，以致計(jì)算的時(shí)間復(fù)雜度和空間復(fù)雜度都很高[10]。因CN-FDTD法特別適合模擬具有寬頻帶、非線性特征的精細(xì)結(jié)構(gòu)，而PRDGS本身結(jié)構(gòu)較為簡(jiǎn)單，又兼具寬頻帶特性，故本文應(yīng)用CN-FDTDM方法直接處理PRDGS的寬頻帶特性。

1 等效電路模型法和CN-FDTD法

1.1 等效電路模型法

根據(jù)微波電路理論，可以得出缺陷接地結(jié)構(gòu)(Defected Ground Structure，DGS)單元的等效電路模型。文獻(xiàn)[11]給出了通常DGS單元的等效電路模型，具有該模型結(jié)構(gòu)的微波電路在特定頻率上具有衰減極點(diǎn)和截止頻率，可以近似等效為并聯(lián)LC電路，如圖1所示。

圖1中的LC等效電路具有電路參數(shù)ZLC，XLC，其求取方法如下式，其中，DGS單元的電抗ZLC可表示為：

(1)

(2)

(3)

由于DGS單元在特定頻率上具有衰減極點(diǎn)、截止頻率，所以可等效為1階的Butterworth低通濾波電路，如圖2所示。

依據(jù)微波電路理論，可得圖1中的等效電容表達(dá)式和等效電感表達(dá)式為：

(4)

(5)

通過以上分析，可以得到DGS單元的等效電路模型和相關(guān)電路參數(shù)，再通過計(jì)算可得相關(guān)微波電路的傳輸特性。

1.2 CN-FDTD法

CN-FDTD法采用Yee離散格式，如圖3所示，其中電場(chǎng)和磁場(chǎng)在空間中相互錯(cuò)置，即每個(gè)電場(chǎng)分量均環(huán)繞四個(gè)磁場(chǎng)分量，同理每個(gè)磁場(chǎng)分量環(huán)繞四個(gè)電場(chǎng)分量，同時(shí)電場(chǎng)與磁場(chǎng)間在時(shí)間上錯(cuò)開半個(gè)步長(zhǎng)，即時(shí)間上相差半個(gè)步長(zhǎng)。因而，當(dāng)待解電磁場(chǎng)的初始值與邊界條件存在時(shí)，可由CN-FDTD法單步遞推來求解以后各時(shí)刻的空間內(nèi)的電磁場(chǎng)。

通過分析可知，由于CN-FDTD法是一種無條件穩(wěn)定的FDTD方法，因而對(duì)于在CN-FDTD法中的時(shí)間步長(zhǎng)的增長(zhǎng)因子，始終可保證||≤1成立，由離散的Maxwell方程可推得一組離散方程組，將聯(lián)立離散后的方程組通過整理，可得關(guān)于E分量的線性方程組，此處直接給出線性方程組如式(6)所示，其它類似的線性方程組可同理得到，具體論證詳見相關(guān)文獻(xiàn)[9]，然后，求解該線性方程組可解得E分量，最后再解得H分量。在CN-FDTD法的計(jì)算過程中可引入高斯脈沖源，對(duì)CN-FDTD過程中的稀疏矩陣方程組可以采用共軛梯度法進(jìn)行處理。

(6)

2 數(shù)值結(jié)果及分析

2.1 PRDGS結(jié)構(gòu)

PRDGS結(jié)構(gòu)示意圖和相關(guān)結(jié)構(gòu)參數(shù)，如圖4和圖5所示；完全匹配層(Perfectly Matched Layer，PML)以及吸收邊界條件的設(shè)置如圖6所示。

2.2 PRDGS的傳輸特性

采用等效電路模型法對(duì)圖示PRDGS的傳輸特性進(jìn)行計(jì)算，可得其阻帶中心頻率f0：

(7)

式中，c為光速，εr為介質(zhì)的相對(duì)介電常數(shù)，d為周期單元的間距。

采用CN-FDTD法對(duì)圖示PRDGS的傳輸特性進(jìn)行計(jì)算。相關(guān)結(jié)構(gòu)參數(shù)和計(jì)算參數(shù)的取值為l1=120 mm，l2=30 mm，w=3 mm，d=20 mm，a=7 mm，b=7 mm，x=0 mm、y=0 mm，z=0 mm，介質(zhì)板厚度設(shè)為1 mm，周期單元設(shè)為5，離散網(wǎng)格為4×60×480，介質(zhì)的相對(duì)介電常數(shù)r分別取2.65，2.90和3.30。PML吸收層數(shù)設(shè)為15層，傳統(tǒng)FDTD法時(shí)間步長(zhǎng)設(shè)為Courant-Friedrich-Levy(CFL)時(shí)間步長(zhǎng)，即tFDTD=0.42 ps，推進(jìn)步數(shù)為18 000 步，物理時(shí)間為7 560 ps。而CN-FDTD法的時(shí)間步長(zhǎng)分別設(shè)為CFL時(shí)間步長(zhǎng)的2倍、6倍、10倍、14倍、18倍、22倍，即2tFDTD=0.84 ps，6tFDTD=2.52 ps，10tFDTD=4.2 ps，14tFDTD=5.88 ps，18tFDTD=7.56 ps，22tFDTD=9.24 ps，其它如物理時(shí)間等設(shè)置不變。

經(jīng)過CN-FDTD法計(jì)算，在矩形單元尺寸一定時(shí)，介電常數(shù)取值不同情況下，PRDGS的傳輸特性，如圖7所示。圖中橫軸為頻率，縱軸為幅頻特性，由圖可知，隨著介質(zhì)的介電常數(shù)不斷增大，阻帶中心頻率整體向左偏移，而阻帶的帶寬和深度基本不變。將等效電路模型法和CN-FDTD法所計(jì)算出的阻帶中心頻率f0進(jìn)行比較，如表1所示。由表可知，兩種方法所計(jì)算出的阻帶中心頻率一致性很好，平均相對(duì)誤差不超過1.72 %。其中方法一為等效電路模型法，方法二為CN-FDTD法，平均相對(duì)誤差以CN-FDTD計(jì)算結(jié)果為標(biāo)準(zhǔn)。

進(jìn)一步，在介電常數(shù)一定時(shí)，采用CN-FDTD對(duì)不同矩形單元尺寸情況下的PRDGS傳輸特性進(jìn)行計(jì)算，計(jì)算結(jié)果分別如圖8和圖9所示。此時(shí)介質(zhì)的相對(duì)介電常數(shù)r=2.65，其它結(jié)構(gòu)參數(shù)和計(jì)算參數(shù)取值不變。

表1 兩種方法所計(jì)算出的阻帶中心頻率

圖8給出了PRDGS在保持b=7 mm不變，a分別等于6 mm，7 mm，8 mm時(shí)，阻帶特性的變化。圖9給出了PRDGS在保持a=7 mm不變，b分別為6 mm，7 mm，8 mm時(shí)，阻帶特性的變化。由圖可知，PRDGS具有良好的阻帶特性，其阻帶中心頻率在4.6 GHz左右，符合等效電路模型法的計(jì)算結(jié)果，且阻帶寬度和阻帶深度會(huì)隨著矩形單元尺寸的增大而增大。

2.3 CN-FDTD時(shí)間步長(zhǎng)與計(jì)算效率、計(jì)算精度

為衡量CN-FDTD法相較于傳統(tǒng)FDTD法在計(jì)算效率和計(jì)算精度上的優(yōu)勢(shì)，引入時(shí)間節(jié)省率rTS和平均相對(duì)誤差eMR兩參量，同時(shí)將CN-FDTD法的時(shí)間步長(zhǎng)與CFL時(shí)間步長(zhǎng)的比值作為變參量，用nCFLN表示。時(shí)間節(jié)省率rTS，為將傳統(tǒng)FDTD法計(jì)算所需時(shí)間與采用CN-FDTD法計(jì)算所需時(shí)間做差，再除以傳統(tǒng)FDTD法計(jì)算所需時(shí)間后所得結(jié)果的百分比。平均相對(duì)誤差eMR，為將CN-FDTD法計(jì)算結(jié)果與傳統(tǒng)FDTD法計(jì)算結(jié)果做差，再除以傳統(tǒng)FDTD法計(jì)算結(jié)果，將所得結(jié)果的絕對(duì)值取百分比，然后再取其平均值。

由上可做出采用CN-FDTD法后衡量計(jì)算效率的時(shí)間節(jié)省率rTS曲線圖，如圖10所示，其橫軸為時(shí)間步長(zhǎng)倍率nCFLN。當(dāng)nCFLN取值為2，6，10，14，18，22，即表示CN-FDTD法計(jì)算的時(shí)間步長(zhǎng)取為CFL時(shí)間步長(zhǎng)的2倍、6倍、10倍、14倍、18倍、22倍時(shí)，縱軸時(shí)間節(jié)省率rTS分別對(duì)應(yīng)9.1%，37.2%，62.8%，77.2%，82.8%，86.5%。觀察圖中rTS曲線可知，隨著nCFLN倍數(shù)增大，衡量計(jì)算效率的時(shí)間節(jié)省率rTS參量不斷提高，即采用CN-FDTD法時(shí)所設(shè)時(shí)間步長(zhǎng)越長(zhǎng)總體效率越高。

同理可做出采用CN-FDTD法后衡量計(jì)算精度的平均相對(duì)誤差eMR曲線圖，如圖11所示，其橫軸亦為nCFLN。當(dāng)nCFLN時(shí)間步長(zhǎng)倍率逐漸增大時(shí)，平均相對(duì)誤差eMR越亦逐漸增大時(shí)，即計(jì)算精度越低。綜上，在采用CN-FDTD法進(jìn)行數(shù)據(jù)處理時(shí)，需根據(jù)實(shí)際要求平衡計(jì)算效率與計(jì)算精度兩者的權(quán)重，一般情況下，應(yīng)根據(jù)實(shí)際需要選擇能夠同時(shí)兼顧計(jì)算效率和計(jì)算精度的時(shí)間步長(zhǎng)，而當(dāng)需要較高精度時(shí)，則需要在允許的范圍內(nèi)適當(dāng)降低計(jì)算效率，反之亦然。

3 結(jié)論

采用等效電路模型法和CN-FDTD法，研究了不同介電常數(shù)對(duì)PRDGS傳輸特性的影響。在矩形單元尺寸一定時(shí)，隨著介質(zhì)的介電常數(shù)不斷增大，阻帶中心頻率會(huì)整體向左偏移，而阻帶的寬度和深度基本保持不變。在介電常數(shù)不變的情況下，PRDGS的阻帶寬度和阻帶深度，會(huì)隨著矩形單元尺寸a，b的增大而增大，同時(shí)阻帶中心頻率基本保持不變。將兩種方法所計(jì)算出的阻帶中心頻率進(jìn)行比較，一致性很好，平均相對(duì)誤差小于1.72 %。通過分析CN-FDTD時(shí)間步長(zhǎng)與效率和精度的關(guān)系，得出當(dāng)計(jì)算時(shí)間步長(zhǎng)遠(yuǎn)大于CFL法的時(shí)間步長(zhǎng)時(shí)，其計(jì)算結(jié)果與傳統(tǒng)FDTD計(jì)算結(jié)果仍然吻合，同時(shí)計(jì)算效率能提高77.2 %。這些有益的結(jié)論，能為PRDGS的精確計(jì)算提供一定指導(dǎo)。