CN114491992B - 一种基于等效偶极矩和物理光学法的高效电磁散射方法 - Google Patents
一种基于等效偶极矩和物理光学法的高效电磁散射方法 Download PDFInfo
- Publication number
- CN114491992B CN114491992B CN202210026307.3A CN202210026307A CN114491992B CN 114491992 B CN114491992 B CN 114491992B CN 202210026307 A CN202210026307 A CN 202210026307A CN 114491992 B CN114491992 B CN 114491992B
- Authority
- CN
- China
- Prior art keywords
- equivalent
- current
- area
- magnetic
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 90
- 230000003287 optical effect Effects 0.000 title claims abstract description 31
- 238000000790 scattering method Methods 0.000 title claims abstract description 9
- 230000005684 electric field Effects 0.000 claims abstract description 61
- 230000008878 coupling Effects 0.000 claims abstract description 25
- 238000010168 coupling process Methods 0.000 claims abstract description 25
- 238000005859 coupling reaction Methods 0.000 claims abstract description 25
- 230000008569 process Effects 0.000 claims abstract description 8
- 230000006870 function Effects 0.000 claims description 61
- 239000013598 vector Substances 0.000 claims description 41
- 239000011159 matrix material Substances 0.000 claims description 28
- 230000006698 induction Effects 0.000 claims description 14
- 230000001808 coupling effect Effects 0.000 claims description 4
- 238000000638 solvent extraction Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 19
- 230000005672 electromagnetic field Effects 0.000 abstract description 4
- 239000011248 coating agent Substances 0.000 abstract description 2
- 238000000576 coating method Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 5
- 239000002131 composite material Substances 0.000 description 4
- 239000004020 conductor Substances 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000005405 multipole Effects 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measuring Magnetic Variables (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基函数离散:
其中,J(r')为等效电流,M(r')为等效磁流,r'为散射体模型的表面任意点,1≤n≤N,N为基函数个数,fn(r')为第n个RWG基函数,an和bn分别为等效电流和等效磁流中第n个RWG基函数的展开系数;表示第n个三角面元对的正负三角形,/>表示第n个三角面元对正负三角形的面积,ln表示第n个三角面元对的公共边长度,/>表示自由顶点和正负三角形中心所连成的向量;
步骤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区域构建散射体模型表面的电场积分方程:
其中,j为虚部符号,η0、k0和G分别表示波阻抗、波数和自由空间的格林函数,Et为散射体模型表面电场的切向分量,Einc表示入射场;s'为r'点所属三角面元对的区域,S为散射体模型表面区域,·表示对J(r')求散度,/>表示对G求梯度,/>×表示求旋度;
将用RWG基函数展开后的等效电流和等效磁流代入电场积分方程,得到离散后的电场积分方程:
进一步的,步骤4-2的电场积分方程组为:
其中,fm(r)表示检验RWG基函数,r为fm(r)对应三角面元对的中心位置矢量,1≤m≤N。
进一步的,步骤4-3中,表面电场和磁场的关系式为:
等效电流和等效磁流之间的关系式:
L(Q)=LD(Q)-LR(Q)
其中,Q表示向量变量,k'表示介质的波数,a'0、a'1和b'分别为由高阶阻抗边界条件确定的系数,表示散射体模型表面的法向单位向量;下标t表示切向分量。
进一步的,步骤4-5中,利用等效偶极矩方法计算步骤4-4的阻抗矩阵的具体过程为:
将步骤4-4中的阻抗矩阵表示为:
其中,[ZJ]表示等效电流贡献的部分,[ZM]表示等效磁流的贡献,[ZE]表示Et经过伽辽金检验后形成的矩阵;
其中,为矩阵[ZJ]第m行第n列元素,rmn表示由第n个RWG基函数中心指向第m个RWG基函数中心的向量,||表示取模,/>和/>分别表示第m个和第n个RWG基函数的等效电偶极矩,C=(1+1/jk0|rmn|)/|rmn|2;/>为矩阵[ZM]第m行第n列元素,对高阶阻抗边界条件的系数进行泰勒近似,得到/>
进一步的,修正电压向量表示为:
其中,表示步骤5中PO区产生的感应电流中第i个基函数的展开系数,/>表示步骤6-1中耦合感应电流中第i个基函数的展开系数,Np表示步骤5中基函数的数目,1≤i≤Np,/>第i个RWG基函数的等效电偶极矩,rmi由MoM区第m个RWG基函数中心指向PO区第i个RWG基函数中心的向量,C'=(1+1/jk0|rmi|)/|rmi|2。
进一步的,对于各向异性介质,由高阶阻抗边界条件,步骤4-3中,等效电流和等效磁流的关系式改写为:
其中,和/>分别是沿着设定各向异性方向的单位向量,kp是各向异性介质的波数,Jp(r')和Mp(r')分别表示沿各向异性方向的等效电流和等效磁流,/>和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基函数离散:
其中,J(r')为等效电流,M(r')为等效磁流,r'为散射体模型的表面任意点,1≤n≤N,N为基函数个数(与三角对的个数相同),fn(r')为第n个RWG基函数,an和bn分别为等效电流和等效磁流中第n个RWG基函数的展开系数;表示第n个三角面元对的正负三角形,表示第n个三角面元对正负三角形的面积,ln表示第n个三角面元对的公共边长度,/>表示自由顶点和正负三角形中心所连成的向量;RWG基函数示意图如图1所示。
步骤4-1.在MoM区域构建散射体模型表面的电场积分方程:
其中,j为虚部符号,η0、k0和G分别表示波阻抗、波数和自由空间的格林函数,Et为散射体模型表面电场的切向分量,Einc表示入射场;s'为r'点所属三角面元对的区域,S为散射体模型表面区域,·表示对J(r')求散度,/>表示对G求梯度,/>×表示求旋度;
将用RWG基函数展开后的等效电流和等效磁流代入电场积分方程,得到离散后的电场积分方程:
步骤4-2.对离散后的电场积分方程等式两端采用伽辽金检验,即对方程两端同时用N个RWG基函数作内积,将方程个数扩展为N个,得到电场积分方程组:
其中,fm(r)表示检验RWG基函数,r为fm(r)对应三角面元对的中心位置矢量,1≤m≤N。
步骤4-3.利用高阶阻抗边界条件(HOIBC),构造表面电场和磁场的关系式:
进而得到等效电流和等效磁流之间的关系式:
L(Q)=LD(Q)-LR(Q)
其中,Q表示向量变量,k'表示介质的波数,a'0、a'1和b'分别为由高阶阻抗边界条件确定的系数,表示散射体模型表面的法向单位向量;下标t表示切向分量;
将等效磁流用等效电流表示,代入电场积分方程组中;
步骤4-4.通过整理步骤4-3得到的方程组,可将方程组中未知数个数减半,并将等效电流系数提取出来,进而方程组改写为如下形式:
[Z][I]=[V]
其中,[I]表示等效电流的展开系数向量,[I]=[a1,a2,…,aN]T,[V]表示入射电场经过伽辽金检验后得到的电压向量;[Z]为离散后的电场积分方程左端经过伽辽金检验后得到的阻抗矩阵,表示为:
其中,[ZJ]表示等效电流贡献的部分,[ZM]表示等效磁流的贡献,[ZE]表示Et经过伽辽金检验后形成的矩阵;
步骤4-5.当伽辽金检验中的源点和场点间距大于0.15个波长时,利用等效偶极矩方法计算步骤4-4的阻抗矩阵和电压向量,可以避免高斯积分计算中求积阶数的影响;
其中,为矩阵[ZJ]第m行第n列元素,rmn表示由第n个RWG基函数中心指向第m个RWG基函数中心的向量,||表示取模,/>和/>分别表示第m个和第n个RWG基函数的等效电偶极矩,C=(1+1/jk0|rmn|)/|rmn|2;/>为矩阵[ZM]第m行第n列元素,对高阶阻抗边界条件的系数进行泰勒近似,得到/>等效电偶极矩和等效磁偶极矩示意图如图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]:
其中,表示步骤5中PO区产生的感应电流中第i个基函数的展开系数,/>表示步骤6-1中耦合感应电流中第i个基函数的展开系数,Np表示步骤5中基函数的数目,1≤i≤Np,/>第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中,等效电流和等效磁流的关系式改写为:
其中,和/>分别是沿着设定各向异性方向的单位向量,kp是各向异性介质的波数,Jp(r')和Mp(r')分别表示沿各向异性方向的等效电流和等效磁流,/>和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基函数离散:
其中,J(r')为等效电流,M(r')为等效磁流,r'为散射体模型的表面任意点,1≤n≤N,N为基函数个数,fn(r')为第n个RWG基函数,an和bn分别为等效电流和等效磁流中第n个RWG基函数的展开系数;表示第n个三角面元对的正负三角形,/>表示第n个三角面元对正负三角形的面积,ln表示第n个三角面元对的公共边长度,/>表示自由顶点和正负三角形中心所连成的向量;
步骤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区域构建散射体模型表面的电场积分方程:
其中,j为虚部符号,η0、k0和G分别表示波阻抗、波数和自由空间的格林函数,Et为散射体模型表面电场的切向分量,Einc表示入射场;s'为r'点所属三角面元对的区域,S为散射体模型表面区域,表示对J(r')求散度,/>表示对G求梯度,/>表示求旋度;
将用RWG基函数展开后的等效电流和等效磁流代入电场积分方程,得到离散后的电场积分方程:
5.根据权利要求4所述的基于等效偶极矩和物理光学法的高效电磁散射方法,其特征在于,步骤4-2的电场积分方程组为:
其中,fm(r)表示检验RWG基函数,r为fm(r)对应三角面元对的中心位置矢量,1≤m≤N。
6.根据权利要求5所述的基于等效偶极矩和物理光学法的高效电磁散射方法,其特征在于,步骤4-3中,表面电场和磁场的关系式为:
等效电流和等效磁流之间的关系式:
L(Q)=LD(Q)-LR(Q)
其中,Q表示向量变量,k'表示介质的波数,a'0、a'1和b'分别为由高阶阻抗边界条件确定的系数,表示散射体模型表面的法向单位向量;下标t表示切向分量。
7.根据权利要求6所述的基于等效偶极矩和物理光学法的高效电磁散射方法,其特征在于,步骤4-5中,利用等效偶极矩方法计算步骤4-4的阻抗矩阵的具体过程为:
将步骤4-4中的阻抗矩阵表示为:
其中,[ZJ]表示等效电流贡献的部分,[ZM]表示等效磁流的贡献,[ZE]表示Et经过伽辽金检验后形成的矩阵;
其中,为矩阵[ZJ]第m行第n列元素,rmn表示由第n个RWG基函数中心指向第m个RWG基函数中心的向量,| |表示取模,/>和/>分别表示第m个和第n个RWG基函数的等效电偶极矩,C=(1+1/jk0|rmn|)/|rmn|2;/>为矩阵[ZM]第m行第n列元素,对高阶阻抗边界条件的系数进行泰勒近似,得到/>
8.根据权利要求7所述的基于等效偶极矩和物理光学法的高效电磁散射方法,其特征在于,修正电压向量表示为:
其中, 表示步骤5中PO区产生的感应电流中第i个基函数的展开系数,/>表示步骤6-1中耦合感应电流中第i个基函数的展开系数,Np表示步骤5中基函数的数目,/> 第i个RWG基函数的等效电偶极矩,rmi由MoM区第m个RWG基函数中心指向PO区第i个RWG基函数中心的向量,C'=(1+1/jk0|rmi|)/|rmi|2。
9.根据权利要求8所述的基于等效偶极矩和物理光学法的高效电磁散射方法,其特征在于,对于各向异性介质,由高阶阻抗边界条件,步骤4-3中,等效电流和等效磁流的关系式改写为:
其中,和/>分别是沿着设定各向异性方向的单位向量,kp是各向异性介质的波数,Jp(r')和Mp(r')分别表示沿各向异性方向的等效电流和等效磁流,/>和b'p分别表示由高阶阻抗边界条件确定的各向异性介质系数;
对于p=1和p=2,分别执行以下步骤:
将等效磁流用等效电流表示,代入电场积分方程组中,执行步骤4-4到步骤6-3。
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 CN114491992A (zh) | 2022-05-13 |
CN114491992B true 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) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115184689B (zh) * | 2022-05-26 | 2024-03-19 | 西安电子科技大学 | 一种计算蒸发波导环境中目标电磁散射的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1404415A (zh) * | 2000-02-22 | 2003-03-19 | 基因谱公司 | 微阵列制造技术及设备 |
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 | 广州大学 | 一种反激式变换器电路板的电磁辐射仿真预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9125581B2 (en) * | 2012-05-23 | 2015-09-08 | Seiko Epson Corporation | Continuous modeling for dipole localization from 2D MCG images with unknown depth |
-
2022
- 2022-01-11 CN CN202210026307.3A patent/CN114491992B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1404415A (zh) * | 2000-02-22 | 2003-03-19 | 基因谱公司 | 微阵列制造技术及设备 |
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)
Title |
---|
An Efficient Hybrid Method of Iterative MoM-PO and Equivalent Dipole-Moment for Scattering From Electrically Large Objects;Li Xiao 等;《IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》;20170215;第16卷;第1723-1726页 * |
Efficient EDM-PO Method for the Scattering From Electrically Large Objects With the High-Order Impedance Boundary Condition;Hui Zhou 等;《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》;20220530;第70卷(第9期);第8242-8249页 * |
介质体电磁散射分析中的等效偶极矩法;郁超 等;《微波学报》;20100415;第26卷(第2期);第26-29页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114491992A (zh) | 2022-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106407723B (zh) | 面向低副瓣的稀疏排布阵列天线激励电流幅度的确定方法 | |
CN109765538B (zh) | 非均匀介质目标体的散射场确定方法 | |
CN114491992B (zh) | 一种基于等效偶极矩和物理光学法的高效电磁散射方法 | |
CN110737873B (zh) | 一种大规模阵列天线散射的快速分析方法 | |
CN107515956A (zh) | 一种基于hfss单元法的大型有限平面阵列分析方法 | |
CN112329204A (zh) | 考虑载体平台耦合的重复性结构电磁特性特征模快速分析方法 | |
CN110276086A (zh) | 基于空间映射的隐身超表面雷达散射截面减缩方法 | |
CN105630740B (zh) | 基于矩阵泰勒级数展开的电磁分析方法 | |
Zhao et al. | Parallel hybrid method of HOMoM–MLFMA for analysis of large antenna arrays on an electrically large platform | |
CN110287549B (zh) | 具有不确定源的多处薄涂敷飞行器rcs预测方法 | |
Zhou et al. | Efficient EDM-PO method for the scattering from electrically large objects with the high-order impedance boundary condition | |
CN113691294B (zh) | 一种近场稀疏阵列天线波束建立方法及装置 | |
CN113567941B (zh) | 一种目标近场多普勒回波快速仿真方法 | |
CN106156394B (zh) | 基于显式差分格式的电磁特性提取方法 | |
Chen et al. | Transformation optics-based finite difference time domain algorithm for scattering from object with thin dielectric coating | |
Tian et al. | Accelerated hybrid method for electromagnetic scattering from multiple complex targets above a rough surface | |
CN109657196B (zh) | 一种基于空域强耦合路径加速的mom-ka混合方法 | |
Guo et al. | Analysis of antenna array with arbitrarily shaped radomes using fast algorithm based on VSIE | |
Hou et al. | Acceleration of Trimmed MLFMA Using the PO Method for Large-Scale Scattering Problems | |
Wang et al. | Novel reduced matrix equation constructing method accelerates iterative solution of characteristic basis function method | |
CN111931353A (zh) | 一种应用于仿真fss结构的散射场求解方法 | |
Enriquez et al. | optIFnet: A Capacitive Antenna Dipole Indention-Flexure Predictive Model Optimized Using Hybrid Lichtenberg Algorithm and Neural Network | |
Gu et al. | Parallel Integral Equation Based Non-overlapping DDM for Fast Solving Electromagnetic Scattering Problems with Changeable Parts | |
CN114417557B (zh) | 超宽带rcs缩减超表面阵列的优化方法、系统及介质 | |
KR102380267B1 (ko) | 표적 탐지 방법 |
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 |