CN105183939A - 基于等效电路的带孔阵腔体屏蔽效能计算方法 - Google Patents
基于等效电路的带孔阵腔体屏蔽效能计算方法 Download PDFInfo
- Publication number
- CN105183939A CN105183939A CN201510428920.8A CN201510428920A CN105183939A CN 105183939 A CN105183939 A CN 105183939A CN 201510428920 A CN201510428920 A CN 201510428920A CN 105183939 A CN105183939 A CN 105183939A
- Authority
- CN
- China
- Prior art keywords
- battle array
- hole
- equivalent
- cavity
- holes
- 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.)
- Pending
Links
Landscapes
- Shielding Devices Or Components To Electric Or Magnetic Fields (AREA)
Abstract
本发明公开了一种基于等效电路的带孔阵腔体屏蔽效能计算方法;其步骤包括构建等效单孔模型,构建等效电路模型,计算等效单孔阻抗并代替孔阵阻抗及计算带孔阵腔体屏蔽效能。本发明通过构建等效单孔模型,在保证仿真精度的前提下降低了建模复杂度,显著减少了网格数量,提高了仿真效率,简化了带孔阵复杂腔体的建模过程,大幅度的提高了等效电路方法在计算带孔阵腔体屏蔽效能时的精度,实现带孔阵腔体屏蔽效能的快速计算,使电磁兼容和散热的协同设计成为可能,有效提高了产品的设计周期。
Description
技术领域
本发明属于电子设备电磁兼容特性技术领域,尤其涉及一种基于等效电路的带孔阵腔体屏蔽效能计算方法。
背景技术
电磁屏蔽用于阻挡电磁能量的传输,是抑制电磁干扰、防止能量泄漏的重要手段。由于输入/输出,控制面板及通风散热等需要,屏蔽腔体上通常会存在各种形状和尺寸的孔及孔阵。这些孔在适当的电磁频率下,会成为高效的辐射天线,并成为EMI(ElectromagneticInterference)问题的干扰源,是影响腔体屏蔽特性的关键因素。
解析方法和数值方法可用于带孔腔体屏蔽效能的计算。解析方法计算速度快,结果较准,但只能用于简单腔体。解析方法中,Bethe提出的小孔衍射理论将无限大零厚度平面上的小孔看成是电偶极子和磁偶极子的组合,Mendez对其理论进行了改进,但只适用于低频范围。近年来等效电路方法得到广泛关注,在计算单孔腔体的屏蔽效能时具有很高的精度及效率,但处理孔阵时误差较大。数值方法包括矩量法、有限元法、时域有限差分法、传输线矩阵法等,广泛用于宽频带上带孔阵复杂腔体的屏蔽效能计算,但是需要牺牲大量的内存及计算时间。在对带孔阵腔体屏蔽效能进行数值计算时,由于组成孔阵的小孔数量众多且与腔体总体尺寸相差很大,导致腔体模型的网格数量巨大,特别是在大型复杂腔体(如机舱)的电磁特性分析中,往往由于存在大量小孔而导致整机的仿真无法进行。
发明内容
本发明的发明目的是:为了解决现有技术中孔阵建模困难及等效电路方法在计算带孔阵腔体屏蔽效能精度不高等问题,本发明提出了一种基于等效电路的带孔阵腔体屏蔽效能计算方法。
本发明的技术方案是:本发明的基于等效电路的带孔阵腔体屏蔽效能计算方法,包括以下步骤:
A、根据带孔阵腔体的孔阵参数构建等效单孔模型;
B、根据步骤A中构建的等效单孔模型构建其等效电路模型;
C、计算步骤A中构建的等效单孔模型中等效单孔的阻抗,并代替孔阵的阻抗;
D、根据步骤B中构建的等效电路模型计算带孔阵腔体的屏蔽效能。
进一步地,所述步骤A根据带孔阵腔体的孔阵参数构建等效单孔模型,具体包括以下分步骤:
A1、根据带孔阵腔体中孔阵面积及组成孔阵的小孔长宽比,构建与孔阵面积及组成孔阵的小孔长宽比相同的单孔模型;
A2、根据孔阵中小孔的数量参数计算孔阵缩放比例,再将单孔模型的长宽尺寸按照孔阵缩放比例进行缩放,构建等效单孔模型。
进一步地,所述孔阵缩放比例的计算公式,具体为:
k=1.283×n-0.1407-0.2829
其中,k为缩放比例,n为组成孔阵的小孔数量。
进一步地,所述步骤B根据步骤A中构建的等效单孔模型构建其等效电路模型,具体为:将等效单孔等效成共面的带状传输线,将带孔阵腔体等效成终端短路的波导。
本发明具有以下有益效果:
(1)本发明通过构建等效单孔模型,在保证仿真精度的前提下降低了建模复杂度,显著减少了网格数量,提高了仿真效率,简化了带孔阵复杂腔体的建模过程;
(2)本发明构建的等效单孔模型,大幅度的提高了等效电路方法在计算带孔阵腔体屏蔽效能时的精度;与孔阵导纳模型对比,可以不受小孔形状及孔阵形式的限制;与广义模矩量法(GMMOM)及专业的电磁仿真软件CST的结果对比,可以在保持计算精度的前提下减少计算时间提高仿真效率,实现带孔阵腔体屏蔽效能的快速计算;
(3)本发明的基于等效电路的带孔阵腔体屏蔽效能计算方法具有快速性和准确性的特定,使电磁兼容和散热的协同设计成为可能,有效提高了产品的设计周期。
附图说明
图1是本发明的基于等效电路的带孔阵腔体屏蔽效能计算方法流程示意图。
图2-a是孔阵等效前带孔阵腔体内的电场分布结构示意图。
图2-b是孔阵等效后带孔阵腔体内的电场分布结构示意图。
图3-a是孔阵等效前的网格结构示意图。
图3-b是孔阵等效后的网格结构示意图。
图4-a是现有的孔阵模型结构示意图。
图4-b是等效单孔模型结构示意图。
图4-c是等效单孔模型的等效电路模型结构示意图。
图5-a是本发明实施例的孔阵模型结构示意图。
图5-b是本发明实施例的孔阵等效单孔模型结构示意图。
图5-c是本发明与现有技术的计算结果对比示意图。
图6-a是本发明与GMMOM方法的计算结果对比示意图。
图6-b是本发明与电磁仿真软件CST的计算结果对比示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1所示,为本发明的基于等效电路的带孔阵腔体屏蔽效能计算方法流程示意图。本发明的基于等效电路的带孔阵腔体屏蔽效能计算方法,包括以下步骤:
A、根据带孔阵腔体的孔阵参数构建等效单孔模型;
B、根据步骤A中构建的等效单孔模型构建其等效电路模型;
C、计算步骤A中构建的等效单孔模型中等效单孔的阻抗,并代替孔阵的阻抗;
D、根据步骤B中构建的等效电路模型计算带孔阵腔体的屏蔽效能。
在步骤A中,由于在对带孔阵腔体的电磁屏蔽特性进行研究时,发现同面积下带孔阵腔体屏蔽效能高于单孔,且孔阵的小孔数量越多(尺寸越小),腔体的屏蔽效能越高;其次,在腔体第4个主谐振点频率以下,随着孔数量的增多,屏蔽效能只是幅值的提高,而谐振点不受影响;而产生这种现象的原因主要是小孔的尺寸效应。如图2-a所示,为孔阵等效前带孔阵腔体内的电场分布结构示意图。如图2-b所示,为孔阵等效后带孔阵腔体内的电场分布结构示意图。本发明在对带孔阵腔体的电磁屏蔽特性进行数值计算时,为了减少模型网格数量并降低孔阵建模复杂度,提出一种孔阵等效单孔建模方法,具体包括以下分步骤:
A1、根据带孔阵腔体中孔阵面积及组成孔阵的小孔长宽比,构建与孔阵面积及组成孔阵的小孔长宽比相同的单孔模型;
A2、根据孔阵中小孔的数量参数计算孔阵缩放比例,再将单孔模型的长宽尺寸按照孔阵缩放比例进行缩放,构建等效单孔模型。
如图4-a所示,为现有的孔阵模型结构示意图;其中,a,b,c分别为带孔阵腔体机箱长、高、宽尺寸,x,y分别为组成孔阵的小孔长宽尺寸,p为监测点P距离前面板的距离。这里的单孔面积与孔阵面积相同,单孔长宽比与组成孔阵的小孔长宽比相同,与孔阵的外包络形状无关。
如图5-a所示,为本发明实施例的孔阵模型结构示意图,机箱尺寸为300mm×120mm×300mm,单孔直径为12mm,水平间距为23mm,垂直间距为20mm,孔数量为20,孔阵中心位于机箱面板的中心位置。构建与孔阵同面积的单孔模型,通过计算得到单孔的直径为53.666mm,并保持单孔中心位置与孔阵中心位置相同。
在步骤A2中,孔阵缩放比例的计算公式,具体为:
k=1.283×n-0 . 1407-0.2829
其中,k为缩放比例,n为组成孔阵的小孔数量,1.283、0.1407及0.2829为根据大量仿真实验数据拟合得到的数据。本发明的孔阵缩放比例计算公式适应性好,精确度高。
根据组成孔阵的小孔数量结合孔阵缩放比例的计算公式即可计算得到孔阵缩放比例。再对步骤A中构建的单孔模型的长宽尺寸根据孔阵缩放比例进行等比例缩放,从而构建孔阵的等效单孔模型。如图4-b所示,为等效单孔模型结构示意图,l,w分别为等效单孔的长宽尺寸。构建孔阵的等效单孔模型,根据孔阵缩放比例的计算公式计算得到孔阵的缩放比例为0.5588,对同面积的等效单孔进行缩放,得到孔阵等效单孔的直径为30mm,如图5-b所示,为本发明实施例的孔阵等效单孔模型结构示意图。
本发明所构建的等效单孔模型,可以直接用于电磁场的仿真软件(如CST、HFSS、FEKO等)计算,即将等效单孔模型导入电磁学仿真软件中计算带孔阵腔体内部监测点的屏蔽效能,可以在保证等效前后腔体场分布相同的前提下,大幅度的降低网格数量,提高计算效率。
在步骤C中,具体为:
C1、根据步骤B中构建的等效单孔模型构建其等效电路模型;
C2、计算等效单孔的阻抗,并代替孔阵的阻抗;
C3、计算带孔阵腔体的屏蔽效能。
在步骤B中,根据步骤A中构建的等效单孔模型构建其等效电路模型,具体为将矩形机箱等效成一端短路的波导,等效单孔等效成共面的带状传输线。如图4-c所示,为等效单孔模型的等效电路模型结构示意图。
在步骤C中,计算等效单孔的阻抗的过程具体为:
孔缝的有效宽度we为
其中,w为等效单孔宽度,t为带孔阵腔体机箱厚度。
根据特性阻抗计算公式,可得屏蔽体特性阻抗Zos为
根据传输线理论可以计算得到孔缝的特性阻抗Zap为
其中,j为复数单位,k0=2π/λ,k0为自由空间的传输常数,λ为自由空间中电磁波波长,l为单孔的长度,a为屏蔽体长度。
将计算得到的等效单孔的阻抗代替孔阵的阻抗,突破了小孔形状及孔阵形式的限制,提高了等效电路方法在计算带孔阵腔体屏蔽效能时的精度,减少计算时间提高仿真效率。
在步骤D中,根据戴维宁定理,分别计算有无机箱时,监测点处的电压值,并利用屏蔽效能计算公式来计算得到监测点处的屏蔽效能值。屏蔽效能SE的计算公式具体为:
其中,vp为监测点P的电压值,v'p为没有屏蔽机箱时P点的电压值,v0为辐射源的电压值。
如图5-c所示,为本发明与现有技术的计算结果对比示意图,即利用等效电路方法计算带孔腔体的屏蔽效能,其结果与等效电路原模型、孔阵导纳模型及实验测试的对比结果。可以发现,本发明与孔阵导纳模型及实验测试结果非常的相近,而原来的等效电路方法由于只是将小孔阻抗线性相加,没有考虑小孔的尺寸效应及孔间耦合作用,导致计算结果与实验测试的结果相差较大。值得注意的是,孔阵导纳模型只能处理圆形孔和正方形孔正交或60°交叉排列的孔阵,且没有将机箱厚度的影响考虑进去,而本发明的方法可以处理任意的小孔形状及孔阵形式问题,并且考虑了机箱厚度的影响,应用范围更广。
本发明还进行了与广义模矩量法(GMMOM)及专业的电磁仿真软件CST求解结果的对比,如图6-a所示,为本发明与GMMOM方法的计算结果对比示意图。如图6-b所示,为本发明与电磁仿真软件CST的计算结果对比示意图。可以发现,利用本发明的计算方法的求解结果与GMMOM及CST的计算结果非常接近,而计算时间却相差很大,如表1所示。
表1.本发明与GMMOM及CST计算时间比较表
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (4)
1.一种基于等效电路的带孔阵腔体屏蔽效能计算方法,其特征在于,包括以下步骤:
A、根据带孔阵腔体的孔阵参数构建等效单孔模型;
B、根据步骤A中构建的等效单孔模型构建其等效电路模型;
C、计算步骤A中构建的等效单孔模型中等效单孔的阻抗,并代替孔阵的阻抗;
D、根据步骤B中构建的等效电路模型计算带孔阵腔体的屏蔽效能。
2.如权利要求1所述的基于等效电路的带孔阵腔体屏蔽效能计算方法,其特征在于,所述步骤A根据带孔阵腔体的孔阵参数构建等效单孔模型,具体包括以下分步骤:
A1、根据带孔阵腔体中孔阵面积及组成孔阵的小孔长宽比,构建与孔阵面积及组成孔阵的小孔长宽比相同的单孔模型;
A2、根据孔阵中小孔的数量参数计算孔阵缩放比例,再将单孔模型的长宽尺寸按照孔阵缩放比例进行缩放,构建等效单孔模型。
3.如权利要求2所述的基于等效电路的带孔阵腔体屏蔽效能计算方法,其特征在于,所述孔阵缩放比例的计算公式,具体为:
k=1.283×n-0.1407-0.2829
其中,k为缩放比例,n为组成孔阵的小孔数量。
4.如权利要求1所述的基于等效电路的带孔阵腔体屏蔽效能计算方法,其特征在于,所述步骤B根据步骤A中构建的等效单孔模型构建其等效电路模型,具体为:将等效单孔等效成共面的带状传输线,将带孔阵腔体等效成终端短路的波导。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510428920.8A CN105183939A (zh) | 2015-07-21 | 2015-07-21 | 基于等效电路的带孔阵腔体屏蔽效能计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510428920.8A CN105183939A (zh) | 2015-07-21 | 2015-07-21 | 基于等效电路的带孔阵腔体屏蔽效能计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105183939A true CN105183939A (zh) | 2015-12-23 |
Family
ID=54906018
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510428920.8A Pending CN105183939A (zh) | 2015-07-21 | 2015-07-21 | 基于等效电路的带孔阵腔体屏蔽效能计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105183939A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107194159A (zh) * | 2017-05-04 | 2017-09-22 | 电子科技大学 | 外部激励下带孔阵腔体电磁谐振的解析方法 |
CN107255755A (zh) * | 2017-06-27 | 2017-10-17 | 华北电力大学 | 一种计算开孔矩形屏蔽体内部场分布的分区解析模型 |
CN110058183A (zh) * | 2019-05-16 | 2019-07-26 | 华北电力大学 | 一种小孔磁极化系数的提取方法 |
CN110489866A (zh) * | 2019-08-20 | 2019-11-22 | 全球能源互联网研究院有限公司 | 一种计算封闭腔体磁场屏蔽效能的预测方法 |
CN110502864A (zh) * | 2019-08-29 | 2019-11-26 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机机体及其电磁屏蔽效能评估方法 |
CN111259521A (zh) * | 2020-01-09 | 2020-06-09 | 华北电力大学 | 一种矩形金属腔体屏蔽盒上开孔等效阻抗提取方法及系统 |
CN112464602A (zh) * | 2020-12-03 | 2021-03-09 | 电子科技大学 | 一种多谐振点电阻及电感的高频spice模型建立方法 |
CN113056184A (zh) * | 2021-03-15 | 2021-06-29 | 上海应用技术大学 | 一种基于优化干扰信号的屏蔽网仿真设计方法及整体架构 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060095889A1 (en) * | 2004-10-29 | 2006-05-04 | Synopsys, Inc. | Silicon tolerance specification using shapes as design intent markers |
CN102411659A (zh) * | 2011-11-25 | 2012-04-11 | 上海华虹Nec电子有限公司 | 硅片通孔等效电路模型及模型参数提取方法 |
-
2015
- 2015-07-21 CN CN201510428920.8A patent/CN105183939A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060095889A1 (en) * | 2004-10-29 | 2006-05-04 | Synopsys, Inc. | Silicon tolerance specification using shapes as design intent markers |
CN102411659A (zh) * | 2011-11-25 | 2012-04-11 | 上海华虹Nec电子有限公司 | 硅片通孔等效电路模型及模型参数提取方法 |
Non-Patent Citations (2)
Title |
---|
任丹等: ""一种考虑小孔尺寸效应的孔阵等效建模方法"", 《物理学报》 * |
罗静雯等: ""一种基于BLT方程的孔缝箱体屏蔽效能计算方法"", 《物理学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107194159A (zh) * | 2017-05-04 | 2017-09-22 | 电子科技大学 | 外部激励下带孔阵腔体电磁谐振的解析方法 |
CN107194159B (zh) * | 2017-05-04 | 2020-09-15 | 电子科技大学 | 外部激励下带孔阵腔体电磁谐振的解析方法 |
CN107255755A (zh) * | 2017-06-27 | 2017-10-17 | 华北电力大学 | 一种计算开孔矩形屏蔽体内部场分布的分区解析模型 |
CN107255755B (zh) * | 2017-06-27 | 2019-09-20 | 华北电力大学 | 一种计算开孔矩形屏蔽体内部场分布的分区解析模型的方法 |
CN110058183A (zh) * | 2019-05-16 | 2019-07-26 | 华北电力大学 | 一种小孔磁极化系数的提取方法 |
CN110058183B (zh) * | 2019-05-16 | 2020-06-19 | 华北电力大学 | 一种小孔磁化系数的提取方法 |
CN110489866A (zh) * | 2019-08-20 | 2019-11-22 | 全球能源互联网研究院有限公司 | 一种计算封闭腔体磁场屏蔽效能的预测方法 |
CN110502864A (zh) * | 2019-08-29 | 2019-11-26 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机机体及其电磁屏蔽效能评估方法 |
CN111259521A (zh) * | 2020-01-09 | 2020-06-09 | 华北电力大学 | 一种矩形金属腔体屏蔽盒上开孔等效阻抗提取方法及系统 |
CN111259521B (zh) * | 2020-01-09 | 2022-03-04 | 华北电力大学 | 一种矩形金属腔体屏蔽盒上开孔等效阻抗提取方法及系统 |
CN112464602A (zh) * | 2020-12-03 | 2021-03-09 | 电子科技大学 | 一种多谐振点电阻及电感的高频spice模型建立方法 |
CN113056184A (zh) * | 2021-03-15 | 2021-06-29 | 上海应用技术大学 | 一种基于优化干扰信号的屏蔽网仿真设计方法及整体架构 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105183939A (zh) | 基于等效电路的带孔阵腔体屏蔽效能计算方法 | |
Wang et al. | Multisolver domain decomposition method for modeling EMC effects of multiple antennas on a large air platform | |
CN104850738B (zh) | 架空电力线路杆塔处雷电感应电压计算方法 | |
CN102770009B (zh) | 一种吸波超材料 | |
CN104537261B (zh) | 一种基于场路耦合的含贯穿线缆机箱内电路干扰预测方法 | |
CN104200074A (zh) | 快速获取目标电磁散射特性的多层复波束方法 | |
CN104573271B (zh) | 一种参数可变的精确等离子体天线辐射性能仿真方法 | |
CN107255755A (zh) | 一种计算开孔矩形屏蔽体内部场分布的分区解析模型 | |
CN106446470B (zh) | 一种高效并行的非均匀介质频域有限差分方法 | |
CN105844019B (zh) | 一种基于车辆表面电流分布规律的接地点选择方法 | |
CN104915326A (zh) | 基于等效原理的区域分解阶数步进时域积分方法 | |
Zhang et al. | Effect of frequency-dependent soil on the propagation of electromagnetic fields radiated by subsequent lightning strike to tall objects | |
CN104750960A (zh) | 一种快速提取金属桁架式天线罩电性能参数的方法 | |
Kavanagh et al. | A method of moments based indoor propagation model | |
CN105205299B (zh) | 电大目标电磁散射特性快速降维分析方法 | |
CN114298433A (zh) | 一种变电站内5g基站选址方法、装置、设备及介质 | |
Asadpourahmadchali et al. | Hybrid continuous circuit‐trapezoidal integration method analysis of multi‐cross structure of grounding system | |
Xu | Computational electrodynamics and simulation in high speed circuit using finite difference time domain (FDTD) method | |
CN105977119A (zh) | 一种栅控电子枪虚阳极位置的确定方法 | |
He et al. | Analysis of Shielding Effectiveness of Aperture Based on Multilayer Metal and Ferrite Composite Structure | |
Xie et al. | The electromagnetic equivalent modeling and simulation analysis of large size metal structure | |
Wei et al. | Analysis of magnetic fields radiated by lightning strikes to a building | |
Fallahpour et al. | Numerical modeling of electromgantic interference (EMI) and radiation problems for large multiscale structures | |
Liu et al. | Indoor environmental signal strength calculation based on finite difference time domain method | |
Liu et al. | Calculation of shielding coefficient of grounded metal mesh under high voltage transmission lines by iterative charge simulation method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151223 |