LA COMMANDE INDIVIDUELLE DU CALAGE DES PALES
LA MODรLISATION AรROรLASTIQUE
L’environnement donne lieu aux oscillations de รฎa structure d’une รฉoรฎienne ร cause des variations du vent, de l’action des vagues et des effets de la gravitรฉ. Ce chapitre explique la physique des รฉoliennes dans leur environnement et les modรจles mathรฉmatiques qui reprรฉsentent les interactions de la structure, l’aรฉrodynamique et le systรจme de contrรดle, dont l’ensemble est connu sous le nom de la dynamique aรฉroรฉlastique. On ne tient compte ni de la dynamique des interactions hydrodynamiques, ni des effets dynamiques dans l’รฉvolution du couple de la gรฉnรฉratrice รฉlectrique. Le chapitre commence en introduisant la dynamique aรฉroรฉlastique qui s’applique aux basses frรฉquences, puis dรฉcrit les modรจles mathรฉmatiques dans FAST et HAWC2. Le NREL Offshore Baseline Turbine^’85, d’une puissance nominale de 5MW, sert comme cas d’รฉtude dans ce travail. Son rotor a un diamรจtre de 126 m, la hauteur de sa tour est de 90 m et elle est contrรดlรฉe en vitesse variable et en calage. En principe, รฎe modรจle reprรฉsente une machine avec les propriรฉtรฉs semblables ร la machine REpower 5M.
La physique des รฉoliennes
Les รฉoรฎiennes transforment l’รฉnergie cinรฉtique du vent en รฉnergie รฉlectrique. Cette section rรฉsume les รฉtapes du processus, qui est une dynamique continue entre l’aรฉrodynamique, la structure et le systรจme de contrรดle.
L’aรฉrodynamique
Les pales des รฉoliennes sont des profils aรฉrodynamiques tournants conรงues pour produire de รฎa portรฉe et donc du couple, ร partir de l’รฉcoulement de Fair qui les entoure. Dans la perspective d’une moyenne temporelle, le rotor agit comme un disque sustentateur, qui donne lieu ร une chute de pression ร travers le rotor et ร un ombrage qui traรฎne dans la direction de l’รฉcoulement. ร cause de l’extraction de l’รฉnergie, l’รฉcoulement est ralenti par la prรฉsence de l’รฉolienne et, puisque les forces doivent รชtre finies, ce ralentissement est distribuรฉ dans la direction principale de l’รฉcoulement. La forme des profils produit une circulation autour des pales et cette circula tion gรฉnรจre un tourbillon continu au bout et ร la racine de chaque pale. Ces tourbillons traรฎnent dans l’ombrage et ont pour effet de diminuer la portรฉe produite aux extrรฉmitรฉs des pales.
L’induction dynamique
On appelle induction le ralentissement de l’รฉcoulement qui se produit par la perte de l’รฉnergie dans l’รฉcoulement. Quand la vitesse moyenne du vent ou les angles de calage des pales changent dans le temps, l’รฉtablissement de l’ombrage et donc l’induction, crรฉรฉ une dynamique. On peut comprendre cette dynamique comme intermรฉdiaire entre deux cas idรฉaux : le premier est celui oรน la vitesse du vent dans le plan du rotor prend immรฉdiatement la valeur qu’elle atteindrait en รฉquilibre (l’ombrage รฉquilibrรฉ) et le second est celui oรน le rotor n’a pas d’effet sur la perturbation de la vitesse du vent (l’ombrage fixe). Dans les cas oรน l’induction est รฉlevรฉe, normalement aux basses vitesses du vent, la situation tend plus vers l’ombrage fixe, puisque les vitesses du vent dans l’ombrage ne sont pas assez รฉlevรฉes pour atteindre un รฉquilibre pendant les variations de l’รฉcoulement. Quand l’induction est moins รฉlevรฉe, ce qui arrive quand l’angle de calage et donc la vitesse du vent, est plus รฉlevรฉe, l’รฉquilibre s’รฉtablit plus rapidement et la situation tend plutรดt vers l’ombrage รฉquilibrรฉ.
L’aรฉrodynamique transitoire
Les effets de l’induction dynamique sont importants dans un temps caractรฉristique de l’ordre du rayon du rotor divisรฉ par la vitesse nominale du vent, tandis que souvent, les pales vivent des conditions beaucoup plus transitoires, celles-ci ayant un temps caractรฉristique de l’ordre de la corde du profil divisรฉe par la vitesse locale du vent relative ร la pale. Normalement, cette derniรจre durรฉe de temps est d’un ordre plus rapide que le temps caractรฉristique d’intรฉrรชt dans l’induction dynamique^101. Afin de mieux comprendre cette aรฉrodynamique transitoire, on pose les variations dans l’รฉcoulement comme les dรฉplacements du profil autour de l’axe qui longe la pale, dans l’axe parallรจle ร la corde du profil et dans l’axe perpendiculaire ร la corde (voir la figure 1.2). Si l’รฉcoulement reste attachรฉ au profil, la circulation s’ajuste constamment aux changements dans l’angle d’attaque, ce qui donne lieu ร un tourbillon qui se dรฉclenche de l’arriรจre du profil ร cause des forces produites par l’accรฉlรฉration de l’air autour du profil.
Les dรฉplacements du profil sont aussi liรฉs aune dynamique dans le processus de la sรฉparation de l’รฉcoulement du profil. Les dรฉplacements du profil influencent la distribution de la pression autour de celui-ci et c’est justement cette distribution qui donne lieu et qui maintient la sรฉparation de l’รฉcoulement. Les profils des pales des รฉoliennes sont relativement รฉpais (>15%), donc la sรฉparation ne commence guรจre au bord d’attaque, mais plutรดt de l’arriรจre du profil. Ce dernier type de sรฉparation commence par une migration du point arrรชt sur la face arriรจre du pro fil (dessus) ร cause d’une basse pression vers l’arriรจre du profil. L’รฉcoulement se sรฉpare quand un tourbillon se forme au bord de fuite et couvre le dessus du profil. Ce tourbillon est transportรฉ lentement dans le sens de l’รฉcoulement et part du profil quand l’angle d’attaque prend une cible critique (figure 1.3). Pour les dรฉplacements angulaires dans l’axe qui longe la pale et ceux dans l’axe perpendiculaire ร la corde du profil, ce n’est que la rรฉgion de basse pression qui cause la sรฉparation ; dans le cas des dรฉplacements dans l’axe de la corde du profil, il y a aussi un ralentissement relatif de l’รฉcoulement qui encourage la migration du point arrรชt.
ย La dynamique structurelle
Les charges aรฉrodynamiques, gravitationnelles et inertielles sont la source de la plupart des vibrations dans les composantes flexibles des รฉoliennes. Les frรฉquences naturelles de ces vibrations dans les รฉoliennes multi-mรฉgawatt, telles que la machine NREL Offshore Baseline Turbine, sont de l’ordre de 1,5 ร 10 fois la frรฉquence de rotation du rotor pour les premiers modes de flexion et de l’ordre de 10 ร 15 fois la frรฉquence de rotation pour les deuxiรจme modes de flexion (cf. Rรฉf. 8). La plupart des modes ont un faible amortissement de la part de la structure, tandis que l’amortissement dominant provient des effets de l’aรฉrodynamique. En pratique, il est commun d’augmenter l’amortissement de ia tour au moyen d’un amortisseur liquide dans le haut de la tour’-12′ 13-1. Les charges de ia gravitรฉ dominent les charges parallรจles au plan du rotor et donc elles sont normalement associรฉes aux dรฉflexions des pales dans ie plan du rotor et aux charges dans รฎ’axe de rotation des enroulements de calage. La gravitรฉ gรฉnรจre aussi les couples qui tendent
ร fausser la pale quand celle-ci est flรฉchie hors le plan du rotor. Les charges qui surgissent ร cause de l’inertie, par contre, sont produites par les accรฉlรฉrations relatives des composantes de la structure. Durant l’opรฉration de ia machine, ces charges font ie lien entre les vibrations des diverses composantes de la structure et dรฉterminent donc les couplages entre les composantes dans la dynamique de la machine. L’inertie s’avรจre trรจs importante aussi dans les รฉvรฉnements extrรชmes, tels que les arrรชts soudains, oรน les accรฉlรฉrations sont sรฉvรจres.
|
Table des matiรจres
RรSUMร
LISTE DES FIGURES
LISTE DES TABLEAUX
LISTE DES ABRรVIATIONS
LISTE DES SYMBOLES
INTRODUCTION
CHAPITRE I LA MODรLISATION AรROรLASTIQUE
1.1 La physique des รฉoliennes
1.1.1 L’aรฉrodynamique
1.1.2 La dynamique structurelle
1.1.3 La commande
1.2 Les propriรฉtรฉs de la machine NREL Offshore Baseline Turbine
1.3 La simulation del’aรฉroรฉlasticitรฉ
1.3.1 HAWC2
1.3.2 FAST
1.3.3 Les modรจles linรฉaires
1.4 Les conclusions
CHAPITRE II LE CHOIX DES CAPTEURS
2.1 L’analyse de l’observabilitรฉ
2.1.1 La distance ร l’inobservabilitรฉ : definition et calcul
2.1.2 Le conditionnement de l’estimateur d’รฉtats
2.1.3 Les effets de la dimensionnalitรฉ dans la comparaison des sorties
2.1.4 Le gramian d’observabilitรฉ
2.1.5 Les aspects numรฉriques du choix des capteurs
2.2 La conception des rรฉseaux de capteurs
2.2.1 Les stratรฉgies de recherche
2.3 Les conclusions
CHAPITRE III L’ESTIMATION DES รTATS
3.1 La thรฉorie de Testimation des รฉtats
3.1.1 L’estimation par la mรฉthode des moindres carrรฉs
3.1.2 Le filtre de Kalman
3.1.3 L’estimation ร rรฉponse impulsionnelle finie
3.1.4 L’observateur optimal LTP
3.1.5 L’estimation des modรจles incertains
3.2 L’estimation des รฉtats dans l’environnement HAWC2
3.2.1 L’initialisation
3.2.2 Les รฉvรฉnements en temps simulation
3.3 Les rรฉsultats 56
3.3.รฎ L’optimisation des paramรจtres de l’estimateur
3.3.2 Les รฉcoulements dรฉterministes
3.3.3 Les รฉcoulements turbulents
3.3.4 L’analyse des suppositions
3.4 Les conclusions
CHAPITRE IV LA COMMANDE INDIVIDUELLE DU CALAGE DES PALES
4.1 La rรฉgulation linรฉaire et quadratique avec mitigation des perturbations
4.2 La rรฉgulation de l’angle d’attaque
4.3 Les rรฉsultats
4.4 Les conclusions
CONCLUSION
RรFรRENCES
ANNEXEA ANALYSE DES SYSTรMES LPV
A.l La thรฉorie FIoquet-Ljapunov
ANNEXEB LES SORTIES FAST
ANNEXEC
LA DESCRIPTION DES CODES
Cl La modรฉlisation aรฉroรฉlastique
C.2 L’analyse des modรจles linรฉaires
C.3 L’estimation
C.4 La commande
ANNEXED
LES FICHIERS D’ENTRรE FAST ET HAWC2
D.l Les fichiers FAST
D.l.l Lefichierist
D.l.2 Les donnรฉes structurelles de la tour
D.1.3 Les donnรฉes structurelles des pales
D.1.4 Les commandes aรฉrodynamique
D.1.5 Les entrรฉes du vent
D.1.6 Les commandes de linรฉarisation
D.2 Les fichiers d’entrรฉe HAWC2
Tรฉlรฉcharger le rapport complet