段佳艷， 鄭勤紅， 翟實慶， 姚斌， 向泰， 曹湘琪， 楊繼孔

（云南師范大學 物理與電子信息學院，云南 昆明650092）

自1947年S.B.Cohn［1］提出脊波導的概念以來，人們對脊波導的研究不斷深入.與普通波導相比，脊波導具有單模帶寬更寬、主模截止波長更長等優(yōu)點，因此脊波導被廣泛應用在微波和毫米波器件中.近年來，隨著對脊波導不斷的深入研究，對更多不同形狀脊波導的研究不斷增加.文獻［2］用積分方程法計算了非對稱單脊波導的各種參數(shù)；文獻［3］用FDTD法分析了非對稱單脊波導，獲得了模場分布圖、波導波長和輻射特性；文獻［4］用多極理論分析了脊波導族的本征問題；文獻［5］用有限元法分析了倒梯形和三角形雙脊波導的TE模的傳輸特性；文獻［6］用差分法分析了不對稱雙脊波導的截止特性；文獻［7］用差分法分析了山字形脊波導的傳輸特性.近年來，新的計算方法和現(xiàn)代高性能計算機的發(fā)展為波導的研究提供了新的平臺，本文用FDTD法分析山字形脊波導的截止特性，并與文獻報道結(jié)果進行比較，分析了部分介質(zhì)填充山字形脊波導的截止特性.

1 基本理論與方法

電磁場所滿足的麥克斯韋旋度方程為：

對于二維問題，?／?z＝0，對于TE模，對 （1）式作時域差分離散，假設(shè)＝0，＝0，電磁參數(shù)?。剑娇傻?/p>

（3）、（4）、（5）三式即為各向同性線性介質(zhì)填充波導的FDTD迭代公式.

對于各向異性介質(zhì)，其電磁參數(shù)為：

對 （1）、（2）式作時域差分離散，可得如下差分公式［8－9］

由 （6）、（7）可獲得電場分量和磁場分量的FDTD迭代計算結(jié)果，詳見文獻［9］、［10］.

2 計算實例

2.1 山字形脊波導截止特性

圖1為山字形脊波導的截面示意圖，主脊位于波導中間，左右為副脊，且主脊和副脊的寬度均為s.為了驗證FDTD法的有效性，波導采用和文獻［7］相同的尺寸，即b／a＝ 0.8，（其中，d1≥d2，d2／b＝0.1，s／a＝0.1） ，并采用如下所示的激勵源：

（8）式中t0為脈沖到達峰值的時間，T決定高斯脈沖的寬度，T＝15dt，t0＝140dt，dt為時間步長.表1列出了用FDTD法獲得的各種幾何尺寸條件下的主模歸一化截止波長，并與文獻［7］用差分法（FDM）獲得的計算結(jié)果進行了比較，由表1可見兩種方法獲得的結(jié)果吻合較好.

表1 山字形脊波導的歸一化截止波長Table 1 The normalized cut－off wavelength of SHAN ZI ridge waveguide

2.2 部分介質(zhì)填充山字形脊波導

圖2為介質(zhì)填充山子形脊波導的截面示意圖，陰影部分為介質(zhì)填充部分.s為介質(zhì)厚度，＝εrε0，d2／b＝0.4.將波導劃分為200×160的網(wǎng)格，計算結(jié)果見表2.

圖1 山字形脊波導截面圖Fig.1 The section view of SHAN ZI ridge waveguide

圖2 部分介質(zhì)填充山字形脊波導截面圖Fig.2 The section view of the SHAN ZI ridge waveguide partially filled with medium

表2 εr取不同值時山子形脊波導的歸一化截止波長Table 2 The normalized cut－off wavelength of the SHAN ZI ridge waveguide filled with different dielectric（εr）

由表2可知：當εr保持不變時，歸一化截止波長隨著d1／b的增大而增大；當d1／b保持不變時，歸一化截止波長隨著εr的增大而增大.圖3給出了山子形脊波導的脊寬度變化時對歸一化截止波長的影響，由圖可知：d1／b不變，s／a越大，即脊寬度越大，歸一化截止波長越大；s／a不變，歸一化截止波長隨著d1／b的增大而增大.

圖3 TE波主模截止波長與s／a的關(guān)系Fig.3 The relation of normalized cut－off wavelength of TE－mode vs s／a

圖4 TE波主模截止波長與加載不同介質(zhì)的關(guān)系Fig.4 The relation of normalized cut－off wavelength of TE－mode vs different loaded medium

圖4 給出了山字形脊波導在填充不同介質(zhì)時，波導TE波主模的截止特性.其中，d2／b＝0.4，s／a＝0.1.各向同性線性均勻介質(zhì)取值為：εr＝9.6，單軸各向異性介質(zhì)的介電參數(shù)張量取值為：，非單軸各向異性介質(zhì)的介電參數(shù)張量取值為

由圖4可知：

（1）填充介質(zhì)后，歸一化截止波長較未填充介質(zhì)時大.

（2）無論有無填充介質(zhì)，λc／a隨著d1／b的增大而增大.

（3）當d2／b，d1／b保持不變，介質(zhì)填充山字形脊波導時，單軸各向異性介質(zhì)填充波導的歸一化截止波長最大，各向同性均勻介質(zhì)最小，非單軸各向異性介質(zhì)介于單軸各向異性介質(zhì)和各向同性線性均勻介質(zhì)之間.

3 結(jié) 論

本文用FDTD法分析了山字形脊波導的截止特性，并與其他方法計算的結(jié)果做了比較，驗證了本文方法的可行性.其次，計算部分介質(zhì)填充山字形脊波導的截止特性，計算和分析了脊寬度改變時對波導截止特性的影響.在此基礎(chǔ)上，文章還分析了波導在εr取不同的值時和填充不同的介質(zhì)時，山字形脊波導的截止特性.

［1］COHN S B.Properties of ridge waveguide［J］.Proc IRE，1974，35（8）：783－788.

［2］魯加國，樊德森.非對稱單脊波導的積分方程法分析［J］.微波學報，1998，14（2）：108－115.

［3］韓一平，聶在平.非對稱單脊波導的FDTD分析［J］.電子學報，1999，27（3）：126－128.

［4］鄭勤紅，曾華，解福瑤，等.脊波導族的多極理論分析［J］.微波學報，2001，17（3）：24－28.

［5］陳小強，李明，逯邁，等.兩種新型雙脊波導傳輸特性的研究［J］.西安電子科技大學學報，2007，34（3）：495－499.

［6］章寶歌，李根，王梓丞.不對稱雙脊波導的特性計算［J］.蘭州交通大學學報，2009，28（4）：102－104.

［7］韓國軍，王偉，張黎明，等.山字脊波導傳輸特性分析［J］.重慶工學院學報：自然科學版，2008，22（7）：91－95.

［8］楊利霞.復雜介質(zhì)電磁散射的FDTD算法及相關(guān)技術(shù)研究［D］.西安：西安電子科技大學，2006.

［9］王建敏，鄭勤紅，姚斌，等.用FDTD法分析各向異性介質(zhì)填充波導和微帶線的傳輸特性［J］.云南師范大學學報：自然科學版，2011，31（2）：60－64.

［10］張方濤，鄭勤紅，向泰，等.各向異性介質(zhì)填充矩形波導的FDTD分析［J］.云南師范大學學報：自然科學版，2012，32（3）：54－57.