Dimensionnement D Un Rotor D Eolienne
Dimensionnement D Un Rotor D Eolienne
Dimensionnement D Un Rotor D Eolienne
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
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
• Bernouilli et Newton
( )
.
m v1 − v 2
2 2
P=
2
P=
ρ
4
(v
1
2
− v2
2
)(v + v )S
1 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
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
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
• Portance
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
• 4 points remarquables
C
• Finesse f = z
C x
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!
µ υ 20 °C]
µ est la valeur de la viscosité [1,81 10-5 m2s-1 à
20 °C].
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
• 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
DECOLLEMENT
DE LA COUCHE LIMITE
• 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 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
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
• Graphique: Recritique
pour profil FX-M2 =
80.000
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
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
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
• 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
– 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
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
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
• 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
ρ
• 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
ω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
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
• acritique = 0,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
• Freinage à la 300
0,8
(Pnom) 200
0,6
0,5
0,4
• Rendement vs. 150
0,3
50
0,1
• Rendement maximal à 0
0 5 10 15 20 25
0
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
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
CHOIX 2:
TYPE 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
• Φ tan φ = (1 − a) ⋅ (1 − a) (1 − 3 ⋅ a )
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
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
50
VÉRIFICATION
DES RESULTATS
Fortis - Montana
51
CONCLUSIONS ET
POSSIBLES AMÉLIORATIONS
• Changement de logiciel
– temps de résolution
– qualité du dessin
– Exemples: Matlab
52
CONCLUSIONS ET
POSSIBLES AMÉLIORATIONS
• La théorie mathématique
• Étude économique