bilinearquadelementme.m

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

function y=BilinearQuadElementME(ND,NEI,NEJ)
for J=1:NEJ
    for I=1:NEI
        NEE=(J-1)*NEI+I;
        ME(1,NEE)=I+(J-1)*(NEI+1);
        ME(2,NEE)=I+(J-1)*(NEI+1)+1;
        ME(3,NEE)=I+J*(NEI+1)+1;
        ME(4,NEE)=I+J*(NEI+1);
    end
end
y=ME;