第4章 结构刚度矩阵存贮和组集.pptVIP

对于大型结构，应该特别注意选取恰当的结点编号方式。这是因为总存贮量大体上与带宽成正比的原故（对等带宽存贮，则准确地与最大半带宽成正比）。 4、程序中应通过试算验证 编好的程序，必须经过试算，以检查刚度矩阵是否被正确地形成。通常，可采用简单但具有典型性的算例进行计试算，输出形成后的刚度矩阵。也可以在组集刚度矩阵的程序段内设置一些输出语句，输出少量信息以便检查。导致刚度矩阵未被正确形成的因素是多方面的，可能是子程序编写有错，也可能计算程序完全正确而输入数据本身有错，例如支承约束条件填写有错等等。 总之，计算机铁面无情，不允许在任何环节上发生任何差错。另一方面，它也是很公正的，对于所有认真而细心的使用者，都迅速准确地提供所要求的计算结果。 4.2 结构刚度矩阵的等带宽存贮 1、结构刚度矩阵的稀疏性 在第3章里，我们曾经指出，结构刚度矩阵是具有大量零元素的稀疏矩阵。 现在对这一性质作进一步的考察。 在对结构系统进行平衡分析时知，如果i节点发生单位位移，只有与节点i有直接联系的那些节点才会受到影响，产生节点力。反之， i节点是否产生节点力也只会受与它直接相连的节点的影响。 设结构中共有100个节点，其中第i（i = 4）号节点及其与i号节点直接联系的节点如图4-4所示。 图4-4 平面框架中的节点 2 5 j=7 i=4 则结构刚度矩阵中的3i-2, 3i-1, 3i 行的形式为： 2 10 11 12 1 ? ? ? 100 10 21 12 上图中每个小方块是一个三阶子块阵，所示出的与节点i相应的100个三阶子块中，只有4个非零块，其余96个皆为零元素组成的三阶子块。通常，对于任何一个节点，与它直接联系的节点个数和问题中的节点总数比较起来总是小得多。因而，呈现出结构刚度矩阵的高度稀疏性。 图4-5 4 由于刚度矩阵为对称矩阵，只要对节点适当编号，就能使刚度矩阵的非零元素聚集在主对角线两侧，呈现为图4-6所示带状矩阵。 UBW UBW 0 0 对 称 图4-6 一般，刚度矩阵各行的非零元素个数和排列是不同的。对第i行而言，从主对角线元素Kii算起，直到它最右边的一个非零元素为止，其间所有元素（包括可能有的零元素）的个数，称为该行的右带宽。以图4-5为例，第10行的右带宽为12，即包括 K10，10，K10，11，…，K10，21。 这12个元素。 把相互连接的两个节点i和j的号码差∣j – i∣称为相邻点号差，当ij时从主对角线子块[Kii]到最右边的非零块[Kij] ，共有∣j – i∣+ 1个三阶子块；当ij时，从[Kjj]到[Kji]也共有∣j – i∣+ 1个三阶子块。因而： 图中所示三行中最大的右带宽为： （∣j – i∣+ 1 ）* 3 2 10 11 12 1 ? ? ? 100 10 21 12 图4-7 4 刚度矩阵各行右带宽中之最大者，称为最大右带宽。因刚度矩阵为对称矩阵，最大右带宽和最大左带宽是相同的，故简称最大右带宽为最大半带宽，记为UBW（Upper Band Width ）。 当结构的节点编号确定后，刚度矩阵的最大半带宽随之确定。找出所有相邻点号差中的最大者，记为rmax，则刚度矩阵最大半带宽为 UBW =（1 + rmax ）* NDF （4-9） 其中，NDF为节点自由度。对平面刚架，NDF = 3。对弹性平面问题，NDF=2。 2、结构刚度矩阵的等带宽存贮 鉴于工程结构都是几何不变结构，从现在起，我们只讨论受约束结构的刚度矩阵。为书写简便起见，用[K]代表受约束结构的刚度矩阵，即以前的[Kff]。 因为刚度矩阵是对称的，只须存贮它的一半即可。约定：只存贮[K]的上三角部分。带区之外全是零元素，自然不必存贮，因此只需存贮包括主对角元素在内的上半带宽区域内的元素。这样，大大节约了需占用的计算机内存。称这种存贮方式为等带宽存贮。 设N为独立结点位移分量总数，则原始形式的刚度矩阵[K]是一个N行N列的方阵。所谓等带宽存贮，就是把[K]的上半带宽区域内全部元素用一个N行UBW列的矩形数组予以存贮： UBW UBW 行 号 1 → IR → N→ 1 列 号 JC 行 号 1 → IR → N→ 1 JC-(IR-1) 方阵形式 等带宽形式 图4-8 对比上图知： （1）等带宽存贮实际上是矩形数组存贮，矩形数组的行数和方阵行数