一维杆单元的组装：
function z=Bar1D2Node_Assembly(KK,k,i,j)
%input:单元刚度矩阵k,单元的节点编号i，j
%output:总刚度矩阵
DOF(1)=i;
DOF(2)=j;
for n1=1:2
for n2=1:2
KK(DOF(n1),DOF(n2))=KK(DOF(n1),DOF(n2))+k(n1,n2);
end
end
z=KK;
二维杆单元组装文章来源：https://www.toymoban.com/news/detail-561854.html
function z=Bar2D2Node_Assembly(KK,k,i,j)
%input:单元刚度矩阵k,单元的节点编号i，j
%output:总刚度矩阵
%A点，两个自由度
DOF(1)=2*i-1;
DOF(2)=2*i;
%B点，两个自由度
DOF(3)=2*j-1;
DOF(4)=2*j;
for n1=1:4
for n2=1:4
KK(DOF(n1),DOF(n2))=KK(DOF(n1),DOF(n2))+k(n1,n2);
end
end
z=KK;
三维情况下以此类推。文章来源地址https://www.toymoban.com/news/detail-561854.html
到了这里，关于Matlab-杆单元整体刚度矩阵组装（有限元基础-曾攀）的文章就介绍完了。如果您还想了解更多内容，请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章，希望大家以后多多支持TOY模板网！
