平面波照射下无限大导体板上周期孔阵屏蔽效能的解析研究*

2019-06-04白婉欣李天乐郭安琪成睿琦焦重庆

物理学报 2019年10期

白婉欣李天乐郭安琪成睿琦焦重庆

(华北电力大学，新能源电力系统国家重点实验室，北京 102206)

1 引言

平面波对无限大导体板上开孔阵列的电磁透射属于经典电磁问题,在电磁屏蔽、频率选择表面等方面有实际应用[1-11].理论上,该问题可以通过数值计算方法求解[12-16],例如,文献[16]采用矩量法计算了近场辐射源透过无限大导体板的开孔阵列时的电磁耦合特性.然而数值方法实施相对复杂,且耗费CPU和存储资源较多,特别是在需要实时仿真或反复计算以便进行大量参数优化的场合.解析方法求解虽然适用范围窄,但具有唯一性、物理意义清晰、易于计算,特别是具有可作为数值方法的检验基准等优点.因此,对于给定的模型,得到其解析解不仅具有理论价值,也具有实际意义[17-20].

本文结合电磁屏蔽应用实际,在电磁波波长明显大于孔阵周期单元边长,且孔阵周期单元边长明显大于开孔尺寸的前提下,研究透射场的解析表达式.目前,对于平面波垂直入射圆形孔阵情况,已有解析结果[21].其中孔阵的影响通过表面阻抗描述,而表面阻抗需要运用分离变量法和平均场方法求解,过程复杂.在电磁波入射角的影响方面,一般简单对照平面波对均匀平板材料电磁透射的入射角依赖规律,缺少严格的理论分析.对于平面波入射金属网情况,文献[22]进行了分析.网线的影响同样通过表面阻抗描述,考虑了平面波入射角的影响.但是该方法普适性有限,无法推广到其他形状开孔板的情况.文献[13]针对开单个孔缝的无限大平面板的磁场屏蔽效能进行了分析.而在实际应用中,孔缝通常以阵列形式出现,该方法无法直接推广至多个孔缝的情况.

本文首先基于Bethe理论[23],将透射场表示成偶极子阵列的辐射场.然后对磁偶极子阵列进行平均化,即将其偶极矩分布于板上,由此得到均匀磁化强度.进而,将磁化强度用平板两侧的磁化面电流表征,然后可以容易地求出面电流所产生电磁场的解析表达式.对于电偶极子阵列(仅存在于TM极化斜入射情况),通过研究电偶极辐射和磁偶极辐射的对应关系,从而依据磁偶极阵列辐射场导出电偶极阵列辐射场,最终通过二者叠加得出总的透射场.所导出的解析解理论上适用于所有形状磁/电极化系数已知的开孔.对于诸如椭圆、圆形孔,其极化系数是已知的.对于其他形状的开孔,通过数值计算可以获得极化系数,如矩形开孔等已有相关结果发表.

2 理论

2.1 垂直入射

图1为平面波垂直入射导体板示意图,入射磁场Hi=-Hie¯jkzey,入射电场Ei=-Eie¯jkzex,平面波波矢量k=ke,其中,相邻z两孔x方向距离为d1,y方向距离为d2,板厚度为t.假设坐标原点位于板上表面.

图1 平面波垂直入射开孔导体板示意图Fig.1.Conductor plate with holes illuminated by a plane wave with vertical polarization.

根据Bethe小孔理论[23],从开孔耦合到z＞ 0区域的电磁场,可借助位于开孔中心的等效偶极矩为p,m的电偶和磁偶极子来描述：

式中αe为电极化系数,αmx为小孔沿x方向磁极化系数,αmy为小孔沿y方向磁极化系数.En为小孔作短路处理(即小孔不存在或被导体封堵)时开孔面上电场的法向分量,Htx,Hty分别为短路处理下开孔面上切向磁场的x,y分量.对于垂直入射情况,Htx=0,Hty=2Hi,En=0 .此时,只有磁偶极子作用,z＞ 0区域的电磁场近似为沿导体板上分布的磁偶极阵列的辐射场.依据镜像原理[24],导体板可以去除,同时磁偶极矩加倍.

通过平均化的方法,将每一个磁偶极子的磁偶极矩均匀分布在一个周期单元内,从而得到板的均匀磁化模型,其等效磁化强度为

其中S代表一个周期单元的面积,ΔV=St.对于图 1所示情况,S=d1d2.所以板上、下表面磁化面电流密度为：

当导体板不开孔时,板表面传导面电流密度为

对比(4)式—(6)式可得

不难理解,z＞ 0区域,(6)式激发的电磁场与入射场完全抵消.进而依据(7)式可得,(4)式和(5)式所述电流产生的磁场分别为：

其中z1,z2分别是场点到板上表面和下表面的距离.假设场点P的坐标为 (xp,yp,zp)且zp＞0,则z1=zp,z2=zp-t.场点P处透射场的磁场强度为

由于kt趋近于0,所以 1-ejkt≈-jkt,(10)式可写为

就此得出,平面波垂直入射开孔金属板情况下,屏蔽效能公式为

2.2 斜入射-TE极化

如图 2所示,TE极化平面波入射角为θ,入射磁场Hi=(-Hicosθex+Hisinθez)e¯jk·r,入射电场Ei=Eie¯jk·rey,波矢量k=ksinθex+kcosθez.

图2 TE极化平面波入射均匀开孔导体板示意图Fig.2.Conductor plate with holes illuminated by a plane wave of TE polarization.

根据TE极化斜入射情况,En=0,Ht=-2Hicosθe¯jkxsinθex.此时,只有磁偶极子作用.考虑镜像加倍效应后,等效磁偶极矩为

同理,可得上、下表面磁化面电流密度为

当导体板不开孔时,板表面传导面电流密度为

同理,通过分别比较(14)式、(15)式和(16)式的比例关系,可得(14)式和(15)式所述电流在z＞ 0区域产生的磁场分别为：

同理,叠加可得场点P处磁场强度为

磁场屏蔽效能为

由电场、磁场和传播方向三者关系可得,场点P处电场强度为

2.3 斜入射-TM极化

如图 3所示,TM极化平面波入射角为θ,入射电场Ei=(-Eicosθex+Eisinθez)e¯jk·r,入射磁场Hi=-Hie¯jk·rey.

图3 TM极化平面波入射均匀开孔金属板Fig.3.Conductor plate with holes illuminated by a plane wave of TM polarization.

根据TM极化斜入射情况,En=Eisinθe¯jk·rez,Ht=-2Hie¯jkxsinθey,故电偶极子、磁偶极子同时作用.考虑镜像加倍效应后,等效磁偶极矩、电偶极矩分别为：

(24)式表明,求得电偶极子、磁偶极子产生的x方向电场,即可求出相应的y方向磁场.另一方面,TM极化下,磁场只有y分量,因此y方向磁场即代表总磁场.

令TE极化中磁偶极子激发的z方向磁场、y方向电场分别为Hz(mTE),Ey(mTE),TM极化中电偶极子、磁偶极子激发的x方向电场分别为Ex(pTM),Ex(mTM).可以证明(证明过程见附录)

根据(19)式和(21)式,TE极化下磁偶极子产生的电场Ey和Hz为：

根据(25)式—(28)式可得：

由(24)式、(29)式和(30)式得,TM极化下,电偶极子、磁偶极子产生的磁场分别为：

由(31)式和(32)式得TM极化下透射场的磁场强度为

所以屏蔽效能为

比较(20)式和(34)式看出,TE极化下透射场幅值和入射角的关系比较简单,即与入射角的余弦成正比,因此屏蔽效能随入射角的增大而增加；TM极化下,透射场幅值与入射角的关系相对复杂,随着入射角的增大,屏蔽效能有减小的趋势.当电极化系数相对磁化系数很小时,附加项可以忽略,但在电极化系数和磁化系数接近的情况下,附加项将会对屏蔽效能产生显著影响.

3 结果与讨论

3.1 开孔的极化系数

对于圆孔和椭圆孔,其电、磁极化系数有解析公式,如表 1所列[25].

表1中K为第一类完全椭圆积分,E为第二类完全椭圆积分：

3.2 垂直入射的分析

为了验证上述公式的准确性,本文选定几种情况进行分析.由(12)式可以看出,屏蔽效能的大小与周期单元边长比(d1/d2)无关,为了验证该结论,在周期单元面积S= 16 cm2,开孔为半径r= 0.25 cm圆孔阵列的情况下,选择d1/d2= 1,4,16进行计算,CST软件全波模拟仿真(微波工作室,采用unitcell边界)结果如图 4所示.随着频率的增大,屏蔽效能减小.对于d1/d2= 1和d1/d2= 4两种情况,仿真结果几乎一致,对于d1/d2= 16,仿真结果与其余两种情况相差约1 dB.所以在d1,d2大于最小开孔尺寸的前提下,屏蔽效能与周期边长比关系不大,后续分析均取d1=d2= 4 cm.

在电磁屏蔽领域,实际应用中平面板厚度为mm级.当厚度相对小孔尺寸的比值(以圆孔为例,t/r)增大时,小孔的波导效应(以圆孔为例,SE波导≈16t/r)随之增大,会带来一个附加项的衰减.图 5展示了开孔半径为1 cm,厚度t= 0.001 mm和t= 1 mm两种情况下的公式计算结果和全波仿真结果.根据波导效应的公式,本文考察的算例中,波导效应带来的衰减并不大,两种厚度下的屏蔽效能基本一致.后续分析中,取t= 0.001 mm.

表1 常见开孔形状极化系数Table 1.Polarization coefficients of typical opening shapes.

图4 d1/d2对屏蔽效能的影响Fig.4.Dependence of the SE on frequency for different d1/d2.

图5 板厚度对屏蔽效能的影响Fig.5.Dependence of the SE on frequency for different thicknesses of the plane.

图6展示了圆形开孔下,屏蔽效能随频率的变化情况.可以看出,随着开孔半径的增大,屏蔽效能下降.r= 0.25,0.5 cm时,公式计算结果和仿真结果相差均约为4 dB,r= 1 cm时,两者相差约为2 dB,且频率越高,一致性越好.

对于图 7所示l/w= 10的椭圆形开孔,图 8展示了长轴l分别为0.5,1 和2 cm时屏蔽效能随频率的变化,其中入射波的电场沿x方向.对于长轴l= 0.5 cm和l= 1 cm公式计算结果与全波仿真结果一致性较好,对于长轴l= 2 cm,两者差值随着频率的增大而增大.

图6 圆形孔开孔大小对屏蔽效能的影响Fig.6.Dependence of the SE on frequency for circular apertures with different radius.

图7 椭圆开孔示意图Fig.7.Diagram of elliptical opening.

图8 椭圆孔开孔大小对屏蔽效能的影响Fig.8.Dependence of SE on frequency of different elliptical apertures with different size.

针对如图 9所示的方形(图 9(a))、十字形(图9(b))开孔,文献[26]分别给出了其磁化系数,将磁化系数代入(12)式可得其屏蔽效能.图 10展示了两种开孔形状板的屏蔽效能随频率的变化.对于边长a= 1 cm的方形开孔,全波仿真结果与公式计算结果一致性较好.对于十字形开孔(l= 30 cm,w/l= 0.33),本文比较了入射波频率为0—5 GHz的情况,在频率低于3 GHz时公式计算结果与全波仿真结果比较一致,在3 GHz之后,随着屏蔽效能变差,两种方法的一致性也稍有下降.

图9 方形、十字形开孔示意图Fig.9.Diagram of square and cross opening.

图10 不同形状开孔板的屏蔽效能随频率的变化Fig.10.Dependence of SE on frequency of different elliptical apertures with different size.

对于混合形状开孔的情况,满足各个开孔形状磁化系数相同即可计算其屏蔽效能.由圆孔、椭圆孔磁化系数的表达式可以得到：r= 0.25 cm的圆孔和l= 24.97 mm,l/w= 10的椭圆孔y方向磁化系数相等,所以他们在磁场沿y方向电磁波照射下的屏蔽效能可由(12)式得出.图 11所示为混合形状开孔周期单元示意图,图 12展示了周期单元如图 11所示的开孔板的屏蔽效能随频率的变化,公式计算结果和仿真结果比较一致.

图11 混合形状开孔周期单元示意图Fig.11.Diagram of the unit cell of mix shape opening.

图12 混合形状开孔的屏蔽效能随频率的变化Fig.12.Dependence of SE on frequency of apertures with mix shape.

3.3 斜入射的分析

取开孔为半径0.5 cm的圆孔,波频率3 GHz,图 13给出了屏蔽效能随入射角θ的变化.对于TE极化入射波,公式计算结果和全波仿真结果的差值约为4 dB；对于TM极化入射波,差值约为3 dB.开孔为l= 2 cm的椭圆,屏蔽效能随入射角变化情况如图 14所示.入射波为TE极化时,公式计算结果与全波仿真的差值约为2 dB,入射波为TM模时,公式计算结果和仿真结果差值大部分在5 dB内,且随着入射角的增大而增加.

图13 圆形孔屏蔽效能随入射角的变化（r = 0.5 cm）Fig.13.Dependence of SE on angle of incidence of different polarization of circular apertures (r = 0.5 cm).

图14所示的TM极化中,磁场方向沿y方向,与椭圆开孔方向垂直,(34)式中电极化系数和磁化系数比较接近,附加项影响比较显著,故随着入射角的增大,屏蔽效能先增大后减小.

图14 椭圆孔屏蔽效能随入射角的变化（l = 2 cm）Fig.14.Dependence of SE on angle of incidence of different polarization of elliptical apertures (l = 2 cm).

4 结论

本文提出了一种求解平面波对无限大导体板上周期孔阵的电磁透射场的解析方法.该方法适用的前提条件是：波长大于周期单元尺寸；周期单元尺寸大于开孔尺寸.该方法首先基于Bethe理论将透射场表示为偶极子阵列的辐射场,然后通过平均化将偶极阵列近似成均匀磁化/极化板,进而通过等效磁化面电流求得磁偶极阵列的透射场,以及通过电偶极和磁偶极的对应关系求得电偶极阵列的透射场,最后二者叠加得出总的透射场.依据该方法分别导出了TE和TM两种极化方式下不同入射角时透射场的表达式.通过对几个案例的计算以及和全波仿真的比较验证了公式的有效性,并得到了以下结论.

1)透射场强幅值与波频率、孔的磁化系数成正比,与一个周期的面积成反比.

2)TE极化下,由于仅磁偶极阵列起作用,透射场幅值与入射角的关系十分简单,即与入射角的余弦成正比.此时,屏蔽效能随入射角的增加而变大.

3)TM极化下,磁偶极阵列和电偶极阵列同时起作用.透射场幅值与入射角的关系相对复杂,不仅与入射角的余弦成反比,还需要乘以一个包括了电极化系数以及入射角正弦值的附加项.比较而言,磁偶极比电偶极对透射场的贡献要大.总体上,此时屏蔽效能随入射角的增加而变小.

附录A

对于TE极化,只存在x方向的磁偶极子,一个位于x=na,y=pb,z= 0处磁偶极子m=mTEex在场点(xp,yp,zp)处产生的电场为(只考虑远场)