利用有限元法分析EMC暗室的NSA_图文

1、利用有限元法分析EMC 暗室的NSA孟东林,沙斐(北京交通大学电子与信息工程学院电磁兼容实验室,北京100044摘要:EMC 暗室是最常用的电磁兼容测试场地,常受限于低频性能。为了在建造EMC 暗室之前,对其低频性能作准确的预测分析,利用平面分层等效模型来模拟暗室中的角锥吸波体,研究了暗室的不同墙面铺设不同吸波材料时其性能的变化情况。仿真结果表明,适当增加钝口吸波体的高度能够显著改善暗室性能;在暗室4个墙位置中,处于第一菲涅尔区的侧墙(与电波传播方向平行位置的吸波材料更加关键。实测数据的验证,表明上述方法分析暗室低频性能时准确可靠。关键词:电波暗室;射线追踪法;归一化场地衰减;吸波体;有限元法

2、;电磁兼容中图分类号:TM935;U285文献标志码:A 文章编号:100326520(20070920041205基金资助项目:北京交通大学“十五”重大科技专项基金项目资助(DX J 04002。Project Supported by Key Science and Technology Foundation of Beijing Jiaotong University During t he 10t hFive 2yearPlan(DX J04002.Analysis of NSA of Anechoic Chamber with FEMM EN G Dong 2lin ,SHA Fei

3、(EMC Lab ,School of Electronic and Information Engineering ,Beijing Jiaotong U niversity ,Beijing 100044,China Abstract :An accurate method was introduced to simulate and analyze normalized site attenuation (NSA of anechoic chambers (AC in the low f requency range (30100M Hz with the finite element

4、method (FEM .The practical pyramid absorbers were modeled by way of effective permittivity and planar layers ,in order to reduce the mesh size of simulation.Since FEM was good at analyzing high loss ,f requency 2dependent and anisotropic media ,the FEM was adopted to simulate NSA of ACs ,because NSA

5、 is a very important criterion for measuring the performance of a semi 2AC.The principle ,details of simulation ,as well as programming flow chart was provided.The reflectivity of two heights of pyramid absorbers was also provided for convenience.Next the effects of different heights of absorb 2ers

6、lined in different side walls of an AC were studied with this method in order to find the right configuration of ab 2sorbers before constructing an AC.The simulation results show that NSA of ACs can be improved by lengthening pyramid absorbers.The absorbers in the first Fresnel zone are key to ACs ;

7、thus its performance can be improved just by optimize the absorbers in this zone.Finally ,the simulation results are verified by practical measurement da 2ta.In short ,the performance of an AC at low f requency can be simulated with this method in a personal computer.K ey w ords :anechoic chamber ;o

8、pen area test site (OA TS ;NSA ;absorbers ;finite element method ;electromagnetic compatibility0引言电磁兼容(EMC 暗室是最重要的电磁兼容测试场地之一,其重要性体现为:EMC 暗室不仅是国际标准和中国国家标准规定的电磁兼容测试的基本场地之一,而且建造EMC 暗室的费用比较昂贵123。事实上,EMC 暗室集中体现了多种电磁兼容技术。EMC 暗室属于修正型半电波暗室,即电磁屏蔽壳体的四周墙壁和天花板均安装吸波材料,地面铺设良导体时,构成半电波暗室,本质上相当于半自由空间;目的是模拟开阔试验场,用于30

9、M Hz 18GHz 的辐射发射测试;相应的指标为归一化场地衰减(Normalized Site Attenuatio n 2NSA ,标准有对于暗室的设计和建造527,已经由早期的先试验再整改的方法逐步过渡到采用射线追踪法在建造暗室之前对其性能作预测分析8。射线追踪法是一种高频近似方法,在低频(100M Hz 时分析暗室的性能误差较大9,因此人们采用各种全波分析方法,如矩量法(Met hod of Moment s 2MoM 10,11,有14第33卷第9期2007年9月高电压技术High Voltage EngineeringSep.2007限元法(Finite Element Met h

10、od2FEM12,时域有限差分法(Finite2Difference Time2Domain2 FD TD13215等。其中文10,11推导出均值法计算公式;文12仅对铺设铁氧体瓦的情况作了仿真; Ignacio根据模分析提出了一种安装了铁氧体瓦时的快速算法16,但是忽略了铁氧体瓦上层的聚氨酯角锥吸波材料的效果,对于没有铁氧体瓦仅安装角锥的暗室无法仿真。本文介绍了采用有限元法(FEM来分析半电波暗室的归一化场地衰减(NSA,并用该方法研究暗室不同墙面的吸波体对NSA的影响情况。1暗室全波分析方法分层等效模型v=(1-g0+gah=0(1+g2(a-0(1+g0+(1-gav=(1-g0+gah

11、=0(1+g2(a-0(1+g0+(1-ga,(1其中0是真空中的介电常数,而a=ar-jal是吸波材料的介电常数6(ar、al分别为a实部、虚部;0是真空中的磁导率,a= ar-jal是磁性吸波材料的磁导率(ar、al分别为a实部 、虚部;g代表吸波材料的占空比,v、h分别表示法向和切向的磁导率,v、h分别表示法向和切向的介电常数。法向为角锥方向,切向为与吸波材料角锥方向相垂直的方向。式(1中法向公式是精确的,切向公式是近似的,误差为5%。公式(1的本质是当频率较低时,角锥结构可以采用各向异性的多层平面来近似,即每一平面分层的材质都是各向异性的,且各层材质介电常数都不相同。原因是,沿角锥方向

12、是慢变化,而沿与角锥相切的方向,是周期性的快变化。根据该公式,利用有限元法算得的实际吸波材料和等效分层以后的反射率比较见文18,关于场均匀性的有限元法仿真情况见文19。分层等效模型见图1。显然,上述等效分层模型是高色散、高损耗和各图1暗室仿真的等效分层方法示意图Fig.1Schem atic of effective layers for modelinganechoic chamber图2电波暗室有限元法仿真流程Fig.2F low ch art of simulating anechoic ch amb er by FE M向异性的媒质。划分层数较多时比较接近实际吸波材料,但是计算时所需内

13、存也越大。通常需要根据材料高度分成1020层,误差分析见文18。仿真流程见图2。2仿真实例和NSA影响因素研究半电波暗室内部共有5个面安装吸波体;因此需要对比研究,哪一个墙面上的吸波最重要。众所周知,吸波材料的作用是吸收电磁波使得反射出去的电磁波能量密度尽量小;吸波材料越长,其反射率越低;且正入射反射率优于斜入射反射率20。然而当吸波材料变长以后(即图3中的L+ D值变大,入射角将变大(因为发射天线T x和接收天线R x的间距必须恒定不变,见图3。入射角变大将导致吸波材料的反射率R变差,f为频率,见图4。其中P TM为平行极化波,指入射波的电场极化方向平行于入射面(平面波的波矢量方向与分界面法

14、向构成的平面称为入射面,P TE为垂直极化波,指入射波的电场极化方向垂直于入射面。EMC164吸波体尺寸为:底座厚度D=20cm、底座24 图3侧墙吸波材料长与入射角的关系Fig.3R elationship betw een absorber sheight and incident angle图4EMC164和EMC122反射率的比较Fig.4Comparison of reflectivity betw eenEMC 164and EMC122边长A =30cm 、角锥顶部钝口边长a =10cm ,角锥高度L =144cm ;EMC122吸波体尺寸为:D =20cm 、A =30cm 、

15、a =10cm ,L =102cm 。增加吸波体高度,一方面其反射率下降,吸波性能变好;另一方面暗室中电磁波的入射角将变大,导致吸波体的性能下降。故需研究屏蔽壳体尺寸不变时,增加侧墙上吸波体的高度是否有用。针对屏蔽壳体尺寸为1214m 718m 6m 的半电波暗室作仿真,逐一分析影响暗室NSA 的因素。半电波暗室NSA 的实际测量方法可参阅文4。本文仅给出半电波暗室两种情况下长度方向的中轴线位置上的NSA 值。1暗室地面为理想电导体(PEC ,侧墙安装吸波体(分别为EMC164或EMC122,其余3个面均为理想吸波体。此时的NSA 偏差D 见图5,D 为实际的NSA 值与理论的NSA 值之差。

16、其中1m 、TM 指发射天线距离金属地面高度为1m ,水平极化方式;1m 、TE 指发射天线距离地面1m ,垂直极化方式。其余类推。对比图5(a 和(b 可见,当吸波体加长以后,暗室中的NSA 性能变好,即增加侧墙上吸波体的高度有积极作用。2当端墙安装吸波体,地面为PEC ,侧墙为理想吸波体时的NSA 偏差D 见图6。图6中的(a 和图5侧墙吸波材料性能对暗室NSA 的影响Fig.5R ole of side 2w all absorbers on NSA图6端墙吸波材料性能对暗室NSA 的影响Fig.6E ffect of end 2w all absorbers on NSA(b 的结果再

17、次说明:当吸波体加长以后,暗室的NSA 性能会变好,即吸波体高度增加对NSA 有积极作用。分别对比图5和图6中的(a 以及(b 可见,在安装同样的吸波体条件下,显然侧墙的吸波体性能更关键,即改善侧墙吸波体的性能,能更加显著的改善半电波暗室的NSA 性能。地面为PEC ,其余5个面均安装EMC164吸波体以后,得到的NSA 偏差,见图7。显然,该性能不达标。本文仿真中相当于把吸波体安装在由理想导体(PEC 构成的屏蔽壳体中,但是实际的吸波体总是342007年9月高电压技术第33卷第9期安装在屏蔽铁壳上的,屏蔽铁壳在低频时对磁场具有一定吸收效果,因此铺设上述吸波体后,半电波暗室的实际效果会比仿真效

18、果稍好一点,当然,为了提高半电波暗室性能,可采用更长的吸波体,或采用复合吸波体。Holloway等采用FD TD方法在时域作了仿真10。本文的结果是针对1214m718m6m尺寸的半电波暗室(显然,该屏蔽壳体比一般3m法电波暗室的屏蔽壳体大很多,是在32位PC机的限制下完成的,单个线程占用最高内存达到2G B。图8是屏蔽体1260cm780cm600cm的某电波暗室的实测数据。图8(a是刚建造完的NSA 偏差数据,发现在低频(37M Hz不达标,后侧墙采用反射率降低23dB的新吸波材料,其余条件不变所得实测数据,发现NSA性能明显好转(图8 (b。该实测结果也证明,改善侧墙吸波材料的性能,能够

19、显著提高半电波暗室的性能。3计算方法比较电波暗室通常采用的预测分析方法是射线追踪法,但是在低频时射线追踪法等高频近似法不再有效,导致采用该法分析暗室低频性能的误差很大。针对该问题,开始研究采用全波分析方法来预测暗室的低频性能。然而,对电波暗室直接做全波分析非常困难,因为吸波体是色散、有耗的,而且电波暗室中的吸波体数目庞大结构复杂,无论对计算能力还是内存容量的要求都很高。由于FEM非常适合计算高损耗和高色散媒质,且FEM采用四面体来离散划分网格,使得FEM能够比较准确地处理大型、复杂的结构。因此,本文在采用平面分层等效模型的基础上,采用FEM来计算电波暗室的低频性能。但需注意在低频时采用有限元法

20、来分析电波暗室时,边界条件、模型中的网格划分优良情况等,都会影响计算精度。如果网格划分不合理,将产生1 2dB的误差。此外,对于暗室的NSA仿真分析,相比场均匀性(FU的仿真要困难一些,因此文14216均研究了此问题。本文除了阐释如何采用有限元法来预测分析电波暗室外,还重点研究了暗室的不同墙面上的吸波体的效果,目的是优化暗室吸波体的配置,节省成本。4结论a通常电波暗室的低频性能较差,但是采用射线追踪法等方法预测的误差较大。 采用本文介绍的图7高度164cm的吸波材料安装后所得NSAFig.7NSA of an anechoic chamber linedwith absorber EMC164

21、图8某电波暗室的实测NSAFig.8Measured NSA of an anechoic chamber方法能够准确的预测半电波暗室的低频性能。b在半电波暗室的各个墙面上吸波体的重要程度并不完全相同,侧墙上的吸波体更加重要。c目前出现了许多新型吸波体,如空心吸波体,甚至空壳状吸波体,而本文主要分析聚氨酯角锥吸波体。如何采用FEM来仿真分析安装了这些新型吸波体的电波暗室的性能,有待进一步研究。参考文献1Lin G,Alvarado M J.Construction of a102m semianechoicchamberJ.Compliance Engineering,2001,17(7:62

22、.2Ray P D,Green J.T rials and tribulations of building a modern EMCfacilityC.IEEE International Symposium on Electromagnetic Compatibility:1.Washington,DC,USA,2000:2752280.3陈嵩,沙斐,王国栋,等.3m法半电波暗室改造施工案例J.安全与电磁兼容,2005(增刊:34236.CH EN Song,SHA Fei,WAN G Guo2dong,et al.Reconstruc2 tion of a3meter semi2anec

23、hoic chamberJ.Safety&EMC, 2005(S:34236.4ANSI C63.422003American national standard for met hods of44measurement of radio2noise emissions from low2voltage electrical and electronic equipment in t he range of9k Hz to40GHzS, 2003.5Kurihara H,Yanagawa M,Suzuki Y,et al.Realization of3msemi anechoic chambe

24、r by using hybrid em wave absorber con2 sisting of composite magnetic materialJ.IEICE Transactions on Electronics,2006,E892C(1:69271.6Suzuki H,Nishikata A,Higashida Y,et al.Free space meas2urement met hod wit h parallel electromagnetic wave beam by u2 sing dielectric lenses and horn antennas for ref

25、lectivity of elec2 tromagnetic absorbers millimeter wavesJ.IEICE Transactions on Electronics,2006,E892C(1:24229.7Y oon J,Chatterjee I,McPherson D,et al.Design,character2ization,and optimization of a broadband mini exposure chamber for studying catecholamine release from chromaffin cells exposed to m

26、icrowave radiation:Finite2difference time2domain technique J.IEEE Trans on Plasma Science,2006,34(4:145521469.8Anzai H,K ishimoto T,Y amazaki T,et al.Analysis of semi2ane2choic chamber using ray2tracing techniqueC.Electromagnetic Compatibility Symposium Record.Santa Clara,1996:1432145. 9Kawabata M,I

27、shida Y,Shimada K,et al.Comparison of siteattenuation analysis result s between FD TD and ray2t racing met hod using compact anechoic chamberJ.IEICE Transac2 tions on Communications,2005,E882B(8:315223157. 10Kuester E,Holloway C.A low2frequency model for wedgepyramid absorber arrays2I:t heoryJ.IEEE

28、Transactions on EMC,1994,36(4:3002306.11Holloway C,Kuester E.A low2frequecy model for wedge orpyramid absorber arrays II2computed and measured result sJ.IEEE Transactions on EMC,1994,36(4:3072313.12Baker A E,Sitzia R M,Preston I,et al.Characterisation of e2lectromagnetic anechoic chamber by finite e

29、lement met hodJ.IEEE Transactions on Magnetics,1996,32(3:151321516. 13Askri A,Vollaire C,Nicolas L,et al.Normalized site attenua2tion standard correction from numerical computing,magnetics J.IEEE Transactions on Magnetics,2002,38(2:6932696.14Gradoni G,Moglie F,Pastore A P,et al.Numerical and ex2peri

30、mental analysis of t he field to enclosure coupling in reverber2ation chamber and comparison wit h anechoic chamberJ.IEEE Trans on Electromagnetic Compatibility,2006,48(1:2032211.15Christopher L Holloway,Paul M Mc Kenna,Roger A Dalke,et al.Time2doman,modeling,characterization,and measure2 ment s of

31、anechoic and semi2anechoic elect romagnetic test cham2 bersJ.IEEE Trans on EMC,2002,44(1:1022118.16Prez I M,Nuo L,Pereira F D Q,et al.Low2frequency modelfor rectangular semi2anechoic chambersJ.IEEE Transactions on Electromagnetic Compatibility,2006,48(4:7252733. 17Holloway C L,DeLyser R R,G erman R

32、F,et al.Electromagneticabsorbers used in anechoic and semi2anechoic chambers for emissions and immunity testing of digital devicesJ.IEEE T ransactions on E2 lectromagnetic Compatibility,1997,39(1:33247.18MEN G Donglin,SHA Fei,WAN G Guodong.Simulation ofanechoic chambers t hrough effective electromag

33、netic propertiesC.7t h Internatioal Conference on Electronic Measurement&Inst rument s.Beijing,2005,1:3892392.19孟东林,沙斐,王国栋.近距离观察电波暗室的场均匀性J.微波学报,2007,23(3:66270.MEN G Dong2lin,SHA Fei,WAN G Guo2dong.A closer look at t he field uniform of anechoic chambers2FEM simulationJ.Journal of Microwave,2007,23(

34、3:66270.20孟东林,沙斐,王国栋,等.电波暗室用吸波材料的性能优化和选用J.安全与电磁兼容,2005(增刊:59262.MEN G Dong2lin,SHA Fei,WAN G Guo2dong,et al.Optimi2 zation design and choice of absorbers in anechoic chambersJ.Safety&EMC,2005(S:592 62.孟东林孟东林1973,男,博士生,研究方向为计算电磁学、吸波材料、电磁兼容暗室和混响室等。沙斐1946,男,教授,博导,研究方向为复杂大系统间的电磁干扰及故障诊断与排除。收稿日期2006212219编辑陈蔓(上接第27页12舒立春,蒋兴良,田