王宏建 劉世華,2 郝齊焱,2 易 敏 陳 雪 劉 廣
(1.中國科學(xué)院空間科學(xué)與應(yīng)用研究中心微波遙感技術(shù)重點實驗室,北京 100190;2.中國科學(xué)院研究生院,北京 100049)
波導(dǎo)縫隙陣列天線以其結(jié)構(gòu)緊湊、重量輕、增益高、副瓣低、口徑利用率高、口徑場分布易于控制等優(yōu)點,廣泛應(yīng)用于雷達(dá)、通信等領(lǐng)域。通常縫隙陣列的設(shè)計可分為三個步驟[1]:① 按天線電性能要求,確定波導(dǎo)尺寸、面陣元數(shù)、各縫上的電壓;② 結(jié)合縫隙特性曲線,考慮各個縫間互耦,由Elliott設(shè)計方程調(diào)整縫隙長度及偏置量;③ 重復(fù)步驟②直到陣中縫處于諧振狀態(tài)并滿足所要求的電壓分布。關(guān)于縫隙特性的描述,已有眾多的學(xué)者采用不同的方法[2]研究過,近年時域有限差分(FDTD)等技術(shù)[3-4]也被用來精確描述縫隙特性。
商業(yè)電磁軟件以其快速、高效率的特點,被越來越多地用于波導(dǎo)縫隙天線的設(shè)計中。以基于有限元法的高頻電磁仿真軟件Ansoft HFSS為例,齊美清等[5]采用它得到開縫波導(dǎo)的散射參數(shù),提取了輻射縫隙的自導(dǎo)納和有源導(dǎo)納,有效地提高了縫隙天線的設(shè)計效率;金劍等[6]利用HFSS的S參數(shù)仿真結(jié)果,采用等效網(wǎng)絡(luò)法對互耦環(huán)境下縫隙的導(dǎo)納參數(shù)進(jìn)行了計算;趙懷成等[7]采用數(shù)值逼近法直接計算陣列中的縫隙歸一化導(dǎo)納,然后用HFSS驗證了設(shè)計方法的正確性;還有人用HFSS的參數(shù)掃描及優(yōu)化功能設(shè)計了波導(dǎo)窄邊縫隙陣天線[8-9],并通過試驗樣機(jī)作了驗證。而Elliott提出的有源導(dǎo)納法為波導(dǎo)裂縫陣天線的輻射陣面設(shè)計提供了理論基礎(chǔ),該方法考慮了輻射縫隙間互耦的影響[10]和波導(dǎo)內(nèi)部高次模的影響[11],可用于包括電介質(zhì)填充的波導(dǎo)等陣列天線的設(shè)計[12]。
在星載Ku波段縫隙陣列天線設(shè)計中,由于電性能要求和結(jié)構(gòu)空間限制,選定半高的非標(biāo)準(zhǔn)波導(dǎo)。考慮到半高波導(dǎo)中更容易產(chǎn)生高次模影響的特點,以及設(shè)計的成本和效率問題,文章將HFSS與Fortran編程相結(jié)合應(yīng)用于該類天線的設(shè)計。在對HFSS進(jìn)行準(zhǔn)確性驗證的基礎(chǔ)上,通過HFSS仿真分析縫隙模型,獲得了輻射縫隙和耦合縫隙的特性。然后結(jié)合Elliott提出的設(shè)計方程,用Fortran編程計算了8×10波導(dǎo)縫隙陣列天線的縫隙參數(shù),并用HFSS優(yōu)化天線陣整體結(jié)構(gòu)參數(shù),最后加工了8×10波導(dǎo)縫隙陣列天線。測試結(jié)果與仿真結(jié)果基本吻合,均滿足要求,良好的測試結(jié)果驗證了用HFSS獲取縫隙特性曲線的正確性及設(shè)計方法的有效性。
采用諧振式波導(dǎo)縫隙陣列天線,每個縫隙的中心間距取λg/2(λg代表波導(dǎo)波長),為了實現(xiàn)高效率,分支波導(dǎo)終端短路,形成駐波陣。各個縫隙的理想工作狀態(tài)是諧振,其長度參數(shù)均在λ0/2(λ0代表自由空間波長)左右。在實際工作前,特別針對WR90型標(biāo)準(zhǔn)波導(dǎo),研究了f0=9.375GHz時的縫隙特性,將HFSS仿真結(jié)果與Stegen[13]及L.G.Josefsson[14]的實測結(jié)果進(jìn)行了比較,良好的結(jié)果驗證了HFSS在該類應(yīng)用中的準(zhǔn)確性。
1.1.1 寬邊縱向輻射縫隙
寬邊縱向輻射縫隙模型如圖1所示,縫隙的電參數(shù)包括:縫隙寬度w,縫隙高度t,縫隙長度lr,縫隙中心偏離波導(dǎo)中心線的距離x.綜合考慮到加工難度、天線工作帶寬及微放電等因素,在設(shè)計初已分析并確定了縫隙寬度w和高度t的具體值。該段波導(dǎo)用空氣填充,總長為3λg/4,距離縫隙中心λg/4處是波導(dǎo)的終端短路板,另一端距離縫隙中心λg/2處是激勵波端口。
圖1 寬邊縱向輻射縫隙模型
圖2 縫隙諧振長度與偏置的關(guān)系
在波端口處的Y參數(shù)可以等效為距端口λg/2處的縫隙中心的Y參數(shù),根據(jù)波導(dǎo)縫隙的基本設(shè)計理論[2],在諧振時縱向縫隙的導(dǎo)納為實數(shù),即Im(Y)=0.給定縫隙偏置,通過計算不同縫隙長度下的參數(shù)ActiveY(有源導(dǎo)納),可分析出縫隙的諧振長度,再結(jié)合波阻抗Z0,就可計算出不同偏置下諧振時的歸一化輸入電導(dǎo)。圖2是縫隙諧振長度與偏置的關(guān)系曲線,實線表示縫隙終端為矩形頭時的變化關(guān)系,虛線表示縫隙終端為圓形頭時的變化關(guān)系,分析結(jié)果發(fā)現(xiàn)二者在數(shù)值上近似為平移關(guān)系,在最后加工天線時,可根據(jù)該結(jié)果把矩形頭輻射縫隙修正為圓形頭縫隙。另外,與文獻(xiàn)[3,13]中的結(jié)果相比,圖2中曲線變化更陡一些,這是由于半高非標(biāo)準(zhǔn)波導(dǎo)的高度小,更容易激勵起高次模,縫隙的諧振長度對其偏置變化更加敏感??p隙諧振長度與偏置的關(guān)系可表示為[15]
圖3是縫隙的歸一化諧振電導(dǎo)與偏置的關(guān)系,將HFSS仿真結(jié)果與Stevenson經(jīng)典公式[2]的計算結(jié)果進(jìn)行比較,在偏置較小時,二者吻合較好,隨著偏置的增大,二者逐漸偏離,因為Stevenson的公式在縫隙偏置為0時的情況下結(jié)果準(zhǔn)確,偏置越大,誤差就會越大[15]??p隙歸一化諧振電導(dǎo)與偏置的關(guān)系[15]表示為
在縫隙陣列中,由于輻射縫隙之間的互耦存在,如果使用單個縫隙諧振長度時的值,縫隙的有源阻抗并不是純實數(shù),這就涉及到縫隙長度在非諧振時的導(dǎo)納。Stegen[13]在對波導(dǎo)寬邊縱向縫隙的研究中發(fā)現(xiàn),以縫隙電導(dǎo)和電納與諧振電導(dǎo)之比為縱坐標(biāo),縫隙長度與諧振長度之比為橫坐標(biāo),在縫隙不同偏置條件下得到的歸一化曲線形狀十分相似。為了獲得半高非標(biāo)準(zhǔn)波導(dǎo)中縫隙的Stegen曲線,在不同偏置下,結(jié)合圖2結(jié)果,對0.9lr≤l≤1.1lr范圍內(nèi)的縫隙長度進(jìn)行計算,得到其有源導(dǎo)納并歸一化處理,結(jié)果如圖4所示,左右縱軸分別是歸一化電導(dǎo)和歸一化電納,這里歸一化對象是縫隙的諧振電導(dǎo),而圖3的歸一化對象是波導(dǎo)的波阻抗。圖中擬合曲線由不同偏置下的結(jié)果取平均得到。該曲線可用函數(shù)[15]h和變量y 表示為
圖3 縫隙歸一化諧振電導(dǎo)與偏置的關(guān)系
式中:橫坐標(biāo)y=l/lr,表示縫隙長度和諧振縫隙長度之比;實部h1(y)=G(x,y)/Gr(x,y=1);虛部h2(y)=B(x,y)/Gr(x,y=1).
圖4 縫隙歸一化導(dǎo)納與歸一化長度的關(guān)系
1.1.2 中心傾斜耦合縫隙
寬邊中心傾斜耦合縫隙模型如圖5所示,包括兩個長度均為λg的正交矩形波導(dǎo),中間由耦合縫隙相連接,該模型將波導(dǎo)壁厚度近鄰為零,只取波導(dǎo)內(nèi)腔結(jié)構(gòu),有益于節(jié)省計算時間。耦合縫隙的電參數(shù)包括:縫隙寬度w,縫隙高度t,縫隙長度lr,縫隙長度方向中心線偏離波導(dǎo)中心線的角度θ.用HFSS分析時,兩根波導(dǎo)共四個端口均設(shè)定為激勵端口,以實現(xiàn)終端阻抗匹配。
圖5 寬邊中心傾斜耦合縫隙模型
S.R.Rengarajan定義耦合縫隙的諧振條件[16-17]為縫隙處前向散射波與入射TE10波反相,該條件等同于[3]:在頻率和傾角一定時,以縫隙長度為變量,諧振時耦合波導(dǎo)向分支波導(dǎo)耦合的能量最大。在圖5模型中,可以通過傳輸參數(shù)S31和S41[3]來判斷耦合能量的多少。在45°≤θ≤90°范圍內(nèi)的傳輸參數(shù)與0°≤θ≤45°內(nèi)的傳輸參數(shù)幅度相同,相位相反[18],因此對0°到45°內(nèi)的縫隙特性作分析即可。圖6、圖7分別是諧振長度和傳輸參數(shù)S31、S41、S11與傾斜角的關(guān)系。
圖6 諧振長度隨傾斜角的變化曲線
從圖6可見,耦合縫隙諧振長度隨傾斜角的變化很小,由于是半高非標(biāo)準(zhǔn)波導(dǎo),與文獻(xiàn)[13]中的標(biāo)準(zhǔn)波導(dǎo)相比,該曲線相對陡一些。圖7中,隨著傾斜角的增大,從端口1耦合到端口3、4的能量逐漸增多,在θ=45°時最大,實際耦合波導(dǎo)中,多個耦合縫隙以λg/2為間距排列,傳輸參數(shù)還會受縫隙間的互耦影響。
圖7 傳輸參數(shù)隨傾斜角的變化曲線
Elliott提出的平面縫隙陣列設(shè)計方程[8-10]為
式中:Yna為單個縫隙的有源導(dǎo)納;Vns為縫隙電壓;Vn為模式電壓;gmn為第m 個和第n個縫隙間的互耦電導(dǎo)計算公式[9]。該設(shè)計方法考慮了外部縫隙間的互耦影響。在應(yīng)用中,第一設(shè)計方程寫為
第二個設(shè)計方程,考慮了波導(dǎo)內(nèi)部高次模的影響,特別有益于非標(biāo)準(zhǔn)波導(dǎo)中的縫隙陣設(shè)計。分母中第一項是自諧振分量,第二項為外耦合分量,要滿足各縫隙的有源導(dǎo)納諧振,即為純實數(shù),應(yīng)該滿足
還有輸入導(dǎo)納匹配方程,由于采用了波導(dǎo)中心耦合饋電,歸一化有源導(dǎo)納之和為2,即
根據(jù)Elliott設(shè)計方程(4)~(8)和縫隙特性曲線(圖2~4、6、7),通過Fortran編程,設(shè)計了波導(dǎo)饋電的8×10平面縫隙陣列。該陣列天線包括3部分:輻射陣面(由若干條分支波導(dǎo)組成)、與分支波導(dǎo)交叉放置的耦合波導(dǎo)、為耦合波導(dǎo)饋電的饋電網(wǎng)絡(luò)。該陣列天線的輻射陣面分成四個區(qū),用四條耦合波導(dǎo)耦合饋電,通過1~4等幅同相波導(dǎo)功分網(wǎng)絡(luò)對耦合波導(dǎo)饋電。天線實物如圖8所示。
圖8 8×10天線實物圖
用HFSS仿真了8×10陣列天線的電性能。方向圖如圖9所示,E面和H面的副瓣電平分別是-24.2dB和-29.2dB;主波束有良好的對稱性,3 dB波束寬度為10.9°.駐波如圖10所示,中心頻點處的駐波為1.04,由于設(shè)計中采用了頻帶展寬措施,中心頻率附近的駐波曲線變化平坦,1.3以下帶寬為420MHz,具有較寬的頻帶特性。
圖9 8×10陣列天線方向圖仿真結(jié)果
在中科院空間科學(xué)與應(yīng)用研究中心的天線近場測試系統(tǒng)中完成了電性能測試。方向圖如圖11所示,E面和H面副瓣電平分別是-29dB和-27.5 dB,與仿真結(jié)果相比,E面實測副瓣電平較好,H面實測副瓣電平提升了1.2dB;主波束對稱性較好,3 dB波束寬度為10°.駐波如圖12所示,中心頻點處的駐波為1.1,駐波在1.3以下的帶寬為340MHz,與仿真結(jié)果相比,中心頻率處的駐波高了0.06,1.3以下的駐波帶寬比仿真結(jié)果少了80MHz.
副瓣電平、駐波、波束寬度等的實測結(jié)果與仿真結(jié)果有偏差,主要是加工誤差和測試誤差的影響所致。諧振陣列工作在中心頻率時,理論上縫隙的有源導(dǎo)納是實數(shù),天線的口徑場相位面為同相。當(dāng)工作頻率變化時,產(chǎn)生以下效應(yīng):縫隙的有源導(dǎo)納特性產(chǎn)生變化;駐波峰值位置偏離縫隙中心;縫隙間距偏離λg/2;縫隙終端所接的λg/4長的短路板,在工作頻率變化時產(chǎn)生一定的電納;縫隙間的互耦發(fā)生變化。而加工誤差會引起中心頻率的偏移,受以上諸因素影響,天線的匹配條件變壞,口徑場的幅度、相位分布變壞,天線的副瓣電平抬高、增益下降,也會提高天線的輸入駐波比。
整體上實測結(jié)果與仿真結(jié)果基本一致,均滿足設(shè)計指標(biāo)要求。
結(jié)合天線實際需求將HFSS與Fortran編程相結(jié)合用于星載Ku波段波導(dǎo)縫隙陣。通過HFSS給出了半高非標(biāo)準(zhǔn)波導(dǎo)中的寬邊輻射縫隙和耦合縫隙的特性曲線,結(jié)合Elliott提出的平面縫隙陣設(shè)計方程,考慮了輻射縫隙間互耦因素及波導(dǎo)內(nèi)部高次模影響,用Fortran編程計算了8×10波導(dǎo)縫隙陣列天線的縫隙參數(shù)。天線近場測試結(jié)果表明:天線副瓣電平、駐波、主波束寬度等和仿真結(jié)果基本一致,均滿足設(shè)計指標(biāo)。與工程使用的實驗法比較,這種方法克服了傳統(tǒng)方法中為了減少試驗件加工工作量帶來的誤差、節(jié)約了設(shè)計成本、提高了工作效率,同時還考慮了非標(biāo)準(zhǔn)波導(dǎo)易產(chǎn)生高次模影響的特點。良好的測試結(jié)果驗證了用HFSS獲取縫隙特性曲線的正確性及設(shè)計方法的有效性,可用于后期大型陣列的設(shè)計。
[1]王宏建,高本慶,劉瑞祥,等.遺傳算法在波導(dǎo)裂縫陣列天線設(shè)計中的應(yīng)用[J].電子學(xué)報,2002,30(3):309-312.WANG Hongjian,GAO Benqing,LIU Ruixiang,et al.The application of genetic algorithms(gas)on the design of the waveguide slotted array[J].Acta Electronica Sinica,2002,30(3):309-312.(in Chinese)
[2]林昌祿,聶在平.天線工程手冊[M].北京:電子工業(yè)出版社,2002:273-298.
[3]WU Ren,GAO Benqing.Full wave analysis of broad wall slot's characteristics in rectangular waveguides[J].IEEE Transactions on Antennas and Propagation,2004,52(9):2436-2444.
[4]丁曉磊,王 建,林昌祿.耦合縱縫饋電的短路波導(dǎo)寬邊縱縫陣的分析與設(shè)計[J].電波科學(xué)學(xué)報,2001,17(4):422-426.DING Xiaolei,WANG Jian,LIN Changlu.The analysis and synthesis of a longitudinal slot-fed array in a broad wall of a shorted waveguide[J].Chinese Journal of Radio Science,2001,17(4):422-426.(in Chinese)
[5]齊美清,汪 偉,金謀平.基于HFSS的波導(dǎo)裂縫有源導(dǎo)納的計算方法[J].雷達(dá)科學(xué)與技術(shù),2006,4(2):121-124.QI Meiqing,WANG Wei,JIN Mouping.Computation method for active admittance of waveguide slot based on HFSS[J].Radar Science and Technology,2006,4(2):121-124.(in Chinese)
[6]金 劍,萬笑梅,汪 偉,等.波導(dǎo)平板裂縫天線陣的設(shè)計[J].雷達(dá)科學(xué)與技術(shù),2007,5(3):232-235.JIN Jian,WAN Xiaomei,WANG Wei,et al.Design of plannar waveguide slotted antenna arrays[J].Radar Science and Technology,2007,5(3):232-235.(in Chinese)
[7]趙懷成,吳 文,湯一銘,等.小型波導(dǎo)縫隙天線陣的改進(jìn)設(shè)計方法研究[J].電波科學(xué)學(xué)報,2007,22(6):1025-1028.ZHAO Huai cheng,WU Wen,TANG Yi ming,et al.Improvement of design method for small rectangular waveguide slot arrays[J].Chinese Journal of Radio Science,2007,22(6):1025-1028.
[8]武 偉,路志勇.波導(dǎo)窄邊縫隙陣列天線的分析與設(shè)計[J].無線電工程,2009,39(1):42-44.WU Wei,LU Zhiyong.Analysis and design of edge slot waveguide array antennas[J].Radio Engineering of China,2009,39(1):42-44.(in Chinese)
[9]李知新,任 朗.矩形波導(dǎo)窄邊縫隙分析方法述評[J].電波科學(xué)學(xué)報,1999,14(2):235-241.LI Zhixin,REN Lang.Review of analysis method of the inclined slots on the narrow wall of a rectangular waveguide[J].Chinese Journal of Radio Science,1999,14(2):235-241.(in Chinese)
[10]ELLIOTT R S,KURTZ L A.The design of small slot arrays[J].IEEE Transactions on Antennas and Propagation,1978,26(2):214-219.
[11]ELLIOTT R S,O'LOUGHLIN W R.The design of slot arrays including internal mutual coupling[J].IEEE Transactions on Antennas and Propagation,1986,34(9):1149-1154.
[12]Elliott R S.An improved design procedure for small slot arrays of shunt slots[J].IEEE Transactions on Antennas and Propagation,1983,31(1):48-53.
[13]JASIK H.Antenna Engineering Handbook[M].2nd ed.New York:McGraw-Hill,1989.
[14]JOSEFSSON L G.Analysis of longitudinal slots in rectangular waveguides[J].IEEE Transactions on Antennas and Propagation,1987,35(12):1351-1357.
[15]ELLIOTT R S.Antenna Theory and Design[M].Englewood Cliffs,NJ:Prentice-Hall,1981:408-411.
[16]RENGARAJAN S R.Characteristics of a longitudinal/transverse coupling slot in crossed rectangular waveguides[J].IEEE Microwave Theory and Techniques Society,1989,37(8):1171-1177.
[17]RENGARAJAN S R.Compound coupling slots for arbitrary excitation of waveguide-fed planar slot arrays[J].IEEE Microwave Theory and Techniques Society,1990,38(2):276-280.
[18]RENGARAJAN S R.Analysis of a centered-inclined waveguide slot coupler[J].IEEE Microwave Theory and Techniques Society,1989,37(5):884-889.