利用matlab计算多层框架位移源代码.docVIP

下载本文档

60
0
约3.82千字
约 4页
2018-05-23 发布于河南
举报
版权申诉

利用matlab计算多层框架位移源代码.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

利用matlab计算多层框架位移源代码

源代码 091210322 蒋伟 function[y,M] =bent(h,L,EIc,EAc,EIb,EAb,EIo,EAo,q) %变量输入% h=input(单层高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=); %柱子的刚度% a1=EAc/h; a2=12*EIc/(h^3); a3=6*EIc/(h^2); a4=4*EIc/h; a5=2*EIc/h; ka=[a1 0 0 -a1 0 0 0 a2 a3 0 -a2 a3 0 a3 a4 0 -a3 a5 -a1 0 0 a1 0 0 0 -a2 -a3 0 a2 -a3 0 a3 a5 0 -a3 a4]; %梁的刚度% b1=EAb/L; b2=12*EIb/(L^3); b3=6*EIb/(L^2); b4=4*EIb/L; b5=2*EIb/L; kb=[b1 0 0 -b1 0 0 0 b2 b3 0 -b2 b3 0 b3 b4 0 -b3 b5 -b1 0 0 b1 0 0 0 -b2 -b3 0 b2 -b3 0 b3 b5 0 -b3 b4]; %斜杆的刚度% C=sqrt(h^2+L^2); c1=EAo/C; c2=12*EIo/(C^3); c3=6*EIo/(C^2); c4=4*EIo/C; c5=2*EIo/C; kc=[c1 0 0 -c1 0 0 0 c2 c3 0 -c2 c3 0 c3 c4 0 -c3 c5 -c1 0 0 c1 0 0 0 -c2 -c3 0 c2 -c3 0 c3 c5 0 -c3 c4]; %坐标转换% T=[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]; d=L/sqrt(h^2+L^2); D=h/sqrt(h^2+L^2); t=[d D 0 0 0 0 -D d 0 0 0 0 0 0 1 0 0 0 0 0 0 d D 0 0 0 0 -D d 0 0 0 0 0 0 1]; Ka=T*ka*T; Kb=kb; Kc=t*kc*t; %柱的整合% K1=zeros(165,165); K2=zeros(165,165); K3=zeros(165,165); K4=zeros(165,165); K5=zeros(165,165); for i=0:43 K2((1+3*i):(3+3*i),(1+3*i):(3+3*i))=Ka(1:3,1:3); K3((34+3*i):(36+3*i),(34+3*i):(36+3*i))=Ka(4:6,4:6); K4((34+3*i):(36+3*i),(1+3*i):(3+3*i))=Ka(4:6,1:3); K5((1+3*i):(3+3*i),(34+3*i):(36+3*i))=Ka(1:3,4:6); K1=K1+K2+K3+K4+K5; K2=zeros(165,165); K3=zeros(165,165); K4=zeros(165,165); K5=zeros(165,165); end %梁的整合% K2=zeros(165,165); K3=zeros(165,165); K4=zeros(165,165); K5=zeros(165,165); K6=zeros(165,165); for i=12:54 if(rem(i,11)==0) continue else K3((3*i-2):3*i,(3*i-2):3*i)=Kb(1:3,1:3); K4((3*i-2):3*i,(3*i+1):(3*i+3))=Kb(1:3,4:6); K5((3*i+1):(3*i+3),(3*i-2):3