第03章 平面问题的有限元法.ppt

　　⑵ 位移模式必须能包含单元的常应变。每个单元的应变一般都是包含着两个部分：一部分是与该单元中各点的坐标位置有关的应变（即所谓各点的变应变）；另一部分是与位置坐标无关的应变（即所谓的常应变）。从物理意义上看， 　　为了保证解答的收敛性，要求位移模式必须满足以下三个条件，即 　　⑴ 位移模式必须包含单元的刚体位移。也就是说，当节点位移是由某个刚体位移所引起时，弹性体内将不会产生应变。所以，位移模式不但要具有描述单元本身形变的能力，而且还要具有描述由于其它单元形变而通过节点位移引起单元刚体位移的能力。例如，在§3-2节的位移模式(b)中，常数项?1、?4 就是用于提供刚体位移的。? ⑶ 位移模式在单元内要连续、且在相邻单元之间的位 移必须协调。当选择多项式来构成位移模式时，单元内的 连续性要求总是得到满足的，单元间的位移协调性，就是 要求单元之间既不会出现开裂也不会出现重叠的现象。通 常，当单元交界面上的位移取决于该交界面上节点的位移 时，就可以保证位移的协调性。 当单元尺寸无限缩小时，每个单元中的应变应该趋于常量。 因此，在位移模式中必须包含有这些常应变，否则就不可 能使数值解收敛于正确解。很显然，在前面的位移模式 (b)中，与?2、?3、?5、?6 有关的线性项就是提供单元中的 常应变的。 在有限单元法中，把能够满足条件1和2的单元，称为完备单元；满足条件3的单元，叫做协调单元或保续单元。前面讨论过的三角形单元和矩形单元，均能同时满足上述三个条件，因此都属于完备的协调单元。在某些梁、板及壳体分析中，要使单元满足条件3比较困难，所以实践中有时也出现一些只满足条件1和2的单元，其收敛性往往也能够令人满意特别是放松条件3的单元，即完备而不协调的单元，已获得了很多成功的应用。对于不协调单元，其主要的缺点是不能事先确定其刚度与真实刚度之间的大小关系。但值得指出的是，不协调单元一般不象协调单元那样刚硬（即比较柔软），因此有可能会比协调单元收敛得快。 　　在选择多项式作为单元的位移模式时，其阶次的确定，要考虑解答的收敛性，即单元的完备性和协调性要求。实践证明，这两项确实是所要考虑的重要因素，但并不是唯一的 因素。选择多项式位移模式阶次时，需要考虑的另一个因素是，所选的模式应该与局部坐标系的方位无关，这一性质称为几何各向同性。对于线性多项式，各向同性的要求通常就等价于位移模式必须包含常应变状态。对于高次位移模式，就是不应该有一个偏惠的坐标方向，也就是位移形式不应该随局部坐标的更换而改变。经验证明，实现几何各向同性的一种有效方法是，根据图3-10所示的巴斯卡三角形来选择二维多项式的各项。在二维多项式中，如果包含有对称轴一边的某一项，那么就必须同时包含有另一边的对称项。 　　选择多项式位移模式时，还应该要考虑的一个因素就是，多项式中的项数必须等于或稍大于单元边界上的外节点的自由度数。通常是取项数与单元的外节点的自由度数相等，取过多的项数是不恰当的。 ⒋ 刚度矩阵[K]是一个稀疏矩阵。 如果遵守一定的节点 编号规则，就可使矩阵的非零元素都集中在主对角线附 近呈带状。 前面在讨论总刚子矩阵的计算时曾指出，总刚中第 r双行的子矩阵[Krs ]，有很多位置上的元素都等于零， 只有当第二个下标s等于r或者s与r同属于一个单元的节 点号码时才不为零，这就说明，在第r双行中非零子矩 阵的块数，应该等于节点r周围直接相邻的节点数目加 一。可见，[K]的元素一般都不是填满的，而是呈稀疏 状（带状）。 以图3-6a所示的单元网格为例，其整体刚度矩阵中 的非零子块（每个子块为2行2列）的分布情况如图3-6b 所示。 图 3-6 a 图 3-6 b 半带宽B=（相邻节点号的最大差值D+1）*2 若第r双行的第一个非零元素子矩阵是[Krl ]，则从[Krl ] 到[Krr ]共有（r-l+1）个子矩阵，于是[K]的第2r行从第一个 非零元素到对角元共有2 (r-l+1) 个元素。显然，带状刚度矩 阵的带宽取决于单元网格中相邻节点号码的最大差值D。我 们把半个斜带形区域中各行所具有的非零元素的最大个数 叫做刚度矩阵的半带宽（包括主对角元），用B表示，即 B = 2 (D+1)。 通常的有限元程序，一般都利用刚度矩阵的对称和稀 疏带状的特点，在计算求解中，只存储上半带的元素，即 所谓的半带存储。因此，在划分完有限元网格进行节点编 号时，要采用合理的编码方式，使同一单元中相邻两节点 的号码差尽可能小，以便节省存储空间、提高计算效率。 ⒌