République Algérienne Démocratique et Populaire

Ministère de l'Enseignement Supérieur et de la Recherche Scientifique

Ecole Nationale Supérieure Polytechnique

Département de Génie Mécanique

Projet de Fin d'Etudes

Pour l'obtention du diplôme
D'Ingénieur d'Etat en Génie Mécanique

Thème

Dimensionnement et conception dune éolienne à

axe horizontal destinée à être installée sur le toit
du bâtiment vert de l’ENSP

Proposé et dirigé par: Présenté par :

Dr. A. SMAILI HADBI Ramdane Azeddine

Mr [Link]

Promotion 2009
A mes parents et toute ma famille
A mes très chers camarades
A mes enseignants
 Remerciements

Je remercie DIEU tout puissant qui m’a donne le

courage, la force et la volonté pour réaliser ce modeste

travail.

Je tiens à exprimer ma reconnaissance à mon promoteur

Monsieur SMAILI Arezki qui m’a donné l’occasion de

travailler sur un sujet passionnant, et pour ses conseils

judicieux et aide précieuse.

Je tiens egalement à remercier mon co-promoteur,

monsieur RIACHE Lakhdar, pour son encadrement, sa

disponibilité ainsi que ses conseils judicieux.

Je remercie Mr RECHAK et Mr BELHADEF membres du

jury qui ont accepté d’examiner mon travail

Je tiens enfin à remercier également mes camarade

pour leurs disponibilité et encouragements.

 ‫ﻣﻠﺨﺺ‬

‫ ﯾﺘﻢ اﻧﺠﺎز‬. ‫ﯾﮭﺪف ھﺪا اﻟﻌﻤﻞ إﻟﻰ ﺗﺼﻮر ﺗﻮر ﺑﯿﻨﺔ رﯾﺎح ذات ﻣﺤﻮر أﻓﻘﻲ ﺑﺎﻻﺳﺘﻌﺎﻧﺔ ﺑﺎﻟﺒﺮﻧﺎﻣﺞ اﻟﺤﺎﺳﻮﺑﻲ ﺳﻮﻟﯿﺪ وورﻛﺲ‬
‫ اﻷوﻟﻰ ﺗﺘﻤﺜﻞ ﻓﻲ ﺟﻤﻊ ﻣﻌﻠﻮﻣﺎت ﺣﻮل ﻣﺨﺘﻠﻒ اﻟﺘﻘﻨﯿﺎت اﻟﻤﺴﺘﻌﻤﻠﺔ ﻓﻲ اﻻﻧﺠﺎزات اﻟﺴﺎﺑﻘﺔ ﻟﺘﺮﺑﯿﻨﺎت‬.‫اﻟﻌﻤﻞ وﻓﻖ ﺛﻼث ﻣﺮاﺣﻞ‬
‫اﻟﻤﺮﺣﻠﺔ اﻟﺜﺎﻧﯿﺔ ھﻲ دراﺳﺔ آﯾﺮودﯾﻨﺎﻣﯿﻜﯿﺔ ﻟﻠﺪوار اﻟﺘﻲ ﺑﺈﻣﻜﺎﻧﻨﺎ ﻣﻦ ﺧﻼﻟﮭﺎ ﺗﺤﺪﯾﺪ ﺳﺮﻋﺖ اﻟﺮﯾﺎح اﻟﺪﻧﯿﺎ اﻟﺘﻲ ﺗﺒﺪآ ﻗﯿﮭﺎ‬. ‫اﻟﺮﯾﺎح‬
‫ آﺧﺮ ﻣﺮﺣﻠﺔ ﺗﺘﻤﺜﻞ ﻓﻲ ﺣﺴﺎب دﯾﻨﺎﻣﯿﻜﻲ ﻟﻸﺟﺰاء اﻟﻤﻜﻮﻧﺔ ﻟﻠﺘﺮ ﺑﯿﻨﺔ ﺑﺎﻻﺳﺘﻌﺎﻧﺔ ﺑﺒﺮﻧﺎﻣﺞ‬. ‫اﻟﺘﺮﺑﯿﻦ ﺑﺘﻮﻟﯿﺪ اﻟﻄﺎﻗﺔ اﻟﻜﮭﺮﺑﺎﺋﯿﺔ‬
.‫ﻛﻮﺳﻤﻮس وورﻛﺲ ﻟﻠﺘﺤﻘﻖ ﻣﻦ ﺻﺤﺔ اﻟﻔﺮﺿﯿﺎت اﻟﻤﺴﮭﻠﺔ ﻟﻠﺤﺴﺎب اﻟﯿﺪوي‬

‫ﻛﻠﻤﺎت ﻣﻘﺘﺎﺣﯿﺔ‬

‫ ﺣﺴﺎب اﯾﺮودﯾﻨﺎﻣﯿﻜﻲ – دراﺳﺔ رﻗﻤﯿﺔ‬- ‫ﺗﻮر ﺑﯿﻨﺔ رﯾﺎح‬- ‫ﺣﺴﺎب دﯾﻨﺎﻣﯿﻜﻲ‬

Résumé

Le but de ce travail est la conception d’une eolienne à axe horizontale fonctionnel à

l’aide du logiciel SoidWorks. La conception se fera en trois grandes étapes. La
première consiste à rassembler des informations sur les différentes technologies
éoliennes. La seconde est une étude aérodynamique de l’hélice permettant de
déterminer la vitesse du vent minimal de production d’électricité et la courbe de
puissance en fonction de la vitesse du vent. La dernière étant un calcul dynamique
des organes constitutifs de la turbine éolienne et simulations par le logiciel
CosmosWorks pour obtenir la forme final de la conception.

Mots clés

Turbine éolienne, calcul aérodynamique, simulation numérique, calcul dynamique.

Abstract

The purpose of this work is the conception of an horizontal axis wind turbine,using
SoidWorks. The conception will be made in three important steps. The first one,
consists in gathering information on the various wind technologies. The second is an
aerodynamic study of the rotor allowing to determine the munimum wind speed of
electricity production and the curve of power as function of the wind speed. The last
one being a dynamic calculation of the constituent parts of the wind turbine and the
simulations by the software SolidWorks to obtain the final shape of the conception.

Keywords

Wind turbine, Aerodynamic calculation, Dynamic calculation, numerical simulation

 Table des matières
INTRODUCTION....................................................................................................................................... 1
CHAPITRE 1 : GENERALITES .................................................................................................................... 4
1.1 Historique : .................................................................................................................................... 4
1.2 Types et caractéristiques des éoliennes.................................................................................... 5
1.2.1 Eolienne à axe vertical........................................................................................................ 5
[Link] Rotor Savonius............................................................................................................. 5
[Link] Rotor Darrieus ............................................................................................................. 6
[Link] Caractéristiques........................................................................................................... 7
1.2.2 Eolienne à axe horizontal ................................................................................................... 7
[Link] Eolienne face au vent « upwind » ............................................................................... 7
[Link] Eolienne sous le vent « downwind »........................................................................... 8
1.3 Régulation.................................................................................................................................. 8
1.3.1 Eolienne à contrôle passif ( stall control ) ........................................................................ 10
1.3.2 Eolienne à contrôle actif « pitch control » ....................................................................... 11
1.4 Aérodynamique des éoliennes ................................................................................................ 13
1.4.1 Théorie de Betz................................................................................................................. 13
[Link] Hypothèses................................................................................................................ 13
[Link] Facteur d’induction axiale ......................................................................................... 16
[Link] Coefficient de poussée et de puissance .................................................................... 15
[Link] Limite de Betz............................................................................................................ 16
1.4.2 Méthode BEM................................................................................................................... 17
[Link] Théorie de l’élément de pale .................................................................................... 17
[Link] Bilan de quantité de mouvement appliqué à l’élément de pale............................... 20
[Link] Couplage théorie de quantité de mouvement et élément de pale .......................... 20
1.4.3 Pertes en bous de pales.................................................................................................... 21
1.4.4 Model de décrochage de Laino & Hansen........................................................................ 22
CHAPITRE 2 : CONCEPTION AERODYNAMIQUE ................................................................................... 24
2.1 Design du rotor............................................................................................................................ 24
2.1.1 Paramètres liés au rotor....................................................................................................... 24
2.1.2 Géométrie de la pale ............................................................................................................ 27
2.1.3 Calcul des performances et modifications du design ........................................................... 26
 2.3 Résultats et validations des codes de calculs .............................................................................. 30
2.3.1 Programme de géométrie de la pale.................................................................................... 30
2.3.2 Courbes de puissance........................................................................................................... 33
CHAPITRE 3 : LOGICIEL DE SIMULATION SOLIDWORKS....................................................................... 37
3.1 Qu’est ce que le logiciel SolidWorks ?......................................................................................... 37
3.2 Potentiel de SolidWorks .............................................................................................................. 37
3.2.1 Création de formes volumiques ........................................................................................... 37
3.2.2 Effectuer des assemblages ................................................................................................... 36
3.2.3 Conception, analyse et simulation ....................................................................................... 36
3.3 Exemple d’utilisation de solidworks............................................................................................ 38
CHAPITRE 4 : CONCEPTION MECANIQUE............................................................................................. 43
4.1 Schéma cinématique et organes constitutifs .............................................................................. 43
4.2 Etude technologique .................................................................................................................... 44
4.3 Calcul dynamique........................................................................................................................ 48
4.3.1 Les pales ............................................................................................................................... 46
4.3.2 L’arbre................................................................................................................................... 49
4.3.3 Montage pales-moyeu ......................................................................................................... 51
[Link] Calculs du diamètre de la tige ....................................................................................... 51
[Link] Calcul du couple de serrage .......................................................................................... 52
[Link] Résistance à la traction, cisaillement et matage ........................................................... 53
4.3.4 Montage moyeu-arbre ......................................................................................................... 54
[Link] Montages par clavette................................................................................................... 54
[Link] Montages par cannelure ............................................................................................... 55
4.3.5 La tour................................................................................................................................... 57
[Link] Résistances à la flexion-compression............................................................................ 57
[Link] Vérifications au flambement......................................................................................... 58
[Link] Calculs du diamètre des câbles de soutien ................................................................... 59
4.3.6 Les accouplements ............................................................................................................... 60
[Link] Descriptions................................................................................................................... 61
[Link] Fonctionnement ............................................................................................................ 61
[Link] Détermination de l’accouplement ................................................................................ 61
4.3.7 La génératrice....................................................................................................................... 63
[Link] Puissance conventionnelle de la machine..................................................................... 63
[Link] Diamètre interne du stator ........................................................................................... 63
 [Link] Diamètre externe du stator........................................................................................... 63
[Link] Calcul du pas polaire ..................................................................................................... 64
[Link] Longueur de la machine ................................................................................................ 64
4.3.8 Le multiplicateur................................................................................................................... 65
[Link] Calcul des diamètres des arbres.................................................................................... 65
[Link] Détermination des modules des engrenages................................................................ 65
4.3.9 Calcul du frein....................................................................................................................... 66
4.4 Vérification des calculs dynamiques par solidworks ............................................................... 69
4.4.1 La tour............................................................................................................................... 69
4.4.2 L’arbre............................................................................................................................... 70
4.4.3 La pale............................................................................................................................... 72
CONCLUSION......................................................................................................................................... 76
BIBLIOGRAPHIE ..................................................................................................................................... 77
ANNEXES ............................................................................................................................................... 79
INTRODUCTION

CHAPITRE 1 : GENERALITES
1.1 historique :

1.2 types et caractéristiques des éoliennes

1.2.1 éolienne à axe vertical

[Link] rotor Savonius

[Link] Rotor Darrieus

[Link] Caractéristiques

1.2.2 éolienne à axe horizontale

[Link] éolienne face au vent « upwind »

[Link] éolienne sous le vent « downwind »

1.3 Régulation

1.3.1 Éolienne à contrôle passif ( stall control )

1.3.2 Éolienne à contrôle actif « pitch control »

1.4 Aérodynamique des éoliennes

1.4.1 théorie de Betz

[Link] hypothèses

[Link] facteur d’induction axiale

[Link] coefficient de poussée et de puissance

[Link] limite de Betz

1.4.2 Méthode BEM

[Link] Théorie de l’élément de pale

[Link] Bilan de quantité de mouvement appliqué à l’élément de pale

[Link] Couplage théorie de quantité de mouvement et élément de pale

1.4.3 Pertes en bous de pales

1.4.4 Model de décrochage de Laino & Hansen

CHAPITRE 2 : CONCEPTION AERODYNAMIQUE

2.1 Design du rotor
2.1.1 Paramètres liés au rotor

2.1.2 Géométrie de la pale

2.1.3 Calcul des performances et modifications du design

2.2 Organigrammes de calcul

2.2.1 Organigramme de déterminations de la géométrie

2.2.2 Organigramme d’évaluation de la puissance

2.3 Résultats et validations des codes de calculs

2.3.1 Programme de géométrie de la pale

2.3.2 Courbes de puissance

CHAPITRE 3 : LOGICIEL DE SIMULATION SOLIDWORKS

3.1 Qu’est ce que le logiciel SolidWorks ?

3.2 Potentiel de SolidWorks

3.2.1 Création de formes volumiques

3.2.2 Effectuer des assemblages

3.2.3 Conception, analyse et simulation

3.3 Exemple d’utilisation de solidworks

CHAPITRE 4 : conception mécanique

4.1 Schéma cinématique et organes constitutifs

4.2 Etude technologique

4.3 Calcul dynamique

4.3.1 Les pales

4.3.2 L’arbre

4.3.3 Montage pales-moyeu

[Link] calculs du diamètre de la tige

[Link] Calcul du couple de serrage

[Link] résistance à la traction, cisaillement et matage

4.3.4 Montage moyeu-arbre

[Link] montages par clavette

[Link] montages par cannelure

4.3.5 La tour

[Link] résistances à la flexion-compression

[Link] vérifications au flambement

[Link] calculs du diamètre des câbles de soutien

4.3.6 Les accouplements

[Link] descriptions

[Link] fonctionnement

[Link] détermination de l’accouplement

4.3.7 la génératrice

[Link] Puissance conventionnelle de la machine

[Link] Diamètre interne du stator

[Link] Diamètre externe du stator

[Link] Calcul du pas polaire

[Link] Longueur de la machine

4.3.8 le multiplicateur

[Link] Calcul des diamètres des arbres

[Link] Détermination des modules des engrenages

4.3.9 calcul du frein

4.4 vérification des calculs dynamiques par solidworks

4.4.1 La tour

4.4.2 L’arbre

4.4.3 La pale

CONCLUSION

BIBLIOGRAPHIE

ANNEXES
 Liste des figures
Figure 1.1: Effet du vent sur un corps creux ........................................................................................... 5
Figure 1.2: Principe de l’incidence variable............................................................................................. 6
Figure 1.3 : technologies d’éoliennes à axe vertical................................................................................ 6
Figure 1.4: Diagramme de la puissance utile sur l'arbre en fonction de la vitesse du vent.................... 9
Figure 1.5: Ecoulement d’air sur une pale en régime « normal » (à gauche) et en situation de
décrochage (à droite) ............................................................................................................................ 10
Figure 1.6: Comparaison des systèmes pitch à une vitesse de rotation et stall à deux vitesses Fixe... 11
Figure 1.7: Effet de l’angle de pas sur le coefficient aérodynamique de puissance éolienne .............. 12
Figure 1.8: Modélisation du rotor éolien par un disque actuateur....................................................... 13
Figure 1.9: Paramètres opératoires d’une turbine Betz........................................................................ 16
Figure 1.10: Subdivision de la pale ........................................................................................................ 17
Figure 1.11: Profil d’une pale d’éolienne à axe horizontale .définition des variables .......................... 18
Figure 2.1 : Variation de la puissance en fonction du coefficient spécifique de vitesse....................... 25
Figure 2.2: Organigramme de détermination de la géométrie…………………………………………………………..28

Figure 2.3: Organigramme d'évaluation de la puissance …………………………………………………………………29

Figure 2.4 : Evolution de la corde et l’angle relatif en fonction du rayon : comparaison avec les
travaux antérieurs ................................................................................................................................. 32
Figure 2.5 : Courbe de la puissance....................................................................................................... 33
Figure 3.1 : Obtention d’un cylindre par révolution d’un rectangle ..................................................... 35
Figure 3.2 : Obtention d’un cylindre par extrusion d’une esquisse de cercle....................................... 36
Figure 3.3: Représentation des 6 premiers plans correspondants aux 6 premiers profils ................... 38
Figure 3.4 : Esquisse du premier points ................................................................................................ 39
Figure 3.5 : Forme finale du profil à la racine avec angle de vrillage .................................................... 39
Figure3. 6 : Lissage entre les deux premiers profils .............................................................................. 40
Figure 3.7 : Forme finale de la partie active de la pale ......................................................................... 40
Figure 3.8 : Support de la pale .............................................................................................................. 41
Figure 3.9: Forme finale de la pale ........................................................................................................ 42
Figure 4.1: Schéma cinématique de l’éolienne à axe horizontale......................................................... 43
Figure 4.2: Types de support de nacelles .............................................................................................. 44
Figure 4.3: Moyeu du rotor de l'éolienne ............................................................................................. 45
Figure 4.4: Profils de la pale avant la couche de résine ........................................................................ 45
Figure 4.5: Vue 3D de la pale................................................................................................................. 47
Figure 4.6: Modélisation du chargement sur la pale............................................................................. 48
Figure 4.7:Chargement de l'arbre ......................................................................................................... 49
Figure 4.8: Schématisation du chargement sur sur l'arbre ................................................................... 50
Figure 4.9: Vus éclatée et coupe du montage moyeu-pales ................................................................. 51
Figure 4.10: Montage par clavette ........................................................................................................ 55
Figure 4.11: Détermination de s’........................................................................................................... 56
Figure 4.12: Caractéristiques géométriques d’une cannelure .............................................................. 56
Figure 4.13: Schématisation du chargement de la tour ........................................................................ 57
Figure 4.14: Schématisation du chargement de la tour haubanée ....................................................... 59
Figure 4.15: Développement d'un accouplement Rotex sur DAO 3D ................................................... 60
Figure 4.16: Diamètre en fonction de la puissance conventionnelle.................................................... 64
Figure 4.18 : Entrefer et charge linéaire en fonction du pas polaire .................................................... 64
Figure 4.19: Multiplicateur à 2 étages................................................................................................... 66
Figure 4.20: Représentation d'une dent ............................................................................................... 66
Figure 4.21: Etrier de frein .................................................................................................................... 66
Figure 4.22: Zone de contacte garniture/frein...................................................................................... 66
Figure 4.23: Distribution des contraintes (tour simple) ........................................................................ 68
Figure 4.24: Distribution des contraintes (tour haubanée) .................................................................. 68
Figure 4.25: Distribution des déplacements (tour simple).................................................................... 69
Figure 4.26: Distribution des déplacements (tour haubanée) .............................................................. 69
Figure 4.27: Distribution des contraintes.............................................................................................. 70
Figure 4.28: Distribution des déplacements ......................................................................................... 71
Figure 4.29: Distribution des contraintes.............................................................................................. 71
Figure 4.30: Distribution des déplacements ......................................................................................... 72
Figure 4.31: Chargement du premier élément de pale......................................................................... 72
Figure 4.32 : Distribution de contrainte sur le premier élément de pale ............................................. 73
Figure 4.33: Distribution des déplacements sur le premier élément de pale....................................... 74
Figure 4.34 : Chargement sur le dernier élément de pale .................................................................... 74
Figure 4.35 : Distribution de contrainte sur le dernier élément de pale .............................................. 74
Figure 4.36 : Distribution des déplacements sur le dernier élément de pale....................................... 74
 Liste des tableaux
Tableau 2.1: Coefficients spécifiques de vitesse suivant le nombre de pales ...................................... 25
Tableau 2.2: Géométrie de la pale calculé par le code………………………………………………………………………………………30

Tableau 2.3: Géométrie de la pale pour le cas idéal (Cl=1, Cd=0).......................................................... 31

Tableau 4.1 : Valeurs du coefficient de température St ........................................................................ 62
Tableau 4.2: Démarrage Sz .................................................................................................................... 62
Tableau 4.3: A-coup SA/SL ...................................................................................................................... 62
Tableau 4.4 : Diamètre externe en fonction du nombre de pole……………………………………………………….63

Tableau 4.5 : Puissance nominale……………………………………………………………………………………………………..63

Nomenclature
Lettres latines

A : Section [m²]

a : Facteur d’induction axial

a’ : Facteur d’induction angulaire

B : Nombre de pales

c : Corde d’un profil [m]

Cd : Coefficient de trainé

Cl : Coefficient de portance

Cp : Coefficient de puissance

CS : Coefficient de sécurité

D : Diamètre [m]

dFd : Force élémentaire de trainée [N/m]

dFl : Force élémentaire de portance [N/m]

dFN : Force élémentaire normale [N/m]

dFT : Force élémentaire tangentielle [N/m]

F : Facteur de correction de prandlt

Fc : Force centrifuge [N]

Fp : Force de poussé du rotor (force axiale) [N]

I : Moment quadratique

M : Moment [Nm]

݉̇ : Débit massique [kg/s]

N : Nombre de subdivision

U : Vitesse [m/s]

P : Puissance [W]

P : Poids [N]

Q : Couple [Nm]
R : Rayon du rotor [m]

T : Poussée [N]

Lettres grecs

α : Angle d’attaque [°]

β : Angle de calage [°]

Ф : Angle d’écoulement relatif [°]

λ : Vitesse spécifique

σ : Contrainte [N/mm²]

µ : Coefficient de frottement

ηm : Rendement mécanique

τ : Contrainte de cisaillement [N/mm²]

Indices

a : Air

adm : Admissible

c : Centrifuge

D : Démarrage

d : Trainée

f : Flexion

l : Portance

M : Maximale

n : Nominale

p : Poussée

r : Rotor

rel : Relatif

s : Solide

t : Tige
 INTRODUCTION

INTRODUCTION

La consommation mondiale en électricité s’accroit au fil des années. L’essentiel de

cette électricité est produit dans des centrales à énergie fossile. Ces centrales consomment
énormément de combustible et rejettent de très importantes quantités de CO2. L’énergie
nucléaire de part le haut degré de maîtrise technologique qu’elle nécessite et la gestion des
risques qui l’entoure reste problématique. Par ailleurs l’énergie hydraulique représente une
forme d’énergie renouvelable non négligeable. Mais elle reste limitée aux pays dotés de cours
d’eau importants et nécessitent de gros investissements à travers la construction de barrages.
De nos jours et par souci de préserver l’environnement, on se tourne vers la production
d’énergies non polluantes et renouvelables telles que les énergies solaires ou éoliennes. C’est
dans le cadre de la sensibilisation et la contribution à la protection de l’environnement, que
nous avons pensé à élaborer un document qui peut servir de base et de guide pour la
conception d’une éolienne de faible puissance.

Dans un premier temps, nous exposerons un aperçu historique ainsi que l’état de l’art
des différentes technologies, leurs évolutions, avantages et inconvénients. L’état de l’art
consiste à présenter les différents types d’éoliennes et leurs architectures ; préciser les
caractéristiques de chacun des types pour enfin connaître les avantages et inconvénients de
chacune d’elles. Ceci nous permettra d’avoir une première idée sur la technologie qui sera
choisie pour la conception de notre éolienne.

Puis un deuxième temps, nous ferons un rappel des principales théories dont nous
aurons besoin lors de la conception, spécialement celle du dimensionnement du rotor .Nous
commencerons par la théorie de base développée par Betz au début des années 1920. Cette
théorie nous permet d’avoir une approximation de l’énergie maximale récupérable de
l’énergie cinétique du vent. Plusieurs approches ont été formulées après la théorie de Betz
dont la théorie de l’élément de pale qui consiste à subdiviser la pale en un grand nombre
d’éléments pour ensuite, en utilisant les caractéristiques aérodynamiques de la géométrie du
profil, calculer les forces agissant sur ces dernières. Une autre théorie largement utilisée de
nos jours combine la théorie de quantité de mouvement et celle de l’élément de pale appelée
BEM (Blade Element Momentum theory) ou théorie de quantité de mouvement-élément de
pale. Après cela, nous détaillerons la procédure de design du rotor basée sur la méthode BEM.
Un organigramme de calcul de la géométrie et un autre pour la détermination de la courbe de

1
 INTRODUCTION

puissance seront présentés, suivi des résultats de calcul pour une éolienne de puissance
nominale égale à 25 kW pour une vitesse du vent de 15m/s.

Enfin, nous donnerons des indications pour le choix et le dimensionnement des

organes mécaniques qui consiste en premier abord en un calcul analytique approché avec
certaines approches simplificatrices. L’étape suivante est la vérification des dimensions par le
logiciel SolidWorks qui sera suivi par d’illustrations de montage des différents éléments de
l’éolienne.

2
 CHAPITRE 1 : GENERALITES

CHAPITRE 1 : GENERALITES
Dans ce chapitre, nous présenterons des informations qui seront nécessaires lors du
dimensionnement et la conception. Nous commencerons par rappeler l’historique de la
technologie éolienne ainsi que son évolution. Nous passerons en suite au classement des
éoliennes selon leurs configurations (position du rotor) et le mode de fonctionnement. Une
description des caractéristiques nous aidera à formuler notre choix quant à la configuration
qui sera adopté pour la conception. Nous finirons par présenter les théories liées à
l’aérodynamique du rotor sur lesquelles s’appuiera notre démarche pour le calcul de la
géométrie des pales.

1.1 Historique

L’origine des éoliennes est sujet de débat entre les historiens. Certains considèrent que
l’éolienne remonte au 12ème siècle avec l’apparition des moulins à vents « windmills » en
Europe, qui ressemblent fortement aux éoliennes les plus répandus de nos jours. D’autres
historiens affirment que la technologie éolienne remonte bien plus loin dans le passé à
l’époque du règne de la civilisation Perse. Les Chinois, les Perses et les Égyptiens utilisaient
déjà la force du vent il y a des dizaines de siècles. Plus tard les peuples de la région de la mer
Égée utilisèrent le même principe.

L’origine la plus probable, du fait que « Eole » est un nom perse qui signifie dieux des
vents, remonte à la civilisation Perse.

Malgré l’ancienneté de cette technologie, son utilisation pour la production d’énergie

électrique était faible, jusqu’en 1974, l’année de la crise pétrolière, qui relança les études,
recherches, et réalisations. L’expérience californienne a été la première à grande échelle (le
Wind-rush) au début des années 80, notamment avec des turbines de moyenne puissance (55
kW) et grâce à une incitation fiscale très volontariste. On passa ainsi de 144 machines (pour
un total de 7MW) en 1981, à 4687 machines (386 MW) en 1985[1]. Mais c’est vers la fin des
années 1980 que le marché des systèmes raccordés au réseau a réellement décollé en Europe,
dans le reste des USA et également en Asie et en Afrique du Nord.

Depuis ce temps là des études et des recherchent ont été faites permettant ainsi une
évolution rapide de la technologie. On peut classer les éoliennes en deux principaux types :
les éoliennes à axe horizontal et les éoliennes à axe vertical.

4
 CHAPITRE 1 : GENERALITES

1.2 Types et caractéristiques des éoliennes

1.2.1 Eolienne à axe vertical

Les éoliennes à axe vertical sont encore au stade d’essai de prototypes, hormis
quelques réalisations de petite puissance industrialisées par la société Winside [2].

[Link] Rotor Savonius

C’est le même principe que l’anémomètre (la traînée différentielle). Les efforts exercés
par le vent sur chacune des faces d’un corps creux sont d’intensités différentes (figure1.1). Il
en résulte donc un couple moteur, que l’on peut utiliser pour entraîner un générateur
électrique ou un autre dispositif mécanique tel qu’une pompe.

Le rotor Savonius, du nom de son inventeur, est amélioré comparé à l’anémomètre, et ce,
en permettant la circulation de l’air entre les deux cylindres ce qui augmente le couple
moteur.

Figure 1.1: Effet du vent sur un corps creux [2]

5
 CHAPITRE 1 : GENERALITES

[Link] Rotor Darrieus

Tout profil placé dans un écoulement d’air selon différents angles est soumis à des
forces d’intensité et de directions qui changent avec le changement de l’écoulement (principe
de l’incidence variable). La combinaison de ces forces génère alors un couple moteur
Figure1.2. Darrieus, aux années 1930 a breveté son invention qui repose sur ce principe. De
tels rotors peuvent être de forme cylindrique, tronconique, parabolique…

ሬ
ࢂሬሬሬ⃗ ᇱ
૙: ࢜࢏࢚ࢋ࢙࢙ࢋ ࢇ࢞࢏ࢇ࢒ࢋ ࢊ࢛ ࢜ࢋ࢔࢚ࢋ࢔ ࢇ࢓ ࢕࢔࢚ࢊࢋ ࢒ࢋ࢕࢒࢏ࢋ࢔࢔ࢋ
ࢹ: ࢜࢏࢚ࢋ࢙࢙ࢋ ࢇ࢔ࢍ࢛࢒ࢇ࢏࢘ࢋ ࢊࢋ ࢘࢕࢚ࢇ࢚࢏࢕࢔ ࢊࢋ࢙࢖ࢇ࢒ࢋ࢙
R : rayon de la pale
ሬሬ⃗: ࢜࢏࢚ࢋ࢙࢙ࢋ ࢛࢔࢏࢚ࢇ࢏࢘ࢋ ࢚ࢇ࢔ࢍࢋ࢔࢚ࢋ ࢇ ࢒ࢇ ࢖ࢇ࢒ࢋ
ࢁ

Figure 1.2: Principe de l’incidence variable [2]

Figure 1.3 : Technologies d’éoliennes à axe vertical [1]

6
 CHAPITRE 1 : GENERALITES

[Link] Caractéristiques

Les machines de types Savonius simples ne peuvent démarrer spontanément. Il est

clair qu’au démarrage (phase d’établissement du vent), les cylindres sont orientés par rapport
au vent de sorte qu’il en résulte un couple nul (figure 1.1). Pour remédier à ce défaut on peut
superposer deux machine de type Savonius mais avec un décalage angulaire ce qui permettra
le démarrage de l’éolienne quelque soit la direction du vent ou bien prévoir un système de
démarrage (batterie de démarrage).

Les machine de types Darrieus possèdent le même défaut que les machines de types
Savonius .On peut palier à ce défaut en plaçant une éolienne type Savonius sur le rotor de la
machine Darrieus ou utiliser la génératrice en mode moteur en utilisant des batteries de
démarrages.

Un autre désavantage des éoliennes à axe vertical est l’édification des machines près
du sol ce qui signifie que le capteur d’énergie se situe dans une zone peu favorable (gradient
de vent, turbulence due aux accidents du terrain en amont de la machine), ce qui réduit
significativement l’efficacité de la machine. Par ailleurs, le principe même de fonctionnement,
basé sur des variations incessantes de charges aérodynamiques sur les pales, fait que ces
éoliennes sont très sujettes aux problèmes d’aéroélasticité [2]. Enfin, pour des éoliennes de
grande puissance, la surface occupée au sol par le haubanage est très conséquente.

1.2.2 Eolienne à axe horizontal

Le concept de l’éolienne à axe horizontale (HAWT Horizontal Axis Wind Turbine)

consiste en une roue à aubes (hélice) supporté par une tour d’une hauteur équivalente au triple
du rayon, et ceux pour capter le maximum de vent en évitant les obstacles tels que les
bâtiments ou les arbres. Les premières applications de cette technologie remontent au moyen
âge, apparitions des moulins à vents en Europe.

[Link] Eolienne face au vent « upwind »

L’écoulement de vent attaque le rotor par la face. Le principal avantage de cette

configuration est qu’elle permet d’éviter l’effet d’ombrage (nuance) derrière la tour qui est
très important comparé au même effet avant la tour. En effet l’écoulement commence à se
plier loin avant d’attaquer la tour, ce qui provoque un léger effet de nuance. Notons qu’à

7
 CHAPITRE 1 : GENERALITES

chaque fois que l’une des pales du rotor coïncide avec la tour la puissance subit une
perturbation qui entraine une diminution du rendement de l’éolienne.

[Link] Eolienne sous le vent « downwind »

L’écoulement du vent est d’incidence derrière le rotor. Le principale avantage de cette

configuration est la possibilité de suivre la direction du vent sans nécessitée de mécanisme
annexes à condition que le rotor et la nacelle aient été conçu de sorte à permettre cette
opération. Cet avantage permet une construction plus flexible des pales du rotor du fait de la
mobilité de la nacelle. La flexibilité du rotor permet, aux grandes vitesses, d’absorber une
partie des charges sur la nacelle et ce en se pliant tout en restant dans la limite d’élasticité du
matériau. La rotation de la nacelle pour suivre le vent n’est pas sans danger. Le risque se situe
au niveau du câble électrique qui, avec la rotation de la nacelle, subit une torsion. Au delà
d’un certain nombre de tours le câble sera cisaillé ce qui est très dangereux d’où nécessitée
d’un compteur de tours et d’un mécanisme qui permet le redressement de la nacelle.
L’inconvénient majeur est le passage du rotor dans une zone d’ombrage qui provoque des
fluctuations de la puissance et une fatigue des organes mécaniques.

1.3 Régulation

Les objectifs de la régulation sont d’assurer la sécurité de l’éolienne par vents forts et
de limiter la puissance. Une turbine éolienne est dimensionnée pour développer sur son arbre
une puissance dénommée puissance nominale Pn. La puissance Pn est obtenue à partir d'une
vitesse du vent Un , dénommée vitesse nominale. Lorsque la vitesse du vent est supérieure à
Un la turbine éolienne doit modifier ses paramètres afin d'éviter la destruction mécanique, de
sorte que sa vitesse de rotation reste pratiquement constante.
On définit aussi :
 la vitesse de démarrage UD , à partir de laquelle l'éolienne commence à fournir de
l'énergie,
 la vitesse maximale du vent, UM , pour laquelle la turbine ne convertit plus l'énergie
éolienne, pour des raisons de sûreté de fonctionnement.

Les vitesses UD , Un et UM définissent quatre zones sur le diagramme de la puissance utile en

fonction de la vitesse du vent :

8
 CHAPITRE 1 : GENERALITES

UD Un UM
Figure 1.4: Diagramme de la puissance utile sur l'arbre en fonction de la vitesse du vent [1]

 la zone I, où P = 0 (la turbine ne fonctionne pas),

 la zone II, dans laquelle la puissance fournie sur l'arbre dépend de la vitesse du vent U,
 la zone III, où la vitesse de rotation est maintenue constante et où la puissance P
fournie reste égale à Pn,
 la zone IV, dans laquelle le système de sûreté du fonctionnement arrête le transfert de
l'énergie.

La plupart des grandes turbines éoliennes utilisent deux principes de contrôle aérodynamique
pour limiter la puissance extraite à la valeur de la puissance nominale de la génératrice :

 Système à pas ou calage variable « Variable Pitch Contrôle » qui permet d’ajuster la
portance des pales à la vitesse du vent pour maintenir une puissance sensiblement
constante dans la zone III de vitesses
 Système décrochage aérodynamique « Stall Controll », le plus robuste car c’est la
forme des pales qui conduit à une perte de portance au-delà d’une certaine vitesse de
vent, mais la courbe de puissance chute plus vite.
 Système « Stall Actif » dans lequel le décrochage aérodynamique est obtenu
progressivement grâce à une orientation minime des pales nécessitant des moyens de
réglage plus économiques et plus robustes que dans le système à calage variable.

9
 CHAPITRE 1 : GENERALITES

La figure 1.5 montre les écoulements d’air dans une situation « normale » et dans une
situation de « décrochage » :

Figure 1.5: Ecoulement d’air sur une pale en régime « normal » (à gauche) et en situation de
décrochage (à droite) [1]

1.3.1 Eolienne à contrôle passif (Stall control )

Le contrôle à décrochage de charge est un système de contrôle passif qui réagit au-delà
d’une certaine vitesse du vent. Les pales du rotor sont fixes et elles ne peuvent pas tourner
autour de leur axe longitudinal. L’angle de calage est choisi afin de permettre le décrochage
du flux d’air de la surface de la pale pour des vitesses du vent plus grandes que la vitesse
nominale. Cela réduit la portance et augmente la traînée.

En comparaison avec les turbines à angle de calage variable, celles à décrochage ont les
avantages suivants :
 pas de système de contrôle d’angle de calage ;
 construction plus simple et moins coûteuse du rotor ;
 maintenance plus aisée ou meilleure fiabilité (nombre plus réduit de pièces
mécaniques en mouvement).

On peut combiner ce dispositif à une machine à deux vitesses de rotation, ce qui permet
d’étendre la plage de production. La figure 1.6 compare l’évolution du coefficient de
puissance, CP, des systèmes pitch (à vitesse de rotation fixe) et stall à deux vitesses.

10
 CHAPITRE 1 : GENERALITES

Figure 1.6: Comparaison des systèmes pitch à une vitesse de rotation et stall à deux vitesses
Fixe [1]

1.3.2 Eolienne à contrôle actif « pitch control »

L’orientation des pales est fonction de la vitesse du vent, les pales sont face au vent
aux basses vitesses (de la vitesse minimale à la vitesse de base) puis s’inclinent pour atteindre
la position « drapeau » à la vitesse maximale.

Mécaniquement, il s’agit d’un dispositif sophistiqué (commande par vérins

hydrauliques dans les grosses machines), surtout pour les petites machines où il grèverait le
coût et pour les très grosses puissances (au delà du 1 MW) où il devient fragile.

La modification de l'angle de calage de la pale de l'éolienne à quatre objectifs distincts :

 démarrage à une vitesse du vent, UD, réduite,

 régulation de la vitesse de rotation Ω, pour U > Un,
 optimisation du régime de conversion de l'énergie, quand la vitesse du vent évolue
entre les limites [UD, Un],
 protection de l'éolienne contre les grandes vitesses du vent (U > UMax), par la mise en
"Drapeau" des pales de l'hélice

11
 CHAPITRE 1 : GENERALITES

Figure 1.7: Effet de l’angle de pas sur le coefficient aérodynamique de puissance éolienne [1]

Ce système de régulation offre les avantages suivants :

 la possibilité de pouvoir effectuer un contrôle actif de la puissance sous toutes les
conditions de vent (en dehors de la vitesse limite de sécurité), y compris à une
puissance partielle ;
 elles peuvent fournir la puissance nominale même quand la densité de l’air est faible
(des sites sur les hauteurs, grandes températures) ;
 une production d’énergie plus importante que les éoliennes à contrôle à décrochage
dans les mêmes conditions de fonctionnement ;
 démarrage rapide par changement de l’angle de calage ;
 il n’y a pas besoin de freins puissants pour les arrêts d’urgence ;
 les contraintes des pales sont plus faibles quand les puissances sont supérieures à la
valeur nominale ;
 la masse réduite des pales du rotor mène à une masse réduite de la turbine éolienne.

12
 CHAPITRE 1 : GENERALITES

1.4 Aérodynamique des éoliennes

1.4.1 Théorie de Betz

L’une des théories élémentaires utilisées dans la conception des rotors d’éoliennes est
la théorie attribuée à Betz (1929) [1] (modèle du disque actuateur). Les développements de
Betz peuvent être utilisés pour déterminer, et ce pour une éolienne idéale, la poussée
développée par le vent et l’influence du rotor sur l’écoulement du vent.

Modèle du disque actuateur : en mécanique des fluides, le disque actif est défini comme une
surface de discontinuité où des forces de surface agissent sur l’écoulement [3]. Ce modèle est
extrêmement simplifie et repose sur les hypothèses suivantes

[Link] Hypothèses
 Fluide homogène et incompressible ;
 La géométrie du rotor est effacée et ce dernier n’est représenté que par un
disque d’épaisseur nulle de diamètre D du rotor;
 Ecoulement stationnaire, irrotationnel ;
 Poussé uniforme dans tout le disque rotor ;
 La pression statique loin en amont et en aval du rotor est égale à la pression
statique de l’ambiance ;
 Non rotation du disque ;

U1
U4
P+
U3

U2
P-

Figure 1.8: Modélisation du rotor éolien par un disque actuateur [3]

13
 CHAPITRE 1 : GENERALITES

[Link] Facteur d’induction axiale

En appliquant la conservation de quantité de mouvement au volume de contrôle (fig 1.8) on

peut écrire :

ܶ = ܷଵ(ߩ௔ ‫ܷܣ‬ଵ) − ܷସ(ߩ௔ ‫ܷܣ‬ସ) (1.1)

Avec : A section du disque rotor

ߩ௔ : masse volumique de l’air

U : vitesse de l’air

ௗ௠̇
Sachant que l’écoulement est stationnaire = 0 on aura :
ௗ௧

(ߩ௔ ‫ܷܣ‬ଵ) = (ߩ௔ ‫ܷܣ‬ସ) = ݉̇ (1.2)

Il en résulte ܶ = ݉̇ (ܷଵ − ܷସ) (1.3)
La poussée étant positive, la vitesse décroit de l’amont vers l’avale U1 > U4

A présent essayons de déterminer une relation entre les vitesses des différentes zones.

En amont et en aval du disque, nous avons aucune production de travail, ce qui nous permet
d’appliquer l’équation de Bernoulli dans les deux volumes contrôle amont et aval.

1 1 (1.4)
‫݌‬ଵ + ߩ௔ ܷଵଶ = ‫݌‬ଶ + ߩ௔ ܷଶଶ
2 2
1 1
‫݌‬ଷ + ߩ௔ ܷଷଶ = ‫݌‬ସ + ߩ௔ ܷସଶ (1.5)
2 2
La poussée peut être exprimée en fonction des forces agissant de part et d’autre du disque
rotor :

T = Aଶ(‫݌‬ଶ − ‫݌‬ଷ) (1.6)

D’après les hypothèses on a p1=p4, et en admettant que la vitesse axiale au niveau du rotor
reste constante (U2=U3), on obtient en combinant les équations (1.4), (1.5) et (1.6) :

1 (1.7)
ܶ= ߩ ‫ ܷ( ܣ‬ଶ − ܷସଶ)
2 ௔ ଶ ଵ

14
 CHAPITRE 1 : GENERALITES

En égalisant les équations (1.7) et (1.3) avec ݉̇ = ߩ‫ܣ‬ଶܷ :

ܷଵ + ܷସ (1.8)
ܷଶ =
2
Finalement on trouve comme résultat que la vitesse au niveau du rotor est la moyenne des
vitesses amont et aval de l’écoulement.

On définit le coefficient d’induction axiale, a , pour avoir les expressions des vitesses au
niveau du rotor et en aval en fonction de la vitesse amont U1.

ܷଵ − ܷଶ (1.9)
ܽ=
ܷଵ
ܷଶ = ܷଵ(1 − ܽ) (1.10)
ܷସ = ܷଵ(1 − 2ܽ) (1.11)
Remarque :

 Des équations (1.10) et (1.11) on déduit que la vitesse induite au niveau du rotor est
égale a la moitie de la vitesse induite à l’infini aval.
 On peut aussi remarquer que la vitesse en aval du rotor peut être égale à zéro si le
coefficient a=1/2. La théorie de Betz est limitée par le coefficient d’induction axiale
qui ne peut dépasser la valeur limite de 0.5.

[Link] Coefficients de puissance et de poussée

On définit le coefficient de puissance comme suit :

‫ݎ݋ݐ݋ݎ ݑ݀݁ܿ݊ܽݏݏ݅ݑ݌‬ ܲ (1.12)

‫ܥ‬௣ = =
‫ ݐ݊݁ݒ݈݁ݏ݈ܾ݊ܽ݀݁݋݊݋݌ݏ݅݀݁ܿ݊ܽݏݏ݅ݑ݌‬1 ߩ‫ ܷܣ‬ଷ
2
La puissance du rotor est égale a la poussée que multiplie la vitesse axiale U2.

1 1 (1.13)
ܲ = ߩ௔ ‫ܣ‬ଶ(ܷଵଶ − ܷସଶ)ܷଶ = ߩ௔ ‫ܣ‬ଶܷଶ(ܷଵ − ܷସ)(ܷଵ + ܷସ)
2 2
En substituant U4 et U2 par les équations (10) et (11) la puissance du rotor s’exprime par :

1 (1.14)
ܲ= ߩ ‫ ܷܣ‬ଷ 4ܽ (1 − ܽ)²
2 ௔

15
 CHAPITRE 1 : GENERALITES

Alors le coefficient de puissance est :

‫ܥ‬௣ = 4ܽ (1 − ܽ)² (1.15)

On définit le coefficient de poussée par la relation suivante :

ܶ ‫ݏݏݑ݋݌‬é݁ (1.16)
‫= ்ܥ‬ =
1 ଶ ݂‫ݐ݊݁ݒ ݑ݀݁ݑݍ݅ ݉ܽ݊ݕ݀݁ܿݎ݋‬
2 ߩ௔ ‫ܷܣ‬
1 (1.17)
ܶ = ߩ௔ ‫ܷܣ‬ଵଶ4ܽ(1 − ܽ)
2
‫ = ்ܥ‬4ܽ(ܽ − 1) (1.18)
[Link] Limite de Betz

Pour déterminer le coefficient de puissance maximale, il suffit de dériver l’expression

(1.15) et de l’égaliser à zéro, ce qui nous donne :

16 (1.19)
‫ܥ‬௣௠ ௔௫ = ≅ 0,5962
27

Si on arrive à concevoir un rotor permettant d’avoir une vitesse axial égale au 2/3 de la
vitesse amont du vent, alors le rotor nous fournira le maximum de puissance; et d’après les
lois basiques de la physique c’est la puissance maximale possible.

Le coefficient de poussée a pour expression données par eq. (1.18). CT maximum a

pour valeur 1, qui correspond à un facteur d’induction axiale a=0 ,5 et à une vitesse aval égale
à zéro. Pour un maximum de puissance, le coefficient de poussée a la valeur de 8/9. La figure
1.9 illustre les variations des coefficients de puissance et de poussée pour un rotor idéal en
fonction du coefficient d’induction axiale a.

Figure 1.9: Paramètres opératoires d’une turbine Betz [4]

16
 CHAPITRE 1 : GENERALITES

1.4.2 Méthode BEM

La méthode la plus utilisée pour calculer la charge appliquée par l’écoulement sur une
éolienne et le comportement d’une éolienne, consiste à utiliser le bilan de quantité de
mouvement appliqué à un élément de pale. Cette méthode est basée sur la division de
l’écoulement en volumes de contrôle annulaires d’épaisseur dr, auxquels on applique le bilan
de quantité de mouvement et d’énergie. Ces anneaux s’étendent de l’infini amont jusqu’à
l’infini aval par rapport au rotor. Les hypothèses principales de cette méthode sont :
 que la vitesse induite dans le plan de rotation est égale à la moitié de la vitesse induite
à l’infini aval,
 que l’on peut analyser l’écoulement par la division de la pale en nombre d’éléments
indépendants et que la force d’un élément de pale est seule responsable de la variation
de quantité de mouvement de l’air qui passe dans l’anneau balayé par cet élément.

[Link] Théorie de l’élément de pale

Les forces agissant sur la pale peuvent être exprimées en fonction des coefficients de
portance et de trainé ainsi que l’angle d’attaque qui caractérisent un profil bidimensionnel.
Pour l’analyse qui suit nous diviserons la pale en N éléments comme montré sur la figure
1.10.

Figure 1.10:Subdivision de la pale [5]

17
 CHAPITRE 1 : GENERALITES

La figure ci-dessous schématise un élément de pale avec toutes les forces agissant sur
lui et les différents paramètres définissant la géométrie du profil.

dFL

dFD

Urel

Figure 1.11: profil d’une pale d’éolienne à axe horizontale .définition des variables [4]

D’après la figure on peu tirer les relations suivantes :

ܸ(1 − ܽ) 1−ܽ (1.20)

tan(߮) = =
ߗ‫(ݎ‬1 + ܽ' ) (1 + ܽ' )ߣ௥
Avec a’ : facteur d’induction angulaire. ([3] §§ 3.3 page 89)

ܷ(1 − ܽ) (1.21)
ܷ௥௘௟ =
sin (߮)
1 ଶ (1.22)
݀‫ܨ‬௟ = ‫ܥ‬௟ߩ௔ ܷ௥௘௟ ܿ݀‫ݎ‬
2
1 ଶ
݀‫ܨ‬ௗ = ‫ܥ‬ௗ ߩ௔ ܷ௥௘௟ ܿ݀‫ݎ‬ (1.23)
2
݀‫ܨ‬ே = ݀‫ܨ‬௟cos(߮) + ݀‫ܨ‬ௗ sin (߮) (1.24)
݀‫ܨ݀ = ்ܨ‬௅ sin(߮) − ݀‫ܨ‬஽ sos (߮) (1.25)

18
 CHAPITRE 1 : GENERALITES

Urel : vitesse d’écoulement relative ;

Φ : angle d’ecoulement relative ;

dFL : force élémentaire de portance ;

dFD : force élémentaire de traînée ;

Cl : coefficient de portance ;

Cd : coefficient de traînée ;

c : corde ;

dFN : force élémentaire normale ;

dFT : force élémentaire tangentielle;

La force normale totale agissant sur tout le disque est :

1 (1.26)
݀‫ܨ‬ே = ‫ ܷ ߩܤ‬ଶ [‫ ܥ‬cos(߮) + ‫ܥ‬ௗ sin(߮)]ܿ݀‫ݎ‬
2 ௔ ௥௘௟ ௟
Avec B : nombre de pales.

Le couple élémentaire agissant sur un élément dr est donnée par :

݀ܳ = ‫்ܨ݀ݎܤ‬ (1.27)
1 ଶ [‫ܥ‬ (1.28)
݀ܳ = ‫ߩܤ‬௔ ܷ௥௘௟ ௅ sin(߮) − ‫ܥ‬ௗ sin(߮)]ܿ‫ݎ݀ݎ‬
2
La traînée influence négativement le couple donc la puissance. Plus la traînée est
élevée plus la puissance chute. Donc pour avoir le maximum de puissance, il nous faut
optimiser le profil pour obtenir le minimum de traînée.

La théorie de l’élément de pale nous permet de définir les forces normales et

tangentielles agissant sur le disque rotor en fonction des caractéristiques bidimensionnelles du
profil et des angles de l’écoulement.

19
 CHAPITRE 1 : GENERALITES

[Link] Bilan de quantité de mouvement appliqué à l’élément de pale

Le bilan de quantité de mouvement sur un élément de pale nous permet de déterminer

la poussée et le couple élémentaire. Un développement similaire à celui de la théorie de
Betz en considérant un élément annulaire nous permet d’avoir comme résultat ce qui suit :

݀ܶ = 4ߩ௔ ܷ ଶܽ(1 − ܽ)ߨ݀‫ݎ‬ (1.29)

݀ܳ = 4ܽ' (1 − ܽ)ߩ௔ ܷߨ‫ݎ‬ଷߗ݀‫ݎ‬ (1.30)

[Link] Couplage théorie de quantité de mouvement et élément de pale

஻௖
On définit la solidité ߪ' =
ଶగ௥

Les équations (1.24) et (1.26) s’écrivent alors :

ܷ ଶ(1 − ܽ)ଶ (1.31)

݀‫ܨ‬ே = ߪ ߨߩ௔
' [‫ܥ‬௟cos(߮) + ‫ܥ‬ௗ sin(߮)]‫ݎ݀ݎ‬
‫݊݅ݏ‬²(߮)
ܷ ଶ(1 − ܽ)ଶ (1.32)
݀ܳ = ߪ' ߨߩ௔ [‫ܥ‬௟sin (φ) − ‫ܥ‬ௗ cos(߮)]‫ݎ‬²݀‫ݎ‬
‫݊݅ݏ‬²(߮)
La méthode BEM consiste en le calcul des facteurs d’induction axiale et angulaire afin
de pouvoir déterminer, dans un premier temps, la géométrie de la pale (variation de la corde
et angle de vrillage) pour ensuite prédire les performances de cette pale.

La méthode BEM accepte comme hypothèse simplificatrice que le coefficient de

traînée est négligeable et peut être pris égale à zéro sans entraîner des erreurs importantes [4]
donc en égalisant les équations (1.28) et (1.30) on obtient :

ܽ' ߪ'‫ܥ‬௅ (1.33)

=
(1 − ܽ) 4ߣ௥sin (߮)
En égalisant les équations (1.28) et (1.31) on obtient :

ܽ ߪ'‫ܥ‬௟cos (߮) (1.34)

=
(1 − ܽ) 4‫݊݅ݏ‬²(߮)
Après quelques manipulations algébriques et en tenant compte des considérations
géométriques, nous obtenons les relations suivantes :

ܿ‫ ߮ݏ݋‬− ߣ௥‫߮݊݅ݏ‬ (1.35)

‫ܥ‬௟ = 4sin (߮)
ߪ' (‫߮݊݅ݏ‬
+ ߣ௥ܿ‫)߮ݏ݋‬

20
 CHAPITRE 1 : GENERALITES

ܽ' ߪ'‫ܥ‬௟ (1.36)

=
1 + ܽ' 4ܿ‫߮ݏ݋‬
ܽ ߣ௥ (1.37)
=
ܽ' ‫߮݊ܽݐ‬
1
ܽ= (1.38)
4‫݊݅ݏ‬ଶ߮
൤1 + ' ൨
ߪ ‫ܥ‬௟ܿ‫߮ݏ݋‬
1 (1.39)
ܽ' =
4ܿ‫߮ݏ݋‬
൤ ' − 1൨
ߪ ‫ܥ‬௟

Pour déterminer les forces agissant sur le disque rotor ainsi que les conditions
d’écoulement il faut résoudre les équations ci-dessus et ce en utilisant une méthode itérative
sur les facteurs a et a’ ou bien sur le coefficient de portance Cl et l’angle d’attaque α.

1.4.3 Pertes en bous de pales

Dans un écoulement autour d’un profil d’aube l’air est accéléré sur l’extrados et
décéléré sur l’intrados. Ceci se traduit par une différence de pression entre les deux parties.
Les particules de fluides parcourant l’intrados ont tendance à s’enrouler en bout de pale du
fait de la dépression. Il en résulte des pertes de performances qui se répercutent sur la
puissance.

Un certain nombre de méthodes ont été suggérées pour inclure l'influence des pertes
en bout de pales sur la puissance. L'approche la plus utilisée est celle développée par Prandtl,
qui consiste à introduire un facteur F, dans les développements de la méthode BEM. Le
coefficient de correction de Prandtl dépend de l’angle relatif d’écoulement ߮, du nombre de
pale B et du rayon de la station. Son expression est donnée par

‫ܤ‬ ‫ݎ‬
2 ቀ1 − ܴቁ+
2
‫ = ܨ‬acos exp ቎− ቐ ‫ݎ‬ ቑ቏
ߨ (1.40)
sin ߮
ܴ
Le résultat de l’inverse de la fonction cosinus doit être considéré en radian et notons que le
coefficient est toujours inferieur à 1.

Les différentes équations affectées par le facteur de Prandtl sont les suivantes :

݀ܶ = ‫ߩܨ‬௔ ܷ ଶ4ܽ(1 − ܽ)ߨ݀‫ݎ‬ (1.41)

21
 CHAPITRE 1 : GENERALITES

݀ܳ = 4‫ܽ ܨ‬ᇱ(1 − ܽ)ߩ௔ ܷߨ‫ݎ‬ଷߗ݀‫ݎ‬ (1.42)

cos ߮ − ߣ௥ sin ߮ (1.43)
‫ܥ‬௟ = 4‫߮ ݊݅ݏ ܨ‬
ߪ (sin ߮ + ߣ௥ cos ߮)
1 (1.44)
ܽ=
4‫݊݅ݏܨ‬ଶ߮
൤1 + ' ൨
ߪ ‫ܥ‬௟ܿ‫߮ݏ݋‬
1
ܽ' =
4‫߮ݏ݋ܿܨ‬ (1.45)
൤ ' − 1൨
ߪ ‫ܥ‬௟

1.4.4 Model de décrochage de Laino & Hansen

Le modèle consiste à utiliser les coefficients de portance et trainée qui sont développés en
écoulement sur trois dimensions contrairement à notre première approche qui utilise les
données en deux dimensions. Ce modèle est déduit par les relations suivantes :

‫ܥ‬௟,ଷௗ = ‫ܥ‬௟,ଶௗ + ݂ఈ ∆‫ܥ‬௟ (1.46)

‫ܥ‬ௗ,ଷௗ = ‫ܥ‬ௗ,ଶௗ + ݂ఈ ∆‫ܥ‬ௗ (1.47)

∆‫ܥ‬௟ = ݂௟൫‫ܥ‬௟,௣ − ‫ܥ‬௟,ଶௗ ൯ (1.48)

‫ ߙ ݊݅ݏ‬− ݇ܿ‫ߙݏ݋‬
∆‫ܥ‬ௗ = ൬ ൰∆‫ܥ‬௟
݇‫ ߙ݊݅ݏ‬+ ܿ‫ߙݏ݋‬ (1.49)
‫ܥ‬௟,௣ = 2ߨ(ߙ − ߙ଴) (1.50)
ௗோ
⎡ ܿ ܿ ௸௥ ⎤
1 ⎢ 1.6 ‫ ܽ⎛ ݎ‬− ቀ‫ݎ‬ቁ ⎞ ⎥ (1.51)
݂௟ = ⎢ ⎜ ௗ ோ ⎟ − 1⎥
2ߨ 0.1267 ܿ ௸௥
⎢ ܾ − ቀ‫ݎ‬ቁ ⎥
⎣ ⎝ ⎠ ⎦
9ߙ
݂ఈ = max ൬min ൬1; 5 − ൰; 0൰ (1.52)
ߨ
ߗܴ
߉=
ඥܷ ଶ + (ߗܴ)ଶ (1.53)

k = 0.12 , 0.8 ≤ a ≤ 1.2 , 0.8 ≤ b ≤ 1.2 , 0.8 ≤ c ≤ 1.2

22
 CHAPITRE 2 : CONCEPTION AERODYNAMIQUE

CHAPITRE 2 : CONCEPTION AERODYNAMIQUE

La première étape de la conception d’une éolienne est la détermination de la géométrie

de l’organe principale à savoir le rotor plus précisément les pales. Pour cela et dans ce
chapitre nous essayerons, selon les conditions de fonctionnement optimal et en utilisant les
théories décrites au chap.1, de déterminer la conception et les dimensions du rotor.

2.1 Design du rotor

Dans cette partie nous nous intéresserons à la détermination de la géométrie de la pale

qui permettra de transformer le maximum de l’énergie contenu dans le vent en énergie
mécanique avec le meilleur rendement possible.

Pour cela nous nous baserons sur la théorie élémentaire de Betz et la méthode BEM.
Globalement il sera questions d’utiliser un profil dans les caractéristiques sont connus, (S809
dans notre cas) [6] et de se donner des paramètres de fonctionnement (vitesse nominale du
vent, coefficient spécifique de vitesse optimale…) et d’essayer d’optimiser la conversion de
l’énergie cinétique du vent sur toute la longueur de la pale. Ceci nous permettra d’avoir la
forme de la pale de notre éolienne.

2.1.1 Paramètres liés au rotor

La puissance désirée est choisie pour une vitesse nominale du vent. A cela s’ajoute la
connaissance approximative du rendement mécanique de la chaine de conversion, et avec un
coefficient de puissance qu’on se donne pour commencer le calcul nous pourrons déterminer
le premier paramètre principal de notre optimisation qui est le rayon du rotor. L’expression de
la puissance est ainsi donnée par la relation

1 (2.1)
ܲ= ߩ ‫ܴߨ ߟ ܥ‬ଶܷ ଷ
2 ௔ ௣ ௠
Le rayon étant connu nous passerons au choix de la vitesse spécifique λ. Selon le type
d’application et le nombre de pales du rotor, la vitesse spécifique λ change. Le tableau ci-
dessous propose des valeurs suivant le nombre de pales.

24
 CHAPITRE 2 : CONCEPTION AERODYNAMIQUE
Tableau 2.1: Coefficients spécifiques de vitesse suivant le nombre de pales [4]

Nombre de pale B λ
8-24 1
6-12 2
3-6 3
3-4 4
1-3 >4

Nos éoliennes possédants trois pales, notre choix s’est porté sur λ=8. Les courbes de
variations du coefficient de puissance en fonction de λ mettent en évidence que pour 6 < λ < 9
le coefficient de puissance est maximale, comme illustre la (fig. 2.1)

Figure 2.1: Variation de la puissance en fonction de la vitesse spécifique [1]

2.1.2 Géométrie de la pale

On divise la pale en N éléments pour calculer les paramètres géométriques de chaque

section.

En considérant les caractéristiques du profil S809 (variation des coefficients de

portance et de traînée en fonction de l’angle d’attaque), nous pourrons choisir les paramètres
optimaux de design. Le fonctionnement de la pale serait optimal quand la traînée est

25
 CHAPITRE 2 : CONCEPTION AERODYNAMIQUE

minimale. On déduit que les paramètres de design Cddesign , Cldesign , αdesign doivent vérifier
la condition Cd/Cl est au minimum et ce pour chaque section de la pale (cas ou Cl dépend de
l’angle d’attaque et du nombre de Reynolds). Les relations ci-dessous permettent d’effectuer
les calculs appropriés ; à savoir

‫ݎ‬௜ (2.2)
ߣ௥௜ = ߣ
ܴ
2 1
߮௜ = atan ൬ ൰ (2.3)
3 ߣ௥௜
8ߨ‫ݎ‬௜
ܿ௜ = (1 − cos(߮௜))
‫ܥܤ‬௟ௗ௘௦௜௚௡,௜ (2.4)
ߚ௜ = ߮௜ − ߙௗ௘௦௜௚௡ (2.5)

2.1.3 Calcul des performances et modifications du design

Le calcul des performances consiste en la détermination de la courbe de puissance en

fonction de la variation de la vitesse et ce, en utilisant les paramètres réels de fonctionnement
(Cl , Cd , λ variant, α) .Ceci nous permettra d’opérer des modifications sur la géométrie de la
pale pour obtenir la forme optimale.

La détermination de la courbe de puissance, connaissant les caractéristiques du profil

et la géométrie de la pale déterminée précédemment, passe par la résolution des équations de
la méthode BEM en procédant comme suit :

Pour une station i

 On se donne une plage de variation de l’angle d’attaque α,

 on calcul pour chaque angle d’attaque le coefficient de portance en utilisant les
formules empiriques,
 On calcul Cl de la formule donnée par la méthode BEM, en tenant compte des
pertes en bout de pale eq.1.40 ;

cos ߮ − ߣ௜sin ߮ (2.6)

‫ܥ‬௟ = 4‫ ܨ‬sin(߮)
ߪ௜(sin ߮ + ߣ௜cos ߮)
Avec ߮ = ߙ + ߚ + ߚ଴

 on cherche les valeurs absolue des différences, entre les Cl calculé par l’équation (2.6)
et ceux calculé par les relations empiriques, soient inferieurs à 0.01, et on prend
l’angle d’attaque αoptimum correspondant au maximum des différences,

26
 CHAPITRE 2 : CONCEPTION AERODYNAMIQUE

 on fait une interpolation, dans le tableau donnant les coefficients de portance et de

traînée en fonction de l’angle d’attaque, pour la valeur de αoptimum,
 on calcul ߮௜ de la relation :

߮௜ = ߚ௜ + ߙ௢௣௧௜௠ ௨௠ (2.7)
 on applique le model de décrochage de Laino & hansen, pour déterminer les
coefficients Cl et Cd tridimensionnel, d’une pale en rotation
 on peut alors calculer :

8ߨ‫ݎ‬௜ (2.8)
ܿ௜ = (1 − cos ߮௜)
‫ܥ ܤ‬௟ଷௗ
‫ܥ‬௬
ܽ௜ = (2.9)
4 (sin ߮)ଶ + ‫ܥ‬௬
‫ܥ‬௬ = ‫ܥ‬௟ଷௗ sin ߮ − ‫ܥ‬ௗଷௗ cos ߮
(2.10)
ܽ௜ tan ߮௜
ܽᇱ௜ = (2.11)
ߣ௜
 calculons le couple :

ܳ௜ = 4ߗ ‫( ܨ‬1 − ܽ)ܽᇱߩ௔ ܷߨ‫ݎ‬ସ (2.12)

 Calcul de la puissance :

ܲ௜ = ܳ ∗ ߗ (2.13)

Les figures 2.2 et 2.3 résument les principales étapes de calcul sous forme d’organigramme.

27
CHAPITRE 2 : CONCEPTION AERODYNAMIQUE

Entrer les données P, λ, Cp, ߟm, U, B, N

Entrer les caractéristiques du profil

Calculer le rayon R équation (2.1)

Subdiviser la pale en N éléments et calculer les rayons ri des

stations

Calculer pour chaque stations les paramètres λri, Φi, ci, βi

équations (2.1)(2.2)(2.3)(2.4)

Calculer ai,1 a’i,1 équations (2.14)(2.15)

Calculer Φi,j , Fi,j équations (2.5)(2.6)

ri < R

Figure 2.2: organigramme de détermination de la géométrie

28
 CHAPITRE 2 : CONCEPTION AERODYNAMIQUE

Introduire les données c, r, tableau1 [α, Cl(α), Cd(α)], R, Ω, ߩ௔

Proposer un domaine de variation de l’angle d’attaque α

Calculer, avec les formules empiriques propres au profil, Cl1 (α)

dans le domaine proposé ci-dessus

i=1 à dimension (c)

Calculer F, Cl2 équations (1.39) (2.6)

Trouver le vecteur H=α(k), avec k positions des éléments du vecteur |∆‫ܥ‬௟| inferieur à
0.01

αoptimum= max(H)

Déterminer du tableau 1 les valeurs de Cl et Cd correspondant à αoptimum

i ≤ dimension (c)

Utiliser le model de décrochage de Laino & Hansen pour déterminer Cl3d

et Cd3d

i=1 à dimension (c)

Calculer a et a’ des équations (2.9) (2.11)

Calcul du couple de l’équation (2.12)

Calculer la puissance locale

i ≤ dimension (c)

Calculer la puissance totale pour la vitesse choisie

Figure 2.3: organigramme d'évaluation de la puissance

29
 CHAPITRE 2 : CONCEPTION AERODYNAMIQUE

2.3 Résultats et validations des codes de calculs

2.3.1 Programme de géométrie de la pale

Les résultats de la variation de la corde et de l’angle de vrillage sont présentés dans le tableau
suivant :

Tableau 2.2: Géométrie de la pale calculé par le code

r (m) c (m) β (°) Φ (°)

0.4780 0.8011 26.4512 32.6512
0.6264 0.7507 21.3047 27.5047
0.7747 0.6850 17.3633 23.5633
0.9231 0.6204 14.3085 20.5085
1.0714 0.5622 11.8998 18.0998
1.2197 0.5115 9.9661 16.1661
1.3681 0.4677 8.3871 14.5871
1.5164 0.4299 7.0776 13.2776
1.6648 0.3972 5.9764 12.1764
1.8131 0.3687 5.0389 11.2389
1.9615 0.3438 4.2321 10.4321
2.1098 0.3219 3.5311 9.7311
2.2582 0.3025 2.9166 9.1166
2.4065 0.2852 2.3739 8.5739
2.5549 0.2697 1.8912 8.0912
2.7032 0.2558 1.4593 7.6593
2.8516 0.2431 1.0707 7.2707
2.9999 0.2317 0.7191 6.9191
3.1483 0.2212 0.3997 6.5997
3.2966 0.2117 0.1082 6.3082

Pour valider le code nous avons comparé les résultats dans le cas d’utilisation des paramètres
de design d’une pale parfaite à savoir : Cl=1, Cd=0, λ=6 et F=1 avec les résultats tirés de la
référence [3].

Les résultats de notre code pour une subdivision de N=9 sont présentés par le tableau 2.3
30
 CHAPITRE 2 : CONCEPTION AERODYNAMIQUE
Tableau 2.3: Géométrie de la pale pour le cas idéal (Cl=1, Cd=0)

r/R c/R ߮ (݀݁݃°)

0.2000 0.1765 26.5370
0.3000 0.1423 19.3697
0.4000 0.1154 15.0799
0.5000 0.0960 12.2900
0.6000 0.0818 10.3494
0.7000 0.0711 8.9283
0.8000 0.0627 7.8455
0.9000 0.0561 6.9943
1.0000 0.0507 6.3082

31
 CHAPITRE 2 : CONCEPTION AERODYNAMIQUE

C/R (a)
. Présent travail

Autre travail antérieur [4]

r/R

߮ (݀݁݃°) (b)

. Présent travail

Autre Travail [4]

r/R

Figure 2.4 : Évolution de la corde et l’angle relatif en fonction du rayon : comparaison avec les travaux
antérieurs [4]

Comme illustre les figures 2.2 (a et b) les résultats sont assez proches, et donc le code est
validé.

32
 CHAPITRE 2 : CONCEPTION AERODYNAMIQUE

2.3.2 Courbe de puissance

P (W)

V (m/s)

Figure 2.5: Courbe de la puissance

La figure 2.3 présente les résultats de la courbe de puissance. L’éolienne a été

dimensionnée pour une puissance de 25 kW. Nous remarquons de la courbe de puissance que
le maximum est atteint autour de 22.5 kW pour une vitesse de 15 m/s soit une perte de 2.5
kW. Ceci est du aux hypothèses utilisées lors de la détermination de la géométrie. Notons que
nous n’avons pas considéré l’effet de la viscosité de l’air qui se traduit par un frottement sur
les pales.

Nous déduisons du graphe que la production de puissance commence des que la

vitesse du vent atteint les 6 m/s.

33
 CHAPITRE 3 : LOGICIEL DE SIMULATION SOLIDWORKS

CHAPITRE 3 : LOGICIEL DE SIMULATION SolidWorks

3.1 Qu’est ce que le logiciel SolidWorks ?

SolidWorks est un logiciel de dessin et conception assisté par ordinateur DAO et

CAO. Il est considéré comme un outil essentiel dans le domaine de la conception mécanique,
et nombreuses entreprises et bureaux d’études en mécanique de renommé l’utilisent. Il
permet, avec ses diverses fonctionnalités, d’aborder des problèmes de dimensionnement
complexes et de les résoudre dans des délais très courts, comparés aux méthodes
conventionnelles. Il est apprécié par des ingénieurs de conception pour sa souplesse et son
potentiel en constante progression grâce au programme de développement.

3.2 Potentiel de SolidWorks

3.2.1 Création de formes volumiques

L’interface DAO de SolidWorks permet la création de formes en 3D de la manière la

pus facile et la plus rapide. Il fait appel à des fonctions comme la révolution, l’extrude,
enlèvement de matière ou la symétrie pour obtenir des formes 3D à partir de dessin plan. Pour
une même forme il peut y avoir différentes manières de procéder pour sa création. Citons
l’exemple simple d’un cylindre qui peut être obtenu par révolution ou extrude :

 Révolution : obtenir un cylindre à partir d’une esquisse de rectangle en le faisant

tourner autour d’un axe, figure 3.1

Figure 3.1 : Obtention d’un cylindre par révolution d’un rectangle

35
 CHAPITRE 3 : LOGICIEL DE SIMULATION SOLIDWORKS

 Extrude : obtenir un cylindre par extrusion de la surface d’un cercle, Figure 3.2

Figure 3.2 : Obtention d’un cylindre par extrusion d’une esquisse de cercle

3.2.2 Effectuer des assemblages

A partir des formes volumiques, on peut créer des assemblages de mécanismes ou

structure en introduisant les contraintes appropriées.

Exemple d’assemblage : assemblage de deux pièces par boulons.

 Contrainte 1 : coaxialité entre le boulon et les perçages des pièces ;

 Contrainte 2 : Coïncidences des deux surfaces de contact des pièces ;
 Contrainte 3 : Coïncidence de la surface basse de la tète du boulon et la surface de
la pièce

3.2.3 Conception, analyse et simulation

SolidWorks était à la base un logiciel de modélisation 3D, mais avec le temps et le

développement de l’outil numérique, il a été doté de plusieurs compléments dont on citera :
CosmosWorks, MotionWorks, FlowWorks, SolidWorks-Piping, SolidWorks Animator…

 SolidWorks-Piping : interface de conception de tuyauteries doté d’une

bibliothèque d’accessoires normalisés ;
 SolidWorks Animator : interface permettant la création d’animation de
mouvement des produits de conceptions ;

36
 CHAPITRE 3 : LOGICIEL DE SIMULATION SOLIDWORKS

 FlowWorks : interface permettant la simulation d’écoulements et des

phénomènes de transfert de chaleur ;
 MotionWorks : interface d’étude des mouvements d’assemblages ;
 CosmosWorks : interface de calcul de structures intégré à SolidWorks. Il
permet différents types d’analyse tels que l’analyse statique (déplacements,
déformations et contraintes), modale (fréquences propres et déformées
modales), du flambage(charges critiques de flambage) ou de l’optimisation
dimensionnelle.
SolidWorksS offre la possibilité de développer des produits plus rapidement, pour un
coût inférieur et de manière plus efficace qu'avec des prototypes, aussi précis soient ils, tout
en permettant d'étudier les performances obtenues avec un bien meilleur niveau de détail.

Quelle que soit l'application industrielle, de l'aéronautique au secteur médical,

SolidWorks offre des avantages significatifs en termes de qualité des produits, car le logiciel
permet aux ingénieurs d'aller au-delà des simples calculs manuels et de vérifier la faisabilité
de leurs conceptions. Une analyse rapide et peu coûteuse aboutit souvent à des solutions qui
n'avaient pas été envisagées initialement et permet aux ingénieurs de mieux maîtriser les
différentes caractéristiques des produits.
Les éventuels dysfonctionnements ou défaillances des produits peuvent entraîner des
rappels coûteux et des problèmes de responsabilité. CosmosWorks présente le fonctionnement
d'un modèle dans des conditions réelles, avant sa fabrication. Les ingénieurs peuvent ainsi
prévoir le comportement réel des pièces et des assemblages dans toutes les conditions de
chargement possibles.
SloidWorks permet aux concepteurs de répondre aux spécifications fonctionnelles
requises sans sur dimensionner et donc sans gaspillage de matière. L'économie de petites
quantités de matière en apparence insignifiantes multipliée par des dizaines de composants
peut aboutir à des réductions de coûts considérables au niveau de la production, de
l'emballage et de l'expédition.
Il permet d’éviter l’étape du prototype. Près de 80% des coûts de fabrication d'un
produit sont liés à la validation de la conception. La possibilité de réaliser des itérations de
conceptions rapides et peu coûteuses avant l'approbation finale du modèle offre donc un
avantage concurrentiel déterminant. SolidWorks permet de réaliser des itérations de
conceptions rapides et peu coûteuses sur des modèles virtuels plutôt que sur des prototypes
réels coûteux.

37
 CHAPITRE 3 : LOGICIEL DE SIMULATION SOLIDWORKS

3.3 Exemple d’utilisation de SolidWorks

Nous prenons pour exemple la pale dont les dimensions ont été déterminées au
chapitre 2. La première étape est le dessin qui peut être fait en 4 étapes.

 Étape 1 : l’étape 1 consiste à préparer les plans sur lesquelles nous allons dessiner les
20 profils obtenus par le code de calcul fig.3.3 (le profil S809 peut être dessiné avec
40 points dont les abscisses et ordonnées dépendent de la corde tableau annexe).
Pour obtenir les plans il suffit de maintenir la touche contrôle du clavier enfoncée et
faire glisser le plans de face et indiquer la distance voulue. On fera de même pour tout
les plans

Figure 3.3: Représentation des 6 premiers plans correspondants aux 6 premiers profils

 Etape 2 : une fois les plans dessinés nous commençons à esquisser les points
déterminants le contour du profil S809 correspondant à la corde à la racine. Fig 3.4 et

38
CHAPITRE 3 : LOGICIEL DE SIMULATION SOLIDWORKS

3.5. Après avoir dessiné les 40 points nous les relions à l’aide du spline et nous
effectuons une rotation d’une valeur correspondante à l’angle de vrillage.

Figure 3.4: Esquisse du premier point

Figure 3.5: Forme finale du profil à la racine avec angle de vrillage

39
CHAPITRE 3 : LOGICIEL DE SIMULATION SOLIDWORKS

 Etape 3 : cette étape nous permettra d’obtenir la forme de la partie active de la pale
(calculée par le code au chapitre 2) fig.3.7 à l’aide de l’option lissage de surface. Il
suffit de sélectionner les 20 profils en essayant d’avoir un alignement des points
apparaissant en bleu sur la figure 3.6.

Figure 3.6: Lissage entre les deux premiers profils

40
 CHAPITRE 3 : LOGICIEL DE SIMULATION SOLIDWORKS

Figure 3.7 : Forme finale de la partie active de la pale

Etape 4 : Nous allons maintenant dessiner la partie qui servira de liaison entre la partie active
de la pale et le moyeu. Nous commencerons par dessiner une esquisse carré sur un plan
distant du plan 1 et lui faire une extrude fig 3.7. Nous utiliserons après l’outil perçage pour
effectuer des alésages permettant le passage des boulons de fixations. Pour lier ces deux
volume nous utilisons la fonction lissage

Figure 3.8: support de la pale

41
CHAPITRE 3 : LOGICIEL DE SIMULATION SOLIDWORKS

Figure 3.9: Forme final de la pale

42
 CHAPITRE 4 : CONCEPTION MECANIQUE

CHAPITRE 4 : CONCEPTION MECANIQUE

4.1 Schéma cinématique et organes constitutifs

La figure 4.1 est une représentation du schéma cinématique d’une éolienne à axe
horizontal.

Arbre lent
Accouplement 2
Moyeux + pales

M G

Frein Accouplement 3

M : multiplicateur
M : multiplicateur
Accouplement 1
Paliers G : génératrice
G : génératrice

Figure 4.1:Schéma cinématique de l’éolienne à axe horizontale

 Pales : au nombre de 2 ou 3 selon la configuration. La forme aérodynamique spéciale

permet la transformation de l’énergie cinétique contenue dans le vent en une énergie
mécanique de rotation. Le matériau utilisé dans les constructions actuelles est un
matériau composite à base de fibre de verre ;
 Moyeu : de formes différentes selon la conception, permet de supporter les pales et de
transmettre la rotation à l’arbre. Les moyeux sont généralement en fonte ;
 Arbre : de forme cylindrique, transmet la rotation au multiplicateur ;
 Génératrice : les deux types de génératrices sont utilisés (synchrone et asynchrone).
L’utilisation de l’une ou de l’autre dépend du choix du concepteur;

43
 CHAPITRE 4 : CONCEPTION MECANIQUE

 Multiplicateur : l’utilisation d’une génératrice asynchrone entraîne l’utilisation d’un

multiplicateur. En effet les génératrices de types asynchrone fonctionnent à des
fréquences de rotation élevées tandis que le rotor tourne à une vitesse assez faible ;
 Paliers : permettent le guidage en rotation ;
 Accouplements : permettent la liaison entre les différents arbres ;
 Frein : le rotor ayant une vitesse de rotation limite pour éviter l’emballement de ce
dernier, l’utilisation du frein permet de contrôler la vitesse de rotation du rotor par le
biais d’un compteur de tour ;
 La nacelle : réceptacle contenant les différents organes cités précédemment pour les
protéger ;
 La tour : supporte la nacelle, de forme différente. La figure 4.2 ci-dessous illustre les
différents configurations ;

Figure 4.2: Types de support de nacelles- tours ou mats

 Éléments de fixation et assemblages : ensemble de boulonnerie, clavetage, cannelures

4.2 Etude technologique

 Rotor : on appelle rotor l’ensemble moyeu-pales.

Le moyeu : en générale obtenu par moulage. La forme la plus simple est cylindrique
avec des dégagements ou méplats permettant de recevoir les pales. La forme
cylindrique du moyeu peut être optimisée pour obtenir une large variété de modèles.
La figure 4.3 illustre la forme d’un moyeu;

44
 CHAPITRE 4 : CONCEPTION MECANIQUE

Figure 4.3: Moyeu du rotor de l'éolienne

Les pales : la forme aérodynamique complexe rend la réalisation difficile. La
réalisation de nos pales se fera par moulage de 20 pièces de la forme du profil S809
dont les cordes ont été calculées au « chapitre 2 ». Ces 20 pièces seront reliées entre
elle par des tiges et le tout recouvert d’une couche de résine permettant d’obtenir la
forme aérodynamique et l’état de surface souhaitée, l’épaisseur de la couche sera
déterminée par le calcul dynamique ;

Figure 4.4:Profils de la pale avant la couche de résine

 Les paliers : l’arbre étant soumis à une force axiale provenant de la poussée
développée par le rotor, il est préférable d’utiliser des paliers à roulement conique. Il
existe deux possibilités de montage: montage en X et montage en O ;

45
 CHAPITRE 4 : CONCEPTION MECANIQUE

 Montage arbre-rotor : deux solutions peuvent êtres adoptées : montage par clavette
et montage par cannelures. Le montage par cannelures est utilisé lorsque la condition
limite d’utilisation des clavette n’est pas satisfaite (l > 2Da l : longueur de la
clavette, Da : diamètre de l’arbre). Pour notre montage nous pourrons, avec les calculs
dynamique, décider du montage qui sera adopté ;
 Accouplement : Les accouplements sont utilisés pour transmettre la vitesse et le
couple, ou la puissance, entre deux arbres de transmission en prolongement l'un de
l’autre. Il existe plusieurs types d’accouplement. Notre choix s’est porté sur un
accouplement élastique du fabricant Rotex [9] afin d’assurer une protection de la
génératrice et du multiplicateur, en cas de fluctuation de la charge ou présence de
vibrations. Le choix de l’accouplement dépend du couple à transmettre et de la vitesse
de rotation.
 Frein : sert à contrôler la vitesse de rotation du rotor et éviter l’emballement. Il peut
être disposé au niveau de l’arbre ou après le multiplicateur. Le couple après le
multiplicateur étant inferieur à celui au niveau de l’arbre lent, le frein sera disposé au
niveau de l’arbre de sortie du multiplicateur ;
 L’arbre : il doit contenir une partie fileté pour recevoir le nez du rotor qui permettra
le blocage axiale du moyeu. Cette partie sera suivie par un logement de clavette ou
des cannelures. L’arbre contiendra éventuellement des épaulements pour le blocage
des roulements ;
 La tour : en considérant le critère d’esthétique sachant que l’éolienne est destinée à
être installé sur un toit, nous avons opté pour l’utilisation d’une tour tubulaire
haubanée ;

4.3 Calcul dynamique

Le calcul que nous allons entreprendre permettra d’avoir les dimensions approximatives
que nous utiliserons pour notre conception. Les simulations par le logiciel apporterons les
corrections et permettrons une optimisation quand c’est nécessaire.

4.3.1 Les pales

La géométrie de la pale obtenue par le calcul aérodynamique est assez complexe et un

calcul analytique s’avère difficile.

46
 CHAPITRE 4 : CONCEPTION MECANIQUE

Figure 4.5: Vue 3D de la pale

Nous adopterons une approche simplificatrice dont nous vérifierons les résultats en
simulation par SolidWorks.

Nous considérerons que la pale est en forme de tronc de cône et d’épaisseur variable
dépendant des contraintes au niveau de chaque section. Les charges appliquées sur la pale
sont les forces aérodynamiques et la force gyroscopique.

De l’algorithme de détermination de la courbe de puissance nous connaissons les

forces appliquées aux vingt stations ainsi que la distance par rapport au centre de l’hélice.
Nous pouvons alors calculer la contrainte de flexion et la contrainte normale appliquée a
chaque section « i » caractérisée par :

 di diamètre du cercle égal à la corde de la station i ;

 ei épaisseur de l’anneau ;
 ∆r : espacement entre section ;

La condition de résistance d’une section au chargement est :

ߪ௘ (4.1)
ߪ௫ <
ܿ‫ݏ‬
σe : limite élastique en extension du matériau CS : coefficient de sécurité.

‫ ܯ‬௙݀ ‫ܨ‬௖ (4.2)

ߪ௫ = ߪ௙ + ߪ௡ = +
2‫ܣ ܫ‬

47
 CHAPITRE 4 : CONCEPTION MECANIQUE

Mf : moment de flexion ; Fc : force centrifuge ; I : moment quadratique de la section

ds

∆r ∆r

Di D
Fc
x

F1 F2 Fi F20
Figure 4.6: modélisation du chargement sur la pale

Fi sont les résultantes des forces tangentielle et axiales pour chaque section « i ».

(4.3)
‫ܨ‬௜ = ට‫ܨ‬௧ଶ௜ + ݂௔௜
ଶ

(ܴ − ‫ݎ‬௠ − ݀‫)ݏ‬ (4.4)

∆‫= ݎ‬
20
rm rayon du moyeu ;

a : longueur du support de la pale ;

 Calcul du moment fléchissant

Pour une section donnée « i » :

௜ିଵ

‫ݎ∆݅[ ∈ ݔ‬, (݅+ 1)∆‫]ݎ‬ ‫ ܯ‬௙௜ = ቌ෍ (݅− ݇)‫ܨ‬௞ቍ ∆‫ݎ‬ (4.5)

௞ୀଵ

Sachant que, pour une station i:

ߨ ସ (4.6)
‫ܫ‬௜ = (݀ − (݀௜ − 2݁௜)ସ)
64 ௜
ߨ(݀௜ଶ − (݀௜ − 2݁௜)ଶ) (4.7)
‫ܣ‬௜ =
4
L’expression de la contrainte de flexion s’exprime par :

48
 CHAPITRE 4 : CONCEPTION MECANIQUE

൫∑௜௞ୀଵ
ିଵ
(݅− ݇)‫ܨ‬௞൯∆‫݀(ݎ‬௜ − 2݁௜) (4.8)
ߪ௙ = ߨ ସ
(݀ (݀ )ସ)
32 ௜ − ௜ − 2݁௜

La condition de résistance s’écrit alors :

൫∑௜௞ୀଵ
ିଵ
(݅− ݇)‫ܨ‬௞൯∆‫݀(ݎ‬௜ − 2݁௜) ‫ܨ‬௖ ߪ௘ (4.9)
ߨ ସ +ߨ <
(݀ (݀ )ସ) (݀ଶ (݀ )ଶ) ‫ܵܥ‬
32 ௜ − ௜ − 2݁௜ 4 ௜ − ௜ − 2݁௜
La détermination des valeurs de l’épaisseur « e » se fait par calcul itératif.

Les données nécessaires au calcul ainsi que les résultats sont détaillés en annexe A.

Remarque : Le diamètre D du pied de la pale doit être vérifié au cisaillement.

∑ଶ଴
௜ୀଵ ‫ܨ‬௜ ߬௔ௗ௠ (4.10)
߬= ߨ <
‫ܦ‬ଶ ‫ܵܥ‬
4
τadm contrainte admissible en cisaillement.

4.3.2 L’arbre

L’arbre est soumis aux contraintes de compression (force de poussée du rotor),

contrainte de flexion (poids du rotor et poids propre de l’arbre), et une contrainte de torsion
due au couple de l’hélice. La figure (4.7) schématise le chargement de l’arbre.

Pr
Palier2 Palier1
Fp
C

Pa
Figure 4.7:Chargement de l'arbre
 Calcul de la contrainte normale due à la compression

‫ܨ‬௣ (4.11)
ߪ௡ =
‫ܣ‬
Fp : force de poussée du rotor (résultante des forces axiales) A : section de l’arbre,

 Calcul de la contrainte de flexion

49
 CHAPITRE 4 : CONCEPTION MECANIQUE

Pr R1y Pa R2y
Fa

Q
b a a b

Figure 4.8:Schématisation du chargement sur l'arbre

෍ ‫=ܨ‬0 , ෍ ‫ ܯ‬/௢ = ܱ (4.12)

ܴଵ௬ + ܴଶ௬ = ܲ௔ + ܲ௥
ቊ

ܴଵ௬ . ܿ+ ܴଶ௬ . (2ܽ + ܾ) = ܲ(ܽ + ܾ) (4.13)
Avec 2a+2b=L, a=2b

La résolution du système d’équation nous donne :

ܲ 5
ܴଵ௬ = + ܲ
൞ 2 4 ௥
(4.14)
1 1
ܴଶ௬ = ൬ܲ − ܲ௥൰
2 2

Les moments de flexion sont :

‫ =ݔ‬0 , ‫ܯ‬௙ = 0
൜

‫ܾ =ݔ‬ , ‫ ܯ‬௙ = −ܲ௥. ܿ

ܲ௥
‫ ܾ = ݔ‬+ ܽ , ‫ ܯ‬௙ = ൬ܲ − ൰ܿ

ቐ 2
‫ ܾ = ݔ‬+ 2ܽ , ‫ܯ‬௙ = 0

Le diagramme des moments fléchissant indique que la section dangereuse se situe au niveau
du palier 1. De ce fait la condition de résistance d’après le critère de Tresca s’écrit :

4‫ܨ‬௣ 32‫ ܯ‬௙ ଶ 32ܳ ଶ ߪ௘ ଶ (4.15)

൬ ଶ+ ൰ + ൬ ൰ < ቀ ቁ
ߨ݀ ߨ݀ଷ ߨ݀ଷ ‫ܵܥ‬

Le calcul du diamètre se fera par un calcul itératif.

50
 CHAPITRE 4 : CONCEPTION MECANIQUE

4.3.3 Montage pales-moyeu

Comme l’indique la figure 4.9, chaque pale est reliée au moyeu par le biais de deux
boulons. Nous indiquerons dans ce qui suit la procédure à adopter pour ce calcul.

Figure 4.9: Vus éclatée et coupe du montage moyeu-pales

[Link] Calculs du diamètre de la tige

Les quatres tiges sont soumises en cisaillement dû aux forces Fa , Pp , Fc .

Fa : force axiale due à l’écoulement autour de la pale ;

Fc : force centrifuge ;

Pp : poids de la pale ;

Toutes ces forces se trouvent sur un même plan. On en déduit la résultante

ଶ (4.16)
‫ = ܨ‬ට ൫ܲ௣ + ‫ܨ‬௖൯ + ‫ܨ‬௔ଶ

Comme nous avons quatre boulons, chaque tige doit alors supporter la moitié de la force F.

La contrainte de cisaillement est donnée par la formule suivante :

51
 CHAPITRE 4 : CONCEPTION MECANIQUE

‫ܨ‬ ߨ݀ (4.17)
߬= ‫ܣ‬௧ = ∶ ‫݁݃݅ݐ ݈ܽ݁݀ ݊݋݅ݐܿ݁ݏ‬
4‫ܣ‬௧ 4

Pour que la tige résiste au cisaillement, il faut que la contrainte de cisaillement soit inferieur à
la contrainte pratique du matériau moyennant un coefficient de sécurité.

Le diamètre est alors déduit de l’équation ci-dessous :

4‫ܨ ܵܥ‬
݀ ≥ ඨ
ߤߪ௘௚ (4.18)

σeg : limite élastique en cisaillement (glissement).

Le diamètre étant déterminé, nous pouvons choisir dans les tables normalisé le boulons que
nous allons utiliser et connaître toutes ces dimensions.

[Link] Calcul du couple de serrage

Dans les assemblages boulonnés, on évite de contraindre les boulons en cisaillement.

Pour cela nous allons calculer le couple de serrage permettant une adhérence suffisante entre
les surfaces en contacte (Fad) qui soit supérieur à la force provocant la contrainte de
cisaillement, à savoir la force F.

H H 3
8 2
P= pas ; H = p = 0,866P
2
H
H
D1 = d1 = d2 - = d - 1,0825P
2
D H
Vi d 3H
D2 = d2 = d - = d - 0,6495P
D2 s 4

D1 P
d
d d
d3 = d - 2 H = d - 1,2268P
H
3
4

Le couple de serrage est donné par la formule suivante :

‫ܳ = ܥ‬ൣ݂‫ݎ‬௠ ௢௬ + (݂ᇱ + tan ߙ)‫ݎ‬௠ᇱ௢௬ ൧ (4.19)

f : coefficient de frottement entre l’écrou et la pièce

52
 CHAPITRE 4 : CONCEPTION MECANIQUE

f’ : coefficient de frottement entre les filets ( vis \ écrou )

rmoy rayon moyen de la surface plane de l’écrou r2
r’moy rayon de l’hélice moyenne
r1

r23  r13 d
rmoy = r’moy = – H (H : hauteur du pas)
r22  r12 2
L’adhérence entre les surfaces est assurée par la force Q.

Sachant que µ est le coefficient de frottement entre les pièces

en contacte, la force d’adhérence s’exprime par :

‫ܨ‬௔ௗ = ߤܳ > ‫ܨ‬ (4.20)

[Link] Résistance à la traction, cisaillement et matage

 Traction au niveau de la tige :

ܳ ߪ௘ (4.21)
ߪ௡ = <
ܵ௧ ‫ܵܥ‬
 Cisaillement de la tete du boulon

Contrainte de cisaillement :
ܳ ߪ௘௚ (4.22)
ߪ= <
ܵ ‫ܵܥ‬
S = dh
σeg : résistance elastique au cisaillement.
ܳ. ‫ܵܥ‬ (4.23)
ℎ≥
ߨ݀ ߪ௘௚
C’est la valeur la plus grande de h qui est prise en considération.
 matage des filets

La vérification de la résistance au matage permet de déduire le nombre de filet nécessaire.

ܳ (4.24)
ܲ௖ = < ܲ௔ௗ௠
ܵ௣
Pc : pression de contacte ;

53
 CHAPITRE 4 : CONCEPTION MECANIQUE

Sf : surface totale de filet en prise ;

݀² − ݀ଶଶ (4.25)
ܵ௣ = ݊ߨ
4
n : nombre de filets

Padm : pression admissible ;

 cisaillement des filets

ܳ ߪ௘௚ (4.26)
߬= <
ܵ௖ ‫ܵܥ‬
Sc : surface cisaillée ܵ௖ = ݊ߨ݀ଶ‫( ݌‬p : pas de filetage)

Le nombre des filets boulons est : N  n + a ; (écrou : a = 2 vis : a = 4)

La longueur filetée de la tige est : L = N.p

4.3.4 Montage moyeu-arbre

[Link] Montages par clavette

Le calcul de clavette consiste à déterminer la longueur nécessaire permettant de

supporter la pression due aux forces de contactes. Les dimensions des sections transversales
des clavettes sont normalisées et dépendent du diamètre de l’arbre.

Figure 4.10: montage par clavette [6]

La résultante des force de contacte peut être déterminé du couple moteur et du rayon
de l’arbre.

54
 CHAPITRE 4 : CONCEPTION MECANIQUE

‫ܥ‬ெ ‫ݎݑ݁ݐ݋ ݈݉݁݌ݑ݋ܥ‬ (4.26)

ܶ=
ܴ ∶ ‫݈݁݀ ݊݋ݕܽݎ‬ᇱܽ‫݁ݎܾݎ‬
où ‫ܥ‬ெ : ‫ݎݑ݁ݐ݋ ݈݉݁݌ݑ݋ܥ‬et ܴ ∶ ‫݈݁݀ ݊݋ݕܽݎ‬ᇱܽ‫݁ݎܾݎ‬
Nous pouvons alors calculer la pression de contacte p de la formule :

ܶ ܶ (4.27)
‫=݌‬ =
ܵ ݈∗ ݁

Avec e : épaisseur du flanc de la clavette e ≈ b/2

Sachant que la pression de contacte doit être inferieur a la pression admissible (Padm)
du matériau nous pouvons déduire la longueur de la clavette :

‫ܥ‬ (4.28)
݈≥
ܴ ∗ ݁∗ ܲ௔ௗ௠

[Link] Montages par cannelure

Le montage par cannelure est utilisé pour transmettre des couples importants. Le
nombre de cannelure et la section des flancs dépendent du diamètre de l’arbre et de la série de
cannelure dont on dénombre trois : la série faible, moyenne et forte. La figure 4.11 indique le
nombre et la valeur de la surface s’ (surface réelle d’appui par mm de longueur de contacte)
en fonction du diamètre et de la série.

Figure 4.11: Détermination de s’[7]

La pression de contacte s’exprime par la relation :

55
 CHAPITRE 4 : CONCEPTION MECANIQUE

‫ܥ‬ (4.29)
ܲ=
‫ݏ‬ᇱ݈ܴ௠ ௢௬

Rmoy : rayon moyen mesuré a mi-hauteur d’une dent.

Figure 4.12: Caractéristiques géométriques [7]

Remarque : Le calcul des clavette et cannelure consiste en la détermination de la langueur l

qui doit vérifier la condition suivante : l < 2,5 D. dans le cas de calcul de clavette si cette
condition n’est pas vérifié, le calcul doit être tourné vers un calcul de montage par cannelure.
Si la condition n’est toujours pas vérifiée il faudra alors changer la série ou le nombre de
cannelures.

4.3.5 La tour

Dans les constructions d’éolienne, la hauteur du mat qui supporte la nacelle est
équivalente à 1-1,5 le diamètre du rotor. Pour notre conception nous avons approximé la
hauteur à 10 m. Pour faciliter la réalisation et l’installation, la tour sera former de deux partie
de hauteur h1 et h2 tel que :

4
ℎଵ = ℎ ݁‫ݐ‬ℎଵ + ℎଶ = 10݉
5 ଶ

La tour est soumise au poids Pn de la nacelle, la force de poussé Fp du rotor et à la

force du vent. La figure ci-dessous modélise le chargement de la tour.

56
 CHAPITRE 4 : CONCEPTION MECANIQUE

y
h1 h2

D1
x
Pn D2

Fp
Force du vent

Figure 4.13: Schématisation du chargement de la tour

[Link] résistances à la flexion-compression

 Contrainte normale

ܲ௡ ܲ
⎧ ߪ௡ଵ = =ߨ ‫ܮ ݎݑ݁݅ݎ݁݌ݑݏ݁݅ݐݎܽ݌‬ଵ
⎪ ܵଵ (‫ܦ‬ ଶ
− ‫ܦ‬ ଶ
)
4 ଵ ௜


⎨ߪ = ܲ௡ + ߩଵ ݃ܵଵ‫ܮ‬ ଵ ܲ௡ + ߩ ଵ݃ܵଵ‫ܮ‬ ଵ
⎪ ௡ଶ = ߨ ‫ܮݎݑ݁݅ݎ݂݁݊݅݁݅ݐݎܽ݌‬ଶ
ܵଶ (‫ܦ‬ ଶ
− ‫ܦ‬ ଶ
)
⎩ 4 ଶ ௜

ρ1 : masse volumique du matériau de la tour

Di : diamètre de passage du câblage électrique.

 Contrainte de flexion

Le vent applique une pression Pv=ρ V²/2 sur la tour

ρ : masse volumique de l’air

La contrainte de flexion s’exprime par :

Partie supérieure h1

ߩߨܸଶℎଵ‫ܦ‬ଵ (4.30)
32 ‫ܮ‬ଵ ൬‫ܨ‬௣ + ൰
4
ߪ௙ଵ =
ߨ‫ܦ‬ଵଷ

57
 CHAPITRE 4 : CONCEPTION MECANIQUE

Partie inférieure h2

9 7 ߩߨܸଶℎଵ‫ܦ‬ଵ ߩߨܸଶℎଶ‫ܦ‬ଶ
32 ‫ܮ‬ଶ ൬ ‫ܨ‬௣ + 20 + ൰
5 4 4
ߪ௙ଶ = (4.31)
ߨ‫ܦ‬ଶଷ

Un calcul numérique nous permettra de déterminer les diamètres D1 et D2

[Link] Vérifications au flambement

Le problème de calcul des tubes en acier est très complexe. Dans la pratique, il peut, sauf
rares exceptons, être ramené à quelques cas simples. Ceux qui sont étudiés dans les règles CM
1956 suffisent ordinairement. Nous transposons ici les données de ce règlement pour leur
application limitée à des tubes isolés comme il est le cas pour notre conception. [8]

Le calcul consiste à vérifier que la somme des taux de fatigues en compression simple et
en flexion simple amplifiés chacun de son coté par des coefficients différents kn et kf, est
inferieur à la fatigue admise, elle-même réduite par u coefficient (1-0.2K0).les tableaux joints
en annexe (A) permettent de calculer ces Kn , Kf et K0 .[8]

La formule de vérification est :

‫ܭ‬௡ ߪ௡ + ‫ܭ‬௙ߪ௙ < (1 − 0.2‫ܭ‬଴)ߪ (4.32)

 a) on cherche dans le tableau a[8] le cas de flexion correspondant à la configuration de

chargement et on note les valeurs de l, A et M0.
 b) On détermine les valeurs de :

ܲ (4.33)
ߪ௡ =
ܵ
‫ ܯ‬଴ܴ (4.34)
ߪ௙ =
‫ܫ‬
‫ ܯ‬଴ܴ (4.35)
ߪ௙ =
‫ܫ‬
γ : élancement ; r : rayon de giration ;

En fonction de l’élancement on détermine du tableau K1a[8] la contrainte critique d’Euler nc.

CS étant le coefficient de sécurité on détermine la quantité :

58
 CHAPITRE 4 : CONCEPTION MECANIQUE
݊௖ (4.36)
ߤ=
‫ߪ ∗ ܵܥ‬௡

Du tableau K2b[8], en fonction de μ on tire les valeurs de Kn et Kf et le coefficient K0 en

fonction de μ et de A.

[Link] Calculs du diamètre des câbles de soutien

Comme la tour sera installée sur un toit de

bâtiment, nous ne pouvons pas utiliser un mat
autoporteur car il nécessite une fondation robuste.
Pour cela nous fixerons la tour au sol par boulons
et nous ajouterons trois câbles disposés à 120°
dont nous allons présenter la méthode de calcul
du diamètre.

La determination du diametre des cables passe par

le calcul des tensions T1 et T2 qui sont dus aux
forces :

Fp : la force de poussé du rotor ;

Figure 4.14: Schématisation du chargement de la
tour haubanée
F1 : la resultante de la force du vent sur la partie
superieur de la tour ;

F2 : la la resultante de la force du vent sur la partie inferieur de la tour ;

ߩߨܸ²‫ܦ‬ଵ‫ܮ‬ଵ
⎧ ‫ܨ‬ଵ =
2
(4.37)
⎨ ߩߨܸ²‫ܦ‬ଶ‫ܮ‬ଶ
⎩‫ܨ‬ଶ = 2

Pour déterminer les tensions T1, T2 il faut résoudre le système d’équations suivant :

෍ ‫ =ܨ‬0, ‫ܨ‬௔ + ‫ܨ‬ଵ + ‫ܨ‬ଶ = ܶଵ sin ߙ + ܶଶ sin ߚ

൞

‫ܮ‬ଵ ‫ܨ‬ଶ‫ܮ‬ଶ
෍ ‫ ܯ‬/௢ = 0 , ‫ܨ‬௔ (‫ܮ‬ଵ + ‫ܮ‬ଶ) + ‫ܨ‬ଵ ൬ + ‫ܮ‬ଶ൰+ = ܶଵ sin ߙ (‫ܮ‬ଵ + ‫ܮ‬ଶ) + ܶଶ sin ߚ ‫ܮ‬ଶ
2 2

59
 CHAPITRE 4 : CONCEPTION MECANIQUE

Les diamètres des câbles seront déterminés par le calcul de la contrainte normal qui doit rester
inférieure à la limite élastique du matériau.

(4.38)
4 ܶ௜ ‫ܵܥ‬
݀௜ ≥ ඨ
ߨߪ௘

4.3.6 Les accouplements

Pour notre conception nous avons suivi la procédure que préconise le fabricant
ROTEX dans son [Link] 4.15

Figure 4.15:Développement d'un accouplement Rotex sur DAO 3D [9]

[Link] Descriptions

Les accouplements Rotex se caractérisent par un encombrement réduit une masse peu
importante et une faible inertie pour une transmission de couples élevés. Performances
techniques et durée de vie sont accrues par un usinage précis de tout l’[Link] assurent
une transmissions des efforts en amortissant les vibrations trigonnelles et absorbent les à-
coups provoqués par les cycles irréguliers des organes moteurs.

[Link] Fonctionnement

Contrairement aux accouplements élastiques ordinaires dont l’élément central est

soumis a la torsion et s’use rapidement, les dents élastiques des accouplements Rotex ne sont

60
 CHAPITRE 4 : CONCEPTION MECANIQUE

soumises qu’à un effort de compression. Ainsi, l’effort supporté par chacune d’elles peut être
plus important. Les élastomères se déforment à l’effort et lors des couples élevés.

L’angle maximal de torsion des accouplements Rotex toute taille confondue est de 5°.

[Link] Détermination de l’accouplement

La détermination des accouplements Rotex se fait d’après la norme DIN 740-2. Les
accouplements Rotex doivent être déterminés de telles sortes qu’en tenant compte des
différents facteurs de services, le couple qu’ils peuvent transmettre ne soit pas dépassé. Il faut
donc procéder à une comparaison des charge de produisant avec les caractéristiques
préconisées [9].

 Charge occasionnée par le couple nominal

Le couple nominal TKN autorisé de l’accouplement doit être en tenant compte de la

température ambiante, au moins aussi grande que le couple nominal de l’installation. Soit

ܶ௄ே ≥ ܶே ∗ ܵ௧ (4.32)
 traversé de la resonnance

Le couple de rotation de pointe TS se produisant lors de la traversé de la resonnace ne doit pas

etre plus grand que le couple de rotation maximal TKmax.

ܶ௄ ௠ ௔௫ ≥ ܶௌܵ௧ (4.39)

 charge causée par les couples alternés impulsionnels

En tenant compte de la temperature ambiante, le couple alterné TKW autorisé ne doit pas
etre depassé par le couple alterné TW pendant la phase en vitesse de regime.

ܶ௄ௐ ≥ ܶௐ ܵ௧ (4.40)
Pour des fréquences de vibrations f > 10 Hz, il est nécessaire de tenir compte de la
puissance à dissiper PW. La puissance à dissiper provient de l’échauffement de l’élastomère
lors de son travail d’amortissement.

61
 CHAPITRE 4 : CONCEPTION MECANIQUE

La puissance d’amortissement PKW autorisé de l’accouplement dépend de la température

ambiante et ne doit pas être dépassée par la puissance d’amortissement nécessaire à
l’application.

ܲ௄ௐ ≥ ܲௐ (4.41)
Tableau 4.1 : valeurs du coefficient de température St [9]

T -30°C à +40°C +60°C +80°C

+30°C
St 1 1.2 1.4 1.8

Tableau 4.2:démarrage Sz [9]

Démarrage/h 100 200 400 800

Sz 1 1.2 1.4 1.6

Tableau 4.3: A-coup SA/SL [9]

SA/SL 1.5 1.8 2.5

A-coup leger A-coup moyen A-coup fort

4.3.7 La génératrice
[Link] Puissance conventionnelle de la machine

ܲே ‫ܭ‬௘ (4.42)
ܲே' =
cos (߮௡ )ߟ௡

Ke : facteur qui tien compte de la chute dans l’enroulement statorique

Ke=0.985-0.004p avec p=2 nombre de paires de pôles machine à 1500 tr/mn

[Link] Diamètre interne du stator

Il doit être déterminé à partir du tabeau 4.16 qui donne le diamètre interne en fonction de la
puissance conventionnelle et du nombre de paires de pôles.

62
 CHAPITRE 4 : CONCEPTION MECANIQUE

[Link] Diamètre externe du stator

Tableau 4.4 : Diamètre externe en fonction du nombre de pole

2p 2 4 6 8 à 10
D ext (1,85 à 1, 95)D1 (1,56 à 1,61)D1 (1,41 à 1,44)D1 (1,34 à 1,39)D1

Le diamètre doit être réduit et choisit selon l’encombrement standard.

Tableau 4.5 : Puissance nominale

Figure 4.16 : Diamètre en fonction de la puissance conventionnelle

63
 CHAPITRE 4 : CONCEPTION MECANIQUE

[Link] Calcul du pas polaire

ߨ ‫ܦ‬ଵ (4.43)
߬௖௠ =
2‫݌‬

De la figure (4.18) on peut déterminer l’induction dans l’entrefer Bs et la charge linéaire du

stator en A/cm

Figure 4.18 : entrefer et charge linéaire en fonction du pas polaire

[Link] Longueur de la machine

6,1 ∗ 10ଵଵܲே' (4.44)

‫ܮ‬௦ =
ߙఋ‫ܭ‬ఉ ‫ܭ‬௕௢௕‫ܣ‬ଵߚఋ݊ଵ‫ܦ‬ଵ²

Pour les machines ordinaires le dimensionnement dépendra De A1, βδ

ߙఋ = 0,715 , ‫ܭ‬ఉ = 1,09 , ‫ܭ‬௕௢௕ = 0,92

4.3.8 le multiplicateur

La fréquence de rotation de la génératrice est de 1500 tr/mn ce

qui est équivalent à 20 fois la vitesse nominal de rotation du
rotor. Pour cela nous disposerons un multiplicateur entre le
rotor et la génératrice. La figure ci-dessous décrit le type de
multiplicateur choisit.
Figure 4.19:multiplicateur à 2
étages [10]

64
 CHAPITRE 4 : CONCEPTION MECANIQUE

[Link] Calcul des diamètres des arbres

Les arbres sont soumis essentiellement en torsion.

la condition de résistance s’écrit :

32 ‫ܥ‬ (4.45)
߬= < ߬௔ௗ௠
ߨ݀ଷ
On en déduit les diamètres Di :

(4 .46)
య32 ‫ܥ‬
‫ܦ‬௜ ≥ ඨ
ߨ߬௔ௗ௠

Ci : couple sur l’arbre i

[Link] Détermination des modules des engrenages

Les valeurs courantes des largeurs de dentures virent entre 10m et 12m.

Figure 4.20: représentation d'une dent[7]

Le calcul de résistance d’une dent est modélisé par le calcul d’une poutre encastré soumise à
un effort F à l’extrémité libre. Ce calcul nous permet d’écrire la formule approximative qui
relie le module à l’effort F (F=2C/D=2C/mZ) du au couple transmis par l’engrenage :

‫ܨ‬ (4.47)
݇‫ ݉ݕ‬ଶ ≥
ߪ௘௣
k= 10 à 12 ; y : facteur de forme dépend du nombre de dent.

En utilisant les résultats de l’étude d’une éolienne de 5 kW [8], le nombre de dents des
différents engrenages permettant un rapport de multiplication de 4x5 est :

Z1 = 44 ; Z2 = 11, Z3 = 85 ; Z4 = 17.

En utilisant ces résultats sachant que m1=m2 et m3=m4 (condition d’engrènement) on obtient :

65
 CHAPITRE 4 : CONCEPTION MECANIQUE

Z1=44, y=0.391, k=11, C1=3184 Nm, m1=m2=5.5, D1=242 mm, D2=60,5 mm, l=60.5 mm

Z3=85, y=0,430, k=11, C2=796 Nm, m3=m4=3 ; D3=255 mm, D4=51 mm, l=33 mm

4.3.9 calcul du frein

Le dispositif de freinage retenu est un frein à disque composé d'un disque (d),
représenté Figure 4.20 et de deux étriers. Le frein est placé à la sortie du multiplicateur ou le
couple est inferieur comparé a l’entré du multiplicateur. Le freinage est réalisé par pression
des garnitures d'usure assurant le serrage de part et d'autre du disque. Chaque étrier supporte 2
garnitures pour lesquelles la surface de contact est représentée Figure 4.21. Le dispositif de
freinage comporte deux étriers.

Figure 4.21 : étrier de frein figure 4.22 : zone de contacte garniture/frein

Il nous faut trouver la relation liant le couple de freinage à la pression que doit exercer l’étrier
sur le disque.

Le couple de freinage doit être supérieur au couple de l’arbre de sortie du multiplicateur.

Sachant que le couple de freinage est donné par la formule suivante :

‫ܦ‬ (4.48)
‫ܥ‬௙௥ = ‫ܨ‬
2
F : force de frottement du a la pression exercée par l’étrier

‫ܵ݌ߤ = ܨ‬௖ (4.49)

µ : coefficient de frottement.

66
 CHAPITRE 4 : CONCEPTION MECANIQUE

4.4 vérification des calculs dynamiques par solidworks

L’interface COSMOS Works nous permettra de valider les dimensions des organes
principaux (pales, arbre, tour) et d’opérer des optimisations ou correction quand c’est
nécessaire.

4.4.1 La tour

Le calcul analytique nous donne les dimensions suivantes :

D1= 90 mm ; e1=18 mm ; h1=4500 mm

D2=110 mm ; e2=22 mm ; h2=5500 mm

Les forces appliquées sont :

Fa=1200 N ; Pn=30 000 N ; F1=450 N ; F2= 671 N

Nous avons effectué deux études statique, l’une sans utiliser les câbles et la seconde avec afin
de justifier son importance.

Figure 4.23:distribution des contraintes (tour simple)

67
CHAPITRE 4 : CONCEPTION MECANIQUE

Figure 4.24: Distribution des contraintes (tour haubanée)

Figure 4.25:Distribution des déplacements (tour simple)

68
 CHAPITRE 4 : CONCEPTION MECANIQUE

Figure 4.26:Distribution des déplacements (tour haubanée)

De la figure 4.23 nous remarquons la présence de zone coloré en orange au niveau de

la partie inferieur de la tour prés du [Link] contrainte avoisine les 135 MPA ce qui dépasse la
contrainte pratique du matériau utilisé 110 MPa. Par contre en utilisant les câbles de soutien,
et d’après les résultats de simulations illustrés figure 4.24, la contrainte maximale que subit la
tour est de l’ordre de 25 Mpa.

Pour les déplacements, nous constatons en analysant les résultats de simulations pour
une tour sans haubanage fig.4.25, que le sommet de la tour se déplace de 500 mm ce qui est
dangereux pour la structure et inacceptable. Pour le cas d’une tour haubanée les déplacements
sont minimes et se situent au niveau des mi-hauteurs des deux parties du mat de l’ordre de 0.5
mm.

4.4.2 L’arbre

Les dimensions de l’arbre sont les suivantes :

D= 35 mm diamètre du noyau ;

L=1500 mm ;

69
 CHAPITRE 4 : CONCEPTION MECANIQUE

Longueur des cannelures l=60 mm ;

Figure 4.27: Distribution des contraintes

Figure 4.28: Distribution des déplacements

L’étude statique indique des contraintes et des déplacements dépassant les limites du
matériau. Pour palier à ce problème, et vus la géométrie non uniforme de l’arbre nous avons
opérer des modifications aux dimensions à savoir :

Augmentation du diamètre du noyau ;

Ajouter des congés au niveau des sections droites ;

70
 CHAPITRE 4 : CONCEPTION MECANIQUE

Les résultats de l’analyse après optimisation sont les suivants :

Figure 4.29: Distribution des contraintes

Figure 4.30:Distribution des déplacements

Les résultats des simulations après modifications sur les dimensions et géométrie de l’arbre
sont satisfaisants. De la figure 4.29 nous déduisons que la contrainte maximale que subit
l’arbre est de l’ordre de 101 MPa (zone colorée en jaune) ce qui est inferieur à la contrainte
pratique du matériau. Nous remarquons de la figure 4.30 des déplacements de l’ordre du 1/10
mm.

4.4.3 La pale

Les dimensions des éléments de la pale sont obtenues par le calcul aérodynamique. Le calcul
dynamique permet la détermination de la variation de l’épaisseur en fonction du rayon.

71
 CHAPITRE 4 : CONCEPTION MECANIQUE

La pale contenant 20 éléments de longueur égale nous avons opté pour une simulation
élément par élément sachant que la force appliqué en une section i est la somme des forces
appliquées aux (i-1) forces précédentes en partant du bout de la pale.

Figure 4.31: Chargement du premier élément de pale

72
CHAPITRE 4 : CONCEPTION MECANIQUE

Figure 4.32 : Distribution de déplacements sur l’élément en bout de pale

Figure 4.33: Distribution des contraintes sur l’élément en bout de pale

73
 CHAPITRE 4 : CONCEPTION MECANIQUE

Figure 4.34 : distribution des déplacements sur le dernier élément de pale

Les résultats des simulations indiquent que la pale n’est pas soumise à des contraintes
dépassants la limite du matériau. Les déplacements relative de chaque éliment par rapport à
l’élément adjacent sont insignifiants de l’ordre de 2/100 mm. Le déplacement en bout de pale
est la somme des déplacements relatifs qui avoisine le ½ mm

74
 CONCLUSION

CONCLUSION

Le but du présent projet était le dimensionnement et la conception d’une turbine

éolienne de faible puissance en utilisant un outil informatique, à savoir le logiciel de
conception SolidWorks. Etant donné que le sujet traitant de la protection de l’environnement
et l’utilisation de l’énergie renouvelable est d’actualité dans le monde, nous avons dirigé notre
choix de mécanisme dans ce sens et choisis de concevoir une éolienne à axe horizontale de
faible puissance (25 kW) pour être installé sur le toit du futur bâtiment vert de l’ENSP.

L’étude aérodynamique nous a permis de déterminer le diamètre du rotor permettant la

production de l’énergie électrique requise, ainsi que les dimensions de la géométrie de la pale.
Nous avons pu établir une courbe de puissance en fonction de la vitesse du vent dont nous
avons déduit la vitesse minimale (V=6 m/s) ou l’éolienne commence à être productive.

A l’aide des calculs de résistance et l’étude technologique nous avons pu déterminer le

mode d’assemblage entre les organes constitutifs, le type de matériau, la géométrie permettant
l’économie de matière et un fonctionnement sécurisé. Nous avons aussi cherché à faciliter
l’opération de maintenance et de démontage-montage des différents organes. Les simulations
par le logiciel nous ont permis d’effectuer les corrections nécessaires du fait des hypothèses
simplificatrices utilisées lors du calcul analytique. Nous citons à titre d’exemple la tour.
L’étude effectuée par COSMOSWorks sans utilisations de câbles de soutien indique un
déplacement important au niveau de la nacelle ce qui est dangereux pour la structure. Par
contre en utilisant la même étude avec les câbles de soutien, la structure ne subit que des
déplacements minimes de l’ordre du 0.5 mm.

En résumé nous avons pu effectuer une ébauche de conception d’une éolienne à axe
horizontale. En effet cette conception peut être perfectionnée. Nous proposons une étude
aérodynamique plus poussée notamment en utilisant l’interface FLOWWorks de solidWorks
qui peut remplacer les essais en souffleries. Il est possible d’augmenter les performances
aérodynamiques de l’éolienne en utilisant une régulation du pas des pales (pitch control).
Concernant le coté mécanique, un choix différent d’organes, des simulations de vibrations,
étude de flambage et fatigue, résistance des assemblages permettraient d’affiner la conception.

76
Bibliographie

[1] B. MULTON, [Link], [Link] AHMED, [Link], [Link], [Link],

[Link] ÉTAT DE L’ART DANS LES AÉROGÉNÉRATEURS ÉLECTRIQUES-Extrait
du rapport de synthèse ECRIN « L’Electronique de Puissance Vecteur d’Optimisation Pour
les Energies Renouvelables », paru mai en 2002
[2] [Link], [Link], [Link], Éoliennes – technique de l’ingénieur-BM
4 640
[3] [Link]- Développement d’un modèle représentatif d’une éolienne afin d’´etudier
l’implantation de plusieurs machines sur un parc éolien. Thèse pour obtenir le grade de
Docteur De l’Ecole Nationale Supérieure d’Arts et Métiers. 20 décembre 2007
[4] J.F MANWELL, J.G McGOAN, A.L ROGERS wind enrergy explained-theory,design
and application, Amherts, USA
[5] [Link]- L’aérodynamique de la turbine éolienne et ses composantes mécaniques-
laboratoire de recherche en énergie éolienne LREE-université du Quebec à Rimouski-22-25
octobre 2007
[6] R. R. Ramsay, [Link], [Link]- Wind Tunnel Testing of Three S809
Aileron Configurations for use on Horizontal Axis Wind Turbines- The Ohio State University
Columbus, Ohio- July 1996
[7] [Link]- guide du calcul en mécanique- HACHETTE Technique
[8] le tube d’acier-profil de construction -edition 1989- chambre syndicale des fabricants de
tubes en fer et en acier
[9] [Link]
[10] M. Khechana*, M. Nadjah, L. Laiche†, T. Ouksel et C. Mahfoudi- Etude de la nacelle et
du mat d’une éolienne de faible puissance- Revue des Energies Renouvelables CISM’08 Oum
El Bouaghi (2008) 195 – 203
[11] N. Jovičić, D. Gordić, [Link]ć, V. Šušteršič- NUMERI^KI ALGORITAM ZA
PRORA^UN AEROTURBINA MALE SNAG- Univerzitet u Kragujevcu
[12] [Link]- guide du dessinateur industriel-HACHETTE Technique édition 2007
[13] M.A. HAMDAOUI-Conception et Dimensionnent D’une éolienne tripale-École
Nationale Supérieure d’Arts et Métiers MEKNES-2007/2008
[14] [Link], [Link]- Résistance des matériaux-technique de l’ingénieur-
Form.C 2 060

77
[15] [Link]- ÉTUDE DU COMPORTEMENT DES FIXATIONS PAR VIS DES
COURONNES DE GUIDAGE DE GRAND DIAMETRE : ELABORATION D’UN OUTIL DE
DIMENSIONNEMENT- INSAT-05 février 2008
[16] M. Khechana*, M. Nadjah, L. Laiche†, T. Ouksel et C. Mahfoudi- Etude de l’hélice
d’une éolienne de 5 kW- Revue des Energies Renouvelables CISM’08 Oum El Bouaghi
(2008) 195 – 203
[17] Marcel BOYER- assemblage par éléments filetés- technique de l‘ingénieur- B 5 565
[18] R. R. Ramsay, [Link], [Link]- Wind Tunnel Testing of Three S809
Aileron Configurations for use on Horizontal Axis Wind Turbines- The Ohio State University
Columbus, Ohio- July 1996

78
Annexes
A : données et résultats:

A.1 caractéristiques du profil

Tableau 1:coordonnées du profil S809

79
Tableau A.2:variation des coefficients de portance et trainé pour le profil S809

α Cl Cd
-180.00 . 000 .1748
-170.00 . 230 0.2116
-160.00 . 460 0.3172
-150.00 . 494 0.4784
-140.00 . 510 0.6743
-130.00 . 486 0.8799
-120.00 . 415 1.0684
-110.00 . 302 1.2148
-100.00 . 159 1.2989
-90.00 . 000 1.308
-80.00 -. 159 1.2989
-70.00 -. 302 1.2148
-60.00 -. 415 1.0684
-50.00 -. 486 0.8799
-40.00 -. 510 0.6743
-30.00 -. 494 0.4784
-20.10 -. 560 0.3027
-18.10 -. 670 0.3069
-16.10 -. 790 0.1928
-14.20 -. 840 0.0898
-12.20 -. 700 0.0553
-10.10 -. 630 0.039
-8.20 -. 560 0.0233
-6.10 -. 640 0.0131
-4.10 -. 420 0.0134
-2.10 -. 210 0.0119
. 10 . 050 0.0122
2.00 . 300 0.0116
4.10 . 540 0.0144
6.20 . 790 0.0146
8.10 . 900 0.0162
10.20 . 930 0.0274
11.30 . 920 0.0303
12.10 . 950 0.0369
13.20 . 990 0.0509
14.20 1.010 0.0648
15.30 1.020 0.0776

80
 16.30 1.000 0.0917
17.10 . 940 0.0994
18.10 . 850 0.2306
19.10 . 700 0.3142
20.10 . 660 0.3186
30.00 . 705 0.4784
40.00 . 729 0.6743
50.00 . 694 0.8799
60.00 . 593 1.0684
70.00 . 432 1.2148
80.00 . 227 1.2989
90.00 . 000 1.308
100.00 -. 159 1.2989
110.00 -. 302 1.2148
120.00 -. 415 1.0684
130.00 -. 486 0.8799
140.00 -. 510 0.6743
150.00 -. 494 0.4784
160.00 -. 460 0.3172
170.00 -. 230 0.2116
180.00 . 000 0.1748

A.2 force aérodynamiques sur une pale :

Tableau 3 : distribution de la corde et des forces aérodynamiques suivant le rayon

r/R C F f r/R C F f
0,165 0,530 21,427 11,054 0,605 0,195 8,762 3,453
0,209 0,466 23,370 20,398 0,649 0,182 3,801 1,320
0,253 0,411 22,139 56,734 0,692 0,171 2,871 7,811
0,297 0,364 22,272 18,341 0,736 0,162 8,818 2,369
0,341 0,326 28,157 22,231 0,780 0,153 15,566 5,108
0,385 0,294 27,654 20,230 0,824 0,145 22,757 7,770
0,429 0,267 24,305 16,276 0,868 0,138 30,417 10,370
0,473 0,245 21,025 12,911 0,912 0,131 38,535 38,731
0,517 0,226 17,563 9,819 0,956 0,126 47,175 46,258
0,561 0,209 13,403 6,611 1,000 0,120 56,278 17,880

81
B : accouplement Rotex
Tableau B.1 : données techniques

Figure B.1 représentation des dimensions de l’accouplement

82