MLD結構下電磁兼容混響室內場均勻性分析

劉重陽，黃 華，張 行

（1.海裝航訂部，北京100841；2.船舶重工集團公司723所，揚州225001）

0 引 言

混響室作為當前主要的電磁兼容測試設備，越來越廣泛地運用于輻射發射和輻射免疫性測試，以及封閉和屏蔽電纜的屏蔽效力等的測試。然而傳統混響室的每次測試都是要伴隨攪拌器足夠多的步進數才能獲得統計均勻的場，這樣便使測試時間大大增加。如果采用散射體來代替攪拌器的作用，那么1次測試時間就僅是攪拌器步進1次的間歇時間，這樣就在相當大的程度上縮短了測試時間，使得混響室更具有經濟實用價值。因此，為了獲得均勻的場而合理安排散射體在混響室中的位置就顯得非常重要。

一般而言，幾乎所有的混響室都是長方體的，利用模攪拌和模調諧來改變混響室內電磁場的邊界條件，從而產生統計均勻分布的場［1］。國內外研究者在攪拌器的基礎上提出了改變混響室的形狀［2］、用散射體［3］以及散射體和攪 拌器相結合等方法［4－5］來優化混響室的性能。其中，混響室的形狀可以改變成三角形和不對稱的長方形；散射體類型最常用的有最大長度散射體（MLD）、施羅德二次剩余散射體（QRD）和原根剩余散射體（PRD）。

本文將從MLD的散射原理出發，具體給出它的設計步驟，最后分析了混響室內MLD不同位置下的場均勻性。

1 MLD散射原理與設計

1.1 MLD散射原理

散射體與相位反射柵的原理很相似，如圖1所示。

圖1 MLD的結構和工作原理示意圖

不同的阱深dn（n＝1，2，…）導致反射波與入射波的相位差為［6］：

式中：λ為電磁波在自由空間的波長。

實際上，散射體由在同一基準面上設置1組高度不同的長方體“脊”構成，習慣上，人們將長方體“脊”稱為“阱”。散射體能夠使混響室內某一區域的場成為統計均勻的原理如圖1所示。具體說明如下：混響室內被關注區域（即企圖使場為統計均勻的區域，一般在混響室的中部）任意一個點上的場，都是入射到散射體上的電磁波在散射體上各點產生的反射波的疊加，而由于散射體的反射面由1組寬度相同而深度不同的阱有序排列組成，因而來自于不同反射點而到達同一點后，各自相位是不同的（因各反射波的波程不同）。由于被關注區域內各點的場均由這樣1組相位不同的反射波疊加而成，因而各點的場的統計效果相同。

本文采用MLD是因為它是上述3種散射體中結構形式最簡單的。MLD結構依賴于一個偽隨機序列，而且這一偽隨機序列是伽羅瓦域（Galois Field）上1個數學群的基［7］。為了更好地對 MLD的設置過程予以理解，此處有必要對伽羅瓦域予以簡單介紹。

伽羅瓦域是個有限域，通常用G（p）表示，其中p為質數且稱之為該域的階。有限域的一個重要性質是每個有限域G（p）至少要包含一個叫做α的本原元素，它能生成該域中的每個元素。實際中可以將G（p）延伸為一個含有pm個元素的域，這稱為G（p）的擴展域，表示為G（pm），m 是一個非零的正整數。G（p）是G（pm）的子集。此處取p＝2，則G（2）為二進制域，除了數字0和1，在擴展域中還有本原元素α，則該域中的元素由｛0，1，α｝形成，根據文獻［7］知擴展域G（2m）的元素由下面序列構成：

式中：該有限域中的元素個數為2m，非零元素個數為2m－1。

以G（24）［8］為例，即 m＝4，其編碼規則如表1所示，其中α－power表示擴展域G（24）中元素α的指數，4－tuple表示與α－power相對應的4位二進制編碼，從表中可知α15＝α0，這表明擴展域G（24）中僅有15個線性無關的非零元素，其它的非零元素可用這15個線性無關的非零元素構造而成，即αn＝

1.2 MLD結構設計

圖1給出了MLD的結構和工作原理示意圖。其中各參數的具體設置如下述：

在MLD結構中，所有阱的深度dn是一致的：

MLD結構中各個阱的寬度不同，其最大的阱寬度為：

表1 G（24）對應編碼

而MLD中其余各阱寬wa則由伽羅瓦域上的基決定，在MLD阱寬wa確定后，MLD阱之間的距離wb也相應地被確定下來。

上面給出了MLD各個參數的約束條件，下面將給出工作頻率范圍為300MHz～2GHz、各散射體的總寬度相同（WT＝97.5cm）、長度一致（L＝100cm）的具體參數設計步驟：

（1）根據文獻［8］可知，MLD的有限頻率范圍為1倍頻，由于此處的上限頻率為2GHz，故只需設計工作在1GHz頻點上的MLD參數便可滿足300MHz～2GHz頻率范圍要求。由于1GHz工作頻率的波長λ0＝30cm，從而由前述方法可確定MLD阱的最大寬度wamax＝0.25λ0＝7.5cm；又由于總寬度WT已經確定為195cm，故在WT中所含編碼的個數N＝WT／wamax＝26。

（2）根據前面對伽羅瓦域的介紹及表1可知，G（24）是采用4位二進制進行編碼的，則每位編碼所對應的阱寬 wn＝0.25×wamax＝0.25×7.5cm＝1.875cm。再由最大的阱寬度不能超過

0.25 λ0，即式（4）的約束條件可得，隨機數的產生應滿足一定要求，將其對應到G（24）編碼上則要滿足其不能同時出現5個1，即不能是11111，與之相對應的元素α的指數則滿足如下8個約束條件：

（a）3后面不能為12；

（b）6后面不能為10、12；

（c）7后面不能為12；

（d）9后面不能為12；

（e）11后面不能為4、7、10、11、12；

（f）12后面不能為0、4、7、8、10、12、13；

（g）13后面不能為10、12；

（h）14后面不能為12。

（3）生成1組滿足條件（1）和（2）的26個隨機正整數，通過前面分析可知，所有的整數均可由G（24）中的15個線性無關的非零元素構成，則產生的這26個隨機整數的編碼可通過對15取模后再將其映射到伽羅瓦域G（24）編碼上而得（表2給出了這么1組偽隨機數及其對應的元素α的指數和編碼），最后將這1組編碼對應成MLD各個阱寬，現將該隨機數編碼下的MLD實際結構示意于圖2。

圖2 MLD實際結構

2 FDTD建模

在此，以尺寸7.6m×5.2m×3.3m的混響室為例進行仿真分析。圖3為仿真計算中的3種類型混響室的結構示意圖：圖3（a）為裝有攪拌器的混響室結構示意圖；圖3（b）為在6個壁面上設置有散射體的混響室的結構示意圖；圖3（c）為在8個頂角上均設置有散射體的混響室結構示意圖。

本文采用時域有限差分（FDTD）方法對混響室內的場進行仿真分析。為了獲得較高的計算精度，網格大小取為Δx＝Δy＝Δz＝2.5cm；同時，為了滿足數值穩定性的要求，時間步長取為Δt＝4.81×10－11s；為使場達到穩定，所選取的總時間步數則根據散射體及其所處的位置而定。混響室壁和散射體選用相同的良導體材料，其電導率σ＝9×106S／m，相對介電常數εr＝1.0，混響室腔壁厚1cm。分析時采用完美匹配層吸收邊界條件（UPML），計算區域沿x、y、z方向的網格數為324×228×152，其中，每個方向上包括10個網格厚度吸收邊界層；在預設精度下的均勻區的尺寸為100Δx×148Δy×72Δz。波源輻射出正弦平面電磁波，其幅值為1V／m，波源的位置設置在點（44，44，50）處。

表2 滿足條件的1組偽隨機數及其G（24）編碼

3 仿真計算及結果分析

為了進行對比性實驗，分析了f＝915MHz，以下3種情況時混響室內場的特性（如圖3所示）：情況1，只設置有攪拌器的混響室；情況2，散射體位于混響室的6個面上；情況3，散射體位于混響室的8個角落處。

3.1 電場強度的歸一化標準偏差

仿真分析中計算了3種情況下混響室內8個點處的電場強度Ex、Ey、Ez和｜E｜的值，這8個點為：（1）（0.75m，0.75m，0.75m）；（2）（0.75m，4.475m，0.75m）；（3）（0.75m，0.75m，2.55m）；（4）（0.75m，4.475m，2.55m）；（5）（3.25m，0.75m，0.75m）；（6）（3.25m，4.475m，0.75m）；（7）（3.25m，0.75m，2.55m）；（8）（3.25m，4.475m，2.55m）。這8個點是依據IEC61000－4－21測試標準要求選取的，其中，每2個點相距0.4m，且離地面有0.7m的距離。歸一化標準偏差如表3所示。從表中可以看出在只設置有攪拌器情況下混響室內的歸一化標準偏差略小于IEC61000－4－21標準所規定的3dB要求，亦即表明其內部場均勻滿足標準要求；散射體位于6個面和8個角落處時，混響室內部的場都是均勻的，且散射體位于8個角落時混響室內的場比散射體位于6個面時混響室內場的均勻性要好。

表3 混響室內電場強度的歸一化標準偏差σ（dB）

3.2 電場強度的累積概率分布

電場強度Ex、Ey和Ez3個分量的累積概率分布函數（CDF）如圖4所示。所分析計算的這3種情況下的混響室的內部空間都是一致的，均為7.6m×5.2m×3.3m，且電場強度Ex、Ey和Ez3個分量的值取的都是z＝82時平面x－y內的值。從圖中可以明顯看出，圖4（b）與圖4（c）中Ex、Ey和Ez的CDF曲線更接近理論值，說明其均勻性比圖4（a）好，圖4（c）的CDF比圖4（b）的均勻性要稍微好一點。這就表明本文所設置散射體的位置可引起更均勻的場分布，且比文獻［2］中的場均勻性要好。

圖3 仿真計算中3種類型混響室的結構圖

圖4 電場強度的累積概率分布

4 結束語

本文提出了混響室中MLD散射體的一種新的位置，并采用FDTD對其所產生的場均勻性進行了仿真分析，將其與設置攪拌器的混響室以及文獻［2］中的散射體位于混響室的6個面的場均勻性進行了對比，驗證了該位置的可行性。散射體放置于混響室的角落處不僅有效提高了混響室中的場均勻性，還大大減少了利用混響室進行測試的時間，從而為混響室的設計和測試提供了可靠的理論指導。

［1］IEC－61000－4－21，Testing and Measurement Techniques Reverberation Chamber Testethods［S］.

［2］Rhee J G.Field uniformity characteristics of an asymmetric structure reverberation chamber by FDTD method［A］.Asia－Pacific Conference on Environmental Electromagnetics［C］，2003：426－429.

［3］Amaut L R.Operation of electromagnetic reverberation chambers with wave diffractors at relatively low frequencies［J］.IEEE Transactions on Electromagnetic Compatibility，2001，43（4）：637－653.

［4］Petisch M，Schwab J S.Investigation of the field uniformity of a mode－stirred chamber using diffusors based on acoustic theory［J］.IEEE Transactions on E－lectromagnetic Compatibility，1999，41（4）：446－451.

［5］何鵬，蔣全興，周香，等.聲學散射體對混波室場均勻性的影響［J］.微波學報，2009，25（2）：38－41.

［6］Antred M，Chroeder R S.Progressin architectural acoustics and artificial reverberation：concert hall acousticsand number theory［J］.Journal of The Audio Engineer Society，1984，32（4）：94－203.

［7］Petirsch M，Schwab A J.Investigation of the field uniformity of a mode－stirred chamber using diffusers based on acoustic theory［J］.IEEE Transactions on E－lectromagnetic Compatibility，1999，41（4）：446－451.

［8］Schroeder M R.Number Theory in Science and Communication［M］.Berlin，Germany：Springer－Verlag，1944.