1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79
| clc,clear all,close all
%......parametre de l'eolienne ....%
%......Cp parameters....%
f=50; %.....frequence (HZ)....%
c1=0.5176;
c2=116;
c3=0.4;
c4=5;
c5=21;
c6=0.068;
Beta=[0;5;20]
Beta1=0; %....L'angle de calage (Degres)....%
Beta2=5;
Beta3=10;
Beta4=15;
Beta5=20;
lambda =[0:0.5:20]; %la vitesse reduite....%
vwind=[4:1:12]; %.....vitesse de vent (m/s)....%
rho=1.229;%.....la masse volumique de l'air (Kg/m^3).......%
R=42.5; %.....rayon du rotor (m).......%
omega=53; %nominal speed [rpm]
f=50; %frequency [HZ]
p=2; %pole pairs
%....calcul du vitesse reduite....%
%....lambda=omega_n*R/vwind....%;
lambda_i=1./(lambda+0.08.*Beta1)-0.035./(Beta1.^3+1);
lambda_i=1./(lambda+0.08.*Beta2)-0.035./(Beta2.^3+1);
lambda_i=1./(lambda+0.08.*Beta3)-0.035./(Beta3.^3+1);
lambda_i=1./(lambda+0.08.*Beta4)-0.035./(Beta4.^3+1);
lambda_i=1./(lambda+0.08.*Beta5)-0.035./(Beta5.^3+1);
%....calcul du coefficient de puissance....%
cp1 = c1*(c2*lambda_i - c3*Beta1 - c4 ).* exp( -c5 *lambda_i)+c6.*lambda;
cp2 = c1*(c2*lambda_i - c3*Beta2 - c4 ).* exp( -c5 *lambda_i)+c6.*lambda;
cp3 = c1*(c2*lambda_i - c3*Beta3 - c4 ).* exp( -c5 *lambda_i)+c6.*lambda;
cp4 = c1*(c2*lambda_i - c3*Beta4 - c4 ).* exp( -c5 *lambda_i)+c6.*lambda;
cp5 = c1*(c2*lambda_i - c3*Beta5 - c4 ).* exp( -c5 *lambda_i)+c6.*lambda;
cp1=1/2*(cp1+abs(cp1))%....uniquement pour les valeurs positif ....%
cp2=1/2*(cp2+abs(cp2))%....uniquement pour les valeurs positif ....%
cp3=1/2*(cp3+abs(cp3))%....uniquement pour les valeurs positif ....%
cp4=1/2*(cp4+abs(cp4))%....uniquement pour les valeurs positif ....%
cp5=1/2*(cp5+abs(cp5))%....uniquement pour les valeurs positif ....%
hold on
figure(1)
axis([0,20,0,1])
plot(lambda,cp1,'*r',lambda,cp2,'*m',lambda,cp3,'*g',lambda,cp4,'*k',lambda,cp5,'*b')
title('Cp-Lambda')
xlabel('Lambda')
ylabel('Cp')
legend('beta1= 0','beta2= 5','beta3= 10','beta4= 15','beta5= 20')
grid on
hold off
%....calcul de puissance aerodynamique....%
for k=1:length(vwind)
Paer = (1/2*rho*pi*R^2).*vwind(k).^3.*cp1;
%.....calcul du coefficient aerodynamique......%
cq1=cp1/omega;
cq2=cp2/omega;
cq3=cp3/omega;
cq4=cp4/omega;
cq5=cp5/omega;
%.....calcul du couple aerodynamique......%
Ta = 1/2*rho*pi*R^3.*vwind(k).^2.*cq1;
cq1=1/2*(cq1+abs(cq1))%....uniquement pour les valeurs positif ....%
cq2=1/2*(cq2+abs(cq2))%....uniquement pour les valeurs positif ....%
cq3=1/2*(cq3+abs(cq3))%....uniquement pour les valeurs positif ....%
cq4=1/2*(cq4+abs(cq4))%....uniquement pour les valeurs positif ....%
cq5=1/2*(cq5+abs(cq5))%....uniquement pour les valeurs positif ....%
end
hold on
figure(2)
axis([0,20,0,0.1])
plot(lambda,cq1,'*m',lambda,cq2,'*k',lambda,cq3,'*g',lambda,cq4,'*b',lambda,cq5,'*c')
title('Cq-Lambda ')
xlabel('la vitesse reduite')
ylabel('coefficient de couple')
legend('beta1= 0','beta2= 5','beta3= 10','beta4= 15','beta5= 20')
grid on
hold off
hold on |
Partager