CN103226644A - 基于柱面等效源区域分解的电磁散射特性仿真方法 - Google Patents
基于柱面等效源区域分解的电磁散射特性仿真方法 Download PDFInfo
- Publication number
- CN103226644A CN103226644A CN2013101866524A CN201310186652A CN103226644A CN 103226644 A CN103226644 A CN 103226644A CN 2013101866524 A CN2013101866524 A CN 2013101866524A CN 201310186652 A CN201310186652 A CN 201310186652A CN 103226644 A CN103226644 A CN 103226644A
- Authority
- CN
- China
- Prior art keywords
- sub
- rwg
- basis function
- cylinder
- equivalent
- 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
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于柱面等效源区域分解的电磁散射特性仿真方法,步骤如下:建立模型,划分子区域并进行剖分;展开各个子等效柱面和子散射体上的电磁流,并确定每个子等效柱面的等效入射电磁流基函数系数;将非旋转对称体的子等效柱面等效入射电磁流基函数系数转换成RWG基函数系数;确定各个子等效柱面的等效入射电磁流在对应子散射体表面产生的入射电场,确定各个子区域内子散射体表面的散射电流;确定子等效柱面上的散射电磁流,将其作为其它各子区域的入射流,更新各子等效柱面上的散射电磁流系数直至达到平衡;确定远场雷达散射截面积,完成电磁散射特性仿真。本发明提供了一种稳定高效、适用于任意形状金属目标的电磁仿真方法。
Description
技术领域
本发明涉及电磁仿真技术领域,特别是一种基于柱面等效源区域分解的电磁散射特性仿真方法。
背景技术
旋转对称体是指绕一轴线旋转对称的三维目标,由于它的结构特性,可以由局部全域基函数展开,使用旋转对称体矩量法快速分析计算(M.Andreasen,"Scattering frombodies of revolution,"Antennas and Propagation,IEEE Transactions on,vol.13,pp.303-310,1965.),相比较传统的分域RWG基函数的矩量法(Rao,S.M.;Wilton,D.;Glisson,A.W.,Electromagnetic scattering by surfaces of arbitrary shape,"Antennas and Propagation,IEEE Transactions on,vol.30,no.3,pp.409,418,May1982)大大节约了计算资源。但实际目标中有很多目标结构不具有或者仅部分具有旋转对称体结构,不适合单独使用旋转对称体矩量法分析。Durham,T.E.,Christodoulou,C.G使用旋转对称基函数分析旋转对称部分,RWG基函数分析非旋转对称部分,旋转对称部分和非旋转对称部分的相互耦合使用的时两者基函数相互直接耦合的方法仿真(Durham,T.E.;Christodoulou,C.G.,"Electromagnetic radiation from structures consisting of combined body of revolution andarbitrary surfaces,"Antennas and Propagation,IEEE Transactions on,vol.40,no.9,pp.1061,1067,Sep1992),该方法的缺点是两种基函数直接相互作用计算时间和计算内存的消耗大严重的制约了该方法仿真问题的规模,对于完全不具有旋转对称体结构的目标仍然不能使用旋转对称体矩量法。
区域分解方法是近年来电磁仿真技术中热门的研究方向。W.-D.Li,W.Hong,andH.-X.Zhou提出了积分方程的重叠型区域分解方法(W.-D.Li,W.Hong,and H.-X.Zhou,“Integral equation-based overlapped domain decomposition method for the analysis ofelectromagnetic scattering of3D conducting objects,”Microwave and Optical TechnologyLetters,vol.49,no.2,pp.265-274,2007.),M.K.Li,and W.C.Chew提出了基于等效原理的区域分解方法(M.K.Li,and W.C.Chew,“Wave-Field Interaction With ComplexStructures Using Equivalence Principle Algorithm,”Antennas and Propagation,IEEETransactions on,vol.55,no.1,pp.130-138,2007.),P.Zhen,W.Xiao-Chuan,and L.Jin-Fa提出了非重叠型区域分解方法(P.Zhen,W.Xiao-Chuan,and L.Jin-Fa,“Integral EquationBased Domain Decomposition Method for Solving Electromagnetic Wave Scattering FromNon-Penetrable Objects,”Antennas and Propagation,IEEE Transactions on,vol.59,no.9,pp.3328-3338,2011.),这些方法都是建立在积分方程方法的基础上的,它们将原先的大问题分解成一个个小问题来求解,但是由于问题的维数较多,增加了计算资源的开销,且操作步骤比较复杂。
发明内容
本发明的目的在于提供一种高效稳定的基于柱面等效源区域分解的电磁散射特性仿真方法,能够对任意形状目标的电磁散射特性进行快速的仿真。
实现本发明的技术解决方案是:一种基于柱面等效源区域分解的电磁散射特性仿真方法,步骤如下:
第1步,建立散射体模型及等效柱面模型,划分子区域,对每个子区域内部的子散射体和等效柱面进行剖分;
第2步,根据第1步中子区域的剖分方式,定义相应的基函数,展开各个子等效柱面和子散射体上的电磁流,并确定每个子等效柱面的等效入射电磁流旋转对称体基函数系数;
第3步,对于非旋转对称体子区域的子等效柱面,将对应的等效入射电磁流旋转对称体基函数系数转换成RWG基函数系数;
第4步,确定各个子等效柱面的等效入射电磁流在对应子散射体表面产生的入射电场;
第5步,由第4步求得的入射电场,确定各个子区域内子散射体表面的散射电流;
第6步,由第5步得到的各个子区域内子散射体上的散射电流向外辐射,确定对应子等效柱面上的散射电磁流;
第7步,将第6步的各个子区域内子等效柱面上的散射电磁流,作为除自身之外其它各子区域的入射流,更新各个子区域内子等效柱面上的入射电磁流系数,并重新计算第3~7步,更新各子区域内子等效柱面上的散射电磁流系数,直至散射电磁流系数达到平衡,进入下一步;
第8步,由第7步所得等效柱面的表面散射电磁流,确定远场雷达散射截面积,完成电磁散射特性仿真。
本发明与现有技术相比其显著效果是:(1)利用了目标结构特性进行区域分解,很好的结合了传统的旋转对称体方法和多层快速多极子方法的优越性;(2)节约了计算资源的开销,对于各种类导弹结构的散射特性的分析本发明方法尤为适合;(3)整个方法理论可靠,易于实现,对于任意形状金属目标的电磁散射特性的快速仿真着重要的意义。
附图说明
图1为散射体柱面等效源区域分解示意图,其中(a)为任意散射体柱面等效源区域分解示意图,(b)为含部分旋转对称结构的散射体柱面等效源区域分解示意图。
图2为等效柱面剖分示意图,其中(a)为各子散射体块及等效柱面剖分示意图,(b)为非旋转对称部分三角形面片剖分示意图,(c)为旋转对称部分母线剖分示意图。
图3为RWG基函数示意图。
图4为实施例中导弹模型示意图,(a)为侧视图,(b)为主视图。
图5为实施例中导弹模型柱面等效源区域划分示意图。
图6为实施例中双站雷达散射截面积。
具体实施方式
下面结合附图,以一部分含有旋转对称体结构的散射体电磁散射特性的仿真为例,对本发明作进一步详细说明。
本发明基于柱面等效源区域分解的电磁散射特性仿真方法,包含以下步骤:
第1步,建立散射体模型及其等效柱面,划分子区域,对每个子区域内部的子散射体和子等效柱面进行剖分。
(1.1)建立一个可以完全包围待求散射体的封闭等效柱面,若散射体结构中部分为旋转对称体,那么等效柱面的轴向和散射体旋转对称部分的轴向一致;
(1.2)划分子区域:若模型不包含旋转对称体时,将等效柱面分解成尺寸完全相同的子等效柱面,每个子等效柱面及其包围的子散射体形成对应的子区域,如图1(a)所示;若模型是由旋转对称体部分和非旋转对称部分组合而成的散射体时,如图1(b)所示,则将等效圆柱面等分成两种子等效柱面:一种子等效柱面包围旋转对称体部分,为旋转对称体子区域;另一种子等效柱面包围非旋转对称部分,为非旋转对称体子区域。
(1.3)对每个子区域进行剖分,如图2(a)所示;对于非旋转对称体子区域的子等效柱面,分别进行三角形面片剖分和母线线段剖分,对于非旋转对称体子区域的子散射体进行三角形面片剖分,图2(b)为三角形面片剖分示意图;对于旋转对称体子区域的子等效柱面和子散射体的母线进行线段剖分,图2(c)为母线线段剖分示意图;两个子等效柱面相交的面使用一致的线段剖分,保证两个子区域电磁流的连续性,记录各个子区域的剖分信息。
第2步,根据第1步中子区域的剖分方式,定义相应的基函数,展开各个子等效柱面和子散射体上的电磁流,并确定每个子等效柱面的等效入射电磁流旋转对称体基函数系数。
(2.1)定义基函数:
其中n为RWG基函数编号,即为三角形剖分单元中第n条边,ln为第n条边的长度,表示第n条边对应的上三角形的面积,表示第n条边对应的下三角形的面积,表示第n条边对应的上三角形中r点指向自由顶点(即除去边n两个顶点剩下的三角形顶点)的方向向量,表示第n条边对应的下三角形中自由顶点指向r点的方向向量,如图3(a)所示。
则散射体表面的电磁流用RWG基函数展开为:
其中J(r)表示散射体表面任意一点r点的电流,M(r)表示r点的磁流,表示电流J(r)使用RWG基函数展开对应于第n个基函数的展开系数,表示磁流M(r)使用RWG基函数展开对应于第n个基函数的展开系数,Nrwg表示等效柱面使用RWG基函数展开后未知量的个数。
其中上标BoR表示旋转对称体,Tn'(t)表示三角基函数,是一维局部基函数,Tn'(t)定义在两条相连接的剖分线段上,我们分别称这两条线段为前段和后段,其表达式为Tn'(t)表示三角基函数,其表达式为:
其中n'为BoR三角基函数编号,t表示r点的母线方向分量;ρ(r)表示r点到旋转对称体旋转轴的垂直距离;φ表示r点的周向角;ejαφ表示傅里叶级数展开对应于第α个模式的指数项;表示r点的母线方向;表示r点的周向方向;NBoR表示旋转对称体基函数对应的未知量的个数;表示第n个三角基函数对应的前段的起点切向分量,表示第n个三角基函数对应的前段的终点切向分量即后段的起点切向分量,表示第n个三角基函数对应的后段的终点切向分量;Δn表示前段的长度,Δn+1表示后段的长度;
其中表示电流J(r)使用旋转对称体基函数展开对应于第α个模式的第n'个基函数母线方向的展开系数;表示电流J(r)使用旋转对称体基函数展开对应于第α个模式的第n'个基函数周向方向的展开系数;表示磁流M(r)使用旋转对称体基函数展开对应于第α个模式的第n'个基函数母线方向的展开系数;表示磁流M(r)使用旋转对称体基函数展开对应于第α个模式的第n'个基函数周向方向的展开系数。
(2.2)计算每个子等效柱面的等效入射电磁流旋转对称体基函数系数:
均匀平面波照射在子等效柱面:
(6)
Einc为入射电场,Hinc为入射磁场,为平面波传播方向单位矢量,k为传播常数,表示电场方向单位矢量,η为平面波阻抗。下面使用c标记子等效柱面,p标记子散射体,则柱面i上等效入射电流等效入射磁流分别为:
定义旋转对称体测试函数为基函数的共轭:
其中下标m'表示旋转对称体测试函数编号,β表示旋转对称体测试函数模式数编号,表示对应于切向方向第β个模式的第m'个测试函数,表示对应于周向方向的第β个模式第m'个测试函数,表示第i个子等效柱面上旋转对称体基函数对应的未知量的个数,Mod表示需要的总模式数;
使用公式(8)对公式(7)测试,生成子等效柱面的激励向量
Vc,i=[Vc,i,-Mod,…,,Vc,i,β…,Vc,i,Mod-1,Vc,i,Mod]T:
其中i=1,…,N,N表示子区域总数,β=-Mod,…,Mod,确定所有子区域的子等效柱面的等效入射电磁流。
确定第i个子区域的子等效柱面的等效入射电磁流旋转对称体基函数系数求解矩阵(UBoR i)-1,其中UBoR i的元素由下式给出:
其中Tm'q,i、Tn'q,i为三角基函数,由公式(4)给出,Tm'q,i表示第i个子区域对应于第q条线段的第m'个测试函数、Tn'q,i表示第i个子区域对应于第q条线段的第n'个三角基函数,Mq表示一个基函数包含的剖分线段数,这里取2,Δq,i表示第i个子区域对应的的第q条剖分线段的长度,由于三角基函数只和相邻的两个基函数有重叠,因此矩阵是一个三带条型矩阵,可以通过追赶法求解得到电流系数求解矩阵
则子区域i的子等效柱面模式β对应的等效入射电磁流系数向量为:
其中表示第i个子等效柱面c上的母线方向t对应的模式β等效入射电流系数向量;表示第i个子等效柱面c上的周向φ对应的模式β等效入射电流系数向量;表示第i个子等效柱面c上的母线方向t对应的模式β等效入射磁流系数向量;表示第i个子等效柱面c上的周向φ对应的模式β等效入射磁流系数向量。
第3步,对于非旋转对称体子区域的子等效柱面,将对应的等效入射电磁流旋转对称体基函数系数转换成RWG基函数系数;
(3.1)对于非旋转对称体子区域的子等效柱面,柱面等效入射电磁流使用RWG基函数表示,通过旋转对称体基函数系数可以确定第i个子等效柱面c上的任意一点r处的电流Jc,i(r)和磁流Mc,i(r):
其中表示旋转对称体基函数母线方向方向的转对称体基函数,表示旋转对称体基函数周向方向的转对称体基函数,表示第i个区域柱面c上的电流Jc,i(r)使用旋转对称体基函数展开对应于第β个模式的第n'个基函数母线方向的展开系数,表示表示第i个区域柱面c上的电流Jc,i(r)使用旋转对称体基函数展开对应于第β个模式的第n'个基函数周向方向的展开系数,表示第i个区域柱面c上的磁流Mc,i(r)使用旋转对称体基函数展开对应于第β个模式的第n'个基函数母线方向的展开系数,表示表示第i个区域柱面c上的磁流Mc,i(r)使用旋转对称体基函数展开对应于第β个模式的第n'个基函数周向方向的展开系数。
(3.2)由于无论使用哪种基函数表示,任意一点的电磁流的值是固定的,RWG基函数展开电磁流的表达式是:
其中Jc,i(r)、Mc,i(r)与公式(13)中含义相同,表示第i个子区域的子等效柱面c上电流RWG基函数展开系数,表示第i个子区域的子等效柱面c上磁流RWG基函数展开系数,表示第i个子等效柱面使用RWG基函数展开后未知量的个数。
(3.3)选择RWG基函数为测试函数对公式(14)两边进行测试,得到如下方程:
将公式(1)带入方程(15)并写出矩阵方程形式如下:
Ic,i为待求RWG基函数系数向量、为已知激励向量,(Urwg c,i)-1为RWG基函数系数求解矩阵,Urwg c,i具体元素为:
m表示RWG基函数测试函数编号,n表示RWG基函数编号,矩阵Uc,i,m,n为稀疏阵,每行每列包含5个非零元素,但对于未知量很大的情况,计算时间长和内存消耗大。为此我们利用RWG基函数是分域基函数,分别在每个三角形面片上建立方程(17),逐个三角形求解RWG基函数的系数,则Uc,i -1被划分成Ntriangle个3×3的子矩阵求逆的问题,Ntriangle为三角形个数。将这个求逆问题计算复杂度降低到了O(N)。则RWG基函数系数可以由求得。
同理,我们也可以通过同样的理论由RWG基函数的系数计算出旋转对称体基函数的系数,使用旋转对称体测试函数测试公式(13),并将公式(14)带入即能得到系数求解方程,其系数求解矩阵在第2步中已经给出,不再赘述。
第4步,确定各个子等效柱面的等效入射电磁流在对应子散射体表面产生的入射电场。
电流和磁流在空间中产生的电场E(r)和磁场H(r)可以由下式计算得到,其中k0为自由空间的传播常数:
其中
其中x代表电流J或磁流M,A(x)表示矢量位,ψ(x)表示标量位,C(x)表示旋度场,pv表示主值积分。g表示源点r'到场点r的格林函数,其表达式为:
定义微积分算子L、K:
K(x)=C(x)
则:
根据公式(22),我们可以计算出子等效柱面的等效入射电磁流在其内部子散射体上产生的电场:
EP,i表示第i个子区域的子等效柱面等效入射电流Jc,i和等效入射磁流Mc,i在对应内部子散射体P上产生的电场。
旋转对称体子区域的子散射体表面使用旋转对称体测试函数离散公式(23),得到由旋转对称体基函数表示的入射场EPBoR,i;非旋转对称体子区域的子散射体表面使用RWG测试函数确定入射场EPRWG,j。
第5步,由第4步求得的入射电场,确定各个子区域内子散射体表面的散射电流。
由电场积分公式:
Lpp,i(J)|t=-Ep,i|t (24)
Lpp,i(J)具体形式由公式(21)给出,下标pp表示积分算子L作用在第i个子区域内部的子散射体p上。
若内部为旋转对称体,则使用旋转对称体矩量法建立阻抗矩阵方程,母线剖分段数在500以内可以使用LU分解直接计算各个模式阻抗矩阵方程的逆,大于500则可以使用矩阵行列抽取技术加速迭代求解方程,求得旋转对称子散射体表面的散射电流系数散射电流系数由Mod个子模式电流组成β表示模式数;
若内部为非旋转对称散射体,则使用基于RWG基函数的矩量法求解,对每个子区域以柱面大小建立八叉树,利用多层快速多极子技术快速迭代求解子散射体表面的电流系数。多层快速多极子求解散射体电磁散射技术是一个较为成熟的技术,这里不再给出详细过程,求解到的非旋转对称子散射体表面的散射电流系数
第6步,由第5步得到的各个子区域内子散射体上的散射电流向外辐射,确定对应子等效柱面上的等效散射电磁流。
由公式(7)和公式(18)可以得到如下关系:
(1)确定具有旋转对称结构子散射体对应的子等效柱面上的散射电磁流:
使用旋转对称体测试函数对(25)两边测试:
则对应子等效柱面上的等效散射电磁流系数为:
其中系数求解矩阵由第二步中求得。表示第i个柱面c上的母线方向对应的模式β等效散射电流系数向量;表示第i个柱面c上的周向对应的模式β等效散射电流系数向量;表示第i个柱面c上的母线方向对应的模式β等效散射磁流系数向量;表示第i个柱面c上的周向对应的模式β等效散射磁流系数向量
(2)确定具有非旋转对称结构子散射体对应的子等效柱面上的散射电磁流:
使用RWG测试函数对(25)两边测试:
其中n±表示第j个区域的第n条边对应于上、下三角形各定义一个系数,分别表示上、下三角形面片上的电磁流系数。
按照三角形循环,使用公式(17)各个三角形求解表示电磁流系数:
其中表示三角形ntr对应的电磁流系数求解矩阵,表示子等效柱面上三角形ntr对应的电流、磁流散射向量,由(28)计算得到,表示子等效柱面上三角形ntr对应的等效散射电流系数,表示子等效柱面上三角形ntr对应的等效散射磁流系数。
接下来使用电磁流不同基函数系数转换算法将(29)式求得的RWG基函数电磁流系数转换为旋转对称体基函数系数,方法参考第3步。
第7步,将第6步的各个子区域内子等效柱面上的散射电磁流,作为除自身之外其它各子区域的入射流,更新各个子区域内子等效柱面上的入射电磁流系数,并重新计算第3~7步,更新各子区域内子等效柱面上的散射电磁流系数,直至散射电磁流系数达到平衡,进入下一步;
根据公式(18),各个子区域的子等效柱面上的散射电磁流会在除自身之外的其它子区域表面形成新的入射场:
其中 分别表示第6步中求得的第i个子区域柱面上的等效散射电、磁流, 分别表示第j个子区域更新的等效入射电、磁流, 叠加上第2步得到的外部均匀平面波产生的等效入射电磁流。重复第3~7步,更新各个子区域内子等效柱面上的散射电磁流系数,直至各个子区域内子等效柱面表面的散射电磁流达到平衡,不再跟新,进入下一步。当两个子区域相接时,要保证两边剖分线段一致且两边电磁流相同。
第8步,由第7步所得等效柱面的表面散射电磁流,确定远场雷达散射截面积,完成电磁散射特性仿真。
下面结合具体的实施例对本发明做进一步详细描述。
实施例1
根据本发明所述方法对模型号导弹模型使用本发明方法建模仿真其电磁散射特性。
某型号导弹模型,外形和尺寸如图4所示。弹径2m,弹长24m,弹身部分旋转对称体结构,方向舵部分为任意结构部分。工作频率为:300MHz,入射角度为弹头方向入射(00,00) ,垂直极化。建立一个半径为4.1m,高度为25m的柱面将弹体全部包围。以5.0m间隔横切柱面成5等分,如图5所示。每个子圆柱具有相同的结构,对其中一个建立母线的线段剖分和三角形基函数剖分,并将信息公用给其它柱面。内部子散射体根据内部结构选择线段剖分或三角形贴片剖分,区域1内子散射体使用RWG基函数计算,区域2到区域5使用旋转对称体基函数计算。按照本发明的计算步骤实施,最终计算得到双站雷达散射截面积,如图6中所示,结果很好的跟商用仿真软件FEKO吻合,证明本发明方法的正确性。
综上可知,本发明利用了目标结构特性进行区域分解,很好的结合了传统的旋转对称体方法和多层快速多极子方法的优越性;节约了计算资源的开销,对于各种类导弹结构的散射特性的分析本发明方法尤为适合;整个方法理论可靠,易于实现,对于任意形状金属目标的电磁散射特性的快速仿真着重要的意义。
Claims (3)
1.一种基于柱面等效源区域分解的电磁散射特性仿真方法,其特征在于,步骤如下:
第1步,建立散射体模型及等效柱面模型,划分子区域,对每个子区域内部的子散射体和等效柱面进行剖分;
第2步,根据第1步中子区域的剖分方式,定义相应的基函数,展开各个子等效柱面和子散射体上的电磁流,并确定每个子等效柱面的等效入射电磁流旋转对称体基函数系数;
第3步,对于非旋转对称体子区域的子等效柱面,将对应的等效入射电磁流旋转对称体基函数系数转换成RWG基函数系数;
第4步,确定各个子等效柱面的等效入射电磁流在对应子散射体表面产生的入射电场;
第5步,由第4步求得的入射电场,确定各个子区域内子散射体表面的散射电流;
第6步,由第5步得到的各个子区域内子散射体上的散射电流向外辐射,确定对应子等效柱面上的散射电磁流;
第7步,将第6步的各个子区域内子等效柱面上的散射电磁流,作为除自身之外其它各子区域的入射流,更新各个子区域内子等效柱面上的入射电磁流系数,并重新计算第3~7步,更新各子区域内子等效柱面上的散射电磁流系数,直至散射电磁流系数达到平衡,进入下一步;
第8步,由第7步所得等效柱面的表面散射电磁流,确定远场雷达散射截面积,完成电磁散射特性仿真。
2.根据权利要求1所述的基于柱面等效源区域分解的电磁散射特性仿真方法,其特征在于:第1步所述建立散射体模型及等效柱面模型,划分子区域,对每个子区域内部的子散射体和等效柱面进行剖分,具体为:
(1.1)建立一个可以完全包围待求散射体的封闭圆柱面,若散射体结构中部分为旋转对称体,那么圆柱面的轴向和旋转对称部分的轴向一致;
(1.2)划分子区域:若模型不包含旋转对称体时,将等效柱面分解成尺寸完全相同的子等效柱面,每个子等效柱面及其包围的子散射体形成对应的子区域;若模型是由旋转对称体部分和非旋转对称部分组合而成的散射体时,则将等效圆柱面等分成两种子等效柱面:一种子等效柱面包围旋转对称体部分,为旋转对称体子区域;另一种子等效柱面包围非旋转对称部分,为非旋转对称体子区域;
(1.3)对每个子区域进行剖分:对于非旋转对称体子区域的子等效柱面,分别进行三角形面片剖分和母线线段剖分,对于非旋转对称体子区域的子散射体进行三角形面片剖分;对于旋转对称体子区域的子等效柱面和子散射体的母线进行线段剖分。
3.根据权利要求1所述的基于柱面等效源区域分解的电磁散射特性仿真方法,其特征在于:第3步所述对于非旋转对称体子区域的子等效柱面,将对应的等效入射电磁流旋转对称体基函数系数转换成RWG基函数系数,具体过程如下:
(3.1)对于非旋转对称体子区域的子等效柱面,柱面等效入射电磁流使用RWG基函数表示,通过旋转对称体基函数系数可以确定第i个子等效柱面c上的任意一点r处的电流Jc,i(r)和磁流Mc,i(r):
其中表示旋转对称体基函数母线方向方向的转对称体基函数,表示旋转对称体基函数周向方向的转对称体基函数,表示第i个区域柱面c上的电流Jc,i(r)使用旋转对称体基函数展开对应于第β个模式的第n'个基函数母线方向的展开系数,表示表示第i个区域柱面c上的电流Jc,i(r)使用旋转对称体基函数展开对应于第β个模式的第n'个基函数周向方向的展开系数,表示第i个区域柱面c上的磁流Mc,i(r)使用旋转对称体基函数展开对应于第β个模式的第n'个基函数母线方向的展开系数,表示表示第i个区域柱面c上的磁流Mc,i(r)使用旋转对称体基函数展开对应于第β个模式的第n'个基函数周向方向的展开系数;Ni BoR表示第i个子等效柱面上旋转对称体基函数对应的未知量的个数,n'为BoR三角基函数编号,β表示模式数编号;
(3.2)RWG基函数展开电磁流的表达式是:
其中Jc,i(r)、Mc,i(r)与公式(1)中含义相同,表示第i个子区域的子等效柱面c上电流RWG基函数展开系数,表示第i个子区域的子等效柱面c上磁流RWG基函数展开系数;表示第i个子等效柱面使用RWG基函数展开后未知量的个数,n为RWG基函数编号;
(3.3)选择RWG基函数为测试函数对公式(2)两边进行测试,得到如下方程:
将RWG基函数公式带入方程(3)并写出矩阵方程形式如下:
m表示RWG基函数测试函数编号,n表示RWG基函数编号,i=1,…,N,N表示子区域总数,r'为源点,r为场点;矩阵Uc,i,m,n为稀疏阵,分别在每个三角形面片上建立方程如式(5),逐个三角形求解RWG基函数的系数,则Uc,i -1被划分成Ntriangle个子矩阵求逆的问题,Ntriangle为三角形个数,则RWG基函数系数可以由得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310186652.4A CN103226644B (zh) | 2013-05-17 | 2013-05-17 | 基于柱面等效源区域分解的电磁散射特性仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310186652.4A CN103226644B (zh) | 2013-05-17 | 2013-05-17 | 基于柱面等效源区域分解的电磁散射特性仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103226644A true CN103226644A (zh) | 2013-07-31 |
CN103226644B CN103226644B (zh) | 2016-01-20 |
Family
ID=48837089
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310186652.4A Active CN103226644B (zh) | 2013-05-17 | 2013-05-17 | 基于柱面等效源区域分解的电磁散射特性仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103226644B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104133376A (zh) * | 2014-06-25 | 2014-11-05 | 电子科技大学 | 一种用于射频仿真中目标的矢量控制方法 |
CN104699870A (zh) * | 2013-12-06 | 2015-06-10 | 南京理工大学 | 电大复杂有耗介质目标电磁散射抛物线快速仿真方法 |
CN104699879A (zh) * | 2013-12-10 | 2015-06-10 | 南京理工大学 | 复杂多目标电磁散射的多次旋转等效仿真方法 |
CN104778286A (zh) * | 2014-01-09 | 2015-07-15 | 南京理工大学 | 掠海飞行器电磁散射特性快速仿真方法 |
CN105955924A (zh) * | 2016-04-20 | 2016-09-21 | 南京航空航天大学 | 一种快速求解局部变化目标的电磁散射特性的方法 |
CN106202594A (zh) * | 2015-05-07 | 2016-12-07 | 南京理工大学 | 分析混合目标瞬态电磁散射特性的时域不连续伽辽金方法 |
TWI574001B (zh) * | 2013-11-26 | 2017-03-11 | Asml荷蘭公司 | 用於計算結構之電磁散射屬性及估計其幾何與材料參數之方法與裝置 |
CN106503349A (zh) * | 2016-10-25 | 2017-03-15 | 中国人民解放军国防科学技术大学 | 一种类周期结构目标电磁散射特性快速计算方法 |
CN108445303A (zh) * | 2018-03-08 | 2018-08-24 | 北京环境特性研究所 | 一种近场电磁散射特性模拟方法 |
CN109063386A (zh) * | 2018-09-26 | 2018-12-21 | 中国人民解放军陆军工程大学 | 一种密集旋转对称体群电磁散射仿真方法 |
CN113673098A (zh) * | 2021-08-12 | 2021-11-19 | 广州广电计量检测股份有限公司 | 一种无线设备的电磁辐射比吸收率仿真检测方法及装置 |
CN114065559A (zh) * | 2022-01-18 | 2022-02-18 | 常州星宇车灯股份有限公司 | 一种基板集成波导天线的设计方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508220A (zh) * | 2011-10-24 | 2012-06-20 | 西瑞克斯(北京)通信设备有限公司 | 均匀双各向同性媒质物体的雷达散射截面获取方法 |
-
2013
- 2013-05-17 CN CN201310186652.4A patent/CN103226644B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508220A (zh) * | 2011-10-24 | 2012-06-20 | 西瑞克斯(北京)通信设备有限公司 | 均匀双各向同性媒质物体的雷达散射截面获取方法 |
Non-Patent Citations (4)
Title |
---|
WANG QUANQUAN等: "analysis of transient electromagnetic scattering using uv method enhanced time-domain integral equations with laguerre polynomials", 《MICROWAVE AND OPTICAL TECHNOLOGY LETTERS》, vol. 53, no. 1, 31 January 2011 (2011-01-31), pages 158 - 163 * |
周世岗: "目标电磁散射特征建模", 《中国优秀硕士学位论文全文数据库 基础科学辑》, no. 11, 15 November 2009 (2009-11-15) * |
王湃等: "多层快速多极子分析三维复杂目标的谐振区电磁散射特性", 《北京大学学报(自然科学版)》, vol. 42, no. 3, 31 May 2006 (2006-05-31), pages 395 - 400 * |
阳晓丽: "旋转对称体散射特性的高效分析", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 11, 15 November 2008 (2008-11-15) * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9939250B2 (en) | 2013-11-26 | 2018-04-10 | Asml Netherlands B.V. | Methods and apparatus for calculating electromagnetic scattering properties of a structure and for estimation of geometrical and material parameters thereof |
TWI574001B (zh) * | 2013-11-26 | 2017-03-11 | Asml荷蘭公司 | 用於計算結構之電磁散射屬性及估計其幾何與材料參數之方法與裝置 |
CN104699870A (zh) * | 2013-12-06 | 2015-06-10 | 南京理工大学 | 电大复杂有耗介质目标电磁散射抛物线快速仿真方法 |
CN104699870B (zh) * | 2013-12-06 | 2017-12-12 | 南京理工大学 | 电大复杂有耗介质目标电磁散射抛物线快速仿真方法 |
CN104699879A (zh) * | 2013-12-10 | 2015-06-10 | 南京理工大学 | 复杂多目标电磁散射的多次旋转等效仿真方法 |
CN104699879B (zh) * | 2013-12-10 | 2018-04-03 | 南京理工大学 | 复杂多目标电磁散射的多次旋转等效仿真方法 |
CN104778286A (zh) * | 2014-01-09 | 2015-07-15 | 南京理工大学 | 掠海飞行器电磁散射特性快速仿真方法 |
CN104778286B (zh) * | 2014-01-09 | 2019-01-15 | 南京理工大学 | 掠海飞行器电磁散射特性快速仿真方法 |
CN104133376A (zh) * | 2014-06-25 | 2014-11-05 | 电子科技大学 | 一种用于射频仿真中目标的矢量控制方法 |
CN106202594A (zh) * | 2015-05-07 | 2016-12-07 | 南京理工大学 | 分析混合目标瞬态电磁散射特性的时域不连续伽辽金方法 |
CN105955924A (zh) * | 2016-04-20 | 2016-09-21 | 南京航空航天大学 | 一种快速求解局部变化目标的电磁散射特性的方法 |
CN106503349A (zh) * | 2016-10-25 | 2017-03-15 | 中国人民解放军国防科学技术大学 | 一种类周期结构目标电磁散射特性快速计算方法 |
CN106503349B (zh) * | 2016-10-25 | 2019-08-13 | 中国人民解放军国防科学技术大学 | 一种类周期结构目标电磁散射特性快速计算方法 |
CN108445303A (zh) * | 2018-03-08 | 2018-08-24 | 北京环境特性研究所 | 一种近场电磁散射特性模拟方法 |
CN109063386A (zh) * | 2018-09-26 | 2018-12-21 | 中国人民解放军陆军工程大学 | 一种密集旋转对称体群电磁散射仿真方法 |
CN109063386B (zh) * | 2018-09-26 | 2023-01-03 | 中国人民解放军陆军工程大学 | 一种密集旋转对称体群电磁散射仿真方法 |
CN113673098A (zh) * | 2021-08-12 | 2021-11-19 | 广州广电计量检测股份有限公司 | 一种无线设备的电磁辐射比吸收率仿真检测方法及装置 |
CN114065559A (zh) * | 2022-01-18 | 2022-02-18 | 常州星宇车灯股份有限公司 | 一种基板集成波导天线的设计方法 |
CN114065559B (zh) * | 2022-01-18 | 2022-04-15 | 常州星宇车灯股份有限公司 | 一种基板集成波导天线的设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103226644B (zh) | 2016-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103226644A (zh) | 基于柱面等效源区域分解的电磁散射特性仿真方法 | |
CN103323845B (zh) | 一种非均匀采样综合孔径辐射计的图像反演方法 | |
CN103218487B (zh) | 旋转对称天线罩和抛物面天线一体化电磁散射仿真方法 | |
CN101944145B (zh) | 可去除伪直流模式的微波管高频电路有限元仿真方法 | |
CN103279589B (zh) | 基于矩阵嵌套压缩的旋转对称体电磁散射特性仿真方法 | |
CN103870654A (zh) | 基于并行矩量法与物理光学混合的电磁散射仿真方法 | |
CN103425816B (zh) | 快速获取金属旋转对称体电磁散射特性的矩阵抽取方法 | |
CN102087359A (zh) | 一维镜像综合孔径辐射成像方法 | |
CN104112051B (zh) | 基于单元形函数的反射面天线机电集成设计方法 | |
CN103279601A (zh) | 导体目标宽带电磁散射特性的仿真方法 | |
CN107992676A (zh) | 一种运动目标时域散射回波的快速仿真建模方法 | |
CN104535998A (zh) | 旋转镜像综合孔径辐射计及测量方法 | |
CN104750944A (zh) | 一种快速的平面稀疏阵列综合方法 | |
CN104573376A (zh) | 一种时域有限差分法计算电磁散射的瞬态场远场外推方法 | |
CN102054094A (zh) | 一种平面微带电路的快速多层方向性仿真方法 | |
Sherman et al. | Structure factors of the kagome-lattice Heisenberg antiferromagnets at finite temperatures | |
CN103412988A (zh) | 基于相移降阶模型周期结构的三维电磁场仿真模拟方法 | |
CN104699879B (zh) | 复杂多目标电磁散射的多次旋转等效仿真方法 | |
CN104915326A (zh) | 基于等效原理的区域分解阶数步进时域积分方法 | |
Yang et al. | FDTD/TDPO hybrid approach for analysis of the EM scattering of combinative objects | |
CN105277927A (zh) | 飞行器编队瞬态电磁特性时域阶数步进分析方法 | |
CN106156475A (zh) | 电大尺寸目标的瞬态电磁特性快速提取方法 | |
CN106021766A (zh) | 基于机电集成优化的圆极化网状天线高指向精度设计方法 | |
CN103983948A (zh) | 基于稀疏表示的目标角度距离联合估计方法 | |
CN104346488B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |