%************************************************************************** % % Plot malha: Plot de Tensões em uma treliça % % %numbering == 1: Plota a malha em preto com os nós e elementos % 0: Plota a malha em preto somente % % % inci(tmat, tgeo, no1, no2 no3,...) % %coord( x1,y1; x2,y2 ; ...) % % Author: Renato Pavanello % Copyright (c) 2016 by The DMC/FEM/UNICAMP. % $Revision: 1.0 $ % %************************************************************************** function plotSxx(coord,inci,Sxx,tit) % Plota a malha de linhas com cores proporcial a Sxx nel=length(inci(:,1)); %numero de elementos nnos=length(coord(:,1));%numero de nós nnel=size(inci(:,5:end),2); %Numero de ntoseleme por elemento dimension=size(coord,2); %Dimensao da malha colormap(jet); %Initializacao das matrices X = zeros(nnel,1) ; Y = zeros(nnel,1) ; Z=zeros(nnel,1); if dimension==2 %For 2D plots for iel=1:nel for i=1:nnel X(i,1)=coord(inci(iel,i+4),1); Y(i,1)=coord(inci(iel,i+4),2); C(i,1)=Sxx(iel); end surf([X X],[Y Y],[Z Z],[C C],'EdgeColor','interp','FaceColor','none','Marker','o','LineWidth',2) ;hold on; end hold off; shading flat %// so each line segment has a plain color view(2) %// view(0,90) %// set view in X-Y plane colorbar(); title(tit) ; axis off; axis equal; end