billinezk.m

来自「有限元平面问题的求解相关程序。 包括质量矩阵」· M 代码 · 共 23 行

function ZK=billinezk(ar,line,nodnum,tnod,k1,k2)
tnum=zeros(8,1); %定义当前组装单元节点对应行号存储程序
ttnum=zeros(2,2); %定义当前组装单元节点号存储程序
ZK=zeros(nodnum*2,nodnum*2);    %组装后全部单元的总刚度矩阵
for j=1:ar
    for i=1:line
    ttnum=tnod(i:i+1,j:j+1);
    tnum(1,1)=2*ttnum(2,1)-1;
    tnum(2,1)=2*ttnum(2,1);
    tnum(3,1)=2*ttnum(2,2)-1;
    tnum(4,1)=2*ttnum(2,2);
    tnum(5,1)=2*ttnum(1,2)-1;
    tnum(6,1)=2*ttnum(1,2);
    tnum(7,1)=2*ttnum(1,1)-1;
    tnum(8,1)=2*ttnum(1,1);
    if (rem(i,3)==rem(j,3) || (rem(i,3)==0 && rem(j,3)==1 || (rem(i,3)==1 && rem(j,3)==0)))
        ZK=Asube(ZK,k1,tnum); %总刚度矩阵
    end
    if (rem(i,3)==2 && rem(j,3)==1) ||(rem(i,3)==0 && rem(j,3)==2) || (rem(i,3)==1 && rem(j,3)==2) || (rem(i,3)==2 && rem(j,3)==0)
        ZK=Asube(ZK,k2,tnum); %总刚度矩阵
    end
    end
end