CN103412988B - 基于相移降阶模型周期结构的三维电磁场仿真模拟方法 - Google Patents
基于相移降阶模型周期结构的三维电磁场仿真模拟方法 Download PDFInfo
- Publication number
- CN103412988B CN103412988B CN201310330253.0A CN201310330253A CN103412988B CN 103412988 B CN103412988 B CN 103412988B CN 201310330253 A CN201310330253 A CN 201310330253A CN 103412988 B CN103412988 B CN 103412988B
- Authority
- CN
- China
- Prior art keywords
- formula
- periodic structure
- phase shift
- equation
- electromagnetic 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
Landscapes
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于相移降阶模型周期结构的三维电磁场仿真模拟方法,通过在单个相移点下将有限元法获得的原本征方程进行Taylor级数展开,缩减为维数更低的本征方程,这一过程也就是建立降阶模型。一旦维数更低的本征方程即降阶模型建立后,就能非常迅速地计算多个相移点下周期结构的本征方程,从而获得宽频带内周期结构的高频特性。利用本发明提出的基于相移降阶模型周期结构的三维电磁场仿真模拟方法可以在一个相移点的情形下对周期结构进行有限元本征分析,能够快速精确地获得周期结构在整个带宽内的高频特性。
Description
技术领域
本发明属于三维电磁场数值求解技术领域,涉及一种基于相移降阶模型周期结构的三维电磁场仿真模拟方法。
背景技术
周期结构在微波管中应用非常广泛,包括波导截面的周期性变化,波导周期加载膜片,周期填充介质等。在微波管中,普遍采用周期结构作为器件的高频电路,形成电子注与高频场相互作用进行能量交换以实现微波振荡或放大的场所。周期结构的高频特性:包括色散特性、阻抗特性与衰减特性,直接影响器件的工作频率、频带宽度、换能效率和输出功率,以及其他一系列整管性能。高精度地获得周期结构的高频特性有着极其重要的意义。
随着计算机技术的发展,利用有限元法对微波管周期结构进行三维本征分析,从而获得其高频特性,已经成为微波管周期结构研究和设计中最常用和最有效的方法之一。这类方法的主要过程是用准周期边界条件将微波管周期结构进行截断得到一个周期长度的计算区域,然后对该计算区域进行三维网格离散,采用有限元法可以将周期结构的电磁问题转换成一个大型广义本征值方程,对求解该本征值方程所得到的本征值和本征向量进行一系列的后处理就能最终获得周期结构的高频特性。
目前,周期结构有限元本征分析算法通常在多个不同的相移点下对周期结构进行本征分析,得到各个相移点对应的高频特性。然后利用插值或拟合等数值技术得到整个带宽内的高频特性曲线。对具有宽频带特性的周期结构,这些离散相移点的数量必须足够多,才能保证插值或拟合的精度。这直接导致了周期结构高频特性的整体仿真效率很低,从而严重影响高性能的周期结构的设计。
发明内容
本发明的目的是克服现有技术中存在的缺陷,提供一种基于相移降阶模型周期结构的三维电磁场仿真模拟方法,该方法通过在单个相移点下将有限元法获得的原本征方程进行Taylor级数展开,缩减为维数更低的本征方程,这一过程也就是建立降阶模型。一旦维数更低的本征方程即降阶模型建立后,就能非常迅速地计算多个相移点下周期结构的本征方程,从而获得宽频带内周期结构的高频特性。
其技术方案为:
一种基于相移降阶模型周期结构的三维电磁场仿真模拟方法,包括以下步骤:
A.建立周期结构内的电磁场边值问题,通过有限元法的标准变分原理得到电磁场边值问题的泛函方程;
B.采用四面体网格剖分求解域,并保证周期边界主面和从面上的网格匹配;
C.选择基函数,将电场强度矢量在所有网格内用基函数展开,并运用里兹方法得到周期结构有限元本征分析中的广义本征值方程;
D.将步骤C所得到的广义本征值方程中的所有变量用Taylor级数展开;
E.获得步骤D中本征向量的Taylor展开系数构成的子空间K;
F.根据步骤E获得的子空间K,建立周期结构有限元本征分析中的相移降阶模型;
G.利用步骤F中获得的周期结构的相移降阶模型,可以快速计算周期结构在给定相移范围内的频率以及对应的电场,通过进一步后处理可得周期结构高频特性。
进一步优选,所述步骤A中,首先,根据麦克斯韦方程组与周期结构的边界条件与导体属性得到周期结构内电磁场的边值问题,如下:
公式(1)中第一个式子为频域矢量波动方程,是周期结构有限元仿真中的主方程;其中,Ω为周期结构的仿真区域空间范围,即为公式(1)的求解域,是矢性偏微分算子符号,μr为求解域Ω中介质的相对磁导率,E为求解域Ω的电场强度矢量,k0为自由空间波数,εr为求解域Ω中介质的相对介电常数;
公式(1)中第二个式子为准周期边界条件,其中,ΓPBC表示准周期边界,由主面及从面组成,其中主面定义为沿着周期结构周期性的方向上仿真区域的初始端面,从面定义为沿着周期结构周期性的方向上仿真区域的最后端面,主面与从面的距离为一个空间周期,为边界的外法向单位矢量;Em和Es分别表示周期边界主面和从面上的电场;j为虚数单位符号;β为相位常数;L为周期长度;βL即为一个空间周期对应的相移Φ,准周期边界条件的物理意义是在周期边界从面上的电磁场和主面上的电磁场,相差一个复数相位系数e-jβL;
公式(1)中第三个式子为理想导体的电壁边界条件,其中,ΓPEC表示电壁边界,理想导体的电壁边界条件的物理意义是理想导体的切向电场为零;
周期边界ΓPBC和电壁边界ΓPEC组成了求解域Ω的外边界;
从周期结构内电磁场的边值问题,即从公式(1)出发,通过有限元法的标准变分原理得到电磁场边值问题的泛函方程F(E),即公式(2):
公式(2)中,上标*表示对物理量取共轭,dΩ表示三维体积分的微元,使泛函方程即公式(2)取极小值,并且满足公式(1)中的第二个方程的电场函数E即为周期结构内电磁场边值问题即公式(1)的解。
进一步优选,所述步骤B中采用四面体网格剖分求解域是有限元方法,剖分后的求解域被人为分割为多个三维四面体网格,从而将连续的几何结构空间转化为离散的网格空间。
进一步优选,所述步骤C中选择合适的基函数Nm,将公式(2)和公式(1)中第二个方程中的电场E在所有网格内用基函数Nm展开,即
公式(3)中,下标m取值从0到MT,MT为网格中所有基函数的个数,上标t将电场展开系数xm和基函数Nm按照所在区域进行区分,t∈{I,M,S},I表示基函数Ni在除周期边界外的计算区域内,M表示基函数Nm在周期边界主面上,S表示基函数Nm在周期边界从面上;
将所有网格内的电场E用基函数Nm展开后,代入泛函方程,即公式(2),运用里兹方法,并利用周期边界条件将展开系数用替代,得到周期结构有限元本征分析中的广义本征值方程;
最后得到周期结构有限元本征分析中的广义本征值方程为:
公式(4)中,为广义本征值方程的本征值,k0=ω/c为自由空间波数,ω为角频率,c为自由空间光速;x为广义本征值方程的本征向量,也即公式(3)中插值系数组成的向量,A和B为大型稀疏M×M维有限元矩阵,M等于所采用基函数Nm的总数,有限元矩阵A和B的每一项Aij和Bij可以分别由公式(5)和公式(6)计算可得:
公式(5)和公式(6)中,Ni和Nj为插值基函数,dΩ表示三维体积分的微元,Ω为周期结构的仿真区域空间范围,∫Ω(·)dΩ表示在仿真区域空间范围内的体积分,Wi与Wj为周期边界条件因子,如果插值基函数为从面上的基函数,则对应的周期边界条件因子W=e-jβL,否则W=1;
从公式(5)和公式(6)可以知道,有限元矩阵A和B是相移Φ=βL的函数。
进一步优选,所述步骤D中从步骤C得到了周期结构有限元本征分析的广义本征值方程,
Ax=λBx(7)
将公式(7)中的所有变量分别采用Taylor级数展开,得到:
公式(8)中的A0、A1、A2…,B0、B1、B2…,λ0、λ1、λ2…,x0、x1、x2…为采用Taylor级数展开后的多项式的展开系数,当对相移Φ=βL进行Taylor级数展开时,u=Φ=βL,ue为展开相移点;
有限元矩阵A和B中的元素分别由公式(5)和公式(6)给出,是相移u=βL的函数,因此有限元矩阵A和B关于相移u=βL的Taylor级数展开系数A0、A1、A2…,B0、B1、B2…可以直接求得;
将公式(8)代入至公式(7)中展开后,对比系数即可以得到:
A0x0=λ0B0x0(9)
(A0-λ0B0)x1=λ1B0x0-(A1-λ0B1)x0(10)
(A0-λ0B0)x2=λ2B0x0-(A1-λ0B1)x1-(A2-λ0B2)x0+λ1(B0x1+B1x0)(11)。
进一步优选,所述步骤E中在展开点处求解广义本征方程公式(9),可以获得x0以及λ0,然后在公式(10)中,两边同时乘以x0的厄密向量后,公式(10)的左端就变为零,于是通过计算就可以得到λ1,再将λ1代入公式(10)求解关于矩阵(A0-λ0B0)的方程可以得到x1,接着针对公式(11)按照相同的方式递归求解,最终可以获得由本征向量展开系数构成的子空间K=[x0,x1,…,xN]。
进一步优选,所述步骤F中对K矩阵进行奇异值分解获得正交归一化矩阵Q,将公式(7)两端同时乘以Q的厄密矩阵QH,并利用正交归一化矩阵Q满足的QHQ=I,I为单位矩阵,可以得到相移降阶模型
其中
由于矩阵Q为M×N维矩阵,QH为N×M矩阵,A(u)和B(u)都为M×M维矩阵,由公式(13)可知,和为N×N矩阵,N是公式(8)中Taylor级数的阶数,而M为步骤C中所采用基函数Nm的总数,N远远小于M,广义本征值方程的相移降阶模型即公式(12)的维数,和原广义本征值方程即公式(7)的维数相比较非常小,利用广义本征值方程的相移降阶模型就可以在改变相移点u下快速获得相应的本征值以及本征向量了,实现周期结构的快速扫相移本征分析。
进一步优选,所述步骤G中在指定的相移范围内取一个确定的相位值u=Φ=βL,代入相移降阶模型即公式(12)并求解该广义本征方程,可得本征值λ(u)和本征向量,利用公式(14)(15)得到周期结构在指定相移下的波数k0和对应的电场基函数插值系数组成的本征向量x,
根据指定相移对应的波数k0即可获得角频率ω,根据对应的电场基函数插值系数组成的本征向量x,利用公式(3)得到电场,最后通过后处理可以进一步得到周期结构的高频特性。
与现有技术相比,本发明的有益效果:
利用本发明提出的基于相移降阶模型周期结构的三维电磁场仿真模拟方法可以在一个相移点的情形下对周期结构进行有限元本征分析,能够快速精确地获得周期结构在整个带宽内的高频特性。
附图说明
图1是本发明基于相移降阶模型周期结构的三维电磁场仿真模拟方法的流程图。
具体实施方式
下面结合附图和具体实施例来详细描述本发明的技术方案。
参照图1,一种基于相移降阶模型周期结构的三维电磁场仿真模拟方法,包括以下步骤:
A.建立周期结构内的电磁场边值问题,通过有限元法的标准变分原理得到电磁场边值问题的泛函方程。
首先,根据麦克斯韦方程组与周期结构的边界条件与导体属性得到周期结构内电磁场的边值问题,如下:
公式(1)中第一个式子为频域矢量波动方程,是周期结构有限元仿真中的主方程;其中,Ω为周期结构的仿真区域空间范围,即为公式(1)的求解域。是矢性偏微分算子符号,μr为求解域Ω中介质的相对磁导率,E为求解域Ω的电场强度矢量,k0为自由空间波数,εr为求解域Ω中介质的相对介电常数。
公式(1)中第二个式子为准周期边界条件,其中,ΓPBC表示准周期边界,由主面及从面组成,其中主面定义为沿着周期结构周期性的方向上仿真区域的初始端面,从面定义为沿着周期结构周期性的方向上仿真区域的最后端面。主面与从面的距离为一个空间周期。为边界的外法向单位矢量;Em和Es分别表示周期边界主面和从面上的电场;j为虚数单位符号;β为相位常数;L为周期长度;βL即为一个空间周期对应的相移Φ。准周期边界条件的物理意义是在周期边界从面上的电磁场和主面上的电磁场,相差一个复数相位系数e-jβL;
公式(1)中第三个式子为理想导体的电壁边界条件,其中,ΓPEC表示电壁边界。理想导体的电壁边界条件的物理意义是理想导体的切向电场为零。
周期边界ΓPBC和电壁边界ΓPEC组成了求解域Ω的外边界;
从周期结构内电磁场的边值问题,即从公式(1)出发,通过有限元法的标准变分原理得到电磁场边值问题的泛函方程F(E),即公式(2):
公式(2)中,上标*表示对物理量取共轭,dΩ表示三维体积分的微元。使泛函方程即公式(2)取极小值,并且满足公式(1)中的第二个方程的电场函数E即为周期结构内电磁场边值问题即公式(1)的解。
B.采用四面体网格剖分求解域,并保证周期边界主面和从面上的网格匹配。
采用四面体网格剖分求解域是有限元方法中的一种公知过程,因此本步骤不再详细描述。剖分后的求解域被人为分割为多个三维四面体网格,从而将连续的几何结构空间转化为离散的网格空间。
C.选择基函数,将电场强度矢量在所有网格内用基函数展开,并运用里兹方法得到周期结构有限元本征分析中的广义本征值方程。
选择合适的基函数Nm,将公式(2)和公式(1)中第二个方程中的电场E在所有网格内用基函数Nm展开,即
公式(3)中,下标m取值从0到MT,MT为网格中所有基函数的个数。上标t将电场展开系数xm和基函数Nm按照所在区域进行区分。t∈{I,M,S},I表示基函数Ni在除周期边界外的计算区域内,M表示基函数Nm在周期边界主面上,S表示基函数Nm在周期边界从面上。
将所有网格内的电场E用基函数Nm展开后,代入泛函方程,即公式(2),运用里兹方法,并利用周期边界条件将展开系数用替代,得到周期结构有限元本征分析中的广义本征值方程。该过程为电磁场有限元分析中众所周知的过程,这里不再赘述。
最后得到周期结构有限元本征分析中的广义本征值方程为:
公式(4)中,为广义本征值方程的本征值,k0=ω/c为自由空间波数,ω为角频率,c为自由空间光速;x为广义本征值方程的本征向量,也即公式(3)中插值系数组成的向量。A和B为大型稀疏M×M维有限元矩阵,M等于所采用基函数Nm的总数。有限元矩阵A和B的每一项Aij和Bij可以分别由公式(5)和公式(6)计算可得:
公式(5)和公式(6)中,Ni和Nj为插值基函数,dΩ表示三维体积分的微元,Ω为周期结构的仿真区域空间范围,∫Ω(·)dΩ表示在仿真区域空间范围内的体积分。Wi与Wj为周期边界条件因子。如果插值基函数为从面上的基函数,则对应的周期边界条件因子W=e-jβL,否则W=1。其他符号的含义与公式(1)同。
从公式(5)和公式(6)可以知道,有限元矩阵A和B是相移Φ=βL的函数。
D.将步骤C所得到的广义本征值方程中的所有变量用Taylor级数展开。
步骤C得到了周期结构有限元本征分析的广义本征值方程,
Ax=λBx(7)
将公式(7)中的所有变量分别采用Taylor级数展开,得到:
公式(8)中的A0、A1、A2…,B0、B1、B2…,λ0、λ1、λ2…,x0、x1、x2…为采用Taylor级数展开后的多项式的展开系数。当对相移Φ=βL进行Taylor级数展开时,u=Φ=βL,ue为展开相移点。
有限元矩阵A和B中的元素分别由公式(5)和公式(6)给出,是相移u=βL的函数,因此有限元矩阵A和B关于相移u=βL的Taylor级数展开系数A0、A1、A2…,B0、B1、B2…可以直接求得。
将公式(8)代入至公式(7)中展开后,对比系数即可以得到:
A0x0=λ0B0x0(9)
(A0-λ0B0)x1=λ1B0x0-(A1-λ0B1)x0(10)
(A0-λ0B0)x2=λ2B0x0-(A1-λ0B1)x1-(A2-λ0B2)x0+λ1(B0x1+B1x0)(11)
E.获得步骤D中本征向量的Taylor展开系数构成的子空间K。
在展开点处求解广义本征方程公式(9),可以获得x0以及λ0。然后在公式(10)中,两边同时乘以x0的厄密向量后,公式(10)的左端就变为零,于是通过计算就可以得到λ1。再将λ1代入公式(10)求解关于矩阵(A0-λ0B0)的方程可以得到x1。接着针对公式(11)按照相同的方式递归求解,最终可以获得由本征向量展开系数构成的子空间K=[x0,x1,…,xN]。
F.根据步骤E获得的子空间K,建立周期结构有限元本征分析中的相移降阶模型。
对K矩阵进行奇异值分解获得正交归一化矩阵Q。将公式(7)两端同时乘以Q的厄密矩阵QH,并利用正交归一化矩阵Q满足的QHQ=I,I为单位矩阵,可以得到相移降阶模型
其中
由于矩阵Q为M×N维矩阵,QH为N×M矩阵,A(u)和B(u)都为M×M维矩阵,由公式(13)可知,和为N×N矩阵。N是公式(8)中Taylor级数的阶数,而M为步骤C中所采用基函数Nm的总数。N远远小于M。因此广义本征值方程的相移降阶模型即公式(12)的维数,和原广义本征值方程即公式(7)的维数相比较非常小。利用广义本征值方程的相移降阶模型就可以在改变相移点u下快速获得相应的本征值以及本征向量了,实现周期结构的快速扫相移本征分析。
G.利用步骤G中获得的周期结构的相移降阶模型,可以快速计算周期结构在给定相移范围内的波数k0以及对应的电场基函数插值系数组成的本征向量x,通过进一步后处理可得周期结构的高频特性。
在指定的相移范围内取一个确定的相位值u=Φ=βL,代入相移降阶模型即公式(12)并求解该广义本征方程,可得本征值λ(u)和本征向量,利用公式(14)(15)得到周期结构在指定相移下的波数k0和对应的电场基函数插值系数组成的本征向量x。
根据指定相移对应的波数k0即可获得角频率ω。根据对应的电场基函数插值系数组成的本征向量x,利用公式(3)得到电场。最后通过后处理可以进一步得到周期结构的高频特性。该过程为本领域的公知过程,因此不再详细描述。
以上所述,仅为本发明较佳的具体实施方式,本发明的保护范围不限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可显而易见地得到的技术方案的简单变化或等效替换均落入本发明的保护范围内。
Claims (8)
1.一种基于相移降阶模型周期结构的三维电磁场仿真模拟方法,其特征在于,包括以下步骤:
A.建立周期结构内的电磁场边值问题,通过有限元法的标准变分原理得到电磁场边值问题的泛函方程;
B.采用四面体网格剖分求解域,并保证周期边界主面和从面上的网格匹配;
C.选择基函数,将电场强度矢量在所有网格内用基函数展开,并运用里兹方法得到周期结构有限元本征分析中的广义本征值方程;
D.将步骤C所得到的广义本征值方程中的所有变量用Taylor级数展开;
E.获得步骤D中本征向量的Taylor展开系数构成的子空间K;
F.根据步骤E获得的子空间K,建立周期结构有限元本征分析中的相移降阶模型;
G.利用步骤F中获得的周期结构的相移降阶模型,可以快速计算周期结构在给定相移范围内的频率以及对应的电场,通过进一步后处理可得周期结构高频特性。
2.根据权利要求1所述的基于相移降阶模型周期结构的三维电磁场仿真模拟方法,其特征在于,所述步骤A中,首先,根据麦克斯韦方程组与周期结构的边界条件与导体属性得到周期结构内电磁场的边值问题,如下:
公式(1)中第一个式子为频域矢量波动方程,是周期结构有限元仿真中的主方程;其中,Ω为周期结构的仿真区域空间范围,即为公式(1)的求解域,是矢性偏微分算子符号,μr为求解域Ω中介质的相对磁导率,E为求解域Ω的电场强度矢量,k0为自由空间波数,εr为求解域Ω中介质的相对介电常数;
公式(1)中第二个式子为准周期边界条件,其中,ΓPBC表示准周期边界,由主面及从面组成,其中主面定义为沿着周期结构周期性的方向上仿真区域的初始端面,从面定义为沿着周期结构周期性的方向上仿真区域的最后端面,主面与从面的距离为一个空间周期,为边界的外法向单位矢量;Em和Es分别表示周期边界主面和从面上的电场;j为虚数单位符号;β为相位常数;L为周期长度;βL即为一个空间周期对应的相移Φ,准周期边界条件的物理意义是在周期边界从面上的电磁场和主面上的电磁场,相差一个复数相位系数e-jβL;
公式(1)中第三个式子为理想导体的电壁边界条件,其中,ΓPEC表示电壁边界,理想导体的电壁边界条件的物理意义是理想导体的切向电场为零;
周期边界ΓPBC和电壁边界ΓPEC组成了求解域Ω的外边界;
从周期结构内电磁场的边值问题,即从公式(1)出发,通过有限元法的标准变分原理得到电磁场边值问题的泛函方程F(E),即公式(2):
公式(2)中,上标*表示对物理量取共轭,dΩ表示三维体积分的微元,使泛函方程即公式(2)取极小值,并且满足公式(1)中的第二个方程的电场函数E即为周期结构内电磁场边值问题即公式(1)的解。
3.根据权利要求1所述的基于相移降阶模型周期结构的三维电磁场仿真模拟方法,其特征在于,所述步骤B中采用四面体网格剖分求解域是有限元方法,剖分后的求解域被人为分割为多个三维四面体网格,从而将连续的几何结构空间转化为离散的网格空间。
4.根据权利要求2所述的基于相移降阶模型周期结构的三维电磁场仿真模拟方法,其特征在于,所述步骤C中选择合适的基函数Nm,将公式(2)和公式(1)中第二个方程中的电场E在所有网格内用基函数Nm展开,即
公式(3)中,下标m取值从0到MT,MT为网格中所有基函数的个数,上标t将电场展开系数xm和基函数Nm按照所在区域进行区分,t∈{I,M,S},I表示基函数Nm在除周期边界外的计算区域内,M表示基函数Nm在周期边界主面上,S表示基函数Nm在周期边界从面上;
将所有网格内的电场E用基函数Nm展开后,代入泛函方程,即公式(2),运用里兹方法,并利用周期边界条件将展开系数用替代,得到周期结构有限元本征分析中的广义本征值方程;
最后得到周期结构有限元本征分析中的广义本征值方程为:
公式(4)中,为广义本征值方程的本征值,k0=ω/c为自由空间波数,ω为角频率,c为自由空间光速;x为广义本征值方程的本征向量,也即公式(3)中插值系数组成的向量,A和B为大型稀疏M×M维有限元矩阵,M等于所采用基函数Nm的总数,有限元矩阵A和B的每一项Aij和Bij可以分别由公式(5)和公式(6)计算可得:
公式(5)和公式(6)中,Ni和Nj为插值基函数,dΩ表示三维体积分的微元,Ω为周期结构的仿真区域空间范围,∫Ω(·)dΩ表示在仿真区域空间范围内的体积分,Wi与Wj为周期边界条件因子,如果插值基函数为从面上的基函数,则对应的周期边界条件因子W=e-jβL,否则W=1;
从公式(5)和公式(6)可以知道,有限元矩阵A和B是相移Φ=βL的函数。
5.根据权利要求4所述的基于相移降阶模型周期结构的三维电磁场仿真模拟方法,其特征在于,所述步骤D中从步骤C得到了周期结构有限元本征分析的广义本征值方程,
Ax=λBx(7)
将公式(7)中的所有变量分别采用Taylor级数展开,得到:
公式(8)中的A0、A1、A2…,B0、B1、B2…,λ0、λ1、λ2…,x0、x1、x2…为采用Taylor级数展开后的多项式的展开系数,当对相移Φ=βL进行Taylor级数展开时,u=Φ=βL,ue为展开相移点;
有限元矩阵A和B中的元素分别由公式(5)和公式(6)给出,是相移u=βL的函数,因此有限元矩阵A和B关于相移u=βL的Taylor级数展开系数A0、A1、A2…,B0、B1、B2…可以直接求得;
将公式(8)代入至公式(7)中展开后,对比系数即可以得到:
A0x0=λ0B0x0(9)
(A0-λ0B0)x1=λ1B0x0-(A1-λ0B1)x0(10)
(A0-λ0B0)x2=λ2B0x0-(A1-λ0B1)x1-(A2-λ0B2)x0+λ1(B0x1+B1x0)(11)。
6.根据权利要求5所述的基于相移降阶模型周期结构的三维电磁场仿真模拟方法,其特征在于,所述步骤E中在展开点处求解广义本征方程公式(9),可以获得x0以及λ0,然后在公式(10)中,两边同时乘以x0的厄密向量后,公式(10)的左端就变为零,于是通过计算就可以得到λ1,再将λ1代入公式(10)求解关于矩阵(A0-λ0B0)的方程可以得到x1,接着针对公式(11)按照相同的方式递归求解,最终可以获得由本征向量展开系数构成的子空间K=[x0,x1,…,xN]。
7.根据权利要求6所述的基于相移降阶模型周期结构的三维电磁场仿真模拟方法,其特征在于,所述步骤F中对K矩阵进行奇异值分解获得正交归一化矩阵Q,将公式(7)两端同时乘以Q的厄密矩阵QH,并利用正交归一化矩阵Q满足的QHQ=I,I为单位矩阵,可以得到相移降阶模型
其中
由于矩阵Q为M×N维矩阵,QH为N×M矩阵,A(u)和B(u)都为M×M维矩阵,由公式(13)可知,和为N×N矩阵,N是公式(8)中Taylor级数的阶数,而M为步骤C中所采用基函数Nm的总数,N远远小于M,广义本征值方程的相移降阶模型即公式(12)的维数,和原广义本征值方程即公式(7)的维数相比较非常小,利用广义本征值方程的相移降阶模型就可以在改变相移点u下快速获得相应的本征值以及本征向量了,实现周期结构的快速扫相移本征分析。
8.根据权利要求7所述的基于相移降阶模型周期结构的三维电磁场仿真模拟方法,其特征在于,所述步骤G中在指定的相移范围内取一个确定的相位值u=Φ=βL,代入相移降阶模型即公式(12)并求解该广义本征方程,可得本征值λ(u)和本征向量利用公式(14)(15)得到周期结构在指定相移下的波数k0和对应的电场基函数插值系数组成的本征向量x,
根据指定相移对应的波数k0即可获得角频率ω,根据对应的电场基函数插值系数组成的本征向量x,利用公式(3)得到电场,最后通过后处理可以进一步得到周期结构的高频特性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310330253.0A CN103412988B (zh) | 2013-08-01 | 2013-08-01 | 基于相移降阶模型周期结构的三维电磁场仿真模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310330253.0A CN103412988B (zh) | 2013-08-01 | 2013-08-01 | 基于相移降阶模型周期结构的三维电磁场仿真模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103412988A CN103412988A (zh) | 2013-11-27 |
CN103412988B true CN103412988B (zh) | 2016-07-06 |
Family
ID=49606000
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310330253.0A Active CN103412988B (zh) | 2013-08-01 | 2013-08-01 | 基于相移降阶模型周期结构的三维电磁场仿真模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103412988B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106294975B (zh) * | 2016-08-05 | 2019-01-01 | 大连理工大学 | 一种基于降阶模型的梁式结构自由振动分析方法 |
CN107526869B (zh) * | 2017-07-18 | 2020-03-27 | 电子科技大学 | 一种基于函数逼近自适应三维微波管输入输出窗模型降阶的数值方法 |
CN107992696B (zh) * | 2017-12-13 | 2020-12-29 | 电子科技大学 | 一种复杂色散媒质中的改进的指数时间积分构造方法 |
CN108446429B (zh) * | 2018-02-05 | 2021-07-06 | 电子科技大学 | 一种应用于pic静电模型的粒子受力有限元求解算法 |
CN108873084B (zh) * | 2018-05-10 | 2019-10-08 | 中南大学 | 一种基于单位分解积分的直流电阻率无单元正演方法 |
CN111931458B (zh) * | 2020-10-09 | 2021-01-15 | 北京智芯仿真科技有限公司 | 基于混合阶有限元的三维集成电路电磁场计算方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101944145A (zh) * | 2010-08-31 | 2011-01-12 | 电子科技大学 | 可去除伪直流模式的微波管高频电路有限元仿真方法 |
CN102592057A (zh) * | 2012-01-17 | 2012-07-18 | 电子科技大学 | 周期结构指定频率的本征分析方法 |
CN102930071A (zh) * | 2012-08-29 | 2013-02-13 | 电子科技大学 | 基于非匹配网格的周期结构的三维电磁场仿真模拟方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050021319A1 (en) * | 2003-06-03 | 2005-01-27 | Peng Li | Methods, systems, and computer program products for modeling nonlinear systems |
-
2013
- 2013-08-01 CN CN201310330253.0A patent/CN103412988B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101944145A (zh) * | 2010-08-31 | 2011-01-12 | 电子科技大学 | 可去除伪直流模式的微波管高频电路有限元仿真方法 |
CN102592057A (zh) * | 2012-01-17 | 2012-07-18 | 电子科技大学 | 周期结构指定频率的本征分析方法 |
CN102930071A (zh) * | 2012-08-29 | 2013-02-13 | 电子科技大学 | 基于非匹配网格的周期结构的三维电磁场仿真模拟方法 |
Non-Patent Citations (2)
Title |
---|
High-Frequency Circuit Simulator: An Advanced Three-Dimensional Finite;Li Xu,et al;《IEEE TRANSACTIONS ON ELECTRON DEVICES》;20090531;第56卷(第5期);1141-1150 * |
微波管高频电路快速有限元理论与CAD技术研究;徐立;《中国博士学位论文全文数据库 信息科技辑》;20121215;第I135-8页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103412988A (zh) | 2013-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103412989B (zh) | 基于参数化降阶模型周期结构的三维电磁场仿真模拟方法 | |
CN103412988B (zh) | 基于相移降阶模型周期结构的三维电磁场仿真模拟方法 | |
CN102930071B (zh) | 基于非匹配网格的周期结构的三维电磁场仿真模拟方法 | |
CN103279601B (zh) | 导体目标宽带电磁散射特性的仿真方法 | |
CN102592057B (zh) | 周期结构指定频率的本征分析方法 | |
Huang et al. | Topological strings and quantum spectral problems | |
CN101944145B (zh) | 可去除伪直流模式的微波管高频电路有限元仿真方法 | |
Slone et al. | Multipoint Galerkin asymptotic waveform evaluation for model order reduction of frequency domain FEM electromagnetic radiation problems | |
CN102609598B (zh) | 一种对大型电力系统进行电磁暂态仿真的方法 | |
Zhou et al. | Node-based smoothed radial point interpolation method for electromagnetic-thermal coupled analysis | |
CN102156764A (zh) | 一种分析天线辐射和电磁散射的多分辨预条件方法 | |
CN103970717A (zh) | 基于Associated Hermite 正交函数的无条件稳定FDTD算法 | |
Christophe et al. | An implicit hybridized discontinuous Galerkin method for the 3D time-domain Maxwell equations | |
CN103425816B (zh) | 快速获取金属旋转对称体电磁散射特性的矩阵抽取方法 | |
CN104679920B (zh) | 波导器件微波气体放电数值模拟方法 | |
CN106446470B (zh) | 一种高效并行的非均匀介质频域有限差分方法 | |
CN104050307A (zh) | 基于时域谱元法的微波无源电路电磁热一体化分析方法 | |
CN107515955A (zh) | 基于eb连续‑不连续伽辽金混合的时域有限元方法 | |
CN105468864A (zh) | 基于増维精细积分的高压输电线路电磁暂态数值计算方法 | |
CN102298658A (zh) | 一种行波管注波互作用的模拟方法 | |
CN106777472A (zh) | 基于拉盖尔多项式的减少分裂误差的完全匹配层实现方法 | |
CN104752245A (zh) | 高功率脉冲对场效应管放大器性能影响的数值分析方法 | |
CN104731996A (zh) | 一种快速提取电大尺寸金属腔体目标瞬态散射信号的仿真方法 | |
CN105022850A (zh) | 一种声表面波器件的建模方法 | |
CN105303022A (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 |