桁架杆有限元计算.rar_有限元 杆_杆 桁架_杆有限元_桁架_桁架 刚度矩阵

function y=SpaceTrussAssemble(K,k,i,j) K(3*i-2,3*i-2)=K(3*i-2,3*i-2)+k(1,1) K(3*i-2,3*i-1)=K(3*i-2,3*i-1)+k(1,2) K(3*i-2,3*i)=K(3*i-2,3*i)+k(1,3) K(3*i-2,3*j-2)=K(3*i-2,3*j-2)+k(1,4) K(3*i-2,3*j-1)=K(3*i-2,3*j-1)+k(1,5) K(3*i-2,3*j)=K(3*i-2,3*j)+k(1,6) K(3*i-1,3*i-2)=K(3*j-2,3*i-2)+k(2,1) K(3*i-1,3*i-1)=K(3*j-2,3*i-1)+k(2,2) K(3*i-1,3*i)=K(3*j-2,3*i)+k(2,3) K(3*i-1,3*j-2)=K(3*j-2,3*j-2)+k(2,4) K(3*i-1,3*j-1)=K(3*j-2,3*j-1)+k(2,5) K(3*i-1,3*j)=K(3*j-2,3*j)+k(2,6) K(3*i,3*i-2)=K(3*i-2,3*i-2)+k(3,1) K(3*i,3*i-1)=K(3*i-2,3*i-1)+k(3,2) K(3*i,3*i)=K(3*i-2,3*i)+k(3,3) K(3*i,3*j-2)=K(3*i-2,3*j-2)+k(3,4) K(3*i-2,3*j-1)=K(3*i-2,3*j-1)+k(3,5) K(3*i-2,3*j)=K(3*i-2,3*j)+k(3,6) K(3*j-2,3*i-2)=K(3*j-2,3*i-2)+k(4,1) K(3*j-2,3*i-1)=K(3*j-2,3*i-1)+k(4,2) K(3*j-2,3*i)=K(3*j-2,3*i)+k(4,3) K(3*j-2,3*j-2)=K(3*j-2,3*j-2)+k(4,4) K(3*j-2,3*j-1)=K(3*j-2,3*j-1)+k(4,5) K(3*j-2,3*j)=K(3*j-2,3*j)+k(4,6) K(3*j-1,3*i-2)=K(3*j-1,3*i-2)+k(5,1) K(3*j-1,3*i-1)=K(3*j-1,3*i-1)+k(5,2) K(3*j-1,3*i)=K(3*j-1,3*i)+k(5,3) K(3*j-1,3*j-2)=K(3*j-1,3*j-2)+k(5,4) K(3*j-1,3*j-1)=K(3*j-1,3*j-1)+k(5,5) K(3*j-1,3*j)=K(3*j-1,3*j)+k(5,6) K(3*j,3*i-2)=K(3*j,3*i-2)+k(6,1) K(3*j,3*i-1)=K(3*j,3*i-1)+k(6,2) K(3*j,3*i)=K(3*j,3*i)+k(6,3) K(3*j,3*j-2)=K(3*j,3*j-2)+k(6,4) K(3*j,3*j-1)=K(3*j,3*j-1)+k(6,5) K(3*j,3*j)=K(3*j,3*j)+k(6,6) y=K; end