%**************************************************************************
%
%
%                 Volume Constraint 
%       
%
%**************************************************************************
%   Authors:Renato Pavanello
%	Copyright (c) Outubro 2023  by The LTM/FEM/UNICAMP.	
%   $ Revision: 1.0 $					
%**************************************************************************
function [c,ceq]=nlcon(Tgeo)

global inci coord nel nnos Tmat fg FreeDofs ug ite vComp vVol Wlim

We=zeros(nel,1);
for e=1:nel
     noi=inci(e,5);    
     noj=inci(e,6);
     xi=coord(noi,1);
     yi=coord(noi,2);
     xj=coord(noj,1);
     yj=coord(noj,2);
     l=sqrt((xj-xi)^2+ (yj-yi)^2);
     rho=Tmat(inci(e,3) ,4);
     A=Tgeo(inci(e,4) ,1);
     We(e)=A*l*rho;
end
 c=sum(We)-Wlim;
 ceq=0;