Compensation active des vibrations électriques
Compensation active des vibrations électriques
machines électriques
Pierre Granjon
THESE
DOCTEUR DE L’INPG
par
Pierre GRANJON
le 19 décembre 2000
Titre :
JURY
Merci aussi à la crème de la recherche grenobloise, JIM et POA, qui ont démontré
expérimentalement que mesurer plus de 180 cm, faire du foot, du ski, des enfants, des concours
de boules privés, quelques excès et être chercheur pouvait très bien être des caractéristiques
d’une seule et même personne. Cela ne va tout de même pas sans une ou deux blessures dues
principalement à la fatigue : fracture du coude causée par un ballon récalcitrant, déchirure
musculaire causée par un survêtement rebelle, etc.
i
REMERCIEMENTS
dernières fantaisies du Grenoble «by night» (et viva Argentina !), même s’il avait déjà
rencontré pas mal de monde qui comprenne et partage sa passion de la nuit,
– Cédric, pour avoir apporté beaucoup de finesse et de sagesse dans ce labo de brutes,
– Steph, éminent spécialiste des langues (étrangères ou non, d’ailleurs), qui a survolé le
«championnat du placement incongru de mots complexes dans les conversations les plus
banales»,
– Seb, pour avoir courageusement infiltrés les rangs ennemis pendant plusieurs mois et au
retour de cette mission, ramené des informations capitales pour notre lutte perpétuelle
contre les américains,
– Steeve, pour m’avoir fait découvrir que je n’aimais décidément pas le zouc,
– Antoine, avec plaisir,
– Cyrille, Cyril, Syril ou Sirryhl, je sais plus,
– aux inséparables compères Jacky et Jean-Jacques pour leur science innée des véhicules
motorisés,
– à tous les autres que je n’ai pas cité mais qui se reconnaı̂tront. . .
Je n’oublierai pas bien sûr de citer mes chers et dignes géniteurs ainsi que ma grand-mère
Virginie de qui j’ai hérité «la fibre».
Pour finir, merci aux sourires et à la patience de Linda, Romain et Mathis, et félicitations
pour avoir réussi à me supporter et à me remettre les pieds sur terre durant les sombres mois
de la rédaction.
ii
Table des matières
Remerciements i
Introduction 1
iii
TABLE DES MATIÈRES
iv
TABLE DES MATIÈRES
Bibliographie 123
v
Table des figures
vii
TABLE DES FIGURES
3.3 Exemple de réponse impulsionnelle g(r, k) d’un système LVPT causal avec P = 2 48
3.4 Structure à modulation d’entrée d’un système LVPT P -périodique . . . . . . 49
3.5 Structure à modulation de sortie d’un système LVPT P -périodique . . . . . . 50
3.6 Exemple de relation entre les spectres d’entrée et de sortie d’un système linéaire
stationnaire et causal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.7 Exemple de module de réponse bifréquentielle pour un système LVPT P -
périodique causal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
3.8 Exemple de relation entre les spectres d’entrée et de sortie pour un système
LVPT P -périodique causal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.9 Représentation vectorielle de la transformée de Fourier d’un signal s(n) pour
P =3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
3.10 Propriétés de base du gain complexe matriciel . . . . . . . . . . . . . . . . . . 62
3.11 Modulation d’amplitude sinusoı̈dale . . . . . . . . . . . . . . . . . . . . . . . 64
ˆ
3.12 ĥmod (α, β) correspondant au système de la figure 3.11 . . . . . . . . . . . . . 65
3.13 Représentation graphique de la décomposition en valeurs singulières du gain
complexe matriciel Ĝ(α) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
3.14 Structure interne d’un système LVPT H . . . . . . . . . . . . . . . . . . . . . 68
3.15 Structure interne du système inverse Hinv . . . . . . . . . . . . . . . . . . . . 69
3.16 Modèle du système LVPT à temps continu à analyser . . . . . . . . . . . . . 73
3.17 Caractéristiques du système G . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
3.18 Exemple schématique de module de la réponse bifréquentielle (3.61) . . . . . 76
3.19 Fraction de la réponse impulsionnelle (3.70) en niveaux de gris . . . . . . . . 78
3.20 Structure à modulation d’entrée du système LVPT réel . . . . . . . . . . . . . 79
3.21 Fraction du module de la réponse bifréquentielle (3.74) . . . . . . . . . . . . . 80
3.22 Gain complexe matriciel du système de la figure 3.16 . . . . . . . . . . . . . . 81
3.23 Module du gain complexe matriciel Ĝ(0) . . . . . . . . . . . . . . . . . . . . . 83
3.24 Valeurs singulières du gain complexe matriciel Ĝ(0) . . . . . . . . . . . . . . 84
3.25 Module du gain complexe matriciel Ĝinv (0) . . . . . . . . . . . . . . . . . . . 84
3.26 Valeurs singulières du gain complexe matriciel Ĝinv (0) . . . . . . . . . . . . . 85
viii
LISTE DES FIGURES
ix
Introduction
De nos jours, les machines tournantes électriques sont utilisées dans des domaines aussi
variés que le transport (trains, véhicules motorisés, etc), la production électrique (alter-
nateurs), l’industrie de production, ou encore l’électroménager. Les vibrations qu’elles en-
gendrent peuvent être gênantes à plus d’un titre :
– elles sont à l’origine d’une partie du bruit rayonné par cette machine et sont donc
indésirables pour les utilisateurs situés à leur proximité,
– elles peuvent être transmises aux structures avoisinantes et accélérer leur détérioration
ou leur vieillissement,
– elles peuvent également endommager les machines tournantes elles-mêmes.
Il est donc intéressant de développer des méthodes visant à réduire le niveau vibratoire de
ces machines, ainsi que celui de leur environnement direct.
Ces méthodes de compensation peuvent en général être qualifiées de passives ou d’actives.
Comme nous le verrons au chapitre 1, les premières n’utilisent aucun apport supplémentaire
d’énergie. Au contraire, les méthodes actives minimisent les vibrations naturelles à l’aide de
«contre-vibrations» engendrées par des actionneurs nécessitant une alimentation extérieure.
L’objectif de ce travail de recherche est d’élaborer une méthode active de compensation des
vibrations d’une machine tournante électrique, dont le but est de réduire le niveau vibratoire
d’une zone de sa carcasse statorique. Cette zone pourra alors être utilisée pour la fixation de
cette machine et permettre ainsi de minimiser la quantité de vibrations qu’elle transmet à
son environnement direct. L’originalité de cette méthode tient au fait qu’elle n’utilise aucun
actionneur externe, ceux-ci étant remplacés par les enroulements statoriques de la machine
électrique parcourus par des courants de commande additifs. Le but de ce système de contrôle
actif est donc d’obtenir la valeur «optimale» de ces courants, minimisant la puissance des
signaux de vibrations mesurés par des capteurs accéléromètriques fixés sur le stator de la
machine considérée.
Le premier chapitre a pour vocation de définir entièrement le problème à résoudre. L’ob-
jectif général ainsi que le principe physique utilisé par la méthode de compensation à élaborer
sont tout d’abord exposés à la section 1.2. Les vibrations devant être atténuées sont ensuite
caractérisées à la section 1.3. Enfin, le système sur lequel est appliqué le contrôle actif, qui
correspond au transfert situé entre les courants de commande parcourant les enroulements
statoriques de la machine et les contre-vibrations qu’ils engendrent, est déterminé à la sec-
tion 1.4. Nous verrons que sous certaines hypothèses, il peut être considéré comme linéaire
1
INTRODUCTION
2
Chapitre 1
1.1 Introduction
Le but de ce premier chapitre est de spécifier le sujet du travail de recherche présenté,
c’est à dire la compensation des vibrations d’une machine tournante électrique, en utilisant
les forces magnétiques comme moyen d’action.
La section 1.2 expose l’objectif du travail, c’est à dire une méthode d’isolation vibratoire
d’une machine synchrone. Le principe physique employé pour atteindre ce but, ainsi que la
machine sur laquelle des tests seront effectués, y sont présentés.
Dans la section 1.3, nous détaillons l’origine des vibrations d’une telle machine, et présen-
tons un rapide panorama des méthodes d’analyse vibratoire. Nous donnons ensuite l’allure
théorique de son spectre vibratoire à partir des sources de vibrations définies précédemment,
et vérifions ces hypothèses sur un signal expérimental provenant de la machine décrite à la
section précédente. Le paragraphe 1.3.4 permet aussi de déterminer les composantes les plus
gênantes de ce signal vibratoire, qui doivent être éliminées en priorité pour atteindre le but
fixé.
Enfin, le modèle mathématique du transfert à travers lequel les forces magnétiques agissent
sur les vibrations est établi à la section 1.4. Nous verrons qu’il s’agit d’un système linéaire
variant périodiquement dans le temps (LVPT), nécessitant l’emploi de techniques de com-
pensation appropriées. À la fin de ce chapitre, le problème physique est finalement modélisé
en termes de traitement de signal, comme la génération de signaux de commande permettant
d’atténuer «au mieux» un signal périodique à travers un système LVPT.
1.2.1 Objectif
Un problème fréquemment rencontré en analyse vibratoire est celui où une structure
vibrante «contamine» son environnement direct, en lui transmettant tout ou partie de ses
vibrations par l’intermédiaire de ses points de fixation. La structure vibrante constitue alors
3
CHAPITRE 1. Définition du problème traité
la source de vibrations, tandis que les structures recevants ces vibrations sont appelées les
«structures réceptrices». Un exemple simple est représenté sur la figure 1.1, où les vibrations
d’une machine tournante sont transmises à une structure réceptrice par l’intermédiaire du
socle de fixation de cette machine.
L’action d’isoler les structures réceptrices des vibrations de la source est appelée «isolation
machine tournante
=
source vibratoire
vibrations transmises
vibratoire» [Ful96, pages 186-187]. Elle consiste à minimiser la transmission des vibrations
de la source vers les structures réceptrices, même si cette minimisation peut conduire à
augmenter le niveau vibratoire global (c’est à dire spatialement moyenné) de la source [Ber97].
En ce sens, elle est à différencier de l’atténuation vibratoire globale [Ful96, pages 153-154],
qui a pour but de minimiser le niveau vibratoire global de la source, mais qui peut provoquer
l’augmentation de la quantité de vibrations transmises à son environnement direct.
Le but du travail de recherche exposé dans ce mémoire est de réaliser l’isolation vibratoire
d’une machine synchrone, pour la configuration présentée sur la figure 1.1. Le niveau vibra-
toire global de la machine n’est donc pas mis en cause, seule la transmission de ses vibrations
doit être atténuée.
La section suivante expose le principe physique utilisé pour agir sur les vibrations de la
machine et réaliser l’isolation vibratoire recherchée.
4
1.2. Point de départ
constituées de matériaux absorbant les vibrations. La forme de la source vibratoire peut aussi
être optimisée pour la rendre plus discrète et donc minimiser son niveau vibratoire global (ce
principe a par exemple été employé dans [Der92] et [Far95] sur une machine synchrone).
Les méthodes utilisant ce principe ne nécessitent aucun apport supplémentaire d’énergie
au système, d’où leur dénomination de «passives». Malheureusement, leur coût de mise en
œuvre peut rapidement devenir prohibitif (par exemple, changer les matériaux constitutifs ou
la forme d’une machine existante). De plus, elles sont peu efficaces sur les vibrations situées
dans les basses fréquences (typiquement en dessous de 1 à 2 kHz) [Ful95].
Ces inconvénients peuvent être contournés en utilisant des actionneurs extérieurs au système
comme moyens d’action sur les vibrations. Parmi les actionneurs les plus couramment utilisés,
on peut remarquer les électroaimants [Ber97], les vérins hydrauliques ou pneumatiques [Ful96,
chapitre 7], les actionneurs piézoélectriques [Ful96, chapitre 5], etc. Pour pouvoir être mis en
œuvre, ils nécessitent un apport d’énergie extérieure, ce qui explique la dénomination de
«méthode active». La bande passante de ces actionneurs permet un contrôle efficace pour les
basses fréquences où les méthodes passives n’ont que peu d’effet, et leur coût de revient est
directement lié au type et au nombre d’actionneurs utilisés.
Afin de réaliser l’isolation vibratoire d’une machine tournante électrique classique, on
propose une méthode active «sans actionneur». Une telle machine est constituée d’un rotor
et d’un stator, chacun générant son propre champ d’induction. Leur somme à l’instant t et
−
→
à la position angulaire θ forme le champ total b (t, θ). Si on suppose ce champ radial, on
peut montrer [Juf95, pages 41-47], [Car67, pages 202-209] que l’élément de surface ds de la
carcasse statorique situé à la position angulaire θ, subit à l’instant t une force f (t, θ) de
direction radiale et de module :
−
→ 2
k b (t, θ)k
f (t, θ) = ds (1.1)
2µ0
−
→ −
→
où µ0 est la permittivité du vide et k b (t, θ)k le module du vecteur b (t, θ).
Cette force, dénommée «force de Maxwell», est donc à tout instant et en tout point de la
surface du stator proportionnelle au carré de l’induction magnétique. C’est en utilisant ce
principe que l’on peut espérer pouvoir agir sur les vibrations naturelles de la machine, pro-
voquées par sa rotation. En effet, supposons que l’on alimente les bobines statoriques avec des
courants additionnels appelés «courants de compensation». Ces courants créent un champ
magnétique additionnel dans l’entrefer, qui, d’après la relation (1.1), engendre une force en
tout point de la surface du stator. Celui-ci réagit par des vibrations supplémentaires qui,
ajoutées à celles générées par la rotation de la machine, forment les vibrations totales de la
carcasse statorique. Ce principe physique est résumé par la figure 1.2. Ces vibrations addi-
tionnelles peuvent être utilisées comme des «contre-vibrations» par rapport aux vibrations
naturelles, afin de minimiser les vibrations totales du stator au voisinage d’une position angu-
laire donnée. On pourra alors utiliser cette «zone de calme» pour fixer la machine et réaliser
son isolation vibratoire vis-à-vis de la structure réceptrice. Il faut toutefois garder à l’esprit
que la force donnée par (1.1) suit une direction radiale. Par conséquent, les seules vibrations
5
CHAPITRE 1. Définition du problème traité
vibrations
naturelles
du stator
que l’on peut atténuer par cette méthode, et donc les seules vibrations considérées dans la
suite de ce document, sont les vibrations statoriques radiales.
Le principe physique décrit dans cette section permet d’avoir une idée précise de l’im-
plantation pratique d’une telle méthode. Elle est représentée par la figure 1.3, où une coupe
transversale de la machine tournante de test est située sur la gauche. Un ou plusieurs cap-
vibrations
amplificateur
de
capteurs puissance
de processeur
vibrations
DSP
teurs, localisés autour d’une position angulaire donnée, sont utilisés pour mesurer les vibra-
tions statoriques à atténuer. Les signaux ainsi obtenus sont discrétisés, puis un processeur de
traitement du signal est chargé d’élaborer les signaux de commande nécessaires à leur mi-
nimisation. Les sorties de ce processeur alimentent, par l’intermédiaire d’un amplificateur de
puissance, une ou plusieurs des bobines statoriques de la machine. Les actionneurs employés
(les bobines statoriques de la machine) font donc partie intégrante du système cible, ce qui
constitue un intérêt économique et pratique considérable, et donne un sens à l’expression de
méthode «sans actionneur».
Ce principe a été employé pour la première fois dans [Mae87] sur des machines à courant
continu. Il a ensuite été adapté aux machines synchrones et implanté sous forme analogique
dans [Der92] et [Far95], puis sous forme numérique dans [Jar95]. Nous verrons dans le second
chapitre, qu’une non-stationnarité présente dans le transfert entre les courants de commande
et les vibrations qu’ils engendrent, amoindrit considérablement les performances obtenues
6
1.2. Point de départ
dans cette application. Le but de ce travail de recherche est donc d’élaborer un algorithme
de commande prenant en compte cette non-stationnarité, et minimisant «au mieux» les
vibrations du stator autour d’une position angulaire donnée.
Après avoir présenté le principe physique utilisé par la méthode de compensation à
élaborer, la section suivante décrit les caractéristiques et les particularités du banc d’essais
sur lequel cette technique sera testée.
Afin de tester les performances de cette technique d’isolation vibratoire sur une machine
réelle, on dispose d’un banc d’essais constitué de deux machines. La première est la machine
de test, c’est à dire celle sur laquelle on doit mettre en œuvre l’isolation vibratoire. C’est une
machine synchrone triphasée à deux paires de pôles, d’une puissance de 10 kW. Sa struc-
ture est le résultat des travaux de recherche menés dans [Der92] et [Far95], visant à atténuer
son niveau vibratoire global à l’aide de techniques de contrôle passif. Sa coupe transver-
sale est représentée sur la figure 1.4. Elle est constituée d’un stator à 36 encoches inclinées
carcasse
rotor lisse statorique
à aimants
permanents
bobinage
statorique
continûment d’un pas dentaire, avec un bobinage à pas diamétral. Cette inclinaison contribue
fortement à la discrétion vibratoire de la machine [Der92, Far95]. On peut noter sur cette fi-
gure l’absence d’ailettes de refroidissement afin de pouvoir disposer les capteurs de vibrations
au plus près de la carcasse statorique, et donc au plus près des phénomènes étudiés. Le rotor
est à aimants permanents, dont la forme et les aimantations ont été optimisées de manière à
générer un champ d’induction quasiment sinusoı̈dal. Le très faible taux d’harmoniques de ce
champ magnétique permet de limiter les vibrations d’origine électromagnétiques [Far95]. La
deuxième machine est à courant continu, et permet d’entraı̂ner la machine de test si celle-ci
fonctionne en alternateur, ou de constituer sa charge si elle fonctionne en moteur. La figure
1.5 montre la structure générale de ce banc d’essais. L’accouplement, ainsi que les diverses
suspensions, ont pour rôle d’isoler la machine de test des vibrations du monde extérieur et de
celles engendrées par la machine à courant continu. Les mesures vibratoires effectuées sur la
7
CHAPITRE 1. Définition du problème traité
accouplement
souple
suspensions
machine synchrone ne prendront finalement en compte que ses propres vibrations, comme il
a été montré dans [Der92] et [Far95]. Toutefois, les fréquences de résonance des diverses sus-
pensions pouvant atteindre 16 Hertz [Far95], la fréquence de rotation des machines doit être
supérieure à 20 Hertz pour ne pas exciter une résonance du banc, et risquer sa détérioration.
L’objectif de ce travail, le principe physique utilisé pour l’atteindre, ainsi que la machine
de test sont maintenant connus. La section suivante cherche à mieux comprendre la nature des
vibrations dont on veut minimiser la transmission, c’est à dire celles des machines tournantes.
La vibration d’une structure, même complexe, est communément définie (et mesurée)
en un point de cette structure. Il s’agit alors du mouvement oscillatoire de ce point de la
structure autour d’une position d’équilibre, en réponse à une ou plusieurs forces excitatrices.
Or les grandeurs physiques usuelles permettant de caractériser un mouvement, et donc une
vibration, sont au nombre de trois :
– le déplacement,
– la vitesse,
– l’accélération.
De plus, le passage de l’une à l’autre s’effectue par simple dérivation ou intégration. On peut
donc utiliser indifféremment une de ces trois grandeurs comme mesure représentative d’une
vibration. Toutefois, l’accélération est celle qui permet de mettre en évidence les phénomènes
vibratoires dont les fréquences sont les plus élevées [Mor92, pages 85-86]. C’est donc cette
8
1.3. Vibrations d’une machine tournante électrique
grandeur qui est choisie pour représenter les vibrations de la machine synchrone, afin de
pouvoir caractériser les effets de l’isolation vibratoire sur la plage fréquentielle la plus large
possible.
Pour pouvoir effectuer ce type de mesure, nous disposons d’accéléromètres piézoélectriques
[Bou98, pages 53-57], dont la sortie fournit une quantité de charges électriques proportionnelle
à l’accélération subie. Cette information est transformée en tension par un amplificateur de
charge [Bou98, pages 59-61], dont la sortie est un signal analogique directement exploitable
par un oscilloscope, un analyseur de spectre, ou tout autre appareil de mesure. La bande
passante de ces capteurs s’étend de 5 Hz à 12 kHz, et nous verrons par la suite qu’elle est
largement suffisante pour notre application.
D’après la définition donnée à la section précédente, seules les forces excitatrices appliquées
à une structure sont à l’origine de ses vibrations. Le transfert entre ces forces et les vibrations
qu’elles engendrent est la réponse mécanique de la structure, illustrée par la figure 1.6.
Pour une machine tournante électrique, ces forces sont de trois natures différentes :
– aérodynamique,
– mécanique,
– électromagnétique.
La contribution de chacune de ces sources est fonction du type de machine étudiée, mais
quelques comportements généraux peuvent être soulignés.
Les forces d’origine aérodynamique proviennent de l’écoulement du fluide de refroidis-
sement dans les différentes parties du moteur, comme l’entrefer. Cet écoulement provoque
des variations aléatoires de pression dans le moteur, ce qui engendre des forces ayant sur-
tout un caractère large bande. Elles peuvent donc exciter les résonances mécaniques de la
machine, et génèrent des vibrations qui sont elles aussi aléatoires et large bande. Toutefois,
dans le cas de machines avec ventilateur de refroidissement, des vibrations sinusoı̈dales dont
la fréquence est le produit du nombre d’ailettes du ventilateur par la fréquence de rotation
peuvent apparaı̂tre. La machine de test décrite dans la section précédente étant refroidie
à l’air et ne comportant pas de ventilateur, ses vibrations d’origine aérodynamique ont un
caractère purement aléatoire large bande.
Les forces d’origine mécanique sont dues à des défauts de fabrication, de jeu ou d’usure.
On peut par exemple citer le phénomène de balourd dynamique provenant d’un mauvais
équilibrage du rotor, ou l’usure de l’accouplement reliant deux machines différentes. Ces
9
CHAPITRE 1. Définition du problème traité
phénomènes génèrent en général des forces périodiques, dont les fréquences sont liées à la
fréquence de rotation de la machine. Les vibrations qu’elles engendrent seront donc elles
aussi périodiques, et leur transformée de Fourier consistera en un ensemble de raies spec-
trales, de fréquences harmoniques ou sous-harmoniques de la fréquence de rotation. En ce qui
concerne la machine utilisée pour notre application, ces vibrations se limitent principalement
au balourd mécanique, générant des raies spectrales de fréquence égale à la fréquence de
rotation ainsi qu’à ses premiers harmoniques.
Les forces d’origine électromagnétique sont principalement des forces de Maxwell, donc
directement proportionnelles au carré du champ d’induction présent dans l’entrefer de la ma-
chine (cf. équation (1.1)). Ce champ n’est jamais à répartition purement sinusoı̈dale, et ceci
pour un grand nombre de raisons (distribution spatiale imparfaite des forces magnétomotrices
créées par les enroulements et/ou les aimants, entrefer non constant, courants d’alimentation
non sinusoı̈daux, etc). Les harmoniques de l’induction engendrent donc un grand nombre de
champs de force tournants à répartition périodique, qui excitent la structure de la machine. À
ces forces viennent aussi s’ajouter celles dues à des défauts de fabrication comme le balourd
magnétique, traduisant une excentricité du champ d’induction de la machine. Pour une ma-
chine synchrone, elles génèrent toutes des vibrations périodiques, de fréquences harmoniques
à la fréquence de rotation. Une étude très complète de ces différentes fréquences est faite
dans [Tim89, chapitre 3]. Pour la machine décrite au paragraphe 1.2.3, il a été montré [Far95]
que l’ensemble des raies spectrales les plus importantes provoquées par ce phénomène est
constitué de la raie d’encoche et de ses harmoniques. Elles ont pour origine la présence d’un
entrefer denté. En effet, la variation de perméance due aux encoches provoque une variation
de l’induction au droit de ces encoches. La machine de test ayant un rotor lisse et un stator à
36 encoches, la fréquence des raies vibratoires engendrées se situe autour des multiples entiers
de 36 fois la fréquence de rotation.
Bien que cette section soit relativement succincte, il faut être conscient de l’extrême
complexité du processus de génération des vibrations d’une machine tournante électrique. Ce
rapide exposé sur l’origine des vibrations d’une telle structure, largement inspiré du premier
chapitre de [Der92] traitant la question de manière plus approfondie, n’a pas la prétention
d’être exhaustif. Son but est seulement de donner une idée des différents types de vibrations
qui peuvent être engendrées dans une machine tournante, et des différents processus physiques
qui peuvent leur donner naissance. Pour la machine synchrone utilisée, on peut retenir que
les vibrations mesurées seront de deux types différents :
– les vibrations d’origine aérodynamiques auront un caractère aléatoire large bande,
– les vibrations d’origine mécanique ou électromagnétiques auront un caractère périodique,
donc bande étroite.
Maintenant que les caractéristiques des vibrations générées par la machine de test sont
connues, la section suivante permet de faire un choix parmi les nombreuses méthodes d’ana-
lyse existant pour un tel signal.
10
1.3. Vibrations d’une machine tournante électrique
Préliminaires
La plupart des méthodes d’analyse étant implantées sur calculateur numérique, le premier
conditionnement à réaliser sur un signal vibratoire est son échantillonnage. De plus, beaucoup
de ces méthodes font intervenir la notion de fréquence. Il est donc également nécessaire de
définir la transformation permettant de passer du domaine temporel au domaine fréquentiel.
Cette opération est réalisée, sur des signaux échantillonnés, par la transformation de Fou-
rier en fréquence réduite (TFr). De sa définition découle la condition d’échantillonnage de
Shannon, suivant laquelle l’échantillonnage d’un signal se fait sans perte d’information. Nous
allons voir comment la TFr et cette condition s’expriment pour des signaux mono ou bidi-
mensionnels, c’est à dire dépendant d’une ou de deux variables. En effet, ces derniers seront
nécessaires dans le chapitre 3, pour l’étude des systèmes variant périodiquement dans le
temps.
signaux monodimensionnels
Les TFr directe et inverse d’un signal discret s(n) sont respectivement définies par les
relations :
+∞
X
ŝ(α) = s(n)e−j2παn (1.2)
n=−∞
Z 1
2
s(n) = ŝ(α)ej2παn dα (1.3)
− 12
Pour que cette transformation existe, le signal étudié doit être de module sommable. En effet,
c’est une condition suffisante pour que la série exprimée en (1.2) converge absolument vers
une fonction continue [Kun91].
La fréquence réduite sans dimension α est liée à la fréquence réelle ν et à la fréquence
d’échantillonnage νe , exprimées toutes deux en Hertz, par :
¸ ¸
ν 1 1
α= ∈ − ; (1.4)
νe 2 2
11
CHAPITRE 1. Définition du problème traité
signaux bidimensionnels
De plus, les fréquences réduites α et β ont le même intervalle de définition qu’en (1.4). La
condition de Shannon bidimensionnelle impose donc au signal analogique original de ne pas
avoir de composantes fréquentielles supérieures à ν2e Hertz, et ceci suivant ses deux variables
fréquentielles [Tek95].
Les différentes méthodes d’analyse d’un signal vibratoire peuvent être classées en trois
catégories différentes.
La première regroupe les méthodes dites classiques, permettant surtout d’analyser des si-
gnaux stationnaires. Dans le domaine temporel, on retrouve par exemple l’analyse par in-
dicateurs, tels que l’amplitude maximale du signal, ou son facteur de crête [Big95, pages
157-161]. Des outils statistiques tels que l’histogramme, les moments d’ordre 1 ou 2, ou la
recherche de corrélations peuvent aussi être utilisés [Mor92, pages 131-133]. Dans le domaine
fréquentiel, la méthode la plus employée actuellement reste l’analyse spectrale, donnant accès
à la répartition de la puissance ou de l’énergie d’un signal dans le domaine fréquentiel. Elle
permet aussi de calculer les fonctions de transfert fréquentielles éventuelles reliant différents
signaux vibratoires. Enfin, une de ses variantes, l’analyse cepstrale introduite dans [Bog63],
permet de mettre facilement en évidence les modulations d’amplitude, même aux basses
fréquences [Bou98, pages 99-102].
Les méthodes classées dans la deuxième catégorie peuvent être qualifiées de modernes, et per-
mettent l’analyse de signaux non-stationnaires. Les principales sont les transformées temps-
fréquence (spectrogramme, Wigner-ville, scalogramme, etc) [Fla93, Whi00, Bas92]. Elles ont
été utilisées plus tardivement que les précédentes du fait de leur complexité théorique, et de
la puissance de calcul qu’elles nécessitent.
Enfin, la troisième catégorie regroupe les méthodes émergentes et novatrices, d’un usage en-
core peu courant dans l’industrie. Ces outils concernent par exemple les techniques d’analyse
cyclostationnaire [Gar91, Mac98, Ran00, Cap00] comme la corrélation spectrale, représentant
12
1.3. Vibrations d’une machine tournante électrique
les liens statistiques entre différentes fréquences d’un même signal. Récemment, des statis-
tiques d’ordre supérieur à deux ont aussi été utilisées [Lac97]. En effet, on a défini dans
le domaine temporel de nouveaux indicateurs, comme des moments d’ordre supérieur à 2
(skewness, kurtosis) qui permettent de mieux caractériser la loi de répartition de l’amplitude
d’un signal. Dans le domaine fréquentiel, l’emploi des multispectres (bispectre, trispectre,
etc) peut aider à la recherche de couplages non-linéraires [Gra99b]. On peut enfin noter l’em-
ploi de modélisations non-linéaires comme les filtres de Volterra, afin de mieux approcher les
transferts mécaniques complexes [Mio96].
Pour le présent problème d’isolation vibratoire, un outil permettant non seulement de
caractériser, mais aussi de quantifier les effets de la méthode employée est nécessaire. La
densité spectrale de puissance répond à ces exigences, puisqu’elle permettra de localiser les
fréquences des vibrations affectées, et de mesurer la variation de puissance (globale et/ou par
bande de fréquence) provoquée. Le prochain paragraphe expose donc les définitions nécessaires
à son utilisation.
13
CHAPITRE 1. Définition du problème traité
La grandeur obtenue a bien une dimension de puissance (énergie divisée par temps), et est
une fonction réelle de la fréquence. Il faut noter que la α est une fréquence réduite sans
dimension, dont l’intervalle de définition et le lien avec la fréquence réelle ν (en Hertz) sont
donnés par (1.4). En généralisant (1.6), on peut déduire la densité interspectrale de puissance
entre deux signaux discrets de puissance moyenne finie x(n) et y(n) :
· ¸
x̂N (α)ŷN (α)∗
γxy (α) = lim E (1.7)
N →∞ N
Sous les hypothèses conjointes de stationnarité (invariance temporelle des propriétés sta-
tistiques) et d’ergodisme (les moyennes temporelles tendent vers des valeurs finies et certaines)
que nous supposerons toujours vérifiées par les signaux vibratoires analysés, l’espérance
mathématique de la relation (1.6) peut être remplacée par une moyenne temporelle. On
obtient alors γs (α) en découpant le signal discret s(n) en M tranches de L échantillons, et
en effectuant la moyenne suivante :
M −1
1 X |ŝLi (α)|2
γs (α) = lim
M →+∞ M L
i=0
où ŝLi (α) est la transformée de Fourier de la ième tranche de L échantillons notée sLi (n).
Cette relation n’est toujours pas utilisable en pratique à cause du nombre infini de tranches.
On construit donc une estimée de la densité spectrale de puissance notée γ˜s (α) en prenant
M fini mais suffisamment élevé :
M −1
1 X |ŝLi (α)|2
γ˜s (α) = (1.8)
M L
i=0
Cette méthode est appelée estimation par périodogramme moyenné, et on peut montrer que
c’est un estimateur consistant de γs (α). Les différentes densités spectrales présentées dans la
suite de ce rapport ont été calculées par cette méthode.
D’après le paragraphe sur l’origine des vibrations d’une machine tournante électrique, les
forces excitant leur structure sont soit périodiques, soit aléatoires large bande. Les signaux
14
1.3. Vibrations d’une machine tournante électrique
vibratoires qu’elles peuvent engendrer ont donc les mêmes caractéristiques, et l’allure générale
de leur densité spectrale de puissance est montrée à la figure 1.7, duale de la figure 1.6, mais
en adoptant un point de vue fréquentiel.
spectres de puissance fonction de transfert spectre de puissance
des forces excitatrices de la réponse mécanique d’un signal de vibration
γa (α)
γ(α)
α
0 1
2
γp (α)
α
0 1 α
2 0 1
2
α
0 1
2
Fig. 1.7: Allure théorique de la densité spectrale de puissance d’un signal vibratoire de ma-
chine tournante électrique
La densité spectrale de puissance γp (α) des forces périodiques est un ensemble de raies
d’amplitude variable. Celle γa (α) des forces aléatoires a un caractère large bande. Par l’in-
termédiaire de la réponse mécanique dont la fonction de transfert présente plusieurs résonances
et anti-résonances, ces forces génèrent un signal vibratoire dont la densité spectrale de puis-
sance γ(α) présente à la fois des raies spectrales, et un caractère large bande. La section
suivante est dédiée à vérifier ce résultat en estimant la densité spectrale de puissance d’un
signal vibratoire mesuré sur la machine synchrone de test.
Domaine temporel
La figure 1.8 donne l’allure temporelle du signal de vibrations mesuré, ainsi que celle du si-
gnal de top-tour, dont une période représente une révolution complète du rotor de la machine.
15
CHAPITRE 1. Définition du problème traité
On peut noter visuellement la présence d’un déphasage constant entre le top-tour et certaines
800 vibration
top−tour
600
200
−200
−400
−600
400 400
200 200
0 0
−200 −200
−400 −400
−600 −600
16
1.3. Vibrations d’une machine tournante électrique
par :
N −1
1 X
PsN = |sN (k)|2 (1.9)
N
k=0
Dans le cas présent, on obtient, pour les composantes synchrones, une puissance 16 fois
supérieure à celle des composantes non synchrones. Ceci montre bien leur prédominance dans
ce signal vibratoire.
Le paragraphe suivant donne une représentation de ces signaux dans le domaine fréquentiel,
grâce à leurs densités spectrales de puissance.
Domaine fréquentiel
80
raie d’encoche
70 et ses harmoniques
60
decibels
50
40
30
20
résonances
importantes 10
0 500 1000 1500 2000 2500 3000 3500
frequence (Hertz)
large bande et de raies spectrales, comme prévu par la figure 1.7. On peut clairement iden-
tifier les raies spectrales dues aux balourds mécanique et magnétique, situées aux premières
fréquences harmoniques de νr . On trouve ensuite, autour de 36 × νr ≃ 857 Hertz, un groupe
de raies d’encoche, ainsi que leurs harmoniques autour des multiples entiers de cette valeur.
Pour le côté large bande, on peut noter la présence de deux résonances importantes (autour
de 450 et 1700 Hertz), contribuant à l’amplification de certaines raies spectrales, dont le pre-
mier harmonique de la raie d’encoche. Les densités spectrales de puissance des composantes
synchrones et non synchrones sont données respectivement par les figures 1.11(a) et 1.11(b).
On retrouve, comme principaux phénomènes synchrones, les raies de balourd ainsi que les
raies d’encoche, visibles jusqu’à 3400 Hz, et toutes situées à des fréquences harmoniques de νr .
Les phénomènes non synchrones sont, quant à eux, principalement à caractère large bande,
et couvrent toute la bande du filtre passe-bas utilisé dans la chaı̂ne d’acquisition. On note la
présence des deux résonances importantes repérées sur la densité spectrale (1.10), amplifiant
17
CHAPITRE 1. Définition du problème traité
80 80
70 70
60 60
decibels
50
decibels
50
40 40
30 30
20 20
10 10
0 500 1000 1500 2000 2500 3000 3500 0 500 1000 1500 2000 2500 3000 3500
frequence (Hertz) frequence (Hertz)
Cette section a donc permis de donner un rapide aperçu sur l’origine et la nature des
vibrations d’une machine tournante électrique, ainsi que sur les principales méthodes de
traitement du signal utilisées pour les analyser. De plus, un exemple de signal vibratoire
réel mesuré sur la machine de test a montré que pour cette machine, les vibrations radiales
prédominantes sont synchrones à la vitesse de rotation et que la méthode d’isolation vibratoire
envisagée doit atténuer en priorité leur transmission. Elles sont représentées, dans le domaine
fréquentiel, par un ensemble de raies spectrales situées aux fréquences harmoniques de νr .
Nous verrons au chapitre 4 que leur synchronisme et leur périodicité permettent de développer
des algorithmes beaucoup plus simples que pour des vibrations à caractère aléatoire. Ce point
est important dans le sens où la méthode élaborée doit être implantable dans un processeur
numérique, afin de pouvoir fonctionner en temps-réel. La nature des vibrations à atténuer
étant maintenant connue, le modèle du transfert grâce auquel on pourra agir sur celles-ci va
être déterminé dans la section suivante.
18
1.4. Modélisation physique du transfert
vn (t, θ)
19
CHAPITRE 1. Définition du problème traité
machine, quel que soit son état de fonctionnement. Le champ d’induction est donc toujours
proportionnel au champ magnétique, quelle que soit sa valeur.
Après ces premières simplifications, nous allons «progresser» le long du schéma de la
figure 1.12, en partant des courants de commande pour finalement en déduire l’expression de
vc (t, θ).
= bobine statorique 2
= bobine statorique 3
lié au stator fixant la position angulaire θ est l’axe de la bobine statorique numéro 1. En
supposant ces trois enroulements identiques et idéalement répartis, et en négligeant l’effet
de saillance dû aux faibles variations de l’entrefer lisse, le champ d’induction créé par les
courants de commande prend la valeur suivante :
· ¸
2π 4π
bc (t, θ) = K cos(pθ)i1c (t) + cos(pθ − )i2c (t) + cos(pθ − )i3c (t) (1.10)
3 3
20
1.4. Modélisation physique du transfert
la fois par le rotor et le stator. Afin de déterminer l’expression de b(t, θ), on se place dans le
cas où la machine synchrone de test est en fonctionnement normal. Elle est donc utilisée en
moteur alimenté par un système de courants triphasés équilibré, ou en alternateur débitant
sur une charge triphasée équilibrée, sans aucun courant de compensation parcourant ses
enroulements. Dans [Far95], la structure de cette machine a été optimisée de sorte que pour
un tel fonctionnement, son champ d’induction ait une répartition temporelle et fréquentielle
quasiment sinusoı̈dale. On peut alors exprimer b(t, θ) sous la forme :
b(t, θ) = B cos(2πpνr t − pθ + ϕ)
où l’expression de bc (t, θ) en fonction des courants de compensation est donnée par la relation
(1.10).
L’équation (1.1), donnant l’expression de la force due au champ d’induction de la machine
et subie par l’élément de surface statorique ds situé en θ peut s’écrire, en utilisant (1.11) :
ds
f (t, θ) = [B cos(2πpνr t − pθ + ϕ) + bc (t, θ)]2
2µ0
B 2 ds Bds ds
= cos2 (2πpνr t − pθ + ϕ) + cos(2πpνr t − pθ + ϕ)bc (t, θ) + bc (t, θ)2
2µ0 µ0 2µ0
Le premier terme de cette somme est engendré par le seul champ d’induction dû au fonction-
nement normal de la machine. Il provoque l’apparition d’un champ de force périodique de
fréquence 2pνr , et se retrouve dans tous les cas de figure, même en l’absence des courants de
compensation. Ce n’est pas le cas des deux derniers termes, qui apparaissent avec bc (t, θ), et
donc avec i1c (t), i2c (t) et i3c (t). Ils donnent l’expression de la force de compensation fc (t, θ)
créée par ces courants, que l’on peut mettre sous la forme :
Le transfert entre le champ de compensation et la force qu’il engendre est donc, dans le
cas général, à la fois non-linéaire et variant dans le temps. En effet, ce dernier caractère est
contenu dans le premier terme de la somme (1.12), où l’entrée (bc (t, θ)) est multipliée par
un gain fonction du temps. Quant à la non-linéarité, elle provient du second terme de cette
somme qui est quadratique.
21
CHAPITRE 1. Définition du problème traité
Une dernière hypothèse, portant sur les courants de compensation, va nous permettre de
simplifier cette relation. En effet, le rôle de ces courants est d’atténuer les vibrations de
la machine de test, mais ils ne doivent en aucun cas perturber son fonctionnement. En ce
sens, ils doivent rester faibles devant les courants statoriques présents lorsque la machine est
utilisée normalement. Il en découle que le champ de compensation bc (t, θ) qu’ils engendrent
sera petit devant celui de la machine en fonctionnement normal. Donc pour des courants de
compensation assez faibles, le terme quadratique de la relation (1.12) peut être négligé devant
l’autre terme, d’où l’approximation suivante :
Le dernier transfert qu’il reste à déterminer sur la figure 1.12, relie la force de com-
pensation appliquée à la carcasse statorique en θ, à la vibration additionnelle vc (t, θ) que
cette dernière donne en réponse. Il s’agit donc de la réponse mécanique du stator, qui est
fonction de la position angulaire considérée. On imagine bien la complexité d’un tel trans-
fert, c’est pourtant sur lui que nous allons faire l’hypothèse la plus forte : il est supposé
linéaire et invariant dans le temps. Il peut donc être entièrement modélisé par une fonction
de transfert fréquentielle au sens classique du terme, qui est à la fois fonction de la fréquence
réelle ν et de la position angulaire θ et notée ĝ(ν, θ). Cette approche a déjà été adoptée
dans [Mae87], [Der92] et [Far95], avec de bons résultats de modélisation. De plus, quelques
éléments de réponse peuvent être trouvés dans [Ful96, chapitre2], où l’auteur montre qu’un
cylindre de longueur finie répond à divers types de champ de force sinusoı̈daux en temps, par
un déplacement de même fréquence. Ceci conforte donc l’hypothèse de linéarité et d’inva-
riance dans le temps pour la réponse mécanique du stator, puisque sa structure est bien celle
d’un cylindre de longueur finie. Si son comportement s’avérait malgré tout plus compliqué,
ĝ(ν, θ) permettrait tout de même de modéliser la partie linéaire de sa réponse mécanique.
De plus, on rappelle qu’en pratique, les signaux vibratoires sont fournis par des accéléro-
mètres, dont le transfert est linéaire et invariant dans le temps dans la bande fréquentielle
considérée. Le transfert de ce capteur devrait donc apparaı̂tre sur la figure 1.12 après le
signal vibratoire v(t, θ). On peut en fait considérer que la fonction de transfert fréquentielle
du stator englobe celle du capteur, puisqu’elles sont de même type. ĝ(ν, θ) servira donc de
modèle à la fois à la carcasse statorique et à l’accéléromètre.
22
1.4. Modélisation physique du transfert
Tous les transferts de la figure 1.12 ont donc été déterminés. En utilisant les équations
(1.10) et (1.13) ainsi que les remarques du paragraphe précédent, on peut représenter le
modèle du système global comme sur la figure 1.14.
cos(pθ)
ic1 (t)
2π
m(t, θ) = F cos(2πpνr t − pθ + ϕ) vn (t, θ)
cos(pθ − 3 )
ic3 (t)
Le système situé entre les courants de commande et la vibration additionnelle qu’ils en-
gendrent est bien linéaire. Toutefois, la présence du signal sinusoı̈dal modulant en amplitude
le champ d’induction bc (t, θ) fait varier ses caractéristiques périodiquement dans le temps.
Le terme «caractéristiques» désigne ici les grandeurs représentatives du système, comme sa
réponse impulsionnelle ou sa fonction de transfert fréquentielle. Le système global peut donc
être vu comme un système linéaire variant périodiquement dans le temps (LVPT), à trois
entrées et une sortie. Il est à noter que ce transfert a été validé grâce à des signaux réels dans
[Gra98] et [Gra99b] à l’aide des statistiques d’ordre supérieur (bispectre et bicohérence). Il
est sera donc considéré comme parfaitement valable, dans la suite de ce mémoire.
Enfin, on peut remarquer que chaque courant de commande ne peut agir sur toute la cir-
conférence de la carcasse statorique. En effet, si on se limite au courant ic1 (t), le gain cos(pθ)
qui lui est appliqué en entrée lui interdit de pouvoir engendrer une quelconque vibration
π
additionnelle pour des positions angulaires vérifiant θ = (2k + 1) 2p , où k ∈ {0, 1, . . . , 2p − 1}.
Par exemple, pour la machine de test le nombre de paires de pôles p étant égal à deux, cette
entrée ne peut engendrer de vibrations aux positions angulaires π4 , 3π 5π 7π
4 , 4 et 4 radians. Par
contre, les autres entrées peuvent prendre le relais puisque leurs gains respectifs ne sont alors
pas nuls. On peut donc d’ores et déjà voir qu’en se limitant à un système mono-entrée, on ne
peut avoir d’action sur tout le stator, contrairement à un système multi-entrées.
23
CHAPITRE 1. Définition du problème traité
1.5 Conclusion
Ce premier chapitre a permis de bien spécifier le sujet du travail de recherche présenté
dans ce mémoire.
La première section définit tout d’abord l’objectif général de ce travail : réaliser l’isolation
vibratoire de la machine synchrone de test présentée au paragraphe 1.2.3. Il a ensuite été
montré qu’une manière d’atteindre cet objectif est de minimiser les vibrations de la carcasse
statorique autour d’une position angulaire, et de se servir de cette «zone de calme» comme
zone de fixation. On a vu enfin que le principe d’action sur les vibrations de la machine est
l’utilisation des forces radiales d’origine magnétique, générées par des courants de compen-
sation dans les enroulements de cette machine.
Dans la seconde section, il a été montré que les vibrations synchrones au phénomène de ro-
tation (donc périodiques et de fréquence harmonique de la fréquence de rotation νr ) sont
prédominantes dans les vibrations radiales de la machine de test. La méthode à élaborer doit
donc chercher à les atténuer en priorité, et doit pouvoir couvrir tout le spectre vibratoire
fréquentiel jusqu’à 3 kHz.
Enfin, la modélisation du système mis en jeu, réalisée dans la troisième et dernière section, à
mis en évidence qu’il pouvait être considéré comme linéaire et variant périodiquement dans le
temps, sous l’hypothèse de courants de compensation assez faibles. Le but à atteindre est donc
d’atténuer «au mieux» un signal formé de composantes spectrales de fréquences harmoniques
de νr , à travers un système linéaire et variant périodiquement dans le temps de fréquence pνr
connue, tout en gardant la possibilité de limiter le niveau des signaux de commande.
Avant d’essayer de résoudre ce problème, nous allons nous pencher, dans le prochain
chapitre, sur l’étude des travaux qui ont déjà été effectués sur ce sujet.
24
Chapitre 2
2.1 Introduction
Ce chapitre a pour but de présenter les méthodes «classiques» employées dans le domaine
du contrôle actif de vibrations, et de montrer leur insuffisance pour la présente application.
En effet, ces méthodes supposent la présence de transferts linéaires et invariants dans le
temps au sein du système global de contrôle. Elles ne sont donc pas adaptées au système
linéaire et variant périodiquement dans le temps dont le modèle a été déterminé au chapitre
précédent. Nous verrons que cette particularité diminue fortement leurs performances, et
qu’elles deviennent alors difficilement exploitables.
La section 2.2 est une introduction au concept de base du contrôle actif, et décrit les
différentes stratégies d’implantation ainsi que les algorithmes les plus couramment utilisés.
Les résultats obtenus par ce type de méthodes dans le cadre de la compensation de
vibrations envisagée sont ensuite présentés dans la section 2.3. Les divers travaux de recherche
dont ce sujet a précédemment fait l’objet, ainsi que les résultats obtenus y sont analysés.
2.2.1 Généralités
Historiquement, les techniques de contrôle actif ont été élaborées afin de lutter contre des
ondes sonores indésirables [Lue36, Ols53, Con56]. La propriété fondamentale qu’elles utilisent
est le principe de superposition. En effet, une source appelée «source primaire» est supposée
générer un bruit que l’on veut éliminer. Un système de contrôle actif est chargé d’émettre,
par l’intermédiaire d’une «source secondaire», un son de même amplitude que ce bruit, mais
de phase opposée. Idéalement, ces deux ondes s’ajoutent l’une à l’autre et forment ainsi une
erreur résiduelle nulle, comme indiqué sur la figure 2.1. Toutefois, les erreurs modélisation
pouvant se produire sur les différents transferts présents dans le système, ainsi que le caractère
25
CHAPITRE 2. Méthodes classiques de contrôle actif
bruit à éliminer
source primaire
de bruit
✂✁✁✁
capteur
✂✁✂
✂✁✁
✁✁✁✂✁✂ erreur résiduelle
✂✁✂✁✂ nulle
source secondaire
souvent aléatoire des signaux de perturbation, font que dans le cas général, on ne peut annuler
complètement l’erreur résiduelle. Par contre, on peut chercher à l’atténuer «au mieux» en
calculant le bruit secondaire «optimal» qui minimise une fonction de coût dépendant de
l’erreur résiduelle. Ce coût est choisi quadratique dans l’immense majorité des cas, d’une part
parce que cette grandeur, homogène à une puissance sonore, caractérise bien les performances
d’atténuation du système de contrôle, et d’autre part parce que ce type de coût permet une
grande simplification des calculs. Il est à remarquer que le bruit élaboré par le système de
contrôle actif ne sera optimal qu’au sens de cet indice de performance.
Plus généralement, le principe du contrôle actif peut être appliqué à divers types de per-
turbations, dont les vibrations. La figure 2.2 reprend la figure 2.1, mais dans un cadre plus
général, où les perturbations peuvent être d’une nature physique quelconque. La perturbation
commande i
à élaborer source
secondaire
générée uniquement par la source primaire au niveau du capteur est notée d. Celle générée
au même point par la source secondaire, et donc par la commande i, est notée v, et peut être
considérée comme une «contre-perturbation». Le système de contrôle actif doit être capable
d’élaborer la commande iopt , optimale au sens où elle doit minimiser un critère quadratique
dépendant de l’erreur e = d + v. Une fois ce minimum atteint, on peut dire que la puissance
de l’erreur résiduelle a été diminuée, et la perturbation compensée grâce à l’action de la
26
2.2. Introduction au contrôle actif classique
source primaire x
P
de perturbations
d
v
i S e
On définit tout d’abord un signal x, que l’on nomme «signal de référence». Il est dû exclusive-
ment à la source primaire, et permet de modéliser le trajet parcouru par la grandeur physique
perturbatrice de cette source au capteur. En effet, le «transfert primaire» noté P, relie la
référence x à la perturbation mesurée d. Il est à noter que le signal de référence, s’il est acces-
sible à la mesure, permettra au système de contrôle actif de prédire certaines caractéristiques
des perturbations à éliminer au niveau du capteur. Ensuite, le transfert donnant le lien entre
le signal de commande i de la source secondaire, et la contre-perturbation v qu’elle engendre
au niveau du capteur, est noté S. Finalement, on voit qu’on ne peut atténuer d que par
l’intermédiaire de ce transfert secondaire. Il est donc naturel de penser que les performances
du système de contrôle actif seront très dépendantes des caractéristiques de S.
Afin d’illustrer les différents signaux et systèmes de la figure 2.3, nous allons décrire
quelques exemples applicatifs de contrôle actif de vibrations. Le premier exemple, dont le
schéma de principe est donné par la figure 2.4, concerne le contrôle actif des vibrations d’une
poutrelle métallique de longueur infinie, soumise à des ondes de flexion qui se propagent de la
gauche vers la droite. Un premier accéléromètre fournit le signal de référence x, représentant
capteur
=
accéléromètre
x e
onde de flexion
✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁
✂✁ ✂✁
✂✁✂✁
✂✁✂✁
✂✁✂✁
✂✁✂✁
✂✁✂✁
✂✁✂✁
✂✁✂✁
✂✁✂✁
✂✁✂✁
✂✁✂✁
✂✁✂✁
✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✁✂✂
✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁✁
✂✁
source
secondaire
poutrelle métallique
=
vérin i
les vibrations de la poutrelle en une position amont par rapport à un second accéléromètre.
Ce dernier mesure le signal d’erreur e à minimiser, ou, en l’absence de signal de commande i,
27
CHAPITRE 2. Méthodes classiques de contrôle actif
les vibrations d à éliminer. Le transfert P de la figure 2.3 est donc celui reliant les signaux x
et e, en l’absence de commande i. Cette dernière agit additivement sur d en créant la contre-
vibration v grâce à une source secondaire, constituée par un vérin. Le transfert S relie quant
à lui la commande i du vérin, au signal d’erreur e, en l’absence des perturbations engendrées
par les ondes de flexion. Pour cette application, le signal de référence peut être considéré
comme une version des perturbations qui vont affecter le second accéléromètre, et peut donc
être utilisé afin de les prédire.
Le second exemple auquel on peut s’intéresser est celui de l’application décrite au chapitre
précédent, concernant le contrôle actif des vibrations d’une machine tournante électrique.
Dans ce cas, la difficulté revient à déterminer un signal de référence représentatif des vibra-
tions à éliminer, dont la seule source est la rotation de la machine. Intuitivement, on peut
penser que le signal de sortie d’un capteur angulaire fixé sur l’arbre de rotation (signal de
top-tour) peut remplir cette condition. En effet, il fournit une sinusoı̈de dont la fréquence est
celle de rotation νr . Ce point est confirmé par les résultats obtenus à la section 1.3.4, où il a
été montré que les vibrations de la machine de test étaient principalement constituées d’un
ensemble de sinusoı̈des, dont les fréquences sont les harmoniques de νr . Dans la suite, c’est
donc ce choix qui sera fait pour le signal de référence, qui fournira le fondamental des vibra-
tions périodiques à éliminer. Le schéma correspondant à cette configuration est celui de la
figure 2.5. Le signal d’erreur est mesuré par un accéléromètre placé sur la carcasse statorique
i1 , i2 , i3
capteur angulaire
Fig. 2.5: Contrôle actif des vibrations d’une machine tournante électrique
de la machine, alors que le signal de référence est fourni par un capteur angulaire placé en
bout d’arbre rotorique. Les signaux de commande sont au nombre de trois, et correspondent
aux courants injectés dans chacune des bobines statoriques de la machine. Le transfert P
relie le signal de top-tour x aux vibrations à éliminer d, que l’on peut mesurer en sortie de
l’accéléromètre lorsqu’il n’y a aucun signal de commande. Le transfert secondaire S reliant
les signaux de commande aux contre-vibrations qu’elles engendrent a déjà été étudié à la sec-
tion 1.4, et sa structure est représentée à la figure 1.14. Cette fois, la référence ne représente
pas une version future des perturbations à éliminer, mais permet de déterminer certaines de
leurs caractéristiques, comme leurs fréquences.
Il existe beaucoup d’autres applications qui ont été étudiées dans le cadre du contrôle
28
2.2. Introduction au contrôle actif classique
actif des vibrations (voir en particulier les exemples donnés dans [Ber97, Ful95, Ful96]). On
peut néanmoins noter qu’elles reposent toutes sur le modèle de la figure 2.3.
Après avoir établi le modèle général d’un système de contrôle actif, nous allons déterminer
celui du contrôleur à utiliser, qui est fonction de la disponibilité d’un signal de référence x
représentatif des vibrations à éliminer.
Stratégie «feedback»
Dans un premier temps, on suppose ne pas disposer d’un signal de référence correct. Ceci
peut arriver par exemple si on ne possède pas le capteur correspondant (le capteur angulaire
dans notre cas de figure), ou si ce signal est physiquement impossible à mesurer (le bruit am-
biant dans une pièce contenant de nombreuses personnes, par exemple). La seule information
disponible sur la perturbation à éliminer et sur l’atténuation réalisée par le système est donc
le signal d’erreur. Il constitue la seule entrée du contrôleur, alors que sa sortie est le signal de
commande à élaborer. On obtient une structure à boucle de rétroaction, plus communément
dénommée «structure feedback», dont le principe et les transferts mis en jeu sont représentés
sur la figure 2.6, et où le transfert du contrôleur est noté C. Cette structure a été intensément
étudiée en automatique, et une multitude de méthodes ont été élaborées pour trouver un
correcteur satisfaisant, qui ne seront pourtant pas utilisées par la suite. En effet, la plupart
d’entre elles supposent le transfert S linéaire et invariant dans le temps alors qu’il a été
montré linéaire et variant périodiquement dans le temps (voir section 1.4). D’autre part, elles
imposent toutes d’étudier la stabilité du système bouclé, ce qui peut devenir difficile avec un
tel transfert. Enfin, nous disposons d’un signal de référence fourni par le capteur angulaire
décrit précédemment. Comme nous allons le voir dans le paragraphe suivant, la présence de
ce signal permet d’adopter une structure différente pour le contrôleur.
29
CHAPITRE 2. Méthodes classiques de contrôle actif
Stratégie «feedforward»
capteur
source primaire
de perturbations e
source secondaire
i
contrôleur
(a) principe
source primaire x
P
de perturbations
d
i v
C S
e
référence constitue l’entrée de ce contrôleur, qui utilise les informations qu’il contient sur les
perturbations à éliminer pour élaborer la meilleure commande i. Le signal d’erreur, quant
à lui, est utilisé pour mesurer les performances atteintes par le système de contrôle actif,
et modifier les caractéristiques du contrôleur si elles ne sont pas satisfaisantes. Grâce à ce
principe, le contrôleur et la commande i, vont converger vers leur valeur optimale, et la
fonction de coût vers son minimum.
D’après la figure 2.5, nous disposons, pour réaliser le système de contrôle actif désiré, d’un
signal de référence grâce au capteur angulaire monté sur l’arbre de rotation de la machine. Le
contrôleur utilisé aura donc une structure feedforward, et le prochain paragraphe a pour but
30
2.2. Introduction au contrôle actif classique
Cas général
On peut donner, dans le cas général de la figure 2.7(b), l’expression du transfert optimal
Copt du contrôleur feedforward ne vérifiant aucune contrainte de causalité ni de stabilité. Pour
cela, il faut exprimer le signal d’erreur e en fonction des signaux et des transferts du système
global. Dans un premier temps, les transferts P et S ne sont pas supposés invariants dans le
temps, et le symbole mathématique «◦» se lisant «appliqué à» doit être utilisé pour relier
leur entrée à leur sortie. Par exemple, le signal de commande i est la sortie du transfert C
d’entrée x, ce qui se note :
i=C◦x (2.1)
e = d + v = P ◦ x + S ◦ C ◦ x = (P + S ◦ C) ◦ x (2.2)
Pour annuler l’erreur, le transfert du contrôleur optimal doit donc vérifier l’équation :
S ◦ Copt = −P
On peut alors définir Sinv , le transfert inverse de S qui vérifie, pour un signal y quelconque :
Cette équation montre qu’il est nécessaire d’inverser le transfert secondaire S pour obtenir
une compensation parfaite de la perturbation d. Les caractéristiques de ce transfert auront
donc une grande influence sur les performances obtenues par le système de contrôle actif
feedforward, ce qui peut même poser des problèmes de réalisation si Sinv est instable ou
non-causal.
Si les transferts P et S sont linéaires et invariants dans le temps (LIT), ils sont entièrement
décrits par leur réponse impulsionnelle, ou par sa transformée de Fourier. Ce point sera étudié
31
CHAPITRE 2. Méthodes classiques de contrôle actif
plus en détail dans le chapitre 3, portant sur l’étude théorique des systèmes linéaires variant
périodiquement dans le temps. On peut toutefois définir dès maintenant les quelques notations
qui nous serons nécessaires dans ce chapitre, en prenant l’exemple du transfert C. Sa réponse
impulsionnelle est notée c(r), où r représente le retard en nombre d’échantillons par rapport
à celui de l’impulsion d’excitation. Dans le domaine temporel, le signal de sortie est alors
donné par la convolution du signal d’entrée par c(r). Pour des signaux représentés sous leur
forme temporelle, on a donc une relation entrée-sortie équivalente à (2.1), où l’opérateur «◦»
devient l’opérateur de convolution. La transformée de Fourier ĉ(α) de c(r) est appelée «gain
complexe» ou «fonction de transfert en fréquence» du système C. Elle permet quant à elle de
donner la relation entrée-sortie du système dans le domaine fréquentiel, où l’opérateur «◦»
devient alors un simple produit. Pour simplifier les calculs qui vont suivrent et faciliter la
compréhension des résultats, on suppose que tous les signaux mis en jeu dans la figure 2.7(b)
sont déterministes et de module sommable. Dans ce cas, on peut calculer leur transformée de
Fourier et la relation générale (2.1) transposée en fréquence devient :
î(α) = ĉ(α)x̂(α)
où le domaine de définition de la fréquence réduite α est donnée par (1.4), et la définition de
la transformée de Fourier par (1.2). Les fonctions complexes î(α) et x̂(α) sont respectivement
les transformées de Fourier des signaux i et x.
En appliquant ce résultat à la relation générale (2.2), on obtient l’expression du signal d’erreur
en fréquence :
ˆ + ĉ(α)ŝ(α)x̂(α)
ê(α) = d(α) (2.5)
À partir de cette équation, on peut définir la fonction de coût que doit minimiser le
contrôleur optimal. Le choix de son expression est primordial, car elle doit représenter à elle
seule les performances que l’on attend du système de contrôle actif. Elle peut prendre une
forme quelconque, mais sa minimisation sera rendue infiniment plus simple si elle ne comporte
pas de minima locaux. Une forme très largement utilisée et vérifiant cette propriété est la
fonction de coût des moindres carrés, ou coût quadratique. Sa forme la plus simple dans le
domaine fréquentiel est le module carré de la transformée de Fourier du signal d’erreur :
C = |ê(α)|2 (2.6)
32
2.2. Introduction au contrôle actif classique
C est donc bien une fonction quadratique du gain complexe du contrôleur ĉ(α). Si on décompose
ce gain complexe en sa partie réelle ℜ [ĉ(α)] et sa partie imaginaire ℑ [ĉ(α)], la fonction C
peut alors être vue comme une fonction réelle de ces deux variables réelles. Tracée dans le
plan orthonormé dont l’axe des abscisses représente ℜ [ĉ(α)] et celui des ordonnées ℑ [ĉ(α)],
33
CHAPITRE 2. Méthodes classiques de contrôle actif
elle forme un paraboloı̈de. Cette surface est convexe, et comporte donc un unique minimum
seulement si la valeur réelle |ŝ(α)|2 |x̂(α)|2 + ǫ est positive, ce qui peut toujours être imposé
grâce à ǫ. La figure 2.8 montre un exemple d’une telle surface.
500
400
valeur de 300
C
200
100
0
1
0 1
ℑ [ĉ(α)] −1 0
−1 −2
ℜ [ĉ(α)]
Les coordonnées du point minimum sont les parties réelles et imaginaires du gain complexe
ĉopt (α) du contrôleur optimal Copt . Son expression analytique peut être facilement déterminée
en remarquant que pour ĉ(α) = ĉopt (α), le gradient de C s’annule dans toutes les directions.
On peut alors définir un gradient complexe ∇C dont le module et l’argument donnent respec-
tivement, au point ĉ(α) de coordonnées (ℜ [ĉ(α)] , ℑ [ĉ(α)]), la valeur de la plus grande pente
de C ainsi que sa direction [Nel92, pp 416-420] :
∂C ∂C
∇C = ∇C [ĉ(α)] = +j (2.10)
∂ℜ [ĉ(α)] ∂ℑ [ĉ(α)]
Le gain complexe optimal ĉopt (α) minimisant la fonction de coût (2.7) est donc la valeur
complexe annulant le gradient ci-dessus. En calculant (2.10) dans le cas du coût quadratique
régularisé (2.9), on obtient :
£ ¤
ˆ
∇C = 2ŝ(α)∗ d(α)x̂(α)∗
+ 2ĉ(α) |ŝ(α)|2 |x̂(α)|2 + ǫ (2.11)
Cette solution est souvent appelée «solution de Wiener», et on peut noter qu’elle est direc-
tement proportionnelle à l’interspectre entre les signaux d et x, donc à leur intercorrélation.
Reportée dans l’équation (2.9), cette expression donne la valeur minimale que peut prendre
34
2.2. Introduction au contrôle actif classique
la fonction de coût :
ˆ
|ŝ(α)|2 |x̂(α)|2 |d(α)|2
ˆ
C [ĉopt (α)] = Cmin = |d(α)|2
− 2 2
|ŝ(α)| |x̂(α)| + ǫ
ˆ
ǫ|d(α)| 2
= (2.13)
|ŝ(α)| |x̂(α)|2 + ǫ
2
p̂(α)
copt
ˆ (α) = − (2.14)
ŝ(α)
Cette expression, qui est un cas particulier de l’équation générale (2.4), est donc celle du
contrôleur annulant complètement l’énergie du signal d’erreur à la fréquence α. Toutefois, il
n’y a alors aucune contrainte sur le signal de commande, dont l’énergie peut atteindre des
valeurs très élevées. De plus, si la fréquence à laquelle on cherche à minimiser l’énergie de
l’erreur est un zéro du transfert secondaire (c’est à dire ŝ(α) = 0), copt ˆ (α) ne peut alors
pas être calculé. Par contre, si ǫ est choisie positive et différente de zéro, non seulement le
gain complexe optimal existe quel que soit ŝ(α), mais elle fournit aussi un moyen de limiter
l’énergie de la commande. Malgré tout, la relation (2.13) montre que la valeur minimum de
la fonction de coût croı̂t avec ǫ, jusqu’à ce que le contrôleur ait un gain complexe nul et que
ˆ
le coût reste figé à |d(α)|2 . Le réglage de la valeur de ǫ est donc une tâche très délicate, où il
faut déterminer le juste milieu entre un coût minimal suffisamment faible et une commande
optimale d’énergie peu importante.
L’expression du correcteur optimal permettant de minimiser la fonction de coût qua-
dratique régularisée (2.7), ainsi que ses performances sont donc connues dans le cas où les
transferts primaire et secondaire sont LIT. Tout le problème revient maintenant à déterminer
la manière de faire converger le contrôleur C vers cette valeur. Le paragraphe suivant donne
le principe d’un algorithme classiquement utilisé en contrôle actif pour réaliser cette tâche.
35
CHAPITRE 2. Méthodes classiques de contrôle actif
cas d’un changement de son régime de fonctionnement, d’un désserrement progressif de son
socle de fixation, d’un changement du couple de charge, etc.
L’algorithme recherché doit minimiser une fonction de coût quadratique C [ĉ(α)], ne com-
portant qu’un seul minimum. Dans ce cas, on peut utiliser l’algorithme du gradient pour
rechercher la valeur optimale de ĉ(α), dont la forme générale est :
où ĉ(α)l est la valeur du gain complexe à l’itération courante l, et ĉ(α)l+1 sa prédiction à
l’itération suivante.
La valeur de ĉ(α)l+1 est donc déduite de ĉ(α)l en lui enlevant une fraction du gradient de
la fonction de coût calculé au point ĉ(α)l . Nous avons vu précédemment que ce gradient
donne la valeur et la direction de la plus grande pente de la surface C au point considéré.
Le gain complexe va donc «descendre» progressivement le long de cette courbe, de sa valeur
originale fixée par l’utilisateur à sa valeur optimale ĉopt (α), en suivant le chemin de la plus
grande pente. Le coefficient réel positif µ, appelé «pas d’adaptation» et fixé par l’utilisateur,
quantifie l’importance de la réactualisation, et donc le temps de convergence de cet algorithme.
Afin d’obtenir une relation itérative utilisable en pratique, on peut écrire le gradient (2.11)
de la fonction de coût choisie sous la forme :
La première propriété à étudier est celle de la convergence. Pour cela, on définit une
variable El , traduisant l’erreur entre le gain complexe courant et sa valeur optimale :
Ce vecteur utilisé conjointement avec les équations (2.12) et (2.5) permet d’écrire la relation
(2.17) sous la forme :
£ ¡ ¢¤ £ ¡ ¢¤l
El+1 = 1 − 2µ |ŝ(α)|2 |x̂(α)|2 + ǫ El = 1 − 2µ |ŝ(α)|2 |x̂(α)|2 + ǫ E0
Si cette valeur tend vers 0 lorsque l tend vers l’infini, l’algorithme est dit convergent. Cette
condition est réalisée quelle que soit la valeur d’initialisation E0 , si et seulement si |1 −
¡ ¢
2µ |ŝ(α)|2 |x̂(α)|2 + ǫ | < 1, ce qui est équivalent à :
1
0<µ< (2.18)
ŝ(α)|2 |x̂(α)|2 + ǫ
Donc quelle que soit la valeur d’initialisation de l’algorithme (2.17), si le pas d’adaptation
vérifie cette relation, l’algorithme convergera vers le gain complexe optimal recherché.
36
2.3. Résultats antérieurs
Pour implanter d’une manière pratique l’algorithme précédent, il est nécessaire de réaliser
des traitements par blocs. En effet, afin de connaı̂tre les transformées de Fourier des signaux
d’erreur e et de référence x directement mesurés par les deux capteurs présents dans le
système de contrôle, la solution la plus simple est d’employer un algorithme de transformée
de Fourier rapide, portant sur des blocs d’échantillons. Ce type d’implantation permet de
lever l’hypothèse très contraignante suivant laquelle les signaux utilisés doivent être de mo-
dule sommable pour assurer l’existence des transformées de Fourier. Cette condition est alors
vérifiée quelle que soit leur nature puisqu’ils sont observés sur un intervalle de temps fini
correspondant à la longueur d’un bloc. Il est aussi nécessaire de connaı̂tre la valeur du gain
complexe ŝ(α) du transfert secondaire, qui devra être estimé durant une étape préliminaire
à l’étape de contrôle. En surmontant les grandeurs estimées d’un tilde (e), on obtient l’algo-
rithme pratique suivant :
g ∗ x̂(α)
ĉ(α)l+1 = (1 − 2µǫ)ĉ(α)l − 2µŝ(α) g ∗ ê(α)
g (2.19)
l l
Contrairement à (2.17), ce n’est plus le gradient théorique qui est utilisé dans cette relation
mais son estimée, d’où le nom d’«algorithme du gradient stochastique». Il est aussi appelé
«algorithme des moindres carrés x filtré pondéré» [Ell87], et constitue une variante de «l’al-
gorithme des moindres carrés x filtré» (FxLMS), originellement développé dans le domaine
temporel [Wid85, Chap. 11] pour le contrôle actif.
Cet algorithme fréquentiel a été généralisé au cas multi-entrées multi-sorties dans [Ell92],
et employé avec succès pour le contrôle de perturbations périodiques. Il faut toutefois remar-
quer que le traitement par bloc génère un retard dans le signal de commande proportionnel
au nombre d’échantillons contenus dans un bloc, et qu’il ne peut être employé dans le cas
de perturbations aléatoires. Pour ce type de perturbations, un algorithme récursif temporel
réactualisant le contrôleur échantillon par échantillon est préférable. Ces algorithmes sont la
plupart du temps dérivés de celui du gradient, et sont l’équivalent temporel de l’équation
(2.19). Ils ne sont pas exposés dans ce mémoire, le but étant de donner leur principe général
étudié à l’aide de l’exemple précédent, et non une bibliothèque de traitements déjà existante
dans plusieurs ouvrages [Nel92, Ful96, Kuo96].
Après avoir présenté les méthodes courantes de contrôle actif, la prochaine section décrit
les résultats qu’elles ont permis d’obtenir dans le cadre de l’application décrite au chapitre
précédent.
37
CHAPITRE 2. Méthodes classiques de contrôle actif
Le premier contrôleur issu des travaux décrits dans [Der92, Far95] est analogique. Il a
pour but de minimiser les vibrations du stator à une seule fréquence et pour une seule
position angulaire. Il s’agit d’un contrôleur de structure feedback, dont l’entrée est le signal
de vibrations statoriques totales, correspondant à l’erreur. Son principe est décrit sur la
figure 2.10, où les spectres des différents signaux sont schématisés en gris. Cette figure est à
rapprocher de la figure 2.6, décrivant un tel système dans le cas général.
En sortie de l’accéléromètre mesurant les vibrations totales de la carcasse statorique, la
composante spectrale vibratoire de fréquence νd est isolée par un filtre passe-bande sélectif.
Ce signal bande étroite est ensuite modulé en amplitude par une sinusoı̈de de fréquence pνr ,
ce qui le décale aux fréquences νd ±pνr . La fréquence de cette modulation est tirée du signal de
top-tour, fourni par un capteur angulaire monté en bout d’arbre rotorique. Un filtre réjecteur
rejette alors la composante de plus petite fréquence. Il reste donc, à l’entrée du contrôleur,
une seule composante bande étroite, de fréquence νd + pνr . Le contrôleur est linéaire et
invariant dans le temps, et applique une amplification et un déphasage sur son entrée, que
l’on peut régler manuellement. On a donc une composante de fréquence νd + pνr en entrée de
l’enroulement de compensation qui, une fois modulée et filtrée par le transfert de la machine,
génère des contre-vibrations à deux fréquences différentes : νd et νd + 2pνr . Le déphasage et
l’amplification imposée par le contrôleur permettent de minimiser la composante de fréquence
38
2.3. Résultats antérieurs
MACHINE SYNCHRONE
vibrations
statoriques
F cos(2πpνr t + φ) naturelles
vibrations
statoriques
νd + pνr νd νd + 2pνr système νd νd + 2pνr totales
LIT
cos(2πpνr t)
contrôleur
analogique νd − pνr νd
CONTRÔLEUR FEEDBACK
Dans un second temps, une version numérique adaptative de ce contrôleur a été réalisée
dans le but d’augmenter le nombre de composantes spectrales vibratoires touchées par ce
traitement, et de suivre les éventuelles variations de vitesse de la machine [Jar95]. Une étape
importante a alors été franchie, puisque le contrôleur prend une structure feedforward afin
d’éviter tout problème de stabilité. Son principe est donné sur la figure 2.11, qui est un cas
particulier de la figure générale 2.7.
Le signal d’erreur est filtré par un filtre passe-bande sélectif centré sur la fréquence N νr .
C’est donc cette composante spectrale du spectre vibratoire de la machine qui doit être
atténuée et qui constitue la fonction de coût. Ce signal réactualise dans ce cas le contrôleur
adaptatif en utilisant l’algorithme décrit au paragraphe 2.2.4. Le signal de référence, qui
doit être de même fréquence que l’erreur, consiste en une sinusoı̈de dont la fréquence N νr
peut est déduite du signal de top-tour grâce à un multiplieur fréquentiel. Le même principe
de prémodulation et de filtrage permet alors à ce contrôleur d’éliminer du spectre vibra-
toire de la machine la composante visée, tout en générant une raie parasite incontrôlable
de fréquence (N + 2p)νr . On obtient donc les mêmes performances de contrôle qu’avec la
39
CHAPITRE 2. Méthodes classiques de contrôle actif
vibrations
statoriques
cos(2πpνr t) F cos(2πpνr t + φ) naturelles
vibrations
statoriques
contrôleur (N ± p)νr (N + p)νr N νr (N + 2p)νr système N νr (N + 2p)νr totales
N νr adaptatif (N − p)νr LIT
N νr
N νr
40
2.4. Conclusion
spectre
sans
compensation
(échelle linéaire)
∆ν = 100 Hz
spectre
avec composante spectrale
compensation parasite
(échelle linéaire)
fréquence (Hertz)
Fig. 2.12: Résultat de compensation active utilisant le contrôleur feedforward (issu de [Jar95])
2.4 Conclusion
Ce chapitre, composé de deux sections, a permis de montrer l’insuffisance des méthodes
classiques de contrôle actif, pour l’atténuation des vibrations d’une machine tournante en
utilisant le principe décrit au chapitre 1.
La première section est une introduction aux méthodes de contrôle actif, utilisées pour les
systèmes linéaires et invariants dans le temps. L’étude du contrôleur de structure feedforward
y est particulièrement mise en avant, puisque c’est ce choix qui sera fait pour l’application
étudiée. La définition d’un coût quadratique fonction de ce contrôleur, permet de déterminer
l’expression de son transfert optimal, ainsi qu’un algorithme récursif convergent vers cette
solution.
Une étude des travaux de recherche déjà effectués sur cette application de contrôle actif
nous permet de voir, dans la seconde section, que l’utilisation de ces méthodes ne donne
41
CHAPITRE 2. Méthodes classiques de contrôle actif
42
Chapitre 3
3.1 Introduction
Ce chapitre livre les outils théoriques nécessaires à l’étude du type de système rencontré
dans le problème de compensation active posé au chapitre 1.
La section 3.2 est consacrée à l’étude théorique des systèmes linéaires variant périodiquement
dans le temps. Les propriétés de leur réponse impulsionnelle, ainsi que leurs propriétés
fréquentielles, permettent de caractériser les effets qu’ils produisent sur le signal d’entrée.
Finalement, on aboutit à la définition d’une matrice, jouant le même rôle que la fonction
de transfert des systèmes linéaires invariants dans le temps. L’étude de cette matrice peut
fournir des informations importantes sur le système.
Enfin, dans la section 3.3, les outils théoriques précédemment définis sont utilisés pour ana-
lyser un système linéaire variant périodiquement dans le temps de même structure que celui
exposé à la section 1.4.5. Les résultats obtenus dans cette section seront donc utilisés au cha-
pitre suivant, et permettront d’élaborer un algorithme de contrôle actif de vibrations efficace
pour ce type de système.
43
CHAPITRE 3. Étude théorique du système
plus simple, et des résultats facilement utilisables peuvent être mis à jour. De plus, ce type
de systèmes se rencontre couramment dans la pratique. En effet, un échantillonneur est typi-
quement un système LVPT [Cla82], tout comme une modulation d’amplitude sinusoı̈dale, ou
un filtre à capacités commutées [Lux90]. Un système de sur- ou de sous-échantillonnage peut
lui aussi être considéré comme un système LVPT [Cro81]. Enfin, certains systèmes physiques
sont naturellement LVPT comme celui présenté au chapitre 1. Ce type de systèmes n’est donc
pas dénué d’intérêt, et leur étude détaillée va être réalisée dans cette section.
Après avoir argumenté le choix du modèle utilisé pour représenter ces systèmes, leur étude
est d’abord réalisée en temps, puis en fréquence. On aboutit alors à un opérateur de type
«gain complexe» que l’on peut mettre sous une forme matricielle, et qui caractérise leur effet
d’un point de vue fréquentiel. L’étude de ce «gain complexe matriciel» grâce à des techniques
d’algèbre matricielle (décomposition en valeurs singulières et calcul de normes) fournit des
informations intéressantes sur le système. Ces techniques permettent, en outre, de donner les
conditions d’inversibilité d’un tel opérateur.
Deux grands types de modèles sont utilisés dans la littérature, pour décrire les systèmes
LVPT.
Le premier est le modèle d’état, régi par des matrices périodiques en temps [Bit86, Kit88].
Son étude peut être réalisée grâce au théorème de Floquet [Wer90]. Il montre l’existence d’une
transformation linéaire variant périodiquement dans le temps, permettant de transformer le
système de départ en un système invariant. Ce dernier contient beaucoup d’informations sur
le système périodique original, qu’il soit continu ou discret [Cor79]. Pour les systèmes discrets,
on peut aussi utiliser une transformation par blocs permettant d’obtenir des matrices d’état
invariantes [Mey75].
Le second modèle est celui de la réponse impulsionnelle. L’étude de ce signal, dans le domaine
temporel aussi bien que dans le domaine fréquentiel, permet de caractériser complètement
les effets des systèmes LVPT sur leur entrée, qu’ils soient continus ou discrets.
Pour le premier modèle, il est nécessaire de choisir les variables d’état composant le vecteur
d’état du système. Ce choix n’est jamais facile, et peut même poser des problèmes si ces
variables ne résument pas parfaitement l’état du système à un instant donné. Par contre, la
réponse impulsionnelle ne peut être choisie, puisqu’elle est fournie par le système lui-même
en réponse à une impulsion. Cette deuxième solution écarte donc tout risque d’erreur de
modélisation, lié au choix des variables d’état. De plus, ce modèle permet une interprétation
fréquentielle quasi immédiate des effets du système étudié, puisque la transformée de Fourier
de la réponse impulsionnelle fournit directement la fonction de transfert en fréquence. Pour
ces différentes raisons, le modèle de la réponse impulsionnelle sera utilisé dans la suite de ce
mémoire pour représenter et étudier les systèmes LVPT.
La dernière question à se poser sur le modèle utilisé est de savoir s’il sera discret ou continu.
Le but de ce travail étant d’élaborer une stratégie de contrôle implantée, à terme, dans un
44
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
processeur numérique de signal, le système présenté au chapitre 1 à travers lequel les pertur-
bations doivent être compensées, sera vu par le processeur comme un système discret. Seuls
des systèmes LVPT discrets seront donc considérés ici. De plus, si leur réponse impulsion-
nelle provient de l’échantillonnage de celle d’un système continu, elle sera supposée avoir été
correctement échantillonnée, et vérifier la condition de Shannon énoncée au paragraphe 1.3.3,
afin d’éviter tout problème de recouvrement spectral.
3.2.2 Définition
Soit le système H discret, linéaire, d’entrée x(n) et de sortie y(n) à temps discret n,
représenté à la figure 3.1. où l’opérateur «◦» peut se lire «appliqué à» et a été défini au
x(n) H y(n)
paragraphe 2.2.3.
Les signaux x(n) et y(n) sont discrets, et la relation d’entrée-sortie qu’ils vérifient s’exprime
par :
y(n) = H ◦ x(n) (3.1)
Afin de définir correctement le concept de système discret variant ou invariant dans le temps,
il est nécessaire d’introduire l’opérateur retard de N échantillons. Noté RN et appliqué à un
signal s(n), il vérifie :
RN ◦ s(n) = s(n − N ) (3.2)
Le système linéaire H est dit invariant dans le temps (LIT) si, pour une même entrée
retardée de N échantillons, on obtient la même sortie, retardée de N échantillons. Sachant
que (3.1) est vérifiée, ceci revient à écrire :
Pour un tel système, les opérateurs H et RN peuvent donc commuter quelle que soit la valeur
du retard. En effet, en utilisant les relations (3.2) et (3.3), on obtient :
Si cette commutation n’est pas réalisée pour tout N , le système linéaire est dit variant
dans le temps (LVT).
Si la relation (3.4) est vérifiée pour une valeur particulière de N , elle l’est aussi pour tous
ses multiples entiers. Le système linéaire H est alors dit variant périodiquement dans le temps
(LVPT) et sa période, notée P , est la plus petite valeur de N pour laquelle la commutation
45
CHAPITRE 3. Étude théorique du système
(3.4) est possible. Ce type de systèmes est donc un cas particulier des systèmes LVT. De
plus, on remarque que les systèmes invariant dans le temps peuvent être considérés comme
périodiques de période P = 1, et constituent eux-mêmes un cas particulier des systèmes
LVPT.
réponse impulsionnelle
On se place dans le cas général où le système H est LVT. Il n’y a donc aucune hypothèse
sur la loi suivant laquelle l’équation (3.4) est vérifiée en fonction de N . Néanmoins, la pro-
priété de linéarité de H permet à ce système d’être entièrement caractérisé par sa réponse
impulsionnelle [Cla82]. Elle peut être exprimée en temps absolu, c’est à dire sans utiliser la
notion de retard ou de différence entre deux instants. Elle est alors notée h(n, k) et représente
la réponse du système observée à l’échantillon n, à une impulsion envoyée à l’échantillon k.
On peut définir deux autres formes de cette réponse, en introduisant par changement de
variables le retard r = n − k qui est le laps de temps entre l’instant d’observation de la sortie
et l’impulsion d’excitation [Cla82, Pra92] :
g(r, k) est donc la réponse du système, observée r échantillons après l’impulsion d’excitation
envoyée en k. Quant à f (n, r), elle peut s’interpréter comme étant la réponse du système ob-
servée à l’instant n, à l’impulsion d’excitation envoyée r échantillons auparavant. L’équation
reliant ces trois formes de la réponse impulsionnelle peut s’écrire, en éliminant r de (3.5) et
(3.6) :
h(n, k) = g(n − k, k) = f (n, n − k) (3.7)
Enfin, si le système considéré est causal, sa réponse impulsionnelle est nulle pour r = n−k < 0.
La figure 3.2 donne un exemple d’allure temporelle de h(n, k) et de g(r, k) pour un système
LVT causal. f (n, r) n’a pas été représentée d’une part parce qu’elle contient exactement la
même information, et d’autre part parce qu’elle est moins explicite d’un point de vue visuel.
On voit que la forme du signal obtenu en sortie du système varie en fonction de l’instant
d’excitation k. Si le système avait été invariant, toutes les réponses auraient été de forme et
d’amplitude identiques, et la réponse impulsionnelle n’aurait été fonction que du retard r.
relation d’entrée-sortie
La relation temporelle d’entrée-sortie d’un système LVT est donnée par [Hua80] :
+∞
X
y(n) = H ◦ x(n) = h(n, k)x(k) (3.8)
k=−∞
46
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
k k
instant instant
d’excitation n=k d’excitation
n
instant r
d’observation retard
L’équation (3.7) permet d’exprimer la sortie du système en fonction des deux autres formes
de la réponse impulsionnelle :
+∞
X
y(n) = g(n − k, k)x(k) (3.9)
k=−∞
+∞
X
y(n) = f (n, n − k)x(k) (3.10)
k=−∞
Ces équations montrent que la réponse impulsionnelle permet de relier complètement l’entrée
x(n) à la sortie y(n) du système, et qu’elle résume donc bien son comportement entrée-sortie.
La définition de f et g ne simplifie en rien la notion de réponse impulsionnelle, mais elle
permet de revenir rapidement au cas d’un système invariant. En effet, ces signaux discrets
ne sont alors plus fonction du moment d’excitation k ou d’observation n, mais seulement du
retard r = n − k. Les relations (3.9) et (3.10) se simplifient donc directement en une équation
de convolution, caractéristique des systèmes LIT.
réponse impulsionnelle
47
CHAPITRE 3. Étude théorique du système
r
retard
Fig. 3.3: Exemple de réponse impulsionnelle g(r, k) d’un système LVPT causal avec P = 2
relation d’entrée-sortie
Les relations de périodicité (3.13) et (3.14), ainsi que les relations d’entrée-sortie (3.9)
et (3.10) vont nous permettre de déterminer la structure interne d’un système LVPT de
période P , et nous aider ainsi à comprendre son effet sur le signal d’entrée. Par exemple, le
signal bidimensionnel g(r, k) est P -périodique suivant k, et peut être développé en série de
Fourier discrète (ou d’une manière équivalente, en transformée de Fourier discrète) suivant
cette variable [Gas90] :
P
X −1
l
g(r, k) = gl (r)ej2π P k (3.15)
l=0
avec
P −1
1 X l
gl (r) = g(r, k)e−j2π P k (3.16)
P
k=0
La relation d’entrée-sortie (3.9) peut donc s’écrire, en utilisant (3.15) :
+∞
ÃP −1 !
X X l
y(n) = gl (n − k)ej2π P k x(k)
k=−∞ l=0
P −1
à +∞
!
X X
j2π Pl k
= gl (n − k)x(k)e (3.17)
l=0 k=−∞
48
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
La sortie y(n) du système considéré est donc constituée d’une somme de P éléments. Chacun
de ces éléments est la sortie d’un système LIT de réponse impulsionnelle gl (r) dont l’entrée
est une version modulée de x(n) par une exponentielle complexe de fréquence Pl . En effet, on
l
reconnaı̂t en (3.17) la somme de P équations de convolutions entre gl (n) et x(n)ej2π P n . La
figure 3.4 représente la traduction graphique de l’équation (3.17), et montre qu’un système
LVPT P -périodique peut être considéré comme un banc de P filtres parallèles à entrées
modulées. Cette représentation est dénommée «structure à modulation d’entrée».
e0 = 1
g0 (r)
1
ej2π P n
P −1
ej2π P n
gP −1 (r)
Une structure tout à fait équivalente peut être déterminée en décomposant f (n, r) P -
périodique suivant n en série de Fourier discrète :
P
X −1
l
f (n, r) = fl (r)ej2π P n (3.18)
l=0
avec
P −1
1 X l
fl (r) = f (n, r)e−j2π P n (3.19)
P
n=0
+∞
ÃP −1 !
X X
j2π Pl n
y(n) = fl (n − k)e x(k)
k=−∞ l=0
P −1
à +∞
!
X X
j2π Pl n
= e fl (n − k)x(k) (3.20)
l=0 k=−∞
Le signal d’entrée passe donc d’abord dans P systèmes LIT de réponse impulsionnelle fl (r).
Les sorties de ces filtres sont ensuite modulées en amplitude par des exponentielles complexes
de fréquence Pl . La structure correspondante est représentée sur la figure 3.5. Par analogie
avec le cas précédent, cette représentation sera appelée «structure à modulation de sortie».
49
CHAPITRE 3. Étude théorique du système
e0 = 1
f0 (r)
1
ej2π P n
P −1
ej2π P
n
fP −1 (r)
Afin d’établir un lien entre les systèmes LIT présents dans ces deux structures, l’équation
(3.17) peut s’écrire :
P −1
à +∞ !
l l X X l
y(n) = ej2π P n e−j2π P n gl (n − k)x(k)ej2π P k
l=0 k=−∞
P −1
à +∞
!
X X
j2π Pl n −j2π Pl (n−k)
= e gl (n − k)e x(k) (3.21)
l=0 k=−∞
Par identification entre (3.20) et (3.21), la relation entre fl (r) et gl (r) est donnée par :
l
fl (r) = gl (r)e−j2π P r (3.22)
Les réponses impulsionnelles fl (r) et gl (r) sont donc liées par une modulation d’amplitude
sinusoı̈dale complexe de fréquence Pl , ce qui équivaut à un décalage fréquentiel de même
fréquence. Cette relation permet de confirmer l’équivalence entre les représentations à modu-
lation d’entrée et de sortie.
La décomposition en série de Fourier de la réponse impulsionnelle d’un système LVPT P -
périodique conduit donc aux structures décrites par les figures 3.4 et 3.5, qui permettent
d’analyser facilement l’effet d’un tel système. Les différentes modulations provoquent des
décalages fréquentiels de Pl (l ∈ {0, . . . , P − 1}), alors que le banc de P systèmes LIT à
un effet de filtrage linéaire et invariant sur les signaux. On retrouve finalement en sortie P
versions fréquentiellement décalées et filtrées de l’entrée. L’effet caractéristique d’un système
LVPT P -périodique sera donc à la fois de filtrer, mais aussi de décaler le signal d’entrée en
fréquence, ce que ne fait pas un système LIT.
Cette caractéristique permet d’ores et déjà d’affirmer que les exponentielles complexes ne
sont pas les signaux propres des systèmes LVPT. En effet, si l’entrée est une exponentielle
complexe, on ne retrouve pas en sortie la même exponentielle amplifiée et déphasée comme
pour les systèmes LIT [Lac96, page 16], mais une somme de P exponentielles de fréquences
50
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
différentes. Cette particularité compromet donc l’existence d’un gain complexe du même type
que celui des systèmes LIT. Ce point sera discuté plus longuement à la section 3.2.4.
Enfin, ces deux structures peuvent être vues comme une généralisation du cas invariant. Si
le système analysé est LIT, il peut être considéré comme LVPT de période P = 1, et seule
la première branche des représentations existe. Le système est alors caractérisé par la seule
réponse impulsionnelle g0 (r) = f0 (r), et sa relation d’entrée-sortie est donnée par :
+∞
X
y(n) = g0 (n − k)x(k) (3.23)
k=−∞
Après avoir réalisé une étude temporelle des systèmes LVPT en analysant leurs effets
grâce aux propriétés de leurs réponses impulsionnelles, nous allons nous attacher à les décrire
d’une manière purement fréquentielle.
réponse fréquentielle
La réponse fréquentielle d’un système LIT, appelée aussi gain complexe, est définie comme
la transformée de Fourier de sa réponse impulsionnelle [Lac96, pages 16-17]. Pour de tels
systèmes, une transformée de Fourier monodimensionnelle suffit, la réponse impulsionnelle
étant fonction de la seule variable r. Le gain complexe représente alors les effets du système
sur chacune des composantes fréquentielles du signal d’entrée, en terme d’amplification et de
déphasage.
Dans le cas d’un système LVT, les diverses formes de la réponse impulsionnelle sont
toutes fonctions de deux variables temporelles. Il est donc nécessaire, dans ce cas, d’utiliser la
transformation de Fourier bidimensionnelle. Une telle opération, définie par (1.5) et appliquée
à l’équation (3.7), permet d’écrire les relations :
ˆ ˆ β + α) = fˆˆ(α + β, −β)
ĥ(α, β) = ĝ(α, (3.24)
Les fonctions définies en (3.24) sont complexes, et dépendent des deux variables fréquentielles
α et β, associées respectivement aux variables temporelles n et k. Elles sont appelées «réponses
en fréquence bidimensionnelles» ou «réponses bifréquentielles». Dans la suite, on se limitera
ˆ
à l’étude de ĥ(α, β), puisque ces trois fonctions sont toutes équivalentes et contiennent la
même information, à un changement de variable près.
Pour représenter graphiquement les caractéristiques d’une telle fonction, comme son mo-
dule ou sa phase, il est nécessaire de définir un repère bidimensionnel, nommé plan bifréquence
[Zad50], dont les abscisses et les ordonnées sont des fréquences réduites. La localisation de
ˆ
ĥ(α, β) dans ce plan pour les systèmes linéaires, causaux et stationnaires a été étudié dans
[Cla82]. Avant d’énoncer les conclusions de cet article, il est intéressant de développer le
concept de système stationnaire. Il découle du théorème de Wold [Pic93, pages 516-522] qui
51
CHAPITRE 3. Étude théorique du système
statue que toute réalisation d’un signal aléatoire stationnaire peut se décomposer en deux
parties :
– une partie purement stochastique, pouvant être générée par la sortie d’un système LIT
excité par un bruit blanc,
– une partie purement déterministe, consistant en une somme de sinusoı̈des, parfaitement
définies par leur fréquence, leur amplitude et leur phase.
Un signal déterministe est alors dit stationnaire s’il est exclusivement constitué d’une somme
de signaux sinusoı̈daux. De même, un système est dit stationnaire si, excité par une entrée
déterministe stationnaire, il délivre un signal de sortie déterministe stationnaire. Les auteurs
ˆ
montrent alors que pour un système linéaire, stationnaire et causal, la fonction ĥ(α, β) ne peut
être localisée que sur un réseau de droites parallèles, et de pente −1 dans le plan d’abscisses
α et d’ordonnées β.
Il faut remarquer l’indépendance des propriétés de stationnarité et d’invariance dans le temps.
Un exemple frappant est celui d’un système LVPT P -périodique. D’après la section 3.2.3, ses
effets se limitent à des décalages fréquentiels proportionnels à P1 , et à des filtrages LIT. À
une entrée déterministe stationnaire, un tel système fait donc correspondre une sortie elle
aussi déterministe stationnaire. En effet, même si les sinusoı̈des de l’entrée sont décalées en
fréquence et filtrées par des systèmes LIT, le signal de sortie reste constitué d’une somme de
sinusoı̈des. Un tel système est donc stationnaire, tout en étant variant dans le temps. S’il est
de plus causal, sa réponse bifréquentielle restera limitée à des droites de pente −1 dans le
plan bifréquence.
relation entrée-sortie
52
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
ˆ
La réponse bifréquentielle ĥ(α, β) est donc bien une extension bidimensionnelle du gain com-
plexe d’un système LIT. En effet, comme le montre cette équation, elle permet de relier le
spectre du signal de sortie à celui du signal d’entrée. La figure 3.6 illustre la relation (3.26) en
montrant graphiquement ce lien pour un système linéaire, stationnaire et causal. La réponse
bifréquentielle reste alors confinée sur des droites parallèles de pente −1, avec des ordonnées à
l’origine quelconques. Le spectre d’amplitude du signal d’entrée, représenté sur la gauche de la
figure, est celui d’un signal monofréquentiel réel. Ses composantes fréquentielles «traversent»
ˆ
le plan bifréquence horizontalement pour aller se «refléter» sur les droites constituant ĥ(α, β).
Ceci donne alors naissance à de nouvelles composantes dans le spectre d’amplitude de la sor-
tie, représenté en bas de cette figure. À des fins de simplification, seules les droites où est
ˆ
localisée ĥ(α, β) sont représentées, et non son module.
β β
+ 12 + 12
|x̂(β)| α
− 12 + 12
− 12 − 12
|ŷ(α)|
− 12 + 12
α
Fig. 3.6: Exemple de relation entre les spectres d’entrée et de sortie d’un système linéaire
stationnaire et causal
réponse fréquentielle
Si le système considéré est LVPT de période P , les équations (3.7), (3.15) et (3.18) per-
mettent d’écrire la réponse impulsionnelle sous deux formes :
P
X −1 P
X −1
j2π Pl k l
h(n, k) = gl (n − k)e = fl (n − k)ej2π P n (3.27)
l=0 l=0
53
CHAPITRE 3. Étude théorique du système
P
X −1 P
X −1
ˆ l l
ĥ(α, β) = gˆl (α)δ(α + β − ) = f̂l (−β)δ(α + β − ) (3.28)
P P
l=0 l=0
D’après (3.28), la réponse bifréquentielle d’un système LVPT reste confinée sur des droites
parallèles de pente −1, comme pour un système LVT. Ceci est dû à la présence des distribu-
ˆ
tions de Dirac dans l’expression de ĥ(α, β). Toutefois, la contrainte de périodicité impose à
ces droites d’être régulièrement espacées de P1 , fréquence de variation du système. Cette pro-
priété montre bien que les systèmes LVPT sont un cas particulier des systèmes LVT, et que
la périodicité de leurs variations se retrouve directement dans le plan bifréquence à travers
cette contrainte d’espacement.
Enfin, on retrouve sur chacune de ces droites les gains complexes gˆl (α) ou f̂l (−β), corres-
pondant aux P systèmes LIT constituant le système LVPT, et qui apparaissent sur les figures
3.4 et 3.5. On peut vérifier l’égalité entre gˆl (α) et f̂l (−β) en appliquant une transformée de
Fourier à (3.22), grâce à laquelle on obtient :
l
f̂l (λ) = gˆl (λ + ) (3.29)
P
l
f̂l (−β) = f̂l (α − ) = gˆl (α) (3.30)
P
ˆ
Les expressions de ĥ(α, β) données par (3.28) sont donc parfaitement équivalentes et contien-
nent la même information. Dans la suite, la réponse bifréquentielle sera préférentiellement
exprimée en fonction des gains complexes gˆl (α). ¯ ¯
¯ˆ ¯
La figure 3.7 représente en perspective ¯ĥ(α, β)¯ pour un système LVPT P -périodique
causal. Les seuls points non nuls sont situés sur les droites de pente −1 et d’ordonnée à l’origine
proportionnelle à P1 . Cette réponse bifréquentielle pourrait être obtenue, par exemple, par
échantillonnage puis par transformée de Fourier de la réponse impulsionnelle d’un système
LVPT continu.¯ Puisqu’on
¯ suppose vérifier la condition d’échantillonnage de Shannon, il est
¯ˆ ¯
normal que ¯ĥ(α, β)¯ soit nul lorsque |α| ou |β| sont égaux à 12 . L’emplacement des différents
ˆ
gains complexes gˆl (α) est déterminé grâce à la propriété de périodicité de la fonction ĥ(α, β),
qui est 1-périodique suivant α et β.
54
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
ˆ
ĥ(α, β)
β |g0 |
|gP −1 |
|gP −3 |
α
− 12
− 12
Fig. 3.7: Exemple de module de réponse bifréquentielle pour un système LVPT P -périodique
causal
relation entrée-sortie
Les équations (3.28) et (3.26) permettent d’établir la relation d’entrée-sortie d’un système
LVPT P -périodique dans le domaine fréquentiel :
P −1
" Z 1 #
X 2 l
ŷ(α) = gˆl (α) x̂∗ (β)δ(α + β − )dβ
−1 P
l=0 2
P
X −1
l
= gˆl (α)x̂∗ (−α + ) (3.31)
P
l=0
Ceci montre que gˆ0 (α) représente la partie invariante du système LVPT, alors que les gains
complexes gˆl (α) pour l ≥ 1 caractérisent sa partie variante dans le temps.
La figure 3.8 représente schématiquement, pour le même système que celui de la figure 3.7,
les liens existant entre le spectre du signal d’entrée
¯ et celui
¯ du signal de sortie. Les droites de
¯ˆ ¯
pente −1 en pointillés représentent les points de ¯ĥ(α, β)¯ pouvant être non nuls, alors que les
55
CHAPITRE 3. Étude théorique du système
mêmes droites en traits pleins représentent les points réellement non nuls de cette fonction.
Elle peut être lue de deux manières différentes :
– sens entrée → sortie : la composante spectrale de l’entrée située en β = β0 traverse la
réponse bifréquentielle
¯ ¯ horizontalement (lignes tirets). Elle peut influencer, suivant les
¯ˆ ¯
valeurs de ¯ĥ(α, β)¯, le spectre de la sortie aux fréquences −β0 + Pl , avec l ∈ {0, . . . , P −
1}.
– sens sortie → entrée : on remonte cette fois la réponse bifréquentielle verticalement
(lignes point-point-tiret). la composante spectrale
¯ de¯ la sortie à la fréquence α0 peut
¯ˆ ¯
donc être influencée, suivant les valeurs de ¯ĥ(α, β)¯, par le spectre de l’entrée aux
l
fréquences −α0 + P, avec l ∈ {0, . . . , P − 1}. C’est l’illustration directe de l’équation
(3.31).
β β
gˆ0 gˆ1 gˆ2 ĝP −1
+ 12 + 12
1
ĝP −1
−α0 + P 1
ĝP −2 P
−α0
1
−α0 − P
|x̂(β)| α
− 12 + 12
β0
gˆ1
− 12 − 12
|ŷ(α)|
− 12 α0 −β0 + 12
α
1
−β0 − P
Fig. 3.8: Exemple de relation entre les spectres d’entrée et de sortie pour un système LVPT
P -périodique causal
P −1
Si le système analysé est LIT, alors la somme Pl=0 ||ĝl (α)||2 se réduit à ||gˆ0 (α)||2 et ∆ prend
la valeur 0. Au contraire, si le système ne comporte pas de partie invariante, alors ||gˆ0 (α)||2
devient nul et ∆ vaut 1.
56
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
Les résultats obtenus dans le domaine fréquentiel confirment un résultat déjà ennoncé à la
section 3.2.3, où il avait été montré que les exponentielles complexes n’étaient pas les signaux
propres des systèmes
¯ LVPT.
¯ En effet, une exponentielle complexe en entrée produira, suivant
¯ˆ ¯
les valeurs de ¯ĥ(α, β)¯, non pas une, mais un ensemble d’exponentielles complexes en sortie.
La notion de fonction de transfert est donc à revoir pour les systèmes LVPT, et la section
suivante est dédiée à définir un opérateur du même type utilisable pour ce genre de systèmes.
systèmes LIT
Il est bien connu que les exponentielles complexes sont les signaux propres des systèmes
LIT. Cette propriété peut être démontrée à partir de leur relation d’entrée-sortie en temps,
donnée par (3.23). Si le signal d’entrée x(k) est une exponentielle complexe d’amplitude
unitaire, de phase nulle et de fréquence réduite α, cette relation devient :
+∞
X
y(n) = g0 (n − k)ej2παk
k=−∞
L’équation (1.2) permet d’identifier le terme entre crochets comme la valeur de la transformée
de Fourier de la réponse impulsionnelle du système LIT, à la fréquence α. Cette équation peut
donc s’écrire :
y(n) = gˆ0 (α)ej2παn (3.35)
Finalement, on retrouve en sortie la même exponentielle complexe qu’en entrée pondérée par
le nombre complexe gˆ0 (α), d’où le terme de «signaux propres». Ce nombre, dépendant de la
fréquence, est appelé réponse fréquentielle ou gain complexe. Il résume, par son module et sa
phase, l’effet de filtrage qu’a ce type de système sur son entrée. On comprend alors qu’il soit
important de pouvoir décomposer un signal sur une base d’exponentielles complexes, ce qui
est réalisé par la transformation de Fourier.
La définition des signaux propres des systèmes LIT nous a conduit à la notion de gain
complexe, la même démarche va donc être adoptée pour les systèmes LVPT.
systèmes LVPT
Pour de tels systèmes, les signaux propres ne sont plus de simples exponentielles com-
plexes. En effet, supposons que le signal d’entrée x(k) d’un système LVPT P -périodique soit
57
CHAPITRE 3. Étude théorique du système
Le double effet de filtrage et de décalage fréquentiel de ces système, mis en évidence dans les
sections précédentes, est confirmé par cette équation. En effet, on ne retrouve pas en sortie
la seule exponentielle complexe d’entrée, mais la somme de P exponentielles complexes de
¡ ¢
fréquence α + Pl , pondérées respectivement par gˆl α + Pl pour l = 1, . . . , P .
Une exponentielle complexe en entrée engendre donc P exponentielles complexes en sortie,
décalées d’un multiple de P1 en fréquence. La TFr étant une fonction 1-périodique, on peut
espérer qu’un signal d’entrée constitué de P exponentielles régulièrement espacées de P1 en
fréquence, redonne le même type de signal en sortie. Afin de vérifier cette proposition, prenons
un signal d’entrée de la forme :
P
X −1
m
x(k) = ej2π(α+ P )n (3.36)
m=0
58
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
¡ ¢ q
Les fonctions gˆl α + Pq et ej2π(α+ P )n étant 1-périodique suivant α, l’indice q du deuxième
terme de Sl peut être exprimé modulo P , ce qui conduit à :
P
X −1 ³ q ´ j2π(α+ q )n X ³
l−1
q ´ j2π(α+ q )n
Sl = ĝl α + e P + ĝl α + e P
P P
q=l q=0
P
X −1 ³ q ´ j2π(α+ q )n
= ĝl α + e P
P
q=0
La réponse du système à une somme d’exponentielles complexes de la forme (3.36) est donc
une somme pondérée des mêmes exponentielles complexes, aux mêmes fréquences. Ce résultat
montre bien que ces signaux constituent les signaux propres des systèmes LVPT.
On peut donc espérer trouver une relation d’entrée-sortie simple en décomposant les
signaux sur une base de P exponentielles complexes, décalées en fréquence de P1 . Pour
un signal s(n), ceci revient à «éclater» sa transformée de Fourier ŝ(α) en P intervalles
fréquentiels différents de longueur P1 . On sélectionne alors, à un même instant, une compo-
sante fréquentielle dans chaque intervalle. Rangées sous forme vectorielles, ces composantes
forment le vecteur ŝ(α) suivant :
ŝ (α)
¡ ¢ ¸ ¸
ŝ α − P1 1 1
ŝ(α) =
.. avec α ∈ − 2P ; 2P (3.39)
.
¡ ¢
ŝ α − PP−1
Ceci constitue une description fréquentielle équivalente à la transformée de Fourier, mais
en observant P composantes fréquentielles à la fois. Cette décomposition de la réponse
fréquentielle est appelée «décomposition en composantes de recouvrement» dans [She94],
«réponse fréquentielle augmentée» dans [Zha97], ou encore «décomposition à modulations»
dans [Vet87] où la même étude a été réalisée dans le domaine de la transformée en z. Le
domaine de définition adopté pour α permet à ŝ(α) de décrire sans redondance toute la
transformée de Fourier du signal analysé. Un exemple graphique du lien entre la transformée
de Fourier s(α) et le vecteur ŝ(α) est donné par la figure 3.9 pour P = 3. Enfin, on peut
remarquer que si P = 1, ŝ(α) se réduit à ŝ(α).
59
CHAPITRE 3. Étude théorique du système
|ŝ(α)|
α
|ŝ(α)| − 16 1
6
|ŝ(α − 13 )|
|ŝ(α)| =
α
− 16 1
6
α |ŝ(α − 23 )|
− 12 − 16 0 1
6
1
2
α
− 16 1
6
Fig. 3.9: Représentation vectorielle de la transformée de Fourier d’un signal s(n) pour P = 3
Afin d’étendre la notion de gain complexe initialement définie pour les systèmes LIT aux
systèmes LVPT P -périodiques, prenons un tel système, d’entrée x(k) réelle et de sortie y(k).
Le but étant de donner la relation liant les vecteurs x̂(α) et ŷ(α), une composante quelconque
¡ ¢
du vecteur de sortie ŷ α − m P doit être exprimée en fonction de celles du vecteur d’entrée.
La relation d’entrée-sortie fréquentielle (3.32) permet d’écrire :
³ µ ¶
m´ X ³ m´
P −1
l+m
ŷ α − = ĝl α − x̂ α −
P P P
l=0
³ m´
P −1+m
X ³ m´ ³ q´
ŷ α − = ĝq−m α − x̂ α −
P q=m
P P
P
X −1 ³ m´ ³ q´
P −1+m
X ³ m´ ³ q´
= ĝq−m α − x̂ α − + ĝq−m α − x̂ α −
q=m
P P P P
q=P
¡ ¢
La TFr étant 1-périodique, l’indice q présent dans la fonction x̂ α − Pq peut s’exprimer
modulo P . Il en est de même pour la fonction ĝq−m (α − m P ) puisque l’équation (3.16) implique
l’égalité des fonctions gl (r) et gl+P (r). Le deuxième terme de cette expression peut donc être
¡ ¢
exprimé sous la forme d’une somme allant de 0 à m − 1, et la composante ŷ α − m P devient
finalement fonction de x̂(α) :
³ m´
P
X −1 ³ m´ ³ q´
ŷ α − = ĝq−m α − x̂ α −
P P P
q=0
h ³ m´ ³ m´ ³ m ´i
= ĝ−m α − ĝ1−m α − ··· ĝP −1−m α − · x̂(α)
P P P
60
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
En utilisant cette dernière expression, et le fait que la fonction gl (k) soit P -périodique suivant
l’indice l, on peut lier les vecteurs x̂(α) et ŷ(α) par la relation matricielle suivante :
¸ ¸
1 1
ŷ(α) = Ĝ(α) · x̂(α) avec α ∈ − ; (3.40)
2P 2P
où Ĝ(α) est une matrice carrée de dimension P × P , fonction de la fréquence et de la forme :
ĝ0 (α) ĝ1 (α) ··· ĝP −1 (α)
1
ĝP −1 (α − P) ĝ0 (α − P1 ) ··· ĝP −2 (α − P1 )
Ĝ(α) =
.. .. .. ..
(3.41)
. . . .
P −1 P −1 P −1
ĝ1 (α − P ) ĝ2 (α − P ) ··· ĝ0 (α − P )
Les équations (3.40) et (3.41) fournissent donc une relation simple entre les vecteurs x̂(α) et
ŷ(α), décrivant dans le domaine fréquentiel les signaux d’entrée et de sortie du système. La
matrice Ĝ(α) sera appelée «gain complexe matriciel» ou «réponse fréquentielle matricielle»
dans la suite de ce document. La section suivante est dédiée à son interprétation, et à l’étude
de quelques unes de ses propriétés.
D’après les équations (3.39), (3.40) et (3.41), l’élément situé sur la ligne m et la colonne
n de la matrice Ĝ(α) représente le gain complexe existant entre la composante fréquentielle
de fréquence α − Pn du signal d’entrée, et celle de fréquence α − m P du signal de sortie. Le
gain complexe matriciel peut donc prendre en compte les relations éventuelles existant entre
les signaux d’entrée et de sortie à des fréquences différentes. De telles relations n’existant pas
pour un système LIT, son gain complexe matriciel est diagonal, et seul ĝ0 (α) peut être non
nul.
De plus, un lien peut être établi entre la réponse bifréquentielle d’un système LVPT et
son gain complexe matriciel. h(α, β), définie à la section 3.2.4, représente le lien entre la
transformée de Fourier du signal d’entrée et celle du signal de sortie du système. Le gain
complexe matriciel joue le même rôle entre les réponses fréquentielles augmentées de ces
deux signaux. Or ce type de représentation, définie par (3.39), peut être vue comme une
version échantillonnée de la transformée de Fourier du signal. La matrice Ĝ(α) n’est donc
rien d’autre qu’une version échantillonnée de la réponse bifréquentielle, l’emplacement des
échantillons variant avec la fréquence α.
Enfin, on peut montrer que le gain complexe matriciel vérifie les mêmes propriétés de
base que le gain complexe classique. Pour cela, supposons deux systèmes LVPT H1 et H2 de
même fréquence. Leur association en parallèle est représentée à la figure 3.10(a). Le signal de
sortie global y est dans ce cas :
y = y1 + y2 = H1 ◦ x + H2 ◦ x
61
CHAPITRE 3. Étude théorique du système
H1 y1 x
H Hinv
y
x y x H1 H2 z x
y2 Hinv H x
H2
En représentant les signaux par leur réponse fréquentielle augmentée, et les systèmes par leur
gain complexe matriciel, on obtient :
³ ´
ŷ(α) = Ĝ1 (α) · x̂(α) + Ĝ2 (α) · x̂(α) = Ĝ1 (α) + Ĝ2 (α) · x̂(α)
Le gain complexe matriciel global de plusieurs systèmes LVPT de même fréquence en parallèle
est donc la somme des gains complexes matriciels de chaque système.
La figure 3.10(b) représente l’association en série des deux systèmes précédents. Le signal de
sortie global est cette fois donné par :
z = H2 ◦ H1 ◦ x
Le gain complexe matriciel global de plusieurs systèmes LVPT de même fréquence mis en série
est donc le produit des gains complexes matriciels de chaque système. Toutefois, contrairement
aux systèmes LIT, l’ordre des gains complexes matriciels ne peut être changé puisque le produit
matriciel n’est pas commutatif. En effet, dans le cas général, les systèmes LVT, et donc les
systèmes LVPT, ne peuvent être intervertis sans changer la forme de leurs signaux de sortie.
Le système inverse d’un système LVPT H, noté Hinv , est défini par la figure 3.10(c). Il vérifie
donc les équations :
H ◦ Hinv ◦ x = Hinv ◦ H ◦ x = x
Cette équation, traduite dans le domaine fréquentiel à l’aide des outils précédemment définis,
devient :
Ĝ(α) · Ĝinv (α) · x̂(α) = Ĝinv (α) · Ĝ(α) · x̂(α) = IP · x̂(α)
où IP est la matrice identité de dimension P × P . Alors, si le vecteur x̂(α) est différent du
vecteur nul, Ĝinv (α) vérifie les équations suivantes :
62
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
1 1
Si la matrice Ĝ(α) est régulière quelle que soit la fréquence α comprise entre − 2P et 2P , ces
équations sont résolues par :
¸ ¸
−1 1 1
Ĝinv (α) = Ĝ(α) avec α ∈ − ; (3.43)
2P 2P
−1
où Ĝ(α) est la matrice inverse de Ĝ(α). Le gain complexe matriciel de l’inverse d’un
système est donc la matrice inverse du gain complexe matriciel de ce système.
Le gain complexe matriciel associé aux systèmes LVPT vérifie donc, comme le gain com-
plexe classique, les propriétés énoncées ci-dessus. De plus, il prend en compte d’éventuels
liens entre des fréquences de valeurs différentes des signaux d’entrée et de sortie. Il peut donc
être interprété comme une extension de la notion de gain complexe, appliquée aux systèmes
LVPT. Des exemples de gain complexe matriciel de systèmes simples sont donnés dans la
suite.
système LIT
Si nous sommes en présence d’un système LIT (P = 1), les vecteurs et les matrices définis
précédemment redeviennent tous des scalaires. x̂(α) et ŷ(α) se réduisent en x̂(α) et ŷ(α), et
Ĝ(α) en gˆ0 (α). L’équation matricielle (3.40) redevient le simple produit de fonctions (3.33),
l’intervalle de définition de α étant alors donné par (1.4). Malgré tout, il peut être utile de
calculer le gain complexe matriciel de ce système avec une dimension supérieure à 1. En effet,
cette matrice sera nécessaire s’il est associé à des systèmes LVPT de période P > 1. Comme
il a été vu précédemment, les gains complexes ĝi (α) d’un système LIT sont nuls si i 6= 0. On
obtient donc, un gain complexe matriciel diagonal de la forme :
ĝ0 (α) 0
..
ĜLIT (α) =
.
(3.44)
P −1
0 ĝ0 (α − P )
Ce résultat confirme le fait que le gain complexe matriciel est une version échantillonnée de
la réponse bifréquentielle. En effet, dans le cas d’un système LIT, les seuls points de h(α, β)
non nuls sont ceux de sa diagonale principale.
Le deuxième exemple est une modulation d’amplitude sinusoı̈dale. Le système, formé d’un
simple multiplieur, est donné par la figure suivante :
Le signal modulateur est un cosinus d’amplitude A, de phase φ et de fréquence réduite P1 .
Ce système est donc LVPT P -périodique, et son signal de sortie est donné par :
µ ¶
2π
y(n) = A cos n + φ x(n)
P
63
CHAPITRE 3. Étude théorique du système
¡ 2π ¢
A cos P
n+φ
x(n) y(n)
x(n) = δn−k
ˆ Aejφ 1 Ae−jφ 1
ĥmod (α, β) = δ(α + β − ) + δ(α + β + )
2 P 2 P
Cette fonction étant 1-périodique en α et β, le deuxième terme peut s’écrire :
ˆ Aejφ 1 Ae−jφ P −1
ĥmod (α, β) = δ(α + β − ) + δ(α + β + )
2 P 2 P
Par identification avec (3.28), on peut voir que ce système LVPT est composé de deux
systèmes LIT de transformées de Fourier constantes :
A jφ
ĝ1 (α) = e
2
A −jφ
ĝP −1 (α) = e
2
Cette réponse bifréquentielle est donc confinée sur deux droites différentes, d’ordonnées à
l’origine ± P1 et de pente −1. Les deux morceaux de droite supplémentaires apparaissant dans
64
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
β
+ 12
1
+P
A ejφ
2
− 12 + 12
α
A e−jφ
2
1
−P
− 12
ˆ
Fig. 3.12: ĥmod (α, β) correspondant au système de la figure 3.11
0 ejφ e−jφ
e−jφ 0
A
Ĝmod (α) = 2
ejφ
0
ejφ e−jφ 0
les coins haut droit et bas gauche de la figure 3.12 appartiennent aux droites d’équation
ˆ
α + β ± PP−1 . Leur présence est due à la périodicité de ĥmod (α, β).
Finalement, la définition du gain complexe donnée en (3.41) montre que pour ce système,on
obtient une matrice P × P de la forme :
Sa structure est bien identique à celle de la réponse bifréquentielle de la figure 3.12.
Après avoir interprété le gain complexe matriciel comme une extension du gain complexe
classique, la prochaine section est dédiée à l’emploi de techniques d’algèbre matriciel afin
d’extraire des informations supplémentaires sur le système.
L’intérêt d’obtenir un gain complexe sous une forme matricielle est de pouvoir lui appli-
quer des techniques d’algèbre linéaire, susceptibles d’apporter des informations supplémen-
taires sur le système correspondant. Dans cette section, la décomposition en valeurs sin-
gulières, ainsi que le calcul de normes matricielles ont été employés dans ce but. Dans la
suite, l’opérateur (·)H représente le transposé complexe conjugué d’un vecteur ou d’une ma-
trice.
65
CHAPITRE 3. Étude théorique du système
définition
où :
– Û(α) est la matrice P × P unitaire formée des vecteurs propres de Ĝ(α) · Ĝ(α)H ,
– V̂(α) est la matrice P × P unitaire formée des vecteurs propres de Ĝ(α)H · Ĝ(α),
– Σ̂(α) est la matrice P × P diagonale dont la diagonale principale est formée de la racine
carrée des valeurs propres de Ĝ(α) · Ĝ(α)H , appelées «valeurs singulières». Elles sont
notées σi (α) et rangées par ordre décroissant :
σ̂1 (α) ≥ σ̂2 (α) ≥ · · · σ̂r (α) > σ̂r+1 (α) = · · · σ̂P (α) = 0
Seules les r premières valeurs singulières sont donc supposées non nulles, r pouvant
aller de 1 à P .
La figure 3.13 représente l’équation (3.46) sous une forme graphique. De plus, elle montre
σ̂1 (α)
r r
que l’on peut obtenir une expression réduite de l’équation (3.46), en négligeant les matrices
Û0 (α) et V̂0 (α) :
H
Ĝ(α) = Ûr (α) · Σ̂r (α) · V̂r (α) (3.47)
La matrice Σ̂r (α) est alors une matrice diagonale de dimension r × r, ne comportant aucune
valeur singulière nulle sur sa diagonale principale. Les matrices Ûr (α) et V̂r (α) sont, quant
à elles, de dimension P × r.
Interprétation
66
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
Ces deux vecteurs peuvent être considérés comme les représentations fréquentielles aug-
mentées des signaux temporels xp (n) et yp (n). L’équation (3.49) montre que ces nouvelles
coordonnées vérifient la relation :
La différence entre les relations (3.40) et (3.52) vient du fait que la matrice Σ̂(α) est toujours
diagonale. Une composante de x̂p (α) n’agit alors plus que sur une seule composante de ŷp (α)
par l’intermédiaire de la valeur singulière correspondante. Toutes les composantes de ces
vecteurs sont donc bien découplées. Le même type de démarche a été effectuée dans [Cab99]
en appliquant la décomposition en valeurs singulières sur la matrice de gains complexes d’un
système LIT multi-entrées multi-sorties, pour une fréquence fixe.
Les matrices Û(α) et V̂(α) peuvent être vues comme des matrices de changement de base,
H H
il en est donc de même pour leurs inverses Û(α) et V̂(α) . Les transformations décrites
par les équations (3.51) et (3.50) peuvent donc être vues comme deux changements de base
différents, permettant de passer des composantes de départ x̂(α) et ŷ(α) aux composantes
principales x̂p (α) et ŷp (α). Dans le cas présent, ces matrices peuvent aussi être vues comme
les gains complexes matriciels correspondant à différents systèmes :
H
– le gain complexe matriciel V̂(α) correspond à un système LVPT noté P1 ,
– le gain complexe matriciel Û(α) correspond à un système LVPT noté P2 ,
– le gain complexe matriciel Σ̂(α) étant diagonal, il correspond à un système LIT noté
I.
L’équation (3.48) permet alors de décomposer le système LVPT d’origine H en trois systèmes
associés en série, dont la sortie est reliée à l’entrée par la relation :
y = H ◦ x = P2 ◦ I ◦ P1 ◦ x
67
CHAPITRE 3. Étude théorique du système
Comme le montre la figure 3.14, H peut donc être vu comme deux systèmes LVPT différents
entourant un système LIT. Le premier système LVPT P1 permet de passer de l’entrée x(n)
H
xp (n) yp (n)
x(n) P1 I P2 y(n)
système
LIT
systèmes
LVPT
au signal xp (n). Ce signal passe ensuite par un système LIT I qui applique, à chaque compo-
sante fréquentielle principale, un gain égal à la valeur singulière correspondante. Le deuxième
système LVPT P2 permet alors de repasser de yp (n) au signal de sortie y(n) du système
global.
Ce modèle, mis à jour grâce à la décomposition en valeurs singulières du gain complexe
matriciel, permet donc de mieux comprendre la structure interne des systèmes LVPT. Nous
allons voir qu’il permet aussi de donner leurs conditions d’inversion.
Les équations (3.43) et (3.46) montrent que le gain complexe du système inverse de H
peut s’exprimer par :
¸ ¸
−1 −1 H 1 1
Ĝinv (α) = Ĝ(α) = V̂(α) · Σ̂(α) · Û(α) avec α ∈ − ;
2P 2P
Ce système peut donc être vu comme la mise en série des trois systèmes suivants :
H
– le système de gain complexe matriciel Û(α) , inverse de Û(α), et donc noté P2 inv ,
−1
– le système LIT de gain complexe matriciel diagonal Σ̂(α) , inverse de Σ̂(α), et noté
Iinv ,
H
– le système de gain complexe matriciel V̂(α), inverse de V̂(α) , et noté P1inv .
La figure 3.15 montre la structure interne de ce système. On peut noter qu’il est constitué
de deux systèmes LVPT et d’un système LIT, qui sont les inverses de ceux décrits par la
figure 3.14. L’existence des gains complexes matriciels de P1inv et P2inv est toujours assurée
−1
puisque ce sont les transposés hermitiens de ceux de P1 et P2 . La matrice Σ̂(α) , quant
à elle, n’existe pas forcément. En effet, sa diagonale principale est constituée des éléments
inverses de celle de Σ̂(α), et elle ne peut être calculée si l’un d’eux est nul. Un système sera
donc non-inversible en α si une ou plusieurs des valeurs singulières de son gain complexe
matriciel s’annulent pour cette fréquence.
Pour palier à ce problème et essayer d’inverser «au mieux» un système non-inversible, on
peut utiliser la notion de matrice inverse généralisée [Sch90, page 49]. Pour le gain complexe
68
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
Hinv
yp (n) xp (n)
y(n) P2 inv Iinv P1 inv x(n)
−g
matriciel Ĝ(α), elle est notée Ĝ(α) , et vérifie les propriétés suivantes :
−g
Ĝ(α) · Ĝ(α) · Ĝ(α) = Ĝ(α) (3.53)
−g −g −g
Ĝ(α) · Ĝ(α) · Ĝ(α) = Ĝ(α) (3.54)
−g −g
Ĝ(α) · Ĝ(α) et Ĝ(α) · Ĝ(α) sont hermitiennes (3.55)
−g −1 H
Ĝ(α) = V̂r (α) · Σ̂r (α) · Ûr (α) (3.56)
où V̂r (α), Σ̂r (α) et Ûr (α) sont définies par l’équation (3.47) et la figure 3.13. L’avantage
d’utiliser l’inverse généralisée par rapport à l’inverse exacte réside dans le fait qu’elle existe
toujours, même si Ĝ(α) n’est pas inversible. En effet, la matrice diagonale Σ̂r (α) ne comporte
−1
aucun élément nul sur sa diagonale principale, de telle sorte que Σ̂r (α) est toujours définie.
−g −1
De plus, on peut remarquer que si Ĝ(α) est inversible, alors Ĝ(α) est égale à Ĝ(α) . Si, par
contre, Ĝ(α) est singulière, alors son inverse généralisée tend à avoir le même comportement
que l’inverse exacte (relations (3.53), (3.54), et (3.55)), sans pour autant vérifier la relation
stricte (3.42). En fait, l’équation (3.56) indique que si Ĝ(α) a une valeur singulière nulle, la
valeur singulière correspondante de son inverse généralisée restera nulle, alors qu’elle devrait
prendre une valeur infinie pour l’inverse exacte. Cette remarque peut être illustrée grâce à
un exemple simple. Soit la matrice G complexe :
1 0 0 0
j 0 0 0
G=
0 0 1 0
0 0 0 1
69
CHAPITRE 3. Étude théorique du système
Cette matrice est singulière, et ne comporte pas d’inverse exacte. Le calcul de son inverse
généralisée nécessite sa décomposition en valeurs singulières :
√
√1
2
0 0 √j2 2 0 0 0 1 0 0 0
j √
√ 0 0 2
G = 2 · 0 1 0 0 · 0 0 1 0
0 0 1 0 0 0 0 1
0 1 0 0
0 0 1 0 0 0 0 0 0 1 0 0
| {z } | {z } | {z }
U Σ VH
√1 0 0
√
2
2 0 0 1 0 0 0
√j
0 0
= 2·
0 1 0 · 0 0 1 0
0 1 0
0 0 1 0 0 0 1
0 0 1 | {z } | {z }
| {z } Σr Vr H
Ur
Cette décomposition, ainsi que l’équation (3.56) donnent l’inverse généralisée recherchée :
1 √ 1 j
1 0 0 √ − 0 0
√ 0 0 2 −j 2 0 0 2 2
0 0 0 2
0 0 0 0
G−g =
0 1 0 · 0 1 0 · 0 0 1 0 =
0 0 1 0
0 0 1 0 0 0 1
0 0 1 0 0 0 1
Les matrices G et G−g vérifient bien les relations (3.53), (3.54) et (3.55), et on a de plus les
deux relations :
1 0 0 0
0 0 0 0
−g
G ·G =
0 0 1 0
0 0 0 1
1 j
2 −2 0 0
j 1
2 0 0
G · G−g = 2
0 0 1 0
0 0 0 1
Ces produits donnent donc «presque» la matrice identité, mais l’inverse généralisée ne fait
qu’approcher le comportement de l’inverse exacte vérifiant (3.42).
Cette méthode d’inversion, appliquée à un système LVPT H de gain complexe Ĝ(α), donne
les résultats suivants :
– Si on se trouve à une fréquence α pour laquelle Ĝ(α) n’a aucune valeur singulière nulle,
elle est alors inversible. Le gain complexe matriciel du système inverse Hinv est donc,
pour cette fréquence, la matrice inverse exacte de Ĝ(α). Le gain complexe matriciel du
système constitué par la mise en série de H et Hinv est alors la matrice identité.
– Si on se trouve à une fréquence α pour laquelle Ĝ(α) a une ou plusieurs de ses valeurs
singulières qui s’annulent, elles seront conservées dans le gain complexe matriciel du
70
3.2. Étude des systèmes linéaires variant périodiquement dans le temps
système inverse. L’association en série des deux systèmes aura alors un gain complexe
matriciel global approchant «au mieux» la matrice identité. On peut montrer [Sch90,
pages 393-398] que le critère permettant de dire que cette matrice est approchée «au
mieux» est quadratique. En effet, l’inverse généralisée minimise la norme quadratique
−g
de la différence entre la matrice identité Ip et le produit Ĝ(α) · Ĝ(α) .
Ce résultat est très facilement interprétable pour un système LIT HLIT . Pour les fréquences
où son gain complexe ĝ0 (α) est différent de 0, le système inverse HLIT inv a alors comme gain
complexe ĝ01(α) . Par contre, pour les fréquences où ĝ0 (α) est nul, le gain complexe de HLIT inv
l’est aussi. L’association en série de ces deux systèmes donne donc un gain complexe global
égal à 1, sauf pour les zéros de ĝ0 (α) où il est également nul.
La décomposition en valeurs singulières du gain complexe matriciel a donc permis de
préciser la structure interne générale d’un système LVPT, et d’établir ses conditions d’in-
versibilité en fonction de ses valeurs singulières. Nous allons maintenant explorer quelles
informations supplémentaires peut fournir le calcul de la norme de son gain complexe matri-
ciel.
Calcul de normes
Afin d’introduire correctement le concept de norme, il est nécessaire de mieux définir les
objets sur lesquels vont porter les calculs, c’est à dire les systèmes linéaires discrets, et les
signaux présents en entrée et en sortie de ces systèmes. Afin de vérifier la condition d’existence
de la transformée de Fourier (1.2) donnée au paragraphe 1.3.3, les signaux numériques qui
suivent sont supposés de module sommable. Pour un tel signal s, on peut montrer que l’on a
également :
+∞
X
|s(k)|2 < ∞
k=−∞
P
On parle alors de signal à «énergie finie», puisque la quantité +∞ 2
k=−∞ |s(k)| est homogène
à l’énergie contenue dans s. Si, de plus, on suppose que ces signaux sont à valeurs complexes,
alors ils peuvent être considérés comme les éléments d’un espace vectoriel sur le corps C des
complexes. En ajoutant à cet ensemble la notion de distance entre deux de ses éléments, on
obtient un espace vectoriel normé, que l’on peut montrer complet (il est alors appelé «espace
de Banach»). Enfin, en définissant le produit scalaire entre s1 et s2 par :
+∞
X Z + 12
∗
hs1 , s2 i = s1 (n) s2 (n) = ŝ1 (α)∗ ŝ2 (α)dα (3.57)
k=−∞ − 12
on obtient un ensemble dénommé espace de Hilbert des signaux à module sommable, noté H1 .
Dans la suite de ce paragraphe, on suppose que tous les signaux rencontrés sont des éléments
de H1 . Ils peuvent donc être vus comme des vecteurs, et représentés sous différentes formes
– la forme temporelle, constituée de l’ensemble de leurs échantillons,
– la forme fréquentielle classique, c’est à dire leur transformée de Fourier en fréquence
réduite,
71
CHAPITRE 3. Étude théorique du système
Les signaux étant équivalents à des vecteurs, leur «mesure» va pouvoir être quantifiée
grâce à une norme. Il existe une multitude de normes vectorielles, mais la plus utilisée en
traitement du signal reste la norme introduite par le produit scalaire (3.57). La norme au
carré d’un élément s de H1 est alors définie par :
+∞
X Z + 12
ksk22 = hs, si = 2
|s(k)| = |ŝ(α)|2 dα (3.58)
k=−∞ − 12
Cette valeur a une signification physique importante puisqu’elle représente l’énergie Es conte-
nue dans le signal. Elle peut encore s’exprimer en fonction de la représentation fréquentielle
augmentée de s par :
Z +1
2P
2
ksk2 = ŝ(α)H · ŝ(α)dα = kŝ(α)k22
1
− 2P
La norme de l’opérateur linéaire qu’est un système LVPT H est donnée par celle de son
gain complexe matriciel Ĝ(α). Une fois de plus, il existe une multitude de normes matri-
cielles, mais les normes d’opérateurs sont le plus souvent induites par des normes vectorielles.
La norme quadratique d’un système LVPT d’entrée x et de sortie y, induite par la norme
vectorielle (3.58), est définie par [Gol96, pages 14-16] :
Pour calculer cette norme, il faut donc chercher le signal d’entrée non nul, qui engendre le
«gain en norme» maximal entre l’entrée et la sortie du système. La valeur de cette norme
élevée au carré est équivalente au «gain en énergie» maximal que peut présenter le système.
En effet, il a été vu que la norme au carré d’un signal représente l’énergie qu’il contient. Il
est montré, dans [Mir96], que cette norme matricielle est donnée par :
h i
kHk2 = kĜ(α)k2 = sup σmax Ĝ(α)
− 12 <α≤ 12
h i
où σmax Ĝ(α) représente la plus grande valeur singulière de la matrice Ĝ(α).
La norme d’un système LVPT est donc la plus grande valeur singulière que puisse avoir le gain
complexe matriciel, sur l’ensemble de définition de α. Cette valeur représente l’amplification
72
3.3. Application
maximum que peut provoquer le système de son entrée vers sa sortie. Par exemple, pour un
système LIT, cette norme est égale au maximum du module de son gain complexe |ĝ0 (α)|.
Ce résultat permet donc de définir un gain d’entrée-sortie pour les systèmes LVPT, mais il
permet aussi de donner un sens plus physique aux valeurs singulières. En effet, leurs valeurs,
élevées au carré, peuvent être vues comme de véritables gains d’entrée-sortie en énergie,
appliqués sur les composantes fréquentielles principales définies au paragraphe précédent.
Les outils théoriques introduits dans cette section vont pouvoir être utilisés afin d’étudier
les caractéristiques du système LVPT présent dans le transfert décrit à la fin du chapitre 1.
3.3 Application
Cette section est consacrée à l’utilisation des outils théoriques précédemment définis, dans
le but d’étudier un système LVPT de même structure que celui décrit au chapitre 1, situé
entre les courants de commande et les vibrations qu’ils engendrent.
Tout d’abord, le paragraphe 3.3.1 rappelle la structure de ce système qui est à temps continu.
Sa réponse impulsionnelle ainsi que sa réponse bifréquentielle y sont aussi étudiées.
Ensuite, les contraintes que doit vérifier la fréquence d’échantillonnage, pour que l’échantil-
lonnage de ce système permette d’obtenir un système LVPT discret sont données au para-
graphe 3.3.2.
Enfin, les trois derniers paragraphes (3.3.3 à 3.3.5) sont consacrés à l’étude des caractéristiques
de ce système discret, en reprenant le plan utilisé dans la section précédente : domaine tem-
porel, domaine fréquentiel, puis détermination et analyse du gain complexe matriciel.
fc (t) vc (t)
ic (t) G v(t)
système LVPT
courant de commande, est modulée par un signal sinusoı̈dal m(t) = A cos (2πν0 t + φ) de
73
CHAPITRE 3. Étude théorique du système
La structure du système à étudier étant connue, on doit fixer les valeurs des différentes
variables qu’il comporte pour détailler l’allure de sa réponse impulsionnelle. Dans un cas
pratique, il faudrait donc passer par une phase d’identification afin de déterminer les ca-
ractéristiques de m(t) et la réponse impulsionnelle g(τ ). Dans notre cas, nous allons fixer ces
valeurs numériquement pour illustrer les caractéristiques d’un tel système :
– L’amplitude et la phase du signal de modulation sont choisies d’une manière arbitraire.
Sa fréquence, quant à elle, est imposée égale à la fréquence statorique nominale de la
machine de test, c’est à dire ν0 = 50 Hz, ce qui correspond à une fréquence de rotation
de νr = 25 Hz. Ce signal est alors donné par :
µ ¶
2π
m(t) = 2, 5 cos 2π × 50 × t + (3.59)
3
– Afin d’avoir une réponse fréquentielle réaliste, tout en étant assez simple pour faciliter
l’interprétation des résultats, on choisit un système G dont le module du gain complexe
est représenté sur la figure 3.17(a). Sa réponse en fréquence présente deux résonances à
0,5 et 1,5 kHz, ainsi qu’une anti-résonance à 1 kHz, puis a tendance à s’annuler pour les
fréquences supérieures à 4 kHz. La réponse impulsionnelle correspondante est donnée
à la figure 3.17(b).
Une grandeur intéressante à calculer afin de caractériser le système LVPT de la figure
3.16, est sa réponse impulsionnelle h(t, s). Il faut tout d’abord noter que le modèle employé
est un modèle à temps continu, ce qui impose aux variables temporelles d’être des variables
continues, contrairement à ce qui a été décrit dans la section précédente. On peut néanmoins
définir h(t, s) de la même manière que dans la partie 3.2.3 : c’est la réponse du système
observée au temps t, à une impulsion envoyée au temps s. Une impulsion à temps continu
étant modélisée par une distribution de Dirac, l’entrée du système est alors :
ic (t) = δ(t − s)
74
3.3. Application
5.5 0.2
5
0.15
4.5
0.1
4
3.5
échelle linéaire
échelle linéaire
0.05
3
0
2.5
2 −0.05
1.5
−0.1
1
−0.15
0.5
0 −0.2
0 1000 2000 3000 4000 5000 0 0.005 0.01 0.015 0.02 0.025 0.03
fréquence réelle (Hertz) temps (Secondes)
Cette réponse impulsionnelle est bien celle d’un système LVPT, puisqu’elle vérifie la relation
h(t, s) = h(t + T0 , s + T0 ) avec T0 = ν10 , relation identique à (3.12), mais pour les systèmes
continus.
À partir de cette réponse impulsionnelle, on peut définir, comme en 3.2.4, la réponse bifréquen
ˆ
-tielle ĥ(ξ, ζ). Elle correspond à la transformée de Fourier à deux dimensions de h(t, s), mais
ne possède pas de propriété de périodicité, comme son homologue définie en (3.28) pour les
systèmes discrets. En partant de l’équation (3.60), on obtient :
· ¸ · ¸
ˆ A −jφ A jφ
ĥ(ξ, ζ) = e ĝ(ξ) δ(ξ + ζ + ν0 ) + e ĝ(ξ) δ(ξ + ζ − ν0 ) (3.61)
2 2
Cette fonction ne peut être non nulle que sur les deux droites de pente −1 et d’ordonnées
à l’origine ±ν0 . Un exemple schématique est donné sur la figure 3.18. Le module du gain
complexe de G est représenté sur la partie gauche de la figure. Étant de largeur de bande B,
£ ¤
il est donc nul si la fréquence est à l’extérieur de l’intervalle − B2 ; + B2 . Le résultat obtenu
pour le module de la réponse bifréquentielle est montré sur la partie droite de cette figure.
ˆ
On peut voir que |ĥ(ξ, ζ)| s’annule si |ξ| > B2 ou |ζ| > B2 + ν0 . Le domaine où la réponse
bifréquentielle est non nulle est donc fonction de la largeur de bande de G et de la fréquence
75
CHAPITRE 3. Étude théorique du système
ζ ¯ ¯
¯ˆ ¯
¯ĝ(ξ, ζ)¯
B
2
+ ν0
A
2
|ĝ(ξ)|
|ĝ(ξ)|
ν0
+ B2
ξ
ξ − B2
− B2 + B2 −ν0
− B2 − ν0
3.3.2 Échantillonnage
76
3.3. Application
Pour que cette réponse impulsionnelle soit celle d’un système linéaire discret ℘-périodique,
elle doit vérifier la relation (3.12), c’est à dire :
νe = ℘ν0 (3.64)
Pour que le système discrétisé soit périodique, il faut que la fréquence d’échantillonnage soit
égale à un nombre entier de fois celle du système LVPT continu. Il est donc nécessaire de syn-
chroniser l’échantillonnage des signaux d’entrée et de sortie du système sur la fréquence de va-
riation de ce système. Cette méthode est appelée «échantillonnage synchrone». Dans le cadre
de l’application décrite au chapitre 1, le signal de top-tour fournit directement une sinusoı̈de
de fréquence νr . En utilisant ce signal, il est donc plus facile de synchroniser l’échantillonnage
des signaux sur νr , plutôt que sur ν0 . Au lieu d’utiliser la relation (3.64), on va donc partir
de l’hypothèse suivant laquelle la fréquence d’échantillonnage est un multiple d’ordre P de
νr :
νe = P νr (3.65)
77
CHAPITRE 3. Étude théorique du système
ce système, une valeur de 10 kHz convient parfaitement pour νe . En effet, elle est à la fois
largement supérieure à sa limite inférieure 8000+2×50 = 8100 Hz, et reste un multiple entier
de νr = 25 Hz. On obtient une période de P = 10000 25 = 400 et une réponse impulsionnelle de
la forme : µ ¶
π 2π
h(n, k) = 2, 5 cos k+ g(n − k) (3.68)
100 3
Les valeurs de g(n − k) sont connues par discrétisation de la réponse impulsionnelle continue
donnée à la figure 3.17(b).
On suppose dans cette partie que la fréquence d’échantillonnage a une valeur de 10 kHz,
et que la réponse impulsionnelle du système analysé est donnée par (3.68). En utilisant la
relation (3.5) et en définissant le retard r = n − k, on peut exprimer la réponse impulsionnelle
sous une autre forme :
³ p ´
g(r, k) = A cos 2π k + φ g(r) (3.69)
µP ¶
π 2π
= 2, 5 cos k+ g(r) (3.70)
100 3
La figure 3.19 représente une partie de cette réponse impulsionnelle, représentée en niveaux
de gris. On voit apparaı̂tre une période de P = 400 échantillons suivant l’axe des k, alors que
500
0.4
450
0.3
400
0.2
k (numéro d’échantillon)
350
300 0.1
250 0
200 −0.1
150
−0.2
100
−0.3
50
−0.4
0
0 50 100 150 200 250
r (numéro d’échantillon)
les réponses impulsionnelles du type de celle de la figure 3.17(b) sont situées suivant l’axe
des retards r. La valeur Pp = 400 2 = 200 étant un nombre entier, ce système à également
une période de 200 échantillons suivant l’axe des k, mais il reste considéré comme étant 400-
périodique. Comme prévu par (3.13), g(r, k) est P -périodique en k, et donc décomposable en
série de Fourier discrète par (3.15) et (3.16). En écrivant le cosinus présent dans (3.69) sous
sa forme exponentielle, on obtient immédiatement la série de Fourier recherchée :
A jφ p A P −p
g(r, k) = e g(r)ej2π P k + e−jφ g(r)ej2π P k
2 2
78
3.3. Application
Cette équation, par identification avec (3.15), permet de définir les réponses impulsionnelles
des deux systèmes LIT composant le système LVPT analysé :
A jφ
gp (r) = e g(r) (3.71)
2
A −jφ
gP −p (r) = e g(r) (3.72)
2
Il peut donc être représenté par une structure à modulation d’entrée à deux branches, comme
le montre la figure 3.20. On pourrait en déduire sa structure à modulation de sortie par la
p
ej2π P n
A jφ
2
e g(r)
ic (n) P −p
vc (n)
ej2π P n
A −jφ
2
e g(r)
79
CHAPITRE 3. Étude théorique du système
valeurs numériques précédemment définies. La réponse bifréquentielle est bien formée des
deux droites d’équation β = −α± Pp = −α± 200 1
, sur lesquelles sont situés les gains complexes
A jφ A −jφ
2 e ĝ(α) et 2 e ĝ(α). On retrouve, sur chacune de ces droites, l’allure en niveaux de gris
du module du gain complexe de la figure 3.17(a).
Finalement, la relation entrée-sortie en fréquence dans le cas de signaux réels est déduite de
−0.05
5
−0.1
4
−0.15
β
−0.2
2
−0.25 1
p P −p
ŷ(α) = ĝp (α)x̂(α − ) + ĝP −p (α)x̂(α − )
P P
A jφ p A P −p
= e ĝ(α)x̂(α − ) + e−jφ ĝ(α)x̂(α − )
2 P 2 P
A jφ p A p
= e ĝ(α)x̂(α − ) + e−jφ ĝ(α)x̂(α + )
2 P 2 P
Cette expression montre explicitement que la transformée de Fourier de la sortie ŷ(α) est la
somme de deux versions décalées et pondérées (ou filtrées) de celle de l’entrée. Les mêmes
conclusions peuvent être déduites de l’étude du module de la réponse bifréquentielle, en
80
3.3. Application
Calcul
p−1
ĝP −p (α − P
)
p
ĝP −p (α − P
)
Ĝ(α) =
P −p−1
ĝp (α − P
)
P −p
ĝp (α − P
)
P −1 P −1
ĝp (α − P
) ĝP −p (α − P
)
¤ 1 1 ¤
Cette matrice contient des fonctions de la variable α ∈ − 2P ; 2P . Les traitements décrits
à la section précédente peuvent lui être appliqués, mais les résultats sont alors eux-mêmes
fonctions de α (matrice inverse, décomposition en valeurs singulières). Afin d’obtenir des
81
CHAPITRE 3. Étude théorique du système
résultats facilement interprétables, on impose à la fréquence une valeur fixe. Une valeur
particulièrement intéressante est la fréquence nulle, pour laquelle la relation (3.40) entre les
réponses fréquentielles augmentées de l’entrée et de la sortie devient :
Or, en se basant sur la définition de la réponse fréquentielle augmentée d’un signal donnée
par (3.39), et en utilisant la périodicité de la TFr, on peut écrire, pour α = 0 :
x̂ (0) x̂ (0)
¡ ¢ ¡ ¢
x̂ − P1 x̂ PP−1
x̂(0) =
.. =
..
. .
¡ P −1 ¢ ¡1¢
x̂ − P x̂ P
Si le signal x est un signal discret P -périodique, ce vecteur regroupe les composantes non
nulles de sa transformée de Fourier (ou les coefficients de sa série de Fourier discrète) et
le caractérise entièrement. En effet, il permet de retrouver ses caractéristiques temporelles
aussi bien que fréquentielles. L’équation (3.75) donne donc la réponse d’un système LVPT
P -périodique de gain complexe matriciel Ĝ(α), à un signal d’entrée x P -périodique. y est
alors lui-même un signal P -périodique, caractérisé par le vecteur ŷ(0).
On retrouve exactement le même cas dans le problème de contrôle actif original. D’une part,
on vient de voir que le système à travers lequel on peut agir sur les vibrations de la machine est
P -périodique. D’autre part, on peut montrer (voir section 4.2) que pour atténuer le type de
vibrations envisagé dans le chapitre 1, l’entrée de ce système est nécessairement P -périodique.
Les résultats obtenus en analysant le gain complexe matriciel du système de la figure 3.16
pour α = 0 nous seront donc très utiles par la suite.
L’application numérique, à l’aide des valeurs précédemment définies, fournit une matrice Ĝ(0)
de taille 400 × 400 dont le module est représenté sur la figure 3.23. Comme on pouvait s’y
attendre, seules les sur- et sous-diagonales d’ordre p = 2 et P − p = 398 sont non nulles
et portent les valeurs des gains complexes ĝ2 (α) et ĝ398 (α) à des fréquences multiples de
1
400 , calculés par (3.71) et (3.72). Toutefois, les diagonales les plus éloignées de la diagonale
principale ne sont pas visibles du fait de leur très faible amplitude. L’équation (3.75) et la
forme du gain complexe matriciel obtenu montrent que chaque composante spectrale du signal
de sortie y est fonction de deux composantes spectrales de x. Ceci est dû à la modulation
d’amplitude sinusoı̈dale présente dans le système LVPT global, qui décale en fréquence le
signal d’entrée de deux valeurs différentes. On peut remarquer que si le signal modulant m(t)
avait eut une série de Fourier comprenant plus de deux coefficients, d’autres droites parallèles
seraient apparues dans la réponse bifréquentielle, ainsi que dans le gain complexe matriciel.
Le paragraphe suivant est consacré à l’analyse de ce gain complexe matriciel particulier, à
l’aide des outils théoriques introduits à la section 3.2.7.
82
3.3. Application
50
5
100
150 4
ligne
200
3
250
300
1
350
400 0
50 100 150 200 250 300 350 400
colonne
Analyse
L’information la plus importante est concentrée sur la diagonale principale de Σ̂(0) contenant
les valeurs singulières du gain complexe matriciel, rangées par ordre décroissant. Elles sont
représentées sur la figure 3.24, et on peut noter que la valeur singulière maximale est σ1 (0) ≃
9, 1 alors que la minimale est de l’ordre de σ400 (0) ≃ 1, 6 10−4 . L’interprétation de la norme
d’un gain complexe matriciel a montré que chaque valeur singulière élevée au carré, est
équivalente à un gain en énergie. Ce système LVPT peut donc appliquer un gain en énergie
dont l’ordre de grandeur varie de 10−8 à 100, suivant le signal d’entrée.
Inversion
Aucune valeur singulière n’étant nulle, le gain complexe matriciel du système inverse est
−1
donné par l’inverse exacte de Ĝ(0). On obtient alors la matrice Ĝinv (0) = Ĝ(0) , dont le
module est représenté sur la figure 3.25. Ses valeurs singulières sont les inverses de celles de
83
CHAPITRE 3. Étude théorique du système
10
0
0 50 100 150 200 250 300 350 400
50
500
100
400
150
ligne
200 300
250
200
300
100
350
400 0
50 100 150 200 250 300 350 400
colonne
Ĝ(0), et sont représentées sur la figure 3.26. Certaines prennent des valeurs très élevées, tout
comme le gain en énergie du système inverse, qui peut approximativement varier entre 10−2
et 108 .
Le système inverse exact est donc susceptible de présenter des gains en énergie très élevés,
ce qui peut être gênant pour certaines applications. Comme dans le cas des systèmes non-
inversibles, l’inverse généralisée peut fournir une solution à ce problème. En effet, le gain
complexe matriciel du système inverse peut être calculé en prenant l’inverse généralisée de
84
3.4. Conclusion
7000
6000
5000
4000
3000
2000
1000
0
0 50 100 150 200 250 300 350 400
celui du système direct, où on aura annulé les valeurs singulières situées en dessous d’un seuil
choisi. On obtient alors une approximation de Ĝinv (α) dont la valeur singulière maximum
est moins élevée que celle de l’inverse exacte, ce nouveau maximum étant fixé par l’inverse
du seuil précédent.
3.4 Conclusion
Ce chapitre est donc dédié à l’étude des systèmes LVPT, ce type de système étant présent
dans le problème de contrôle actif posé au chapitre 1.
Les outils théoriques nécessaires à leur étude sont fournis dans la section 3.2. L’étude de
leur réponse impulsionnelle (domaine temporel), ainsi que de leur réponse bifréquentielle
(domaine fréquentiel), permet de mettre à jour les effets qu’ils provoquent sur leur signal
d’entrée. Ils sont non seulement capables de le filtrer comme le ferait un système LIT, mais
aussi de le décaler en fréquence. Ces deux effets sont particulièrement visibles lorsque ces
systèmes sont représentés sous leur forme à modulation d’entrée ou de sortie, comme sur
les figures 3.4 et 3.5. Leur sortie peut donc être considérée comme une somme de versions
filtrées et décalées en fréquence du signal d’entrée, ce qui est résumé mathématiquement
par la relation fréquentielle (3.32). De plus, une matrice fonction de la fréquence appelée
«gain complexe matriciel», que l’on peut considérer comme l’extension du gain complexe
des systèmes LIT, apparaı̂t en étudiant les signaux propres de ces systèmes. Elle relie les
représentations fréquentielles «augmentées» des signaux d’entrée et de sortie, tout comme la
réponse en fréquence des systèmes LIT relie leur transformée de Fourier. Son analyse, réalisée
à l’aide de techniques d’algèbre matriciel comme la décomposition en valeurs singulières et
le calcul de normes, permet de définir des notions importantes pour le système. En effet, il
apparaı̂t que le carré de chaque valeur singulière de la matrice précédente représente un gain
d’entrée-sortie en énergie du système correspondant. De plus, si une de ses valeurs singulières
est nulle, le système n’est pas inversible.
L’application de ces outils théoriques sur un système analogue à celui rencontré au chapitre 1
85
CHAPITRE 3. Étude théorique du système
est réalisée dans la section 3.3. Son étude dans les domaines temporels et fréquentiels montre
qu’il provoque seulement deux décalages fréquentiels sur son signal d’entrée, qui sont dûs
à la modulation d’amplitude sinusoı̈dale présente dans le système global. De plus, dans le
cadre de notre application, le signal d’entrée peut être montré périodique de même période
que le système (voir section 4.2). L’étude du gain complexe matriciel correspondant est alors
limitée à la fréquence nulle, et la relation d’entrée-sortie en fréquence devient (3.75), qui sera
intensivement utilisée dans le chapitre 4. Enfin, pour l’exemple numérique choisi, le système
est montré inversible, mais le système inverse exact présente des gains entrée-sortie très élevés.
Une partie de ces résultats sera utilisée dans le chapitre suivant afin d’élaborer une solution
au problème de contrôle actif original.
86
Chapitre 4
4.1 Introduction
Le but de ce chapitre est d’élaborer les solutions théoriques au problème original de
contrôle actif, puis de tester leurs performances sur le système LVPT considéré.
Dans le chapitre 1, le transfert entre un courant de commande et un signal de vibrations a
été modélisé par un système LVPT. Le gain complexe matriciel de ce système a été étudié et la
relation qu’il permet d’écrire entre les représentations fréquentielles augmentées de ses signaux
d’entrée et de sortie à été développée à la section 3.3. La section 4.2 utilise cette relation pour
modéliser le système global à temps discret entre les trois courants statoriques et un ou
plusieurs capteurs accéléromètriques placés sur la carcasse statorique de la machine de test.
L’expression fréquentielle augmentée des signaux de vibration à compenser est finalement
déterminée.
La section 4.3 propose une fonction de coût quadratique pour minimiser la puissance
globale sur un ou plusieurs capteurs de vibrations. Elle permet de déterminer la commande
dite optimale minimisant cet indice de performance. Un algorithme récursif est finalement
développé afin de faire converger les valeurs de la commande vers cet optimal. Il permet aussi
à la commande de s’adapter à d’éventuelles variations des vibrations à minimiser au cours
du temps.
Enfin, la section 4.4 décrit les performances atteintes par cet algorithme sur le système
considéré, ainsi que la manière pratique de les implanter.
Le modèle du système global, représenté sur la figure 1.14, a été obtenu à la section 1.4
sous la forme d’un système LVPT à temps continu. Le système de contrôle actif à élaborer
87
CHAPITRE 4. Détermination et mise en œuvre des solutions
devant être implanté sur un processeur numérique, les signaux d’entrée et de sortie doivent
être échantillonnés. Le problème revient donc à déterminer les conditions que doit respecter
la fréquence d’échantillonnage νe , pour que la version échantillonnée de ce transfert soit bien
LVPT, afin de pouvoir lui appliquer les outils théoriques développés au chapitre 3.
On peut remarquer que la structure du transfert entre l’une quelconque des entrées ic1 (t),
ic2 (t) ou ic3 (t) et la sortie vc (t, θ) est la même que celui de la figure 3.16, étudié à la section
3.3. Les conditions d’échantillonnage précédemment obtenues peuvent donc être étendues au
système global. On rappelle que la fréquence réelle des variations de ce système est de pνr
Hertz, et que le gain complexe ĝ(ν, θ) du système LIT Gθ est supposé s’annuler après une
fréquence limite de + B2 Hertz. On a alors les conditions d’échantillonnage suivantes :
– si νe est supérieure à B +2pνr Hz, il n’y a pas de recouvrement spectral dans les signaux
échantillonnés,
– si νe est un multiple d’ordre P de la fréquence de rotation de la machine νr , le
système LVPT obtenu peut être considéré comme P -périodique, c’est à dire de fréquence
réduite P1 .
Il est donc nécessaire, comme à la section 3.3, d’avoir une fréquence d’échantillonnage suffi-
samment élevée, et un échantillonnage synchrone à la rotation de la machine, ce qui peut être
facilement réalisé en utilisant le signal de top-tour. Dans la suite, on suppose que νe vérifie
ces propriétés, et que le modèle numérique des transferts entre l’une quelconque des entrées
échantillonnées ic1 (n), ic2 (n) ou ic3 (n) et la sortie vc (n, θ), est celui d’un système LVPT dis-
cret P -périodique. On obtient le transfert à trois entrées et une sortie dont la structure est
représentée à la figure 4.1.
SYSTÈME LVPT
cos(pθ)
ic1 (n)
m(n, θ) = F cos(2π Pp n − pθ + ϕ) vn (n, θ)
2π
cos(pθ − 3 )
bc (n, θ) vc (n, θ)
ic2 (n) Gθ v(n, θ)
4π
cos(pθ − 3 )
ic3 (n)
Fig. 4.1: Modèle à temps discret du système global, discrétisé à la fréquence d’échantillonnage
νe = P νr
88
4.2. Modélisation mathématique de l’erreur
89
CHAPITRE 4. Détermination et mise en œuvre des solutions
la figure 3.22. Les seuls éléments non nuls qu’il contient sont situés sur ses pième et (P − p)ième
sur- et sous-diagonales, et sont les transformées de Fourier en fréquence réduite des signaux :
F j(−pθ+ϕ)
gp (r, θ) = e g(r, θ)
2
F −j(−pθ+ϕ)
gP −p (r, θ) = e g(r, θ)
2
c’est à dire :
F j(−pθ+ϕ)
ĝp (α, θ) = e ĝ(α, θ)
2
F −j(−pθ+ϕ)
ĝP −p (α, θ) = e ĝ(α, θ)
2
La représentation fréquentielle augmentée du signal de sortie v(n, θ) est donc donnée par :
v̂c (α, θ) = Ĝ(α, θ) · C(θ) · îc (α) = T̂(α, θ) · îc (α) (4.3)
où le gain complexe matriciel T̂(α, θ) regroupe les caractéristiques fréquentielles du système.
Cette relation constitue bien la relation entrée-sortie en fréquence du système LVPT discret
décrit à la figure 4.1. Elle peut être étendue au cas multicapteurs en utilisant de nouveau
la notation en longs vecteurs. En effet, si M capteurs sont placés sur la carcasse statorique
de la machine aux positions angulaires θ1 , . . ., θM , M signaux de sortie vérifiants (4.3) sont
alors disponibles. En regroupant les représentations fréquentielles augmentées de chacun de
ces signaux dans un même vecteur v̂c (α, Θ), où Θ = {θ1 , . . . , θM }, on obtient :
v̂c (α, θ1 ) Ĝ(α, θ1 ) · C(θ1 )
.. ..
v̂c (α, Θ) =
. =
. · îc (α) = T̂(α, Θ) · îc (α)
(4.4)
v̂c (α, θM ) Ĝ(α, θM ) · C(θM )
Cette nouvelle équation constitue bien une généralisation de (4.3), et relie les représentations
fréquentielles augmentées de tous les signaux de sortie du système LVPT considéré à celles
de ses entrées.
90
4.2. Modélisation mathématique de l’erreur
est nécessaire pour établir celle de la commande optimale, et peut facilement être déduite
dans le cas mono et multicapteurs, des relations (4.3), (4.4) et de la figure 4.1. En effet, en
se plaçant dans le cas monocapteur, les propriétés de la transformée de Fourier et la relation
(4.3) permettent d’écrire :
Cette relation est importante dans le sens où elle montre que les seules composantes fréquentielles-
des signaux de commande capables d’agir sur v̂(0, θ), et donc de le minimiser, sont conte-
nues dans le vecteur îc (0). Les signaux de commande optimaux ne contiennent donc des
composantes fréquentielles qu’aux fréquences multiples de P1 , ou d’une manière équivalente,
sont nécessairement P -périodiques. Le problème revient maintenant à déterminer le vecteur
îcopt (0), optimal dans le sens où il minimise un critère quadratique à déterminer, fonction de
l’erreur v̂(0, θ).
Ce résultat est directement généralisable au cas multicapteurs, pour lequel la relation (4.5)
devient :
v̂(α, Θ) = T̂(α, Θ) · îc (α) + v̂n (α, Θ)
avec
h i
v̂(α, Θ)T = v̂(α, θ1 )T ··· v̂(α, θM )T
h i
v̂n (α, Θ)T = v̂n (α, θ1 )T ··· v̂n (α, θM )T
91
CHAPITRE 4. Détermination et mise en œuvre des solutions
Cette dernière équation est tout à fait générale, et se réduit au cas monocapteur (4.6) en
fixant M = 1.
La section suivante utilise ce modèle mathématique de l’erreur afin de déterminer l’ex-
pression de la commande optimale qui, avec une puissance bornée, permet de minimiser le
vecteur (4.7).
Cette relation est l’extension vectorielle de (2.6), où le coût C devient proportionnel à la puis-
sance contenue dans le signal d’erreur aux fréquences multiples de P1 . En effet, en développant
cette expression en fonction des composantes du vecteur v̂(0, θ), on obtient :
P
X −1 ¯ µ ¶¯2
¯ l ¯
C= ¯v̂ , θ ¯¯
¯ P
l=0
92
4.3. Résolution théorique : commande optimale
Minimiser ce coût remplit donc le premier objectif correspondant à minimiser les composantes
vibratoires synchrones au phénomène de rotation. Par contre, cela ne donne aucun moyen
d’action sur les composantes de la commande îc (0). En s’inspirant des résultats du paragraphe
2.2.3, on peut définir une fonction de coût quadratique qui permet de contourner ce problème :
H
C = v̂(0, θ)H · v̂(0, θ) + ǫ îc (0) · îc (0) (4.9)
On retrouve une expression du même type que (2.7) qui peut être développée en fonction du
vecteur îc (0) à optimiser :
Comme au paragraphe 2.2.3, on peut interpréter physiquement cette expression, où les gran-
deurs ne sont considérées que pour un ensemble de fréquences constitué des seuls harmoniques
de P1 . Le premier terme correspond au coût en l’absence de signal de commande, sa valeur
étant proportionnelle à la puissance contenue dans le signal de perturbation à atténuer. En
utilisant l’équation (4.3), le terme situé sur la seconde ligne peut être vu comme la puissance
des contre-vibrations vc (n, θ) additionnée d’une fraction de la puissance de la commande.
Enfin, les termes de la troisième ligne correspondent à la puissance d’interaction entre les
contre-vibrations engendrées par la commande et les perturbations à atténuer. Ce sont les
seuls à pouvoir prendre des valeurs négatives et seront donc utilisés pour minimiser C.
Cette fonction de coût est une version vectorielle de la forme quadratique complexe (2.9).
Elle est appelée forme quadratique hermitienne [Nel92, pages 416-420], et peut être exprimée
comme une fonction réelle des 2P variables hréellesi définies
h pariles parties réelle et imaginaire
du vecteur îc (0), respectivement notées ℜ îc (0) et ℑ îc (0) . Elle prend alors une forme
de paraboloı̈de, donth les propriétés de courbure
i sont fixées par les valeurs propres de la ma-
H
trice hermitienne T̂(0, θ) · T̂(0, θ) + ǫ I3P . En effet, on peut montrer que la fonction C
est convexe et comporte donc un unique minimum si et seulement cette matrice est définie
positive, ce qui peut être imposé grâce à la constante réelle positive ǫ.
La minimisation d’une fonction de coût de la forme (4.10) dans le cadre d’un système de
contrôle actif a déjà été étudiée en détail dans [Nel92]. En effet, on obtient un coût identique
si le transfert secondaire sur lequel est appliqué le contrôle actif est un système LIT, multi-
entrées et multi-sorties. Les résultats exposés dans les sections suivantes étant directement
inspirés de ceux donnés dans ce papier, certaines démonstrations seront omises.
Les coordonnées du point où C prend sa valeur minimale définissent le vecteur com-
plexe îcopt (0) regroupant les valeurs optimales des composantes fréquentielles du signal de
93
CHAPITRE 4. Détermination et mise en œuvre des solutions
Cette solution est l’équivalent de la solution de Wiener obtenue en filtrage optimal LIT clas-
sique. Elle est fonction des caractéristiques fréquentielles du système contenues dans T̂(0, θ),
et de celles du signal de perturbation contenues dans v̂n (0, θ).
L’influence du paramètre ǫ sur la valeur de cet optimal peut être étudiée en exprimant
une grandeur proportionnelle à la puissance du signal de commande :
H
h H
i−2 H
îcopt (0) · îcopt (0) = v̂n (0, θ)H · T̂(0, θ) · T̂(0, θ) · T̂(0, θ) + ǫ I3P · T̂(0, θ) · v̂n (0, θ)
(4.13)
Cette grandeur est inversement proportionnelle au carré de ǫ, et diminue lorsque ce paramètre
augmente. Pour que la puissance de la commande optimale reste inférieure à un seuil fixé par
l’utilisateur, ǫ doit donc être choisi suffisamment élevé.
Ce même paramètre joue également sur la valeur minimale de la fonction de coût, qui
peut être déterminée en reportant l’expression du vecteur optimal (4.12) dans (4.10) :
· h i−1 ¸
H H H
Cmin = v̂n (0, θ) · IP − T̂(0, θ) · T̂(0, θ) · T̂(0, θ) + ǫ I3P · T̂(0, θ) · v̂n (0, θ) (4.14)
H
On se place dans le cas où la matrice T̂(0, θ) · T̂(0, θ) inversible. Si ǫ est choisi nul, la
commande optimale annule alors complètement la fonction de coût, au prix d’une puissance
de commande éventuellement très élevée. Si le système expérimental ne peut supporter une
telle puissance, le paramètre ǫ est imposé strictement positif. La valeur minimale du coût
n’est alors plus nulle, et augmente pour atteindre v̂n (0, θ)H · v̂n (0, θ) lorsque ǫ est infini.
94
4.3. Résolution théorique : commande optimale
Les grandeurs indicées l correspondent aux valeurs de ces grandeurs pour la lième itération.
Cette relation permet donc de prédire la valeur du vecteur îc (0) à la prochaine itération en
fonction de sa valeur courante îc (0)l , et a le même comportement que la relation récursive
du paragraphe 2.2.4. Une expression plus pratique est obtenue en notant que le gradient de
la fonction de coût (4.11) peut s’écrire, à l’aide de (4.6) :
h i H
∇C îc (0) = 2T̂(0, θ) · v̂(0, θ) + 2ǫ îc (0) (4.16)
95
CHAPITRE 4. Détermination et mise en œuvre des solutions
Les propriétés de convergence de cet algorithme peuvent être étudiées de la même manière
qu’au paragraphe 2.2.4. On définit un vecteur E l représentant, à l’itération l, l’erreur entre
le vecteur commande courant et sa valeur optimale :
On obtient alors une relation de récurrence sur cette erreur en reportant (4.18), (4.6) et (4.12)
dans (4.17) :
h ³ H
´i
E l+1 = I3P − 2µ T̂(0, θ) · T̂(0, θ) + ǫ I3P · E l
h ³ H
´il
= I3P − 2µ T̂(0, θ) · T̂(0, θ) + ǫ I3P · E0
La condition à vérifier par le pas d’adaptation µ pour que ce vecteur converge vers zéro,
ainsi que l’influence que peut avoir le paramètre ǫ sur cette convergence ont été étudiées en
détail dans [Nel92] à l’aide de la décomposition en valeurs singulières de la matrice T̂(0, θ).
Il apparaı̂t tout d’abord que le pas d’adaptation µ doit avoir une valeur comprise dans
l’intervalle :
1
0<µ< (4.19)
σ̂1 (0, θ)2 + ǫ
où σ̂1 (0, θ) est la valeur singulière la plus élevée de T̂(0, θ).
Cette condition assure la convergence du vecteur de commande îc (0)l vers sa valeur optimale
îc (0)opt , quelle que soit sa valeur initiale îc (0)0 .
De plus, il est montré que la convergence globale du vecteur E l esth due à celle de différents i
H
modes indépendants, tous associés à une valeur propre de la matrice T̂(0, θ) · T̂(0, θ) + ǫ I3P .
Plus cette valeur propre est faible, plus le mode qui lui est associé convergera lentement. Or,
chacune de ces valeurs propres est formée par la somme d’une valeur propre de la matrice
H
T̂(0, θ) · T̂(0, θ) et du paramètre ǫ. Ce dernier, s’il est correctement choisi, permet donc de
donner une borne inférieure positive aux valeurs propres globales les plus faibles, sans avoir
d’effet notable sur les plus élevées. Le temps de convergence de l’algorithme est ainsi amélioré
par ǫ, puisqu’il accélère les modes les plus lents sans trop influencer les modes rapides.
Les divers résultats qui viennent d’être obtenus dans le cas monocapteur vont être généra-
lisés au cas multicapteurs dans le paragraphe suivant.
96
4.4. Implantation des solutions et étude des performances obtenues
Le second terme de cette somme reste équivalent à celui du cas monocapteur (4.9). Par
contre, le premier terme est maintenant formé de la somme des puissances contenues dans les
M signaux d’erreur, aux fréquences multiples de P1 .
Par la même méthode qu’au paragraphe précédent, on peut montrer que le vecteur optimal
minimisant ce coût prend la forme :
h H
i−1 H
îcopt (0) = − T̂(0, Θ) · T̂(0, Θ) + ǫ I3P · T̂(0, Θ) · v̂n (0, Θ) (4.21)
Les signaux de commande dont les caractéristiques fréquentielles vérifient la relation précédente
minimisent donc en même temps la puissance des M signaux d’erreur issus des M capteurs.
La relation de récurrence permettant de faire converger le vecteur îc (0) vers sa valeur
optimale devient alors :
H
îc (0)l+1 = (1 − 2µǫ)îc (0)l − 2µT̂(0, Θ) · v̂(0, Θ)l (4.22)
Les expressions trouvées dans le cas multicapteurs sont donc tout à fait équivalentes à
(4.9), (4.12) et (4.17), obtenues dans le cas monocapteur. Il faut toutefois utiliser les matrices
et les vecteur «longs», définis à la fin du paragraphe (4.2.3). De plus, la valeur du pas d’adap-
tation vérifie le même type de contrainte que celle définie parh (4.19), où σ̂1 (0, θ) est remplacée
i
H
par σ̂1 (0, Θ), la plus grande valeur singulière de la matrice T̂(0, Θ) · T̂(0, Θ) + ǫ I3P . En-
fin, le paramètre ǫ garde les mêmes propriétés, et permet de limiter la puissance des signaux
de commande et d’accélérer la convergence de l’algorithme (4.22), tout en limitant les per-
formances optimales du système de contrôle actif.
L’équation générale (4.22) donne donc l’algorithme récursif permettant de faire converger
le vecteur îc (0)l , contenant les composantes fréquentielles de la commande, vers sa valeur
optimale (4.21). La section suivante est dédiée à l’étude de l’implantation pratique d’un tel
algorithme, et des performances qu’il atteint lorsqu’il est appliqué à des signaux synthétiques
ou réels de vibrations.
97
CHAPITRE 4. Détermination et mise en œuvre des solutions
vibrations
amplificateur
capteur
de
angulaire
puissance
capteurs
de processeur
multiplieur
vibrations
de fréquence
par P DSP
fréquence réduite P1 donné par la figure 4.1. À chaque instant d’échantillonnage, on dispose
finalement dans la mémoire du processeur numérique, des échantillons courants et précédents
des M signaux d’erreur v(n, θ1 ), . . . , v(n, θM ), échantillonnés d’une manière synchrone à la
rotation de la machine.
La programmation de l’algorithme général (4.22) peut maintenant être abordée. Cette
récurrence s’effectue dans le domaine fréquentiel, et nécessite le calcul de transformées de Fou-
rier discrètes (TFD), portant sur des blocs d’échantillons. La longueur optimale de ces blocs
peut être déterminée en notant qu’il est nécessaire, afin de réactualiser le vecteur commande
îc (0)l , de connaı̂tre le vecteur v̂(0, Θ)l . Celui-ci contient les P composantes fréquentielles de
chaque signal d’erreur, aux fréquences réduites multiples de P1 . Ces coefficients peuvent être
calculés avec une grande précision par une TFD portant sur une longueur de bloc propor-
tionnelle à P . Les blocs de traitement sont donc choisis de longueur égale à la période du
système LVPT, c’est à dire avec P échantillons.
L’algorithme (4.22) consiste donc à acquérir P échantillons de chaque signal d’erreur pour
former M blocs de traitement de longueur P et d’indice l. Une fois ces blocs complets, on leur
applique une TFD, dont le résultat permet de construire le vecteur v̂(0, Θ)l . On en déduit
ensuite le vecteur îc (0)l+1 par l’équation de récurrence (4.22). Les P échantillons des trois
98
4.4. Implantation des solutions et étude des performances obtenues
signaux de commande correspondants peuvent alors être déterminés par TFD inverse sur ce
vecteur. Ils sont envoyés à l’entrée du système LVPT en même temps que l’acquisition des
blocs suivants, d’indice l + 1. Le chronogramme de la figure 4.3 résume ce traitement, où
chaque instant d’échantillonnage est symbolisé par un cercle de couleur différente suivant le
bloc auquel il appartient.
bloc l bloc l + 1
instants
1 2 P 1 2 P d’échantillonnage
99
CHAPITRE 4. Détermination et mise en œuvre des solutions
pensable pour déterminer l’expression exacte de T̂(0, Θ) mais il est difficile à connaı̂tre pour
tout point de fonctionnement. En pratique, il sera donc nécessaire de l’estimer en étudiant le
champ magnétique naturel de la machine pour le point de fonctionnement considéré, et sera
supposé connu pour les différentes simulations qui suivent.
Nous allons tout d’abord nous placer dans le cas où le seul courant de commande ic1 (t)
parcourant la bobine statorique numérotée 1 à la figure 1.13, est utilisé pour atténuer les
vibrations mesurées par un seul accéléromètre fixé sur le stator. Sa position angulaire est
choisie telle qu’il soit situé dans l’axe de cette bobine, c’est à dire θ = 0. On a alors un
système LVPT à une entrée et une sortie. Afin de pouvoir utiliser les résultats numériques
obtenus à la section 3.3, nous nous plaçons dans les mêmes conditions de fonctionnement :
– la machine synchrone est supposée avoir p = 2 paires de pôles, et tourner à une fréquence
de νr = 25 Hz,
– le signal de modulation m(t, θ) = m(t, 0) est donné par la relation (3.59),
– le gain complexe ĝ(ν, θ) = ĝ(ν, 0) modélisant la réponse statorique en θ = 0 est supposé
avoir été identifié au préalable, et ses caractéristiques sont données par la figure 3.17,
– La fréquence d’échantillonnage est de 10 kHz.
Le modèle du transfert LVPT discret obtenu est alors exactement celui de la section 3.3.
La matrice T̂(0, Θ) = T̂(0, 0) correspondante à ce système se réduit donc au gain complexe
matriciel déterminé au paragraphe 3.3.5, dont le module est représenté à la figure 3.23.
De plus, on se trouve maintenant dans le domaine discret, et toutes les fréquences sont
données sous leur forme réduite, le retour à la fréquence réelle exprimée en Hertz se faisant
par une simple multiplication par νe = 10000. La perturbation à atténuer est choisie la plus
simple possible, et comporte une seule sinusoı̈de d’amplitude unitaire et de phase nulle. Sa
fréquence réduite est fixée à 0,15, qui est l’harmonique d’ordre 60 de celle du système LVPT
1
(0, 15 = 60 × 400 ). Elle correspondrait à la présence d’une raie de fréquence 1,5 kHz dans le
spectre vibratoire de la machine. Ces caractéristiques permettent de déterminer la valeur du
vecteur v̂n (0, Θ) = v̂n (0, 0), dont une seule composante est alors non nulle et égale à 12 .
Performances optimales
Nous allons tout d’abord étudier les performances optimales que peut atteindre ce système
de contrôle actif, c’est à dire celles obtenues lorsque l’algorithme récursif élaboré précédemment
a convergé. Pour cela, on utilise directement les équations du paragraphe 4.3.2. La relation
(4.12) permet de calculer le vecteur optimal îcopt (0) qui, une fois placé dans (4.6), donne la va-
leur de v̂(0, 0) obtenue pour une telle commande en θ = 0. On peut tirer de ces deux vecteurs
les valeurs complexes des composantes fréquentielles optimales des signaux de commande et
d’erreur, aux fréquences réduites multiples de P1 = 400 1
. Les modules de ces valeurs calculés
pour ǫ = 0, 05 sont représentés sur la figure 4.4 en fonction de la fréquence réduite. Le module
100
4.4. Implantation des solutions et étude des performances obtenues
0.6 0.6
perturbation perturbation
commande erreur
0.5 0.5
amplitude linéaire
amplitude linéaire
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0 0
0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.1 0.12 0.14 0.16 0.18 0.2 0.22
fréquence réduite fréquence réduite
(a) |v̂n (0, 0)|=tirets, |îcopt (0)|=traits pleins (b) |v̂n (0, 0)|=tirets, |v̂(0, 0)|=traits pleins
des composantes fréquentielles du signal de perturbation est représenté en tirets sur les fi-
gures 4.4(a) et 4.4(b), où l’on peut vérifier que sa seule valeur non nulle se situe à la fréquence
réduite 0,15. Celui du signal de commande optimale est donné en traits pleins sur la figure
4.4(a). Il s’agit d’un ensemble de sinusoı̈des réparties autour de 0,15 et dont les fréquences
4 2
sont espacées de 400 . Elles vont toutes être décalées de ± 400 par la modulation d’amplitude
sinusoı̈dale que comporte le système LVPT. Certaines composantes pourront alors interagir
destructivement avec celle de la perturbation, mais il y aura également création de raies pa-
4
rasites espacées de 400 autour de la fréquence visée dans le signal d’erreur. C’est le résultat
que montre la figure 4.4(b), où l’on a représenté le module des composantes fréquentielles du
signal d’erreur. On peut remarquer que la raie spectrale perturbatrice a bien été atténuée,
au prix de l’apparition de plusieurs raies parasites.
Les allures temporelles de ces différents signaux peuvent être obtenues en appliquant une
TFD inverse sur les vecteurs précédents. Les résultats sont exposés sur la figure 4.5 où une
sinusoı̈de de même fréquence et de même phase que le signal sinusoı̈dal de modulation m(n, 0)
est représentée. Cette courbe est donc l’image des variations du système LVPT global et per-
met une interprétation plus aisée des résultats. La figure 4.5(a) montre que la commande
essaie de «contrer» l’effet de la modulation d’amplitude sinusoı̈dale en prenant une ampli-
tude inversement proportionnelle à m(n, 0). Toutefois, lorsque la valeur de la modulation est
très proche de zéro, la commande ne peut devenir infinie, la valeur positive du paramètre ǫ
l’obligeant à garder une puissance bornée. Elle a alors tendance à s’annuler, et ne cherche
plus à atténuer la perturbation. Ce phénomène se retrouve dans le signal d’erreur tracé sur la
figure 4.5(b). En effet, celui-ci est faible lorsque la modulation d’amplitude prend des valeurs
importantes, mais il devient beaucoup plus grand autour des zones où m(n, 0) est proche de
zéro. Ce sont ces «bouffées» espacées de 100 échantillons dans le signal d’erreur, qui corres-
1
pondent aux composantes fréquentielles parasites espacées d’une fréquence réduite de 100 à la
101
CHAPITRE 4. Détermination et mise en œuvre des solutions
2
3 variations système variations système
commande erreur
1.5
perturbation
2
1
amplitude linéaire
amplitude linéaire
1 0.5
0 0
−0.5
−1
−1
−2
−1.5
0 50 100 150 200 250 300 0 50 100 150 200 250 300
numéro d’échantillon numéro d’échantillon
(a) ic1 opt (n) (b) v(n, 0)=traits pleins, vn (n, 0)=pointillés
Fig. 4.5: Allures temporelle des signaux optimaux de commande et d’erreur calculées pour
ǫ = 0, 05
figure 4.4(b). Elles seront présentes dans tous les résultats ultérieurs, et nous verrons qu’elles
ne peuvent être évitées qu’à l’aide d’une commande de puissance très importante, inadaptée
à la présente application. Ceci est bien entendu lié au fait que les valeurs singulières du gain
complexe matriciel T̂(0, 0) de ce système LVPT sont très faibles. Celles du système inverse qui
intervient dans l’expression de la commande optimale annulant complètement l’erreur sont
donc très importantes, ce qui engendre une commande optimale de puissance trop élevée pour
être utilisée.
On peut finalement quantifier les performances de cette méthode en calculant la valeur mi-
nimale de la fonction de coût Cmin , ainsi que la puissance du signal d’erreur optimal Pv .
Lorsqu’il n’y a aucun contrôle, ces grandeurs sont toutes deux égales à la puissance du signal
de perturbation notée Pvn , et peuvent donc être exprimées en pourcentage relativement à
cette valeur. Pour ǫ = 0, 05, on obtient Cmin ≃ 25 % de Pvn et Pv ≃ 11 % de Pvn . La puis-
sance de la commande optimale Pic1 opt nécessaire peut elle-même être calculée relativement à
la puissance de la perturbation, ce qui donne Pic1 opt ≃ 270 % de Pvn . Ces résultats montrent
que pour ce cas particulier, il est nécessaire d’avoir une commande dont la puissance est 2,7
fois celle de la perturbation afin d’atténuer 90 % de cette dernière.
Afin de vérifier l’influence du paramètre ǫ sur les performances obtenues, on peut calculer,
pour les mêmes conditions que précédemment, les grandeurs Cmin , Pv et Pic1 opt en fonction de
ǫ. Les résultats obtenus sont représentés sur la figure 4.6, où 4.6(a) montre le coût minimum
ainsi que la puissance du signal d’erreur et 4.6(b) la puissance de la commande optimale pour
ǫ variant de 0 à 2. Ces courbes montrent que si la valeur de ce paramètre est proche de zéro,
le coût minimal et la puissance de l’erreur seront eux-mêmes très faibles, au détriment d’une
puissance de commande élevée. Inversement, si ǫ devient trop important, la puissance de la
commande restera faible, mais la perturbation ne sera alors quasiment pas atténuée, ce qui
correspond bien aux conclusions du paragraphe 4.3.2.
102
4.4. Implantation des solutions et étude des performances obtenues
1 3
cout minimum
0.9 puissance erreur
2.5
0.8
0.7
2
echelle lineaire
echelle lineaire
0.6
0.5 1.5
0.4
1
0.3
0.2
0.5
0.1
0 0
0 0.5 1 1.5 2 0 0.5 1 1.5 2
ε ε
(a) Cmin et puissance du signal d’erreur Pv (b) Puissance de la commande optimale Pic1 opt
30
variations systeme
25 commande
20
15
amplitude lineaire
10
5
−5
−10
−15
−20
0 50 100 150 200 250 300
numero d’echantillon
différences se situent autour des zones où la modulation d’amplitude est proche de zéro. En
effet, pour ǫ nul, la commande optimale peut devenir aussi importante que nécessaire et tend
alors à prendre de très grandes valeurs. C’est donc principalement en ces endroits que le
paramètre ǫ limite l’amplitude de la commande, et permet d’obtenir un signal de puissance
bornée.
On rappelle que le module |ĝ(ν, 0)| du gain complexe du transfert LIT présent dans le
système LVPT global est donné par la figure 3.17(a). Il comporte une résonance importante
à 500 Hz, une antirésonance à 1 kHz et une autre résonance plus faible à 1,5 kHz. Pour
comprendre l’influence de sa forme sur les performances optimales de compensation, on définit
un signal de perturbation formé de quatre sinusoı̈des de fréquences différentes. La première
est située juste au niveau de la plus grande résonance, avec une fréquence réduite de 0,05.
103
CHAPITRE 4. Détermination et mise en œuvre des solutions
La seconde sinusoı̈de de fréquence réduite 0,1 se trouve dans l’antirésonance. Enfin, les deux
dernières ont des fréquences réduites très rapprochées (0,15 et 0,1575) et se trouvent au niveau
de la seconde résonance. La figure 4.8 expose les résultats obtenus pour une telle perturbation,
où les graphiques du haut donnent les composantes fréquentielles des signaux et ceux du bas
leurs allures temporelles. Les figures de gauche concernent le signal optimal de commande et
celles de droite sont dédiées à celui de l’erreur. On peut noter sur 4.8(b) que les composantes
0.7 0.7
perturbation perturbation
commande erreur
0.6 0.6
0.5 0.5
amplitude lineaire
amplitude lineaire
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0 0
0 0.05 0.1 0.15 0.2 0.25 0 0.05 0.1 0.15 0.2 0.25
frequence reduite frequence reduite
(a) |v̂n (0, 0)|=tirets, |îcopt (0)|=traits pleins (b) |v̂n (0, 0)|=tirets, |v̂(0, 0)|=traits pleins
8 6
variations systeme variations systeme
commande 5 erreur
6
perturbation
4
4
amplitude lineaire
amplitude lineaire
3
2 2
0 1
0
−2
−1
−4
−2
−6 −3
0 50 100 150 200 250 300 0 50 100 150 200 250 300
numero d’echantillon numero d’echantillon
(c) ic1 opt (n) (d) v(n, 0)=traits pleins, vn (n, 0)=pointillés
Fig. 4.8: Caractéristiques des signaux optimaux de commande et de erreur calculées pour
ǫ = 0, 05
fréquentielles de la perturbation situées sur des résonances sont fortement atténuées, alors
que celle située dans l’antirésonance ne diminue que très peu. On la retrouve présente dans
4.8(d), où le signal d’erreur est principalement constitué d’une sinusoı̈de de période égale à 10
échantillons. En effet, cette composante ne peut être diminuée qu’au prix d’une commande de
forte puissance, ce qui n’est pas autorisé par la valeur choisie pour le paramètre ǫ. De plus, on
remarque en 4.8(c) que la commande se comporte de la même manière que précédemment,
et tend à s’annuler lorsque la modulation d’amplitude devient faible. On peut enfin noter
104
4.4. Implantation des solutions et étude des performances obtenues
que ce type de commande peut agir sur plusieurs composantes fréquentielles à la fois, pourvu
qu’elles soient harmoniques de la fréquence fondamentale des variations du système LVPT.
La dernière variable dont on peut étudier l’influence sur les performances de cette méthode
de compensation est la position angulaire du capteur θ. Pour simplifier l’interprétation des
résultats, la perturbation est de nouveau supposée sinusoı̈dale de fréquence réduite 0,15 et
le transfert ĝ(ν, θ) est choisi identique au précédent, quelle que soit la position angulaire θ.
Les courbes de la figure 4.9 sont tracées en coordonnées polaires et ont toutes été calculées
pour ǫ = 0, 05. Elles représentent le coût et la puissance du signal d’erreur en fonction de θ
(figure 4.9(a)) ainsi que la puissance de la commande en fonction de la même variable (figure
4.9(b)). Elles peuvent être interprétées comme une coupe transversale de la machine, sur
0 0
30 330 30 330
60 300 60 300
(a) Cmin (trait plein) et puissance du signal d’er- (b) Puissance de la commande optimale Pic1 opt
reur Pv (tirets)
105
CHAPITRE 4. Détermination et mise en œuvre des solutions
0 0
30 330 30 330
60 300 60 300
(a) Cmin (trait plein) et puissance du signal d’er- (b) Pic1 opt (trait plein), Pic2 opt (tirets),Pic3 opt
reur Pv (tirets) (pointillés), Pic opt (tirets-points)
voit que les performances obtenues deviennent constantes quelle que soit la position angulaire
θ, et que cette méthode de compensation ne comporte alors plus de limitations spatiales. La
figure 4.10(b) représente, en fonction de θ, les puissances respectives de chaque signal de
commande Pic1 opt , Pic2 opt et Pic3 opt , ainsi que leur somme Pic opt exprimées elles aussi en
pourcentage relativement à Pvn . Picopt reste constante quel que soit l’endroit où se trouve le
capteur mesurant les vibrations à minimiser, mais les autres courbes mettent en évidence le
rôle joué par chaque signal de commande. Par exemple, pour θ = 0 c’est Pic1 opt qui est la
plus importante et qui génère la plupart des contre-vibrations nécessaires à l’atténuation des
perturbations. Au contraire, pour θ = π4 , Pic1 opt s’annule et ce sont les deux autres signaux de
commande qui prennent le relais afin de minimiser le signal d’erreur. C’est donc le déphasage
spatial des trois bobines statoriques qui leur permet d’agir à des endroits différents sur le
stator et de garder des performances de compensation constantes. Dans la suite, les trois
commandes seront toujours utilisées en même temps, et les performances obtenues seront
indépendantes de la position angulaire du capteur θ.
Après avoir étudié les performances optimales que peut atteindre cette méthode de contrôle
actif, nous allons nous attarder sur le comportement de l’algorithme récursif décrit au para-
graphe 4.3.3, qui permet au signal de commande de converger vers la solution désirée.
106
4.4. Implantation des solutions et étude des performances obtenues
Algorithme récursif
Le comportement de l’algorithme récursif va tout d’abord être étudié dans le cas où il
réactualise les composantes fréquentielles de trois signaux de commande afin d’atténuer, à
la position angulaire θ = 0, une perturbation sinusoı̈dale de fréquence réduite 0,15. Le gain
complexe ĝ(ν, 0) garde la même valeur et le paramètre ǫ est fixé à 0,05. Les performances
globales doivent donc être les mêmes que celles de la figures 4.10 pour la position angulaire
considérée. Le pas d’adaptation µ de l’algorithme (4.17) a été choisi de telle sorte qu’il y
ait convergence et prend la valeur 0,08. Les courbes de convergence, généralement appelées
«courbes d’apprentissage», sont données sur la figure 4.11. Elles représentes la valeur du coût
C, celle de la puissance du signal d’erreur Pv ainsi que la somme des puissances des commandes
Pic opt non pas en fonction des échantillons, mais en fonction de la période du système LVPT
sur laquelle ces grandeurs ont été calculées. Chaque unité de l’axe des abscisses représente
alors une période, ce qui correspond à 400 échantillons temporels. La forme de ces courbes,
200
150
n
% de Pv
100 cout
puissance erreur
puissance commandes
50
0
0 50 100 150 200 250
numero de periode du systeme LVPT
ainsi que les valeurs vers lesquelles ces grandeurs convergent correspondent à ce qui était
attendu. De plus, la figure 4.12 montre les composantes fréquentielles de l’erreur ainsi que
celles des trois signaux de commande obtenues pour la période numérotée 250, c’est à dire une
fois que l’algorithme a convergé. Le vecteur îcopt (0) regroupant les composantes fréquentielles
des trois commandes converge donc bien vers sa valeur optimale. De plus, on peut vérifier
que le signal de commande le plus important pour cette position angulaire est ic1 (n), comme
prévu par la courbe 4.10(b).
L’influence du pas d’adaptation µ sur la vitesse de convergence est illustré par la figure
4.13. L’évolution de la valeur du coût en fonction de la période du système LVPT y est
représentée pour trois valeurs différentes de µ, et on peut vérifier que le temps de convergence
de l’algorithme étudié est inversement proportionnel à ce paramètre.
L’algorithme récursif élaboré au paragraphe 4.3.3 se comporte donc bien comme prévu dans
le cas de perturbations sinusoı̈dales synchrones aux variations du système LVPT. Nous allons
maintenant étudier ce qu’il se passe en présence de bruit large bande et de sinusoı̈des non
synchrones dans le signal de perturbation à atténuer.
107
CHAPITRE 4. Détermination et mise en œuvre des solutions
0.6 0.6
perturbation perturbation
erreur commande 1
0.5 0.5
amplitude lineaire
amplitude lineaire
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0 0
0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.1 0.12 0.14 0.16 0.18 0.2 0.22
frequence reduite frequence reduite
(a) |v̂n (0, 0)|=tirets, |v̂(0, 0)|=traits pleins (b) |v̂n (0, 0)|=tirets, |îc1opt (0)|=traits pleins
0.6 0.6
perturbation perturbation
commande 2 commande 3
0.5 0.5
amplitude lineaire
amplitude lineaire
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0 0
0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.1 0.12 0.14 0.16 0.18 0.2 0.22
frequence reduite frequence reduite
(c) |v̂n (0, 0)|=tirets, |îc2opt (0)|=traits pleins (d) |v̂n (0, 0)|=tirets, |îc3opt (0)|=traits pleins
100
µ=0.005
90 µ=0.02
µ=0.08
80
70
n
% de Pv
60
50
40
30
20
10
0 50 100 150 200 250
numero de periode du systeme LVPT
108
4.4. Implantation des solutions et étude des performances obtenues
Les conditions de simulation restent complètement identiques aux précédentes, mis à part
le signal de perturbation. Il est constitué d’un bruit aléatoire gaussien large bande ajouté à
deux sinusoı̈des pures, l’une de fréquence 0,15 qui est un multiple de la fréquence du système
LVPT, l’autre non synchrone de fréquence 0,131. Le rapport entre la puissance de la sinusoı̈de
représentant les composantes synchrones et celle des autres composantes de ce signal (le bruit
et la sinusoı̈de non synchrone) peut être interprété comme le rapport signal à bruit (RSB).
Les puissances de chaque sinusoı̈de ainsi que celle du bruit étant choisies égales, ce rapport
est de 12 ou, ce qui est équivalent, de -3 décibels. La figure 4.14 représente les courbes d’ap-
prentissage de l’algorithme pour µ = 0, 08, qui sont du même type que celles de la figure
4.11. On remarque que la puissance du signal d’erreur, une fois que l’algorithme a convergé,
100
90
80
70
60
n
% de Pv
50
40
30
20 cout
puissance erreur
10
puissances commandes
0
0 50 100 150 200 250
numero de periode du systeme LVPT
Fig. 4.14: Courbes d’apprentissage C (trait plein), Pv (tirets) et Picopt (pointillés) pour une
perturbation sinusoı̈dale bruitée
représente environ 70 % de sa valeur initiale. Cette méthode de contrôle n’agissant que sur
les composantes synchrones, les 30 % qui ont été atténués doivent normalement apparte-
nir exclusivement à la sinusoı̈de de fréquence 0,15. Ceci peut être vérifié sur la figure 4.15,
représentant suivant une échelle logarithmique les densités spectrales de puissance des si-
gnaux de perturbation, de commande et d’erreur. En effet, la présence d’une composante
aléatoire dans la perturbation impose de calculer les densités spectrales de puissance comme
elles ont été définies au paragraphe 1.3.3. On remarque sur la figure 4.15(a) que la commande
n’est quasiment pas perturbée par la présence du bruit et de la sinusoı̈de non synchrone.
En effet, ses composantes fréquentielles sont réparties autour de la seule fréquence 0,15 et
ne prennent en compte ni le bruit, ni la sinusoı̈de non synchrone. La figure 4.15(b) prouve
que cette dernière se retrouve inchangée dans le spectre du signal d’erreur, et que ce traite-
ment n’atténue que les composantes synchrones aux variations du système LVPT. Une légère
amplification du bruit se produit toutefois au niveau de la forte résonance que comporte le
système LIT ĝ(ν, 0). Pour éviter ce genre de problèmes en pratique, il faudra donc limiter la
gamme de fréquence sur laquelle porte le contrôle actif, en filtrant le signal d’erreur par un
filtre dont la bande passante ne correspond pas aux fréquences de forte résonance du système.
Après avoir vu que la présence de bruit aléatoire large bande et de sinusoı̈des non syn-
109
CHAPITRE 4. Détermination et mise en œuvre des solutions
0 0
10 10
spectre perturbation spectre perturbation
spectre commande 1 spectre erreur
−1 −1
10 10
echelle logarithmique
−2
−2
10 10
−3
−3 10
10
−4
−4 10
10 0.05 0.1 0.15 0.2
0.05 0.1 0.15 0.2 frequence reduite
(a) Spectre du signal de commande (trait plein) et (b) Densité spectrale de puissance du signal d’er-
de perturbation (pointillés) reur (trait plein) et de perturbation (pointillés)
Fig. 4.15: Densité spectrale de puissance des signaux de commande et d’erreur pour une
perturbation sinusoı̈dale bruitée
chrones aux variations du système LVPT dans les vibrations naturelles de la machine ne gênait
par outre mesure le calcul de la commande optimale, l’algorithme (4.17) va être appliqué à
des signaux vibratoires mesurés sur une machine synchrone en rotation.
110
4.4. Implantation des solutions et étude des performances obtenues
250
cout
puissance erreur
puissance commandes
200
150
% de Pv
100
50
0
0 20 40 60 80 100
numero de periode systeme LVPT
Fig. 4.16: Courbes d’apprentissage C (trait plein), Pv (tirets) et Pic opt (pointillés) pour une
perturbation vibratoire réelle
La figure 4.17 montre les densités spectrales de puissance des signaux de perturbation, de
commande et d’erreur une fois que l’algorithme récursif a convergé. Comme prévu, les com-
30 30
perturbation perturbation
25 commande 1 25 erreur
20 20
15 15
10 10
decibels
decibels
5 5
0 0
−5 −5
−10 −10
−15 −15
−20 −20
1200 1400 1600 1800 2000 2200 1200 1400 1600 1800 2000 2200
frequence reelle (Hertz) frequence reelle (Hertz)
(a) Densité spectrale de puissance du signal de com- (b) Densité spectrale de puissance du signal d’er-
mande (trait plein) et de perturbation (pointillés) reur (trait plein) et de perturbation (pointillés)
Fig. 4.17: Densités spectrales de puissance des signaux de commande et d’erreur pour une
perturbation vibratoire réelle
111
CHAPITRE 4. Détermination et mise en œuvre des solutions
1.5 1.5
amplitude lineaire
amplitude lineaire
1 1
0.5 0.5
0 0
−0.5 −0.5
−1 −1
−1.5 −1.5
3.46 3.47 3.48 3.49 3.5 3.51 3.52 3.46 3.47 3.48 3.49 3.5 3.51 3.52
temps (secondes) temps (secondes)
Fig. 4.18: Allures temporelles des composantes synchrones et non synchrones des signaux
d’erreur (traits pleins) et de perturbation (pointillés) filtrés pour une perturbation
vibratoire réelle
30 30
erreur erreur
25 perturbation 25 perturbation
20 20
15 15
10 10
decibels
decibels
5 5
0 0
−5 −5
−10 −10
−15 −15
−20 −20
1200 1400 1600 1800 2000 2200 1200 1400 1600 1800 2000 2200
frequence reelle (Hertz) frequence reelle (Hertz)
Fig. 4.19: Densités spectrales de puissance des composantes synchrones et non synchrones
des signaux d’erreur (traits pleins) et de perturbation (pointillés) filtrés pour une
perturbation vibratoire réelle
ment que la commande optimale agit presque exclusivement sur les composantes synchrones
du signal de perturbation. En effet, les composantes non synchrones (4.18(b) et 4.19(b)) ne
sont pratiquement pas modifiées, alors que les composantes synchrones (4.18(a) et 4.19(a)),
principalement constituées par l’harmonique de la raie d’encoche, ont été fortement atténuées.
Finalement, une fois que l’algorithme récursif a fait converger le signal de commande vers sa
forme optimale, la puissance des vibrations totales du stator en θ = 0 est réduite à environ
17 % de sa valeur initiale. Ceci démontre l’efficacité d’une telle méthode de contrôle actif
pour les signaux vibratoires générés par la machine synchrone de test, dans le cas où un seul
accéléromètre est utilisé.
112
4.4. Implantation des solutions et étude des performances obtenues
À des fins de comparaison, la figure 4.20 montre les résultats obtenus en appliquant au même
signal vibratoire la méthode de compensation développée dans [Far95] et [Jar95] et décrite
au paragraphe 2.3. La raie spectrale correspondant à l’harmonique 2 de la raie d’encoche est
30
erreur
25 perturbation
20
15
10
decibels 5
−5
−10
−15
−20
1200 1400 1600 1800 2000 2200
frequence reelle (Hertz)
Fig. 4.20: Densité spectrale de puissance du signal d’erreur (trait plein) et de perturbation
(pointillés)
atténuée de 10 dB, et une composante parasite de forte amplitude est créée comme prévu
en 1666 + 4 × 23.8 = 1761, 2 Hz. La puissance du signal d’erreur représente, après conver-
gence, environ 40 % de celle de la perturbation à éliminer. On voit donc que par rapport
à la méthode précédente, les performances obtenues sont moins bonnes, même si le spectre
vibratoire est modifié sur une plage fréquentielle plus restreinte.
Il faut maintenant étudier la possibilité de minimiser les vibrations mesurées par plusieurs
capteurs à la fois, pour que le contrôle actif n’agisse plus seulement en une seule position
angulaire, mais sur une zone de la carcasse statorique.
113
CHAPITRE 4. Détermination et mise en œuvre des solutions
utilisées pour calculer la matrice T̂(0, Θ) où Θ = {θ1 , θ2 , θ3}, nécessaire à l’équation (4.22).
Les transferts LVPT situés entre les signaux de commande et les contre-vibrations qu’ils
engendrent aux positions angulaires θ1 , θ2 et θ3 peuvent donc être simulés, alors que les
vibrations traitées sont des signaux réels. Différentes courbes d’apprentissage obtenues en
fixant µ = 0.02 et ǫ = 0.001 sont tracées sur la figure 4.21. 4.21(a) représente l’évolution de la
140 100
cout puissance erreur 1
puissance commandes 90 puissance erreur 2
120
puissance erreur 3
80
100 70
60
n
80
v
v
% de P
% de P
50
60
40
40 30
20
20
10
0 0
0 20 40 60 80 100 0 20 40 60 80 100
numero de periode du systeme LVPT numero de periode du systeme LVPT
(a) C (trait plein) et Pic opt (pointillés) (b) Puissance des signaux d’erreur en sortie de
chaque capteur
Fig. 4.21: Courbes d’apprentissage pour un système de contrôle actif à trois commandes et
trois capteurs, portant sur des signaux de vibrations réels
somme des puissances de chaque commande Pic opt , ainsi que celle du coût global C calculé par
(4.20), ces grandeurs étant exprimées en pourcentage relativement à la somme des puissances
des trois perturbations. Le comportement de ce système de contrôle actif est le même que dans
le cas monocapteur, où le coût global converge vers un minimum qui devient proche de 30 %
de sa valeur initiale. La figure 4.21(b) permet de connaı̂tre dans quelle mesure les vibrations
sont atténuées pour chaque capteur puisqu’elle représente l’évolution de la puissance du signal
en sortie de chaque capteur, exprimée en pourcentage de sa valeur initiale. On remarque que
ces valeurs sont toutes réduites à approximativement 30 % de leur valeur initiale, et quel que
soit le capteur.
La figure 4.22 représente les densités spectrales de puissance de chaque signal de commande
(colonne de gauche) et de chaque signal d’erreur (colonne de droite), obtenues après la phase
de convergence de l’algorithme. On voit que ces courbes ont toutes l’allure de celles obtenues
dans le cas monocapteur, et qu’il y a bien compensation des vibrations sur chacun des trois
capteurs. On peut donc penser qu’il est possible d’atténuer les vibrations statoriques d’une
machine synchrone sur une zone continue plutôt qu’en une seule position angulaire.
114
4.4. Implantation des solutions et étude des performances obtenues
30 commande 1 30 erreur 1
perturbation 1 perturbation 1
20 20
10 10
decibels
decibels
0 0
−10 −10
−20 −20
1200 1400 1600 1800 2000 2200 1200 1400 1600 1800 2000 2200
frequence reelle (Hertz) frequence reelle (Hertz)
(a) Densité spectrale de puissance de la commande (b) Densités spectrales de puissance de la pertur-
numéro 1 bation numéro 1 (tirets) et de l’erreur numéro 1
(traits pleins)
30 commande 2 30 erreur 2
perturbation 2 perturbation 2
20 20
10 10
decibels
decibels
0 0
−10 −10
−20 −20
1200 1400 1600 1800 2000 2200 1200 1400 1600 1800 2000 2200
frequence reelle (Hertz) frequence reelle (Hertz)
(c) Densité spectrale de puissance de la commande (d) Densités spectrales de puissance de la pertur-
numéro 2 bation numéro 2 (tirets) et de l’erreur numéro 2
(traits pleins)
30 commande 3 30 erreur 3
perturbation 3 perturbation 3
20 20
10 10
decibels
decibels
0 0
−10 −10
−20 −20
1200 1400 1600 1800 2000 2200 1200 1400 1600 1800 2000 2200
frequence reelle (Hertz) frequence reelle (Hertz)
(e) Densité spectrale de puissance de la commande (f) Densités spectrales de puissance de la pertur-
numéro 3 bation numéro 3 (tirets) et de l’erreur numéro 3
(traits pleins)
Fig. 4.22: Densités spectrales de puissance des signaux de commande et d’erreur pour trois
perturbations vibratoires réelles 115
CHAPITRE 4. Détermination et mise en œuvre des solutions
L’implantation de l’algorithme (4.22) a déjà été discutée au paragraphe 4.4.1. Ceci a per-
mis de déterminer la méthode donnée par la figure 4.3, qui a été utilisée pour réaliser les
simulations présentées lors des paragraphes 4.4.2 et 4.4.3. Nous allons reprendre exactement
le même schéma, mais quelques modifications sont nécessaires afin de l’implanter en temps
réel sur un processeur numérique de signal. Les diverses TFD doivent être réalisées à l’aide
d’un algorithme de transformée de Fourier rapide (FFT = «Fast Fourier Transform»), qui
nécessite des blocs de traitement contenant un nombre d’échantillons égal à une puissance
de 2. La valeur de P , qui est celle du multiplieur fréquentiel de la figure 4.2, doit donc être
une puissance de 2, mais elle fixe aussi celle de la fréquence d’échantillonnage par la relation
νe = P νr . P doit donc être choisi de telle sorte qu’il n’y ait pas de recouvrement spec-
tral dans les signaux échantillonnés. Un exemple concret permet de déterminer une valeur
possible pour ce paramètre. Supposons que la machine synchrone sur laquelle on désire im-
planter cette méthode de contrôle actif tourne à une vitesse de rotation de νr = 25 Hz, et
que l’atténuation de ses vibrations soit nécessaire jusqu’à 3 kHz. Des filtres anti-repliement
peuvent donc être insérés dans la chaı̂ne d’acquisition afin de supprimer les composantes
fréquentielles des signaux vibratoire au-delà de 3 kHz, et νe doit alors être choisie supérieure
à 6 kHz pour vérifier la condition d’échantillonnage de Shannon. La valeur de P doit fina-
lement vérifier 25 × P > 6000, ce qui est équivalent à P > 240. Pour cet exemple précis, P
prendra la puissance de 2 immédiatement supérieure à 240, c’est à dire 256. Les différentes
TFD directes et inverses apparaissant dans la figure 4.3 porteront dans ce cas sur des blocs de
256 échantillons. On peut alors se rendre compte qu’il va être très difficile, voire impossible de
réaliser tous les calculs nécessaires à la réactualisation (4.22) en un instant d’échantillonnage.
Pour remédier à cela, on peut modifier légèrement la figure 4.3 en calculant, à partir des
mesures effectuées au bloc l, non pas les échantillons des signaux de commande du bloc l + 1
mais ceux du bloc l + 2 suivant. Ce principe, représenté à la figure 4.23, laisse au processeur
le temps d’effectuer les calculs nécessaires, mais il provoque un retard d’un bloc pour l’envoie
de la commande. Toutefois, nous avons vu que dans le cas de perturbations synchrones aux
variations du système LVPT, la commande était nécessairement périodique, de période P
correspondant à celle du système, mais aussi à la longueur d’un bloc. Ce retard sera donc
complètement transparent au niveau de la convergence de l’algorithme, ou au niveau des
performances optimales d’atténuation.
Cet algorithme a été implanté sur un processeur numérique de type TMS320C32 fabriqué
par Texas instruments, intégré sur une carte électronique fabriquée par Blue Wave Systems.
La version programmée comporte une entrée et une sortie, et calcule donc en temps-réel les
échantillons d’une seule commande, afin de minimiser le signal mesuré par un seul capteur
de vibrations. Elle a été testée sur des signaux issus de générateurs et donne des résultats
analogues à ceux du paragraphe 4.4.2. Par manque de temps, elle n’a malheureusement pas
pu être mise en œuvre sur la machine de test.
116
4.5. Conclusion
4.5 Conclusion
Ce chapitre a permis d’élaborer des solutions au problème de contrôle actif posé au cha-
pitre 1, et qui peuvent être implantées en temps-réel à l’aide d’un processeur numérique.
Dans la section 4.2, l’utilisation des notions théoriques du chapitre 3 sur les systèmes
LVPT a permis de donner l’expression (4.7), qui représente les composantes fréquentielles du
signal d’erreur à minimiser en fonction des autres caractéristiques du système.
Cette expression a été utilisée tout au long de la section 4.3 afin de développer la relation
(4.21) donnant les composantes fréquentielles optimales des signaux de commande minimisant
la fonction de coût globale (4.20). Un algorithme récursif permettant de suivre les éventuelles
variations des perturbations à minimiser a également été déterminé et est donné par l’équation
(4.22).
L’implantation pratique de ces résultats, ainsi que les performances qu’ils permettent
d’obtenir sont ensuite exposées à la section 4.4. Que ce soit avec un ou plusieurs capteurs,
les diverses simulations des paragraphes 4.4.2 et 4.4.3, et en particulier celles portant sur des
signaux réels de vibrations, montrent que l’algorithme proposé permet d’atténuer fortement
les vibrations en une zone donnée de la carcasse statorique. La plupart des résultats présentés
dans cette section ont également été exposés dans [Gra99c, Gra99a, Gra00]. Finalement, le
paragraphe 4.4.4 décrit le principe de programmation qui a été utilisé pour implanter cet
algorithme dans un processeur numérique du signal. Les résultats alors obtenus correspondent
à ceux des simulations précédentes.
Les performances des solutions présentées dans ce chapitres atteignent donc les objectifs
de compensation fixés, et peuvent de plus fonctionner en temps-réel.
117
Conclusion et perspectives
Le travail de recherche présenté dans ce mémoire a permis d’élaborer une solution pratique
au problème d’isolation vibratoire d’une machine synchrone posé au chapitre 1.
L’originalité de cette méthode est d’engendrer des contre-vibrations à l’aide des forces
magnétiques générées par des courants additionnels parcourant les enroulements statoriques
de la machine considérée. Ces contres-vibrations doivent minimiser les vibrations totales de
la carcasse statorique, mesurées par des accéléromètres fixés en plusieurs points de celle-ci.
En faisant abstraction de la chaı̂ne de commande (processeur et amplificateur de puissance),
aucun actionneur externe à la machine n’est donc nécessaire, ceux-ci étant remplacés par ses
bobines statoriques.
Cet avantage matériel est acquis au prix d’une plus grande complexité théorique puisque le
transfert entre les courants de commande et les contre-vibrations qu’ils engendrent n’est plus,
comme dans la plupart des cas pratiques, invariant dans le temps. En effet, la modélisation
effectuée à la section 1.4 du premier chapitre permet de montrer que les caractéristiques de
ce système linéaire varient périodiquement dans le temps, avec une fréquence proportionnelle
à celle de la fréquence de rotation de la machine. Il fait donc partie du groupe des systèmes
linéaires variant périodiquement dans le temps (LVPT), qui constituent une généralisation des
systèmes linéaires invariants dans le temps (LIT). De plus, l’analyse d’un signal de vibrations
mesuré sur une machine synchrone en rotation est effectuée au paragraphe 1.3.3. Elle montre
que les composantes synchrones au phénomène de rotation sont particulièrement importantes
dans le spectre vibratoire obtenu. Ce sera donc ce type de composantes que le système de
contrôle cherchera à atténuer.
Les travaux précédemment effectués sur le sujet, dont les principaux résultats sont exposés
au chapitre 2, ont abouti à l’élaboration d’un premier système de contrôle. Celui-ci peut
utiliser les techniques classiques de contrôle actif malgré la variabilité temporelle du transfert,
grâce à un ensemble de filtrages sélectifs et d’une prémodulation. Ce système est capable
d’annuler une composante fréquentielle du spectre vibratoire de la machine, mais provoque
alors l’apparition d’une raie parasite, dont l’amplitude ne peut être contrôlée par l’utilisateur.
Cet inconvénient est exclusivement dû à la présence du système LVPT, et peut conduire à
des performances de compensation inacceptables.
Il s’avère donc nécessaire de faire une étude plus approfondie des systèmes LVPT, afin
d’élaborer une stratégie de contrôle prenant en compte leurs particularités. Ceci est réalisé
au chapitre 3, où le modèle choisi est celui de la réponse impulsionnelle de ces systèmes. Une
119
CONCLUSION ET PERSPECTIVES
étude dans le domaine temporel est réalisée au paragraphe 3.2.3, et permet de mettre en
évidence le double effet de décalage fréquentiel et de filtrage LIT qu’ils provoquent sur leur
entrée. Ce double effet apparaı̂t également dans le domaine fréquentiel, où les liens existant
entre les composantes fréquentielles des signaux d’entrée et de sortie sont donnés par la
réponse bifréquentielle définie au paragraphe 3.2.4. Enfin, l’étude des signaux propres de ces
systèmes, réalisée au paragraphe 3.2.5, conduit à la définition d’un «gain complexe matriciel»
qui constitue l’extension du gain complexe classique des systèmes LIT. En effet, cette matrice
permet de relier les «représentations fréquentielles augmentées» des signaux d’entrée et de
sortie à l’aide d’un simple produit matriciel (équations (3.40) et (3.41)). Comme il est montré
au paragraphe 3.2.7, son analyse peut être conduite à l’aide d’outils classiques de l’algèbre
linéaire (décomposition en valeurs singulières, calcul de normes), et fournit des renseignements
utiles sur le système considéré, comme son inversibilité, ou les gains entrée-sortie minimaux
et maximaux qu’il peut présenter. Un exemple de système LVPT est alors analysé en utilisant
les différents outils théoriques qui viennent d’être présentés. Ce système a la même structure
que celui dont le modèle a été élaboré au chapitre 1, et sa réponse impulsionnelle, sa réponse
bifréquentielle, ainsi que son gain complexe sont calculés et interprétés à la section 3.3.
Après avoir présenté les outils théoriques permettant de caractériser les systèmes LVPT,
une solution pratique au problème initialement posé peut être déterminée. Ceci est l’objet
du chapitre 4, où il est tout d’abord montré que pour obtenir un système LVPT discret,
l’échantillonnage des différents signaux doit être synchronisé sur la fréquence de rotation de
la machine. On utilise alors le gain complexe matriciel de ce système pour pouvoir exprimer les
signaux de vibrations statoriques à minimiser, appelés «signaux d’erreur». Cette expression,
donnée par la relation (4.7), est valable pour plusieurs capteurs fixés sur le stator de la
machine et fournit les composantes fréquentielles des signaux d’erreur en fonction des autres
paramètres du système. Elle est utilisée dans l’équation (4.20) pour construire un critère
quadratique appelé «fonction de coût». Cette grandeur peut être physiquement interprétée
comme la somme de la puissance des signaux d’erreur, avec une fraction de celle des signaux
de commande. La commande minimisant cette fonction est alors considérée comme optimale
vis-à-vis de ce critère, et peut être calculée soit d’une manière directe (équation (4.21)), soit
par un algorithme récursif (équation (4.22)). Ce dernier est plus avantageux que la forme
directe car il autorise le suivi d’éventuelles variations du spectre des vibrations à atténuer.
La section 4.4 étudie l’implantation pratique matérielle et logicielle de cet algorithme (figu-
res 4.2 et 4.3) pour aboutir, au paragraphe 4.4.4, à l’implantation de sa version la plus simple
(un seul courant de commande minimisant les vibrations mesurées par un seul capteur de
vibrations) dans un processeur numérique. Il faut remarquer que cette mise en œuvre nécessite
la connaissance du gain complexe matriciel du système LVPT considéré, ce qui réclame :
– une phase d’identification effectuée machine à l’arrêt (avant le contrôle), qui n’est va-
lable que pour une seule fréquence d’échantillonnage et donc une seule vitesse de rota-
tion,
– la connaissance de la phase du champ magnétique généré dans la machine synchrone,
120
CONCLUSION ET PERSPECTIVES
Le problème qu’il reste à résoudre pour parvenir à une méthode de contrôle facilement im-
plantable sur une machine synchrone est donc principalement l’identification du gain complexe
matriciel nécessaire à l’algorithme de calcul de la commande. En effet, si cette identification
était réalisée en ligne, elle permettrait de suivre les éventuelles fluctuations de cette matrice,
mais aussi d’adapter automatiquement le système de contrôle à un changement de vitesse
de rotation, et donc de fréquence d’échantillonnage. Il semble que le domaine temporel soit
plus propice à l’élaboration d’un algorithme d’identification automatique pour les systèmes
LVPT, en utilisant les techniques de gradients adaptatifs.
De plus, il serait intéressant de pouvoir adapter ce système à d’autre types de machines tour-
nantes, en particulier les machines de type asynchrone. Le problème est que pour de telles
machines, les fréquences mécaniques et électriques ne sont pas identiques, ce qui génère des
vibrations comprenant deux familles d’harmoniques, dont l’une est liée au phénomène de
rotation, et l’autre au fondamental des courants statoriques. Une solution possible serait de
générer des signaux de commande constitués par la somme de deux signaux différents, cal-
121
CONCLUSION ET PERSPECTIVES
culés par un algorithme synchronisé pour l’un sur la rotation de la machine, et pour l’autre
sur ses courants statoriques.
Enfin, l’étude théorique des systèmes LVPT réalisée au chapitre 3 permet d’envisager leur
application dans la modélisation de signaux vibratoires de machines. En effet, nous avons vu
qu’un tel système provoque sur son entrée à la fois des décalages fréquentiels, et des filtrages
LIT. S’il est correctement synchronisé et alimenté par un bruit blanc ajouté au fondamental
de la rotation de la machine considérée, sa sortie est alors formée d’une composante aléatoire
stationnaire, d’une autre cyclostationnaire, et d’un ensemble de sinusoı̈des de fréquences mul-
tiples à la fréquence de rotation. Ces caractéristiques correspondent à celles d’un signal de
vibrations d’une machine synchrone, qui peut donc être modélisé par le système LVPT en
question. Ce modèle pourrait être utilisé à des fins de diagnostic ou de surveillance pour de
telles machines. En effet, un dysfonctionnement modifiant le signal vibratoire ferait alors de
même sur les paramètres du système LVPT, et pourrait donc être détecté.
122
Bibliographie
123
BIBLIOGRAPHIE
124
BIBLIOGRAPHIE
125
BIBLIOGRAPHIE
126
BIBLIOGRAPHIE
127
RÉSUMÉ
Ce travail de recherche est consacré à l’élaboration d’une méthode de compensation active des vibrations d’une
machine tournante électrique. Son originalité tient au fait que les enroulements statoriques sont alimentés par
des courants de commande additionnels afin d’engendrer des forces radiales sur le stator. Celui-ci répond alors
par des vibrations additionnelles qui interagissent avec les vibrations naturelles de la machine. Le but de ce
système de contrôle actif est donc de calculer la valeur ”optimale” de ces courants, permettant de minimiser
la puissance des signaux vibratoires au niveau de capteurs accéléromètriques fixés à la périphérie du stator.
Dans un premier temps, la modélisation du transfert situé entre les commandes et les contre-vibrations en-
gendrées conduit à un système linéaire et variant périodiquement dans le temps (LVPT). La fréquence fonda-
mentale de ses variations est alors proportionnelle à la fréquence de rotation de la machine.
Après avoir montré l’insuffisance des méthodes classiques de contrôle actif pour cette application, une étude
théorique détaillée des systèmes LVPT est réalisée. Elle conduit à la définition d’une matrice de transfert,
jouant le même rôle que la fonction de transfert classique employée pour les systèmes linéaires et invariants
dans le temps. Cette matrice permet d’écrire simplement la relation entre les entrées et les sorties du système
considéré dans le domaine fréquentiel.
Finalement, les résultats précédents sont utilisés afin de déterminer l’expression optimale des courants de com-
mande minimisant la puissance des signaux de vibrations mesurés. Un algorithme récursif permet également de
converger vers cet optimum, et de prendre en compte d’éventuelles variations des perturbations vibratoires à
éliminer. Divers résultats, obtenus sur des signaux synthétiques puis sur des signaux vibratoires réels, illustrent
les performances obtenues par ce système de contrôle actif. Il permet une réduction significative des vibrations
synchrones au phénomène de rotation, sans pour autant modifier les caractéristiques des autres. Enfin, son
implantation en temps-réel dans un processeur numérique de traitement de signal est discutée et réalisée.
ABSTRACT
This work is devoted to elaborate an active control system of rotating machine vibrations. It is based on
additional currents supplying the stator coils of the machine. They generate radial forces on the stator frame,
and finally create additional vibrations which interact with the machine ones. Therefore, the aim of this system
is to process the optimal value of the input currents, in order to minimize the vibration signals power measured
on the stator frame by several accelerometers.
First, the transfer function between the input currrents and the engendered vibrations is modelized by a linear
and periodically time-varying (LPTV) system. Moreover, its fondamental frequency is shown proportional to
the machine rotating frequency.
Applied to the present problem, classical active control methods reach middling performance, a theoritical
study of LPTV systems is thus realized. It leads to the definition of a transfer matrix, having the same
properties than the classical transfer function used for linear and time-invariant systems. By using this matrix,
a simple frequency relation can be obtained linking inputs and outputs of a LPTV system.
Then, previous results are used to determine the optimal value of the input currents minimizing the power
of the measured vibration signals. Moreover, an adaptive algorithm permits to reach this optimum and to
take into account some potential variations of the vibratory disturbances. The performance obtained by this
processing are illustrated thanks to simulations on synthetic and real data. The synchronous vibrations with
the machine rotation are significantly reduced, unlike the asynchronous ones, which are left unchanged. The
real-time implementation of this algorithm is finally described and realized thanks to a digital signal processor.
MOTS-CLÉS
contrôle actif ; systèmes non stationnaires ; vibrations ; machine synchrone ; traitements adaptatifs ; systèmes
variant périodiquement dans le temps ; machines tournantes ; temps-réel
BIBLIOGRAPHIE
129