beamelementstiffness6.m

有限元程序的matlab工具包

function y = BeamElementStiffness6(E,G,A,Ix,Iy,Iz,L)

k=zeros(12,12);
k(1,1)=E*A/L;
k(2,2)=12*E*Iz/L^3;
k(3,3)=12*E*Iy/L^3;
k(4,4)=G*Ix/L;
k(5,3)=-6*E*Iy/L^2;
k(5,5)=4*E*Iy/L;
k(6,2)=6*E*Iz/L^2;
k(6,6)=4*E*Iz/L;
k(7,1)=-E*A/L;
k(7,7)=E*A/L;
k(8,2)=-12*E*Iz/L^3;
k(8,6)=-6*E*Iz/L^2;
k(8,8)=12*E*Iz/L^3;
k(9,3)=-12*E*Iy/L^3;
k(9,5)=6*E*Iy/L^2;
k(9,9)=12*E*Iy/L^3;
k(10,4)=-G*Ix/L;
k(10,10)=G*Ix/L;
k(11,3)=-6*E*Iy/L^2;
k(11,5)=2*E*Iy/L;
k(11,9)=6*E*Iy/L^2;
k(11,11)=4*E*Iy/L;
k(12,2)=6*E*Iy/L^2;
k(12,6)=2*E*Iz/L;
k(12,8)=-6*E*Iz/L^2;
k(12,12)=4*E*Iz/L;
k(3,5)=k(5,3);
k(2,6)=k(6,2);
k(1,7)=k(7,1);
k(2,8)=k(8,2);
k(6,8)=k(8,6);
k(3,9)=k(9,3);
k(5,9)=k(9,5);
k(4,10)=k(10,4);
k(3,11)=k(11,3);
k(5,11)=k(11,5);
k(9,11)=k(11,9);
k(2,12)=k(12,2);
k(6,12)=k(12,6);
k(8,12)=k(12,8);
y=k;