cummulatedrestrainedlist3.m

来自「平面刚架（带弹簧支座）有限元程序」· M 代码 · 共 23 行

function [n,ne,nr,crl]=CummulatedRestrainedList3(nnj,rl)
kk1=0;
for jx=1:nnj
    if rl(jx)==0
        kk1=kk1+1;
        crl(jx)=kk1;
    end
end
n=kk1;%非约束位移数
for jx=1:nnj
    if rl(jx)==1
        kk1=kk1+1;
        crl(jx)=kk1;
    end
end
ne=kk1-n;%弹性约束位移数
for jx=1:nnj
    if rl(jx)==2
        kk1=kk1+1;
        crl(jx)=kk1;
    end
end
nr=kk1-n-ne;%弹性约束位移数