Scribd
Importer
37 vues2 pages

Coefficient de puissance Cp en MATLAB

Le code MATLAB simule le coefficient de puissance Cp en fonction du ratio de vitesse lambda pour différentes orientations de pales (beta). Il utilise des coefficients définis pour calculer Cp et trace les résultats sur un graphique. Les valeurs complexes sont supprimées pour assurer la validité des résultats affichés.

Transféré par

Lynda Meziane
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
37 vues2 pages

Coefficient de puissance Cp en MATLAB

Le code MATLAB simule le coefficient de puissance Cp en fonction du ratio de vitesse lambda pour différentes orientations de pales (beta). Il utilise des coefficients définis pour calculer Cp et trace les résultats sur un graphique. Les valeurs complexes sont supprimées pour assurer la validité des résultats affichés.

Transféré par

Lynda Meziane
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

Code matlab :

clc; clear; close all;

% Définition des coefficients


C1 = 0.5176;
C2 = 116;
C3 = 0.4;
C4 = 5;
C5 = 21;
C6 = 0.0068;

% Plage de valeurs pour lambda (ratio de vitesse)


lambda = linspace(0, 16, 200);

% Angles d'orientation des pales


beta_values = [0, 5, 10, 15, 20];

% Initialisation de la figure
figure;
hold on;
grid on;

% Boucle pour tracer Cp en fonction de lambda pour différentes valeurs de beta


for beta = beta_values
lambda_1 = 1 ./ (1 ./ (lambda + 0.08 .* beta) - 0.035 ./ (beta.^3 + 1));
Cp = C1 * (C2 ./ lambda_1 - C3 .* beta - C4) .* exp(-C5 ./ lambda_1);

% Correction : suppression des valeurs complexes


Cp(imag(Cp) ~= 0) = NaN;

% Tracé des courbes


plot(lambda, Cp, 'DisplayName', sprintf('\\beta = %d°', beta), 'LineWidth', 1.5);
end

% Mise en forme du graphique


xlabel('Ratio de vitesse \lambda', 'FontSize', 12);
ylabel('Coefficient de puissance C_p', 'FontSize', 12);
title('Coefficient de puissance en fonction du ratio de vitesse et de l’angle d’orientation',
'FontSize', 13);
legend('show', 'Location', 'Best');
ylim([0 0.5]);
hold off;

Résultats de la simulation :

Vous aimerez peut-être aussi