偏微分方程的有限元方法SetinArtile.pdf

16.901 讲义笔记 2002 年 4 月 17 日 一维高阶有限元方法 重点： ﹡节点基的扩展 ﹡参考元中的高阶检验函数 ﹡残差积分法和刚度矩阵的组成 节点基的扩展 直到现在，我们已经考虑了一些方程的解，并允许这些 解在跨越单元时(至多)有线性的变化： 在一个单元内 为常数 在归纳高阶情形时，我们考虑在一个单元中的二次变分： 在各单元为常数 注意： v(x)只有在节点处是连续的。 为了定义一维单元中一般的二次函数，我们需要三个自由 度： 2 V( x) c =+c x+c x quad 0 1 2 有三个自由度（c ，c ，c ）的一般二次函数。 0 1 2 回忆一下，对于线性元我们有两个自由度： V( x) c =+c x linear 0 1 在这个例子中，我们用了“节点基”这个概念使得未知量在 节点对应V( x) 的值： V( x) Vζ x =+V ζ x 1 1( ) 2 2 ( ) element 这里， ， 是两个自由度 V V 1 2 上图中，ζ x ，ζ x 是两个线性内插式 1( ) 2 ( ) ⎧ 0, if i ≠ j 已知： ⎪ ζ (x ) ⎨ i j 1, if i j ⎪⎩ 对于二次的情况，我们通过在单元内部增加一个节点来扩展 节点基，通常的选择是将增加的节点选在单元的中点处： 新的内插节点选在 单元区域 所以，在一个单元内： V( x) Vζ x +V ζ x +Vζ x element 1 1( ) 2 2 ( ) 3 3 ( ) 这里，ζ (x) 是二次函数，通常满足 i ⎧ 0, if i ≠ j ζ (x ) ⎨⎪ i j 1, if i j ⎪⎩ 下面给出两倍数量的节点集合： 可见，对于二次元而言，N 个单元对应 2N+1 个未知量，而 线性单元只有 N+1 个未知量。 参考元中的二次内插式 回忆前面我们学过的将 x 映射为ξ 的参考元，对于二次参考 元，就像下图所示： 也就是：ξ =−1 ，ξ 0 ，ξ 1 1 2 3 下一步就需要我们来找出ζ (ξ) ，ζ (ξ) 和ζ (ξ) 。 1 2 3 经过很多步工作，我们求得： 画出 为了求 x ，我们先给定ξ ，在一个单元中利用等参数映射内 插x(ξ) ，得： x(ξ) xζ (ξ) =+x ζ (ξ) +x ζ (ξ) 1 1 2 2