MATLAB求解总体刚度矩阵.docx

clearclcE=2*1011;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: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);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 Endk2=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);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(8,8)=1;Q=0;2;3;4;5;6;0;0;delta=KKQFunction 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