Bonjour,

voila je trouve les courbes ccomme je cherche le probleme quand je cherche trouver le courbe de Cp en fonction de lambda et beta en 3D avec plot3 ou bien meshgrid je ne trouve pas la solution

merci

Code : Sélectionner tout - Visualiser dans une fenêtre à part
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