Comparaison des Convertisseurs CUK et SEPIC
Comparaison des Convertisseurs CUK et SEPIC
كلية التكنولوجيا
Faculté de Technologie
قسم اإللكترونيـك
Département d’Électronique
BOUAMRA ASMA
HAMIDI OUSSAMA
Pour l’obtention du diplôme MASTER en Électronique option system de vision et
robotique
Thème
Bouamra Asma
REMERCEMENT
Je remercie en premier lieu Dieu qui m’a donné ce bien là et pour que
je vie ce jour et la force et la patience pour terminer ce travail.
Je tiens à exprimer mes remerciements à mon encadreur Prof.
[Link] qui a proposé et a dirigé ce travail.
Je remercie monsieur le président de jury, ainsi que les membres de jury
d’avoir accepté de juger ce travail.
Je tiens à remercier les responsables et tout le personnel du département
d'Electronique de Blida pour les facilités qu’ils m’ont accordés pour
terminer ce travail.
Le travail de mémoire que j’ai effectué doit beaucoup à certaines
personnes que je tiens à les remercier sincèrement.
Enfin, je remercie tous ceux qui ont contribué
de près ou de loin à ma formation et à l’élaboration de ce modeste
mémoire.
OUSSAMA HAMIDI
Taleb de matière
Introduction générale 01
Chapitre 01 Généralité sur l’énergie Solaire et Système PV
1.1 Introduction 05
1.2 L’énergie solaire 06
1.3 Le gisement solaire en Algérie 07
1.4 Avantages et Inconvénients de l’énergie Photovoltaïque 07
1.5 Cellule photovoltaïque 08
1.5.1 Caractéristiques électriques d’une cellule photovoltaïque réelle 08
1.5.2Paramètres d’une cellule photovoltaïque 10
1.1 1.6 Générateur photovoltaïque 11
1.6.1 Constitution d'un générateur photovoltaïque 11
1.6.2 Caractéristique Courant-Tension 11
1.6.3 Influence de l’éclairement et de la température sur les courbes I=f(V) et P=f(V) 12
1.7 Introduction d’un étage d’adaptation 15
1.8 conclusion 16
Chapitre 02 Les convertisseurs et la commande MPPT
2.1 Introduction 17
2.2 Convertisseurs DC-DC pour les systèmes d'énergie solaire 17
2.2.1 Topologie de convertisseurs DC/DC 18
2.2.2Convertisseur SEPIC 19
[Link] Synoptique d’hacheur SEPIC 20
[Link] Principe de fonctionnement 21
2.3 Dimensionnement d’un convertisseur SEPIC 23
2.3.1 Le rapport cyclique 23
2.3.2 Les inductances 24
2.3.3 Le condensateur de couplage 24
2.3.4 Condensateur de filtrage 24
2.3.5 Choix du MOSFET(Transistor) 24
2.3.6 Choix de la diode de sortie 25
2.3.7 Application numérique 25
2-4 Convertisseur CUK 26
2.4.1 Principe de fonctionnement 26
2.4.2 Dimensionnement des composants 28
2.4.3 Application numérique 29
59
Chapitre01 Généralité surl’énergie Solaire etSystème PV
1.1 Introduction
Le soleil est une source énergétique quasiment illimitée, cette énergie solaire est disponible
en abondance sur toute la surface terrestre, et malgré une atténuation importante lorsqu'elle
traverse l'atmosphère, la quantité qui reste est encore assez importante quand elle arrive au sol.
On peut ainsi compter sur 10 000 w/m² crête dans les zones tempérées et jusqu'à 14 000 W/m²
lorsque l'atmosphère est faiblement polluée. Malheureusement, toute cette énergie optique n’est
pas convertible en énergie électrique: en laboratoire, les meilleurs rendements de conversion
atteints par des cellules solaires (multicouches en l’occurrence) ne dépassent pas 40%.
L'énergie solaire photovoltaïque provient de la transformation directe d'une partie du
rayonnement solaire en énergie électrique aux moyens des cellules généralement à base de
silicium cristallin qui reste la filière la plus avancées sur le plan technologiques et industriel, en
effet le silicium et l’un des éléments les plus abondants sur terre sous forme de silice non toxique.
Cette conversion d'énergie s'effectue par le biais d'une cellule dite photovoltaïque (PV) basée sur
un phénomène physique appelé effet photovoltaïque.
L’énergie photovoltaïque consiste à produire une force électromotrice lorsque la surface de
cette cellule est exposée à la lumière. La tension générée peut varier en fonction du matériau
utilisé pour la fabrication de la cellule. L’association de plusieurs cellules PV en série/parallèle
donnent lieu à un générateur photovoltaïque (GPV) qui a une caractéristique courant-tension (I-V)
non linéaire présentant un point de puissance maximale(MPPT)[8]dépendant du niveau
d’éclairement, de la température de la cellule ainsi que du vieillissement de l’ensemble. Pour
pouvoir « récolter » et utiliser cette énergie, la connexion à une charge est nécessaire. Il en résulte
un point de fonctionnement correspondant au point d’intersection des caractéristiques
électriques du GPV et de la charge. Selon les choix des charges envisagées, le régime de
fonctionnement obtenu est optimisé ou non.
Ce chapitre sera consacré au rappel du principe et la performance de l'effet photovoltaïque
surla cellule photovoltaïque, dans le but de comprendre le fonctionnement et la production de
l’énergie par le générateur photovoltaïque et les principes de son optimisation.
5
Chapitre01 Généralité surl’énergie Solaire etSystème PV
Le nombre d’heures de puissance crête est obtenu par intégration de l’irradiance sur toutes
les heures de clarté (figure 1.1) ci-dessus.
6
Chapitre01 Généralité surl’énergie Solaire etSystème PV
La connaissance de cette valeur permet de quantifier rapidement les possibilités offertes par
le générateur solaire. Un module de 50W fournira 150WhH dans un site correspondant à un
ensoleillement de 3 heures de puissance crête [9][10][11].
7
Chapitre01 Généralité surl’énergie Solaire etSystème PV
- Sa structure fixe.
- Son coup de maintenance bas.
- Sa flexibilité (dimensionnement selon les besoins, modularité).
- L’installation ne produit aucun bruit.
- Son potentiel illimité, 5% de la surface des déserts suffiraient pour alimenter la planète entière.
Les inconvénients de l’énergie photovoltaïque sont :
- Le coût élevé de l’installation.
- Le rendement relativement bas de l’effet photovoltaïque.
- La puissance est réduite lorsque les conditions climatiques sont défavorables (nuages).
- Le stockage de l’énergie électrique sous forme chimique (batterie) est nécessaire pour une
installation autonome.
- Même si l’électricité produite par une installation photovoltaïque est sans pollution, la
fabrication, l’installation et l’élimination des panneaux ont un impact sur l’environnement.
Figure 1.2 Caractéristique I(V) d'une cellule photovoltaïque pour un éclairement et une température
donnés
Les résistances Rs et Rsh permettent de tenir en compte des pertes liées aux défauts de
fabrication.
8
Chapitre01 Généralité surl’énergie Solaire etSystème PV
• La résistance série (Rs) est la résistance interne de la cellule ; elle dépend principalement de la
résistance du semi-conducteur utilisé, de la résistance de contact des grilles collectrices et de la
résistivité de ces grilles.
• La résistance shunt (R sh) est due à un courant de fuite au niveau de la jonction; elle dépend de
la façon dont celle-ci a été réalisée.
I=Iph-Isat-Irsh (I.1)
Vpv+(Ipv∗Rs)
Irsh= (I.2)
Rsh
e(Vpv+(Ipv∗Rs))
Id=Io ∗ [exp( )− 1] (I.3)
nKT
e(Vpv+(Ipv∗Rs)) Vpv+(Ipv∗Rs)
Ipv = Iph − Isat [exp( ) − 1] − (I.4)
nKT Rsh
Où :
e(Icc∗Rs) (Icc∗Rs)
Icc = Iph − Isat [exp( ) − 1] − (I.5)
nKT Rs
Pour la plupart des cellules (dont la résistance série est faible), on peut négliger le terme
e(Icc∗Rs)
Isat [exp( ) − 1] devant Iph.L’expression approchée du courant de court-circuit est alors:
nKT
Iph
Icc ≅ Rs (I.6)
(1+Rsh)
eVpv (Vpv)
0 = Icc − Isat [exp( ) − 1] − (I.7)
nKT Rs
nKT Iph
Vco= ln[ Isat + 1] (I.8)
E
10
Chapitre01 Généralité surl’énergie Solaire etSystème PV
11
Chapitre01 Généralité surl’énergie Solaire etSystème PV
12
Chapitre01 Généralité surl’énergie Solaire etSystème PV
13
Chapitre01 Généralité surl’énergie Solaire etSystème PV
la tension de circuit ouvert. Quand la température augmente la tension de circuit ouvert diminue.
Par conséquent la puissance maximale du générateur subit une diminution [18].
14
Chapitre01 Généralité surl’énergie Solaire etSystème PV
15
Chapitre01 Généralité surl’énergie Solaire etSystème PV
1.8 Conclusion
Ce chapitre s’est consacré à présenter les différents composants d’un système
photovoltaïque, l’étude du principe de l’effet photovoltaïque, la cellule PV. Le rappel sur les
systèmes PV et leurs performancesbasés sur le modèle électrique équivalentainsi les différents
facteurs associés à la production et à l’adaptation d’énergie produite sont étudiéspour différents
variations météorologique.
16
Chapitre02 Les convertisseurs et la commande MPPT
2.1 Introduction
Comme on a vu dans le chapitre précédent, les caractéristiques d’un module solaire
montrent bien que la puissance maximale générée dépend fortement de l’intensité des radiations
solaires ainsi que de la température. En plus de ces dépendances, le module solaire(GPV) ne peut
générer sa puissance maximale que pour une certaine tension et courant de fonctionnement, la
position de ce point de fonctionnement dépendant à son tour de l’ensoleillement et de la
température ainsi que de la charge. Pour que le module fournisse sa puissance maximale
disponible il faut une adaptation permanente de la charge avec le générateur photovoltaïque.
Cette adaptation pourra être réalisée par l’insertion d’un convertisseur DC-DC (hacheur)
contrôlé par un mécanisme de poursuite «Maximum Power Point Tracking» (MPPT)[22][23] [24].
Le problème du couplage parfait entre un générateur photovoltaïque et une charge (surtout
de type continue) n’est pas encore réellement résolu. Un des verrous technologiques qui existe
dans ce type de couplage est le transfert de la puissance maximale du GPV à la charge qui souffre
souvent d’une mauvaise adaptation. Le point de fonctionnement qui en découle est alors parfois
très éloigné du point de puissance maximale PPM [25].
Le but de ce chapitre est l’étude des deux convertisseurs Ćuk et SEPIC qui sont utilisées dans
les systèmes photovoltaïques, on décrit aussi la commande MPPT de ces deux convertisseurs on
spécifie les deux algorithmes (perturbation et observation (P&O) et incrémentation de
l’inductance).
La conversion DC/DC (continue-continue) a pour but de contrôler l'écoulement de l'énergie
entre le panneau solaire et la charge. Pour ce type de convertisseur, la gamme de puissance que
l'on peut traiter s'étend de quelques watts, à des centaines de kilowatts. Cependant, l'apparition
de composants de puissance sans cesse plus performants, permet d'envisager le traitement de
puissances plus élevées.
17
Chapitre02 Les convertisseurs et la commande MPPT
18
Chapitre02 Les convertisseurs et la commande MPPT
Boost 1 NON
1−α
Buck-Boost −α NON
1−α
Ćuk −α NON
1−α
SEPIC α NON
1−α
Flyback α OUI
𝐾
1−α
Push-pull Kα OUI
Forward Kα OUI
Tableau II.1 Rapport de transformation des principaux convertisseurs DC-DC
19
Chapitre02 Les convertisseurs et la commande MPPT
20
Chapitre02 Les convertisseurs et la commande MPPT
a- De 0 à α T : Phase de conduction
Pendant la phase de conduction, S est en état passant (ON), la tension aux bornes de
l'inductance L1 est égale à Vin, le condensateur C1 est connecté en parallèle avec L2et la tension
aux bornes de L2 est la même que la tension du condensateur = -Vin. Voir Figure 2.4. La diode D
est en polarisation inverse donc bloquée. Le courant de la charge est fourni par le condensateur
C2. Pendant cette période le courant est stocké en L1 par l'entrée et en L2 par C1.
b- De α T à T :
État de récupération lors de l’ouverture de l’interrupteur S, la diode D se met à conduire.
Pendant la phase de récupération, S est en état bloqué (OFF), le courant dans L1 continue à
circuler à travers C1 et la diode D puis dans le condensateur C2, durant cet intervalle le
condensateur C1 est chargé et prêt pour le prochain cycle. Aussi le courant emmagasiné dans L2
circule dans C2 et la charge, rendant C2 prêt pour le prochain cycle Figure 2.5 [12].
21
Chapitre02 Les convertisseurs et la commande MPPT
22
Chapitre02 Les convertisseurs et la commande MPPT
Avec : V𝐼𝑁 est la tension d’entrée, Vout est la tension de sortie et 𝑉𝐷 la tension de seuil de la
diode. Avec un rapport cyclique proche de 50%, la tension d’entrée est égale à la tension de sortie,
le gain est unitaire. Le rapport cyclique varie en fonction de la tension d’entrée appliquée afin
d’avoir une tension de sortie constante. Ainsi le rapport cyclique maximal est:
Vout+VD
αmax = V (II.2)
IN(min) +Vout+VD
23
Chapitre02 Les convertisseurs et la commande MPPT
Avec Iout est le courant dans la charge, Vripple est l’ondulation de la tension de sortie.
24
Chapitre02 Les convertisseurs et la commande MPPT
supporter une tension 𝑉𝑆 = 𝑉𝐼𝑁 +𝑉𝑂𝑈𝑇 , et doit supporter un courant crête égal àIS(peak) =
IL1(peak)+IL2(peak) .
(Vout+VIN(min) +VD)×(Vout+VD)
IS(rms) = Iout × √ (II.8)
V2IN(min)
VIN = 26.3 V
IOUT = 3.125 A
Rload = 15.5Ω
F =100 KHZ
A- Rapport cyclique
VIN(min) 26.3
L1 = L2 = L = × αmax = 0.57∗105 × 0.646 =0.298 mH
∆I L ∗f
25
Chapitre02 Les convertisseurs et la commande MPPT
C- Le condensateur de couplage
Le courant efficace dans C1 sera calculé par :
Vout + VD 48
IC1(rms) = Iout × √ = 3.125 × √ = 2.29A
VIN(min) 26.3
En prenant une ondulation de tension 0.19, la valeur du condensateur de couplage s’obtient par:
Iout×αmax 3.125×0.646
C1 = = =106.25µF
∆V C1 ×f 0.19×105
D- Le condensateur de filtrage
La tension de sortie accepte une ondulation de 1V.
Iout×αmax 3.125×0.646
C2 = = =40.37µf
Vripple×0.5×f 1×0.5×105
26
Chapitre02 Les convertisseurs et la commande MPPT
IL1 α
= 1−α (II.13)
I L2
VIN I
= I L1 (II.16)
V OUT L2
27
Chapitre02 Les convertisseurs et la commande MPPT
• α=0.5 V OUT = V IN
VIN=26.3 V
IIN =7.61 A
Vout=48 V
P=150W
IOUT = 8.2 A
Rload = 15.36 Ω
a- Rapport cyclique
Vin α VOUT
= Donc : α=
VOUT 1−α VOUT +VIN
b- Les inductances
Les valeurs des inductances sont choisies de telles sorte que ∆IL1 n’est pas supérieur à 5% de
IIN et ∆IL2 n’est pas supérieur à 5% de IL2 .
On calcule L par la relation :
∆IL1 =0.05× IL1 =0.05×7.61=0.3805 A
V ×α
Donc : L1 =∆Iin (II.18)
L1 ×f
V ×α
L2 = ∆Iin ×f (II.19)
L2
c- Les condensateurs
Le critère de choix des condensateurs est que la tension d’oscillation à travers eux doit être
moins de5%.La tension moyenne à travers le condensateur (C1) :
28
Chapitre02 Les convertisseurs et la commande MPPT
b- les inductances
V ×α 26.3×0.646
L1 =∆Iin ×f=L1 = 0.3805×105=0.44 mH
L1
V ×α 26.3×0.646
L2 = ∆Iin = = 0.41 mH
L2 ×f 0.41×105
c- les condensateurs
VOUT ×α 48×0.646
C1 = = =5.434 µF
∆VC1×Rload ×f 3.715×15.36×105
1−α 1−0.646
C2 = ∆V = 0.05 =10.36µF
8×( OUT )L2f2 8×(
48
)×0.41×10−3 ×(105 )2
VOUT
d- choix de diode
Dans ce cas, nous choisissons une diode de faible seuil à l’état passant pour minimiser la
perte de puissance par conduction et faible temps de commutation(en général 5 à10nS), la tension
inverse maximale de la diode est la même que la tension moyenne de condensateur
VC1 égale à74.3V,le courant direct moyen de la diode est la combinaison des courants de sortie et
d’entrée lorsque le Transistor de commutation est fermé, alors il est :
ID = IL1 + IL1 (II.22)
e- choix du Transistor de commutation
29
Chapitre02 Les convertisseurs et la commande MPPT
Les transistors MOSFET sont les plus utilisés pour les domaines de faible et moyenne
puissance .la tension crête à du commutateur (S) est obtenue par l’application de la loi de
Kirchhoff lorsque le commutateur est fermé :
dIL1
VOUT = IIN − (II.23)
dt
30
Chapitre02 Les convertisseurs et la commande MPPT
autour de la tension idéaleVmp de fonctionnement. Ceci cause une perte de puissance qui dépend
de la largeur du pas d'une perturbation simple Le processus est répété périodiquement jusqu’à
ce que le MPP soit atteint. Le système oscille alors autour du MPP, ce qui provoque des pertes de
puissance. L’oscillation peut être minimisée en diminuant la taille de la perturbation. Cependant,
une taille de perturbation trop petite ralentit considérablement la poursuite du MPP. Il existe alors
un compromis entre précision et rapidité [34].
Si la largeur du pas est grande, l'algorithme MPPT répondra rapidement aux changements
soudains des conditions de fonctionnement, mais les pertes seront accrues dans les conditions
stables ou lentement changeantes [35].
Si la largeur du pas est très petite, les pertes dans les conditions stables ou lentement
changeantes seront réduites, mais le système ne pourra plus suivre les changements rapides de la
température ou de l'insolation.
L’inconvénient de la méthode P&O est lorsqu’une augmentation brutale de l’ensoleillement
est produite on aura une augmentation de la puissance du panneau, l’algorithme précédent réagit
comme si cette augmentation est produite par l’effet de perturbation précédente, alors il continue
dans la même direction qui est une mauvaise direction, ce qu’il éloigne du vrai point de puissance
maximale.
Ce processus continu jusqu’à la stabilité de l’ensoleillement où il revient au vrai point de
puissance maximale. Ceci cause un retard de réponse lors des changements soudains des
conditions de fonctionnement et des pertes de puissance.
31
Chapitre02 Les convertisseurs et la commande MPPT
dP/dV=0 au MMP
dP d(IV) dI ∆I
Puisque : = = I+V ≃ I+V
dV dV dV ∆I
32
Chapitre02 Les convertisseurs et la commande MPPT
Alors:
∆I I
=− au MPP
∆V V
∆I I
> −V à Gauche de MPP
∆V
∆I I
<− V à droite de MPP
∆V
point de fonctionnement correspondant est maintenu, à moins q'un changement dans∆Iest noté,
indiquant un changement des conditions atmosphériqueset donc du MPP.L'algorithme augmente
ou diminueVref pour suivre le nouveau MPP.
La taille de l’incrément détermine la rapidité de la poursuite du MPP. Une poursuite rapide
peut être obtenue avec un incrément plus grand, mais le système ne pourrait pas fonctionner
exactement au MPP et oscille autour de celui-ci. Il y a donc, comme pour la méthode P&O, un
compromis entre rapidité et précision. Cette méthode peut être améliorée en amenant le point de
33
Chapitre02 Les convertisseurs et la commande MPPT
fonctionnement près du MPP dans une première étape, puis d’utiliser l’algorithme conductance
incrémentale pour suivre exactement le MPP dans une deuxième étape.
2.6 Conclusion
Dans ce chapitre l’étude et le dimensionnement du convertisseur SEPIC et CUK, leurs
fonctionnements et ses performances ont été présentés. Les directives d'optimisation de la
conception, nous ont aidés à donner un bon dimensionnement aux convertisseurs.
Nous avons aussi présenté deux méthodes de la commande MPPT pour les deux
convertisseurs SEPIC et CUK pour la rechercher du point maximale de puissance du générateur
photovoltaïque, à savoir la méthode perturbe et observe et l’incrémentale conductance.
34
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
3.1 Introduction
Dans ce chapitre nous abordons deux parties distinctes, la première partie sera
consacrée à la simulation des deux convertisseurs SEPIC et CUK adapté par deux commandes
MPPT numérique «Perturbation et Observation, Incrémentale Conductance» appliquées à un
système photovoltaïque. Dans le but de comparer les performances de ces deux convertisseurs
contrôlés chacun par les deux algorithmes développés dans le chapitre précédent, les
architectures sont implémentées dans l’environnement SIMULINK. Puis les performances des deux
convertisseurs contrôlés par un seul MPPT, à savoir Incrémentale Conductance sera implémenté
sur cible FPGA en utilisant XILINX System Generator. Et enfin une comparaison des deux
performances pour les deux environnements sera discutée.
35
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
Description Estimation
Inductance(L1) 298 µH
Inductance(L2) 298 µH
Condensateur(C1) 106.25 µF
Condensateur(C2) 47 µF
MOSFET(S) STW40NF20
Charge Résistive (R load) 15.5 Ω
Diode de puissance BYT30-1000
Fréquence de coupure 100 khz
Rapport cyclique α 0.646
Tableau III.1 Valeur des composants pour convertisseur SEPIC
36
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
37
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
[Link].2 Etage d’adaptation convertisseur CUK
Le système photovoltaïque composé de panneau photovoltaïque, un MPPT, et un étage
d’adaptation représentée par un convertisseur CUK.
38
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
[Link] La comparaison des deux Algorithmes MPPT pour les deux
convertisseurs
On a écrit les deux Algorithmes dans Embedded Matlab Function BLOC comme suit (voir
annexe I) :
[Link].1 Algorithme P&O (perturb and observe) appliqué au convertisseur
Ćuk
39
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
de code matlab qui donne a sa sortie le rapport cyclique α ;et un bloc PWM, nous intéresse de
rapport α qui soit injectée directement au interrupteur (MOSFET) de convertisseur SEPIC
[Link].3La comparaison des deux algorithmes pour le convertisseur Ćuk
Dans ce cas on a choisis le convertisseur Ćuk qui est présente dans la figure avec les deux
Algorithme P&O et InCn et un générateur GPV KC200GT, est on obtient le graphe comme suit:
Figure3.7 Système photovoltaïque avec convertisseur Ćuk contrôlé par deux algorithmes MPPT
sous SIMULINK
40
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
Figure 3.8 Rapport cyclique des deux MPPT pour convertisseur Ćuk
Les deux figures A et B présentent la variation de rapport cyclique des deux Algorithme
A présente le rapport cyclique de MPPT perturbe and observe.
B présente le rapport cyclique de MPPT conductance incrémentale.
Figure3.9 Puissance à la sortie de convertisseur Ćuk contrôlé par deux algorithmes (P&O et InCn)
41
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
Les deux graphes de figure (3.9)présentent les résultats de la simulation.
On remarque qu’avec l’MPPT P&O les ondulations de puissance dans l’état transitoire est inferieur
par rapport àl’MPPT de InCn, par contre les puissances de sortie à partir de l’instant t=0.02 S
devient légèrement différente pour les deux algorithmes, il est clair que la puissance pour le cas de
InCn est supérieur à celle de cas P&O.
[Link].4 La comparaison des deux algorithmes pour le convertisseur SEPIC
Dans ce cas on a choisit le convertisseur SEPIC qui est présenté dans la figure 3.10 avec les
deux Algorithme P&O et InCn et un générateur électrique GPV catégorie KC200GT qui donne à se
sortie une tension de 26.3V, est on obtient les résultats ci-dessous:
42
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
Figure 3.11 Système photovoltaïque avec convertisseur SEPIC contrôlé par deux algorithmes
MPPT sous SIMULINK
43
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
Figure 3.12 Rapport cyclique des deux MPPT pour le convertisseur SEPIC
Les deux figures A et B présentent la variation de rapport cyclique des 2 Algorithme pour un
convertisseur SEPIC
- La figure A présente le rapport cyclique de MPPT Conductance Incrémentale.
- La figure B présente le rapport cyclique de MPPT Perturb and Observ.
Après la simulation on obtient le graphe suivant :
44
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
Figure3.13 Puissance à la sortie de convertisseur SEPIC contrôlé par deux algorithmes (P&O et
InCn)
On remarque que les ondulations dans l’algorithme MPPT P&O sont inférieures par rapport
au MPPT InCn, on remarque aussi que la puissance délivrée avec l’Algorithme InCn est plus
grande.
Jusqu'à t=0.02s les deux puissances sont presque les mêmes, à partir de l’instant t=0.02 la
puissance devient légèrement différente pour les deux algorithmes, il est clair que la puissance
pour le cas de InCn est supérieure à celle du cas P&O.
Dans ce cas on peut dire que le meilleur Algorithme de MPPT pour notre système est
l’algorithme InCn (conductance incrémentale), parce qu’il nous donne une puissance optimale.
45
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
Figure 3.14 Convertisseurs SEPIC et Ćuk contrôlé par l’algorithme InCn sous SIMULINK
Les résultats de simulation obtenues sont représentés par les figures suivantes:
46
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
Figure 3.15 Rapport cyclique de l’algorithme InCn pour les deux convertisseurs
La figure suivante représente la puissance à la sortie des deux convertisseurs:
Figure 3.16 Puissances à la sortie des convertisseurs Ćuk et SEPIC contrôlés par l’algorithme InCn
Les figures présentent la puissance à la sortie des deux convertisseurs. En remarque que :
- les ondulations dans le cas de convertisseur Ćuk sont plus atténuées que dans le cas de
SEPIC .
- Les deux puissances à la sortie des convertisseurs Ćuk et SEPIC sont presque les mêmes
- l’état transitoire dans le cas de convertisseur Ćuk atteint son régime de stabilisation après
celui de convertisseur SEPIC. Le convertisseur Ćuk se stabilise entre 0.002 S jusqu’à 0.003 S par
contre le convertisseur SEPIC se stabilise après 0.01 S.
Dans ce cas on peut dire que le convertisseur Ćuk est meilleur que le convertisseur SEPIC.
47
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
[Link].2 Comparaison des deux convertisseurs contrôlés par P&O
Dans ce cas on a choisit un MPPT avec l’Algorithme P&O (perturb and observ) qui est
représenté dans la figure ci-dessous avec les deux convertisseurs
Figure 3.17 Convertisseurs SEPIC et Ćuk contrôlés par l’algorithme P&O sous SIMULINK
48
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
Figure 3.19 Puissance à la sortie des deux convertisseurs contrôlés par l’algorithme P&O
Les graphes présentent la puissance à la sortie des deux convertisseurs. En remarque que :
- les ondulations dans le cas du convertisseur SEPIC sont importants par rapport au cas du
convertisseur Ćuk.
- La valeur de la puissance délivrée par le convertisseur Ćuk dépasse légèrement celle
délivrée par le convertisseur SEPIC dans l’intervalle 0.01s jusqu’à 0.06s.
- l’état transitoire de convertisseur Ćuk est inférieur à celui du convertisseur (SEPIC), le
convertisseur Ćuk atteint la stabilité après un temps max de 0.002s, par contre le convertisseur
SEPIC dépasse un temps de 0 .01s pour être stable.
Dans ce cas on confirme que le convertisseur CUK est meilleur que convertisseur SEPIC.
49
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
3.3 Implémentation du système sur XSG :
3.3.1INTRODUCTION
FPGA (field-programmable gatearray, réseau de portes programmables in situ)est un circuit
logique programmable, ou réseau logique programmable, est un circuit intégré logique qui peut
êtrereprogrammé après sa fabrication. Notons qu'il serait impropre de parler de programmation
au sens logiciel (contrairement à un microprocesseur, il n'exécute aucune ligne de code). Ici, mieux
vaudrait parler de « reconfiguration » plutôt que de reprogrammation (on modifie des connexions
ou le comportement du composant, on connecte des portes logiques entre elles, etc.). Le verbe
programmer est toutefois le plus fréquent, mais au sens de personnaliser. Il s'agit bel et bien de
réseaux logiques reconfigurables et modifiables.
Il est composé de nombreuses cellules logiques élémentaires et bascules logiques librement
connectables (c'est justement la reconfiguration, ou programmation, du composant qui définit les
connexions faites entre portes logiques)[38].
50
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
spécifier la manière dont la génération de code et de simulation doivent être manipulée figure
4.16 [44].
La précision et les passerelles Xilinx SG fonctionne avec les modèles Simulink standard. Deux
blocs appelés "Gateway In" et "Gateway Out" définissent la limite du FPGA du modèle simulé en
Simulink voir la figure 4.17. La passerelle dans le bloc convertit l’entrée à un nombre à virgule fixe
[40].
« Gateway In »block est utilisé pour définir la limite d’entrée FPGA et convertit la virgule
flottante à pointfixe.
Le « Gateway Out » block est utilisé pour définir la limite de sortie FPGA et convertit à la double
précision.
51
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
A B
Figure3.23 Tableauxdu Choix de nombre de bit de 𝐕𝐩𝐯 𝐞𝐭 𝐈𝐩𝐯
La figure présente le choix des nombres de bits pour les entrés I et V de notre MPPT.
On a choisit un nombre de bit dans les deux tableaux A et B
Dans le tableau A on a une valeur de tension V=32.9 V lorsqu’on a convertit en binaire on
trouve 25 =32 pour la partie entière ça veut dire qu’on doit réserver 5 bits pour cette partie, la
partie fractionnelle doit être le plus précise possible à 0.9 donc on a choisit un nombre de bit de 16
bits, par ce que la valeur de la partie fractionnelle doit être le plus proche possible de la valeur
calculée précédemment.
Dans le tableau B on a une valeur de courant I=8.2 A lorsqu’on a convertit en binaire on
trouve 24 =8 pour la partie entière ça veut dire qu’on doit réserver 4 bits pour cette partie, la
partie fractionnelle doit être la plus précise possible à 0.2 donc on a choisis un nombre de bit de
16 bits, par ce que la valeur de la partie fractionnelle doit être la plus proche possible que valeur
calculée.
52
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
La figure suivante illustre L’architecture de l’algorithme InCn sous XSG.
53
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
A
Figure 3.25 Comparaison de la puissance de sortie POUT pour les deux convertisseurs contrôlés par
un MPPT InCn sous XSG et Simulink
54
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
Le graphe A représente les résultats obtenus sous XSG et le graphe B ceux obtenus sous
SIMULINK. On analysons ces résultats on constate une bonne concordance entre les deux cas pour
différents niveaux d’irradiation. Ceci justifie amplement le bon choix de notre architecture sous
XSG. En conséquence la matérialisation de la commande devient plus aisée en ciblant une carte
FPGA.
55
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
56
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
Les deux figures représentent le circuit de connexion de notre MPPT sur FPGA.
57
Chapitre03 Simulation et implémentation sous SIMULINK et XSG
3.3.3 Conclusion :
Dans la première partie de ce chapitre nous avons simulé notre MPPT sous
Simulink/MATLAB pour différents niveaux d’irradiation avec deux convertisseurs. Les résultats
obtenus sont acceptables.
Cette étude présenté une comparaison de deux MPPT plus populaire contrôleur, perturb and
Observ Contrôleur avec incrémental conductance Controller. Cette orientation du papier sur la
comparaison des deux convertisseurs différent qui sera relié avec le dispositif de commande. Un
panneau solaire simple qui a la norme la valeur de l'insolation et la température a été inclus dans
la simulation circuit. De tous les cas, le meilleur contrôleur de MPPT est incrémentalecontrôleur de
conductance. Ce contrôleur donne une meilleure valeur de sortie pour son argent, stimuler et les
deux convertisseurs Cuk et SEPIC. Ainsi ce contrôleur donnera genre différent de courbes pour
l'ensemble des convertisseurs.
Les résultats montrent que tous les systèmes détecte la puissance maximale Point bien SEPIC
est moine stable avec moins d'entraînement de puissance par rapport au convertisseur Ćuk à la
puissance maximale sur D'autre part, l'avantage de convertisseur Ćuk est la réduction du
paramètres de circuit (inductance et de condensateurs) par rapport à Convertisseur SEPIC et le
coût réduit le rendement. La sélection de l'une de deux convertisseurs dépend de l’exigence du
système et le budget. Pour une bonne qualité de la production la puissance, le convertisseur Ćuk
est favorable.
Dans la deuxième partie la description et la simulation sous XSG du MPPTInCnNous avons
simulé l’algorithme InCn sous MATLAB, Simulink, puis générer du code pour l’ FPGA Xilinx. Nous
avons utilisé des blocs spécifiques. Les résultats obtenus sont comparables avec ceux obtenus par
58
Liste des tableaux
1
Introduction général
2
Introduction général
3
Introduction général
4
Bibliographie
[14] Bahrame Khadabakkchian, Fouis.A Dessaint et Michel lavoie, "simulation en temps réel
complément numérique de grand raison de transport EHT, en vue d’analyse et de tester leurs
circuits de commande et de protection", IEE Canadian preview-summer/été 1999.
[15] C. Lerouge, et al.."Recherche & Industrie Photovoltaïque (PV) aux Etats-Unis 2006", SCIENCES
PHYSIQUES ETATS-UNIS.
[16] MULLER, J.-C., "Électricité photovoltaïque, in BE 8578. 2007", Techniques d'ingenieurs.
[17] Antonio Luque and Steven Hegedus,"Handbook of Photovoltaic Science and Engineering",
John Wiley & Sons Ltd, 2003.
[18] [Link], "A Modiffied Adaptative Hill Climbing Maximum Power Point Tracking (MPPT) Control
Method For Photovoltaic Power Systems", The University of British Columbia, 2003.
[19] Akihiro Oi, "Design and Simulation of Photovoltaic Water Pumping System", Faculty of
California Polytechnic State University, 2005.
[20] S. Issaadi, "Commande d’une poursuite du point de puissance maximum (MPPT) par les
Réseaux de Neurones ", Mémoire de magister, Ecole Nationale Polytechnique, Alger, Algérie,
2006.
[21] M. Adouane, "Etude et conception d’une stratégie de commande d’un onduleur connecté au
réseau électrique", Mémoire de magistère, Ecole Nationale Polytechnique, Alger, Algérie, 2008.
[22] S. Bellakehal, "Conception & Commande des Machines à Aimants Permanents Dédiées aux
Energies Renouvelables", thèse de Doctorat, Université de Constantine, Algérie, 2010.
[23] M. Makhlouf, "Etude et optimisation d'un modèle de conversion d'énergie photovoltaïque
application au pompage", mémoire de magister, université Mentouri Constantine, Algérie, 2006.
[24] Cedric Cabal," Optimisation énergétique de l’étage d’adaptation électronique dédié à la
conversion photovoltaïque.", UNIVERSITE TOULOUSE III –juin 2008.
[25] S. M. Ait-Cheikh, "Etude, Investigation et conception d’algorithmes de commande appliqués
aux systèmes photovoltaïques", Thèse de Doctorat d’état, Ecole Nationale Polytechnique, Alger,
Algérie, 2007.
[26] N. Abouchabana, " Etude d’une nouvelle topologie buck-boost appliquée à un MPPT ",
Mémoire de Magistère, Ecole Nationale Polytechnique, Alger, Algérie, 2009.
[27] M. Moulay-Amar et M. Loghouini, "Etude et réalisation d’un système de poursuite de point de
puissance maximale à base de microcontrôleur destiné a une installationphotovoltaïque ",
Mémoire d’ingéniorat, université de Ouargla, Algérie, 2005.
Bibliographie
[28] Angel Cid Pastor, "Conception et Réalisation de Modules Photovoltaïques Electroniques, thèse
de doctorat, soutenue le 29 septembre 2006.
[29] Fang Lin Luo, Hong Ye , Muhammad Rashid . "Digital Power Electronics and Applications"
,Elsevier Academic press,2005.
[30] Muckund R,Ph.D.,P.E. "Wind AND Solar Power Systems".CRC Press LLC,1999.
[31] [Link]
[32] Joris Pallier, ‘Dimensionnement des composants pour convertisseur SEPIC’, Projet Polytech-
Clermont-Ferrand.
[33] M. Djarallah : ‘contribution a l’étude des systèmes photovoltaïques résidentiels couples au
réseau électrique’, Thèse de doctorat université de Batna 2008.
[34] L. Protin et al., Convertisseurs photovoltaïques, Technique de l’ingénieur, traité génie
électrique, pp 3-360, 1992.
[35] S. Masoum et al., Microprocessor-controlled new class of optimal battery chargers for
photovoltaic application, Vol. 19, September 2004.
[36] F. Chekired, Etude et implémentation d’une commande MPPT neuro-floue sur FPGA,
Mémoire de Magister de l’Ecole Nationale Polytechnique, Alger, 2008.
[37] Trishan Esram, Student Member, IEEE, and Patrick L. Chapman, Senior Member, IEEE,
"Comparison of Photovoltaic Array Maximum Power Point Tracking Techniques", IEEE
TRANSACTIONS ON ENERGY CONVERSION, 2006.
[38] J. Fadat, "Chaîne de conversion d’énergie photovoltaïque, Commande MPPT", La Revue 3EI n°
43 – Décembre 2005, ESIEE-Amiens 14 Quai de la Somme – BP 100 – 80 082 Amiens.
[39] A. Oi, "Design and simulation of photovoltaïque water pumping système", Thèse de Magister
de l’université California, USA, Septembre 2005.
[40] A. Meflah, « Modélisation et commande d’une chine de pompage photovoltaïque», Mémoire
de Magistère, Université de Tlemcen, Algérie, 2011.
[41] [Link]
[42] htt:/xilinx system generator and HDL coder-FPGA.
[43] JEAN-GABRIEL MAILLOUX, ’ PROTOTYPAGE RAPIDE DE LA COMMANDE VECTORIELLE SUR
FPGA À L'AIDE DES OUTILS SIMULINK - SYSTEM GENERATOR’, Mémoire dans le cadre du
programme de maîtrise d’ingénierie, l’université du QUÉBEC à CHICOUTIMI, MARS 2008.
Annexes
Photovoltaic Module KC200GT
Le panneau photovoltaïque
function i = fcn(v,t,g,ns)
i=15;
tn=298;
gn=1000;
ki=0.0032;
ipvn=8.214;
k=1.380*10^-23;
q=1.602*10^-19;
kv=-0.1230;
vocn=32.9;
a=1.3;
rs=0.221;
rp=412.405;
iscn=8.21;
ipv=(g/gn)*(ipvn+(ki*(t-tn)));
vt=(ns*k*t)/q;
i0=(iscn+(ki*(t-tn)))/(exp((vocn+(kv*(t-tn)))/(vt*a))-1);
i=ipv-(i0*(exp((v+rs*i)/(vt*a))-1))-((v+rs*i)/rp);
Annexes
3.2.1 Algorithme P&O( perturbe and observe)
function D = PO(V,I,T)
persistent P2 P1 dP d dd n;
if isempty(V)
V=29.53; %V=30;
end
if isempty(I)
I=0;
end
if isempty(P2)
P2=0;
end
if isempty(P1)
P1=0;
end
if isempty(dP)
dP=0;
end
if isempty(d)
d=1;
end
if isempty(dd)
dd=0;
end
if isempty(n)
n=1;
end
%%%%%%%%%%%%%%%%%%%%%
if (T > n*0.02) %chu ki lay mau 0.02s
n = n + 1;
P1=P2;
P2=V*I;
dP=P2-P1;
if (dd==0)
if dP>1
dd=0.01;
d=d+dd;
else
if dP<-1
dd=-0.01;
d=d+dd;
else
dd=0;
Annexes
end
end
else
if ((dP<1)&&(dP>-1))
dd=0;
d=d+dd;
else
if ((dP/dd)>0)
dd=0.01;
d=d+dd;
else
dd=-0.01;
d=d+dd;
end
end
end
end
D=d;
if D~=0.68
d=0.68;
end
end
3.2.2 Algorithme incremental conductance (InCn)
function y= MPPT(V,I,V0,I0,D)
m=0;
dV=V-V0;dI=I-I0;
if dV==0
if dI==0,m=D;
else
if dI>0,m=D-0.01;
else
m=D+0.01;
end
end
else if dI/dV==-(I/V)
else
if dI/dV>-(I/V),m=D-0.01;
else
m=D+0.01;
end
end
end
y=m;