有限元法(2007-07-06).ppt

有限元和本征函数展开 因此，可以确定常数amn和bmn 为了导出S1上的边界条件，我们对E（x,y,z）取旋度， 有限元和本征函数展开 然后用 叉乘它，得到 将amn和bmn代入上式，得到 在z=z1处 其中 有限元和本征函数展开 类似地，在虚构面S2上的场可展开为 其中 有限元和本征函数展开 由此可导出边界条件 其中P(E)仍由前式给出，z2表示平面S2的位置。注意， 如果波导截面在S2处与S1处的尺寸不同，应使用a和b的 新参数。 虚构面S1与S2处的边界条件导出后，由波导壁、 S1和S2 围成区域中场的边值问题就唯一地确定了。该问题的 泛函由下式给出 在z=z2处 有限元和本征函数展开 将P(E)和Uinc代入上式，我们得到只用E表示的泛函。 然后，使用矢量单元离散上式给出的泛函，即可求得 所需要的数值场解。一旦得到S1和S2处的场，我们可以 确定amn和bmn等常数，从而求出每个反射模和传输模的 幅度，由此可确定反射系数和传输系数。 有限元和本征函数展开 开放区域散射 前述本征函数展开方法的基本思想也可用于开放区域 的散射分析。现以图示二维问题为例说明该方法的具 体应用。 散射体 ? 二维散射问题 有限元和本征函数展开 设?是一个大得足以包括散射体的虚构圆形边界。在? 上及其外边，外行散射场可以展开为 式中， 表示第二类n阶汉克尔函数，a是?的半径。 在展开式中，每个本征函数 都是齐次亥 姆霍兹方程的解，且满足索末菲辐射条件。它们形成 一套完备基函数，能表达任何外行场。 有限元和本征函数展开 由正交关系 得 由上式可得 有限元和本征函数展开 上式可认为是第三类广义边界条件，它可写为 其中 它可直接应用在散射场公式中。对于总场公式，我们 将?sc＝? -? inc代入边界条件得到 有限元和本征函数展开 上式边界条件唯一地确定了?内场的边值问题，它的解 可通过寻找下式给出的泛函的驻点而获得 有限元和本征函数展开 其中，对Ez极化，?＝Ez，u=?r，v=?r；对Ez极化，?＝Hz ，u=?r ，v=?r。该泛函可通过有限元方法离散，结果产 生一个部分满秩、部分稀疏的矩阵方程，由该矩阵方程 可解出?上和?内的场。一旦求出?上的场，利用展开系 数公式，我们可计算出?外的散射场。 用有限元法求解波导问题 如果单元et的i，j，m相应的总体编号分别为p，g，r，则 对其它以总体编号p为顶点的单元也可得出类似的公式，当然求出的局部 系数第二下标将随单元的不同而不同。这样， 上式就可写成下面的形式 用有限元法求解波导问题 式中系数 是由各个单元求出的同样下标的局部B系数组合而成。最后得到 式中 用有限元法求解波导问题 由于K，B两矩阵中各元是根据各个三角元求出的局部K，B系数矩阵中相同下标的元组成，而局部K，B系数矩阵是对称的，所以K，B矩阵也是对称的，也可以证明K，B矩阵是正定阵。 编制用有限元法计算波导问题的程序 输入节点总数为L1，待求场值的节点数为L0，单元总数为JE0，各节点坐标值，各个三角形i，j，m相应的整体编号 用有限元法求解波导问题 对TM波将边界节点置零 NE=1 （NE为某个三角元的编号） 求NE三角形的各个b，c值 当i(NE)=p， j(NE)=q， m(NE)=r时，求NE三角形的面积和九个局部K系数和九个局部B系数 用有限元法求解波导问题 将同样下标的局部K，B系数连加形成总体K，B系数 NEJE0? 调用广义特征值问题各个子程序求各模式波的kkp，u 是 打印输出 否 以高次多项式为插值函数的有限元法 三角形上的三点线性插值是一种简单的插值方法，只是当单元取得足够小时才有较好的计算精度。为使精度提高，可以采用高次多项式作为插值函数，下面就以拉普拉斯方程为例来说明它的求解过程。 · 拉普拉斯方程的等价修正的泛函为 这时试探函数u的选取不受边界条件的限制。 以高次多项式为插值函数的有限元法 当采用有限元法求解时，仍先将讨论的场域剖分成很多三角形单元。但现在括值函数采用高次多项式。因为试探函数u不需要满足边界，所以可以来用完整的N阶多项式比较方便，即令 以高次多项式为插值函数的有限元法 为了求上式中的M个ci常数，就需要在三角形单元上选M个节点，根据每个节点，例如j点的坐标(xj，yj)和电位uj ，代入式后就得到M个方程，根据这M个方程就能解出M个ci值。M个节点的配置一般是在三角形每个边上等间隔配N+1个点(包括三角形的顶点)，其余M -3N个点处于三角形的内部。下图示出了在N=2和N=3时三角形上M个节