MEFtave/Dessin/plot2D.m

function plot2D
% trace du maillage d'une structure plane 
% avec les numero des noeuds et des elements
% les conditions aux limites (deplacements imposes) 
% et les Charges nodales
%
% appel plot2D
%
%  H.Oudin 
global nddln nnod nddlt nelt nnode ndim
global Coord Connec Ncl F

hold on
title('maillage de la structure')
axis equal

XY = Coord;
Cu = Ncl(1:nddln:nddlt); Cv = Ncl(2:nddln:nddlt);
Fx = F(1:nddln:nddlt);   Fy = F(2:nddln:nddlt);
if ndim == 1 
    XY(:,2) = zeros(nnod,1); 
    if (nddln == 1) Cv = zeros(nnod,1); Fy = zeros(nnod,1); end
end
  
for inod=1:nnod     % ----- Noeuds , conditions aux limites et charges nodales
    text(XY(inod,1),XY(inod,2),[' ',int2str(inod)],'color','m','FontSize',14)
    if Cu(inod)==1 plot(XY(inod,1),XY(inod,2),'>','color','k','MarkerSize',10);end
    if Cv(inod)==1 plot(XY(inod,1),XY(inod,2),'^','color','k','MarkerSize',10);end
    if Fx(inod)~= 0  
        text(XY(inod,1),XY(inod,2),['  \rightarrow',num2str(Fx(inod))],...
            'HorizontalAlignment','left','color','g','FontSize',14);
    end
    if Fy(inod)~= 0  
        text(XY(inod,1),XY(inod,2),['\uparrow',num2str(Fy(inod))],...
            'VerticalAlignment','bottom','color','g','FontSize',14);
    end   
end
plot(XY(:,1),XY(:,2),'.','color','m','MarkerSize',12);

for iel = 1:nelt        %----- visualisation du maillage
    loce=[];            %----- table de localisation pour l'element
    for i=1:nnode 
		if Connec(iel,i) > 0 loce=[loce,Connec(iel,i)]; end
	end;  
    Pos  = XY(loce,:);
    X = [[Pos(:,1)]; Pos(1,1)];Y = [[Pos(:,2)]; Pos(1,2)];
    line(X,Y,'color','b')
    text(mean(Pos(:,1)),mean(Pos(:,2)),[' ',int2str(iel)],'color','b','FontSize',14)
end 
return