有限元法基础通常将有限元法分为两大类变分法和加权余量.doc

PAGE PAGE 13 有限元法基础 通常将有限元法分为两大类：变分法和加权余量法。两种方法的出发点不同，但最后都归结为：①离散化：用若干个子区域（即单元）代替整个连续区域，②算子解析方程，即偏微分方程转化为代数方程组：区域的物理性质可以用节点上有限个自由度来描述，再应用离散系统分析方法将其汇集在一起。 §3-1 算子方程及变分原理 3.1.1 算子的概念 （1）静电场中，泊松方程 可以写为 ，其中称为算子。 （2）稳态磁场中，双旋度方程 其中算子是 （3）时变场中，波动方程 其中算子 3.1.2 泛函 1、泛函的概念 泛函是函数空间H中，函数到数的映像，如 也可以说泛函是函数的函数，函数空间中的某一函数有一个I值与之对应，变量I就是D空间的函数的泛函。 例如 求所表示的曲线长度及所围面积。 曲线长度 曲线所围面积 不同的，有不同的I与之对应，不同的 图3-1 求曲线长度及所围面积 构成了函数空间H。 2、泛函连续 若对于的微小改变，有泛函的微小改变与之对应，就称泛函是连续的。 3、线性泛函 若泛函满足 c为常数 或 则称其为线性泛函。 4、函数的变分 泛函的宗量的变分是的微小增量 5、泛函的变分 对于宗量的变分，泛函的增量为 式中，是对的线性泛函，是的主要部分，称为一阶（或一次）变分 是误差项。 与dy的区别： 当自变量x的增量充分小时，可用dx来表示，dx称为x的微分。相应地，函数y的增量 当充分小时，可用dy来表示，dy称为y的微分，dy 是x的变化引起的微分，是函数增量 的线性主要部分，即可记为 在泛函中，当宗量的增量足够小，即有变分时，泛函的增量 其中，对而言为线性，称为一阶变分。 图3-2 函数的增量 图3-3 泛函的增量 6、泛函的极值 设时泛函取得极值，那么，泛函在极值函数上的变分等于0，即 当泛函是多元函数的泛函，泛函在上有极值时，变分。因此，泛函取得极值的必要条件是使变分。 3.1.3 算子（微分、积分、矩阵方程）方程的变分原理 各种类型电磁场的微分方程都可对应于D空间中的算子方程 它可以转化为与之等价的变分问题，即泛函求极值问题。 定理：若L为正算子，而在D上有解，则此解必然使泛函 取极小值。反之，在D上使泛函I取得极小值的函数，必是方程的解。（证明略，参见颜威利《电气工程电磁场数值分析》P24-25. 也就是说，当L为正算子时，求解算子方程的问题与求泛函的极小值问题等价，即与泛函的变分问题等价。 3.1.4 算子方程的泛函公式 1、静态场 对于静电场和恒定磁场，泛函I有明确的物理意义，它代表场域中的总位能，即当总位能最小时，场是稳定的（汤姆逊定理），因此，对应于无界空间中的算子方程 的泛函形式应该为 泊松方程的变分公式 泊松方程 为了得到正算子L，改写上式 对应的算子方程 式中，，若材料为均匀，为常数。 边界条件： 相应的泛函为 有内积的定义（在单元中可以认为是常数） 根据格林定理 泛函可以写为 由于算子方程与变分等价，最后一项是在泛函的H空间的边界г上积分，因此，D空间（即空间）边界条件不能直接代入。应该将泛函的被积函数写成D空间的积分形式，再代入边界条件，即 因此泛函可以写为 也可以写成能量泛函 由于第二、三类边界条件已包含在泛函中，其极值问题就只需要满足第一类边界条件。（强调①从偏微分方程要求函数二阶连续偏导，降低为一阶连续偏导，得到弱解，适用范围更广；②为不同媒质分界面上自动满足场的切向分量连续或法向分量连续的证明打下基础；③因为泛函具有能量的意义和量纲，故又称能量泛函，描述静电现象的“最小作用原理”，即汤姆逊定理。） 恒定磁场方程的变分公式 矢量泊松方程 若材料是线性的 边界条件 （场的切向分量，位函数的法向分量） 相应的泛函为 根据矢量恒等式 有 泛函可以写为 有高斯散度定理 泛函的一般表达式为 根据矢量恒等式 因此 同上理由，边界条件不能直接代入泛函，将被积函数写成积分形式后再代入边界条件，则有 泛函可