第四节结构刚度矩阵的一维变带宽存储.ppt

第四节 结构刚度矩阵的一维变带宽存储 以平面刚架（图a）为例 一. 主元地址数组KD(N) （二）用途 一维数组TK中的地址JI的计算公式: （三）计算机自动生成主元地址数组KD(N) 2.利用单元定位向量求矩阵[K]的各列半带宽 单刚[K](2)的元素在总刚[K]中的位置，如图（b）所示→ (1)总刚[K]由于单独叠加单刚[K](e) 而产生的第J列半带宽的计算公式 例： 3.累加各列半带宽形成主元地址数组 4.形成主元地址数组KD(N)的程序框图（图3-20） 与程序框图（图3-20）相应的子程序ADE 二. 先处理法形成一维存储的结构刚度矩阵TK(NN) 2.先处理法形成一维存储的总刚TK(NN)的程序框图 与程序框图（图3-21）相应的子程序TSM 例： 解： 习题3-5 * * 为了比二维等带宽存储方式更有效地节省 计算机内存，按照对称带状的结构刚度矩阵[K] 各列（或行）的实际半带宽，将 [K]的上半带 内元素一列接一列（或将 [K]的下半带内元素 一行接一行）地存储在一个一维数组TK(NN)中， 而不存储各列（或行）的带外零元素。这种存 储方式称为结构刚度矩阵的一维变带宽存储。 上一节 返回上级菜单 本节结合先处理法程序PFF，介绍将结构总刚[K]中上半带变带宽内的元素按列一维存储的方式和相应的程序框图。 由节点位移统一编号知结构总刚[K]为10×10阶方阵，由总刚[K]的组成规律形成总刚，如图（b）所示， 0──零元素（无单刚元素送入）； ×──非零元素（有单刚元素送入） 。 将结构总刚[K]的上半带变带宽按列一维存储在TK(NN)数组中，图（C）中数字为一维数组TK(NN)中各元素的地址。 用计算机确定[K]的上半带内任一元素在一维数组TK(NN) 中的地址时，需首先引入主元地址数组。 主元──指总刚[K]N×N中每列的主对角线元素。 主元地址──指主元在一维变带宽存储数组TK(NN)中的序号。 主元地址数组──指存放主元地址的数组，用 KD(N)表示。其中KD(J)是主元kjj的地址。 (一)定义 KD(N)={1，3，6，10，14，19，25，29，34，40} 如图（C）所示。 例如：对于图（b） （2）根据主元地址数组KD(N)， 确定[K]上半带的非对角元素kij(ij)在一维数组TK中的地址JI。 分析图(b)可以看出： （1）主元地址数组KD(N)的最后一个元素表示了一维数组TK中存储的元素个数NN。 如NN=KD(10)=40 例如图(b)中k67在一维数组TK中的地址为： 依此类推，得： 1. 列半带宽（或列高）的定义 在结构刚度矩阵[K]的上三角部分，从某列的第一个非零元素起到该列主对角元素为止所具有的元素个数称为该列的半带宽（或称为列高）。 例如：5列的列高为4， 6列的列高为5。 以图(a)所示结构的单元(2)为例， 该单元定位向量为 {λ}(2)={5,6,7,9,0,10} 其中最小非零分量为5。 (1) (2) (3) (4) 总刚[K]由于单独叠 加单刚[K](2) 的元素 所产生的各列半带宽aj为： a5=5-5+1=1 a6=6-5+1=2 a7=7-5+1=3 a9=9-5+1=5 a10=10-5+1=6 由以上计算可归纳出： 其中 MIN──单元(e)定位向量{λ} (e)中的最小非零分量； J ──单元(e)定位向量{λ} (e)中的某个非零分量。 (2)总刚[K]的第J列真正半带宽： （1） （2） （3） （4） 对于图(a)所示结构，按式（3-12） 利用各单元定位向量求总刚[K]的各列半带宽。 解： 否 否 M=1,NE 数组KD置零 J=0 是 CALL ELV MIN=N J=LV(I) NW=J-MIN+1 KD(J)=NN NN=NN+KD(J) KD(J)=NW NN=1 J=LV(I) MIN=J I=1,6 J=2,N I=1,6 NW≤KD(J) J=0 J≥MIN 是 是 是 否 否 对单元循环 形成单元定位向量初值 对定位向量的六个分量循环 取定位向量的分量 排除零分量 求定位向量的最小零分量 对定位向量的六个分量循环 取定位向量的分量 排除零分量 列半带宽 第J列的真正半带宽 初值 总刚[K]的列循环 累加各列半带宽， 形成KD数组 见128-129页（自学） 1.步骤 （1）利用单元定位向量确定单元刚度矩阵[K](e)的第L行和第K 列元素在结构刚度矩阵[K]N×N 中的行码I和列码J； （2）按式（3-11）求得总刚[K]N×N的元素kij在一维数组TK(NN)中的地址JI，JI也是单刚元素EK(L,K