矩阵位移法连续梁计算程序分析报告.docxVIP

结构力学矩阵位移法连续梁计算程序——分析报告王识水工13班 2011010284程序框架程序源代码及注释*********************************************************************“main.m”function main()%各种输入n = input(请输入杆的数量n，必须为整数\n); l = input1(请输入各杆的杆长l,用向量表示，如[2 3 4 5],长度为n\n,n);EI = input1(请输入各杆的抗弯刚度EI,用向量表示，如[2 3 4 5],长度为n\n,n);link = input1(请输入两端的连接情况，用向量表示，如[1 0],0为铰接，1为刚接，长度为2\n,2);f = input1(请输入各杆的集中力荷载，用向量表示，如[2 0 4 0],长度为n\n,n);q = input1(请输入各杆的均布力荷载，用向量表示，如[2 0 4 0],长度为n\n,n);M = input1(请输入各节点的集中力矩，用向量表示，如[2 0 4 0 4],长度为n+1\n,n+1);i = zeros(1,n);for k = 1:n(1); i(k) = EI(k)./l(k);%各杆的刚度end K = Ki(n,i,link);%求出整体刚度矩阵Fp = Fi(n,l,f,q);%求出荷载