第三章FEM2014讲解.ppt

利用豪斯豪尔法 ：对称三角矩阵 经过多次相似变换 满足指定精度的对角阵 对角阵的每个元素就是其特征值 反变换 求得的一系列特征值中非负的最小非零特征值 就给出相应波导中最低型（主模）的截止波长 截 止 波 长 理论值 相对误差 波型 (cm) 数值解 TE10 4.572 4.512 4.569 4.574 0.066 TE20 2.286 2.193 2.262 2.275 1.05 TE01 2.032 1.847 1.979 2.61 TE11 1.857 1.635 1.790 3.61 TE30 1.524 1.410 1.488 2.36 TE21 1.519 1.250 1.440 5.20 TE31 1.219 1.021 1.141 6.40 TE40 1.143 0.943 1.093 4.37 TE02 1.016 0.922 0.923 9.15 TE41 0.996 0.859 0.916 8.03 TE22 0.928 0.795 0.891 3.99 * * * 强加边界条件的处理 设第m个节点上具有强加边界条件 =1 =0 处理强加边界条件小结 将对角线元素 置1 该行的右端项改为强加电位值 第m行与m列的其他元素全部置零 除第m行外，其他各行右端项为原右端项减去强加电位置 与对应第m列未变换前的系数的乘积 若强加边界条件节点共有N0个，则如法处理N0次 N阶方阵 与强加边界节点电位值及有关量组成的列阵 系数矩阵的存储 利用其对称性 、稀疏性 等带宽存储 变带宽存储 具体作法 将[K]矩阵的下三角部分带内元素按行的顺序， 从各行的第一个非零元素起， 至主对角元素（包括带内的零元素及对角元素）为止的元素 （共有MD+1），依次存入一维数组KL(n)中 ★ ★ ★ 说明KL(n)数组中各元素在原矩阵[K]中的行、列位置 指针数组L(n) ﹢ 存储原[K]阵中主对角元素 在KL(n)中的位置 需要计算的量 [K]的下三角阵各行的半带宽MD 第i行的半带宽MD=i-min(j) 主对角元素在KL(n)的地址JO JO=L(i)=L(i-1)+MD+1 非主对角元素在KL(n)的地址JO JO=L(i)-(i-j) 非主对角元素的地址为该行主对角元素的地址向前推移(i-j)个位置 所有以节点i为顶点的各单元上节点总体编号的最小值 前一行主对角元的地址 行号i MD 在KL(n)中 位置（即L(n)值） 在KL(n) 中位置 第一行 1 第二行 1 3 的JO＝2 第三行 2 6 的JO＝4 的JO＝5 第四行 2 9 的JO＝7 的JO＝8 第五行 2 12 的JO＝10 的JO＝11 第六行 3 16 的JO＝13 的JO＝14 的JO＝15 [K]阵中各元素在一维数组KL(n)中的位置为 KL(1) KL(2) KL(3) KL(4) KL(5) KL(6) KL(7) KL(8) KL(9) KL(10) KL(11) KL(12) KL(13) KL(14) KL(15) KL(16) 高斯消元法求解线性代数方程组 消元 将原方程的系数矩阵化为对角元素为1的上三角阵 回代 从最末一个方程（只包含一个未知数）求出一个未知数，再代入上一个方程，求出另一个未知数 消元公式 第一个方程的系数 及右端项的处理 回代公式 3.10轴对称场的有限元方程 轴对称场的有限元方程 二维直角坐标系 是具有非齐次自然边界条件的边界单元的贡献 对称性条件的利用及处理 左右对称 强加边界条件左右对称 左右反对称 多个对称轴 场域 介质系数 边界条件 关于几条直线对称 3.11电场强度的计算 一单元中E 的数值解 与坐标x、y无关， 仅与单元节点的坐标 及电位有关 单元的节点一般与相邻几个单元相联系 节点场强必须将有关单元对场强综合贡献 分别总和有关单元的 ，然后取算术平均值 分别将有关单元的 加权后总和，以单元面积 为权因子，然后取其平均值 场强的方向 剖分细 对称 分界面处细剖分 高精度： 3.12波导场的有限元解 用有限元法求解的优越 将截面逐步细分，将导致本征值向极限单调减小 可保证有一个更快的收敛速度趋于本征值 对难处理的边界形状更容易处理，而不导致非对称矩阵 波导中的奇点无需特别处理 假设 波导壁为完纯导体 波导内的介质系均匀、线性且各向同性的理想介质 波导中无自由电荷和传导电流 波导工作在匹配状态，具有均匀截面，所以在分析时只考虑入射波，无反射波 波导中传播的电磁波可分为TE波或TM波 求解相应的场纵向分量所描述的定解问题