第五章 有限元法-5-时谐场波导本征值问题

1、5.6 有限元的应用问题有限元的应用问题 确定性问题与本征模问题确定性问题与本征模问题 p简而言之，有限元方法就是离散泛函数变分数学表达形式的离散化方简而言之，有限元方法就是离散泛函数变分数学表达形式的离散化方 法。法。 p有限元法在电磁学的演进主要围绕两条主线：有限元法在电磁学的演进主要围绕两条主线： n一是如何求解一是如何求解本征模问题本征模问题；(粗略地说，就是谐振问题）粗略地说，就是谐振问题） n一条是如何求解一条是如何求解确定性问题确定性问题。(前面提到的都是确定性问题，主要是传输前面提到的都是确定性问题，主要是传输 传播问题）传播问题） n本征模问题可细分为本征模问题可细分为波导本

2、征模波导本征模和和谐振腔本征模谐振腔本征模。这两类进而又可细分。这两类进而又可细分 为空、介质填充、介质三种。为空、介质填充、介质三种。 n确定性问题可分为确定性问题可分为闭域传输问题闭域传输问题和和开域散射、辐射问题开域散射、辐射问题。 p确定性确定性问题中，控制微分方程和边界条件其中之问题中，控制微分方程和边界条件其中之 一、或它们两者，都是一、或它们两者，都是非齐次非齐次的。的。 p不同于确定性问题，不同于确定性问题，本征值本征值问题的控制微分方程问题的控制微分方程 和边界条件均是和边界条件均是齐次齐次的。的。 n从物理学的观点看，在本征值问题中从物理学的观点看，在本征值问题中不存在任何

3、形式不存在任何形式 的源或激励的源或激励。 有源，由激励产生 无源，是固有的 p本征值问题用有限元方法处理导出的方程组具有本征值问题用有限元方法处理导出的方程组具有 下列下列广义本征值方程广义本征值方程的形式的形式 p在电磁学中，经常出现的在电磁学中，经常出现的本征值问题包括本征值问题包括腔体谐振腔体谐振问题以问题以 及波在封闭和开放结构及波在封闭和开放结构(波导波导)中的中的传播模式传播模式问题。问题。 n封闭和开放结构封闭和开放结构(波导波导)包括金属波导、开放和屏蔽微带传输线、包括金属波导、开放和屏蔽微带传输线、 光波导以及光纤等。光波导以及光纤等。 n在这些问题中，须确定在这些问题中，

4、须确定对应于本征值的谐振频率或传播常数对应于本征值的谐振频率或传播常数，和，和 确定确定对应于本征向量的有关谐振模式或传输模式对应于本征向量的有关谐振模式或传输模式。 n可先考虑填充均匀各向同性媒质的金属波导可先考虑填充均匀各向同性媒质的金属波导，再考虑填充不均匀，再考虑填充不均匀 媒质的金属波导。然后考虑填充单轴各向异性媒质和双轴各向异媒质的金属波导。然后考虑填充单轴各向异性媒质和双轴各向异 性媒质的波导问题。之后，再考虑能够处理一般各向异性波导的性媒质的波导问题。之后，再考虑能够处理一般各向异性波导的 矢量场公式，并讨论开放导波结构的分析。矢量场公式，并讨论开放导波结构的分析。 p介质填充

5、波导本征模介质填充波导本征模 n广义上说，波导是指一切用来引导电磁波的传输线。广义上说，波导是指一切用来引导电磁波的传输线。 n对于规则波导，由波导理论可知，其内存在一系列可单独存在的对于规则波导，由波导理论可知，其内存在一系列可单独存在的 模式，这些模式在横截面上的场分布和纵向上的传播常数都不同，模式，这些模式在横截面上的场分布和纵向上的传播常数都不同， 如何确定它们就是本征值问题。如何确定它们就是本征值问题。 n由于有限元方法能将问题转化成数学上标准的矩阵本征值或广义由于有限元方法能将问题转化成数学上标准的矩阵本征值或广义 本征值问题，因而有限元法比矩量法或时域有限差分法更适于解本征值问题

6、，因而有限元法比矩量法或时域有限差分法更适于解 决本征值问题。决本征值问题。 n它们属于时谐场问题。它们属于时谐场问题。 5.6.1 时谐电磁场中的有限元法时谐电磁场中的有限元法 p对于任一时谐电磁场，可应用不同的物理量进行对于任一时谐电磁场，可应用不同的物理量进行 研究。研究。 n例如，既可采用场量例如，既可采用场量 ，也可采用时变场中的位函，也可采用时变场中的位函 数，即动态数，即动态 来着手分析。来着手分析。 p关于分析对象的选择，完全取决于实际问题的性关于分析对象的选择，完全取决于实际问题的性 质，取决于能否简明地给出问题的定解条件。质，取决于能否简明地给出问题的定解条件。 H、E 、

7、 A p在时谐电磁场问题的研究中，为了在时谐电磁场问题的研究中，为了简化简化分析，通分析，通 常引入以下假设：常引入以下假设： n媒质系线性且各向同性，并忽略时变场角频率媒质系线性且各向同性，并忽略时变场角频率w w对于媒对于媒 质的影响；质的影响； n忽略铁磁材料的磁滞效应与导体电阻率的温度效应；忽略铁磁材料的磁滞效应与导体电阻率的温度效应； n实际三维场问题被理想化为二维平行平面场或轴对称实际三维场问题被理想化为二维平行平面场或轴对称 场。场。 5.6.2 波导场问题的有限元方程波导场问题的有限元方程 p当分析波导中电磁波传播问题时，常进一步假设：当分析波导中电磁波传播问题时，常进一步假设

8、： n波导壁由完纯导体（波导壁由完纯导体（g g）构成；）构成； n波导中无自由电荷和传导电流（波导中无自由电荷和传导电流（r r=0，J=0），也就），也就 是说，波导是远离激励源的；是说，波导是远离激励源的； n波导工作在匹配状态，具有均匀的截面。因此，在分波导工作在匹配状态，具有均匀的截面。因此，在分 析中只考虑向前传播的入射波，无反射波。这样，由析中只考虑向前传播的入射波，无反射波。这样，由 于激励源激发的不同状况，波导中传播的电磁波可分于激励源激发的不同状况，波导中传播的电磁波可分 为横电波（为横电波（TE波）或横磁波（波）或横磁波（TM波）两种类型。波）两种类型。 p电磁场理论中有

9、关波导问题的分析表明，波导中电磁场理论中有关波导问题的分析表明，波导中 场量在随时间作正弦变化的同时，也在波导空间场量在随时间作正弦变化的同时，也在波导空间 中沿波导方向（设为中沿波导方向（设为z方向）呈行波特征的变化，方向）呈行波特征的变化， 而而一旦得知场量的纵向分量，便可求出相应的横一旦得知场量的纵向分量，便可求出相应的横 向分量向分量。 p因此，对于因此，对于TE波波，应取，应取 为着手分析的对象，根为着手分析的对象，根 据向量波动方程，可知应满足如下的标量波动方据向量波动方程，可知应满足如下的标量波动方 程（令程（令k2 2=w=w2 2meme）：）： z H z E 同理，同理，

10、TM波波则归结为求解其纵向分量则归结为求解其纵向分量 所对应的波动方程所对应的波动方程 因此，若以因此，若以f f标记相应的纵向分量标记相应的纵向分量Hz或或Ez，则波导场的分析，则波导场的分析 可归结为如下定义于波导横截面（可归结为如下定义于波导横截面（x，y）平面内的二维标）平面内的二维标 量波动方程（亥姆霍兹方程）的解答，即量波动方程（亥姆霍兹方程）的解答，即 式中式中l l2 = k2 kz2；kz描述场量沿描述场量沿z方向每单位长度中相位的方向每单位长度中相位的 变化，被称为相位系数（亦称波数）。应注意，上式中变化，被称为相位系数（亦称波数）。应注意，上式中f f （Hz或或Ez）以

11、及相应的各横向分量均仅是）以及相应的各横向分量均仅是x，y的函数。的函数。 亥姆霍兹方程是椭圆型方程亥姆霍兹方程是椭圆型方程 对于对于TE波波，应令，应令f f = Hz。由于前述假设。由于前述假设(1)，故导壁内部磁场，故导壁内部磁场 为零，此时在导壁表面将呈现面电流为零，此时在导壁表面将呈现面电流JS，所以作为定解条件，所以作为定解条件， 即有即有 ；0= n f 对于对于TM波波，f f = Ez。同理，因导壁内没有电场，故电场。同理，因导壁内没有电场，故电场 强度的切向分量应等于零，即有定解条件强度的切向分量应等于零，即有定解条件f f = 0。 由此可知，描述波导场的定解问题为 式中

12、，边界式中，边界L即为波导壁。根据变分原理，容易证明，即为波导壁。根据变分原理，容易证明， 与上述边值问题的泛定方程（亥姆霍兹方程）对应的泛与上述边值问题的泛定方程（亥姆霍兹方程）对应的泛 函为函为 因此，与因此，与TE波波的波导场定解问题（的波导场定解问题（5-135a）、（）、（5-135b） 等价的无条件变分问题为等价的无条件变分问题为 而与而与TM波波的波导场定解问题（的波导场定解问题（5-135a）、（）、（5-135c）等）等 价的则为条件变分问题价的则为条件变分问题（5-137）、（）、（5-135c）。 式中，对应于所选取的有限单元式中，对应于所选取的有限单元e，其各个单元矩阵

13、的，其各个单元矩阵的 元素分别为元素分别为 对泛函（对泛函（5-137）取极值，并经有限元的离散化处理，可以）取极值，并经有限元的离散化处理，可以 导得如下的有限元方程导得如下的有限元方程 和和 式（式（5-138）即是所谓）即是所谓广义代数特征值问题广义代数特征值问题. 广义代数特征值问题（广义代数特征值问题（5-138）将给出）将给出波导场定解问题波导场定解问题 （即亥姆霍兹方程的特征值问题）（即亥姆霍兹方程的特征值问题）式（式（5-135a）、（）、（5- 135b）或式（）或式（5-135a）、（）、（5-135c）的特征值）的特征值l l2和特征和特征 向量向量f f的解答，然后，由

14、此可得给定波导中截止频率：的解答，然后，由此可得给定波导中截止频率： 截止波长：截止波长： 及各种可能波型的场分布图等。及各种可能波型的场分布图等。 5.6.3、有限元的三个方面、有限元的三个方面 p有限元法的演进，主要在三个方面：有限元法的演进，主要在三个方面： n一是泛函变分表达式各种形式的研究；一是泛函变分表达式各种形式的研究； n二是基函数的恰当构造二是基函数的恰当构造（形状和插值参量选取）；（形状和插值参量选取）； n三是广义本征值方程的快速求解。三是广义本征值方程的快速求解。 （1）泛函变分表达式）泛函变分表达式 p一般说来，泛函变分表达式是由描述问题的偏微分方程推导而来。一般说来

15、，泛函变分表达式是由描述问题的偏微分方程推导而来。 p而描述问题的偏微分方程通常有多种形式，这样描述问题的泛函变分而描述问题的偏微分方程通常有多种形式，这样描述问题的泛函变分 表达式一般也就多种多样。表达式一般也就多种多样。 p这些泛函数变分表达式往往有适用范围窄宽之分，效率高低之别。这些泛函数变分表达式往往有适用范围窄宽之分，效率高低之别。 n就介质填充波导本征模问题而论，就有纵向电场和磁场共同构建的，有就介质填充波导本征模问题而论，就有纵向电场和磁场共同构建的，有 横向电场和磁场共同构建的，还有矢量位和标量位共同构建的等多种泛横向电场和磁场共同构建的，还有矢量位和标量位共同构建的等多种泛

16、函数变分表达式。函数变分表达式。 n就适用范围和效率综合而论，从矢量波动方程导出的全电场或磁场泛函就适用范围和效率综合而论，从矢量波动方程导出的全电场或磁场泛函 变分表达式较为恰当。这种泛函变分表达式适用于任意填充介质，泛函变分表达式较为恰当。这种泛函变分表达式适用于任意填充介质，泛函 中的未知量只有电场或磁场。中的未知量只有电场或磁场。 （2）基函数的选取）基函数的选取 p与矩量法一样，离散泛函，首先要选取基函数，与矩量法一样，离散泛函，首先要选取基函数， 将未知变量用基函数线性组合的方式表达出来。将未知变量用基函数线性组合的方式表达出来。 p基函数的选取有如矩量法，关键在两点：基函数的选取

17、有如矩量法，关键在两点：形状和形状和 插值参量插值参量。 p当求解域为面，很明显基函数的形状要以当求解域为面，很明显基函数的形状要以三角形三角形 最为灵活方便。最为灵活方便。 （3）插值参量的选取）插值参量的选取 p至于插值参量的选取如何方为恰当，并非那么明显。至于插值参量的选取如何方为恰当，并非那么明显。 n如用三角形三个顶点的矢量电场作为插值参量，则求解结果中可能含有如用三角形三个顶点的矢量电场作为插值参量，则求解结果中可能含有 很多伪解（很多伪解（spurious solutions)。 n其原因是以三角形三个顶点的矢量电场作为插值参量，得到的基函数不其原因是以三角形三个顶点的矢量电场作

18、为插值参量，得到的基函数不 仅保证了相邻单元切向电场连续，同时也仅保证了相邻单元切向电场连续，同时也额外强加了法向电场连续额外强加了法向电场连续，这，这 是不符合物理意义的。是不符合物理意义的。 n解决的办法是：除了纵向电场解决的办法是：除了纵向电场Ez仍选用三顶点处值作为插值参量，横向仍选用三顶点处值作为插值参量，横向 电场分量电场分量Et改用三角形各边中点处的切向电场改用三角形各边中点处的切向电场Eti用为插值参量。这样构用为插值参量。这样构 建的基函数称之为建的基函数称之为边缘元（边缘元（edge-element)基函数。基函数。 n三角形三顶点按逆时针方向分别标为三角形三顶点按逆时针方

19、向分别标为1、2、3，面积有，面积有 33 22 11 1 1 1 2 1 yx yx yx = 棱边元 5.6.4、三维波导不连续性问题、三维波导不连续性问题 p很多微波器件如滤波器、定向耦合器、环行器等很多微波器件如滤波器、定向耦合器、环行器等 的理论分析都可归结为波导不连续性问题的分析。的理论分析都可归结为波导不连续性问题的分析。 p波导不连续性问题可细分为波导不连续性问题可细分为开波导不连续性开波导不连续性和和闭闭 波导不连续性波导不连续性两类，两类， n像光纤中的不连续性就属于开波导不连续性问题，像光纤中的不连续性就属于开波导不连续性问题， n而金属波导中不连续性就属于闭波导的不连续

20、性问题。而金属波导中不连续性就属于闭波导的不连续性问题。 p这两类问题的分析有本质不同，因为这两类问题的分析有本质不同，因为 n开波导开波导不连续性产生的散射场会向无限大空间辐射，其求解域外为无限不连续性产生的散射场会向无限大空间辐射，其求解域外为无限 大；大； n而而闭波导闭波导不连续性产生的散射场将被波导壁封在波导内，其解域只限在不连续性产生的散射场将被波导壁封在波导内，其解域只限在 波导之内。波导之内。 p对于对于开波导开波导不连续性问题，不连续性问题，如辐射现象较弱，可以忽略如辐射现象较弱，可以忽略，则其分析过，则其分析过 程与闭波导不连续性分析并无二致；程与闭波导不连续性分析并无二致

21、；如果辐射现象严重，有限元分析如果辐射现象严重，有限元分析 就较为困难，实际中一般也就不用，转而采用矩量法分析或合元极技就较为困难，实际中一般也就不用，转而采用矩量法分析或合元极技 术。术。 p可见，就三维波导不连续性问题的有限元求解而言，可见，就三维波导不连续性问题的有限元求解而言，闭波导闭波导不连续性不连续性 问题的求解更问题的求解更典型典型。 5.6.5、本征模解法对比、本征模解法对比 p本征模的成熟方法可分频域法和时域法两大类。本征模的成熟方法可分频域法和时域法两大类。 n频域法的基本思想是通过寻找频率域的本征模和本征频域法的基本思想是通过寻找频率域的本征模和本征 值获得特定边值问题的

22、解；值获得特定边值问题的解； p频域法包括：矩量法、有限元法、有限差分法。频域法包括：矩量法、有限元法、有限差分法。 n而时域法则与电磁波的瞬态特性相联系，通过模拟电而时域法则与电磁波的瞬态特性相联系，通过模拟电 磁波与目标作用的瞬态特性，得到时间相关的四维数磁波与目标作用的瞬态特性，得到时间相关的四维数 值解。值解。 p矩量法矩量法，也常称为，也常称为“表面积分技术表面积分技术”。这个方法。这个方法 通常在分析无边界辐射问题时可得到很好结果。通常在分析无边界辐射问题时可得到很好结果。 它也它也擅长于分析纯导体结构和单一介质材料结构擅长于分析纯导体结构和单一介质材料结构 的电磁场问题。的电磁场

23、问题。但这个方法不太适用于多种材料但这个方法不太适用于多种材料 的复合结构分析和复杂结构的计算分析。的复合结构分析和复杂结构的计算分析。 p有限元法有限元法要求将整个被分析区进行网格化，而不要求将整个被分析区进行网格化，而不 像矩量法只网格化物体表面。每个网格可以是不像矩量法只网格化物体表面。每个网格可以是不 同的材料，因此同的材料，因此可以用于多种材料的复合结构分可以用于多种材料的复合结构分 析。但它在分析无边界辐射问题时不如矩量法析。但它在分析无边界辐射问题时不如矩量法。 5.6.7 新方法新方法 p有限元法的主要特点是可以用任意形状的小单元来描述问有限元法的主要特点是可以用任意形状的小单

24、元来描述问 题的复杂的几何形状和复杂的媒质分布。题的复杂的几何形状和复杂的媒质分布。 p除频域应用，还有时域有限元法（除频域应用，还有时域有限元法（TDFEM）和谐波平衡）和谐波平衡 有限元法（有限元法（BH（harmonic balance)-FEM）。）。 n谐波平衡有限元法谐波平衡有限元法以解决姆霍兹方程为重点，因时谐波动方程出以解决姆霍兹方程为重点，因时谐波动方程出 现在开放空间的电磁辐射问题中。现在开放空间的电磁辐射问题中。 n时域有限元时域有限元的矩阵方程可由最常用的的矩阵方程可由最常用的Newmark-b b法来解。法来解。 n时域有限差分法中的时域有限差分法中的吸收边界条件和理

25、想匹配层（吸收边界条件和理想匹配层（PML）条件同条件同 样被用于有限元的分析中。样被用于有限元的分析中。 1、时域有限元法和谐波平衡有限元法、时域有限元法和谐波平衡有限元法 2、节点元法和棱边元法、节点元法和棱边元法 p早期有限元方法，用插值节点数值而获得的节点基单元来早期有限元方法，用插值节点数值而获得的节点基单元来 表示矢量电场或磁场，但会遇到几个严重问题。表示矢量电场或磁场，但会遇到几个严重问题。 n首先，可能会有非物理的赝解出现，这通常是由于未强加散度条首先，可能会有非物理的赝解出现，这通常是由于未强加散度条 件引起的；件引起的； n其次，在材料界面和导体表面强加边界条件不方便；其次，在材料界