cummulatedrestrainedlist3.m

平面刚架（带弹簧支座）有限元程序

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;%弹性约束位移数