CN114491992A - 一种基于等效偶极矩和物理光学法的高效电磁散射方法 - Google Patents

一种基于等效偶极矩和物理光学法的高效电磁散射方法 Download PDF

Info

Publication number
CN114491992A
CN114491992A CN202210026307.3A CN202210026307A CN114491992A CN 114491992 A CN114491992 A CN 114491992A CN 202210026307 A CN202210026307 A CN 202210026307A CN 114491992 A CN114491992 A CN 114491992A
Authority
CN
China
Prior art keywords
equivalent
current
magnetic
electric field
area
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.)
Granted
Application number
CN202210026307.3A
Other languages
English (en)
Other versions
CN114491992B (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202210026307.3A priority Critical patent/CN114491992B/zh
Publication of CN114491992A publication Critical patent/CN114491992A/zh
Application granted granted Critical
Publication of CN114491992B publication Critical patent/CN114491992B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Magnetic Variables (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于等效偶极矩和物理光学法的高效电磁散射方法，属于电磁场建模仿真技术领域。本发明所述方法将目标划分为矩量区(MoM区)和物理光学区(PO区)。MoM区满足表面电场积分方程，利用高阶阻抗边界条件构造表面等效电流和等效磁流的关系式，引入等效偶极矩方法加速计算；PO区利用物理光学法计算，同时将MoM区的等效电磁流以等效偶极矩的形式作为额外入射波；MoM区和PO区的互耦作用也利用等效偶极矩计算。本发明克服了积分过程中高斯求积阶数对计算速度的影响，采用等效偶极矩方法和高阶阻抗边界条件加速计算，适用范围广，可以用来分析各向同性或各向异性介质涂敷目标的电磁特性，具有高效、高精度、内存占用少等优点。

Description

一种基于等效偶极矩和物理光学法的高效电磁散射方法
技术领域
本发明属于电磁场建模仿真技术领域，具体涉及一种基于等效偶极矩和物理光学法的高效电磁散射方法。
背景技术
电磁场数值计算方法已经发展的相当成熟，然而，计算精度和效率之间的权衡仍然没有得到很好的解决。电磁场数值方法可分为低频方法和高频方法，低频方法如：矩量法(Method ofMoment，MoM)和有限元法(Finite Element Method，FEM)等，在处理包含精细结构的目标体时，可以保证相当高的精度；高频方法如：物理光学法(Physical Optics，PO)和几何光学法(Geometric Optics，GO)等，以牺牲计算精度为代价，可以极大提高计算效率。
随着目标电尺寸的增大和计算频点的增多，现有的数值方法所需的计算时间非常长，不能满足实际工程领域许多电磁问题高效分析的要求。因此需要提出新的改进方法，目前对电磁数值计算方法的改进主要体现在两个方面：一是改进算法本身，采用高低频混合方法，如：MoM-PO等；二是针对计算过程进行改进，如等效偶极矩法(Equivalent DipoleMoment，EDM)和快速多极子方法(Fast Multipole Method，FMM)等。
文献“An efficient hybrid method of iterative MoM-PO and equivalentdipole-moment for scattering from electrically large objects”公开了一种将等效偶极矩方法应用于迭代形式的MoM-PO中的方法，分析了PEC复合物体的电磁散射特性；文献“Efficient electromagnetic analysis for complex planar thin-layer compositeobjects by a hybrid method”公开了一种利用等效偶极矩方法加速贴片天线阵列电磁特性分析的方法。但目前所提出的等效偶极矩方法只分析了PEC物体的散射，对于各向同性和各向异性物体特性的分析，现有算法如MoM-PO也可以处理，但效率不高，尤其是电大尺寸目标，计算时间往往难以接受。
发明内容
本发明的目的是克服上述现有技术的缺陷，提供一种基于等效偶极矩和物理光学法的高效电磁散射方法，全面分析了PEC、各向同性和各向异性介质涂敷PEC物体的电磁散射特性，并且利用EDM方法大大减少CPU计算时间，利用高阶阻抗边界条件保证了较好的精度，并减少了计算机内存占用。
本发明所提出的技术问题是这样解决的：
一种基于等效偶极矩和物理光学法的高效电磁散射方法，利用电磁仿真软件将目标表面以三角形网格划分，根据目标表面曲率半径将目标表面划分为MoM区域和PO区域，该方法计算过程分为三个部分，第一部分用于处理MoM区域，第二部分用于计算PO区域的电磁流，第三部分用于计算MoM区域和PO区域的耦合作用。
本发明所述方法具体包括以下步骤：
步骤1.在电磁仿真软件中建立散射体模型，设置电磁波频率和三角形边长，利用电磁仿真软件对模型表面进行三角形网格划分，导出每个三角形序号和顶点坐标；
步骤2.记录每个三角形的中心坐标及面积，将相邻三角形组成一个三角面元对，编号并记录三角对的个数N，记录每对三角面元对的公共边长度；
步骤3.对模型进行分区，将曲率半径大于10个波长的区域划分为物理光学区即PO区，将其他区域划分为矩量区即MoM区；
步骤4.MoM区的处理采用矩量法计算；
步骤4-0.将散射体模型的表面等效电流和等效磁流用RWG基函数离散：
Figure BDA0003464049870000021
Figure BDA0003464049870000022
Figure BDA0003464049870000023
其中，J(r')为等效电流，M(r')为等效磁流，r'为散射体模型的表面任意点，1≤n≤N，N为基函数个数，fn(r')为第n个RWG基函数，an和bn分别为等效电流和等效磁流中第n个RWG基函数的展开系数；
Figure BDA0003464049870000024
表示第n个三角面元对的正负三角形，
Figure BDA0003464049870000025
表示第n个三角面元对正负三角形的面积，ln表示第n个三角面元对的公共边长度，
Figure BDA0003464049870000026
表示自由顶点和正负三角形中心所连成的向量；
步骤4-1.在MoM区域构建散射体模型表面的电场积分方程，并将用RWG基函数展开后的等效电流和等效磁流代入电场积分方程，得到离散后的电场积分方程；
步骤4-2.对离散后的电场积分方程等式两端采用伽辽金检验，即对方程两端同时用N个RWG基函数作内积，将方程个数扩展为N个，得到电场积分方程组；
步骤4-3.利用高阶阻抗边界条件，构造表面电场和磁场的关系式，进而得到等效电流和等效磁流之间的关系式，并将等效磁流用等效电流表示，代入电场积分方程组中；
步骤4-4.将电场积分方程组中的等效电流展开系数提取出来，改写为如下矩阵方程：
[Z][I]＝[V]
其中，[Z]表示离散后的电场积分方程左端经过伽辽金检验后得到的阻抗矩阵，[I]表示等效电流的展开系数向量，[I]＝[a1,a2,…,aN]T，[V]表示入射电场经过伽辽金检验后得到的电压向量；
步骤4-5.当伽辽金检验中的源点和场点间距大于0.15个波长时，利用等效偶极矩方法计算步骤4-4的阻抗矩阵和电压向量，当源点和场点间距小于0.15个波长时，则用矩量法计算步骤4-4的阻抗矩阵和电压向量；进而求解步骤4-4得到的矩阵方程，得到等效电流的展开系数；
步骤4-6.利用等效电流的展开系数计算等效电流，根据步骤4-3的等效电流和等效磁流之间的关系式，得到等效磁流及展开系数；
步骤5.物理光学区的处理采用物理光学法计算；
对PO区域采用物理光学法计算，在入射电场照射下，计算PO区产生的感应电流和感应磁流，并用RWG基函数展开；
步骤6.利用等效偶极矩方法计算MoM区和PO区的互耦作用；
步骤6-1.将步骤4-5得到的MoM区的等效电流和等效磁流看作等效电偶极子和等效磁偶极子，利用等效偶极矩方法分别计算等效电偶极子和等效磁偶极子产生的电场和磁场，并将电场和磁场作为PO区域耦合的入射场，利用物理光学法计算在该入射场下产生的耦合感应电流和耦合感应磁流，并用RWG基函数展开；
步骤6-2.将步骤6-1得到的耦合感应电流和步骤5得到的感应电流看作等效电偶极子，将步骤6-1得到的耦合感应磁流和步骤5得到的感应磁流看作等效磁偶极子，利用等效偶极矩方法分别计算等效电偶极子和等效磁偶极子产生的电场和磁场，将电场和磁场作为MoM区的修正入射场，对修正入射场进行伽辽金检验得到修正电压向量[ΔV]，并对步骤4-4的电压向量进行修正，即令：
[V]＝[V]+[ΔV]
执行步骤4-5至步骤4-6得到修正后的MoM区的等效电流和等效磁流的展开系数；
步骤6-3.判断由步骤6-2得到的修正后MoM区的等效电流展开系数和修正前的等效电流展开系数的相对误差，若相对误差小于阈值10-5，则判定修正后的等效电流和等效磁流即为MoM区的真实电流和真实磁流，步骤5得到的PO区产生的感应电流与步骤6-1得到的耦合感应电流之和为PO区的真实电流，步骤5得到的PO区产生的感应磁流与步骤6-1得到的耦合感应磁流之和为PO区的真实磁流；若相对误差大于10-5，则重复步骤6-1至步骤6-3。
进一步的，步骤1中，电磁波频率设置为300MHz。
进一步的，步骤1中，三角形边长设置为0.1个波长。
进一步的，步骤4-1的具体过程为：
在MoM区域构建散射体模型表面的电场积分方程：
Figure BDA0003464049870000041
Figure BDA0003464049870000042
Figure BDA0003464049870000043
其中，j为虚部符号，η0、k0和G分别表示波阻抗、波数和自由空间的格林函数，Et为散射体模型表面电场的切向分量，Einc表示入射场；s'为r'点所属三角面元对的区域，S为散射体模型表面区域，
Figure BDA0003464049870000044
·表示对J(r')求散度，
Figure BDA0003464049870000045
表示对G求梯度，
Figure BDA0003464049870000046
×表示求旋度；
将用RWG基函数展开后的等效电流和等效磁流代入电场积分方程，得到离散后的电场积分方程：
Figure BDA0003464049870000047
进一步的，步骤4-2的电场积分方程组为：
Figure BDA0003464049870000048
其中，fm(r)表示检验RWG基函数，r为fm(r)对应三角面元对的中心位置矢量，1≤m≤N。
进一步的，步骤4-3中，表面电场和磁场的关系式为：
Figure BDA0003464049870000049
Figure BDA00034640498700000410
Figure BDA00034640498700000411
等效电流和等效磁流之间的关系式：
Figure BDA0003464049870000051
L(Q)＝LD(Q)-LR(Q)
Figure BDA0003464049870000052
Figure BDA0003464049870000053
其中，Q表示向量变量，k'表示介质的波数，a'0、a'1和b'分别为由高阶阻抗边界条件确定的系数，
Figure BDA00034640498700000513
表示散射体模型表面的法向单位向量；下标t表示切向分量。
进一步的，步骤4-5中，利用等效偶极矩方法计算步骤4-4的阻抗矩阵的具体过程为：
将步骤4-4中的阻抗矩阵表示为：
Figure BDA0003464049870000054
其中，[ZJ]表示等效电流贡献的部分，[ZM]表示等效磁流的贡献，[ZE]表示Et经过伽辽金检验后形成的矩阵；
Figure BDA0003464049870000055
Figure BDA0003464049870000056
Figure BDA0003464049870000057
其中，
Figure BDA0003464049870000058
为矩阵[ZJ]第m行第n列元素，rmn表示由第n个RWG基函数中心指向第m个RWG基函数中心的向量，||表示取模，
Figure BDA0003464049870000059
和
Figure BDA00034640498700000510
分别表示第m个和第n个RWG基函数的等效电偶极矩，C＝(1+1/jk0|rmn|)/|rmn|2；
Figure BDA00034640498700000512
为矩阵[ZM]第m行第n列元素，对高阶阻抗边界条件的系数进行泰勒近似，得到
Figure BDA00034640498700000511
进一步的，修正电压向量表示为：
Figure BDA0003464049870000061
其中，
Figure BDA0003464049870000062
表示步骤5中PO区产生的感应电流中第i个基函数的展开系数，
Figure BDA0003464049870000063
表示步骤6-1中耦合感应电流中第i个基函数的展开系数，Np表示步骤5中基函数的数目，1≤i≤Np，
Figure BDA0003464049870000064
第i个RWG基函数的等效电偶极矩，rmi由MoM区第m个RWG基函数中心指向PO区第i个RWG基函数中心的向量，C'＝(1+1/jk0|rmi|)/|rmi|2。
进一步的，对于各向异性介质，由高阶阻抗边界条件，步骤4-3中，等效电流和等效磁流的关系式改写为：
Figure BDA0003464049870000065
Figure BDA0003464049870000066
其中，
Figure BDA0003464049870000067
和
Figure BDA0003464049870000068
分别是沿着设定各向异性方向的单位向量，kp是各向异性介质的波数，Jp(r')和Mp(r')分别表示沿各向异性方向的等效电流和等效磁流，
Figure BDA0003464049870000069
和b'p分别表示由高阶阻抗边界条件确定的各向异性介质系数；
对于p＝1和p＝2，分别执行以下步骤：
将等效磁流用等效电流表示，代入电场积分方程组中，执行步骤4-4到步骤6-3。
本发明的有益效果是：
本发明所述方法利用等效偶极矩方法，不仅加速了MoM区阻抗矩阵的填充，还加速了MoM和PO区的互耦计算，避免了高斯积分中求积阶数对计算速度的影响，大大提高了计算效率；等效偶极矩和物理光学法相结合，通过对目标表面的区域划分，使未知数个数大大减少，同时，引入高阶阻抗边界条件，构造目标表面电磁流的联系，进一步使未知数个数减半，在大幅提高计算效率的同时，还可以极大降低计算机内存占用。并且，该方法应用范围广，适用于PEC、各向同性或各向异性介质，对介质涂敷导体目标也同样适用。
附图说明
图1为RWG基函数示意图；
图2为等效电偶极矩和等效磁偶极矩示意图；
图3为MoM和PO区域划分示意图；
图4为介质涂敷导体球示意图；
图5为复合目标示意图。
具体实施方式
下面结合附图和实施例对本发明进行进一步的说明。
本实施例提供一种基于等效偶极矩和物理光学法的高效电磁散射方法，包括以下步骤：
步骤1.在电磁仿真软件中建立散射体模型，设置电磁波频率为300MHz，设置三角形边长为0.1个波长，利用电磁仿真软件对模型表面进行三角形网格划分，导出每个三角形序号和顶点坐标；
如图4和图5所示，本实施例的散射体模型选用各向异性介质涂敷导体球与各向异性正方形平板复合目标。PEC球体半径为0.48m，表面涂敷一层厚度为0.02m的各向异性介质，介质的相对介电系数为εuu＝4和εvv＝2-j。平板表面的阻抗特性为Zuu＝2η0和Zvv＝5η0，主轴偏转角α＝35°。
步骤2.记录每个三角形的中心坐标及面积，将相邻三角形组成一个三角面元对，编号并记录三角对的个数N，记录每对三角面元对的公共边长度；
步骤3.对模型进行分区，将曲率半径大于10个波长的区域划分为物理光学区(PO区)，将其他区域划分为矩量区(MoM区)；MoM和PO区域划分示意图如图3所示。
步骤4.MoM区的处理采用矩量法(MoM)计算；
将散射体模型的表面等效电流和磁流用RWG基函数离散：
Figure BDA0003464049870000071
Figure BDA0003464049870000072
Figure BDA0003464049870000073
其中，J(r')为等效电流，M(r')为等效磁流，r'为散射体模型的表面任意点，1≤n≤N，N为基函数个数(与三角对的个数相同)，fn(r')为第n个RWG基函数，an和bn分别为等效电流和等效磁流中第n个RWG基函数的展开系数；
Figure BDA0003464049870000081
表示第n个三角面元对的正负三角形，
Figure BDA0003464049870000082
表示第n个三角面元对正负三角形的面积，ln表示第n个三角面元对的公共边长度，
Figure BDA0003464049870000083
表示自由顶点和正负三角形中心所连成的向量；RWG基函数示意图如图1所示。
步骤4-1.在MoM区域构建散射体模型表面的电场积分方程：
Figure BDA0003464049870000084
Figure BDA0003464049870000085
Figure BDA0003464049870000086
其中，j为虚部符号，η0、k0和G分别表示波阻抗、波数和自由空间的格林函数，Et为散射体模型表面电场的切向分量，Einc表示入射场；s'为r'点所属三角面元对的区域，S为散射体模型表面区域，
Figure BDA0003464049870000087
·表示对J(r')求散度，
Figure BDA0003464049870000088
表示对G求梯度，
Figure BDA0003464049870000089
×表示求旋度；
将用RWG基函数展开后的等效电流和等效磁流代入电场积分方程，得到离散后的电场积分方程：
Figure BDA00034640498700000810
步骤4-2.对离散后的电场积分方程等式两端采用伽辽金检验，即对方程两端同时用N个RWG基函数作内积，将方程个数扩展为N个，得到电场积分方程组：
Figure BDA00034640498700000811
其中，fm(r)表示检验RWG基函数，r为fm(r)对应三角面元对的中心位置矢量，1≤m≤N。
步骤4-3.利用高阶阻抗边界条件(HOIBC)，构造表面电场和磁场的关系式：
Figure BDA00034640498700000812
Figure BDA00034640498700000813
Figure BDA00034640498700000814
进而得到等效电流和等效磁流之间的关系式：
Figure BDA0003464049870000091
L(Q)＝LD(Q)-LR(Q)
Figure BDA0003464049870000092
Figure BDA0003464049870000093
其中，Q表示向量变量，k'表示介质的波数，a'0、a'1和b'分别为由高阶阻抗边界条件确定的系数，
Figure BDA0003464049870000094
表示散射体模型表面的法向单位向量；下标t表示切向分量；
将等效磁流用等效电流表示，代入电场积分方程组中；
步骤4-4.通过整理步骤4-3得到的方程组，可将方程组中未知数个数减半，并将等效电流系数提取出来，进而方程组改写为如下形式：
[Z][I]＝[V]
其中，[I]表示等效电流的展开系数向量，[I]＝[a1,a2,…,aN]T，[V]表示入射电场经过伽辽金检验后得到的电压向量；[Z]为离散后的电场积分方程左端经过伽辽金检验后得到的阻抗矩阵，表示为：
Figure BDA0003464049870000095
其中，[ZJ]表示等效电流贡献的部分，[ZM]表示等效磁流的贡献，[ZE]表示Et经过伽辽金检验后形成的矩阵；
步骤4-5.当伽辽金检验中的源点和场点间距大于0.15个波长时，利用等效偶极矩方法计算步骤4-4的阻抗矩阵和电压向量，可以避免高斯积分计算中求积阶数的影响；
Figure BDA0003464049870000096
Figure BDA0003464049870000097
Figure BDA0003464049870000098
其中，
Figure BDA0003464049870000099
为矩阵[ZJ]第m行第n列元素，rmn表示由第n个RWG基函数中心指向第m个RWG基函数中心的向量，||表示取模，
Figure BDA0003464049870000101
和
Figure BDA0003464049870000102
分别表示第m个和第n个RWG基函数的等效电偶极矩，C＝(1+1/jk0|rmn|)/|rmn|2；
Figure BDA0003464049870000103
为矩阵[ZM]第m行第n列元素，对高阶阻抗边界条件的系数进行泰勒近似，得到
Figure BDA0003464049870000104
等效电偶极矩和等效磁偶极矩示意图如图2所示。
当源点和场点间距小于0.15个波长时，则用矩量法计算步骤4-4的阻抗矩阵和电压向量；进而求解步骤4-4得到的矩阵方程，得到等效电流的展开系数；
步骤4-6.利用等效电流的展开系数计算等效电流，根据步骤4-3的等效电流和等效磁流之间的关系式，得到等效磁流及展开系数。
步骤5.物理光学区的处理采用物理光学法(PO)计算；
对PO区域采用物理光学法计算，在入射电场照射下，计算PO区产生的感应电流和感应磁流，并用RWG基函数展开；
步骤6.利用等效偶极矩方法(EDM)计算MoM区和PO区的互耦作用；
步骤6-1.将步骤4-5得到的MoM区的等效电流和等效磁流看作等效电偶极子和等效磁偶极子，利用等效偶极矩方法分别计算等效电偶极子和等效磁偶极子产生的电场和磁场，并将电场和磁场作为PO区域耦合的入射场，利用物理光学法计算在该入射场下产生的耦合感应电流和耦合感应磁流，并用RWG基函数展开；
步骤6-2.将步骤6-1得到的耦合感应电流和步骤5得到的感应电流看作等效电偶极子，将步骤6-1得到的耦合感应磁流和步骤5得到的感应磁流看作等效磁偶极子，利用等效偶极矩方法分别计算等效电偶极子和等效磁偶极子产生的电场和磁场，将电场和磁场作为MoM区的修正入射场，对修正入射场进行伽辽金检验得到修正电压向量[ΔV]：
Figure BDA0003464049870000105
其中，
Figure BDA0003464049870000106
表示步骤5中PO区产生的感应电流中第i个基函数的展开系数，
Figure BDA0003464049870000107
表示步骤6-1中耦合感应电流中第i个基函数的展开系数，Np表示步骤5中基函数的数目，1≤i≤Np，
Figure BDA0003464049870000108
第i个RWG基函数的等效电偶极矩，rmi由MoM区第m个RWG基函数中心指向PO区第i个RWG基函数中心的向量，C'＝(1+1/jk0|rmi|)/|rmi|2；
对步骤4-4的电压向量进行修正，即令：
[V]＝[V]+[ΔV]
重复步骤4-5至步骤4-6得到修正后的MoM区的等效电流和等效磁流的展开系数。
步骤6-3.判断由步骤6-2得到的修正后MoM区的等效电流展开系数和修正前的等效电流展开系数的相对误差，若相对误差小于设定阈值10-5，则判定修正后的等效电流和等效磁流即为MoM区的真实电流和磁流，步骤5得到的PO区产生的感应电流与步骤6-1得到的耦合感应电流之和为PO区的真实电流，步骤5得到的PO区产生的感应磁流与步骤6-1得到的耦合感应磁流之和为PO区的真实磁流。若相对误差大于10-5，则重复步骤6-1至步骤6-3。通常，只需要通过5-10次迭代就能使相对误差小于阈值，从而终止循环。
对于各向异性介质，由高阶阻抗边界条件，步骤4-3中，等效电流和等效磁流的关系式改写为：
Figure BDA0003464049870000111
Figure BDA0003464049870000112
其中，
Figure BDA0003464049870000113
和
Figure BDA0003464049870000114
分别是沿着设定各向异性方向的单位向量，kp是各向异性介质的波数，Jp(r')和Mp(r')分别表示沿各向异性方向的等效电流和等效磁流，
Figure BDA0003464049870000115
和b'p分别表示由高阶阻抗边界条件确定的各向异性介质系数。
对于p＝1和p＝2，分别执行以下步骤：
将等效磁流用等效电流表示，代入电场积分方程组中，执行步骤4-4到步骤6-3。
以上所述，仅为本发明的具体实施方式，本说明书中所公开的任一特征，除非特别叙述，均可被其他等效或具有类似目的的替代特征加以替换；所公开的所有特征、或所有方法或过程中的步骤，除了互相排斥的特征和/或步骤以外，均可以任何方式组合；本领域的技术人员根据本发明技术方案的技术特征所做出的任何非本质的添加、替换，均属于本发明的保护范围。

Claims (9)

1.一种基于等效偶极矩和物理光学法的高效电磁散射方法，其特征在于，包括以下步骤：
步骤1.在电磁仿真软件中建立散射体模型，设置电磁波频率和三角形边长，利用电磁仿真软件对模型表面进行三角形网格划分，导出每个三角形序号和顶点坐标；
步骤2.记录每个三角形的中心坐标及面积，将相邻三角形组成一个三角面元对，编号并记录三角对的个数N，记录每对三角面元对的公共边长度；
步骤3.对模型进行分区，将曲率半径大于10个波长的区域划分为物理光学区即PO区，将其他区域划分为矩量区即MoM区；
步骤4.MoM区的处理采用矩量法计算；
步骤4-0.将散射体模型的表面等效电流和等效磁流用RWG基函数离散：
Figure FDA0003464049860000011
Figure FDA0003464049860000012
Figure FDA0003464049860000013
其中，J(r')为等效电流，M(r')为等效磁流，r'为散射体模型的表面任意点，1≤n≤N，N为基函数个数，fn(r')为第n个RWG基函数，an和bn分别为等效电流和等效磁流中第n个RWG基函数的展开系数；
Figure FDA0003464049860000014
表示第n个三角面元对的正负三角形，
Figure FDA0003464049860000015
表示第n个三角面元对正负三角形的面积，ln表示第n个三角面元对的公共边长度，
Figure FDA0003464049860000016
表示自由顶点和正负三角形中心所连成的向量；
步骤4-1.在MoM区域构建散射体模型表面的电场积分方程，并将用RWG基函数展开后的等效电流和等效磁流代入电场积分方程，得到离散后的电场积分方程；
步骤4-2.对离散后的电场积分方程等式两端采用伽辽金检验，即对方程两端同时用N个RWG基函数作内积，将方程个数扩展为N个，得到电场积分方程组；
步骤4-3.利用高阶阻抗边界条件，构造表面电场和磁场的关系式，进而得到等效电流和等效磁流之间的关系式，并将等效磁流用等效电流表示，代入电场积分方程组中；
步骤4-4.将电场积分方程组中的等效电流展开系数提取出来，改写为如下矩阵方程：
[Z][I]＝[V]
其中，[Z]表示离散后的电场积分方程左端经过伽辽金检验后得到的阻抗矩阵，[I]表示等效电流的展开系数向量，[I]＝[a1,a2,…,aN]T，[V]表示入射电场经过伽辽金检验后得到的电压向量；
步骤4-5.当伽辽金检验中的源点和场点间距大于0.15个波长时，利用等效偶极矩方法计算步骤4-4的阻抗矩阵和电压向量，当源点和场点间距小于0.15个波长时，则用矩量法计算步骤4-4的阻抗矩阵和电压向量；进而求解步骤4-4得到的矩阵方程，得到等效电流的展开系数；
步骤4-6.利用等效电流的展开系数计算等效电流，根据步骤4-3的等效电流和等效磁流之间的关系式，得到等效磁流及展开系数；
步骤5.物理光学区的处理采用物理光学法计算；
对PO区域采用物理光学法计算，在入射电场照射下，计算PO区产生的感应电流和感应磁流，并用RWG基函数展开；
步骤6.利用等效偶极矩方法计算MoM区和PO区的互耦作用；
步骤6-1.将步骤4-5得到的MoM区的等效电流和等效磁流看作等效电偶极子和等效磁偶极子，利用等效偶极矩方法分别计算等效电偶极子和等效磁偶极子产生的电场和磁场，并将电场和磁场作为PO区域耦合的入射场，利用物理光学法计算在该入射场下产生的耦合感应电流和耦合感应磁流，并用RWG基函数展开；
步骤6-2.将步骤6-1得到的耦合感应电流和步骤5得到的感应电流看作等效电偶极子，将步骤6-1得到的耦合感应磁流和步骤5得到的感应磁流看作等效磁偶极子，利用等效偶极矩方法分别计算等效电偶极子和等效磁偶极子产生的电场和磁场，将电场和磁场作为MoM区的修正入射场，对修正入射场进行伽辽金检验得到修正电压向量[ΔV]，并对步骤4-4的电压向量进行修正，即令：
[V]＝[V]+[ΔV]
执行步骤4-5至步骤4-6得到修正后的MoM区的等效电流和等效磁流的展开系数；
步骤6-3.判断由步骤6-2得到的修正后MoM区的等效电流展开系数和修正前的等效电流展开系数的相对误差，若相对误差小于阈值10-5，则判定修正后的等效电流和等效磁流即为MoM区的真实电流和真实磁流，步骤5得到的PO区产生的感应电流与步骤6-1得到的耦合感应电流之和为PO区的真实电流，步骤5得到的PO区产生的感应磁流与步骤6-1得到的耦合感应磁流之和为PO区的真实磁流；若相对误差大于10-5，则重复步骤6-1至步骤6-3。
2.根据权利要求1所述的基于等效偶极矩和物理光学法的高效电磁散射方法，其特征在于，步骤1中，电磁波频率设置为300MHz。
3.根据权利要求1所述的基于等效偶极矩和物理光学法的高效电磁散射方法，其特征在于，步骤1中，三角形边长设置为0.1个波长。
4.根据权利要求1所述的基于等效偶极矩和物理光学法的高效电磁散射方法，其特征在于，步骤4-1的具体过程为：
在MoM区域构建散射体模型表面的电场积分方程：
Figure FDA0003464049860000031
Figure FDA0003464049860000032
Figure FDA0003464049860000033
其中，j为虚部符号，η0、k0和G分别表示波阻抗、波数和自由空间的格林函数，Et为散射体模型表面电场的切向分量，Einc表示入射场；s'为r'点所属三角面元对的区域，S为散射体模型表面区域，
Figure FDA0003464049860000034
表示对J(r')求散度，
Figure FDA0003464049860000035
表示对G求梯度，
Figure FDA0003464049860000036
表示求旋度；
将用RWG基函数展开后的等效电流和等效磁流代入电场积分方程，得到离散后的电场积分方程：
Figure FDA0003464049860000037
5.根据权利要求4所述的基于等效偶极矩和物理光学法的高效电磁散射方法，其特征在于，步骤4-2的电场积分方程组为：
Figure FDA0003464049860000038
其中，fm(r)表示检验RWG基函数，r为fm(r)对应三角面元对的中心位置矢量，1≤m≤N。
6.根据权利要求5所述的基于等效偶极矩和物理光学法的高效电磁散射方法，其特征在于，步骤4-3中，表面电场和磁场的关系式为：
Figure FDA0003464049860000039
Figure FDA0003464049860000041
Figure FDA0003464049860000042
等效电流和等效磁流之间的关系式：
Figure FDA0003464049860000043
L(Q)＝LD(Q)-LR(Q)
Figure FDA0003464049860000044
Figure FDA0003464049860000045
其中，Q表示向量变量，k'表示介质的波数，a'0、a'1和b'分别为由高阶阻抗边界条件确定的系数，
Figure FDA0003464049860000046
表示散射体模型表面的法向单位向量；下标t表示切向分量。
7.根据权利要求6所述的基于等效偶极矩和物理光学法的高效电磁散射方法，其特征在于，步骤4-5中，利用等效偶极矩方法计算步骤4-4的阻抗矩阵的具体过程为：
将步骤4-4中的阻抗矩阵表示为：
Figure FDA0003464049860000047
其中，[ZJ]表示等效电流贡献的部分，[ZM]表示等效磁流的贡献，[ZE]表示Et经过伽辽金检验后形成的矩阵；
Figure FDA0003464049860000048
Figure FDA0003464049860000049
Figure FDA00034640498600000410
其中，
Figure FDA00034640498600000411
为矩阵[ZJ]第m行第n列元素，rmn表示由第n个RWG基函数中心指向第m个RWG基函数中心的向量，| |表示取模，
Figure FDA00034640498600000412
和
Figure FDA00034640498600000413
分别表示第m个和第n个RWG基函数的等效电偶极矩，C＝(1+1/jk0|rmn|)/|rmn|2；
Figure FDA00034640498600000414
为矩阵[ZM]第m行第n列元素，对高阶阻抗边界条件的系数进行泰勒近似，得到
Figure FDA0003464049860000051
8.根据权利要求7所述的基于等效偶极矩和物理光学法的高效电磁散射方法，其特征在于，修正电压向量表示为：
Figure FDA0003464049860000052
其中，
Figure FDA0003464049860000053
Figure FDA0003464049860000054
表示步骤5中PO区产生的感应电流中第i个基函数的展开系数，
Figure FDA0003464049860000055
表示步骤6-1中耦合感应电流中第i个基函数的展开系数，Np表示步骤5中基函数的数目，
Figure FDA0003464049860000056
Figure FDA0003464049860000057
第i个RWG基函数的等效电偶极矩，rmi由MoM区第m个RWG基函数中心指向PO区第i个RWG基函数中心的向量，C'＝(1+1/jk0|rmi|)/|rmi|2。
9.根据权利要求8所述的基于等效偶极矩和物理光学法的高效电磁散射方法，其特征在于，对于各向异性介质，由高阶阻抗边界条件，步骤4-3中，等效电流和等效磁流的关系式改写为：
Figure FDA0003464049860000058
Figure FDA0003464049860000059
其中，
Figure FDA00034640498600000510
和
Figure FDA00034640498600000511
分别是沿着设定各向异性方向的单位向量，kp是各向异性介质的波数，Jp(r')和Mp(r')分别表示沿各向异性方向的等效电流和等效磁流，
Figure FDA00034640498600000512
和b'p分别表示由高阶阻抗边界条件确定的各向异性介质系数；
对于p＝1和p＝2，分别执行以下步骤：
将等效磁流用等效电流表示，代入电场积分方程组中，执行步骤4-4到步骤6-3。
CN202210026307.3A 2022-01-11 2022-01-11 一种基于等效偶极矩和物理光学法的高效电磁散射方法 Active CN114491992B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210026307.3A CN114491992B (zh) 2022-01-11 2022-01-11 一种基于等效偶极矩和物理光学法的高效电磁散射方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210026307.3A CN114491992B (zh) 2022-01-11 2022-01-11 一种基于等效偶极矩和物理光学法的高效电磁散射方法

Publications (2)

Publication Number Publication Date
CN114491992A true CN114491992A (zh) 2022-05-13
CN114491992B CN114491992B (zh) 2024-05-24

Family

ID=81509792

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210026307.3A Active CN114491992B (zh) 2022-01-11 2022-01-11 一种基于等效偶极矩和物理光学法的高效电磁散射方法

Country Status (1)

Country Link
CN (1) CN114491992B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115184689A (zh) * 2022-05-26 2022-10-14 西安电子科技大学 一种计算蒸发波导环境中目标电磁散射的方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1404415A (zh) * 2000-02-22 2003-03-19 基因谱公司 微阵列制造技术及设备
US20130317337A1 (en) * 2012-05-23 2013-11-28 Chenyu Wu Continuous Modeling for Dipole Localization from 2D MCG Images with Unknown Depth
CN104915326A (zh) * 2014-03-14 2015-09-16 南京理工大学 基于等效原理的区域分解阶数步进时域积分方法
CN105929253A (zh) * 2016-06-15 2016-09-07 北京航空航天大学 采用偶极矩模型对pcb电路电磁场进行反演的方法
CN110765412A (zh) * 2019-09-26 2020-02-07 哈尔滨工程大学 质心切分的小波矩量法求电大物体电磁散射的方法
CN111029760A (zh) * 2019-12-07 2020-04-17 复旦大学 基片集成波导可编程超材料天线方向图估计方法
CN112528580A (zh) * 2020-11-19 2021-03-19 广州大学 一种反激式变换器电路板的电磁辐射仿真预测方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1404415A (zh) * 2000-02-22 2003-03-19 基因谱公司 微阵列制造技术及设备
US20130317337A1 (en) * 2012-05-23 2013-11-28 Chenyu Wu Continuous Modeling for Dipole Localization from 2D MCG Images with Unknown Depth
CN104915326A (zh) * 2014-03-14 2015-09-16 南京理工大学 基于等效原理的区域分解阶数步进时域积分方法
CN105929253A (zh) * 2016-06-15 2016-09-07 北京航空航天大学 采用偶极矩模型对pcb电路电磁场进行反演的方法
CN110765412A (zh) * 2019-09-26 2020-02-07 哈尔滨工程大学 质心切分的小波矩量法求电大物体电磁散射的方法
CN111029760A (zh) * 2019-12-07 2020-04-17 复旦大学 基片集成波导可编程超材料天线方向图估计方法
CN112528580A (zh) * 2020-11-19 2021-03-19 广州大学 一种反激式变换器电路板的电磁辐射仿真预测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HUI ZHOU 等: "Efficient EDM-PO Method for the Scattering From Electrically Large Objects With the High-Order Impedance Boundary Condition", 《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》, vol. 70, no. 9, 30 May 2022 (2022-05-30), pages 8242 - 8249, XP011922678, DOI: 10.1109/TAP.2022.3177560 *
LI XIAO 等: "An Efficient Hybrid Method of Iterative MoM-PO and Equivalent Dipole-Moment for Scattering From Electrically Large Objects", 《IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》, vol. 16, 15 February 2017 (2017-02-15), pages 1723 - 1726, XP011653113, DOI: 10.1109/LAWP.2017.2669910 *
郁超 等: "介质体电磁散射分析中的等效偶极矩法", 《微波学报》, vol. 26, no. 2, 15 April 2010 (2010-04-15), pages 26 - 29 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115184689A (zh) * 2022-05-26 2022-10-14 西安电子科技大学 一种计算蒸发波导环境中目标电磁散射的方法
CN115184689B (zh) * 2022-05-26 2024-03-19 西安电子科技大学 一种计算蒸发波导环境中目标电磁散射的方法

Also Published As

Publication number Publication date
CN114491992B (zh) 2024-05-24

Similar Documents

Publication Publication Date Title
CN109765538B (zh) 非均匀介质目标体的散射场确定方法
CN113158527B (zh) 一种基于隐式fvfd计算频域电磁场的方法
CN111046603A (zh) 基于gpu并行加速特征基函数算法的电磁散射特性分析方法
Cai et al. Volume surface integral equation method based on higher order hierarchical vector basis functions for EM scattering and radiation from composite metallic and dielectric structures
CN107515956A (zh) 一种基于hfss单元法的大型有限平面阵列分析方法
CN114491992A (zh) 一种基于等效偶极矩和物理光学法的高效电磁散射方法
CN111144013B (zh) 高精度介质体目标散射的仿真方法
CN110737873A (zh) 一种大规模阵列天线散射的快速分析方法
CN106294894A (zh) 快速分析非均匀目标电磁散射特性的有限元边界积分方法
He et al. A multihybrid FE-BI-KA technique for 3-D electromagnetic scattering from a coated object above a conductive rough surface
Cai et al. Singularity treatment and high‐order RWG basis functions for integral equations of electromagnetic scattering
Shi et al. Comparison of interpolating functions and interpolating points in full-wave multilevel Green's function interpolation method
Rong et al. Fast direct surface integral equation solution for electromagnetic scattering analysis with skeletonization factorization
CN116401921B (zh) 一种各项异性磁化等离子体媒质处理方法及系统
US7359929B2 (en) Fast solution of integral equations representing wave propagation
Zhou et al. Efficient EDM-PO method for the scattering from electrically large objects with the high-order impedance boundary condition
CN115017782A (zh) 考虑介质各向异性的三维天然源电磁场计算方法
CN109657196B (zh) 一种基于空域强耦合路径加速的mom-ka混合方法
Guo et al. Analysis of antenna array with arbitrarily shaped radomes using fast algorithm based on VSIE
Wang et al. Novel reduced matrix equation constructing method accelerates iterative solution of characteristic basis function method
CN110765412A (zh) 质心切分的小波矩量法求电大物体电磁散射的方法
Xiang et al. The application of Barycentric subdivision method for numerical integration in method of moments
Xiang et al. Application of barycentric subdivision method for singularity integration in method of moments
Gu et al. Parallel Integral Equation Based Non-overlapping DDM for Fast Solving Electromagnetic Scattering Problems with Changeable Parts
Wei et al. Domain decomposition WLP-FDTD for efficient electromagnetic scattering analysis

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