0% ont trouvé ce document utile (0 vote)
223 vues54 pages

Dimensionnement D Un Rotor D Eolienne

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1/ 54

Angela Valle Guimaré

Tuteur: André Ferrand

DIMENSIONNEMENT
D’UN ROTOR
D’ÉOLIENNE
Théorie mathématique du dimensionnement d’un rotor
d’éolienne à axe horizontal.
Création d’un logiciel de dimensionnement.
Programmation avec VBA.
TABLE DES MATIÈRES
 Objectif

 CHAPITRE I: L’éolienne et l’énergie éolienne

 CHAPITRE II: Introduction à l’aérodynamique

 CHAPITRE III: Introduction au logiciel

 CHAPITRE IV: Base théorique du logiciel

 CHAPITRE V: Utilisation du logiciel

 Vérification des résultats

 Conclusions et possibles améliorations

2
Corrections
Corrections

La
La structure
structure Finalisation
Finalisation
du
du rapport
rapport

OBJECTIF
OBJECTIF

Vérification
Vérification
des Complémentation
des résultats
résultats

L’implantation
L’implantation Recherche
Recherche
dans
dans le
le logiciel
logiciel

3
CHAPITRE I
L’éolienne et l’énergie éolienne

4
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel
CHAPITRE V: Utilisation du logiciel

ENERGIE ET PUISSANCE
D’UN FLUX LIBRE

1 · 2 1 Vol 2 1 · 2
1) E = mV = ρ V = ρ QV
2 2 t 2

m 2 m 
3
Vol
2) = V ·A =  m  =  
t s   s 

1 1
P = ρV (V · A) = ρπR 2V 3
2

2 2 5
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel LE TUB DE
CHAPITRE V: Utilisation du logiciel
COURANT

• Ec  Emec
• V2<V1
.
Q = VA = cnt

A2>A1
Transformation d’énergie plus efficace
• V2≠0
• V2≠V1 ratio de vitesses

6
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel LA LOI
CHAPITRE V: Utilisation du logiciel
DE BETZ

• Théorie du Disque Actuateur


.
m = ρS
(v1 + v2 ) S = surface _ balayée _ rotor
2

• Bernouilli et Newton

( )
.
m v1 − v 2
2 2
P=
2
P=
ρ
4
(v
1
2
− v2
2
)(v + v )S
1 2

• Puissance d’un flux d’air non perturbé


ρ
P0 =
3
v1 S v2 1
2 =
v1 3

• Ratio: Puissance Récupérée vs.


P 16
Puissance du Vent = 0.59 =
P0 max 27
P 1   v2    v2  
2

= 1−   1+    16  1 8
P0 2   v1    v1   PBETZ =   ⋅ ⋅ ρ ⋅ π ⋅ R 2 ⋅V 3 = ⋅ ρ ⋅ π ⋅ R 2 ⋅V 3
   27  2 27

7
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel PUISSANCE
CHAPITRE V: Utilisation du logiciel
RECUPEREE

 16  1 8
PBETZ =   ⋅ ⋅ ρ ⋅ π ⋅ R 2 ⋅V 3 = ⋅ ρ ⋅ π ⋅ R 2 ⋅V 3
 27  2 27

ηhelice = 0,7 − 0,85


η mécanique = 0 ,98
η global
η alternater = 0,96
η transmissi on = 0,92
8
Préc = η global ⋅ ⋅ ρ ⋅ π ⋅ R 2 ⋅V 3
27
8
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel
CHAPITRE V: Utilisation du logiciel

SYSTÈMES DE
RÉGULATION
• Pitch  Pas variable OU FREINAGE
– Variation de l’angle de calage
de les pales
– Freinage: système
électronique
– Profit maximal du vent

• Stall  Pas fixe


– Pales fixes par rapport au
moyeu
– Freinage: décrochage
aérodynamique
– Pertes d’efficacité

9
CHAPITRE II
Introduction à l’aérodynamique
10
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel PORTANCE
CHAPITRE V: Utilisation du logiciel
TRAÎNEE

• Aérodynamique de l’aviation

• Aile d’avion ≈ Pale d’éolienne

• Portance

• Traînée  résistance à l’avancée

11
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel
CHAPITRE V: Utilisation du logiciel
RESULTANT
AERODYNAMIQUE
Rx
Cd = C x =
ρ Z
⋅ l ⋅V 2
2
Rz
Cl = C z =
ρ
⋅ l ⋅V 2
2

Cz = f(α) et Cx = g(α)

Vent relatif
Ra =2 RX2 + RZ2
X

φ : angle du vent relatif avec le plan de rotation [rad].


β: angle de calage du profil [rad].
α: angle d’incidence [rad].
12
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel LA POLAIRE
CHAPITRE V: Utilisation du logiciel
DE L’AILE

• 4 points remarquables

C
• Finesse f = z
C x

• f exprime la capacité M1:portance nulle

pour pénétrer à l’air M2:traînée minimale


M3:finesse maximale
M4:portance maximale

13
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel
CHAPITRE V: Utilisation du logiciel LE NOMBRE
DE REYNOLDS

• Alors, il existe une polaire d’aile unique pour chaque aile? NON!

polaire _ d ' aile = f (R e)


V est la vitesse relative (ms-1) de la pale par
rapport à l’air.

lVρ lV l est la longueur de la corde moyenne [m]

Re = = ρ est la masse volumique de l’air [1,225 kg/m3 à

µ υ 20 °C]
µ est la valeur de la viscosité [1,81 10-5 m2s-1 à
20 °C].

• Re = rapport entre les forces liées à la vitesse et les forces visqueuses

m2  Re ≈ 68 000 Vl
• Conditions normales: υ =1.45⋅10  
−5

s
14
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel
CHAPITRE V: Utilisation du logiciel
IMPORTANCE
DU NOMBRE DE REYNOLDS

• Re faible  les forces de


résistance dominent

• Re croît  portance et finesse


maximales plus grandes

• Une polaire pour chaque Re

• Préférable  Re élevé
Re=1E6 Re=0.5E6

1,8 1,6

1,6 1,4
1,4
1,2
1,2
1
1
Cl

Serie1
Cl

0,8 Serie1
0,8
0,6
0,6

0,4 0,4

0,2 0,2

0 0
0 5 10 15 20
Cd*1000
25 30 35 40 45 0 10 20 30 40 50 60
15
Cd * 1000
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel
LA COUCHE
CHAPITRE V: Utilisation du logiciel LIMITE

• Air ≈ lames infiniment minces


superposées
• Vitesse de l’écoulement libre
(Vel)
• Vitesses croissantes
• Couche limite  des lames se
déplaçant à des vitesses ≠ Vel
• Loi de variation  dépend de
Re
• ↓ Re  égalise les vitesses
• ↑ Re  lames + indépendants
16
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel
CHAPITRE V: Utilisation du logiciel

DECOLLEMENT
DE LA COUCHE LIMITE

• Couche limite laminaire

• Couche limite turbulente

• Décollement de la Couche Limite


 perte de portance

• Régime laminaire  décollement


plus rapide et probable

• Aéronautique  turbulence

Régime turbulente
17
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel
CHAPITRE V: Utilisation du logiciel

INFLUENCE DU
NOMBRE DE REYNOLDS
• Petits Re
• Décollement à couche
limite décollée de
l’extrados
A point de décollement
B point de transition
C zone morte
L écoulement laminaire
Régime sous-critique T écoulement turbulent

• Portance ↓
• Traînée ↑

18
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel
CHAPITRE V: Utilisation du logiciel

• Re élevés (>106) INFLUENCE DU


NOMBRE DE REYNOLDS
• Couche limite décollée 
atteint le bord de fuite

Re critique
A point de décollement
• Traînée ↓ B point de transition
D bulle décollé laminaire
L écoulement laminaire
• Portance ↑ T écoulement turbulent

• Re ↑ c.l. turbulente plus fine


 traînée ↓↓

Régime sur-critique
19
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel
CHAPITRE V: Utilisation du logiciel

LE NOMBRE DE
REYNOLDS CRITIQUE

• Re < Recritique  ↓ des


performances de la
pale

• Graphique: Recritique
pour profil FX-M2 =
80.000

Finesse maximal vs. Reynolds

20
CHAPITRE III
Base théorique du logiciel
21
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel

LA THÉORIE DU
MOMENT AXIAL DE FROUDE

• Modèle le plus simple: Actuator Disk Model

• But : démonstration de Vitesse rotor ≠ Vitesse de l’écoulement libre

• Disque actuateur  dispositif de traînée

22
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel
FACTEUR
D’INDUCTION AXIAL
• Bernouillli
• Équation du moment axial
• Conservation de la masse

U - u1 = a U
U - u2 = 2a U

• a = facteur d’induction axial


• Vmin écoulement libre = 0 
amax = 0.5
23
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel
COEFFICIENT
DE PUISSANCE
• 1ère loi thermodynamique

P = 0,5 ρ A (U 2- u22) u1 = 0,5 ρ A u1 (U + u2) (U - u2)

U - u1 = a U et U - u2 = 2a U

P 4U 3 a(1 − a ) 2
Cp = = 3 C p = 4a(1 − a) 2
1 U
ρU A
3

• Cp max= (16/27)=0,593

24
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel

CONSIDÉRATIONS
DU ROTATIONNEL
• Rotation dans le
sillage

• Equation d’Euler

• aΩ (b)= facteur
d’induction
rotationnel

• Graphique:

25
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel

LE TIP SPEED RATIO (TSR)


Ωr
• λ =
U ∞

• TSR  efficience
• TSR  niveau de bruit
• TSR  nécessité de
pales plus résistantes
• Courbe verte  éolienne
à axe horizontale
• TSR = 6/7  meilleur
rendement

26
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel

LA THÉORIE DE
L’ELEMENT PALE
• Complète l’étude: dessin de la
pale  forme, section, calage

• Modélisation bidimensionnelle
du rotor

• Intégration  la poussée, le
couple et la puissance

• Discrétisation du V.C  N
éléments d’épaisseur dr

– Éléments aillées
– Nombre de pales infini

27
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel
TRIANGLE
DES VITESSES

• Triangle des vitesses


d’un élément pale
– Vitesse axiale: V(1-a)

– Vitesse tangentielle:
Ωr(1+b)
(1 − a) V
tan φ = ⋅
(1 + b) Ωr

φ = β +α
Φ : angle du vent relatif avec le plan de

W = V +V 2
ax
2
tan
rotation [rad].
β: angle de calage du profil [rad]
α: angle d’incidence [rad].
28
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel
COEFFICIENTS
ET FORCES
• Coefficients de portance et traînée
Rz
Cd =
ρ
⋅ l ⋅V 2 ρ: densité du fluide
2
Rz: force en direction du flux
Rx Rx: force en direction normal au flux
Cl =
ρ V: vitesse du flux
⋅ l ⋅V 2
2

• Les composantes poussée, CT et tangentielle, CM


CT = Cl ⋅ cos φ + Cd ⋅ sin φ

C M = Cl ⋅ sin φ − C d ⋅ cos φ

29
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel DYNAMIQUE
DE LA PALE

• Les forces Intégration en long de la pale


ρ
dRz = Cl ⋅ ⋅ l ⋅ W 2 ⋅ dr R

2 FT = p ⋅ ∫ dFT
• Poussée axiale
ρ 0
dRx = Cd ⋅ ⋅ l ⋅ W 2 ⋅ dr
2
R
QT = p ⋅ ∫ r ⋅ dFT
• Couple de poussée
• Poussée 0

dFT = dRz ⋅ cos φ + dRx ⋅ sin φ FM = p ⋅ ∫ dFM


• Force motrice 0

• Force • Couple de moteur QM = p ⋅ ∫ r ⋅ dFM


0

dFM = dRz ⋅ sin φ − dRx ⋅ cos φ R


P = Ω ⋅ QM = p ⋅ Ω ⋅ ∫ r ⋅ dFM
• Puissance 0

• Couple
• Moment de flexion M f = QT2 + QM2
dQM = r ⋅ dM 30
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel
CALCUL DES
PARAMETRES a ET b
• Débit de masse  ρ ⋅ V (1 − a ) ⋅ 2πr ⋅ dr

• Variation de la quantité de mouvement 


ρ ⋅V (1 − a) ⋅ 2πr ⋅ dr ⋅ [V − V (1 − 2a)]

ρ
• Force appliquée par les p pales  p⋅
2
⋅ W 2 ⋅ l ⋅ CT ⋅ dr

V ⋅ (1 − a )
• Triangle des vitesses  W
= sin φ

a p ⋅ l ⋅ Ct
=
(1 − a ) 8πr ⋅ sin 2 φ

31
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel
CALCUL DES
PARAMETRES a ET b

• Moment de la quantité de mouvement 


ρ
ρ ⋅ V (1 − a ) ⋅ r ⋅ 2 π r ⋅ dr ⋅ 2 b Ω = ⋅ p ⋅W 2
⋅l ⋅C M ⋅ dr
2

ωr (1 + b)
= cos φ
• Triangle des vitesses  W

b p ⋅ l ⋅ CM
=
(1 + b) 4πr ⋅ sin 2φ

32
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel

CORRECTIONS DES
PARAMETRES a ET b
• Correction du nombre infini de
pales  Prandtl
• Nappes tourbillonnaires ≈
lignes parallèles à la distance d

2π ( r )
d= R ⋅ sin φ
t
p
avec φt = φ pour
r =1
R

• Facteur de réduction de la Visualisation du vortex crée en bout d’une


pale de 10 dans une soufflerie de la NASA
circulation de Prandtl
π ⋅ (1´−( r R) p 1 − ( R)
r
cos −1 [exp( − f )]
2 f = = ⋅
F= avec
π d 2 sin φt

33
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel

CORRECTIONS DES
PARAMETRES a ET b

• 2ème correction  contrôler les hauts


facteurs d’induction axiale

• acritique = 0,2

• Correction empirique si l’on est au-dessus


de acritique 1  a= 2 + K (1 − 2a c ) − ( K (1 − 2a c ) + 2) 2 + 4( Ka c − 1)
2

2  

4 F sin φ 2
K=
σCt

34
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel
CHAPITRE V: Utilisation du logiciel

FACTEURS a ET b CORRIGÉS
p ⋅l
• σ = solidité  σ =
2πr

• Prandtl aF ⋅ (1 − aF ) =
σ ⋅ CT
⋅ (1 − a ) 2
b σ ⋅ CM
=
4 ⋅ sin φ 2
(1 + b) 2 F ⋅ sin 2φ

−1
 4 F sin 2 φ 
• Nouvelles valeurs a corr = 1 +
σCT 
si a ≤ 0,2

1
acorr = 2 + K (1 − 2ac ) − ( K (1 − 2a c ) + 2) 2 + 4( Ka c − 1)  >
2 si a 0,2
2  
−1
 4 F sin φ ⋅ cos φ 
bcorr = − 1
 σ C M  35
CHAPITRE IV
Introduction au logiciel
36
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel

INTRODUCTION
AU LOGICIEL
• Le but du logiciel:
courbes de puissance
et rendement
• Données initiales

• Choix 1: type de
régulation
• Choix 2: type de profil
• Choix 3: longueur de
corde
• Choix 4: angles de
calage

37
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel

LE BUT DU LOGICIEL:
COURBES DE PUISSANCE
ET RENDEMENT

• Puissance récupérée Puissance en kW Puissance et Rendement


Coefficient de
puissance
vs. Vitesse du vent 400 1
Puissance totale
Coeficient de puissance 0,9
350

• Freinage à la 300
0,8

Puissance Nominale 250


0,7

(Pnom) 200
0,6

0,5

0,4
• Rendement vs. 150

0,3

Vitesse du vent 100


0,2

50
0,1

• Rendement maximal à 0
0 5 10 15 20 25
0

la vitesse moyenne Vitesse du Vent en m/s

38
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel DONNEES
CHAPITRE V: Utilisation du logiciel
INITIALS

• p : nombre de pales de
l’éolienne
• ρ : densité de l’air
• V : vitesse du vent
moyenne au lieu choisi
• Pnom : Puissance nominale
de la génératrice de l’éolienne
(ou R : rayon de l’éolienne)
• λ : rapport entre vitesse du
vent et le productif de vitesse
angulaire et rayon
• η : rendement global de
l’éolienne

39
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel
CALCUL DE LA
PUISSANCE ET LE RAYON

• Vdem, Vnom et Pnom  assimilation de la courbe à une droite


• Prec à la vitesse moyenne (approximation)

8 P
P = η· ·ρ ·π ·V 3 ·R 2 R=
8
27 η· ·ρ ·π ·V 2
27
• Puissance récupérée à partir du rayon
• R  Préc
• Vdem, Vnom et Prec  Pnom

40
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel

CHOIX 1:
TYPE DE RÉGULATION

• Objectif: freiner à la PNOM de la


génératrice
• Pitch: profit maximal de
l’énergie  control externe
• ηhélice = maximal = 0,85
• Stall: décrochage
aérodynamique
• ηhélice = non maximal
41
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel

CHOIX 2:
TYPE DE PROFIL

• Profil NACA 4412


– On connaît les coefficients
de portance, Cz, et traînée,
Cx, pour chaque angle
d’incidence, α

• Profil NACA 63215


– On ne connaît que les Cx
et Cz pour les α jusqu’à
environ 30º Viterna-
Corrigan
• Coefficients introduits par
l’utilisateur
42
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel CHOIX
CHAPITRE V: Utilisation du logiciel
DE PROFIL

• NACA 63215
• Problème: essais insuffisants de soufflerie
• Objectif: Cl e Cd à la zone post stall Viterna
Corrigan
Cl/Cd  cos 2 α  Cd = B1 sin 2 α + B2 cos α
C l = A1 ⋅ sin α + A2 
2

 sin α 

A1/B1 B1 B1 = Cdmáx
A1 =
2 C dmáx = 1,11 + 0,018 ⋅ AR
A2/B2  sinα   1 
A2 = (Cls − Cdmáx ⋅ sinαs ⋅ cosαs ) ⋅  2 s  B2 =   ⋅ (Cds − C dmáx ⋅ sin 2 α s )
 cos αs   cosα s 
αs =15° Cds = 0.075
43
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel
CHAPITRE V: Utilisation du logiciel

CHOIX 3:
LONGUEUR DE CORDE
1. Longueurs optimales  la
méthode de Glauert
2. Longueurs introduites par
l’utilisateur

CHOIX 4: ANGLES DE CALAGE

1. Calage optimal  la méthode


de Glauert

2. Calage introduit par l’utilisateur

3. Calage à partir d’une droite 


β= AR+B 44
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel
CALAGES ET
LONGUEURS OPTIMAUX
(1 − 3 ⋅ a )
b=
(4 ⋅ a − 1)
x = (4 ⋅ a − 1) ⋅ (1 − a) (1 − 3 ⋅ a)
b ⋅ (1 + b) V2 1
= 2 2 = 2
a ⋅ (1 − a) Ω ⋅ r x

• Miller  série  1   2   10   418 


a ≅  − + − + .....
2  4  6 
approximative  3   81 ⋅ x   729 ⋅ x   59049 ⋅ x 

• Φ tan φ = (1 − a) ⋅ (1 − a) (1 − 3 ⋅ a )

• Corde optimal p⋅Ω  (4 ⋅ a − 1)


⋅ (l ⋅ Cl ) =   ⋅ (1 − a ) ⋅ (1 − 3 ⋅ a )
8 ⋅ π ⋅V  (1 − 2 ⋅ a )
45
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Base théorique du logiciel
CHAPITRE IV: Introduction au logiciel
CHAPITRE V: Utilisation du logiciel
CALAGES ET
LONGUEURS OPTIMAUX
• Calcul du rotor
optimal  Glauert
120

• Supposition: traînée 100

nulle  Cd=0 80

• Calage optimal 60

Finesse
Serie1
40

20

β = φ −α -10
0
0 10 20 30 40 50 60 70

-20
angle d'incidence
α dont la finesse est
maximale
• Corde optimal
p⋅Ω  (4 ⋅ a − 1)
⋅ (l ⋅ Cl ) =   ⋅ (1 − a) ⋅ (1 − 3 ⋅ a)
• Cl  coefficient de 8 ⋅ π ⋅V  (1 − 2 ⋅ a )
portance à α optimal

46
ALGORITHME
MATEMATIQUE

1)Initialisation de
2) Angle du vent relatif Φ
a et b  a=b=0
 (1 − a ) ⋅ V0 
φ = tan −1  
 (1 + a ) ⋅ ω ⋅ r 

8) Test de convergence
3) Angle d’incidence α
anouveau − aancien < ξ
α =φ −β
bnouveau − bancien < ξ

7) Nouveaux a et b 4) Cl et Cd

6) Calcul du facteur F 5) Coefficients des forces


Cn = Cl .cos φ + Cd .sin φ
Ct = Cl .sin φ − Cd .cos φ
47
CHAPITRE V
Utilisation du logiciel
48
CHAPITRE I: L’éolienne et l’énergie éolienne
CHAPITRE II: Introduction à l’aérodynamique
CHAPITRE III: Introduction au logiciel
CHAPITRE IV: Base théorique du logiciel
CHAPITRE V: Utilisation du logiciel

EXEMPLE DE
PROGRAMMATION
Sub iteration()

Dim i As Integer
'Itérations
i=0 'stockage des valeurs de a et b
While (i < 100 And Range("L20:M29").Select
Sheets("Feuille1").Range("J13 Selection.Copy
").Value <> 0) Range("H3").Select
'copie des valeurs des profils Selection.PasteSpecial
(feuille Polaires ou Paste:=xlValues,
Polaires_ext en fonction du
choix) Operation:=xlNone,
SkipBlanks:=False,
If(Sheets("Polaires_ext").Rang Transpose:=False
e("M6").Value <> 1) Then
'polaires viterna i=i+1
Call Coeffs_viterna Wend
Else 'polaires ext End Sub
Call coeffs_ext
End If
49
VÉRIFICATION
DES RESULTATS
Enercon E32/300

• Type régulation : Stall


• Longueur de la pale : 16 m
• Puissance nominale : 300 kW
• Vitesse moyenne : 8 m/s
• Vitesse démarrage : 4 m/s
• Vitesse nominal : 14 m/s
• Rendement = 0.7
• λ=7

50
VÉRIFICATION
DES RESULTATS
Fortis - Montana

• Type régulation : Stall


• Longueur de la pale : 16 m
• Puissance nominale : 5800 W
• Vitesse moyenne : 8 m/s
• Vitesse démarrage : 3 m/s
• Vitesse nominal : 14 m/s
• Rendement = 0.7
• λ=7

51
CONCLUSIONS ET
POSSIBLES AMÉLIORATIONS

• Choix de plus de profils

• Changement de logiciel

– temps de résolution

– qualité du dessin

– possibilité de apprendre un logiciel plus nouveau

– Exemples: Matlab

• Approfondir la recherche pour la régulation pitch

52
CONCLUSIONS ET
POSSIBLES AMÉLIORATIONS

• La recherche: Reynolds, Mach

• La théorie mathématique

• L’expérience avec la petite éolienne

• Étude économique

Ce programme n’a pas la vocation de


dimensionner précisément une éolienne, il s’agit
d’un outil pratique de pré dimensionnement et
d’analyse des interactions entre les différents
paramètres de modélisation
53
MERCI PAR
VOTRE ATTENTION!
54

Vous aimerez peut-être aussi