11-有限元若干问题讨论-2.pptVIP

有限元分析中的若干问题讨论 提高计算精度的h方法和p方法 二、提高计算精度的p方法(p-version或p-method)： 保持有限元的网格剖分固定不变，增加各单元上基底函数的阶次，从而改善计算精度。大量的实践表明：p方法的收敛性大大优于h方法；p方法的收敛性可根据Weierstrass定理来论证；由于p方法使用高阶多项式作为基底函数，会出现数值稳定性问题，另外，由于计算机容量和速度的限制，多项式的阶次不能太高（一般情况下多项式函数的最高阶次p9）。 有限元分析中的若干问题讨论 单元位移函数构造与收敛性要求 　　单元中的位移模式(displacement model)一般采用设有待定系数的有限多项式作为近似函数，有限多项式选取的原则应考虑以下几点： 　　(1) 待定系数是由节点位移条件确定的，因此它的个数应与节点位移DOF数相等。如平面3节点三角形单元共有六个节点位移DOF，待定系数的总个数应取6个，则两个方向的位移u和v应各取三项多项式。对于平面4节点的矩形单元，待定系数为8，位移函数可取四项多项式作为位移模式。 一、单元位移函数的选择 1、选择单元位移函数的一般原则 有限元分析中的若干问题讨论 单元位移函数构造与收敛性要求 　　(2) 在选取多项式时，必须要选择常数项和完备的一次项。位移模式中的常数项和一次项可以反映单元刚体位移和常应变的特性。这是因为当划分的单元数趋于无穷时，即单元缩小趋于一点，此时单元应变应趋于常数，3节点三角形单元的位移模式正好满足这个基本要求。 　　(3) 选择多项式应由低阶到高阶，尽量选取完全多项式以提高单元的精度。一般情况下对于每边具有2个端节点的单元应选取一次完全多项式的位移模式，每边有3个节点时应取二次完全多项式；若由于项数限制不能选取完全多项式时，选择的多项式应具有坐标的对称性。并且一个坐标方向的次数不应超过完全多项式的次数。 1、选择单元位移函数的一般原则 有限元分析中的若干问题讨论 单元位移函数构造与收敛性要求 一、选择单元位移函数的一般原则 2、多项式函数构成的帕斯卡三角形(Pascal triangle) 图5-8 图5-9 有限元分析中的若干问题讨论 单元位移函数构造与收敛性要求 一、选择单元位移函数的一般原则 3、常用的二维单元和三维单元 有限元分析中的若干问题讨论 单元位移函数构造与收敛性要求 二、关于收敛性问题 １、收敛(convergence) 　　在有限元分析中，当节点数目趋于无穷大时（即当单元尺寸趋近于零时）或单元插值位移的项数趋于无穷大时，最后的解答如果能够无限地逼近准确解，那么这样的位移函数（或形状函数）是逼近于真解的，这就称为收敛。 有限元分析中的若干问题讨论 单元位移函数构造与收敛性要求 二、关于收敛性问题 为使有限元分析的解答收敛，位移函数必须满足一些收敛准则(convergence criterion)： 当单元尺寸趋于零时，其位移函数及其应变总是趋向于某一常数。 　　对位移函数的基本要求应当是：函数本身应在单元上连续，还要包括使得位移函数及对应于应变的导数都为常数的项，即常位移项和常应变项。 要保证单元的收敛性，还要考虑单元之间的位移协调。不仅节点处的位移应协调，沿整个单元边界上的位移都应当是协调的（或相容的） 有限元分析中的若干问题讨论 单元位移函数构造与收敛性要求 二、关于收敛性问题 收敛性准则1：完备性要求(针对单元内部） 如果在（势能）泛函中所出现位移函数的最高阶导数是m阶，则有限元解答收敛的条件之一是选取单元内的位移场函数至少是m阶完全多项式。 二维问题和三维问题的完全多项式参见图5-8和图5-9。可以看出所要求的m阶完全多项式已包含了刚体位移和常应变项。 收敛性准则2：协调性要求(针对单元之间） 如果在（势能）泛函中位移函数出现的最高阶导数是m阶，则位移函数在单元交界面上必须具有直至(m-1)阶的连续导数，即Cm-1连续性。 下面就一般的平面问题进行讨论。 有限元分析中的若干问题讨论 单元位移函数构造与收敛性要求 三、平面单元位移函数选取的要求 在平面问题中，势能泛函为 可以看出所出现的物理量关于位移的最高阶导数是1，因此m=1。 由准则1，形状函数至少应包含完整的一次多项式，即 （5-97） （5-98） 有限元分析中的若干问题讨论 单元位移函数构造与收敛性要求 三、平面单元位移函数选取的要求 （5-97） （5-98） 这代表刚体位移和常应变的位移模式。在平面3节点单元和平面4节点单元中，其位移模式都包含了(5-97)和(5-98)的多项式。 由准则2，平面3节