偏微分方程的有限元方法分析.pptVIP

只利用插值节点处函数值得信息构造出来的插值函数称作是Lagrange型插值函数。 除利用插值节点处的函数值信息，还利用插值节点处的导数信息作为插值条件构造出来的插值函数称作是Hermite型插值函数。 就整体而言，Hermite型插值函数的光滑程度会有所提高。 一般情况下，高次元有限元解的精度比线性有限元解的精度高。但在形成和求解有限元方程组的过程中所花费的工作量也随之猛增。通常除特殊需要，一般不必采用高次元，而是采用线性元、网格细分、外推等措施。 除了三角形元和矩形元，还可考虑任意(甚至是曲边)的四边形单元及曲边三角形单元。 (b) 从虚功原理及节点基函数出发形成有限元方程。 5 用有限元方法解题主要步骤 ② 对求解域作网格。 ③ 构造基函数（或单元形状函数）。 ⑤ 求解线性代数方程组。 ④ 形成有限元方程。 (a) 从变分原理及单元形状函数出发，先形成单元刚度矩阵，再由单元刚度矩阵叠加成总刚矩阵。 ① 把边值问题化成等价的变分问题（或建立等价的变分方程）。 线性无关，它们可组成试探函数空间的基，常称为节点基函数。 几何形状如图 a b 任一试探函数 可表示为 用这类插值型基函数，可以构造出适合各种边界条件的试探函数。 若借助前述放射变换 节点基函数可用变量?表示为 ① 直接形成有限元方程 (a) 把表达式 代入泛函； （4）从Ritz方法出发形成有限元方程 (b) 将泛函表达式中积分区间[a,b]变到[0,1]； (c) 由 达到极小值的条件 得到含 的有限元方程 这儿 (d) 解出有限元方程的数值解 ，就得到使二次泛函取极小的近似函数(有限元解) 有限元方程可用矩阵表示为 其中 称为总刚矩阵。 工程中形成有限元方程时，通常先在每个单元上形成单元矩阵（称为单元刚度矩阵），然后由单元刚度矩阵形成总刚度矩阵（称为总体合成）。 ② 用单元刚度分析形成有限元方程 (a) 把 按单元组织，则在第i个单元上，令 其中 称为单元刚度矩阵。各元素可计算得到。 再把 扩展成n?n矩阵，使其第i?1行、第i行和第i?1列、第i列交叉位置的元素就是单元刚度矩阵的四个元素，其余全为零（ 只是第一行，第一列元素非零）。即 记 则 其中 称为总刚矩阵。 (b) 由 达到极小值的条件 (c) 解出有限元方程的数值解 ，就得到使二次泛函取极小的近似函数(有限元解) 得到有限元方程 。 （5）从Galerkin方法出发形成有限元方程 把表达式 代入变分方程 对前面的两点边值问题，变分方程变为 其中 与Ritz方法相比， Galerkin方法形成的有限元方程其系数矩阵就是总刚矩阵。 该方程即为Galerkin法形成的有限元方程。 由Galerkin方法推导有限元方程更加方便直接，且适用面广。 若希望在每个单元上提高逼近的精确度，则可通过提高插值多项式次数来实现， 在单元 上可构造一、二、三及高次插值多项式，其方法有两种： 2 一维问题的高次元 整个问题计算的全过程除分析单元插值外，均与前面框架类似。 ① Lagrange型：在单元内部增加一些插值节点。 ② Hermite型：在节点引进一阶、二阶乃至更高阶导数。 ① 线性元(Lagrange型) 要求：在每一个单元上是一次多项式，在单元节点处连续。 插值条件：在单元的两个端点取指定值。 ② 二次元 (Lagrange型) 要求：在每一个单元上是二次多项式，在单元节点处连续。 插值条件：在单元的两个端点及单元中点取指定值。 ③ 三次元（Hermite型） 要求：在每一个单元上是三次多项式，在单元节点处连续。 插值条件：在两个端点取指定的函数值和一阶导数值。 采用高次元，有限元方程形成的方法和线性元类似，但工作量增加。一是计算积分的复杂性增加，二是矩阵的带宽增加。 高次元的主要优点是收敛阶高，且提高了函数逼近的光滑性。 假定区域G可以分