鞠文靜 周忠元 蔣全興 景莘慧 周 香

(東南大學機械工程學院,南京 211189)

?

基于有限積分法的電磁兼容吸波材料反射率的建模仿真

鞠文靜 周忠元 蔣全興 景莘慧 周 香

(東南大學機械工程學院,南京 211189)

針對電磁兼容吸波材料反射率試驗測試成本較高、理論方法難以精確預(yù)測的問題,提出了基于有限積分法的吸波材料反射率的建模仿真方法.使用該方法計算了矩形同軸測試裝置空載反射系數(shù),以及分別加載鐵氧體瓦和角錐泡沫吸波材料后的反射系數(shù).吸波材料的介電常數(shù)和磁導(dǎo)率使用二階通用色散模型進行擬合.對于鐵氧體瓦,仿真得到的反射率略優(yōu)于實測值,與產(chǎn)品提供的反射率曲線相比，在諧振點處相差10 dB.對于角錐泡沫吸波材料,仿真結(jié)果與實測值在有效測試頻率范圍內(nèi)相差7 dB.對仿真結(jié)果與實驗結(jié)果的差別進行分析,證明了該方法的可行性.

吸波材料;電磁兼容;有限積分法

吸波材料在電子通訊等領(lǐng)域具有非常廣泛的應(yīng)用,如在電磁兼容暗室、橫電磁波傳輸室等測試場地和設(shè)備中,吸波材料常被置于金屬壁面或作為終端裝置以減少電磁波的反射[1-3].吸波材料最重要的特性是其吸波性能,反射率是吸波材料吸波性能的重要指標,其大小取決于吸波材料的幾何結(jié)構(gòu)和電磁參數(shù).目前電磁兼容吸波材料主要有鐵氧體瓦吸波材料和角錐狀泡沫吸波材料,這些材料都屬于色散介質(zhì),其電磁參數(shù)與頻率有關(guān),從而使得從理論上預(yù)測其反射率較為困難.工程上主要是通過同軸反射法、拱形法等對其進行測試[4],但是搭建這些測試平臺成本較高,而且每種測試方法都有一些不足之處.因此,采用仿真方法得到吸波材料的反射率成為一種重要途徑[5-6].

目前建模仿真中對吸波材料的處理主要有2種方法:① 采用波導(dǎo)端口替代橫電磁波傳輸室終端的角錐泡沫吸波材料,并在端口上設(shè)置足夠多的模式.周小俠[7]使用這種方法對GTEM小室的場強進行了研究,由于端口的模式吸收是全頻帶的,而在實際GTEM小室的終端使用的角錐狀泡沫吸波材料只對100 MHz以上的高頻微波有較好的吸收效果,因而這種方法計算的場強與實驗測試的結(jié)果存在差異.② 將整個終端的吸波材料整體簡化為多層片狀結(jié)構(gòu),并在每層代入不同的電磁參數(shù).Thye等[8]使用這種方法對GTEM小室終端的吸波材料進行了建模,然而這種方法無法處理不同入射方向電磁波在角錐狀吸波材料間的折射問題.

本文采用有限積分法[9-10]對電磁兼容吸波材料進行實體建模,計算吸波材料的反射率.通過仿真結(jié)果與實測數(shù)據(jù)的對比,發(fā)現(xiàn)此建模方法能夠較為準確地反映吸波材料的吸波性能,有助于電磁兼容暗室和橫電磁波傳輸室的優(yōu)化設(shè)計.

1 建模仿真的基本原理與流程

本文基于有限積分法對電磁兼容吸波材料反射率進行仿真.有限積分法最早由Weiland[9]提出,目前發(fā)展較為完善,可用于解決多種電磁場問題.有限積分法的基本思想是,首先將一個有限計算區(qū)域分割為許多小的網(wǎng)格單元,然后將麥克斯韋方程組在每個網(wǎng)格面上離散,得到完全離散化的麥克斯韋網(wǎng)格方程組.最后把與材料特性相聯(lián)系的介質(zhì)方程代入完全離散化的麥克斯韋網(wǎng)格方程組中,通過數(shù)值方法對所有的電學量和磁學量進行求解.在得到端口的場強值后,通過反射率的基本定義進行計算.

在劃分網(wǎng)格時,本文使用了正交六面體網(wǎng)格,并采用時域求解器進行數(shù)值求解,求解過程中需要考慮方程的穩(wěn)定性,其穩(wěn)定條件為

(1)

式中,Δt為時間步長;Δx,Δy,Δz分別為3個方向的網(wǎng)格步長;ε和μ為材料的介電常數(shù)和磁導(dǎo)率.

圖1給出了本文的仿真建模流程.首先建立一個矩形同軸測試裝置模型,通過計算其空載反射系數(shù)來檢驗?zāi)Ｐ偷恼_性.然后在測試裝置中加載吸波材料,并給定吸波材料的電磁參數(shù).最后通過有限積分法計算端口的反射系數(shù).將仿真建模得到的反射系數(shù)與實驗測試結(jié)果進行比較,用以驗證本文提出的仿真模型的正確性.

圖1 吸波材料反射系數(shù)仿真流程圖

2 吸波材料的性能測試技術(shù)

吸波材料建模仿真的關(guān)鍵是所建模型能夠真實體現(xiàn)材料的吸波性能,仿真結(jié)果能夠與實測值一致.因此必須對吸波材料的反射率進行實驗測試.本文使用東南大學實驗室自行研制的矩形同軸測試裝置(見圖2(a))測試吸波材料在低頻段內(nèi)的反射率[11].測試裝置制作參照IEEE STD 1128—1998[4]標準,總高8 m,配合Agilent公司的5061B矢量網(wǎng)絡(luò)分析儀,測試有效頻段為30～500 MHz.測試時,用同軸電纜將矢量網(wǎng)絡(luò)分析儀連接至同軸測試裝置的輸入端口,并采用時域門技術(shù)濾除由阻抗不連續(xù)引起的反射,然后分別測得同軸測試裝置空載時的反射系數(shù)以及裝載吸波材料后的反射系數(shù).采用這些測試數(shù)據(jù)對本文提出的仿真建模方法進行檢驗.

(a) 實物圖

(b) 仿真模型

3 吸波材料的仿真計算及結(jié)果分析

對矩形同軸測試裝置進行幾何建模,然后將鐵氧體瓦吸波材料和角錐泡沫吸波材料模型分別置于建好的測試裝置模型中進行仿真計算,其中吸波材料的幾何參數(shù)可以任意調(diào)節(jié),材料的電磁參數(shù)以離散點的形式代入模型,作為基本的輸入?yún)?shù).

3.1 矩形同軸線的建模仿真

矩形同軸測試裝置的頂端為圓形同軸接頭,其特性阻抗為50 Ω,而主體段為內(nèi)、外均為矩形的同軸線,其特性阻抗為60 Ω,故在圓形同軸接頭與主體段之間存在一過渡段,且此段結(jié)構(gòu)非常復(fù)雜.但是在實際測試過程中,過渡段中由阻抗不連續(xù)引起的反射可通過時域門技術(shù)濾除,測試中起主要作用的仍是主體段.為了縮短計算時間,將矩形同軸測試裝置簡化為僅有主體段,模型見圖2(b).建模時,設(shè)置矩形同軸測試裝置內(nèi)、外導(dǎo)體橫截面邊長分別為0.6和1.8 m.為了檢驗試驗裝置高度對仿真結(jié)果的影響,仿真時分別設(shè)置測試裝置的高度為4和1 m,仿真結(jié)果顯示測試裝置的高度對反射系數(shù)的影響僅表現(xiàn)在波紋震蕩程度上,并不會改變反射系數(shù)極值點的位置.為了減少計算量,仿真時統(tǒng)一設(shè)置測試裝置的高度為1 m.輸入端口設(shè)置為60 Ω矩形同軸輸入端,然后計算空載反射系數(shù)Γ.

如圖3所示,矩形同軸測試裝置空載時,在0～600 MHz內(nèi),反射系數(shù)接近理論值0 dB,最大谷點處僅為-0.067 dB,證實了測試裝置建模的正確性.

圖3 矩形同軸測試裝置空載反射系數(shù)

3.2 鐵氧體瓦的建模仿真

(2)

式中,f為某一頻率點的相對介電常數(shù)或相對磁導(dǎo)率;f∞為無窮頻率時的相對介電常數(shù)或相對磁導(dǎo)率;α0,α1,β0和β1為擬合參數(shù).采用模型(2)可以有效并準確地描述材料的主要物理特性.

圖4 鐵氧體瓦相對磁導(dǎo)率的輸入值及擬合曲線

圖5中實線是仿真計算出的鐵氧體瓦反射率曲線,與產(chǎn)品提供的反射率曲線相比,兩曲線吻合較好,僅在諧振頻率點處有10 dB的偏差,而此偏差可能與仿真頻率步長有關(guān).而與矩形同軸線法的實測結(jié)果相比,由于測試過程中測試樣品的數(shù)量為228片,片與片之間存在不規(guī)則氣隙,使得實測結(jié)果略差于仿真結(jié)果.

圖5 鐵氧體瓦反射率的比較

3.3 角錐泡沫吸波材料的建模仿真

首先建立角錐泡沫吸波材料的幾何模型,其中角錐材料總高為600 mm,底座高為230 mm,底邊邊長為200 mm.然后將角錐材料放入矩形同軸測試裝置模型中,并代入相應(yīng)電磁參數(shù).

圖6 角錐泡沫吸波材料相對介電常數(shù)輸入值及擬合曲線

圖7為角錐泡沫吸波材料實測結(jié)果與仿真結(jié)果對比.仿真曲線在低頻段出現(xiàn)的谷點頻率恰好對應(yīng)擬合模型虛部的峰值點頻率.2條曲線在100～350 MHz頻率段重合較好,但在350～600 MHz頻率段最大偏差約為7 dB,這與仿真模型的簡化和電磁參數(shù)測試的誤差有關(guān).由于在500 MHz以下的頻率段中,吸波材料的反射率小于-25 dB已滿足設(shè)計需要,因而仿真結(jié)果和測試結(jié)果之間的偏差可以接受.

圖7 角錐泡沫吸波材料反射率的比較

4 結(jié)語

本文通過采用有限積分法求解麥克斯韋電磁方程組,較精確地解決了常用吸波材料的建模仿真問題.首先對測試裝置進行建模,然后根據(jù)所研究吸波材料的實際形狀建立幾何模型,接著代入相應(yīng)的電磁參數(shù),最后通過求解電磁方程得到所研究吸波材料的反射率.使用該方法對電磁兼容吸波材料鐵氧體瓦和角錐泡沫吸波材料的反射率進行了仿真研究.對于鐵氧體瓦,仿真結(jié)果得到的反射率與產(chǎn)品提供的反射率相比僅在諧振點處相差10 dB,且略優(yōu)于試驗測試結(jié)果;對于角錐泡沫吸波材料,仿真值與實測值在有效測試頻率范圍內(nèi)最大偏差僅為7 dB.證明本文提出的方法是可行的.本文提出的仿真方法和仿真流程對具有復(fù)雜幾何結(jié)構(gòu)的吸波材料的建模仿真非常有效,可將其應(yīng)用于暗室和GTEM室的優(yōu)化設(shè)計.

References)

[1]Pouhe D, Monich G. On the role of essential higher order modes in a GTEM cell [C]//IEEEInternationalSymposiumonElectromagneticCompatibility. Austin, TX, USA, 2009: 13-17.

[2]Pouhe D, Sandstede O, Monich G. A system theoretical approach for investigating the performance of GTEM cell hybrid termination [C]//InternationalConferenceonElectromagneticsinAdvancedApplications. Torino, Italy, 2013: 49-53.

[3]Pouhe D, Sandstede O, Monich G. Optimization of the GTEM cell wideband termination [C]//InternationalConferenceonElectromagneticsinAdvancedApplications. Torino, Italy, 2011:1285-1289.

[4]IEEE Electromagnetic Compatibility Society. IEEE STD 1128—1998: IEEE recommended practice for RF evaluation in the range of 30 MHz to 5 GHz [S]. New York, USA: The Institute of Electrical and Electronics Engineers, Inc., 1998.

[5]Lubkowski G, Ihsan Z, Adami C, et al. Modelling of pyramid absorbers used in EMC facilities [C]//IEEEInternationalSymposiumonElectromagneticCompatibility. Rome, Italy, 2012:1-4.

[6]Xu J, Koledintseva M Y, De S, et al. FDTD modeling of absorbing materials for EMI applications [C]//Asia-PacificSymposiumonElectromagneticCompatibility. Beijing, China, 2010:173-176.

[7]周小俠. 電大尺寸電磁結(jié)構(gòu)的時域仿真實踐[D]. 成都:電子科技大學物理電子學院,2005:20-40.

[8]Thye H, Armbrecht G, Koch M. Pulse propagation in gigahertz transverse electromagnetic cells [J].IEEETransactionsonElectromagneticCompatibility, 2009, 51(3): 592-603.

[9]Weiland T. A discretization method for the solution of Maxwell’s equations for six-component fields [J].ElectronicsandCommunication, 1977, 31: 116-120.

[10]Clemens M, Weiland T. Discrete electromagnetism with the finite integration technique [J].ProgressinElectromagneticsResearch, 2001, 32: 65-87.

[11]彭鵬,周忠元,蔣全興.射頻吸波材料同軸測試裝置改造及其性能比較[J]. 電子機械工程,2009, 25(4):4-6. Peng Peng, Zhou Zhongyuan, Jiang Quanxing. Performance analysis of improved coaxial test equipment for absorber testing [J].Electro-MechanicalEngineering, 2009, 25(4): 4-6. (in Chinese)

[12]Jing S H, Jiang Q X. Transmission/reflection method based on coaxial line for RF materials characterization measurement [J].JournalofAstronautics, 2005, 26(5): 630-634.

FIT-based simulation of reflectivity of absorbing materials for electromagnetic compatibility

Ju Wenjing Zhou Zhongyuan Jiang Quanxing Jing Shenhui Zhou Xiang

(School of Mechanical Engineering, Southeast University, Nanjing 211189, China)

Considering the high cost in experimental measurement and the difficulty of accurate prediction by theoretical method, a simulation method based on finite integration technique (FIT) is proposed to obtain the reflectivity of absorbing materials for electromagnetic compatibility. By using this method, the reflectivity of rectangle coaxial device without loading, and the reflectivity with loading either ferrite tiles or pyramid foam absorbers are calculated. The general second-order dispersion model is used to fit the permittivity and permeability of the absorbing materials. For ferrite tiles, the reflectivity from simulation results is slightly superior to experimental results, and the difference of reflectivity around the resonance frequency is 10 dB between the simulation results and the reflectivity curve provided by the producer. For pyramid foam absorbers, the difference is only 7 dB within effective testing frequency range. The difference of the results from simulation and experiment is analyzed. The practicability of this method is verified.

absorbing materials; electromagnetic compatibility; finite integration technique

10.3969/j.issn.1001-0505.2015.03.011

2014-11-04. 作者簡介: 鞠文靜(1985—),女,博士生;周忠元(聯(lián)系人),男,博士,副教授,博士生導(dǎo)師,zyzhou@seu.edu.cn.

總裝預(yù)研基金資助項目(9140A33040114JW06001).

鞠文靜,周忠元,蔣全興,等.基于有限積分法的電磁兼容吸波材料反射率的建模仿真[J].東南大學學報:自然科學版,2015,45(3):474-477.

10.3969/j.issn.1001-0505.2015.03.011

TN813

A

1001-0505(2015)03-0474-04