matlab编程四层七跨刚架的计算.doc

矩阵位移法编程大作业 根据题目的要求，我要编的程序是处理四层七跨且第二跨带有斜杆的刚架结构。 编程原理 处理这种刚架的问题，主要采用的是矩阵位移法。由单刚组建成整个刚架的总刚利用公式K*Δ=F求解位移，在利用位移求解力。 程序说明 程序仅限于四层七跨，第二跨带斜杆的刚架。先输入所需变量单层高h、单跨度l、柱子抗弯刚度EIc、柱子抗压刚度EAc、梁的抗弯刚度EIb、梁的抗压刚度EAb、斜杆抗弯刚度EIo、斜杆抗压刚度EAo、横向的均布荷载q，最后运行程序后将得出每一层最右端节点的侧向位移和弯矩值。 算法流程 首先对每个节点进行编号如上图，每个节点都对应三个位移并编号，以及对杆件的编号。接着建立柱子、梁和斜杆的单刚K1e、K2、K3e，由于柱子和斜杆都是有转角的所以需要转化分别建立转化矩阵T1、T2最后得出他们的单刚为K1、K2、K3，然后对单刚进行整合。然后是将每个杆件的刚度整合到总刚里，先建立一个120阶的零矩阵。柱子的整合是按照编号找到规律，利用每个柱子的两个节点差值为8的特点，以柱子的下端为基准从节点1到32，将柱子的单刚分成四部分用循环有规律的对应编号的位置放入总刚中，有的位置还需要叠加。梁的整合是根据梁的两端节点值相差1建立的规律，以梁的左端为基准节点从1到40（除去8的倍数，因为8的倍数的节点在刚架的最右端没有以他们为基准存在的梁），将梁的单刚分成四个部分用循环有规律的对应编号的位置放入总刚，有的位置还需要叠加。斜杆是根据两端节点值相差9的规律，一斜杆的左端为基准节点为2、10、18、26，将斜杆的单刚分成四个部分用循环有规律的对应编号的位置放入总刚，有的位置还需要叠加。建立好总刚后，由于节点1到8是固定端所以采用前处理法，先将总刚中的k(1:24,1:120)和k(1:120,1:24)变成零矩阵，再将k(1:24,1:24)变为单位矩阵。接着剪力外力矩阵，将横向荷载化成等效节点荷载对应相应的位移放在相应的位置其余位置为零荷载，因为采用了前处理法所以荷载矩阵的前24个元素一定为零。接着求解对应的位移阵，取出其中要求的位置的位移，再根据要求解除最右端节点对应不同杆件所产生的弯矩。 程序的源代码 %091210225张刚 四层七跨 第二跨有斜杆的桁架计算程序 h=input(单层高h\n); l=input(单跨度l\n); E=input(弹性模量E\n); EIc=input(柱子抗弯刚度EIc\n); EAc=input(柱子抗压刚度EAc\n); EIb=input(梁的抗弯刚度EIb\n); EAb=input(梁的抗压刚度EAb\n); EIo=input(斜杆抗弯刚度EIo\n); EAo=input(斜杆抗压刚度EAo\n); q=input(横向的均布荷载q\n); %柱子的刚度 K1e=[EAc/h,0,0,-EAc/h,0,0; 0,12*EIc/h^3,6*EIc/h^2,0,-12*EIc/h^3,6*EIc/h^2; 0,6*EIc/h^2,4*EIc/h,0,-6*EIc/h^2,2*EIc/h; -EAc/h,0,0,EAc/h,0,0; 0,-12*EIc/h^3,-6*EIc/h^2,0,12*EIc/h^3,-6*EIc/h^2; 0,6*EIc/h^2,2*EIc/h,0,-6*EIc/h^2,4*EIc/h]; n1=[0,1,0;-1,0,0;0,0,1]; n2=zeros(size(n1)); T1=[n1,n2;n2,n1]; K1=T1.*K1e*T1; %梁的刚度 K2=[EAb/l,0,0,-EAb/l,0,0; 0,12*EIb/l^3,6*EIb/l^2,0,-12*EIb/l^3,6*EIb/l^2; 0,6*EIb/l^2,4*EIb/l,0,-6*EIb/l^2,2*EIb/l; -EAb/l,0,0,EAb/l,0,0; 0,-12*EIb/l^3,-6*EIb/l^2,0,12*EIb/l^3,-6*EIb/l^2; 0,6*EIb/l^2,2*EIb/l,0,-6*EIb/l^2,4*EIb/l]; %斜杆的刚度 m=sqrt(h*h+l*l); K3e=[EAo/m,0,0,-EAo/m,0,0; 0,12*EIo/m^3,6*EIo/m^2,0,-12*EIo/m^3,6*EIo/m^2; 0,6*EIo/m^2,4*EIo/m,0,-6*EIo/m^2,2*EIo/m; -EAo/m,0,0,EAo/m,0,0; 0,-12*EIo/m^3,-6*EIo/m^2,0,12*EIo/m^3,-6*EIo/m^2; 0,6*EIo/m^2,2*EIo/m,0,-6*EIo/m^2,4*EIo/m]; a1=[l/m,h/m,0;-