萬 博，白嶸峰，張 龍，楊 洋，崔耀中

(1. 總裝西安軍事代表局駐華縣地區(qū)軍事代表室，西安 710071；2. 西安軍事代表局政治部，西安 710071；3. 軍械工程學(xué)院 靜電與電磁防護(hù)研究所，石家莊 050003)

0 引 言

混響室是一種能在屏蔽腔體內(nèi)產(chǎn)生空間統(tǒng)計(jì)均勻、各向同性、隨機(jī)極化電磁環(huán)境的模擬設(shè)備。通過對混響室進(jìn)行數(shù)值仿真計(jì)算，對于了解其內(nèi)部電磁場的分布，減少相關(guān)的調(diào)試工作具有重要的作用。建立一個能夠真實(shí)反映實(shí)際混響室場性能的混響室數(shù)值仿真模型，將會對實(shí)際混響室內(nèi)的電磁兼容測試具有重大指導(dǎo)意義。雖然國內(nèi)外開展了一些建立混響室仿真模型的研究工作，但在模型建立方面大多是采用二維結(jié)構(gòu)建模、將混響室面壁設(shè)為理想導(dǎo)體，以及激勵源簡化為偶極子天線模型［1，2］，以致無法真實(shí)準(zhǔn)確反映實(shí)際混響室的場性能特點(diǎn);而在模型驗(yàn)證方面或者僅進(jìn)行電場強(qiáng)度大小的比較［3］，或者僅進(jìn)行混響室統(tǒng)計(jì)規(guī)律上的驗(yàn)證［4］，以致無法對混響室性能進(jìn)行全面驗(yàn)證。本文采用基于矩量法的電磁仿真軟件FEKO 建立混響室仿真模型，按照混響室實(shí)際材料構(gòu)成設(shè)置面壁材料參數(shù)，以實(shí)際測試中的對數(shù)周期天線為激勵源，并從電場強(qiáng)度大及場均勻性兩方面對混響室模型進(jìn)行驗(yàn)證。

1 混響室仿真模型的建立

以實(shí)驗(yàn)室現(xiàn)有的混響室為依據(jù)，建立1 ∶1 的混響室仿真模型，其內(nèi)部結(jié)構(gòu)分別如圖1(a)、1(b)所示?；祉懯夷Ｐ桶ㄆ帘吻惑w、攪拌器、發(fā)射天線三部分。

圖1 混響室內(nèi)部結(jié)構(gòu)

1.1 屏蔽腔體建模

建立10.5 m(長)×8 m(寬)×4.3 m(高)的矩形屏蔽腔體。根據(jù)諧振腔理論，其諧振頻率fmnp為

式中，L、W、H 分別為混響室的長、寬、高，單位為m;m，n，p 為0 或正整數(shù);c 為光速。由式(1)可知，此屏蔽腔體最低諧振頻率f110為23.6 MHz。為滿足混響室場均勻性要求，最低工作頻率一般取最低諧振頻率的3 倍，即其工作頻率應(yīng)大于70.7 MHz，本文選擇工作頻率80 MHz、150 MHz 進(jìn)行研究。腔體材料選擇鍍鋅鋼板，其電導(dǎo)率與相對磁導(dǎo)率之比σ/μr為4.4 ×106S·m－1。

1.2 攪拌器建模

攪拌器是混響室內(nèi)的重要組成部分。通過攪拌器的旋轉(zhuǎn)，混響室邊界條件不斷發(fā)生改變，腔體內(nèi)駐波分布不斷變化，最終形成統(tǒng)計(jì)均勻的電磁環(huán)境。

本文建立了兩個結(jié)構(gòu)不同的攪拌器，如圖2 所示。其中橫向放置的攪拌器為Z 型結(jié)構(gòu)，由四個葉片組成，葉片夾角φ 為128.3°;豎直放置的攪拌器由六個V 型槳葉組成，槳葉夾角θ 為90°，為增大攪拌效率，槳葉采用V 型切口設(shè)計(jì)。

1.3 發(fā)射天線建模

設(shè)計(jì)工作頻率60 MHz ～1.9 GHz，對數(shù)周期偶極子陣列作為發(fā)射天線，如圖3 所示。其中最短振子長0.08 m，最長振子長2. 52 m，增長因子τ 為0.88，天線全長1.6 m。根據(jù)IEC61000 －4 －21 標(biāo)準(zhǔn)［5］，發(fā)射天線位置應(yīng)對準(zhǔn)墻角，天線輸入功率100 W。

2 仿真計(jì)算

混響室場強(qiáng)直接反映混響室產(chǎn)生的電磁環(huán)境場強(qiáng)情況，而場均勻性則反映了混響室的統(tǒng)計(jì)規(guī)律特性，在驗(yàn)證混響室仿真模型時這兩個參數(shù)不可偏廢。采用IEC61000 －4 －21 標(biāo)準(zhǔn)給出的方法來判定混響室工作區(qū)域的場均勻性，即在攪拌器步進(jìn)一周后，求解混響室工作區(qū)域8 個頂點(diǎn)位置處對輸入功率歸一化的最大電場值的標(biāo)準(zhǔn)偏差，標(biāo)準(zhǔn)偏差越小，說明工作區(qū)域的場均勻性越好。工作區(qū)域距離腔室6 個面、天線及攪拌器的最短距離應(yīng)大于λ/4 或者1 m，λ 為最低工作頻率波長。在此選擇矩形工作區(qū)域坐標(biāo)區(qū)間x∈［2.9，7.0］;y∈［4.3，9.5］;z∈［1.0，3.3］。

IEC61000 －4 －21 標(biāo)準(zhǔn)介紹了兩種攪拌方式，即連續(xù)攪拌和步進(jìn)攪拌，本文選擇步進(jìn)攪拌。通過編寫Matlab 程序控制FEKO，實(shí)現(xiàn)仿真模型中攪拌器的自動步進(jìn)，共步進(jìn)50 個位置。在攪拌器步進(jìn)一周內(nèi)，記錄測試區(qū)域8 個頂點(diǎn)位置在正交軸上的最大電場強(qiáng)度Emaxik(i 代表正交軸x、y、z 方向，k 代表頂點(diǎn)位置1 到8)及平均輸入功率Pinput，求出各頂點(diǎn)處Emaxik對Pinput的歸一化最大電場值，為

各正交軸歸一化最大電場值標(biāo)準(zhǔn)偏差

8 個頂點(diǎn)，每個頂點(diǎn)3 個正交軸的歸一化電場最大值Eik標(biāo)準(zhǔn)偏差

標(biāo)準(zhǔn)偏差σ 相對于平均值可以表示為單位為dB 的量

隨著頻率的增加，混響室內(nèi)模式數(shù)增多，各個模式相互疊加，將使得混響室場均勻性越來越好。仿真得到的混響室工作區(qū)域標(biāo)準(zhǔn)偏差如圖4 所示反映了這一點(diǎn)。由圖4 可知，隨著頻率的增大，雖然混響室的標(biāo)準(zhǔn)偏差σr(r =x，y，z，24)會出現(xiàn)后者大于前者的情況，但是整體的趨勢是越來越小，尤其體現(xiàn)在表征整體均勻性的σ24上，從而驗(yàn)證了混響室的場均勻性隨著頻率的增大而越來越好。

圖4 混響室標(biāo)準(zhǔn)偏差

3 混響室模型的實(shí)驗(yàn)驗(yàn)證

本文從場強(qiáng)及反映混響室統(tǒng)計(jì)規(guī)律的場均勻性兩個指標(biāo)來驗(yàn)證混響室仿真模型。為驗(yàn)證優(yōu)化仿真的正確性，需要對實(shí)際混響室進(jìn)行測試，分別對80 MHz、150 MHz 這兩個頻點(diǎn)進(jìn)行了測試。采用Lab Windows/CVI 軟件平臺開發(fā)了混響室場均勻性測試程序，通過GPIB 卡對信號源、功率計(jì)、攪拌器以及場強(qiáng)計(jì)進(jìn)行控制，自動采集數(shù)據(jù)并進(jìn)行處理，最后輸出電場強(qiáng)度及電場標(biāo)準(zhǔn)偏差數(shù)據(jù)文件，測試流程如圖5 所示。

3.1 電場強(qiáng)度驗(yàn)證

圖6、圖7 分別為80 MHz、150 MHz 時，仿真得到的工作區(qū)域各頂點(diǎn)歸一化最大電場強(qiáng)度Eik與實(shí)際測試的比較。從圖6、圖7 中可知，仿真計(jì)算得到的各頂點(diǎn)Ezk歸一化最大電場強(qiáng)度趨勢與實(shí)測結(jié)果相似性較好，Exk、Eyk變化趨勢相似性稍差的原因是Exk、Eyk的方向具有沿軸向旋轉(zhuǎn)的任意性，與仿真值不具有一一對應(yīng)性。另外，仿真計(jì)算得到的歸一化場強(qiáng)值Eik高于實(shí)際測試的結(jié)果，出現(xiàn)此現(xiàn)象的原因是由于:仿真計(jì)算忽略了實(shí)際混響室中通風(fēng)口、電纜口的電磁損耗;仿真計(jì)算未考慮混響室實(shí)際測試中支架、場強(qiáng)計(jì)等的加載效應(yīng)及電纜的損耗。

圖5 測試流程圖

圖7 頻率為150 MHz 時Eik仿真結(jié)果與實(shí)測結(jié)果的比較

3.2 工作區(qū)域場均勻性驗(yàn)證

仿真得到的工作區(qū)域標(biāo)準(zhǔn)偏差與實(shí)際測試結(jié)果的比較見表2。根據(jù)表2，仿真與實(shí)測標(biāo)準(zhǔn)偏差的相對誤差值見表3，其中相對誤差ε 計(jì)算表達(dá)式為

從表3 可知，仿真計(jì)算得到的標(biāo)準(zhǔn)偏差值與實(shí)測結(jié)果的相對誤差平均值保持在12%左右，說明仿真模型能真實(shí)反映混響室的場分布情況。

表2 仿真得到的標(biāo)準(zhǔn)偏差與實(shí)測結(jié)果比較

4 結(jié) 語

本文依據(jù)實(shí)驗(yàn)室現(xiàn)有的混響室建立了1 ∶1 的混響室仿真模型，在頻率80 MHz、150 MHz 時對混響室進(jìn)行仿真計(jì)算，從場強(qiáng)大小和場均勻性兩個方面對混響室模型進(jìn)行了驗(yàn)證。結(jié)果表明:仿真計(jì)算得到的各頂點(diǎn)歸一化最大電場強(qiáng)度趨勢與實(shí)測結(jié)果具有較好的一致性，仿真計(jì)算得到的標(biāo)準(zhǔn)偏差值與實(shí)測結(jié)果的相對誤差平均值保持在12%左右，說明仿真模型能真實(shí)反映混響室的場分布情況，從而為混響室的設(shè)計(jì)測試工作奠定一定的基礎(chǔ)。

［1］DAVID A HILL. Electronic Mode Stirring for Reverberation Chamber［J］. IEEE transactions on electromagnetic compatibility，1994，36(4):295-299.

［2］ JANET CLEGG，ANDREW C MARVIN，JOHN F DAWSON，et al.Optimization of Stirrer Designs in a Reverberation Chamber［J］. IEEE Transactions on Electromagnetic Compatibility，2005，40(7):824-831.

［3］CHRISTIAN BRUNS，RUEDIGER VAHLDIECK.A Closer Look at Reverberation Chambers-3-D Simulation and Experimental Verification［J］. IEEE transactions on electromagnetic compatibility，2005，47(5):612-626.

［4］ALYSE COATES，HUGH G SASSE，DAWN E COLEBY，et al.Validation of a Three-Dimensional Transmission Line Matrix(TLM)Model Implementation of a Mode-Stirred Reverberation Chamber［J］. IEEE Transactions on Electromagnetic Compatibility，2007，49(4):734-744.

［5］IEC61000-4-21，Testing and measurement techniques Reverberation Chamber Test Methods［S］.2003.