Chapter6有限元分析中及单元性质特征及误差处理.pptVIP

6.1单元节点编号与带宽存储 6.2形状函数矩阵与刚度矩阵的性质 6.3边界条件的处理与支反力计算 6.4单元刚度矩阵的缩聚 6.5为以函数构造与收敛性要求 6.6C0型单元与C1型单元 6.7单元的拼片试验 6.8有限元分析数值解的精度与性质 6.9单元应力的计算结果的误差与平均处理 6.10控制误差和提高精度的h方法和p方法 6.1单元节点编号与带宽存储 计算机进行有限元分析时，需要存储所有单元和节点信息，随着所求解问题自由度的增大，计算规模的增大，整体刚度矩阵的规模非常巨大。 由于整体刚度矩阵中显现出相邻单元之间的关联性，因此矩阵中的大部分数据都为零，反映非零数据的一个指标就是带宽。 由于刚度矩阵是对称的，可以看出，若节点的自由度数目为m，则每一个单元在整体刚度矩阵的半带宽为 di=（第i个单元中节点编号的最大差值+1）＊m d=max（ di ） （i=1，2……n） 其中n为整个结构系统的单元数。显然 对于二维问题，m=2 对于三维问题，m=3 6.2形状函数矩阵与刚度矩阵的性质 以一维杆单元为例，杆单元的位移场为 6.2形状函数矩阵与刚度矩阵的性质 仍然以一维杆单元为例，它的刚度方程为 1、考虑单元左端发生单位位移，右端固定情况 2、考虑单元右端发生单位位移，左端固定情况 3、考察刚体位移 性质1：单元刚度矩阵的对角元素kii表示要使单元的第i个节点产生单位位移，而其它的节点位移为0时，需要在i点施加的节点力。 性质2：单元刚度矩阵的对角元素kij（i≠j）表示要使单元的第j个节点产生单位位移，而其它的节点位移为0时，需要在i点施加的节点力。 性质3：单元刚度矩阵是对称的。这可以由功的互等定理得到。对于线弹性体，力所做的功跟加载次序无关，这可以利用上面的性质1和2得到。 性质4：单元刚度矩阵是半正定的。 性质5：单元刚度矩阵是奇异的。 性质6：单元刚度矩阵的任意行或列代表一个平衡力系，当节点位移全部为线位移时，任意行或列的代数和应该为0。 同样，由单元刚度矩阵所组装的整体刚度矩阵也有以下性质： 1）对称性 2）奇异性 3）半正定性 4）稀疏性 5）非零元素呈现带状分布 6.3边界条件的处理与支反力的计算 位移边界条件在大多数情况下有两种类型。 1、零位移边界条件 2、给定具体数值的位移边界条件 根据上述两类边界条件，刚度方程的求解有以下几种方法： 1、直接法 2、置“1”法 3、乘大数法 4、罚函数法 直接法 置“1”法 乘大数法 罚函数法 罚函数法的最大好处是可以直接求出位移边界上的支反力。 支反力的计算： 除了罚函数法能够求出支反力以外，其它的方法都需要求解一定的方程得到。 6.4单元刚度矩阵的缩聚 采用高次位移函数的单元也常被称为高阶单元。对于高次单元来说，除了几何端点以外，其余的那些节点可能与其它的单元不发生关系，当中间的节点与其它单元无关时，我们称作是内部节点。而其余的节点是外部节点。既然内部节点与其他单元无关，那么在组成整体刚度之前，就可以把他们消去，也就是把内部节点的位移用外部节点的位移来表示。 以一维三节点杆单元为例 以一维三节点杆单元为例 6.5位移函数构造与收敛性要求 单元中的位移模式一般采用设有待定系数的有限多项式作为近似函数，优先多项式的选取原则应该考虑以下几个方面： 1、待定系数是由节点位移条件确定的，因此它的个数应该与节点位移DOF个数相等。 2、在选取多项式时，必须选择常数项和完备的一次项。单元位移模式中的常数项和一次项可以反映单元的刚体位移合唱应变的特性。这是因为当划分的单元数趋于无穷时，即单元缩小趋于一点，此时单元应变趋于常数。 3、选择多项式应该由低到高，尽量选取完全多项式以提高单元的精度。 因此，在构造一个单元的位移函数时，应该参考由多项式函数构成的Pascal三角形和上述原则进行函数项次的选取与构造。 收敛性问题 在有限元分析中，当节点数目或单元插值函数的项数趋于无穷大时，即单元尺寸趋于零时，最后的解答如果能够无线的逼近准确解，那么这样的位移函数或形函数是逼近于真实的，这就称为收敛。 为使有限元分析的解答收敛，位移函数必须满足一些收敛准则，这些准则都经过过严密的理论验证。主要包括以下三个方面。 收敛性准则 定义：当单元尺寸趋于零时，有限元的解趋于真实解。