CN105302992A - 一种用于混响室内搅拌器的快速优化仿真设计方法 - Google Patents

一种用于混响室内搅拌器的快速优化仿真设计方法 Download PDF

Info

Publication number
CN105302992A
CN105302992A CN201510802353.8A CN201510802353A CN105302992A CN 105302992 A CN105302992 A CN 105302992A CN 201510802353 A CN201510802353 A CN 201510802353A CN 105302992 A CN105302992 A CN 105302992A
Authority
CN
China
Prior art keywords
stirrer
incident
reverberation chamber
random
electric 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.)
Pending
Application number
CN201510802353.8A
Other languages
English (en)
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.)
Shanghai Radio Equipment Research Institute
Original Assignee
Shanghai Radio Equipment Research Institute
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 Shanghai Radio Equipment Research Institute filed Critical Shanghai Radio Equipment Research Institute
Priority to CN201510802353.8A priority Critical patent/CN105302992A/zh
Publication of CN105302992A publication Critical patent/CN105302992A/zh
Pending legal-status Critical Current

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

本发明公开的一种用于混响室内搅拌器的快速优化仿真设计方法,首先设计该搅拌器的外形结构,并确定要求优化设计的搅拌器参数;其次选取多个等效入射源,确定每个所述等效入射源的随机入射方向和极化角;采用混响室测试标准中的混响室均匀性统计方法,对搅拌器的所有观察点的电场矢量平均结果进行统计分析;根据不同的优化设计的搅拌器参数,选择搅拌器的最佳优化设计的搅拌器参数组形成的设计模型。本发明能够提高搅拌器的优化设计效率,避免搅拌器传统设计方法中存在的设计难点和设计局限性;并且本发明公开的快速优化仿真设计方法设计效率提高的量级与其应用到的混响室具体尺寸相关,混响室尺寸越大,设计效率越高。

Description

一种用于混响室内搅拌器的快速优化仿真设计方法
技术领域
本发明涉及电磁兼容性测试领域,具体涉及一种用于混响室内搅拌器的快速优化仿真设计方法。
背景技术
随着军用产品、汽车、航空工业产品的辐射抗扰度对场强强度、以及测试均匀性和入射辐射场具有方向、极化随机等要求,混响室相比于传统的测试场地(如开阔场地、半电波暗室、全电波暗室等),具有测试频带宽、重复性好、测试时间短以及在测试过程中采用较小的功率就能产生很高的场强等特点,许多相关学者开始对混响室及其设计进行了深入研究。此外包括CISPR、美国军用标准以及我国的国家标准等在内的标准化机构也已逐渐将混响室测试内容纳入标准中,并不断修正和完善。
混响室内部电场均匀特性是衡量混响室性能的重要指标。混响室腔体结构与尺寸确定以后,搅拌器的材料、形状、大小、个数都是影响混响室内部电场均匀性的重要因素。因此搅拌器的优化设计就成为进一步改善混响室场均匀性的重要手段。目前国内外学者已经在搅拌器的外形、叶片夹角、场均匀性区域的选择等方面展开了大量的研究,但这些研究大部分是对给定的几种搅拌器方案进行比较选择,得出了各自因素对混响室场均匀性独立影响的结论。
但在搅拌器优化设计方面,仍然需要对搅拌器和混响室的联立模型进行仿真。由于混响室一般均为超电大尺寸,且为高谐振结构,目前仿真分析混响室内的电场分布过程中计算量巨大、耗费时间漫长,且对计算机硬件配置要求高等难点,特别在某些高频段,有时用计算平台实现仿真分析几乎是不可能的,因此更加无法实现搅拌器快速优化设计。
发明内容
本发明的目的在于提供一种用于混响室内搅拌器的快速优化仿真设计方法,首先设计该搅拌器的外形结构,并确定要求优化设计的搅拌器参数;其次选取多个等效入射源,确定每个所述等效入射源的随机入射方向和极化角;采用混响室测试标准中的混响室均匀性统计方法,对搅拌器的所有观察点的电场矢量平均结果进行统计分析;根据不同的优化设计的搅拌器参数,选择搅拌器的最佳优化设计的搅拌器参数组形成的设计模型。本发明能够提高搅拌器的优化设计效率,避免搅拌器传统设计方法中存在的设计难点和设计局限性;并且本发明公开的快速优化仿真设计方法设计效率提高的量级与其应用到的混响室具体尺寸相关,混响室尺寸越大,设计效率越高。
为了达到上述目的,本发明通过以下技术方案实现:
一种用于混响室内搅拌器的快速优化仿真设计方法,该仿真设计方法包含:
S1,根据搅拌器的设计要求参数,设计该搅拌器的外形结构,并确定要求优化设计的搅拌器参数;
S2,选取多个等效入射源,在确保所有等效入射源能够模拟混响室内部随机电磁环境的前提下,确定每个所述等效入射源的随机入射方向和极化角;
S3,针对每个所述等效入射源进行仿真计算,获取以所述搅拌器为中心的球面上各个观察点的电场矢量分布计算结果;
S4,对所有所述等效入射源的电场矢量分布计算结果进行矢量平均计算,采用混响室测试标准中的混响室均匀性统计方法,对所述搅拌器的所有观察点的电场矢量平均结果进行统计分析;
S5,修改所述步骤S1中的所述要求优化设计的搅拌器参数,并重复执行所述步骤S3-S4,获取在多组不同优化设计的搅拌器参数下,对应的多组搅拌器参数的所有观察点的电场矢量平均结果的累加标准偏差;
S6,以电场矢量累加标准偏差越小、所述搅拌器的搅拌性能越好为评判原则,筛选获得所述搅拌器的最佳优化设计的搅拌器参数组形成的设计模型。
优选地,所述步骤S1包含:
所述搅拌器的设计要求参数包含:搅拌器的旋转空间参数、搅拌器的回转直径、搅拌器的最大尺寸及搅拌器的表面材料;
根据上述搅拌器的设计要求参数设计该搅拌器的外形结构,并确定要求优化设计的搅拌器参数为:所述搅拌器的桨叶长度、桨叶宽度及相邻桨叶之间的夹角。
优选地,所述步骤S2包含:
S2.1,设定第m个所述等效入射源的随机入射方向角为θm第m个所述等效入射源的随机极化角为ηm;其中,每个所述等效入射源的随机入射方向角、极化角均要求满足均匀分布统计特性的要求;
S2.2,要求第m个所述等效入射源的随机入射方向角θm的随机取值应在[0°,180°]之间;要求第m个所述等效入射源的随机入射方向角的随机取值应在[0°,360°]之间;
S2.3,要求第m个所述等效入射源的随机极化角ηx的随机取值应在[0°,360°]之间。
优选地,所述步骤S2.2包含:
S2.2.1,对于第m个所述等效入射源的随机入射方向角θm,在[0,1]内抽取随机数作为随机入射方向角θm的系数A,则θm=A·180°;
S2.2.2,对于第m个所述等效入射源的随机入射方向角的计算方法为:其中为随机取值,当的某一随机值使得与其对应的值落在[0°,360°]之间,则保留该值。
优选地,所述步骤S2.3包含:
对于第m个所述等效入射源的随机极化角ηm,在[0,1]内抽取随机数作为随机入射方向角ηm的系数B,则ηm=B·360°。
优选地,所述步骤S3包含:
分别针对每个所述等效入射源进行仿真计算,其中每次仿真计算对象仅为搅拌器模型,从而获取以所述搅拌器为中心的球面上所有观察点的电场矢量分布计算结果。
优选地,所述步骤S4包含:
S4.1,对所有所述等效入射源分别激励下的电场矢量计算结果进行矢量平均;
S4.2,采用混响室测试标准中的混响室均匀性统计方法,对所述搅拌器的所有观察点x、y、z三个极化方向进行计算获得三个极化方向的电场矢量标准偏差Dx、Dy及Dz,并计算上述三个极化方向的标准偏差Dtotal
优选地,所述步骤S5包含:
用户根据所述搅拌器的桨叶长度、桨叶宽度及相邻桨叶之间的夹角的不同,重复执行所述步骤S3-S4;获取多组不同优化设计的搅拌器参数情况下,对应的多个电场矢量平均结果的x、y、z三个极化方向的电场矢量标准偏差Dx、Dy、Dz及上述三个极化方向的标准偏差Dtotal
优选地,所述步骤S6包含:
S6.1,将所述步骤S5获取的每组优化设计的搅拌器参数对应的电场矢量平均结果的x、y、z三个极化方向的电场矢量标准偏差Dx、Dy、Dz及上述三个极化方向的标准偏差Dtotal进行如下处理计算,获取搅拌器电场矢量累加标准偏差D:
D=Dx+Dy+Dtotal
S6.2,通过比较所述步骤S5获取多组不同优化设计的搅拌器参数对应的搅拌器电场矢量累加标准偏差D,根据电场矢量累加标准偏差越小、所述搅拌器的搅拌性能越好为评判原则,筛选获得所述搅拌器的最佳优化设计的搅拌器参数组形成的设计模型。
本发明与现有技术相比具有以下优点:
本发明公开的一种用于混响室内搅拌器的快速优化仿真设计方法,首先设计该搅拌器的外形结构,并确定要求优化设计的搅拌器参数;其次选取多个等效入射源,确定每个所述等效入射源的随机入射方向和极化角;采用混响室测试标准中的混响室均匀性统计方法,对搅拌器的所有观察点的电场矢量平均结果进行统计分析;根据不同的优化设计的搅拌器参数,选择搅拌器的最佳优化设计的搅拌器参数组形成的设计模型。本发明能够提高搅拌器的优化设计效率,避免搅拌器传统设计方法中存在的设计难点和设计局限性;并且本发明公开的快速优化仿真设计方法设计效率提高的量级与其应用到的混响室具体尺寸相关,混响室尺寸越大,设计效率越高。本发明具有很强的实用性、操作性以及良好的应用前景,特别在电磁兼容性测试、航空航天器高强辐射场测试中舱室搅拌器设计、混响室内部搅拌器设计等领域中。
附图说明
图1为本发明一种用于混响室内搅拌器的快速优化仿真设计方法的整体流程示意图。
图2为本发明一种用于混响室内搅拌器的快速优化仿真设计方法的第m个等效入射源入射方向、极化角的示意图。
图3a为本发明一种用于混响室内搅拌器的快速优化仿真设计方法的现有技术等效入射源随机分布示意图。
图3b为本发明一种用于混响室内搅拌器的快速优化仿真设计方法的等效入射源随机分布示意图。
图4a为本发明一种用于混响室内搅拌器的快速优化仿真设计方法中根据搅拌器设计原则初步设计的搅拌器结构之一。
图4b为本发明一种用于混响室内搅拌器的快速优化仿真设计方法中根据搅拌器设计原则初步设计的搅拌器结构之二。
图4c为本发明一种用于混响室内搅拌器的快速优化仿真设计方法中根据搅拌器设计原则初步设计的搅拌器结构之三。
图4d为本发明一种用于混响室内搅拌器的快速优化仿真设计方法中根据搅拌器设计原则初步设计的搅拌器结构之四。
图5为本发明一种用于混响室内搅拌器的快速优化仿真设计方法的多组优化设计的搅拌器参数对应的累加标准偏差统计分析示意图。
图6a为本发明一种用于混响室内搅拌器的快速优化仿真设计方法的仿真示意图之一。
图6b为本发明一种用于混响室内搅拌器的快速优化仿真设计方法的仿真示意图之二。
具体实施方式
以下结合附图,通过详细说明一个较佳的具体实施例,对本发明做进一步阐述。
如图1所示,一种用于混响室内搅拌器的快速优化仿真设计方法,该仿真设计方法包含:
S1,根据搅拌器的设计要求参数,设计该搅拌器的外形结构,并确定要求优化设计的搅拌器参数。该步骤S1包含:
搅拌器的设计要求参数包含:搅拌器的旋转空间参数、搅拌器的回转直径、搅拌器的最大尺寸及搅拌器的表面材料。
根据上述搅拌器的设计要求参数设计该搅拌器的外形结构,并确定要求优化设计的搅拌器参数为:搅拌器的桨叶长度、桨叶宽度及相邻桨叶之间的夹角。
本实施例中,搅拌器的旋转空间应该足够大。搅拌器的回转直径不应小于最低使用频率对应波长的1/4倍,一般选用的回转直径范围为最低使用频率对应波长的1/3~1/2倍。搅拌器的最大尺寸不应小于对应混响室最小尺寸的75%。搅拌器的表面材料应为导电率高的金属材料。同时,本实施例要求搅拌器应具有非对称的外形结构。
如图4a-图4d所示,由于目前Z字形的搅拌器的桨叶是最有效且最简单的设计形式,因此本发明中,设计了四种不同桨叶长度、桨叶宽度以及不同相邻桨叶之间的夹角的搅拌器。
S2,选取多个等效入射源,在确保所有等效入射源能够模拟混响室内部随机电磁环境的前提下,确定每个等效入射源的随机入射方向和极化角。该步骤S2包含:
S2.1,设定第m个等效入射源的随机入射方向角为θm第m个等效入射源的随机极化角为ηm;其中,每个等效入射源的随机入射方向角、极化角均要求满足均匀分布统计特性的要求。
本实施例中,设定共有N个等效入射源能够模拟混响室内部随机电磁环境,其中0<m≤N。
本发明中,设置在球坐标系中的等效入射源的分布方位,(由于入射方向均指向球心,因此分布方位随机即代表入射方向也是随机的)、极化角也必须满足均匀分布统计特性。
本发明中,由于θm的随机性要受到球坐标中坐标变量取值范围的约束,因此在方位角θm极化角ηm随机生成时需要设置如以下步骤S2.2、S2.3的约束条件。
S2.2,要求第m个等效入射源的随机入射方向角θm的随机取值应在[0°,180°]之间;要求第m个等效入射源的随机入射方向角的随机取值应在[0°,360°]之间。该步骤S2.2包含:
S2.2.1,对于第m个等效入射源的随机入射方向角θm,在[0,1]内抽取随机数作为随机入射方向角θm的系数A,则θm=A·180°。
S2.2.2,对于第m个等效入射源的随机入射方向角的计算方法为:其中为随机取值,当的某一随机值使得与其对应的值落在[0°,360°]之间,则保留该值。
的随机取值中应特别值得注意的问题是,有时随机取值的结果有可能出现某一区域的取值样本过于密集,这种情况最有可能发生在高纬度即θm很小或θm很大的区域(即两极区域),对于发生这种情况的原因可作这样的解释:当θm很小或θm很大时,同球面上与之对应的纬度圈的周长2πrsinθm将很小,因而在取同样数量的样本个数的情况下,其样本点的密度便远大于其它纬度圈上的密度,因而若不加约束地按此种情况随机设置等效入射源,其结果将严重失实。因此,采用步骤S2.2.2对进行约束。
如图3a、图3b所示,随机等效源N=1000是,对应的等效源在球面上的分布情况。如图3a所示,若在步骤S2.2.2中不对进行约束,则球面两极附近区域等效源过于密集。如图图3b所示,经过在步骤S2.2.2中对进行约束处理后,两极附近区域等效源与其它区域随机分布保持一致。
S2.3,要求第m个等效入射源的随机极化角ηx的随机取值应在[0°,360°]之间。该步骤S2.3包含:
对于第m个等效入射源的随机极化角ηm,在[0,1]内抽取随机数作为随机入射方向角ηm的系数B,则ηm=B·360°。
如图2所示,在步骤S2中,第m个等效入射源入射方向、极化角如图中所示,方位角θm与随机等效源一一对应。由于等效源入射方向均指向球心,因此等效源分布方位随机即代表入射方向也是随机的。
S3,针对每个等效入射源进行仿真计算,获取以搅拌器为中心的球面上各个观察点的电场矢量分布计算结果。该步骤S3包含:
分别针对每个等效入射源进行仿真计算,其中每次仿真计算对象仅为搅拌器模型,从而获取以搅拌器为中心的球面上所有观察点的电场矢量分布计算结果。
本发明中,将突破传统的设计方法,通过随机等效入射源的引入,而仅仅需要针对搅拌器模型进行仿真,即类似“自由空间”下的仿真计算。一方面大大降低仿真计算的网格数,另一方面由于不存在像混响室这样的高谐振结构。因此其对计算机硬件资源需求大大降低、仿真时间呈几个数量级别的减少,从而极大提高了设计效率。
此外,在以搅拌器为中心的球面上,均匀选取电场观察点,获取所有观察点的电场矢量数据。
S4,对所有等效入射源的电场矢量分布计算结果进行矢量平均计算,采用混响室测试标准中的混响室均匀性统计方法,对搅拌器的所有观察点的电场矢量平均结果进行统计分析。该步骤S4包含:
S4.1,对所有等效入射源分别激励下的电场矢量计算结果进行矢量平均。
S4.2,采用混响室测试标准中的混响室均匀性统计方法,对搅拌器的所有观察点x、y、z三个极化方向进行计算获得三个极化方向的电场矢量标准偏差Dx、Dy及Dz,并计算上述三个极化方向的标准偏差Dtotal
本实施例中,采用IEC61000-4-21的国家标准中混响室均匀性统计方法,对所有观察点的电场矢量平均结果进行统计分析。
S5,修改步骤S1中的要求优化设计的搅拌器参数,针对每组优化设计的搅拌器参数重复执行步骤S3-S4,当完成获取在所有组不同优化设计的搅拌器参数时,计算对应的多组搅拌器参数的所有观察点的电场矢量平均结果的累加标准偏差。该步骤S5包含:
用户根据搅拌器的桨叶长度、桨叶宽度及相邻桨叶之间的夹角的不同,重复执行步骤S3-S4;获取多组不同优化设计的搅拌器参数情况下,对应的多个电场矢量平均结果的x、y、z三个极化方向的电场矢量标准偏差Dx、Dy、Dz及上述三个极化方向的标准偏差Dtotal
如图4a-图4d所示,图4a中搅拌器一的桨叶长度:0.75m,桨叶宽度:0.5m;图4b中搅拌器二的桨叶长度:1m,桨叶宽度:0.5m;图4c中搅拌器三的桨叶长度:0.75m,桨叶宽度:0.6m;图4d中搅拌器二的桨叶长度:0.75m,桨叶宽度:0.5m,并采用折叠结构。其中选用五十个随机等效入射源对上述四种搅拌器进行仿真计算。
S6,以电场矢量累加标准偏差越小、搅拌器的搅拌性能越好为评判原则,筛选获得搅拌器的最佳优化设计的搅拌器参数组形成的设计模型。该步骤S6包含:
S6.1,将步骤S5获取的每组优化设计的搅拌器参数对应的电场矢量平均结果的x、y、z三个极化方向的电场矢量标准偏差Dx、Dy、Dz及上述三个极化方向的标准偏差Dtotal进行如下处理计算,获取搅拌器电场矢量累加标准偏差D:
D=Dx+Dy+Dtotal
S6.2,通过比较步骤S5获取多组不同优化设计的搅拌器参数对应的搅拌器电场矢量累加标准偏差D,根据电场矢量累加标准偏差越小、搅拌器的搅拌性能越好为评判原则,筛选获得搅拌器的最佳优化设计的搅拌器参数组形成的设计模型。
图5为四种搅拌器累加标准偏差D统计分析示意图。从搅拌器的设计要求可知,搅拌器尺寸越大,其搅拌性能更好,从图中也可以明显看出该结论:搅拌器一和二相比较,搅拌器二的累计标准偏差相对较小,其性能更好,这主要是因为搅拌器二的旋转直径相对较大;搅拌器一和三相比较,除了个别频点外,搅拌器三的累计标准偏差也相对较小,其性能更好,这主要是因为搅拌器三的桨叶高度相对较大;搅拌器一和四相比较,低频段搅拌器四性能显著提高,而在高频段搅拌器四性能逐渐下降,比搅拌器一相对更差。这是因为搅拌器四是在搅拌器一的基础上进行了折叠处理,但是其四个桨叶变形比较对称,因此在某些频段可能无法显著提高搅拌性能。
为了进一步验证本发明中提出的搅拌器快速优化仿真设计方法,将搅拌器一和搅拌器二分别放置在尺寸为4.9m(长)×3m(宽)×2.3m(高)混响室中,按照传统的仿真分析方法进一步确定两者的优劣性。
如图6a、图6b所示,搅拌器二的搅拌性能要优于搅拌器一,因为在搅拌器二的搅拌作用下,场均匀性统计的标准偏差明显小于搅拌器一对应的标准偏差。因此搅拌器二的搅拌性能优于搅拌器一,这与步骤S6中得到的结论相一致。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。

Claims (9)

1.一种用于混响室内搅拌器的快速优化仿真设计方法,其特征在于,该仿真设计方法包含:
S1,根据搅拌器的设计要求参数,设计该搅拌器的外形结构,并确定要求优化设计的搅拌器参数;
S2,选取多个等效入射源,在确保所有等效入射源能够模拟混响室内部随机电磁环境的前提下,确定每个所述等效入射源的随机入射方向和极化角;
S3,针对每个所述等效入射源进行仿真计算,获取以所述搅拌器为中心的球面上各个观察点的电场矢量分布计算结果;
S4,对所有所述等效入射源的电场矢量分布计算结果进行矢量平均计算,采用混响室测试标准中的混响室均匀性统计方法,对所述搅拌器的所有观察点的电场矢量平均结果进行统计分析;
S5,修改所述步骤S1中的所述要求优化设计的搅拌器参数,并重复执行所述步骤S3-S4,获取在多组不同优化设计的搅拌器参数下,对应的多组搅拌器参数的所有观察点的电场矢量平均结果的累加标准偏差;
S6,以电场矢量累加标准偏差越小、所述搅拌器的搅拌性能越好为评判原则,筛选获得所述搅拌器的最佳优化设计的搅拌器参数组形成的设计模型。
2.如权利要求1所述的用于混响室内搅拌器的快速优化仿真设计方法,其特征在于,所述步骤S1包含:
所述搅拌器的设计要求参数包含:搅拌器的旋转空间参数、搅拌器的回转直径、搅拌器的最大尺寸及搅拌器的表面材料;
根据上述搅拌器的设计要求参数设计该搅拌器的外形结构,并确定要求优化设计的搅拌器参数为:所述搅拌器的桨叶长度、桨叶宽度及相邻桨叶之间的夹角。
3.如权利要求1所述的用于混响室内搅拌器的快速优化仿真设计方法,其特征在于,所述步骤S2包含:
S2.1,设定第m个所述等效入射源的随机入射方向角为θm第m个所述等效入射源的随机极化角为ηm;其中,每个所述等效入射源的随机入射方向角、极化角均要求满足均匀分布统计特性的要求;
S2.2,要求第m个所述等效入射源的随机入射方向角θm的随机取值应在[0°,180°]之间;要求第m个所述等效入射源的随机入射方向角的随机取值应在[0°,360°]之间;
S2.3,要求第m个所述等效入射源的随机极化角ηx的随机取值应在[0°,360°]之间。
4.如权利要求3所述的用于混响室内搅拌器的快速优化仿真设计方法,其特征在于,所述步骤S2.2包含:
S2.2.1,对于第m个所述等效入射源的随机入射方向角θm,在[0,1]内抽取随机数作为随机入射方向角θm的系数A,则θm=A·180°;
S2.2.2,对于第m个所述等效入射源的随机入射方向角的计算方法为:其中为随机取值,当的某一随机值使得与其对应的值落在[0°,360°]之间,则保留该值。
5.如权利要求3所述的用于混响室内搅拌器的快速优化仿真设计方法,其特征在于,所述步骤S2.3包含:
对于第m个所述等效入射源的随机极化角ηm,在[0,1]内抽取随机数作为随机入射方向角ηm的系数B,则ηm=B·360°。
6.如权利要求1所述的用于混响室内搅拌器的快速优化仿真设计方法,其特征在于,所述步骤S3包含:
分别针对每个所述等效入射源进行仿真计算,其中每次仿真计算对象仅为搅拌器模型,从而获取以所述搅拌器为中心的球面上所有观察点的电场矢量分布计算结果。
7.如权利要求2所述的用于混响室内搅拌器的快速优化仿真设计方法,其特征在于,所述步骤S4包含:
S4.1,对所有所述等效入射源分别激励下的电场矢量计算结果进行矢量平均;
S4.2,采用混响室测试标准中的混响室均匀性统计方法,对所述搅拌器的所有观察点x、y、z三个极化方向进行计算获得三个极化方向的电场矢量标准偏差Dx、Dy及Dz,并计算上述三个极化方向的标准偏差Dtotal
8.如权利要求7所述的用于混响室内搅拌器的快速优化仿真设计方法,其特征在于,所述步骤S5包含:
用户根据所述搅拌器的桨叶长度、桨叶宽度及相邻桨叶之间的夹角的不同,重复执行所述步骤S3-S4;获取多组不同优化设计的搅拌器参数情况下,对应的多个电场矢量平均结果的x、y、z三个极化方向的电场矢量标准偏差Dx、Dy、Dz及上述三个极化方向的标准偏差Dtotal
9.如权利要求8所述的用于混响室内搅拌器的快速优化仿真设计方法,其特征在于,所述步骤S6包含:
S6.1,将所述步骤S5获取的每组优化设计的搅拌器参数对应的电场矢量平均结果的x、y、z三个极化方向的电场矢量标准偏差Dx、Dy、Dz及上述三个极化方向的标准偏差Dtotal进行如下处理计算,获取搅拌器电场矢量累加标准偏差D:
D=Dx+Dy+Dtotal
S6.2,通过比较所述步骤S5获取多组不同优化设计的搅拌器参数对应的搅拌器电场矢量累加标准偏差D,根据电场矢量累加标准偏差越小、所述搅拌器的搅拌性能越好为评判原则,筛选获得所述搅拌器的最佳优化设计的搅拌器参数组形成的设计模型。
CN201510802353.8A 2015-11-19 2015-11-19 一种用于混响室内搅拌器的快速优化仿真设计方法 Pending CN105302992A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510802353.8A CN105302992A (zh) 2015-11-19 2015-11-19 一种用于混响室内搅拌器的快速优化仿真设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510802353.8A CN105302992A (zh) 2015-11-19 2015-11-19 一种用于混响室内搅拌器的快速优化仿真设计方法

Publications (1)

Publication Number Publication Date
CN105302992A true CN105302992A (zh) 2016-02-03

Family

ID=55200260

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510802353.8A Pending CN105302992A (zh) 2015-11-19 2015-11-19 一种用于混响室内搅拌器的快速优化仿真设计方法

Country Status (1)

Country Link
CN (1) CN105302992A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106294998A (zh) * 2016-08-09 2017-01-04 中国计量科学研究院 一种电磁混响室腔体的尺寸优化方法及装置
CN108710035A (zh) * 2018-02-28 2018-10-26 北京无线电计量测试研究所 一种改善混响室均匀性的方法和装置
CN111337756A (zh) * 2020-04-24 2020-06-26 中国人民解放军63892部队 一种评价脉冲激励混响室测试区域电场均匀性方法
CN111337786A (zh) * 2020-04-24 2020-06-26 中国人民解放军63892部队 一种提高混响室场均匀性的线搅拌器设计方法
CN114113810A (zh) * 2021-11-24 2022-03-01 中国人民解放军陆军工程大学 边界形变混响室空间电场统计均匀性测试方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120268340A1 (en) * 2009-09-16 2012-10-25 Agence Spatiale Europeenne Aperiodic and Non-Planar Array of Electromagnetic Scatterers, and Reflectarray Antenna Comprising the Same
CN104133126A (zh) * 2014-07-01 2014-11-05 陕西海泰电子有限责任公司 一种源搅拌电磁混响室及其搅拌方法
CN104407248A (zh) * 2014-11-27 2015-03-11 中国人民解放军军械工程学院 一种基于混响室平台的电子系统电磁环境效应试验方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120268340A1 (en) * 2009-09-16 2012-10-25 Agence Spatiale Europeenne Aperiodic and Non-Planar Array of Electromagnetic Scatterers, and Reflectarray Antenna Comprising the Same
CN104133126A (zh) * 2014-07-01 2014-11-05 陕西海泰电子有限责任公司 一种源搅拌电磁混响室及其搅拌方法
CN104407248A (zh) * 2014-11-27 2015-03-11 中国人民解放军军械工程学院 一种基于混响室平台的电子系统电磁环境效应试验方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
崔耀中 等: "混响室关键配置对场分布影响分析", 《电波科学学报》 *
汪宗福 等: "基于双因素方差分析的混响室搅拌器优化研究", 《宇航计测技术》 *
黄华 等: "时域有限差分法分析混响室中场的均匀性", 《电波科学学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106294998A (zh) * 2016-08-09 2017-01-04 中国计量科学研究院 一种电磁混响室腔体的尺寸优化方法及装置
CN106294998B (zh) * 2016-08-09 2019-09-17 中国计量科学研究院 一种电磁混响室腔体的尺寸优化方法及装置
CN108710035A (zh) * 2018-02-28 2018-10-26 北京无线电计量测试研究所 一种改善混响室均匀性的方法和装置
CN111337756A (zh) * 2020-04-24 2020-06-26 中国人民解放军63892部队 一种评价脉冲激励混响室测试区域电场均匀性方法
CN111337786A (zh) * 2020-04-24 2020-06-26 中国人民解放军63892部队 一种提高混响室场均匀性的线搅拌器设计方法
CN111337786B (zh) * 2020-04-24 2023-02-24 中国人民解放军63892部队 一种提高混响室场均匀性的线搅拌器设计方法
CN114113810A (zh) * 2021-11-24 2022-03-01 中国人民解放军陆军工程大学 边界形变混响室空间电场统计均匀性测试方法及装置

Similar Documents

Publication Publication Date Title
CN105302992A (zh) 一种用于混响室内搅拌器的快速优化仿真设计方法
CN106099393A (zh) 基于遗传算法的子阵稀布与阵元稀疏天线阵列综合布阵方法
CN105243280B (zh) 基于cpu与gpu混合异步并行方式的时域物理光学计算方法
CN108959772B (zh) 一种大规模有限周期阵列结构特征模式分析方法
CN103425816B (zh) 快速获取金属旋转对称体电磁散射特性的矩阵抽取方法
CN106252873A (zh) 一种共形承载天线远场功率方向图的区间分析方法
CN107515956A (zh) 一种基于hfss单元法的大型有限平面阵列分析方法
CN110472356B (zh) 一种电磁波多入射角度下的复合材料电磁参数计算方法
Sherman et al. Structure factors of the kagome-lattice heisenberg antiferromagnets at finite temperatures
CN110362877B (zh) 不确定因素的环境电磁散射特性分析方法
CN112285435B (zh) 一种大功率磁场辐射源的等效模拟方法
CN107305536B (zh) 混合阶时域不连续伽略金方法
Xu et al. Improved synthetic basis functions method for nonperiodic scaling structures with arbitrary spatial attitudes
CN106503349B (zh) 一种类周期结构目标电磁散射特性快速计算方法
Zhou et al. An adaptive proper orthogonal decomposition method for evaluating variability bounds of antenna responses
CN113343423A (zh) 基于强度空间变异性的随机裂隙网络生成方法
Chen et al. Parallel DGTD method for transient electromagnetic problems
CN103136401B (zh) 一种人造微结构的试验设计表获取方法及装置
Kim et al. Performance evaluation of a mode-stirred reverberation chamber using the finite difference time domain (FDTD) simulation
Wang et al. Edge diffraction in NURBS-UTD method
CN111339650B (zh) 正交维数缩减的电磁基函数构造方法及装置
Hara et al. Numerical derivation of design guidelines for tightness and shaking amplitude of vibrating intrinsic reverberation chamber by method of moment
Clegg et al. Optimisation of stirrer designs in a mode stirred chamber using TLM
Zhang et al. Fast Wideband RCS Analysis of Large-Scale Finite Periodic Structures Using MB-CBF Method and Rational Interpolation
CN111832157B (zh) 基于子全域基函数方法的大规模准周期结构电磁散射特性分析方法

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

Application publication date: 20160203

RJ01 Rejection of invention patent application after publication