plot_fem_data.m

FEM tools for caculation of nonlinear problems

 maxZ = max(in_data.ND(:,4));  minZ = min(in_data.ND(:,4));
 labx = (maxX / 9); laby = (maxY / 9); labz = (maxZ / 9);
 labx = min([labx laby labz]); laby = labx;  labz=labx;

for i=1:size(in_data.LOAD_)
   node_i = find(in_data.ND(:,1)==in_data.LOAD_(i,1));
   if in_data.LOAD_(i,2)~=0 & in_data.LOAD_(i,2)>0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)+labx],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'r-','LineWidth',3); hold on;
   end;
   if in_data.LOAD_(i,2)~=0 & in_data.LOAD_(i,2)<0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)-labx],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'r-','LineWidth',3); hold on;
   end;
   if in_data.LOAD_(i,3)~=0 & in_data.LOAD_(i,3)>0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)+laby],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'r-','LineWidth',3); hold on;
   end;
   if in_data.LOAD_(i,3)~=0 & in_data.LOAD_(i,3)<0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)-laby],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'r-','LineWidth',3); hold on;
   end;
   if in_data.LOAD_(i,4)~=0 & in_data.LOAD_(i,4)>0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)+labz],'r-','LineWidth',3); hold on;
   end;
   if in_data.LOAD_(i,4)~=0 & in_data.LOAD_(i,4)<0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)-labz],'r-','LineWidth',3); hold on;
   end;
end;

for i=1:size(in_data.CON)
   node_i = find(in_data.ND(:,1)==in_data.CON(i,1));
   if in_data.CON(i,2)==0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)-labx],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'g-','LineWidth',2); hold on;
   end;
   if in_data.CON(i,3)==0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)-laby],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'g-','LineWidth',2); hold on;
   end;
   if in_data.CON(i,4)==0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)-labz],'g-','LineWidth',2); hold on;
   end;
   if in_data.CON(i,4)==0 | in_data.CON(i,5)==0 | in_data.CON(i,6)==0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'gs','LineWidth',2); hold on;
   end;
end;

rotate3d(gca); set(gcf,'Pointer','arrow');
axis([(minX-labx) (maxX+labx) (minY-laby) (maxY+laby) (minZ-labz) (maxZ+labz)]);
hold off




function plot_input_tetrah ( in_data, dof_,PND) 

hold off;

axis equal;
if 0
for i=1:size(in_data.EL)
    node1 = find(in_data.ND(:,1)==in_data.EL(i,3));
    node2 = find(in_data.ND(:,1)==in_data.EL(i,4));
    node3 = find(in_data.ND(:,1)==in_data.EL(i,5));
    node4 = find(in_data.ND(:,1)==in_data.EL(i,6));

    plot3([in_data.ND(node1,2) in_data.ND(node2,2) in_data.ND(node3,2) in_data.ND(node1,2)  in_data.ND(node4,2)],...
       [in_data.ND(node1,3) in_data.ND(node2,3) in_data.ND(node3,3)    in_data.ND(node1,3)  in_data.ND(node4,3)],...
       [in_data.ND(node1,4) in_data.ND(node2,4) in_data.ND(node3,4)    in_data.ND(node1,4)  in_data.ND(node4,4)],...
          'Color',[0.4 0.1 0.7],'LineWidth',1);    

    plot3([in_data.ND(node2,2) in_data.ND(node4,2) in_data.ND(node3,2) ],...
       [in_data.ND(node2,3) in_data.ND(node4,3) in_data.ND(node3,3)],...
       [in_data.ND(node2,4) in_data.ND(node4,4) in_data.ND(node3,4)],...
         'Color',[0.4 0.1 0.7],'LineWidth',1);
     if PND == 1
         h=text(in_data.ND(node1,2) , in_data.ND(node1,3), in_data.ND(node1,4), num2str(node1)); set(h,'FontSize',8); 
         set(h,'Color','k');
         h=text(in_data.ND(node2,2) , in_data.ND(node2,3), in_data.ND(node2,4), num2str(node2)); set(h,'FontSize',8); 
         set(h,'Color','k');
         h=text(in_data.ND(node3,2) , in_data.ND(node3,3), in_data.ND(node3,4), num2str(node3)); set(h,'FontSize',8); 
         set(h,'Color','k');
         h=text(in_data.ND(node4,2) , in_data.ND(node4,3), in_data.ND(node4,4), num2str(node4)); set(h,'FontSize',8); 
         set(h,'Color','k');

             x1= in_data.ND(node1,2)/4  +  in_data.ND(node2,2)/4 +  in_data.ND(node3,2)/4  + in_data.ND(node4,2)/4;
             y1= in_data.ND(node1,3)/4  +  in_data.ND(node2,3)/4 +  in_data.ND(node3,3)/4  + in_data.ND(node4,3)/4;
             z1= in_data.ND(node1,4)/4  +  in_data.ND(node2,4)/4 +  in_data.ND(node3,4)/4  + in_data.ND(node4,4)/4;
             h=text(x1 , y1, z1, num2str(in_data.EL(i,1))); set(h,'FontSize',6); set(h,'Color','m');
     end;
 end;
end % ZERO ---
tri2=surftri(in_data.ND(:,2:4),in_data.EL(:,3:6));
trimesh(tri2,in_data.ND(:,2),in_data.ND(:,3),in_data.ND(:,4),'facecolor','none','edgecolor','k'); axis equal
axis off; view(3);
h=title('Input model'); set(h,'FontSize',12);  hold on;

 maxX = max(in_data.ND(:,2));  minX = min(in_data.ND(:,2));
 maxY = max(in_data.ND(:,3));  minY = min(in_data.ND(:,3));
 maxZ = max(in_data.ND(:,4));  minZ = min(in_data.ND(:,4));
 labx = (maxX / 9); laby = (maxY / 9); labz = (maxZ / 9);
 labx = min([labx laby labz]); laby = labx;  labz=labx;

for i=1:size(in_data.LOAD_,1)
   node_i = find(in_data.ND(:,1)==in_data.LOAD_(i,1));
   if in_data.LOAD_(i,2)~=0 & in_data.LOAD_(i,2)>0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)+labx],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'r-','LineWidth',3); hold on;
   end;
   if in_data.LOAD_(i,2)~=0 & in_data.LOAD_(i,2)<0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)-labx],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'r-','LineWidth',3); hold on;
   end;
   if in_data.LOAD_(i,3)~=0 & in_data.LOAD_(i,3)>0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)+laby],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'r-','LineWidth',3); hold on;
   end;
   if in_data.LOAD_(i,3)~=0 & in_data.LOAD_(i,3)<0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)-laby],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'r-','LineWidth',3); hold on;
   end;
   if in_data.LOAD_(i,4)~=0 & in_data.LOAD_(i,4)>0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)+labz],'r-','LineWidth',3); hold on;
   end;
   if in_data.LOAD_(i,4)~=0 & in_data.LOAD_(i,4)<0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)-labz],'r-','LineWidth',3); hold on;
   end;
end;

for i=1:size(in_data.CON,1)
   node_i = find(in_data.ND(:,1)==in_data.CON(i,1));
   if in_data.CON(i,2)==0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)-labx],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'g-','LineWidth',2); hold on;
   end;
   if in_data.CON(i,3)==0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)-laby],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'g-','LineWidth',2); hold on;
   end;
   if in_data.CON(i,4)==0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)-labz],'g-','LineWidth',2); hold on;
   end;
   if in_data.CON(i,4)==0 | in_data.CON(i,5)==0 | in_data.CON(i,6)==0
      plot3([in_data.ND(node_i,2) in_data.ND(node_i,2)],[in_data.ND(node_i,3) in_data.ND(node_i,3)],...
         [in_data.ND(node_i,4) in_data.ND(node_i,4)],'gs','LineWidth',2); hold on;
   end;
end;

axis([(minX-labx) (maxX+labx) (minY-laby) (maxY+laby) (minZ-labz) (maxZ+labz)]);
rotate3d(gca); set(gcf,'Pointer','arrow');

hold off;




function plot_heat_CST_input ( in_data, dof, NF, PND) 

hold off;
if 1
   maxX = max(in_data.ND(:,2)); minX = min(in_data.ND(:,2));
   maxY = max(in_data.ND(:,3)); minY = min(in_data.ND(:,3));
   labx = abs(maxX / 15); laby = abs(maxY / 15);
   labx = max([labx laby]); laby = labx;


   if 0

       % hold on
       for i=1:size(in_data.LOAD_)
           node_i = find(in_data.ND(:,1)==in_data.LOAD_(i,1));
           if in_data.LOAD_(i,2)~=0 & in_data.LOAD_(i,2)>0
               plot([in_data.ND(node_i,2) in_data.ND(node_i,2)+labx],...
                   [in_data.ND(node_i,3) in_data.ND(node_i,3)],'r-','LineWidth',3); hold on;
           end;
           if in_data.LOAD_(i,2)~=0 & in_data.LOAD_(i,2)<0
               plot([in_data.ND(node_i,2) in_data.ND(node_i,2)-labx],...
                   [in_data.ND(node_i,3) in_data.ND(node_i,3)],'r-','LineWidth',3);hold on;
           end;
           if in_data.LOAD_(i,3)~=0 & in_data.LOAD_(i,3)>0
               plot([in_data.ND(node_i,2) in_data.ND(node_i,2)],...
                   [in_data.ND(node_i,3) in_data.ND(node_i,3)+laby],'r-','LineWidth',3);hold on;
           end;
           if in_data.LOAD_(i,3)~=0 & in_data.LOAD_(i,3)<0
               plot([in_data.ND(node_i,2) in_data.ND(node_i,2)],...
                   [in_data.ND(node_i,3) in_data.ND(node_i,3)-laby],'r-','LineWidth',3);hold on;
           end;
       end;
  

       for i=1:size(in_data.CON)
           node_i = find(in_data.ND(:,1)==in_data.CON(i,1));
           if in_data.CON(i,2)==0
               plot([in_data.ND(node_i,2) in_data.ND(node_i,2)-labx],...
                   [in_data.ND(node_i,3) in_data.ND(node_i,3)],'g-','LineWidth',2.5);
           end;
           if in_data.CON(i,3)==0
               plot([in_data.ND(node_i,2) in_data.ND(node_i,2)],...
                   [in_data.ND(node_i,3) in_data.ND(node_i,3)-laby],'g-','LineWidth',2.5);
           end;
       end;
   end; % % -- end if NF > 100
   
   
   plot(in_data.ND(:,2),in_data.ND(:,3),'y.','MarkerSize',2);
   axis equal; axis off; view(2); hold on;
   for i=1:size(in_data.EL)
      yy = 0.9; % color scale: white to black 
      if in_data.EL(i,2)==44   
         node1 = find(in_data.ND(:,1)==in_data.EL(i,3));
         node2 = find(in_data.ND(:,1)==in_data.EL(i,4));
         node3 = find(in_data.ND(:,1)==in_data.EL(i,5));
         plot([in_data.ND(node1,2) in_data.ND(node2,2) in_data.ND(node3,2) in_data.ND(node1,2)], ...
            [in_data.ND(node1,3) in_data.ND(node2,3) in_data.ND(node3,3) in_data.ND(node1,3)],...
            'Color',[0.4 0.1 0.7],'LineWidth',.1);
        if PND == 1
            h=text(in_data.ND(node1,2) , in_data.ND(node1,3), 0, num2str(node1)); set(h,'FontSize',8); 
            set(h,'Color','k');
            h=text(in_data.ND(node2,2) , in_data.ND(node2,3), 0, num2str(node2)); set(h,'FontSize',8); 
            set(h,'Color','k');
            h=text(in_data.ND(node3,2) , in_data.ND(node3,3), 0, num2str(node3)); set(h,'FontSize',8); 
            set(h,'Color','k');
            x1= in_data.ND(node1,2)/3  +  in_data.ND(node2,2)/3 +  in_data.ND(node3,2)/3;
            y1= in_data.ND(node1,3)/3  +  in_data.ND(node2,3)/3 +  in_data.ND(node3,3)/3;
            h=text(x1 , y1, num2str(in_data.EL(i,1))); set(h,'FontSize',6); set(h,'Color','m');
        end;
      end;
 end;

 h=title('Input model'); set(h,'FontSize',12);
 axis([(minX-labx) (maxX+labx) (minY-laby) (maxY+laby)]); drawnow
 rotate3d(gca); set(gcf,'Pointer','arrow');
 hold off;
 
 % =======END
end