MATLAB求解总体刚度矩阵
《MATLAB求解总体刚度矩阵》由会员分享,可在线阅读,更多相关《MATLAB求解总体刚度矩阵(5页珍藏版)》请在装配图网上搜索。
1、clearclcE=2*1OT1;J=1000;L=5;NE=3;NP=4;K=zeros(2*NP,2*NP);ElementInf=112EJL2 23EJL3 34EJL;Fore=1:1:NEEe=ElementInf(e,4);Je=ElementInf(e,5);Le=ElementInf(e,6);ie=ElementInf(e,2);je=ElementInf(e,3);ke=Element_stiffness_matrix(Ee,Je,Le);K=BeamZongGang(K,ke,ie,je);endBC=110200300411;ForI=1:NPIfBC(I,2)=1K
2、(2*I-1,:)=0;K(:,2*I-1)=0;K(2*I-1,2*I-1)=1;EndIfBC(I,3)=1K(2*I,:)=0;K(:,2*I)=0;K(2*I,2*I)=1;EndEndk2=k1;k3=k1;fori=1:1:4;j=1:1:4;K(i,j)=K(i,j)+k1(i,j);K(i+2,j+2)=K(i+2,j+2)+k2(i,j);K(i+4,j+4)=K(i+4,j+4)+k3(i,j);endKK=K;KK(1,:)=0;KK(:,1)=0;KK(1,1)=1;KK(7,:)=0;KK(:,7)=0;KK(7,7)=1;KK(8,:)=0;KK(:,8)=0;KK
3、(8,8)=1;Q=0;2;3;4;5;6;0;0;delta=KKQFunctionK=BeamZongGang(K,ie,je)K(2*ie-1,2*ie-1)=K(2*ie-1,2*ie-1)+ke(1,1);K(2*ie-1,2*ie)=K(2*ie-1,2*ie)+ke(1,2);K(2*ie-1,2*je-1)=K(2*ie-1,2*je-1)+ke(1,3);K(2*ie-1,2*je)=K(2*ie-1,2*je)+ke(1,4);K(2*ie,2*ie-1)=K(2*ie,2*ie-1)+ke(2,1);K(2*ie,2*ie)=K(2*ie,2*ie)+ke(2,2);K(
4、2*ie,2*je-1)=K(2*ie,2*je-1)+ke(2,3);K(2*ie,2*je)=K(2*ie,2*je)+ke(2,4);K(2*je-1,2*ie-1)=K(2*je-1,2*ie-1)+ke(3,1);K(2*je-1,2*ie)=K(2*je-1,2*ie)+ke(3,2);K(2*je-1,2*je-1)=K(2*je-1,2*je-1)+ke(3,3);K(2*je-1,2*je)=K(2*je-1,2*je)+ke(3,4);K(2*je,2*ie-1)=K(2*je,2*ie-1)+ke(4,1);K(2*je,2*ie)=K(2*je,2*ie)+ke(4,2);K(2*je,2*je-1)=K(2*je,2*je-1)+ke(4,3);K(2*je,2*je)=K(2*je,2*je)+ke(4,4);K(2*i-1,2*i,2*j-1,2*j,2*i-1,2*i,2*j-1,2*j)=K(2*i-1,2*i,2*j-1,2*j,2*i-1,2*i,2*j-1,2*j)+ke
- 温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。