SHARE
TWEET

Untitled

a guest Sep 20th, 2019 108 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. % Cálculo de Potência: Dener A. L. Brandão - UFMG 2019
  2. R = 37.5; % Raio das pás: 37,5 m
  3. P_turb = 1.5e6; % Potência da turbina: 1,5 MW
  4. p = 1.225; % Densidade do ar (1,225 kg/m^3 a 15ºC e 1 atm
  5. A = pi*R^2;
  6. beta = 0; % Ângulo de ataque igual a 0
  7. i = 0;
  8. j = 0;
  9. for V_ven = 5.5:1:14.5 % Variação da velocidade do vento;
  10.     j = j+1;
  11.     fator_rad = pi/30;
  12.     w_turb = 0:3:30; % Variação de w_turb de 0.1 a 30 rpm
  13.     w_turb1 = fator_rad.*w_turb;
  14.     lambda = w_turb1.*R/V_ven;
  15.     lambda_i = 1./(1./(lambda + 0.08*beta) - 0.035/(beta^3 + 1)); % Fórmula de lambda_i
  16.     Cp = 0.22.*(116./lambda_i - 0.4*beta - 5).*exp(-12.5./lambda_i); %Fórmula de Cp.
  17.     P = p*A*V_ven^3.*Cp/2;
  18.     handle = plot(w_turb, P/1000); %Plota cada vetor de potência
  19.     ylabel('Potência (kW)','FontSize',15);
  20.     xlabel('V_{Turbina} (RPM)','FontSize',15);
  21.     title('Gráfico Potência x Velocidade da Turbina','FontSize',15);
  22.     set(gca,'FontSize',15); % Define tamanho da fonte dos eixos
  23.     set(handle,'LineWidth',[2]);
  24.     grid on;
  25.     hold on;
  26.     legend on
  27.     legend('5,5 m/s', '6,5 m/s','7,5 m/s','8,5 m/s','9,5 m/s','10,5 m/s','11,5 m/s','12,5 m/s','13,5 m/s','14,5 m/s','Max_P');
  28.     axis([0 30 0 4e3]);
  29. end
  30. for i = 1:1:length(P)
  31.         Max_P(i) = max(P(i));
  32.         if (Max_P(i) >= 1.5e6)
  33.             Max_P(i) = 1.5e6;
  34.         end
  35.     end
  36. handle2 = plot(w_turb, Max_P/1000,'m');
  37. set(handle2,'LineWidth',[2]);
  38. hold on
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