版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、矩阵位移法编程大作业(091210211)一、编制原理本程序的原理是基于结构力学矩阵位移法原理,以结构结点位移作基本未知量,将要分析的结构拆成已知节点力结点力位移关系的单跨梁集合,通过强令结构发生待定的基本未知位移,在各个单跨梁受力分析结果的基础上通过保证结构平衡建立位移法的线性方程组,从而求得基本未知量。二、程序说明本程序是计算3层11跨框架右侧结点的位移和弯矩的程序,编译过程是按照矩阵位移法的先处理法进行的。首先将结构杆件的交汇点作为结点,共有36个结点和108个位移编号,然后根据梁、柱、斜杆的不同分别建立单元刚度矩阵,然后转换为整体坐标系下的刚度矩阵,然后将所有杆件的单元刚度矩阵整合成为
2、总体刚度矩阵,在进行整合时连续运用for函数,最终形成108阶的总体刚度矩阵。然后通过对荷载的分析自己确定出荷载矩阵,直接写进程序。这样就可以把36个结点的108个位移求得,然后再利用各个单元的单元刚度矩阵和所得的位移求得单元杆件的内力。离散化编号如下图:三、算法流程利用杆件单元刚度矩阵和所求位移求内力建立方程,求解位移确定综合结点荷载矩阵先处理法把各个杆件的单元刚度矩阵整合成总体刚度矩阵建立梁、柱、斜杆的局部坐标系下的单元刚度矩阵确定梁、柱、斜杆在整体坐标系下的刚度矩阵单元分析结构离散化编号结束输出结果四、源代码%结构力学大作业 3层11跨框架矩阵位移法编程 王贝 091210211h=in
3、put('输入单层高h:');L=input('输入单跨度L:');EIc=input('输入柱子的抗弯刚度EIc:');EAc=input('输入柱子的抗压刚度EAc:');EIb=input('输入梁的抗弯刚度EIb:');EAb=input('输入梁的抗压刚度EAb:');EIo=input('输入斜杆的抗弯刚度EIo:');EAo=input('输入斜杆的抗压刚度EAo:');q=input('输入侧向均布荷载集度q:');T1=1,0,0,0
4、,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°的转换矩阵x=atan(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(
5、x),0; 0,0,0,0,0,1;%斜杆的转换矩阵T3=T;%梁的单元刚度矩阵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
6、;%柱子的单元刚度矩阵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;%斜杆的单元刚度矩阵H=sqrt(h*h+L*L);ko
7、0=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;kb=T1'*kb0*T1;%总体坐标下梁的单元刚度矩阵kc=T2'
8、;*kc0*T2;%总体坐标下柱子的单元刚度矩阵ko=T3'*ko0*T3;%总体坐标斜杆的单元刚度矩阵X=zeros(108,108);Y=zeros(108,108);Z=zeros(108,108);%定义108阶0矩阵K1=zeros(108,108);K2=zeros(108,108);K3=zeros(108,108);K4=zeros(108,108);K5=zeros(108,108);K6=zeros(108,108);K7=zeros(108,108);K8=zeros(108,108);K9=zeros(108,108);%把梁杆单元矩阵整合到总体刚度矩阵的循环语
9、句for ii=1:11 X(3*ii-2:3*ii+3,3*ii-2:3*ii+3)=kb; K1=K1+X;X=zeros(108,108);endfor ii=13:23 Y(3*ii-2:3*ii+3,3*ii-2:3*ii+3)=kb; K1=K1+Y;Y=zeros(108,108);endfor ii=25:35 Z(3*ii-2:3*ii+3,3*ii-2:3*ii+3)=kb; K1=K1+Z;Z=zeros(108,108);end%把柱杆单元矩阵整合到总体刚度矩阵的循环语句for jj=1:36 K2(3*jj-2:3*jj,3*jj-2:3*jj)=kc(4:6,4:6
10、);endfor jj=1:24 K3(3*jj-2:3*jj,3*jj-2:3*jj)=kc(1:3,1:3);endfor jj=1:24 K4(3*jj-2:3*jj,3*jj+34:3*jj+36)=kc(1:3,4:6);endfor jj=1:24 K5(3*jj+34:3*jj+36,3*jj-2:3*jj)=kc(4:6,1:3);end%把斜杆杆单元矩阵整合到总体刚度矩阵的循环语句for gg=3:12:27 K6(3*gg-2:3*gg,3*gg-2:3*gg)=ko(4:6,4:6);endfor gg=2:12:14 K7(3*gg-2:3*gg,3*gg-2:3*gg
11、)=ko(1:3,1:3);endfor gg=2:12:14 K8(3*gg-2:3*gg,3*gg+37:3*gg+39)=ko(1:3,4:6);endfor gg=2:12:14 K9(3*gg+37:3*gg+39,3*gg-2:3*gg)=ko(4:6,1:3);endK=K1+K2+K3+K4+K5+K6+K7+K8+K9;%总体刚度矩阵P=zeros(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;A=KP;%结构位移B1=kb*A(103:108,1);B2=kb*A(67:72,1);B3=kb*A(
12、31:36,1);D1=zeros(6,1);D1(1:3,1)=A(70:72,1);D1(4:6,1)=A(106:108,1);D2=zeros(6,1);D2(1:3,1)=A(34:36,1);D2(4:6,1)=A(70:72,1);C1=kc*D1;C2=kc*D2;C3=kc(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
13、i=1:3 m=36*i-2;fprintf('第%d层最右侧节点的位移是%dn',i,A(m,1)endfprintf('第1层最右侧节点左侧杆的弯矩是%fn',M6)fprintf('第1层最右侧节点下侧杆的弯矩是%fn',M7)fprintf('第1层最右侧节点上侧杆的弯矩是%fn',M8)fprintf('第2层最右侧节点左侧杆的弯矩是%fn',M3)fprintf('第2层最右侧节点下侧杆的弯矩是%fn',M4)fprintf('第2层最右侧节点上侧杆的弯矩是%fn',M5)fprintf('第3层最右侧节点左侧杆的弯矩是%fn',M1)fprintf('第3层最右侧节点下侧杆的弯矩是%fn',M2)五、试算算例输入数据:输入单层高h:1输入单跨度L:1输入柱子的抗弯刚度EIc:1输入柱子的抗压刚度EAc:1输入梁的抗弯刚度EIb:1输入梁的抗压刚度EAb:1输入斜杆的抗弯刚度EIo:1输入斜杆的抗压刚度EAo:1输入侧向均布荷载集度q:1计算结果:第1层最右侧节点的位移是-6.219850e-003第2层最右侧节点的位移是-2.152659e-002第3层最右侧节点的位移是-4.1
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 窗帘装饰合同三篇
- 《评价大纲写法案例》课件
- 《证券投资模拟实训》课件
- 《信息传输技术基础》课件
- 《计算机数据安全》课件
- 函数应用题专题复习课件
- 堆放材料整改报告范文
- 《电工电子技术》课件-第9章
- 旅游市场旅游学概论课件全集之
- 2025年资阳货运从业资格仿真考题
- 2025年中考英语复习:热点话题作文与范文汇编
- 2024-2030年中国潮牌衣服行业市场发展趋势与前景展望战略分析报告
- 中国传统文化漫谈智慧树知到答案2024年西华大学
- 初中语文八年级上册 20《人民英雄永垂不朽》公开课一等奖创新教学设计
- 生猪屠宰兽医卫生检验人员理论考试题库及答案
- DL-T 2680-2023 电力建设施工企业安全生产标准化实施规范
- 液体容积式油流量计在线校准规范
- 互联网文案写作智慧树知到期末考试答案章节答案2024年山东管理学院
- 51个行业领域重大事故隐患判定标准和重点检查事项汇编
- 江苏开放大学本科行政管理专业060191西方行政学说期末试卷
- 服装厂生产制度
评论
0/150
提交评论