% Curva Cp
% Código: Dener A. L. Brandão - UFMG 2019
clc; % Limpa linha de comandos
clear; % Limpa variáveis do Workspace

R = 37.5; % Raio das pás: 37,5 m
P_turb = 1.5e6; % Potência da turbina: 1,5 MW
p = 1.225; % Densidade do ar (1,225 kg/m^3 a 15ºC e 1 atm
A = pi*R^2;
beta = 0:5:30; % Variação de beta de 0º até 30º, passo 5º
fator_rad = pi/30;
w_turb = 0:3:30; % Variação de w_turb de 0 a 30 rpm
w_turb = fator_rad.*w_turb;
k = 0;
for h = 1:1:length(beta)
    figure(h)
 for V_ven = 5.5:1:14.5 % Variação da velocidade do vento;
    for i = 1:1:length(w_turb)
        lambda(i) = w_turb(i)*R/V_ven;
        % A divisão ./ é utilizada para dividir cada ponto do vetor.
        lambda_i(i) = 1/(1/(0.08*beta(h) + lambda(i)) - 0.035/(beta(h)^3 + 1)); % Fórmula de lambda_i

        % Cada vetor da matriz Cp corresponde a um valor de Beta. Ex: i = 1 @ beta
        % = 2º
        Cp(i) = 0.22*(116/lambda_i(i) - 0.4*beta(h) - 5)*exp(-12.5/lambda_i(i)); %Fórmula de Cp.
    end
  %  handle = plot(w_turb*30/pi,Cp); %Plota cada vetor de Cp
    handle = plot(lambda,Cp); %Plota cada vetor de Cp
    set(handle,'LineWidth',2);
    hold on
 end
grid on
axis([0 9 0 0.44]);
ylabel('Cp','FontSize',15);
xlabel('\lambda (Relação de Velocidades)','FontSize',15);
title('Curva Cp','FontSize',15);
set(gca,'FontSize',15)
hold off
end