bilinearquadelementshowmesh.m

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

function BilinearQuadElementShowmesh(X,Y,ME,NE,ND)
%    LinearTriangleElement的作用是画出有限元网格、显示各个结点和单元的编号。
%    在结点上标出的是结点号，为了防止单元和结点的编号混淆，我把单元的编号放在单元中心位置，且用方括号表示，如[2]，表示第二个单元，
plot(0,0);
hold on;
for I=1:NE
       for J=1:ND
        MEE=ME(J,I);
        x(J)=X(MEE);
        y(J)=Y(MEE);
        text(x(J),y(J),int2str(MEE));
    end
    x(ND+1)=x(1);
    y(ND+1)=y(1);
    plot(x,y,'r');
    s=0;
    t=0;
    for m=1:ND
       s=s+x(m);
       t=t+y(m);
   end
   dx=x(1)-x(2);
   dy=y(1)-y(3);
   xx=s/ND;
   yy=t/ND;
   text(xx+dx/6,yy-dy/8,['[',int2str(I),']']);
end