CN107255755A - 一种计算开孔矩形屏蔽体内部场分布的分区解析模型 - Google Patents

一种计算开孔矩形屏蔽体内部场分布的分区解析模型 Download PDF

Info

Publication number
CN107255755A
CN107255755A CN201710530340.9A CN201710530340A CN107255755A CN 107255755 A CN107255755 A CN 107255755A CN 201710530340 A CN201710530340 A CN 201710530340A CN 107255755 A CN107255755 A CN 107255755A
Authority
CN
China
Prior art keywords
mrow
msub
msup
mfrac
perforate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710530340.9A
Other languages
English (en)
Other versions
CN107255755B (zh
Inventor
焦重庆
何日
胡逸帆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Original Assignee
North China Electric Power University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Priority to CN201710530340.9A priority Critical patent/CN107255755B/zh
Publication of CN107255755A publication Critical patent/CN107255755A/zh
Application granted granted Critical
Publication of CN107255755B publication Critical patent/CN107255755B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08Measuring electromagnetic field characteristics
    • G01R29/0807Measuring electromagnetic field characteristics characterised by the application
    • G01R29/0814Field measurements related to measuring influence on or from apparatus, components or humans, e.g. in ESD, EMI, EMC, EMP testing, measuring radiation leakage; detecting presence of micro- or radiowave emitters; dosimetry; testing shielding; measurements related to lightning
    • G01R29/0835Testing shielding, e.g. for efficiency
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Electromagnetism (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Measuring Magnetic Variables (AREA)

Abstract

本发明公开了属于电磁屏蔽技术领域的一种快速计算开孔矩形屏蔽体内部场分布的分区解析模型。首先将开孔矩形屏蔽体内部的区域划分为小孔区、偶极区、波导区,在不同的区域分别建立解析模型,各区域解析模型分别为:小孔区采用带开孔的无限大金属板等效、偶极区采用偶极辐射模型近似、波导区采用波导模式进行等效。利用上述解析模型可以计算出开孔矩形屏蔽体内部的场分布,进而得出电磁场屏蔽效能,通过与全波仿真的结果进行比较,验证该解析模型的正确性。本发明物理意义清晰,不仅计算速度快、效率高,而且能够得较为准确的描述开孔腔体内部各个分区的场分布规律,对进一步研究开孔屏蔽体内部场分布具有指导意义。

Description

一种计算开孔矩形屏蔽体内部场分布的分区解析模型
技术领域
本发明属于电磁屏蔽技术领域,特别涉及一种计算开孔矩形屏蔽体内部场分布的分区解析模型。
背景技术
电磁屏蔽是抑制电磁干扰的主要技术措施之一,应用范围广泛,常见的方法为将金属腔体应用于设备的外壳,以隔离骚扰源和敏感设备,达到电磁屏蔽的目的。完全封闭的金属腔体具有很好的电磁屏蔽效能,但是在实际应用中,出于散热、通风等的考虑和技术水平的限制,腔体上不可避免地存在着开孔和缝隙,使得屏蔽效能下降,因此,开孔腔体的电磁场屏蔽效能问题是电磁兼容领域的研究热点。
目前关于开孔屏蔽体屏蔽效能的研究主要集中在腔体中部,很少关注开孔附近的屏蔽效能。目前开孔屏蔽体屏蔽效能的计算方法主要有解析理论和数值计算两类,解析理论有已提出的如Bethe小孔耦合理论和等效电路法等,但这些方法只在腔体内部某一部分计算结果较为准确,不适用于整个腔体内部电磁屏蔽效能的计算。而数值方法虽然适用范围广,计算精度高,得到的结果准确,但计算速度慢、时间长,需要装备专业软件,所以发展计算效率高、结果准确,物理意义清晰的解析模型具有实际意义。
发明内容
本发明的目的是提出一种计算开孔矩形屏蔽体内部场分布的分区解析模型,其特征在于,计算平面波垂直于小孔表面入射时,开孔中心轴线上的屏蔽体内部的场分布,将开孔屏蔽体内部空间分解为三个区域:围绕开孔附近的小孔区、靠近屏蔽体背板的波导区、以及小孔区和波导区之间的偶极区,各个区域的场分布模型分别为:小孔区采用带开孔的无限大金属平板模型替代、偶极区采用偶极辐射模型近似、波导区采用波导模式场分布等效。对开孔屏蔽体内部进行分区,各个区域建立场分布计算模型,主要包括:
主要包括:
1)小孔区:在开孔矩形屏蔽体内部,距小孔面1倍直径内的场分布可以等效成无限大金属板附近的场分布,故运用计算无限大金属板开孔附近场分布的方法建立的计算模型如下:
磁场:在平面波垂直于小孔表面入射的条件下,推导出在小孔区,开孔中心轴线上的磁场强度为:
其中H1为当无限大金属板上没有开孔时金属板入射波侧的磁场切向分量,数值上等于入射平面波磁场沿y轴分量H0的两倍,R为小孔的半径,z为观测点到小孔中心的距离;ey指的是y轴的单位矢量
电场:在平面波的垂直入射下,推导出在小孔区,开孔中心轴线上的电场强度为:
其中E0为入射平面波沿x轴分量的电场强度,λ为波长,ρ2=R2+z2,k为波数,ex指的是x轴的单位矢量、ejkρ指的是e的jkρ次方,j为虚数单位,efks指的是e的jkz次方;
2)偶极区:开孔矩形屏蔽体内部,距小孔面倍小孔直径的范围内的电磁场,依据Bethe小孔耦合理论,并结合上述模型及入射平面波的性质,用位于开孔中心的磁偶极矩m的激发场来描述,m则由入射平面波决定,
m=-αmyH0,yey (3)
式中,H0,y为入射磁场HO在开孔面切向的y分量,极化系数αm与开孔的形状、尺寸有关,对于直径为dc的圆形孔的极化系数表达式为
由此得出磁偶极矩:
式中,E0为入射电场沿x轴分量的振幅,η0为自由空间波阻抗,η0=377。
可以得到在偶极区,开孔中心轴线上的磁场强度为:
电场强度为:
式中,SI=2m,为波数,z为观测点到孔的距离;
3)波导区:在开孔矩形屏蔽体内部,距小孔面4倍直径之后的范围内,根据传输线波导理论,将开孔金属腔体看成到开孔金属板距离为d的金属板短路的传输线,由此可得开孔面中心轴线上的电场强度和磁场强度分别为:
E=C1sinkg(d-z)ex (8)
H=C2coskg(d-z)ey (9)
由于偶极区与波导区有一部分区域是重叠的,这部分区域用偶极区计算方法算出的电磁场强度与波导区计算方法算出的相同,因此利用这一点,当z在距小孔倍直径内取值时,令z=z1,联立得到
解得
综上,在波导区,开孔面中心轴线上的电场强度为
磁场强度为
其中,E0为入射平面波沿x轴分量的电场强度振幅,dc为圆形开孔的直径,λ为波长,z为观测点到孔的距离,z1在距小孔3~4倍直径内取值,a为开孔屏蔽体的高度,kg为传播常数,为波数,η0为自由空间波阻抗,η0=377。C1、C2为计算过程设定的两个未知数;
至此,已给出了计算开孔矩形屏蔽体内部场分布的全套方程,按惯例,腔内任一点处的电场屏蔽效能定义为
SE=20log(E0/E) (17)
式中,E0为入射平面波沿x轴方向的电场强度振幅,E为腔内任一点沿x轴方向的电场强度振幅;
腔内任一点处的磁场屏蔽效能定义为:
SH=20log(H0/H) (18)
式中,H0为入射平面波沿y轴方向的磁场强度振幅,H为腔内任一点沿y轴方向的磁场强度振幅。
本发明的有益效果是本发明物理意义清晰,具有计算速度快,计算效率高的特点,并且可以准确地描述开孔矩形屏蔽体内部各个分区的场分布规律,为进一步研究开孔腔体内部场分布提供了新思路。
附图说明
图1是带孔矩形屏蔽体示意图。
图2-5开孔屏蔽体内部磁场屏蔽效能随观测点位置变化情况示意图。
图6-9波导区的磁场屏蔽效能的CST计算结果的示意图/
具体实施方式
本发明提出一种计算开孔矩形屏蔽体内部场分布的分区解析模型,具体计算平面波垂直于小孔表面入射时,开孔中心轴线上的屏蔽体内部的场分布,将开孔屏蔽体内部的区域划分为小孔区、偶极区、波导区,在不同的区域建立场分布计算模型;以该区域的解析模型对其场分布进行计算;按惯例,
腔内任一点处的电场屏蔽效能定义为
SE=20log(E0/E)
式中,E0为入射平面波沿x轴方向的电场强度振幅,E为腔内任一点沿x轴方向的电场强度振幅;
腔内任一点处的磁场屏蔽效能定义为:
SH=20log(H0/H)
式中,H0为入射平面波沿y轴方向的磁场强度振幅,H为腔内任一点沿y轴方向的磁场强度振幅。
下面结合附图,对优选实施例作详细说明。
实施例1,具体实施步骤
步骤1:根据开孔矩形屏蔽体内部到开孔面的距离,可将开孔矩形屏蔽体内部分为三个区域:其中1倍小孔直径内的空间为小孔区;倍小孔直径内的空间为偶极区;4倍小孔直径外的部分为波导区。
步骤2:本优选实施例所用的模型如图1所示,是一个高=120mm,宽=300mm,深度=300mm;开孔半径为r=10mm,直径为dc=20mm,材料为理想导体的金属腔体。施加的入射平面波垂直入射并沿x轴极化,且极化电场强度为E0=1V/m。观测点位于开孔中心轴线上,且观测点在轴线上均匀分布,分别计算在500MHz和1GHz下各个区域的场分布。
步骤3:用公式(1)、(2)计算小孔区的场分布,并将模型数据代入公式,可得在小孔区,开孔中心轴线上的电场强度和磁场强度的值。
步骤4:用公式(3)、(4)、(5)计算得磁偶极矩m,将m和模型数据带入(6)、(7),求得偶极区开孔中心轴线上的电场强度和磁场强度。
步骤5:将模型数据带入公式(15)、(16),可计算得到波导区开孔中心轴线上的电场强度和磁场强度。
步骤6:将以上得到的各个区域的电场强度和磁场强度,带入公式(17)、(18),计算得到腔体开孔中心轴线上的电场屏蔽效能和磁场屏蔽效能。
步骤7:根据步骤1所述,1倍小孔直径20mm内的空间为小孔区,倍小孔直径内的空间为偶极区,4倍小孔直径80mm外的部分为波导区,将得到的各区域电场屏蔽效能和磁场屏蔽效能以观测点到开孔中心的距离为横坐标,电场屏蔽效能或磁场屏蔽效能为纵坐标,将各区域电场屏蔽效能和磁场屏蔽效能画在同一图上,便于比较。
由于波导区的计算公式中的C1、C2是由偶极区的计算公式(6)、(7)与波导区的计算公式(8)、(9)联立求得,其中C1、C2的大小随观测点距开孔的距离而变化,而能联立求解的条件是z在偶极区与波导区重叠的区域内取值,此时令z=z1,故z1的取值范围为
实施例2:
1GHz和500MHz时,开孔屏蔽体内部磁场屏蔽效能随观测点位置变化的情况如图2、3所示:
实施例3
1GHz和500MHz时,开孔屏蔽体内部电场屏蔽效能随观测点位置变化的情况如图4、5所示:
实施例4
如图6、7所示,1GHz和500MHz时,波导区的磁场屏蔽效能的计算公式中的r取55mm、76mm、85mm时的结果与CST计算结果的示意图。
实施例5
如图8、9所示,1GHz和500MHz时,波导区的电场屏蔽效能的计算公式中的r取55mm、76mm、85mm时的结果与CST计算结果的示意图。
由图可以看出,利用由z1之间取值得到的波导区计算公式,计算得到的结果与用CST软件计算得到的波导区的电磁屏蔽效能结果一致的比较好,因此z1之间取值得到的波导区计算公式能较为准确的反映波导区电磁场的变化规律。
综上所述,将开孔屏蔽体内部分为小孔区、偶极区、波导区三个区域,在不同区域建立其场分布解析模型,以该解析模型计算得到的结果与基于全波仿真软件CST的频域求解器得到的结果一致性较好,能准确反映腔体内部电磁场的变化规律,因此该解析模型是合理可行的。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (1)

1.一种计算开孔矩形屏蔽体内部场分布的分区解析模型,其特征在于,计算平面波垂直于小孔表面入射时,开孔中心轴线上的屏蔽体内部的场分布,将开孔矩形屏蔽体内部空间分解为三个区域:围绕开孔附近的小孔区、靠近屏蔽体背板的波导区、以及小孔区和波导区之间的偶极区,各个区域的场分布模型分别为:小孔区采用带开孔的无限大金属平板模型替代、偶极区采用偶极辐射模型近似、波导区采用波导模式场分布等效。对开孔屏蔽体内部进行分区,各个区域建立场分布计算模型,主要包括:
1)小孔区:在开孔矩形屏蔽体内部,距小孔面1倍直径内的场分布可以等效成无限大金属板附近的场分布,故运用计算无限大金属板开孔附近场分布的方法建立的计算模型如下:
磁场:在平面波垂直于小孔表面入射的条件下,推导出在小孔区,开孔中心轴线上的磁场强度为:
<mrow> <mi>H</mi> <mo>=</mo> <mfrac> <msub> <mi>H</mi> <mn>1</mn> </msub> <mi>&amp;pi;</mi> </mfrac> <mrow> <mo>(</mo> <msup> <mi>tan</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mfrac> <mi>R</mi> <mi>z</mi> </mfrac> <mo>-</mo> <mfrac> <mrow> <mi>R</mi> <mi>z</mi> </mrow> <mrow> <msup> <mi>R</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>z</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> <msub> <mi>e</mi> <mi>y</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
其中H1为当无限大金属板上没有开孔时金属板入射波侧的磁场切向分量,数值上等于入射平面波磁场沿y轴分量H0的两倍,R为小孔的半径,z为观测点到小孔中心的距离;ey指的是y轴的单位矢量
电场:在平面波的垂直入射下,推导出在小孔区,开孔中心轴线上的电场强度为:
<mrow> <mi>E</mi> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <msub> <mi>E</mi> <mi>o</mi> </msub> <mo>&amp;CenterDot;</mo> <mn>16</mn> <mi>R</mi> </mrow> <mrow> <mn>3</mn> <mi>&amp;lambda;</mi> <mo>&amp;lsqb;</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mfrac> <mn>1</mn> <mrow> <mi>j</mi> <mi>k</mi> <mi>R</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>cos</mi> <mo>(</mo> <mrow> <mi>k</mi> <mi>R</mi> </mrow> <mo>)</mo> <mo>+</mo> <mi>j</mi> <mi>sin</mi> <mo>(</mo> <mrow> <mi>k</mi> <mi>R</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mi>k</mi> <mi>z</mi> </mrow> </msup> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mfrac> <msup> <mi>z</mi> <mn>2</mn> </msup> <msup> <mi>&amp;rho;</mi> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <msup> <mi>R</mi> <mn>2</mn> </msup> <mrow> <msup> <mi>jk&amp;rho;</mi> <mn>3</mn> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mi>k</mi> <mi>&amp;rho;</mi> </mrow> </msup> <mo>&amp;rsqb;</mo> <msub> <mi>e</mi> <mi>x</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
其中E0为入射平面波沿x轴分量的电场强度,λ为波长,ρ2=R2+z2,k为波数,ex指的是x轴的单位矢量、ejkρ指的是e的jkρ次方,j为虚数单位,ejkz指的是e的jkz次方;
2)偶极区:开孔矩形屏蔽体内部,距小孔面倍小孔直径的范围内的电磁场,依据Bethe小孔耦合理论,并结合上述模型及入射平面波的性质,用位于开孔中心的磁偶极矩m的激发场来描述,m则由入射平面波决定,
m=-αmyH0,yey (3)
式中,H0,y为入射磁场H0在开孔面切向的y分量,极化系数αm与开孔的形状、尺寸有关,对于直径为dc的圆形孔的极化系数表达式为
<mrow> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>m</mi> <mi>y</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>d</mi> <mi>c</mi> <mn>3</mn> </msubsup> <mo>/</mo> <mn>3</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
由此得出磁偶极矩:
<mrow> <mi>m</mi> <mo>=</mo> <msubsup> <mi>&amp;eta;</mi> <mn>0</mn> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <msub> <mi>E</mi> <mn>0</mn> </msub> <msub> <mi>e</mi> <mi>y</mi> </msub> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>m</mi> <mi>y</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
式中,E0为入射电场沿x轴分量的振幅,η0为自由空间波阻抗,η0=377。可以得到在偶极区,开孔中心轴线上的磁场强度为:
<mrow> <mi>H</mi> <mo>=</mo> <mfrac> <mrow> <msup> <mi>k</mi> <mn>3</mn> </msup> <mi>S</mi> <mi>I</mi> </mrow> <mrow> <mn>4</mn> <mi>&amp;pi;</mi> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mi>k</mi> <mi>z</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mi>j</mi> <msup> <mrow> <mo>(</mo> <mi>k</mi> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msup> <mrow> <mo>(</mo> <mi>k</mi> <mi>z</mi> <mo>)</mo> </mrow> <mn>3</mn> </msup> </mfrac> <mo>&amp;rsqb;</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mi>k</mi> <mi>z</mi> </mrow> </msup> <msub> <mi>e</mi> <mi>y</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
电场强度为:
<mrow> <mi>E</mi> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <msup> <mi>jk</mi> <mn>3</mn> </msup> <mi>S</mi> <mi>I</mi> </mrow> <mrow> <mn>4</mn> <mi>&amp;pi;</mi> </mrow> </mfrac> <msub> <mi>&amp;eta;</mi> <mn>0</mn> </msub> <mo>&amp;lsqb;</mo> <mfrac> <mi>j</mi> <mrow> <mi>k</mi> <mi>z</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msup> <mrow> <mo>(</mo> <mi>k</mi> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mfrac> <mo>&amp;rsqb;</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mi>k</mi> <mi>z</mi> </mrow> </msup> <msub> <mi>e</mi> <mi>x</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> 1
式中,SI=2m,为波数,z为观测点到孔的距离;
3)波导区:在开孔矩形屏蔽体内部,距小孔面4倍直径之后的范围内,根据传输线波导理论,将开孔金属腔体看成到开孔金属板距离为d的金属板短路的传输线,由此可得开孔面中心轴线上的电场强度和磁场强度分别为:
E=C1sin kg(d-z)ex (8)
H=C2cos kg(d-z)ey (9)
由于偶极区与波导区有一部分区域是重叠的,这部分区域用偶极区计算方法算出的电磁场强度与波导区计算方法算出的相同,因此利用这一点,当z在距小孔倍直径内取值时,令z=z1,联立得到
<mrow> <msub> <mi>C</mi> <mn>1</mn> </msub> <mi>sin</mi> <mi> </mi> <msub> <mi>k</mi> <mi>g</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mo>-</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <msup> <mi>jk</mi> <mn>3</mn> </msup> <mi>S</mi> <mi>I</mi> </mrow> <mrow> <mn>4</mn> <mi>&amp;pi;</mi> </mrow> </mfrac> <msub> <mi>&amp;eta;</mi> <mn>0</mn> </msub> <mo>&amp;lsqb;</mo> <mfrac> <mi>j</mi> <mrow> <msub> <mi>kz</mi> <mn>1</mn> </msub> </mrow> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msup> <mrow> <mo>(</mo> <msub> <mi>kz</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mfrac> <mo>&amp;rsqb;</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>jkz</mi> <mn>1</mn> </msub> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>C</mi> <mn>2</mn> </msub> <mi>cos</mi> <mi> </mi> <msub> <mi>k</mi> <mi>g</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mo>-</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msup> <mi>k</mi> <mn>3</mn> </msup> <mi>S</mi> <mi>I</mi> </mrow> <mrow> <mn>4</mn> <mi>&amp;pi;</mi> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <msub> <mi>kz</mi> <mn>1</mn> </msub> </mrow> </mfrac> <mo>+</mo> <mfrac> <mi>j</mi> <msup> <mrow> <mo>(</mo> <msub> <mi>kz</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msup> <mrow> <mo>(</mo> <msub> <mi>kz</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>3</mn> </msup> </mfrac> <mo>&amp;rsqb;</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>jkz</mi> <mn>1</mn> </msub> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
解得
<mrow> <msub> <mi>C</mi> <mn>1</mn> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mi>j</mi> <mn>2</mn> <msup> <mi>&amp;pi;</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>0</mn> </msub> <msubsup> <mi>d</mi> <mi>C</mi> <mn>3</mn> </msubsup> </mrow> <mrow> <mn>3</mn> <msup> <mi>&amp;lambda;</mi> <mn>3</mn> </msup> <mi>sin</mi> <mo>&amp;lsqb;</mo> <msub> <mi>k</mi> <mi>g</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mo>-</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <mfrac> <mrow> <mi>j</mi> <mi>&amp;lambda;</mi> </mrow> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;z</mi> <mn>1</mn> </msub> </mrow> </mfrac> <mo>+</mo> <mfrac> <msup> <mi>&amp;lambda;</mi> <mn>2</mn> </msup> <mrow> <mn>4</mn> <msup> <mi>&amp;pi;</mi> <mn>2</mn> </msup> <msup> <msub> <mi>z</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> </mrow> </mfrac> <mo>&amp;rsqb;</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;z</mi> <mn>1</mn> </msub> </mrow> <mi>&amp;lambda;</mi> </mfrac> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>C</mi> <mn>2</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <msup> <mi>&amp;pi;</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>0</mn> </msub> <msubsup> <mi>d</mi> <mi>C</mi> <mn>3</mn> </msubsup> </mrow> <mrow> <mn>3</mn> <msub> <mi>&amp;eta;</mi> <mn>0</mn> </msub> <msup> <mi>&amp;lambda;</mi> <mn>3</mn> </msup> <mi>cos</mi> <mo>&amp;lsqb;</mo> <msub> <mi>k</mi> <mi>g</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mo>-</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <mo>-</mo> <mfrac> <mi>&amp;lambda;</mi> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;z</mi> <mn>1</mn> </msub> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mi>j&amp;lambda;</mi> <mn>2</mn> </msup> </mrow> <mrow> <mn>4</mn> <msup> <mi>&amp;pi;</mi> <mn>2</mn> </msup> <msup> <msub> <mi>z</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <mfrac> <msup> <mi>&amp;lambda;</mi> <mn>3</mn> </msup> <mrow> <mn>8</mn> <msup> <mi>&amp;pi;</mi> <mn>3</mn> </msup> <msup> <msub> <mi>z</mi> <mn>1</mn> </msub> <mn>3</mn> </msup> </mrow> </mfrac> <mo>&amp;rsqb;</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;z</mi> <mn>1</mn> </msub> </mrow> <mi>&amp;lambda;</mi> </mfrac> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>k</mi> <mi>g</mi> </msub> <mo>=</mo> <mi>k</mi> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mi>&amp;lambda;</mi> <mrow> <mn>2</mn> <mi>a</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
综上,在波导区,开孔面中心轴线上的电场强度为
<mrow> <mi>E</mi> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mi>j</mi> <mn>2</mn> <msup> <mi>&amp;pi;</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>0</mn> </msub> <msubsup> <mi>d</mi> <mi>C</mi> <mn>3</mn> </msubsup> </mrow> <mrow> <mn>3</mn> <msup> <mi>&amp;lambda;</mi> <mn>3</mn> </msup> <mi>sin</mi> <mo>&amp;lsqb;</mo> <msub> <mi>k</mi> <mi>g</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mo>-</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <mfrac> <mrow> <mi>j</mi> <mi>&amp;lambda;</mi> </mrow> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;z</mi> <mn>1</mn> </msub> </mrow> </mfrac> <mo>+</mo> <mfrac> <msup> <mi>&amp;lambda;</mi> <mn>2</mn> </msup> <mrow> <mn>4</mn> <msup> <mi>&amp;pi;</mi> <mn>2</mn> </msup> <msup> <msub> <mi>z</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> </mrow> </mfrac> <mo>&amp;rsqb;</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;z</mi> <mn>1</mn> </msub> </mrow> <mi>&amp;lambda;</mi> </mfrac> </mrow> </msup> <mi>sin</mi> <mi> </mi> <msub> <mi>k</mi> <mi>g</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mo>-</mo> <mi>z</mi> <mo>)</mo> </mrow> <msub> <mi>e</mi> <mi>x</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
磁场强度为
<mrow> <mi>H</mi> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <msup> <mi>&amp;pi;</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>0</mn> </msub> <msubsup> <mi>d</mi> <mi>C</mi> <mn>3</mn> </msubsup> </mrow> <mrow> <mn>3</mn> <msub> <mi>&amp;eta;</mi> <mn>0</mn> </msub> <msup> <mi>&amp;lambda;</mi> <mn>3</mn> </msup> <mi>cos</mi> <mo>&amp;lsqb;</mo> <msub> <mi>k</mi> <mi>g</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mo>-</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <mo>-</mo> <mfrac> <mi>&amp;lambda;</mi> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;z</mi> <mn>1</mn> </msub> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mi>j&amp;lambda;</mi> <mn>2</mn> </msup> </mrow> <mrow> <mn>4</mn> <msup> <mi>&amp;pi;</mi> <mn>2</mn> </msup> <msup> <msub> <mi>z</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <mfrac> <msup> <mi>&amp;lambda;</mi> <mn>3</mn> </msup> <mrow> <mn>8</mn> <msup> <mi>&amp;pi;</mi> <mn>3</mn> </msup> <msup> <msub> <mi>z</mi> <mn>1</mn> </msub> <mn>3</mn> </msup> </mrow> </mfrac> <mo>&amp;rsqb;</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;z</mi> <mn>1</mn> </msub> </mrow> <mi>&amp;lambda;</mi> </mfrac> </mrow> </msup> <mi>cos</mi> <mi> </mi> <msub> <mi>k</mi> <mi>g</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mo>-</mo> <mi>z</mi> <mo>)</mo> </mrow> <msub> <mi>e</mi> <mi>y</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
其中,E0为入射平面波沿x轴分量的电场强度振幅,dc为圆形开孔的直径,λ为波长,z为观测点到孔的距离,z1在距小孔3~4倍直径内取值,a为开孔屏蔽体的高度,kg为传播常数,为波数,η0为自由空间波阻抗,η0=377。C1、C2为计算过程设定的两个未知数;
至此,已给出了计算开孔矩形屏蔽体内部场分布的全套方程,按惯例,腔内任一点处的电场屏蔽效能定义为
SE=20log(E0/E) (17)
式中,E0为入射平面波沿x轴方向的电场强度振幅,E为腔内任一点沿x轴方向的电场强度振幅;
腔内任一点处的磁场屏蔽效能定义为:
SH=20log(H0/H) (18)
式中,H0为入射平面波沿y轴方向的磁场强度振幅,H为腔内任一点沿y轴方向的磁场强度振幅。
CN201710530340.9A 2017-06-27 2017-06-27 一种计算开孔矩形屏蔽体内部场分布的分区解析模型的方法 Expired - Fee Related CN107255755B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710530340.9A CN107255755B (zh) 2017-06-27 2017-06-27 一种计算开孔矩形屏蔽体内部场分布的分区解析模型的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710530340.9A CN107255755B (zh) 2017-06-27 2017-06-27 一种计算开孔矩形屏蔽体内部场分布的分区解析模型的方法

Publications (2)

Publication Number Publication Date
CN107255755A true CN107255755A (zh) 2017-10-17
CN107255755B CN107255755B (zh) 2019-09-20

Family

ID=60024807

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710530340.9A Expired - Fee Related CN107255755B (zh) 2017-06-27 2017-06-27 一种计算开孔矩形屏蔽体内部场分布的分区解析模型的方法

Country Status (1)

Country Link
CN (1) CN107255755B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108536929A (zh) * 2018-03-23 2018-09-14 电子科技大学 一种应用arpack求解波导结构色散特性的方法
CN109581096A (zh) * 2018-10-18 2019-04-05 上海无线电设备研究所 屏蔽室下金属屏蔽材料屏蔽效能分析模型及其测量方法
CN110058183A (zh) * 2019-05-16 2019-07-26 华北电力大学 一种小孔磁极化系数的提取方法
CN110489866A (zh) * 2019-08-20 2019-11-22 全球能源互联网研究院有限公司 一种计算封闭腔体磁场屏蔽效能的预测方法
CN110489808A (zh) * 2019-07-23 2019-11-22 国网内蒙古东部电力有限公司电力科学研究院 一种提取金属网表面阻抗的方法及系统
CN111046524A (zh) * 2019-11-12 2020-04-21 天津博迈科海洋工程有限公司 一种海洋电气间散热优化方法
CN111259521A (zh) * 2020-01-09 2020-06-09 华北电力大学 一种矩形金属腔体屏蔽盒上开孔等效阻抗提取方法及系统
CN111965458A (zh) * 2020-08-20 2020-11-20 国网经济技术研究院有限公司 一种平面状周期开孔屏蔽体的屏蔽效能确定方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103743985A (zh) * 2014-01-24 2014-04-23 中国人民解放军军械工程学院 一种测试复杂电磁环境下织物材料屏蔽效能的方法
CN104915499A (zh) * 2015-06-10 2015-09-16 电子科技大学 一种预测开孔腔体谐振点频率的快速算法
CN105183939A (zh) * 2015-07-21 2015-12-23 电子科技大学 基于等效电路的带孔阵腔体屏蔽效能计算方法
CN105259450A (zh) * 2015-11-05 2016-01-20 北京东方计量测试研究所 一种磁屏蔽效能估算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103743985A (zh) * 2014-01-24 2014-04-23 中国人民解放军军械工程学院 一种测试复杂电磁环境下织物材料屏蔽效能的方法
CN104915499A (zh) * 2015-06-10 2015-09-16 电子科技大学 一种预测开孔腔体谐振点频率的快速算法
CN105183939A (zh) * 2015-07-21 2015-12-23 电子科技大学 基于等效电路的带孔阵腔体屏蔽效能计算方法
CN105259450A (zh) * 2015-11-05 2016-01-20 北京东方计量测试研究所 一种磁屏蔽效能估算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
程二威等: "屏蔽体内部电磁场分布特性分析", 《河北科技大学学报》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108536929A (zh) * 2018-03-23 2018-09-14 电子科技大学 一种应用arpack求解波导结构色散特性的方法
CN109581096A (zh) * 2018-10-18 2019-04-05 上海无线电设备研究所 屏蔽室下金属屏蔽材料屏蔽效能分析模型及其测量方法
CN110058183A (zh) * 2019-05-16 2019-07-26 华北电力大学 一种小孔磁极化系数的提取方法
CN110058183B (zh) * 2019-05-16 2020-06-19 华北电力大学 一种小孔磁化系数的提取方法
CN110489808A (zh) * 2019-07-23 2019-11-22 国网内蒙古东部电力有限公司电力科学研究院 一种提取金属网表面阻抗的方法及系统
CN110489808B (zh) * 2019-07-23 2023-05-12 国网内蒙古东部电力有限公司电力科学研究院 一种提取金属网表面阻抗的方法及系统
CN110489866A (zh) * 2019-08-20 2019-11-22 全球能源互联网研究院有限公司 一种计算封闭腔体磁场屏蔽效能的预测方法
CN111046524A (zh) * 2019-11-12 2020-04-21 天津博迈科海洋工程有限公司 一种海洋电气间散热优化方法
CN111259521A (zh) * 2020-01-09 2020-06-09 华北电力大学 一种矩形金属腔体屏蔽盒上开孔等效阻抗提取方法及系统
CN111259521B (zh) * 2020-01-09 2022-03-04 华北电力大学 一种矩形金属腔体屏蔽盒上开孔等效阻抗提取方法及系统
CN111965458A (zh) * 2020-08-20 2020-11-20 国网经济技术研究院有限公司 一种平面状周期开孔屏蔽体的屏蔽效能确定方法及系统
CN111965458B (zh) * 2020-08-20 2023-05-23 国网经济技术研究院有限公司 一种平面状周期开孔屏蔽体的屏蔽效能确定方法及系统

Also Published As

Publication number Publication date
CN107255755B (zh) 2019-09-20

Similar Documents

Publication Publication Date Title
CN107255755A (zh) 一种计算开孔矩形屏蔽体内部场分布的分区解析模型
Klinkenbusch On the shielding effectiveness of enclosures
Nikita et al. A study of uncertainties in modeling antenna performance and power absorption in the head of a cellular phone user
Azadifar et al. Partial discharge localization using electromagnetic time reversal: A performance analysis
Kong et al. A novel dielectric conformal FDTD method for computing SAR distribution of the human body in a metallic cabin illuminated by an intentional electromagnetic pulse (IEMP)
Jiao et al. Low-frequency magnetic shielding of planar shields: A unified wave impedance formula for the transmission line analogy
Xiong et al. A simple local approximation FDTD model of short apertures with a finite thickness
Štumpf et al. Pulsed electromagnetic field signal transfer across a thin magneto-dielectric sheet
Bayati et al. Analyzing the near and far field using finite difference and finite element method
Park et al. Analysis of quasistatic magnetic field penetration into multiple slits in a conducting plane loaded with a ferrite sheet
Belokour et al. A 2D-transmission line model for the EM field estimation inside enclosures with apertures
Fan et al. FDTD and PSTD simulations for plasma applications
Milutinović et al. Circuital and Numerical Models for Calculation of Shielding Effectiveness of Enclosure with Apertures and Monitoring Dipole Antenna Inside
Yang et al. Aperture coupling analysis on shielding enclosure excited by lightning electromagnetic pulse
SU et al. An analytical method for calculating the shielding effectiveness of an enclosure with an oblique rectangular aperture
Nicolas et al. Specific absorption rate assessments based on a selective isotropic measuring system for electromagnetic fields
Li et al. Analytical Formulation for Electromagnetic Leakage from an Apertured Rectangular Cavity.
Zhang et al. Specific absorption rate assessment of fifth generation mobile phones with specific anthropomorphic mannequin model and high‐resolution anatomical head model
Wei et al. Efficient EMC simulation of enclosures with apertures residing in an electrically large platform using the MM-UTD method
Park et al. Electromagnetic shielding analysis of multiple slits on a metal plate coated with a ferrite sheet
Bogdanov et al. Comparison of MoM implementations of surface impedance boundary conditions for material sheets with analytical and measurement data
Panagamuwa et al. A study of the validation of RF energy Specific Absorption Rates for simulations of anatomically correct head FDTD simulations and truncated DASY4 standard equipment measurements
Newman et al. Moment method analysis of the electric shielding factor of a conducting TM shield at ELF
Smitha et al. Increase in Specific Absorption Rate (SAR) Of Mobile Phone-A Threat to Human
Sun et al. Approximate simulation of low-frequency magnetic shielding of a rectangular shielding box with all walls perforated periodical holes

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190920