CN113295941B - 一种基于干扰概率的混响室环境电磁辐射敏感度测试方法 - Google Patents
一种基于干扰概率的混响室环境电磁辐射敏感度测试方法 Download PDFInfo
- Publication number
- CN113295941B CN113295941B CN202110450645.5A CN202110450645A CN113295941B CN 113295941 B CN113295941 B CN 113295941B CN 202110450645 A CN202110450645 A CN 202110450645A CN 113295941 B CN113295941 B CN 113295941B
- Authority
- CN
- China
- Prior art keywords
- eut
- electromagnetic radiation
- test
- effect
- cable
- 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
- 238000012360 testing method Methods 0.000 title claims abstract description 208
- 230000005670 electromagnetic radiation Effects 0.000 title claims abstract description 110
- 230000035945 sensitivity Effects 0.000 title claims abstract description 45
- 230000000694 effects Effects 0.000 claims abstract description 107
- 238000000034 method Methods 0.000 claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 230000007613 environmental effect Effects 0.000 claims abstract description 14
- 230000008878 coupling Effects 0.000 claims description 65
- 238000010168 coupling process Methods 0.000 claims description 65
- 238000005859 coupling reaction Methods 0.000 claims description 65
- 230000005684 electric field Effects 0.000 claims description 31
- 238000003756 stirring Methods 0.000 claims description 26
- 230000005855 radiation Effects 0.000 claims description 25
- 238000005259 measurement Methods 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 12
- 230000000903 blocking effect Effects 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 4
- 230000002452 interceptive effect Effects 0.000 claims description 3
- 230000009916 joint effect Effects 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 230000004064 dysfunction Effects 0.000 claims description 2
- 230000007257 malfunction Effects 0.000 claims description 2
- 239000002904 solvent Substances 0.000 claims 1
- 230000001965 increasing effect Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 230000005672 electromagnetic field Effects 0.000 description 8
- 238000010998 test method Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 238000002347 injection Methods 0.000 description 5
- 239000007924 injection Substances 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 230000010287 polarization Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000013019 agitation Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- YBMRDBCBODYGJE-UHFFFAOYSA-N germanium dioxide Chemical compound O=[Ge]=O YBMRDBCBODYGJE-UHFFFAOYSA-N 0.000 description 2
- 230000009191 jumping Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000004132 cross linking Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000005669 field effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012113 quantitative test Methods 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/001—Measuring interference from external sources to, or emission from, the device under test, e.g. EMC, EMI, EMP or ESD testing
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Testing Electric Properties And Detecting Electric Faults (AREA)
Abstract
本发明公开了一种基于干扰概率的混响室环境电磁辐射敏感度测试方法，包括以下步骤：步骤一，规划搅拌器步进规则；步骤二，混响室馈入功率预调节；步骤三，EUT电磁辐射效应测试；步骤四，确定EUT等效天线方向性系数最大值；步骤五，EUT临界干扰场强的定量计算；步骤六，测定EUT的临界干扰场强曲线；本发明的基于干扰概率的混响室环境电磁辐射敏感度测试方法，结合工程实际，形成简便易行、便于第三方应用的测试方法与测试流程，提高了测试准确度和测试效率。
Description
技术领域
本发明涉及一种基于干扰概率的混响室环境电磁辐射敏感度测试方法，属于电磁技术领域。
背景技术
随着信息技术的高速发展，电子信息设备辐射功率不断攀升，电磁环境日趋恶劣，电磁干扰现象层出不穷；电磁辐射敏感性测试在评价设备抗电磁干扰能力、揭示电磁辐射敏感通道与作用机理、对症下药进行电磁防护加固等方面发挥了越来越大的作用。目前，传统的全电平电磁辐射效应测试方法已经成熟，但仍难以解决大尺寸受试设备、互联系统的强场电磁辐射效应测试难题，发展电磁环境高效模拟技术、电磁辐射效应等效测试方法，解决等效测试方法与传统全电平电磁辐照法测试结果的一致性问题是该领域的研究重点和热点。为此，相继发展了大电流注入、直接电流注入等共模电流注入测试方法和差模注入-线性外推、接地电流注入等差模电流注入测试方法，一定程度上解决了线缆束、天线与同轴线缆、地线耦合导致的强场电磁辐射效应测试技术难题，但难以全面反映受试设备的多径电磁辐射耦合交联问题，导致测试结果与全电平电磁辐照法的相关性变差。
混响室是一种内部具有搅拌器的高品质因数屏蔽腔体，利用金属表面的高反射特性，使电磁波在腔体内部反复反射，达到电磁能量高效利用、以较低输入功率激发强场电磁环境的目的；搅拌器通过不断改变电磁场的边界条件，在测试区域内形成各向同性、随机极化、统计均匀的电磁环境，特别适合进行设备强场电磁辐射效应测试。但是，混响室测试环境中不同位置电磁场的传播方向、极化方向随机变化且不可控，场强大小各处不同，难以直接用场强计度量受试设备所处测试环境的场强大小。因此，在混响室环境中对受试系统进行电磁辐射敏感性测试需要解决以下两个关键问题：一是临界干扰场强如何度量？二是敏感性测试结果与均匀场环境测试结果的相关性、等价性。经过多年的系统研究，基本解决了上述理论基础和关键技术，但如何结合工程实际，形成简便易行、便于第三方应用的测试方法与测试流程，提高测试准确度，尚有许多问题需要破解。
国际电工委员会标准IEC 61000-4-21给出了混响室条件下的辐射敏感度测试方法：通过测量混响室的空载归一化场强和被测设备(Equipment Under Test，EUT)加载前后的参考天线校准系数，确定混响室的加载系数，进而给出加载后的归一化场强；根据混响室的输入功率、加载归一化场强，计算EUT加载后的混响室最大场强；基于EUT受扰与场强最大值相关，通过调节混响室的输入功率来控制测试区域场强的最大值，完成EUT辐射敏感度的测试。由于混响室内的归一化场强基于其立方体测试区域8个顶点处电场强度的最大值统计得出，混响室内电场统计均匀、各向同性的前提是搅拌位置数和每一位置的场分布模式数足够多，而按标准进行校准和效应测试时搅拌位置数仅有12～50，难以完全满足电场分布统计均匀的要求。在搅拌位置数有限的情况下，不同位置的电场强度最大值存在不确定性且随着搅拌步数N变化，校准位置出现的电场强度最大值，EUT所在位置未必能够经历。如图1所示，随着搅拌位置数的增加，电场分量的期望值逐步提高，但分布范围逐步缩小，搅拌位置数高达100时，电场最大值的变化范围仍高达7dB以上；搅拌器步进一周，不同位置处的电场最大值差别很大；另外，混响室不同位置出现场强最大值时的极化方向不同，最大场强值难以与EUT的敏感方向对应，从而导致测试结果存在较大差异。
上述两方面问题导致混响室环境电磁辐射敏感度测试结果重复性差、且与均匀场环境测试结果相关性差，制约了标准的推广应用；电磁辐射敏感性作为受试设备的固有属性，不应随测试场地的改变而出现较大差异。为此，本发明提出了一种基于干扰概率的混响室环境电磁辐射敏感度测试方法，并通过试验验证了方法的准确性。
发明内容
为解决上述问题，本发明提出了一种基于干扰概率的混响室环境电磁辐射敏感度测试方法，结合工程实际，形成简便易行、便于第三方应用的测试方法与测试流程，提高了测试准确度和测试效率。
本发明的基于干扰概率的混响室环境电磁辐射敏感度测试方法，包括以下步骤：
步骤一，规划搅拌器步进规则，
根据混响室的搅拌器数量，规划搅拌器步进规则，使搅拌器步进一周的混响室搅拌位置数达到30左右，最好在30～50之间且使各个搅拌器等角度步进的位置数大致相等；
步骤二，混响室馈入功率预调节，
根据EUT测试要求，选定测试频率，设定混响室激励天线的初始馈入功率，对EUT进行电磁辐射效应快速测试，监视或检测确定EUT的效应类型；参照变步长升降法调节混响室激励天线的馈入功率，在搅拌器步进一周的过程中，使EUT出现某一效应或多种效应的概率分别达到20％～60％；
步骤三，EUT电磁辐射效应测试，
固定混响室激励天线的馈入功率，按搅拌器步进规则步进搅拌器一周，记录EUT出现每种效应的概率估计值
测试过程中，在测试区域远离EUT的某一固定位置测试搅拌器处于每一搅拌位置时的电场强度有效值E0n，按以下公式计算混响室测试环境电场强度有效均值E0：
步骤四，确定EUT等效天线方向性系数最大值，
根据EUT出现的电磁辐射效应及EUT、线缆电尺寸，按EUT电磁辐射耦合等效天线方向性系数的选取原则，分别采用用频设备天线系数或按数学拟合估值公式计算，确定EUT等效天线方向性系数最大值的估计值；
步骤五，EUT临界干扰场强的定量计算，
根据EUT出现每种效应的概率估计值
混响室测试环境电场强度有效均值E0和每种效应对应的EUT等效天线方向性系数最大值的估计值
利用以下公式：
计算EUT在测试频点的临界干扰场强估值
则临界干扰场强估值如以下公式所示：
步骤六，测定EUT的临界干扰场强曲线，
改变测试频率，重复上述测试步骤二至步骤五，测试确定不同频率对应的EUT临界干扰场强，给出不同效应对应的EUT临界干扰场强曲线。
进一步地，所述步骤四的具体操作步骤如下：
不同电磁辐射效应下等效天线方向性系数的选取和计算，
(1)当EUT的电磁辐射效应为用频设备阻塞、虚警与损伤效应，则采用混响室测试用频设备的电磁辐射敏感度时，直接用天线系数代替Dmax即可；
(2)当EUT的电磁辐射效应为强场电磁辐射效应，则采用混响室测试该类设备的电磁辐射敏感度时，按EUT接地线或共地线缆的电长度，采用线缆共模耦合等效天线方向性系数代替Dmax，计算受试系统的临界辐射干扰场强Es；当线缆长度L＜3.5λ时，采用线缆共模耦合等效天线方向性系数对数中值的平滑线代替Dmax，当L≥3.5λ时Dmax取整不再变化；为方便计算，给出Dmax估值的数学拟合估值公式：
(3)当EUT的电磁辐射效应为测控设备的扰动效应，采用混响室进行效应试验时，若线缆差模耦合和孔缝耦合单独作用，则低频段以线缆耦合为主，高频段以孔缝耦合为主；当受试测控设备在线缆差模耦合和孔缝耦合联合作用下出现扰动效应时，则等效天线方向性系数最大值Dmax用统一的拟合公式来估值：
当线缆长度L＜3.5λ时，a＝L；否则a为包围EUT的最小球体半径。
再进一步地，所述强场电磁辐射效应包括死机、重启和功能紊乱效应；所述功能紊乱效应包括显示乱码和功能失效。
再进一步地，当EUT的电磁辐射效应为强场电磁辐射效应，采用混响室测试该类设备的电磁辐射敏感度时，若EUT具有多根不同长度的接地线或共地线缆时，则根据测试频率选择与半波长较接近的线缆长度计算Dmax的估值。
再进一步地，当EUT的电磁辐射效应为测控设备的扰动效应，若EUT采用屏蔽线缆，其差模干扰信号源于线缆屏蔽层的共模感应电流，当线缆长度L＜3.5λ时，采用公式(4)计算Dmax的估值；否则，按公式(5)进行计算。
再进一步地，所述强场电磁辐射效应还包括射频前端硬损伤效应。
本发明与现有技术相比较，本发明的基于干扰概率的混响室环境电磁辐射敏感度测试方法，结合工程实际，形成简便易行、便于第三方应用的测试方法与测试流程，提高了测试准确度和测试效率。
附图说明
图1是混响室中固定位置某一电场分量最大值相对于测试空间电场振幅平均值的标准概率密度分布曲线。
图2是孔缝耦合Dmax的取值范围示意图。
图3是线缆差模耦合Dmax的取值范围示意图。
图4是线缆共模耦合Dmax的取值范围示意图。
图5是1m长线缆的感应电流幅频曲线。
图6是干扰概率
和搅拌次数N对ES测试误差的影响示意图。
图7是本发明的受试电台结构示意图。
图8是本发明的屏蔽室试验布置示意图。
图9是本发明的混响室试验布置示意图。
图10是本发明的混响室试验和屏蔽室试验的试验流程示意图。
图11是本发明的电压变化数值与场强之间的关系示意图；
其中，图11(a)是本发明的0.55GHz对电台的某一方向进行辐照时的测试结果示意图；图11(b)是本发明的3.2GHz对电台的某一方向进行辐照时的测试结果示意图。
图12是本发明的受试电台临界干扰场强测试结果对比示意图。
具体实施方式
本发明的基于干扰概率的混响室环境电磁辐射敏感度测试方法，研究过程如下：
首先，针对混响室环境电磁辐射敏感度定量测试的技术需求，分析了现有测试方法的局限性，指出了存在的问题，从混响室电磁场统计分布规律和电磁辐射效应机理出发，建立了EUT临界干扰场强与混响室测试环境电场强度有效均值、EUT受扰概率和等效天线方向性系数最大值之间的关系。
其次，统筹考虑效应类型、等效天线耦合效率等因素，从用频设备干扰与损伤效应、强场电磁辐射效应和测控设备扰动效应三方面，给出了EUT等效天线方向性系数的选取原则和拟合公式；具体如下：
由于混响室中电磁场分布特性与传统均匀场测试环境有显著的差异，要保证在混响室中获取的EUT电磁辐射敏感度测试结果与均匀场环境一致，需要建立混响室与均匀场之间的关联关系；
第一，推导临界辐射干扰场强计算公式，
现有技术中Amador以混响室场强直角分量的幅值服从瑞利分布、而场强直角分量的幅值大于等于受试设备的临界辐射干扰场强Es时EUT受到干扰为依据，从统计分析入手，提出了通过EUT的干扰概率来计算混响室测试环境EUT临界辐射干扰场强的方法，建立了混响室环境EUT临界干扰场强的计算公式：
式中P为EUT在搅拌器步进一周的N个位置中受干扰的概率，σ可由混响室测试区域某一固定位置在N个搅拌状态的场强直角分量Exn来确定；具体地，根据混响室中场强直角分量服从瑞利分布，设第n个搅拌位置混响室区域内一点场强x方向分量为Exn，σ可以通过最大似然估计法进行估计：
上述方法通过干扰概率来计算EUT的临界辐射干扰场强，为混响室环境EUT临界辐射干扰场强测试提供了一种新颖的思路，但其推导过程中未能考虑EUT电磁耦合敏感方向与电场极化方向之间的差异，公式(6)不具有普适性；
由于混响室中的电磁环境比较复杂，根据EUT外部的电场很难对EUT所面临的辐射场强进行准确度量；按照现有技术中Amador处理问题的思路，假定EUT场路耦合的等效天线耦合特性与电磁场分布无关、EUT的敏感元件接收功率/电压达到其临界干扰功率/电压时EUT即可受到干扰；在此基础上，将EUT的敏感元件等效为接收天线负载，令混响室与均匀场中敏感元件的临界干扰功率/电压相等，则现有技术中(胡德洲等人提出的混响室与均匀场辐射敏感度测试相关性研究)推导了混响室条件下设备临界辐射干扰场强瞬时值的普适计算公式：
式中P、σ的含义与式(6)相同，Dmax是EUT场路耦合等效天线的方向性系数最大值；
若临界干扰场强用有效值ESE表示，测试区域某一固定位置的电场强度在混响室搅拌器步进一周的N个位置的有效均值为E0(第n个搅拌位置时的有效值为E0n)，则：
若搅拌器步进一周的N个位置中EUT受扰NS次，干扰概率P可用其估计值代替：
E0通过测量确定后，只要知道EUT场路耦合等效天线的方向性系数最大值Dmax，即可根据式(3)计算EUT的临界干扰场强；
第二，等效天线方向性系数最大值的近似估计，
EUT的等效天线方向性系数最大值Dmax与其有效电磁耦合途径紧密相关，耦合途径不同，Dmax的变化规律不同；天线的方向性系数容易测量，现有技术中(胡德洲等人提出的混响室条件下临界辐射干扰场强测试的影响因素分析)利用MATLAB编程计算的方法，系统研究了孔缝、线缆耦合的等效天线方向性系数统计特性，给出了其取值范围的变化规律；图2所示为EUT电尺寸ka(a为能够包围EUT的最小球体半径)对孔缝耦合Dmax取值范围的影响规律，图中同时给出了95％置信度对应的置信区间上、下限及其对数中值；除非EUT的电尺寸太小(此时孔缝电尺寸太小，难以成为EUT受扰的主要电磁辐射耦合通道)，以其对数中值作为Dmax的估值
时，误差可以控制在3dB之内；
线缆耦合分为差模耦合与共模耦合两种情况，Dmax的分布规律大不相同，线缆差模耦合时，95％置信度对应的Dmax的置信区间及其对数中值如图3所示，在线缆长度大于半个波长后，Dmax几乎随线缆长度增加而周期性变化；而线缆共模耦合时，Dmax随线缆长度增加而波动性增加，图4同时给出了95％置信度对应的线缆共模耦合Dmax置信区间的MATLAB编程计算结果和FEKO软件仿真结果，从图中可以看出：FEKO软件仿真结果相当于MATLAB编程计算给出的Dmax线性比例中值；
第三，根据等效天线方向性系数的选取原则，给出了EUT等效天线方向性系数拟合公式，
EUT的电磁辐射效应不同，如阻塞、虚警、死机、重启、功能紊乱等，其有效电磁辐射耦合通道随之改变；同时，不同耦合通道的电磁辐射耦合效率大不相同，确定EUT的等效方向性系数必须统筹考虑效应类型、耦合效率和线缆、EUT电尺寸等多方面因素；
(1)针对用频设备阻塞、虚警与损伤效应，
阻塞效应源于电路的非线性，干扰信号增强导致有用信号增益降低所致；虚警效应源于干扰信号对用频设备的直接作用，将干扰信号作为有用信号进行处理所致；损伤效应源于高电压、大功率信号输入，超过射频前端的损伤阈值所致；尽管天线、线缆、孔缝耦合均能够导致上述3种效应，但由于天线耦合效率远高于线缆和孔缝，用频设备普遍具备接收天线，接收天线才是用频设备出现上述多种效应的有效电磁辐射耦合通道；因此，采用混响室测试用频设备的电磁辐射敏感度时，直接用天线系数代替Dmax即可；对射频前端之外的电磁辐射损伤效应，归入下述强场电磁辐射效应进行处理；
(2)针对强场电磁辐射效应，
强场电磁辐射效应主要包括：死机、重启、功能紊乱(显示乱码、功能失效)和硬损伤等4种类型；除了前述射频前端硬损伤外，导致其它强场电磁辐射效应的根本原因一般为地电位波动，其电磁辐射耦合的主要通道是接地线或输入、输出共地线缆；采用混响室测试该类设备的电磁辐射敏感度时，按EUT接地线或共地线缆的电长度，采用线缆共模耦合等效天线方向性系数代替Dmax，计算受试系统的临界辐射干扰场强Es；
从图4仿真结果可见，线缆共模耦合等效天线方向性系数随线缆电长度的增加波动增加，其原因是线缆长度大于半波长时，线缆中不同线元的电磁辐射在空间不同位置同相或反相迭加，增强了电磁辐射的方向性；当线缆长度较长时，一方面线缆弯曲导致其等效天线方向性系数下降；另一方面，线缆越长其电磁辐射接收能力越弱，图5所示为平面波均匀辐照下1m长同轴电缆外导体感应电流幅频曲线的计算结果，当线缆长度接近半波长的奇数倍时感应电流出现极大值，辐射频率越高，归一化电流幅值越低，线缆长度大于3.5个波长时，其电磁辐射选频特性已不明显；因此，综合考虑线缆弯曲导致等效天线方向性系数下降和线缆长度L大于3.5λ时线缆选频特性已不明显两方面因素，当线缆长度L＜3.5λ时，采用线缆共模耦合等效天线方向性系数对数中值的平滑线代替Dmax，当L≥3.5λ时Dmax取整不再变化；为方便计算，给出Dmax估值的数学拟合估值公式：
当EUT具有多根不同长度的接地线或共地线缆时，根据测试频率选择与半波长较接近的线缆长度计算Dmax的估值；
(3)针对测控设备的扰动效应，
测控设备的扰动效应源于干扰信号与有用信号的迭加，属于差模干扰范畴，采用混响室进行效应试验时需要考虑线缆差模耦合、共模耦合和壳体孔缝耦合哪个效率更高；图5说明线缆耦合在低频段效率较高，而壳体孔缝耦合低频效率很低，频率越高耦合效率越高；因此，导致受试测控设备出现扰动效应时，低频段以线缆耦合为主，高频段以孔缝耦合为主；
比较图2、图3孔缝耦合、线缆差模耦合Dmax取值范围的对数中值可以发现：在低频段线缆差模耦合与孔缝耦合的Dmax取值相差较小，高频段孔缝耦合的Dmax明显高于线缆差模耦合的Dmax；因此，当受试测控设备在线缆差模耦合、孔缝耦合联合作用下出现扰动效应时，等效天线方向性系数最大值Dmax可以用统一的拟合公式来估值：
当线缆长度L＜3.5λ时，a＝L；否则a为包围EUT的最小球体半径；
若EUT采用屏蔽线缆，其差模干扰信号源于线缆屏蔽层的共模感应电流，当线缆长度L＜3.5λ时，采用式(4)计算Dmax的估值；否则，按式(5)进行计算；
第四，搅拌次数和干扰概率的选取，
混响室内电磁环境统计均匀、各向同性是电磁场多模式、多状态统计分析的结果，混响室模式数密度再高，离开搅拌器多状态的贡献，其电磁环境也难以达到统计均匀；因此，为实现混响室与均匀场环境电磁辐射敏感度测试的等效性，在混响室环境进行电磁辐射效应测试时，搅拌次数必须足够多，出现干扰的概率也需适当，不能过大或过小，以降低统计误差；
针对上述问题，现有技术中(胡德洲等人提出的混响室条件下临界辐射干扰场强测试的影响因素分析)通过系统的仿真分析和试验验证，结果表明：混响室的搅拌次数、测试时的干扰概率是影响测试结果准确性的主要因素；测试误差随搅拌次数的增加单调降低，搅拌次数小于20时，测试误差随搅拌次数增加迅速降低；搅拌次数在30左右时即可对混响室内电场强度有效均值E0进行较为准确的估计，由此带来的测试误差控制在10％之内；搅拌次数再增加，测试误差的下降速率减缓，对测试误差降低的贡献不再显著；无论混响室搅拌次数大小，测试过程中受试系统的干扰概率过高或过低都会导致临界干扰场强测试误差的增加，干扰概率取10％～60％比较适宜，与均匀场测试相比，误差能够控制在3dB左右，如图6所示；干扰概率再增加，不仅会造成混响室电磁激励能量的浪费，同时还会导致测试误差增加；
综上所述，只要根据EUT的效应类型合理选择其等效天线方向性系数最大值Dmax，效应测试过程中使混响室的搅拌位置数达到30左右、调节混响室馈入功率使EUT的干扰概率达到10％～60％，在测试区域某一固定位置测试电场强度的有效值，按式(1)计算混响室测试环境的有效场强均值E0，最后按式(2)计算EUT的临界干扰场强有效值ESE，即可得到与均匀场测试环境等效的测试辐射敏感度。
接着，进行电磁辐射敏感度混响室测试，由前述可知：在混响室环境测试EUT的电磁辐射敏感度时，若采用基于干扰概率的方法计算EUT的临界干扰场强，则无需事先确定混响室加载前后的归一化场强，只要在测试过程中测试某一固定位置的辐射场强即可；因此，只要通过定期校准保证混响室电磁环境的各向同性、统计均匀，在每次对EUT进行电磁辐射敏感度测试时，不用单独进行加载前后的混响室校准，大幅度降低了测试工作量，测试流程具体如下：
本发明的基于干扰概率的混响室环境电磁辐射敏感度测试方法，包括以下步骤：
步骤一，规划搅拌器步进规则，
根据混响室的搅拌器数量，规划搅拌器步进规则，使搅拌器步进一周的混响室搅拌位置数达到30左右，最好在30～50之间且使各个搅拌器等角度步进的位置数大致相等；
步骤二，混响室馈入功率预调节，
根据EUT测试要求，选定测试频率，设定混响室激励天线的初始馈入功率，对EUT进行电磁辐射效应快速测试，监视或检测确定EUT的效应类型；参照变步长升降法调节混响室激励天线的馈入功率，在搅拌器步进一周的过程中，使EUT出现某一效应或多种效应的概率分别达到20％～60％；
步骤三，EUT电磁辐射效应测试，
固定混响室激励天线的馈入功率，按搅拌器步进规则步进搅拌器一周，记录EUT出现每种效应的概率估计值
测试过程中，在测试区域远离EUT的某一固定位置测试搅拌器处于每一搅拌位置时的电场强度有效值E0n，按以下公式计算混响室测试环境电场强度有效均值E0：
步骤四，确定EUT等效天线方向性系数最大值，
根据EUT出现的电磁辐射效应及EUT、线缆电尺寸，按EUT电磁辐射耦合等效天线方向性系数的选取原则，分别采用用频设备天线系数或按数学拟合估值公式式(4)和式(5)计算，确定EUT等效天线方向性系数最大值的估计值；
步骤五，EUT临界干扰场强的定量计算，
根据EUT出现每种效应的概率估计值
混响室测试环境电场强度有效均值E0和每种效应对应的EUT等效天线方向性系数最大值的估计值
利用以下公式：
计算EUT在测试频点的临界干扰场强估值
则临界干扰场强估值如以下公式所示：
步骤六，测定EUT的临界干扰场强曲线，
改变测试频率，重复上述测试步骤二至步骤五，测试确定不同频率对应的EUT临界干扰场强，给出不同效应对应的EUT临界干扰场强曲线。
试验验证，
为了验证电磁辐射敏感度混响室测试方法的准确性，选择一种军用超短波通信电台为EUT，去掉天线，并将天线接口屏蔽，进行了带外干扰试验；
如图7所示，受试电台由电源模块1、功放模块2、收发模块3构成，受试电台主机外形尺寸为24.6cm×26.0cm×22.5cm，受试电台共设3根外露线缆：线1A为电源模块与市电电源间的连接线，其采用三芯非屏蔽电源线，长1.37m；、线2B为电源模块与功放模块间的连接线，其采用非屏蔽平行双线直流输出线，长0.83m，线3C为功放模块与收发模块间的连接线，其采用多芯屏蔽信号线，长0.053m；电源模块尺寸24.6cm×20.5cm×6.5cm，功放模块尺寸24.6cm×26.0cm×8.0cm，收发模块尺寸24.6cm×13.5cm×8.0cm；受试电台带有多种孔缝，包括通电源模块、收发模块的通风孔以及显示窗口；最大线度的孔缝为电源模块的电压显示窗口4，几何尺寸为5.7cm×2.8cm；如图8和9所示，分别在均匀场和混响室电磁环境进行带外临界干扰场强测试，试验流程如图10所示；受试电台正常工作电压为20V，在电磁辐射作用下，电源模块输出电压随辐射强度增加出现小幅单调下降或升高(与辐射频率有关)，继续增加辐射强度，电压显示发生较大跳变甚至出现显示错误，电源模块掉电，受试电台关机而无法工作；若停止电磁辐射，电压显示恢复正常，受试电台自动重启，恢复正常工作状态；
根据电台的干扰效应，测试时有两种可供选择的指标作为干扰判据：第一种是电台重启，第二种是电源模块输出电压发生一定的幅度的跳变；由于在较大场强下电台才会发生重启，测试时在较为敏感的频点以电台重启作为干扰判据，在相对钝感的频点以电压跳变2V作为干扰判据；采用第二种判据时需满足电压跳变的幅度与干扰场强呈正比，即电场强度越大，电压变化的幅度也越大；为此，在屏蔽室中选取部分频点测试了电压变化幅度与电场强度之间的关系，如图11所示，给出了0.55GHz和3.2GHz对电台的某一方向进行辐照时的测试结果；可以看到，虽然在不同频率的电磁波辐照下电源模块的显示电压出现升高或降低，但跳变的幅度与辐照场强均满足正比关系；因此，选择电压跳变作为干扰判据是可行的；
试验时，控制试验频率范围0.17GHz～3.3GHz；根据不同测试频段电台的敏感程度，选择了两种效应作为电台的干扰判据：在0.17GHz～2.3GHz，电台对于外界电磁辐射较为敏感，则以电台重启为干扰判据；在2.4GHz～3.3GHz，电台很难出现重启效应，但电压跳变依然存在，则以电源模块输出电压跳变2V为干扰判据；
在混响室环境中按上述测试流程进行测试，部分频点干扰概率P＜10％，用于检验干扰概率对测试误差的影响；均匀场环境测试时，通过受试电台底部转台的转动来改变被照射方向，在EUT转动一周内均匀选取36个方向进行临界干扰场强测试，取其最小值作为测试频点的EU；试验结果如表1所列，ESE、EU分别是混响室和均匀场环境受试电台临界干扰场强测试结果，
表1受试电台临界干扰场强测试结果
由于试验过程中受试电台主要出现关机-重启，属于强场电磁辐射效应，其本质原因是电磁辐射导致的电源模块输出电压急剧下降，由地电位波动引起，电磁辐射耦合通道是共地线缆的共模耦合；分析受试电台的3根外露线缆可以看出：线1-电源线需通过变压、整流滤波等电路与电源输出电路相连，二者没有共地连接；线3-信号线虽然与电源输出电路共地，但距离较远、影响较弱；线2-电源模块输出线，电磁感应信号直接影响电源模块的正常工作；因此，将线2的长度L＝0.83m代入式(4)计算受试电台电磁辐射耦合等效天线系数最大值Dmax的估值
见表1第2列，借助试验过程中测试确定的环境电场有效均值E0、干扰概率
按式(3)计算受试电台的临界干扰场强有效值ESE，与均匀场环境测试确定的受试电台的临界干扰场强有效值EU比较，计算二者之间的测试误差，如表1所列，均小于5dB；根据表1所列数据，绘制受试电台的临界干扰场强测试结果对比曲线，如图12所示；在2.4GHz～3.3GHz频段，受试电台难以出现关机-重启效应，以电源模块输出电压跳变2V为干扰判据，因此，该效应属于测控设备的扰动效应，应按式(5)计算受试电台的等效天线系数Dmax的估值；在此频段，线2的长度L＝0.83m＞3.5λ，将包围受试电台主机的最小球体半径a＝(24.62+26.02+22.52)1/2/2＝21.1cm代入式(5)计算
利用表1所列数据按式(3)计算受试电台的临界干扰场强有效值ESE，将其随频率的变化曲线用不带数据标记的实线绘于图12右侧，其与均匀场测试结果的误差普遍小于按强场效应计算的结果，但二者差距较小，说明了按效应类型选择EUT等效天线系数
的可行性和准确性；
从图12测试结果对比可以看出：除了频率最低的2个测试频点外，混响室环境受试电台的临界干扰场强测试值普遍低于均匀场环境中的测试值，且二者之间的误差均小于5dB；混响室与均匀场环境测试确定的受试电台临界干扰场强随频率的变化规律基本一致，但0.80GHz频点均匀场测试值明显偏高；说明混响室环境能够得到与均匀场环境一致的电磁辐射敏感度测试结果，且更容易找到最敏感测试状态，较好地反映了EUT的电磁辐射敏感特性；0.17GHz、0.20GHz两个频点混响室中EUT临界干扰场强测试值较高的原因可能是混响室电磁场模式数较少，EUT的敏感状态与场强最大值相遇的概率较低，可通过增加搅拌位置数加以改善。
综上可知，
混响室电磁环境各项同性、统计均匀，与电波暗室、开阔场等传统均匀场电磁环境具有本质的区别，在这两类不同的电磁环境中对EUT开展电磁辐射敏感度测试，一般难以得到等价的测试结果；但是，从混响室环境中EUT电磁辐射耦合的统计特性出发，通过建立EUT电磁辐射临界干扰场强与混响室测试环境电场强度有效均值、EUT受扰概率和等效天线方向性系数最大值之间的关系，能够解决这一技术难题，主要研究结论如下：
(1)混响室电磁环境的分布特点导致EUT位置的电场强度、极化方向难以度量和控制，采用临界干扰场强与EUT效应一一对应的方式，难以表征和测试EUT的电磁辐射敏感度；但是，混响室测试区域内各处电磁环境统计均匀、各向同性，在搅拌器步进一周的过程中，可以认为测试区域任何一点的电场强度有效均值与EUT处相同；从统计规律出发，基于混响室与均匀场中EUT敏感元件的临界干扰功率/电压相等，建立了混响室环境测试EUT临界干扰场强的计算公式，解决了混响室与均匀场环境EUT辐射敏感度测试结果等价的技术难题；
(2)确定EUT的等效天线方向性系数最大值是制约混响室环境电磁辐射敏感度测试的关键因素；从EUT的电磁辐射效应类型和等效天线的接收效率出发，给出了EUT等效天线方向性系数最大值的选取原则和量化计算公式，克服了基于干扰概率的电磁辐射敏感度混响室测试方法的工程运用障碍；
(3)测试过程中EUT出现干扰的概率、混响室的步进搅拌数是制约混响室环境EUT临界干扰场强测试准确度的主要影响因素；为保证足够的测试准确度并提高电磁场能量的利用率，干扰概率应控制在10％～60％之间，步进位置数一般应控制在30左右，在5倍的混响室最低可用频率以下频段，步进位置数宜选取50左右；
(4)一种军用超短波通信电台带外强场辐射关机-重启效应、电源模块测控扰动效应临界干扰场强测试结果表明：混响室与均匀场环境EUT的电磁辐射临界干扰场强测试值随辐射频率的变化规律一致，二者之间的最大测试误差为5dB，能够满足工程测试需求。
上述实施例，仅是本发明的较佳实施方式，故凡依本发明专利申请范围所述的构造、特征及原理所做的等效变化或修饰，均包括于本发明专利申请范围内。
Claims (5)
1.一种基于干扰概率的混响室环境电磁辐射敏感度测试方法，其特征在于，包括以下步骤：
步骤一，规划搅拌器步进规则，
根据混响室的搅拌器数量，规划搅拌器步进规则，使搅拌器步进一周的混响室搅拌位置数在30～50之间且使各个搅拌器等角度步进的位置数大致相等；
步骤二，混响室馈入功率预调节，
根据EUT测试要求，选定测试频率，设定混响室激励天线的初始馈入功率，对EUT进行电磁辐射效应快速测试，监视或检测确定EUT的效应类型；参照变步长升降法调节混响室激励天线的馈入功率，在搅拌器步进一周的过程中，使EUT出现某一效应或多种效应的概率分别达到20％～60％；
步骤三，EUT电磁辐射效应测试，
固定混响室激励天线的馈入功率，按搅拌器步进规则步进搅拌器一周，记录EUT出现每种效应的概率估计值
测试过程中，在测试区域远离EUT的某一固定位置测试搅拌器处于每一搅拌位置时的电场强度有效值E0n，按以下公式计算混响室测试环境电场强度有效均值E0：
步骤四，确定EUT等效天线方向性系数最大值，
根据EUT出现的电磁辐射效应及EUT、线缆电尺寸，按EUT电磁辐射耦合等效天线方向性系数的选取原则，分别采用用频设备天线系数或按数学拟合估值公式计算，确定EUT等效天线方向性系数最大值的估计值；具体地，
不同电磁辐射效应下等效天线方向性系数的选取和计算，
(1)当EUT的电磁辐射效应为用频设备阻塞、虚警与损伤效应，则采用混响室测试用频设备的电磁辐射敏感度时，直接用天线系数代替Dmax即可；
(2)当EUT的电磁辐射效应为强场电磁辐射效应，则采用混响室测试该类设备的电磁辐射敏感度时，按EUT接地线或共地线缆的电长度，采用线缆共模耦合等效天线方向性系数代替Dmax，计算受试系统的临界辐射干扰场强Es；当线缆长度L＜3.5λ时，其中，λ为波长，采用线缆共模耦合等效天线方向性系数对数中值的平滑线代替Dmax，当L≥3.5λ时Dmax取整不再变化；为方便计算，给出Dmax估值的数学拟合估值公式：
(3)当EUT的电磁辐射效应为测控设备的扰动效应，采用混响室进行效应试验时，若线缆差模耦合和孔缝耦合单独作用，则低频段以线缆耦合为主，高频段以孔缝耦合为主；当受试测控设备在线缆差模耦合和孔缝耦合联合作用下出现扰动效应时，则等效天线方向性系数最大值Dmax用统一的拟合公式来估值：
当线缆长度L＜3.5λ时，a＝L；否则a为包围EUT的最小球体半径；
步骤五，EUT临界干扰场强的定量计算，
根据EUT出现每种效应的概率估计值
混响室测试环境电场强度有效均值E0和每种效应对应的EUT等效天线方向性系数最大值的估计值
利用以下公式：
计算EUT在测试频点的临界干扰场强估值
则临界干扰场强估值如以下公式所示：
步骤六，测定EUT的临界干扰场强曲线，
改变测试频率，重复上述测试步骤二至步骤五，测试确定不同频率对应的EUT临界干扰场强，给出不同效应对应的EUT临界干扰场强曲线。
2.根据权利要求1所述的基于干扰概率的混响室环境电磁辐射敏感度测试方法，其特征在于，所述强场电磁辐射效应包括死机、重启和功能紊乱效应；所述功能紊乱效应包括显示乱码和功能失效。
3.根据权利要求1所述的基于干扰概率的混响室环境电磁辐射敏感度测试方法，其特征在于，当EUT的电磁辐射效应为强场电磁辐射效应，采用混响室测试该类设备的电磁辐射敏感度时，若EUT具有多根不同长度的接地线或共地线缆时，则根据测试频率选择与半波长较接近的线缆长度计算Dmax的估值。
4.根据权利要求1所述的基于干扰概率的混响室环境电磁辐射敏感度测试方法，其特征在于，当EUT的电磁辐射效应为测控设备的扰动效应，若EUT采用屏蔽线缆，其差模干扰信号源于线缆屏蔽层的共模感应电流，当线缆长度L＜3.5λ时，采用公式(2)计算Dmax的估值；否则，按公式(3)进行计算。
5.根据权利要求1所述的基于干扰概率的混响室环境电磁辐射敏感度测试方法，其特征在于，所述强场电磁辐射效应还包括射频前端硬损伤效应。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110450645.5A CN113295941B (zh) | 2021-04-25 | 2021-04-25 | 一种基于干扰概率的混响室环境电磁辐射敏感度测试方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110450645.5A CN113295941B (zh) | 2021-04-25 | 2021-04-25 | 一种基于干扰概率的混响室环境电磁辐射敏感度测试方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113295941A CN113295941A (zh) | 2021-08-24 |
| CN113295941B true CN113295941B (zh) | 2022-10-25 |
Family
ID=77320326
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110450645.5A Active CN113295941B (zh) | 2021-04-25 | 2021-04-25 | 一种基于干扰概率的混响室环境电磁辐射敏感度测试方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113295941B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113687276B (zh) * | 2021-08-30 | 2022-02-22 | 北京航空航天大学 | 一种基于瞬态电磁环境的设备敏感指标设计方法 |
| TWI771232B (zh) * | 2021-11-23 | 2022-07-11 | 英業達股份有限公司 | 測試輻射敏感度之系統及方法 |
| CN114460382A (zh) * | 2022-01-25 | 2022-05-10 | 中国船舶重工集团公司第七二四研究所 | 电机控制器辐射敏感度快速评估方法 |
| CN114884586B (zh) * | 2022-04-02 | 2024-04-30 | 深圳市通用测试系统有限公司 | 用于混响室的搅拌器和具有该搅拌器的混响室 |
| CN117634187B (zh) * | 2023-11-28 | 2024-06-07 | 北京航空航天大学 | 一种电磁场耦合放电对燃油危害的安全阈值评估方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108445313A (zh) * | 2018-01-31 | 2018-08-24 | 中国人民解放军陆军工程大学 | 用频装备电磁辐射效应综合预测方法及终端设备 |
| KR20190115272A (ko) * | 2018-04-02 | 2019-10-11 | 한국전자통신연구원 | 전자파잔향실에서의 복사성 내성 평가를 위한 전기장 형성 장치 및 방법 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107037283B (zh) * | 2017-03-06 | 2019-08-20 | 中国人民解放军军械工程学院 | 混响室条件下辐射敏感度测试方法及装置 |
| CN108490280B (zh) * | 2018-01-31 | 2020-06-16 | 中国人民解放军陆军工程大学 | 用频装备带内多频电磁辐射效应预测方法及终端设备 |
| CN108647393B (zh) * | 2018-04-11 | 2022-03-04 | 北京空间飞行器总体设计部 | 一种航天器设备电磁辐射敏感度控制量级设计方法 |
| CN112485538B (zh) * | 2020-11-03 | 2022-12-20 | 飞特质科(北京)计量检测技术有限公司 | 一种提高混响室射频辐射敏感度校验精度的系统 |
-
2021
- 2021-04-25 CN CN202110450645.5A patent/CN113295941B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108445313A (zh) * | 2018-01-31 | 2018-08-24 | 中国人民解放军陆军工程大学 | 用频装备电磁辐射效应综合预测方法及终端设备 |
| KR20190115272A (ko) * | 2018-04-02 | 2019-10-11 | 한국전자통신연구원 | 전자파잔향실에서의 복사성 내성 평가를 위한 전기장 형성 장치 및 방법 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113295941A (zh) | 2021-08-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113295941B (zh) | 一种基于干扰概率的混响室环境电磁辐射敏感度测试方法 | |
| KR101360280B1 (ko) | 흡수재를 구비하지 않은 다중채널 근접장 측정 시스템 | |
| US20110193566A1 (en) | Multichannel absorberless near field measurement system | |
| CN102981086B (zh) | 一种电压驱动型辐射源电磁辐射分析测量方法 | |
| KR20190022960A (ko) | 전자파 잔향실 기반의 복사성 방출 측정 방법 및 복사성 방출 측정 시스템 | |
| CN103063939B (zh) | 一种地面级联系统外部射频电磁环境的试验方法 | |
| Hong et al. | Optimization of stirrer with various parameters in reverberation chamber | |
| Sarkar et al. | Low cost RF sensor for partial discharge detection of high voltage apparatus | |
| US11631927B2 (en) | Probe calibration system and method for electromagnetic compatibility testing | |
| Silfverskiold et al. | Microwave field-to-wire coupling measurements in anechoic and reverberation chambers | |
| García-Pérez et al. | Time-domain shielding effectiveness of enclosures against a plane wave excitation | |
| Paez et al. | A proposed method for quantifying uncertainty in RF immunity testing due to EUT presence | |
| Murata et al. | Susceptibility of notebook computers to HPM | |
| Ali et al. | FDTD analysis of rectangular waveguide in receiving mode as EMI sensors | |
| Merfeldas et al. | Characterization of radiated susceptibility directivity for microstrip and coplanar lines | |
| CN113295940B (zh) | 一种等效天线方向性系数最大值的近似估计方法 | |
| Boozari et al. | Investigation on transient response of linear dipole antennas in the presence of HEMP based on equivalent circuit | |
| Coetzer et al. | Investigating PV Module Far-Field Radiation: A Measurement and Circuit-CEM Hybrid Modelling Approach | |
| Tremola et al. | An interpolation method to calibrate electromagnetic probes in semi-anechoic chambers | |
| Poodineh et al. | Standard Deviation Characterization of a Small Size Reverberation Chamber by Using Full wave Simulation and E Field Probe | |
| US20240329108A1 (en) | Measuring system and measuring method of antenna pattern based on near field to far field transformation | |
| Hanmantrao et al. | Far field prediction of a PCB using simulation and validation | |
| Abdulhameed et al. | Using of LPDA Antenna for UFA Evaluation in Radiated Immunity Test | |
| Choeysakul et al. | Reducing the LUF of a Reverberation Chamber based on the concept of MIMO for Electromagnetic Emission Measurements for Radio Astronomy Applications | |
| MUSIAL et al. | Threats associated with electromagnetic field immunity tests of electronic devices according to selected standards |
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 |