MATLAB中基于结构动力学原理生成梁单元刚度矩阵，并且可以对单元刚度矩阵进行组装（支持一点程度的材料变化）

%矩阵组装中，不同结构的矩阵连接（两者的连接）% function K=AssemblytoolConnection(Ke1,n1,Ke2,n2) %Ke1,Ke2为待组装的两个半成品矩阵，由part1,part2生成% %Ne1,Ne2为两个组成两个半成品矩阵对应的单位矩阵阶数% %n1,n2为识别参数，取值范围：Ne=3,n=123; Ne=2,n=12,n=13,n=23; Ne=1,n=1,n=2,n=3% %单元矩阵规则为：轴向，水平，转动% if n1==123&&n2==123 Ne1=size(Ke1); Ne1=Ne1(1); Ne2=size(Ke2); Ne2=Ne2(1); K=zeros(Ne1+Ne2-3,Ne1+Ne2-3); for i=1:Ne1 for j=1:Ne1 a=Ke1(i,j); K(i,j)=K(i,j)+a; end end for i=1:Ne2 for j=1:Ne2 b=Ke2(i,j); K(Ne1-3+i,Ne1-3+j)=K(Ne1-3+i,Ne1-3+j)+b; end end else K=0; end