PDF

UNIVERSITÉ DU QUÉBEC À TROIS-RIVIÈRES
MÉMOIRE PRÉSENTÉ À
L'UNIVERSITÉ DU QUÉBEC À TROIS-RIVIÈRES
COMME EXIGENCE PARTIELLE POUR L'OBTENTION D'UN

DIPLÔME DE MAÎTRE ES SCIENCES APPLIQUÉES EN
ÉLECTRONIQUE INDUSTRIELLE
PAR
SYLVAIN LAHAIE
ÉTUDE ET CONCEPTION D'UN BANC D'ESSAI:

APPLICATION À L'ALIMENTATION DES
ENTRAÎNEMENTS À VITESSE VARIABLE
NOVEMBRE 1994
Université du Québec à Trois-Rivières
Service de la bibliothèque
Avertissement
L’auteur de ce mémoire ou de cette thèse a autorisé l’Université du Québec

à Trois-Rivières à diffuser, à des fins non lucratives, une copie de son
mémoire ou de sa thèse.
Cette diffusion n’entraîne pas une renonciation de la part de l’auteur à ses

droits de propriété intellectuelle, incluant le droit d’auteur, sur ce mémoire
ou cette thèse. Notamment, la reproduction ou la publication de la totalité
ou d’une partie importante de ce mémoire ou de cette thèse requiert son
autorisation.
" "
RESUME
Ce mémoire de Maîtrise présente la conception d'un banc d'essai pour

les Entraînements à Vitesse Variable (EVV) réalisé pour le Laboratoire
des Technologies Électrochimiques et des Électrotechnologies (LTEE)
d'Hydro-Québec. Il s'agit là d'un projet industriel de grande envergure
échelonné sur deux années dont l'étude préliminaire a été menée à
bien grâce à une collaboration étroite entre le LTEE et le Groupe de
Recherche en Électronique Industrielle de l'Université du Québec à
Trois-Rivières.
La conception prend comme point de départ la description des essaIS

requis dans le document SC-610 de "l'Electric Power Research
Institute" (EPRI). Elle repose également sur l'analyse des informations
contenues dans la littérature scientifique et sur des travaux de
simulation.
Ce travail couvre la partie du projet se rapportant aux simulations du

banc d'essai appliquées aux tests décrits dans le document SC-610. Le
simulateur est bâti dans l'environnement du logiciel MA TL AB ® et de son
interface graphique SIMULINK®. Les éléments constituant le banc
® SIMULINK et MATLAB sont les logiciels (droits réservés) de "The MATH WORKS Inc."
1
d'essai pour l'étude des perturbations sur les EVV ont été modélisés
avec ce simulateur.
La partie critique de la conception du banc d'essai se retrouve dans la

section décrivant les tests d'immunité des EVV aux perturbations sur
la forme d'onde de la source de tension. Pour cela, nous proposons un
montage qUI a pour objet de générer une forme d'onde de tension
perturbée, à la fréquence du réseau, tout en ne nécessitant qu'une
source électronique ayant une pUIssance réduite. Ce montage (appelé
"générateur d'harmoniques de tension" (GRT» et son circuit de
commande constituent la principale innovation technique des travaux
décrits de ce projet.
En plus de la conception du GRT, ce mémoire contribue au

développement des modèles de simulation numérique qui seront, dans
un premier temps, utiles à la conception du banc d'essai et par la suite,
utilisés pour les tests de compatibilité des systèmes.
Ce travail présente également l'étude d'un filtre actif dont la

configuration a inspiré la conception du GRT ainsi qu'un programme
MA TLAB sur la conception de filtres passifs accordés.
II
REMERCIEMENTS
Mes remerciements vont vers mon directeur de recherche, le

Professeur Venkatachari Rajagopalan. Son suivi, ses conseils éclairés et
ses encouragements ont été pour moi des facteurs déterminants dans
ma motivation de compléter ces études de deuxième cycle. Je lui suis
infiniment reconnaissant pour sa confiance et ses encouragements
soutenus. Au cours des années où j'ai eu le prévilège de le côtoyer j'ai
découvert en lui un chercheur émérite d'une clairvoyance
exceptionnelle.
Je remercie le chercheur post-doctoral Yao Ziwen, pour ses critiques et

ses conseils prodigués au cours de ces derniers mois de travail. Sa
gentillesse, sa disponibilité et sa compétence remarquable ont rendu le
travail avec lui très agréable et enrichissant. Qu'il soit assuré de mon
estime sincère.
Je remerCIe aussi les membres du corps professoral de la maîtrise en

Électronique Industrielle. Qu'ils soient assurés de mon respect et de
mon appréciation: Andrej Barwicz, Ahmed Cheriti, Jean-Luc Dion, Yves
Dubé, André Jacob et Pierre Sicard.
III
Un grand merCI à Madame Nicole Montour et Madame Marlène
Tremblay, secrétaires au LTEE, pour leur gentillesse et pour n'avoir
JamaIs refusé de me venir en aide tout au long de mes études de
second cycle.
Je remerCIe la direction d'Hydro-Québec; spécialement M. Antoine

Duchesne, chef du service d'Électrotechnologie au LTEE, de m'avoir
permis et encouragé à faire mes études de maîtrise en électronique
industrielle.
Finalement Je désire remerCIer tout particulièrement M. André

Charette, chercheur au LTEE, pour m'avoir constamment soutenu et
encouragé durant ces années d'étude. Nos discussions ont permIs,
grâce à ses vastes connaissances en électronique industrielle et son
esprit critique remarquable , de surmonter biens des problèmes
rencontrés au cours de mon travail.
IV
TABLE DES MATIÈRES
RÉSUMÉ .......................................................................................................................... 1
REMERCIEMENTS .................................................................................................... 111
TABLE DES MATIÈRES ........................................................................................... V
LISTE DES FIGURES ................................................................................................. XI
LISTE DES TABLEAUX ...................................................................................... XVII
LISTE DES SYMBOLES ET DES ABRÉVIATIONS ..................................... XVIII
CHAPITRE 1:
INTRODUCTION GÉNÉRALE .................................................................................... 1
1.1 Cadre et motivation du projet.. .............................................................. 1
1.2 Simulation numérique du banc d'essai .............................................. 2
1 .3 Conception du banc d'essai.. .................................................................... 2
. t'Ion des sUJe

1 . 4 Descnp . ts traI't'es ................................................................ .. 3
v
CHAPITRE 2:
ÉTUDE DES PERTURBATIONS ÉLECTRIQUES PRODUITES PAR LES
EVV ET CONCEPTION DES FILTRES PASSIFS ACCORDÉS ............................ 5
2.1 Introduction ................................................................................................... 5
2.2 Notions de base sur les harmoniques et sur le calcul des

filtres assisté par ordinateur .................................................................. 6
2.2.1 Perturbations de la forme d'onde de tension ......................... 6
2.2.2 Simulation des configurations les plus communes

d'EVV ..................................................................................................... 12
2.3 Réduction des harmoniques injectés sur le réseau .................... 1 7
2.3.1 Généralités ............................................................................................. 1 7
2.3.2 Conception des filtres passifs accordés ................................... l 9
2.3.2.1 Calcul des filtres dans le domaine fréquenciel.. ........... l 9
2.3.2.2 Description du programme de calcul de filtre .............. 20
2.3.2.3 Critères de conception des filtres ....................................... 2 1
2.3.2.4 Étapes de design des filtres .................................................. 22
2.3.2.5 Organigramme du programme de calcul des filtres ... 36
2.4 Exemple de calcul de filtres accordés .............................................. 3 9
2.4.1 Critères de conception spécifiés ............................................... 39
2.4.2 Résultats du calcul de filtre ......................................................... 39
2.4.3 Analyse des performances des filtres dans le

domaine fréquenciel ....................................................................... 40
VI
2.4.4 Analyse des performances des filtres dans le
domaine temporeL ........................................................................... 4 2
2.5 Conclusion .................................................................................................... 43
CHAPITRE 3:
ÉTUDE D'UN FILTRE ACTIF .................................................................................. 4 4
3.1 Introduction ................................................................................................ 44
3.2 Notions de base sur les filtres actifs ................................................ .4 5
3.2.1 Configuration classique de filtres actifs ................................ .4 5
3.2.2 Combinaison de filtres actifs et passifs ................................. .4 6
3.3 Description d'un filtre actif relié en série avec des filtres
passifs ............................................................................................................ 48
3.3.1 Portée de l'étude ................................................................................. 4 8
3.3.2 Description de la configuration du circuit.. .............................. 49
3.3.3 Description du circuit de contrôle ............................................... 50
3.3.4 Analyse des performances du filtre actif ................................ 53
3.3.4.1 Problème de résonnance parallèle .................................... 53
3.3.4.2 Problèmes de résonnance série .......................................... 56
3.3.5 Simulateur du filtre actif de Fujita et Akagi.. ...................... 59
3.3.5.2 Exemple de simulation d'un problème de résonnance

parallèle ............................................................................................. 5 9
VII
série ..................................................................................................... 61
3.4 Conclusion ...................................................................................................... 63
CHAPITRE 4:
CONCEPTION D'UN GÉNÉRATEUR D'HARMONIQUES DE TENSION ......... 64
4.1 Introduction .................................................................................................. 64
4.2 Étude des perturbations de la tension de ligne ............................. 65
. 4.2.1 Fluctuation de la fréquence de 60 Rz ...................................... 6 6
4.2.2 Variation de l'amplitude: .............................................................. 66
4.2.3 Variation de la forme d'onde de tension: .............................. 6 9
4.2.4 Dissymétrie du système triphasé (déséquilibre

des phases) ......................................................................................... 7 1
4.3 Description du générateur d'harmoniques de tension .............. 7 1
4.3.1 Tests d'immunité et utilité du GRT ............................................. 71
4.3.2 Principe de fonctionnement du GRT .......................................... 73
4.3.2.1 Transformation des harmoniques de courant en

harmoniques de tension ............................................................. 73
4.3.2.2 Réglage de la tension harmonique de ligne ................... 7 7
4.3.2.3 Utilité des filtres passifs accordés ..................................... 78
4.4 Conception et analyse du GRT dans le domaine

fréquenciel. .................................................................................................. 7 8
VIII
4.4.1 Conception des filtres accordés utilisé par le GHT ............... 7 8
4.4.1.1 Critères de conception des filtres du GHT: ..................... 7 9
4.4.1.2 Résultats du calcul des filtres du GHT .............................. 7 9
4.4.2 Analyse des performances du GHT dans le

domaine fréquenciel ....................................................................... 80
4.5 Simulation numérique du GHT dans le domaine temporeL ..... 8 5
4.5.1 Simulation du GHT appliqué aux tests d'immunité ............. 85
4.5.2 Simulation des tests d'immunité aux harmoniques

de tension ............................................................................................ 85
4.5.3 Simulation des tests d'immunité aux encoches

dans la forme de l'onde de tension ......................................... 96
4.5.4 Simulation des tests d'immunité aux perturbations

causées par la commutation de condensateurs ............. 1 00
4.6 Mode de réglage des harmoniques de tension ............................ 101
4.7 Conclusion ................................................................................................... 105
CHAPITRE 5:
CONCLUSION ........................................................................................................... 1 07
RÉFÉREN"CES ............................................................................................................ 1 09
ANNEXE A:
MODÈLES DE SIMULATIONS DES EVV ........................................................ 112
IX
ANNEXEB:
FICHIERS DU PROGRAMME DE CALCUL ET DE CONCEPTION DES

FILTRES PASSIFS ACCORDÉS ........................................................................... 116
ANNEXEC:
SIMULATEUR DU FILTRE ACTIF DE FUJITA ET AKAGI.. ...................... 1 28
ANNEXED:
SIMULATEUR DU GHT ........................................................................................ 1 29
ANNEXEE:
DÉTAIL DES SIMULATIONS DU GHT ............................................................. 1 34
x
LISTE DES FIGURES
Figure 2.1 Convertisseur CC/CA du point de vue des

harmoniques de courant .............................................................. 7
Figure 2.2a Courant de ligne théorique d'un redresseur triphasé

de 6 impulsions à charge inductive ....................................... 8
Figure 2.2b Valeurs théoriques et typiques des harmoniques de

courant produits par différents montages
redresseurs ..................................................................................... 11
Figure 2.3 Schéma fonctionnel d'un EVV ayant uniquement un

filtre capacitif comme lien en tension continue .............. 1 3
Figure 2.4 Courant de ligne du circuit de la figure 2.3 ..................... 1 3
Figure 2.5 Schéma fonctionnel d'un système d'alimentation

d'un EVV dans lequel on utilise un filtre Lf-Cf.. ............. 1 5
Figure 2.6 Schéma fonctionnel d'un système d'alimentation

d'un EVV avec un filtre inductif du côté du réseau ...... 1 5
Figure 2.7 Courant de ligne du circuit présenté à la figure 2.5 .... 1 6
Figure 2.8 Courant de ligne du circuit présenté à la figure 2.6 .... 1 7
Figure 2.9 Exemple de filtres passifs accordés ...................................... 1 8
Figure 2.10 Représentation symbolique monophasée d'un

convertisseur avec des filtres passifs accordés .............. 2 2
XI
Figure 2.11 Impédance caractéristique du réseau modèle a) [9] .... 2 4
Figure 2.12 Impédance caractéristique du réseau modèle b) [9] .... 24
Figure 2.13 a Organigramme du programme de calcul des

filtres .............................................................................................. 37
Figure 2.13 b Organigramme du programme de calcul des filtres

(suite) ............................................................................................. 38
Figure 2.14 Impédance individuelle des filtres ..................................... .4 0
Figure 2.15 Impédance du système dans le plan complexe ............. .4 1
Figure 2.16 Formes d'onde de tension et courant de ligne ................ .4 1
Figure 3.1 Topologie d'un filtre actif relié directement au

réseau ................................................................................................ 46
Figure 3.2 Combinaison d'un filtre actif relié en série avec la

ligne et de filtres passifs accordés reliés en
parallèle ........................................................................................... 47
Figure 3.3 Combinaison d'un filtre actif relié en série avec des
filtres passifs, l'ensemble étant relié en parallèle
avec la ligne .................................................................................... 4 8
Figure 3.4 Filtres actifs et passifs accordés reliés en série l'un

avec l'autre ..................................................................................... 50
Figure 3.5 Circuit monophasé équivalent de la figure 3.4 ............... 54
Figure 3.6 Circuit équivalent du circuit monophasé pour ILh ....... 55
XII
Figure 3.7 Caractéristiques du filtre pour les courants
harmoniques .................................................................................. 55
Figure 3.8 Circuit équivalent du filtre actif pour Vsh ........................ 57
Figure 3.9 Caractéristiques de filtration en fonction de Vsh .......... 58
Figure 3.10 Simulation d'un problème de résonnace parallèle ........ 6 0
;
Figure 3.11 Simulation d'un problème de resonnance série .............. 62
Figure 4.1 Exemple de test d'immunité aux creux de tension

(Test 16) ........................................................................................... 67
Figure 4.2 Exemple de test d'immunité aux crêtes de tension

(Test 17) ........................................................................................... 68
Figure 4.3 Schéma de base pour les tests d'immunité 16 et I7 ..... 68
Figure 4.4 Exemple de pollution de la tension de ligne par des

harmoniques (Test 13, 14, 15) ................................................. 6 9
Figure 4.5 Exemple de test d'immunité aux encoches dans la

forme d'onde de tension de ligne (Test 19) ...................... 70
Figure 4.6 Exemple de test d'immunité à la commutation des

condensateurs (Test 18) ............................................................ 70
Figure 4.7 Représentation d'un problème de résonnance

parallèle ........................................................................................... 74
Figure 4.8 Module de l'impédance Zsh vu du convertisseur ........... 7 4
XIII
Figure 4.9 Forme d'onde de tension de ligne résultant du
produit Zs * Ih tel que présenté à l'équation 4.1. ........ 7 4
Figure 4.10 Schéma fonctionnel du GHT ..................................................... 7 5
Figure 4.11 Circuit équivalent du GHT pour un harmonique ............ 76
Figure 4.12 Module de l'impédance des composantes du GHT ......... 8 0
Figure 4.13 Circuit équivalent monophasé du GHT ................................ 82
Figure 4.14 Fonction de transfert de la source électronique et du

transformateur de distribution en fonction de f ............ 84
Figure 4.15 Génération d'harmoniques: combinaison #1 .................... 88
Figure 4.21 Circuit de réglage d'un harmonique de tension .............. 91
Figure 4.22 Réponse dynamique du circuit de réglage pour

l'exemple de la figure 4.15 (h5=33 volts et 0°) .............. 91
Figure 4.23 Formes d'onde pour Ve = 0 volt pour une charge de

600 kV A ......................................................................................... 93
XIV
Figure 4.24 Forme d'onde de courant et tension du système
pour l'essai #6 de génération d'harmoniques de
tension .............................................................................................. 9 4
Figure 4.25 Génération d'encoche de tension: essaI #1 ........................ 97
Figure 4.31 Forme d'onde de courant et tension du système

pour l'essai #6 de génération d'une encoche de
tension .............................................................................................. 9 9
Figure 4.32 Exemple de forme d'onde de tension produite pour

les tests d'immunité aux commutations de
condensateurs ............................................................................. 100
Figure 4.33 Schéma fonctionnel du contrôleur utilisé pour le

réglage des harmoniques de tension ................................ 104
Figure A.l Schéma SIMULINK d'un système d'alimentation

d'un EVV avec onduleur "6 steps" dans lequel on
utilise un filtre Lf-Cf .............................................................. 1 1 3
Figure A.2 Schéma SIMULINK d'un EVV à modulation MLI

ayant uniquement un filtre capacitif comme lien en
tension continue ........................................................................ Il 4
xv
Figure A.3 Schéma SIMULINK d'un système d'alimentation
d'un EVV à modulation MLI avec un filtre inductif
du côté du réseau ...................................................................... Il 5
Figure B.1 Schéma SIMULINK du modèle simplifié d'un EVV

vue du réseau de distribution ............................................. 118
Figure B.2 Schéma SIMULINK du système d'alimentation

triphasé avec des filtres passifs accordés et un EVV
comme charge ............................................................................. Il 9
Figure C.1 Schéma SIMULINK du système comprenant des

filtres passifs accordés et un EVV comme charge ....... 1 28
Figure D.1 Schéma de sortie du simulateur ......................................... 130
Figure D.2 Schéma principal du simulateur ......................................... 13 0
Figure D.3 Schéma principal du circuit de contrôle des

harmoniques ................................................................................ 1 3 1
Figure DA Schéma du circuit de contrôle pour un harmonique .. 131
Figure E.1 Résultats de la superposition d'un harmonique ........... 1 35
XVI
LISTE DES TABLEAUX
Tableau 2.1 Courant harmonique typique des redresseurs

triphasés ...................................................................................... 9
Tableau 4.1 Combinaisons d'harmoniques utilisées pour les

tests d'immunité ................................................................... 8 1
Tableau 4.2 Résultats des essais du GHT pour des exemples

de génération d'harmoniques de tension ................... 8 2
Tableau 4.3 Combinaisons d'encoches de tension utilisées

pour les tests d'immunité 19 ........................................... 89
Tableau E.l Résultats de simulation de superposition d'un

harmonique .......................................................................... 1 3 4
XVII
LISTE DES SYMBOLES ET DES ABRÉVIATIONS
Constante intermédiaire relative au calcul de

QtN
BEEPSEVV Banc d'Essai pour l'Étude des Perturbations Sur
les Entraînements à Vitesse Variable
CA/CC transformation d'un courant alternatif vers un

courant continu
condensateur du filtre
Condensateur du filtre N
Valeur total des condensateurs à utiliser pour

obtenir un coût minimum sur chacun des filtres
DHTIs Distorsion harmonique totale sur le courant de

la source
Distorsion Harmonique Totale sur la tension

ligne-neutre
EPRI Electric Power Research Institute
EVV Entraînements à Vitesse Variable
fI fréquence fondamentale du réseau
fn fréquence de résonnance du filtre n
FP Facteur de Puissance
Gill Générateur d'Harmoniques de Tension
XVIII
ho harmonique choisi pour le calcul de la
fréquence de résonnance du filtre passe-haut
HP horsepower
ishu, ishv, ishw composantes harmoniques du courant de la

source des phases u, v et w
ISU, ISV, ISW courants de la source des phases u, v et w
source de courant de l'harmonique h
Courant total traversant le filtre passe-haut.
valeur efficace du courant de ligne
courant harmonique h pénétrant dans la source
Valeur efficace du courant à la fréquence

fondamentale
15, 17 Courant traversant les filtre L5-C5 et L7-C7
jXCh impédance complexe du condensateur Ch
jXLh impédance complexe de l'inductance Lh
Je Coût des condensateurs en $/kVAR
JL Coût des inductances en $/kVAR
kVA kilovoltampère
kVAR kilovoltampère réactif
K gain du circuit de réglage du filtre actif
MLI Modulation de Largeur d'Impulsion
pu par unité "per unit"
QiN Capacité du filtre N en VAR
XIX
Facteur de qualité du filtre n
Qr puissance réactive reqUIse
Re X lm plan complexe réel et imaginaire
partie résistive de l'impédance de la source
Sn puissance nominale du transformateur
V L-N eff tension efficace ligne-neutre
composante fondamentale de la tension aux

bornes du condensateur
composante fondamentale de la tension aux

bornes de l'inductance
inductance de commutation
partie inductive de l'impédance de la source
YfhN admittance du filtre N
impédance caractéristique du filtre
impédance du filtre accordé a" la 5e harmonique
impédance du filtre accordé a" la 7e harmonique
li impédance totale des filtres passifs
Zph impédance du filtre passe-haut
impédance de la source
Zsh impédance de la source à la fréquence de

l'harmonique h
facteur d'incertitude sur la fréquence de

résonnance
xx
angle de l'impédance de la source
Pulsation de fréquence fondamentale (radis)
XXI
CHAPITRE 1
INTRODUCTION GÉNÉRALE
1.1 CADRE ET MOTIVATlON DU PROJET
Les travaux de recherches décrits dans ce mémoire de Maîtrise

s'inscrivent dans le cadre d'un projet de conception d'un banc d'essai
pour les Entraînements à Vitesse Variable (EVV) pour le Laboratoire
des Technologies Électrochimiques et des Électrotechnologies d'Hydro-
Québec (le LTEE d'H-Q).
Le laboratoire possède un banc d'essai pour la mesure du rendement

des moteurs à induction et désire étendre sa fonctionnalité aux tests
des EVV d'une puissance de 1 à 600 kV A (1 à 800 HP). Le banc d'essai
a pour objet d'étudier les problèmes d'émission des courants
harmoniques et de défaillance des EVV . faces aux perturbations
produites sur une ligne de distribution (test d'émission et d'immunité) .
La conception prend comme point de départ la description des essaIS

requis dans le document SC-610 de l'Electric Power Research lnstitute
(EPRl). Ce document décrit une série de protocoles de tests sur les EVV
de faible puissancd 1] . Elle repose également sur l'analyse des
informations contenues dans la littérature scientifique et sur des
travaux de simulation.
1
1 .2 SIMULATION NUMÉRIQUE DU BANC D'ESSAI
Ce travail couvre la partie du projet se rapportant aux simulations du

banc d'essai appliquées aux tests décrits dans le document SC-610 [1].
Le simulateur utilisé est le logiciel S IMULINK [2], qui est l'interface
graphique du programme MA TLAB [3]. Le Groupe de Recherche en
Électronique Industrielle (GREI) de l'Université du Québec à Trois-
Rivières (UQTR) a mis au point deux librairies de données dans
SIMULINK appelé SIMUSEC [4] et SEMALINK [5]. Ces librairies sont la base
du simulateur SIMUSEP (Simulation des Systèmes en Électronique de

~uissance) utilisé pour fin d'analyse temporelle de systèmes
électroniques en commutation. Partant de ce simulateur nous nous
proposons de développer les nouveaux modèles pour les éléments
constituant le banc d'essai "BEEPSEVV" (Banc d'Essai pour l'Étude des
Perturbations S ur les Entraînements à Vitesse Variable).
1.3 CONCEPTION DU BANC D'ESSAI
La partie critique de la conception et de la réalisation du banc d'essai

se retrouve dans la section décrivant les tests d'immunité des EVV aux
perturbations sur la forme d'onde de la source de tension. Le
document SC-610 suggère d'utiliser une source électronique
programmable raccordée directement à l'entrée de l'EVV pour
effectuer ce type de test. En pratique, dû au coût des sources de
puissance, cette façon de faire, bien que simple et élégante, n'est pas
économiquement réalisable sur des EVV de plus de 100 kVA. Nous
avons donc dû étudier des approches différentes pour créer un
2
système d'alimentation ayant les caractéristiques désirées malS pour
un coût bien inférieur au montage suggéré par le document SC-610.
Pour cela, nous proposons le montage présenté au chapitre 4 qUI a

pour objet de générer une forme d'onde de tension perturbée, à la
fréquence du réseau, tout en ne nécessitant qu'une source
électronique d'une puissance réduite. Ce montage et son circuit de
commande constituent la principale innovation technique des travaux
décrits dans ce projet. En raison de ce mérite, ce "générateur
d'harmoniques de tension" (GRT) fera l'objet d'une analyse plus
approfondie.
En plus de la conception du GRT, ce mémoire a également comme

objectif de développer des modèles de simulation numérique des
composantes constituant le banc d'essai. Ces modèles seront, dans un
premier temps, utiles à la conception du banc d'essai et par la suite,
utilisés pour les tests de compatibilité des systèmes (tests se limitant à
des simulations numériques d'une installation réelle[1]) .
1 .4 DESCRIPTION DES SUJETS TRAITÉS
Le chapitre 2 présente les calculs pertinents à la conception de filtres

passifs accordés par le développement d'un programme MA TLAB. Ce
type de filtre est utilisé dans le montage du GRT et dans un montage
classique afin de limiter la distorsion sur l'onde de tension pour les
tests de rendement et d'émission des EVV.
3
Le chapitre 3 présente un type de filtre actif dont la configuration a
inspiré la conception du GHT présenté au chapitre 4. On introduit dans
ce chapitre la méthode de calcul utilisée par le circuit de réglage du
GHT qui est le calcul des "puissances imaginaires instantanées"
("Instantaneous Imaginary Power" [6]). Nos résultats de simulation de
ce dispositif sont confrontés aux résultats de mesure publiés par le
concepteur de ce filtre. Ceci nous permet de valider notre modèle de
simulation avec des mesures expérimentales.
Le chapitre 4 présente la conception du générateur d'harmoniques de

tension qui sera utilisé pour les tests d'immunité des EVV. On y
retrouve les modèles de simulation développés pour ce projet avec les
résultats obtenus pour différentes perturbations typiques de l'onde de
tension. On y présente également les solutions suggérées pour
effectuer ou simuler les autres tests décrits dans le document SC-610.
Le chapitre 5 est une conclusion générale du mémoire. Les points

importants de la conception du banc d'essai y sont résumés. Les
simulations démontrent la faisabilité du projet de conception du GHT.
4
CHAPITRE 2
ÉTUDE DES PERTURBATIONS ÉLECTRIQUES PRODUITES PAR

LES EVV ET CONCEPTION DES FILTRES PASSIFS ACCORDÉS
2.1 INTRODUCTION
Les charges électriques non linéaire, telles que les entraînements à

vitesse variable, qui utilisent des convertisseurs de fréquence à semi-
conducteurs, sont la cause de la production de courants harmoniques.
Ces courants peuvent se propager à l'intérieur du réseau électrique
d'une usine et créer des perturbations à l'onde de tension. Ils sont
également susceptibles de créer des interférences électromagnétiques
par rayonnement ou d'engendrer d'autres problèmes reliés à la
conduction directe de ces courants. Le banc d'essai qui sera implanté
au L TEE couvre cet aspect du problème qu'est l'émission des courants
harmoniques produits par les EVV.
L'étude des émissions de courants harmoniques a deux utilités pour le

banc d'essai du L TEE. La première s'applique à la caractérisation des
émissions de courants harmoniques, selon une procédure standard [1],
pour évaluer les perturbations produites par différents convertisseurs.
La seconde consiste à utiliser les résultats de cette étude afin de
concevoir des filtres permettant de réduire l'amplitude de ces
courants harmoniques en deçà des limites acceptables [7, 8]. Des filtres
5
accordés sont utilisés dans les deux cas et c'est pourquoI nous avons
développé un programme de conception de filtres. Ce chapitre traite
de la conception de ce programme.
Le chapitre débute par un bref rappel des notions de base sur les
émissions de courants harmoniques avec les techniques envisageables
pour leurs réductions. Nous abordons ensuite la conception des filtres
passifs accordés à l'aide du calcul d'impédance dans le domaine
fréquenciel. Nous verrons de quelle façon jumeler les modèles de
simulation, dans le domaine temporel, à ce programme de calcul de
filtres de manière à vérifier les performances du système en régime
permanent ou transitoire avec des charges non linéaire et variable en
amplitude.
2.2 NOTIONS DE BASE SUR LES HARMONIOUES ET SUR LE CALCUL DES

FILTRES ASSISTÉ PAR ORDINATEUR
2.2.1 Perturbations de la forme d'onde de tension
Dans ce travail, on s'intéresse aux perturbations causées par les

charges électroniques qui peuvent modifier la forme d'onde de tension
d'alimentation du réseau électrique de distribution. Ces charges
peuvent être assimilées à des générateurs de courants harmoniques
dont les amplitudes sont pratiquement indépendantes de l'impédance
du réseau d'alimentation. La figure 2.1 montre un convertisseur CA/CC
représenté uniquement du point de vue des harmoniques de courant
produits par le convertisseur du côté du réseau électrique.
6
Transformateur Convertisseur Transformateur
~~-----f"-=--71
Rése:J~ I...--~_I_h~
Injection d'harmoniques
Figure 2.1: Convertisseur CC/CA du point de vue des harmoniques de

courant
Les perturbations typiques des convertisseurs classiques sont bien

connues. On les représente souvent comme des sources d'harmoniques
de courant. Cette représentation provient d'une décomposition
mathématique de forme d'onde, non sinusoïdale, en une somme de
signaux sinusoïdaux appelée série de Fourier telle que présentée à
l'équation 2.1.
00
l(t) = Co + L Ch cos(hW + eh)

h=J
J
(2.1 )
où: roI = 21tf l , fI est la fréquence fondamentale
h rang de l'harmonique
Co est la composante continue du signal

Ch est l'amplitude maximale de l'harmonique h
Sh est le déphasage de l'harmonique h
7
Par exemple, la figure 2.2a qUI présente la forme d'onde idéale d'un
convertisseur de courant à 6 impulsions peut être décomposée en une
série de Fourier tel que défini par l'équation 2.2. Ces valeurs peuvent
être utilisées pour la conception des filtres. Toutefois, dû à l'inductance
de commutation Xc ou au déséquilibre des phases, il peut y avoir des
différences importantes entre les valeurs idéalisées de ces formes
d'onde et les valeurs mesurées sur un réseau complexe. Ainsi il est
préférable d'utiliser des valeurs typiques de courant harmonique
lorsque celles-ci sont disponibles. La figure 2.2b présente un
graphique des valeurs typiques des courants de différents types de
redresseurs aInSI que les courants harmoniques théoriques du
redresseur à 6 impulsions. Ces données présentées au tableau 2.1,
proviennent des normes IEEE 519-1992 [8] qui se veut un guide sur le
contrôle des harmoniques et sur la compensation de la pUlssance
réactive consommée par les convertisseurs de pUIssance.
1.2
-
5 o 20 25 30
· 1.2
temps (ms)
Figure 2.2a: Courant de ligne théorique d'un redresseur triphasé de

6 impulsions à charge inductive
8
Les données présentées au tableau 2.1 correspondent à des
convertisseurs où le courant de charge est parfaitement lissé. Ces
données sont utiles mais elles ne sont pas représentatives des EVV de
faible ou de moyenne puissance qui utilisent généralement des filtres
capacitifs sur le lien à tension continue. Dans ce cas, on fera appel à la
mesure des courants ou à un simulateur fonctionnant dans le domaine
temporel, afin d'évaluer de façon plus précise les amplitudes des
courants harmoniques produits par l'EVV.
Tableau 2.1: Courant harmonique typique des redresseurs triphasés
Nombre Rang de l'harmonique
d'Impulsions 5 7 11 13 17 19 23 25
Théorique .200 .143 .091 .077 .059 .053 .043 .040

6
Typique 175 .110 .045 .029 .015 .010 .009 .008
Théorique .000 .000 .091 .077 .059 .053 .043 .040

12
Typique .021 .014 .075 .059 .0 Il .009 .015 .011
Théorique .000 .000 .000 .000 .059 .053 .043 .040

18
Typique .026 .016 .007 .004 .015 .010 .001 .001
Théorique .000 .000 .000 .000 .000 .000 .043 .040

24
Typique .026 .016 .007 .004 .002 .001 .009 .008
9
· sin(mt+ 4Jt) - (~Sin(5mt + 4J5)) - (~Sin(7mt+ 4J7)) +
l(t)={ilt (2.2)
(1\ sin(llmt + 4J1I)) + (1~ sin(13mt + 4J13)) - (...)
où: i(t) courant instantané théorique d'un redresseur idéal de
6 impulsions (Voir figure 2.2a)
00 00= 21tf 1, fI est la fréquence fondamentale de 60 Hz
h 5, 7, 11, 13 : rang de l'harmonique
Il Valeur efficace du courant à la fréquence fondamentale
On peut noter que les amplitudes typiques des harmoniques sont en

général plus faibles que leurs valeurs théoriques. Par contre, il n'y a
pas d'annulation complète des harmoniques de rang plus faible lorsque
l'on utilise des configurations de redresseur de 12, 18 ou
24 impulsions. Par exemple, pour un redresseur de 12 impulsions la
5e harmonique est théoriquement nulle alors que sa valeur typique
est de 2.1 % du fondamental.
10
Théo. 6 imp.
Typ. 6 imp.
5e 7e Typ. 12 imp.
11 e 13 e Typ. 18 imp.
17 e 19 e Typ. 24 imp.
Rang des harmoniques 23 e 25 e
• Typ. 24 imp . • Typ. 18 imp. • Typ. 12 imp.
o Typ. 6 imp. 0 Théo . 6 imp.
Figure 2.2b: Valeurs théoriques et typiques des harmoniques de

courant produits par différents montages redresseurs
11
2.2.2 Simulation des configurations les plus communes d'EVV
Dans la gamme de puissances d'EVV étudiée pour ce projet, les

convertisseurs les plus répandus utilisent un onduleur à commande de
type Modulation de Largeur d'Impulsion (MLI) alimenté par un
redresseur à diodes de 6 impulsions. Les différences les plus
marquées entre les différents convertisseurs commerCIaux se
retrouvent au nIveau du patron de commutation de l'onduleur
(fréquence, algorithme de réglage, etc.) et selon le type de filtre utilisé
pour produire le lièn en tension continue.
Dans ce chapitre, on s'intéresse aux harmoniques produits sur le

réseau; comme celles-ci dépendent principalement du circuit
redresseur nous étudierons et ferons la simulation numérique de
différentes configurations de filtres utilisés par le circuit redresseur.
La charge raccordée sera un EVV de 50 ou de 500 HP. Les simulations
de l'EVV de 50 HP ont pour objet de vérifier les modèles utilisés par
des mesures au laboratoire alors que les simulations de l'EVV de
500 HP permettent de vérifier les performances des filtres à la
puissance nominale du banc d'essai. Cette démarche permet de
montrer la méthode de calcul de filtre utilisé aInSI que de présenter
les différents modèles d'EVV développés à cette fin. Ces modèles sont
également utiles pour la conception du générateur d'harmoniques et
pour les tests de compatibilité des systèmes.
12
La figure 2.3 présente le premier type d'EVV étudié. Il s'agit du cas où
le filtre sur le lien à tension continue est complètement capacitif. Dans
ce cas seule l'impédance du transformateur d'alimentation électrique
limite l'appel de courant du condensateur Cf. Le modèle de simulation
numérique SIMUSEC [14] de cette configuration est présenté en annexe.
Les résultats de simulation sont présentés à la figure 2.4
Les principaux paramètres de simulation sont les suivants:
Source de tension: Transformateur 575 V, 3Ph, Impédance de fuite de 5 %

Filtre sur le lien en tension continue: Cf=2000 pF
Charge: Moteur à induction de 50 HP.
*
Lien à
Redresseur Ondule ur Moteur à induction
Source de
tension
avec
hnpédance I--_--I_~
Cf
-Ô 1-_ _--1
interne 1---+---1___ ---1
Figure 2.3: Schéma fonctionnel d'un EVV ayant uniquement un filtre

capacitif comme lien en tension continue.
Reponse temporel: courant de ligne Analyse frequenciel du courant de ligne
80~----~----~----~~
;;0.8 . .... :. ' " . '. ' ..... . .
20 ........ al
~0. 6 .. . . ... . .'. .. .. ..' ..... :. ....
: .
o E
Ci
-20 .... ... . .... . .. . . . . . .. . ~ 0.4
'0
E
-40 li
~ 0.2 . · · ···la DH1.-sur Il est'de 60.9 % ..
-60 ..
-80
0.4 0.405 0.41 0.415 5 10 15 20 25 30
temps (sec) harmonique
Figure 2.4: Courant de ligne du circuit de la figure 2.3
13
Pour cette topologie de circuit la forme d'onde du courant de ligne
dépend de la valeur de Cf et de Lf ainsi que de la charge. Dans ce cas
la façon la plus adéquate de calculer le contenu harmonique du
courant de ligne est d'effectuer une simulation numérique du circuit
dans le domaine temporel.
Les figures 2.5 et 2.6 présentent les deux autres configurations de

redresseurs étudiés. Dans chacun des cas, on ajoute une inductance
afin de limiter la crête du courant d'appel du condensateur. Cette
inductance peut être placée soit sur le lien à tension continue, soit à
l'entrée du redresseur du côté courant alternatif. Les nouveaux
modèles SIMUSEP de ces configurations sont présentés respectivement
aux figures A.2 et A.3 de l'annexe A. Les formes d'onde de courant,
tirées des résultats de la simulation numérique, sont présentées aux
figures 2.7 et 2.8. Dans ces configurations, le contenu harmonique de
courant dépend du choix d'inductances utilisées: une faible valeur
d'inductance nous rapproche des résultats de la figure 2.2 alors qu'une
valeur élevée d'inductance procure le contenu harmonique d'un
convertisseur de courant tel que présenté au tableau 2.1.
On utilise surtout la topologie capacitive (figure 2.2) pour les

convertisseurs de faible puissance «100 H.P.). Cette configuration, plus
économique, permet d'obtenir la tension moyenne la plus élevée et
requiert mOIns d'espace que les deux autres topologies. Les
convertisseurs de grande puissance utilisent plûtot la configuration
présentée à la figure 2.5 (avec filtre passe-bas LC). L'ajout de
14
l'inductance série sur le lien à courant continu réduit la distorsion sur
la forme d'onde de courant. La topologie de la figure 2.6 correspond au
cas où l'on ajoute un transformateur d'isolation à l'entrée de l'EVV à la
topologie de la figure 2.3. Cet ajout devient nécessaire si l'EVV est relié
à un transformateur présentant une faible impédance afin de limiter
le courant tiré par le condensateur Cf. Le choix de la topologie à

utiliser varie d'un fabricant à l'autre car les deux configurations
procurent le même effet bénéfique sur la réduction des harmoniques
de courant produites par l'EVV.
lien à tension
Redresseur continue Onduleur Moteur à induction
Source de
tension
avec
impédance
interne
L...-_ _----" L l ' - - - - -......
Filtre passe bas
Figure 2.5: Schéma fonctionnel d'un système d'alimentation d'un

EVV dans lequel on utilise un filtre Lf-Cf.
Redresseur Onduleur Moteur à induction

Source de LI
tension
LI
avec
impédance
interne
LI
* Cf
-Ô-
Figure 2.6: Schéma fonctionnel d'un système d'alimentation d'un
EVV avec un filtre inductif du côté du réseau
15
L'analyse spectrale des formes d'onde de courant produites par les
trois configurations étudiées et simulées, sera utilisée pour le calcul de
filtres présenté à la section 2.3. La définition du calcul du taux de
distorsion utilisé dans l'analyse spectrale est présentée à l'équation
2.3. Elle provient du guide IEEE Std 519-1992[8] qui aborde les notions
importantes à la compréhension et à l'analyse des convertisseurs
électroniques.
50
L Amplitude des harmoniques 2
Facteur de Distorsion = h=2 (2.3)

Amplitude de la fondarnentale 2
Les paramètres de simulation des EVV présentés aux figures 2.5 et 2.6
sont les mêmes que l'exemple précédent sauf pour l'ajout des
inductances Ll dans chaque phase.
Filtre sur le lien en tension continue [13]: Cj=2000 pF, L]=2 mH

60r------,.---~--~,........,
n O.8 . .
.... "., . . . . ... .. . . , .. .. ' ...... . .. .
a:
E
5. 20 . :s
Q)
Q)
c: ~0. 6 . .. .. . :. . . . .. .. . . .. ... .........
' '
. .
Cl
ca
E
~ 0 o
'0
~ 0.4 ........ . .. . . .... . .. .. ............ .
'E
~ -20 ......... , .. . ... . '0
:::> ~
8 C.
-40 .... .. ...... ... .. ~ 0.2 .. . . : ... · La DHTsur Il est'de 34.1 % ..
_60L----~-~~-~<...--J OLL
: n :nn : n n: n~ n
0.4 0.405 0.41 0.415 o 5 10 15 20 25 30

Figure 2.7: Courant de ligne du circuit présenté à la figure 2.5
16
60r-----~----~----~~
c:
40
li 0.8 - - - - - - - . - - --- - ..", " . .. . .. .
E
$ 20 al
al
C
Cl
~0.6
~ 0 §
"0 o
'E ~ 0.4
~ -20 "0
:J ~
8 ëi
-40 ------- . -;- ... . -.-.-:- ~ 0.2 - - - - -La DHT.-sur 1L est-dEt 27'. 0 % - .
-60 L -_ _ _ _----'-_ _ _ _ ~ _ _ _ _~___._l

°o~~~~~~~~--~~
0.4 0.405 0.41 0.415 5 10 15 20 25 30
Figure 2.8: Courant de ligne du circuit présenté à la figure 2.6
2.3 RÉDUCTION DES HARMONIOUES INJECTÉs SUR LE RÉSEAU
2.3.1 Généralités
Le taux d'harmoniques de courants injectés sur le réseau électrique de

distribution peut être réduit de plusieurs façons.
La première méthode consiste à réduire les courants harmoniques à la

source; par exemple, en utilisant des configurations du redresseur
(groupement série ou parallèle) permettant l'annulation de certains
courants harmoniques. On a vu par exemple, au tableau 2.1, qu'un
redresseur de 12 impulsions permet de réduire considérablement les
harmoniques de rangs 5 et 7. L'utilisation d'un transformateur
d'isolation, dans certaines configurations, peut également avoir le
même effet.
17
La seconde méthode consiste à utiliser des filtres afin de "trapper" les
courants harmoniques (filtres passifs accordés) ou à les compenser par
l'utilisation d'une seconde source de courants harmoniques jouant le
rôle d'un filtre actif (Voir chapitre 3).
La figure 2.9 présente un montage utilisant une série de filtres passifs

accordés disposés à l'entrée d'un redresseur triphasé. Ces filtres LC
augmentent la capacité de court-circuit entre les lignes (à certaines
fréquences) de manière à diminuer l'absorption des courants
harmoniques par le réseau.
Vs = S7S V 60 Hz. Zs ~ a ~~~m

~________-r~~~=!"~______VT
__________~=!~~----_
Filtres
passifs
Filtres accordés au Se et 7e harmonique + fitre passe-haut
Figure 2.9: Exemple de filtres passifs accordés
18
Le principal inconvénient de ces filtres est que leurs caractéristiques
sont fortement influencées par l'impédance de la source, Zs.
- Les valeurs des éléments réactifs des filtres peuvent entrer en
résonnance parallèle avec l'inductance de source et ainsi créer une
amplification des courants harmoniques se transformant en pointe

de tension.
Les filtres peuvent entrer en résonnance série, avec les
harmoniques produits par la source de tension et ainsi drainer des
courants beaucoup plus grands que ceux prévus pour compenser les
courants harmoniques produits par la charge non linéaire.
Les simulations numériques, dans le domaine temporel, permettent de
déceler ce type de problèmes. L'utilisation de filtres actifs raccordés
en série avec les filtres passifs élimine ces inconvénients. Ces filtres
sont présentés au chapitre 3 avec un exemple de problème de
résonnance série et de résonnance parallèle.
2.3.2 Conception des filtres passifs accordés
2.3.2.1 Calcul des filtres dans le domaine fréquenciel
Habituellement, le calcul des filtres se fait uniquement dans le

domaine fréquenciel, en régime établi et stable. On associe alors
l'impédance des inductances et des condensateurs aux représentations

complexes jXL et -jXc selon les équations 2.4 et 2.5. En général, on
19
représente les résistances R par une impédance réelle et constante sur
la plage de fréquence étudiée.
(2.4 )
1
-jX c = (2.5)
h j 2 n( h1;) C
où: h est le rang de l'harmonique étudié de 1 à 50,

est la fréquence fondamentale de l'onde de tension,
-jXc h et jXL h sont les impédances complexes du condensateur C
et de l'inductance L en régime sinusoïdal établi à la

fréquence harmonique h *f 1.
Le calcul des filtres utilise la théorie de superposition. Pour

chacun des 50 premIers harmoniques, on calcule un circuit dans le
domaine fréquenciel et on en déduit les tensions et les courants
produits par une source de courant harmonique. L'onde de tension ou
de courant peut ensuite être reconstituée dans le domaine temporel
par la sommation des harmoniques en utilisant l'équation 2.1.
2.3.2.2 Description du programme de calcul de filtre
Nous avons développé un programme pour la conception de

filtre dans le langage de programmation MATLAB fourni à l'annexe B. Ce
programme est simple car le langage MATLAB manipule directement les
nombres complexes utilisés par les équations 2.4 et 2.5. De plus, il
20
utilise les données de simulation du système provenant du logiciel
SIMULINK qui est l'interface graphique de MA TLAB. Le logiciel SIMULINK
calcul les courants harmoniques produits par l'EVV, selon les modèles
présentés aux annexes A et B, puis le logiciel MA TLAB fait le calcul des
filtres. Finalement, on vérifie les performances de ces filtres dans le
domaine temporel à l'aide du modèle présenté à la figure BI. Les
données transitent d'une façon naturelle d'un programme à l'autre
puisqu'elles proviennent du même lanaaae
1:> 1:>
de base. L'utilisation
conjointe de ces deux programmes, dans le domaine fréquenciel et
temporel, constitue un avantage important par rapport à la méthode
classique de conception de filtres.
2.3.2.3 Critères de conception des filtres
Pour ce projet les critères de conception des filtres dépendent du

but qu'on définit. Nous avons mentionné à la section 2.1 que ce
programme servirait à deux applications distinctes. La première va
aux tests d'émission et de rendement des EVV sur le banc d'essai. Dans
ce cas, le document SC-610 précise que l'on doit respecter une
distorsion harmonique totale (DHT) maximale de 1 % pour l'onde de
tension à l'entrée du redresseur. Le calcul des filtres tient compte de
ce critère qui est en fonction de la charge mise sous essai et de
l'impédance du transformateur d'alimentation électrique utilisé. La
seconde application est de fournir des recommandations aux clients
d'Hydro-Québec sur le choix de filtres à utiliser afin de réduire
l'émission des harmoniques de courant. Cela demande une conception
21
qUI mlmmlse le coût des installations tout en limitant l'amplitude des
courants harmoniques à un niveau acceptable selon la norme IEEE Std
5 19- 1992 [8]. De plus, ce type de filtre est souvent utilisé afin de
corriger le facteur de puissance de l'usine; ainsi le programme tient
également compte de ce critère de conception.
2.3.2.4 Étapes de design des filtres
Les étapes de conception des filtres sont énumérées aux points 1

à 10 présentés dans cette sous-section. Elles consistent à calculer, en
valeurs réelles, les éléments du circuit équivalent monophasé de la
source de courant harmonique de la figure 2.9 en fonction de la plage
de fréquences étudiées afin de répondre aux critères de design.
La figure 2.10 présente l'ajout de filtres accordés au circuit de la

figure 2.1. Ces filtres permettent de minimiser le courant Is h
pénétrant dans la source Vs et ainsi diminuer la distorsion de la
tension de ligne Vth (Vth = Ish Zsh).
Zsh
Vth
1fh 5 t1fh? tlfh hp

r- t r- I -1
I r-
1 IT 1 1 1
Vs ru
1
1Z5 1
1 1
1
1Z? 1
1 1
1
1Zhp
1
1
•Ih
1- --1 1- --1
Figure 2.10: Représentation symbolique monophasée d'un

convertisseur avec des filtres passifs accordés.
22
Étape 1 - Déterminer l'impédance de la source
On représente l'impédance de la source Zs sous forme d'un angle em

constant. Ceci correspond à une inductance pure reliée en série avec
une résistance dont la valeur augmente de façon linéaire avec la
fréquence f (Zs = R + j XL s' où R = k f RS). Cette valeur n'est
habituellement pas disponible mais, en première approximation, pour

les transformateurs de capacité supérieure à 500 kV A, la valeur de
XL s est au moins 5 fois plus grande que celle de la résistance Rs; soit
6 m >78°. La figure 2.11 montre cette représentation simplifiée (modèle

a) de l'impédance du réseau. La figure 2.12 présente un modèle (b)
plus complet malS plus difficile à implanter dans un simulateur, car il
nécessite de connaître RI, R2 et l'angle em. Nous avons décidé de
choisir le modèle de la figure 2.11 dans notre programme car il
permet de calculer facilement un facteur de qualité optimal des filtres
en fonction des variations des paramètres du système (Voir l'étape 8)
Dans cette étape du calcul, on doit également spécifier la puissance du

transformateur Snom en kV A, son impédance Zs en p.u., sa tension
ligne-neutre V l-n et la fréquence fondamentale du réseau f1 . Ces
données sont disponibles sur la plaque signalétique du

transformateur.
23
tx 181<8m
R ~
Figure 2.11: Impédance caractéristique du réseau modèle a) [9]
t
x
181<8m
R R2
Figure 2.12: Impédance caractéristique du réseau modèle b) [9]
Étape 2 - Caractériser la charge en terme de courants

harmoniques
La caractérisation de la charge peut se faire des deux façons sui vantes:
La première est de fournir au programme un vecteur avec le contenu

normalisé des courants harmoniques mesurés (amplitude et phase), la
valeur du courant de ligne, et le facteur de puissance de la charge. Le
24
vecteur des courants harmoniques peut également être calculé par le
programme dans le cas général de l'équation 2.2.
La seconde façon de procéder consiste à utiliser les données provenant
de la simulation du système (EVV + Transformateur+autres composants
du système) dans le domaine temporel selon les modèles fournis à
l'annexe A (modèle A.I, A.2 ou A.3). Dans ce cas, le programme lit les
résultats de simulation sauvegardés dans un fichier MA TLAB et fait
l'analyse spectrale nécessaire au calcul des filtres. Les données en
entrée proviennent de la simulation numérique du modèle S IMU SEP
correspondant à l'EVV et au moteur raccordé au banc d'essai. Un bon
choix de configuration d'EVV est disponible et présenté à l'annexe A.
Étape 3 - Calculer la pUIssance réactive à corriger
À partir du facteur de pUIssance calculé ou défini par
l'utilisateur, on calcule la puissance réactive qui corrige le facteur de
puissance du convertisseur à l'aide de l'équation (2.6).
(2.6)
où: FP est le facteur de pUIssance,
FP varie entre 0.9 et 0.95 pour un redresseur
6 impulsions à diode.
FP vaut 0.82 pour un redresseur 6 impulsions à thyristor

avec a. = 30°.
s= V L-N eff* IL eff, est la puissance apparente
25
monophasée en VA
V L-N eff est la tension efficace ligne-neutre
IL eff est le courant de ligne

Qr est la puissance réactive reqUIse en V AR
On utilise cette valeur de Qr dans le cas où l'on s'intéresse à la
correction du facteur de pUIssance en priorité sur la diminution des
courants ou tensions harmoniques.
Étape 4- _.- Calculer la valeur totale du condensateur à utiliser
Comme première approximation, on calcule la valeur de condensateur

CT à utiliser pour fournir la puissance réactive capacitive Qr à la
fréquence fondamentale fI selon l'équation (2.7).
(2.7)
où: Qr est la puissance réactive capacitive reqUIse en V AR
VIC est la tension aux bornes du condensateur
Étape 5 - Distribuer cette valeur de CT sur les filtres requis
Pour cette application nous avons choisi de distribuer cette valeur
totale de capacité sur trois filtres parallèles accordés aux 5 e et

7 e harmoniques et à un filtre passe-haut centré sur le
Ile harmonique.
(2.8)
26
L'équation (2.8) comporte une infinité de combinaisons possibles. Nous
avons choisi de distribuer la pUIssance réactive totale des
condensateurs en fonction de l'amplitude des courants harmoniques
traversant chacun des filtres et selon une procédure de minimisation
du coût des filtres [12] présentée aux équations 2.9 à 2.19.
La puissance réactive capacitive totale du filtre (tenant compte de la

composante fondamentale de la tension de ligne et de l'harmonique de
courant traversant le filtre) est donnée par l'équation 2.9
[2 v.2 [2
IÎ
~
= v.le2 m1C + _h _ = Q + _I_h
hm C 1 hQ
(2 .9)
1 1
où:
Qc = Puissance réactive capacitive totale du filtre (V AR)
VIC = Tension de fréquence fondamentale aux bornes du
condensateur (Volt)
(01 = Pulsation de fréquence fondamentale = 27t*60 (radIs)
C= Valeur du condensateur (Farad)
Ih = Courant de l'harmonique 'h' traversant le filtre
h= Rang de l'harmonique auquel le filtre est accordé
QI = ~2m1C = Capacité du filtre pour la fréquence fondamentale fI
(VAR)
On détermine la puissance réactive inductive totale du filtre (QL) en
fonction de la valeur du filtre Sf selon les équations 2.10 à 2.14.
(2.10)
27
À la fréquence de résonnance (h *f1) les valeurs de L et de C sont liées
par l'équation 2.11:
(2.11)
En remplaçant la valeur de L de l'équation 2.11 dans l'équation 2.10,

on détermine l'équation 2.12:
(2.12)
La tension VIL due à la fréquence fondamentale exprimée en fonction

de la tension VIC est donnée par l'équation 2.13:
(2.13)
Finalement la pUIssance réactive inductive totale du filtre en fonction

de la valeur du filtre Qf est déterminée en remplaçant la tension VIL
de l'équation 2.13 dans l'équation 2.12.
(2.15)
Dans les équations (2.10) à (2.15):

Qr. = Puissance réactive inductive totale du filtre (V AR)
V IL = Tension de fréquence fondamentale aux bornes de
l'inductance (Volt)
L= Valeur de l'inductance du filtre (Henry)
28
Le coût total du filtre est donné par l'équation 2.16 (si on néglige le
coût de la résistance du filtre).
(2.16)
où:
Je = Le coût unitaire du condensateur du filtre ($/kV AR)
JL = Le coût unitaire de l'inductance du filtre ($/kVAR)
Les équations 2.9, 2.15 et 2.16 permettent de calculer le coût total du

filtre selon l'équation 2.17.
"
Cout/oral J [I;~~](Je + J)
L
-- Sf ( Je +/1 ) + hS L
(2.17)
f
(2.18 )
où:
A=(Je + h2
JL )
Le coût minimal du filtre est obtenu lorsque la dérivée de l'équation

2.18 est égale à zéro. On détermine alors la valeur du filtre Qf avec
l'équation 2.19.
(2.19)
29
On applique les équations 2.18 et 2.19 pour chacun des filtres
sélectionnés afin de minimiser leur coût individuel. Afin de simplifier
les calculs nous avons traité chacun des filtres de façon individuelle en
supposant que l'harmonique de courant Ih5 passe uniquement par le
filtre LsCs, que Ih7 passent par L7C7 et que tous les autres courants de
haute fréquence passent par LhpChp.
(2.20)
(2.21 )
(2.22)
(2.23)
(2.24 )
(2.25)
(2.26)
(2.27)
(2.28 )
(2.29)
30
B =
hp
1121~~
11
(J +J )
C L (2.30)
(2.31)
(2.32)
(2.33)
(2.34)
où:
Je Coût des condensateurs en $/kVAR
h Coût des inductances en $/kVAR
~ Condensateur du filtre N
Qm Capacité du filtre N en VAR
AN, BN Constante intermédiaire relative pour le calcul de SfN
15, h Courant traversant les filtre L5-C5 et L7-C7
1h p Courant total traversant le filtre passe-haut.
CT$ Valeur total des condensateurs à utiliser pour obtenir un
coût minimum sur chacun des filtres (farad).
Si on a pas de contingence sur la correction du facteur de puisance de

l'installation on utilisera ces valeurs de C5, C7 et Chp comme valeurs
initiales de filtre, sinon on augmentera la capacité total de l'installation
de façon proportionnelle sur chacun des filtres selon les équation 2.35
à 2.37.
31
(2.35)
Sn
C7 =c·-
T S
(2.36)
fT
chp
= C . S jhp
T S (2.37)
fT
Étape 6 - Calculer les inductances accordées à ces valeurs de

condensateur
(2.34 )
(2.35)
(2 .36)
Étape 7 - Calculer le facteur d'incertitude [9] sur la fréquence de

résonnance des filtres
Ce calcul présenté à l'équation 2.37 dépend:
a) des dérives thermiques,

b) des fluctuations de la fréquence fondamentale fI ,
c) de la précision des spécifications des valeurs de L et C.
(2.37)
32
où: f n Est la fréquence de résonnance du filtre n. (In = 2nJr:c:);
LnCn
Est le condensateur du filtre n·,

Est l'inductance du filtre n;
Est le facteur d'incertitude sur la fréquence de résonnance
qui représente la variation de la fréquence de coupure
Étape 8 - Calculer le facteur de qualité des filtres
Ce calcul présenté à l'équation (2.38) optimise l'admittance des filtres

de manière à minimiser le taux d'harmonique sur la tension [9].
(2.38)
où: Qn Est le facteur de qualité du filtre n (Q= XL)

R
es n Est l'angle de l'impédance de la source (voir figure 2.11)
on Est le facteur d'incertitude sur la fréquence de résonnance
n = 5, 7 ou 11.
Normalement un facteur de qualité élevé est utilisé afin de réduire les

pertes. Dans ce cas, l'utilisateur peut imposer pour le calcul une valeur
supérieure à celle suggérée par l'équation 2.27 dans l'étape d'analyse
des performances des filtres.
L'utilisation d'un filtre actif, en série avec les filtres passifs accordés,
permet également d'augmenter le facteur de qualité des filtres.
L'analyse de cette configuration particulière fait l'objet du chapitre 3.
33
Étape 9 - Calculer la valeur de résistance parallèle du filtre
passe-hau t
Un filtre passe-haut est utilisé pour réduire l'émission des

harmoniques d'ordre supérieur à Il. Les valeurs de l'inductance et de
la capacité du filtre ont été calculées à l'étape 7; la valeur de la
résistance parallèle est donnée par les équations 2.39 et 2.40.
(2 .39)
R = cr Zo (2.40)
Les valeurs suivantes peuvent être choisies
0.5 < cr < 2

et ho 5: -J2 hrnin
où: Zo est l'impédance caractéristique du filtre
hm i n est la plus petite valeur d'harmonique à être filtrée
(ici hmin= Il)

ho est l'harmonique choisi pour le calcul de la fréquence
de résonnance du filtre LhpChp
Étape 10 - Vérification du respect des critères de conception
Le premIer critère de conception, qui est la correction du facteur de

puissance, a été rempli à l'étape 3. Si l'utilisateur désire plutôt limiter
le taux d'harmonique de courant ou de tension sur la ligne il est
possible que la capacité de filtres installés soit insuffisante. On calcule
34
alors les courants harmoniques pénétrant dans la source et les
tensions harmoniques de ligne à l'aide des équations 2.41, 2.42 et 2.43
selon le circuit équivalent monophasé présenté à la figure 2.10.
Impédance totale des filtres à la fréquence de l'harmonique h:
Iz 1--Iy
jl1r
jll~
1
+ yfh, + yfh hp 1
(2 .41 )
Courant harmonique pénétrant dans la source à la fréquence de

l'harmonique h:
(2.42)
Tension harmonique de ligne VTh à la fréquence de l'harmonique h:
(2.43)
Source de courant de l'harmonique h

y fh N Admittance du filtre N (5 , 7 ou hp) à la fréquence de
l'harmonique h (Y fhN= 1/ZfhN)
Z sh Impédance de la source à la fréquence de l'harmonique h
lsh Courant harmonique h pénétrant dans la source
35
On calcul les distorsions harmoniques totales de courant et de la
tension avec les équations 2.44 et 2.45 tirées de l'équation 2.2.
l'fIs "
FD = ~.:. :. . 2h=-=--:-'_' (2.44)
Is Is 2 1
50
I.Vr/
h=2
FD y T = (2.45)
VS 2
L'utilisateur du programme peut spécifier un taux de distorsion

maximum sur 1s et sur VT. Si ces conditions ne sont pas respectées le
programme augmente la valeur de CT calculée à l'étape 4.
L'augmentation de CT a pour effet de diminuer l'impédance des filtres
et ainsi de réduire les courants harmoniques pénétrant dans la source
(eq. 2.42 ) et la distorsion sur la tension de ligne (eq. 2.43). On refait
les étapes 5 à 10 jusqu'à ce que les courants pénétrant dans la source
I shT et la distorsion sur l'onde de tension VT respectent les critères de
design.
2.3.2.5 Organigramme du programme de calcul des filtres
La figure 2.13 présente l'organigramme du programme de calcul de

filtre. On y retrouve les étapes de conception de 1 à 10 énumérées à la
section précédente ainsi que les étapes d'analyse qui seront abordées
aux sections 2.3.2.7 et 2.3.2.8.
36
DÉBUT
•
Entrer les caractéristiques de la
source et de la charge:
S, Zs, VL, IL, f
•
Entrer les harmoniques de
courant produit par la charge:
Valeurs typiques
ou
Valeurs mesurées
ou
Valeurs provenant d'une
simulation numérique
•
Entrer la déviation maximum
sur les composantes L et C des
filtres et sur la fréquence
fondamentale de la source.
+
Entrer le coût des condensateurs
et des inductances des filtres en
$/kVAR
•
Calcul du facteur de qualité des
filtres accordés.
•
Calcul de Qr et CT requis pour
la correction du F.P.
cb
Figure 2.13 a Organigramme du programme de calcul des filtres
37
Distribution de CT sur les
filtre accordés au Se et 7e
harmoniques et sur un filtre
passe-haut.
Augmenter la
valeur de Ct de 5%
Calcul de la DHT sur le courant

et la tension de ligne.
OUI
Afficher les résultats
Module de l'impédance Impédance des filtres Analyse fréquenciel et

des filtres. dans le plan complexe reconstitution du signal
de IL et Vt
FIN
Figure 2.13 b: Organigramme du programme de calcul des filtres

(suite)
38
2.4 EXEMPLE DE CALCUL DE FILTRES ACCORDÉS
2.4.1 Critères de conception spécifiés
Puissance du transformateur (Snom): 500 kVA

Impédance du transformateur (Zs): 0.05 pu (0.0086+jO.0492 pu)
Tension de phase (Vp): 575/-./3 Volts
Courant de charge (IL): 400 Ampères
Fréquence fondamentale (fI ): 60 Hertz
Facteur de puissance (FP): 0.9
Distorsion maximum sur VL (DHTvd 5 %
Distorsion maximum sur Is (DHTIs) 15 %
Le courant de charge est représenté par le contenu harmonique

typique d'un convertisseur de courant de 6 impulsions.
2.4.2 Résultats du calcul de filtre
Filtre accordé au 5e harmonique: Cs=354 pF, Ls = 0.8 mH, Qs=15
Filtre accordé au 7e harmonique: C7=180 pF, L7 = 0.8 mH, Q7=15

Filtre pas se -haut: C p h=270 pF. Lph = 0.22 mH, Rph=1.1 .Q
Di storsion harmonique totale:

- sur la tension ligne-neutre : DHTVL _N = 4 %,
- sur le courant de la source : DHTIs = 14 %,

Puissance réactive totale installée: QTfilrre = 100 kVAR
L'annexe E donne le listage du programme FILTRE.m
39
2.4.3 Analyse des performances des filtres dans le domaine
fréquenciel
Le programme de calcul de filtre présente l'analyse des performances

des filtres, dans le domaine fréquenciel, avec trois séries de courbes.
La figure 2.14 montre l'amplitude de l'impédance individuelle des
filtres. La seconde courbe (figure 2.15) présente l'impédance totale du
système (filtre + transformateur) dans un plan complexe (Re X lm) et
en fonction de la fréquence. Cette courbe permet de visualiser les
pôles du circuit correspondant aux fréquences de résonnance série et
parallèle du système. La figure 2.16 présente une reconstitution des
formes d'onde de tension (VL-N) et de courant de ligne (IL) à partir de
l'analyse fréquencielle obtenue par le programme de calcul de filtre.
On présente également à la figure 2.16 une comparaIson de cette
reconstitution du signal avec les résultats de simulation numérique
provenant du simulateur SIMUSEP qui effectue ses calculs directement
dans le domaine temporel (Ce simulateur est présenté à l'annexe A) .
...... . . .. . ..- ....... ', -. -

... .. ..... .. ... , .................
,' '
. ... .... -. . ..... - .... . .. .. ,.
,' ' '
r · · · · · · · · · · · '" .. ... ' " . -. ','" '". ','. ' ,' - .. -

••••••••••• ".' •••••
.
'" ••• '"
,
••
.
'. ' . '". '," l' l'
.
Figure 2.14: Impédance individuelle des filtres
40
~ -1
O ,(G ~~
1 3~
· .
~ /......... ..: ... . .. .... ':' .... ..... .
02 : .
~-2
.§
.. . .. :. . . . :. :. ... .
" " " "
.. -3
Rang de I:harmonique
10-2,---~ _ _'- ' _._._._.~
. _____ ~ __...J
-4~------~--------~------~
1if 1~ 1~ o 0.5 1 1.5
Rang de l'harmonique Re(Zt) (Ohms)
Figure 2.15: Impédance du système dans le plan complexe
Comparaison des methodes de simulation: tension ligne-neutre

500
~
!!!
ic
J,
c
g
.§
~ simulation nùmerique dans lé :
. .
reconstitution du signal provena
-500
0 0 .005 0.01 0.Q15 0.02 0.025 0.03 0.035
Temps (sec)
Comparaison des methodes de simulation: courant de ligne

500
Ci:
.[
al
c
~
al 0
"0
'E
e!
::J
0
()
-500
0 0 .005 0.01 0.015 0.02 0.025 0.03 0.035
Temps (sec)
Figure 2.16: Formes d'onde de tension et courant de ligne
41
2.4.4 Analyse des performances des filtres dans le domaine temporel
La dernière étape du calcul de filtre consiste à vérifier ses

performances en utilisant le simulateur S lM U S EP dans le domaine
temporel. À partir des modèles numériques fournis en annexe on
introduit les valeurs calculées du filtre et on fait une simulation
jusqu'à ce que le régime permanent soit atteint. On récupère ensuite
les données dans MA TL AB afin d'en faire une analyse spectrale. À ce
niveau du travail, il est possible de vérifier l'effet d'une multitude de
perturbations au système telles que:
- déviation en fréquence de la source (frnoy +/- 1 %),

- variation de l'impédance de la source (Zrnoy +/- 0.02 pu),
- variation brusque de la charge (0-100 %),
- déviation des valeurs de C et L des filtres (C rnoy et Lmoy +/-3 %),
- ajout d'un deuxième EVV sur le transformateur,
etc ...
Si les filtres répondent correctement à ces différentes variables du

système, les valeurs recommandées par le programme seront
retenues. Sinon , on modifie les critères d'entrée du programme de
manière à cornger la situation.
42
2.5 CONCLUSION
Nous avons vu, dans ce chapitre, comment les EVV perturbent la

forme d'onde de tension de ligne du réseau électrique de distribution
par l'injection d'harmoniques de courant sur celui-ci. Des modèles de
simulation ont été développés afin de reproduire ces effets avec un
calculateur numérique. Nous avons également développé un
programme de calcul de filtres qui a comme particularité d'être couplé
à un programme de simulation, dans le domaine temporel, pour
obtenir les paramètres électriques de calcul et afin de vérifier ses
performances en régime dynamique.
Un exemple de calcul de filtre a été fait et des mesures sur un EVV de

60 HP ont permis de valider les résultats de simulation des modèles
SIMUSEP. Cette mesure fera partie des essais de mise en service du
futur BEEPSEVV.
43
CHAPITRE 3
ÉTUDE D'UN FILTRE ACTIF
3.1 INTRODUCTION
Les filtres actifs sont des dispositifs électroniques utilisés pour réduire
l'émission des harmoniques de courants produits par les charges non
linéaires.
Leur étude est motivée par le projet de conception du générateur

d'harmonique de tension (GHT) pour le banc d'essai pour EVV du LTEE.
Ce projet prend comme base bibliographique la littérature scientifique
sur les filtres actifs [10 -14 ] qui, de par leur configuration et leur
principe de fonctionnement, sont les dispositifs qui se rapprochent le
plus du GHT. L'étude des filtres actifs intéresse également le LTEE, en
tant que projet de recherche potentiel, puisque des essais pourraient
être menés à bien à des frais raisonnables, en utilisant les mêmes
équipements qui serviront au GHT.
Le chapitre débute par un bref rappel des notions de base sur les
filtres actifs. Puis, nous étudions un système composé de filtres actifs
et passifs reliés en série l'un avec l'autre. Cette configuration est
particulièrement intéressante puisqu'elle ressemble beaucoup à celle
du GHT. De plus, elle permet d'éliminer les inconvénients inhérents
44
aux filtres passifs accordés utilisés seuls. Finalement, on présente les
résultats de simulation de cas types de problèmes de filtres passifs
entrant en résonnance série ou parallèle avec l'impédance de la
source.
3.2 NOTIONS DE BASE SUR LES FILTRES ACTIFS
3.2.1 Configuration classique de filtres actifs [10]
Les filtres actifs injectent sur la ligne des harmoniques de

courant de même amplitude que la charge non linéaire mais en
opposition de phase avec ceux-ci. Les effets perturbants de ces deux
sources d'harmoniques s'annulent l'un l'autre et le courant tiré par le
réseau de distribution devient sinusoïdal.
Contrairement aux filtres passifs accordés ces filtres ont

l'avantage d'avoir des performances indépendantes de l'impédance du
réseau de distribution. Dans certains cas ils sont également
.... utilisés
pour corriger le facteur de puissance de l'installation.
Il est toutefois difficile de construire des filtre s actifs de grande

capacité en courant ayant une bande passante élevée (bande passante
requise à l'élimination des harmoniques de rang élevé). De plus, les
coûts d'installation et d'opération des filtres actifs sont élevés.
La figure 3.1 présente un exemple de filtre actif.
45
Zs Is IL Redresseur
~
VT ~
Charge
produissant
des
harmoniques
J CO
Filtre actif
ill
Figure 3.1: Topologie d'un filtre actif relié directement au réseau
3.2.2 Combinaison de filtres actifs et passifs [11]
Ce dispositif se veut un compromis entre les filtres passifs accordés et

les filtres actifs. On a vu au chapitre 2 que les performances des filtres
passifs, utilisés pour réduire l'émission des harmoniques de courant,
sont dépendantes de l'impédance du réseau. L'objet de ce nouveau
type de filtre est d'éliminer cet inconvénient en utilisant un circuit
électronique de puissance réduite (typiquement le circuit électronique
ne fournit que 1 à 2 % de la puissance réactive de la charge car 98 %
de la puissance réactive est fournie par les filtres passifs).
Les auteurs de ce nouveau dispositif[ll] ont étudié deux configurations

différentes de filtres. Dans sa première version le filtre actif est
46
disposé en série avec la ligne (figure 3.2) alors que dans sa dernière
publication l'auteur place le filtre actif en série avec les filtres passifs
accordés (figure 3.3). Les deux configurations procurent les
performances recherchées mais la configuration de la figure 3.3
possède certains avantages du point de vue de l'isolation électrique et
de la puissance nécessaire du filtre actif.
Ce système possède les avantages des filtres actifs classiques malS ne

nécessite qu'une faible partie de la capacité en courant de ces
derniers. L'analyse et la modélisation de ce dispositif, présentées à la
section 3.3, ont été utilisées lors de la conception du GHT. Cette étude
devrait mener à un projet de recherche au LTEE utlisant les
équipements du banc d'essai EVV.
Filtre actif
Impédance~zs ~
du transfo. 1-----i'V )----r-----r,'+-...,
t ~vv
Courant harmonique
produit par l'EVV
Vs
EVV
Tension du
réseau
Filtres accordés aux

harmoni ues 5, 7 et P-H
Figure 3.2: Combinaison d'un filtre actif relié en série avec la ligne et
de filtres passifs accordés reliés en parallèle.
47
Impédance~zs
du transfo.
1h
t EVV
Courant harmonique
produit par l'EVV
EVV
Tension du
réseau
Filtres accordés aux

~----+-- Filtre actif
harmoni ues S, 7 et P-H
Figure 3.3: Combinaison d'un filtre actif relié en série avec des filtres
passifs, l'ensemble étant relié en parallèle avec la ligne.
3.3 DESCRIPTION D'UN FILTRE ACTIF RELIÉ EN SÉRIE AVEC DES

FILTRES PASSIFS
3.3.1 Portée de l'étude
Dans cette section nous étudions le filtre actif développé par Fujita et
Akagi [6] et nous le comparons au système utilisant uniquement des
filtres passifs. Les expressions analytiques des courants harmoniques
circulant dans la source et dans les filtres sont développées dans
chacun des cas afin de démontrer les caractéristiques de filtration des
deux systèmes. De plus, le circuit de contrôle, la source de tension et la
charge sont simulés à l'aide du simulateur SIMUSEP. Les résultats de
simulation sont comparés à ceux présentés dans la publication de
Fujita et d'Akagi afin de vérifier notre méthode de simulation.
48
Les résultats sont présentés en grandeur réelle de manière à faciliter
leur comparaIson avec les résultats expérimentaux publiés.
L'impédance de base du système est de 2 ohms et la fréquence
fondamentale utilisée est de 50 Hz.
3.3.2 Description de la configuration du circuit
Le système étudié est présenté à la figure 3.4. Il est composé de filtres

actifs et passifs reliés en série l'un avec l'autre et ayant comme objet
de filtrer les harmoniques produits par une charge non linéaire.
Des filtres passifs accordés au 5e et 7e harmonique amSI qu'un filtre

passe-haut sont installés sur chacune des trois phases.
Le filtre actif est composé de six MOSFET commutés à une fréquence de

20 kHz en modulation de largeur d'impulsion. Un filtre passe-bas
(LR, CR) permet de filtrer les hautes fréquences venant de l'onduleur
produisant amSI une source de tension raccordée en série avec les
filtres passifs.
Typiquement le filtre actif ne requiert que 1 à 2 % de la pUIssance

apparente de la charge perturbante.
49
VI 200 V 50 Hz.
- Il
VT
- L Redresseur
Filtro posaif (10 kVA)
J LR
Vcc
J
Clrcun de puissance
Clrcun de contrOle
Calcul de ish calcul de P et q
GainK
Figure 3.4: Filtres actifs et passifs accordés reliés en série l'un avec
l'autre.
3.3.3 Description du circuit de contrôle
Le circuit de contrôle est présenté à la figure 3.4. Les courants de ligne

de la source, isu, isv, isw, sont mesurés et le circuit de contrôle en isole
les composantes harmoniques, ishu, ishv, ishw en appliquant un calcul
permettant de déduire les valeurs instantanées de puissance active p
et puissance réactive q.
50
Les composantes continues de p et de q (puissances active et réactive
moyennes totales) sont filtrées pour ne conserver que les composantes
alternatives soit P et q. De la transformation inverse de P et q il en
résulte les courants harmoniques ish u , ishv et ishw injectés dans la
source .
Les équations (3.1), (3.2) et (3.3) présentent la transformation des

courants et tensions de lignes en pUIssances actives et réactives
instantanées [6]
[::]=~[~ ~
1
2
2
-~
~
-2
rJ
ev
ew
(3 .1)
1
2 (3.2)
~
2
(3.3)
où: ea est la composante directe de la tension de ligne,
ef3 est la composante en quadrature de la tension de ligne,
iSa est la composante directe du courant de source,
iSf3 est la composante en quadrature du courant de source.
51
L'équation (3.4) présente la transformation inverse de p et q en
courants harmoniques de source après avoir filtré leur composante
continue avec un filtre passe-haut ayant une fréquence de coupure de
10 Hz.
iSh ] 11
~shv =~-Z-
• u 2
(3.4 )
[
lshw 1
2 2
Les courants harmoniques sont ensuite amplifiés d'un facteur K et sont

appliqués comme signal de commande de l'onduleur à modulation de
largeur d'impulsion (MLI). Après démodulation du filtre LR -CR une
tension Vc* = K ISh est appliquée en série avec les filtres passifs.
En plus, notons que le circuit de commande permet de bâtir la tension

continue aux bornes du condensateur du circuit onduleur. Un circuit
sommateur agit sur la composante q du signal, permettant d'introduire
une composante continue de pUIssance réactive (à la fréquence
fondamentale) sur le circuit de commande du modulateur d'impulsion.
Les courants provenant de la source de tension peuvent ainsi passer
de façon contrôlée à travers les filtres passifs et bâtir la tension
continue désirée (V cc *).
Cette tension est bâtie VIa les diodes antiparallèles sur les MOSFET et
asservie par la composante dq = Vcc*-Vcc à l'entrée du module de
calcul des courants harmoniques ish.
52
3.3.4 Analyse des performances du filtre actif
Le but de cette configuration de filtre actif est d'éliminer les

problèmes inhérents aux filtres passifs utilisés seuls. Dans cette
section on analyse les performances du circuit pour le système de
filtre passif avec et sans circuit actif pour des conditions de
résonnance série et parallèle avec l'impédance du réseau. Une
simulation, dans le domaine temporel, reproduit les résultats
expérimentaux du concepteur de ce filtre pour un problème de
résonnance série.
3.3.4.1 Problème de résonnance parallèle
Le circuit équivalent monophasé du circuit de la figure 3.2 est

présenté à la figure 3.3. Lorsque le filtre actif n'est pas activé (K=O),
les courants harmoniques de charge (IL h) sont distribués entre
l'impédance totale des filtres passifs Zf et l'impédance de source Zs
selon l'équation (3.5).
Zf
ISh = 1Lh (3 .5)
Zs+zf
Si la source a une impédance faible ( Zs"",O ) ou si IZfl > 1Zsl une grande
partie des courants harmoniques seront absorbés par la source. De
plus une résonance parallèle, à certaines fréquences entre Zs et Zf,
peut survenir et créer un phénomène d'amplification des courants
harmoniques de source.
53
Lorsque le filtre actif est mIS en opération et asservi comme une
source de tension, V c=K*Ish, celui-ci bloque le passage des courants
harmoniques à l'intérieur de la source Vs.
Is = Zf IL _ KIs h
h Zs+Zf 11 Zs+Zf
1s
h
(1 + Zs +K Zf J
=( Zf JIL
Zs + Zf h
Is -( zf JI4, (3.6)
h - K+Zs+Zf
On observe que, du point de vue analytique, le filtre agit comme une

résistance pure de valeur K. Le circuit de la figure 3.4 peut être
redessiné comme le circuit équivalent de la figure 3.5 que décrit
l'équation 3.6.
Zs
Vth
Ifh
l
TI
1 1Zf
1
1
L ILh
Figure 3.5: Circuit monophasé équivalent de la figure 3.4
54
Zs K
Vth
\r--
ï ~h
111
Ish 1 .) Zf
)1 , ILh
L~
Figure 3.6: Circuit équivalent du circuit monophasé pour ILh
La figure 3.7 montre les caractéristiques de filtration du système avec

et sans le filtre actif pour Zs=O.02 pu On observe que sans le filtre
actif il existe une résonnance parallèle, près du 4 e harmonique, entre
le filtre et l'impédance de source. Ainsi, SI la charge produit un
quatrième harmonique de courant, dû à un déséquilibre des phases ou
autre, ce courant sera amplifié d'un facteur important et produira des
perturbations de la forme d'onde de tension (résonnance parallèle).
A vec le filtre actif (K=2), ce problème est éliminé.
~ Or-~~~~==~~~~
:!2.
.J::
..J
-
~
<Il
- -10 .. :.. ~ . ~. \ ' ." ..
-20 L--'----'---"-'--'-_ _ _"--_"------'-----'------'---'---'----L--I..._ _ _-'--_"------'-----'--'--'--J........J

1
10° 10
Rang de l'harmonique (1111)
Figure 3.7: Caractéristiques du filtre pour les courants harmoniques
55
3.3.4.2 Problèmes de résonnance série
Un autre problème communément rencontré dans l'utilisation

des filtres passifs est l'apparition de résonnance série entre
l'impédance de source et l'impédance du filtre passif. Examinons par
exemple le cas où la charge ne produirait pas de courant harmonique
ou serait débranchée mais que la source de tension Vs ne serait pas
parfaitement sinusoïdale. Dans ce cas, celle-ci produirait des courants
harmoniques selon l'équation (3.7)
Is = VS h - KIsh
h Zs+Zf
Is
h
(1 + Zs+Zf
K J- _VS-",-h-
- Zs+Zf
Is = VSh
h Zs+Zf+K (3.7)
Le filtre actif est alors perçu par la source de tension harmonique

comme une résistance pure de valeur K en série avec les filtres
accordés. La figure 3.8 présente le circuit équivalent monophasé
décrivant l'équation 3.7.
56
Zs
Vth
ÎSii r ~Ifh
ITI
1Zf
Vsh '\,
U
K
Figure 3.8: Circuit équivalent du filtre actif pour Vsh
La figure 3.9 montre les caractéristiques de filtration du système avec

et sans le filtre actif en fonction des harmoniques de tensions V sh
produites par la source. On observe que sans le filtre actif il existe une
résonnance série près du 4e harmonique, entre le filtre et l'impédance
de source (pour Zs=0.02 pu). Ainsi, si la source de tension produit une
quatrième harmonique de tension de seulement 1 %, cette tension
résultera en une augmentation du courant harmonique à l'intérieur du
filtre d'environ 20 % si le filtre actif est absent (K=O) et de seulement
1 % si celui-ci est branché (K=2). Le filtre actif permet de régler les
problèmes de résonnance série créés par les filtres passifs seuls.
On présente à la section 3.3.6 deux exemples de simulation numérique

du filtre actif dans le domaine temporel. Le premier exemple
reproduit des conditions de résonnance série et le deuxième exemple,
un problème de résonnance parallèle. Les résultats obtenus
confirment les performances du filtre telles qu'elles ont été analysées
dans le domaine fréquenciel.
57
Les caractéristiques idéales du filtre sont atteintes lorsque l'on choisit
K égal l'infini. On obtient alors:
Ish = 0,
VTh = 0,
Vc = ZF*ILh + Vsh.
Cependant, en réalité, un gam K trop élevé peut causer des problèmes
d'instabilité. Les simulations numériques ont démontré que pour le
circuit mis à l'essai (figure 3.4) avec un gain supérieur à 2.5 le filtre
actif devient instable et engendre ou amplifie les perturbations de
haute fréquence.
Les caractéristiques du système présentées aux figures 3.7 et 3.9 ont

été simulées selon la formulation analytique des équations 3.6 et 3.7
dans le domaine fréquenciel avec le logiciel de calcul MA TLAB. Le
fichier de simulation des filtres est fourni à l'annexe B. La valeur de
base de l'impédance pour cet exemple est de 2 n.
. .. " .'. , .. ,., ... " .. , .... ,." ,. " " . . " . . .. . . . ,' , '
10'
: : : : : : : ! : : : : : : : ~ : .: :: ::':: ~ ~ ~ ~ ~:~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ 1 ~ 1 ~ ; ~ ~ ~ ~ ~ ~ ~ ~ ~ 11 ~ 1~ ~ ~ 1:111 ~ 11 ~ ~
" , . : : : : :' : : : : : : : ,
:~~.;....;...;..;....
;....;.
"
...;... ..
.. ;....; . .;...;
.. .. . .......
l ' · · ' · ' · ."
...;..
~ . ................
• •• ,. _ . " .",.,
.
"., .. . . ' , . , . , ... ".' ... '", ..... '. .". ..., .. . .
, ." ..' " .,. , .,
.::k~b : ~ ~ ~
...........
:! ~ ~ ~ ~~
.... - ., ....
-,
! ! ~ : ~ : ::: : : : : : : ::: : : : : : : ! : : : : : : : : : : : : ! : : : : ~ ~ :
-,., ..... . , .. .. , ' . " . , . . ............... .
,
!~!~
,
.. : ,!l~l~~
,
.. , . , !~
.. ' . ,
:~ ~!~
.
.",. ,
, , , . ,' , .. ,
...
,. ,
::: K=~ :::: :: ::::: : ::::: ::::::::: ::: :::::::::::::::::::::::::::

10-2 --'-1_ _--'-1_ _---1.1_ _--'
2 4 6 8 10 12 14 16 18 20
Rang de l'harmonique (1111)
Figure 3.9: Caractéristiques de filtration en fonction de Vsh
58
3.3.5 Simulateur du filtre actif de Fujita et Akagi [6J
Le circuit de contrôle du système, tel que décrit par les équations 3.1 à
3.4, a été simulé dans le domaine temporel avec le logiciel SIMULINK.
Le schéma fonctionnel de ce simulateur et les paramètres de
simulation sont présentés à l'annexe C.

parallèle
Les problèmes de résonnance parallèle surviennent lorsque

l'impédance de la source entre en résonnance parallèle avec
l'impédance des filtres accordés tels que présentés à la figure 3.10. Les
harmoniques de courant produits par la charge à cette fréquence de
résonnance sont alors amplifiés par les filtres au lieu d'être atténués.
Cela peut même aller jusqu'à perturber la forme de l'onde de tension
et générer des problèmes de résonnance série sur les transformateurs
reliés sur le même lien en tension.
Nous avons reproduit le problème de résonnance parallèle décrit dans

l'article de Fujita et Akagi [6] afin de démontrer les performances du
filtre actif dans le domaine temporel avec un exemple réel.
La figure 3.10 présente les résultats de simulation du circuit de la

figure 3.4.
59
Courant de source (Is)
80,-----,-----~----,_----,_----,_----_,----_r----_,----_.----~
ê
~ 0
G)
." -20
'E
8'" -60
~ -40
. .. "
~ ....... '. ... .... '. ' ...... ,
', ' ...
-80~~~~--~~--~~--~~--~~--~~--~~--~----~----~
M ~ ~ ~ ~ M ~ ~ ~ 0,58
Temps (sec)
(a)
Is sans fittre act~ Is avec filtre actif

60r-~------~------r-----~--~ 60.--,------~----~----~__,
Ci: 40 Ci: 40
~
E
~
ê
i
.g
:::l
5l
G)
."
~-20 ~ -20
g 5
o
() -40 () - 40
_60L-~----~----~----~~ -60~~----~~----~----~~
0,42 0.43 0.44 0.45 0.56 0.57 0.58 0,59
Temps (sec) Temps (sec)
(b) (c)
Analyse frequenciel du courant de ligne

w,----,-----,-----,----,-----,-----.-----r----,-----.----~
OHT = 21 . ~ %
40 .... .. . , ...... . . ..: • .. • , .. .OI:1T.=: 7.2 % .....
Ci:30 ... . .. : ... . ....' .. . .. . .. ' ... .... . '. .

.[
~20 ...... ~ ....... " , . ..... ' ." .... .. "-
,
.. ... . .. . .. .. ...... . . . . .. ...... . -.. ... ... . ,, . . . . .. .
.
............. .
.... .. .. ... . ... .
Avant la mise en marche du filtre
10 .. ' !ipres liùi1i~e 'eiùi1Mche (fu·filtre ':' ...... ':' . . .
O~~~~~~~=i~~-u~~~~~l-~----~--~~~~----~
o 2 6 8 10 12 14 16 18 20
Rang de l'harmonique
(d)
Figure 3.10: Simulation d'un problème de résonnace parallèle
60
3.3.5.3 Exemple de simulation d'un problème de résonnance série
Les problèmes de résonnance série surviennent lorsque la source de

tension n'est pas parfaitement sinusoïdale et qu'un harmonique de
tension entre en résonnance série avec l'impédance d'un filtre accordé.
Il apparaît alors un courant harmonique de source même si la charge
ne produit pas d'harmonique de courant.
La figure 3.11 présente les résultats de simulation lorsque la source de

tension possède une 4e harmonique de tension ayant 1% de
l'amplitude de la fondamentale. Il apparaît alors une distorsion
harmonique totale de 9 % du courant tiré par la source. Lorsque le
filtre actif entre en action cette distorsion est réduite à 1 %.
Les conditions de simulation pour le problème de résonnance parallèle

présenté à la figure 3.10 et pour le problème de résonnance série
présenté à la figure 3.11 sont les suivantes:
Source de tension: Transformateur 200 V 3Ph, 50 Hz 20 kVA 0.02 pu

Filtre au Se, 7e et passe-haut:
C5=340 )JP, L5=1.2 mH, Q5 = 15
C7=170)JF L7=1.2 mH, Q7 = 15
Cph=300 )JP, Lph=0.26 mH, Rph = 3 Q.
Charge pour la figure 3.10:

Moteur de 20 HP allimenté par un EVV en MU à 20 Hz.
Charge pour la figure 3.11:

Résistance triphasée de 2 Ohms relié en étoile (20 kW )
61
Courant de source (Is)
~ 0
~
~ -50
ë
f!
8 -100 . .
()
Avant la :mise en mElrche du filtl'\3
-150
~~~ _ _ _ _- L_ _ _ _~_ _ _ _L -_ _~_ _ _ _- L____~_ __ _L -__~____~
0.4 0.44 0.46 0.48 0.5 0.52 0.54 0.56 0.58
Temps (sec)
(a)
Is sans fi~re actW Is avec filtre actif
0:
E
~
~
::>
5l
al
'0
ë
f!
::>
0
()
0.42 0.43 0.44 0.45 0.56 0.57 0.58 0.59

(b) (c)
Analyse frequenciel du courant de ligne

1oor-----,------,-----,-----,,-----,-----,------,-----,------,-----,
80
DHT = 9.3 %
0:60 ....... : ........' . .. ... . . ' . . . . ... -;........ : .. . . . . . : ...... . .:- . . .... ·OHT·= ·1'.2 % " ..
!
:ïii 40
20 .. : ..... .. .' .. ...... ' ...... .. ', ... . .. . " ....... : ...... . ,' ........ ' ....... . . . . .. . .
Avant la mise en marche du filtre :
2 4 6 8 10 12 14 16 18 20
(d)
Figure 3.11: Simulation d'un problème de résonnance série
62
3.4 CONCLUSION
Il est difficile de juger des performances réelles de ce système

uniquement avec les données fournies dans la publication. Les
simulations numériques du système, avec un circuit de réglage
utilisant strictement les équations développées dans la publication, ont
fait ressortir certains problèmes affectant les performances du filtre.
Lorsque la tension de ligne n'est pas parfaitement sinusoïdale, la
reconstitution des sIgnaux de commande par la méthode des
pUIssances imaginaires instantanées passe par des points de calcul
singuliers (lors d'une inversion matricielle). Cela nuit aux
performances du filtre mais n'empêche pas le fonctionnement car ces
discontinuités sont aténuées par la bande passante limitée du filtre.
L'auteur de l'article n'aborde pas ce problème et ne spécifie pas SI un
traitement numérique est appliqué pour rendre la commande plus
robuste.
L'idée de base, qui est de rendre les performances des filtres

passifs accordés indépendantes de l'impédance du réseau et de la
distorsion existante sur la forme d'onde de tension, est excellente. Une
étude plus poussée, par des simulations numériques et avec la source
électronique qui sera utilisée par le ORT, nous permettrait d'évaluer le
potentiel réel de ce type de filtre.
63
CHAPITRE 4
CONCEPTION D'UN GÉNÉRATEUR D'HARMONIQUES DE TENSION
4.1 INTRODUCTION
La qualité de l'onde de tension distribuée par le réseau n'est pas

toujours parfaite. Certains de ces défauts de qualité sont susceptibles
d'influencer le fonctionnement d'appareils électroniques, tels que les
EVV, qui y sont particulièrement sensibles. La susceptibilité de ces
équipements nécessite de prendre des dispositions pour limiter soit ce
que l'on peut définir comme étant une source de pollution, soit les
inconvénients dus à cette pollution.
Dans ce chapitre on aborde le problème de compatibilité

électromagnétique entre le réseau d'alimentation électrique et les
EVV. Le traitement du sujet comportera trois parties distinctes. Soit:
1 - L'étude des perturbations susceptibles d'affecter la forme

d'onde fournie par le réseau (harmoniques, surtensions, sous-
tensions, fautes, etc.).
2 - La conception d'une source permettant de produire des formes
d'onde de tension harmonique du 60 Hz appelé "le Générateur
d'Harmoniques de Tension" (GHT).
64
3 - La simulation numérique du fonctionnement du GHT utilisé
afin de reproduire des perturbations typiques de la forme
d'onde de tension.
Tout comme le filtre actif présenté au chapitre 3, le GHT utilise un

calcul des puissances imaginaires instantanées pour la transformation
des tensions de ligne en leurs composantes harmoniques réelles et
imaginaires. Toutefois, la similitude apparente entre les deux systèmes
s'arrête là. Le pricipe de fonctionnement, l'utilité et la mise en oeuvre
des deux dispositifs sont totalement différents et ne doivent pas être
confondus. Le GHT n'est pas un filtre actif et n'est pas utilisé pour
réduire l'émission des harmoniques de courant ou la correction du
facteur de puissance. Il agit sur la forme de l'onde de tension qui peut
être réglée avec un taux d'harmonique quelconque; incluant un THD de
o% correspondant à une onde de tension sinusoïdale.
4.2 ÉTUDE DES PERTURBATIONS DE LA TENSION DE LIGNE
Les défauts à la qualité de la tension peuvent se manifester sous

différentes formes et affecter les quatre paramètres définissant les
caractéristiques du système de distribution électrique qui sont: la
fréquence, l'amplitude , la forme d'onde (sinusoïdale) et la symétrie
des systèmes triphasés (égalité des modules des tensions et de leurs
déphasages relatifs). Il existe donc quatre possibilités distinctes de
perturbations et leurs combinaisons éventuelles. Elles seront
développées aux paragraphes 4.2.1 à 4.2.4 ci-dessous en fonction des
tests d'immunité prévus dans le document SC-610 [1].
65
4.2.1 Fluctuation de la fréquence de 60 Hz
Elles sont rares et ne sont observées que lors de circonstances

exceptionnelles telles que les défauts graves en réseau au niveau de la
production et du transport. L'étude de ce type de perturbation ne
comporte donc pas d'intérêt particulier et n'est pas couvert par les
tests décrits dans le document SC-610 [iL
4.2.2 Variation de l'amplitude:
Il ne s'agit pas de variations lentes de tension qUI sont corrigées par

les régleurs de charge des transformateurs, mais de variations rapides
de tension, ou de creux de tension se présentant souvent sous forme
d'à-coups brusques ou isolés soit plus ou mOlfiS répétitifs, de forme
régulière ou non.
Le document SC-610 présente quatre tests d'immunité se rapportant

aux variations de l'amplitude de l'onde de tension. Les deux premiers,
Il et 12 vérifient le fonctionnement de l'EVV lorsque l'amplitude de
l'onde de tension, en régime établi, est inférieure ou supérieure à la
tension nominale du réseau. Ces tests peuvent être réalisés en utilisant
un variac triphasé déjà existant sur le banc d'essai moteur. Les deux
autres tests 16 et 17 correspondent à des baisses ou des montées de
tension momentanées telles que présentées aux figures 4.1 et 4.2. Ces
variations de tension sont trop rapides pour être exécutées par le
variac qui possède une constante de temps mécanique assez grande
(environ 6 volts/sec). Pour reproduire ce type de perturbation, on
66
compte utiliser une source de tension électronique triphasée reliée en
série avec la tension de ligne fournie par le transformateur. La figure
4.3 présente le schéma fonctionnel suggéré pour ces tests. La
réalisation de ces tests ne présente pas de problème technique
particulier et ne fera pas l'objet d'une analyse plus approfondie dans
ce travail.
T3
200
J T
Temps (msec)
Figure 4.1 Exemple de test d'immunité aux creux de tension (Test I6)
67
T.1. ~-' ;111111 T2 ~I T3 ~I
1.5 '{\
Crête de""
A A 1
- --n
1
1
tension 1 :
n n
~
~~ A
-
:::l
Cl..
0.5
1
-
1
'-"
0
tv
> 0 E0 1 )0 2 00
-0.5
-1 V V V
V V
-1.5 V V
Temps (msec)
Figure 4.2 Exemple de test d'immunité aux crêtes de tension (Test I7)
Réseau de
Moteur à
distribution
EVV induction
Vs Zs
Source électronique
avec trois sorties
monophasées
Figure 4.3: Schéma de base pour les tests d'immunité I6 et I7
68
4.2.3 Variation de la forme d'onde de tension:
Cette onde n'est plus sinusoïdale et elle peut être représentée par une
onde fondamentale à 60 Hz associée souvent à des harmoniques de
fréquences multiples du 60 Hz, mais aussi, parfois, à des ondes de
fréquences quelconques de valeurs supérieures ou inférieures à 60 Hz.
Le document SC-6l0 présente cmq tests d'immunité de ce type. Trois

d'entre eux traitent des harmoniques de tension d'entrée en régime
établi et deux autres tests en régime transitoire sur les perturbations
créées par la commutation des condensateurs et sur des encoches dans
la forme d'onde de tension de ligne. Ces cinq tests font 1'objet de la
section 4.3 .
--
:J
0-
-
>
0
-0.5
0
T
30
~ 40
vJ
1
1
(\ '
~ !\
-1 - 1 v j
-1.5 -
Temps (msec)
Figure 4.4 Exemple de pollution de la tension de ligne par des

harmoniques (Test 13, 14, 15)
69
1
0 .8
0.6
0.4
--
::l
Cl.
0.2
0
:> -0.2 20 40
-0.4
-0.6
-0.8
-1
Temps (msec)
Figure 4.5 Exemple de test d'immunité aux encoches dans la forme

d'onde de tension de ligne (Test 19)
1.5
--
: :l
Cl.
0.5
0
:>
-0.5
-1
-1.5
0
t
.1
1
1 20 30 40
Temps (msec)
Figure 4 .6 Exemple de test d'immunité à ·la commutation des

condensateurs (Test 18)
70
4.2.4 Dissymétrie du système triphasé (déséquilibre des phases)
Ce type de problème peut être reproduit à l'aide du montage de la

figure 4.3. Toutefois, ce type de défaut ne fait pas parti des tests
d'immunité décrits par le document SC-610 et ne sera donc pas étudié
dans ce travail.
4.3 DESCRIPTION DU GÉNÉRATEUR D'HARMONIQUES DE TENSION (GHT)
4.3.1 Tests d'immunité et utilité du GHT
Les défauts de la forme d'onde de tension décrits aux paragraphes

4.2.1 à 4.2.4 sont causés par différents types de perturbations du
réseau de distribution. Par exemple:
Un court-circuit ligne-neutre sur le réseau produira une baisse de

tension momentanée jusqu'à ce que les contacteurs d'isolation
entrent en opération.
La commutation de groupes de condensateurs perturbera la forme

d'onde de tension pour quelques cycles.
Un problème de résonnance parallèle, entre les condensateurs de

correction du facteur de puissance et l'impédance du réseau,
produira une amplification des courants harmoniques générant
des harmoniques de tension sur la ligne de distribution.
71
Ces problèmes peuvent créer des défaillances des EVV pouvant
résulter en des manoeuvres intempestives du système de commande
électronique (pouvant être dues au décalage du passage par zéro ou à
la crête de l'onde résultante de tension) et des claquages dus à des
surtensions de la valeur crête.
Les tests d'immunité consistent à vérifier en laboratoire le

comportement des EVV lors de ces perturbations du réseau. Il n'est
évidemment pas réaliste de tenter de recréer les causes des
perturbations malS il est possible de recréer leurs effets, c'est-à-dire
la déformation de l'onde de tension de ligne.
Une source électronique, ayant une bande passante suffisante, peut

être utilisée à cette fin, toutefois dû au coût élevé des sources cette
solution n'est pas économiquement viable pour tester des EVV d'une
puissance supérieure à 150 kVA. Pour les tests d'EVV de 150 à
600 kV A nous avons développé un système consistant en une mise en
parallèle d'une source électronique et du réseau de distribution. Cette
configuration prend avantage du fait que toutes les formes d'onde qui
doivent être reproduites ont comme point commun d'avoir une
fréquence fondamentale fixée à celle du réseau de 60 Hz. De plus,
l'amplitude de la fondamentale de la plupart de ces formes d'onde
demeure constante lors du test. Nous avons appelé ce dispositif le
"générateur d'harmonique de tension" (GHT). Celui-ci permet de régler
l'amplitude et la phase des 60 premiers harmoniques de tension (sauf
la composante continue). Il est donc possible de reproduire les formes
72
d'onde reqUIses par les tests d'immunité en passant par une analyse
spectrale des formes d'ondes et une reconstitution du signal dans le
domaine temporel par une sommation des harmoniques selon
l'équation générale 2.1.
4.3.2 Principe de fonctionnement du GHT
4.3.2.1 Transformation des harmoniques de courant e n

harmoniques de tension
Afin de comprendre le fonctionnement du GHT il peut s'avérer utile de

rappeler comment les harmoniques de tension peuvent être produites
de façon naturelle sur le réseau de distribution. La figure 4.7 présente
le schéma fonctionnel monophasé d'un système composé d'une charge
non linéaire reliée à un transformateur de distribution ayant un
condensateur de correction du facteur de puissance.
La figure 4.8 présente l'impédance Zs en fonction de la fréquence. On a

choisi, dans cet exemple, une valeur de condensateur qui entre en
résonnance parallèle avec l'impédance interne de la source. À cette
fréquence ( Ile harmonique ) l'impédance Zs devient très grande et le
pro duit "h armonique de courant" et "impédance de source" donnera
naIssance à des harmoniques de tension de forte amplitude
(V t Il . = Il 1 * Zs Il). La figure 4.9 illustre la tension de ligne de ce
système lorsqu'il alimente un convertisseur de courant typique (ayant
une Ile harmonique d'amplitude Al/11).
73
~rn:Jnsforf --.:~ ::;o~:c~sdel 1
Reseau 1_ => "'- 1
= z~ Vth
Condensateur de ~:- 1
correction du F.P.
Injection d'hannoniques
Figure 4.7: Représentation d'un problème de résonnance parallèle .
.... .... ...... : ........ .. .... ...... .. .-. .. .. . . . . " . .. . .. .. .... ':' .. ..... : ..... -:' ... :' .. : .. -:: ' . :'. .
.... ...... ... ... ..... _......... ..-.. .... .. . . .. ... . . ... . ...... _. . . .. . . . .. .. .. ... ... ... . .. . ... ..... .
~
"
0" 0 ::: : : : : : : :
~- ::-, !!:, ",:,!:!!,1!',!,::!!!::!:,::,:1:::,!,:::/! rI::!!!!!!!!!:!!::i,:::!::!1:::!!l'!!: r" i!' 1':l!' J

~ ~ ~~ ~~ ~~ ~~~~ ~ ~~~. ~ ~ ~~ ~ ~ ~~. ~ ~ ~~ ~~. ~ ~ ~ ~~~. ~ ~~~. ~~~. ~ ~~. ~ . ~~~. ~ ~ ~ ~ ~~ ~ ~ ~ ~ ~ ~ ~ ~~~o
~ ~~~~.
~ ~ ~~ ~ ~ ~ ~ r~::
. ~~. ~ ~ ~~. ~ ~~. ~.~~ ~ ~
10-2
• • • •
.. .. .. .. .. .. .. .. .
o • •• o • •• • • f • • • o • • • '.' ••••• ~ • •• ' .' •• ' , ' • • , •• ,
.
.
.
.. .
. ',' ••••• •
..
. .. .. ..
•• o' • • • 0.' . . . . . . . .. . . . . . 0.' . . . . . . . . . . . . '. ' • ... .
. .
•
10° 10'
Rang de l'harmonique h
Figure 4.8: Module de l'impédance Zsh vu du convertisseur.
Tension harmonique ligne-neutre

2'---~----~----'----T----~----~--~----~--~
-s
~
co
c
,g'
CI) .., .. . . ... . . ,
"0
c
o
eii .. .... . .. . . .. .. . .. ... . . . . ...
cCI) -1
1-
-2~--~----~----~----~----~----~----~----~--~
o 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018
Temps (sec)
Figure 4.9: Forme d'onde de tension de ligne résultant du produit

Zs * rh tel que présenté à l'équation 4.1.
74
00
Vt(t) = L.,Ih cos(hav+ eh) * ZSh (4.1)

h=l
Le GHT utilise un principe similaire. Un filtre accordé à la fréquence

fondamentale est ajouté à la sortie du transformateur de distribution
de manière à produire une impédance élevée pour tous les
harmoniques de courant (ZSh» Zb as e). Une source de tension
électronique reliée à la charge au travers de filtres passifs permet de
régler l'amplitude et la phase des harmoniques de courant qui
produiront la tension harmonique de ligne désirée Vth (Vth= Ih * ZSh)·
La figure 4.10 présente le circuit complet du GHT.
A u t~
transformateur à
tension de sortie
variable
0-700 Valls
Filtre accordé à 60 Hz.
10--
L_J-----T----t---j
vr
-- IL
ASD SOO HP
Moteur à induction
60 Hz. 500 KYA

Fr=60 Hz.
Filtres passifs accordés

Se,7e, Ue, 13e et passe-haut
(200 kVA lotal)
Source
électronique
programmable 1--_ _---'
0-120 Volts
100kVA
~In:u~ puissance_
Cln:uit de contrOle
Vthu'
.:i'3d
~~a..-----.:.:::..::...j Calcul de Vth par
Vthv' TransfonnaHon ,.. transformation ae ParX et
LHf+------oi~ inve~e de VIh d el
~--+_O«:~....:v.:::.
t3 ~q filtre passe-bas
Vthw VIh q
'" d*
vt3
Vthu'
'"tN~d Calcul de VIh par
~~a.._----.:V:::.:
Vthv' TransfonnaHon
'--I---~:....j inv ~e de Vth d et
-1N q transfonnation de Pari< et
~--+-f~~"":"::':~ filtre passe-bas
Vthw VIh q
Figure 4.10: Schéma fonctionnel du GHT
75
La figure 4.11 présente le circuit équivalent monophasé simplifié pour
la fréquence fondamentale et pour un harmonique h. L'équation 4.2
montre le principe de fonctionnement simplifié du GHT.
Impédance du Z1
transformateu~Zs r • • • • • • • • • • • •,
~-r'~ ~~'~~------~~~
1 ~~V
............,
1 r.
~~
Ish Zs
Courant harmonique
produit par l'EVV
Filtre série
accordé à 60 Hz
Vs
EVV
Tension du
réseau
Filtre accordé à
l'harmonique h
Tension
harmonique
de ligne Vth Vth=lsh*Zsh
(Ih EVV - Ih s.elect.)

Ih
Ih Ih
Zsh T Sélect. EVV
Im~ce totale de la
source ZsT = ZS+Z1
\
Courant harmonique
produit la source
pour l'harmonique h électronique
Figure 4.11: Circuit équivalent du GHT pour un harmonique.
76
50
Vt(t) = ~ sin( mIt) + I(lchnsin(nuv + <PCh n
) - Isensin(nmlt + <Pse.))*IZsnl (4.2)
n=2
Le réglage des courants harmoniques injectés sur la ligne est

complètement contrôlé par la source de tension et cela avec environ
25 % des kVA requis par l'EVV. Ainsi une source de 150 kVA possède
suffisamment de puissance pour tester des EVV jusqu'à 600 kV A.
4.3.2.2 Réglage de la tension harmonique de ligne
Le réglage de la tension de sortie de la source électronique est fait en

fonction des harmoniques de tension à l'entrée de l'EVV. On utilise la
théorie des puissances imaginaires instantanées[6] comme instrument
de calcul des tensions harmoniques et pour le calcul de déphasage des
différents signaux de tension. Cette technique de calcul est rapide et
peut facilement être adaptée à une boucle d'asservissement de type
Proportionnel et Intégral (PI).
Nous avons sélectionné un contrôleur qui permet d'imposer un poids

élevé sur la consigne dans le calcul de son signal de sortie. De plus une
fenêtre de saturation du signal de sortie a été ajoutée afin d'éviter des
dépassements excessifs de la tension de sortie. Les détails de ce circuit
de réglage sont abordés à la section 4.5.
77
4.3.2.3 Utilité des filtres passifs accordés
Vus du réseau de distribution, les filtres accordés aux fréquences

harmoniques reliés en série avec la source électronique sont perçus
comme une impédance élevée. Cela a l'avantage de limiter le passage
du courant de fréquence fondamentale dans la source électronique.
Une source électronique de 140 volts ligne-neutre est alors suffisante
pour produire le nIveau de tension harmonique que l'on doit
superposer à l'onde de fréquence fondamentale. De plus ces filtres ont
une impédance capacitive ayant comme effet bénéfique non
négligeable de corriger le facteur de puissance de l'installation.
4.4 CONCEPTION ET ANALYSE DU GIIT DANS LE DOMAINE FRÉOUENCIEL
4.4.1 Conception des filtres accordés utilisés par le GHT
Nous avons utilisé une version modifiée du programme de conception

de filtres passifs accordés qui a été présentée au chapitre 2. Un filtre
accordé à 60 Hz a été ajouté en série avec l'impédance de la source Zs
du transformateur afin de bloquer le passage des harmoniques de
courant qui seront déviés vers le GHT VIa les filtres parallèles
accordés. Les critères de conception de ces filtres sont présentés à la
section 4.4.1.1 suivante et les résultats de calcul sont présentés à la
section 4.4.1.2.
78
4.4.1.1 Critères de conception des filtres du GHT
Puissance du transformateur (Sn): 2500 kVA
Impédance du transformateur (Zs) : 0.1 pu
Filtre en série avec Zs: C1=JOOO pP L1=7.04 mH, Q1=30
Tension de phase (V p): 575/ ~ Volts

Courant de charge (IL): 600 Ampères
Fréquence fondamentale (fI): 60 Hertz
Facteur de puissance (FP): 0.9
Distorsion maximum sur VL (DHTvd 3 % @ 600 Amp.
Le courant de charge est représenté par le contenu harmonique

typique d'un convertisseur de courant de 6 impulsions.
4.4.1.2 Résultats du calcul des filtres du GHT
Filtre accordé au Se harmonique: CS=780 JlF, LS=0.36 rnH, RS=0.03 .Q
Filtre accordé au 7e harmonique: C7=400 JlF, L7=0.36 rnH, R7=0.05 .Q
Filtre accordé au Ile harmonique: Cl1=200 JlF, Ll1=0.29 rnH, Rl1=0.06 .Q
Filtre accordé au 13e harmonique: C13=300 JlF, L13=0.09 rnH, R13=0.03 .Q
Filtre passe-haut: C ph=260 JlF, Lph = 0.094 rnH, R ph=0.8 .Q
Distorsion harmonique totale:

- sur la tension ligne-neutre: DHTIL =2.9 %,
- sur le courant de la source: DHTV L _N = 0.2 %,
Puissance réactive totale installée: QTfiitre = 174 kVAR
Filtre accordé à la fondamentale:
QTfiltre = 3* XC I *IL2 =425 kV AR
79
4.4.2 Analyse des performances du GHT dans le domaine fréquenciel
La figure 4.12 présente l'amplitude des impédances du GHT (définies

par les équations 4.13 à 4.16) en fonction de la fréquence.
L'impédance du système, aux fréquences harmoniques
caractéristiques d'un redresseur triphasé, est inférieure à celle du
transformateur de distribution utilisé seul. On observe toutefois une
résonnance parallèle très prononcée à une fréquence de 105 Hz. Cette
fréquence ne fait pas partie de la zone d'intérêt du GHT et sera évitée
afin de ne pas créer de surtension de ligne. Des résonnances parallèles
apparaissent également autour du 6e harmonique et entre le 8e et ge
harmonique. Ces résonnances parallèles ne causent pas de réel
problème puisqu'elles ne sont pas à des fréquences caractéristiques du
redresseur et ne sont pas très prononcées.
. . . . . .'. . . .'.. . . . . .'.. . . . ... . . . . . . . . . ~ :. ; ~ ~ :' :~ :' ~~ : ' :~ :' ....
r
. . . . . . . . . . . . . .'.
,
....... ~
,
..... : '. :
!:: ~::: 1::: ~ 1: 11: U1U!.[l:_:~!- :~-H 1hl [ë:-0:rTDln-:

101
0E __ _
~100
~ . ...................... , . .............. ',' .... , .. . .:.-' -
. ~~~~~~~Li~
.. . . ........ , .... .
,
~ 1::: UUW-Wl~L~\~Jtlt.~: ~u~~llW~tœ n:!: : : 1- ::1::: 1- -!: ~::: ~ 111::

1
j 10-
1
:V _!!
lmpedance du filtre accotde a Ia:frequence fondamentale :
!::::::: 11;:: !:!·~·!·~ ~~ ~·~ l~=~::::g::~:~::M~~r~~~fFri~1~~1~~~I?~~ :; l:::;:
-2 - - : - - - - -
10
10'
Figure 4.12: Module de l'impédance des composantes du GHT
80
(4.3)
(4.4 )
1
Zh(f)=~l~--~l----~l----~l----~l-- ( 4.5)
--+ + + +----
Zfs (f) Z1; (f) Zft 1(f) Zft 3 (f) Zfp_h (f)
1
ZCHT(f) = 1 1 (4.6)
-----+---
ZI (f) + Zs(f) ZfT (f)
Zs(f): est l'impédance du transformateur de distribution en fonction

de la fréquence.
Zl(f) est l'impédance du filtre accordé à 60 Hz en série avec

l'impédance du transformateur.
ZfT(f) est l'impédance totale des filtres passifs accordés aux Se, 7e,
11e,13e harmonique et passe-haut
ZGHT est l'impédance totale du GHT vue de la charge.
Z,, (f) = RII + j(27if L n - 2":C

"'./ n
J Zp_h(f ) = -_--;-l~+----;--- + -j27ifC
---
p _h
Rp _ 1t j27ifL />_h
Rn est la résistance en série avec le filtre Ln-Cn pour n=5, 7, Il, 13.
Les équations 4.3 à 4.8 font référence à la figure 4.13.
81
Comme nous l'avons déjà mentionné le GHT a pour objet de faire la
sommation de la tension provenant du réseau de distribution (Vs) et
la tension d'une source électronique (Ve) afin de produire une tension
de ligne (Vt) dont le contenu harmonique est programmable.
La figure 4.13 présente le circuit équivalent du GHT pour la fréquence

fondamentale et pour les fréquences harmoniques. Les rapports V s/Vt
et V s/Ve présentées aux équations 4.7 et 4.8 sont déduites de ce
circuit. À la fréquence fondamentale et pour sa charge nominale l'EVV
est remplacé par une impédance Zcharge=0.52+jO .19 ohms, soit une
charge triphasée consommant une puissance de 600 kV A avec un
facteur de puissance de 0.9. Pour les harmoniques, l'EVV est considéré
comme une source de courant avec une très grande impédance
d'entrée.
Impédance du Zs Zl
transformateu~ ,. .... __ ..... .,
...........................
r-T:~ ~~'4-~--~ __-+~
1
Filtre série
.... .... .. ....
:
accordé à 60 Hz :
Vs
/
Tension du
réseau
IV
/T!
Filtre accordé à .
Vr EVV
l'harmonique h
Figure 4.13: Circuit équivalent monophasé du GHT
82
1 1 1 1
V (60) = Y:(60) * ZfT(60) + Zcharge
[ Il
T ( 4.7)
ZST(60) + 1 1
ZfT(60) + Zcharge
v, Cf) - v Cf) ZST (f) (4.8)

T - e ZsTCf) + ZfT(f)
où: ZST(f) est l'impédance de source totale ZI (f)+Zs(f) (eq. 4.3 et 4.4).
ZfT(f) est l'impédance totale des filtres en fonction de f (eq. 4.5) .
Zcharge = 0.52+jO.19 ohms
À l'équation 4.7 si (Zf(60) Il Zcharge) » ZST alors: Vr =:: 1 et

~
V,
à l'équation 4.8 si ZST(f) » Zh(f) alors: ; =:: 1;
e
qui sont les caractéristiques idéales recherchées pour le GHT.
La figure 4.15 présente les rapports Vt( 60)/V s( 60) et Vt(f)/Ve (f)
obtenues avec le choix de filtres présentés au paragraphe 4.4.1.2.
Notons que le transformateur de distribution a une fréquence fixe et
que la source électronique fonctionne à une fréquence variable.
83
Le rapport Vt(60)/Vs(60) a un module de 0.97 avec un angle de phase
de -0.4 0 • Ceci corespond à une chute de tension de 3% avec un
déphasage nul. La fonction de transfert Vt(f)/Ve(f) varie de 0.92 à
1.08 dans la plage de fréquences de la source électronique avec des
variations du déphasage compris entre 0 et 10 o. Si on examine cette
fonction de transfert pour les fréquences harmoniques d'intérêt (5, 7,
11, 13, 17, 19, ... ) on note des variations de moins de 2 % autour de la
valeur unité et de moins de IOde déphasage.
: : : : ~ : : ~ : ~ : -: ~ :. :- ! ~ : : : : ~ : : ! : : : ': ~ ~ ~ : r'\ :
~150 ..
g>
• •
!:
::
• r • •
:!:::: ~ : .:
:!::!:
~ • .. • "
~ : .:
.. " " J \
~ ~ ~
~ :. :-
!
!
~
~
• • • •
: : : : : : :
\ • • \
: : : : ! : : ! : ! :.: !
~
• 1 • '.
: : :.: !
1 1
~
~
1
: i
! !
1
-
e.
~
....... .
100 .. .: .. :.:.::::::., . . : . . :.: . :.:: ::
> . . : ::: ::: , . ,: ::::
15 . . . ...... , . ......
~ 50 . . .. : .. : ,:.;: :::: .. . . : .. : . : . :. : :::
. . . : ' . : . : . ~ : : : : : . .. . G VtNs : . :. : : : :
0 : ~ : ~: : ,::: : : : :
l l l ~ i l l il; l ~ ; r. ; i ~ l ~ l 11 YVY~ i 1:1 i ; ; i g> .. .
Cl)
>-1'
: : : : ; : : ;: : ~ : .: ~ :. :. ;: ~ : : : : , : : ; : j :.: ; ~ ;: ;
e( <;I angle(Vt1Vs) : :
.;. angle(Vti'Ve) : :
1O-2'--_ _ _ _ ~ _ _ _ _ _...J
-50'------~----_-.1
1~ 1~ li 1~ 1~ li
Rang de l'harmonique Rang de l'harmonique
1.1 r--r---,---------..., 10
a;
~ 1.05 . !
~ .. .. ..,.. " ,. , ..
. ..... . ' . . ,
Qi 1 · .•. . , .•. ' .• l ' '

i" 5 ... ,
, ...'.', .'
'"
~
0
~
"3 . . ~Qi"
~0.95 . " : " : ' : ' : . : :: ... ' GVtNs : . :. :::: Oang le(Vti'Vs) : :
~ , . . , , , ..
:" VtNe : : :::: e(
7 angle(Vt1Ve) : :
0 . 9L...J....----~----_...J
1~ 1~ 10° 10
1
Rang de l'harmonique Rang de l'harmonique
Figure 4.14: Fonction de transfert de la source électronique et du

transformateur de distribution en fonction de f.
84
Ces résultats démontrent que le GHT se comporte réellement comme
deux sources de tension indépendantes produisant une tension de
ligne VT(t)=Vs(t)+Ve(t) si on assume un fonctionnement à 60 Hz pour
Vs et à plus de 180 Hz pour Ve. Une simulation de ce dispositif, dans

le domaine temporel, est toutefois nécessaire afin de valider ces
résultats car la charge (l'EVV) est non linéaire et l'analyse
fréquencielle peut être inadéquate dans ces conditions.
4.5 SIMULATION NUMÉRIQUE DU GHT DANS LE DOMAINE TEMPOREL
4.5.1 Simulation du GHT appliqué aux tests d'immunité
Le GHT a pour objet principal d'effectuer les tests d'immunité des EVV
décrits dans le document SC-610 d'EPRI. On présente, dans cette
section, l'analyse des performances du GHT lorsqu'il est utilisé pour
reproduire chacun de ces tests. Dans chacun des exemples de tests
nous comparerons son comportement en boucle ouverte avec celui en
mode de régulation de type proportionnel et intégral (PI). Le
simulateur utilisé est présenté à l'annexe D
4.5.2 Simu1ation des tests d'immunité aux harmoniques de tension
Ces tests sont décrits dans le document SC-610 comme les tests
d'immunité de premier niveau (les plus importants) identifiés par les
abréviations 13, 14 et 15. Le document suggère un test où le client
spécifie l'amplitude et la phase de différentes combinaisons
d'harmoniques de tension procurant une distorsion harmonique totale
85
(DHT) de 10%. Le tableau 4.1 présente les combinaisons d'harmoniques
qui ont été choisies afin de vérifier le comportement du GHT.
Tableau 4.1: Combinaisons d'harmoniques utilisées pour les tests

d'immunité [1]
Combinaison # 5e 7e Il e 13e 17e DHT
Harm. 1 33 V, 0° OV OV OV OV 10 %
Harm . 2 OV 33 V,180° OV OV OV 10 %
Harm. 3 OV OV 33 V,180° OV OV 10 %
Harm. 4 OV OV OV 33 V, 0° OV 10%
Harm. 5 OV OV OV OV 33 V, 0° 10 %
Harm. 6 66 V, 0° 47 V,180° 30 V,180° 25 V, 0° OV 27%
Les résultats de ces tests sont présentés aux figures 4.15 à 4.21 sous la
forme d'une analyse fréquencielle des signaux de tension obtenus
pour le test en boucle ouverte et en mode de régulation, comparés aux
valeurs désirées. On présente également les formes d'ondes désirées et
obtenues dans chacun des cas et les formes d'ondes de tension et
courant du système pour la combinaison #6 avec comme charge un
EVV de 600 kV A. Les résultats de ces simulations ainsi que d'autres
combinaisons d'harmoniques sont présentés en détail à l'annexe F. Le
tableau 4.2 résume les résultats obtenus pour différentes
combinaisons mises à l'essai.
86
Tableau 4.2: Résultats de simulation du GHT pour des exemples de
génération d'harmoniques de tension
Combinaison # VT* VT bo Err VTbo VTbf Err VTbf
Harm. 1 (HS) 33 V, 0° 33.9 V, 5° 2.5%, 5° 32.4 V,9° -1.7 %, 9°
Harm. 2 (H7) 33 V, 180° 33.5 V, 194 0 1.5 %,14° 32.7 V, 181 0 -0.9%, 1°
Harm. 3 (Hl1) 33 V, 180° 31.8 V, 198 0 -3.6%, 18° 32.8 V, 201 0 -0 .6 %, 21°
Harm. 4 (H13) 33 V, 0° 32.1 V, 18° -2.7 %, 18° 33.0 V, 6° 0.0 %, 6°
Harm. 5 (H17) 33 V, 0° 29.5 V,17° -11 %, 17° 32.9 V,32° -0.3 %, 32°
Harm. 6 (HS + 66 V, 0° 64.5 V, 4° 2.3 %,4° 66.4 V, 8° 0.6 %, 8°
H7+ 47 V, 180° 43 .2 V, 186 0 8.1 %, 6° 46.9 V, 181 0 0.2 %, 1°
Hll+ 30 V, 180° 28 .5 V, 191 0 5.0 %, 11° 29.9 V, 197 0 0.3 %, 17°
H 13) 25 V, 0° 22.2 V, 12° 11.2%, 12° 24.5 V, 3° 2.0 %, 3°
Où: VT*: Valeur désirée de tension de ligne de l'hamonique h.
VT bo: Valeur de tension de ligne de l'hamonique h obtenue

avec le test en boucle ouverte.
ErrVTbo: Erreur sur l'amplitude et la phase de l'harmonique que

l'on désire Imposer en boucle ouverte.
VTbf: Valeur de tension de ligne de l'harmonique h obtenue

avec le test en boucle fermée (mode de régulation
proportionnel et intégral).
Err VTbf: Erreur sur l'amplitude et la phase de l'harmonique que

l'on désire Imposer en boucle fermée.
87
Test en boucle ouverte et en mode de regulation
~r-------~--------'---------'--------r--------.--------.-----.
.~ -200
'c:"
CI) Valeur desiree : h 5=33.0 vOlts / 0.0 degres
f- - 400 Valèùr 'obièriuii:éri ooùclè oùVèri~ .. .. . . :' ... h'5=3!nhi~Hs "FM ilégres' ':' .. .. . .
_ ~L-
Valeur
______obtenue :en mode
~~
de reguJation
______ ~
h 5=33.0 -VQ1ts
________L-______ L____
/ 3.4 degres
__ ~ ~ ______ ~ ____ ~
o 0.005 0.01 0.015 0.02 0.025 0.03

Temps (sec)
Figure 4.15: Génération d'harmoniques de tension: combinaison #1

~'--------r--------.--------''--------'--------'--------'-----'
_ ~L-
Valeur
______obtenue 'en mode de reguJation
- L________ ~ ______ ~~ ______- L______ ~~ ______ ~ ____ ~
o 0.005 0.01 0.03

~.--------.--------.---------,--------.--------.--------.-----.
h11 =33.0 v01ts /180.0 degreS

. . :' . .. h11 ';3f.8 v~fts 1':'1Elf.5 oegres' ..
_ ~L- ______- L______ ~~
h11 =32.8 v Ils /-158.9 degrès ______
________L-______- L______ ~~ ~ ____ ~
o 0.005 0.01 0.Q15 0 .02 0.025 0.03

Temps (sec)
88
~r--------.--------'---------.--------'--------'--------'-----'
h13=33.0 volts 10.0 degres

. :' . .. h1~32 . 1 V9rts l1é.O ·deiiès · . ....
_~L- ______ ~~ ______ ~
h13=33.0-L
________L -______ volts /5.9 degres
____ __ ~ ~ ______ ~ ____ ~
o 0.005 0.01 0.Q15 0.02 0.025 0.03

Temps (sec)

~r---------.---------.---------.---------'---------'---------'-----~
0.005 0.Q1 0.Q15 0.02 0.025 0.03

Temps (sec)
Figure 4.19 : Génération d'harmoniques de tension: combinaison #5

~r--------.--------'---------.--------'--------'--------'-----'
Valeur desiree .
_~L- ______ ~ ______ ~~ ________L -______-L________ ~ ______ ~ ____ ~
o 0.005 0.01 0.Q15 0.02 0.025 0.03

Temps (sec)
89
La figure 4.21 présente le circuit de réglage du GHT pour un
harmonique de tension. Le réglage de l'amplitude et de la phase de
l'harmonique nécessite une transformation de coordonnées polaires
(module et phase) vers un système de coordonnées rectangulaires (axe
direct et en quadrature de la tension selon une transformation de
Park). C'est pourquoi nous présenterons les résultats de simulation
sous ces deux formes.
La figure 4.22 présente la dynamique de ce circuit de régulation pour

la génération d'un harmonique de tension (combinaison #1 de la figure
4.15). Ces caractéristiques dynamiques dépendent essentiellement des
paramètres choisis pour les deux régulateurs de type PI (régulateurs
utilisés pour régler les composantes désirées, directe et en quadrature,
de la tension harmonique). Elles dépendent également des
caractéristiques du filtre passe-bas utilisé (filtre Butterworth du 4e
ordre dont la fréquence de coupure est de 20 Hz) afin de séparer
l'harmonique désiré des autres harmoniques et de la fondamentale.
Pour cet exemple on observe un dépassement de 17% sur la consigne

d'amplitude et cela pour une durée de 0.03 seconde. De plus, il faut
0.3 seconde au circuit afin de réduire à zéro l'erreur sur la consigne
d'amplitude. Si un tel dépassement s'avérait inacceptable pour
certains essais, il pourrait facilement être évité en modifiant les
paramètres des régulateurs PI qui comprennent chacun un circuit
d'antidépassement (voir figure 4.33) qui n'était pas actif pour cet
exemple.
90
Valeure meeurees de Vd etVq
Fi~re Anal. Butt.
o alpha
Tenson obota
V llipha, bole Ca/rut do Vd 01 Vq Transformation inverae
de VdA et VqA -> Vu .. , Vv- , Vw-
L---+---.la~=~~~==~--
l Polaire a Ca_on l ___ -l
Conaigne en COOfdMnMII Consl""e tranafonnee en eoordomeea

Polaire: Mo(h~e et phase Rectangulaire: Vd ot Vq
Figure 4.21: Circuit de réglage d'un harmonique de tension

ëcAO
:;35 ......
rL .
't-'-'-=-.;-:..:'::::'='=':;=: ='~'-''!-:-' ---1
~ 5r-~--'-----'------r-----r-----'
0-
1--
> 4
.: .. . ...... .... .. .: .. .. . .
j
: -_
. =
' • • ••• 1 ••
~
1-
>30 ....
~ ~ 3
f
Ir '
>25
~
~
::l
'§ 2
, ,.
,
1
i 20
ï5 15 .... . . ~ .. .... ~ .. ...

~
::l
0-
••• '\ ' • 1 •
.cil c:
Q)
œlO Valeur de consigne sur Je Se har nonique

.cil
~5 . 1\ ...
0 1 \
. . Valeur mesuree du 5e·harmoniqb e œ0
8. "
o O u---~'~----~----~~----~--~ E
8 -10
o 0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5
40r----,,-~~----_.----_.----~ 100r----,,---~----_,----_,--__.
c-
35 ......... L 50 .,., .
,-
~30 'ê' 1
~ 0~1~
/ ~~~~~~--~----~--~
~25
~20 ~ -50 ,.. . . .. ; ... ... ;.
~ 15
'5
~
Q)
, • •
g>-100 ... .. . : .. . .. . : . ... . . :.. .. . . ....... .

_. -. -
~ 10 "\
..... ,.. . . ...
Valeur de consigne sur Je Se har flonique
<:
- 150 . . . . . . . . . . . . . . . . . . . . .'. . . . . .
J ..
~
5 -,' . \ . Valeur mesurée du 5e· harmoniq~ e
o -200L---~----~----~----~--~
o
o 0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5
Figure 4.22: Réponse dynamique du circuit de réglage pour l'exemple

de la figure 4.15 (h5=33 volts et 0°)
91
La figure 4.23 présente les formes d'onde de tension ligne-neutre et
de courant de ligne lorsque la source électronique n'est pas en action
pour une charge de 600 kV A. Cette charge correspond à la pUIssance
maXImum que le frein moteur du banc d'essai peut absorber et qUI
deviendra ainsi la puissance maximum du banc d'essai pour EVV.
Les courbes des figures 4.15 à 4.20 sont des résultats de simulation
pour un EVV de 500 HP, ce qui correspond à la limite maXImum
actuelle du banc d'essai moteur lorsqu'il est raccordé à un
autotransformateur de 500 kV A. Pour le futur banc d'essai pour EVV,
on remplacera l'autotransformateur par un transformateur de
distribution de 2.5 MV A ce qui permettra de tester des EVV de
600 kVA.
La figure 4.24 présente les formes d'onde de tension et courant pour

l'exemple de la combinaison #6 de production d'harmoniques sur un
EVV de 600 kVA. Dans cet essai on impose une DHT de près de 30 %;
ce qui est trois fois supérieur aux recommandations contenues dans le
document SC-610. Notons également que dans cet exemple les
harmoniques de courant qui doivent être fournies par la source
électronique sont en opposition de phase avec les courants
harmoniques produits par la charge. Les puissances mises en jeu dans
cet exemple correspondent donc au pire cas (puissance maximum
nécessaire, DTH maximum et phase des harmoniques en opposition
avec celles produites par la charge) que l'on pense rencontrer lors des
essais réels.
92
Ci:
E Tension et courant de ligne pour une phase pour Ve=O vott
~ 1()()()
Cl>
c:
g
-8
~::::1
500
i
~
!!!
'Sc:Cl>
-500
~
Tension de ligne = 310 volts
Courant de ligne =625 amperes :
.~-1()()()~------~--------~--------~--------~--------~--------~--~
~ 0
0.005 0.01 0.D15
Temps (sec)
0.02 0.025 0.03
Puissance apparente triphasee tiree par la charge
600
[500
Cl>
-400
~
<CI
~300
<CI
... . ................... . .' ..
' .. ", .. , . . ....
' , ....... ,' . . .... . ..... ' ......... . . . ..... .
',
~ . . .. . .....
C:200 ..................
. . ......
. . .... . . ... . ...... . ... , .. . ... ',' , ,'
.Œ
::::1
c.. 100 s ~ Y.Tu~ilu .+. VTv·jLv. + Iflw··iLW ... . ....... , .... .
Puissance moyenne totale: 555.7 kW
0.005 0.01 0.D15 0.02 0.025 0.03

Temps (sec)
Ci: Tension et courant de la source electronique pour Ve=O volt

~ 600~------~----~~----~------~------~------~------~
=
!!!
'5-200
~
cl, -400
c:
~-600 .. . ip!l .d.e. llg.ri!'l:-(l~W~ ~ !é\ ~lJrç~ ~e.ctr9(li~ . ..=. q Y. . ....... . ....... ,' ...... .
c: Courant de la source electronique =Courant totàle des fittres passifs =285 amperes
.2~_600L- ________ ________-L________
~ ~ __________L -________ ~ ________-L________ ~
~ 0 0.005 0.01 0.D15 0.02 0.025 0.03 0.035

Temps (sec)
Figure 4.23: Formes d'onde pour Ve o volt pour une charge de

600 kVA
93
0:
E Tension et courant de ligne pour l'essai #6
~
CIl
c:
g
~
~
i '5
Il! -200
CIl
........ .
c 400 ... . ...... .' . . . .. .

cI, -
c:
0>-600 "T:e.n~!<>fl .qe. l!g~.=: ~?~ Y<?I!s.
'8 Courant de lig~ =572 ampe
._ _~~____~~____~~~~~~~____~~______-L~~~~____~
~ 0 0 .005 0 .01 0.D15 0 .02 0 .025 0.03
~ ~~~
Puissance apparente triphase pour l'essai #6
0.D15 0 .02 0 .025 0 .03

Temps (sec)
Ci: Tension et courant de la source electronique pour l'essai #6

E
~
=
ë
~
:::l
0
~
""~0
~
Il!
'S
~
cl,
c:
g
0 0 .005 0 .01 0.015 0 .02 0 .025 0.03

Temps (sec)
Figure 4.24: Forme d'onde de courant et tension du système pour

l'essai #6 de génération d'harmoniques de tension.
94
Les résultats présentés au tableau 4.2 confirment que le
fonctionnement du système s'approche des caractéristiques idéales
recherchées. L'erreur sur l'amplitude des harmoniques, en boucle
ouverte, est inférieure à 4 % jusqu'à la 13e harmonique et augmente à
Il % pour les harmoniques de rang supérieur à 17. Le système de
correction en boucle fermée, avec un réglage de type proportionnel et
intégral, diminue l'erreur sur l'amplitude des harmoniques de tension
imposées à moins de 2% sur toute la plage de fréquence utile de 180 à
3000 Hz (Voir l'annexe F pour le détail des résultats).
L'erreur sur le déphasage des harmoniques dépend principalement de

la précision de la synchronisation entre l'onde de tension sortant du
transformateur de distribution et le signal envoyé à la source
électronique. Une erreur de déphasage temporel se traduit par une
erreur de phase amplifiée par le rang de l'harmonique qui est imposé.
Ce problème est aussi valable pour le circuit de réglage avec
rétroaction car dans les deux cas nous avons noté une erreur de phase
de 5 à 20 0 jusqu'à la 13e harmonique et de 10 à 32 0 à partir de la
17 e harmonique.
Le courant efficace traversant la source électronique est de

370 ampères avec des crêtes de 1000 ampères. La tension efficace de
la source est de 99 volts avec des crêtes à 200 volts. Si l'on tient
compte que celle-ci peut fournir des crêtes de courant de 2.5 fois son
courant nominal, la puissance minimum requise pour cet exemple sera
de 160 kV A @ 140 volts l-n et 380 ampères.
95
4.5.3 Simulation des tests d'immunité aux encoches dans la forme de
l'onde de tension (ttResponse to Voltage-Notching tt )
Ce test de deuxième niveau (test souhaitable mais non reconnu comme

test de base indispensable) est décrit dans le document SC-610 comme
le test d'immunité 19. Le document suggère d'imposer des encoches de
tension d'une largeur de 10 à 30 degrés centrées autour de 60 degrés
de phase et dont la profondeur est augmentée progressivement
jusqu'à ce qu'il y ait défaillance de l'EVV sous test. Le tableau 4.3
présente les combinaisons qui ont été choisies afin de vérifier le
comportement du ORT. Les résultats de ces tests sont présentés aux
figures 4.24 à 4.30 sous la forme de courbes comparatives, dans le
domaine temporel, entre les valeurs désirées et les valeurs obtenues.
Tableau 4.3: Combinaisons d'encoches de tension utilisées pour les

tests d'immunité 19.
Largeur de Profondeur de Debut de l'encoche

Combinaison #
l'encoche (deg) l'encoche (volt) (deg )
Enc 1 10 ° 100 V 50 °
Enc 2 10 0 200 V 50 0
Enc 3 10 0 300 V 50 0
Enc 4 20 ° 100 V 40 0
Enc 5 20° 200 V 40 0
Enc 6 20° 300 V 40 °
96
Test en boucle ouverte: profondeur de 100 volts et largeur de 10 degres
500 .. ". " ..... . ... " ,"
l~
'S
2!
d,
c:
J!
.~<lJ
c:
Valeur desiree
~
Valeur obtenu~ par simulation :nume
-500
0 100 200 300 400 500 600 700
Angle (deg)
Figure 4.25: Génération d'encoche de tension: essaI #1
500
~
~
'S
CI)
c:
d,
c:
g
.~<lJ
!
f-
Valeur desiree
-500
0 100 200 300 400 500 600 700
Angle (deg)
500
=-~
~
'S
2!
g
~
c:
0
.~
Valeurdesire6
~
-500
0 100 200 300 400 500 600 700
Angle (deg)
97
500
~
~
~
'5
2!
cl,
c::
g
c::
.2
<Il
c:: Valeur desiroo
~
Valeur obtenu!! par simulation :nume
-500
0 100 200 300 400 500 600 700
Angle (deg)
500
~
~
~
OS
Q)
c::
1
Q)
c::
g
8
"Cn
!ii
1-
Valeur desiroo
-500
0 100 200 300 400 500 600 700
Angle (deg)

500 ... . .. .
-5OOL-______~~------~~------~~------~~·-·-·-·-··-·-·-·-·L
·-··-·-·-·----~________~
o 100 200 300 400 500 600 700
Angle (deg)
98
La figure 4.31 présente les formes d'ondes de tension et courant pour
l'essai #6 de production d'une encoche dans l'onde de tension pour une
charge de 350 kVA. Le courant efficace de la source électronique est
dans ce cas de 274 ampères avec des crêtes de 500 ampères et la
tension crête de la source électronique est de 190 volts . La puissance
mIll1mUm requIse pour cet exemple sera donc de 126 kV A
@ 140 volts ligne-neutre et 300 ampères .
Ci:
E Tension et courant de ligne pour l'essai #6
~
Q)
~,--------,-------,--------.--------.--------~-------.----~
c:
g 600 ............: ....... . . .. ....... ... . . . . . ... . ....... , ... .
-8
~400
§200
~
2! -200
~
~ -400
Cl-600 re.n~l<>fI.d.e. I!9.~. ,,: ~?~ y9l!s. . . . . . . . . . . .... . ..... . . ... ........... .. ... .
'8 Courant de ligne = 327 am peres :
._-~~------~------~--------~------~--------~------~----~
~ 0 0.005 0.01 0.Q15 0.02 0.025 0.03
Temps (sec)
Ci: Tension el courant de la source electronique pour l'essai #6

i 600.--------.-------,,--------,-------,--------,--------,-----,
~400
o
Q
c 200
ï
~
c:
.2
~ -600
~ 0 0.005 0.01 0.015 0.02 0.025 0.03
Temps (sec)
Figure 4.31: Forme d'onde de courant et tension du système pour

l'essai #6 de génération d'une encoche de tension.
99
4.5.4 Simulation des tests d'immunité aux perturbations causées par la
commutation de condensateurs ("Response to transients from
capacitor switching")
Ce test de deuxième niveau est décrit dans le document SC-610

comme le test d'immunité 18. Le document est incomplet en ce qui
concerne la définition de la forme d'onde de tension à imposer. Nous
avons alors choisi de superposer à l'onde de tension de ligne une
sinusoïde de lIe harmonique atténué sur une période de 33 msec car
la forme d'onde résultante est similaire à celle présentée dans le
document SC-610 . Les résultats de simulation obtenus sont présentés
à la figure 4.32.
Test de commutation de condensateur

~.--.-------.------.------.-------.------,------.------~
~ 400 ... . . : . . . . . . . .
~
~200
!i
ê
~
g
. .
Valeur obtenue par simula~on numerique :
-~
0.08 0.09 0.1 0.11 0.12 0.13 0.14
Temps (sec)
Figure 4.32: Exemple de forme d'onde de tension produite pour les

tests d'immunité aux commutations de condensateurs
100
4.6 MODE DE RÉGLAGE DES HARMONIQUES DE TENSION
Il a été démontré aux sections 4.5 et 4.6 que le GHT agit comme une
source de tension harmonique. On pourrait régler l'amplitude et la
phase des harmoniques de tension en boucle ouverte en utilisant un
simple générateur d'onde arbitraire comme signal de commande
toutefois cette façon de faire comporte certains inconvénients:
• L'erreur sur le réglage de l'amplitude et de la phase des

harmoniques d'ordre supérieur à 17 est importante (>10%);
• la réponse des filtres au régime transitoire peut entraîner des

dépassements si la consigne est appliquée trop brusquement;
• ce système requiert un calculateur ou un circuit électronique

pour synchroniser le signal du générateur avec la tension
fournie par le transformateur.
Ainsi il est recommandé d'utiliser une boucle de régulation pour

chacun des harmoniques qui doivent être réglés. Nous proposons de
mesurer ces harmoniques en transformant l'onde de tension triphasé
dans un système de coordonné biphasé selon les équations 4.9 à 4.12.
On obtient alors une composante directe et en quadrature de VT dont
les parties continues sont les valeurs réelles et imaginaires de l'onde
de tension harmonique de V T. On utilise un filtre passe-bas afin
d'éliminer la partie alternative du signal de mesure et on compare ces
valeurs mesurées avec les valeurs cibles VT*d et VT*q. Les valeurs
cibles proviennent d'un changement de coordonné polaire (Amplitude*
101
et phase*) à un système de coordonné rectangulaire (VT réel * ,
VTimag.*).
v: réf. ] -fi2 [ sin(h * * t) ] (O.
[Vv"ririff.. = 3 sin(h
w
sin(h * * 21rj3)
(O.
* * + 21rj3)
(O.
t-
t
(4.9)
~ réf.] f2[1 ~ --J3~ ][V:ri

[ ~réf.
f
= ~3 0 -J3 v"ri/. .] (4.10)

2 -2 Vwri/.
1
2 (4.11)
-J3
2
d] = [: ri f. ~rif.][~] (4.12)
rq ~ri/. ~rif. ~
où: h est le rang de l'harmonique à mesurer,

roI est la pulsation fondamentale de l'onde de tension
(ro 1=21t*60),
"
d est la composante continue de d,
q" est la composante continue de q,
W 1 est la pulsation fondamentale de l'onde de tension (21t*60),

VT u ,v ,w
L'équation (4.13) présente la transformation inverse de d et q en une

tension harmonique de ligne après avoir filtré leur composante
102
alternative avec un filtre passe-bas ayant une fréquence de coupure
de 10Hz pour obtenir J et q.
o
~ [~réf. ~réf.]-I[~] (4 . 13 )
2
~
-Vqréf· Vdréf. q
2 2
Étant donné que le système à régler possède une fonction de transfert

proche de l'unité nous avons utilisé une circuit de réglage qui
comporte un poids important (-0.8-0.9) sur la consigne ainsi qu'un
circuit de limitation sur le signal de sortie du contrôleur. Ceci permet
d'obtenir un réglage rapide et évite les dépassements excessifs qui
pourraient survenir en régime transitoire. Les valeurs de chaque
harmonique sont règlées individuellement par deux contrôleurs PI de
ce type . Le schéma fontionnel du régulateur est présenté à la
figure 4.33. et sa fonction de transfert est donnée par les équations
4.14 à 4.16.
e = (ypc - y) (4.14)
P = K*(b * ypc - y) (4.15)
dI/dt = (KlTi) * e + lI(T t (v-u)) (4.16)
v = saturation( P + l ) (4.17)
où: e: erreur entre la valeur mesurée "y" et la consIgne "y pc " ,
103
P: est la partie "Proportionnel" du contrôleur
K: est le gain proportionnel du contrôleur
b: est le poids sur la consigne de la composante proportionnelle,
I: est la partie "Intégral" du contrôleur,
Ti: est le temps d'intégration de la partie "Intégral" du

contrôleur,
Tt: est le temps d'intégration du circuit d'anti-dépassement du

contrôleur,
u: est la sortie du contrôleur écrêté à une valeur mlmmum et

maximum de saturation
v: est la sortie du contrôleur
Controleur PI avec circuit d'anti-depassement
Point de consigne ypc.
2r-------,--t====~
Valeur mesure y
Poids sur la consigne
+
Sortie du controleur
Somme P+1 Proportionel
Figure 4.33: Schéma fonctionnel du contrôleur utilisé pour le réglage

des harmoniques de tension.
104
Ce circuit a été utilisé pour le réglage des harmoniques de tension
présenté à la section 4.5.2. Il permet de régler rapidement l'amplitude
de l'harmonique désiré avec une précision acceptable (err<l 0%) après
16 msec. et avec une erreur quasi nulle «2% ) après 0.2 sec et cela sur
toute la plage de fréquence étudiée (300-3000 Hz.)
4.7 CONCLUSION
Nous avons vu dans ce chapitre le pnncipe de fonctionnement, l'utilité

et la conception du GHT. Ses performances ont été établies à partir
d'un modèle simplifié dans le domaine fréquenciel puis vérifié par des
exemples concrets avec des simulations numériques dans le domaine
temporel. Nous avons choisi et calculé les principales parties
constituant le GHT et suggéré l'installation de filtres passifs accordés
série et parallèle. La puissance suggérée de la source électronique est
de 150 kV A @ 140 volts avec un facteur de crête répétitif (cress
factor) de 2.5. La somme de ces composantes constitue un GHT devant
pouvoir tester l'immunité des EVV d'une puissance de 600 kV A et
moins . La faisabilité du GHT a été confirmée par des résultats de
simulation qUI sont fort encourageants tant au niveau des
performances en mode de régulation qu'en boucle ouverte.
La suite des travaux consiste à optimiser le design du GHT. On devrait

pOUVOIr diminuer d'environ de moitié le courant de fréquence
fondamentale pénétrant dans les filtres passifs accordés. Ce résultat
pourait être obtenu soit en diminuant le nombre de kVAR des filtres
105
passifs soit en utilisant le GHT afin de limiter le passage de ces
courants de fréquence fondamentale. Dans le premier cas, la distorsion
harmonique totale, lorsque le GHT ne fonctionne pas, s'en trouverait
augmentée mais il serait facile de compenser cet effet avec le circuit
de réglage du GHT ajusté au point de consigne 0 volts pour toutes les
harmoniques désirés. Dans le deuxième cas, il serait possible d'utiliser
une version modifiée du filtre actif présentée au chapitre 3 afin
d'obtenir cet effet en réglant la valeur de ~q qui minimise le courant
de 60 Hz pénétrant dans les filtres passifs et la source électronique.
Il serait également souhaitable de diminuer le nombre de kV AR du

filtre disposé en série avec le transformateur de distribution de
manière à diminuer le coût du GHT. Nous devrons également étudier
plus en détail, par des simulations numériques, les effets des
fluctuations de la fréquence fondamentale du réseau et des
perturbations sur la lecture des signaux de tension de ligne. Les
travaux en ce sens et sur d'autres détails connexes vont se poursuivre
jusqu'à la mise à l'essai d'une maquette de 200 kVA prévue pour la
fin de l'année 1994.
106
CHAPITRE 5
CONCLUSIONS
Ce mémoire présente la conception d'un banc d'essai permettant de

vérifier l'immunité aux perturbations de la forme d'onde de tension de
ligne. Il s'agit d'un projet industriel d'envergure réparti sur deux ans.
Dans ce cas une étude préliminaire approfondie, par des simulations
numériques du système était nécessaire puisque ce type de banc
d'essai n'a JamaIs été réalisé pour des niveaux de puissance aUSSI
élevée (600 kVA). Cette étude a été réalisée et menée à bien grâce à
une collaboration étroite entre le L TEE et le GREI de l'UQTR.
Dans ce travail, nous avons élaboré la partie du projet se rapportant

aux simulations du banc d'essai appliquées aux tests décrits dans le
document SC-610 d'EPRI. Nous avons également présenté une étude
des filtres actifs et passifs par des simulations numériques et le
développement d'un programme de conception de filtres passifs
accordés. Des modèles de simulation des EVV ont été développés. Ces
modèles qui ont été forts utiles à la conception du GHT, pouront être
utilisés pour les tests de compatibilité des systèmes décrits dans le
document SC-610.
107
Le programme de calcul de filtres qUI a été présenté au chapitre 2 a
comme particularité d'être couplé à un programme de simulation, dans
le domaine temporel, pour obtenir les paramètres électriques de calcul
et afin de vérifier ses performances en régime dynamique. Cette
technique de conception est très avantageuse par rapport au calcul
dans le domaine fréquenciel puisqu'elle permet de vérifier les
performances des filtres en régime transitoire ou déséquilibré sur des
systèmes pouvant être très complexes à étudier par d'autres méthodes
de calcul.
L'étude du filtre actif présenté au chapitre 3 a été utile afin de vérifier

notre technique de simulation dont les résultats ont pu être comparés
à ceux publiés par le concepteur de ce filtre [10]. Nous avons détecté
certaines faiblesses du circuit de réglage qui devront être éliminé si on
décide d'incorporer ce filtre au ORT afin de minimiser les courants de
fréquence fondamentale pénétrant dans la source électronique. Ces
travaux font partie de la suite du projet qui se terminera avec la mIse
en service d'une maquette du banc d'essai prévue pour la fin de
l'année 1994. D'ici là, nous tenterons d'améliorer le design actuel du
ORT de manière à minimiser le coût des équipements et de la mise en
service.
Les performances du ORT ont été établies dans le domaine fréquenciel

pUIS vérifiés par des exemples concrets avec des simulations
numériques dans le domaine temporel. Les simulations ont démontré
la faisabilité du projet ainsi que les limites de fonctionnement du ORT.
108
RÉFÉRENCES
[1]. Alexander Domijan and Ernesto Embriz-Santander, "SYSTEM

COMPATIBILITY SC-610 ADJUSTABLE SPEED DRIVES USED IN LOW-
VOLTAGE AC POWER SYSTEMS", Florida Power Affiliate & Power
Electronics Consortium University of Florida, Gainesville FL 32611,
January. 1994.
[2]. SIMULINK®, Dynamic System Simulation Software, "SIMULINK

User's Guide for the X Window System", MATHWORKS Inc.24 Prime
Park Way, Natrick, Mass., Dec. 1993
[3]. MAT L A B ® , High-Performance Numeric Computation and

Visualization Software, "MATLAB Reference Guide", MATHWORKS
Inc.24 Prime Park Way, Natrick, Mass., Oct. 1992
[4]. A. Ba-Razzouk, K. Debebe, V. Rajagopalan, Conception assistée par

ordinateur des systèmes d'entraînement des machines électriques, à
l'aide du "logiciel SEMALINK", Rapport de recherche, UQTR, décembre
1992.
[5]. A. Ba-Razzouk, A. Pite t, V. Rajagopalan, Technique de conception

assistée par ordinateur des systèmes électroniques de puissance à
l'aide du logiciel SIMULINK (logiciel SIMUSEC), Rapport de recherche,
109
UQTR, mars 1993.SEMALINK, Simulation des système en électronique
de ,
[6]. Hirofumi Akagi, Yoshira Kanazawa, and Akira Nabae,

"Instantaneous Reactive Power Compensators Comprising Switching
Deviees without Energy Storage Components", IEEE Transactions on
Industry Applications, vol. IA-20, no. 3, pp.625-630, May.lJune.
1984, .
[7]. Wilsun W.Xu, "A Practical Harmonie Guideline for Adjustable

Speed Drive Applications", IEEE Transactions on Power Delivery, vol. 7,
no. 1, pp.427-433, January. 1992,.
[8]. IEEE Guide for Harmonie Control and Reactive Compensation for
Statie Power Converters, IEEE Std. 519-1981, New York, 1981.
[9]. HVDC Power Transmission Systems, Technology and System

Interactions, K.R. Padiyar, John Wiley & Sons, ppI45-158,.1990.
[10]. C.E. Lin, M.T. Tsai, Y.S. Shiao and C.L. Huang, "An Active Filter for
Reactive and Harmonie Conpensation Using Voltage Source Inverter",
IEEE International Conference on Advences in Power System Control
Operation and Management, November 1991, Hong Kong, pp.588-593.
110
[11]. Hideaki Fujita and Hirofumi Akagi, "A Practical Approach to
Harmonic Compensation in Power Systems-Series Connection of
Passive and Active Filters", IEEE Transactions on Industry
Applications, vol. 27, no. 6, pp.l020-1025, Nov.lDec. 1991,.
[12]. Ali Moshref, "SYSTEM MODELING FOR HARMONIC STUDIES",

CYME International Inc., CEA-CYME Workshop on Power System
Harmonics, Montréal, 1989.
[13]. Richard W. Y. Cheung,"VARIABLE SPEED DRIVE AND POWER

SYSTEM COMPATIBILITY", Association Canadienne de l'électricité,
Rapport # 9133U860, ONTARIO HYDRO TECHNOLOGIES, Toronto,
Ontario, Décembre 1993.
[14]. S. Lahaie, Z. Yao, V. Rajagopalan, "Guide d'utilisation du logiciel

SIMUSEP", Chaire de recherche industrielle Hydro-Québec-CRSNG,
UQTR, décembre 1994
111
ANNEXE A
MODÈLES DE SIMULATION DES EVV
On présente dans cette annexe une description sommaIre des modèles

de simulation des EVV qui ont été utilisés pour la conception du banc
d'essai. Les modèles ont été développés avec le logiciel de simulation
Simulink version 1.3 sur une station de travail de Hewlett Packard.
Bien que cela n'ait pas été vérifié les fichiers, fournis sur disquette,
devraient être compatibles sur PC pour la même version de Simulink.
On présente aux figures A.l, A.2 et A.3 trois exemples de modèles

d'EVV qui ont été utilisés ou développés pour ce projet. Dans chacun
des cas on peut distinguer six parties différentes:
1 Une source de tension avec son impédance interne
2 Un redresseur à diode triphasé
3 Le circuit de filtre sur le lien à tension continue ou du côté

réseau
4 Un onduleur triphasé
5 Le circuit de commande de l'onduleur qui utilisera une

modulation par largeur d'impulsion (MLI) ou une commande
"six steps"
6 Un modèle dq d'un moteur à induction
112
Il est possible de vaner les paramètres de chacune de ces parties du
circuit et de faire différentes combinaison de configuration de manière
à simuler de nombreux types d'EVV raccordés sur des alimentations
électriques différentes.
Source de tension
~y
Circuit de commande
Onduleur-3PH
Moteur a InducHo~
Figure A.1: Schéma SIMULINK d'un système d'alimentation d'un EVV

avec onduleur "6 steps" dans lequel on utilise un filtre
Lf-Cf.
Il 3
svy 'lw
Sou rce de tension
avec imp. interne
Ponl de diodes
triphas.
Va:
• ed
Tm""
Onduleur-3PH
Irb
Moteur a inductio~
Figure A.2: Schéma SIMULINK d'un EVV à modulation MLI ayant

uniquement un filtre capacitif comme lien en tension
continue.
114
Source de tension ILw
avec imp. interne 1
OndJteur-3PH
ire
Moteur a induction
Figure A.3: Schéma SIMULINK d'un système d'alimentation d'un EVV

à modulation MLI avec un filtre inductif du côté du
réseau.
115
ANNEXE B
FICHIERS DU PROGRAMME DE CALCUL ET DE CONCEPTION DES

FILTRES PASSIFS ACCORDÉS
B.1 PRÉSENTATION DU PROGRAMME DE CALCUL DE FILTRES PASSIFS
On présente dans cette annexe le programme de calcul de filtres passif

accordés. Le programme a été développé avec le logiciel de calcul
Matlab version 4.2 sur une station de travail de Hewlett Packard série
700. Ce programme est fournit sur une disquette compatibles sur PC
pour la même version de Matlab.
L'entrée des données ne se fait pas de façon interactive; on doit éditer

le fichier "fpassif.m" et modifier les paramètres du système aux
endroits qui ont étés idendifiés clairement dans le programme par des
lignes de commentaires avec points d'intérogation "%?????????????,'.
Ces paramètres sont:
1- Les paramètres électrique du système
la puissance du transformateur de distribution (SnomT) en kV A,
l'impédance du transformateur de distribution (Imp) en "pu",
la tension efficace ligne-neutre du réseau (Vp) en volt,
Le courant de ligne à pleine charge (IL) en ampère,
116
La fréquence du réseau fI (60 Hz),
2- Les spécifications désirées
Le facteur de puissance désiré (Fpr) sans unité de 0 à 1,
La distortion harmonique maXImum sur la tension (DHTv) en %,
La distortion harmonique maXImum sur le courant (DHTi) en %,
3 - Les déviations possibles des valeurs des filtre ou du systèmec
ausés par des dérives thermiques ou par la tolérence des
spécifications du fabriquant
La déviation maximum de la fréquence du réseau (dfI) en %,
La déviation maximum des valeurs des condensateurs utilisé par

les filtres accordés (dC) en % ,
La déviation maximum des valeurs des condensateurs utilisé par

les filtres accordés (dC) en % ,
Et de façon très optionnelle le coût relatif des condensateurs et
inductances des filtres
Le coût des condensateurs des filtres (JC) en $/kV AR,
Le coût des inductances des filtres (JL) en $/kV AR,
Note: Les paramètres électrique du système n'ont pas à être entrés par
l'utilisateur si on procède tout d'abord à une simulation numérique
dans le domaine temporel. Cette simulation est faite à partir des

modèles SIMULINK présentés à l'annexe A. Elle doit être suivit d'une
analyse fréquenciel des résultats avec le programme "analyse.m"
Il 7
fourni sur la même disquette que le programme de calcul de filtre. Ces
paramètres font alors partie du modèle et doivent être spécifiés dans
le bloc "source de tension avec impédance interne" et pour le modèle
de l'EVV complet (Moteur de la bonne puissance et de charge
adéquate) . Il est également possible d'utiliser un modèle simplifié de
la charge produissannt à peu près les mêmes formes d'onde de courant
de ligne. Ce modèle présenté à la figure B.1 permet de modifier
facilement l'amplitude de la charge en jouant sur la valeur de
l'impédance de charge "Zcharge RL (VI)"
Modele simplifie de l'onduleur de tension avec filtre Le

(produit un courant de ligne similaire a celui d'un EVV)
~ 1
Source de tension Y ILw
8
Vcc
Zcharge RL (UI)
Sous- systeme de: "mod_s_oUc.m"
Figure B.1: Schéma SIMVLINK du modèle simplifié d'un EVV vue du

réseau de distribution
118
Les valeur obtenues après l'exécution du programme peuvent être
réintroduites dans le modèle du système incluant des filtres passifs
accordés tel que présenté à la figure B.2. Il est alors possible de
vérifier l'effet d'une multitude de perturbations au système telles que:
- déviation en fréquence de la source,

- variation de l'impédance de la source,
- variation brusque de la charge,
- déviation des valeurs de C et L des filtres,
- ajout d'un deuxième EVV sur le transformateur,
- etc.
r---+i 4
Courants des filtres
Tensions
Filtres accordes (UI) au transfo
H5,H7,HP
Source de tension
avec imp. interne et
ajustement de tension
t---~3
L...-_ _ _ _ _ _ _ + 2
Courants d'entree de
Courants de source
l'onduleur
au transfo
Courants et tensions
sur le circuit CC
el
Ose. flottant
Sous-systeme de: "filtresJlass.m"
Figure B.2: Schéma SIMULINK du système d'alimentation triphasé

avec des filtres passifs accordés et un EVV comme
charge.
119
B.2 LISTAGE DU PROGRAMME DE CALCUL DE FILTRES PASSIFS
Le listage complet du programme de calcul de filtres passifs est

présenté ci-dessous. Le programe est également disponible sur une
disquette fourni en annexe.
% Programme Matlab de calcul de filtres accordes
% Par Sylvain Lahaie, mars 94
%
% Exemple pour trois filtres accordes sur h5, h7 et hll+HPF
% Reference
% 1- "POWER SYSTEM HARMONICS", J.Arrillaga,D.A. Brandley, P.S.Bodger,
% Edition WILEY, Chap. 10, pp. 296-304
% 2- "HVDC Power Transmission Systems", K.R. Padiyar, Edition WILEY, pp.146-156
%
% Les donne es necessaire au programme doivent etre incorpore a l'interieur du code matlab
% aux endroits ou il y a des "???"
%
% Etape 1
% Entreer la frequence fondamentale et la puissance du transformateur pour calculer
% l'impedance de la source a cette frequence
%
% Note: Habituellement ont associe l'impedance du transfo comme ayant un angle
% de phase constant en fonction de la frequence .
% Pour les transfo de 500 kV A et plus Xp est au moins 5 fois plus
% grand que Rp (78 deg)
%
% Exemple: S transfo= 20 kVA, 5 pu., Vp=200 volts 3ph, phi=60 degre
%
% Zsec = 2001\2120000 * 0.05;
% = 0.1 ohms avec un angle de 60 degres.
%
% *********************************************************-******************************-****
% Entrez les donnees suivantes:
% -Puissance nominale du transformateur en KV A (Snom),
% -son impedance en pu. (Imp)
% -angle de phase du transfo en degre, (phi_d)
% -tension de ligne-neutre en volt, (Vp )
% -frequence fondamentale en Hz. (fl_h)
% -taux de distortion maximum acceptable sur la tension en %
% -taux de distortion maximum acceptable sur le courant en %
% ????????????????????
SnomT = 500;
Imp = 0.05;
phi_d = 80;
Vp = 575/(31\0.5);
IL = 400;
fI = 60;
120
DHTv = 5;
DHTi = 15;
9& ????????????????????
Zs = (Vp"2 * Imp * 3 / (SnomT*lOOO) ) * ( cos( phi_d * pi/180) + i*(sin( phi_d * pi/180»);
9& **-********************************-*********************--**************************
9&
9& Les filtres accordees peuvent etre utilise afin de corriger le facteur de puissance
9& de la charge dont ont veut minimiser les harmoniques c'est pourquoi ont peut specifier
9& le nombre de kV A inductif que l'on veut corriger. Ce calcul suppose que l'on connait
9& bien la nature de la charge.
9&
9& Exemple pour un redresseur a diode de 6 impulsions le facteur de puissance est de 0.955
9& et les courants harmoniques sont les suivants
if exis t('ffti ')== 1

Snom=S/1000 ;
Ih(l: 101 )=ffti(2: 102);
vvh(1: 101)=fftv(2: 102);
else
p = 6;
Snom = 3 * Vp*IL/3"0.5/(1000) ;
9& Valeurs typique normalisee des harmoniques

9& d'un convertisseur de courant
Ih(lOl)=O;
Ih( 1)= -0 . 156446710 -0 .987686400*i

Ih(5)= 0.148064310 +0.140184420*i
Ih(7)= 0.104431590 +0.061761586*i
Ih(11)= -0.074166240 +0.011678924*i
Ih(l3)= -0.050942974 +0.024151776*i
Ih(17)= 0.015149245 -0.033635810*i
Ih(19)= 0.003102993 -0.028864241 *i
Ih(23)= 0.011650275 +0.014390335*i
Ih(25)= 0.012733232 +0 .007425433 *i
Ih(29)= -0.008581318 +0. 003990703 *i
Ih(31 )= -0 .005080063 +0. 005990499*i
Ih(33)= 0.000006860 -0.000051163*i
Ih(35)= -0.002413420 -0 .005663730*i
Ih(37)= -0.004146078 -0.003746264*i
Ih(41)= 0.004986067 -0.001296043*i
Ih(43)= 0.003816915 -0.002830746*i
Ih(47)= -0 .000177053 +0.004197943*i
Ih( 49)= 0.001228466 +0.003579645*i
9& Valeurs en grandeur reel

k=O;
121
for n=2:50
k = (abs(Ih(n)))A2 + k;
end
k=l+k"0.5;
Ih( 1: 50)=ILlk*Ih( 1:50);
end
% Calcul du courant RMS dans le condensateur du filtre passe-haut
It_hp=O;
for n=II:50
IChp = IChp + (abs(Ih(n)))"2;
end
It_hp=It_hp"0.5;
%
% Entrez le facteur de puissance de la charge et le facteur de puissance requis (Fpr)
% ou impose le nombre total de kV A requis (Qr) des filtres a fl
% 11111111111111111????1????????
if exist('Fp') == 0
Fp = 0.9;
end
Fpr = 1.0;
Qr = Snom* 1000*(1.0 -Fp"2)"0.5;
% 1????1??1??111??????1?????????
Qrr = Qr;
%********************-*********************************-***********************
%
% Le calcul du facteur de qualite des filtres depend de la stabilite de la frequence
% fondamentale et des valeurs de C et L en fonction du temps et des derive thermiques.
% On doit specifier ces parametres
% Exemple:
% Deviation de la frequence = 1 % = delta FI/FI
% Derive thermique caussant une variation de la valeur de C de 2 % = delta C / CO
% Derive thermique caussant une variation de la valeur de L de 3 % = delta L / LO
%
%**************************************************************************************
% Entrez la variation de FI, C et L en %
% ??1?1?11??1??11?11???11?
dfl = 2.0;
dC = 2.0;
dL = 2.0;
% ???11111?11?1111?1????1?
%*******--***-**---************--***********-******-**-
%
% Le calcul des filtres donnent une infinitee de combinaisons possibles.
122
% On peut tenter de faire le design afin de mInImiSer le cout des filtres
% Entrez le cout des condensateurs utilises en $/kVar (JC) et
% le cout de inductance en $/kV AR (JL)
% ??????????????????????????????????
JC = 0.6;
JL = 0.4;
% ???????????????????????????????????
%
% -*******-*********-**--****************************--*********--
% Debut des calculs
%
% Calcul de C total pour fournir l'energie reactive voulut a la frequence FI pour une phase.
% et distribuer sur les trois filtres.
Wl = 2*pi*f1;
Ls = imag(Zs)/WI;
Rs =real(Zs);
Ct = Qr/(3*Wl *Vp"2);
Drss=IOO; Drssi=100;
while (Drss > DHTv/lOO 1 Drssi > DHTi/lOO ),
% Clacul des derives thermiques et frequencielle
dT = (dfl + 0.5 * ( dC + dL » / 100 ;
% Clacul du facteur de qualite optimum.
phi = phi_d * pi/180;

Q = (1 +cos(phi))/(2*dT*sin(phi));
% calcul iteratif pour tenter de distribuer la capacite total selon un critere d'economie
% ce calcul debute avec les valeurs initial du courant de ligne.
I5=Ih (5); I7=Ih(7) ; III=It_hp;

Zs5 = Zs*5; Zs7=Zs*7 ; Zs 1 1=Zs * Il ;
A5 = JC + JLl25;
B5 = «abs(I5»"2 * (abs(Vp»"2)/5 * ( JC + JL );
Qf5 = ( B5 / A5 )"0.5;
A7 = JC + JLl49;
B7 = «abs(I7»"2 * (abs(Vp»"2)/7 * ( JC + IL );
Qt7 = ( B7 / A7 )"0.5;
123
Qfll=O;
for n=11 :50
AlI = JC + JLln"2;
B Il = «abs(Ih(n»)"2 * (abs(Vp»"2)/n * ( JC + JL );
Qfll = ( B Il / All )"0.5 + Qfll;
end
Qt = Qf5 + Qt7 + Qfll;

CC5= Qf5/(Vp*Vp*120*pi);
CC7= Qf7/(Vp*Vp*120*pi);
CClI= Qfll/(Vp*Vp* l20*pi);
CCt=CC5+CC7+CCll;
if CCt > Ct
Ct=CCt;
end
C5 = Qf5 /Qt * Ct;

C7 = Qt7 /Qt * Ct;
CIl = Qfll/Qt * Ct;
% Calculer l'impedance correspondant a ces capacite pour les harmoniques 5,7 et Il.
L5 = 1 / «W1*5)"2 * C5 );
L 7 = 1 / ( (W1 *7)"2 * C7 );
% Pour le filtre passe-haut ont choisit hO < 1.41 *hmin (ici hmin=ll hOmax=15.5)
hO = Il;
LI1 = 1 / ( (WI *hO)"2 * Cll);
% Calcul des resistances serie total des branches 5 et 7
R5 = WI *5*L5 / Q;
R7 = WI*7*L7 / Q;
Rll = W1 *hO*LIl / 30;
% Calcul de la resistance parallele du filtre passe-haut

% m = L / ( R"2 * C) ou: "m" est compris entre 0.5 et 2 ([1] p.308)
m = 0.5;
Rhp = ( L Il/(C II *m) )"0.5;
% Calcul du taux de distortion sur la tension
Ifht=O; Drss=O; Drssi=O;

for h=1:50
Z5 = R5 + i*( WI *h*L5 - ( l/(WI *h*C5) ) );
Z7 = R7 + i*( Wl*h*L7 - (1/(Wl*h*C7»);
Zll = 1I( 1IRhp + II(RIl+(i*W1 *h*LIl» ) + 1I( i*W1 *h*Cll );
Zfh = 1/ (l/Z5 + IIZ7 + IIZ11);
Zsh=Zs*h;
124
if h==1
Vh( 1)=-Vp*i-Ih( 1)*Zs;
Ifh(l) = -Vh(I)/Zfh;
%fh(l) = Vp*i/Zfh
else
Ifh(h) = Ih(h) * Zsh / ( Zsh + Zfh );
Vh(h) = -Ih(h) / ( l/Zsh + l/Zfh );
Drss = Drss + abs(Vh(h»"2;
Drssi= Drssi+ abs( Ih(h)-Ifh(h) )"2;
end
Ifht = Ifht + abs(Ifh(h»"2;
end
Ifht = Ifht"0.5;
DF = (Drss/Vp"2)"0.5;
Drss = Drss"0.5 / Vp;
Drssi = Drssi"0.5 / abs( Ih(I)-Ifh(1) );
Ct = Ct*1.05;
end
% affichage des resultats
f5= l/(2*pi*(L5*C5)"0.5);
f7= l/(2*pi *(L 7*C7)"0.5);
fll=1I(2*pi*(LlI *CII)"0.5);
R5 = WI *5*L5/Q;
R7 = Wl *7*L7/Q;
Qr = (C5+C7+Cll) * (3*Wl *Vp"2);
sI = sprintf(' L5= %5 .3f mH C5= %5.lf uf R5 = %5.3f ohms (Q=%2d) Fr= %d Hz\n',
L5*1000, C5*le6, R5, round(Q), round(f5»;
s2 = sprintf(' L7= %5.3f mH C7= %5.1f uf R7 = %5.3f ohms (Q=%2d) Fr: %d Hz\n',
L7*1000, C7*le6, R7, round(Q), round(f7» ;
if(LlI > 0.le-3),
s3 = sprintf(' Lph=%5.3f mH Cph=%5.lf uf Rph = %5.lf ohms Fr= %d Hz\n\n',
Lll*1000,Cll*le6, Rhp, round(fll»;
e Ise
s3 = sprintfC Lph=%5.2f uH Cph=%5.1f uf Rph = %5.1f ohms Fr= %d Hz\n\n',
Lli * le6,Cli * le6, Rhp, round(fll»;
end
s-+ = sprintf(' Le nombre de kVAR requis pour obtenir Fp = %2. lf est de %3.lf kVar\n', Fpr,
round(QrrIlOO)/lO );
s5 = sprintfC Le nombre de kVAR intalle est de %2.1f kVar\n', round(Qr/l00)/10 );
s6 = sprintf(' La distortion harmonique sur la tension est: DHTv=%3.2f\n', Drss*100 );
s7 = sprintf(' Le distorsion harmonique sur le courant est: DHTi=%3.2f\n', Drssi* 100 );
s8 = sprintf(' Le courant total dans les filtres est de %4.lf amp. RMS\n', Ifht );
Resultat = sprintf('%s%s%s%s%s%s%s%s' ,s l,s2,s3,s4,s5 ,s6,s7 ,s8)
% Calcul de l'impedance des filtres
125
ind = 1;
for f=fl: 1:3000
Z5(ind) = R5 + i*( 2*pi*f*L5 - ( l/(2*pi*f*C5» );
Z7(ind) = R7 + i*( 2*pi*f*L7 - ( l/(2*pi*f*C7» );
Zh(ind) = 1I( lIRhp + 1I(Rl1+(i*2*pi*f*L11») + 1I( i*2*pi*f*Cll );
Zt(ind) = 11 ( lIZ5(ind) + lIZ7(ind) + lIZh(ind) );

ff(ind)=f;
ffh(ind)=flfl ;
ind = ind+l;
end
figure(l)
clf
subplot(2, 1, 1);
loglog(ffh,abs(Z5),'r', ffh, abs(Z7),'g', ffh, abs(Zh),'b' )
xlabeI('Rang de l"harmonique')
yIabeI('ModuIe de l"impedance (Ohms)')
%titIe('lmpedance des filtres accordes 5, 7 et passe-haut')
grid
subpIot(2,I,2);
axis([2, 100,2, 10])
axis off
text(2,9,sl);
text(2,8,s2);
text(2,7,s3);
text(2,6,s4);
tex t(2,5 ,s5);
text(2,4,s6);
text(2,3,s7);
drawnow
% press any key to continue

figure(2)
clf
subpIot(2,2,1 );
logIog(ffh,abs(Zt),'b' )
yIabeI('ModuIe(IZtl) (Ohms)')
xIabeI('Rang de l"harmonique')
%title('lmpedance vue par la charge en fonction de la frequence ')
grid
subpIot(2,2,2);
plot(real (Zt),imag(Zt),' g'
xlabel('Re(Zt) (Ohms)')
ylabel('Im(Zt) U Ohms)')
% title('Impedance vue par la charge dans le plan complexe')
grid
hold on;
pos=axis;
indice = 1;
126
for tesCbode=1 : 1:50
if (abs(Ih(tescbode» > abs(O.Ol *Ih(1» 1 tesCbode < 13 )
while (-(ffh(indice) == tescbode» & (indice < size(ffh,2»
indice = indice+ 1;
end
end
if(ffh(indice) == test_bode)
plot( real(Zt(indice) ),imag(Zt(indice», '0'
if tesCbode < 23
text( real(Zt(indice) )+pos(2)/30,imag(Zt(indice) ),sprintf('%d' ,tescbode)
end
end
end
pos=axis;
plot( pos( 1)+pos(2)/1 0,pos(3 )+(pos( 4 )-pos(3»/ 10, '0' )
text( pos( 1)+pos(2)/1 0+pos(2)/40,pos(3 )+(pos( 4 )-pos(3 »/1 0, 'Rang de l"harmonique' )
t=O: l/(fl *4000) :2/fl;
h=I:1 :50;
Vh(1 )=-Vp*i+Ih(l)*Zs;
for x=I:1:8001
v_t(x)= 1.4142 *sum(abs(Vh( 1:50». *cos(W 1*h*t(x)+angle(Vh( 1 :50»»;
i_t(x)= 1.4142 *sum(abs(Ih( 1:50)-Ifh( 1 :50» . *cos(W 1*h *t(x)+angle(Ih( 1 :50)-Ifh( 1 :50»»;
end
hold on
%subplot(2,1,2)
%title('Impedance vue par la charge en fonction de la frequence ou dans le plan compexe')
figure(3)
clf
subplot(2,2 , 1)
plot(t , v _t,'g')
title('tension harmonique ligne-neutre ')
xlabel('temps (sec)')
ylabel('tensiion (volt)')
s=sprintf('Le THD V l-n est de %4.2f %%',Drss*100) ;
text(O.1 *max(t) ,O. 7*max(v _t) ,s)
subplot(2,2,2 )
plot(t, i_t,'g' )
title('courant harmonique de ligne')
xlabel('temps (sec)')
ylabel('courant (amp)')
s=sprintf('Le THD IL est de %4.2f %%',Drssi* 100);
text(O . 1*max(t),0.7*max(i_t) ,s)
subplot(2 ,2,3)
hv _t=fft(v _t( 1 :4001»;
bar(abs(h v _t(2: 50) )/abs(h v _t(2», 'r')
subplot(2,2,4 )
hi_t=fft(i_t(l :4001»;
bar( abs(hi_t(2: 5 0) )/abs(hi_t(2», 'r')
127
ANNEXE C
SIMULATEUR DU FILTRE ACTIF DE FUJITA ET AKAGI[ll]
On présente dans cette annexe le simulateur que nous avons

développé du filtre actif de Fujita et AkagHll]. La figure C.I présente le
schéma SIMULINK du simulateur. Dans cet exemple la charge du
redresseur triphasé est très inductive. On obtient alors les mêmes
formes d'onde de courant que celles présentées dans l'article de Fujita.
19
Fioating Scope
Isw Vw '----+1 utw

w5
CictJit de controle
Source de tension
avec imp. interne
eTv
(9,
Clock
Figure C.I: Schéma SIMULINK du système comprenant des filtres

passifs accordés et un EVV comme charge.
128
ANNEXED
SIMULATEUR DU GHT
On présente dans cette annexe une description sommaue du

simulateur du GHT. Le simulateur a été développé sur une station de
travail avec la version Matlab 4.2 et Simulink 1.3. Les fichiers
devraient être compatibles sur PC pour la même version de Simulink
bien que les temps de simulation risque d'être excessif pour ce type
de machine puisque le modèle qui a été développé est très complexe.
On présente à la figure D.1 le schéma de sortie du simulateur. On peut

observer que tous les paramètres importants du système sont
sauvegardés dans un fichier pour analyse ultérieure . Si on ouvre le
bloc "Simulateur" on découvre le schéma de la figure D.2. Celui-ci est
composé de cinq blocs principaux:
1 Le circuit de contrôle (figure D.3 et D.4)
2 Les filtres accordés au 5e , 7e et passe-haut
3 Les filtres accordés à la fréquence fondamentale
4 La source de tension qui est un auto transformateur triphasé

motorisé
5 La charge qui est le modèle d'EVV présenté à l'annexe A ou

un modèle simplifié comme celui de la figure B 1.
129
Simulateur du banc d'essai d'entrainemnt a vitesse variable du LTEE: Generation de cinq harmoniques de tension au choix.
Simulateur du banc r---l-__~____~____+-____~________________~MUX

d'essai d'E.V.V. du I--....__ ~~ __-+I
L T.E.E. - H.a.
Simulateur parametres de revv
FIChier: "gh'-blcombi.m"
Figure D.I: Schéma de sortie du simulateur
Tensions du filtre actif
Source de tension
avec imp. interne et
ajustement de tension
Cieuit de controle
2
Courants de source
au transfo
Courants et tensions
sur le circuit CC
el
Ose. flottant
~------------------~3
Tension de ligne
cote charge Sous-systeme de: "ght_bfcombi,m"
Figure D.2: Schéma principal du simulateur
130
C ircuit de CXlntrQte
Tension de ligna
cota transfo.
IUII.IV'I 'ond.
Generateur d'harmoniques 2
Figure D.3: Schéma principal du circuit de contrôle des harmoniques
pid
CirOJit de contra'- - Generateur d"harmon~.
paramaltndu
pid
o alpha
Tension abeta
V alpha, bota
Cajeu'''''po'q Trans. inv
el
Osc. t
Polar '0Canes_
Sous-systame da: -ghCbfcombl.m-
Figure DA: Schéma du circuit de contrôle pour un harmonique
131
Il est à noter que les lignes à trait larges désignes plusieurs lignes
multiplexés (habituellement 3 lignes pour le réseau triphasé).
Chacun de ces blocs peuvent être modifiés de manière à simuler

différents systèmes. Une librairie de modèle d'EVV a été faite afin de
pOUVOIr modifier rapidement notre modèle de simulation en "coupant"
l'ancien bloc et en "collant" le nouveau bloc provenant de la librairie.
Le même principe s'applique au circuit de contrôle qui sera choisi en
fonction du test d'immunité que l'on désire exécuter.
Les modèles d'EVV disponible sont:
1- Modèle simplifié avec:
filtre Lf-Cf sur le lien à tension continue:

50, 100, 200, 400, 800 HP
filtre Cf sur le lien à tension continue:

50, 100, 200, 400, 800 HP
filtre Le à l'entrée et condensateur Cf sur le lien à tension continue:

50, 100, 200, 400, 800 HP
2- Modèle complet avec:
filtre Lf-Cf sur le lien à tension continue:

Modèle MLI et "six steps" 50, 100, 200, 400, 800 HP
filtre Cf suc le lien à tension continue:

filtre Le à l'entrée et condensateur Cf sur le lien à tension continue:

132
Les modèles disponibles de circuit de commande sont:
Modèle pour les tests de commutation de condensateurs
Modèle pour les tests de génération d'encoche de tension
Modèle pour les tests de génération d'harmonique:

En boucle ouverte pour une séquence de 5 harmoniques
appliqués successivement
En boucle ouverte pour l'addition de 5 harmoniques différentes
En mode de régulation pour l'addition de 1, 4 et 9 harmoniques

différentes
Des fichiers d'analyse, écrit en langage Matlab, ont également été

développés afin de faire l'analyse fréquenciel des sIgnaux de sortie.
133
ANNEXE E
DÉTAIL DES SIMULATIONS DU GHT
On présente dans cette annexe quelques résultats suplémentaires des

simulations présentés au chapitre 4 sur la superposition
d'harmoniques de tension
Le tableau E.1 présente les essais de superposition d'un harmonique

de tension procurant une DHT de 10 % sur la tension.
Tableau E.1: Résultats de simulation de superposition d'un harmonique
H V* Vbo Vreg err bo err bt

5 33 33.9 33.0 -2.7 0.0
7 33 33.5 32.7 -1.5 0.9
11 33 31.8 32.8 3.6 0.6
13 33 32.1 33 2.7 0
17 33 29.5 32.9 10.6 0.3
19 33 29.9 32.9 9.4 0.3
23 33 25.9 32.7 21.5 0.9
31 33 28.3 14.2
35 33 27.9 33.02 15.5 -0.1
41 33 27.7 33.01 16.1 0
43 33 27.9 33.28 15.5 -0.8
47 33 28.5 33 13.6 0
49 33 32.98 0.1
134
La figure E.l présente les résultats du tableau El de façon graphique.
35 25
*=* - t • •
30
20
25
15
...
.--.
Q 20
.--.
~
:>
'-'
'-'
"-
.•=
~
CIl
10
=
Q,/
"-
"-
~ 15 r.:l
~
5
10
o
5
/
o -5
o 10 20 30 40 50
• V* ~ Vbo - -. - Vreg
- ------- e r r ---0--- e rr
bo bf
Figure E.l: Résultats de la superposition d'un harmonique
135

PDF

Transféré par

Droits d'auteur :

Formats disponibles

PDF

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

PDF

Transféré par

Droits d'auteur :

Formats disponibles

UNIVERSITÉ DU QUÉBEC À TROIS-RIVIÈRES

COMME EXIGENCE PARTIELLE POUR L'OBTENTION D'UN

ÉTUDE ET CONCEPTION D'UN BANC D'ESSAI:

L’auteur de ce mémoire ou de cette thèse a autorisé l’Université du Québec

Cette diffusion n’entraîne pas une renonciation de la part de l’auteur à ses

Ce mémoire de Maîtrise présente la conception d'un banc d'essai pour

La conception prend comme point de départ la description des essaIS

Ce travail couvre la partie du projet se rapportant aux simulations du

La partie critique de la conception du banc d'essai se retrouve dans la

En plus de la conception du GRT, ce mémoire contribue au

Ce travail présente également l'étude d'un filtre actif dont la

Mes remerciements vont vers mon directeur de recherche, le

Je remercie le chercheur post-doctoral Yao Ziwen, pour ses critiques et

Je remerCIe aussi les membres du corps professoral de la maîtrise en

Je remerCIe la direction d'Hydro-Québec; spécialement M. Antoine

Finalement Je désire remerCIer tout particulièrement M. André

REMERCIEMENTS .................................................................................................... 111

TABLE DES MATIÈRES ........................................................................................... V

LISTE DES FIGURES ................................................................................................. XI

LISTE DES TABLEAUX ...................................................................................... XVII

LISTE DES SYMBOLES ET DES ABRÉVIATIONS ..................................... XVIII

1.1 Cadre et motivation du projet.. .............................................................. 1

1.2 Simulation numérique du banc d'essai .............................................. 2

1 .3 Conception du banc d'essai.. .................................................................... 2

. t'Ion des sUJe

2.1 Introduction ................................................................................................... 5

2.2 Notions de base sur les harmoniques et sur le calcul des

2.2.1 Perturbations de la forme d'onde de tension ......................... 6

2.2.2 Simulation des configurations les plus communes

2.3 Réduction des harmoniques injectés sur le réseau .................... 1 7

2.3.1 Généralités ............................................................................................. 1 7

2.3.2 Conception des filtres passifs accordés ................................... l 9

2.3.2.1 Calcul des filtres dans le domaine fréquenciel.. ........... l 9

2.3.2.2 Description du programme de calcul de filtre .............. 20

2.3.2.3 Critères de conception des filtres ....................................... 2 1

2.3.2.4 Étapes de design des filtres .................................................. 22

2.3.2.5 Organigramme du programme de calcul des filtres ... 36

2.4 Exemple de calcul de filtres accordés .............................................. 3 9

2.4.1 Critères de conception spécifiés ............................................... 39

2.4.2 Résultats du calcul de filtre ......................................................... 39

2.4.3 Analyse des performances des filtres dans le

2.5 Conclusion .................................................................................................... 43

3.1 Introduction ................................................................................................ 44

3.2 Notions de base sur les filtres actifs ................................................ .4 5

3.2.1 Configuration classique de filtres actifs ................................ .4 5

3.2.2 Combinaison de filtres actifs et passifs ................................. .4 6

3.3.1 Portée de l'étude ................................................................................. 4 8

3.3.2 Description de la configuration du circuit.. .............................. 49

3.3.3 Description du circuit de contrôle ............................................... 50

3.3.4 Analyse des performances du filtre actif ................................ 53

3.3.4.1 Problème de résonnance parallèle .................................... 53

3.3.4.2 Problèmes de résonnance série .......................................... 56

3.3.5 Simulateur du filtre actif de Fujita et Akagi.. ...................... 59

3.3.5.2 Exemple de simulation d'un problème de résonnance

3.4 Conclusion ...................................................................................................... 63

4.1 Introduction .................................................................................................. 64

4.2 Étude des perturbations de la tension de ligne ............................. 65

. 4.2.1 Fluctuation de la fréquence de 60 Rz ...................................... 6 6

4.2.2 Variation de l'amplitude: .............................................................. 66