Matlab-杆单元整体刚度矩阵组装（有限元基础-曾攀）

Smith Eric

已于 2022-05-05 10:53:11 修改

阅读量1w

点赞数 20

分类专栏：有限元程序文章标签： matlab 线性代数性能优化学习力学

于 2022-04-29 17:31:17 首次发布

本文链接：https://blog.csdn.net/qq_39642029/article/details/124491784

版权

有限元程序专栏收录该内容

8 篇文章

订阅专栏

一维杆单元的组装：

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;

二维杆单元组装

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;

三维情况下以此类推。