These Buck Boost
These Buck Boost
THESE
Prsente Par
Youssef TOUZANI
Ingnieur Polytech' Marseille
DOCTEUR DUNIVERSITE
Spcialit : Electronique et systmes
Que Monsieur Michel DHOME, Professeur lUniversit Blaise Pascal, trouve ici ma
profonde reconnaissance pour avoir accept de faire partie de ce jury.
Je risque de ne jamais trouver les mots justes pour exprimer ma reconnaissance mon
directeur de thse Monsieur Khalil EL KHAMLICHI DRISSI, Professeur lUniversit
Blaise Pascal. Je le remercie pour son aide constante dans la phase finale de ce travail, pour
ses marques de soutien des moments importants, pour ses conseils prcieux et pour son
amiti.
Je remercie mes ex collgues de bureau Nadia, Fulbert et David qui mont soutenu
pendant le droulement de ma thse.
Enfin, je ne saurais oublier mon pouse et mes parents qui mont apport soutien est
encouragement durant toutes ces annes.
Table de matires
Table de matires
Introduction gnrale............................................................................................................................ 5
Introduction........................................................................................................................................................... 9
Conclusion ........................................................................................................................................................... 37
Introduction......................................................................................................................................................... 39
Conclusion ........................................................................................................................................................... 73
Introduction......................................................................................................................................................... 75
Introduction....................................................................................................................................................... 103
Introduction....................................................................................................................................................... 132
Bibliographie...................................................................................................................................... 171
Introduction gnrale
- 5 / 188 -
Introduction gnrale
Dans le premier chapitre, aprs un bref rappel de la pollution harmonique, de ses effets ainsi
que des normes en vigueur, les solutions existantes pour minimiser cette pollution sont
prsentes. En outre, une tude comparative des diffrentes familles des redresseurs
prlvement sinusodal, des contraintes lies leur fiabilit et des diffrentes techniques de
modulation susceptibles de rduire les harmoniques et amliorer leur qualit de point de vue
de la compatibilit lectromagntique sont dveloppes.
- 6 / 188 -
Introduction gnrale
total des capteurs utiliss dans le systme de commande. Cette architecture est caractrise
par des algorithmes de reconstruction des diffrents signaux dans un environnement temps
rel, et avec des contraintes lies aux variations alatoires de la modulation. Cet
environnement et ces contraintes imposent la prcision et la robustesse des mthodes et
techniques utilises dune part, et la rapidit dexcution des instructions dautre part. Ce
compromis est le critre principal du choix des mthodes retenues dans cette architecture de
commande.
- 7 / 188 -
Chapitre 1 : Etat de lart
Conclusion ........................................................................................................................................................... 37
- 8 / 188 -
Chapitre 1 : Etat de lart
Introduction
Au cours des dix dernires annes, on a assist une multiplication des convertisseurs
dnergie lectrique raccords au rseau de distribution. Ces charges non linaires absorbent
des courants non sinusodaux et les courants perturbateurs (harmoniques) de telles charges
sont lorigine de distorsions sur londe de tension au travers de limpdance de court-circuit
au point de raccordement. Cest cette pollution harmonique des rseaux, qui remet donc en
cause leur aptitude dlivrer une tension sinusodale quelle que soit la forme du courant
absorb par la charge.
Plusieurs solutions existent pour lutter contre cette pollution harmonique. Le redressement
correction de facteur de puissance dispos en srie avec le convertisseur incrimin et le
filtrage actif plac en parallle avec le convertisseur polluant. Ces deux solutions permettent
de limiter la pollution harmonique la source.
Aprs une brve description de la pollution harmonique, des nuisances occasionnes par le
raccordement des charges non linaires au rseau, et des solutions pour lutter contre cette
pollution, nous prsenterons dans ce chapitre les diffrentes topologies des convertisseurs
prlvement sinusodal du courant. Nous voquerons ensuite les contraintes de fiabilit lies
ce type de convertisseurs. Nous rappellerons les diffrentes techniques de modulation
susceptibles dtre utilises par les convertisseurs prlvement sinusodal, ainsi que les
possibilits quelles offrent pour amliorer la qualit de ces convertisseurs du point de vue de
la compatibilit lectromagntique. Nous terminerons par une brve prsentation des objectifs
de notre tude ainsi que des outils utiliss.
I. La pollution harmonique
I.1. Principales charges polluantes
Outre les machines et les transformateurs, les sources principales de la pollution harmonique
ont pour origine les charges non linaires, principalement les fours arc (en ce qui concerne
les fortes puissances) et les convertisseurs d'nergie. On qualifie de charges non linaires des
- 9 / 188 -
Chapitre 1 : Etat de lart
charges qui absorbent un courant non sinusodal, donc des courants harmoniques, mme
lorsquelles sont alimentes par des tensions parfaitement sinusodales.
Les convertisseurs statiques sont les plus reprsentatifs de ces charges en terme de puissance
dformante gnre :
- Les montages gradateurs, que ce soit pour les entranements faible puissance, pour les
systmes dclairage et de chauffage ou pour les systmes de conduite des rseaux.
- Les redresseurs de tension diodes ou thyristors, prsents dans la mtallurgie ou encore
dans les variateurs de vitesse.
- Les redresseurs de courant diodes avec filtrage capacitif, prsents dans les alimentations
des ordinateurs ou des tlviseurs.
A titre dexemple, la forme donde et le spectre du courant absorb par un redresseur triphas
diodes sont prsents sur la figure I.1.1.
15,00
20
12,50
10
10,00
0
7,50
-10
5,00
-20
2,50
-30 0,00
0 10 20 30 40 50 60 0,00 0,50 1,00 1,50 2,00
- 10 / 188 -
Chapitre 1 : Etat de lart
fait que la section du fil du neutre est gnralement plus petite que celle des fils de phase, car
en rgime sinusodal le neutre ne conduit que les courants dus au dsquilibre entre les
phases.
Plus gnralement et indpendamment de leffet Joule, les effets des courants harmoniques
doivent tre pris en considration au niveau du dimensionnement de chaque appareil, voire de
chaque composant [TOL-96]. Prenons, titre dexemple le redresseur diodes de la figure
I.1.1, le choix des diodes constituant le pont ne peut tre guid par de simple considrations
sur la puissance active dlivre par le redresseur. En effet, ce type de convertisseur est connu
pour prlever des pics de courant sur le rseau lectrique. Une contrainte forte lors du choix
des diodes : elles doivent tre dimensionnes pour supporter ces pics de courant.
- 11 / 188 -
Chapitre 1 : Etat de lart
Le niveau des perturbations lectriques engendres par les appareils lectriques sur le rseau
basse tension peut tre qualifi laide du facteur de puissance, du spectre du courant absorb
et du spectre des ondes lectromagntiques rayonnes. Le facteur de puissance permet
destimer la puissance apparente transite par lappareil, qui peut tre plus grande que la
puissance active consomme (Annexe -A-).
- 12 / 188 -
Chapitre 1 : Etat de lart
Classe B : Outils portatifs (outils lectriques tenus la main pendant leur fonctionnement
normal et utiliss pendant une courte priode uniquement).
Tableau 1.1 : Limites pour les missions des courants harmoniques (EN 61000-3-2)
Bien au-del des contraintes normatives, nous avons vu que la minimisation des perturbations
harmoniques au sein dun dispositif lectrique vite un surdimensionnement, permet une
augmentation de sa fiabilit, une rduction de sa fatigue mcanique et une discrtion
- 13 / 188 -
Chapitre 1 : Etat de lart
acoustique Pour pallier les problmes gnrs par les convertisseurs et notamment la
production dharmoniques, plusieurs solutions peuvent tre envisages : les solutions
traditionnelles reposent sur des composants passifs (inductances, condensateurs,
transformateurs) [SO-96], et les nouvelles solutions reposent sur des dispositifs actifs (filtres
actifs, redresseurs absorption sinusodale des courants) [SO-96] [SAN-03].
Lutilisation dun couplage zig-zag au secondaire (figure I.1.2) permet dviter le transfert de
la composante homopolaire (harmonique de rang 3n, n entier naturel) au primaire. Elle est
bien adapte au cas de charges non linaires du type redresseurs de courant et gradateurs, dont
le contenu harmonique homopolaire est statistiquement trs lev [SO-96].
- 14 / 188 -
Chapitre 1 : Etat de lart
sinusode que le courant obtenu avec un seul redresseur. Grce ce montage, on repousse le
rang du premier harmonique prsent sur le rseau, ce qui constitue une solution la rduction
du taux de distorsion harmonique.
Il est possible dtendre cette proprit : en augmentant le nombre de redresseurs et le nombre
de secondaires du transformateur avec des dphasages relatifs appropris, dautres
harmoniques de courant peuvent tre limins. Ainsi, lassociation de 3 redresseurs connects
avec des phases dcales deux deux de 20 donne un courant 18 niveaux. C'est--dire que
seuls restent les harmoniques de rang 18n 1. Lassociation de 4 redresseurs connects avec
des phases dcales deux deux de 15 donne un courant 24 niveaux.
1
Is1 /Io [A]
0
-1
1
Is ' 1 / Io [A]
0
-1
1
Ip1 /Io [A]
0
-1
0.04 0.05 0.06 0.07 0.08 0.09 0.1
Temps [s]
Cette solution est intressante dans la mesure o elle permet de rduire le contenu spectral des
courants de charge et donc dobtenir des taux de distorsion harmonique en courant
relativement faibles. En revanche, ces montages sont complexes et onreux. Leur utilisation
- 15 / 188 -
Chapitre 1 : Etat de lart
se cantonne donc aux redresseurs de trs fortes puissances. A titre dexemple, llectrolyse de
laluminium qui fait appel des puissances de plusieurs MW utilise des montages allant
jusqu 72 phases [GUF-00].
Il sagit gnralement de circuits LC placs en parallle sur le rseau, comme illustr sur la
figure I.1.5, et accords sur la pulsation de lharmonique attnuer. Le filtre, qui prsente une
impdance minimale pour sa frquence daccord, se comporte comme un court-circuit pour
lharmonique considr. Simple, ce mode de filtrage prsente toutefois plusieurs
inconvnients. Il est slectif, ainsi si on souhaite attnuer plusieurs rangs dharmoniques,
autant de filtres sont ncessaires, ce qui peut poser des problmes en terme dencombrement.
Par ailleurs, le filtre est raccord limpdance variable du rseau et il se peut que cette
dernire forme avec lui un circuit rsonant : les harmoniques aux frquences voisines de la
frquence de rsonance risquent alors dtre amplifis.
Enfin dans le cas dune installation, du type rseau embarqu par exemple, soumise des
variations importantes de la frquence (de lordre de 10%), de la puissance de court-circuit et
o les frquences de certains harmoniques sont variables (fonction de la vitesse), les filtres
passifs sont mal adapts : pour tre toujours efficaces, ils doivent tre large bande, ce qui
rend leur ralisation dlicate et entrane une augmentation des pertes, de leur puissance, de
leur volume et de leur cot.
Le principe du filtrage actif est de compenser les harmoniques prsents sur le rseau en
injectant des harmoniques de mme amplitude mais en opposition de phase. On peut agir sur
les courants harmoniques (figure I.1.6) injects par une charge polluante (filtrage actif shunt
- 16 / 188 -
Chapitre 1 : Etat de lart
ou parallle) mais aussi sur les tensions harmoniques prsentes sur le rseau (figure I.1.7) afin
de protger une installation sensible (filtrage actif srie).
Le filtrage srie est destin protger les installations sensibles aux perturbations provenant
du rseau telles que les harmoniques en tension, les surtensions,. En revanche, le filtrage
srie ne permet pas de compenser les courants harmoniques consomms par la charge [GUF-
00].
charge
rseau
polluante
Le filtre actif parallle est un convertisseur statique connect au rseau en parallle sur la
charge polluante (figure I.1.8). Il est command pour absorber un courant compensant les
harmoniques prlevs au rseau par la charge polluante. Globalement, cest un courant
sinusodal qui circule sur le rseau. Plusieurs solutions technologiques base dlectronique
de puissance sont envisageables. Selon que la source de courant commande est construite
partir dun onduleur de tension (figure I.1.8a) ou dun commutateur de courant (figure I.1.8b),
le filtre actif est dit "structure tension" ou "structure courant". Dans tous les cas les filtres
actifs nchangent thoriquement aucune nergie active avec le rseau, si lon fait abstraction
des pertes dans les composants.
- 17 / 188 -
Chapitre 1 : Etat de lart
Les filtres actifs ont fait lobjet de divers travaux de recherche concernant les structures
possibles de convertisseurs [LOT-95], les commandes envisageables [GUF-00] ou
lamlioration de la robustesse et des performances dynamiques [GUF-00]. Des filtres actifs
sont dores et dj proposs sur le march par les industriels. Cest en effet une solution trs
confortable pour lutilisateur : le fait de connecter une "bote noire" sur un appareillage
lectrique le rend immdiatement conforme aux normes. Toutefois, ces filtres sont
dimensionns pour une puissance donne, et doivent tre rtudis si la charge varie.
rseau rseau
Dans le cas du filtrage actif parallle, pour compenser exactement les harmoniques de courant
gnrs par un redresseur de 6 kW un filtre actif de 6 kVA est ncessaire. Ce filtre peut tre
constitu de nimporte quel convertisseur statique dont le courant dentre est contrlable.
Thoriquement un filtre actif ne consomme que les puissances ractive et dformante.
Toutefois, vue les dispositifs de stockage dnergie (condensateurs ou bobines) et les pics en
courant (dus aux harmoniques compenser), pratiquement le filtre actif est dimensionn pour
une puissance quivalente celle de la charge compenser.
On peut donc se demander sil nest pas plus avantageux dutiliser directement un
convertisseur AC/DC dont le courant dentre est asservi une sinusode [CHE-00] [BEN-
98]. Ces convertisseurs statiques absorption quasi sinusodale de courant sont appels PFC
(Power Factor Corrector) : le courant dentre est asservi une sinusode de rfrence par une
technique de dcoupage. Lamplitude de la rfrence est donne par une boucle de rgulation
- 18 / 188 -
Chapitre 1 : Etat de lart
La figure I.2.1 prsente la structure la plus simple des redresseurs prlvement sinusodal,
ne ncessitant quun seul interrupteur commandable.
- 19 / 188 -
Chapitre 1 : Etat de lart
Ce montage conomique permet le dcoupage du courant dentre de sorte quil ny ait pas
dharmoniques en basse frquence [PRA-89] [KOL-93]. Le courant absorb prsente
cependant des ondulations leves la frquence de dcoupage du hacheur, et ncessite donc
la prsence dun filtre passif (LC) encombrant en entre. La tension de ltage continu doit
tre asservie une valeur suprieure au maximum de la tension compose en entre afin
dassurer le contrle du courant absorb sur le rseau. Ce convertisseur ne convient que pour
une gamme de puissance infrieure 8 kW [JAN-98].
Cette structure (figure I.2.2) fonctionne dans des gammes de puissance plus importantes que
la prcdente [MOH-92]. Le principe de fonctionnement est diffrent, car la commande est
inspire du principe du filtrage actif. La commande du hacheur assure deux fonctions :
lasservissement de la tension de sortie et linjection dharmoniques en opposition de phase
avec ceux gnrs par le pont redresseur PD3 dans le rseau, via des filtres LC [MOH-92].
Ceci permet de rduire la puissance dformante. Les filtres LC sont accords pour
lharmonique de rang 3. Dans [SAL-96] une approche utilisant la commutation des
interrupteurs au lieu du filtre LC, est prsente.
Le circuit de puissance de cette structure est illustr sur la figure I.2.3. Cette structure contient
trois interrupteurs commandables. Avec un systme de commande relativement simple [KOL-
94] [KOL-96] on obtient un facteur de puissance proche de 1. La difficult de cette structure
- 20 / 188 -
Chapitre 1 : Etat de lart
rside dans lasservissement de la tension continue, car lquilibre entre les tensions aux
bornes des deux condensateurs doit tre soigneusement maintenu [OOI-85] [WU-88] [WU-
91]. En plus des deux sources de tension continues, cette structure offre aussi la possibilit de
contrler le dphasage courant tension, sans dpasser toutefois 30 [KOL-94], ce qui reste
insuffisant pour assurer la bidirectionnalit en puissance. Par contre si la bidirectionnalit en
courant nest pas ncessaire, ce convertisseur est trs avantageux.
- 21 / 188 -
Chapitre 1 : Etat de lart
Dans ce convertisseur (figure I.2.6), le courant dentre ne peut tre prlev la source que si
la tension dentre est suprieure la tension de sortie. En consquence, le courant dentre
prlev aura de fortes composantes harmoniques de rangs faibles. De plus, le courant dentre
est hach ce qui explique lutilisation du filtre LC du ct rseau [ISM-92]. Si on dsire un
fonctionnement continu tout au long de la priode rseau, ce convertisseur se trouve cart.
- 22 / 188 -
Chapitre 1 : Etat de lart
La figure I.2.7 illustre le circuit de puissance dun redresseur Buck. Les travaux de [MAL-87]
ont dmontr que labsorption sinusodale du courant peut tre obtenue facilement par ce
redresseur, et que cet tat de fonctionnement nest assur que lorsque le dphasage courant
tension dentre est compris entre -30 et +30. Pour = 0, la tension de sortie est infrieure
de 10% la tension dalimentation, et ce pourcentage augmente lorsque le dphasage
sloigne de 0. Baumann [BAU-01] a approfondi ltude de ce convertisseur en valuant
exprimentalement les pertes dues aux commutations ainsi que sa fiabilit.
Les circuits de puissance correspondant cette structure sont illustrs sur la figure I.2.8. Ce
convertisseur peut tre considr comme une extension du redresseur thyristors [KAT-79]
[BUS-82] [WIE-84] [ZIO-85]. Chaque thyristor est remplac par un IGBT avec une diode en
srie. Cette structure assure le prlvement sinusodal. Cependant, le courant absorb prsente
des ondulations leves la frquence de dcoupage, et ncessite la prsence dun filtre passif
(LC) en entre. Linconvnient majeur de ce convertisseur est celui des pertes en conduction
dues la prsence des diodes [WAL-01].
- 23 / 188 -
Chapitre 1 : Etat de lart
Les figures (I.2.9 et I.2.10) prsentent deux structures amliores du redresseur lvateur
VIENNA. Afin dassurer un fonctionnement abaisseur/lvateur, des circuits de stockage
intermdiaires sont ajouts [PAN-99] [SHI-00]. Comme le redresseur VIENNA, ces deux
structures ne permettent pas davoir une rversibilit en courant.
- 24 / 188 -
Chapitre 1 : Etat de lart
Le cas de la figure I.2.9 est inspir du concept du hacheur accumulation capacitive [PAN-
99]. Le second convertisseur (figure I.2.10), est inspir de la structure SEPIC (Single Ended
Primary Inductance Converter) [KOL-97] [SHI-00].
Dans [KOL-97] le principe de fonctionnement ainsi que des rsultats de simulation ont t
prsents. Toutefois les auteurs nont pas approfondi leur travail avec une tude
exprimentale et/ou une tude de la fiabilit comme dans le cas du redresseur Buck. Cela peut
sexpliquer, notre avis, par la complexit engendre par le nombre important de composants
passifs ajouts la structure de base.
- 25 / 188 -
Chapitre 1 : Etat de lart
Les structures rversibles SEPIC (figure I.2.11) et Cuk (figure I.2.12) ont t tudies
respectivement dans les travaux de [SHI-02] et [KIK-02]. Le principe de la commande est
exactement le mme dans les deux structures. La seule diffrence significative concerne les
circuits logiques permettant de gnrer la commande des sept interrupteurs partir de la MLI
vectorielle standard (six interrupteurs). Dans le cadre de ce mmoire, notre travail porte sur la
conception dune nouvelle technique de commande de la structure Cuk rversible. Le principe
de notre mthode pourra facilement tre adapt la structure SEPIC rversible.
- 26 / 188 -
Chapitre 1 : Etat de lart
dalimentation. Cette commande tant ralise dans la majorit des cas avec des circuits
numriques sensibles, les perturbations provenant du rseau peuvent entraner le
dysfonctionnement de ces circuits, voire mme leur destruction.
En plus de lasservissement des courants alternatifs, la commande doit aussi assurer une
tension continue en sortie quelle que soit la nature de la charge. Cette fonction ncessite
lasservissement en temps rel de la tension de sortie en tenant compte de deux paramtres : la
tension de sortie et le courant dans la charge [LEE-00]. La connaissance de ces deux
paramtres exige une augmentation du nombre des capteurs afin de mesurer ces deux signaux.
En consquence, cette augmentation du nombre de capteurs engendre une augmentation
considrable du cot de la ralisation. De plus, daprs [B-E380], le nombre de capteurs et de
composants utiliss est inversement proportionnel la fiabilit des convertisseurs statiques.
Cest pour cette raison que depuis 1997, plusieurs travaux ciblent la rduction du nombre de
capteurs dans les convertisseurs prlvement sinusodal, plus particulirement la structure
Boost rversible en courant. Les travaux les plus rcents [LEE-02] ont permis de raliser un
nouveau systme de commande, pour la structure Boost, nutilisant aucun capteur du ct
alternatif.
- 27 / 188 -
Chapitre 1 : Etat de lart
La modulation MLI est fortement inspire des techniques utilises dans la transmission de
linformation. Elle consiste moduler une onde porteuse (issue dun dcoupage frquence
Fd) par une onde modulante (frquence Fs) frquence beaucoup plus basse.
Il existe nanmoins une diffrence fondamentale puisque linverse de la transposition
spectrale recherche en transmission, une composante basse frquence est gnre partir de
londe porteuse. Pour obtenir ce rsultat, cest la valeur moyenne de cette onde porteuse
dfinie sur la priode Td qui est module tandis quelle reste nulle en transmission. Le
principe et les rsultats purement qualitatifs de cette technique de modulation sont donns
dans la figure I.3.1.
Td Ts
Fs
Fd
La modulation MLI peut tre dcompose en deux grandes familles. La premire famille
rassemble les stratgies autorisant des modifications de la largeur dimpulsion imposes par
les variations dune grandeur instantane gnre par les organes de contrle, lchelle de la
priode de dcoupage Td. On parle alors de MLI instantanes, qui peuvent tre locales
(commande dune cellule ou interrupteur) ou globales (gestion vectorielle de la commande de
plusieurs cellules). Elles correspondent gnralement des cas pour lesquels la frquence de
dcoupage Fd est grande devant la frquence fondamentale Fs [F-D3176].
La seconde famille inclut les stratgies dans lesquelles les formes de londe dcoupe sur une
priode fondamentale sont calcules a priori pour minimiser les composantes harmoniques
parasites, valeur donne de la composante fondamentale. Les rsultats de ces calculs sont
mmoriss puis utiliss selon les besoins en cours de fonctionnement. On ne peut changer de
configuration qu lchelle de la priode de modulation Ts (Ts/4 en monophas, Ts/6 en
triphas). Le principe de cette stratgie empche d'agir en temps rel l'chelle de la priode
- 28 / 188 -
Chapitre 1 : Etat de lart
de dcoupage et l'on perd donc le caractre quasi-instantan ncessaire pour la commande des
convertisseurs prlvement sinusodal des courants.
Comme lobjectif est de contrler en temps rel le courant dentre, seule la premire famille
de modulation prsente un intrt dans la commande des convertisseurs prlvement
sinusodal. En particulier en triphas, la MLI vectorielle reste une stratgie incontournable
dans la conception de ce type de convertisseur [HOL-92].
E 2/3
v 3 v 2
(010) (110)
II
III I
v 4 v 1
(011) (100)
VI
IV
V
v 5 v 6
(001) (101)
Si lon sappuie sur le schma de la figure I.3.2a, il apparat clairement que lobjectif final est
de contrler le systme triphas {vaN,vbN,vcN} aux bornes dune charge triphase quilibre.
On peut alors se ramener un vecteur diphas, image du systme triphas, travers la
transformation de Concordia inverse.
La commande en MLI vectorielle de londuleur consiste le piloter dans le repre , , dans
lequel il peut gnrer six vecteurs dont la reprsentation est donne sur la figure I.3.2b. Si lon
utilise directement ces tats, on se place dans le cadre limit de la commande pleine onde.
Pour contrler le vecteur (v, v) en phase et en amplitude, il faut ajouter un aspect temporel,
cest dire la modulation de largeur dimpulsion. On fabrique un vecteur moyen l'chelle de
- 29 / 188 -
Chapitre 1 : Etat de lart
la priode de dcoupage qui rsulte du moyennage temporel des deux vecteurs les plus
proches du diagramme (figure I.3.3).
v 2
(100) (110)
(110)
fm1
fm2 v
(R2 Td) v2
R2 Td/2
fm3
v 1
(R1 Td) v1 (100)
R1 Td/2 Td/2 Td
La modulation MLI dcoupage alatoire RPWM (Random Pulse Witdh Modulation) est
lune des solutions pour minimiser le problme des harmoniques en haute frquence. Dans la
littrature, cette technique a t utilise dans la commande des hacheurs [TS1-00] [TS2-00]
ainsi que dans les onduleurs triphass [LIA-00]. Le dcoupage alatoire permet dattnuer la
densit spectrale en haute frquence du signal modul [TSE-02] [BEC-99]. Au lieu davoir
des composantes spectrales en chaque multiple de la frquence de dcoupage, on obtient des
composantes tales sur tout le spectre de frquence avec des amplitudes aux frquences
centrales attnues.
Pour expliquer le principe de la modulation dcoupage alatoire, nous allons revenir sur la
notion du rapport cyclique. Rappelons que la fonction de modulation dune commande scrit
sous la forme suivante :
- 30 / 188 -
Chapitre 1 : Etat de lart
1, si Interrupteur ferm
f m (t) =
fm(t) d.T
0, si Interrupteur ouvert
(3)
t0 t0+T
et le rapport cyclique au sein dune priode de dcoupage T est :
T t 0
d=
1 t0 +T
f m (t).dt (4)
Daprs (4) nous remarquons que si la priode de dcoupage T varie, nous pouvons obtenir
le mme rapport cyclique en faisant varier la dure de fermeture de la cellule de commutation.
Dans le mme contexte, dans le cas o la priode T est constante, on obtient ce mme
rapport cyclique en changeant la position de limpulsion de commande sur lintervalle
[t0, t0+T], [BEC-99].
Donc le principe de la RPWM, consiste faire varier alatoirement lun de ces deux
paramtres (la priode de dcoupage ou la position de limpulsion) tout en adaptant le second
paramtre afin dobtenir le mme rapport cyclique d.
Pour llaboration de la modulation MLI commande alatoire, trois concepts existent dans la
littrature : Random Switching (RS), Random Pulse Position (RPP) et Randomized Carrier
Frequency Modulation (RCFM).
- 31 / 188 -
Chapitre 1 : Etat de lart
Valeur moyenne
MLI
t
0 T/2 T
Valeur moyenne
RS
t
0 T/5 2T/5 3T/5 4T/5 T
Dans le deuxime concept RPP [TRZ-94] [BEC-99], la priode de dcoupage est fixe et la
position de limpulsion de commande se dplace alatoirement sur chaque priode (figure
I.3.5). Le fait de varier la position de limpulsion de commande attnue le spectre rsultant de
cette technique de modulation. Toutefois, son effet sur les multiples de la frquence de
dcoupage reste trs limit.
MLI
d1.T d2.T
RPP
t
0 T 2T
Figure I.3.5 : Random Pulse Position (RPP)
Les deux techniques RS et RPP sont utiliss dans les convertisseurs DC/DC pour la
commande des machines courant continu [TS1-00] afin dinjecter un dcoupage alatoire
dans la tension de sortie de ces convertisseurs. Cela permet de rduire le bruit acoustique des
machines.
Le troisime concept RCFM [TRZ-94] [BEC-99] consiste utiliser une porteuse frquence
alatoire. Le signal triangulaire (la porteuse) est frquence fixe pendant une priode de
dcoupage et seule la priode de dcoupage varie alatoirement (figure I.3.6b). Le principe de
cette mthode consiste faire varier alatoirement la priode de dcoupage tout en maintenant
le mme rapport cyclique di de la MLI (figure I.3.6a).
- 32 / 188 -
Chapitre 1 : Etat de lart
MLI
t
T 2.T
RCFM
d1 d2
t
0 Tr1 Tr1+Tr2
La figure I.3.6 illustre le principe de cette mthode, le signal de modulation est un signal
triangulaire frquence alatoire variant entre fmin = 500 Hz et fmax = 1,5 kHz.
10 10
5 5
[V]
[V]
0 0
-5 -5
-10 -10
0.2 0.202 0.204 0.206 0.208 0.21 0.212 0.214 0.216 0.218 0.22 0.2 0.202 0.204 0.206 0.208 0.21 0.212 0.214 0.216 0.218 0.22
10 10
5 5
[V]
[V]
0 0
-5 -5
-10 -10
0.2 0.202 0.204 0.206 0.208 0.21 0.212 0.214 0.216 0.218 0.22 0.2 0.202 0.204 0.206 0.208 0.21 0.212 0.214 0.216 0.218 0.22
[s] [s]
8 8
7 7
6 6
5 5
[V]
[V]
4 4
3 3
2 2
1 1
0 0
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
La figure I.3.7 compare la modulation RCFM avec la modulation MLI. Nous constatons que
dans le cas de la RCFM, la densit spectrale du signal modul, est attnue et tale sur tout le
- 33 / 188 -
Chapitre 1 : Etat de lart
spectre de frquence. Cette technique est la plus adquate pour notre tude car, dune part, le
systme garde exactement les mmes vecteurs de commande et, dautre part, parce que pour
un fonctionnement en temps rel, la RCFM se comporte de la mme manire que la MLI. Une
tude plus dtaille de cette technique de modulation sera labore dans le chapitre III.
Notre tude porte sur le convertisseur Cuk rversible dcrit dans le paragraphe II.3.2. Les
fonctionnalits de ce convertisseur sont les suivantes : abaisseur/lvateur de la tension de
sortie, absorption sinusodale du courant rseau et bidirectionnalit en puissance. Le schma
de principe du convertisseur est reprsent dans la figure I.4.1.
Circuit de dcouplage
g1 g3 g5
g7
g4 g6 g2
Notre tude sinscrit dans la continuit du travail conu par J. KIKUCHI [KIK-02] de
lUniversit de Madison (USA). Daprs [KIK-02], quatre capteurs sont utiliss dans la
structure existante (figure I.4.1). Deux capteurs mesurent les courants de phase afin que la
commande puisse corriger la dformation du courant absorb. Le troisime capteur permet de
mesurer la tension de sortie afin que le systme dasservissement assure sa stabilit. Enfin, le
rle du dernier capteur est de mesurer le courant consomm par la charge. Lobjectif de ce
capteur est de tenir compte de la nature de la charge, ce qui permet au systme de commande
de suivre rapidement les variations du courant de charge. Ce dernier capteur est galement
utilis dans le convertisseur Boost triphas pour les mmes objectifs [SIN-00].
Le circuit de stockage intermdiaire appel aussi circuit de dcouplage, assure deux
fonctionnalits : le fonctionnement en abaisseur/lvateur de la tension de sortie, et la
rversibilit en puissance. La prsentation de toutes les fonctionnalits du convertisseur sera
labore dans le chapitre 2.
- 34 / 188 -
Chapitre 1 : Etat de lart
g1 g3 g5
g7
g4 g6 g2
Pour atteindre ces objectifs, nous avons amlior la commande de faon ne ncessiter que
deux capteurs de mesure (figure I.4.2), et utilisant une technique de modulation dcoupage
alatoire permettant dattnuer la densit spectrale en haute frquence sans utiliser les circuits
de filtrage (paragraphe III.3). Ce nouveau positionnement des capteurs de mesure engendre
- 35 / 188 -
Chapitre 1 : Etat de lart
une reconstruction interne des courants de phase et de la tension de sortie partir des
diffrents signaux mesurs.
La reconstruction des signaux partir des vecteurs de modulation (paragraphe III.2), tout
en tenant compte des variations alatoires de la frquence de dcoupage.
Lutilisation dobservateur des variables dtats afin de tenir compte des retards dus
lacquisition (conversion analogique/numrique) et aux temps de calcul ncessaires
lexcution des instructions, pour la reconstruction des signaux de faon synchrone.
Afin de raliser cette nouvelle commande, la premire tape consiste en une phase de
simulation. Lobjectif est de vrifier tout dabord la faisabilit et les limites de cette approche
avant de se lancer dans la ralisation du convertisseur. Pour les simulations, nous utilisons
MATLAB/SIMULINK avec la bibliothque Power System Toolbox. Le choix de cet outil de
simulation est justifi par sa flexibilit ainsi que par les grands avantages quil offre dans les
simulations des systmes dynamiques grce des outils de programmation adapts tels que
les S-Functions.
Pour la ralisation du nouveau systme de commande, nous utilisons le processeur de signal
TMS320LF2407A. Ce DSP (Digital Signal Processor) est adapt pour la commande des
convertisseurs de puissance en particulier les convertisseurs utilisant la modulation de largeur
dimpulsion. La construction des signaux de commande des IGBT partir des fonctions de
modulation est ralise par des circuits logiques programmables de type ISPGAL.
- 36 / 188 -
Chapitre 1 : Etat de lart
Conclusion
- 37 / 188 -
Chapitre 2 : Prsentation du convertisseur
Conclusion ........................................................................................................................................................... 73
- 38 / 188 -
Chapitre 2 : Prsentation du convertisseur
Introduction
Dans le premier chapitre, nous avons prsent le convertisseur Cuk rversible qui fera l'objet
de notre tude. La conception dun tel convertisseur ncessite une connaissance dtaille de
son modle thorique qui fait appel deux aspects :
- L'tat statique des interrupteurs : il dfinit lvolution dans le temps de la
configuration ON ou OFF de chaque semi-conducteur. Cette volution est reprsente
par des fonctions de modulation associes chaque interrupteur.
- La mise en quation de chaque configuration du montage de puissance : elle permet de
rsoudre squentiellement plusieurs quations d'tat et de dduire l'volution
temporelle des variables d'tat associes. Le dimensionnement des diffrents lments
du convertisseur s'en dduit aisment.
Le modle du convertisseur sera notre base de travail pour l'tude et l'analyse des variables
d'tat tant pour la commande MLI classique qu'alatoire.
Ce chapitre consacr la modlisation est scind en trois parties : La premire partie concerne
la description du principe de fonctionnement du convertisseur Cuk rversible et llaboration
des fonctions de modulation rgissant ce convertisseur. Dans la deuxime partie, nous
prsenterons le dimensionnement des circuits de stockage dnergie et de filtrage du circuit de
puissance tandis que la stratgie de commande sera traite dans la dernire partie de ce
chapitre.
I. Description du convertisseur
La figure 2.1.1 prsente la structure de puissance du redresseur Cuk rversible. linverse
des structures traditionnelles, chaque IGBT est command indpendamment pour assurer le
redressement de la tension de sortie, le fonctionnement en abaisseur lvateur de la tension de
sortie, le prlvement sinusodal du courant de phase et la rversibilit en puissance.
- 39 / 188 -
Chapitre 2 : Prsentation du convertisseur
De plus, durant certaines phases de fonctionnement, les bras du pont triphas sont
volontairement court-circuits afin dassurer le transfert dnergie du circuit de stockage
intermdiaire vers la charge et vice-versa. Dans la suite de ce mmoire, nous allons appeler ce
circuit de stockage intermdiaire circuit de dcouplage. Ce dernier est constitu de
composants ractifs (capacit Cc et inductance Ldc) et dune diode monte en antiparallle
avec un IGBT. Dans ce paragraphe, nous prsentons le rle de ce circuit et ltat du
convertisseur dans les modes de fonctionnement redresseur MLI et onduleur non autonome.
Pour prsenter le fonctionnement en mode redresseur MLI, nous avons choisi une priode de
dcoupage qui dbute l'un instant t1 o ia > 0 > ib > ic et vaca > vacb > vacc. La figure 2.1.2
illustre les impulsions de commande des sept IGBT, rsultantes de la comparaison des
tensions vac(a,b,c) avec la porteuse. Ces impulsions sont caractrises par trois zones de
modulation x, y et z. Les zones x et y correspondent aux vecteurs de commande actifs
au sein de chaque priode de dcoupage, et les zones z correspondent aux vecteurs nuls.
Dans x et y la commande du convertisseur Cuk est identique celle du redresseur Boost
(figure 2.1.2). Dans ces deux intervalles, la diode du circuit de dcouplage est passante, ce qui
permet la charge du condensateur de dcouplage Cc et le maintien de la tension dans Cdc par le
courant iLdc (figure 2.1.3).
Dans la zone de commande z lnergie stocke dans Cc est transfre vers la charge, ce qui
ncessite de court-circuiter au minimum un bras du pont triphas. Trois cas de commande
sont possibles (figure 2.1.4). Dans les figures 2.1.2(a) et 2.1.4(a) le bras du pont
correspondant au courant de phase intermdiaire ib (ic < ib < ia) est court-circuit. Dans ce cas,
les IGBT 3 et 6 du bras de la phase b doivent supporter respectivement les courants
- 40 / 188 -
Chapitre 2 : Prsentation du convertisseur
(ia + iLdc) et (ic + iLdc), et la frquence de commutation moyenne des cellules du pont triphas
est infrieure la frquence de dcoupage (2.1).
2 2
f = .2.f dcoupage = .f dcoupage (2.1)
6 3
Vdcoupage
Vaca.sin(wt - )
Vaca Vacb Vacc
Vacb.sin(wt - - 2/3)
Vacc.sin(wt - + 2/3)
z x y z y x z
g1
g4
g3
Un seul bras du pont triphas
g6
court-circuit dans la zone z
g5 (a)
g2
g7
g1
g4
g3
Deux bras du pont triphas
g6
court-circuits dans la zone z
g5 (b)
g2
g7
g1
g4
g3
Trois bras du pont triphas
g6
court-circuits dans la zone z
g5 (c)
g2
g7
Le deuxime cas est le court-circuit de deux bras du pont triphas (figures 2.1.4(b) et
2.1.2(b)). Ces deux bras correspondent aux phases dont les courants sont de mme signe.
Dans lexemple de la figure 2.1.4(b), les deux bras concerns sont b et c. Les deux IGBT 3
et 5 sont parcourus par les courants respectifs (ib + iLdc/2) et (ic + iLdc/2) et les IGBT 2 et 6
- 41 / 188 -
Chapitre 2 : Prsentation du convertisseur
sont parcourus par le courant iLdc/2. La frquence de commutation moyenne des cellules du
pont triphas est gale la frquence de dcoupage (2.2).
3
f = .2.f dcoupage = f dcoupage (2.2)
6
iLdc
(a) Intervalle x
iLdc
(b) Intervalle y
Figure 2.1.3 : Etat du convertisseur Cuk dans les zones x et y (mode redresseur MLI)
Dans les figures 2.1.2(c) et 2.1.4(c), trois bras du pont triphas sont court-circuits. Dans ce
cas, le courant supplmentaire iLdc est support par les trois bras du pont, ce qui reprsente le
cas le plus avantageux de ce point de vu par rapport aux deux premiers cas. Cependant, vu
que quatre cellules commutent deux fois pendant une priode de dcoupage, la frquence de
commutation moyenne du pont triphas augmente (2.3).
4 4
f = .2.f dcoupage = .f dcoupage (2.3)
6 3
iLdc
- 42 / 188 -
Chapitre 2 : Prsentation du convertisseur
iLdc
iLdc
En rsum, le meilleur compromis est de court-circuiter deux bras du pont triphas, car dune
part le courant supplmentaire (iLdc) qui se superpose aux courant de phase est support par
deux bras du pont, et dautre part la frquence de commutation moyenne du pont est gale la
frquence de dcoupage. La symtrie du circuit de dcouplage et des cellules de commutation
identiques dans le pont triphas sont indispensables pour le bon fonctionnement du
convertisseur. Nous utilisons deux condensateurs de capacit 2Cc et deux bobines
dinductance Ldc/2 pour obtenir lquivalent de la capacit et de linductance de dcouplage
Cc et Ldc ncessaires.
- 43 / 188 -
Chapitre 2 : Prsentation du convertisseur
avec la porteuse. Ces impulsions sont caractrises par quatre zones de modulation z0, x,
y et z1.
Vdcoupage
Vaca.sin(wt - )
Vacb.sin(wt - - 2/3)
Vacc.sin(wt - + 2/3)
z0 x y z1 y x z0
g1
g4
g3
g6
g5
g2
g7
Dans les zones x et y, linterrupteur g7 est ferm, et la commande du pont triphas est
similaire celle du convertisseur Boost (figure 2.1.7). Pendant ces deux phases de
fonctionnement, le condensateur de dcouplage fournit lnergie au rseau (nergie
accumule durant les intervalles z0 et z1).
(a) Intervalle x
(b) Intervalle y
Figure 2.1.6 : Etat du convertisseur Cuk dans les zones x et y (mode onduleur non
autonome)
- 44 / 188 -
Chapitre 2 : Prsentation du convertisseur
LIGBT g7 nest command que dans le mode de fonctionnement onduleur non autonome, car
contrairement au mode redresseur, la conduction naturelle de la diode D7 nest pas assure. La
charge du condensateur de dcouplage est effectue dans les intervalles de commande z0 et
z1 via le pont triphas. Dans ce cas le pont triphas est court-circuit naturellement par la
conduction des diodes en antiparallle dans le pont triphas. La conduction des diodes dpend
de la comparaison du courant iLdc avec le(s) courant(s) de phase de mme signe.
iLdc
(a) Intervalle z0
iLdc
(b) Intervalle z1
Par exemple, lorsque iLdc < ia (figure 2.1.7), dans lintervalle de commande z0 la conduction
des deux diodes D3 et D5 engendre le court-circuit de deux bras du pont triphas. Dans
lintervalle z1 la conduction de la diode D4 permet la charge du condensateur Cc via le bras
de la phase a. Le deuxime exemple de la figure 2.1.8 concerne le cas iLdc > ia, o les trois
bras sont court-circuits par la conduction du pont triphas diodes dans les deux intervalles
z0 et z1.
- 45 / 188 -
Chapitre 2 : Prsentation du convertisseur
iLdc
Aprs cette description du principe de fonctionnement, nous allons laborer la commande des
interrupteurs du convertisseur Cuk rversible par une approche gnrique . Le paragraphe qui
suit est consacr la prsentation des diffrents variables rgissant les signaux de commande
des quatre bras du convertisseur.
v aca - v bn = H a v cc
v - v = H v
acb bn
v acc - v bn = H c v cc
b cc
(2.4)
v 7 = H z v cc
En se basant sur la comparaison directe des tensions de rfrence vac(a,b,c) avec la porteuse,
nous avons dfini les quatre fonctions de modulation du convertisseur. Nous notons Ha_bst,
Hb_bst et Hc_bst les fonctions de modulation rsultantes de cette comparaison directe
(figure 2.2.1), et Hmax_bst et Hmin_bst la fonction H(a,b,c)_bst correspondante la tension de
rfrence vac(a,b,c) respectivement la plus grande et la plus petite. Les relations dans (2.5)
rcapitulent la correspondance entre les fonctions Ha_bst, Hb_bst et Hc_bst avec Ha, Hb, Hc et Hz.
La figure 2.2.1 prsente lallure des quatre fonctions de modulation au sein dune priode de
dcoupage lorsque vaca > vacb > vacc.
- 46 / 188 -
Chapitre 2 : Prsentation du convertisseur
H z
= H max_bst + H min_bst
H a
= H a_bst - H min_bst
H b
(2.5)
H c
= H b_bst - H min_bst
= H c_bst - H min_bst
Vaca.sin(wt - )
Vacb.sin(wt - - 2/3)
Vacc.sin(wt - + 2/3)
Ha_bst
Fonctions de modulation
Hb_bst du convertisseur Boost
triphas
Hc_bst
Ha
Hb Fonctions de modulation
du convertisseur Ck
Hc rversible
Hz
Figure 2.2.1 : Comparaison des fonctions de modulation du convertisseur Cuk avec les
fonctions de modulation du convertisseur Boost triphas
Dans le convertisseur Boost, Lorsque le vecteur de commande (Ha_bst, Hb_bst, Hc_bst)t est
ltat (1, 1, 1)t ou (0, 0, 0)t il est appel vecteur de commande zro, et dans les autres cas il est
appel vecteur de commande actif. Par correspondance avec le convertisseur Boost, le vecteur
de commande du convertisseur Cuk (Ha, Hb, Hc, Hz) est actif lorsque Hz = 0 quelque soit les
valeurs des autres fonctions de commande (figure 2.2.1). Le vecteur de commande zro
correspond Hz = 1.
- 47 / 188 -
Chapitre 2 : Prsentation du convertisseur
OUI
signe_i > 0
g 1 := H a_bst
g 2 := H c_bst
g 3 := H b_bst H a := H a_bst - H min
g 4 := H a_bst H b := H b_bst - H min
g 5 := H c_bst H c := H c_bst - H min
Afin de gnrer les impulsions de commande partir des fonctions de modulation, nous avons
labor une procdure de fonctionnement gnrique (figure 2.2.2). Les paramtres dentre de
lalgorithme sont les suivants :
- 48 / 188 -
Chapitre 2 : Prsentation du convertisseur
Les signaux de sortie sont les impulsions de commande des sept IGBT correspondant
g1g7. Les fonctions intermdiaires calcules, sont dans un premier temps les deux fonctions
Hmax_bst et Hmin_bst, ensuite les quatre fonctions de modulation Ha, Hb, Hc et Hz, pour enfin
gnrer les impulsions de commande spcifiques chaque tat du convertisseur.
Lalgorithme de la figure 2.2.2 est implant dans des circuits logiques programmables de type
ISPGAL en utilisant le langage de description matriel ABEL-HDL.
Dans cette modlisation, nous considrons les cellules de commutation parfaites, tout comme
les inductances, les capacits et la tension rseau. Nous ngligeons les variations du courant
de phase dues au dcoupage par rapport aux variations du fondamental (50 Hz), et nous
considrons que le convertisseur dbite sur une charge purement rsistive Rcharge. Le
fonctionnement du convertisseur Cuk est rgi par les six quations de (2.6).
- 49 / 188 -
Chapitre 2 : Prsentation du convertisseur
di a
L ac dt = -R aci a - H a v cc + v a - v bn
L di b = -R i - H v + v - v
ac dt
ac b b cc b bn
L ac di c = -R aci b - H c v cc + v c - v bn
dv
dt
Cc cc = H a i a + H bi b + H ci c + H z i Ldc
(2.6)
di L
dt
L dc dc = H z v cc - R dci Ldc - v dc
dv dc
dt
Cdc
v
= i Ldc - dc
dt R charge
avec Rac et Rdc sont respectivement les rsistances internes des inductances Lac et Ldc.
Le rseau triphas tant quilibr (ia + ib + ic= 0 et va + vb + vc= 0), les trois premires
quations dans (2.6) permettent de dterminer la tension vbn en fonction de Ha, Hb, Hc et vcc.
v bn = - vcc ( H a + H b + H c )
1
(2.7)
3
Ainsi, en remplaant vbn dans (2.6), les quations de cette dernire peuvent tre exprimes
sous la forme matricielle suivante :
R ac
- L ( -2H a + H b + H c )
va
L
1
ac
0 0 0 0
3L ac
vb
( H a - 2H b + H c )
ac
ia 0 ia
R ac 1
L ac
- 0 0 0
L ac 3L ac
ib ib vc
0 ( H a + H b - 2H c )
R ac 1
(2.8)
d ic i c + L ac
0 - 0 0
dt v cc H a v
L ac 3L ac
=
cc 0
Hb Hc Hz
i L dc C c i L dc
0 - 0
Cc Cc Cc
v dc 0 v dc 0
Hz R dc 1
0 0 - -
L dc L dc L dc
0 0
1 1
0 0 0 -
C dc C dc R charge
Dans la suite de cette tude nous considrons le cas particulier vaca > vacb > vacc, c'est dire
que vaca est la tension la plus positive, ce test est vrifi sur un sixime de la priode
fondamentale. L'extrapolation toute la priode ncessite une adaptation des indices des
fonctions de modulation :
H z = H a_bst + H c_bst
H a = H a_bst - H c_bst
H b = H b_bst - H c_bst
(2.9)
H c = 0
- 50 / 188 -
Chapitre 2 : Prsentation du convertisseur
Dans lhypothse o les variations dues au dcoupage (Fd de l'ordre de quelques kHz) sont
ngligeables devant les variations du signal basse frquence (50 Hz), et en ngligeant les
ondulations de la tension vcc, nous pouvons remplacer les fonctions Ha_bst, Hb_bst et Hc_bst par
les rapports cycliques moyens da, db et dc. Ainsi :
H a_bst = d a = m 1+ sin ( wt - )
1
2
1 2
H b_bst = d b = m 1+ sin wt - -
2 3
(2.10)
1 2
H c_bst = d c = m 1+ sin wt - +
2 3
avec m = 2.Vac_amplitude / vcc et le dphasage entre la tension secteur v(a,b,c) et les tensions
vac(a,b,c). Daprs (2.5), les fonctions Ha, Hb, Hc et Hz deviennent sur un sixime de priode :
H a = m.sin wt - -
3
6
2
H b = m.sin wt - -
3
2
H = 0
2 (2.11)
c
H z = 1- m.sin wt - -
3
2 6
La figure 2.2.3 illustre lvolution de ces fonctions de modulation durant une priode de
20 ms.
Ha
Hb
Hc
Hz
Temps [s]
Figure 2.2.3 : Evolution des fonctions de modulation du convertisseur Cuk sur 20ms
- 51 / 188 -
Chapitre 2 : Prsentation du convertisseur
En appliquant la transformation de Park (ia, ib, ic), (va, vb, vc) et (Ha, Hb, Hc) nous exprimons
la relation (2.8) dans le rfrentiel synchrone d-q-0. La matrice de transformation T dun
rfrentiel stationnaire (a, b, c) un rfrentiel d-q-0 est exprime dans (2.12).
2
sin ( wt ) sin wt - sin wt +
2
3 3
d a a
2 2
q = T b = 3 cos ( wt ) cos wt - 3 cos wt + 3 b
2
0 c
(2.12)
c 1
1 1
2 2 2
Lquation (2.8) exprime dans le rfrentiel d-q-0 scrit alors sous la forme :
R ac vd
-L
Hd
ac Lac
-w 0 - 0 0
vq
Lac
d w id
R ac Hq
i - 0 - 0 0
iq q v0
L ac Lac L ac
0
R ac i
d i0 i 0 Lac
0 - 0 0 0
dt v cc 3 H d v +
Lac
cc 0
= (2.13)
3 Hq
iL
3 H0 H
i Ldc 2 Cc
0 - z 0
dc
2 Cc 2 Cc Cc
v dc 0 vdc
0
Hz R dc 1
0 0 - -
Ldc Ldc Ldc
0
0
1 1
0 0 0 -
Cdc Cdc R charge
avec
H d = 2 m.cos ( )
1
H = 1 m.sin ( )
(2.14)
q 2
- 52 / 188 -
Chapitre 2 : Prsentation du convertisseur
R ac vd
-L L
Hd
ac
-w - 0 0
ac
L ac
id w id 0
R ac Hq
i
- - 0 0
d 3 H d q
L ac L ac
iq
v cc = vcc + 0
3 Hq Hz
dt 2 Cc i
0 - 0 (2.15)
i Ldc Ldc
2 Cc Cc
v 0 v dc 0
dc
Hz R dc 1
0 - -
L dc Ldc L dc
0 0
1 1
0 0 -
Cdc Cdc R charge
La relation (2.15) caractrise le convertisseur Cuk rversible. Elle est la base thorique de
toute la suite de notre tude (dimensionnement du circuit de puissance et du systme
dasservissement). Le seul paramtre variable dans la matrice carre de (2.15) est la fonction
Hz (figure 2.2.3). Le dveloppement en srie de Fourier de cette fonction conduit :
m
3 3 1 1
cos 6nwt - 6n + = D z + d z_6n
3 3
2 n=1 6n -1 6n +1 6
H z = 1- m+ - (2.16)
2
Pour attnuer leffet des variations de cette fonction, le dimensionnement des lments
ractifs du circuit de puissance doit viter toute rsonance avec les frquences 6nf1 (multiples
de 300 Hz).
Linductance Lac tant identique pour les trois phases, la dmarche de dimensionnement est
labore pour une seule phase (phase a). Le rle de cette inductance est de filtrer le courant de
phase et non pas de stocker lnergie. Nous allons dimensionner cette inductance pour une
puissance de 6 kW, en considrant la frquence de dcoupage fixe et en autorisant une
ondulation maximale du courant de phase I = 1 A.
Le courant de phase crot lorsque la fonction de modulation associe est gale 1 et dcrot
lorsquelle est gale zro (figure 2.3.1). Daprs (2.8), le courant de phase ia sexprime sous
la forme d'une quation diffrentielle du 1er ordre :
= -R aci a + Va sin ( wt ) - v cc ( 2H a - H b - H c )
di a 1
Lac (2.17)
dt 3
- 53 / 188 -
Chapitre 2 : Prsentation du convertisseur
Daprs la figure 2.3.1, nous avons I = Imax - Imin = ia(t2) - ia(t1). En exprimant (2.17) dans
lintervalle [t1 ; t2], nous obtenons :
- v cc 2H a ( t1 ) - H b ( t1 ) - H c ( t1 )
2
1
3
Ha ( t ) =
m.sin wt - -
3
6
2
2H a ( t ) - H b ( t ) - H c ( t ) = m.sin ( wt - )
3
i a_M = I M .sin ( wt )
2
m = Vac_amplitude
avec
1
Vcc
2
Lac w.I M
= arctan V
a
ia
Imax
ia_M
Imin
Ha
t t1 t2 t + Td
Td/2
Figure 2.3.1 : Evolution du courant de phase ia pendant une demi priode de dcoupage
Ha ( t )
La frquence de dcoupage Fd peut ainsi scrire :
-R ac i M - sin ( wt ) + Va sin ( wt )
I
1
Td 2.Lac .I 2
Fd = =
- v cc 2H a ( t ) - H b ( t ) - H c ( t )
(2.19)
1
3
- 54 / 188 -
Chapitre 2 : Prsentation du convertisseur
d H a (t )
Si lon se met dans le cas o I est maximal, cela correspond = 0 , et la valeur de
dwt
2
wt satisfaisant cette galit est wt = + . En remplaant cette valeur de wt dans (2.19),
3
on obtient les courbes de la figure 2.3.2.
I = 0,25 A
I = 0,5 A
I = 0,75 A
I = 1 A
I = 1,5 A
Daprs la figure 2.3.2, nous constatons quune inductance de 2,5 mH autorisera une
enveloppe maximale du courant de largeur 1 A pour une frquence de dcoupage 10 kHz.
tant donn que la modulation RCFM utilise des frquences infrieures ou gales 10 kHz,
nous allons adopter cette valeur dinductance dans toute la suite de ce travail.
Le rle des composants du circuit de dcouplage (Cdc, Ldc) est le stockage intermdiaire de
lnergie, tout en attnuant leffet des ondulations rsultantes de la composante dz_6n. Alors
que le condensateur Cdc stabilise la tension de sortie vdc. Nous dimensionnons Cdc pour une
fluctuation maximale de 1% de la tension de sortie.
Dans le cas du redresseur Boost, le dimensionnement du condensateur de sortie est
relativement simple. Car lquation (2.15) permet de dterminer facilement le comportement
de la tension de sortie en fonction de la capacit du condensateur de sortie. Alors que dans le
cas du redresseur Cuk, la prsence de la fonction Hz et des deux variables supplmentaires iLdc
- 55 / 188 -
Chapitre 2 : Prsentation du convertisseur
et vcc augmente le degr de complexit du systme. Ceci engendre une grande diffrence en
terme danalyse entre les deux structures Cuk et Boost.
Dans la suite de cette section, nous posons x = (id, iq, vcc, iLdc, vdc)t vecteur dtat. Pour
rsoudre lquation (2.15), nous nous sommes bass sur le travail de Kikuchi [KIK-02] qui
consiste valuer chaque variable du vecteur dtat en rgime dynamique. Dans un premier
temps, nous dcomposons cette quation en deux parties : statique (2.20) et dynamique (2.21),
ensuite dans un second temps, nous valuons le comportement des variables dtat associs
chaque lment ractif.
R ac Vd
-L L
Hd
ac
-w - 0 0
ac
L ac
I
0 w d 0
R ac Hq
Iq
- - 0 0
0 3 H
L ac L ac
0 = Vcc + 0
3 Hq Dz
2 Cc
d
0 - 0 (2.20)
0 I Ldc
2 Cc Cc
0 0 0
Dz R dc 1
Vdc
0 - -
Ldc Ldc Ldc
0 0
1 1
0 0 -
Cdc Cdc R charge
R ac
- L 0
Hd
ac
-w - 0 0
Lac
I d_6n w I d_6n
0
R ac Hq
I q_6n I q_6n
- - 0 0
I Ldc
Lac Lac
3 Hd
j 6nw Vcc_6n = cc_6n + - .d z_6n
3 Hq Dz
c
V
I
0 - 0 (2.21)
Ldc_6n Ldc_6n V
2 Cc 2 Cc Cc C
cc
I
Vdc_6n 0
Ldc Vdc_6n Ldc
Dz R dc 1
0 - -
Ldc Ldc
0 0
1 1
Cdc R charge
0 0 -
Cdc
La premire partie (2.20) dcrit le comportement en rgime statique du vecteur dtat, elle
permet de dterminer les composantes continues des signaux ainsi que la charge rsistive
utiliser pour que le convertisseur puisse fonctionner dans la gamme de puissance souhaite
(6 kW). La deuxime partie (2.21) prsente les quations rgissant le rgime dynamique du
vecteur dtat. Son rle est dvaluer les valeurs de (Cc, Ldc, Cdc) permettant dempcher la
- 56 / 188 -
Chapitre 2 : Prsentation du convertisseur
rsonance par rapport aux frquences 6nf avec une attnuation maximale des oscillations dues
dz_6n. Lquation (2.21) peut tre rsolue pour chaque frquence 6nf, et le vecteur dtat
rsultant est le suivant :
R ac
- L 0
Hd
ac
-w - 0 0
L ac
w 0
R ac Hq
- - 0 0
3 H I
Lac L ac
A= , B = - dc
3 Hq
Cc
Dz
2 Cc
d L
0 - 0
2 Cc Cc
0 Vcc
L
Dz R dc 1
Ldc
0 - -
dc
Ldc Ldc
0
Cdc R charge
1 1
0
0 0 -
Cdc
Le dveloppement de (2.21) nous a permis de dterminer lvolution des signaux Vcc-6n, ILdc_6n
et Vdc_6n respectivement en fonction de Cc, Ldc et Cdc. La figure 2.3.3 illustre lvolution de la
tension Vcc_6n pour les cinq premires composantes de 6nf.
|Vcc_6n| [V]
Cc [F]
Figure 2.3.3 : Evolution de |Vcc_6n| en fonction de Cc pour les 5 premires composantes de 6nf
(Rac=0,33 ; Lac=2.5mH ; Rcharge=15 ; ILdc=20A ; Rdc=0.24 ; Ldc=2,5mH ; m=0,75 ; =3,5 ; Vcc=810V ; Vdc=300V ; Vd=311V)
- 57 / 188 -
Chapitre 2 : Prsentation du convertisseur
Nous constatons que la valeur de la capacit Cc doit tre suprieure 34 F qui correspond
une rsonance suprieure 350 V pour n = 1. Pour avoir une ondulation de 1% de la tension
vcc (8 V), la capacit utiliser est de lordre de 470 F.
|ILdc_6n| [V]
Ldc [H]
Figure 2.3.4 : Evolution de |ILdc_6n| en fonction de Ldc pour les 5 premires composantes de
6nf (Rac=0,33 ; Lac=2.5mH ; Rcharge=15 ; ILdc=20A ; Rdc=0.24 ; Cc=470F ; Vdc=300V ; Vd=311V ; m=0,75 ; =3,5 ; Vcc=810V)
- 58 / 188 -
Chapitre 2 : Prsentation du convertisseur
|Vdc_6n| [V]
Cdc [F]
Figure 2.3.5 : Evolution de |Vdc_6n| en fonction de Cdc pour les 5 premires composantes de
6nf (Rac=0,33 ; Lac=2.5mH ; Rcharge=15 ; ILdc=20A ; Rdc=0.24 ; Ldc=2,5mH ; m=0,75 ; =3,5 ; Vcc=810V ; Vdc=300V ; Cc=470F)
Dans une tude complmentaire, nous avons dtermin lvolution de la tension de sortie en
fonction de Cc et Ldc. Avec les valeurs choisies ci-dessus, nous navons relev aucun
problme de rsonance sur la tension de sortie.
Aprs le dimensionnement du circuit de puissance du convertisseur, la deuxime tape de
notre tude consiste dimensionner les boucles dasservissement du courant de phase et de la
tension de sortie.
- 59 / 188 -
Chapitre 2 : Prsentation du convertisseur
sans ce symbole correspondent aux grandeurs physiques mesures, et une variable dtat en
majuscule spcifie la partie continue du signal associ.
Dans la boucle de courant, la commande consiste asservir les projections dans le rfrentiel
de Park des courants de phase. Cette structure, qui met profit que le systme est quilibr, ne
comporte que deux boucles de courant couples (composantes id et iq) [LEE-00]. Daprs
(2.15), les courants de phase sont rgis par les deux premires quations. Etant donn que le
condensateur de dcouplage Cc est dimensionn pour attnuer les ondulations 6nf et comme la
tension dentre est purement sinusodale, nous pouvons exprimer ces deux quations sous la
forme suivante :
di d
Lac = -R aci d - Lac wi q - H d Vcc + Vd (2.23)
dt
di q
Lac = -R aci q + Lac wi d - H q Vcc (2.24)
dt
En se basant sur (2.14), HdVcc et HqVcc scrivent :
vaca
vacd
avec la transforme de Park du vecteur tension vacb , et Vac _ amplitude lamplitude des
acq v
acc
v
tensions vac( a ,b ,c ) .
En remplaant (2.25) et (2.26) dans (2.23) et (2.24), la transforme de Laplace des deux
Daprs (2.27) et (2.28) le modle du procd dcrivant le courant de phase ne dpend pas de
Hz, ce qui engendre une boucle de courant similaire celle du convertisseur Boost rversible
[LEE-00]. Nous avons labor le modle de la figure 2.4.1 en se basant sur ces deux relations.
La boucle de rgulation du courant de phase est tablie sous lhypothse que v*acd = vacd et
- 60 / 188 -
Chapitre 2 : Prsentation du convertisseur
v*acq = vacq . Hypothse justifie par le fait que les frquences de la modulation RCFM sont
suffisamment leves par rapport 50 Hz. Les consignes de courant ( id* , iq* ) sont issues de la
Vd
- *
vacd +
id* + id
RCFM 1
+ Modulation - + 1
Correcteur
- vacd Lac s + Rac
-
Lac w
Lac w
Lac w
Lac w
Vq = 0
iq* + + + + iq
RCFM 1
+ Modulation - + 1
Correcteur *
vacq vacq Lac s + Rac
-
Si lon considre que la boucle de courant assure un facteur de puissance unitaire et que la
tension dalimentation est purement sinusodale, le modle du procd en rgime permanent
est quivalent une fonction de transfert du premier ordre (2.29).
G (s ) =
1
(2.29)
Lacs + R ac
La synthse du correcteur de la boucle de courant est tablie pour cette fonction de transfert.
Toutefois, le dcoupage du courant engendre des petites variations dues au couplage entre les
deux composantes id et iq. Dans le systme de commande, nous avons ajout Lacwi(d,q) la
sortie des correcteurs afin de tenir compte de ces variations dans les signaux corrigs sans
toutefois les amplifier. Le correcteur utilis dans la boucle dasservissement est de type PI
[LEE-00]. Les coefficients Kpi et Kii de ce convertisseur sont choisis pour avoir un systme
dasservissement stable, ayant une marge de phase suprieure 45 et une frquence de
coupure avoisinant 2 kHz.
- 61 / 188 -
Chapitre 2 : Prsentation du convertisseur
La figure 2.4.2 illustre le diagramme de bode du modle de courant en boucle ouverte avec les
coefficients du correcteur suivant : Kii = 310 et Kpi = 32. Dans ce diagramme, nous ne
relevons aucun problme de stabilit, la marge de phase est de 90 et la frquence de coupure
est de lordre de 2 kHz. Pour tester la validit de ce dimensionnement, nous avons simul le
cas correspondant une transition du mode redresseur MLI vers le mode onduleur non
autonome, tout en introduisant les variations dues une modulation RCFM sur lintervalle de
frquence [6 kHz ; 10 kHz]. La figure 2.4.3 prsente les rsultats de simulation de ce cas de
fonctionnement.
- 62 / 188 -
Chapitre 2 : Prsentation du convertisseur
(A)
(s)
(s)
(b) Sortie de la boucle de courant id
- 63 / 188 -
Chapitre 2 : Prsentation du convertisseur
Daprs (2.15), le ct continu du convertisseur est rgi par les trois dernires quations et
sexprime sous la forme du systme (2.30). Contrairement au convertisseur Boost triphas, il
ny a pas de relation liant directement la tension de sortie au courant de phase. De plus, la
fonction Hz tant variable dans le temps, elle engendre un procd non linaire et variable
dans le temps. Ceci rend le dimensionnement du systme dasservissement assez complexe
compar au cas du redresseur Boost triphas.
1
0 Hd Hq
1 1
v cc v Cc
- Hz 0
d 1 3 id
Cc Cc
iL = i Ldc + 0 0
dt dc L dc
cc
R dc 1
v 2 iq
Hz - - (2.30)
v dc dc
Ldc L dc
0 0 0
1 1
-
Cdc Cdc R Load
Le modle rsultant de ces trois quations est prsent dans la figure 2.4.4. Dans ce modle la
boucle de tension est tablie sous lhypothse dune boucle de courant parfaite.
iL dc
BC* 1 d 3 H d
vdc* id* i - 1 vdc
+ 1 1
Correcteur Hz
+- 2 + sCc vcc +- Ldc s + Rdc +- sCdc
+
iq* = 0 BC* 1 iq 3 H 1
q
2 R charge
FPB**
(*) : Boucle de courant ; (**) : Filtre Passe Bas (bande passante 30 Hz)
La consigne de courant de phase tant issue de la boucle de tension, cette dernire doit avoir
une bande passante rduite afin de ne pas introduire les ondulations de la tension de sortie
dans id. Cette bande passante doit tre infrieure 50 Hz [LEE-00]. Pour satisfaire cette
contrainte, nous avons insr un filtre passe bas dans la chane de retour du modle de la
figure 2.4.4. La simulation du convertisseur nous a amen choisir un filtre avec une bande
passante de 30 Hz.
- 64 / 188 -
Chapitre 2 : Prsentation du convertisseur
Partant de ce postulat, toutes les oscillations dues dz_6n pourraient tre ngliges. Ceci
conduit modliser un systme non linaire invariant dans le temps, et autorise une sparation
entre les rgimes dynamique et statique. La relation (2.30) devient alors :
- ( D z + d z )
1
Vcc + v cc V + v
0 0
1 cc
C
( Dz + d z )
d
c
1 1
0 -
Cdc Cdc R Load
( D q + H q )
1
( Dd + H d )
(2.31)
1
Cc
3 Id + i d
Cc
+
2 Iq + i d
0 0
0 0
La partie continue de la relation (2.31) sexprime sous la forme (2.32), et spcifie un lien
direct entre la composante continue Id du courant de phase dans le rfrentiel (d ; q) avec les
composantes continues de tension de sortie et du courant de charge (2.33).
1
0 Dd Dq
1 1
0 V Cc
- Dz 0
1 3 Id
Cc Cc
0 = L Dz I Ldc + 0 0
cc
R dc 1
0 dc V 2 Iq
- - (2.32)
dc
L dc L dc
0 0 0
1 1
-
Cdc Cdc R Load
Vdc = 2 D R charge Id
3 Dd
I =
z
(2.33)
d 3 Dd charge
2 D z
I
Cette reprsentation permet davoir un lien direct entre Id et Vdc. Ce qui se traduit par une
boucle de tension qui fonctionnera correctement avec un simple correcteur proportionnel.
Nous avons simul le modle de la figure 2.4.4 sur Matlab/Simulink en introduisant le gain
3Dd R charge
K= lentre de la boucle. La figure 2.4.5 illustre la rponse de ce modle en
2D z
boucle ouverte une consigne chelon de 300 V. Malgr que la sortie se stabilise en rgime
- 65 / 188 -
Chapitre 2 : Prsentation du convertisseur
permanent, les problmes suivant se posent : le temps de rponse est trs lev, le rgime
transitoire est incontrlable, lerreur statique est non nulle et les variations de la tension de
sortie restent inacceptables vue que lun des objectifs de ce convertisseur est davoir une
consigne de tension de sortie variable (abaisseur / lvateur).
(V)
(s)
Figure 2.4.5 : Rponse en boucle ouverte du modle de la boucle de tension
une consigne chelon de 300 V
-1
0 C I Ldc
1
v cc v c
- Dz 0
1 1
Cc
d 1
= i Ldc + Vcc H z
dt dc Ldc
cc
R dc
i
L dc
v dc dc
Dz - -
vdc
L
L dc L
0
0
1 1
Cdc R Load
-
Cdc
1 1 1
(2.34)
Dd Dq Id Iq
1
3 i d 3 H d
C c
C c
C c
C c
+ 0 0 + 0 0
2 i q 2 H q
0 0 0 0
(
sortie ( vdc ) et le vecteur dtat v cc , i Ldc , v dc . )
t
- 66 / 188 -
Chapitre 2 : Prsentation du convertisseur
dans ltude de cette boucle, la fonction Hz peut scrire sous la forme suivante :
H z 1-
3 3
H d2 + H q2 (2.35)
Cette reprsentation de Hz permet dexprimer Hz en fonction de Hd et Hq.
H z H z
H d H q
H z = H d + H q
Dd ,Dq Dd ,Dq
3 3
(2.36)
H q
Dd Dq
Dd2 + Dq2
=- H d +
2 2
Dd + Dq
v cc v Cc
- Dz 0
1 1 3 i d
C Cc
d
c
i Ldc = i Ldc + 0 0
cc
R
dt L dc L dc
v dc
2 i q
Dz - dc -
v dc
1
L dc
0 0 0
1
Cdc R Load
-
Cdc
3 Iq
(2.37)
kDd
I Ldc 3 I d I
Cc 2 Cc
+ kD q Ldc +
H d
Cc 2 Cc
+ -kD d cc
V V
H q
-kDq cc
Ldc L dc
0 0
avec
3 3 1
k= (2.38)
Dd2 + Dq2
Dans lhypothse dune boucle de courant parfaite, les composantes d-q du courant de phase
vues de la boucle de tension restent invariables. Ceci autorise lapproximation suivante :
0 -R aci d - Lac wi q - H d v cc + v d
0 -R aci q + Lac wi d - H q v cc
(2.39)
De plus, vu que la boucle de courant assure un facteur de puissance unitaire (iq = 0) et que la
tension dalimentation est sinusodale, les fonctions Hd et Hq deviennent :
- 67 / 188 -
Chapitre 2 : Prsentation du convertisseur
H d H d
i d vcc
H d = i d + v cc
i d + 2 ( R ac Id - Vd ) v cc
Id ,Vcc Id ,Vcc
(2.40)
R ac 1
=-
Vcc Vcc
H q H q
i d vcc
H q = i d + vcc
Id ,Vcc Id ,Vcc
(2.41)
1 1
= wLac i d - 2 wLac Id v cc
Vcc Vcc
Le remplacement de (3.40) et (3.41) dans (3.37), mne lquation dtat suivante :
x& = Ax + Bu
y = Cx
(2.42)
avec
v cc
i d
C = ( 0 0 1) ,
x = i Ldc , u= , y = v dc
v i q
(2.43)
dc
R ac I d - Vd I Ldc 3 Id wL ac I d
I Ldc 1
Vcc Cc 2 Cc Vcc2
kD d
+ - kD q
- Dz 0
2
Cc Cc
A=
1 R I -V wLac I d R dc 1
D z - kD d ac d d + kD q - - (2.44)
L dc Vcc L dc Vcc L dc L dc Ldc
1 1
0 -
Cdc Cdc R Load
3 1 I Ldc 3 I d wL ac
kDd + Dq
I 3 1
2 Cc Cc 2 Cc Vcc
+ kD q Ldc
Cc 2 Cc
B= 0
R wL ac
kD d ac - kD q (2.45)
Ldc L dc
0
0
- 68 / 188 -
Chapitre 2 : Prsentation du convertisseur
Lobjectif initial de notre travail tant de concevoir une structure de convertisseur Cuk
susceptible de concurrencer la structure Boost triphase, nous avons opt pour un correcteur
simple de type PI comme dans lasservissement de la tension de sortie dans le Boost. En
utilisant un correcteur PI avec les coefficients Kiv = 2.5 et Kpv = 0.25 nous avons obtenu le
diagramme de bode du systme en boucle ouverte de la figure 2.4.6.
Comme nous pouvons le constater dans la figure 2.4.6, la bande passante de la boucle de
tension est de lordre de 20 Hz avec une marge de phase trs satisfaisante. Nous avons simul
dans un premier temps la rponse indicielle en boucle ferme pour une variation de la
consigne de 0 V 300 V (cas reprsentant la mise sous tension du convertisseur), ensuite dans
un second temps nous avons baiss la consigne de 300 V 200 V. Le rsultat de cette
simulation est illustr dans la figure 2.4.7. Nous nobservons aucun dpassement de la tension
- 69 / 188 -
Chapitre 2 : Prsentation du convertisseur
de sortie par rapport la consigne dans les deux cas avec un temps de rponse satisfaisant et
rpondant aux exigences fonctionnelles de la boucle de tension.
Dans la boucle de tension, nous avons introduit un filtre passe bas dans la chane de retour
pour viter une amplification des fluctuations de la tension par le correcteur PI. Cependant, ce
choix dtriore le comportement dynamique de cette boucle lors de la variation du courant de
charge. Dans la figure 2.4.9, nous constatons quune variation du courant de charge de 20 A
10 A engendre un temps de rponse de 45 ms.
- 70 / 188 -
Chapitre 2 : Prsentation du convertisseur
- 71 / 188 -
Chapitre 2 : Prsentation du convertisseur
- ia iq
w.Lac
b - ib Transformation
Observateur abc-dq g1 g7
dtat c ic id
w.Lac
Adaptation des
icoup Estimateur signaux de
de courants
commande
-
V*dc
PI
Im1 +
PI -- V*d
V*aca
iLoad 1 Im2
1 + s / wFF
La figure 2.4.10 illustre ce systme de commande. Les courants de phase sont reconstruits
partir du courant mesur dans le circuit de dcouplage. Alors que la tension de sortie est
reconstruite partir des courants circulant dans la charge et dans le circuit de dcouplage. Les
dmarches et les algorithmes de reconstruction des courants de phase et de la tension de sortie
seront dvelopps dans le chapitre 4 de ce mmoire.
- 72 / 188 -
Chapitre 2 : Prsentation du convertisseur
Conclusion
- 73 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
- 74 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
Introduction
Comme nous lavons signal au chapitre1, la frquence alatoire varie entre deux valeurs fb et
fa reprsentant respectivement les frquences minimale et maximale. Le choix de cet
intervalle ainsi que la loi de variation de la frquence alatoire sont les deux paramtres les
plus importants dans la conception dune telle modulation. En plus, lanalyse spectrale des
signaux obtenus ne peut tre effectue correctement sans la connaissance pralable de ces
deux paramtres. Dans notre cas, nous avons choisi dutiliser cette technique de modulation
pour deux raisons :
Ltude pralable de la technique utilise, la dfinition des diffrents paramtres ainsi que le
choix optimis de lintervalle [fb ; fa] reprsentent trois axes indispensables pour
limplmentation de cette technique ainsi que lanalyse des rsultats.
Dans ce chapitre, aprs la prsentation des lois de variation de la frquence alatoire, nous
effectuerons une comparaison entre ces diffrentes approches. Ensuite, dans la deuxime
partie, les hypothses de ltude ainsi que lanalyse spectrale de la modulation RCFM dans le
convertisseur Cuk, seront prsentes. Enfin le dernier paragraphe sera consacr
linterprtation et l'analyse des rsultats et au choix de lintervalle [fb ; fa].
- 75 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
mieux le fonctionnement de cette technique. Pour cette raison, nous avons choisi dutiliser T
au lieu de f dans la suite de notre tude thorique.
Dans la modulation RCFM, la fonction rgissant les variations de la priode de dcoupage est
le paramtre principal qui dtermine le contenu spectral du signal modul. Cette fonction est
appele la fonction "densit de probabilit (PDF)" dfinie par :
Distribution
P(T) Gaussienne P(T)
Distribution
uniforme
p1
pN
pi
2T p2
Ta T Tb T T1 T2 Ti TN
P(T) = pi .(T Ti )
N
(3.2)
i =1
Avec N le nombre total des priodes entre Ta et Tb, et pi est le poids associ la ime priode
Ti. Ces poids sont strictement positifs et vrifient la relation (3.3).
p =1
N
(3.3)
i =1
i
- 76 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
Toutes les lois de probabilit ont pour origine la loi de distribution uniforme (figure 3.1.1a).
Pour cette raison nous avons choisi dutiliser deux lois de probabilit comme support pour
notre tude : distribution uniforme et non uniforme (Gaussienne).
Dans le cas dune PDF continue et uniforme, la priode alatoire T est dtermine en temps
rel partir de la relation (3.4) o x est un nombre alatoire dans lintervalle [0,1] suivant
cette loi de probabilit.
Si cette PDF est non uniforme, on utilise des mthodes de transformation permettant de passer
de la distribution uniforme la loi dsire. Par exemple dans [BEC-00], un gnrateur de
priodes de dcoupage alatoires utilisant la transformation de la distribution uniforme en une
distribution gaussienne a t labor. Lutilisation de telles distributions exige des calculs
supplmentaires compars au cas de la distribution uniforme. De plus, les travaux de [ELK-
02] ont dmontr que les spectres rsultant des distributions uniformes et non uniformes sont
similaires dans le cas de la RCFM.
Dans le cas de la fonction de densit de probabilit discrte, le nombre de priodes alatoires
Ti est fini, et ces priodes sont dtermines partir de (3.5). Ce type de gnrateur correspond
un choix alatoire de T dans un ensemble {Ti} en effectuant une srie de comparaisons
[BEC-00]. Cette procdure ncessite N-1 oprations de comparaison une fois que x est
connue. Cest une solution trs avantageuse pour une petite valeur de N, mais reste trs lourde
mettre en place pour N >> 10 dans le cas d'un systme temps rel.
T1 si 0 x < p1 ,
T si p x < p + p ,
T= 2
M
1 1 2
(3.5)
TN si p1 + p 2 + L + p N 1 x < 1.
- 77 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
En rsum, on peut considrer quune fonction de densit de probabilit continue ayant une
loi de distribution uniforme est une solution optimale pour le systme de commande de la
structure existante du convertisseur Cuk rversible, dune part par sa simplicit, et dautre part
par sa rapidit d'implmentation.
Les travaux mens sur les convertisseurs DC/DC [TSE1-00] [TSE-02] [TSE3-00] et DC/AC
[LIA-00] [TRZ-94] en utilisant un dcoupage avec une densit de probabilit continue ont
tous dmontr que le spectre de tension est une fonction purement continue, except une
composante discrte correspondant la frquence moyenne et multiples. Cela peut se
gnraliser pour tous les signaux rsultant dune telle modulation [BEC-00].
Cependant dans le cas des distributions discrtes, lattnuation du spectre n'est pas garantie.
Car de nouveaux harmoniques peuvent apparatre des frquences non multiples de la
frquence fondamentale [BEC-99] si lensemble des priodes Ti nest pas choisi
correctement. Aussi, des harmoniques peuvent apparatre la frquence correspondant au plus
petit multiple commun des frquences utilises et ses multiples. Mais les distributions
discrtes peuvent aussi donner exactement les mmes rsultats que les distributions continues,
condition dutiliser un ensemble de priodes dont les frquences associes sont des
multiples de la fondamentale, et que le plus petit multiple commun des frquences utilises se
trouve loign dans le spectre [BEC-00].
Lutilisation dune fonction "densit de probabilit" ayant une distribution continue assure
une attnuation de la densit spectrale en haute frquence sans aucune contrainte sur la loi de
distribution, ce qui renforce nos conclusions dans le dernier paragraphe. Toutefois, cette
garantie ne couvre pas le spectre basse frquence, car ce dernier dpend en plus de la loi de
probabilit de lintervalle de variation des priodes alatoires.
Afin de justifier lintrt de cette technique de modulation et de dterminer lintervalle de
variation de la frquence de dcoupage, les paragraphes (II et III) sont consacrs ltude
thorique de la modulation RCFM dans le convertisseur Cuk rversible.
- 78 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
solutions risquent de ne pas exister. Pour cette raison, nous utiliserons dans notre tude une
PDF discrte mais paramtre de telle sorte que la valeur moyenne de la frquence de
dcoupage concide avec celle du cas continu, et le spectre du courant de phase ne prsente
pas de discontinuits dans la plage de frquence analyse.
Plusieurs tudes menes sur cette technique de modulation ont t publies dans la littrature
[LIA-00] [TSE2-00]. Certes ces travaux ne concernent que les convertisseurs DC/DC et les
onduleurs, mais ils ont tous dmontr que lanalyse spectrale des signaux rsultants de cette
modulation se ramne lanalyse spectrale des fonctions de modulation associes aux cellules
de commutations. Dans le cas de notre convertisseur, la prsence du septime interrupteur
engendre une certaine dpendance entre toutes les fonctions de modulation. Dans la suite du
chapitre, nous analysons leffet de la modulation RCFM sur le courant de phase ia, les
rsultats et conclusions lis ce dernier s'appliquent aux deux autres courants de phase.
Daprs le chapitre 2, et en ngligeant la rsistance interne de linductance Lac nous avons :
di a
Lac dt = -H a .Vcc + v a - v bn
di b
Lac
= -H b .Vcc + v b - v bn (3.6)
dt
Lac dt = -H c .Vcc + v c - v bn
di c
En plus,
ia + ib + ic = 0
va + v b + vc = 0
(3.7)
La somme des trois quations de (3.6) en tenant compte de (3.7) permet de dduire que
- 79 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
Comme nous lavons voqu au chapitre II, les commutations des sept cellules sont labores
partir de quatre fonctions de modulations Ha, Hb, Hc et Hz. Ces fonctions de modulation
correspondent respectivement aux bras des phases a, b, c et au quatrime bras du circuit de
couplage. La figure 3.2.1 rcapitule ces fonctions et les compare aux fonctions de modulation
du redresseur Boost. Dans la figure 3.2.1, nous remarquons que la diffrence avec les
fonctions de modulation du redresseur Boost se situe au niveau des vecteurs zros o les
fonctions Ha,b,c sannulent. En plus, la phase correspondant la tension vac(a,b,c) la plus basse a
une fonction de modulation gale zro. Dans lexemple de la figure 3.2.1, cest la fonction
Hc qui se retrouve gale zro.
Vtriangle
Vaca.sin(wt - )
Vacb.sin(wt - - 2/3)
Vacc.sin(wt - + 2/3)
Ha_bst
Fonctions de modulation
Hb_bst Redresseur Boost rversible
Hc_bst
Ha
Hb Fonctions de modulation
Redresseur Cuk rversible
Hc
Hz
T0 T1 T2 T0 T0 T2 T1 T0
Figure 3.2.1 : Comparaison des fonctions de modulation des deux structures Boost et Cuk
- 80 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
c a b
Ha
t
Figure 3.2.2 : La fonction de modulation de la phase a
Pour que notre tude soit ralisable, nous avons utilis certaines hypothses facilitant la
modlisation des fonctions de modulation du convertisseur. Ces hypothses sont :
(1) Les impulsions de la fonction de modulation sont centres par rapport chaque priode
de dcoupage (figure 3.2.3).
(2) Le signal de rfrence est purement sinusodal (figure 3.2.2).
(3) On suppose que les priodes utilises pour le dcoupage sont connues pralablement.
Nous utilisons J priodes gnres alatoirement suivant une loi de probabilit
uniforme. Lintrt de ce choix est de rduire le temps de calcul et aussi de permettre
lchantillonnage du signal correspondant au rapport cyclique moyen.
Tn,m
- 81 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
H a (t) = H
N K
(t) (3.10)
n = N m =1
a n ,m
Ha
K
(t) reprsente la fonction de modulation correspondant au nme bloc (une priode T1).
m =1
n,m
Han,m reprsente la fonction de modulation au sein dune priode de dcoupage. Elle scrit :
tn,m reprsente linstant initial de la priode de dcoupage Tn,m. n,m dtermine la position de
chaque impulsion par rapport tn,m, et n,m reprsente la largeur des deux impulsions au sein
de chaque priode Tn,m (figure 3.2.3). Ainsi la fonction Ha(t) devient :
(d a n ,m d bn ,m ) si d bn ,m < d cn ,m
a n,m = a n ,m Tn,m avec a n ,m =
(d a n ,m d cn ,m ) si d bn ,m > d cn ,m
(3.14a)
- 82 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
a n,m = a n ,m Tn,m avec a n ,m = (1 da n,m )
1
'a n,m = a n ,m Tn,m = (1 a n ,m 2 a n ,m )
2 (3.14b)
avec a n ,m
Ha Ha Ha
Hb Hb Hb
Hc Hc Hc
H an H an H an
Figure 3.2.5 : Reprsentation de la fonction Han pour toutes les combinaisons possibles dans
lordre des fonctions de modulation (Ha, Hb et Hc)
Afin de modliser la fonction Han, nous lavons paramtre de telle sorte quelle regroupe les
trois cas possibles de la figure 3.2.5. Le modle thorique de cette fonction est prsent dans
la figure 3.2.6.
t n,m 1T 1T t n,m + T
- 83 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
( )
Les diffrents paramtres du modle de la figure 3.2.6 sont les suivants :
1
n,m = 1- max n,m , (3.15)
2
1n,m =
1
2
( )
max n,m - midn,m , (3.16)
2n,m =
1
2
( )
midn,m - min n,m , (3.19)
La figure 3.2.7 illustre la fonction Han sur une dure temporelle de 60ms. La modlisation de
cette fonction sur 2N+1 priodes fondamentales peut scrire sous la faon suivante :
( ) (
H ( t ) = H1 ( t ) + H 2 t - 1n,m Tn,m + H 2 ( t - 'n,m Tn,m ) + H1 t - 'n,m Tn,m - 2n,m Tn,m )
o K est le nombre total des priodes de dcoupage dans une priode fondamentale.
(3.22)
avec n,m
'
= 1n,m + 2n,m + n,m
(3.10) et (3.21) ne sont valables que lorsque le nombre K est le mme dans tous les blocs, ce
qui implique un chantillonnage pralable de tous les blocs en K chantillons.
- 84 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
Han
Nous utilisons dans ltude thorique une loi de distribution uniforme discrte, et paramtre
de telle sorte que la diffrence avec la distribution continue soit minime. En plus des
arguments cits au dbut de cette section (II), ce choix est justifi pour deux raisons :
premirement parce que cest la seule approche qui permet davoir une connaissance
pralable des priodes de modulation, et deuximement car lutilisation de la distribution
continue augmente le nombre des paramtres et des variables, et conduit des calculs
dintgrales trs lourds.
Rappelons que dans le cas des signaux priodiques, la densit spectrale de puissance S(f) est
la transforme de Fourrier de la fonction dautocorrlation (). Lexpression de la fonction
dautocorrelation pour un signal priodique x(t) est la suivante :
( ) = lim x(t).x * (t ) d
1 T0
(3.25)
T0 2T T0
0
S(f ) = ( ) e j2f d
+
(3.26)
Dans le cas des signaux alatoires et daprs [TSE2-00] et [BEC-00], S(f) scrit :
S(f ) = R( ) e j2f d
+
(3.27)
- 85 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
T0 2T T0
0
E {g(T)} = pi .g(Ti )
J
donc
i =1
T0 2T T0
N (2N + 1)T
H(t (t
0
H(t (t
1
+ n,m % % + )) dt
N K
% % Tn,m
n% = N m
% =1
% %
n,m
H(t (t n,m
% % + n,m % % + )) =
+
(3.30)
jw 2 (t n% ,m% + n% ,m% Tn% ,m% + ) jw 2 t
% % Tn,m U n,m
% % (f 2 )e e df 2
U n,m ( f ) = U1n,m ( f ) 1 + e
-j2f ( 'n,m + 2n,m ) Tn,m
+ U 2n,m ( f ) e
-j2f.1n,m Tn,m
(
+e
-j2f. 'n,m Tn,m
) (3.31)
(1- e )
avec
(f ) =
A1n,m -j2f.1n,m Tn,m
U1 (3.32)
n,m
j2f
(f ) = (1- e )
A 2n,m -j2f. 2n,m Tn,m
U2 (3.33)
n,m
j2f
Soit n,m = n,m Tn,m , en utilisant (3.29) et (3.30), la fonction dautocorrlation devient :
- 86 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
( ) = lim
N (2N + 1)KT
K K
1 + N N + jw1 (t n ,m + n ,m ) jw1 t
% =1
U n,m (f1 ) e e
m =1 m n = N n% = N (3.34)
df 2 df1 dt
jw 2 (t n% ,m% + n% ,m% + ) jw 2 t
U n,m
% % (f 2 ) e e
( ) = m,m% ( )
K K
(3.35)
m =1 m
% =1
m,m% ( ) = lim
N (2N + 1)T
1 + N N + jw1 (t n ,m + n ,m ) jw1 t
U n,m (f1 )e e
1 n = N n% = N (3.36)
jw 2 (t n% ,m% + n% ,m% + ) jw 2 t
U n,m
% % (f 2 )e e df 2 df1 dt
e
(3.37)
jw 2 (t n% ,m% + n% ,m% + ) j(w1 + w 2 )t
U n,m
% % (f 2 ) e dt df 2 df1
Comme le terme entre crochs dans (3.37) ne dpend pas de t, nous pouvons simplifier (3.37)
en utilisant la proprit de la distribution de Dirac dcrite dans (3.38) :
e j(w1 + w 2 )t dt = (f1 + f 2 )
+
(3.38)
m,m% ( ) = lim
1 N N + jw1 (t n ,m + n ,m )
N (2N + 1)T
U n,m (f1 )e
1 n = N n% = N (3.39)
(f1 + f 2 ) df 2 df1
jw 2 ( t n% ,m% + n% ,m% + )
U n,m
% % (f 2 )e
- 87 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
m,m% ( ) = lim
1 N N +
N (2N + 1)T
U n,m (f1 )U*n,m
% % (f1 )
1 n = N n% = N (3.41)
jw1 (t n ,m t n% ,m% ) jw1 ( n ,m n% ,m% ) jw1
e e e df1
N (2N + 1)T
Nn
m,m% ( ) = lim
1 + N
U n,m (f )U*n + ,m% (f )
1 n = N = N n (3.42)
jw (t n ,m t n + ,m% ) jw ( n ,m n + ,m% ) jw
e e e df
La relation (3.42) peut tre transforme en (3.43) en introduisant les deux fonctions
y n,m et y n + ,m% :
1, si n < N 1, si n + < N
y n,m = y n + ,m% =
0, sin on 0, sin on
N (2N + 1)T
y
+ +
m,m% ( ) = lim
N
1
y n + ,m% U n,m (f )U*n + ,m% (f )
= n = N
n,m
1 (3.43)
jw ( t n ,m t n + ,m% ) jw ( n ,m n + ,m% ) jw
e e e df
Comme nous lavons dfini dans (3.15), la fonction R associe un signal alatoire est
lesprance mathmatique de sa fonction dautocorrlation. En appliquant cette proprit
(3.43), cette fonction scrit :
N (2N + 1)T
y
+
R m,m% ( ) = lim
1 + N
=
E y n + ,m% U n,m (f )U*n + ,m% (f )
}e
n = N
n,m
1 (3.44)
jw (t n ,m t n + ,m% ) jw ( n ,m n + ,m% ) jw
e e df
En permutant E{.} avec les sommes dans la relation (3.44), cette dernire devient :
N (2N + 1)T
=
*
E
}e
1 n = N (3.45)
jw (t n ,m t n + ,m% ) jw ( n ,m n + ,m% ) jw
e e df
y
+
soit D n, ,m,m% =
jw (t n ,m t n + ,m% ) jw ( n ,m n + ,m% )
y n + ,m% U n,m (f )U*n + ,m% (f )e e (3.46)
=
n,m
En considrant que la loi de distribution uniforme des priodes Ti est un processus alatoire
stationnaire dune priode fondamentale lautre, on peut crire :
- 88 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
alors :
La densit spectrale de puissance de Han est construite partir des densits spectrales de
puissance partielles comme le montre la relation (3.49)
S(f ) = Sm,m% (f )
K K
(3.50)
m =1 m
% =1
M O (3.51)
M M O M
SK,1 (f ) SK,2 (f ) L L SK,K (f )
Cette densit spectrale de puissance est la somme de tous les lments de la matrice carre de
degr K (3.51). Le calcul de chaque lment de cette matrice passe par lvaluation de la
relation (3.49). Cependant cette relation est assez difficile valuer cause de la prsence de
la somme infinie. Dans la section II.3.3 nous prsentons la mthode que nous avons mise en
place pour valuer cette relation.
Lune des solutions permettant de faciliter lvaluation des densits spectrales de puissance
partielles est dchantillonner suivant une priode T gale lesprance mathmatique des
priodes de dcoupage possibles Ti (figure 3.2.8).
- 89 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
dc da db
T = E {T} = pi Ti
J
(3.52)
i =1
- 90 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
da db dc
Nous avons labor un gnrateur de priodes alatoires de telle sorte que chaque priode
fondamentale contienne K priodes de dcoupage, et ce nombre K est le mme dans les
dix priodes. Ces priodes sont gnres suivant une loi de densit de probabilit uniforme (la
fonction rand() de Matlab). Cependant le nombre entier K de priodes de dcoupage doit
satisfaire les contraintes suivantes :
- la somme des K priodes de dcoupage au sein de chaque priode du fondamental
doit tre gale 20 ms.
- La valeur moyenne de la frquence de dcoupage doit correspondre celle rsultante
de lutilisation dune distribution continue.
- Le choix des priodes alatoires dans chaque intervalle est effectu sur un ensemble
de cinq priodes dfinies pralablement afin que les premires composantes
harmoniques haute frquence apparaissent au voisinage de 2,5 MHz.
Il est pratiquement impossible de respecter rigoureusement ces trois contraintes dans tous les
intervalles de frquences alatoires. Afin de rendre notre mthode plus flexible et applicable
des intervalles moins larges que [7kHz ; 10kHz], nous avons allg la premire contrainte. Au
lieu que la somme des K priodes alatoires soit gale exactement 20 ms, nous lui avons
associ une erreur relative maximale de 0.05%. Approximation justifiable par le fait que le
courant de phase varie dans une enveloppe sinusodale reprsentant une erreur de 5% (voir
dimensionnement de linductance Lac). Dans la suite, nous nommerons la mthode qui
respecte rigoureusement les contraintes 1 et 3 "mthode rigoureuse" et la mthode associant
une erreur relative de 0.05% la somme des K priodes alatoires "mthode adapte". La
figure 3.2.10 prsente une comparaison entre les deux mthodes dans le cas de lintervalle des
frquences alatoires [6 kHz ; 10 kHz] pour une distribution uniforme. Pour la mthode
- 91 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
adapte, nous utilisons lensemble des frquences {6, 7, 8, 9, 10} kHz avec K = 154, et pour
la mthode rigoureuse lensemble des frquences est {6, 8, 8.4, 8.75, 10} kHz avec K = 160.
Cette comparaison montre que les deux spectres varient de la mme faon et conservent la
mme composante fondamentale, avec un avantage pour la mthode adapte qui conserve la
mme frquence moyenne de dcoupage par rapport une distribution continue et uniforme.
0
Mthode rigoureuse
Fundamental period(T==20
20ms
ms)
-10 Mthode adapte
Fundamental period 10-3 ms) ms
(T ==2020(+/-0.001)
-20
-30
-40
S(f) [dB]
-50
-60
-70
-80
-90
-100
200 400 600 800 1000 1200 1400 1600 1800 2000
f [Hz]
Figure 3.2.10 : Comparaison des spectres basse frquence de Han rsultant de la mthode
adapte et de la mthode rigoureuse dans le cas de lintervalle [6kHz ; 10kHz].
- 92 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
Dbut
NON
Compteur == K
OUI
NON
Instant final - instant initial == 20ms 0.01ms
OUI
NON
i == 10
OUI
Fin
- 93 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
[5kHz ; 10kHz] {5, 7.5, 8.5, 9, 10} kHz 142 19.9935 [ms]
[8kHz ; 10kHz] {8, 8.5, 9, 9.5, 10} kHz 178 19.9912 [ms]
[9kHz ; 10kHz] {9, 9.25, 9.5, 9.75, 10} kHz 189 19.9954 [ms]
Tableau 3.1 : Caractristiques des gammes de frquence utilises pour gnrer la RCFM
Les figures de 3.3.1 3.3.6 prsentent les densits spectrales de puissance de la fonction Han
sur [50Hz ; 40kHz], respectivement pour les intervalles des frquences alatoires
[4kHz ; 10kHz], [5kHz ; 10kHz], [6kHz ; 10kHz], [7kHz ; 10kHz], [8kHz ; 10kHz] et
[9kHz ; 10kHz]. Ces densits spectrales sont continues, et nous remarquons que plus
lintervalle des frquences alatoires est large plus le spectre est attnu.
- 94 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
Figure 3.3.1 : Spectre de la RCFM utilisant des frquences alatoires dans [4kHz ; 10kHz]
Figure 3.3.2 : Spectre de la RCFM utilisant des frquences alatoires dans [5kHz ; 10kHz]
- 95 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
Figure 3.3.3 : Spectre de la RCFM utilisant des frquences alatoires dans [6kHz ; 10kHz]
Figure 3.3.4 : Spectre de la RCFM utilisant des frquences alatoires dans [7kHz ; 10kHz]
- 96 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
Figure 3.3.5 : Spectre de la RCFM utilisant des frquences alatoires dans [8kHz ; 10kHz]
Figure 3.3.6 : Spectre de la RCFM utilisant des frquences alatoires dans [9kHz ; 10kHz]
- 97 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
Dans la plage de frquence [2kHz ; 40kHz], nous constatons que plus lintervalle des
frquences alatoires est large plus le spectre est attnu. Par exemple, lutilisation de
lintervalle [4kHz ; 10kHz] prsente une attnuation de plus de 10 dB par rapport
lintervalle [9kHz ; 10kHz]. Dans la plage [95kHz ; 125kHz], les remarques concernant
lattnuation et la continuit du spectre restent valable (figure 3.3.7). Cependant, la diffrence
en terme dattnuation du spectre entre les intervalles des frquences alatoires devient
minime. Donc en trs haute frquence (frquence suprieure 100 kHz) tous les intervalles
des frquences de modulation se comportent de la mme faon. En rsum, en haute
frquence les intervalles de modulation les plus intressants sont plus larges que
[7kHz ; 10kHz].
-50 -50
-60 -60
-70 -70
-80 -80
-90 -90
-100 -100
-110 -110
-120 -120
-130 -130
95 100 105 110 115 120 125 95 100 105 110 115 120 125
kHz kHz
-50 -50
-60 -60
-70 -70
-80 -80
-90 -90
-100 -100
-110 -110
-120 -120
-130 -130
95 100 105 110 115 120 125 95 100 105 110 115 120 125
kHz kHz
-50 -50
-60 -60
-70 -70
-80 -80
-90 -90
-100 -100
-110 -110
-120 -120
-130 -130
95 100 105 110 115 120 125 95 100 105 110 115 120 125
kHz kHz
- 98 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
A ce stade de ltude, le comportement des harmoniques basse frquence est le critre qui
reste vrifier pour dterminer lintervalle optimal de la modulation RCFM. En basse
frquence, ce critre est le respect de la norme EN 61000-3-2.
Le paragraphe III.2 est consacr lanalyse du comportement du spectre basse frquence des
diffrents intervalles de modulation de la modulation RCFM.
Daprs (3.25), les harmoniques de courant sont caractriss par le signal Vcc.Hannf / 3Lacwn,
et sa composante fondamentale par (Va 1/3.Vcc.|Hanf1|) / Lacw1. En se basant sur cette
remarque et en considrant que la tension Vcc est constante, nous avons pu dterminer le
spectre de la fonction Han correspondant un courant ia respectant exactement la norme EN
61000-3-2 (tableau 1.1 du chapitre 1). Ainsi, en basse frquence les spectres rsultant de la
RCFM pour les diffrents intervalles de modulation pourront tre compars au spectre
correspondant la limite des harmoniques impose par la norme EN 61000-3-2 (Fig. 3.3.8).
En basse frquence, le comportement du spectre de la modulation RCFM peut se rsumer en
trois points :
- La composante fondamentale nest pas affecte par la modulation RCFM quelque soit
lintervalle des frquences alatoires utilis,
- Plus lintervalle des frquences de modulation est large, plus le spectre des harmoniques
est dtrior. Cela peut sexpliquer par le fait que plus lintervalle est large plus la
frquence de dcoupage moyenne est infrieure 10 kHz. Ce qui est quivalent au spectre
dune MLI frquence de dcoupage infrieure 10 kHz.
- Par rapport la norme EN 61000-3-2, lutilisation de lintervalle [4kHz ; 10kHz]
engendre le non respect de la norme, et le spectre rsultant de lintervalle [5kHz ; 10kHz]
est la limite de cette norme. Ceci carte donc les intervalles plus larges que
[5kHz ; 10kHz] de la liste des intervalles de frquences de modulation susceptibles dtre
utiliss par le convertisseur.
En tenant compte des rsultats danalyse en haute et basse frquence, nous pouvons dduire
que la plage de frquence des intervalles de modulation adquate pour le bon fonctionnement
du convertisseur est situe dans l'intervalle [5,57 kHz ; 10 kHz]. Toutefois, si le cahier des
charges impose une marge plus importante par rapport aux normes basse frquence, des
intervalles moins larges que [7kHz ; 10kHz] pourront tre utiliss. Dans ce cas, la contrainte
significative concerne la dtrioration du spectre dans la plage de frquence [5 kHz ; 100kHz]
- 99 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
0
Standards EN 61000-3-2
6100-3-2
-10 Theoretical results RCFM ([4kHz ; 10kHz])
RCFM ([5kHz ; 10kHz])
RCFM ([6kHz ; 10kHz])
-20 RCFM ([7kHz ; 10kHz])
RCFM ([8kHz ; 10kHz])
-30 RCFM ([9kHz ; 10kHz])
-40
S(f) [dB]
-50
-60
-70
-80
-90
-100
-110
200 400 600 800 1000 1200 1400 1600 1800 2000
f [Hz]
Lobjectif initial tant dattnuer le spectre haute frquence tout en respectant les normes
basse frquence, nous pensons que le meilleur compromis est lintervalle [6kHz ; 10kHz].
Afin de valider cette tude, nous avons labor des simulations [TOU1-04] ainsi quune
maquette exprimentale de l'architecture de commande qui seront prsentes dans le chapitre
5 de ce mmoire.
Conclusion
Dans ce chapitre, nous avons pu analyser thoriquement le comportement de la modulation
RCFM. Cette analyse a concern ltude de leffet de cette modulation sur le courant de phase
du convertisseur. Leffet sur la tension de sortie tant moins contraignant, nous lavons vrifi
par des simulations [TOU1-04] qui seront prsentes au chapitre 5.
Nous avons procd en trois tapes : la premire tape a consist en la dtermination dune
fonction mathmatique dcrivant le courant de phase. Ensuite, nous avons modlis et
paramtr cette fonction afin de pouvoir calculer sa densit spectrale de puissance. Pour le
calcul de cette densit spectrale, nous avons labor une mthode permettant de prvoir le
- 100 / 188 -
Chapitre 3 : Analyse thorique de la RCFM
comportement des harmoniques rsultants de la modulation RCFM. Cette tude nous a permis
de comprendre deux aspects de cette technique de modulation. Le premier aspect est quen
haute frquence, cette technique est trs avantageuse par rapport la MLI et nimpose pas de
contraintes svres sur la largeur de lintervalle de modulation. Ceci explique son utilisation
massive dans les convertisseurs DC/DC qui nimposent pas de contraintes sur les
harmoniques. Le deuxime aspect concerne son comportement en basse frquence.
Lutilisation de la modulation RCFM naffecte pas la composante fondamentale. Mais, pour
utiliser cette mthode dans des convertisseurs AC/DC prlvement sinusodal du courant,
une tude pralable du convertisseur et de la largeur de lintervalle de modulation simpose.
Le comportement des harmoniques rsultant de la RCFM est extrmement li la structure du
convertisseur et aux frquences alatoires utilises.
Nous pensons que la dmarche danalyse que nous avons dvelopp dans ce chapitre est
gnrique et peut tre applicable tous les convertisseurs AC/DC. Toutefois, la premire
tape concernant la modlisation et le paramtrage des fonctions de modulation reste propre
chaque convertisseur.
En utilisant la modulation RCFM, on amliore la qualit du courant prlev tout en rduisant
la frquence de commutation, ce qui rduit les pertes dans les IGBT. Pour continuer sur cette
voie damlioration, nous avons labor une architecture du convertisseur Cuk utilisant moins
de capteurs et amliorant sa fiabilit. Cette architecture est prsente dans le chapitre 4.
- 101 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
- 102 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
Introduction
La relation (4.1) implique que si lon mesure le courant icoup un instant donn, le courant
mesur correspond ncessairement lun des courants ia, ib, ic ou iLdc. Le principe de la
- 103 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
= i Ldc i ch arg e
dvdc
Cdc (4.2)
dt
Le tableau 4.1 rcapitule les sept vecteurs de commande possibles dans le cas du
convertisseur Cuk et le courant associ chaque vecteur. Les vecteurs V1 V6 sont appels
vecteurs actifs, ils permettent destimer les courants de phase. Au sein de chaque priode de
dcoupage, seuls deux vecteurs actifs sont gnrs. Dans lexemple de la figure 4.1.2(a)
(vaca > vacb > vacc), les courants de phase concerns sont ia et ic. Ltat du convertisseur
pendant les dures T1 et T2 respectivement des vecteurs V1(1000) et V2(1100) est illustr dans
les figures 4.1.2(b-c). Lensemble des vecteurs actifs dfinit six zones de modulation (figure
4.1.3) similaires aux zones de modulation dans le redresseur triphas Boost.
Le vecteur V7(0001) est appel vecteur zro, parce que pendant la dure de cet tat, les bras
du pont triphas sont court-circuits (figure 4.1.2(d)) et icoup = -iLdc.
Hz t
-iLdc -iLdc
ia -iLdc
-ic -iLdc ia
-ic
icoup
t
T1 T2 T3 T4 T1 T2 T3 T4
- 104 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
Figure 4.1.2 : Exemple de correspondance entre le courant icoup et les diffrents courants
estims en fonction de la commande
0100 1100
II
III I
0110 1000
IV VI
V
0010 1010
Figure 4.1.3 : Diagramme vectoriel de commande du convertisseur Cuk
(six zones de modulation)
- 105 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
Les tensions de rfrence de la modulation ntant pas constantes pendant chaque priode de
dcoupage, la symtrie des vecteurs de commande dans la figure 4.1.2(a) nest pas toujours
valide. Toutefois, pendant chaque demi priode de dcoupage, deux vecteurs actifs et deux
vecteurs zros sont produits. Donc, nous pouvons reconstruire les courants de phase et la
tension de sortie la fin de chaque demi priode de dcoupage (ou au dbut de la demi
priode suivante) [LEE-02].
i b ( t k ) = ( ia ( t k ) + ic ( t k ) )
T3 Vcc ( t1 )
avec Vaca = la valeur moyenne de la tension Vaca calcule sur [ t1 ; t k ] .
T3 + T4 3
Lune des tensions Vac ( a ,b,c) intervient dans le calcul du courant de phase associ au premier
- 106 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
De plus, les dures des vecteurs de commande doivent tre connues la fin de chaque demi
priode de dcoupage pour pouvoir estimer correctement les tensions Vac ( a ,b,c) partir de Vcc.
Toutefois, malgr son utilisation dans [NOG-98], [LEE-02] et [BOH-97], elle reste une
mthode trs exigeante en terme de temps de calcul.
- 107 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
la premire acquisition associe ce vecteur. Le principe de cette mthode est illustr dans la
figure 4.2.1. Tacq est la priode dacquisition du courant icoup. La gnration des trois
courants de phase estims est effectue la fin de chaque demi priode de dcoupage.
icoup[A] icoup = ia
icoup = -ic
icoup = -iLdc
1 1
V1 V2
0
0
1
0
0 0
tk 1 2 tk + T/2 t
Tacq
Figure 4.2.1 : Acquisition du courant icoup lors de la transition entre deux vecteurs actifs
- 108 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
Lavantage de cette mthode est quelle ne ncessite quun capteur de courant icoup, et ne
ncessite quasiment pas de calcul des courants de phase. Le tableau 4.3 prsente une
comparaison de ces deux mthodes.
Pour les deux mthodes, les trois courants de phase sont calculs la fin de chaque demi
priode de dcoupage, et ces courants estims gardent les valeurs calcules pendant toute la
dure de la demi priode de dcoupage suivante.
Aprs cette brve comparaison des deux mthodes, nous avons choisi dutiliser la mthode
consistant estimer les courants de phase la transition entre deux vecteurs actifs. Lobjectif
tant dlaborer un systme de commande fonctionnant dans un environnement temps rel,
cette mthode est parfaitement adapte au cahier des charges car elle ne ncessite pas de
calculs supplmentaires par rapport au systme de commande (chapitre 2). De plus, elle
nutilise quun seul capteur de mesure, ce qui permet de rduire le cot du convertisseur.
30 30
20 20
10 10
Icoup [A]
Icoup [A]
0 0
-10 -10
-20 -20
A B A B C
-30
C -30
0.3457 0.3457 0.3457 0.3457 0.3457 0.3458 0.3458 0.3458 0.3458 0.3458 0.3457 0.3457 0.3457 0.3457 0.3457 0.3458 0.3458 0.3458 0.3458 0.3458
[s] [s]
Afin de calquer le courant icoup la thorie, nous avons introduit un filtre passe bas de
Butterworth (ordre 2). La frquence de coupure de ce filtre est de 300 kHz. Dans la figure
- 109 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
4.2.2(b), nous constatons que ce filtrage empche les oscillations transitoires du courant icoup.
Ce qui permet damliorer la qualit des courants estims.
De plus, dans le cas du processeur TMS320LF2407, le temps minimal entre deux acquisitions
successives est de lordre de 500 ns. Compte tenu du temps ncessaire pour lexcution de
lalgorithme du systme de commande, le retard des impulsions de commande d aux drivers
des IGBT et le retard engendr par le filtrage du courant icoup, nous estimons la dure Tacq
minimale de lordre de 4 s. La contrainte lie cette priode dacquisition est que le courant
icoup risque de ne pas tre prlev pendant certains vecteurs de commande. Surtout lorsque la
dure du vecteur actif est infrieure Tacq.
0100 1100
Ha
t
II Hb
t
III Hc
I t
Hz t
0110 1000
-iLdc -iLdc
IV VI -iLdc
ia -iLdc ia
icoup
V
T1 T2 T4 T1 T2 T4
0010 t
1010 (tk-1) (tk) (tk+1)
Dans la figure 4.2.3, les parties hachures du diagramme vectoriel de commande illustrent les
zones o risquent de se produire un tel dysfonctionnement. Ces zones sont calcules partir
de la largeur des enveloppes sinusodales des tensions de rfrence de la modulation Vac* ( a ,b,c) .
- 110 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
300
* * *
Vaca Vacb Vacc
200
100
[V]
0
-100
-200
-300
0.5 0.502 0.504 0.506 0.508 0.51 0.512 0.514 0.516 0.518 0.52
T [s]
Figure 4.2.4 : Les tensions de rfrence de la modulation Vac* ( a ,b,c) et les zones susceptibles
dengendrer des estimations errones du courants de phase
Par rapport au courant rel, ces erreurs destimation sont intolrables. Pour contourner cette
anomalie dans lestimation, nous avons tudi et labor deux mthodes de correction des
courant estims. La premire mthode consiste corriger directement les vecteurs de
commande, en agissant sur la largeur du vecteur en fonction de son angle polaire dans le
diagramme de commande vectoriel [TOU2-04]. La deuxime mthode est base sur la
- 111 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
Cette mthode est inspire des travaux de [BLA-97] et [LEE-02]. Son principe consiste fixer
la largeur minimale dun vecteur de commande Tacq. Etant donne que dans les phases de
fonctionnement critiques, deux fonctions de modulation se retrouvent superposes, la
correction consiste laisser lune des deux inchange et largir ou rtrcir lautre en fonction
de langle polaire du vecteur de commande.
4 3
5 2
6 1
7 12
8 11
9 10
Afin de pouvoir suivre les angles des vecteurs de commande et agir sur les fonctions de
modulation, le diagramme vectoriel est dcoup en douze secteurs (figure 4.2.6). La prsence
de langle dans chaque secteur spcifie la fonction de modulation sur laquelle on doit agir.
Si lon se place dans le cas de la figure 4.2.7, un angle de commande dans la zone 1 signifie
que la fonction Hb doit tre plus large dau moins Tacq par rapport Hc_bst. La prsence de cet
angle dans la zone 2 implique que Ha doit tre au moins Tacq plus large que Hc_bst pour
pouvoir estimer le courant de la phase a. Le tableau 4.4 rcapitule le principe de la
correction des fonctions de modulation dans chaque secteur. Une fonction de modulation de
rfrence signifie la rfrence par rapport laquelle on agit (largir ou rtrcir Tacq) sur la
fonction de modulation corriger dans un secteur donn. Le but est que la somme des dures
des deux vecteurs actifs aprs correction reste identique quavant la correction.
- 112 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
= 60
(1100)
2
T2 V = 30
1
T2' V
= 0
T1 T ' (1000)
1
Figure 4.2.7 : La diffrence en terme de correction des vecteurs actifs entre les deux secteurs
(1) et (2), malgr que les deux secteurs appartiennent la mme zone de commande (I)
Fonction Fonction de
Secteur Action
corriger rfrence
1 Hb Hc_bst largir Hb de Tacq
2 Hb Ha_bst rtrcir Hb de Tacq
3 Ha Hb_bst rtrcir Ha de Tacq
4 Ha Hc_bst largir Ha de Tacq
5 Hc Ha_bst largir Hc de Tacq
6 Hc Hb_bst rtrcir Hc de Tacq
7 Hb Hc_bst rtrcir Hb de Tacq
8 Hb Ha_bst largir Hb de Tacq
9 Ha Hb_bst largir Ha de Tacq
10 Ha Hc_bst rtrcir Ha de Tacq
11 Hc Ha_bst rtrcir Hc de Tacq
12 Hc Hb_bst largir Hc de Tacq
Algorithme de dtection
i [1 ; 12]
Dtecter le secteur i dans le diagramme dun secteur dans la MLI
vectoriel. vectorielle.
[D-3643]
Gnrer les trois fonctions Gnrer les trois fonctions Tableau 4.4
de modulation relles de modulation corriges
- 113 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
Lalgorithme simplifi de la figure 4.2.8 donne une vue gnrale de la mthode de correction
des vecteurs actifs.
Lalgorithme de la MLI vectorielle (indispensable dans cette mthode) gnrant les fonctions
de modulation corriges et la connaissance en temps rel de langle polaire sont
relativement complexes et ncessitent lexcution dun nombre important dinstructions [D-
3643]. Ce qui reprsente un handicap majeur de cette mthode. De plus, lorsque lon corrige
un vecteur de commande un angle i donn, le signal modul correspond un angle polaire
antrieur j, ce qui risque dengendrer la dformation des courants de phases estims.
Cette mthode est base sur le principe que le systme triphas ( i a _ est , i b _ est , i c _ est ) est
II.3.2. Correction de la trajectoire du courant
- 114 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
transforme de Park du systme triphas ( i a _ est , i b _ est , i c _ est ) la fin de chaque demi priode
( )
(4.5)
Etant donn que les tensions Vac* ( a ,b,c) ont des enveloppes sinusodales pouvant atteindre jusqu
30 V, nous appliquons cette mthode lorsque la valeur absolue de la diffrence entre deux
tensions de rfrence est infrieure ou gale 30 V. Lalgorithme de la figure 4.2.10 dcrit le
principe de la mthode et son intgration dans la procdure destimation de courants dcrite
prcdemment.
Donnes
La composante id_est(tk-1) de la transforme de
Park du systme triphas (ia_est, ib_est, ic_est)
*
Vaca Vacb
*
30 ou
Vacc 30 ou
* * Non
Vacb
*
Vacc Vaca
*
30
tk = tk +
Tdcoupage Gnrer les fonctions de
2 modulations Ha, Hb, Hc et Hz
Oui
Lavantage de cette mthode par rapport la prcdente est que de part sa simplicit, elle ne
modifie pas les fonctions de modulation, et engendre moins de calcul que lestimation des
courants de phase partir des mesures de icoup. De plus, les dformations du courant estim
sont moins importantes compares la premire mthode (figure 4.2.11).
- 115 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
courant estim, et rduit ainsi son retard par rapport au courant rel (retard d au temps de
calcul et lacquisition) [MOY-91]. Le comportement dynamique de cet observateur doit tre
plus rapide que le processus rel, de manire ce que lerreur dobservation soit aussi faible
que possible, mme dans les phnomnes transitoires.
- 116 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
R ac
w
i d Lac
.
v vacd
= d + 1 d
i
iq w R i Lac v q vacq
ac q
(4.6)
L ac
w 0
L ac Lac
avec
R ac T
1 L T wT L 0
vd ( k ) v acd ( k )
F= , G= , U ( k ) = et T =
vq ( k ) v acq ( k )
Tdcoupage
R ac T
ac ac
wT 1 T 0 Lac
2
Lac
( )
Ainsi lquation de lobservateur scrit :
( k + 1) = F.X
X ( k ) + G.U ( k ) + K Y ( k ) Y
(k) (4.9)
E
e rr ( k + 1) = [ F K E C] e rr ( k ) (4.10)
Nous remarquons que l'erreur est indpendante de lentre du systme et ne dpend que de
matrice [ F K E C] . Si ces valeurs propres sont toutes parties relles ngatives, on aura
l'erreur initiale. La dynamique de lerreur dobservation est rgie par les valeurs propres de la
lim e rr ( t ) = 0 (4.11)
t
En dautres termes, X convergera vers X. De plus, pour que lerreur destimation puisse
converger rapidement vers zro, les valeurs propres doivent tre situes lintrieur du cercle
- 117 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
unitaire de centre 0 [LEE-02] [BH-83]. Aussi, les ples de lobservateur doivent tre plus
rapides que ceux de la boucle de courant sans toutefois augmenter la bande passante du
systme (rendant celui-ci trs sensible au bruit de mesures).
Le vecteur gain KE de lobservateur peut tre dtermin de diffrentes faons (estimation
optimale - filtre de Kalman, placement des ples, compagne observable) pour que
det [sI A + K E C ] = 0
lquation caractristique de (4.10) possde les deux ples dsirs, soit :
(4.12)
10
Courant reconstruit
5
Ia [A]
-5
-10
-15
0.15 0.16 0.17 0.18 0.19 0.2
[s]
- 118 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
= i Ldc i ch arg e
dvdc
Cdc (4.13)
dt
Si lon utilise le mme principe de discrtisation que pour le courant de phase (demi priode
de dcoupage), la relation (4.13) devient :
v dc ( t k +1 ) = v dc ( t k ) + i Ldc ( t k ) i ch arg e ( t k )
Cdc
T
(4.14)
Daprs (4.14), la reconstruction de la tension dpend de sa valeur initiale et des courants iLdc
et icharge. Le courant de charge icharge tant mesurable, la reconstruction de la tension de sortie
dpend de lestimation du courant iLdc. Comme nous lavons signal dans la figure 4.1.2(a), le
courant icoup mesur pendant le vecteur nul (Hz = 1) est gale au courant iLdc. Dans le courant
mesur icoup de la figure 4.3.1, nous constatons des surintensits lors des transitions des
vecteurs actifs aux vecteurs zro. Les endroits o le courant icoup concide exactement avec le
courant iLdc sont les centres des vecteurs nuls. Cette remarque est trs avantageuse pour la
reconstruction de la tension de sortie lorsque la frquence de modulation est fixe. Parce
quavec une frquence dacquisition multiple de la frquence de modulation, le courant icoup
est mesur au voisinage de chaque instant tk. Cependant, dans le cas de la modulation RCFM,
lacquisition du courant icoup correspondant iLdc ncessite la dtection du vecteur zro et
augmente par consquent le temps de calcul.
- 119 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
Figure 4.3.1 : Correspondance entre les courants icoup, iLdc et les fonctions de modulation
1
0 Hd Hq
1 1
v cc v Cc
- Hz 0
d 1 3 id
Cc Cc
iL = i Ldc + 0 0
dt dc L dc
cc
R dc 1
v 2 iq
Hz - - (4.15)
v dc dc
Ldc L dc
0 0 0
1 1
-
Cdc Cdc R Load
- 120 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
Afin de traiter un systme invariant dans le temps, la fonction Hz est approxime par sa
3 3
composante continue 1- m (Chapitre 2). Le courant iLdc tant la seule grandeur
2
mesurable, il est considr l'unique sortie du systme. Lquation de lobservateur dtat
( )
sexprime alors
( k + 1) = F.X
X ( k ) + G .U ( k ) + K Y ( k ) Y
(k)
(k) = C X (k)
V V V V EV V V
(4.16)
Yv v
avec
T
Hd Hq
1 T
Cc
1 - HzT 0
1 3 id ( k )
C Cc
0 , U ( k ) =
c
F= , G= 0
iq ( k )
R T
Ldc 2
H z T 1- dc T -
Ldc Ldc
0 0
T T
0 1-
Cdc Cdc R Load
vcc ( k ) 0 0 0
X V ( k ) = i Ldc ( k ) , Yv ( k ) = i Ldc ( k ) et CV = 0 1 0
v (k) 0 0 0
dc
Lutilisation de lobservateur dtat pour la reconstruction de la tension de sortie conduit des
rsultats semblables la premire mthode en rgime permanent. Toutefois, lors de la mise
sous tension, si lon dmarre le convertisseur avec des condensateurs non chargs
lalgorithme de lobservateur diverge. De plus, le temps de calcul ncessaire pour cette
mthode est trois fois suprieur par rapport la premire mthode. Pour ces raisons, nous
avons opt pour lutilisation de la premire mthode de reconstruction. Thoriquement, la
reconstruction de la tension vdc est relativement simple compare la reconstruction du
courant de phase. Cependant, certaines contraintes imposent des modifications dans la relation
(4.14) et dans lalgorithme de reconstruction.
- 121 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
350
300
Consigne de tension
250
[V]dc
Tension estime
200
150
(a) Tension estime avec une acquisition des courants la fin de chaque
demi priode de dcoupage
303
Tension estime
302
301
[V]dc
300
299
298
Tension mesure
297
(b) Tension estime avec une acquisition des courants lors de chaque vecteur nul
Nous constatons que la tension estime de la figure 4.3.2(a) diverge de la tension consigne
cause dune seule acquisition errone du courant icoup. En effet, la variation des tensions de
rfrence de la modulation Vac* ( a,b,c) a engendr une translation du vecteur zro, et la fin de la
- 122 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
icoup et icharge quivalent au temps sparant deux acquisitions successives dans le processeur
TMS320LF2407. Ceci explique en partie le dcalage entre la tension estime et la tension
relle. Aussi, plus linstant tacq est proche de tk, plus la tension estime converge vers la
tension relle.
Un autre aspect dans la reconstruction de la tension de sortie est si la dure du vecteur nul est
infrieure Tacq, cela ncessite un traitement semblable la reconstruction du courant de
phase. Dans ce cas de figure deux solutions se prsentent : la premire consiste corriger les
trois fonctions de modulation, et la seconde solution consiste affecter vdc(tk-1) directement
vdc(tk). Cette deuxime solution est la plus approprie, et elle est justifie par le fait que seule
la partie purement continue de vdc intervient dans le systme dasservissement. Dans
lestimation de la tension (figure 4.3.2(b)), nous avons pris en compte cette deuxime solution
lorsque Tvecteur_nul < Tacq. Nous constatons que malgr limprcision de lestimateur, la tension
de sortie mesure nest pas dtriore. Ceci implique la convergence de lalgorithme de
reconstruction de la tension vdc.
- 123 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
300
200
100
-100
-200
-300
tn tn+1 tn+2 tn+3 tn+4 tn+5 tn+6 tn+7 tn+8 tn+9 tn+10 tn+11 tn+12 tn+13 tn+14
0.01 0.0101 0.0102 0.0103 0.0104 0.0105 0.0106 0.0107 0.0108 0.0109 0.011
[s]
Figure 4.4.1 : Correspondance entre les instants tk et les priodes de dcoupage alatoires
20
Courant rel
15
10
5
Courant estim
Ia [A]
-5
-10
-15
-20
- 124 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
15
Courant estim
10
5
Ia [A]
Courant rel
-5
-10
-15
- 125 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
302.5
Tension estime
302
301.5
301
300.5
[V]dc
300
299.5
299
298.5
298
Tension relle
297.5
0.33 0.335 0.34 0.345 0.35
[s]
V. Algorithmes de reconstruction
Dans cette tude de la nouvelle architecture du convertisseur Cuk, nous avons trait
sparment les algorithmes de reconstruction des courants alternatifs et de la tension continue.
Toutefois, nous navons intgr ni laspect temporel ni linterdpendance entre ces
algorithmes. Cette section est consacre la prsentation de lalgorithme gnral de la
reconstruction, ainsi qu son excution chaque cycle de calcul.
La figure 4.5.1 illustre lalgorithme gnral de lestimation des signaux non mesurs. Le
Tend correspond au dernier instant dune demi priode de dcoupage. Les trois courants de
phases sont gnrs au dbut de chaque demi priode de dcoupage et gardent les valeurs
calcules pendant toute la dure de cette priode. Le chargement des entres du systme de
commande consiste charger au dbut de chaque cycle de calcul les informations suivantes :
- 126 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
- 127 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
Estimation de la tension de
sortie
*
Vaca *
- Vacb 30 ou
*
Vacb *
- Vacc 30 ou Non 1er vecteur Oui
30
* *
actif ? Acquisition de iLoad et icoup pour
Vacc - Vaca
la reconstruction de la priode
Non suivante
Oui
Pour rduire considrablement le temps de calcul, lacquisition de tous les signaux est
programmable. C'est--dire, que lacquisition dun courant est effectue sur demande de
lalgorithme de calcul. Par exemple, lorsque la relation (4.5) est utilise pour estimer le
courant de phase, aucune acquisition de icoup pendant les deux vecteurs actifs nest effectue.
Lalgorithme de reconstruction est mis jour la fin de chaque demi priode de dcoupage.
Cette mise jour consiste en la dfinition du nouvel instant de fin de la demi priode Tend et
- 128 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
en linitialisation des drapeaux spcifiant lordre dapparition des vecteurs actifs. La figure
4.5.2 illustre cette procdure de mise jour.
Une acquisition de
Tend = Tend +
Non Tdcoupage
iLoad et icoup effectue ?
2
Oui
v dc ( k + 1) = v dc ( k )
Avant dinitialiser les drapeaux et linstant Tend, la premire tape de cette procdure est de
gnrer les courants de phase et la tension de sortie. Les courants de phase tant calculs
durant les vecteurs actifs, leur gnration nengendre pratiquement pas de temps de calcul
(une seule opration daffectation). Cependant la tension de sortie, si elle nest pas calcule
pendant le dernier vecteur nul, elle doit tre estime avant la mise jour. Si aucune
acquisition ne permet cette estimation, le systme gnre cette tension en utilisant la valeur
prcdente de la tension vdc.
La modulation frquence de dcoupage alatoire impose un chantillonnage suivant une
frquence alatoire. Donc au dmarrage du systme, la premire priode de dcoupage est
inconnue. Pour pallier ce problme, nous dmarrons le systme avec une frquence de
dcoupage correspondant la frquence maximale de la gamme de frquence utilise.
Ensuite, une fois la premire priode alatoire connue, le systme prend pour rfrence la fin
de cette priode Tend.
- 129 / 188 -
Chapitre 4 : Reconstruction du courant de phase & de la tension de sortie
Conclusion
Certes, la reconstruction des signaux rduit le cot associ au matriel, mais au dtriment
dune augmentation de la complexit du systme de commande. Cette complexit se traduit
par des algorithmes assez lourds mettre en place, et qui sont souvent inexploitables dans les
systmes temps rels. Dans la conception de ces algorithmes, la grande difficult rside dans
le choix faire entre la prcision et le temps de calcul. Le temps de calcul doit tre le plus
petit possible pour ne pas perturber lalgorithme du systme de commande, et la prcision doit
tre la plus importante possible pour ne pas dtriorer la qualit du convertisseur par rapport
sont tat en prsence des capteurs de mesure. Donc, le compromis temps de calcul / prcision
reprsente le critre principal du choix des mthodes retenues pour la reconstruction.
Dans ce chapitre nous avons tudi et compar plusieurs mthodes destimation et de
correction du courant de phase et de la tension de sortie. Pour la reconstruction du courant,
nous avons prsent le principe de lestimateur, et nous avons tudi diffrentes mthodes
pour pallier les dfauts dacquisition du courant icoup. La solution retenue est lestimation du
courant la transition entre deux vecteurs actifs avec correction de la trajectoire du courant.
Afin damliorer la dynamique de lestimateur du courant, nous avons introduit un
observateur des variables dtat. En plus de lamlioration de la dynamique du systme, cet
observateur amliore la forme du signal reconstruit en introduisant un effet de filtrage.
Dans la reconstruction de la tension de sortie, aprs la description du principe de la mthode
destimation, nous avons tudi et rsolu les contraintes dimplmentation de cette mthode.
Contraintes lies principalement la dtection du vecteur nul et aux dfauts dacquisition.
Nous avons ensuite tudi linfluence du dcoupage alatoire sur les mthodes de
reconstruction. La modulation RCFM distribution discrte permet une reconstruction des
signaux sans effets ngatifs lorsque lensemble des frquences alatoires est judicieusement
choisi.
Afin de complter ltude dveloppe dans ce chapitre ainsi que dans les trois premiers, nous
consacrons le cinquime chapitre aux simulations et ltude exprimentale de cette nouvelle
architecture du convertisseur Cuk.
- 130 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
- 131 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Introduction
Comme dans la majorit des systmes lectriques base de convertisseurs statiques, ltude
analytique des systmes de commande est difficile mener sans lutilisation de plusieurs
hypothses simplificatrices souvent trs rductrices. Etude analytique dautant plus complexe
lorsque les frquences dchantillonnage et de dcoupage sont alatoires. Dans ce domaine, la
simulation numrique constitue un outil dtude, dvaluation et de prconception
incontournable et trs efficace. Surtout lorsque les simulations sont labores en tenant
compte des contraintes exprimentales.
Dans ce chapitre, avant de simuler le fonctionnement du convertisseur, nous proposons de
valider exprimentalement ltude thorique de la modulation RCFM. Cette modulation
engendre des spectres continus avec la possibilit dapparition dharmoniques en basse
frquence. Ceci impose une tude pralable des diffrentes mthodes numriques destimation
de la densit spectrale de puissance, pour aboutir au choix de la mthode la plus adapte
lanalyse de cette modulation. Pour cette raison, nous avons consacr la premire partie de ce
chapitre ltude comparative de ces mthodes, la seconde partie porte sur la prsentation de
la maquette exprimentale validant lanalyse thorique de la RCFM ainsi que les rsultats des
analyses spectrales. La troisime partie de ce chapitre est consacre la simulation des
modles labors dans les chapitres 2, 3 et 4. Enfin la dernire partie, nous valuons les
contraintes exprimentales lies au systme de commande ainsi qu leur effet sur les rsultats
de simulation.
- 132 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
- 133 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Dans cette partie, nous allons tudier et comparer les mthodes non paramtriques, afin de
dterminer la mthode la plus adapte pour lanalyse des fonctions de modulation du
convertisseur Cuk. Mais pour cela, il faut disposer dun signal alatoire dont on connat
parfaitement la densit spectrale thorique, et qui prsente des variations rapides comparables
la fonction Han.
Le signal alatoire variations rapides est celui de la figure 5.1.1. Cest un signal carr
frquence fixe et amplitude variable alatoirement entre -1 et 1. La PSD de ce signal a t
dtermine thoriquement dans les travaux de [FON-00]. Pour une loi de distribution
uniforme, cette PSD thorique sexprime sous la forme (5.1)
S ( f ) = .sin c ( f .) .
2 1
(5.1)
3.T
sin ( .x )
avec la largeur de limpulsion au sein de chaque priode T, et sinc(x) le sinus cardinal de x
0.8
0.6
0.4
0.2
-0.2
-0.4
-0.6
-0.8
-1
0.1 0.1005 0.101 0.1015 0.102 0.1025 0.103 0.1035 0.104 0.1045
[s]
Laspect alatoire et les variations rapides caractrisant ce signal sont similaires aux
caractristiques de la fonction Han. La PSD thorique de ce signal est donne dans la figure
5.1.2 et les spectres rsultants des diffrentes mthodes non paramtriques destimation de la
PSD seront compars ce spectre thorique. Les critres de comparaison des mthodes sont :
la dispersion du spectre autour de sa moyenne, la rsolution et les fuites spectrales de chaque
mthode.
- 134 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
-40
-50
-60
-70
S(f) [dB]
-80
-90
-100
-110
-120
0 1 2 3 4 5 6 7 8 9 10
F [Hz] 4
x 10
Initialement, cette mthode a t utilise pour dtecter une priodicit cache dans un signal
bruit do le nom de "Periodogram". Dans cette mthode, la PSD dun signal numrique
xL(n) de longueur L est estime par la relation (5.2).
Ses ( f ) = F . X L ( f k )
L 2
X ( f ) = 1
L( )
ech
2
(5.2)
L 1 j. .k.n
L k
N
x n .e
N n =0
- 135 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
(a) (b)
(c) (d)
Figure 5.1.3 : Les PSD estimes par la mthode Periodogram suivant quatre fentres de
pondration, et comparaison avec la PSD thorique
Le choix du fentrage de pondration ninfluence pas la dispersion du spectre, qui reste trs
importante et atteint les 20 dB (figure 5.1.3). La rsolution lors de lutilisation dune fentre
rectangulaire (Boxcar window) est approximativement de Fech/L. Dans le cas de la fentre de
Hamming, les fuites spectrales sont minimises par rapport la premire fentre mais la
rsolution passe environ 2.Fech/L. Pour les deux fentres restantes, la rsolution est comprise
entre Fech/L et 2.Fech/L.
Nous avons vu que la mthode Periodogramme conduisait une dispersion trs importante du
spectre. De plus, lutilisation dune fentre de pondration rduit les fuites spectrales au
dtriment de la rsolution.
La dispersion peut tre rduite en moyennant K priodogrammes de K ralisations diffrentes.
Afin de raliser cette opration de moyennage, la mthode de Welch consiste dcouper
lenregistrement de longueur L (200 ms) en K segments de longueur N (20 ms). Ces segments
peuvent ventuellement se chevaucher (figure 5.1.4).
- 136 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
N N N
P
t
(a) (b)
(c) (d)
Figure 5.1.5 : Les PSD estimes par la mthode de Welch dans le cas dun recouvrement
P = 0% pour diffrentes fentres de pondration, et comparaison avec la PSD thorique
Chaque segment est pondr par une fentre wn. Ainsi, cette mthode permet de gagner un
facteur L/N sur la dispersion (non recouvrement) ou par exemple 1,5.L/N en cas de
recouvrement 50%.
- 137 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
NP
Le taux de recouvrement des segments sexprime comme suit : 100.
N
Lorsque le recouvrement est gal 0%, la mthode de Welch consiste moyenner K
priodogrammes de K ralisations distinctes du signal. Ceci permet de rduire la dispersion
du spectre comme nous pouvons le constater dans la figure 5.1.5. Cependant, la fentre de
pondration ne tient pas compte des chantillons situs sur ses bords, ce qui introduit une
perte dinformations. Do lintrt du recouvrement des fentres de pondration. A noter que
la PSD estime concide avec la PSD thorique translate de 3 dB.
(a) (b)
(c) (d)
Figure 5.1.6 : Les PSD estimes par la mthode de Welch dans le cas dun recouvrement
P = 50% pour diffrentes fentres de pondration, et comparaison avec la PSD thorique
Dans la figure 5.1.6, nous avons estim la PSD avec un recouvrement de 50%. En partant
dun recouvrement de 0% jusqu 50%, plus le recouvrement est important plus le spectre est
"liss", et partir de 50% la dispersion augmente. Dans tous les cas de recouvrement,
lutilisation de la mthode de Welch permet destimer une PSD plus lisse et qui colle avec le
spectre thorique translat de 3 dB (due loption de calcul Matlab).
La dispersion est plus importante dans le cas de la fentre rectangulaire par rapport la
fentre de Blackman. Lutilisation des deux fentres Hamming et Bartlett aboutit une
estimation identique de la PSD, et de point de vue dispersion elles sont intermdiaires entre la
- 138 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
La critique que nous pouvons apporter la mthode de Welch est la perte de rsolution. Car
on dispose dun enregistrement de longueur L et on peut esprer une rsolution de lordre de
Fech/L ; mais le dcoupage en K segments de longueur N conduit une rsolution de lordre
de Fech/N (L/N fois moins bonne).
Lide de la mthode "Multitaper" est dutiliser K fentres de pondration diffrentes (notes
w(n), n variant de 0 K-1) de faon obtenir K signaux indpendants (de longueur L) :
y (n ) = x.w (n ) (5.3)
fentres de pondration :
c w ( k ) w ( j) = 0 pour k j
c w ( k ) w ( j) = 1 pour k = j
(5.4)
La solution propose par Matlab pour satisfaire cette condition sont les fentres rsultantes de
la fonction donde sphrodale prolate (Prolate Spheroidal Wave Function). Ces fentres
dpendent de deux paramtres (N : le nombre des points et B : la bande passante) et sont telles
que :
D ( f ).W( ) ( N, B, ) .d = ( ) ( N, B) .W( ) ( N, B, f )
B
N (5.5)
B
k k k
sin ( N..f )
o D N ( f ) = e
N 1
j.2 .
sin ( .f )
.f
2
. le noyau de Dirichlet (transforme de Fourier dune fentre
- 139 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Lexemple de la figure 5.1.7 illustre les fentres de pondration w(k) pour N = 2.104 et k
variant de 0 7. Le nombre de fentres K est li la bande passante B par la relation :
K = 2.N.B 1 (5.6)
Typiquement, on choisit la bande passante (rsolution) B comme un multiple de la rsolution
de la FFT, c'est--dire quon choisit le produit N.B (2, 5/2, 3, 7/2, 4). Dans Matlab, le produit
N.B est not NW.
a (f )
K 1
n n
Ses ( f ) =
.Ses y
k =0 (k)
a
k
K 1
(5.7)
k =0
k
ak sont des cfficients de pondration (ak = (k)). Par exemple le tableau 5.1 illustre les (k) des
fentres de pondration de la figure 5.1.7.
k 0 1 2 3 4 5 6 7
(k) 1 1 1 0.99997 0.99941 0.9925 0.93665 0.69884
Tableau 5.1 : Les valeurs de (k) associes chaque fentre de pondration
Nous navons pas pu utiliser cette mthode pour analyser lenregistrement de longueur 106
points et de frquence dchantillonnage 5 MHz pour insuffisance despace mmoire (out of
Memory dans un Pentium IV ayant 512 Mo de mmoire RAM). Pour tudier le spectre
rsultant de cette mthode, nous avons rduit la frquence dchantillonnage 1 MHz et la
- 140 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
longueur du signal 2.105 chantillons (afin de garder un enregistrement de 200 ms). Les
rsultats de cette mthode pour diffrents NW sont prsents dans la figure 5.1.8.
Laugmentation du nombre de fentres K minimise la dispersion du spectre, mais avec une
lgre augmentation des fuites spectrales (qui restent ngligeables par rapport la mthode de
Welch). Laugmentation de la bande passante B conduit une faible rsolution (compare la
mthode Periodogram), mais le spectre est mieux "liss". Indpendamment des rsultats
danalyse, la mthode Multitaper ncessite un espace mmoire trs important et engendre des
temps de calcul trs levs (45 fois plus suprieur la mthode de Welch).
(a) (b)
(c) (d)
Figure 5.1.8 : Les PSD estimes par la mthode Multitaper pour diffrentes valeurs de NW,
et comparaison avec la PSD thorique
Aprs cette brve description des mthodes non paramtriques, nous pensons que le meilleur
compromis dispersion / rsolution / fuite spectrale est la mthode de Welch avec un
recouvrement de 50% et une fentre de pondration de Blackman. Ce choix nest valable que
pour les signaux alatoires spectre continu (spectre haute frquence de la modulation
RCFM).
Le spectre basse frquence de la fonction Han, est constitu dun spectre continu et dun
spectre discret (les harmoniques basse frquence). Le spectre continu est exprim en V2/Hz,
alors que les harmoniques sont donns en V2. Dans ce cas, la rsolution et les fuites spectrales
- 141 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
sont les paramtres les plus importants dans le choix de la mthode destimation de la PSD.
Pour cette raison, nous avons choisi la mthode Periodogram moyenne sur cinq acquisitions
distinctes (Fech = 5 MHz, L = 106 points) avec un fentrage rectangulaire pour lestimation de
la PSD en basse frquence.
Afin de valider exprimentalement ltude thorique du chapitre 3, nous avons mis en place
une maquette exprimentale (figure 5.1.9) pour gnrer les trois fonctions Han, Hbn et Hcn. En
plus de la validation de la thorie, lobjectif de ce dispositif est de vrifier le comportement de
la modulation RCFM dans le convertisseur Cuk avant son implmentation dans le systme de
commande.
Cette maquette est constitue principalement du processeur de signal TMS320LF2407 et du
circuit logique programmable ISPGAL22LV10. Le processeur gnre les fonctions de
modulation Ha_bst, Hb_bst et Hc_bst (sorties PWM) ainsi que les trois fonctions Cab, Cbc et Cca
(sorties numriques) de lalgorithme dadaptation des signaux de commande (chapitre 2).
Lalgorithme implment dans le circuit logique permet de gnrer et synchroniser les trois
fonctions de modulation Ha, Hb, Hc partir des signaux provenant du processeur. Une fois ces
fonctions de modulation enregistres dans le systme dacquisition, les trois fonctions Han,
Hbn et Hcn sont facilement dtermines dans Matlab en utilisant les relations de (5.8).
Lvaluation de la PSD de ces fonctions est effectue en utilisant les mthodes dcrites dans
la section prcdente.
H an = 2.H a H b H c
H bn = 2.H b H a H c (5.8)
H cn = 2.H c H b H a
Le systme dacquisition utilis, a une rsolution de 8 bits (0,4 %), une profondeur mmoire
de 106 mots et une frquence dchantillonnage de 5 MHz. Avec ces caractristiques du
systme dacquisition, et tant donnes que les fonctions H(a,b,c) sont mesures sur la plage
[0 V ; 6 V], le bruit de quantification moyen est de lordre de -103 dB. La relation dfinissant
la densit spectrale du bruit de quantification moyen est donne dans (5.9).
Sbq = .
1 Q2
(5.9)
N 12
Ce rsultat implique que seules les densits spectrales de puissance suprieure -100 dB
proviennent effectivement de lanalyse spectrale des fonctions de modulation.
- 142 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Evaluation RS232
de la PSD
outputs
PWM
EVM Board
Circuit logique
programmable
Systme
dacquisition
de donnes
outputs
Logic
PC
ISPGAL
TMS320LF2407
Modulation RCFM (Fonction Han)
Lalgorithme de gnration des fonctions de modulation dans le processeur de signal est dcrit
dans la figure 5.1.10. La seule diffrence par rapport lalgorithme de la RCFM dans le
systme de commande est que le signal de rfrence de la modulation est construit en interne
(systme sinusodal triphas quilibr) et ne dpend pas du systme dasservissement.
Linitialisation du processeur consiste :
Initialiser les registres dtat (ncessaires pour lexcution des instructions, stockage et
accs aux donnes),
Initialiser les mmoires des donnes et des programmes,
Initialiser les registres de gestion des vnements (EVA et EVB),
Initialiser le registre de contrle du Timer (base de temps pour la RCFM),
Dclarer et initialiser les interruptions (Interrupt Mask Register et Interrupt Flag
Register),
Activer les sorties numriques et PWM (sans temps mort),
Dclarer les fonctions associes au vecteur interruptions,
Initialiser la premire priode de la modulation RCFM,
Initialiser les variables globales et appel de la fonction "main()".
- 143 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Les trois signaux de rfrence de la modulation (systme triphas) sont calculs suivant un
pas alatoire (dbut de chaque priode de modulation). Afin dviter tout retard entre les
sorties du DSP, nous avons dvelopp une fonction de synchronisation des sorties numriques
avec les sorties PWM. Cette fonction est associe une interruption programme pour suivre
les variations alatoires de la priode de modulation.
Initialisation du DSP
TRCF :=
f min + X rand ( f max f min )
1
Le circuit logique ISPGAL transforme les six sorties du DSP en fonction des modulations Ha,
Hb et Hc. Les caractristiques de ces fonctions sont prsentes dans la section (II.1) du
troisime chapitre de ce mmoire. Les fonctions Han, Hbn et Hcn sont ensuite dduites par les
relations de (5.8). La figure 5.1.11 illustre la fonction Han rsultante de ce dispositif
exprimental.
- 144 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Norme Norme
Norme Norme
Norme Norme
- 145 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Figure 5.1.13 : Les PSD exprimentales dans la plage [2 kHz ; 50 kHz] pour diffrents
intervalles des frquences alatoires
Le spectre basse frquence de la RCFM pour les diffrents intervalles des frquences
alatoires est estim par la mthode Periodogram moyenne sur cinq acquisitions de longueur
106 points chacune et un chantillonnage 5 MHz. La premire remarque sur les spectres de
la figure 5.1.12 concerne la composante fondamentale (50 Hz) de la fonction Han qui nest pas
affecte par le choix de lintervalle des frquences alatoires. De plus, cette composante
fondamentale correspond exactement la valeur thorique.
Lutilisation des intervalles [6 kHz ; 10 kHz] [9 kHz ; 10 kHz] engendre des harmoniques 5
et 7 assez levs par rapport aux autres intervalles des frquences alatoires. Plus
- 146 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Figure 5.1.14 : Les PSD exprimentales de Han en haute frquence dans le cas des deux
intervalles [4 kHz ; 10 kHz] et [9 kHz ; 10 kHz]
- 147 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
reproche que lon peut formuler concernant ltude thorique, porte sur lestimation des
harmoniques 5 et 7 qui napparaissent pas dune faon aussi claire que dans les PSD
exprimentales.
Aprs cette confrontation de ltude thorique avec ltude exprimentale et la validation des
outils destimation de la PSD rsultante des mesures, nous allons vrifier la concordance entre
lanalyse des fonctions de modulation et lanalyse directe des courants de phase. Pour cela,
nous avons dvelopp des modles de simulation dynamique dans Matlab. La section suivante
est consacre la prsentation de ces modles de simulation, aux contraintes de
dveloppement, la viabilit de la nouvelle architecture de commande du convertisseur Cuk
et lanalyse des signaux rsultants des mesures.
II. Simulations
Toutes les simulations du convertisseur Cuk rversible sont effectues dans Matlab/Simulink,
cet outil prsente une grande flexibilit dans la simulation des systmes dynamiques. Grce
aux "S-functions" de Matlab, tout processus physique dont le modle est descriptible par
algorithme est implantable et simulable dans Simulink. Toutefois, lun des points qui restent
optimiser dans Simulink sont les modles des interrupteurs de puissance (Diode, IGBT,
MOSFET, Thyristor). Le paramtrage des modles de la bibliothque "Power System
ToolBox" ne correspond pas aux caractristiques des composants fournies par le constructeur
(figure 5.2.1).
- 148 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Le bon paramtrage des interrupteurs de puissance tant indispensable pour vrifier la validit
de la nouvelle architecture du Cuk, nous avons investigu deux voies :
La premire voie est de dvelopper des modles dIGBT et de diode en se basant sur
les caractristiques relles des composants (les paramtres constructeur). Cependant,
la complexit des algorithmes a engendr des temps de calcul trs considrables et des
dpassements mmoire (out of memory).
Les difficults de la premire dmarche nous ont pouss paramtrer ces composants
en utilisant les modles Matlab, mais en se rfrant loutil de simulation Saber qui
permet de paramtrer les cellules de commutation dune manire plus prcise. Nous
nous sommes bass sur des simulations effectues sur la structure Cuk existante avec
loutil Saber et valides exprimentalement dans [KIK-02]. La premire tape a
consist en la reproduction dun modle de simulation collant exactement avec le
modle de [KIK-02] (commande + composants passifs + frquence fondamentale
60 Hz). Le circuit de puissance est paramtr de telle sorte que le courant de phase, la
tension de sortie et la tension vcc soient identiques aux rsultats obtenus par "Saber".
- 149 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Dans ce mmoire nous nous limitons la description de la dmarche retenue pour valider le
modle du circuit de puissance. Nous avons vrifi le fonctionnement normal en rgime
permanent, et compar le taux de distorsion harmonique du courant de phase. Nous navons
pas pu comparer les spectres du courant de phase rsultant des deux outils de simulation, car
les travaux [KIK-02] nont pas port sur lanalyse spectrale approfondie des courants de
phase. Les caractristiques du modle de test sont les suivantes :
Rdc = 0,24
Puissance 5 kW Inductance de dcouplage Ldc = 2,2 mH
Frquence rseau 60 Hz (cosinusodal) Rsistance de Ldc
Frquence MLI f = 9 kHz Condensateur de sortie Cdc = 2350 F
Inductance dentre Lac = 2,5 mH Charge rsistive Rcharge = 18
Les rsultats de simulation en rgime permanent sont prsents dans la figure 5.2.2. Nous
avons obtenu des rsultats identiques avec Simulink, que ce soit du point de vue respect des
consignes courant et tension, ou du point de vue forme des signaux. Le THD du courant de
phase est le mme dans les deux cas. Une analyse annexe notre tude a consist en la
simulation du cas correspondant un changement instantan du mode redresseur MLI au
mode onduleur non autonome. Lobjectif de ce test est de vrifier leffet du paramtrage du
circuit de puissance sur les diffrentes variables dtat lors de la transition entre ces deux
modes de fonctionnement. Le mme paramtrage des IGBT et des diodes utilis dans le cas
de la figure 5.2.2, a permis dobtenir des rsultats de simulation identiques aux rsultats des
travaux de [KIK-02]. Le paramtrage du circuit de puissance utilis dans ces modles de test
est le mme exploit dans le modle de la nouvelle architecture.
- 150 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
0.795
- 151 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
(5)
(2)
(1)
(4)
(3)
- 152 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
g1
g4
g3
g6
g5
g2
g7
360
340
320 vdc
[V]
300
280
260
240
30
20
Redresseur
10
[A]
0 Icharge
-10 Onduleur non autonome
-20
-30
30
20 Id
10
Iq
[A]
0
-10
-20
-30
300
200 vc va vb
100
[V]
0
-100
-200
-300
30
20 ic ia ib
10
[A]
0
-10
-20
-30
0,47 0,48 0,49 0,50 0,51 0,52 0,53 0,54 0,55
[s]
Figure 5.2.5 : Rsultats de simulation du convertisseur Cuk en rgime permanent
- 153 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
La distribution des frquences alatoires est effectue suivant une loi de probabilit uniforme
et discrte afin que la reconstruction des courants de phase soit la plus proche possible des
signaux rels. La figure 5.2.6 compare le courant de phase ia reconstruit dans le cas de la
RCFM distribution uniforme continue avec celui reconstruit lors de lutilisation de la RCFM
distribution uniforme discrte. Dans cette figure, on distingue clairement que la qualit du
courant de phase rsultant de la reconstruction distribution discrte est meilleure par rapport
au courant rsultant de la reconstruction distribution continue. Ce qui renforce les
conclusion de ltude du chapitre 4.
Les rsultats de simulation de la figure 5.2.5 confirment la qualit de la reconstruction des
signaux, renforce lintrt de la boucle Feed-Forward du courant de charge, et valide le
dimensionnement des boucles dasservissement du courant de phase et de la tension de sortie.
Malgr laugmentation du temps de calcul par rapport au systme de commande existant,
larchitecture de commande dveloppe permet davoir le mme comportement dynamique,
que ce soit en courant ou en tension.
(a) (b)
Figure 5.2.6 : Courant de phase rsultant de la reconstruction en interne
(a) RCFM avec loi de distribution uniforme et continue
(b) RCFM avec loi de distribution uniforme et discrte.
- 154 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
De plus, la dformation des trois courants de phase (non sinusodaux) implique que tout le
dimensionnement labor sous lhypothse dun systme sinusodal triphas quilibr nest
plus approprie. Aussi, pendant la phase de dmarrage lasservissement des deux
composantes du courant de phase dans le rfrentiel de Park id et iq est insuffisant pour
asservir le courant de phase (la composante homopolaire io nest pas nulle). Ceci peut
expliquer en partie la perte de contrle du courant lors de la mise sous tension.
Ce problme de mise sous tension du convertisseur, na pas t trait par Kikuchi [Kik-02]
mme sil remet en cause le grand avantage des convertisseurs abaisseurs / lvateurs qui est
le contrle du courant ds la mise sous tension.
- 155 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
premier pic de la tension de sortie est de lordre de 340 V, reprsentant ainsi une diffrence de
100 V par rapport au cas prcdent.
De plus les trois courants de phase sont sinusodaux et contrls ds la phase de dmarrage.
Certes, cette procdure engendre une certaine lenteur dans la charge des condensateurs de
stockage dnergie, mais le rgime transitoire de la tension de sortie (ncessaire pour la
stabilit de la tension) reste quivalent entre les deux cas et est de lordre de 100 ms.
- 156 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Norme
Norme
Norme
Dans la figure 5.2.10, nous constatons que le spectre basse frquence de la tension Vaca est
image du spectre de la fonction Han. Ceci montre que les hypothses consistant ngliger les
ondulations de la tension Vcc et considrer une rfrence de modulation purement
sinusodale sont fondes et ne dtriorent pas lestimation de la PSD. Si lon revient sur la
relation (3.9) liant la fonction Han au courant de phase, les harmoniques du courant sont
dfinis par (Han_nf.Vcc) / (3.Lac.wnf). La dpendance du dnominateur, de cette relation, de la
frquence explique pourquoi les harmoniques basse frquence du courant sont plus amplifis
par rapport ceux de la fonction Han. Lexemple le plus parlant est celui des harmoniques 5 et
7 qui sont presque visibles dans le spectre de la tension Vaca (figure 5.2.10), et qui deviennent
- 157 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
plus importantes dans le spectre du courant. De plus, dans la relation (3.9) nous avons nglig
la tension Rac.ia afin de trouver une relation directe entre le courant de phase et la fonction
Han. Ceci explique lapparition des harmoniques 11, 13, 17, et 19 dans le cas de lintervalle
[9 kHz ; 10 kHz].
Figure 5.2.10 : La PSD de la tension Vaca en basse frquence dans le cas o la RCFM utilise
des frquences alatoires dans [9 kHz ; 10 kHz]
- 158 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
la PSD doit tre effectue sur un nombre de priodes fondamentales suprieur 20 (car les
calculs sur 20 priodes conduisent des rsultats similaires ceux rsultant de 10 priodes
fondamentales). Cependant, lalgorithme et les outils de calcul que nous avons dvelopp
pour valuer la PSD thorique posent des problmes de saturation mmoire au-del de 20
priodes avec un pas de 5 Hz. Ce travail peut sinscrire dans les perspectives futures de cette
thse.
- 159 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Sorties PWM
Conversion Module "Event Manager"
Entres/Sorties numriques
Analogique/Numrique
- 160 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
3,3
Module de gestion des entres / sorties numriques.
Autres modules : SPI (Serial Peripheral Interface), SCI (Serial Communications
Interface) et CAN (Controller Area Network), ayant le mme degr de priorit que le
module de conversion numrique analogique en terme dinterruptions. Ces trois
modules ne sont pas utiliss dans le cadre de ce travail.
Afin destimer ce temps, nous avons dvelopp et implment le code source du systme de
commande dans le DSP TMS320LF2407. La partie initialisation du DSP nest pas
comptabilise dans le temps dexcution de lalgorithme de commande (car excute une
seule fois au dmarrage). Le temps dexcution de lalgorithme de commande est la somme
des temps dexcution des fonctions suivantes :
Acquisitions des deux courants icoup et icharge.
Gnration de la modulation RCFM.
Construction des sinus et cosinus de rfrence de la transformation de Park.
Les trois correcteurs PI numriques (avec gestion du dmarrage).
Limitation des rfrences de la modulation et de lerreur corrige de lasservissement
de la tension de sortie.
Reconstruction des diffrents signaux (courant de phase et tension de sortie)
Adaptation des impulsions de commande (logique programmable).
- 161 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Comme lexplique la figure 5.3.2, lamplitude de la porteuse dans le DSP est proportionnelle
la priode de dcoupage. Etant donne que les priodes de dcoupage varient alatoirement,
lamplitude de la porteuse suivra cette variation. Afin de maintenir le mme rapport cyclique
par rapport la MLI ( frquence fixe), la rfrence de la modulation doit tre adapte en
consquence. Par exemple si le signal de rfrence de la modulation est damplitude 300, dans
le cas de la MLI, la priode est dfinie par un pas dincrmentation du compteur adquat
(multiple du cycle dhorloge) pour avoir une amplitude de la porteuse de 305 par exemple
(305 pas dincrmentation). Par contre, dans la RCFM cette amplitude variera pour chaque
priode gnre (et ne sera gale 305 que pour une priode correspondant celle de la
MLI), ce qui engendre une variation du rapport cyclique. Pour viter ce problme on doit agir
sur le pas dincrmentation du compteur pour maintenir la mme amplitude de la porteuse.
Cependant, le pas temporel du compteur est dfini dans la phase dinitialisation du DSP est ne
peut tre modifi au cours de lexcution du processeur. La seule solution que lon peut
envisager est de faire varier lamplitude de la rfrence de modulation avec lamplitude de la
porteuse de faon maintenir le mme rapport cyclique par rapport la MLI. Cette adaptation
ncessite deux cycles dhorloge.
- 162 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Pour estimer le temps ncessaire pour ladaptation des impulsions de commande, nous nous
sommes bass sur le retard engendr par le circuit exprimental (figure 5.3.3). Lors de la
conception de cette carte, nous avons optimis lutilisation des circuits programmables en
traduisant lalgorithme de la figure 2.2.2 en plusieurs fonctions logiques excutables en
parallle. Les neuf circuits ISPGAL utiliss, sont rpartis en quatre couches fonctionnelles.
Chaque couche engendre un retard de 10 ns (quivalent au retard dun composant). Ceci
implique un retard de 40 ns entre entre et sortie de ce circuit.
En plus du temps dexcution de lalgorithme de commande, les drivers des IGBT engendrent
aussi un retard entre les ordres du systme de commande et la commande effective des IGBT.
Pour valuer ce retard, nous avons tudi le driver SKH71 de la famille SEMIKRON. Ce
driver est destin commander sept IGBT ou MOSFET dune manire indpendante. Il
engendre un retard maximal lors de lamorage et du blocage de lordre de 0,6 s.
- 163 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Si lon tient compte de tous ces temps que nous avons valu ci-dessus, on peut estimer le
retard engendr par la commande de lordre de 5 s.
- 164 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
La dernire contrainte dont nous avons tenu compte, est le retard engendr par la commande,
nous avons introduit les retards valus dans la section prcdente au niveau des gchettes des
IGBT, des capteurs de mesure, du gnrateur de la RCFM, du reconstructeur des signaux, des
limiteurs et des correcteurs. Les rsultats de simulation du modle tenant compte de toutes ces
contraintes exprimentales, sont illustrs dans la figure 5.3.5.
Ces simulations sont effectues avec une modulation RCFM utilisant des frquences
alatoires dans lintervalle [6 kHz ; 10 kHz]. Dans cette figure, nous constatons que
lenveloppe sinusodale du courant de phase prsente une certaine saturation aux extremums
de chaque sinusode. Cette dformation est directement lie la tension dalimentation. Car
sur le spectre du courant, les harmoniques 5, 7, 9, 11 et 13 sont plus important par rapport aux
rsultats de la figure 5.2.9. De plus, la largeur de lenveloppe sinusodale aux deux
extremums, est suprieure 1 A. Ceci est d la dformation de la tension rseau et aux
variations alatoires de la frquence. En revanche, la composante fondamentale du courant de
phase et la tension de sortie ne sont pas affectes par la dformation de la tension
dalimentation.
- 165 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
- 166 / 188 -
Chapitre 5 : Etude par exprimentation & simulation de larchitecture de commande
Conclusion
Ce chapitre a t consacr ltude exprimentale de la RCFM et aux simulations portant sur
le convertisseur Cuk rversible. Pour valuer les PSD rsultantes des courants et fonctions de
modulation, nous avons dans un premier temps tudi et compar les mthodes numriques
destimation de la PSD, la famille tudie comprend les mthodes non paramtriques. Afin de
tester ladaptabilit de ces mthodes lanalyse de la RCFM, nous avons utilis comme
spectre de rfrence celui dun signal carr dont lamplitude varie alatoirement entre -1 et 1,
et dont la PSD thorique est connue. La PSD en basse frquence est estime par la mthode
Periodogram moyenne sur cinq acquisitions de longueur 106 points chantillonnes 5 MHz.
Alors quen haute frquence, cest la mthode de Welch avec un recouvrement de 50% et un
fentrage de pondration de Blackman qui est utilise.
Concernant lanalyse de la fonction Han, ltude exprimentale a valid la thorie dveloppe
dans le troisime chapitre. De plus, lanalyse de cette fonction prvoit le comportement du
courant de phase. Aussi, cette tude exprimentale a permis de dterminer linsuffisance de la
mthode thorique dans lestimation des harmoniques 5 et 7.
Nous avons pu valider par des simulations, le dimensionnement du convertisseur, sa
dynamique ainsi que la reconstruction des diffrents signaux. Nous avons aussi tenu compte
du retard de la commande, en valuant les temps dexcution des diffrents algorithmes, de
lacquisition et de pilotage des IGBT. Pour valuer ces temps, nous avons implant les
diffrents algorithmes dans le processeur TMS320LF2407, et avons ralis le circuit logique
gnrant les impulsions de commande.
La dformation de la tension dalimentation constate lors des mesures sur le rseau, nous a
amen en tenir compte dans nos modles. Mme si cette dformation ne respecte pas les
hypothses de la modlisation, le courant de phase suit la forme de la tension et le
convertisseur se comporte comme une charge rsistive.
- 167 / 188 -
Conclusion gnrale
Conclusion gnrale
Notre tude sinscrit dans la continuit du travail de J. KIKUCHI [KIK-02]. Nous avons tir
parti des dveloppements rcents des techniques de modulation frquence alatoire et des
techniques "sensorless" labores pour le pilotage des machines lectriques, afin dlargir leur
utilisation aux redresseurs prlvement sinusodal, particulirement le convertisseur Cuk
rversible. Lobjectif est daboutir une commande du Cuk rversible permettant damliorer
la fiabilit et la qualit du convertisseur tout en rduisant le cot de ralisation afin de rendre
notre convertisseur plus comptitif que le convertisseur Boost. Tout au long de ce mmoire,
nous avions pour but dvaluer la faisabilit dune telle structure.
- 168 / 188 -
Conclusion gnrale
montage de puissance suivant l'tat statique des interrupteurs. Lvolution en frquence de ces
modles nous a permis dvaluer les diffrents lments du circuit de puissance. Quant aux
systmes de rgulation, nous avons trait lasservissement du courant de phase dune faon
similaire au convertisseur Boost triphas. Tandis que lasservissement de la tension de sortie a
ncessit une radaptation du modle initial (non linaire et variable dans le temps) en le
transformant en un modle moyen petits signaux sur lequel on a pu appliquer des techniques
danalyse et de commande linaires. Afin damliorer le comportement dynamique de
lasservissement de la tension de sortie vis--vis des variations du courant de charge tout en
rduisant le dpassement par rapport la consigne, nous avons introduit une boucle
"Feedforward" du courant de charge.
Certes le dcoupage en MLI rduit les harmoniques basse frquence, mais au dtriment dune
augmentation des harmoniques au voisinage des multiples de la frquence de dcoupage. Afin
dattnuer ces harmoniques sans dtriorer la composante fondamentale, nous avons opt
pour la modulation RCFM. Modulation qui sintgre au systme de commande dune manire
similaire la MLI. De plus, comme nous lavons signal plus haut, lutilisation des
frquences infrieures la frquence nominale rduit les pertes par commutation dans les
IGBT. Dans le chapitre 3, nous avons valu thoriquement le gabarit de la densit spectrale
de puissance du courant de phase en dveloppant un algorithme de calcul numrique intgrant
les quatre contraintes suivantes : priode fondamentale de 20 ms, nombre entier de priodes
de dcoupage, la moyenne des frquences discrtes du tirage alatoire correspond la
moyenne dune distribution continue et le plus petit multiple commun de ces frquences est
situ au del de 2,5 MHz.
Cette tude nous a amen conclure que lintervalle de frquence [6 kHz ; 10 kHz] est le
meilleur compromis attnuation du spectre haute frquence et dtrioration du spectre basse
frquence. En revanche, la composante fondamentale nest pas dtriore par la RCFM et
reste indpendante de lintervalle des frquences alatoires choisi.
- 169 / 188 -
Conclusion gnrale
RCFM distribution discrte permet une reconstruction des signaux sans effets ngatifs
lorsque le tirage des frquences alatoires est judicieusement contrl.
Dans la dernire partie de ce mmoire, nous avons pu valider par des simulations le
dimensionnement du convertisseur, sa dynamique ainsi que la reconstruction des diffrents
signaux. Nous avons aussi valid ltude thorique de la RCFM par une tude exprimentale
bas niveau. Etude qui a montr certaines limites de la dmarche thorique dans lestimation
de la PSD en basse frquence et pourrait tre complte travers des tudes ultrieures.
Les simulations prsentes dans ce chapitre ont port aussi sur des tests concernant des
contraintes exprimentales telles que la mise sous tension du convertisseur, la quantification
des signaux mesurs, la dformation de la tension dalimentation et lintroduction du retard
rel des diffrents lments du convertisseur. Les rsultats de simulation ont confirm la
faisabilit exprimentale de notre architecture et ont justifi son intrt par rapport la
structure Boost rversible.
Les suites que nous pourrions donner cette tude sont multiples :
Notre tude thorique de la modulation RCFM dans le domaine temporel a point les limites
de cette dmarche. Dautant plus que la modlisation des fonctions de modulation dans le
temps peut varier dun convertisseur lautre. Nous pensons quune approche thorique de
cette modulation dans lespace des vecteurs de commande (Space Vector Modulation SVM)
serait plus gnrique.
- 170 / 188 -
Bibliographie
Bibliographie
Bibliographie -chapitre 1-
[ALL-02] J. M. Aller, A. Bueno and T. Paga, Power system Analysis using space-vector
transformation" IEEE Transactions on Power Systems, vol. 17, no. 4, pp. 957-
965, 2002.
[BUS-82] A. Busse and J. Holtz, Multiloop control of a unity power factor switching ac
to dc converter" IEEE PESC 82, pp. 171-179, 1982.
- 171 / 188 -
Bibliographie
[ISM-92] E. Ismail and R. W. Erickson, A single transistor three phase resonant switch
for high quality rectification" IEEE PESC 92, pp. 1341-1351, 1992.
[JAN-98] Y. Jang and M. M. Jovanovic, A novel robust harmonic injection method for
single-switch three-phase discontinuous-conduction-mode Boost rectifiers"
IEEE Transactions on Power Electronics, vol. 13., pp 824-834, 1998.
[KIK-00] J. Kikuchi, M.D. Manjrekar and T.A. Lipo, Complementary half controlled
Three-Phase PWM Boost rectifier for multi-dc-link applications" IEEE APEC
00, pp. 494-500, 2000.
[KIK-02] Jun Kikuchi and Thomas A. Lipo, Three-Phase PWM Boost-Buck Rectifiers
With Power-Regenerating Capability" IEEE Trans. on Industry Applications,
vol. 38, pp. 1361-1369, 2002.
[KIK-99] J. Kikuchi, M.D. Manjrekar and T.A. Lipo, Performance improvement of half
controlled three phase PWM Boost rectifier" IEEE PESC 99, pp. 319-324,
1999.
[KOL-93] J.W. Kolar, H. Ertl, and F. C. Zach, A comprehensive design approach for a
three-phase high-frequency single-switch discontinuous-mode boost power
factor corrector based on analytically derived normalized converter component
- 172 / 188 -
Bibliographie
[KOL-99] J.W. Kolar and H. Ertl, Status of the techniques of Three-Phase rectifier
systems with low effects on the mains" IEEE INTELEC 99, n 14-1, 1999.
[LEE-00] D. C. Lee, G-M. Lee, and K. D. Lee, DC-Bus voltage control of three-phase
AC/DC PWM converters using feedback linearization" IEEE Trans. on Industry
Applications, vol. 36, no. 3, pp. 826-833, 2000.
[LEE-02] D. C. Lee and D. S. Lim, AC voltage and current sensorless control of three-
phase PWM rectifiers" IEEE Transactions on Power Electronics, vol. 17, pp
883-890, 2002.
[LOT-95] C. Lott, Filtrage actif parallle des harmoniques de courant sur des rseaux
industriels : tude thorique et ralisation dune maquette GTO" thse de
lInstitut National Polytechnique de Lorraine, 17 octobre 1995.
- 173 / 188 -
Bibliographie
rectifiers" IEEE Trans. on Industry Applications, vol. 37, no. 4, pp. 1019-1027,
2001.
[MAL-87] L. Malesani and P. Tenti, Three-phase AC/DC PWM converter with sinusoidal
AC currents and minimum filter requirements," IEEE Trans. on Industry
Applications, vol. IA-23, no. 1, pp. 71-77, 1987.
[MAR-98] D.C. Martins, A.H. de Oliveira and I. Barbi, Three-Phase rectifier using Sepic
dc-dc converter in continues conduction mode for power factor correction"
IEEE INTELEC 98, n 20-4, pp 491-497, 1998.
[PAN-99] C.-T. Pan and J.-J. Shieh, A single-stage three-phase boost-buck AC/DC
converter based on generalized zero-space vectors" IEEE Transactions on
Power Electronics, vol. 14, no. 5, pp. 949-958, 1999.
- 174 / 188 -
Bibliographie
quality problems" ELSEVIER Electric Power Systems Research, vol. 141, no. 1,
pp. 54-64, 2002.
[SHI-00] J.-J. Shieh, SEPIC derived three-phase switching mode rectifier with
sinusoidal input current" IEE Proceedings Part B. -Electrical Power
Applications, vol. 147, no. 4, pp. 286-294, 2000.
[SHI-02] J.-J. Shieh, A reversible Three-Phase step up/down switching mode rectifier"
Journal of Marine Science and Technology, vol. 10, no. 2, pp. 118-127, 2002.
- 175 / 188 -
Bibliographie
[XU-94] J. Xu, Filtrage actif parallle des harmoniques des rseaux de distribution
dlectricit" thse de lInstitut National Polytechnique de Lorraine, 20 janvier
1994.
[ZIO-85] P. D. Ziogas, Y.-G. Kang, and V. R. Stefanovic, PWM control techniques for
rectifier filter minimization" IEEE Trans. on Industry Applications, vol. IA-21,
no. 5, pp. 1206-1214, 1985.
- 176 / 188 -
Bibliographie
Bibliographie -chapitre 2-
[SIN-00] B. N.Singh, P. Jain and G. Joos, Three-phase ac/dc regulated power supplies :
A comparative evaluation of different topologies" Proceeding of the IEEE
APEC Conference, 2000, pp. 513-518.
[LEE-00] D.-C. Lee, Advanced nonlinear control of three-phase PWM rectifiers" IEEE
Proc. Elect. Power Appl., vol. 147, no. 5, pp 361-366, 2000.
Bibliographie -chapitre 3-
[BEC-00] M. M. Bech, Analysis of random pulse-width modulation techniques for power
electronic converters" Ph.D. Thesis, Aalborg University (denmark), August
2000.
- 177 / 188 -
Bibliographie
- 178 / 188 -
Bibliographie
Bibliographie -chapitre 4-
[LEE-02] D-C Lee and D-S Lim, AC voltage and current sensorless control of
threephase PWM rectifiers, IEEE Transactions on Power Electronics, Vol 17,
pp. 883-890, 2002.
[BOH-97] S. Bhowmik, A. Van Zyl, R. SPEE and J. H. Enslin, Sensorless current control
for active rectifiers" IEEE Trans. on Industry Applications, vol. 33, no. 3, pp.
765-773, 1997.
- 179 / 188 -
Bibliographie
Bibliographie -chapitre 5-
[FON-00] R. J. Fontana, A note on power spectral density calculations for jittered pulse
trains" Multispectral Solutions Inc., 2000.
- 180 / 188 -
ANNEXE -A-
ANNEXE -A-
Une grande partie de notre tude reposant sur des grandeurs lectriques sinusodales et non
sinusodales, il est ncessaire de faire un rappel sur les grandeurs lectriques en prsence
dharmoniques.
Toute fonction priodique peut tre dcompose sous la forme dune srie de Fourier dfinie
par :
X(t) = X 0 + X n 2 sin[n(t + n )]
(A.1)
n =1
Avec
X 0 : valeur de la composante continue, gnralement nulle
:pulsation de la frquence fondamentale
n : rang de l ' harmonique
X n :valeur efficace de l ' harmonique de rang n
n.n :phase de l 'harmonique de rang n
Les expressions instantanes du courant i(t) et de la tension v(t) peuvent donc se mettre sous
la forme :
T 0
X eff =
1 T 2
X (t) dt (A.3)
I V
Ieff = 2
et Veff = 2
(A.4)
n =1 n =1
n n
Facteur de puissance
Dans le cas dun rseau quilibr, la puissance apparente S est dfinie par :
- 181 / 188 -
ANNEXE -A-
P = Vn I n cos(n n ) , o n = n n
(A.6)
n =1
V
I n cos(n n )
F= = n =1
n
P
V . I
(A.7)
S
2 2
n =1 n =1
n n
Dans le cas o lon considre une tension rseau parfaitement sinusodale, lexpression du
facteur de puissance de la relation (A.7) se simplifie et devient :
I1cos 1
F=
I
(A.8)
2
n =1
n
Q = V1 I1 sin 1 (A.9)
S = P 2 + Q2 + D2 (A.10)
F=
P
P + Q2 + D2
(A.11)
2
Cette dernire expression fait clairement apparatre le fait que le facteur de puissance est
dgrad par la prsence dharmoniques dune part et par la consommation dnergie ractive
dautre part.
Le taux global de distorsion harmonique (THD) permet dvaluer lcart entre la forme
donde relle et la forme donde sinusodale pour un courant ou une tension. Il reprsente le
rapport de la valeur efficace des harmoniques celle du fondamental :
- 182 / 188 -
ANNEXE -A-
X
2
THD% = n =2
n
.100 (A.12)
X1
D = V1 I1 THD (A.13)
Et le facteur de puissance :
cos 1
F=
1 + THD 2
(A.14)
- 183 / 188 -
ANNEXE -B-
ANNEXE -B-
Pour valuer les PSD partielles de la relation (3.51), lvaluation de lesprance mathmatique
(de la somme infinie des transformes de Fourrier) simpose. Les PSD partielles peuvent tre
dcomposes en trois catgories. Chaque catgorie dpend de la comparaison entre m et m% .
Les trois catgories sont :
Les termes de la diagonale pour ( m = m% )
Les termes de la matrice triangulaire suprieure pour ( m < m% )
Les termes de la matrice triangulaire infrieure pour ( m > m% )
Lorsque ( m = m% )
Tout dabord nous commenons par valuer ( t ,m% t 0,m ). Daprs (B.1) nous remarquons que
....
(T + ... + T + T + ... + T
2,m 2,K 1,1 1,M + T0,1 + ... + T0,m 1 ), = 2
(T1,m + ... + T1,K + T0,1 + ... + T0,m 1 ), = 1
t ,m% t 0,m = 0, =0
T + ... + T + T + ... + T ,
(B.1)
=1
0,m 1,m 1
....
{
+ E U(f ;0,m )U* (f ;1,m% )e jw (T0,m +...+ T1,m1 ) e jw (1,m% T1,m% 0,mT0,m ) }
+ E {U(f ; 0,m )U* (f ; 2,m% )e e } +
jw (T0,m +...+ T2,m1 ) jw ( 2,m% T2,m% 0,m T0,m )
- 184 / 188 -
ANNEXE -B-
{
E {} = U(f ;0,m )e jw0,m T0,m }
{ }
2K-1 termes
{ }
K-1 termes (B.4)
+ E U(f ;0,m )e jw (1 0,m )T0,m
E {U (f ; )e } E {e } E {e jwT }
K-1 termes
+ E {U (f ; )e } E {e
* jw1,m% T1,m jwT0,m +1 1,m 1
} E {e jwT } +
%
1,m
2K-1 termes
De plus, comme on utilise une loi de probabilit uniforme et que toutes les priodes Ti sont
connues pralablement, alors on peut crire :
{ }
En utilisant (B.5), on peut simplifier (B.4) qui devient :
E {} = E U(f ; m )
2
+ E {U(f ; m ) e } E {U (f ; } E {e }
+
jw (1 m )T * jw m% T K 1
jwT
)e
=1
%
m
E {e jwT }
E {e }
+ K 1
jwT K 1 =
1 E {e jwT }
(B.7)
=1
K
Dans (B.6) les lignes 2 et 3 reprsentent chacune des complexes conjugus. Ainsi en utilisant
cette proprit et la relation (B.7), la PSD partielle scrit lorsque m = m
% :
Sm,m% =
1
{
E U(f ; m )
2
}
E {e jwT }
KT
- 185 / 188 -
ANNEXE -B-
Lorsque ( m < m% )
....
(T + ... + T + T + ... + T + T + ... + T = 2
2,m% 2,K 1,1 1,K 0,m 1 ),
t ,m% t 0,m = T0,m + ... + T0,m% 1 , =0
T + ... + T + T + ... + T ,
(B.9)
=1
0,m % 1
....
K + (m
% m) lorsque 0. Le terme ( ,m% 0,m ) sexprime de la mme faon que dans
E {} =
{
+ E U(f ;0,m )U* (f ; 2,m% )e jw (T2,m% +...+ T0,m1 ) e jw (2,m% T2,m% 0,m T0,m ) }
+ E {U(f ; jw (T1,m% +...+ T0,m1 ) jw ( 1,m% T1,m% 0,m T0,m )
}
+ E {U(f ; }
0,m )U* (f ; 1,m% )e e
jw (T0,m +...+ T0,m% 1 ) jw (0 ,m% T0,m% 0,m T0,m ) (B.10)
+ E {U(f ; }
0,m )U* (f ;0,m% )e e
jw (T0,m +...+ T1,m% 1 ) jw (1,m% T1,m% 0,m T0,m )
+ E {U(f ; } +
0,m )U* (f ;1,m% )e e
jw (T0,m +...+ T2,m% 1 ) jw ( 2,m% T2,m% 0,m T0,m )
0,m )U* (f ; 2,m% )e e
{ }
Le dveloppement de (B.10) conduit (B.11)
E {} = U(f ;0,m )e jw0,m T0,m
{ }
% - m) -1 termes
2K - (m
jw (1 0,m )T0,m
E {U (f ; )e } E {e } E {e }
0,m % - m) -1 termes
(m % - m) -1 termes
K - (m (B.11)
+ E {U (f ; )e } E {e } E {e }
* jw0,m% T0,m% jwT0,m +1 jwT0,m% 1
%
0,m % - m) -1 termes
K + (m
+ E {U (f ; )e } E {e } E {e } + .
* jw1,m% T1,m% jwT0,m +1 jwT1,m% 1
%
1,m % - m) -1 termes
2K + (m
- 186 / 188 -
ANNEXE -B-
=1
=0
=
K (m
% m) 1
1 E {e jwT }
(B.13)
=1
K
E { jwT }
E {e } = E {e }
+ K + (m
% m) 1
jwT K + (m m) 1 jwT (m m) 1 + e
1 E {e jwT }
% %
(B.14)
=0
K
jw (1 m% )T E {e }
E {U(f ;m )e m } E {U (f ;m% )e }
1 jwT K (m
% m) 1
= jw T
KT 1 E {e jwT }
*
Sm,m% K
E {e jwT }
+ E {U(f ;m )e } E {U (f ;m% )e } E {e }
(B.15)
K + (m
% m) 1
+
jw (1 m )T % m) 1
1 E {e jwT }
* jw m% T (m
jwT
K
Lorsque ( m > m% )
Dans ce cas :
=0
=1
En crivant :
E {e jwT }
E {e jwT } = E {e jwT }
+ K (m
% m) 1
+
K (m
% m) 1 (m m)
% 1
1 E {e jwT }
(B.18)
=0
K
E { jwT }
E {e }
+ K + (m
% m) 1
jwT K + (m m) 1 = e
1 E {e jwT }
%
(B.19)
=1
K
Ainsi, nous pouvons dterminer la densit spectrale dnergie partielle lorsque m > m
% :
- 187 / 188 -
ANNEXE -B-
E{e jwT }
{ } { }
{ }
K (m
% m) 1
=
(m m)
+
% 1
KT
jwm T jw(1 m% )T
1 E{e jwT }
1
*
Sm,m% E U(f; m )e E U (f; %
m )e E e jwT
K
jwm% T E {e }
+E{U(f;m )e } E{U (f;m% )e }
(B.20)
K + (m
% m) 1
jw(1 m )T
1 E{e jwT }
jwT
*
K
- 188 / 188 -
ANNEXE -B-
Rsum :
Ce mmoire traite de la commande alatoire du convertisseur Buck-Boost triphas
prlvement sinusodal du courant. Cette famille de convertisseurs est peu matrise, et ce
malgr le grand avantage quelle apportent par rapport la structure Boost. Lobjectif de ce
travail est de dvelopper une commande adapte cette structure amliorant la fiabilit et la
qualit du convertisseur tout en rduisant le cot de ralisation. L'objectif tant daboutir un
convertisseur plus comptitif par rapport au convertisseur Boost.
Aprs une prsentation du contexte de ltude dans le premier chapitre, le second chapitre est
consacr au dimensionnement et la mise en uvre de la commande du convertisseur.
Le troisime chapitre dcrit la modulation frquence alatoire RCFM, permettant dattnuer
la densit spectrale des courants de ligne tout en conservant lamplitude de la composante
fondamentale.
Le quatrime chapitre porte sur lamlioration de la fiabilit et la rduction du cot du
convertisseur. Qui consiste en llaboration dune architecture de commande saffranchissant
des capteurs de courant alternatif, et rduisant le nombre total des capteurs utiliss ce qui
amliore sa fiabilit.
Enfin, le dernier chapitre est consacr ltude exprimentale et aux simulations. Ltude
exprimentale de la modulation a permis de valider les rsultats thoriques de lanalyse de la
RCFM. De plus, les simulations (intgrant des contraintes exprimentales) ont permis de
valider le dimensionnement du convertisseur, sa dynamique ainsi que la reconstruction des
diffrents signaux
Abstract :
This thesis deals with the random control of three phase Buck-Boost converter with input
current power factor correction. These structures are less investigated in spite of the great
advantage which they bring compared to the Boost structures. So in this PhD we investigate
new solutions allowing to obtain a Buck-Boost structure more competitive compared to the
Boost (reliability, cost and EMC constraints).
After a short presentation of the study context in the first chapter, the second one deals with
the control system, dimensioning and implementation.
The third chapter describes the randomized carrier frequency modulation. This technique
attenuates the power spectral density of input current and preserves the fundamental
component.
The fourth chapter is dedicated to reliability improvement and reduces of manufacture cost.
This consists to reduce the number of used sensors and only measure DC bus signals. This
sensorless control requires fast algorithms to rebuild all signals (not measured but necessary
for control).
The final chapter is dedicated to the experimental study and simulations. Concerning the
modulation analysis, the experimental study validates theoretical results. Moreover,
simulations (with experimental constraints) allow validating the dynamic of converter, its
dimensioning, as well as the rebuilding of the all signals
Key words : Power Factor Corrector (PFC), Buck-Boost, Random Carrier Frequency
Modulation (RCFM), Power Spectral Density (PSD), Harmonics, Modulation function, Power
electronics, Real time, Digital Signal Processor (DSP).
- 189 / 188 -