CN107657082B - 一种低电平电场照射下非封闭舱室内均匀场区域划定方法 - Google Patents
一种低电平电场照射下非封闭舱室内均匀场区域划定方法 Download PDFInfo
- Publication number
- CN107657082B CN107657082B CN201710773723.9A CN201710773723A CN107657082B CN 107657082 B CN107657082 B CN 107657082B CN 201710773723 A CN201710773723 A CN 201710773723A CN 107657082 B CN107657082 B CN 107657082B
- Authority
- CN
- China
- Prior art keywords
- region
- standard deviation
- cabin
- electric field
- field
- 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.)
- Active
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Medicines That Contain Protein Lipid Enzymes And Other Medicines (AREA)
- Testing Electric Properties And Detecting Electric Faults (AREA)
Abstract
本发明公开了一种低电平电场照射下非封闭舱室内均匀场区域划定方法,包括以下过程:建立非封闭舱室和模式搅拌器的电磁仿真模型,设置舱室内待划定区域;数值分析外部低电平电场照射下舱室内的场强分布,得到待划定区域的场强数据;基于测试标准中的均匀性统计方法,计算区域内观察点的标准偏差值;判断观察点的标准偏差值是否满足均匀性要求,若否,继续缩小待划定区域的体积直至满足均匀性要求。本发明具有解决飞机级低电平扫描场测试时实际舱室非封闭结构带来的难题,不同于封闭混响室结构,考虑了舱室开口对能量密度和场均匀性的影响,为低电平扫描场测试时飞机舱室内接收天线/探头与搅拌器的位置布局提供有效支撑的优点。
Description
技术领域
本发明涉及电磁环境效应技术领域,特别涉及一种低电平电场照射下非封闭舱室内均匀场区域划定方法。
背景技术
飞机上更多的弱电设备取代传统机械或机电控制系统,更多的低屏蔽性能复合材料应用,以及机上开口面、缝隙的存在,使得高强电磁信号更容易耦合进入机内造成干扰、损伤甚至损坏,高强辐射场(HIRF)已成为影响飞机飞行安全的重要因素。欧美、中国等航空局机构在近期颁布了相应条款和专用条件强制规定了民用航空器必须通过HIRF适航符合性试验验证。作为飞机级HIRF试验的重要内容,低电平扫描场测试具有试验便捷、降低待测飞机潜在风险等优点,其测试方法是采用外部电磁波照射飞机舱室,接收天线或场强探头置于舱室内测试内部场强,并通过与校准值对比获得飞机舱室的衰减特性。
低电平扫描场测试过程中由于舱室存在驻波,舱室内不同位置场强变化较大,为了获得舱室内最大场强,以往的做法是在舱室内安装多个的接收天线/探头,或者采用有限接收天线在不同位置开展多次测试,这些方法不仅耗时而且耗力。目前比较公认的测试方法是在舱室内布置搅拌器,通过旋转改变附近的边界条件,使得搅拌器的一个旋转周期内接收天线/探头位置出现场强最大值,达到统计均匀的目的。
专利方面,2015年公开的发明专利“一种用于混响室内搅拌器的快速优化仿真设计方法”(申请号:201510802353.8)通过仿真设计优化搅拌器实现混响室内场强均匀性。2016年公开的发明专利“一种电磁混响室腔体的尺寸优化方法和装置”(申请号:201610648278.9)在设计电磁混响室通过机械搅拌浆边缘切割改善磁场的均匀性。2007年电波科学学报第22卷第4期论文“混响室设计与校准测试”提出了长方体结构混响室设计和均匀性验证方法。目前的文献涉及的搅拌器优化和均匀场分析限定在封闭混响室结构,其方法不能应用于实际的非封闭的飞机舱室的低电平扫描场测试情况。2008年高电压技术第34卷第8期论文“混响室测试区域均匀性分部规律仿真分析”和2015年科学技术与工程第15卷第34期论文“通风孔对混响室场均匀性的影响”研究了舱室内均匀性区域规律和开口对舱室内的局部能量密度和总场均匀性的影响,但是文献仍局限在发射天线置于混响室内部的情况,没有开展外部电磁波照射下舱室内均匀性分析。
另外,实际飞机舱室LLSF测试过程中,很难针对每个待测舱室设计不同尺寸的搅拌器。
因此,需要提出一种低电平电场照射下非封闭舱室内均匀场区域划定方法,解决实际舱室非封闭结构带来的均匀场区域变化的难题,为低电平扫描场测试时飞机舱室内接收天线/探头与搅拌器的位置布局提供有效支撑。
发明内容
本发明的目的是提出一种低电平电场照射下非封闭舱室内均匀场区域划定方法,该方法结合数值分析和均匀性统计方法,通过逐渐缩小待划定区域寻找均匀场有效区域,解决了实际非封闭舱室内场均匀性受舱室开口影响的问题。在考虑外部电磁波照射下,结合数值分析和均匀性统计方法寻找均匀场有效区域的区域划定方法,首先建立非封闭舱室和模式搅拌器的电磁仿真模型,设置舱室内待划定区域,然后数值分析外部低电平电场照射下舱室内的场强分布,得到待划定区域的场强数据,最后基于测试标准中的均匀性统计方法,计算区域内观察点的标准偏差值,验证并指导待划定区域的体积调整以满足均匀性要求,通过考虑外部电磁照射下舱室开口对能量密度和场均匀性的影响,为低电平扫描场测试时飞机舱室内接收天线/探头与搅拌器的位置布局提供有效支撑。
为了实现以上目的,本发明通过以下技术方案实现:
一种低电平电场照射下非封闭舱室内均匀场区域划定方法,其特征在于,包括以下过程:建立非封闭舱室和模式搅拌器的电磁仿真模型,设置舱室内待划定区域。数值分析外部低电平电场照射下舱室内的场强分布,得到待划定区域的场强数据。基于测试标准中的均匀性统计方法,计算区域内观察点的标准偏差值。判断观察点的标准偏差值是否满足均匀性要求,若否,继续缩小待划定区域的体积直至满足均匀性要求。
优选地,所述非封闭舱室和模式搅拌器的电磁仿真模型是将多个模式搅拌器置于非封闭舱室内部,每一模式搅拌器采用步进工作模式。
优选地,所述待划定区域是一个长方体区域,长方体区域的边界与非封闭舱室的距离至少为λ/4,其中λ为最小频率的波长;所述长方体区域的边界与每一模式搅拌器旋转形成的圆柱边界之间的距离至少为λ/4。
优选地,所述数值分析外部低电平电场照射下舱室内的场强分布通过采用矩量法计算1V/m场强照射下非封闭舱室内的场强值分布。
优选地,所述待划定区域的场强数据包含每一模式搅拌器步进状态下区域内观察点的场强数据,每个区域内观察点的场强数据包括笛卡尔坐标系下x、y、z三个方向的场强分量和总场强值。
优选地,所述基于测试标准中的均匀性统计方法为参考RTCA DO-160G标准计算标准偏差值。
优选地,所述区域观察点的标准偏差值的计算包含以下过程:
所述区域内观察点的最大场强值归一化到平均输入功率的平方根:
式中:EMaxx,y,z为每个区域内观察点沿x、y、z轴方向电场分量的最大测量值;为每个区域内观察点沿x、y、z轴方向电场分量归一化的最大测量值;Pinput-empty为平均输入功率;针对每个频率,计算区域内观察点每个分量的电场最大测量值的归一化平均值:
以及所有电场最大测量值的归一化平均值:
式中:<>表示算术平均值。
优选地,所述区域观察点的标准偏差值的计算进一步包含以下过程:在模式搅拌器旋转一周时,以九个所述区域内观察点中每个区域内观察点得到的最大测量值的归一化平均值的标准偏差确定电场均匀性,采用区域内观察点每个分量的单独数据和总数据集计算标准偏差;
标准偏差σ:
用与平均值有关的dB数表示标准偏差:
优选地,所述区域内观察点是长方体区域的八个顶点和一个长方体区域中心点。
优选地,所述标准偏差值是否满足均匀性要求为所有的标准偏差值低于3dB,所述缩小待划定区域的体积为依次缩小长方体区域的长、宽、高,若区域观察点的标准偏差值低于3dB,则划定的区域满足均匀性要求,完成均匀场区域划定,否则,未能满足均匀性要求,则继续缩小所述待划定区域的体积再进行判断,直至当前区域观察点的标准偏差值低于3dB。
本发明与现有技术相比具有以下优点:
采用数值分析的方法得出非封闭舱室内场强,用于划定均匀场区域,可以快速得出结果,指导试验。
均匀场强区域划定过程中,仿真模型参数、外部低电平电场照射条件、待划定区域易设置调整,方便根据实际需求进行分析。
寻找均匀场区域时,采用笛卡尔坐标系三个方向依次缩小,简单方便,结合均匀性统计方法能快速获得均匀场有效区域。
附图说明
图1是本发明的低电平电场照射下非封闭舱室内均匀场区域划定方法的流程图;
图2是实施例中非封闭舱室和搅拌器的仿真模型;
图3是实施例中待划定区域示意;
图4是实施例中待划定区域的标准偏差值计算结果。
具体实施方式
以下结合附图,通过详细说明一个较佳的具体实施例,对本发明做进一步阐述。
如图1所示,本发明一种低电平电场照射下非封闭舱室内均匀场区域划定方法,包含以下过程:
S1、建立非封闭舱室和模式搅拌器的电磁仿真模型,设置舱室内待划定区域。
所述的非封闭舱室和模式搅拌器的电磁仿真模型是模式搅拌器置于非封闭舱室内部,模式搅拌器采用步进工作模式。
所述的待划定区域是一个长方体区域,长方体区域的边界与舱室的距离至少为λ/4,其中λ为最小频率的波长,同时,长方体区域的边界与搅拌器旋转形成的圆柱边界的距离至少为λ/4。
S2、数值分析外部低电平电场照射下舱室内的场强分布,得到待划定区域的场强数据。
所述的数值分析外部低电平电场照射下舱室内的场强分布是采用矩量法计算1V/m场强照射下舱室内的场强值分布。
所述的待划定区域的场强数据包括了每一个搅拌器步进状态下区域内的场强数据,区域内每个点的场强数据包括笛卡尔坐标系下x、y、z三个方向的场强分量和总场强值。
S3、基于测试标准中的均匀性统计方法,计算区域内观察点的标准偏差值。
所述的基于测试标准中的均匀性统计方法是参考RTCADO-160G标准计算标准偏差值。
所述的区域内观察点是长方体区域的八个顶点和一个长方体区中心点。
计算标准偏差值时,将九个观察点的最大场强值归一化到平均输入功率的平方根:
其中:EMaxx,y,z为每个观察点x、y、z分量的最大测量值;为每个观察点x、y、z分量归一化的最大测量值。Pinput-empty表示在搅拌器旋转期间,记录的EMaxx,y,z时舱室内平均输入功率。
针对每个频率,计算观察点每个分量的最大测量值的归一化平均值:
以及所有最大测量值的归一化平均值:
式中:<>表示算术平均值。
在搅拌器旋转一周期间,以九个观察点中每个观察点得到的最大测量值的归一化平均值的标准偏差来确定场均匀性,用观察点每个分量的单独数据和总数据集来计算标准偏差。
标准偏差由以下公式给出:
α在n≤20时为1.06,n>20时为1。
用与平均值有关的dB数表示标准偏差:
S4、判断观察点的标准偏差值是否满足均匀性要求,如不满足,继续缩小待划定区域的体积直至满足均匀性要求。
如果得到观察点的标准偏差值低于3dB,则说明划定的区域满足均匀性要求;否则未能满足均匀性要求,继续缩小待划定区域的体积再进行判断。
下面以一个实施例说明本发明方法在非封闭模拟舱室均匀场区域划定中的应用,按照如下方法进行区域划定:
步骤一:建立非封闭舱室和模式搅拌器的电磁仿真模型,设置舱室内待划定区域。非封闭模拟舱室电磁仿真模型如图2a所示,长宽高三维尺寸为4m×2m×2m,在舱室左右两侧分别有五个0.3m×0.2m的窗口,舱室前侧有一个1m×0.4m的开口,舱室材料为金属。舱室内搅拌器的模型如图2b所示,由两个0.2m×0.2m正方形板折叠形成,模型下方为驱动盒和转轴,搅拌叶片的底座中心置于舱室中心位置,搅拌器以15°步进旋转。
如图3所示外部电磁波e以平面波、垂直极化方式从舱室一侧垂直照射,待分析的均匀区域f,待划定区域与舱室1距离至少为λ/4,与搅拌器2旋转形成的圆柱边界的距离至少为λ/4,根据最低分析频率400MHz,距离设置为0.2m,在直角坐标系中,待划定区域为:x∈(0.5m,1.8m),y∈(-0.8m,0.8m),z∈(-0.8m,0.8m),x、y、z方向步进为0.1m。
步骤二:数值分析外部低电平电场照射下舱室内的场强分布,得到待划定区域的场强数据。
本实施例中,利用FEKO仿真软件中的矩量法,分析了400MHz~1GHz频率范围舱室内的场强,得到待划定区域内的场强数据。
步骤三:基于测试标准中的均匀性统计方法,计算区域内观察点的标准偏差值。
如图4所示,在本实施例中,利用公式(1)~(7)计算了400MHz~1GHz频率范围内九个观察点的标准偏差值,包括x、y、z三个分量Ex、Ey、Ez的标准偏差值和总场Etotal的标准偏差值。
步骤四:判断观察点的标准偏差值是否满足均匀性要求,如不满足,继续缩小待划定区域的体积直至满足均匀性要求。
在本实施例中,800MHz~1GHz,此待划定区域满足均匀性要求,400MHz~800MHz频率范围不满足均匀性要求,继续缩小划定区域:x∈(0.5m,0.8m),y∈(-0.5m,0.5m),z∈(-0.5m,0.5m),采用公式(1)~(7)计算后所有的标准偏差值低于3dB,满足均匀性要求。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (3)
1.一种低电平电场照射下非封闭舱室内均匀场区域划定方法,其特征在于,包括以下过程:
建立非封闭舱室和模式搅拌器的电磁仿真模型,设置舱室内待划定区域;
数值分析外部低电平电场照射下舱室内的场强分布,得到待划定区域的场强数据;
基于测试标准中的均匀性统计方法,计算区域内观察点的标准偏差值;
判断观察点的标准偏差值是否满足均匀性要求,若否,继续缩小待划定区域的体积直至满足均匀性要求;所述非封闭舱室和模式搅拌器的电磁仿真模型是将多个模式搅拌器置于非封闭舱室内部,每一模式搅拌器采用步进工作模式;所述待划定区域是一个长方体区域,长方体区域的边界与非封闭舱室的距离至少为λ/4,其中λ为最小频率的波长;所述长方体区域的边界与每一模式搅拌器旋转形成的圆柱边界之间的距离至少为λ/4;所述数值分析外部低电平电场照射下舱室内的场强分布通过采用矩量法计算1V/m场强照射下非封闭舱室内的场强值分布;所述待划定区域的场强数据包含每一模式搅拌器步进状态下区域内观察点的场强数据,每个区域内观察点的场强数据包括笛卡尔坐标系下x、y、z三个方向的场强分量和总场强值;所述基于测试标准中的均匀性统计方法为参考RTCA DO-160G标准计算区域观察点的电场标准偏差值;所述区域内观察点是长方体区域的八个顶点和一个长方体区域中心点;所述标准偏差值是否满足均匀性要求为所有的标准偏差值低于3dB;
所述缩小待划定区域的体积为依次缩小长方体区域的长、宽、高;
若区域观察点的标准偏差值低于3dB,则划定的区域满足均匀性要求,完成均匀场区域划定;
否则,未能满足均匀性要求,则继续缩小所述待划定区域的体积再进行判断,直至当前区域观察点的标准偏差值低于3dB。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710773723.9A CN107657082B (zh) | 2017-08-31 | 2017-08-31 | 一种低电平电场照射下非封闭舱室内均匀场区域划定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710773723.9A CN107657082B (zh) | 2017-08-31 | 2017-08-31 | 一种低电平电场照射下非封闭舱室内均匀场区域划定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107657082A CN107657082A (zh) | 2018-02-02 |
CN107657082B true CN107657082B (zh) | 2021-04-02 |
Family
ID=61129059
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710773723.9A Active CN107657082B (zh) | 2017-08-31 | 2017-08-31 | 一种低电平电场照射下非封闭舱室内均匀场区域划定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107657082B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109709403A (zh) * | 2018-11-12 | 2019-05-03 | 中航通飞研究院有限公司 | 一种飞机低电平扫描场测试的数据处理方法 |
CN109541512B (zh) * | 2018-12-20 | 2021-03-23 | 北京无线电计量测试研究所 | 一种用于混响室装置的搅拌器及混响室装置 |
CN110728023B (zh) * | 2019-09-06 | 2023-11-21 | 久瓴(江苏)数字智能科技有限公司 | 老虎窗构件节点生成方法、装置、计算机设备和存储介质 |
CN111337756A (zh) * | 2020-04-24 | 2020-06-26 | 中国人民解放军63892部队 | 一种评价脉冲激励混响室测试区域电场均匀性方法 |
CN112114218B (zh) * | 2020-09-22 | 2024-02-23 | 上海无线电设备研究所 | 一种用于hirf试验的高电平试验系统及试验和验证方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104392023A (zh) * | 2014-10-30 | 2015-03-04 | 上海无线电设备研究所 | 一种高强辐射场条件下飞机舱体电磁模型校验方法 |
CN105184108A (zh) * | 2015-10-27 | 2015-12-23 | 上海无线电设备研究所 | 一种hirf条件下复合材料舱室内场强计算方法 |
-
2017
- 2017-08-31 CN CN201710773723.9A patent/CN107657082B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104392023A (zh) * | 2014-10-30 | 2015-03-04 | 上海无线电设备研究所 | 一种高强辐射场条件下飞机舱体电磁模型校验方法 |
CN105184108A (zh) * | 2015-10-27 | 2015-12-23 | 上海无线电设备研究所 | 一种hirf条件下复合材料舱室内场强计算方法 |
Non-Patent Citations (4)
Title |
---|
Environmental Conditions and Test Procedures for Airborne Equipment;RTCA;《RTCA DO-160G标准》;20101208;第20-16-20-41页 * |
混响室测试区场均匀性分布规律仿真分析;张成怀等;《高压电技术》;20080831;第34卷(第8期);第1537-1541页 * |
混响室电磁环境的GPU并行重建算法及应用研究;白冰;《中国博士学位论文全文数据库 信息科技辑》;20160315(第3期);第I135-2页 * |
通风孔对混响室场均匀性的影响;陈亮等;《科学技术与工程》;20151231;第15卷(第34期);第1-6页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107657082A (zh) | 2018-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107657082B (zh) | 一种低电平电场照射下非封闭舱室内均匀场区域划定方法 | |
Bruns et al. | A closer look at reverberation chambers-3-D simulation and experimental verification | |
US6525657B1 (en) | Apparatus and method for production testing of the RF performance of wireless communications devices | |
Georgakopoulos et al. | HIRF penetration through apertures: FDTD versus measurements | |
CN108132390B (zh) | 一种提升混响室内场强和模拟复合场电磁环境的方法及系统 | |
CN104392023B (zh) | 一种高强辐射场条件下飞机舱体电磁模型校验方法 | |
WO2009046516A1 (en) | Multichannel absorberless near field measurement system | |
Dobychina et al. | Unmanned aerial vehicle antenna measurement using anechoic chamber | |
CN109657277A (zh) | 一种舱室结构下复合材料电磁参数提取方法 | |
Bunting | Shielding effectiveness in a two-dimensional reverberation chamber using finite-element techniques | |
CN111487474B (zh) | 数值孪生电磁测量系统 | |
Hamdalla et al. | Characteristic mode analysis prediction and guidance of electromagnetic coupling measurements to a UAV model | |
Brüns et al. | Modeling challenging EMC problems | |
CN115236411A (zh) | 一种屏蔽罩内电子器件的辐射源建模方法 | |
Faust et al. | Analysis, Modeling, and Measurement of Shielding Effectiveness for a Cylindrical Waveguide with a Hexagonal Insert Structure. | |
Gradoni et al. | Determination of the reverberation chamber stirrer uncorrelated positions by means of the spatial and frequency correlation matrix | |
García-Pérez et al. | Time-domain shielding effectiveness of enclosures against a plane wave excitation | |
Le Goff et al. | A novel and innovative near field system for testing radomes of commercial aircrafts | |
Gutierrez et al. | Influence of geometric simplifications on high-intensity radiated field simulations | |
Bakti et al. | Field uniformity and correlation coefficient analysis of KRISS reverberation chamber | |
Scialacqua et al. | Diagnostics on electrically large structures by a nested skeletonization scheme enhancement of the equivalent current technique | |
Pei et al. | Electric field uniformity of cabin containing apertures under low level swept fields | |
Poodineh et al. | Standard Deviation Characterization of a Small Size Reverberation Chamber by Using Full wave Simulation and E Field Probe | |
CN105676006B (zh) | 一种基于实测数据的电子设备危害辐射场预测方法及系统 | |
Řeznicek et al. | TD and FD simulations of internal EM environment in small aircraft and experimental test comparison |
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 | ||
CB02 | Change of applicant information |
Address after: No. 1555, Zhongchun Road, Minhang District, Shanghai, 200233 Applicant after: SHANGHAI RADIO EQUIPMENT Research Institute Address before: 200090 No. 203, Liping Road, Shanghai, Yangpu District Applicant before: SHANGHAI RADIO EQUIPMENT Research Institute |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |