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

global inci coord nel nnos Tmat fg FreeDofs ug ite vComp vVol Wlim EEC ECClim
We=zeros(nel,1);
ECCe=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=(X(e+nel)^3)*(Tmat(1,4)-Tmat(2,4))+Tmat(2,4);
     A=X(inci(e,4),1);
     We(e)=A*l*rho;
     vVol(e)=A*l;
     rhoECC=(X(e+nel)^3)*((Tmat(1,5)*Tmat(1,4))-(Tmat(2,5)*Tmat(2,4)))+(Tmat(2,4)*Tmat(2,5));
     ECCe(e)=A*l*rhoECC;

end
c=sum(We)-Wlim;
%c=sum(ECCe)-ECClim;
ceq=0;