%--------------------------------------------------------------------------------
% Carlos Manuel Ribeiro Cabral  Nº16141

% ROBUCAR - Projecto de Automação/Seminário de processos Tecnológicos (2000/2001)
%--------------------------------------------------------------------------------



% Calculo da reacção ou da força que é necessário exercer na barra transversal para elevar a rampa - Calculado de uma forma analitica
% (ver graficos nº1 do relatório final)
%---------------------------------------------------
clear all
%---------------------------------------------------

syms F beta alpha

P=2649;
c=105;
a=602;


F = (P*c*cos(beta*pi/180)*sin(alpha*pi/180+beta*pi/180)*(sin(beta*pi/180)+cos(beta*pi/180)*tan(alpha*pi/180)))...
   /(a*tan(alpha*pi/180));

domain = [3, 7, eps, 25];
ezsurf(F, domain,65);
title('grafico da força (N)')

T = (P*c*(cos(beta*pi/180))^2*sin(beta*pi/180+alpha*pi/180))/1000;

figure(2)
ezsurf(T, domain,65);
title(' gráfico do binário (N.m)')


F = (P*c*cos(beta*pi/180)*sin(alpha*pi/180+beta*pi/180)*(sin(beta*pi/180)+cos(beta*pi/180)*tan(alpha*pi/180)))...
   /(a*tan(alpha*pi/180));
figure(3)
domain = [6, 15, eps, 25];
ezsurf(F, domain,65);
title('grafico da força (N)')

T = (P*c*(cos(beta*pi/180))^2*sin(beta*pi/180+alpha*pi/180))/1000;

figure(4)
ezsurf(T, domain,65);
title(' gráfico do binário (N.m)')

F = (P*c*cos(beta*pi/180)*sin(alpha*pi/180+beta*pi/180)*(sin(beta*pi/180)+cos(beta*pi/180)*tan(alpha*pi/180)))...
   /(a*tan(alpha*pi/180));
figure(5)
domain = [15, 70, eps, 25];
ezsurf(F, domain,65);
title('grafico da força (N)')

T = (P*c*(cos(beta*pi/180))^2*sin(beta*pi/180+alpha*pi/180))/1000;

figure(6)
ezsurf(T, domain,65);
title(' gráfico do binário (N.m)')