B8、 弹性力学有限元法基本原理(二).ppt

第三单元 弹性力学有限元法基本原理（二） 第一节 有限元解的性质和收敛准则 1、有限元解的收敛准则 有限单元法作为连续问题的数值解法可以看作里兹法的一种特殊形式，即采用分片试探函数（假定位移场）的里兹法。前面通过受轴向力杆的里兹法求解，已经指出里兹解收敛必须满足的条件：除了满足连续性和边界约束之外，试探函数还必须是完备的，即包含完全的低阶多项式。 对于有限元法，解的收敛除了里兹法意义上的收敛外，显然，当单元尺寸趋于零（h→ 0)时，有限元解应该趋于问题的精确解。这就是通常意义上有限元解收敛的涵义（h-收敛）。 事实上，有限元位移法中，在一个单元内用完全多项式逼近实际位移场，当单元尺寸趋于零时，如果整体位移试探函数还满足连续性要求，那么在满足最小势能原理的情况下，整个系统的势能泛函将趋于它的精确值——最小值。在每一单元内位移及其导数将趋于它的精确值（常数），有限元解就趋于精确解，即解是收敛的。 根据以上分析，对弹性力学有限元法，为了使有限元解收敛，单元（一维杆，二、三维实体元）的构造必须满足下列要求： 每个单元的位移模式必须包含完全一次多项式。 位移模式在单元边界之间连续（C0连续）。 单元网格在边界上受到均匀载荷，单元上的有限元解应该具有一致的均匀值。 上述要求可以概括为两个收敛准则： 准则1 ：完备性要求 对弹性力学问题，单元位移模式必须包含一次完全多项式。 满足上述要求的单元称为完备单元。 准则2 ：协调性要求 对弹性力学问题，位移试探函数在单元交界面上必须具有C0连续性（函数值连续）。 满足上述要求的单元称为协调元。 理论上可以证明，同时满足完备性和协调性的单元一定收敛。但协调性不是收敛的必要条件，某些具有非协调位移模式的单元只要满足一定条件也是收敛的。 2、对收敛性和收敛准则的理解 根据前面分析，对于有限元位移法，有两个途径得到不断逼近精确解的有限元解序列：第一，网格不变，不断增加位移模式多项式的阶数；第二，单元位移模式不变，不断增加单元数，即单元尺寸趋于零。通常所指有限元解的收敛性是第二种情况。 关于有限元解的收敛性和收敛准则，数学家已经给出严格的证明。下面以弹性力学问题为例从物理概念上进行理解。 准则2的协调性要求是连续体力学问题的必然要求。它是最小势能原理和里兹法的前提条件。有限元法作为里兹法的特殊形式必然要满足这个要求。有限元的协调性要求在整个弹性体区域上的体现就是试探位移场必须满足的连续性条件。事实上，如果单元尺寸趋于零时，单元交界面上位移不连续，则有限元模型模拟的就不可能是原来的连续结构，获得的有限元解就不可能收敛到问题的真正解。 在有限元法中，一般在粗网格下单元要满足协调性要求。如果某单元在粗网格下不满足协调性，但随着单元尺寸减小，不协调性趋于消失，同时满足完备性，则该单元也能收敛。 不难证明，3节点三角形单元满足完备性准则和协调性准则。 3、收敛速度、精度及其意义 如果单元位移模式满足完备性和协调性，则当单元尺寸趋于零时，有限元解趋于精确解。 根据里兹法的原理，如果单元的位移插值多项式能够精确拟合真正解，则很粗糙的单元划分就能得到精确的解答。比如，假设位移精确解是二次函数，而单元位移模式包含了完全二次多项式，则有限元解一定是精确的。 对于一般的实际位移场，一点附近的位移可以展开为Taylor级数。根据前面结论，在一个单元范围内，有限元解可以拟合实际位移的低阶成分，而忽略的高阶成分就是误差。设单元直径是h，单元位移模式包含p阶完全多项式，则它可以在单元上拟合实际位移Taylor展开中的前p阶。因此有限元位移解的误差是 ，这只是一种量级的估计，不反映误差的绝对数值，但可以反映收敛速度。 例如，3节点三角形单元，位移模式是线性的，所以位移的误差估计是 ，也可以说该单元的收敛速度是 量级。即，如果把有限元网格细化，单元尺寸减半，则有限元位移解的误差大约是前一种网格的(1/2)2 = 1/4。显然，该单元应变、应力解的收敛速度是 量级。 实际工作中，往往需要对误差作出具体估计，对于一般的实际问题，可采取下列办法： （1）用相近的有已知解析解的问题做有限元误差估计，单元类型相同，网格划分相似。则某种网格下其有限元解与解析解的具体误差可以作为实际问题的误差。 （2）根据收敛的含义，可以对网格进行连续多次细化，当两次网格的解相差不大时，可以认为得到的解答足够精确。 （3）利用收敛速度的量级估计精确解。有限元解是单调收敛的，对3节点三角形单元，设第一次网格的位移