SHARE
TWEET

Untitled

a guest Sep 20th, 2019 100 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
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top