成都网站建设 四川冠辰,网页设计大作业模板,wordpress内置函数,移动端网站构成要素矩阵位移法编程大作业(091210211)一、编制原理本程序的原理是基于结构力学矩阵位移法原理#xff0c;以结构结点位移作基本未知量#xff0c;将要分析的结构拆成已知节点力—结点力位移关系的单跨梁集合#xff0c;通过强令结构发生待定的基本未知位移#xff0c;在各个单跨…矩阵位移法编程大作业(091210211)一、编制原理本程序的原理是基于结构力学矩阵位移法原理以结构结点位移作基本未知量将要分析的结构拆成已知节点力—结点力位移关系的单跨梁集合通过强令结构发生待定的基本未知位移在各个单跨梁受力分析结果的基础上通过保证结构平衡建立位移法的线性方程组从而求得基本未知量。二、程序说明本程序是计算3层11跨框架右侧结点的位移和弯矩的程序编译过程是按照矩阵位移法的先处理法进行的。首先将结构杆件的交汇点作为结点共有36个结点和108个位移编号然后根据梁、柱、斜杆的不同分别建立单元刚度矩阵然后转换为整体坐标系下的刚度矩阵然后将所有杆件的单元刚度矩阵整合成为总体刚度矩阵在进行整合时连续运用for函数最终形成108阶的总体刚度矩阵。然后通过对荷载的分析自己确定出荷载矩阵直接写进程序。这样就可以把36个结点的108个位移求得然后再利用各个单元的单元刚度矩阵和所得的位移求得单元杆件的内力。 离散化编号如下图三、算法流程离散结构化编号 确定梁、柱、斜杆在整体 坐标系下的刚度矩阵 先处理法把各个杆件的单元刚度矩阵整合成总体刚度矩阵 利用杆件单元刚度矩阵和所求位移求内力 单元分析 建立梁、柱、斜杆的局部坐标系下的单元刚度矩阵 确定综合结点荷载矩阵 建立方程求解位移输出结果 四、源代码结束 %结构力学大作业 3层11跨框架矩阵位移法编程 王贝 091210211 hinput(输入单层高h); Linput(输入单跨度L);EIcinput(输入柱子的抗弯刚度EIc); EAcinput(输入柱子的抗压刚度EAc); EIbinput(输入梁的抗弯刚度EIb); EAbinput(输入梁的抗压刚度EAb); EIoinput(输入斜杆的抗弯刚度EIo); EAoinput(输入斜杆的抗压刚度EAo); qinput(输入侧向均布荷载集度q:); T1[1,0,0,0,0,0; 0,1,0,0,0,0; 0,0,1,0,0,0; 0,0,0,1,0,0; 0,0,0,0,1,0;0,0,0,0,0,1];%角度为0°的转换矩阵 T2[0,1,0,0,0,0; -1,0,0,0,0,0; 0,0,1,0,0,0; 0,0,0,0,1,0; 0,0,0,-1,0,0;0,0,0,0,0,1];%角度为90°的转换矩阵 xatan(h/L);T[cos(x),sin(x),0,0,0,0; -sin(x),cos(x),0,0,0,0; 0,0,1,0,0,0;0,0,0,cos(x),sin(x),0; 0,0,0,-sin(x),cos(x),0;0,0,0,0,0,1];%斜杆的转换矩阵 T3T;%梁的单元刚度矩阵kb0[EAb/L 0 0 -EAb/L 0 0;0 12*EIb/(L*L*L) 6*EIb/(L*L) 0 -12*EIb/(L*L*L) 6*EIb/(L*L); 0 6*EIb/(L*L) 4*EIb/L 0 -6*EIb/(L*L) 2*EIb/L; -EAb/L 0 0 EAb/L 0 0;0 -12*EIb/(L*L*L) -6*EIb/(L*L) 0 12*EIb/(L*L*L) -6*EIb/(L*L);0 6*EIb/(L*L) 2*EIb/L 0 -6*EIb/(L*L) 4*EIb/L]; %柱子的单元刚度矩阵kc0[EAc/h 0 0 -EAc/h 0 0;0 12*EIc/(h*h*h) 6*EIc/(h*h) 0 -12*EIc/(h*h*h) 6*EIc/(h*h); 0 6*EIc/(h*h) 4*EIc/h 0 -6*EIc/(h*h) 2*EIc/h; -EAc/h 0 0 EAc/h 0 0;0 -12*EIc/(h*h*h) -6*EIc/(h*h) 0 12*EIc/(h*h*h) -6*EIc/(h*h); 0 6*EIc/(h*h) 2*EIc/h 0 -6*EIc/(h*h) 4*EIc/h;]; %斜杆的单元刚度矩阵 Hsqrt(h*hL*L);ko0[EAo/H 0 0 -EAo/H 0 0;0 12*EIo/(H*H*H) 6*EIo/(H*H) 0 -12*EIo/(H*H*H) 6*EIo/(H*H); 0 6*EIo/(H*H) 4*EIo/H 0 -6*EIo/(H*H) 2*EIo/H; -EAo/H 0 0 EAo/H 0 0;0 -12*EIo/(H*H*H) -6*EIo/(H*H) 0 12*EIo/(H*H*H) -6*EIo/(H*H); 0 6*EIo/(H*H) 2*EIo/H 0 -6*EIo/(H*H) 4*EIo/H]; kbT1*kb0*T1;%总体坐标下梁的单元刚度矩阵 kcT2*kc0*T2;%总体坐标下柱子的单元刚度矩阵 koT3*ko0*T3;%总体坐标斜杆的单元刚度矩阵Xzeros(108,108);Yzeros(108,108);Zzeros(108,108);%定义108阶0矩阵 K1zeros(108,108);K2zeros(108,108);K3zeros(108,108); K4zeros(108,108);K5zeros(108,108);K6zeros(108,108); K7zeros(108,108);K8zeros(108,108);K9zeros(108,108); %把梁杆单元矩阵整合到总体刚度矩阵的循环语句 for ii1:11X(3*ii-2:3*ii3,3*ii-2:3*ii3)kb; K1K1X;Xzeros(108,108); endfor ii13:23Y(3*ii-2:3*ii3,3*ii-2:3*ii3)kb; K1K1Y;Yzeros(108,108); endfor ii25:35Z(3*ii-2:3*ii3,3*ii-2:3*ii3)kb; K1K1Z;Zzeros(108,108); end%把柱杆单元矩阵整合到总体刚度矩阵的循环语句 for jj1:36K2(3*jj-2:3*jj,3*jj-2:3*jj)kc(4:6,4:6); endfor jj1:24K3(3*jj-2:3*jj,3*jj-2:3*jj)kc(1:3,1:3); endfor jj1:24K4(3*jj-2:3*jj,3*jj34:3*jj36)kc(1:3,4:6); endfor jj1:24K5(3*jj34:3*jj36,3*jj-2:3*jj)kc(4:6,1:3); end%把斜杆杆单元矩阵整合到总体刚度矩阵的循环语句 for gg3:12:27K6(3*gg-2:3*gg,3*gg-2:3*gg)ko(4:6,4:6); endfor gg2:12:14K7(3*gg-2:3*gg,3*gg-2:3*gg)ko(1:3,1:3); endfor gg2:12:14K8(3*gg-2:3*gg,3*gg37:3*gg39)ko(1:3,4:6); endfor gg2:12:14K9(3*gg37:3*gg39,3*gg-2:3*gg)ko(4:6,1:3); endKK1K2K3K4K5K6K7K8K9;%总体刚度矩阵Pzeros(108,1);P(1,1)h*q;P(37,1)h*q;P(73,1)h*q/2;P(75,1)q*h*h/12; AK\\P;%结构位移 B1kb*A(103:108,1); B2kb*A(67:72,1); B3kb*A(31:36,1);D1zeros(6,1);D1(1:3,1)A(70:72,1);D1(4:6,1)A(106:108,1); D2zeros(6,1);D2(1:3,1)A(34:36,1);D2(4:6,1)A(70:72,1); C1kc*D1;C2kc*D2;C3kc(4:6,4:6)*A(34:36,1);M1(1,1)B1(6,1);M2(1,1)C1(6,1);M3(1,1)B2(6,1);M4(1,1)C2(6,1); M5(1,1)C1(3,1);M6(1,1)B3(6,1);M7(1,1)C3(3,1);M8(1,1)C2(3,1);for i1:3m36*i-2;fprintf(第%d层最右侧节点的位移是%d\\n,i,A(m,1)) endfprintf(第1层最右侧节点左侧杆的弯矩是%f\\n,M6) fprintf(第1层最右侧节点下侧杆的弯矩是%f\\n,M7) fprintf(第1层最右侧节点上侧杆的弯矩是%f\\n,M8) fprintf(第2层最右侧节点左侧杆的弯矩是%f\\n,M3) fprintf(第2层最右侧节点下侧杆的弯矩是%f\\n,M4) fprintf(第2层最右侧节点上侧杆的弯矩是%f\\n,M5) fprintf(第3层最右侧节点左侧杆的弯矩是%f\\n,M1) fprintf(第3层最右侧节点下侧杆的弯矩是%f\\n,M2)五、试算算例
