ROYAUME DU MAROC المملكة المغربية
Ministère de l’Education Nationale, de la
Formation Professionnelle, de l’Enseignement وزارة التربية الوطنية والتكوين المهني
Supérieur et de la Recherche Scientifique والتعليم العالي والبحث العلمي
Université Sidi Mohamed Ben Abdellah جامعة سيدي محمد بن عبد للا
Faculté des Sciences Dhar El كلية العلوم ظهر المهراز فاس
Mahraz - Fès
Année d’étude : Master Année Universitaire : 2023/2024
Module : système de conversion d’énergie éolienne
TP N° 1 : Modélisation, Simulation, et optimisation de la turbine éolienne
Objectifs : analyse des performances d'une éolienne, y compris la création
d'algorithmes, la simulation dans Simulink MATLAB, et l'ajustement de modèles
mathématiques pour représenter le coefficient de puissance Cq.
Aérodynamique
Soient les fonctions Cq et Cp donnée comme suite :
1 𝑐4
𝐶𝑞 = 𝑐1 (1 + 𝑐2 (𝛽 + 𝑐3 )2 ) + (𝑐5 𝜆𝑖 − 𝑐6 𝛽 − 𝑐7 𝛽 𝑐8 − 𝑐9 )𝑒 (−𝑐10𝜆𝑖)
𝜆
1 0.035
avec 𝜆𝑖 = − et 𝐶𝑝 = 𝜆 𝐶𝑞 .
𝜆+0.08∗𝛽 𝑐11 +𝑐12 𝛽 3
Données: c1 = 0.005, c2 = 1.53, c3 = 0.5, c4 = 0.18, c5 = 121, c6 =
27.9, c7 = 198, c8 = 2.36, c9 = 5.74, c10 = 11.35, c11 = 16.1, c12 = 201
1/ Écrire un algorithme permettant de tracer Cq et Cp
2/ Déduire, λopt et Cpmax
3/ Écrire un algorithme permettant de tracer Pv=f(v) et Pt=f(v) pour différentes
vitesses de vent
1 1
𝜌 𝑆 𝑣3,
𝑃𝑣 = 𝑃𝑡 = 𝜌 𝑆 𝑣 3 𝐶𝑝
2 2
Données : 𝑅 = 61, 𝜌 = 1.225.
4/ Ecrire un algorithme permettant de tracer la courbe des puissances optimales
Popt=f(v)
5/ Utilisez l'outil Simulink MATLAB pour retrouver les résultats des calculs de Cq et
Cp que vous avez effectués dans la question n°1.
6/ À partir des résultats de la question n°5, créez un tableau permettant de
représenter la formule théorique de Cq.
7/ À partir du tableau obtenu, pour β=0, trouver les paramètres de la fonction
suivante capable de représenter la fonction Cq.
𝑐2 𝑐
(− 4 )
𝐶𝑞 = 𝑐1 ( − 𝑐3 ) 𝑒 𝜆 + 𝑐5 𝜆
𝜆
8/ Comparer le résultat obtenu avec celui donné dans la question n°1.
Système de calage des pales (pitch system)
1/ système hydraulique
Soit le système d’orientation des pales hydraulique suivant :
𝛽(𝑠) 𝜔𝑛2
= 2
𝛽𝑐 (𝑠) 𝑠 + 2 ∙ 𝜁 ∙ 𝜔𝑛 + 𝜔𝑛2
Ou 𝛽 est l’angle de calage, 𝛽𝑐 est l’angle désirée, 𝜔𝑛 est la pulsation propre du
système, et 𝜁 est le coefficient d’amortissement.
1/ Tracer le schéma équivalent du système sous MATLAB Simulink.
2/ Simuler la réponse indicielle et fréquentielle (boucle ouverte) et commenter les
résultats obtenus.
3/ Proposer une solution pour améliorer les résultats précédents.
4/ Trouver la représentation d'état liée à ce système.
5/ À partir de la question 4, trouvez une commande par retour d'état qui stabilise le
système et tracez les résultats sur MATLAB ou Simulink.
Données : 𝜁 = 0.6, 𝜔𝑛 = 11.11 𝑟𝑎𝑑/𝑠
2/ système électromécanique
Soit le système d’orientation des pales électromécanique suivant :
1 1
𝛽̇ = − 𝛽 + 𝛽𝑐
𝜏 𝜏
1/ Refaire les mêmes questions 1 à 4 de la partie précédente.
2/ Trouver la commande par retour d'état pour ce système.