第4章_平面问题的有限元法-4收敛准则.ppt

* 　　可以证明，对于一个给定的位移模式，其刚度系数的数值比精确值要大。所以，在给定的载荷之下，有限元计算模型的变形将比实际结构的变形小。因而，当单元网格分得越来越细时，位移的近似解将由下方收敛于精确解，即得到真实解的下界。 第七节 收敛准则 　　对于一个数值计算方法，一般总是希望随着网格的逐步细分所得到的解答能够收敛于问题的精确解。根据前面的分析，我们知道，在有限元分析中，一旦确定了单元的形状之后，位移模式的选择将是非常关键的。由于载荷的移置、应力矩阵和刚度矩阵的建立等等，都依赖于单元的位移模式，所以，如果所选择的位移模式与真实的位移分布有很大的差别，那么就很难获得良好的数值解。 　　⑵ 位移模式必须能包含单元的常应变。每个单元的应变一般都是包含着两个部分：一部分是与该单元中各点的坐标位置有关的应变（即所谓各点的变应变）；另一部分是与位置坐标无关的应变（即所谓的常应变）。从物理意义上看， 　　为了保证解答的收敛性，要求位移模式必须满足以下三个条件，即 　　⑴ 位移模式必须包含单元的刚体位移。也就是说，当节点位移是由某个刚体位移所引起时，弹性体内将不会产生应变。所以，位移模式不但要具有描述单元本身形变的能力，而且还要具有描述由于其它单元形变而通过节点位移引起单元刚体位移的能力。 例如，三角形三节点单元位移模式中，常数项?1、?4 就是用于提供刚体位移的。? ? ⑶ 位移模式在单元内要连续、且在相邻单元之间的位移必须协调。当选择多项式来构成位移模式时，单元内的连续性要求总是得到满足的，单元间的位移协调性，就是要求单元之间既不会出现开裂也不会出现重叠的现象。通常，当单元交界面上的位移取决于该交界面上节点的位移时，就可以保证位移的协调性。 当单元尺寸无限缩小时，每个单元中的应变应该趋于常量。 因此，在位移模式中必须包含有这些常应变，否则就不可 能使数值解收敛于正确解。 很显然，三角形三节点单元位移模式中 ，与?2、?3、?5、?6 有关的线性项就是提供单元中的常应变的。 在有限单元法中，把能够满足条件1和2的单元，称为完备单元；满足条件3的单元，叫做协调单元或保续单元。前面讨论过的三角形单元和矩形单元，均能同时满足上述三个条件，因此都属于完备的协调单元。在某些梁、板及壳体分析中，要使单元满足条件3比较困难，所以实践中有时也出现一些只满足条件1和2的单元，其收敛性往往也能够令人满意特别是放松条件3的单元，即完备而不协调的单元，已获得了很多成功的应用。对于不协调单元，其主要的缺点是不能事先确定其刚度与真实刚度之间的大小关系。但值得指出的是，不协调单元一般不象协调单元那样刚硬（即比较柔软），因此有可能会比协调单元收敛得快。 　　在选择多项式作为单元的位移模式时，其阶次的确定，要考虑解答的收敛性，即单元的完备性和协调性要求。实践证明，这两项确实是所要考虑的重要因素，但并不是唯一的 因素。选择多项式位移模式阶次时，需要考虑的另一个因素是，所选的模式应该与局部坐标系的方位无关，这一性质称为几何各向同性。对于线性多项式，各向同性的要求通常就等价于位移模式必须包含常应变状态。对于高次位移模式，就是不应该有一个偏惠的坐标方向，也就是位移形式不应该随局部坐标的更换而改变。经验证明，实现几何各向同性的一种有效方法是，根据如图所示的巴斯卡三角形来选择二维多项式的各项。在二维多项式中，如果包含有对称轴一边的某一项，那么就必须同时包含有另一边的对称项。 　　选择多项式位移模式时，还应该要考虑的一个因素：多项式中的项数必须等于或稍大于单元边界上的外节点的自由度数。通常是取项数与单元的外节点的自由度数相等，取过多的项数是不恰当的。 常数项 线性项 二次项 三次项 四次项 五次项 对称轴 巴斯卡三角形 第八节 有限元分析的步骤 　根据前面的讨论，现以三角形常应变单元为例来说明应用有限元法求解弹性力学平面问题的具体步骤： ①力学模型的确定根据工程实际情况确定问题的力学模型，并按一定比例绘制结构图、注明尺寸、载荷和约束情况等。 ②将计算对象进行离散化，即弹性体划分为许多三角形单元，并对节点进行编号。确定全部节点的坐标值，对单元进行编号，并列出各单元三个节点的节点号。 ③ 计算载荷的等效节点力（要求的输入信息）。 ④ 由各单元的常数bi 、ci 、bj 、cj 、bm 、cm 及行列式2?，计算单元刚度矩阵。 ⑤ 组集整体刚度矩阵，即形成总刚的非零子矩阵。 ⑥ 处理约束，消除刚体位移。 ⑦ 求解线性方程组，得到节点位移。 ⑧ 计算应力矩阵，求得单元应力，并根据需要计算主应力和主方向。 ⑨ 整理计算结果（后处理部分）。 　 为了提高有限元分析计算的效率、达到一定的精度，应该注意以下几个方面。