Matlab-有限单元法-2D梁单元的刚度矩阵组装（曾攀）

Smith Eric

已于 2022-05-08 19:00:33 修改

阅读量2.3k

点赞数 2

分类专栏：有限元程序文章标签：力学

于 2022-04-30 22:52:48 首次发布

本文链接：https://blog.youkuaiyun.com/qq_39642029/article/details/124520532

版权

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

8 篇文章

订阅专栏

function z=Bar2D2Node_Assembly(KK,k,i,j)
%input:单元刚度矩阵k,单元的节点编号i，j
%output:总刚度矩阵
 
%A点，3个自由度
DOF(1)=3*i-2;
DOF(2)=3*i-1;
DOF(3)=3*i;
%B点，3个自由度
DOF(4)=3*j-2;
DOF(5)=3*j-1;
DOF(6)=3*j;
%一个梁单元具有6个自由度
for n1=1:6
    for n2=1:6
        KK(DOF(n1),DOF(n2))=KK(DOF(n1),DOF(n2))+k(n1,n2);
    end
end
z=KK;