第三节结构刚度矩阵的二维等带宽存储.ppt

第三节 结构刚度矩阵的二维等带宽存储 一.结构总刚的最大半带宽NW 如图(a)所示 3.结构总刚的最大半带宽NW的影响因素 4. NW的计算公式和计算框图 ⅱ.NW的计算框图 问题 (2)在先处理法中 实际上，在程序中采用的NW计算公式为 ⅱ.NW的计算框图（图3-11） 例: 二.后处理法形成二维等带宽存储的结构刚度矩阵 例: 习题3-3 (二)后处理法形成等带宽存储的结构原始刚度矩阵[K*]N×NW的框图（图3-13） (三)引入结构的支承条件 (2)对于方阵[K]， (3)对于矩阵[K*]，应将第L行第1个元素KL1*改 为确定450斜线上元素的最大列码Jm，可分两种情况考虑： (4)对于结点荷载{P} 2.乘大数法 三. 先处理法形成二维等带宽存储的结构刚度矩阵 例： 解： 习题3-4 先处理法形成等带宽存储总刚的框图 * * NW N (a) 对称 N NW (b) 上一节 返回上级菜单 在程序PFL和PFF中，总刚[K] 是采用满阵存储，它需要占用计算机的大量内存资源，对于大型结构，满阵存储是不可取的。 为了节省计算机的内存资源，根据总刚的对称性和稀疏性，我们引入了二维等带宽的存储方式──将总刚[K] (如图(a)所示, ×──表示非零元素)的上半带(或下半带)内的元素(包括带内零元素) 存储在一个二维数组(如图(b)所示)中。 1.某行半带宽的定义 在结构总刚[K] 的上三角部分，从某行的主对角线元素到该行最后一个非零元素为止的元素个数，称为该行的半带宽。 2.最大半带宽的定义 在总刚[K] 中各行半带宽的最大值称为结构总刚的最大半带宽，用NW表示。 NW N 对称 (a) 1、2、3、4、5、6、7、8、9、10、11、12行的半带宽分别为6、5、4、6、5、4、6、5、4、3、2、1。 显然，结构总刚的最大半带宽NW=6 ───结点编号顺序 NW=3×3=9 NW=5×3=15 NW=8×3=24 根据结构总刚元素的物理意义可知，对于给定的结构，当采用不同的结点编号顺序时，结构总刚的最大半带宽一般是不同的。 例如图(a)所示刚架，当采用图(b)、(c）、(d)三种编号顺序时，若采用后处理法，则结构原始刚度矩阵的非零子块(用“×”表示) 的分布如图(b’)、(c’）、(d’)所示。 因为每个子块都是3×3阶子矩阵，所以NW如上一张所示。 由此可见，采用图(b)的编号顺序是最好的。 (1)在后处理法中 ⅰ.NW的计算公式 NW=(各单元两端结点号之差的最大值+1)×C 式中：C为结构中一个结点的位移分量个数。 对于平面桁架，C=2； 对于平面刚架、空间桁架、交叉梁系，C=3 。 图3-10 是 NW=0 I=1,NE 初值 对单元循环 L=INT(ABS(IJ(I,1)-IJ(I,2))) L＞NW NW=L NW=3﹡(NW+1) 否 计算单元两端结点号之差 求单元两端结点号之差的最大值 计算最大半带宽 在后处理法中，如何使结构刚度矩阵的最大半带宽NW最小？ 答案： 应使各单元两端结点号之 差的最大值最小。 ⅰ.NW的计算公式 NW=max(NW(e)) NW(e)=MAX－MIN+1 (a) 式中： MAX ──单元(e)的定位向量｛λ｝(e)中的最大分量； MIN ──单元(e)的定位向量｛λ｝(e)中的最小分量 (非零值)； NW(e)──单元(e)的单刚元素在总刚[K]中产生的半带 宽度。 若 NW＜LV(I)－MIN+1 I=1,2,…,6 (b) 则 NW=LV(I)－MIN+1 初值 求[K]的最大半带宽 求单元定位向量中最小非零分量 单元杆端位移号循环 单元定位向量中最小非零分量初值 形成单元定位向量 对单元循环 NW=0 M=1,NE CALL ELV MIN=N+1 I=1,6 LV(I)＜MIN且LV(I)＞0 MIN=LV(I) I=1,6 NW＜LV(I) －MIN+1 NW=LV(I) －MIN+1 是 否 否 是 求下列结构按等带宽存储总刚[K]时所需的半带宽NW 解： (1)先处理法 (0,0,1) (2,3,4) (5,6,7) (5,6,8) (0,0,0) 将刚架的结点位移分量编号 (1) (2) (3) (2)后处理法 NW=(5-2+1)×2=8 (一)矩阵[K]N×N的上半带元素在矩阵[K*]N×NW中的位置 Krs与k*IJ的下标对应关系为： 设方阵[K]N×N中上半带内第r行、s列的元素krs存储在矩阵[K*]N×NW中元素k*IJ的位置上，则这两种