第3章单元和单元插值函数.pdfVIP

第三章 单元和单元插值函数 §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 2 是插值函数 （或称形函数），N 1