Faculté de Technologie
Département de Génie Electrique
Module : Systèmes de conversion de l’énergie éolienne
Niveau : M1-Energies Renouvelables et efficacité énergétique
Année universitaire : 2023 / 2024 – Semestre1
TP N° 2 : Modélisation, Simulation, et optimisation de la turbine éolienne
Objectifs :
- Savoir implanter et simuler un schéma bloc par MATLAB-SIMULINK
- Savoir écrire et exécuter un algorithme sur Matlab
- Savoir visualiser et exploiter les résultats de simulation.
Paramètres de la turbine éolienne de type Savanius (Axe vertical)
Masse volumique de l’air ρ=1,2 (Kg/m3) Rayon de l’éolienne R= 0,5 =(m)
Hauteur de la turbine : H=2 (m) Rapport de multiplicateur de vitesse : G=1
Moment d’inertie total de la machine : J=0,16 Kg.m2; Cœfficient de frottement visqueux f=0.06 ;
Vitesse spécifique optimale : λopt ; Coefficient de puissance maximal : Cp-max
Le coefficient de puissance Cp est donné en fonction de la vitesse spécifique TSR (ou) λ par :
C p 0.21213 0.0856 2 0.2539
Profils de la vitesse de vent utilisés :
Profil N°1 Profil N°2
Time Series Plot:
14
13
V(t)=10+0.2sin(0.1047t)+2*sin(0.2665t)+sin(1.2930t)+0.2sin(3.6645t)
V ite s s e d e v e n t m /s
12
11
10
7
0 10 20 30 40 50 60
Temps s
Travail demandé : PARTIE I
1/ Écrire un algorithme permettant de tracer Cp=f(λ)
2/Déduire, TSR-opt et Cpmax
TSR=0:0.05:1.4;%0.05 est le pas de calcul
Cp=……………………………………………………;
plot(TSR,Cp);grid
axis([0 1.4 0 0.16])
title(’Titre du graphique’);
xlabel(’TSR’);%étiquette de l’axe x
ylabel(’Cp’)%étiquette de l’axe y
Responsable du module: MAYOUF Page 1
3/ Ecrire un algorithme permettant de tracer Paéro=f(Ω) pour
déférentes vitesses de vent
4/ Ecrire un algorithme permettant de tracer la courbe des
puissances optimales : Popt=f(Ω)
PARTIE II
1/ Compléter et simuler le schéma bloc de la turbine éolienne en utilisant le profil du vent N°1
2/ Simuler et visualiser sur le schéma bloc :
- L’allure du profil de vent
- La vitesse de rotation de l’arbre secondaire (à grande vitesse)
- La puissance aérodynamique Paéro
- La vitesse spécifique (λ ), et le coefficient de puissance (Cp)
3/ Traitement des courbes dans l’espace de travail
Reprendre et traiter les allures précédentes dans l’espace de travail en utilisant le bloc (to workspace)
Exemple :
%subplot(2,2,1),la fenètre de visualisation présente 2 lignes et 2 colonnes.
%%subplot(4,1,2), la fenètre de visualisation présente 4 lignes et 1 colonne.
subplot(2,2,1),plot(vent);%la courbe sera visualisée sur le 1er cadran à gauche parmi les 4 cadrans
axis([. . . .])
title(’Profil du vent’);
xlabel(’Temps , s’);
ylabel(’Vitesse de vent, m/s’) ;
4/ Reprendre les étapes 1 et 2 en utilisant le profil du vent N°2
Responsable du module: MAYOUF Page 2
Responsable du module: MAYOUF Page 3