请教：关于板单元的Global stiffness

xurundong

板单元得到的单元刚度矩阵是每个节点3个自由度的Wi,θxi,θyi。

而且这些自由度都是基于板单元单元坐标系----也就是一般重合于总体坐标系。

怎样才能得到可以集合到整体坐标系下的刚度矩阵呢？

或者说单元刚度矩阵是怎样集成到总刚中去的呢？

请指教，代码，公式都行。

topview

一般的有限元书籍上都有具体的做法，建议参考一下。

pengweicai

function y = Assemble(K,ke,Elem,Node,j)
for INODE=1:4               %每个单元节点指针  这里为4节点单元
    NODEI=Elem(j,INODE);               %第 j个单元的节点编号
    for IDOFN=1:3                   %  每个节点的自由度数
        NROWS=(NODEI-1)*3+IDOFN;    %  整体刚度矩阵中的行
        NROWE=(INODE-1)*3+IDOFN;    %  单刚中的行
        for JNODE=1:4
            NODEJ=Elem(j,JNODE);
            for JDOFN=1:3                       %与上面的行相同
                NCOLS=(NODEJ-1)*3+JDOFN;
                NCOLE=(JNODE-1)*3+JDOFN;
                K(NROWS,NCOLS)=K(NROWS,NCOLS)+ ke(NROWE,NCOLE);
            end
        end
    end
end

pengweicai

上面是matlab的代码

lwm11314

:kiss: ;P ;P ;P