(新)第四章平面问题的有限元法——精品.ppt

4.5.2 表面力的等效载荷列阵 4.5.3体积力的等效载荷列阵 { 4.3四结点矩形单元 由（3-41）式可以看出，令节点1在坐标轴x方向的位移u1 =1，而其余的节点位移v1 = u2 = v2 = u3 = v3 = …= u2n = v2n =0，这样就可得到节点载荷列阵等于[K]的第一列元素组成的列阵，即 即表示： 是在j节点有单位位移时，而在I节点所需施加的力。 (s) 返回 ⒉ 刚度矩阵[K]中主对角元素总是正的。 例如，刚度矩阵[K]中的元素k33 是表示节点2在x方向产生单位位移，而其它位移均为零时，在节点2的x方向上必须施加的力，很显然，力的方向应该与位移方向一致，故应为正号。 ⒊刚度矩阵[K]是一个对称矩阵，即[Krs ]= [Ksr ]T。 由（3-32）、（3-36）式得 所以，可以只存储上三角或下三角矩阵。 (t) 返回 ⒋ 刚度矩阵[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)。 通常的有限元程序，一般都利用刚度矩阵的对称和稀 疏带状的特点，在计算求解中，只存储上半带的元素，即 所谓的半带存储。因此，在划分完有限元网格进行节点编 号时，要采用合理的编码方式，使同一单元中相邻两节点 的号码差尽可能小，以便节省存储空间、提高计算效率。 返回 ⒌ 刚度矩阵[K]是一个奇异矩阵，在排除刚体位移后，它是正定阵。 弹性体在{R}的作用下处于平衡，{R}的分量应该满足三个静力平衡方程。这反映在整体刚度矩阵[K]中就意味着存在三个线性相关的行或列，所以[K]是个奇异阵，不存在逆矩阵。 因 代入（3-30）得 (u) 返回 上式左乘{?}T ，并注意到 (3-13) 式， 在集合过程中将[B]扩充到3×2n阶后，有[B]3×2n {?}e2n ×1 = [B]3 ×2n {?}2n ×1 ，故 (v) 由于弹性矩阵[D]是正定的，且t和?都是正的，所以只有当每个单元中都有{? }=0时，才有 否则 也就是说，当排除了弹性体的刚体位移{? }=0之后，若 {? }?0，则二次型{? }T [K]{? }恒大于零，于是[K]必定为正定阵。有关排除整体刚度矩阵奇异性的方法将在后面的章节中予以讨论。 返回 第五节 等效节点力载荷列阵 在上节讨论整体刚度矩阵时已经指出，（3-37）式中的载荷列阵{R}，是由弹性体的全部单元的等效节点力集合而成，而其中单元的等效节点力{R}e 则是由作用在单元上的集中力、表面力和体积力分别移置到节点上，再逐点加以合成求得。根据虚位移原理，等效节点力的大小，应按其所做的功与作用在单元上的三种力在任何虚位移上所做的功相等这一原则来确定。即 上式中等号的左边表示单元的等效节点力{R}e 所做的虚功；等号右边的第一项是集中力{G}所做的虚功、第二项的积分是沿着单元的边界进行，表示面力{q}所 (a) 返回 做的虚功、第三项的积分则是遍及整个单元，表示体积力{p}所做的虚功；t为单元的厚度，假定为常量。 　　将上节中的 (c) 式代入上式，并注意到节点虚位移列阵{?* }e 中的元素都是常量，这就可以把({? *}e ) T 提到积分号的外面，于是有 (b) 上式括号中大第一项与节点虚位移相乘等于集中力所做的虚功