MATLAB求解总体刚度矩阵!10
clear
clc
E=2*10^11;
J=1000;
L=5;
NE=3; NP=4;
K=zeros(2*NP,2*NP);
ElementInf=[1 1 2 E J L
2 2 3 E J L
3 3 4 E J L];
For e=1:1:NE
Ee=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);end
BC=[1 1 0
2 0 0
3 0 0
4 1 1];
For I=1:NP
If BC(I,2)==1
K(2*I-1,:)=0; K(:,2*I-1)=0;
K(2*I-1,2*I-1) =1;
End
If BC(I,3)==1
K(2*I,:)=0; K(:,2*I)=0;
K(2*I,2*I) =1;
End
End
k2=k1;
k3=k1;
for i=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);
end
KK=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(8,8)=1;
Q=[0;2;3;4;5;6;0;0];
delta=KK\Q
Function K=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(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
本文档为【MATLAB求解总体刚度矩阵!10】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。