bilinearquadelementis.m

四节点矩形单元matlab有限元程序

function y=BilinearQuadElementIS(ND,NF,NEE,ME)
for I=1:ND
    for J=1:NF
        II=NF*(I-1)+J;        
        IS(II)=NF*(ME(I,NEE)-1)+J;
    end
end
y=IS;