有限元的性质和收敛性.pdfVIP

有限元的性质和收敛性 一、有限元解的收敛准则 有限单元法作为求解数学微分方程的一种数值方法可以认为是里兹法的一种特殊形式，不同 在于有限单元法的试探函数是定义于单元 （子域）而不是全域。因此有限元解的收敛性可以与里 兹法的收敛性对比进行讨论。里兹法的收敛条件是要求试探函数具有完备性和连续性，也就是说， 如果试探函数满足完备性和连续性要求，当试探函数的项数n∞时，则Ritz 法的近似解将趋 近于数学微分方程的精确解。现在要研究什么是有限元解的收敛性提法？收敛的条件又是什么？ 在有限单元法中，场函数的总体泛函是由单元泛函集成的。如果采用完全多项式作为单元的 插值函数（即试探函数），则有限元解在一个有限尺寸的单元内可以精确地和真正解一致。但是 实际上有限元的试探函数只能取有限项多项式，因此有限元解只能是真正解的一个近似解答。有 限元解的收敛准则需要回答的是，在什么条件下当单元尺寸趋于零时，有限元解趋于真正解。 下面仍以含有一个待求的标量场函数为例，微分方程是： A( φ) = L( φ) + b = 0 (1.1) 相应的泛函是： （1.2） 假定泛函∏中包含φ和它的直至m 阶的各阶导数，若 m 阶导数是非零的，则近似函数 至少必 须是m 次多项式。若取p 次完全多项式为试探函数，则必须满足p ≥m ，这时 及其各阶导数在 一个单元内的表达式如下： （1. ３） 由上式可见，由于 是p 次完全多项式，所以它的直至m 阶导数的表达式中都包含有常数项。 但单元尺寸趋近于零时，在每一单元内 及其直至 m 阶导数将趋近于它的精确值，即趋近于常 数。因此，每一个单元的泛函有可能趋于它的精确值。如果试探函数还满足连续性要求，那么整 个系统的泛函将趋近于它的精确值。有限元解就趋近于精确解，也就是说解是收敛的。 从上述讨论可以得到下列收敛准则： 准则1 完备性要求。如果出现在泛函中场函数的最高阶导数是 m 阶，则有限元解收敛的条件之一 是单元内场函数的试探函数至少是m 次完全多项式。或者说试探函数中必须包括本身和直至m 阶导数为常数的项。 单元的插值函数满足上述要求时，我们称单元是完备的。 1 至于连续性的要求，当试探函数是多项式的情况下，单元内部函数的连续性显然是满足的， 如试探函数是m 次多项式，则单元内部满足 Cm-1 连续性要求。因此需要特别注意的是单元交 界面上的连续性，这就提出另外一个收敛准则。 准则2 协调性要求。如果出现在泛函中的最高阶导数是m 阶，则试探函数在单元交界面上必须具 有Cm-1 连续性，即在相邻单元的交界面上应有函数直至m- 1 阶的连续导数。 当单元的插值函数满足上述要求时，我们称单元是协调的。 简单的说，当选取的单元既完备又协调时，有限元解是收敛的，即但单元尺寸趋于零时，有 限元解趋于真正解。 需要补充指出的是：前面所述有限元解收敛于数学微分方程精确解的进一步含义。因为数学 微分方程的精确解往往不一定能够得到，甚至问题的数学微分方程并未建立（例如对于复杂型式 的结构）。同时有限元解中通常包含多种误差 （例如计算机的截断误差和舍入误差），因此有限元 解收敛于精确解，在更严格意义上说是问题的有限元解的离散误差趋于零。 二、收敛准则的物理意义 为了从物理上加深对收敛准则的理解，我们以平面问题为例加以说明。 在平面问题中，泛函∏p 中出现的是位移u 和v 的一次导数，即εx ，εy ，γxy ，因此m= 1。 收敛准则1 要求插值函数或者位移函数至少是x ，y 的一次完全多项式。我们知道位移及其一阶 导数为常数的项是代表与单元的刚体位移和常应变状态相应的位移模式。实际分析中，各单元的 变形往往包含着刚体位移。因而当单元尺寸趋于无穷小时，各单元的应变也总是趋于常应变。所 以完备性的要求由插值函数所构成的有限元解必须能反映单元的刚体位移和常应变状态。若不能 满足上述要求，那么赋予结点以刚体位移 （零应变）或者常应变的位移值时，在单元内部将产生 非零或非常值的应变，这样有限元解将不可能收敛于真正解。