李 爽 王建國(guó)，2 謝海燕 陸希成

（1.西北核技術(shù)研究所，陜西 西安 710024；2.西安交通大學(xué)電子與信息工程學(xué)院，陜西 西安710049）

引 言

隨著高功率微波（HPM）技術(shù)的快速發(fā)展［1］，現(xiàn)代武器系統(tǒng)所處的電磁環(huán)境更趨復(fù)雜，這給武器裝備的安全性構(gòu)成了嚴(yán)重威脅。因此，有必要對(duì)各類武器系統(tǒng)進(jìn)行HPM輻照效應(yīng)研究。通常對(duì)武器系統(tǒng)進(jìn)行HPM效應(yīng)研究，可以利用外場(chǎng)或者微波暗室來(lái)進(jìn)行，相關(guān)研究也都取得了豐富的成果。但為了得到統(tǒng)計(jì)規(guī)律，這種試驗(yàn)方法需要的統(tǒng)計(jì)樣本較多，試驗(yàn)成本較為昂貴，實(shí)施也比較復(fù)雜。同時(shí)，整系統(tǒng)的效應(yīng)規(guī)律與其內(nèi)部子系統(tǒng)的效應(yīng)特征也密切相關(guān)。對(duì)于復(fù)雜腔體中的子系統(tǒng)設(shè)備，所處的電磁環(huán)境更為復(fù)雜。外部HPM可以通過(guò)多種耦合途徑（天線、腔體上的孔縫等）進(jìn)入腔體內(nèi)部，對(duì)腔體內(nèi)的子系統(tǒng)設(shè)備產(chǎn)生影響，從而嚴(yán)重影響整系統(tǒng)的工作效能。因此，開(kāi)展對(duì)復(fù)雜腔體內(nèi)子系統(tǒng)的HPM效應(yīng)研究也是十分必要的。而傳統(tǒng)試驗(yàn)場(chǎng)地也難以模擬這種復(fù)雜的電磁環(huán)境。這給研究復(fù)雜腔體內(nèi)子系統(tǒng)的HPM效應(yīng)帶來(lái)了很大的困難。

混響室作為一種新型電磁效應(yīng)試驗(yàn)設(shè)備，能夠在腔體內(nèi)產(chǎn)生統(tǒng)計(jì)意義上的各向同性均勻場(chǎng)，可以很好地模擬復(fù)雜電磁環(huán)境，因而對(duì)研究復(fù)雜腔體內(nèi)部子系統(tǒng)的電磁效應(yīng)具有重要意義［2－5］。美國(guó)的Crawford等對(duì)混響室的結(jié)構(gòu)設(shè)計(jì)、性能評(píng)價(jià)方面進(jìn)行了一系列的試驗(yàn)和研究［6］；1992年英國(guó)的Huang教授首次提出了源攪拌混響室的概念［7］。國(guó)內(nèi)方面，北京郵電大學(xué)的高攸綱、沈遠(yuǎn)茂［8－9］以及買望［10］等人均在混響室設(shè)計(jì)及均勻性分析方面做了大量工作。然而，目前的混響室研究主要是基于連續(xù)波激勵(lì)，很少有涉及脈沖波激勵(lì)下的混響室特性研究。本文將從理論分析和數(shù)值模擬兩個(gè)方面，研究HPM激勵(lì)下，不同的腔體形狀對(duì)腔內(nèi)電磁場(chǎng)分布的影響，為該型混響室的設(shè)計(jì)提供技術(shù)參考。

1.理論分析

1.1 混響室的攪拌原理

實(shí)際應(yīng)用中的混響室，一般是通過(guò)其內(nèi)部的源天線激勵(lì)，構(gòu)造電磁場(chǎng)環(huán)境，這相當(dāng)于一個(gè)有源諧振腔。研究混響室內(nèi)的攪拌原理，可以借助有源諧振腔內(nèi)部的電磁場(chǎng)理論進(jìn)行分析。下面以簡(jiǎn)單的二維有源矩形諧振腔為例進(jìn)行研究，如圖1所示。

圖1 二維有源諧振腔

腔體放置在xOy平面上，邊長(zhǎng)分別為a和b，腔體材料為理想導(dǎo)體，其他區(qū)域均為自由空間。激勵(lì)源是沿z軸放置的電振子（I0），位置在（x0，y0）.腔體內(nèi)部的非零電場(chǎng)分量Ez可以表示為［8］

式中：k2＝ω2με.二維模型中的諧振頻率為

式中：m、n分別表示x、y方向的半波數(shù)；a、b分別為矩形腔體的長(zhǎng)、寬。

從式（1）可以看出：在空間介質(zhì)不變的情況下，腔體內(nèi)的場(chǎng)分布取決于以下因素：①腔體的形狀、尺寸；②激勵(lì)源的位置、幅度等。

傳統(tǒng)的混響室是根據(jù)因素①，在腔體內(nèi)設(shè)置機(jī)械攪拌器來(lái)實(shí)現(xiàn)攪拌功能的。當(dāng)攪拌器進(jìn)行旋轉(zhuǎn)時(shí)，即可以不斷地改變腔內(nèi)的邊界條件，進(jìn)而改變腔體內(nèi)模式的諧振狀態(tài)，從而實(shí)現(xiàn)混響室所要求的統(tǒng)計(jì)均勻場(chǎng)。

源攪拌混響室是根據(jù)因素②來(lái)工作的［11－12］。通過(guò)對(duì)激勵(lì)源的擾動(dòng)來(lái)影響有效模式組合中各模式所占的權(quán)重，進(jìn)而對(duì)腔體內(nèi)的場(chǎng)分布產(chǎn)生攪動(dòng)效果。

根據(jù)微波混沌理論，復(fù)雜腔體不同于規(guī)則腔體，其內(nèi)部電磁波的傳播一般具有復(fù)雜性和非周期性，除了一些特殊方向的射線能夠產(chǎn)生規(guī)則軌道之外，在大部分情況下它們都將產(chǎn)生復(fù)雜的不規(guī)則軌道。而且空腔中規(guī)則的短周期軌道具有強(qiáng)烈的不穩(wěn)定性，會(huì)對(duì)本征函數(shù)造成“疤痕”現(xiàn)象［13］。腔體內(nèi)電磁波分布的復(fù)雜度對(duì)電磁場(chǎng)分布的統(tǒng)計(jì)均勻性有重要影響。而在腔體尺寸一定的情況下（參數(shù)a、b保持大致不變），改變腔體內(nèi)部的細(xì)微結(jié)構(gòu)，可以改變腔體內(nèi)部場(chǎng)的邊界條件，進(jìn)而影響場(chǎng)分布的狀態(tài)。因此，腔體形狀對(duì)腔內(nèi)場(chǎng)分布狀態(tài)具有重要的影響。

重點(diǎn)研究HPM激勵(lì)下的混響室，由于脈沖持續(xù)時(shí)間短，遠(yuǎn)小于機(jī)械攪拌中機(jī)械臂運(yùn)動(dòng)的特征時(shí)間，而且腔內(nèi)場(chǎng)的瞬態(tài)特性強(qiáng)，難以在腔體內(nèi)形成類似連續(xù)波那樣的穩(wěn)定狀態(tài)。因此，機(jī)械攪拌的方法不適于用HPM效應(yīng)分析的混響室。這里考慮采用源攪拌方法，通過(guò)多天線激勵(lì)的方式來(lái)實(shí)現(xiàn)對(duì)腔體內(nèi)電場(chǎng)的攪拌，以此為基礎(chǔ)，在相同攪拌條件下研究腔體形狀對(duì)場(chǎng)分布的影響。

1.2 場(chǎng)均勻性分析方法

由于混響室中的場(chǎng)均勻性是統(tǒng)計(jì)意義上的概念，需要借助統(tǒng)計(jì)計(jì)算方法來(lái)驗(yàn)證場(chǎng)的均勻性。根據(jù)IEC 61000－4－21標(biāo)準(zhǔn)［14］，對(duì)于穩(wěn)態(tài)場(chǎng)的混響室校準(zhǔn)，是在測(cè)試區(qū)域的頂點(diǎn)處放置8個(gè)電場(chǎng)探頭，攪拌器以步進(jìn)方式旋轉(zhuǎn)一圈，每個(gè)步進(jìn)位置為一個(gè)采樣點(diǎn)。在每個(gè)采樣點(diǎn)上，用8個(gè)探頭測(cè)量各個(gè)位置三個(gè)正交方向的電場(chǎng)強(qiáng)度。用所有采樣點(diǎn)場(chǎng)值的標(biāo)準(zhǔn)差來(lái)反映場(chǎng)的空間均勻性。但是在脈沖波激勵(lì)下，有耗腔體內(nèi)的場(chǎng)隨時(shí)間劇烈變化，難以形成穩(wěn)定狀態(tài)。本文重點(diǎn)關(guān)注脈沖波對(duì)器件的電磁瞬態(tài)效應(yīng)，所以，不能像穩(wěn)態(tài)場(chǎng)那樣只考慮場(chǎng)的空間均勻性，還要關(guān)注場(chǎng)的時(shí)域特性。在本文中，主要參考IEC 61000－4－21，觀測(cè)一段時(shí)間內(nèi)整個(gè)測(cè)試區(qū)域內(nèi)的場(chǎng)值，分析電場(chǎng)隨時(shí)間的變化特征，并選取若干平面來(lái)分析場(chǎng)分布的空間均勻性。場(chǎng)值分布的標(biāo)準(zhǔn)差越小則說(shuō)明場(chǎng)值分布越均勻，因而也就越有利于在腔體內(nèi)對(duì)器件進(jìn)行相關(guān)的瞬態(tài)效應(yīng)實(shí)驗(yàn)。

為觀測(cè)腔體內(nèi)瞬態(tài)電場(chǎng)的時(shí)域特性，選取一定長(zhǎng)度的時(shí)間窗口作為研究單元。在每個(gè)時(shí)間窗內(nèi)，觀測(cè)整個(gè)測(cè)量空間內(nèi)每個(gè)場(chǎng)點(diǎn)的電場(chǎng)最大值，綜合所有的場(chǎng)點(diǎn)，分析整個(gè)空間內(nèi)的最大值及場(chǎng)值分布標(biāo)準(zhǔn)差，可以得到該時(shí)間窗內(nèi)腔體中的電場(chǎng)分布特征。綜合若干個(gè)時(shí)間窗，還可以得到腔體內(nèi)電場(chǎng)的時(shí)域分布特征。由此，可以得到不同條件下電場(chǎng)的空間分布特征，總結(jié)出改善場(chǎng)分布均勻性的規(guī)律。具體的統(tǒng)計(jì)計(jì)算過(guò)程如下：

腔體內(nèi)某測(cè)試點(diǎn)（xi，yi，zi）在時(shí)刻點(diǎn)t1的電場(chǎng)表示為

該測(cè)試點(diǎn)在某個(gè)時(shí)間窗Tj（時(shí)刻t1到tp）內(nèi)的電場(chǎng)最大值為

單方向上的平均值 〈Eξ〉（ξ＝x，y，z）和3個(gè)方向的聯(lián)合平均值〈Exyz〉可以分別定義為

式中，N為單方向上觀測(cè)的場(chǎng)點(diǎn)個(gè)數(shù)。

單方向上的標(biāo)準(zhǔn)差σξ和3個(gè)方向的聯(lián)合標(biāo)準(zhǔn)差σxyz可以分別定義為

1.3 各向同性

混響室內(nèi)的電磁場(chǎng)分布具有各向同性的特點(diǎn)，腔內(nèi)各個(gè)方向上的能流密度均衡，便于對(duì)受試設(shè)備進(jìn)行全向測(cè)試［15］。腔體內(nèi)電場(chǎng)的各向同性研究，可以通過(guò)比較各頻點(diǎn)下電場(chǎng)分量的平方值得到。為了避免激勵(lì)源參數(shù)（激勵(lì)幅度等）對(duì)該分布的影響，主要研究各方向上的能流密度占總場(chǎng)能流密度的比值理想情況下，三個(gè)比值均應(yīng)為1/3.

2.數(shù)值計(jì)算

2.1 數(shù)值計(jì)算模型

主要計(jì)算并對(duì)比矩形腔體、復(fù)雜腔［16］（在腔體內(nèi)壁上引入若干半球形金屬凸起結(jié)構(gòu)，凸起結(jié)構(gòu)的半徑約為1倍波長(zhǎng)）、凹形腔（凹去部分尺寸為2倍波長(zhǎng)）等形狀腔體內(nèi)的場(chǎng)分布形式。其中，矩形腔代表的是規(guī)則腔體，復(fù)雜腔是根據(jù)聲學(xué)混響室中的返波體引申得到的形狀，可以增加腔體內(nèi)電磁波的反射方向，而凹形腔是電磁統(tǒng)計(jì)學(xué)中典型的Sinai腔模型。研究的腔體尺寸均約為6倍波長(zhǎng)，具體形狀如圖2所示。在腔內(nèi)均采用載頻3GHz，脈寬10ns的脈沖波進(jìn)行激勵(lì)。

對(duì)于每種模型，首先采用時(shí)域有限差分（FDTD）方法分析整個(gè)腔體內(nèi)電場(chǎng)的時(shí)域分布特征［17］，得到場(chǎng)分布結(jié)果；在場(chǎng)值最大的時(shí)間窗內(nèi)，進(jìn)一步分析腔內(nèi)不同平面上的電場(chǎng)分布。選取兩個(gè)z平面進(jìn)行研究，其中z1面靠近腔壁，z2面位于腔體中央。整個(gè)時(shí)域研究長(zhǎng)度選取為10個(gè)時(shí)間窗Ti（i＝0，1，2，……，9），其中T0代表天線開(kāi)始進(jìn)行激勵(lì)的時(shí)間窗，時(shí)間窗的長(zhǎng)度選取與脈寬相同。還分析了腔內(nèi)電場(chǎng)的各向同性分布。

2.2 數(shù)值計(jì)算結(jié)果

2.2.1 整個(gè)空間內(nèi)電場(chǎng)的時(shí)域分布特點(diǎn)

根據(jù)圖3（a）可以看出：復(fù)雜腔、矩形腔和凹形腔在T1、T2時(shí)間窗內(nèi)，腔內(nèi)電場(chǎng)值達(dá)到最大，隨后場(chǎng)值逐漸衰減。這兩個(gè)時(shí)間段是進(jìn)行瞬態(tài)效應(yīng)研究的重點(diǎn)關(guān)注區(qū)域。根據(jù)圖3（b），復(fù)雜腔和凹形腔在T1和T2內(nèi)的場(chǎng)均勻性都比矩形腔好。這是由于加入了復(fù)雜結(jié)構(gòu)，增加了電磁波的反射方向，而各場(chǎng)點(diǎn)上的電場(chǎng)是直射、反射電磁波的疊加，所以，各場(chǎng)點(diǎn)的場(chǎng)值更接近。該結(jié)果表明：在矩形腔結(jié)構(gòu)的基礎(chǔ)上引入若干復(fù)雜結(jié)構(gòu)，可以有效地改善腔體內(nèi)的電場(chǎng)分布，提高場(chǎng)分布的均勻性。不過(guò)，凹形腔的特殊結(jié)構(gòu)會(huì)造成腔內(nèi)測(cè)量空間的減小，影響實(shí)際應(yīng)用范圍。

2.2.2 平面上電磁場(chǎng)的分布特點(diǎn)

從電場(chǎng)的時(shí)域分布圖可以看出，在腔體內(nèi)激勵(lì)剛完成的時(shí)刻，即T1、T2時(shí)間窗內(nèi)，腔體內(nèi)的電場(chǎng)值最大。因此，這兩個(gè)時(shí)間窗是我們開(kāi)展電磁效應(yīng)研究的重點(diǎn)，選取T1下z1和z2兩個(gè)面來(lái)研究腔體內(nèi)部不同位置處的場(chǎng)分布特征，具體結(jié)果見(jiàn)圖4所示。

從圖4可以看出：復(fù)雜腔中的場(chǎng)分布更加均勻，整個(gè)平面上的電場(chǎng)值相差不大。而且，平面z2比平面z1的場(chǎng)分布要更加均勻。而其他兩種模型中，場(chǎng)值的空間分布變化較大，統(tǒng)計(jì)均勻性不如復(fù)雜腔。這說(shuō)明在復(fù)雜腔內(nèi)電場(chǎng)的分布更加均勻，采用該種腔形可以很好地改善腔內(nèi)場(chǎng)分布的均勻性。由于腔壁表面電流的影響，腔內(nèi)中央處的場(chǎng)均勻性也比靠近腔壁處的好。所以，采用類似復(fù)雜腔的結(jié)構(gòu)，在矩形腔的基礎(chǔ)上引入若干復(fù)雜結(jié)構(gòu)，可以增加腔內(nèi)電磁波的反射，更有利于在腔內(nèi)構(gòu)造均勻的電磁環(huán)境。

2.2.3 各向同性

混響室的特征之一就是具有良好的各向同性，即各方向上的能流密度均衡，能夠?qū)Ρ粶y(cè)設(shè)備進(jìn)行全向均勻測(cè)試。從圖5可以明顯看出，復(fù)雜腔和凹形腔比矩形腔內(nèi)的能流分布更均衡。而且頻率越高，各方向上的能量分布就越均衡。因此，利用混響室技術(shù)對(duì)器件進(jìn)行輻照試驗(yàn)時(shí)，可以實(shí)現(xiàn)對(duì)器件的全向輻照，不需要每輻照一個(gè)方向就移動(dòng)一次天線。

圖5 不同模型中的各向同性結(jié)果

3.結(jié) 論

通過(guò)采用源攪拌的方法，研究了在HPM激勵(lì)下腔體形狀對(duì)腔體內(nèi)電場(chǎng)分布的影響效果。首先通過(guò)諧振腔理論和微波混沌理論的分析，證明腔體形狀對(duì)腔體內(nèi)的場(chǎng)分布有重要影響。其次，以多天線激勵(lì)為基礎(chǔ)，通過(guò)數(shù)值模擬計(jì)算驗(yàn)證腔體形狀對(duì)腔內(nèi)場(chǎng)分布的影響。結(jié)果表明：在矩形腔的基礎(chǔ)上引入一定的復(fù)雜結(jié)構(gòu)，可以增強(qiáng)腔內(nèi)電磁波的反射效果，增加腔內(nèi)電磁場(chǎng)的復(fù)雜度，有效地提高場(chǎng)分布的均勻性及能量的均衡性，對(duì)構(gòu)造HPM下的混響室有重要意義。但是，也要注意引入復(fù)雜結(jié)構(gòu)后對(duì)腔體帶來(lái)的其他影響，如本文中的凹形腔體會(huì)影響腔內(nèi)的實(shí)際測(cè)量空間。同時(shí)，由于目前混響室實(shí)驗(yàn)室還處于搭建階段，相關(guān)驗(yàn)證實(shí)驗(yàn)尚未能進(jìn)行。開(kāi)展相關(guān)實(shí)驗(yàn)，對(duì)理論結(jié)果進(jìn)行驗(yàn)證也是下一步的工作重點(diǎn)。

［1］BENFORD J，SWEGLE J，SCHAMILOGLU E.High Power Microwaves［M］.2nd ed.New York：CRC Press，2007：35－37.

［2］HILL D A.Electronic mode stirring for reverberation chamber［J］.IEEE Transactions on Electromagnetic Compatibility，1994，36（4）：294－299.

［3］張林昌.混響室及其進(jìn)展［J］.安全與電磁兼容，2001，（1）：2－8.ZHANG Linchang.The development of reverberation studio［J］.Safety and EMC，2001，（1）：2－8.（in Chinese）

［4］GRADONI G，ARNAUT L R.Generalized extreme－value distributions of power near a boundary inside electromagnetic reverberation chamber ［J］.IEEE Transactions on Electromagnetic Compatibility，2010，52（3）：506－515.

［5］BACKSTROM M G，LOVSTRAND K G.Susceptibility of electronic systems to high－power microwaves：summary of test experience［J］.IEEE Transactions on Electromagnetic Compatibility，2004，46（3）：396－403.

［6］CRAWFORD M L，KOEPKE G H.Design，Evaluation，and Use of a Reverberation Chamber for Performing Electromagnetic Susceptibility/Vulnerability Measurements［S］.Washington：US Department of commerence，National Bureau of standards 1986.

［7］HUANG Y，EDWARDS D J.A novel reverberating chamber：the source－stirred chamber［C］//8th International Conference on Electromagnetic Compatibility.Oxford，UK，1992：120－124.

［8］沈遠(yuǎn)茂，高攸剛.電磁兼容測(cè)試中的源攪拌混響室和電磁干擾接收機(jī)的相關(guān)研究［D］.北京：北京郵電大學(xué)，2006.SHEN Yuanmao，GAO Yougang.Research about Source－Stirred Reverberation Chamber and Electromagnetic Interference Measuring Receiver in the Field of Electromagnetic Compatibility Tests［D］.Beijing：Beijing University of posts and telecommunications，2006.（in Chinese）

［9］沈遠(yuǎn)茂，石 丹，高攸綱，等.利用多天線源攪拌改善混響室場(chǎng)均勻性的分析［J］.電波科學(xué)學(xué)報(bào)，2009，24（4）：682－686.SHEN Yuanmao，SHI Dan，GAO Yougang，et al.Improvement in field uniformity introduced by multiple－antenna in source－stirred reverberation chamber［J］.Chinese Journal of Radio Science，2009，24（4）：682－686.（in Chinese）

［10］買 望，高攸綱.混響室場(chǎng)均勻性的仿真與分析［D］.北京：北京郵電大學(xué)，2006.MAI Wang，GAO Yougang.Simulation and Analyzation of Uniformity inside Reverberation Chamber［D］.Beijing：Beijing University of posts and telecommunications，2006.（in Chinese）

［11］CERRI G，MARIANI V，PENNSI S，et al.Source Stirring Mode for Reverberation Chambers［J］.IEEE Transactions on Electromagnetic Compatibility，2005，47（4）：815－823.

［12］梁雙港，許家棟，劉易勇，等.源攪拌混響室的仿真分析與實(shí)驗(yàn)研究［J］.電波科學(xué)學(xué)報(bào)，2010，25（6）：1058－1063.LIANG Shuanggang，XU Jiadong，LIU Yiyong，et al.Simulation analysis and experimental investigation of SSRC［J］.Chinese Journal of Radio Science，2010，25（6）：1058－1063.（in Chinese）

［13］陸希成，王建國(guó).復(fù)雜腔體電磁耦合效應(yīng)的統(tǒng)計(jì)方法研究［D］.西安：西北核技術(shù)研究所，2010.LU Xicheng，WANG Jianguo.Statistical Method Study in Electromagnetic Coupling of Complex Cavities［D］.Xi’an：Norethwest Institute of Nuclear Technology，2010.（in Chinese）

［14］International Electrotechnical Commission.Electromagnetic compatibility （EMC）－Part 4－21：Testing and measurement techniques－reverberation chamber test methods［R］.Geneva：IEC，2003.

［15］HERLEM Y，SCHAFFAR A，GINESTE P N，et al.Oversized cavity limits assessment by numerical simulation［R］.Florence：ESA Workshop on Aerospace EMC，2009.

［16］PETRISCH M，SCHWAB A J.Investigation of the field uniformity of a mode－stirred chamber using diffusers based on acoustic theory［J］.IEEE Transactions on Electromagnetic Compatibility，1999，41（4）：446－451.

［17］黃 華，牛中奇，白 冰.時(shí)域有限差分法分析混響室中場(chǎng)的均勻性［J］.電波科學(xué)學(xué)報(bào)，2011，26（1）：124－132.HUANG Hua，NIU Zhong－qi，BAI Bing.Using FDTD to analyze the field uniformity in reverberation chamber［J］.Chinese Journal of Radio Science，2011，26（1）：124－132.（in Chinese）