顧長青 王 斌 牛臻弋 李 茁 史劍峰 林志斌 曹麗霞

（南京航空航天大學(xué)電子信息工程學(xué)院，江蘇 南京210016）

引 言

在應(yīng)用功率平衡法的網(wǎng)絡(luò)化公式分析電大復(fù)雜結(jié)構(gòu)混沌腔體的屏蔽性能［1］中，精確預(yù)知孔縫和腔內(nèi)裝載物的耦合截面積是關(guān)鍵.其中腔內(nèi)壁、孔縫、收發(fā)天線等的耦合截面積都有近似解析式或半解析式，而腔內(nèi)一般裝載物的耦合截面積則以混響室實(shí)驗(yàn)測定為主.近年來，在混響室電磁模擬環(huán)境下，提取腔內(nèi)裝載物耦合截面積的研究已取得一定進(jìn)展.

混響室是一種模擬現(xiàn)實(shí)中漫射場電磁環(huán)境的腔體結(jié)構(gòu)，已廣泛應(yīng)用于電磁輻射耦合、抗擾度、屏蔽效能、天線效率和吸收截面等的測試研究.由于是通過模式攪拌器攪拌來實(shí)現(xiàn)多模和場均勻，因此待測量具有統(tǒng)計(jì)性質(zhì).D.A.Hill在20世紀(jì)90年代開展了混響室概率統(tǒng)計(jì)參數(shù)的研究，提出了一種新的模擬混響室中電磁場的方法，即混響室中的場可以看成是自由空間中平面波角譜的積分［2-3］；2003年L.Musso通過平面波角譜的離散化和蒙特卡羅方法數(shù)值模擬了混響室中的場分布，結(jié)合傳輸線理論計(jì)算了傳輸線上的平均耦合電流［4-5］；2009年I.Junqua研究了混響室中電場的空間相關(guān)函數(shù)與傳輸線上耦合電流平方的均值之間的關(guān)系，推導(dǎo)出其半解析式，并計(jì)算了傳輸線的平均耦合截面積［6-7］.2012年M.Magdowski又進(jìn)一步推導(dǎo)出傳輸線的平均耦合截面積的一般閉合式，給出了負(fù)載匹配、開路和短路等特殊情況的平均耦合截面積的簡化式［8］.2012年H.Zhang對(duì)I.Junqua文獻(xiàn)中的半解析式進(jìn)行完善，包含了兩個(gè)正交場分量間的相關(guān)函數(shù)項(xiàng)的計(jì)算［9］.此外，他通過改變I.Junqua文獻(xiàn)中半解析式的積分順序，獲得文獻(xiàn)［6-7］同樣問題的另一種形式的半解析式［10］.

在使用傳輸線理論（BLT方程）研究混響室電磁環(huán)境下場線間的耦合問題中［4-10］，由于BLT方程只考慮了傳輸線的差模響應(yīng)并忽略了高次模影響，因此高頻計(jì)算結(jié)果不再精確.使用全波仿真方法進(jìn)行直接計(jì)算，雖然計(jì)算結(jié)果精確，但是需要計(jì)算大量不同入射方向、極化角的平面波對(duì)傳輸線的響應(yīng)，計(jì)算效率低下.

本文使用平面波角譜概念和蒙特卡羅方法來數(shù)值模擬混響室中的電磁環(huán)境，在此基礎(chǔ)上，將理想混響室內(nèi)靠近理想導(dǎo)電地平面的多導(dǎo)體傳輸線等效成多天線系統(tǒng)，首次應(yīng)用文獻(xiàn)［12］提出的基于互易定理的快速全波分析方法提取了耦合電流模值平方的均值和平均耦合截面積，驗(yàn)證了方法的正確性，并計(jì)算了單根裸線傳輸線、帶包層傳輸線、雙絞傳輸線、傳輸線附近有導(dǎo)體擋板以及多導(dǎo)體傳輸線端口的平均耦合截面積.

1 模擬混響室電磁環(huán)境的蒙特卡羅方法

根據(jù)Hill的理論，理想狀態(tài)下混響室中無源位置處的電場可以由平面波角譜對(duì)空間立體角的積分得到，即

式中：角譜F為復(fù)矢量，每個(gè)分量包含有幅度和相位信息；k為自由空間波矢量；r為觀察點(diǎn)位置矢量.根據(jù)混響室的統(tǒng)計(jì)理論；角譜F是依賴于模式攪拌器旋轉(zhuǎn)位置的隨機(jī)變量.

平面波角譜離散化后，式（1）進(jìn)一步表示為

式（2）表明理想混響室中r處的電場可以選用Ⅰ個(gè)隨機(jī)平面波的疊加進(jìn)行表述，這樣無需再考慮理想混響室的具體結(jié)構(gòu).

在蒙特卡羅方法模擬混響室內(nèi)的場分布時(shí)，需要隨機(jī)地抽取入射平面波電場的相關(guān)物理量，包括入射角（方位角、俯仰角）、極化角和相位.由于混響室是各向同性的環(huán)境，因此對(duì)影響入射平面波的各參數(shù)在其變化范圍內(nèi)可以等概率地選取，也就是說入射平面波的空間立體角在0～4π范圍內(nèi)均勻分布，極化角和相位在0～2π范圍內(nèi)均勻分布.

為了預(yù)測混響室中傳輸線上耦合的電磁能量大小，首先要求出表征混響室攪拌器在某個(gè)位置時(shí)Ⅰ列隨機(jī)平面波在傳輸線負(fù)載上總的耦合電流，然后對(duì)混響室攪拌器J個(gè)不同旋轉(zhuǎn)角度下耦合電流模值進(jìn)行平方后再取均值，得到和能量相關(guān)的量，即耦合電流模值平方的均值為

式中：符號(hào)〈·〉表示系統(tǒng)平均，即在J個(gè)不同旋轉(zhuǎn)角度下耦合電流的算術(shù)平均值；Ⅰi表示第i平面波在傳輸線負(fù)載上產(chǎn)生的耦合電流.

傳輸線平均耦合截面積σ定義為：

式中：RL為傳輸線端口負(fù)載阻抗的實(shí)部——電阻；E0為電場均方根值，一般取作1V／m.

2 互易定理求耦合電流

由第1節(jié)可以看出，在蒙特卡羅方法計(jì)算混響室環(huán)境下傳輸線平均耦合截面積中，總的計(jì)算量主要由〈｜Ⅰ｜2〉的計(jì)算量決定，它需要Ⅰ×J次計(jì)算傳輸線耦合電流Ⅰi，計(jì)算時(shí)間開銷大約為O（Ⅰ×J）.采用傳統(tǒng)的BLT方程計(jì)算Ⅰi，雖然計(jì)算速度較快，但是存在普適性差、高頻時(shí)誤差變大的問題.本文基于文獻(xiàn)［12-13］提出的互易定理方法實(shí)現(xiàn)〈｜Ⅰ｜2〉的快速求解，既避開了空間相關(guān)函數(shù)的求解，又避免了直接做Ⅰ×J次的全波計(jì)算.

對(duì)于理想混響室內(nèi)靠近理想導(dǎo)電地平面的多導(dǎo)體傳輸線，可以用圖1多導(dǎo)體傳輸線的結(jié)構(gòu)示意圖無限大理想導(dǎo)電地平面上多導(dǎo)體傳輸線來表示，其中考慮了地平面的作用.基于多端口網(wǎng)絡(luò)理論，多導(dǎo)體傳輸線可以等效成N根天線系統(tǒng)的多端口網(wǎng)絡(luò)，傳輸線的每一負(fù)載端接處等效為一天線的激勵(lì)端口，這樣所有負(fù)載和激勵(lì)源可以統(tǒng)一地看作是與多天線系統(tǒng)連接的另一多端口網(wǎng)絡(luò).

圖1 多導(dǎo)體傳輸線的結(jié)構(gòu)示意圖

計(jì)算過程可歸納如下：

1）第n天線工作在發(fā)射模式，其它天線端口開路時(shí)，使用全波方法（例如矩量法）計(jì)算輻射電場以及其余第m（m＝1，2，…，N）天線端口上的電壓電流

2）根據(jù)公式

計(jì)算第n天線工作在接收模式時(shí)由蒙特卡羅方法抽取的隨機(jī)入射平面電磁波Einc（θ，φ，ψ）在饋電端口感應(yīng)的開路電壓（θ，φ，ψ）；式中：λ為工作波長；η為自由空間的特性阻抗；（θ，φ）為球坐標(biāo)系中方向角；ψ為極化角.根據(jù)公式

計(jì)算第n天線與第m（n、m＝1，2，…，N）天線之間的端口互阻抗Znm，式中為第n饋電端口上加激勵(lì)時(shí)其端口上的電流.

3）對(duì)步驟1）和步驟2）進(jìn)行n從1至N的循環(huán)運(yùn)算，得到N天線系統(tǒng)的激勵(lì)電壓向量Vopen（θ，φ，ψ）和端口的阻抗矩陣Zant.Vopen向量和Zant矩陣的形式為：

由互易定理知，Zant與天線工作模式無關(guān)，不論工作在接收模式還是工作在發(fā)射模式，其阻抗元素是一樣的，另外Znm＝Zmn.

4）計(jì)算隨機(jī)均勻平面波在傳輸線負(fù)載上的耦合電流響應(yīng).天線端口上的電流向量Ir和電壓向量Vr通過下列式子計(jì)算得到：

5）做兩重循環(huán)運(yùn)算.內(nèi)循環(huán)對(duì)Ⅰ列蒙特卡羅隨機(jī)選取的平面波遍歷，外循環(huán)對(duì)J個(gè)旋轉(zhuǎn)角度遍歷，求出多導(dǎo)體傳輸線所有端口上的耦合電流模值平方的均值，進(jìn)而求得各個(gè)負(fù)載的平均耦合截面積.對(duì)于給定的一組負(fù)載，如果僅僅要求預(yù)測某一個(gè)負(fù)載的平均耦合截面積，只需要做一次全波仿真即可完成.

3 數(shù)值結(jié)果

當(dāng)單導(dǎo)體裸線傳輸線的長度L＝40cm，半徑r＝0.5mm，離地高度h＝5mm時(shí)，傳輸線的特性阻抗ZC為180Ω.圖2給出了傳輸線端口接匹配負(fù)載時(shí)本文方法、BLT方程結(jié)合蒙特卡羅方法和解析式分別計(jì)算得到的〈｜Ⅰ｜2〉曲線.其中，在本文方法中按照文獻(xiàn)［9］選?、瘛罦＝1 000×250，滿足2π／J?1、2π／Ⅰ?1條件.注意的是這里的平面波角譜是在半空間抽?。唤馕鍪剑?］是〈｜Ⅰ｜2〉＝（／2）（h／ZC）2［2-sin（2kL）／（kL）］.從圖中可以看出，三曲線吻合良好，驗(yàn)證了本文方法的正確性.強(qiáng)調(diào)的是：蒙特卡羅方法是一種概率的方法，因此每次計(jì)算所抽取的樣本都不一樣，相應(yīng)的結(jié)果有點(diǎn)差異，但變化不大.

圖3是單導(dǎo)體裸線傳輸線L＝1m，r＝0.5 mm，h＝5cm，＝50Ω，＝50Ω，參數(shù)條件下負(fù)載1的平均耦合截面積曲線（負(fù)載2是一樣的）.與文獻(xiàn)［7］的結(jié)果（圖中帶三角形線）相比，本文方法計(jì)算結(jié)果更接近傳輸線作為天線時(shí)［7］的耦合截面積.

圖2 耦合電流模值平方的均值

圖3 平均耦合截面積

圖4是帶包層單導(dǎo)體傳輸線L＝1m，r＝0.5 mm，h＝5cm，＝50Ω，包層厚度d＝0.3mm，相對(duì)介電常數(shù)εr＝2.5參數(shù)條件下2端口接不同負(fù)載時(shí)的平均耦合截面積曲線.和圖3曲線相比，＝50Ω時(shí)的結(jié)果與裸線傳輸線相比變化不是很大，增大負(fù)載則使端口耦合電流減小，降低了負(fù)載的耦合截面積.

圖4 改變開孔形狀

圖5 改變機(jī)箱大小

圖5給出了利用雙絞傳輸線L＝1m，r＝0.5 mm，h＝5cm，s＝0.9mm，p＝14mm＝50Ω，＝50Ω參數(shù)條件下負(fù)載1的平均耦合截面積曲線（雙絞線上任一點(diǎn)坐標(biāo)可由文獻(xiàn)［14］中公式確定）.從圖中可以看出雙絞線在低頻時(shí)能大大減少輻射干擾，而隨著頻率的升高，平均耦合截面積也增大，到10GHz時(shí)已經(jīng)與普通傳輸線相當(dāng)，甚至不如.因?yàn)樵诘皖l時(shí)外界電磁波在雙絞線上產(chǎn)生的共模信號(hào)受到雙絞線的平衡結(jié)構(gòu)抑制，而到了高頻雙絞線兩根導(dǎo)體之間的等效電容開始增大使其耦合增強(qiáng).

在單導(dǎo)體裸線傳輸線（L＝1m，r＝0.5mm，h＝5cm＝50Ω＝50Ω）附近存在一金屬擋板，擋板長邊與傳輸線平行并垂直放置在地平面上.該金屬擋板長Lm＝1.3m，與傳輸線之間距離dm＝5cm（傳輸線位于擋板正前方）.圖6給出了負(fù)載1的平均耦合截面積隨擋板高度hm的變化情況對(duì)比.hm≤h時(shí)沒有大幅降低傳輸線的耦合截面積，而h≤hm時(shí)，耦合截面積減小很多.

圖6 三個(gè)方孔情況

圖7為雙導(dǎo)體傳輸線，其中單根傳輸線參數(shù)L＝1m，r＝0.5mm，h＝5cm，兩根傳輸線之間距離為2mm，端口阻抗分別為＝50Ω＝100Ω，＝200Ω和＝500Ω時(shí)的各個(gè)端口耦合截面積曲線.

圖7 底面三根傳輸線上最大感應(yīng)電流

4 結(jié) 論

在混響室環(huán)境下，本文提出的傳輸線平均耦合截面積的全波快速提取方法，相比于BLT方法不僅計(jì)算結(jié)果更加精確，而且普適性強(qiáng).既適用于裸線、包層線、雙絞線等一般傳輸線，也適合于處于復(fù)雜環(huán)境中的傳輸線；既適合于單導(dǎo)體傳輸線，也適合于多導(dǎo)體傳輸線；傳輸線既可以靠近混響室的內(nèi)壁，也可以置在混響室的工作區(qū).此外，全波仿真過程是和端接負(fù)載無關(guān)的.

［1］JUNQUA I，PARMANTIER J P，ISSAC F.A network formulation of the power balance method for high-frequency coupling［J］.Electromagnetics，2005，25（7／8）：603-622.

［2］HILL D A.Electromagnetic Fields in Cavities：Deterministic and Statistical Theories［M］.John Wiley ＆Sons，2009.

［3］HILL D A.Plane wave integral representation for fields in reverberation chambers［J］.IEEE Transactions on Electromagnetic Compatibility，1998，40（3）：209-217.

［4］MUSSO L，BERAT V，CANAVERO F，et al.A plane wave Monte Carlo simulation method for reverberation chambers［J］.IEEE Trans on Electromagnetic Compatibility，1998，40（3）：209-217.

［5］MUSSO L.Assessment of Reverberation Chamber Testing for Automotive Applications［D］.Turin：University de Lille，2003.

［6］JUNQUA I，PARMANTIER J P，DEGAUQUE P.Field-to-wire coupling in an electrically large cavity：a semianalytic solution［J］.IEEE Transactions on Elec-tromagnetic Compatibility，2010，52（4）：1034-1040.

［7］JUNQUA I，PARMANTIER J P，DEGAUQUE P.Coupling on cables in an electrically large cavity：a theoretical approach［C］／／20th International Zurich Symposium on Electromagnetic Compatibility.IEEE，2009：89-92.

［8］MAGDOWSKI M，VICK R.Closed-form formulas for the stochastic electromagnetic field coupling to a transmission line with arbitrary loads［J］.IEEE Transactions on Electromagnetic Compatibility，2012，54（5）：1147-1152.

［9］ZHANG H，ZHAO X，YAN L，et al.Some amendments to“field-to-wire coupling in an electrically large cavity：a semianalytic solution”［J］.IEEE Transactions on Electromagnetic Compatibility，2012，54（1）：232-234.

［10］ZHANG H，ZHAO X，YAN L，et al.An alternative semianalytical／analytical solution to field-to-wire coupling in an electrically large cavity［J］.IEEE Transactions on Electromagnetic Compatibility，2012，54（5）：1153-1160.

［11］MAGDOWSKI M，TKACHENKO S V，VICK R.Coupling of stochastic electromagnetic fields to a transmission line in a reverberation chamber［J］.IEEE Transactions on Electromagnetic Compatibility，2011，53（2）：308-317.

［12］VANHEE F，PISSOORT D，CATRYSSE J，et al.Efficient reciprocity-based algorithm to predict worst case induced disturbances on multiconductor transmission lines due to incoming plane waves［J］.IEEE Transactions on Electromagnetic Compatibility，2013，55（1）：208-216.

［13］KILDAL P S.Equivalent circuits of receive antennas in signal processing arrays［J］.Microwave and Optical Technology Letters，1999，21（4）：244-246.

［14］STOLLE R.Electromagnetic coupling of twisted pair cables［J］.IEEE Journal on Selected Areas in Communications，2002，20（5）：883-892.

［15］賈 銳，王慶國，程二威，等.混響室“全向輻照”電磁環(huán)境場線耦合規(guī)律分析［J］.電波科學(xué)學(xué)報(bào)，2014，29（2）：385-390.JIA Rui，WANG Qingguo，CHENG Erwei，et al.Analysis of field-to-wire coupling in full radiation environment in reverberation chamber［J］.Chinese Journal of Radio Science，2014，29（2）：385-390.（in Chinese）

［16］王利萍，周東方，彭 強(qiáng)，等.基于拓?fù)渚W(wǎng)絡(luò)的屏蔽腔體內(nèi)置微帶線響應(yīng)分析［J］.電波科學(xué)學(xué)報(bào)，2014，29（1）：321-327.WANG Liping，ZHOU Dongfang，PENG Qiang，et al.Response analysis of microstrip lines in shielding cavity based on topological network［J］.Chinese Journal of Radio Science，2014，29（1）：321-327.（in Chinese）

［17］劉順坤，聶 鑫，陳向躍.電磁脈沖對(duì)電纜耦合問題的理論研究［J］.電波科學(xué)學(xué)報(bào)，2010，25（2）：348-352.LIU Shunkun，NIE Xin，CHEN Xiangyue.Numerical study on cable coupling effects excited by electromagnetic pulse［J］.Chinese Journal of Radio Science，2010，25（2）：348-352.（in Chinese）