2009Part2浙江大学电气学院电磁场数值计算课件.ppt

December 9, 2009 2）伽辽金有限元法（Galerkin Finite element method) B .伽辽金有限元法（Galerkin Finite element method) 2.有限元法(2D) （3）总体合成Assembly 把节点编号由i、j、m扩展成整体坐标 3）第一类边界条件（essential boundary condition)的处理 A .迭代法求解代数方程组 2.有限元法(2D) 强加边界条件，无需通过有限元方程求解，相反，正是给定这些边界节点上的位函数值的基础上而推求其余节点上的位函数值 理论上可以不处理。实践中处理后收敛速度较快。 3）第一类边界条件（essential boundary condition)的处理 B .消去法 2.有限元法(2D) 强加边界条件，无需通过有限元方程求解，相反，正是给定这些边界节点上的位函数值的基础上而推求其余节点上的位函数值 必须处理，同学自学。 2）伽辽金有限元法（Galerkin Finite element method) B .伽辽金有限元法（Galerkin Finite element method) 2.有限元法(2D) （1）加权余数表示 对于以二维泊松方程描述的低频电磁场的混合边值问题 2）伽辽金有限元法（Galerkin Finite element method) B .伽辽金有限元法（Galerkin Finite element method) 2.有限元法(2D) （1）加权余数表示 2）伽辽金有限元法（Galerkin Finite element method) B .伽辽金有限元法（Galerkin Finite element method) 2.有限元法(2D) （1）加权余数表示 对于第一类边界项，在离散后通过修正（改）方程组，使 在边界?1上精确地等于u0。故无须列入。由此有 2）伽辽金有限元法（Galerkin Finite element method) B .伽辽金有限元法（Galerkin Finite element method) 2.有限元法(2D) （1）加权余数表示 该式称为加权余量的强式（strong form)表述。 2）伽辽金有限元法（Galerkin Finite element method) B .伽辽金有限元法（Galerkin Finite element method) 2.有限元法(2D) （1）加权余数表示 由此，可以通过代入 导出有限元方程 2）伽辽金有限元法（Galerkin Finite element method) B .伽辽金有限元法（Galerkin Finite element method) 2.有限元法(2D) （1）加权余数表示 此有限元方程是基于强式的余量形式导出的，故称为强式有限元方程。 2）伽辽金有限元法（Galerkin Finite element method) B .伽辽金有限元法（Galerkin Finite element method) 2.有限元法(2D) （1）加权余数表示 强式有限元方程的缺点： （i）要求形状函数二阶可导 （ii)相反，对权函数的可导性并无要求 2）伽辽金有限元法（Galerkin Finite element method) B .伽辽金有限元法（Galerkin Finite element method) 2.有限元法(2D) （1）加权余数表示 应用格林公式（Green’s equation），通过部分积分，将二阶偏导数转化为一阶偏导数 2）伽辽金有限元法（Galerkin Finite element method) B .伽辽金有限元法（Galerkin Finite element method) 2.有限元法(2D) （1）加权余数表示 2）伽辽金有限元法（Galerkin Finite element method) B .伽辽金有限元法（Galerkin Finite element method) 2.有限元法(2D) （1）加权余数表示 2）伽辽金有限元法（Galerkin Finite element method) B .伽辽金有限元法（Galerkin Finite element method) 2.有限元法(2D) （1）加权余数表示 原问题的弱（weak form）表述 2）伽辽金有限元法（Galerkin Finite element method) B .伽辽金有限元法（Galerkin Finite element method) 2.有限元法(2D) （1）加权余数表示 在第一类边界上把权函数选为零，则有 2）伽辽金有