Guest User

Untitled

a guest
Sep 20th, 2019
167
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. % Curva Cp
  2. % Código: Dener A. L. Brandão - UFMG 2019
  3. clc; % Limpa linha de comandos
  4. clear; % Limpa variáveis do Workspace
  5.  
  6. R = 37.5; % Raio das pás: 37,5 m
  7. P_turb = 1.5e6; % Potência da turbina: 1,5 MW
  8. p = 1.225; % Densidade do ar (1,225 kg/m^3 a 15ºC e 1 atm
  9. A = pi*R^2;
  10. beta = 0:5:30; % Variação de beta de 0º até 30º, passo 5º
  11. fator_rad = pi/30;
  12. w_turb = 0:3:30; % Variação de w_turb de 0 a 30 rpm
  13. w_turb = fator_rad.*w_turb;
  14. k = 0;
  15. for h = 1:1:length(beta)
  16.     figure(h)
  17.  for V_ven = 5.5:1:14.5 % Variação da velocidade do vento;
  18.     for i = 1:1:length(w_turb)
  19.         lambda(i) = w_turb(i)*R/V_ven;
  20.         % A divisão ./ é utilizada para dividir cada ponto do vetor.
  21.         lambda_i(i) = 1/(1/(0.08*beta(h) + lambda(i)) - 0.035/(beta(h)^3 + 1)); % Fórmula de lambda_i
  22.  
  23.         % Cada vetor da matriz Cp corresponde a um valor de Beta. Ex: i = 1 @ beta
  24.         % = 2º
  25.         Cp(i) = 0.22*(116/lambda_i(i) - 0.4*beta(h) - 5)*exp(-12.5/lambda_i(i)); %Fórmula de Cp.
  26.     end
  27.   %  handle = plot(w_turb*30/pi,Cp); %Plota cada vetor de Cp
  28.     handle = plot(lambda,Cp); %Plota cada vetor de Cp
  29.     set(handle,'LineWidth',2);
  30.     hold on
  31.  end
  32. grid on
  33. axis([0 9 0 0.44]);
  34. ylabel('Cp','FontSize',15);
  35. xlabel('\lambda (Relação de Velocidades)','FontSize',15);
  36. title('Curva Cp','FontSize',15);
  37. set(gca,'FontSize',15)
  38. hold off
  39. end
RAW Paste Data