CN103218487B - 旋转对称天线罩和抛物面天线一体化电磁散射仿真方法 - Google Patents
旋转对称天线罩和抛物面天线一体化电磁散射仿真方法 Download PDFInfo
- Publication number
- CN103218487B CN103218487B CN201310120271.6A CN201310120271A CN103218487B CN 103218487 B CN103218487 B CN 103218487B CN 201310120271 A CN201310120271 A CN 201310120271A CN 103218487 B CN103218487 B CN 103218487B
- Authority
- CN
- China
- Prior art keywords
- alpha
- represent
- rotational
- array antenna
- phi
- 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
Landscapes
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Aerials With Secondary Devices (AREA)
Abstract
本发明公开了一种旋转对称天线罩和抛物面天线一体化电磁散射仿真方法,步骤如下:(1)建立旋转对称天线罩局部直角坐标系、抛物面天线局部直角坐标系、以及虚拟等效球面,然后对各条母线剖分离散;(2)加入激励源,确定激励向量;(3)分别建立两个直角坐标系下的阻抗矩阵及其逆矩阵;(4)确定散射矩阵和传输矩阵;(5)通过坐标系旋转建立两个直角坐标系之间的关系,将虚拟等效球面上的等效电磁流在两个局部直角坐标系之间相互转换;(6)根据第(2)~(5)步的信息建立求解方程组,解算得到旋转对称天线罩和虚拟等效球面上的等效散射电磁流;(7)由互易定理确定雷达散射截面积。本发明提供了一种速度快、内存消耗低、求解精度高的旋转对称天线罩和抛物面天线一体化电磁散射仿真方法。
Description
技术领域
本发明涉及电磁仿真技术,特别是一种旋转对称天线罩和抛物面天线一体化电磁散射仿真方法。
背景技术
利用矩量法分析和研究各类旋转对称体的电磁散射问题一直是电磁散射仿真研究课题中的一个热点。但是随着天线罩尺寸的增大,传统的基于RWG基函数剖分建模的电磁散射仿真方法需要消耗大量的计算机内存和花费很长计算时间,个人计算机已经无法胜任。一般天线罩多为介质材料的旋转对称结构,适用旋转对称体矩量法来进行电磁散射仿真。
Andreasen,M.G在1965年首先提出了旋转对称体矩量法。文中将入射平面波利用傅立叶级数展开为相互正交的柱面波形式,利用各模式间的正交性,分别求解单一模式下的感应电流,然后进行线性叠加,从而求得表面感应电流的分布(M.Andreasen,"Scattering from bodies of revolution,"Antennas and Propagation,IEEE Transactions on,vol.13,pp.303-310,1965.)。1977年T.K.Wu和L.L.Tsai(T.K.Wu and L.L.Tsai,Scatteringfrom arbitrarily-shaped lossy dielectric bodies of revolution,Radio Science,Sep-Oct.1977,vol.12,NO.5)利用Maxwell方程组、等效原理和边界条件建立了介质旋转对称体电磁散射的数学模型,得到了四个积分方程。由于未知量只有两个,所以这是一个超定问题。利用Poggio,Miller,Chang,Harring,Wu的思想将四个积分方程化解为两个等价的积分方程,即PMCHW(以五位学者名字的首字母命名)积分方程,使问题化解为了一个限定问题,并对其应用矩量法求解,将等效电流和等效磁流展开成互相正交的傅立叶级数,再结合伽略金(Galerkin)法,最终解决了旋转介质体的散射问题。但基于PMCHW方程的旋转对称体矩量法方法仅适用于仅有一个旋转轴的旋转对称体电磁散射特性的电磁散射仿真计算。虽然天线罩和抛物面天线都具有旋转对称结构,介质部分可以使用PMCHW方程,金属部分可以使用电场积分方程分别分析(张剑锋,"旋转对称介质壳体的矩量法分析,"2004),但是二者一体化分析时不具有旋转对称特性,旋转对称体矩量法方法不能得到相互正交的独立模式,方法失效。
发明内容
本发明的目的在于提供一种速度快、内存消耗低、求解精度高的旋转对称天线罩和抛物面天线一体化电磁散射仿真方法。
实现本发明的技术解决方案是:一种旋转对称天线罩和抛物面天线一体化电磁散射仿真方法,包括以下步骤:
第一步,建立模型和局部直角坐标系:以旋转对称天线罩对称轴和抛物面天线对称轴的交点为坐标系原点O,在两条对称轴所在的平面分别建立旋转对称天线罩局部直角坐标系XOZ和抛物面天线局部直角坐标系X’OZ’;以该坐标系原点O为球心,建立一个完全包围抛物面天线的虚拟等效球面;在旋转对称天线罩局部直角坐标系XOZ下建立旋转对称天线罩的母线和虚拟等效球面的第一母线,在抛物面天线局部直角坐标系X’OZ’下建立抛物面天线的母线和虚拟等效球面的第二母线;然后对各条母线分别进行剖分离散;
第二步,加入激励源,确定激励向量:根据等效原理确定旋转对称天线罩上的等效入射电磁流和虚拟等效球面上的等效入射电磁流;
第三步,分别建立旋转对称天线罩局部直角坐标系XOZ和抛物面天线局部直角坐标系X’OZ’下的阻抗矩阵,并求解各自阻抗矩阵的逆矩阵;
第四步,确定散射矩阵和传输矩阵:散射矩阵是描述抛物面天线局部直角坐标系下,虚拟等效球面上等效入射电磁流和等效散射电磁流之间的关系;传输矩阵是描述旋转对称天线罩局部直角坐标系下,虚拟等效球面上的感应电流在旋转对称天线罩上产生的二次散射场和旋转对称天线罩上的感应电磁流在虚拟等效球面上产生的二次散射场;
第五步,通过坐标系旋转建立旋转对称天线罩局部直角坐标系XOZ和抛物面天线局部直角坐标系X’OZ’之间的关系,将虚拟等效球面上的等效电磁流在两个局部直角坐标系之间相互转换;
第六步,根据第二步~第五步的信息建立求解方程组,解算得到旋转对称天线罩和虚拟等效球面上的等效散射电磁流;
第七步,由互易定理确定雷达散射截面积。
本发明与现有技术相比其显著效果是:(1)使用了旋转对称体矩量法法、坐标系旋转电磁流转换及基于零场等效原理的区域分解思想,成功仿真了两个不共旋转轴的旋转对称体散射、辐射问题,克服了传统旋转对称体矩量法不能对旋转对称天线罩和抛物面天线一体化高效电磁散射仿真的缺陷;(2)相比传统的矩量法,其很好的利用了模型本来的结构特点,使用傅里叶级数将待求问题降维,降低了电磁仿真对计算机资源的消耗,计算速度快、内存少、精度高;(3)整个方法简单而易于实现,对旋转对称天线罩和抛物面天线的仿真与设计有重要意义。
附图说明
图1是本发明旋转对称天线罩和抛物面天线的结构示意图。
图2是本发明旋转对称天线罩、抛物面天线和虚拟等效球面母线建模示意图。
图3是本发明虚拟等效球面上母线位置随坐标轴旋转示意图。
图4是本发明母线剖分示意图和三角基函数示意图。
图5是本发明散射矩阵描述图。
图6是本发明坐标系旋转示意图。
图7是本发明实施例中旋转对称天线罩和抛物面天线结构示意图。
图8为本发明实施例中双站雷达散射截面积。
具体实施方式
下面结合附图对本发明做进一步详细描述。
以一个旋转对称天线罩和抛物面天线的电磁散射仿真为例(如图1所示),对实现本发明的具体步骤作进一步阐述:
第一步,结合图1,建立模型和局部直角坐标系:以旋转对称天线罩对称轴和抛物面天线对称轴的交点为坐标系原点O,在两条对称轴所在的平面分别建立旋转对称天线罩局部直角坐标系XOZ和抛物面天线局部直角坐标系X’OZ’;以该坐标系原点O为球心,建立一个完全包围抛物面天线的虚拟等效球面;结合图2(a)~图2(c)在旋转对称天线罩局部直角坐标系XOZ下沿Z轴方向建立旋转对称天线罩的母线和虚拟等效球面的第一母线,在抛物面天线局部直角坐标系X’OZ’下沿Z’轴方向建立抛物面天线的母线和虚拟等效球面的第二母线;然后对各条母线分别进行剖分离散;所述旋转对称天线罩局部直角坐标系XOZ以旋转对称天线罩的旋转轴为Z轴,抛物面天线局部直角坐标系X’OZ’以抛物面天线的旋转轴为Z’轴,旋转对称天线罩局部直角坐标系XOZ与抛物面天线局部直角坐标系X’OZ’共平面;旋转对称天线罩完全在虚拟等效球面之外,如图3所示;所述对各条母线分别进行剖分离散,具体方法如下:
使用三角基函数描述切向电磁流,使用傅里叶级数展开的轴向基函数描述周向电磁流,目标表面的电磁流可以展开为:
其中J(r)表示旋转对称体表面任意一点r点的电流,M(r)表示r点的磁流,表示电流J(r)使用旋转对称体基函数展开对应于第α个模式的第n个基函数切向方向的展开系数;表示电流J(r)使用旋转对称体基函数展开对应于第α个模式的第n个基函数周向方向的展开系数;表示磁流M(r)使用旋转对称体基函数展开对应于第α个模式的第n个基函数切向方向的展开系数;表示磁流M(r)使用旋转对称体基函数展开对应于第α个模式的第n个基函数周向方向的展开系数;α表示旋转对称体基函数的模式数;旋转对称体基函数表达式为:
其中Tn(t)表示三角基函数,是一维局部基函数,Tn(t)定义在两条相连接的剖分线段上,我们分别称这两条线段为前段和后段,其表达式为:
t表示r点的切向分量;ρ表示r点到z轴的垂直距离;φ表示r点的周向角;ejαφ表示表示傅里叶级数展开对应于第α个模式的指数项;N表示旋转对称体母线剖分段数,即旋转对称体基函数对应的未知量的个数;表示第n个三角基函数对应的前段的起点切向分量,表示第n个三角基函数对应的前段的终点切向分量即后段的起点切向分量,表示第n个三角基函数对应的后段的终点切向分量;Δn表示前段的长度,Δn+1表示后段的长度;
母线的剖分与三角基函数示意图如图4所示,抛物面天线的母线按照十分之一波长1/(10λ)剖分,旋转对称天线罩的母线按照十分之一介质波长剖分,虚拟等效球面的母线按照七分之一波长1/(7λ)剖分,εr表示介质的相对介电常数,其中λ表示入射电磁波在空气中传播时的波长。
第二步,加入激励源(以平面波入射为例),确定激励向量:根据等效原理确定旋转对称天线罩上的等效入射电磁流和虚拟等效球面上的等效入射电磁流;所述的激励向量包括旋转对称天线罩的激励向量VD和虚拟等效球面的激励向量VH,以平面波入射为例(若电磁散射仿真辐射问题,以偶极子天线的辐射方向图为源入射,这里不做赘述)确定激励向量,具体确定方法如下:
(1)确定旋转对称天线罩的激励向量VD
在旋转对称天线罩表面建立PMCHW方程,激励向量为旋转对称天线罩表面的等效入射电磁流,具体公式如下:
其中下标D表示旋转对称天线罩,上标inc表示入射场,表示旋转对称天线罩表面的入射电场,表示旋转对称天线罩表面的入射磁场,表示旋转对称天线罩表面的等效入射电流,表示旋转对称天线罩表面的等效入射磁流,表示旋转对称天线罩的法向方向;
使用旋转对称基函数的共轭做伽辽金测试:
(5)
其中下标m表示基函数编号,α表示模式数编号,表示对应于第α个模式的第m个基函数切向方向的测试基函数,表示对应于第α个模式的第m个基函数周向方向的测试基函数,N表示旋转对称体母线剖分线段数,Mod表示旋转对称体需要的总模式数;
对公式(4)两边使用公式(5)进行测试得到旋转对称天线罩D的激励向量VD:
VJ D表示旋转对称天线罩D表面上等效入射电流对应的激励向量,VM D表示旋转对称天线罩D表面上等效入射磁流对应的激励向量,其每个元素为:
其中<·>表示内积,ND表示旋转对称天线罩母线剖分线段数,ModD表示旋转对称体天线罩的总模式数;
(2)确定虚拟等效球面的激励向量VH
与公式(1)中确定旋转对称天线罩的激励向量VD的方法相同,以下标H表示虚拟等效球面,那么
表示虚拟等效球面表面的入射电场,表示虚拟等效球面表面的入射磁场,表示虚拟等效球面表面的等效入射电流,表示虚拟等效球面表面的入射磁流,表示虚拟等效球面的法向方向;
使用公式(5)所示同样形式的测试函数对公式(8)两边进行测试,得到虚拟等效面的激励向量VH:
VJ H表示虚拟等效球面H表面上等效J对应的激励向量,VM H表示虚拟等效球面H表面上等效磁流M对应的激励向量,其每个元素为:
其中下标l表示基函数编号,β表示模式数,表示对应于第β个模式的第l个基函数切向方向的测试基函数,表示对应于第β个模式的第l个基函数周向方向的测试基函数,NH表示虚拟等效球面母线剖分线段数,ModH表示虚拟等效球面上入射等效电磁流使用基函数展开需要的最大模式数。
由于基函数使用了傅里叶级数展开,所以需要对对应的模式数进行截断。由于两个局部直角坐标系下的入射波的入射角度不同,以两个模式数截断标准中max(ModD,ModH)做为总的截断标准。
第三步,分别建立旋转对称天线罩局部直角坐标系XOZ和抛物面天线局部直角坐标系X’OZ’下的阻抗矩阵,并求解各自阻抗矩阵的逆矩阵,具体建立过程如下:
(1)建立旋转对称天线罩的阻抗矩阵
在旋转对称天线罩上建立PMCHW方程各个模式阻抗矩阵
表示第α个模式的第m个基函数切向方向对应的电流,表示第α个模式的第m个基函数切向方向对应的磁流,表示第α个模式的第m个基函数周向方向对应的电流,表示第α个模式的第m个基函数周向方向对应的磁流;η1表示空气中的波阻抗,η2表示介质的波阻抗;L、K算子表示积分算子,对应的第一个下标表示场区所在位置,第二个下标表示源区所在位置;而“1”则表示旋转对称天线罩的外表面即与空气接触的表面,“2”表示内表面即指向内部介质的表面,即L1,1、K11则表示外表面上的电磁流在自身产生的场,L2,1,K21表示外表面电磁流在内表面产生的场;算子L、K具体积分表达形式如下:
χ(r')表示等效电流J(r')或等效磁流M(r'),它们将使用公式(1)离散,与各个模式电磁流分量相对应,ω表示角频率,ε表示电导率,μ表示磁导率;G(r,r')表示源点r'到场点r的格林函数其表达式为:
逐个确定出各模式数对应的逆矩阵,使用LU分解算法,表示为
(2)建立抛物面天线的模式阻抗矩阵
在抛物面天线上建立电场积分方程各个模式阻抗矩阵
其中上标P表示抛物面天线,Wβ表示β模式对应的测试函数,ModH表示抛物面天线的总模式数,L(Jβ)表示入射模式电流Jβ在抛物面天线P上产生的激励电场,算子L如公式(12)所示;
逐个确定出各模式数对应的逆矩阵,使用LU分解算法,表示为
第四步,确定散射矩阵和传输矩阵:散射矩阵是描述抛物面天线局部直角坐标系下,虚拟等效球面上等效入射电磁流和等效散射电磁流之间的关系;传输矩阵是描述旋转对称天线罩局部直角坐标系下,虚拟等效球面上的感应电流在旋转对称天线罩上产生的二次散射场和旋转对称天线罩上的感应电磁流在虚拟等效球面上产生的二次散射场;
(1)确定散射矩阵S
在抛物面天线局部直角坐标系x'oz'下,抛物面天线和虚拟等效球面的第二母线是共轴的,如图3中x’oz’坐标平面所示,使用旋转对称矩量法思想确定出各个模式独立的散射矩阵S,散射矩阵S描述的是虚拟等效球面上的等效入射电磁流经过抛物面天线作用后在虚拟等效球面上得到的等效散射电磁流之间的关系:
表示虚拟等效球面表面的等效入射电流,表示虚拟等效球面表面的入射磁流,表示虚拟等效球面表面的等效散射电流,表示虚拟等效球面表面的散射磁流;
图5给出了散射矩阵S描述的物理意义,散射矩阵S由三部分构成:虚拟等效球面上的模式β对应的等效入射电磁流对抛物面天线的作用描述为Zβ PH;抛物面天线自身的作用,描述为第三步中求得的抛物面天线阻抗的逆由抛物面天线自身产生的模式β对应的散射场在虚拟等效球面上得到的模式β对应的等效散射电磁流的作用描述为Zβ HP,则S矩阵可以写成:
其中U-1为电流系数求解矩阵:
其中fαn表示旋转对称体基函数,n表示基函数编号α表示模式数,Wβm表示旋转对称体测试函数,m表示测试函数编号,β表示模式数。代入具体得到:
其中Tmp、Tnp表示对应于m,n第p调线段的三角基函数,Mp表示一个基函数包含的剖分线段数,这里取2。Δp表示对应的剖分线段的长度。由于三角基函数只和相邻的两个基函数有重叠,因此矩阵U是一个三带条型矩阵,可以通过追赶法求解得到电流系数求解矩阵U-1。
L,K积分算子由公式(12)给出,下标H表示虚拟等效球面,P表示抛物面天线表面,即积分算子LPH(Jβ,H)表示虚拟等效球面H外部等效面上的第β模式的电流Jβ,H在抛物面天线P表面产生的β模式的激励电场,积分算子KPH(Mβ,H)表示虚拟等效球面H外部等效面上的第β模式的磁流Mβ,H在抛物面天线P表面产生的β模式的激励磁场;积分算子KHP(Jβ,P)表示抛物面天线P上的β模式的电流Jβ,P在外虚拟等效球面H外部等效面上产生的β模式的激励磁场,积分算子LHP(Jβ,P)表示抛物面天线P上的β模式的电流Jβ,P在外虚拟等效球面H外部等效面上产生的β模式的激励电场;
(2)确定传输矩阵TDH和THD
在旋转对称天线罩局部直角坐标系XOZ下,旋转对称天线罩母线和虚拟等效球面的第二母线是共轴的,使用旋转对称矩量法思想确定出各个模式独立的传输矩阵TDH和THD,TDH表示虚拟等效球面上的感应电流在旋转对称天线罩上产生的二次散射场,THD表示旋转对称天线罩上的感应电磁流在虚拟等效面上产生的二次散射场:
L,K积分算子的积分表达式由公式(12)给出,下标H表示虚拟等效球面,D表示天线罩表面,即积分算子KHD(JD)表示旋转对称天线罩D上的电流J在虚拟等效球面H上的激励磁场,积分算子LHD(MD)表示旋转对称天线罩D上的磁流M在虚拟等效球面H上的激励磁场;积分算子LHD(JD)表示旋转对称天线罩D上的电流J在虚拟等效球面H上的激励电场,积分算子KHD(MD)表示旋转对称天线罩D上的磁流M在虚拟等效球面H上的激励电场;积分算子KDH(JH)表示虚拟等效球面H上的电流JH在旋转对称天线罩D上的激励磁场,积分算子LDH(MH)表示虚拟等效球面H上的电流MH在旋转对称天线罩D上的激励磁场;积分算子LDH(JH)表示虚拟等效球面H上的电流JH在旋转对称天线罩D上的激励电场,积分算子KDH(MH)表示虚拟等效球面H上的电流MH在旋转对称天线罩D上的激励电场。公式(20)和(21)也使用伽辽金测试生成传输矩阵的,其过程和公式(11)完全类似,不在赘述。
第五步,通过坐标系旋转建立旋转对称天线罩局部直角坐标系XOZ和抛物面天线局部直角坐标系X’OZ’之间的关系,将虚拟等效球面上的等效电磁流在两个局部直角坐标系之间相互转换;等效球面上每一点上等效电磁流是唯一的,但在两个局部直角坐标系下表示不同,我们需要通过电磁流的唯一性建立两个坐标系下电磁流表示的转换关系,具体为:
在旋转对称天线罩局部直角坐标系下虚拟等效球面上的电磁流表示记为 在抛物面天线局部直角坐标系下虚拟等效球面上的电磁流表示记为 建立如下关系:
其中ZRP和ZPR表示坐标系旋转矩阵,ZRP表示将在抛物面天线局部直角坐标系下的虚拟等效球面上的等效电磁流转换到旋转对称天线罩局部直角坐标系下虚拟等效球面上的等效电磁流的旋转矩阵,ZPR表示将旋转对称天线罩局部直角坐标系下虚拟等效球面上的等效电磁流转换到在抛物面天线局部直角坐标系下的虚拟等效球面上的等效电磁流的旋转矩阵;
如图6所示,确定将抛物面天线旋转轴旋转到旋转对称天线罩的旋转轴的旋转角度θr,即为z轴和z’轴的夹角,将以抛物面天线旋转轴为z’轴的坐标系下的散射电磁流系数旋转到以旋转对称天线罩旋转轴为z轴的新坐标系下对应的初始散射电磁流系数,J通过第一步中的电磁流系数公式(1)求得,ZR为旋转矩阵,旋转过程是绕y轴旋转,使得z’轴旋转到与z轴重合。
第六步,根据第二步~第五步的信息建立求解方程组,解算得到旋转对称天线罩和虚拟等效球面上的等效散射电磁流;
由二到五步得到的信息建立整体求解方程组如下:
其中α为旋转对称天线罩局部直角坐标系下的旋转对称系统需要的模式数,β为抛物面天线局部直角坐标系下的旋转对称系统需要的模式数,使用迭代算法逐个模式求解上式方程,求解精度小于1e-3,最终得到旋转对称天线罩和虚拟等效球面上的等效散射电磁流。
第七步,由互易定理确定雷达散射截面积(如果确定的是辐射场,则是确定辐射方向图);由第六步确定得到的解向量 和 α=1,……,ModD确定远场,将各个模式下的远场信息线性叠加,得到远场信息,从而得到雷达散射截面积。
下面结合具体实施例对本发明做出进一步详细描述。
实施例1
根据本发明所述方法对一个大型旋转对称天线罩和抛物面天线一体化建模仿真,整个仿真过程可以在个人计算机上完成。
根据本发明所述方法对两个都在Z轴上面的半介质球壳、半金属球壳以及一个喇叭天线组合体,其中半金属球壳在半介质球壳内部,剖面图如图7所示,其中半介质球壳外半径是1m,内半径为0.97m,喇叭天线波导部分半径0.04m,高为0.15m,喇叭口下底面半径为0.15m,上底面半径为0.04m,高为0.15m,金属厚度为0.01m,半金属球壳外半径是0.5m,内半径是0.48m,介质球的介电常数为(1.6,0),工作频率为f=600M,两者的球心距为0.5m,抛物面天线的倾斜角度为15°。平面波的入射角度为θ=0°,金属部分剖分密度为1/10λ,介质部分剖分密度为0.06λ,总得剖分线段为81,旋转对称天线罩使用模式数和抛物面天线使用模式数都为10。图8为VV极化(垂直极化)的双站雷达散射截面积结果,该结果与仿真软件FEKO得到的结果吻合很好,证明本发明方法的正确性。本算例若使用基于三角形剖分RWG基函数矩量法电磁散射仿真未知量为62154,内存消耗为1.2Gb,耗时1.1小时,在同样的测试环境下,本发明消耗内存7Mb,耗时113秒。可见本发明可以对旋转对称天线罩和抛物面天线高效地一体化电磁散射仿真。
Claims (9)
1.一种旋转对称天线罩和抛物面天线一体化电磁散射仿真方法,其特征在于,步骤如下:
第一步,建立模型和局部直角坐标系:以旋转对称天线罩对称轴和抛物面天线对称轴的交点为坐标系原点O,在两条对称轴所在的平面分别建立旋转对称天线罩局部直角坐标系XOZ和抛物面天线局部直角坐标系X’OZ’;以抛物面天线局部直角坐标系X’OZ’的坐标系原点O为球心,建立一个完全包围抛物面天线的虚拟等效球面;在旋转对称天线罩局部直角坐标系XOZ下建立旋转对称天线罩的母线和虚拟等效球面的第一母线,在抛物面天线局部直角坐标系X’OZ’下建立抛物面天线的母线和虚拟等效球面的第二母线;然后对各条母线分别进行剖分离散;
第二步,加入激励源,确定激励向量:根据等效原理确定旋转对称天线罩上的等效入射电磁流和虚拟等效球面上的等效入射电磁流;
第三步,分别建立旋转对称天线罩局部直角坐标系XOZ和抛物面天线局部直角坐标系X’OZ’下的阻抗矩阵,并求解各自阻抗矩阵的逆矩阵;
第四步,确定散射矩阵和传输矩阵:散射矩阵是描述抛物面天线局部直角坐标系下,虚拟等效球面上等效入射电磁流和等效散射电磁流之间的关系;传输矩阵是描述旋转对称天线罩局部直角坐标系下,虚拟等效球面上的感应电流在旋转对称天线罩上产生的二次散射场和旋转对称天线罩上的感应电磁流在虚拟等效球面上产生的二次散射场;
第五步,通过坐标系旋转建立旋转对称天线罩局部直角坐标系XOZ和抛物面天线局部直角坐标系X’OZ’之间的关系,将虚拟等效球面上的等效电磁流在两个局部直角坐标系之间相互转换;
第六步,根据第二步~第五步的信息建立求解方程组,解算得到旋转对称天线罩和虚拟等效球面上的等效散射电磁流;
第七步,由互易定理确定雷达散射截面积。
2.根据权利要求1所述的旋转对称天线罩和抛物面天线一体化电磁散射仿真方法,其特征在于,第一步中所述旋转对称天线罩局部直角坐标系XOZ以旋转对称天线罩的旋转轴为Z轴,抛物面天线局部直角坐标系X’OZ’以抛物面天线的旋转轴为Z’轴;旋转对称天线罩完全在虚拟等效球面之外。
3.根据权利要求1所述的旋转对称天线罩和抛物面天线一体化电磁散射仿真方法,其特征在于,第一步中所述对各条母线分别进行剖分离散,具体方法如下:
使用三角基函数描述切向电磁流,使用傅里叶级数展开的轴向基函数描述周向电磁流,目标表面的电磁流可以展开为:
其中J(r)表示旋转对称体表面任意一点r点的电流,M(r)表示r点的磁流,表示电流J(r)使用旋转对称体基函数展开对应于第α个模式的第n个基函数切向方向的展开系数;表示电流J(r)使用旋转对称体基函数展开对应于第α个模式的第n个基函数周向方向的展开系数;表示磁流M(r)使用旋转对称体基函数展开对应于第α个模式的第n个基函数切向方向的展开系数;表示磁流M(r)使用旋转对称体基函数展开对应于第α个模式的第n个基函数周向方向的展开系数;α表示旋转对称体基函数的模式数;旋转对称体基函数表达式为:
其中Tn(t)表示三角基函数,是一维局部基函数,Tn(t)定义在两条相连接的剖分线段上,我们分别称这两条线段为前段和后段,其表达式为:
t表示r点的切向分量;ρ表示r点到z轴的垂直距离;φ表示r点的周向角;ejαφ表示表示傅里叶级数展开对应于第α个模式的指数项;N表示旋转对称体母线剖分段数,即旋转对称体基函数对应的未知量的个数;表示第n个三角基函数对应的前段的起点切向分量,表示第n个三角基函数对应的前段的终点切向分量即后段的起点切向分量,表示第n个三角基函数对应的后段的终点切向分量;Δn表示前段的长度,Δn+1表示后段的长度;
抛物面天线的母线按照十分之一波长λ/10剖分,旋转对称天线罩的母线按照十分之一介质波长剖分,虚拟等效球面的母线按照七分之一波长λ/7剖分,εr表示介质的相对介电常数,其中λ表示入射电磁波在空气中传播时的波长。
4.根据权利要求1所述的旋转对称天线罩和抛物面天线一体化电磁散射仿真方法,其特征在于,第二步中所述的激励向量包括旋转对称天线罩的激励向量VD和虚拟等效球面的激励向量VH,具体方法如下:
(1)确定旋转对称天线罩的激励向量VD
在旋转对称天线罩表面建立PMCHW方程,激励向量为旋转对称天线罩表面的等效入射电磁流,具体公式如下:
其中下标D表示旋转对称天线罩,上标inc表示入射场,表示旋转对称天线罩表面的入射电场,表示旋转对称天线罩表面的入射磁场,表示旋转对称天线罩表面的等效入射电流,表示旋转对称天线罩表面的等效入射磁流,表示旋转对称天线罩的法向方向;
使用旋转对称基函数的共轭做伽辽金测试:
其中下标m表示基函数编号,α表示模式数编号,表示对应于第α个模式的第m个基函数切向方向的测试基函数,表示对应于第α个模式的第m个基函数周向方向的测试基函数,N表示旋转对称体母线剖分线段数,Mod表示旋转对称体需要的总模式数;
对公式(4)两边使用公式(5)进行测试得到旋转对称天线罩D的激励向量VD:
VJ D表示旋转对称天线罩D表面上等效入射电流对应的激励向量,VM D表示旋转对称天线罩D表面上等效入射磁流对应的激励向量,其每个元素为:
其中〈,〉表示内积,ND表示旋转对称天线罩母线剖分线段数,ModD表示旋转对称体天线罩的总模式数;
(2)确定虚拟等效球面的激励向量VH
与公式(4)中确定旋转对称天线罩的激励向量VD的方法相同,以下标H表示虚拟等效球面,那么
表示虚拟等效球面表面的入射电场,表示虚拟等效球面表面的入射磁场,表示虚拟等效球面表面的等效入射电流,表示虚拟等效球面表面的入射磁流,表示虚拟等效球面的法向方向;
使用公式(5)所示同样形式的测试函数对公式(8)两边进行测试,得到虚拟等效面的激励向量VH:
VJ H表示虚拟等效球面H表面上等效J对应的激励向量,VM H表示虚拟等效球面H表面上等效磁流M对应的激励向量,其每个元素为:
其中下标l表示基函数编号,β表示模式数,表示对应于第β个模式的第l个基函数切向方向的测试基函数,表示对应于第β个模式的第l个基函数周向方向的测试基函数,NH表示虚拟等效球面母线剖分线段数,ModH表示虚拟等效球面上入射等效电磁流使用基函数展开需要的最大模式数。
5.根据权利要求1所述的旋转对称天线罩和抛物面天线一体化电磁散射仿真方法,其特征在于,第三步中所述旋转对称天线罩局部直角坐标系XOZ和抛物面天线局部直角坐标系X’OZ’下的阻抗矩阵、及各自阻抗矩阵的逆矩阵的具体建立过程如下:
(1)建立旋转对称天线罩的阻抗矩阵
在旋转对称天线罩上建立PMCHW方程各个模式阻抗矩阵
其中,Z表示各类阻抗矩阵,其下标L、K算子表示积分算子,对应的第一个下标EJ,HJ表示由电流产生的电磁场,EM,HM表示由磁流产生的电磁场;m表示基函数号;α表示模式数;DD表示场源基函数都是在旋转对称天线罩上;上标tt,tφ,φt,φφ分别表示场源基函数的方向,表示第α个模式的第m个基函数切向方向对应的电流,表示第α个模式的第m个基函数切向方向对应的磁流,表示第α个模式的第m个基函数周向方向对应的电流,表示第α个模式的第m个基函数周向方向对应的磁流;η1表示空气中的波阻抗,η2表示介质的波阻抗;L、K算子表示积分算子,对应的第一个下标表示场区所在位置,第二个下标表示源区所在位置;而“1”则表示旋转对称天线罩的外表面即与空气接触的表面,“2”表示内表面即指向内部介质的表面,即L1,1、K1,1则表示外表面上的电磁流在自身产生的场,L2,1,K2,1表示外表面电磁流在内表面产生的场;算子L、K具体积分表达形式如下:
χ(r')表示等效电流J(r')或等效磁流M(r'),它们将使用公式(1)离散,与各个模式电磁流分量相对应,ω表示角频率,ε表示电导率,μ表示磁导率;G(r,r')表示源点r'到场点r的格林函数其表达式为:
逐个确定出各模式数对应的逆矩阵,使用LU分解算法,表示为
(2)建立抛物面天线的模式阻抗矩阵
在抛物面天线上建立电场积分方程各个模式阻抗矩阵
其中上标P表示抛物面天线,Wβ表示β模式对应的测试函数,ModH表示抛物面天线的总模式数,L(Jβ)表示入射模式电流Jβ在抛物面天线P上产生的激励电场,算子L如公式(12)所示;
逐个确定出各模式数对应的逆矩阵,使用LU分解算法,表示为
6.根据权利要求1所述的旋转对称天线罩和抛物面天线一体化电磁散射仿真方法,其特征在于,第四步所述的散射矩阵和传输矩阵的具体确定方法如下:
(1)确定散射矩阵S
在抛物面天线局部直角坐标系X’OZ’下,抛物面天线和虚拟等效球面的第二母线是共轴的,使用旋转对称矩量法思想确定出各个模式独立的散射矩阵S,散射矩阵S描述的是虚拟等效球面上的等效入射电磁流经过抛物面天线作用后在虚拟等效球面上得到的等效散射电磁流之间的关系:
表示虚拟等效球面表面的等效入射电流,表示虚拟等效球面表面的入射磁流,表示虚拟等效球面表面的等效散射电流,表示虚拟等效球面表面的散射磁流;
散射矩阵S由三部分构成:虚拟等效球面上的模式β对应的等效入射电磁流对抛物面天线的作用描述为抛物面天线自身的作用,描述为第三步中求得的抛物面天线阻抗的逆由抛物面天线自身产生的模式β对应的散射场在虚拟等效球面上得到的模式β对应的等效散射电磁流的作用描述为则S矩阵可以写成:
表示虚拟等效球面的法向量;η表示自由空间中的波阻抗;
其中U-1为电流系数求解矩阵,U矩阵表示如下:
其中fαn表示基函数,n表示基函数编号,α表示模式数,Wβm表示测试函数,m表示测试函数编号,β表示模式数,代入具体得到:
其中Tmp、Tnp分别表示对应于第m个测试函数的第p条线段的三角基函数、第n个基函数的第p条线段的三角基函数,Mp表示一个基函数包含的剖分线段数,这里取2,Δp表示对应的剖分线段的长度,由于三角基函数只和相邻的两个基函数有重叠,因此矩阵U是一个三带条型矩阵,可以通过追赶法求解得到电流系数求解矩阵U-1;
L,K积分算子由公式(12)给出,下标H表示虚拟等效球面,P表示抛物面天线表面,即积分算子LPH(Jβ,H)表示虚拟等效球面H外部等效面上的第β模式的电流Jβ,H在抛物面天线P表面产生的β模式的激励电场,积分算子KPH(Mβ,H)表示虚拟等效球面H外部等效面上的第β模式的磁流Mβ,H在抛物面天线P表面产生的β模式的激励磁场;积分算子KHP(Jβ,P)表示抛物面天线P上的β模式的电流Jβ,P在外虚拟等效球面H外部等效面上产生的β模式的激励磁场,积分算子LHP(Jβ,P)表示抛物面天线P上的β模式的电流Jβ,P在外虚拟等效球面H外部等效面上产生的β模式的激励电场;
(2)确定传输矩阵TDH和THD
在旋转对称天线罩局部直角坐标系XOZ下,旋转对称天线罩母线和虚拟等效球面的第二母线是共轴的,使用旋转对称矩量法思想确定出各个模式独立的传输矩阵TDH和THD,TDH表示虚拟等效球面上的感应电流在旋转对称天线罩上产生的二次散射场,THD表示旋转对称天线罩上的感应电磁流在虚拟等效面上产生的二次散射场:
表示旋转对称天线罩表面的法向量;表示虚拟等效球面表面的法向量;η表示自由空间中的波阻抗,L,K积分算子的积分表达式由公式(12)给出,下标H表示虚拟等效球面,D表示天线罩表面,即积分算子KHD(JD)表示旋转对称天线罩D上的电流J在虚拟等效球面H上的激励磁场,积分算子LHD(MD)表示旋转对称天线罩D上的磁流M在虚拟等效球面H上的激励磁场;积分算子LHD(JD)表示旋转对称天线罩D上的电流J在虚拟等效球面H上的激励电场,积分算子KHD(MD)表示旋转对称天线罩D上的磁流M在虚拟等效球面H上的激励电场;积分算子KDH(JH)表示虚拟等效球面H上的电流JH在旋转对称天线罩D上的激励磁场,积分算子LDH(MH)表示虚拟等效球面H上的电流MH在旋转对称天线罩D上的激励磁场;积分算子LDH(JH)表示虚拟等效球面H上的电流JH在旋转对称天线罩D上的激励电场,积分算子KDH(MH)表示虚拟等效球面H上的电流MH在旋转对称天线罩D上的激励电场。
7.根据权利要求1所述的旋转对称天线罩和抛物面天线一体化电磁散射仿真方法,其特征在于,第四步所述的虚拟等效球面上的等效电磁流在两个局部直角坐标系之间相互转换,具体为:在旋转对称天线罩局部直角坐标系下虚拟等效球面上的电磁流表示记为[JHR MHR]T,在抛物面天线局部直角坐标系下虚拟等效球面上的电磁流表示记为[JHP MHP]T,建立如下关系:
其中ZRP和ZPR表示坐标系旋转矩阵,ZRP表示将在抛物面天线局部直角坐标系下的虚拟等效球面上的等效电磁流转换到旋转对称天线罩局部直角坐标系下虚拟等效球面上的等效电磁流的旋转矩阵,ZPR表示将旋转对称天线罩局部直角坐标系下虚拟等效球面上的等效电磁流转换到在抛物面天线局部直角坐标系下的虚拟等效球面上的等效电磁流的旋转矩阵;
确定将抛物面天线旋转轴旋转到旋转对称天线罩的旋转轴的旋转角度θr,即为z轴和z’轴的夹角,将以抛物面天线旋转轴为z’轴的坐标系下的散射电磁流系数旋转到以旋转对称天线罩旋转轴为z轴的新坐标系下对应的初始散射电磁流系数,JHP、MHP通过第一步中的电磁流系数公式(1)求得,ZRP为旋转矩阵,旋转过程是绕y轴旋转,使得z’轴旋转到与z轴重合。
8.根据权利要求1所述的旋转对称天线罩和抛物面天线一体化电磁散射仿真方法,其特征在于,第六步所述的求解方程组如下:
其中,α为旋转对称天线罩局部直角坐标系下的旋转对称系统需要的模式数,β为抛物面天线局部直角坐标系下的旋转对称系统需要的模式数,表示在模式α下虚拟等效球面上的感应电流在旋转对称天线罩上产生的二次散射场,表示在模式α下旋转对称天线罩上的感应电磁流在虚拟等效面上产生的二次散射场,J,M分别表示电流和磁流,下标HR表示在旋转对称天线罩局部直角坐标系下虚拟等效球面上;下标HP表示抛物面天线局部直角坐标系下虚拟等效球面上;下标D表示天线罩表面;上标sca表示散射;上标inc表示入射:表示在抛物面天线局部直角坐标系下虚拟等效球面上的第β个模式下的散射电流;表示在抛物面天线局部直角坐标系下虚拟等效球面上的第β个模式下的散射磁流;表示在旋转对称天线罩表面的第α个模式下的散射电流;表示在旋转对称天线罩表面的第α个模式下的散射磁流;表示在旋转对称天线罩局部直角坐标系下虚拟等效球面上的第α个模式下的散射电流;表示在旋转对称天线罩局部直角坐标系下虚拟等效球面上的第α个模式下的散射磁流,使用迭代算法逐个模式求解上式方程,求解精度小于1e-3,最终得到旋转对称天线罩和虚拟等效球面上的等效散射电磁流。
9.根据权利要求1所述的旋转对称天线罩和抛物面天线一体化电磁散射仿真方法,其特征在于,第七步所述的由互易定理确定雷达散射截面积,具体为:由第六步得到的解向量 和 α=1,……,ModD确定远场,将各个模式下的远场信息线性叠加,得到远场信息,从而得到雷达散射截面积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310120271.6A CN103218487B (zh) | 2013-04-09 | 2013-04-09 | 旋转对称天线罩和抛物面天线一体化电磁散射仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310120271.6A CN103218487B (zh) | 2013-04-09 | 2013-04-09 | 旋转对称天线罩和抛物面天线一体化电磁散射仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103218487A CN103218487A (zh) | 2013-07-24 |
CN103218487B true CN103218487B (zh) | 2016-08-31 |
Family
ID=48816268
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310120271.6A Active CN103218487B (zh) | 2013-04-09 | 2013-04-09 | 旋转对称天线罩和抛物面天线一体化电磁散射仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103218487B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103559356A (zh) * | 2013-11-11 | 2014-02-05 | 中国电子科技集团公司第五十四研究所 | 一种基于gis的散射链路规划方法 |
CN105095541A (zh) * | 2014-05-13 | 2015-11-25 | 南京理工大学 | 圆柱周期特性的介质目标的频域电磁散射特性分析方法 |
CN103969495B (zh) * | 2014-05-22 | 2016-07-06 | 中国人民解放军国防科学技术大学 | 一种线天线电流的仿真测量方法 |
CN104267381B (zh) * | 2014-08-28 | 2016-09-21 | 北京环境特性研究所 | 组合式鼻锥形电磁散射基准点模型 |
CN105630740B (zh) * | 2014-11-07 | 2018-09-25 | 南京理工大学 | 基于矩阵泰勒级数展开的电磁分析方法 |
CN110008545B (zh) * | 2019-03-22 | 2020-10-09 | 北京航空航天大学 | 一种基于仿真建模的旋转驻极体式机械天线性能评估方法 |
CN111766455B (zh) * | 2020-09-02 | 2020-11-17 | 上海霍莱沃电子系统技术股份有限公司 | 基于口径电流法的相控阵天线方向图预测方法及系统 |
CN113609678B (zh) * | 2021-08-06 | 2024-02-20 | 青岛大学 | 一种改进的永磁体等效面电流计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129523A (zh) * | 2011-03-23 | 2011-07-20 | 南京理工大学 | 基于mda和mlssm的分析复杂目标电磁散射的方法 |
CN102156764A (zh) * | 2010-12-27 | 2011-08-17 | 南京理工大学 | 一种分析天线辐射和电磁散射的多分辨预条件方法 |
CN102608589A (zh) * | 2012-03-16 | 2012-07-25 | 南京理工大学 | 基于仿生模式识别理论的雷达目标识别方法 |
-
2013
- 2013-04-09 CN CN201310120271.6A patent/CN103218487B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102156764A (zh) * | 2010-12-27 | 2011-08-17 | 南京理工大学 | 一种分析天线辐射和电磁散射的多分辨预条件方法 |
CN102129523A (zh) * | 2011-03-23 | 2011-07-20 | 南京理工大学 | 基于mda和mlssm的分析复杂目标电磁散射的方法 |
CN102608589A (zh) * | 2012-03-16 | 2012-07-25 | 南京理工大学 | 基于仿生模式识别理论的雷达目标识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103218487A (zh) | 2013-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103218487B (zh) | 旋转对称天线罩和抛物面天线一体化电磁散射仿真方法 | |
CN103217675B (zh) | 多个不共轴旋转对称体电磁散射特性的仿真方法 | |
Fan et al. | Multidomain pseudospectral time-domain simulations of scattering by objects buried in lossy media | |
Vikram et al. | A novel wideband FMM for fast integral equation solution of multiscale problems in electromagnetics | |
Wang et al. | Multisolver domain decomposition method for modeling EMC effects of multiple antennas on a large air platform | |
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 | |
CN103226644A (zh) | 基于柱面等效源区域分解的电磁散射特性仿真方法 | |
CN103279589A (zh) | 基于矩阵嵌套压缩的旋转对称体电磁散射特性仿真方法 | |
CN103425816B (zh) | 快速获取金属旋转对称体电磁散射特性的矩阵抽取方法 | |
Shao et al. | Analyzing large-scale arrays using tangential equivalence principle algorithm with characteristic basis functions | |
CN113158485A (zh) | 近场条件下电大尺寸目标的电磁散射仿真方法 | |
CN105277927B (zh) | 飞行器编队瞬态电磁特性时域阶数步进分析方法 | |
CN104915326A (zh) | 基于等效原理的区域分解阶数步进时域积分方法 | |
CN109212498B (zh) | 一种飞机编队雷达散射截面的快速算法 | |
Chang et al. | Quadrilateral barycentric basis functions for surface integral equations | |
CN106156475A (zh) | 电大尺寸目标的瞬态电磁特性快速提取方法 | |
Liu et al. | Large-scale electromagnetic computation for modeling and applications | |
CN104778286A (zh) | 掠海飞行器电磁散射特性快速仿真方法 | |
CN103246827B (zh) | 复杂外形金属目标电磁散射的无网格仿真方法 | |
CN105205299B (zh) | 电大目标电磁散射特性快速降维分析方法 | |
Langston et al. | Massively parallel frequency domain electromagnetic simulation codes | |
Dong et al. | Fast and accurate analysis of broadband RCS using method of moments with loop‐tree basis functions | |
Chen et al. | Accelerating the multilevel fast multipole method with parallel preconditioner for large-scale scattering problems | |
Cheng et al. | An efficient high order plane wave time domain algorithm for transient electromagnetic scattering analysis | |
Kireeva et al. | A coupled EFGM–CIE method for acoustic radiation |
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 |