Consideration de L Incertitude A L Etude
Consideration de L Incertitude A L Etude
Consideration de L Incertitude A L Etude
Mouna Hadj Kacem1,2, Abdelkhalak El Hami1, Hassen Trabelsi2, Lassaad Walha2, Mohamed Haddar2
1
Laboratoire de Mécanique de Normandie LMN, Institut Nationale des Sciences Appliquées de Rouen, France
Laboratoire de Mécanique, Modélisation et Productique LA2MP, Ecole National d’Ingénieurs de Sfax, Tunisie
2
RÉSUMÉ. La présente étude porte sur le rendement de transmission dynamique d’un multiplicateur d’éolienne double
étage en considération l’incertitude. Le travail consiste à considérer l’incertitude relative à la détermination du
comportement dynamique de vitesse de déplacement angulaire d’entrée et de la sortie pour aboutir à l’efficacité. Les
incertitudes envisagées dans ce travail sont le module d’engrenage et la rigidité à la torsion, appliquées indépendamment
pour comprendre l’effet de chacun sur le système, et finalement d’une façon couplée pour atteindre des résultats par
interaction.
Afin de s’approcher des valeurs réelles expérimentales, on applique l’incertitude avec la méthode du Monte Carlo dont les
résultats seront la référence de vérification pour l’approche probabiliste Polynôme de Chaos.
ABSTRACT. The present study investigates the dynamic transmission efficiency of a two-stage wind turbine gearbox by
considering the uncertainty. The work consists of considering the uncertainty in determining the dynamic behavior of input
and output angular displacement velocity to achieve the efficiency. The uncertainties considered in this work are the gear
modulus and torsional stiffness, applied independently to understand the effect of each on the system and finally in a
coupled manner to reach results by interaction.
In order to approach the real experimental values, the uncertainty is applied with the Monte Carlo method, whose results
will be the verification reference for the probabilistic Polynomial Chaos approach.
MOTS-CLÉS. Eolienne, Rendement d’engrenage, Analyse d’incertitude, Polynôme de Chaos.
KEYWORDS. Wind Turbine, Gearbox efficiency, Uncertainty analysis, Polynomial Chaos.
1. Introduction
Faisant suite aux études de stabilité dynamique linéaire et non-linéaire [FAK 05], [ABB 11] les
engrenages se présentent comme étant source de bruit résultant des vibrations qui découlent de la
variété des variables incertaines [ELH 13]. Par conséquent, afin de garantir la fiabilité d’un système, il
est manifestement essentiel d’étudier la fiabilité de l’efficacité de ces sous-systèmes.
Dans la littérature, de nombreux outils d’étude de fiabilité sont utilisés, dont la méthode classique
de Monte Carlo (MC), qui est la plus idéale pour tout type de simulation, toutefois plus coûteuse en
termes de temps et du budget [KAL 09]. Par conséquent, diverses méthodes sont utilisées en fonction
du type de simulation, du nombre de variables et de la complexité du calcul. Citons la méthode de
Polynôme de Chaos Généralisé (GPC) [XIU 04] déjà utilisée à plusieurs reprises pour l’étude du
comportement dynamique des systèmes d’engrenages en tenant compte des incertitudes [GUE 15],
[BEY 16].
Le système de conversion d’énergie (SCE) d’une éolienne est un système complexe constitué de
sous-systèmes sujets aux vibrations. Ces dernières résultent de la variation de vitesse du vent au niveau
du couplage fort avec le rotor, ces vibrations sont transmises au multiplicateur, lui-même source
d’autre sortes de vibrations au niveau de la denture d’engrenages, raideur d’engrènement ; pour aboutir
au système le plus sensible dans le SCE, la génératrice. Le fonctionnement de ce troisième sous-
système dépend de la vitesse angulaire lui étant transmise, aussi de la stabilité des performances
dynamique du deuxième sous-système. Afin d’étudier le bon fonctionnement du SCE, il est question
d’étudier le rendement du multiplicateur dont l’objectif de ce papier. Il est présenté par un système
d’engrenages double étage (SEDE) avec un modèle à 12 degrés de liberté soumit à la variation du
couple aérodynamique en fonction de temps et à la fluctuation de la rigidité de l’engrènement.
Le model dynamique du SEDE compte 12 degrés de liberté répartis comme suit : 6 déplacements
linéaires partagés en deux (Xj, Yj) sur chaque bloc (j=1 à 3), tels qu’illustrés dans la figure.1, et 6
déplacements rotationnelles au sein de chaque roue dentée ainsi qu’au niveau du rotor et d’éléments
tournants de la génératrice, conformément au vecteur de coordonnées générales Q(t ) :
Q(t ) x1, y1, x 2, y 2, x 3, y 3, (1,1), (1,2), (2,1), (2,2), (3,2), (3,3)
T
[1]
Comme la montre la figure.1, chaque bloc est maintenus par un palier flexible qui se caractérise par
une rigidité à la flexion K jx et une rigidité à la traction-compression K jy . Nous admettons une certaine
rigidité à la torsion K j aux arbres intermédiaires, étant donné leur masse insignifiante par rapport aux
autres sous-systèmes du SCE.
n t L Q t
n
[2]
Avec L n
dépond de n-contact :
sin 1 ,cos 1 ,sin 1 , cos 1 ,0,0,0, r 12, r 21,0,0,0
b b
L 1 [3]
M Q Cs
L C t Q Ks L
L K t Q Fext
T T
L [5]
La matrice de masse générale est caractérisée par la matrice de masse M L et celle des
inerties M A :
ML 0
M [6]
0 MA
La matrice de rigidité du structure Ks est une matrice constante qui dépende de la rigidité à la
flexion et la rigidité à la traction-compression Kp avec la matrice de rigidité des arbres à la
torsion K :
KP 0
KS [7]
0 K
K t K c K v t [8]
C M Ks K c [9]
Le vecteur des forces extérieures repose sur le couple aérodynamique et le couple du générateur,
variables en fonction du temps :
1 1
Caero t . . . air .R 2 .C p .Vwin
3
t [11]
2 wt
Avec 1 ; 2 4
Caero t
Cgen t [13]
gear
Z12 .Z 22
Avec gear [14]
Z 21.Z 32
Reliant les intrants de la génératrice aux extrants du rotor, ce rapport est exploité cinématiquement,
après avoir envisagé, comme un paramètre à valeur constante, suite à la caractérisation de système à
corps rigide. Mais dans le contexte de notre étude, le système inclut la flexibilité des composants tels
que la flexibilité des paliers, la flexibilité des arbres reliant les sous-systèmes et celle des dentures
d’engrenages. Le rapport adopte une évolution au cours du temps plus concrète fluctuant autour d’une
valeur moyenne qui décrit le mouvement du corps rigide.
32 t gener
Rap t [15]
11 t rotor
L’étude de l’incertitude est aboutie principalement par le Polynôme de Chao Généralisé (GPC) et les
résultats sont validés par la comparaison de ces derniers par rapport à celles achevée par la méthode de
Monte Carlo (MC).
La méthode de MC se traduit par le fait de solver simultanément un modèle déterministe incorporant
des facteurs incertains modélisés par des variables aléatoires [MAK 16]. Cette technique est applicable
à toute sorte de fonction indépendamment de leurs entités probabilistes et des corrélations
correspondantes de leurs distributions de probabilité [KAL 09].
La démarche du Chao polynomial [CHO 16] consiste à projeter des quantités aléatoires dans un
espace de probabilités sous forme de présentation spectrale, en partant de la théorie de Wiener [ZEN
19] relative au chaos homogène. Elle introduit une distinction entre les éléments stochastiques et
déterministes d’une fonction aléatoire [DAM 19].
Toute fonction du second ordre H dans l’espace d’Hilbert L2 , , Pr peut être élaborée comme
une série de fonctions polynomiales convergentes conformément à la logique de Wiener et la notion
généralisée de Cameron Martin [LUO 18] :
U v j j ; ___ 1 , 2 ,... [16]
j 0
U v j j ; [17]
j 0
Avec L0 1; __ L1
4. Simulations et résultat
La simulation numérique réalisée sous MATLAB consiste à déterminer les vecteurs de vitesses de
déplacement des facteurs de liberté angulaire, pour achever l’analyse du RTD. Par la suite, on
procédera au recalcule de ces variables en appliquant la méthode probabiliste, en recourant à la loi
uniforme du polynôme du chaos, en fonction de l’écart-type des paramètres incertains introduits.
Les perturbations de vitesses de déplacement pour les premières fractions de seconde dans la
figure.2 sont plus significatives sur les courbes des engrenages entraînés 21 et 32 , cependant
l’influence de la flexibilité de la structure du système peut également être observée sur le profil de la
courbe finale 33 comparativement à la courbe initiale 11 . Ces irrégularités de vitesses de déplacement
engendrent des perturbations sur le RTD de l’engrenage comme illustré dans la figure.3.
Le présent rapport fonde les simulations d’incertitude par la mise en œuvre du GPC sur le modèle
dynamique approprié et les vérifier avec 1000 simulations de MC. Ces résultats visent la vitesse de
déplacement du rotor 11 et la vitesse de déplacement des éléments rotatifs du générateur 33 , pour
aboutir au comportement du RTD. On considère l’incertitude pour le module (m=16 10-3[m]), qui
décrit généralement les défauts de fabrication [TRA 21], la rigidité à la torsion ( K 5.105 [N.m/rad]).
La figure.4 décrit l’influence de l’incertitude du module par 5 courbes sur chaque figure (a), (b) et
(c) dont on peut trouver la courbe déterministe, les courbes moyennes de GPC et MC et les courbes des
valeurs extrêmes de MC. C’est bien clair que pour une valeur minimale de module la réponse présente
plus de perturbation sur la vitesse de déplacement ainsi que sur l’RTD
L’étude préalable prouve l’influence de chaque paramètre sur le RTD, ce qui entraîne l’étude du
choix de valeurs de ces paramètres avec une considération de l’incertitude sous l’effet de l’interaction
entre eux.
La concordance des réponses GPC et MC est assurée à la suite du recours à un polynôme d’ordre 4
correspondant à une combinaison d’écart-type de m 3%; K 4% . Les illustrations de la
figure.6(a), (b) et (c) décrivent l’influence de l’interaction entre ces deux paramètres au sein de la
considération de l’incertitude associée à la détermination de la vitesse de déplacement d’entrée, de
sortie et du RTD.
La figure.6 met en évidence tout d’abord la correspondance des courbes des valeurs moyennes GPC
et MC par rapport à la réponse déterministe ainsi que 4 autres courbes. Celles-ci présentent la réponse
des valeurs extrêmes du nuage de variables aléatoires de la simulation MC qu’on retrouve dans le
tableau.1 les valeurs de cette simulation {min(m)-min(Kθ) ; min(m)-Max(Kθ) ; Max(m)-min(Kθ) ;
© 2022 ISTE OpenScience – Published by ISTE Ltd. London, UK – openscience.fr Page | 9
Max(m)-Max(Kθ)}, afin d’avoir une évaluation concernant le type de réponse de ces deux variables en
interaction.
Les deux courbes les plus significatives sont initialement, celle du module minimal et du maximum
de rigidité. Cette courbe procure une avance de phase par rapport aux autres courbes. Ensuite, la
courbe du quatrième cas est également importante et prouve une diminution de l’amplitude de vitesse
de déplacement en fonction du le temps, comme illustré à la figure.6(c).
L’évaluation des vibrations des systèmes rotatifs comme les engrenages est procédée,
principalement, avec l’étude de déplacement de chaque élément selon les degrés de liberté associé à
chaque élément. De la même importance la vitesse de ces déplacements est aussi indispensable pour
évaluer son efficacité.
5. Conclusion
Le travail commence par saisir la sensibilisation du système pour chaque variable, indépendamment
à la détermination de la vitesse de déplacement dynamique du rotor de l’éolienne et celui des éléments
tournants du générateur pour aboutir à la sensibilisation du RTD de SEDE. Ensuite, le même travail se
traduit avec la considération de l’interaction entre les deux variables.
Les résultats de la sensibilisation prouvent que le RTD peut être plus stable pour une valeur
maximale de la rigidité à la torsion pendant la considération de l’incertitude par interaction avec le
module.
Bibliographie
[FAK 05] FAKHFAKH, TAHAR, ET AL. "Effect of manufacturing and assembly defects on two-stage gear system
vibration." The international journal of advanced manufacturing technology (2005).
[ABB 11] ABBOUDI, KAMEL, ET AL. "Dynamic behavior of a two-stage gear train used in a fixed-speed wind
turbine." Mechanism and Machine theory 46.12 (2011): 1888-1900.