CN117452081B - 一种电磁干扰计算方法、装置及存储介质、电子终端 - Google Patents

一种电磁干扰计算方法、装置及存储介质、电子终端 Download PDF

Info

Publication number
CN117452081B
CN117452081B CN202311798253.3A CN202311798253A CN117452081B CN 117452081 B CN117452081 B CN 117452081B CN 202311798253 A CN202311798253 A CN 202311798253A CN 117452081 B CN117452081 B CN 117452081B
Authority
CN
China
Prior art keywords
electromagnetic interference
equation set
maxwell equation
working environment
space
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
Application number
CN202311798253.3A
Other languages
English (en)
Other versions
CN117452081A (zh
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.)
State Grid Corp of China SGCC
State Grid Tianjin Electric Power Co Ltd
Marketing Service Center of State Grid Tianjin Electric Power Co Ltd
Marketing Service Center of State Grid Hebei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Tianjin Electric Power Co Ltd
Marketing Service Center of State Grid Tianjin Electric Power Co Ltd
Marketing Service Center of State Grid Hebei Electric Power Co Ltd
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 State Grid Corp of China SGCC, State Grid Tianjin Electric Power Co Ltd, Marketing Service Center of State Grid Tianjin Electric Power Co Ltd, Marketing Service Center of State Grid Hebei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202311798253.3A priority Critical patent/CN117452081B/zh
Publication of CN117452081A publication Critical patent/CN117452081A/zh
Application granted granted Critical
Publication of CN117452081B publication Critical patent/CN117452081B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08Measuring electromagnetic field characteristics
    • G01R29/0807Measuring electromagnetic field characteristics characterised by the application
    • G01R29/0814Field measurements related to measuring influence on or from apparatus, components or humans, e.g. in ESD, EMI, EMC, EMP testing, measuring radiation leakage; detecting presence of micro- or radiowave emitters; dosimetry; testing shielding; measurements related to lightning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141Discrete Fourier transforms
    • G06F17/142Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Electromagnetism (AREA)
  • Discrete Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Testing Electric Properties And Detecting Electric Faults (AREA)

Abstract

本发明公开了一种电磁干扰计算方法、装置及存储介质、电子终端。其中,电磁干扰计算方法,包括如下步骤:对智能电表所处工作环境空间采用六面体进行空间网格剖分,得到离散后的空间网格;基于所述空间网格对麦克斯韦方程组中的电场方程和磁场方程进行差分处理,得到频域下麦克斯韦方程组的离散形式;将加载的智能电表所处工作环境空间中的噪音设备导线中的电流密度作为激励源加入到对应空间网格处的麦克斯韦方程组中的电场方程中;计算麦克斯韦方程组中的可行解,通过频域到时域的傅里叶快速变换得到时域下电磁干扰的电磁场强度信息。利用本发明能够精确计算智能电表周围的电磁干扰强度。

Description

一种电磁干扰计算方法、装置及存储介质、电子终端
技术领域
本发明涉及智能电表应用技术领域,具体涉及了一种电磁干扰计算方法、装置及存储介质、电子终端。
背景技术
智能电表是一种集智能化、数字化、通信化于一体的电表,具有远程抄表、实时监测等功能,是电力行业数字化升级的重要组成部分。目前,由于智能电表中包含许多电子元器件,在存在复杂的电磁环境中,其计量性能很容易受到电磁干扰的影响。然而,在复杂的电磁环境中,麦克斯韦方程连续求解计算量过大,难以精确计算智能电表周围的电磁干扰强度。
发明内容
本发明的目的在于提供一种电磁干扰计算方法、装置及存储介质、电子终端,以解决复杂的电磁环境中,由于麦克斯韦方程连续求解计算量过大,难以精确计算智能电表周围的电磁干扰强度的问题。
为实现本发明的目的,本发明提供的技术方案具体如下:
第一方面
本发明提供了一种电磁干扰计算方法,包括如下步骤:
步骤1:对智能电表所处工作环境空间采用六面体进行空间网格剖分,得到离散后的空间网格;
步骤2:基于所述空间网格对麦克斯韦方程组中的电场方程和磁场方程进行差分处理,得到频域下麦克斯韦方程组的离散形式;
步骤3:将加载的智能电表所处工作环境空间中的噪音设备导线中的电流密度作为激励源加入到对应空间网格处的麦克斯韦方程组中的电场方程中;
步骤4:计算麦克斯韦方程组中的可行解,通过频域到时域的傅里叶快速变换得到时域下电磁干扰的电磁场强度信息。
第二方面
本发明提供了一种电磁干扰计算装置,包括网格剖分单元、差分离散单元、激励源加入单元以及求解单元:
所述网格剖分单元用于对智能电表所处工作环境空间采用六面体进行空间网格剖分,得到离散后的空间网格;
所述差分离散单元用于基于所述空间网格对麦克斯韦方程组中的电场方程和磁场方程进行差分处理,得到频域下麦克斯韦方程组的离散形式;
所述激励源加入单元用于将加载的智能电表所处工作环境空间中的噪音设备导线中的电流密度作为激励源加入到对应空间网格处的麦克斯韦方程组中的电场方程中;
所述求解单元用于计算麦克斯韦方程组中的可行解,通过频域到时域的傅里叶快速变换得到时域下电磁干扰的电磁场强度信息。
第三方面
本发明提供了一种存储介质,所述存储介质中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、所述至少一段程序、所述代码集或指令集由处理器加载并执行以实现如上述的电磁干扰计算方法。
第四方面
本发明提供了一种电子终端,所述电子终端包括处理器和存储器,所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、所述至少一段程序、所述代码集或指令集由所述处理器加载并执行以实现如上述的电磁干扰计算方法。
与现有技术相比,本发明的有益效果为:
本发明提供的方案通过调研智能电能表工作环境,通过多重网格剖分提升电表周围电磁干扰的计算精度,进一步对麦克斯韦方程组中的电场方程和磁场方程进行差分,得到频域下麦克斯韦方程的离散形式,之后将智能电表所处工作环境空间中的噪音设备作为激励源加入到对应网格的电场中,最后通过频域到时域的傅里叶快速变换得到时域下电磁干扰的电磁场强度信息,能够精确计算智能电表周围的电磁干扰强度。
附图说明
图1为本发明实施例中网格剖分示意图;
图2为本发明实施例中电磁场分布示意图;
图3为本发明实施例中智能电表周围电场强度示意图。
具体实施方式
以下结合附图和具体实施例对本发明作进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明的目的是针对复杂的电磁环境中智能电表周围的电磁干扰强度难以精确计算的问题,提出了一种电磁干扰计算方法、装置及存储介质、电子终端。通过调研智能电能表工作环境,通过多重网格剖分提升智能电表周围电磁干扰的计算精度,进一步对麦克斯韦方程组中的电场方程和磁场方程进行差分,得到频域下麦克斯韦方程组的离散形式,之后将智能电表所处工作环境空间中的噪音设备作为激励源加入到对应网格的电场中,最后通过频域到时域的傅里叶快速变换得到时域下电磁干扰的电磁场强度信息。
如图1-图3所示,本发明提供了一种电磁干扰计算方法,包括如下步骤:
步骤1:对智能电表所处工作环境空间采用六面体进行空间网格剖分,得到离散后的空间网格;
步骤2:基于所述空间网格对麦克斯韦方程组中的电场方程和磁场方程进行差分处理,得到频域下麦克斯韦方程组的离散形式;
步骤3:将加载的智能电表所处工作环境空间中的噪音设备导线中的电流密度作为激励源加入到对应空间网格处的麦克斯韦方程组中的电场方程中;
步骤4:计算麦克斯韦方程组中的可行解,通过频域到时域的傅里叶快速变换得到时域下电磁干扰的电磁场强度信息。
优选地,所述步骤1中,包括一次网格剖分和二次网格剖分,其中,一次网格剖分包括在所述智能电表所处工作环境空间的第一范围内采用第一预设步长进行网格剖分,在所述智能电表所处工作环境空间的第二范围内采用第二预设步长进行网格剖分,二次网格剖分为在一次网格剖分结果基础上,对所述智能电表所处工作环境空间采用第三预设步长进行网格剖分;其中,所述第二范围大于所述第一范围,所述第二预设步长大于所述第一预设步长,所述第一预设步长大于所述第三预设步长。
需要说明的是采用40*40*40个网格对智能电表工作环境空间进行剖分,所述第二预设步长为12.8m,所述第一预设步长为1mm,所述第三预设步长为0.5mm。所述第二范围为智能电表工作环境空间外围边界处,所述第一范围为智能电表所处的空间中心处。
优选地,所述步骤2,具体包括如下:
所述麦克斯韦方程组按照所述空间网格在直角坐标系下沿x、y和z三个方向展开后得到:
(1)
(2)
(3)
(4)
(5)
(6)
其中,为电导率,/>为偏微分,/>为真空磁导率,ω为角频率,/>、/>和/>为x、y和z三个方向上的电场分量,/>、/>和/>为x、y和z三个方向上的磁场分量,/>,/>为x、y和z三个方向上噪音设备导线中的电流密度,i为虚数;
对展开后的麦克斯韦方程组进行差分离散,利用差分项代替微分项整理得到频域下麦克斯韦方程组的离散形式:
(7)
(8)
(9)
(10)
(11)
(12)
其中,、/>、/>为直角坐标系下三个方向上的网格步长,i为虚数。
优选地,所述步骤3,具体包括如下:
将加载的智能电表所处工作环境空间中的噪音设备导线通过的电流除以导线的横截面积,得到噪音设备导线中的电流密度,将所述电流密度作为激励源加入到对应空间网格方向上的麦克斯韦方程组中的电场方程中,对应电流密度计算如下:
(13)
(14)
(15)。
优选地,所述步骤4,具体包括如下:将麦克斯韦方程组中的电场方程替换到麦克斯韦方程组中的磁场方程中,并将智能电表所处工作环境空间的第二范围内电场强度和磁场强度全部设置为0向量,通过feasp求解器计算求解,得到直角坐标系下三个方向上的电场分量,通过傅里叶快速变换来获得时域下的电磁干扰的电磁场强度信息。
与上述方法相对应地,本发明还提供了一种电磁干扰计算装置,包括网格剖分单元、差分离散单元、激励源加入单元以及求解单元:
所述网格剖分单元用于对智能电表所处工作环境空间采用六面体进行空间网格剖分,得到离散后的空间网格;
所述差分离散单元用于基于所述空间网格对麦克斯韦方程组中的电场方程和磁场方程进行差分处理,得到频域下麦克斯韦方程组的离散形式;
所述激励源加入单元用于将加载的智能电表所处工作环境空间中的噪音设备导线中的电流密度作为激励源加入到对应空间网格处的麦克斯韦方程组中的电场方程中;
所述求解单元用于计算麦克斯韦方程组中的可行解,通过频域到时域的傅里叶快速变换得到时域下电磁干扰的电磁场强度信息。
优选地,所述差分离散单元,具体用于:
所述麦克斯韦方程组按照所述空间网格在直角坐标系下沿x、y和z三个方向展开后得到:
(16)
(17)
(18)
(19)
(20)
(21)
其中,为电导率,/>为偏微分/>为真空磁导率,ω为角频率,/>、/>和/>为x、y和z三个方向上的电场分量,/>、/>和/>为x、y和z三个方向上的磁场分量,/>,/>,/>为x、y和z三个方向上噪音设备导线中的电流密度,i为虚数;
对展开后的麦克斯韦方程组进行差分离散,利用差分项代替微分项整理得到频域下麦克斯韦方程组的离散形式:
(22)
(23)
(24)
(25)
(26)
(27)
其中,、/>、/>为直角坐标系下三个方向上的网格步长。
优选地,所述激励源加入单元,具体用于:
将加载的智能电表所处工作环境空间中的噪音设备导线通过的电流除以导线的横截面积,得到噪音设备导线中的电流密度,将所述电流密度作为激励源加入到对应空间网格方向上的麦克斯韦方程组中的电场方程中,对应电流密度计算如下:
(28)
(29)
(30)。
另外,本发明还提供了一种存储介质,所述存储介质中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、所述至少一段程序、所述代码集或指令集由处理器加载并执行以实现如上述的电磁干扰计算方法。
另外,本发明还提供了一种电子终端,所述电子终端包括处理器和存储器,所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、所述至少一段程序、所述代码集或指令集由所述处理器加载并执行以实现如上述的电磁干扰计算方法。
以上结合附图详细描述了本发明实施例的可选实施方式,但是,本发明实施例并不限于上述实施方式中的具体细节,在本发明实施例的技术构思范围内,可以对本发明实施例的技术方案进行多种简单变型,这些简单变型均属于本发明实施例的保护范围。

Claims (7)

1.一种电磁干扰计算方法,其特征在于,包括如下步骤:
步骤1:对智能电表所处工作环境空间采用六面体进行空间网格剖分,得到离散后的空间网格;
步骤2:基于所述空间网格对麦克斯韦方程组中的电场方程和磁场方程进行差分处理,得到频域下麦克斯韦方程组的离散形式;
步骤3:将加载的智能电表所处工作环境空间中的噪音设备导线中的电流密度作为激励源加入到对应空间网格处的麦克斯韦方程组中的磁场方程中;
步骤4:计算麦克斯韦方程组中的可行解,通过频域到时域的傅里叶快速变换得到时域下电磁干扰的电磁场强度信息;
所述步骤3,具体包括如下:
将加载的智能电表所处工作环境空间中的噪音设备导线通过的电流除以导线的横截面积,得到噪音设备导线中的电流密度,将所述电流密度作为激励源加入到对应空间网格方向上的麦克斯韦方程组中的磁场方程中,对应电流密度计算如下:
其中,其中,Δx、Δy、Δz为直角坐标系下三个方向上的网格步长;Jx,Jy,Jz为x、y和z三个方向上噪音设备导线中的电流密度;
所述步骤4,具体包括如下:将麦克斯韦方程组中的电场方程替换到麦克斯韦方程组中的磁场方程中,并将智能电表所处工作环境空间的第二范围内电场强度和磁场强度全部设置为0向量,通过feasp求解器计算求解,得到直角坐标系下三个方向上的电场分量,通过傅里叶快速变换来获得时域下的电磁干扰的电磁场强度信息。
2.根据权利要求1所述的一种电磁干扰计算方法,其特征在于,所述步骤1中,包括一次网格剖分和二次网格剖分,其中,一次网格剖分包括在所述智能电表所处工作环境空间的第一范围内采用第一预设步长进行网格剖分,在所述智能电表所处工作环境空间的第二范围内采用第二预设步长进行网格剖分,二次网格剖分为在一次网格剖分结果基础上,对所述智能电表所处工作环境空间采用第三预设步长进行网格剖分;其中,所述第二范围大于所述第一范围,所述第二预设步长大于所述第一预设步长,所述第一预设步长大于所述第三预设步长。
3.根据权利要求2所述的一种电磁干扰计算方法,其特征在于,所述步骤2,具体包括如下:
所述麦克斯韦方程组按照所述空间网格在直角坐标系下沿x、y和z三个方向展开后得到:
其中,σ为电导率,为偏微分,μ0为真空磁导率,ω为角频率,Ex、Ey和Ez为x、y和z三个方向上的电场分量,Hx、Hy和Hz为x、y和z三个方向上的磁场分量,i为虚数;
对展开后的麦克斯韦方程组进行差分离散,利用差分项代替微分项整理得到频域下麦克斯韦方程组的离散形式:
4.一种电磁干扰计算装置,其特征在于,包括网格剖分单元、差分离散单元、激励源加入单元以及求解单元:
所述网格剖分单元用于对智能电表所处工作环境空间采用六面体进行空间网格剖分,得到离散后的空间网格;
所述差分离散单元用于基于所述空间网格对麦克斯韦方程组中的电场方程和磁场方程进行差分处理,得到频域下麦克斯韦方程组的离散形式;
所述激励源加入单元用于将加载的智能电表所处工作环境空间中的噪音设备导线中的电流密度作为激励源加入到对应空间网格处的麦克斯韦方程组中的磁场方程中;
所述求解单元用于计算麦克斯韦方程组中的可行解,通过频域到时域的傅里叶快速变换得到时域下电磁干扰的电磁场强度信息;
步骤3,具体包括如下:
将加载的智能电表所处工作环境空间中的噪音设备导线通过的电流除以导线的横截面积,得到噪音设备导线中的电流密度,将所述电流密度作为激励源加入到对应空间网格方向上的麦克斯韦方程组中的磁场方程中,对应电流密度计算如下:
其中,其中,Δx、Δy、Δz为直角坐标系下三个方向上的网格步长;Jx,Jy,Jz为x、y和z三个方向上噪音设备导线中的电流密度;
步骤4,具体包括如下:将麦克斯韦方程组中的电场方程替换到麦克斯韦方程组中的磁场方程中,并将智能电表所处工作环境空间的第二范围内电场强度和磁场强度全部设置为0向量,通过feasp求解器计算求解,得到直角坐标系下三个方向上的电场分量,通过傅里叶快速变换来获得时域下的电磁干扰的电磁场强度信息。
5.根据权利要求4所述的一种电磁干扰计算装置,其特征在于,所述差分离散单元,具体用于:
所述麦克斯韦方程组按照所述空间网格在直角坐标系下沿x、y和z三个方向展开后得到:
其中,σ为电导率,为偏微分,μ0为真空磁导率,ω为角频率,Ex、Ey和Ez为x、y和z三个方向上的电场分量,Hx、Hy和Hz为x、y和z三个方向上的磁场分量,i为虚数;
对展开后的麦克斯韦方程组进行差分离散,利用差分项代替微分项整理得到频域下麦克斯韦方程组的离散形式:
6.一种存储介质,其特征在于,所述存储介质中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、所述至少一段程序、所述代码集或指令集由处理器加载并执行以实现如权利要求1-3中任一项所述的电磁干扰计算方法。
7.一种电子终端,其特征在于,所述电子终端包括处理器和存储器,所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、所述至少一段程序、所述代码集或指令集由所述处理器加载并执行以实现如权利要求1-3中任一项所述的电磁干扰计算方法。
CN202311798253.3A 2023-12-26 2023-12-26 一种电磁干扰计算方法、装置及存储介质、电子终端 Active CN117452081B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311798253.3A CN117452081B (zh) 2023-12-26 2023-12-26 一种电磁干扰计算方法、装置及存储介质、电子终端

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311798253.3A CN117452081B (zh) 2023-12-26 2023-12-26 一种电磁干扰计算方法、装置及存储介质、电子终端

Publications (2)

Publication Number Publication Date
CN117452081A CN117452081A (zh) 2024-01-26
CN117452081B true CN117452081B (zh) 2024-04-30

Family

ID=89595223

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311798253.3A Active CN117452081B (zh) 2023-12-26 2023-12-26 一种电磁干扰计算方法、装置及存储介质、电子终端

Country Status (1)

Country Link
CN (1) CN117452081B (zh)

Citations (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000003354A (ja) * 1998-06-12 2000-01-07 Matsushita Electric Ind Co Ltd 電磁界の解析方法およびその装置
EP2068262A1 (en) * 2007-11-27 2009-06-10 Fujitsu Limited Numerical simulation with anisotropic local grid refinement
CN102930071A (zh) * 2012-08-29 2013-02-13 电子科技大学 基于非匹配网格的周期结构的三维电磁场仿真模拟方法
CN104833860A (zh) * 2015-03-04 2015-08-12 西北工业大学 一种电磁波跨越海与空界面电磁场的获得方法
CN106055812A (zh) * 2016-06-07 2016-10-26 广东电网有限责任公司电力科学研究院 一种三维瞬态电磁场数值模拟分析方法及装置
CN106443206A (zh) * 2016-10-24 2017-02-22 广州供电局有限公司 高压导线表面电场强度的测量方法及测量装置
CN107422387A (zh) * 2017-06-27 2017-12-01 吉林大学 一种虚拟时域有限差分的瞬变电磁发射源加载方法
CN108287268A (zh) * 2017-12-28 2018-07-17 华中科技大学 一种分布式输电网络工频电磁场分析方法
CN110263463A (zh) * 2019-06-26 2019-09-20 中国人民解放军陆军装甲兵学院 基于Ansys Maxwell软件电磁感应特性仿真分析方法
CN111259533A (zh) * 2020-01-13 2020-06-09 武汉理工大学 Gis隔离开关操作下套管处特快速暂态电磁辐射的计算方法
CN111639447A (zh) * 2020-04-30 2020-09-08 南京理工大学 多级局部时间步进技术的任意高阶混合网格时域不连续伽辽金方法
CN111697415A (zh) * 2020-06-04 2020-09-22 上海理工大学 一种基于Weyl半金属-纳米介孔复合结构的太赫兹增强方法
CN113158527A (zh) * 2021-05-14 2021-07-23 中国空气动力研究与发展中心计算空气动力研究所 一种基于隐式fvfd计算频域电磁场的方法
CN113655733A (zh) * 2021-09-09 2021-11-16 中车长春轨道客车股份有限公司 一种轨道交通车辆计轴器磁场emc半实物仿真方法
CN114137318A (zh) * 2021-11-16 2022-03-04 中国石油大学(华东) 海洋感应电磁场算法、系统、计算机设备和存储介质
CN114528716A (zh) * 2022-03-01 2022-05-24 西安热工研究院有限公司 一种应用于多尺度电磁波问题分析的高效数值仿真方法
CN115017782A (zh) * 2022-08-08 2022-09-06 中南大学 考虑介质各向异性的三维天然源电磁场计算方法
CN115754527A (zh) * 2022-11-13 2023-03-07 重庆理工大学 一种快速计算电动汽车磁场辐射发射量的简易方法
CN115879276A (zh) * 2022-11-03 2023-03-31 中国空气动力研究与发展中心计算空气动力研究所 一种目标对象的电磁特性分析方法、装置及设备和介质
CN116822434A (zh) * 2023-07-19 2023-09-29 经纬恒润(天津)研究开发有限公司 功率开关模块电磁干扰仿真方法、装置、设备及介质
CN117236110A (zh) * 2023-04-13 2023-12-15 安徽大学 基于电热多物理场随机fdtd的电磁辐射强度评估方法

Patent Citations (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000003354A (ja) * 1998-06-12 2000-01-07 Matsushita Electric Ind Co Ltd 電磁界の解析方法およびその装置
EP2068262A1 (en) * 2007-11-27 2009-06-10 Fujitsu Limited Numerical simulation with anisotropic local grid refinement
CN102930071A (zh) * 2012-08-29 2013-02-13 电子科技大学 基于非匹配网格的周期结构的三维电磁场仿真模拟方法
CN104833860A (zh) * 2015-03-04 2015-08-12 西北工业大学 一种电磁波跨越海与空界面电磁场的获得方法
CN106055812A (zh) * 2016-06-07 2016-10-26 广东电网有限责任公司电力科学研究院 一种三维瞬态电磁场数值模拟分析方法及装置
CN106443206A (zh) * 2016-10-24 2017-02-22 广州供电局有限公司 高压导线表面电场强度的测量方法及测量装置
CN107422387A (zh) * 2017-06-27 2017-12-01 吉林大学 一种虚拟时域有限差分的瞬变电磁发射源加载方法
CN108287268A (zh) * 2017-12-28 2018-07-17 华中科技大学 一种分布式输电网络工频电磁场分析方法
CN110263463A (zh) * 2019-06-26 2019-09-20 中国人民解放军陆军装甲兵学院 基于Ansys Maxwell软件电磁感应特性仿真分析方法
CN111259533A (zh) * 2020-01-13 2020-06-09 武汉理工大学 Gis隔离开关操作下套管处特快速暂态电磁辐射的计算方法
CN111639447A (zh) * 2020-04-30 2020-09-08 南京理工大学 多级局部时间步进技术的任意高阶混合网格时域不连续伽辽金方法
CN111697415A (zh) * 2020-06-04 2020-09-22 上海理工大学 一种基于Weyl半金属-纳米介孔复合结构的太赫兹增强方法
CN113158527A (zh) * 2021-05-14 2021-07-23 中国空气动力研究与发展中心计算空气动力研究所 一种基于隐式fvfd计算频域电磁场的方法
CN113655733A (zh) * 2021-09-09 2021-11-16 中车长春轨道客车股份有限公司 一种轨道交通车辆计轴器磁场emc半实物仿真方法
CN114137318A (zh) * 2021-11-16 2022-03-04 中国石油大学(华东) 海洋感应电磁场算法、系统、计算机设备和存储介质
CN114528716A (zh) * 2022-03-01 2022-05-24 西安热工研究院有限公司 一种应用于多尺度电磁波问题分析的高效数值仿真方法
CN115017782A (zh) * 2022-08-08 2022-09-06 中南大学 考虑介质各向异性的三维天然源电磁场计算方法
CN115879276A (zh) * 2022-11-03 2023-03-31 中国空气动力研究与发展中心计算空气动力研究所 一种目标对象的电磁特性分析方法、装置及设备和介质
CN115754527A (zh) * 2022-11-13 2023-03-07 重庆理工大学 一种快速计算电动汽车磁场辐射发射量的简易方法
CN117236110A (zh) * 2023-04-13 2023-12-15 安徽大学 基于电热多物理场随机fdtd的电磁辐射强度评估方法
CN116822434A (zh) * 2023-07-19 2023-09-29 经纬恒润(天津)研究开发有限公司 功率开关模块电磁干扰仿真方法、装置、设备及介质

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
FDTD法建模中激励源的选择与设置;李蓉, 张林昌;铁道学报;20010810(04);全文 *
一种适用于配电开关辐射场仿真分析方法;范闻博;吴燕;赵轶民;关石磊;余勇祥;耿英三;;高压电器;20170916(09);全文 *
三相多功能谐波电能表的EMI分析与EMC设计;滕召胜;罗志坤;唐求;周良璋;李小青;;湖南大学学报(自然科学版);20101125(11);全文 *
条纹变像管时间畸变的分析;惠丹丹;田进寿;卢裕;王俊锋;温文龙;梁玲亮;陈琳;;物理学报;20160629(15);全文 *

Also Published As

Publication number Publication date
CN117452081A (zh) 2024-01-26

Similar Documents

Publication Publication Date Title
CN106295120B (zh) 一种目标电磁散射特性快速计算方法
US6772076B2 (en) Electromagnetic field analysis method based on FDTD method, medium representation method in electromagnetic field analysis, simulation device, and storage medium
CN104123408B (zh) 基于有限元仿真分析的智能电表防外界电磁场干扰方法
CN111062153A (zh) 一种干式空心电抗器磁场分析方法、装置及存储介质
CN104899375A (zh) 一种用于高压电机的主绝缘防晕结构优化计算方法
CN105678002A (zh) 等离子体粒子-场自洽系统长期大规模高保真模拟方法
CN112165179A (zh) 一种地铁接触网与电网的电磁耦合的分析方法及装置
CN117452081B (zh) 一种电磁干扰计算方法、装置及存储介质、电子终端
JP3980811B2 (ja) 電磁界強度算出装置、電磁界強度算出方法およびその方法をコンピュータに実行させるプログラムを記録したコンピュータ読み取り可能な記録媒体
CN110110361B (zh) 柔性直流换流阀阀塔的电场仿真方法、装置及存储介质
CN114676519A (zh) 变压器铁心励磁特性的有限元计算方法、装置及终端
JPH11295365A (ja) 電磁界解析方法
CN115688516A (zh) 一种无线通信设备电磁屏蔽罩性能优化方法及装置
CN112347679B (zh) 一种带有电磁色散的超材料仿真方法及装置
Kong et al. Study on the stability and numerical error of the four-stages split-step FDTD method including lumped inductors
CN106294898B (zh) 一种加速分析介质目标电磁散射特性的复点源求解方法
Golias et al. Adaptive refinement in 2‐D finite element applications
Celuch et al. Properties of the FDTD method relevant to the analysis of microwave power problems
JP2004004054A (ja) Fdtd法を用いた電磁界解析方法、電磁界解析における媒質表現方法、シミュレーション装置、及びプログラム
Pfeiffer et al. Locally Corrected Nyström to Moment Method Conversion for Volume Integral Equations
Zhang et al. Multiquadrics collocation method for transient eddy current problems
Huang et al. Study on Improved CG Algorithm for Reradiation Interference Calculation of Transmission Steel Tower
JP2023156195A (ja) シミュレーションプログラム、シミュレーション方法、およびシミュレーション装置
CN113049998B (zh) 多谐波激励作用下铁磁材料损耗预测方法
CN114048644B (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