岩土工程数值方法有限单元法.pptx

6 单元和插值函数的构造; 一般说来，单元类型和形状的选择依赖于结构或总体求解域的几何特点、方程的类型及求解所希望的精度等因素，而有限元的插值函数则取决于单元的形状，结点的类型和数目等因素。;例如在图6.1上，一个二维域利用一系列三角形或四边形单元进行离散，即将总体求解域理想化为由很 多子域（单元）所组成。 在一般情况下，总体域也可能是一维或三维的，在图6.2上分别给出只具有端结点或角结点的一维、二维和三维单元的几种可能形式。一维单元可以简单地是一直线，二维单元可以是三角形、矩形或四边形，三维单元可以是四面体、五面体、长方体或一般六面体。具有轴对称几何形状和轴对称物理性质的三维域能用二维单元绕对称轴旋转形成的三维单元进行离散。 从结点参数的类型上区别、它们可以是只包含场函数的结点只，也可能同时包含场函数导数的结点值。这主要取决于单元交界面上的连续性要求,;而后者又由泛函中场函数导数的最高阶次所决定。如果泛函中场函数导数的最高阶为一次，则单元交界面上只要求函数值保持连续，即要求单元保持连续性。在次情况下，通常结点参数只包含场函数的结点值。如果泛函中场函数导数最高阶为2次，则要求场函数的一阶导数在交界面上也保持连续，即要求单元保持 连续性，这时结点参数中必须同时包含场函数即其一阶导数的结点值。 关于单元插值函数的形式，有限单元法中几乎全部采用不同阶次幂函数的多项式。如果采用幂函数多项式作为单元的插值函数，对于只满足 连续性的单元（称 型单元），单元内的未知场函数的线性变化能够仅用角（或端）结点之间的边界上适当配置一个边内结点（如图6.3所示）它的三次变化，则必须在;每个边界上配置二个边内结点（如图6.4所示）。配置边内结点的另一原因是常常要求单元的边界是曲线的，沿边界配置适当的边内结点从而可能构成二次或更高次多项式来描述它们。有时为使插值函数保持为一定阶次的完全多项式可能还需要在单元内部配置结点。;6.2 一维单元;;求的全部。因为完备性还要求 型单元场函数的一阶导数应包含常数项???这些将在下一章中讨论。 关于插值函数 的构造，为避免繁琐的推导，不必按前面所述步骤进行，而是直接采用熟知的Lagrange插值多项式。对于n个结点的一维单元， 可采用n-1次Lagrange插值多项式 ，即令;其中 的上标n-1表示Lagrange插值多项式的次数 表示二项式在j的范围内(j=1,2,…,i-1,i+1,…,n)的乘积，n是单元的结点数， 是n个结点的坐标。;如果引入无量钢坐标;;其中 是单元中心的坐标，则对于n=2，有 ;为今后构造其它形式的Lagrange单元方便，在此可将(6.2.6)式改写成;其他形式Lagrange单元的插值函数是有帮助的。 还应指出： 的n-1次完全多项式。它的项数和结点数相同且包含常数项，这样构成的场函数模式是满足收敛准则的。特别地，如令 ，则可从(6.2.1)式得到 ;6.2.2 Hermite单元;其中Hermite多项式具有以下性质;;或;6.3 二维单元;从前面的讨论中已知，对于3结点三角形单元，引入面积坐标;1. 二次单元; 如图6.4所示，二次单元有六个结点，各结点的面积坐标分别标注在括号内。参照6.2.12式，现将需要构造的插值函数表示成 ;;2. 三次单元;对于角结点;对于中心结点;其中p为插值函数的次数。显然，按上式构造的插值函数满足 这一基本要求。 另外，由于 的上述性质，以及结点的数目和配置符合Pascal三角形的要求，可以证明这种单元场函数是满足收敛准则的，当然 这一要求也是恒被满足的。 还可指出，当引入面积坐标后，单元矩阵经常可以表示成下列形式的积分 而此积分可以方便地进行计算，这是三角形面积坐标单元的又一优点。 ;6.3.2 Larange矩形单元和Hermite矩形单元;1. Larange矩形单元;现考虑图6.6所示单元，其中一系列结点布置在的插值函数 ，我们已知lagrange多项式 ; 在结点i上等于1，而在其余所有结点上等于0。这种单元每一边界上的结点数和函数在边界上的变化是协调的。因而也保证了单元之间函数的协调性。 图6.7所示为三种形式的Lagrange矩形单元。虽然构造它们的插值函数是很任意的，但是这种类型的单元存在一定缺点，主要是出现了随插值函数方次增高而增