Projet Eolienne

Rapport Travaux Pratiques Résolution de problèmes énergétiques
Projet d’éolienne
Détermination de la vitesse de rotation d’une éolienne

Filière : Licence 3 « Mécanique génie thermique »
Présentée par : Encadrée par :
EL AZZOUZI Amal Thibaut Menard
1
2
Table des matières
1 Introduction : .................................................................................................... 4
1.1 L’objectif : .................................................................................................. 4
1.2 L’influence de l’éolienne sur l’air : ................................................................ 4
2 Démarche de travail : ......................................................................................... 5
2.1 Détermine toutes les constants du problème ............................................... 5
2.2 Crée les fonctions d’interpolation pour le cx et le cz ...................................... 5
2.3 Créez la fonction principale f(r, w,a) et les fonctions Va(r, ω, a), β(r, ω, a), αv(r),
i(r, ω, a) ................................................................................................................ 5
2.4 Créez la fonction principale h(a , w) et la fonction pour Caero(ω, a) ............... 6
2.5 Calculer la puissance récupérée par l’éolienne ............................................ 7
3 Résultat final : ................................................................................................... 7
3.1 Explication : ............................................................................................... 7
4 Conclusion : ..................................................................................................... 9
3
1 Introduction :
1.1 L’objectif :
L’objectif de ce TP est de créer un programme permettant de déterminer la vitesse de rotation
d’une éolienne dont la forme des pâles est définie par un profil NACA. Ce programme pourra
servir à l’étude de différentes configurations de vrillage des pales et de taille de corde,
influençant la production d’énergie lors de projet. Il a noté que cette modélisation est
simplifiée et prend peu en compte l’influence de l’éolienne sur l’air.
1.2 L’influence de l’éolienne sur l’air :

Tout comme l’énergie géothermique s’appuie sur la chaleur de la Terre, l’énergie éolienne est
obtenue grâce à la puissance de l’air en mouvement aussi appelée force cinétique du vent.
Lorsque le vent souffle, les pales de l’éolienne tournent et entraînent un générateur qui
transforme l’énergie mécanique en énergie électrique.
Avec la seule force du vent, l’hélice de l’éolienne atteint les 10 à 25 tours par minute autour
du rotor, ce qui n’est pas suffisant pour produire de l’électricité sous forme d'énergie éolienne
qui est de fait une énergie renouvelable. Un multiplicateur va alors permettre de multiplier la
vitesse de rotation de l’hélice pour atteindre la vitesse nécessaire de 1000 tours/minute.
• Fp : force potence
• Ft : force de trainé
• F : force qui fait tourne la pale
• Va : vent apparent
Il est important de noter que les effets spécifiques de l'éolienne sur l'air varient en fonction de
nombreux facteurs tels que la localisation, la conception de l'éolienne, les conditions
4
environnementales locales et les efforts de planification pour minimiser les impacts négatifs.
Des études d'impact environnemental et des mesures de mitigation sont souvent menées avant
l'installation d'éoliennes pour minimiser les conséquences potentielles sur l'air et
l'environnement.
2 Démarche de travail :
Dans cette partie en se passe sur Spyder pour déterminer la vitesse de rotation d'une éolienne
en répondant aux questions suivantes qui nous aidons à la programmation :

2.1 Détermine toutes les constants du problème
• Vitesse du vent Vz = 15m.s−1
• Longueur des pales : Lp = 1.5m (= rmax − rmin avec rmin = 0.3m)
• Lp = rmax-rmin , rmax = Lp+rmin
• Longueur de la corde : Lc = 0.2m
• Coefficient de frottement du à l’alternateur A = 15N.m.s.rad−1
• Les polaires de la pale (profil NACA - Cx et Cz en fonction de l’angle d’incidence i)
sont données.
• Angle de callage en début de pale : αdv = 60°
• Angle de callage en fin de pale : αfv = 60°
• Masse volumique de l’air : ρ = 1.2kg.m−3
• Nombre de pâles : Npa = 3
2.2 Crée les fonctions d’interpolation pour le cx et le cz

Application : interpolate : Méthode pour la recherche de zéro
2.3 Créez la fonction principale f(r, w,a) et les fonctions Va(r, ω, a), β(r, ω, a), αv(r), i(r, ω,
a)
f(r, ω, a) = (−Cx(i(r, ω, a))cos(β(r, ω, a)) + Cz (i(r, ω, a))sin(β(r, ω, a))) 1/2 ρ.Lc.Va^2 (r, ω, a
5
Trainé potence
(1−𝑎)𝑉𝑧
✓ 𝛽 = 𝑎𝑟𝑐𝑠𝑖𝑛 ( )
𝑣𝑎
✓ 𝑣𝑎 = √((1 − 𝑎)2 𝑉22 + (𝑟𝑤)2 )

𝑟−𝑟𝑚 𝑛
✓ 𝛼𝑣 = 𝛼𝑣𝑑 + 𝑟 𝑚𝑎𝑥 −𝑟𝑖 ⋅ (𝛼𝑣𝐹 − 𝛼𝑣𝑑 )
𝑚𝑖𝑛
✓ 𝑖 = 𝛽 − 𝛼𝑣
Application :
2.4 Créez la fonction principale h(a , w) et la fonction pour Caero(ω, a)

Chercher cette vitesse finale, revient à se ramener à une recherche de zéro d’une fonction h(a,
ω) définie par:
2𝜌𝑆𝑎 (1 − 𝑎) ⋅ 𝑉23 − 𝐴𝑤𝑓2

h(a, ω) = {
𝑁𝑥 ⋅ 𝐶𝑎𝑒𝑟𝑜(𝑤𝐹, 𝑎) − 𝐴|𝑤𝑓|
La détermination du couple aérodynamique Caero sera fonction des coefficients de portance

et de traînée dépendant de l’angle entre la pâle et le vent apparent (angle qui dépend de la
position r sur la pâle). Cela nécessitera d’intégrer les forces sur toute la pâle pour déterminer
le couple aérodynamique.
𝑟𝑚𝑎𝑥
𝐶𝑎𝑒𝑟𝑜 (𝑤, 𝑎) = ∫ 𝑟𝑓(𝑟, 𝑤, 𝑎) ⅆ𝑟
𝑟𝑚𝑖𝑛
Application :
6
2.5 Calculer la puissance récupérée par l’éolienne
La relation de puissance :
𝑃 = 𝑁𝑝 Caero(𝑤𝑓, 𝑎) × 𝑤𝑓 = 𝐴𝑤 2 𝑓
Application :
3 Résultat final :
3.1 Explication :
• import matplotlib.pyplot as plt :Matplotlib permet de créer des graphiques, des

diagrammes et des visualisations de données de manière efficace et personnalisable.
L'import de cette bibliothèque est généralement réalisé avec l'alias plt pour simplifier
l'utilisation des fonctions de Matplotlib.
• Import timeit : En Python, le module timeit est utilisé pour mesurer le temps d'exécution
de petits fragments de code. Il fournit une méthode simple pour mesurer les performances
relatives de différentes approches pour résoudre un problème, en mesurant le temps
nécessaire à l'exécution d'un bout de code ou d'une fonction.
• Loodtxt : utilise pour présente (texte, tableau …)
• Np.linspace : utilise pour présente les points dans l’intervalle
• ‘x’or’*’ : c’est comment présenté les points dans le graphe
7
8
4 Conclusion :
Ce travail pratique visait à créer un programme utilisant le profil NACA pour déterminer la
vitesse de rotation d'une éolienne. Ce programme est destiné à faciliter l'étude de diverses
configurations de vrillage des pales et de taille de corde, des paramètres cruciaux influençant
la production d'énergie dans des projets éoliens. Il est important de noter que cette
modélisation est simplifiée et ne prend que partiellement en compte l'influence réelle de
l'éolienne sur l'air.
En utilisant la programmation Python, nous avons développé un programme capable de

calculer la vitesse de rotation de l'éolienne en fonction des paramètres du profil NACA, tels
que la puissance récupérée . Ce programme offre la possibilité d'évaluer différentes
configurations de conception pour optimiser la production d'énergie éolienne dans des
conditions idéales.
En somme, ce travail pratique illustre la capacité de la programmation Python à créer des

outils de modélisation et de simulation pour étudier et optimiser les performances des
éoliennes, même si cette modélisation est simplifiée et ne représente qu'une approximation
des résultats réels dans des conditions réelles d'utilisation.

Projet Eolienne

Transféré par

Droits d'auteur :

Formats disponibles

Projet Eolienne

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Projet Eolienne

Transféré par

Droits d'auteur :

Formats disponibles

Rapport Travaux Pratiques Résolution de problèmes énergétiques

Détermination de la vitesse de rotation d’une éolienne

Présentée par : Encadrée par :

EL AZZOUZI Amal Thibaut Menard

1.2 L’influence de l’éolienne sur l’air :

en répondant aux questions suivantes qui nous aidons à la programmation :

• Vitesse du vent Vz = 15m.s−1

• Longueur des pales : Lp = 1.5m (= rmax − rmin avec rmin = 0.3m)

• Lp = rmax-rmin , rmax = Lp+rmin

• Longueur de la corde : Lc = 0.2m

• Coefficient de frottement du à l’alternateur A = 15N.m.s.rad−1

• Les polaires de la pale (profil NACA - Cx et Cz en fonction de l’angle d’incidence i)

• Angle de callage en début de pale : αdv = 60°

• Angle de callage en fin de pale : αfv = 60°

• Masse volumique de l’air : ρ = 1.2kg.m−3

• Nombre de pâles : Npa = 3

2.2 Crée les fonctions d’interpolation pour le cx et le cz

✓ 𝑣𝑎 = √((1 − 𝑎)2 𝑉22 + (𝑟𝑤)2 )

2.4 Créez la fonction principale h(a , w) et la fonction pour Caero(ω, a)

2𝜌𝑆𝑎 (1 − 𝑎) ⋅ 𝑉23 − 𝐴𝑤𝑓2

La détermination du couple aérodynamique Caero sera fonction des coefficients de portance

• import matplotlib.pyplot as plt :Matplotlib permet de créer des graphiques, des

En utilisant la programmation Python, nous avons développé un programme capable de

En somme, ce travail pratique illustre la capacité de la programmation Python à créer des

Vous aimerez peut-être aussi