10第3章单元和单元插值函数.docVIP

第三章 单元和单元插值函数 § 3-1 引言 前面两章较为详细地介绍了基于最小位能原理建立有限元求解方程的途径。我们知道一旦确定了单 元以及相应的插值函数，随后的工作就是遵照完全确定的标准过程进行了，但是，选择适当的单元和单 元插值函数必须动一些脑筋。一般来说，单元的选择依赖于结构或总体求解域的几何特点以及求解所希 望的精度等因素，而单元插值函数则取决于单元的形状、节点的数目和类型等因素。 在有限单元中，单元插值函数几乎全 部采用不同阶次的多项式。这是因为它 具有易于运算和易于满足收敛性要求的 优点，如果采用多项式作为单元插值函 数，对于单元内线性变化的未知场，则 可以仅用角节点的参数表示（如图 3-1 所示）；对于二次变化的未知场，则必须 在角节点之间的边界上配置一个边内节 点（如图 3-2 所示）；对于三次变化的未 图 3 -1、各种形状只有角节点的单元 知场，则必须在每个边界上配置二个边 内节点（如图 3-3 所示）。配置边内节点的另 一个原因是常常要求单元的边界是曲线的， 沿边界配置适当的边内节点，从而可能构成 二次或更高次多项式来描述它们。有时一定 阶次的完全多形式可能还需要在单元内部配 置节点，然而这些内部节点除非是所考虑的 具体情况绝对必需的，否则是不希望的，因 为这些节点的存在将增加有限元表达式和计 图 3-2、二次单元 算上的复杂性。 更具体的考虑，单元的形式可以区分为： （1） 按几何形状：一维、二维和三维； （2） 按节点参数：Lagrange 族（含场函数的 节点值）或 Hermite 族（还含场函数导 数的节点值）； （3） 按插值函数：Lagrange多项式或Hermite 图 3-3、三次或高次单元 多项式； （4） 按单元坐标：笛卡尔坐标和自然坐标。 134 § 3-2 一维单元 如大家熟知的杆单元或弯曲梁单元，实质上是两种形式的一维单元：杆单元的节点参数中只包含 场函数? （表示杆单元的位移u ）的节点值，单元内的场函数可以插值表示为 ? ? N ? N （3-2-1） 1? ? 1 2 2 其中 N ， 1 N 是插值函数（或称形函数）， N1 ?1 ? x l ，N ? x l 2 2 。 另一种弯曲梁单元的节点参数中，除包含场函数（表示弯曲梁单元的位移 v）的节点值外 ，还包含 场函数导数d? dx的节点值。单元内的场函数可以插值表示为 ? ? ? d d ? ? ? ? 1? ? ? ? （3-2-2） ? N N ? ? ? N ? ? N 1 2 3 2 4 ? dx ? ? dx ? 1 2 其中插值函数为 N 1 N 3 ? ? 1? 3 x l 2 3 2 x2 l 2 ? 2 x ? 2 l 3 3 x3 l 3 N N 2 4 ? ? x2 x ? 2 ? l x x ? ? l l 2 3 2 x3 l 2 实际上这两种单元就是以下将讨论的 Lagrange 单元和 Hermite 单元的最简单形式，现对它们的一 般形式进行讨论。 一、Lagrange 单元 对于具有 n 个节点的一维单元，如果它的节点参数中只含有场函数的节点值，则单元内的场函数可 插值表示为 n ? i i ? ? N ? （3-2-3） i?1 其中插值函数具有以下性质 n N ? ? ， ? ? ? ? ? i x ? ? Ni x j ij i?1 ? 1 （3-2-3a） 关于插值函数 N ?x? i 的构造，为避免繁琐的推导，不必按前一章中所述的步骤进行，可以直接采用 Lagrange 插值多项式。对于 n 个节点的一维单元， N ?x? i 可采用 N ?1次 Lagrange 插值多项式，即 x? x ? ?? ? ? ?? ? ? ? x ? x x ? x x ? x x ? x x ? x n ? ? N ? ? ? j 1 2 i ?1 i?1 n 1 2 i ?1 i?1 n ? ?? ? ? ?? ? ? ? i ? x ? x x ? x x ? x x ? x x ? x j 1, j ?i ? x x ? ? i j i 1 i 2 i i ?1 i i?1 i n （3-2-4） ? 表示二项式在 j 的范围内（ j ? 1,2,?,i ?1,i ?1,?n ）的乘积，n 是单元的节点数，x , x , , x 1 ? 是 2 n n 个节点的坐标。 2 ，则插值函数 N ?x? i 就可得到与式（3-2-1）完全一样的结果。 如果n ? 2 ，并且令 x1 ? 0 ， x ? l 如果引入无量纲坐标 x ? x1 ? x ? x ? ? （3-2-5） 1 x ? x l n 1 135 其中l 表示单元的长度，则式（3-2-4）可表