Code covered by the BSD License

# The Plane Stress Problem

### Siva Srinivas Kolukula (view profile)

13 Jun 2011 (Updated )

Plate under uniform tension at its edges is solved using Finite Element Method

PlotFieldonMesh(coordinates,nodes,component)
```function PlotFieldonMesh(coordinates,nodes,component)
%--------------------------------------------------------------------------
% Code written by : Siva Srinivas Kolukula                                |
%                   Senior Research Fellow                                |
%                   Structural Mechanics Laboratory                       |
%                   Indira Gandhi Center for Atomic Research              |
%                   India                                                 |
% E-mail : allwayzitzme@gmail.com                                         |
%--------------------------------------------------------------------------
%--------------------------------------------------------------------------
% Purpose:
%         To plot the profile of a component on mesh
% Synopsis :
%           ProfileonMesh(coordinates,nodes,component)
% Variable Description:
%           coordinates - The nodal coordinates of the mesh
%           -----> coordinates = [node X Y ]
%           nodes - The nodal connectivity of the elements
%           -----> nodes = [node1 node2......]
%           component - The components whose profile to be plotted
%           -----> components = a column vector in the order of node
%                               numbers
%--------------------------------------------------------------------------

nel = length(nodes) ;                  % number of elements
nnode = length(coordinates) ;          % total number of nodes in system
nnel = size(nodes,2);                % number of nodes per element
%
% Initialization of the required matrices
X = zeros(nnel,nel) ;
Y = zeros(nnel,nel) ;
Z = zeros(nnel,nel) ;
profile = zeros(nnel,nel) ;
%
for iel=1:nel
for i=1:nnel
nd(i)=nodes(iel,i);         % extract connected node for (iel)-th element
X(i,iel)=coordinates(nd(i),1);    % extract x value of the node
Y(i,iel)=coordinates(nd(i),2);    % extract y value of the node
end
profile(:,iel) = component(nd') ;         % extract component value of the node
end

% Plotting the FEM mesh and profile of the given component
f3 = figure ;
set(f3,'name','Postprocessing','numbertitle','off') ;
plot(X,Y,'k')
fill(X,Y,profile)
axis off ;
% Colorbar Setting
SetColorbar
end

```