平面问题的有限元法6.pptVIP

图 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 {?}e