Scribd
Importer
5 vues3 pages

Analyse du Profil NACA 4412

Le document présente une analyse du profil aérodynamique NACA 4412, en calculant les caractéristiques telles que la pente de la ligne de cambrure, le coefficient de portance et le moment. Il utilise des méthodes de discrétisation et des coefficients de Fourier pour déterminer les propriétés du profil. Les résultats incluent des graphiques de la cambrure et de sa pente, ainsi que des valeurs numériques pour les coefficients calculés.

Transféré par

Fatimata DIALLO
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
5 vues3 pages

Analyse du Profil NACA 4412

Le document présente une analyse du profil aérodynamique NACA 4412, en calculant les caractéristiques telles que la pente de la ligne de cambrure, le coefficient de portance et le moment. Il utilise des méthodes de discrétisation et des coefficients de Fourier pour déterminer les propriétés du profil. Les résultats incluent des graphiques de la cambrure et de sa pente, ainsi que des valeurs numériques pour les coefficients calculés.

Transféré par

Fatimata DIALLO
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

%% NACA 4-Chiffres - Théorie du Profil Mince (NACA 4412)

% Paramètres du profil

m = 0.04; % cambrure maximale (fraction de la corde)

p = 0.4; % position de la cambrure maximale (fraction de la corde)

t = 0.12; % épaisseur maximale (fraction de la corde)

alpha = 5 * pi/180; % angle d'incidence en radians

% Discrétisation

N = 1000;

theta = linspace(0, pi, N);

xi = (1 - cos(theta)) / 2;

% Calcul de dyc/dx

dyc_dx = zeros(size(xi));

for i = 1:N

if xi(i) <= p

dyc_dx(i) = 2*m/p^2 * (p - xi(i));

else

dyc_dx(i) = 2*m/(1-p)^2 * (p - xi(i));

end

end

% Calcul des coefficients de Fourier

A0 = (2/pi) * trapz(theta, dyc_dx .* cos(0*theta));

A1 = (2/pi) * trapz(theta, dyc_dx .* cos(1*theta));

A2 = (2/pi) * trapz(theta, dyc_dx .* cos(2*theta));


% Pente de portance théorique

dCl_dAlpha = 2*pi;

% Angle de portance nulle

alpha_L0 = A0 / 2;

% Coefficient de portance

Cl = dCl_dAlpha * (alpha - alpha_L0);

% Moments

Cm_c4 = -pi/4 * (A1 - A2);

Cm_LE = Cm_c4 + Cl/4;

% Position du centre de poussée

xCP_c = -Cm_LE / Cl;

% Affichage des résultats

fprintf('A0 = %.4f\n', A0);

fprintf('A1 = %.4f\n', A1);

fprintf('A2 = %.4f\n', A2);

fprintf('dCl/dAlpha = %.4f per rad\n', dCl_dAlpha);

fprintf('Alpha_L0 = %.4f deg\n', alpha_L0*180/pi);

fprintf('Cl = %.4f\n', Cl);

fprintf('Cm_c/4 = %.4f\n', Cm_c4);

fprintf('Cm_LE = %.4f\n', Cm_LE);

fprintf('xCP/c = %.4f\n', xCP_c);


% Tracé de la cambrure et de sa pente

figure;

subplot(2,1,1);

plot(xi, dyc_dx, 'b','LineWidth',1.5);

xlabel('\xi = x/c'); ylabel('dy_c/dx');

title('Pente de la ligne de cambrure'); grid on;

subplot(2,1,2);

% Calcul de yc(x)

yc = zeros(size(xi));

for i = 1:N

if xi(i) <= p

yc(i) = m/p^2 * (2*p*xi(i) - xi(i)^2);

else

yc(i) = m/(1-p)^2 * ((1-2*p) + 2*p*xi(i) - xi(i)^2);

end

end

plot(xi, yc, 'r','LineWidth',1.5);

xlabel('\xi = x/c'); ylabel('y_c/c');

title('Ligne de cambrure du NACA 4412'); grid on;

Vous aimerez peut-être aussi