CN102930071A - 基于非匹配网格的周期结构的三维电磁场仿真模拟方法 - Google Patents

基于非匹配网格的周期结构的三维电磁场仿真模拟方法 Download PDF

Info

Publication number
CN102930071A
CN102930071A CN2012103123593A CN201210312359A CN102930071A CN 102930071 A CN102930071 A CN 102930071A CN 2012103123593 A CN2012103123593 A CN 2012103123593A CN 201210312359 A CN201210312359 A CN 201210312359A CN 102930071 A CN102930071 A CN 102930071A
Authority
CN
China
Prior art keywords
face
rightarrow
interarea
grid
periodic structure
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
CN2012103123593A
Other languages
English (en)
Other versions
CN102930071B (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 CN201210312359.3A priority Critical patent/CN102930071B/zh
Publication of CN102930071A publication Critical patent/CN102930071A/zh
Application granted granted Critical
Publication of CN102930071B publication Critical patent/CN102930071B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明涉及基于非匹配网格的周期结构的三维电磁场仿真模拟方法,包括以下步骤:选取特定的具有周期特性的微波管高频电路;从步骤A中选取的高频电路中截取一个周期长度的结构进行建模,建立该周期长度的高频结构对应的几何结构模型;根据几何结构模型的周期性确定主面和从面,并对所建几何结构模型进行网格划分,将连续的几何结构空间转化为离散的网格空间;根据仿真区域主面面网格与从面面网格生成联合面网格。本发明的有益效果:利用本发明提出的基于非匹配网格的周期结构的三维电磁场仿真模拟方法可以在对网格没有任何限制的条件下,准确快速求解周期结构的高频特性。

Description

基于非匹配网格的周期结构的三维电磁场仿真模拟方法
技术领域
本发明属于三维电磁场数值求解的技术领域,具体涉及一种基于非匹配网格的周期结构本征分析方法。
背景技术
周期结构在微波管中应用非常广泛,包括波导截面的周期性变化,波导周期加载膜片,周期填充介质等。在微波管中,普遍采用周期结构作为器件的高频电路,形成电子注与高频场相互作用进行能量交换以实现微波振荡或放大的场所。周期结构的高频特性(包括色散特性、阻抗特性与衰减特性)直接影响器件的工作频率、频带宽度、换能效率和输出功率,以及其他一系列整管性能。高精度地获得周期结构的高频特性有着极其重要的意义。
目前,利用各种计算电磁学方法对周期结构的高频特性进行仿真分析时,通常利用一定的周期边界条件将仿真区域缩减为一个空间周期,通过对周期结构的一个空间周期进行三维建模和网格划分,在严格要求仿真区域主面和从面上的网格严格匹配的情况下,求解特定边界条件下的麦克斯韦方程组或其等效形式,求出仿真区域内的电磁场分布以及色散、互作用阻抗以及衰减常数等高频特性。
当几何结构比较规则,采用简单的结构化网格(如正六面体单元)就能很好地模拟该结构时,上述的模拟方法是完全可行的。但是,当几何结构日益复杂,日益不规则时,必须采用非结构化的网格单元(如四面体单元,三棱柱等)来精确模拟任意几何形状。通常的做法是采用结构化网格来模拟几何结构的规则部分,而采用四面体单元来模拟不规则部分。当引入非结构化的四面体网格进行空间离散时,仿真区域主面与从面网格严格匹配的强制性要求会导致生成的网格质量非常差,甚至无法完成网格划分。仿真区域主面与从面网格匹配的强制性要求也严重限制了高效率和高可靠性的自适应网格技术的应用。
发明内容
本发明的目的是为了克服现有的周期结构数值分析方法对主面与从面网格严格匹配的强制性要求,提出了基于非匹配网格的周期结构的三维电磁场仿真模拟方法。该方法在对网格划分没有任何特殊限制的条件下,能够精确和高效地求解出周期结构的高频特性。
为了实现上述目的,本发明的技术方案是:基于非匹配网格的周期结构的三维电磁场仿真模拟方法,包括以下步骤:
A.选取特定的具有周期特性的微波管高频电路;
B.从步骤A中选取的高频电路中截取一个周期长度的结构进行建模,建立该周期长度的高频结构对应的几何结构模型;
C.根据几何结构模型的周期性确定主面和从面,并对所建几何结构模型进行网格划分,将连续的几何结构空间转化为离散的网格空间;
D.根据仿真区域主面面网格与从面面网格生成联合面网格;
E.利用有限元法,将麦克斯韦方程组等效的电磁场边值问题在步骤C建立的网格空间进行离散,结合步骤D建立的联合网格,建立考虑介质损耗和有限电导率导体损耗周期结构高频电路的代数本征方程;
F.给定一个频率,求解步骤E所建立的代数本征方程,获得与给定频率相对应的周期结构的相位常数、衰减常数和互作用阻抗;
G.给定不同的频率,重复步骤F,获得周期结构的高频特性。
本发明的有益效果:利用本发明提出的基于非匹配网格的周期结构的三维电磁场仿真模拟方法可以在对网格没有任何限制的条件下,准确快速求解周期结构的高频特性。
附图说明
图1是本发明的主流程图。
图2是联合网格的构造方法示意图。
图3是周期边界主面和从面上的非匹配网格构成的联合网格。
图4是典型螺旋线高频电路的色散特性曲线。
图5是典型螺旋线高频电路的衰减特性曲线。
图6是典型螺旋线高频电路的互作用阻抗特性曲线。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
如图1所示,基于非匹配网格的周期结构的三维电磁场仿真模拟方法包括以下步骤:
A.选取特定的具有周期特性的微波管高频电路;
选取特定的具有周期特性的微波管高频电路,如螺旋线高频电路、耦合腔高频电路、折叠波导高频电路等。
B.从步骤A中选取的高频电路中截取一个周期长度的结构进行建模,建立该周期长度的结构对应的几何结构模型;
根据高频电路的周期性,通常仅建立一个空间周期几何结构,并引入周期边界条件来仿真整个周期结构高频电路的高频特性。具体的结构建模是电磁场数值计算中的一种公知过程,因此本步骤不再详细描述。
C.根据几何结构模型的周期性确定主面和从面,并对所建几何结构模型进行网格划分,将连续的几何结构空间转化为离散的网格空间;
根据几何结构模型的周期性,将沿着周期结构周期性的方向上仿真区域的初始端面定义为主面,沿着周期结构周期性的方向上仿真区域的最后端面定义为从面。主面与从面的距离为一个空间周期。主面上的场与从面上的场服从周期边界条件。
然后,采用四面体网格剖分仿真区域,剖分后的仿真区域被人为分割为多个三维四面体网格,从而将连续的几何结构空间转化为离散的网格空间。此时,主面和从面由一系列离散的面网格单元组成。这里定义附属于主面上的网格单元为主面面网格,附属于从面上的网格单元为从面面网格。由于四面体网格剖分是有限元方法中的一种公知过程,因此本步骤不再详细描述。
与传统的计算电磁学方法分析周期结构要求主面面网格和从面面网格完全匹配不同,这里周期结构的主面面网格和从面面网格没有任何限制,不再强制要求严格匹配。
D.根据仿真区域主面面网格和从面面网格生成联合面网格;
传统的计算电磁学方法分析周期结构的高频特性通常要求主面面网格和从面面网格完全匹配。在此条件下,根据周期边界条件,从面面网格上的场可由主面面网格上的场直接对应得到。反之,主面面网格上的场也可由从面面网格上的场直接对应得到。因此,在模拟过程中,仅需要引入周期边界条件计算主面面网格上的场或者从面面网格上的场即可。
由于在步骤C的网格划分中,没有强制要求主面面网格和从面面网格严格匹配,主面面网格上的场与从面面网格上的场不存在直接对应的关系。为了在后续步骤中能够施加周期边界条件,需要根据主面面网格和从面面网格生成联合面网格。下面给出联合面网格的生成方法与步骤。
1)将主面面网格和从面面网格进行重叠;
一般情况下,主面上的一个面网格和从面上的一个面网格重叠,重叠部分为一不规则的多边形,如图2所示。其中,粗线三角形代表主面上的一个面网格,细线三角形代表一个从面上的面网格,标号1,2,3,4,5,6的六个三角形构成的不规则多边形为主面上的一个面网格和从面上的一个面网格重叠相互重叠的部分。
与主面上的一个面网格和从面上的一个面网格重叠的情形类似,主面上的所有面网格和从面上的所有面网格重叠将形成大量的不规则多边形。
2)将主面和从面上的面网格重叠后形成的各个不规则多边形内部取点,并将该点和相应不规则多边形的各个顶点连接,将各个不规则多边形重新剖分成多个新的三角形,从而形成主面和从面的联合面网格。
在不规则多边形中选取的一点一般尽可能落在多边形的中部,以保证新形成的多个三角形中不出现钝角。一种可取的方案是取相对两个内角的对角线的交点。如图2所示,一个不规则多边形被重新划分为标号1,2,3,4,5,6的六个三角形。
将主面面网格和从面面网格重叠后形成的每一个不规则多边形按照步骤2)用三角形剖分后,得到典型的联合面网格示意图,如图3所示。显然,联合面网格是对主面或从面更为细致的网格划分。此时,将联合网格中的任一三角形网格单元(如图3中的灰色三角形)分别向主面和从面进行投影,就能唯一确定其所属主面上的某个三角形网格和从面上的某个三角形网格。
E.利用有限元法,将麦克斯韦方程组等效的电磁场边值问题在步骤C建立的网格空间进行离散,结合步骤D建立的联合网格,建立考虑介质损耗和有限电导率导体损耗周期结构高频电路的代数本征方程;
有限元法是一种众所周知的近似求解数理边值问题的数值技术,在电磁学中的应用已有40余年的历史,这里不再赘述。
在微波管中,经常采用周期结构作为高频电路,如螺旋线高频电路和耦合腔高频电路。这些周期结构高频电路通常会包含非理想的介质与有限电导率的导体。由于这些非理想的介质和导体的存在,电磁波在周期结构高频电路中传输时会伴随着损耗。这里,直接给出与麦克斯韦方程组等效的,考虑介质损耗和有限电导率导体损耗周期结构高频电路的电磁场边值问题,如式(1)所示。
▿ × μ r - 1 ▿ × E → - k 0 2 ϵ r * E → = 0 inΩ Z s n ^ × μ r - 1 ▿ × E → = j k 0 η 0 n ^ × ( E → × n ^ ) on Γ SIBC n ^ × E → s × n ^ = n ^ × E ^ m × n ^ e - γL on Γ PBC n ^ × μ r - 1 ▿ × E → s = - ( n ^ × μ r - 1 ▿ × E → m ) e - γL on Γ PBC - - - ( 1 )
式(1)中第一个式子为频域矢量波动方程,是周期结构有限元仿真中的主方程;其中,Ω为周期结构的仿真区域空间范围,即为公式(1)的求解域。是矢性偏微分算子符号,μr为求解域Ω中介质的相对磁导率,
Figure BDA00002073601300043
为求解域Ω的电场强度矢量,k0为自由空间波数,εr为求解域Ω中介质的相对介电常数。在考虑介质损耗时,εr=ε′r(1-jtanδ)为复数,ε′r为εr的实部,通常为不随频率变化的常数。j为虚数单位符号。-jε′rtanδ为εr的虚部,tanδ是用以描述介质损耗的损耗角正切,通常随频率升高而增加。
公式(1)中第二个式子为导体的阻抗边界条件。其中,ΓSIBC表示阻抗边界;Zs为良导体的表面阻抗,满足
Figure BDA00002073601300051
其中f为频率,μ和σ分别为导体的磁导率和有限电导率。
Figure BDA00002073601300052
为边界的外法向单位矢量。η0为自由空间波阻抗。
公式(1)中第三和第四式为周期结构的准周期边界条件。ΓPBC表示准周期边界,
Figure BDA00002073601300053
Figure BDA00002073601300054
分别表示主面和从面上的电场强度;传播常数γ=α+jβ,α和β分别为衰减常数和相位常数;L为周期结构一个空间周期的长度。第三个式子表示主面和从面上的切向电场满足弗洛奎定理;第四个式子表示主面和从面上的切向磁场满足弗洛奎定理。弗洛奎定理是该领域公知的研究周期结构高频电路的基础定理,它表明,在给定频率下,对确定的电磁波传播模式,沿着周期结构传播的波在任一横截面上的场分布与离该截面一个周期远处的场分布仅仅相差一个复常数因子e-(α+jβ)L。因为主面和从面的空间距离刚好为周期结构的一个空间周期,符合弗洛奎定理。
周期边界ΓPBC和阻抗边界ΓSIBC共同组成求解域Ω的外边界。
定义矢量任意矢量函数
Figure BDA00002073601300056
Figure BDA00002073601300057
的三维内积为
Figure BDA00002073601300058
对公式(1)中的矢量波动方程和阻抗边界条件选取测试函数
Figure BDA00002073601300059
依照公知的伽略金有限元法推导过程,得到矢量波动方程与导体的阻抗边界条件的联合伽略金弱形式,如
( &dtri; &times; v &RightArrow; , &mu; r - 1 &dtri; &times; E &RightArrow; ) &Omega; - k 0 2 ( v &RightArrow; , &epsiv; r * E &RightArrow; ) &Omega; + k 0 < v &RightArrow; , j &RightArrow; m > &Gamma; m + k 0 < v &RightArrow; , j &RightArrow; s > &Gamma; s + j k 0 &eta; 0 Z s < v &RightArrow; , n ^ &times; E &RightArrow; &times; n ^ > &Gamma; SIBC = 0 - - - ( 2 )
为了在联合面网格的基础上实现周期边界条件,提出了一种等效的二阶TE周期边界条件,通过这种新的施加方式去除周期结构主面和从面网格匹配的限制。首先,引入辅助的表面电场强度矢量
Figure BDA000020736013000511
和表面电流密度矢量
Figure BDA000020736013000512
e &RightArrow; = n ^ &times; ( E &RightArrow; &times; n ^ ) - - - ( 3 )
j &RightArrow; = k 0 - 1 n ^ &times; &mu; r - 1 &dtri; &times; E &RightArrow; - - - ( 4 )
根据公式(1)中电场和磁场的准周期边界条件,经过一系列代数运算,得到其对应的等价形式,如式(5)和式(6)所示。
( k 0 j &RightArrow; s + a e &RightArrow; s + b &dtri; &Gamma; &times; &dtri; &Gamma; &times; e &RightArrow; s ) = e - &gamma;L ( - k 0 j &RightArrow; m + a e &RightArrow; m + b &dtri; &Gamma; &times; &dtri; &Gamma; &times; e &RightArrow; m ) - - - ( 5 )
( k 0 j &RightArrow; m + a e &RightArrow; m + b &dtri; &Gamma; &times; &dtri; &Gamma; &times; e &RightArrow; m ) = e - &gamma;L ( - k 0 j &RightArrow; s + a e &RightArrow; s + b &dtri; &Gamma; &times; &dtri; &Gamma; &times; e &RightArrow; s ) - - - ( 6 )
公式(5)和公式(6)中,
Figure BDA00002073601300061
Figure BDA00002073601300062
分别表示主面与从面上的表面电场强度,
Figure BDA00002073601300064
分别主面与从面上的表面电流密度。
Figure BDA00002073601300065
表示表面旋度算子,a和b为待定系数,这里取为
a=-jk0                            (7)
b = j k 0 + k ~ z - - - ( 8 )
其中,
Figure BDA00002073601300068
为在给定工作频率和网格尺寸下周期边界上能够支持的最大横向波数。
Figure BDA00002073601300069
c为一常数,这里取c=0.5;h为网格最大尺寸,p为有限元法所选取基函数的阶数。当采用边棱元基函数时,p=1;当采用二阶叠层型矢量基函数时,p=2。
公式(5)和公式(6)中包含了横向电场的高阶项
Figure BDA000020736013000611
因而能够显著提高周期边界上存在的TE凋零模式的迭代收敛速度,因此将该周期边界条件称为二阶TE周期边界条件。TE为电磁场中众所周知的横电模的英文缩写。对TE模,电场仅存在于与电磁波传播方向垂直的横向,沿电磁波传播方向没有电场分量。
要施加二阶TE周期边界条件,需要推导对应伽略金弱形式。对(5)式与(6)式,分别选取测试函数
Figure BDA000020736013000612
按照公知的伽略金有限元法推导过程,得到式(5)与式(6)对应的伽略金弱形式,如式(5)与式(9)所示。
k 0 2 < u &RightArrow; s , a - 1 j &RightArrow; s > &Gamma; s + k 0 < u &RightArrow; s , e &RightArrow; s > &Gamma; s + k 0 < &dtri; &Gamma; &times; u &RightArrow; s , b a - 1 &dtri; &Gamma; &times; e &RightArrow; s > &Gamma; s + e - &gamma;L k 0 2 < u &RightArrow; s , a - 1 j &RightArrow; m > &Gamma; p - e - &gamma;L k 0 < u &RightArrow; s , e &RightArrow; m > &Gamma; p - e - &gamma;L k 0 < &dtri; &Gamma; &times; u &RightArrow; s , b a - 1 &dtri; &Gamma; &times; e &RightArrow; m > &Gamma; p = 0 - - - ( 9 )
k 0 2 < u &RightArrow; m , a - 1 j &RightArrow; m > &Gamma; m + k 0 < u &RightArrow; m , e &RightArrow; m > &Gamma; m + k 0 < &dtri; &Gamma; &times; u &RightArrow; m , b a - 1 &dtri; &Gamma; &times; e &RightArrow; m > &Gamma; m + e &gamma;L k 0 2 < u &RightArrow; m , a - 1 j &RightArrow; s > &Gamma; p - e &gamma;L k 0 < u &RightArrow; m , e &RightArrow; s > &Gamma; p - e &gamma;L k 0 < &dtri; &Gamma; &times; u &RightArrow; m , b a - 1 &dtri; &Gamma; &times; e &RightArrow; s > &Gamma; p = 0 - - - ( 10 )
式(9)与式(10)中关于任意两个矢量函数
Figure BDA000020736013000616
Figure BDA000020736013000617
的二维内积定义为
Figure BDA000020736013000618
Γs表示从面,Γm表示主面,Γp表示主面和从面构成的联合面。对于从面上的测试函数与从面上的物理量之间的积分在从面上进行,主面上的测试函数与主面上的物理量之间的积分在主面上进行。涉及从面上的测试函数与主面上的物理量之间的积分,以及主面上的测试函数与从面上的物理量之间的积分则在主面与从面构成的联合面上进行。
方程式(2)、(9)和(10)构成了周期结构电磁场边值问题的伽略金弱形式。将方程式(2)、(9)和(10)在离散后的计算区域Ωh上进行离散,得到对应的离散Galerkin弱形式如下:
( &dtri; &times; v &RightArrow; h , &mu; r - 1 &dtri; &times; E &RightArrow; h ) &Omega; h - k 0 2 ( v &RightArrow; h , &epsiv; r * E &RightArrow; h ) &Omega; h + k 0 < v &RightArrow; h , j &RightArrow; m h > &Gamma; m h + k 0 < v &RightArrow; h , j &RightArrow; s h > &Gamma; s h + j k 0 &eta; 0 Z s < v &RightArrow; h , n ^ &times; E &RightArrow; h &times; n ^ > &Gamma; SIBC h = 0 - - - ( 11 )
k 0 2 < u &RightArrow; s h , a - 1 j &RightArrow; s h > &Gamma; s h + k 0 < u &RightArrow; s h , e &RightArrow; s h > &Gamma; s h + k 0 < &dtri; &Gamma; &times; u &RightArrow; s h , b a - 1 &dtri; &Gamma; &times; e &RightArrow; s h > &Gamma; s h + e - &gamma;L k 0 2 < u &RightArrow; s h , a - 1 j &RightArrow; m h > &Gamma; p h - e - &gamma;L k 0 < u &RightArrow; s h , e &RightArrow; m h > &Gamma; p h - e - &gamma;L k 0 < &dtri; &Gamma; &times; u &RightArrow; s h , b a - 1 &dtri; &Gamma; &times; e &RightArrow; m h > &Gamma; p h = 0 - - - ( 12 )
k 0 2 < u &RightArrow; m h , a - 1 j &RightArrow; m h > &Gamma; m h + k 0 < u &RightArrow; m h , e &RightArrow; m h > &Gamma; m h + k 0 < &dtri; &Gamma; &times; u &RightArrow; m h , b a - 1 &dtri; &Gamma; &times; e &RightArrow; m h > &Gamma; m h + e &gamma;L k 0 2 < u &RightArrow; m h , a - 1 j &RightArrow; s h > &Gamma; p h - e &gamma;L k 0 < u &RightArrow; m h , e &RightArrow; s h > &Gamma; p h - e &gamma;L k 0 < &dtri; &Gamma; &times; u &RightArrow; m h , b a - 1 &dtri; &Gamma; &times; e &RightArrow; s h > &Gamma; p h = 0 - - - ( 13 )
形式上,方程式(2)、(9)和(10)对应的离散Galerkin弱形式与其Galerkin弱形式相似,所有变量通过加上标h来区分。其中,Ωh为由离散四面体单元构成的仿真区域。
Figure BDA00002073601300074
Figure BDA00002073601300076
分别表示网格划分后的主面,从面和联合面。由于没有强制要求主面和从面的网格匹配,
Figure BDA00002073601300078
的网格划分情况不完全一致。
Figure BDA00002073601300079
则为
Figure BDA000020736013000710
Figure BDA000020736013000711
构成的联合网格面。主面,从面与联合面在几何上完全一致,仅网格划分情况不同。对公式(11)、(12)和(13)中涉及的从面上的测试函数与主面上的物理量之间的积分,以及主面上的测试函数与从面上的物理量之间的积分,在联合网格面
Figure BDA000020736013000712
上进行。由于联合网格面上的每一个面网格都能通过投影确定其所属主面上的三角形网格和从面上的三角形网格,因此不需要额外定义联合网格面上的测试函数和物理量,而是通过投影用相对应的主面或从面上的基函数与物理量来表示。
在每一个离散的四面体单元与主面与从面上的三角形单元内,各待求物理量可由基函数展开,求得相应展开系数,即可得到整个求解域内的场。为此,将式(11)中各四面体单元内的电场强度矢量用基函数展开,有
E &RightArrow; h = &Sigma; i E ~ i N &RightArrow; i - - - ( 14 )
其中,
Figure BDA000020736013000715
是四面体基函数,
Figure BDA000020736013000716
为对应展开系数。关于基函数的选择与确定是有限元法中的公知过程,这里不再赘述。
由于主面与从面网格不匹配,主面面网格与从面面网格上的场不存在一一对应的关系,因此需要将式(12)和(13)中主面和从面的表面电场强度同时用基函数
Figure BDA000020736013000717
Figure BDA000020736013000718
展开,即
e &RightArrow; m h = &Sigma; i e ~ i m N &RightArrow; i m e &RightArrow; s h = &Sigma; i e ~ i s N &RightArrow; i s - - - ( 15 )
其中,
Figure BDA00002073601300082
Figure BDA00002073601300083
是对应基函数的展开系数。
Figure BDA00002073601300084
通常取为切向连续的二阶叠层型基函数,
同理,将式(12)和(13)中主面和从面上的表面电流密度同时用基函数
Figure BDA00002073601300086
Figure BDA00002073601300087
展开,即
j &RightArrow; m h = &Sigma; i j ~ i m M &RightArrow; i m j &RightArrow; s h = &Sigma; i j ~ i s M &RightArrow; i s - - - ( 16 )
其中,
Figure BDA00002073601300089
Figure BDA000020736013000810
是对应基函数的展开系数。
Figure BDA000020736013000811
Figure BDA000020736013000812
通常取为切向不连续的二阶叠层型基函数,
显然,由于四面体基函数
Figure BDA000020736013000813
主面基函数
Figure BDA000020736013000814
Figure BDA000020736013000815
以及从面基函数
Figure BDA000020736013000816
Figure BDA000020736013000817
是事先确定的,只需要求解得到对应的展开系数,即可由公式(14)、(15)和(16)得到每个网格单元内的场,进而得到整个求解域内的场分布。
在式(13)两边同时乘以e-γL  ,在每一个离散的四面体网格内,依次用所选基函数
Figure BDA000020736013000818
Figure BDA000020736013000819
Figure BDA000020736013000820
Figure BDA000020736013000821
Figure BDA000020736013000822
替代式(11)、(12)和(13)中的测试函数
Figure BDA000020736013000823
Figure BDA000020736013000824
Figure BDA000020736013000825
Figure BDA000020736013000826
并将式(14)、(15)和(16)代入到式(11)、(12)和(13),整理得到关于e-γL=e-(α+jβ)L为本征值的大型线性广义本征问题
Ax = e - ( &alpha; + j&beta; ) L Bx = &lambda;Bx - - - ( 17 )
其中,λ称为广义本征问题(17)的特征值,x称为广义本征问题(17)的特征向量。广义本征方程(17)是在引入联合网格的基础上建立的,成功克服了周期结构传统的电磁场本征数值模拟要求主面和从面网格匹配的限制。展开系数向量x包含了周期结构所有四面体网格单元以及主面与从面上各面网格单元内各物理量的展开系数,定义如下
x = E ~ e ~ m j ~ m e ~ s j ~ s T - - - ( 18 )
其中,
Figure BDA000020736013000830
为所有四面体网格单元内各基函数展开系数构成的行向量,
Figure BDA000020736013000831
Figure BDA000020736013000832
分别是主面上各三角形单元内表面电场强度和表面电流密度对应展开系数构成的行向量。
Figure BDA000020736013000833
Figure BDA000020736013000834
分别是从面上各三角形单元内表面电场强度和表面电流密度对应展开系数构成的行向量。
式(17)中,A和B矩阵描述为
A = K II + C II K IM 0 K IS 0 K MI K MM T MM ej 0 0 0 0 0 - T MS je - S MS je T MS jj K SI 0 0 K SS T SS ej 0 0 0 T SS je + S SS je T SS jj - - - ( 19 )
B = 0 0 0 0 0 0 0 0 0 0 0 T MM je + S MM je T MM jj 0 0 0 0 0 0 0 0 - T SM je - S SM je T SM jj 0 0 - - - ( 20 )
式(19)和式(20)中,分块矩阵KII,KIM,KMI,KSI,KSS,KIS可统一描述为
( K XY ) ij = ( &dtri; &times; N &RightArrow; i x , &mu; r - 1 &dtri; &times; N &RightArrow; j y ) &Omega; h - k 0 2 ( N &RightArrow; i x , &epsiv; r * N &RightArrow; j y ) &Omega; h - - - ( 21 )
其中,i与j分别表示分块矩阵的行标与列标。下标X,Y∈{I,M,S},基函数上标x,y∈{i,m,s}。I与i均表示基函数不在周期边界的主面与从面上,而在周期结构的内部。M与m均表示基函数在周期边界主面上。S与s则均表示基函数在周期边界从面上。在后续所有表达式中,均遵循这些规定,不再申明。
式(19)和式(20)中,分块矩阵CII是满足阻抗边界条件的三角形面元上各基函数之间的联系矩阵,定义为
( C II ) ij j k 0 &eta; 0 Z s < N &RightArrow; i i , N &RightArrow; j i > &Gamma; SIBC h - - - ( 22 )
式(19)和式(20)中,分块矩阵
Figure BDA00002073601300096
Figure BDA00002073601300097
分别描述了主面与从面上展开电场强度矢量
Figure BDA00002073601300098
的基函数
Figure BDA00002073601300099
与展开表面电流密度矢量
Figure BDA000020736013000910
的基函数
Figure BDA000020736013000911
在相应主面和从面上的积分特性,统一定义为
( T XX ej ) ij = k 0 &CenterDot; < N &RightArrow; i x , M &RightArrow; j x > &Gamma; X h - - - ( 23 )
其中,下标X∈{M,S},基函数上标x,y∈{m,s}。
式(19)和式(20)中
Figure BDA000020736013000913
Figure BDA000020736013000914
Figure BDA000020736013000915
Figure BDA000020736013000916
分别描述了主面,从面与联合面上,展开表面电流密度矢量的基函数
Figure BDA000020736013000918
与展开电场强度矢量
Figure BDA000020736013000919
的基函数
Figure BDA000020736013000920
在相应面上的积分特性,统一定义为
( T XY je ) ij = k 0 &CenterDot; < M &RightArrow; i x , N &RightArrow; j y > &Gamma; z h - - - ( 24 )
式(19)和式(20)中分块矩阵
Figure BDA00002073601300101
Figure BDA00002073601300102
Figure BDA00002073601300103
Figure BDA00002073601300104
为主面,从面或联合面上展开表面电流密度矢量
Figure BDA00002073601300105
的基函数
Figure BDA00002073601300106
Figure BDA00002073601300107
按照式(25)建立的联系矩阵:
( T XY jj ) ij = k 0 2 &CenterDot; < M &RightArrow; i x , a - 1 M &RightArrow; j y > &Gamma; z h - - - ( 25 )
式(19)和式(20)中,分块矩阵
Figure BDA00002073601300109
Figure BDA000020736013001010
Figure BDA000020736013001011
Figure BDA000020736013001012
是主面,从面或联合面上,展开表面电流密度矢量
Figure BDA000020736013001013
的基函数
Figure BDA000020736013001014
与展开电场强度矢量
Figure BDA000020736013001015
的基函数
Figure BDA000020736013001016
按照式(26)建立的联系矩阵:
( S XY je ) ij = k 0 &CenterDot; < &dtri; &Gamma; &times; M &RightArrow; i x , b a - 1 &dtri; &Gamma; &times; N &RightArrow; j y > &Gamma; z h - - - ( 26 )
式(24)、(25)和(26)中,下标X,Y∈{M,S},基函数上标x,y∈{m,s}。当X=Y=M时,
Figure BDA000020736013001018
的下标z=m,表示积分是在主面网格面上进行。当X=Y=S时,z=s,积分在从面网格面上进行。其他情况下,z=p,表示积分是在联合网格面上进行。
F.给定一个频率,求解步骤E所建立的代数本征方程,获得与给定频率相对应的周期结构的相位常数、衰减常数和互作用阻抗;
通过步骤E建立的线性广义本征方程(17)的矩阵A和B包含与频率相关的元素,如有限电导率导体的表面阻抗Zs,所以采用指定频率的本征分析方法。首先给定一个频率,确定大型稀疏矩阵A和B,求解对应矩阵本征方程(17),得到本征方程(17)的特征值λ与特征向量x。根据λ=e-(α+jβ)L,即可得到与给定频率相对应的相位常数β和衰减常数α,即有
&beta; = - 1 L Im [ ln ( &lambda; ) ] - - - ( 27 )
&alpha; = - 1 L Re [ ln ( &lambda; ) ] - - - ( 28 )
式(28)(27)中,Re(·)与Im(·)分别表示取实部和虚部运算。
根据得到的本征方程(17)的特征向量x,即电场展开系数x,结合基函数,由公式(14)、(15)和(16)可以得到求解域内的场分布,进而可以由互作用阻抗的定义得到与指定频率f相对应的互作用阻抗。由电场分布获得互作用阻抗的过程为本领域的公知过程,因此不再详细描述。
G.给定不同的频率,重复步骤F,获得周期结构的高频特性。
指定不同的工作频率,重复步骤F,得到不同频率对应的相位常数β,衰减常数α和互作用阻抗,即可得到周期结构的色散特性、衰减特性与互作用阻抗特性。因为色散特性,衰减特性与互作用阻抗特性是通过不同频率对应的相位常数β,衰减常数α和互作用阻抗的变化趋势来描述的,如图4,图5与图6所示。
周期结构的色散特性、衰减特性与互作用阻抗特性分别描述了电磁波在高频结构内传播时的色散,衰减以及与电子注的互作用强弱的特性,统称为周期结构的高频特性。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。

Claims (4)

1.基于非匹配网格的周期结构的三维电磁场仿真模拟方法,其特征在于,包括以下步骤:A.选取特定的具有周期特性的微波管高频电路;
B.从步骤A中选取的高频电路中截取一个周期长度的结构进行建模,建立该周期长度的高频结构对应的几何结构模型;
C.根据几何结构模型的周期性确定主面和从面,并对所建几何结构模型进行网格划分,将连续的几何结构空间转化为离散的网格空间;
D.根据仿真区域主面面网格与从面面网格生成联合面网格;
E.利用有限元法,将麦克斯韦方程组等效的电磁场边值问题在步骤C建立的网格空间进行离散,结合步骤D建立的联合网格,建立考虑介质损耗和有限电导率导体损耗周期结构高频电路的代数本征方程;
F.给定一个频率,求解步骤E所建立的代数本征方程,获得与给定频率相对应的周期结构的相位常数、衰减常数和互作用阻抗;
G.给定不同的频率,重复步骤F,获得周期结构的高频特性。
2.根据权利要求1所述的基于非匹配网格的周期结构的三维电磁场仿真模拟方法,其特征在于,步骤C的具体过程为:根据几何结构模型的周期性,将沿着周期结构周期性的方向上仿真区域的初始端面定义为主面,沿着周期结构周期性的方向上仿真区域的最后端面定义为从面;主面与从面的距离为一个空间周期;主面上的场与从面上的场服从周期边界条件;然后,采用四面体网格剖分仿真区域,剖分后的仿真区域被人为分割为多个三维四面体网格,从而将连续的几何结构空间转化为离散的网格空间;此时,主面和从面由一系列离散的面网格单元组成;这里定义附属于主面上的网格单元为主面面网格,附属于从面上的网格单元为从面面网格。
3.根据权利要求1所述的基于非匹配网格的周期结构的三维电磁场仿真模拟方法,其特征在于,步骤D中联合面网格的生成方法包括如下步骤:
1)将主面面网格和从面面网格进行重叠;
2)将主面和从面上的面网格重叠后形成的各个不规则多边形内部取点,并将该点和相应不规则多边形的各个顶点连接,将各个不规则多边形重新剖分成多个新的三角形,从而形成主面和从面的联合面网格。
4.根据权利要求1所述的基于非匹配网格的周期结构的三维电磁场仿真模拟方法,其特征在于,步骤E中建立考虑介质损耗和有限电导率导体损耗周期结构高频电路的代数本征方程具体如式(1)所示:
&dtri; &times; &mu; r - 1 &dtri; &times; E &RightArrow; - k 0 2 &epsiv; r * E &RightArrow; = 0 in&Omega; Z s n ^ &times; &mu; r - 1 &dtri; &times; E &RightArrow; = j k 0 &eta; 0 n ^ &times; ( E &RightArrow; &times; n ^ ) on &Gamma; SIBC n ^ &times; E &RightArrow; s &times; n ^ = n ^ &times; E ^ m &times; n ^ e - &gamma;L on &Gamma; PBC n ^ &times; &mu; r - 1 &dtri; &times; E &RightArrow; s = - ( n ^ &times; &mu; r - 1 &dtri; &times; E &RightArrow; m ) e - &gamma;L on &Gamma; PBC - - - ( 1 )
式(1)中第一个式子为频域矢量波动方程,是周期结构有限元仿真中的主方程;其中,Ω为周期结构的仿真区域空间范围,即为公式(1)的求解域。
Figure FDA00002073601200022
是矢性偏微分算子符号,μr为求解域Ω中介质的相对磁导率,
Figure FDA00002073601200023
为求解域Ω的电场强度矢量,k0为自由空间波数,εr为求解域Ω中介质的相对介电常数。在考虑介质损耗时,εr=ε′r(1-jtanδ)为复数,ε′r为εr的实部,通常为不随频率变化的常数。j为虚数单位符号。-jε′rtanδ为εr的虚部,tanδ是用以描述介质损耗的损耗角正切,通常随频率升高而增加。
CN201210312359.3A 2012-08-29 2012-08-29 基于非匹配网格的周期结构的三维电磁场仿真模拟方法 Expired - Fee Related CN102930071B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210312359.3A CN102930071B (zh) 2012-08-29 2012-08-29 基于非匹配网格的周期结构的三维电磁场仿真模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210312359.3A CN102930071B (zh) 2012-08-29 2012-08-29 基于非匹配网格的周期结构的三维电磁场仿真模拟方法

Publications (2)

Publication Number Publication Date
CN102930071A true CN102930071A (zh) 2013-02-13
CN102930071B CN102930071B (zh) 2015-06-17

Family

ID=47644868

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210312359.3A Expired - Fee Related CN102930071B (zh) 2012-08-29 2012-08-29 基于非匹配网格的周期结构的三维电磁场仿真模拟方法

Country Status (1)

Country Link
CN (1) CN102930071B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412988A (zh) * 2013-08-01 2013-11-27 电子科技大学 基于相移降阶模型周期结构的三维电磁场仿真模拟方法
CN103412984A (zh) * 2013-07-22 2013-11-27 西安空间无线电技术研究所 一种规则波导端口微波部件电磁仿真的功率馈电方法
CN103412989A (zh) * 2013-08-01 2013-11-27 电子科技大学 基于参数化降阶模型周期结构的三维电磁场仿真模拟方法
CN105677981A (zh) * 2016-01-08 2016-06-15 西北工业大学 大规模结构有限元模型的多位移边值约束工况处理方法
CN107330162A (zh) * 2017-06-13 2017-11-07 电子科技大学 针对非曲线边界二维模型的网格划分及信息获取方法
CN107515982A (zh) * 2017-08-22 2017-12-26 电子科技大学 一种三维力学有限元模态分析中的接触分析方法
CN109492284A (zh) * 2018-10-30 2019-03-19 电子科技大学 一种波导端口共形卷积完美匹配层吸收边界算法
CN110058315A (zh) * 2019-05-29 2019-07-26 中南大学 一种三维各向异性射频大地电磁自适应有限元正演方法
CN111931458A (zh) * 2020-10-09 2020-11-13 北京智芯仿真科技有限公司 基于混合阶有限元的三维集成电路电磁场计算方法及装置
CN112131774A (zh) * 2020-11-24 2020-12-25 北京智芯仿真科技有限公司 用于集成电路三棱柱网格剖分的混合阶有限元方法及装置
CN114692450A (zh) * 2022-03-23 2022-07-01 电子科技大学 一种基于泰勒展开的表面阻抗边界条件方法
CN117452081A (zh) * 2023-12-26 2024-01-26 国网天津市电力公司营销服务中心 一种电磁干扰计算方法、装置及存储介质、电子终端

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070208547A1 (en) * 2003-09-05 2007-09-06 Graglia Roberto D Numerical Modeling Process and System of Singular Vector Physical Quantities and Corresponding Software Product
CN101770542A (zh) * 2010-02-25 2010-07-07 中国科学院上海光学精密机械研究所 集群计算机模拟电磁波传播的方法
CN102054094A (zh) * 2010-11-24 2011-05-11 南京理工大学 一种平面微带电路的快速多层方向性仿真方法
CN102103650A (zh) * 2011-03-25 2011-06-22 江苏南大先腾信息产业有限公司 输电线路三维仿真沿布搭建方法
CN102207987A (zh) * 2011-05-31 2011-10-05 中国航天标准化研究所 基于OpenCL的GPU加速三维时域有限差分电磁场仿真的方法
US20110251832A1 (en) * 2010-04-13 2011-10-13 CST-Computer Simulation Technology AG Method, device and computer program product for determining an electromagnetic near-field of a field excitation source of an electrical system
CN102592057A (zh) * 2012-01-17 2012-07-18 电子科技大学 周期结构指定频率的本征分析方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070208547A1 (en) * 2003-09-05 2007-09-06 Graglia Roberto D Numerical Modeling Process and System of Singular Vector Physical Quantities and Corresponding Software Product
CN101770542A (zh) * 2010-02-25 2010-07-07 中国科学院上海光学精密机械研究所 集群计算机模拟电磁波传播的方法
US20110251832A1 (en) * 2010-04-13 2011-10-13 CST-Computer Simulation Technology AG Method, device and computer program product for determining an electromagnetic near-field of a field excitation source of an electrical system
CN102054094A (zh) * 2010-11-24 2011-05-11 南京理工大学 一种平面微带电路的快速多层方向性仿真方法
CN102103650A (zh) * 2011-03-25 2011-06-22 江苏南大先腾信息产业有限公司 输电线路三维仿真沿布搭建方法
CN102207987A (zh) * 2011-05-31 2011-10-05 中国航天标准化研究所 基于OpenCL的GPU加速三维时域有限差分电磁场仿真的方法
CN102592057A (zh) * 2012-01-17 2012-07-18 电子科技大学 周期结构指定频率的本征分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
叶珍: "微波管输入输出窗及任意结构腔体有限元理论与CAD技术", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412984A (zh) * 2013-07-22 2013-11-27 西安空间无线电技术研究所 一种规则波导端口微波部件电磁仿真的功率馈电方法
CN103412984B (zh) * 2013-07-22 2014-10-08 西安空间无线电技术研究所 一种规则波导端口微波部件电磁仿真的功率馈电方法
CN103412989A (zh) * 2013-08-01 2013-11-27 电子科技大学 基于参数化降阶模型周期结构的三维电磁场仿真模拟方法
CN103412988B (zh) * 2013-08-01 2016-07-06 电子科技大学 基于相移降阶模型周期结构的三维电磁场仿真模拟方法
CN103412989B (zh) * 2013-08-01 2016-11-16 电子科技大学 基于参数化降阶模型周期结构的三维电磁场仿真模拟方法
CN103412988A (zh) * 2013-08-01 2013-11-27 电子科技大学 基于相移降阶模型周期结构的三维电磁场仿真模拟方法
CN105677981A (zh) * 2016-01-08 2016-06-15 西北工业大学 大规模结构有限元模型的多位移边值约束工况处理方法
CN105677981B (zh) * 2016-01-08 2019-02-22 西北工业大学 大规模结构有限元模型的多位移边值约束工况处理方法
CN107330162A (zh) * 2017-06-13 2017-11-07 电子科技大学 针对非曲线边界二维模型的网格划分及信息获取方法
CN107515982B (zh) * 2017-08-22 2020-08-11 电子科技大学 一种三维力学有限元模态分析中的接触分析方法
CN107515982A (zh) * 2017-08-22 2017-12-26 电子科技大学 一种三维力学有限元模态分析中的接触分析方法
CN109492284A (zh) * 2018-10-30 2019-03-19 电子科技大学 一种波导端口共形卷积完美匹配层吸收边界算法
CN109492284B (zh) * 2018-10-30 2022-05-03 电子科技大学 一种波导端口共形卷积完美匹配层吸收边界算法
CN110058315B (zh) * 2019-05-29 2020-04-14 中南大学 一种三维各向异性射频大地电磁自适应有限元正演方法
CN110058315A (zh) * 2019-05-29 2019-07-26 中南大学 一种三维各向异性射频大地电磁自适应有限元正演方法
CN111931458A (zh) * 2020-10-09 2020-11-13 北京智芯仿真科技有限公司 基于混合阶有限元的三维集成电路电磁场计算方法及装置
CN111931458B (zh) * 2020-10-09 2021-01-15 北京智芯仿真科技有限公司 基于混合阶有限元的三维集成电路电磁场计算方法及装置
CN112131774A (zh) * 2020-11-24 2020-12-25 北京智芯仿真科技有限公司 用于集成电路三棱柱网格剖分的混合阶有限元方法及装置
CN112131774B (zh) * 2020-11-24 2021-02-26 北京智芯仿真科技有限公司 用于集成电路三棱柱网格剖分的混合阶有限元方法及装置
CN114692450A (zh) * 2022-03-23 2022-07-01 电子科技大学 一种基于泰勒展开的表面阻抗边界条件方法
CN114692450B (zh) * 2022-03-23 2023-07-07 电子科技大学 一种基于泰勒展开的表面阻抗边界条件方法
CN117452081A (zh) * 2023-12-26 2024-01-26 国网天津市电力公司营销服务中心 一种电磁干扰计算方法、装置及存储介质、电子终端

Also Published As

Publication number Publication date
CN102930071B (zh) 2015-06-17

Similar Documents

Publication Publication Date Title
CN102930071B (zh) 基于非匹配网格的周期结构的三维电磁场仿真模拟方法
CN103412989B (zh) 基于参数化降阶模型周期结构的三维电磁场仿真模拟方法
CN102156764B (zh) 一种分析天线辐射和电磁散射的多分辨预条件方法
CN103646152B (zh) 一种基于矩阵指数的电力系统电磁暂态仿真方法
CN102081690B (zh) 复杂电路的矩阵分解结合新奇异值分解方法
CN102609598B (zh) 一种对大型电力系统进行电磁暂态仿真的方法
CN101944145B (zh) 可去除伪直流模式的微波管高频电路有限元仿真方法
CN109190169B (zh) 一种三维时域电磁学杂交时域间断伽辽金数值方法
CN102592057B (zh) 周期结构指定频率的本征分析方法
CN103412988B (zh) 基于相移降阶模型周期结构的三维电磁场仿真模拟方法
CN103425816B (zh) 快速获取金属旋转对称体电磁散射特性的矩阵抽取方法
CN102298658B (zh) 一种行波管注波互作用的模拟方法
CN104794289A (zh) 一种扩展直角坐标系下完全匹配吸收边界的实现方法
CN104679920B (zh) 波导器件微波气体放电数值模拟方法
CN104809343A (zh) 一种等离子体中使用电流密度卷积完全匹配层的实现方法
CN105160115A (zh) 基于近似与灵敏度分析的反射面天线机电集成优化设计方法
Bakr et al. Self‐adjoint S‐parameter sensitivities for lossless homogeneous TLM problems
CN105808504A (zh) 一种等离子体中使用辅助微分方程完全匹配层的实现方法
CN106777472A (zh) 基于拉盖尔多项式的减少分裂误差的完全匹配层实现方法
CN105760597A (zh) 基于DG算法的二维色散介质Crank-Nicolson完全匹配层实现算法
CN105205299B (zh) 电大目标电磁散射特性快速降维分析方法
CN104778293A (zh) 非均匀介质目标电磁散射的体积分Nystrom分析方法
CN105760596A (zh) 基于辅助微分方程的二维真空Crank-Nicolson完全匹配层实现算法
Jolly et al. Novel integrated design framework for radio frequency quadrupoles
CN106951607B (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150617

Termination date: 20160829