CN107273627B - 一种快速解析电磁场强度分布的开孔电磁屏蔽体设计方法 - Google Patents
一种快速解析电磁场强度分布的开孔电磁屏蔽体设计方法 Download PDFInfo
- Publication number
- CN107273627B CN107273627B CN201710493738.XA CN201710493738A CN107273627B CN 107273627 B CN107273627 B CN 107273627B CN 201710493738 A CN201710493738 A CN 201710493738A CN 107273627 B CN107273627 B CN 107273627B
- Authority
- CN
- China
- Prior art keywords
- shield
- electric field
- hole
- field intensity
- shielding
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Shielding Devices Or Components To Electric Or Magnetic Fields (AREA)
Abstract
本发明公开了属于电磁屏蔽技术领域的一种快速解析电磁场强度分布的开孔电磁屏蔽体设计方法。利用屏蔽体开孔与无限大开孔平板开孔上的电位分布的近似等效,通过分离变量法得到开孔矩形屏蔽体中的电位与电场强度分布,进而得到开孔电磁屏蔽体的屏蔽效能,并基于此结果指导屏蔽体开孔的大小、形状及位置。本发明提供的方法是具有计算速度较快又不需借助电磁场数值计算软件的解析方法。并且可以作为如何设计开孔屏蔽体的依据。具有重要的实际意义。
Description
技术领域
本发明属于电磁屏蔽技术领域,特别涉及一种快速解析电磁场强度分布的开孔电磁屏蔽体设计方法。
背景技术
随着电网建设的发展,电网输电电压等级不断提高,智能输变电技术不断推广,电磁骚扰问题日益凸显。另一方面,在工业化、信息化发展迅速的当今时代,各类电子产品和电力设备为生活提供了巨大的便利,也带来了各种电磁辐射,不仅影响着人们的正常生活,也可能威胁到身体健康。因此,如何抑制电磁骚扰,是现今备受关注的一个问题。
依靠金属腔体对电磁骚扰进行屏蔽,是常用的技术手段之一。在实际应用中,小到各种电力电子设备,大到各类试验场所等,都常以金属腔体作为屏蔽,隔绝外界的电磁骚扰。一般地,除了对低频磁场的屏蔽效果不太明显,金属腔体可以屏蔽多数电磁骚扰。但事实上,大多数屏蔽体都会出于通讯、通风、散热等需求而带有若干开孔。开孔的存在为屏蔽体内外的电磁场提供了耦合的途径,从而大大削弱了屏蔽体的屏蔽性能。所以,对于开孔屏蔽体的研究具有重要的工程和社会意义。
目前,对于开孔屏蔽体的研究中,电磁波的屏蔽研究较多,但在实际中也会碰到低频电场屏蔽的问题,如变电站、高压设备附近的屏蔽问题。根据电磁环境控制限值标准GB8702-2014,为控制电场所致的公众暴露,环境中电场强度应小于一定的限值,当频率在1Hz25Hz时,电场强度应不大于8000V/m,在常用工频50Hz时,电场强度应不大于4000V/m。因此以静电场(低频电场)为研究对象做出研究,可以为屏蔽体的设计提供指导。
屏蔽效能(SE),是衡量屏蔽体有效性的重要参数。是如何设计开孔屏蔽体的依据;对于开孔屏蔽体屏蔽效能的研究,目前多采用有限元法等数值方法进行仿真计算,但数值方法需要借助专门的计算软件,并建模后进行,不但对工具的要求较高,而且往往计算时间较长。因此,提出一种快速解析开孔屏蔽体的电磁场强度分布的方法,具有重要的实际意义。本发明的目的在于,提供一种用于开孔矩形屏蔽体电场屏蔽效能研究的计算速度较快又不需借助电磁场数值计算软件的解析方法。并且可以作为如何设计开孔屏蔽体的依据。
发明内容
本发明的目的是提出一种快速解析电磁场强度分布的开孔电磁屏蔽体设计方法,其特征在于,利用屏蔽体开孔与无限大开孔平板开孔上的电位分布的近似等效,通过分离变量法得到开孔矩形屏蔽体中的电位与电场强度分布,进而得到开孔电磁屏蔽体的屏蔽效能,并基于此结果指导屏蔽体开孔的大小、形状及位置;具体包括如下步骤:
步骤1:用无限大平板上开孔的电位分布近似等效屏蔽体上开孔的电位分布;
步骤2:通过分离变量法计算得到屏蔽体内开孔的电位分布;
步骤3:对步骤2得到的电位分布求梯度,得到屏蔽体内各个位置的电场强度沿各坐标轴的分量,将屏蔽体内的电位分布转换为电场强度分布;
步骤4:利用屏蔽体内电场强度分布求得屏蔽体内不同位置的电场屏蔽的效能。
步骤5:利用步骤4求得屏蔽体内不同位置的电场屏蔽的效能,并基于此结果指导屏蔽体开孔的大小、形状及位置。
所述步骤1中,利用无限大开孔平板对开孔屏蔽体进行等效时还包括以下步骤:
步骤101:带孔平板所带小孔为圆形或椭圆形时,根据小孔上的电位分布的解析公式,直接通过等效替代得到带圆形及椭圆形孔的屏蔽体开孔上的电位分布;
步骤102:对于带矩形孔的无限大平板,采用椭圆孔进行等效替代;要求等效的椭圆孔与原矩形孔面积相同,且长短轴之比与原矩形孔长宽之比相同;在无限大开孔平板中,椭圆等效替代后得到的结果与矩形孔结果相近。因此,步骤102中屏蔽体带矩形孔时,孔上的电位可以等效为带对应椭圆开孔的无限大平板开孔上的电位分布;因此,所述等效屏蔽体上待计算的小孔仅限圆孔、椭圆孔及矩形孔三种。
所述步骤1,首先将带孔矩形屏蔽体的开孔上的电位分布等效为无限大带孔平板的开孔上的电位分布;因此,当屏蔽体处于一个外加的垂直于孔面的电场E0中时,对于圆形开孔,孔上的电位分布为式(1),
其中,xp和yp是圆孔圆心的x坐标和y坐标,R为圆孔半径;
对于椭圆孔,孔上的电位分布为式(2),
其中,ta和tb是椭圆孔的半长轴和半短轴,e为椭圆离心率,E(e)为第二类完全椭圆积分。
所述步骤2,屏蔽体内的电位表示为式(3),
其中φ(x,y)即为公式1和2中小孔上的电位分布;
所述步骤3,将公式(3)对坐标x、y、z分别求导,得到三个方向上的电场强度分量大小,分别为
则屏蔽体内电场强度大小的分布
所述步骤4,对于电场,屏蔽效能定义
其中,E0为没有屏蔽时空间某个位置的电场强度或外部的垂直电场强度;
E为有屏蔽时该位置的电场强度;利用式(10)即可得到屏蔽体内的屏蔽效能。
本发明的有益效果是本发明采用的计算方法是具有计算速度较快又不需借助电磁场数值计算软件的解析方法。并且可以作为如何设计屏蔽体开孔大小、形状及位置的依据。具有重要的实际意义。
附图说明
图1是带孔矩形屏蔽体示意图。
具体实施方式
本发明提出一种快速解析电磁场强度分布的开孔电磁屏蔽体设计方法,利用屏蔽体开孔与无限大开孔平板开孔上的电位分布的近似等效,通过分离变量法得到开孔矩形屏蔽体中的电位与电场强度分布,进而得到开孔电磁屏蔽体的屏蔽效能,并基于此结果指导屏蔽体开孔的大小、形状及位置。下面结合附图1和具体实施例对本发明作详细说明。
如图1所示的一面带孔的矩形屏蔽体(以圆孔为例,可替换为椭圆孔和矩形孔),采用本发明所述方法进行计算。具体执行步骤如下:
步骤1,首先将带孔矩形屏蔽体的开孔上的电位分布等效为无限大带孔平板的开孔上的电位分布。因此,当屏蔽体处于一个外加的垂直于孔面的电场E0中时,对于圆形开孔,孔上的电位分布为式(1),
其中,xp和yp是圆孔圆心的x坐标和y坐标,R为圆孔半径。
对于椭圆孔,孔上的电位分布为式(2)
其中,ta和tb是椭圆孔的半长轴和半短轴,e为椭圆离心率,E(e)为第二类完全椭圆积分。
步骤2,屏蔽体内的电位可表示为式(3)
其中φ(x,y)即为式(1)和式(2)中小孔上的电位分布。
步骤3,将式(3)对坐标x、y、z分别求导,得到三个方向上的电场强度分量大小,分别为
则屏蔽体内电场强度大小的分布
其中,E0为没有屏蔽时空间某个位置的电场强度,E为有屏蔽时该位置的电场强度。本发明中,E0为附图1中外部的垂直电场强度,是一固定值,且大小已知;E为式(9)中得到的屏蔽体内不同位置的电场强度大小。利用式(10)即可得到屏蔽体内的屏蔽效能。
实施例1
如图1所示,一个一面开圆孔的矩形屏蔽体的屏蔽效能的计算。
令屏蔽体每边长为1m(图1中a=b=c=1m),在上表面的中心(0.5,0.5,1)处开一圆孔,孔半径为0.01m,小孔外部垂直于小孔的均匀电场场强E0大小为2V/m。
按步骤13进行运算,得到屏蔽体内部的电场强度大小。特别地,在步骤2中,利用公式3求解电位分布时,m和n的取值上限无法也无需取至∞,但合理的m、n取值可以使开孔附近的电场强度计算结果更为准确。在实施例1中,令m=n=80。
为验证结果的准确性,利用ANSYS软件对同样模型进行数值仿真,对比结果如下。
然后,为考察屏蔽体内除开孔中轴线上以外的位置的计算结果,设置了三条水平的观测线,分别距小孔0.2m、0.5m、0.8m,对应附图1中的三条线段。提取这三条观测线上场强的计算结果,与ANSYS仿真结果对比于附图5~7中。
最后,进行步骤4,通过计算得到的电场强度以及外部场强E0,可以得到屏蔽体内不同位置的电场屏蔽效能。以开孔中轴线上的屏蔽效能为例,曲线如附图8所示。
实施例2
本发明所述步骤中,对矩形孔的计算采用了椭圆孔近似等效,为验证该近似的准确性,实施例2将实施例1中的圆孔改为矩形孔进行计算。屏蔽体尺寸和位置不变,开孔中心仍为(0.5,0.5,1),矩形孔长0.04m,宽0.02m。外部场强E0不变。
与实施例1相同,先通过步骤1~3得到屏蔽体内的电场强度分布。特别地,在步骤2中,取m=n=55。
虽然采用本发明的方法进行计算时,需要用椭圆形孔近似替代矩形孔,但ANSYS软件进行数值计算时不受该限制,因此,用ANSYS对实施例2中给出的带矩形孔屏蔽体模型进行计算,与本发明中方法做出对比。
开孔中轴线上的电场强度分段对比于图911中,距孔不同距离的三条观测线(观测线同实施例1)上的电场强度对比于图1214中。然后通过步骤4,得到屏蔽体的屏蔽效能。同样地,以开孔中轴线上的屏蔽效能为例,曲线如图15所示。
Claims (3)
1.一种快速解析电磁场强度分布的开孔电磁屏蔽体设计方法,其特征在于,利用屏蔽体开孔与无限大开孔平板开孔上的电位分布的近似等效,通过分离变量法得到开孔矩形屏蔽体中的电位与电场强度分布,进而得到其开孔电磁屏蔽体的屏蔽效能;进而分析屏蔽体开孔的大小、数量,形状及位置;具体包括如下步骤:
步骤1:用无限大平板上开孔的电位分布近似等效屏蔽体上开孔的电位分布;包括以下步骤:
步骤101:带孔平板所带小孔为圆形或椭圆形时,根据小孔上的电位分布的解析公式,直接通过等效替代得到带圆形及椭圆形孔的屏蔽体开孔上的电位分布;
步骤102:对于带矩形孔的无限大平板,采用椭圆孔进行等效替代;要求等效的椭圆孔与原矩形孔面积相同,且长短轴之比与原矩形孔长宽之比相同;在无限大开孔平板中,椭圆等效替代后得到的结果与矩形孔结果相近;因此,步骤102中屏蔽体带矩形孔时,孔上的电位能够等效为带对应椭圆开孔的无限大平板开孔上的电位分布;因此,所述等效屏蔽体上待计算的小孔仅限圆孔、椭圆孔及矩形孔三种;
步骤2:通过分离变量法计算得到屏蔽体内开孔的电位分布;首先将带孔矩形屏蔽体的开孔上的电位分布等效为无限大带孔平板的开孔上的电位分布;因此,当屏蔽体处于一个外加的垂直于孔面的电场E0中时,对于圆形开孔,孔上的电位分布为式(1),
其中,xp和yp是圆孔圆心的x坐标和y坐标,R为圆孔半径;
对于椭圆孔,孔上的电位分布为式(2),
其中,ta和tb是椭圆孔的半长轴和半短轴,e为椭圆离心率,E(e)为第二类完全椭圆积分;
步骤3:对步骤2得到的电位分布求梯度,得到屏蔽体内各个位置的电场强度沿各坐标轴的分量,将屏蔽体内的电位分布转换为电场强度分布;
步骤4:利用屏蔽体内电场强度分布求得屏蔽体内不同位置的电场屏蔽的效能;
步骤5:利用步骤4求得屏蔽体内不同位置的电场屏蔽的效能,并基于此结果指导屏蔽体开孔的大小、形状及位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710493738.XA CN107273627B (zh) | 2017-06-26 | 2017-06-26 | 一种快速解析电磁场强度分布的开孔电磁屏蔽体设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710493738.XA CN107273627B (zh) | 2017-06-26 | 2017-06-26 | 一种快速解析电磁场强度分布的开孔电磁屏蔽体设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107273627A CN107273627A (zh) | 2017-10-20 |
CN107273627B true CN107273627B (zh) | 2020-06-19 |
Family
ID=60068541
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710493738.XA Expired - Fee Related CN107273627B (zh) | 2017-06-26 | 2017-06-26 | 一种快速解析电磁场强度分布的开孔电磁屏蔽体设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107273627B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111122987B (zh) * | 2019-12-26 | 2020-12-04 | 华北电力大学 | 一种磁场屏蔽效能预测方法及系统 |
CN111125610B (zh) * | 2019-12-26 | 2022-02-15 | 华北电力大学 | 一种磁场屏蔽效能预测方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08298391A (ja) * | 1995-04-26 | 1996-11-12 | Nec Ibaraki Ltd | 電子機器の電磁シールド構造 |
CN103635068A (zh) * | 2012-08-27 | 2014-03-12 | 通用电气公司 | 能量转换系统的电磁屏蔽结构和相关方法 |
CN104915499A (zh) * | 2015-06-10 | 2015-09-16 | 电子科技大学 | 一种预测开孔腔体谐振点频率的快速算法 |
CN105911370A (zh) * | 2016-04-28 | 2016-08-31 | 中国科学院电子学研究所 | 一种金属微结构电场传感器 |
CN106199287A (zh) * | 2016-08-19 | 2016-12-07 | 华北电力大学 | 一种基于矩形波导的材料电场屏蔽效能测试系统和方法 |
-
2017
- 2017-06-26 CN CN201710493738.XA patent/CN107273627B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08298391A (ja) * | 1995-04-26 | 1996-11-12 | Nec Ibaraki Ltd | 電子機器の電磁シールド構造 |
CN103635068A (zh) * | 2012-08-27 | 2014-03-12 | 通用电气公司 | 能量转换系统的电磁屏蔽结构和相关方法 |
CN104915499A (zh) * | 2015-06-10 | 2015-09-16 | 电子科技大学 | 一种预测开孔腔体谐振点频率的快速算法 |
CN105911370A (zh) * | 2016-04-28 | 2016-08-31 | 中国科学院电子学研究所 | 一种金属微结构电场传感器 |
CN106199287A (zh) * | 2016-08-19 | 2016-12-07 | 华北电力大学 | 一种基于矩形波导的材料电场屏蔽效能测试系统和方法 |
Non-Patent Citations (3)
Title |
---|
Resonance suppression and electromagnetic shielding effectiveness improvement of an apertured rectangular cavity by using wall losses;Jiao ChongQing;《chin.Phys.B 》;20131231;第22卷(第8期);全文 * |
基于电磁屏蔽的箱体设计及有限元仿真的研究;金鹏;《中国优秀硕士学位论文全文数据库 信息科技辑》;20070915(第3期);全文 * |
开孔矩形腔体的近场电磁屏蔽效能研究;焦重庆 等;《物理学报》;20131231;第62卷(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107273627A (zh) | 2017-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ergul et al. | Linear-linear basis functions for MLFMA solutions of magnetic-field and combined-field integral equations | |
CN107273627B (zh) | 一种快速解析电磁场强度分布的开孔电磁屏蔽体设计方法 | |
CN102385655B (zh) | 一种电子设备辐射耦合电磁敏感性仿真方法 | |
CN107255755A (zh) | 一种计算开孔矩形屏蔽体内部场分布的分区解析模型 | |
Tatematsu | A technique for representing coaxial cables for FDTD-based surge simulations | |
CN105183939A (zh) | 基于等效电路的带孔阵腔体屏蔽效能计算方法 | |
CN104317984B (zh) | 基于分域建模的船舶电磁散射预测方法及系统 | |
Tong et al. | Analysis of electric field influence on buildings under high‐voltage transmission lines | |
CN105844019B (zh) | 一种基于车辆表面电流分布规律的接地点选择方法 | |
CN114417677B (zh) | 一种开关类设备操作强电磁骚扰特性分析方法 | |
CN108037389A (zh) | 一种基于等效场原理的线缆耦合分析方法 | |
CN107992696A (zh) | 一种复杂色散媒质中的改进的指数时间积分构造方法 | |
Andreotti et al. | A new solution for the evaluation of the horizontal electric fields from lightning in presence of a finitely conducting ground | |
Prácser et al. | Two approaches for modeling ELF wave propagation in the earth-ionosphere cavity with day-night asymmetry | |
Hajiaboli et al. | Analysis of grounding systems in the vicinity of hemispheroidal heterogeneities | |
CN112231880B (zh) | 一种超导磁体线圈设计方法 | |
CN115688516A (zh) | 一种无线通信设备电磁屏蔽罩性能优化方法及装置 | |
Lopez-Pena et al. | On the analytic-numeric treatment of weakly singular integrals on arbitrary polygonal domains | |
CN103617372A (zh) | 一种换流阀三维电场计算方法 | |
Jiang et al. | A multiresolution time-domain method for LEMP calculation and comparison with FDTD | |
Ab Ghani et al. | Effects of parameter adjustment on the electromagnetic field of an overhead power transmission line model | |
Bezrukovs et al. | The assessment of wind speed distortions in a simulated flow around a lattice cellular communication mast | |
Zhang et al. | Simulation calculation of electric field protection for live working on Ultra‐high voltage transmission line | |
Hao et al. | Analytical method for load response of a transmission line in a double-enclosure with multiple covered apertures | |
CN107330151B (zh) | 一种小孔电极化系数的提取方法 |
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 |
Granted publication date: 20200619 Termination date: 20210626 |
|
CF01 | Termination of patent right due to non-payment of annual fee |