République Algérienne Démocratique et Populaire

Ministère de l’Enseignement Supérieur et de la Recherche Scientifique

Université de Constantine
Faculté des Sciences de l’Ingénieur
Département d’Electrotechnique

MEMOIRE
Pour l’Obtention du diplôme de
Magister en Electrotechnique
Option: Modélisation et Commande des Machines électriques
Présenté par
BOUCHAREB Ilhem
Ingénieur d’Etat en Electrotechnique de l’Université de Constantine
Thème

Modélisation & Simulation de Défauts

D’une
Machine Synchrone à Réluctance Variable

Soutenu le, /11/2009

Devant le Jury composé de :

Président : Aissa Bouzid Professeur Université de Constantine
Rapporteur : Amar Bentounsi Maître de Conférences Université de Constantine
Examinateur 1 : Hocine Benalla Professeur Université de Constantine
Examinateur 2 : Abdesselam Lebaroud Maître de Conférences Université de Skikda
 REMERCIEMENTS
Je commencerais par exprimer ma profonde gratitude envers mon encadreur

M. Amar BENTOUNSI, enseignant au Dept. d’Electrotechnique de l’U.M.C.,

pour le sujet d’actualité qu’il m’a proposé et les précieux conseils qu’il n’a

cessé de me prodiguer.

Comme je tiens à remercier M. Abdesselam LEBAROUD, enseignant au

Dept. d’Electrotechnique de l’Université de Skikda, pour son aide durant la

mise en œuvre des méthodes de diagnostic dont il a la maitrise !

Je remercierai également tous les membres du jury de soutenance pour leur

intérêt envers ma thématique et leur contribution pour enrichir ce travail, à

savoir MM. Aissa Bouzid (Président), Hocine Benalla et A. Lebaroud

(examinateurs).

Je terminerais par rendre hommage à mes enseignants, tous cycles

confondus ; sans oublier tous ceux et celles qui m’ont aidée et encouragée pour

finaliser ce travail de recherche : MERCI !

_________________
SOMMAIRE

i
 Table des matières

Introduction générale …………………………………………………………………….. 01

Chapitre (I) Le Diagnostic : Etat de l’Art ………………………………………… 03

Introduction ……………………………………………………………………………… 04
I. Généralités sur le diagnostic …………………………………………………………… 05
[Link]érentes notions utilisées en diagnostic……………………………………………. 05
[Link] d’un système………………………………………………………………. 07
[Link]édure de diagnostic………………………………………………………………. 07
[Link] de défauts………………………………………………………………………. 08
[Link]éristiques souhaitables d’un système de diagnostic……………………………. 09
[Link] des mesures dans un système de diagnostic…………………………. 11
[Link] des méthodes de diagnostic…………………………………………….. 13
[Link]érents types de maintenance………………………………………………………. 16
[Link] appliqué aux machines…………………………………………………….. 17
[Link] analysées………………………………………………………………….. 17
II. Différents types de défauts dans les machines électriques…………………………….. 18
[Link] niveau de la machine……………………………………………………………… 18
[Link] défauts au niveau de l’onduleur…………………………………………………. 11
[Link] origines des défauts……………………………………………………………… 21
III. Analyse des défauts par la mise en œuvre de techniques de diagnostic appropriées…. 22
III.1 Analyse fréquentielle………………………………………………………………… 22
Conclusion ………………………………………………………………………………... 26

Chapitre (II) Choix et Etude Théorique d’une

Structure à Réluctance Variable…………………………... 28

Introduction……………………………………………………………………………….. 29
[Link] machine à réluctance variable (MRV)…………………………………………….. 30
[Link] machine synchrone à réluctance (MSR)………………………………………….. 44
Conclusion ……………………………………………………………………………….. 67

Chapitre (III) Analyse par éléments finis (sous FEMM et Flux2D)……………… 70

Introduction…………………………………………………………………… .…………… 71

III.1. Choix et méthodes d’étude…………………………………………………………….. 72

[Link]. Les méthodes numériques ……………………………………………………………. 75
III.3. Le principe du calcul par éléments finis …………………………………………. ….. 75

ii
[Link] logiciel FEMM…………………………………………………………………… 76
[Link] logiciel FLUX 2d…………………………………………………………………. 77
III.4. Modélisation par élément finis………………………………………………………… 78
III.4.1.Régime statique par la mise en œuvre du logiciel FEM……………………………… 78
III.4.1.1.Définition et maillage de la structure d’une MRV 6/4 …………………………….. 78
III.4.1.2.Définition et maillage de la structure d’une MSRV………………………………… 85
III.4.2. Régime dynamique par la mise en œuvre du logiciel FLUX 2d……………………...91
III.4.2.1. LA MRV 6/4………………………………………………………………………. 91
III.4.2.2. La MSRV avec un rotor Axialement Laminé……………………………………... 113
III.5. Les défauts mécaniques………………………………………………………………… 131
III.6. L’analyse spectrale……………………………………………………………………... 133
Conclusion…………………………………………………………………………………… 137

Chapitre (IV) Modélisation Commande et Simulation de Défauts

sous Matlab/Simulink………………………………140
Introduction …………………………………………………………………………… … 141
IV.1. Présentation de l'environnement MATLAB/SIMULINK™…………………… … 141
[Link]èle de simulation…………………………………………………………… ….. 142
[Link] équations électromagnétiques…………………………………………………… 142
[Link] du couple………………………………………………………………… 143
[Link]èle Linéaire de la MRV…………………………………………………………. 144
[Link] stratégie de commande par hysteresis………………………………………….. 145
IV.5.3. Simulation du régime linéaire sain…………………………………………. ……. 148
IV.5.4. Simulation de défauts……………………………………………………………... 151
IV.5.4.1. Le cas d’un court circuit monophasé………………………………………………. 151
IV.5.4.2. Le cas d’un court circuit triphasé……………………………………………………. 154
IV.5.4.3. Le cas d’un circuit ouvert………………………………………………………………. 156
IV.5.5. Commentaire……………………………………………………………………… 158
IV.6. Régime non linéaire………………………………………………………………… 159
IV.6.1. Simulation………………………………………………………………………………………… 162
IV.6.2. La technique de la commande en tension…………………………………………... 163
[Link] de l’onduleur……………………………………………………… …... 163
IV.6.2.2. Le type des interrupteurs constituent l’onduleur ………………………………… 164
IV.6.2.3. Le type de la commande appliqué ……………………………………………….. 165
IV.6.2.4. La commande en tension………………………………………………………….. 166
IV.6.4. Simulation du régime non linéaire sain…………………………………………… 169
IV.6.5. Simulation avec défauts……………………………………………………………. 172
IV.6.5.1. Le cas d’un court circuit monophasé 25% dans la phase a ……………………. 172
IV.6.5.2. Le cas d’un court circuit triphasé………………………………………………... 174
IV.6.5.3. Le cas d’un circuit ouvert…………………………………………………………. 176
IV.6.5.4. Le cas d’un IGTB d’une phase tombe en panne…………………………………. 177
IV.6.6. Simulation sain et en défauts en cas de commande en courant par hystéresis……… 180
IV.6.7. Le cas sain…………………………………………………………………………… 181
IV.6.8. Simulation de défauts………………………………………………………………… 183
IV.6.8.1. Un court circuit monophasé 25 % dans la phase a………………………………….. 184
IV.6.8.2. Le cas d’un court circuit triphasé 25%…………………………………………….. 187

iii
 IV.6.8.3. Le cas d’un circuit ouvert……………………………………………………… 189
IV.6.9. Commentaire…………………………………………………………………… 191
IV.7. Synthèse et Analyse spectrale de défauts par la FFT……………………………. 192
IV.7. 1. Regime non linéaire………………………………………………………….. 192
IV.7. 1.1. Commande en courant le cas d’un court circuit monophasé 25%…………. 192
IV.7. 1.2. Commande en tension le cas d’un court circuit monophasé 25%…………… 193
Conclusion……………………………………………………………………………….. 195

Conclusion générale………………………………………………………………………... 198

Annexe………………………………………………………………………………….201

Annexe A………………………………………………………………………………….
202
Annexe B…………………………………………………………………………………..
204
Annexe C……………………………………………………………………………………. 206
Annexe D……………………………………………………………………………………. 209
Annexe E……………………………………………………………………………………… 212

iiii
 INTRODUCTION
GÉNÉRALE

De part son impact et sa position stratégique dans tout processus industriel, la

maintenance, dans sa conception moderne, ne peut plus se limiter à être corrective ou
préventive mais se doit d’être prédictive. La surveillance des machines est devenue un art qui
permet de diagnostiquer à temps et avec précision les défauts potentiels pour agir rapidement
avant la « casse » !
Les convertisseurs électromécaniques d’énergie occupent une place de plus en plus
importante dans les équipements industriels, surtout avec les nouvelles exigences en matière
de traction électrique ou de conversion d’EnR qui requièrent des structures innovantes de
machines.
Si, de plus, ces dispositifs occupent une place critique dans le processus, ceci génèrera des
contraintes sévères en termes de sûreté de fonctionnement et de disponibilité ; d’où la
nécessité d’un diagnostic robuste d’éventuelles défaillances. Durant la dernière décade, les
recherches dans ce domaine ont fait l’objet de nombreux travaux axés essentiellement sur la
machine asynchrone.
Ainsi, l’objectif du sujet de magister qui nous a été proposé serait d’appliquer ces
méthodes éprouvées à d’autres structures de machines, types machines synchrones à
reluctance (MSR) ou à reluctance variable (MRV). L’approche retenue est une approche
signal permettant une analyse spectrale des signatures par FFT des défauts simulés
(excentricité axe rotation, circuits ouverts/court-circuits des bobinages, des circuits de
puissance des convertisseurs, …).

1
 Après étude théorique et modélisation de 2 topologies de machines (MSR et MRV)
associées à leurs convertisseurs statiques (au début, le cahier de charges ne nous imposait que
l’étude d’une seule machine type MSR !), nous avons simulé divers types de défauts :

• numériquement, par la méthode des éléments finis (sous FEMM et Flux2D) ;

• analytiquement, sous environnement Matlab/Simulink (en dynamique, linéaire et non-

linéaire).

Pour mener à terme ces objectifs, notre mémoire s’articule autour des 4 chapitres :

I. Etat de l’Art (méthodes de diagnostic des machines électriques)

II. Choix et étude théorique de structures de machines réluctantes (MSR et MRV)
III. Analyse de défauts par éléments finis (sous FEMM et Flux 2D)
IV. Simulation de défauts sous Matlab/Simulink (machine et convertisseur statique)

Les résultats obtenus et commentés sont très encourageants, comme nous le montrerons tout
au long des développements qui vont suivre.
_________________________

2
CHAPITRE : I

LE DIAGNOSTIC :

ETAT DE L’ART

3
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

INTRODUCTION

Le diagnostic de défauts des machines électriques est aussi ancien que les machines elles-
mêmes. Chaque famille de machines possédant des caractéristiques différentes, les méthodes
de diagnostic qui leur sont associées sont également différentes. De même, la puissance et
le coût des machines à surveiller jouent un rôle important dans la méthode de diagnostic
choisie.

Grâce à leurs atouts, en termes de coût et de performances, les machines à réluctance

variables trouvent actuellement de multiples applications dans la production d‘énergie
électrique (éoliennes, alternateurs, …), la traction électrique (véhicules électriques, navires,
…), … D’où notre intérêt pour ce type de machines et à leur diagnostic, lequel était jusque là
réservé aux machines asynchrones.

Ce chapitre est essentiellement consacré à l'état de l’art de la thématique de diagnostic.

Nous avons scindé ce chapitre en trois parties: (i) la première concerne les différentes notions
et définitions utilisées en diagnostic et système de diagnostic; (ii) la seconde partie analyse les
différents types de défauts qui affectent les machines électriques ; (iii) la troisième partie est
consacrée à la mise en œuvre des techniques de diagnostic appropriées, en l’occurrence
l’analyse fréquentielle.

4
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

I. Généralités sur le diagnostic [1] [2]

De nos jours, les systèmes industriels modernes deviennent de plus en plus complexes et leurs
commandes de plus en plus sophistiquées. En même temps, la fiabilité, disponibilité et sûreté
de fonctionnement sont devenues très importantes ; elles constituent de véritables enjeux pour
les entreprises actuelles. Le diagnostic des systèmes est apparu dans le but d’améliorer les
points précédents. Terme peu répandu pour les domaines techniques, il y a une trentaine
d’année, le diagnostic est devenu un sujet d’importance stratégique.

Quelle que soit la branche d’industrie concernée, aujourd’hui les procédés industriels sont
couplés à un calculateur numérique qui ne se contente pas de faire l’acquisition des données
mais qui est chargé de la mise en œuvre de l’automatisation. Automatiser peut avoir des
objectifs divers. Les plus fréquents sont d’augmenter les performances du système de
production, de garantir la qualité du produit fabriqué, ou de diminuer les coûts de fabrication.
Mais un autre objectif est de plus en plus pris en compte par les systèmes de contrôle-
commande : améliorer la sécurité de l’installation industrielle, des hommes qui y travaillent et
de son environnement. D’autres objectifs peuvent s’y rajouter tels que garantir la satisfaction
des normes de qualité. On arrive ainsi à une automatisation qui ne s’attache pas seulement au
maintien de quelques variables à des valeurs de référence, mais qui s’intéresse au procédé
dans a globalité, et dans ses divers modes de fonctionnement.

Dans ce qui va suivre, nous allons successivement développer :

• les différentes notions utilisées en diagnostic (défauts, détection, localisation,
identification, résidu, symptôme)
• les caractéristiques souhaitables d’un système de diagnostic
• les transformations des mesures dans un système de diagnostic
• une classification des approches de diagnostic (approches à base de modèle
quantitatif ; de modèle qualitatif ; approches à base d’extraction des caractéristiques
quantitatives ; caractéristiques qualitatives)

I.1. Les différentes notions utilisées en diagnostic : Supervision et diagnostic

La surveillance industrielle joue un rôle primordial dans l’augmentation de la disponibilité
des installations industrielles, réduisant ainsi les coûts directs et indirects de la maintenance
des équipements de production. Les coûts directs de cette maintenance sont ceux relatifs aux

5
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

diverses pièces de rechange, main d’œuvre, etc. par contre, les coûts indirects sont
essentiellement dus au manque à gagner engendré par un arrêt de la production. Par
conséquent, les entreprises soucieuses d’avoir une meilleure maîtrise des coûts de
maintenance ont intérêt à adopter une politique de surveillance.

La sûreté est définie comme un ensemble de propriétés indispensables à une installation

industrielle. La surveillance se réfère à la capacité de reconnaître un comportement anormal et
de le signaler. La supervision est fondée sur la surveillance et la capacité d’agir, pour éviter
des dommages sur l’installation et pour continuer la production, si cela est possible.

Une défaillance est définie comme l’inaptitude d’un composant à remplir toute ou partie de
sa fonction. Un défaut est un écart entre un comportement attendu et un comportement
observé. Constaté grâce à un indicateur de défaut, il révèle la présence d’une défaillance.
L’indicateur de défaut est en pratique élaboré grâce aux variables mesurables de
l’installation.

Un système physique, ou un procédé technologique, est conçu pour réaliser une mission
donnée. Cette mission n’est pas toujours remplie de la même façon, elle peut être totalement
ou partiellement remplie avec différents niveaux de qualité. L’évolution de ce procédé (ou
système) au cours du temps peut être décrite par différents modes de fonctionnement. On
distingue en général :

♦ Les modes normaux sous lesquels la mission est remplie, l’un d’eux est le mode nominal où
la mission est correctement remplie,
♦ Les modes évolutifs qui caractérisent des modes, souvent transitoires, de passage d’un mode
de fonctionnement à un autre, les modes de démarrage et d’arrêt, qui sont évidemment pris en
compte au moment de la conception du système de contrôle-commande.
♦ Les modes anormaux sous lesquels la mission est partiellement ou totalement non remplie,
incluant les modes défaillants, les modes interdits, qui sont les différents états que peut
prendre un procédé sous l’effet d’une panne ou d’un dysfonctionnement. Ces modes
particuliers peuvent nécessiter soit d’arrêter immédiatement l’installation, soit de basculer la
commande vers un mode de repli, ou encore vers un mode dégradé qui consiste à modifier les
consignes des boucles locales ou même carrément l’architecture de la commande ou du
procédé.

6
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

1.2. Le diagnostic d’un système est l’identification du mode de fonctionnement, à chaque

instant, par ses manifestations extérieures (symptômes). Son principe général consiste à
confronter les données relevées au cours du fonctionnement réel du système avec la
connaissance que l’on a de son fonctionnement normal ou défaillant. Si le mode de
fonctionnement identifié est un mode défaillant, le système de diagnostic devrait localiser sa
cause. Il est évident que si l’on veut représenter tous les défauts possibles d’une installation
industrielle, la tâche peut être considérable. En conséquence, la plupart des méthodes que
nous évoquerons ci-dessous se fondent sur les connaissances du comportement normal de
l’installation. Mais il est bien clair que si l’on a une bonne connaissance des anomalies
possibles, il faut l’utiliser pour améliorer la surveillance et le diagnostic.

1.3. La procédure de diagnostic

On distingue généralement, dans la procédure de diagnostic, deux étapes principales :
détection et localisation, plus l’identification pour certaines méthodes. Les définitions
précises de ces notions ainsi que les différents types de défauts et la définition des termes
"résidu" et "symptôme" (souvent utilisés en diagnostic) sont rappelés dans la suite de cette
section.

Détection : La détection (Fault detection) est la première décision que doit prendre un
système de surveillance. C’est le fait de décider qu’un indicateur de défaut est pertinent et
donc de signaler si le système est défaillant ou pas. La détection revient à choisir entre deux
hypothèses la première correspond au fonctionnement normal, la deuxième correspond au
fonctionnement défaillant. Pour détecter les défauts d’un système, il faut donc être capable de
classer les situations observables comme étant normales ou anormales.

Localisation : La localisation (Fault isolation) est le fait de déterminer le(s) sous-ensemble(s)

fonctionnel(s) défaillant(s). Cet ensemble a intérêt à être le plus réduit possible. La
localisation est plus difficile dans le cas de défauts multiples où plusieurs défauts peuvent
avoir lieu en même temps.

Identification : L’identification ou l’estimation (Fault evaluation) vise à fournir des

informations quantitatives sur les caractéristiques du défaut : instant d’apparition, valeur de
l’écart d’un paramètre permettant de le quantifier. Dans le cas de défauts procédé,
l’identification revient à connaître les valeurs de différents paramètres du système.

7
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

I.4. Types de défauts

Un défaut peut être défini comme une modification suffisamment importante et permanente
des caractéristiques physiques d’un système ou d’un composant. Un système défaillant (qui
subit un défaut) est incapable d’accomplir correctement sa ou ses fonctions avec les
performances requises. Comme le montre la Figure I.1, nous distinguons trois types de
défauts :
♦ Les défauts actionneurs
♦ Les défauts procédés
♦ Les défauts capteurs

Il est clair ainsi que pour s’attaquer à ce problème de diagnostic, il faut des connaissances
approfondies sur l’installation : connaissances de son comportement normal ou anormal.
Souvent un défaut est modélisé avec les mêmes outils que ceux utilisés pour représenter le
procédé en état normal.

Figure I.1 – Types de défauts

Résidu : Les résidus (Residual) sont définis comme étant des indicateurs de défauts. Ils
expriment l’écart par rapport au fonctionnement normal dans un espace de représentation
donné. Cet espace dépend de la nature des informations manipulées (variable numérique,
booléenne, variable qualitative, forme). C’est sur la base de leurs valeurs que la décision de
détection est prise ou pas. Souvent, les résidus sont considérés nuls (statistiquement) en
l’absence de défaut.

8
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

Symptôme : Un symptôme (Symptom) est un phénomène (caractère) qui révèle un trouble

fonctionnel. Les alarmes sont des symptômes de comportement anormal souvent utilisés pour
facilité la surveillance et la supervision.

I.5. Caractéristiques souhaitables d’un système de diagnostic

Afin de pouvoir comparer diverses approches de diagnostic, il est utile d'identifier un
ensemble de caractéristiques souhaitables qu'un système de diagnostic devrait posséder. Bien
que ces caractéristiques ne soient habituellement possédées par aucune méthode de diagnostic
particulière, elles sont utiles pour évaluer diverses méthodes en termes d'information a priori
qui doit être fournie, de rapidité, de robustesse, d’adaptabilité, etc. Chaque fois qu'une
anomalie se produit dans un procédé, un système de diagnostic fournit un ensemble
d'hypothèses de défauts qui explique les observations. L’efficacité d'un système de diagnostic
demande que le défaut réel soit un sous-ensemble de l'ensemble des défauts suspectés.
La résolution d'un système de diagnostic demande que l’ensemble des défauts fourni soit
aussi minimal que possible. Les paragraphes suivants présentent un ensemble de
caractéristiques souhaitables qu’un système de diagnostic doit posséder.

I.5.1 Rapidité et sensibilité de détection

Un système de diagnostic devrait répondre rapidement en détectant et en diagnostiquant les
mauvais fonctionnements du procédé. Cependant, les réponses rapides du diagnostic et les
performances tolérables durant un fonctionnement normal sont deux buts différents. Un
système qui est conçu pour détecter un défaut rapidement sera sensible aux influences de
toutes les erreurs (bruit, incertitudes…). Ceci rend le système sensible au bruit et peut
déclencher fréquemment des fausses alarmes pendant un fonctionnement normal.

I.5.2 Robustesse
On voudrait que le système de diagnostic soit robuste aux divers bruits et incertitudes. On
voudrait que ses performances se dégradent lentement au lieu d'échouer totalement et
abruptement. La robustesse est incompatible avec les systèmes très sensibles aux seuils. En
présence du bruit, ces seuils doivent être choisis raisonnablement. Ainsi, les besoins de
robustesse doivent être équilibrés avec ceux de sensibilité.

9
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

I.5.3 Identification de nouveaux défauts

Une des fonctionnalités minimales qu'un système de diagnostic doit posséder, est de pouvoir
décider, étant donné des conditions courantes de fonctionnement, si le procédé fonctionne
normalement ou anormalement, et s’il est anormal, si la cause est un défaut de
fonctionnement connu ou un inconnu. Ce critère est connu comme identification de
nouveauté (nouveaux événements). En général, on dispose de données suffisantes pour
décrire le comportement normal du procédé. Cependant, typiquement on n'a pas de tels
historiques de données disponibles pour modéliser les modes de fonctionnement anormaux
d'une manière satisfaisante. Seuls quelques gabarits de données peuvent être disponibles pour
couvrir une partie du fonctionnement anormal. Ainsi, il est possible qu'une grande partie de la
région de dysfonctionnement ne puisse avoir été modélisée. Ceci posera des défis sérieux
pour l'identification de nouvelles situations. Même dans ces conditions difficiles, on voudrait
que le système de diagnostic puisse identifier l'occurrence des nouveaux défauts et ne les
classifie pas comme d’autres défauts de fonctionnement connus ou comme fonctionnement
normal.

I.5.4 La localisation

la localisation et la capacité du système de diagnostic à distinguer entre différents défauts.

Dans des conditions idéales exemptes de bruit et d'incertitudes de modélisation, ceci revient à
dire que le système de diagnostic devrait pouvoir générer un ensemble d’hypothèses qui ne
contient pas les défauts absents. Naturellement la capacité de concevoir de tels systèmes
dépend en grande partie de la connaissance de la structure et des caractéristiques du procédé.

I.5.5 Adaptabilité
En général les procédés évoluent en raison des changements des entrées externes ou des
changements structurels. Les conditions de fonctionnement du procédé peuvent changer non
seulement en raison des perturbations mais également en raison des changements des
conditions environnementales telles que les changements des quantités de production avec
des demandes variantes, des changements de qualité de matière première, etc. Ainsi le
système de diagnostic devrait être adaptable aux changements. Il devrait être possible de
développer graduellement la portée du système pendant que les nouvelles situations émergent
grâce à la disponibilité d'informations supplémentaires.

10
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

I.5.6 Identification des défauts multiples

La capacité d'identifier les défauts multiples est un objectif important, mais difficile, à cause
de l’influence simultanée de la plupart des défauts. Dans un système non linéaire général, des
interactions se manifestent, et par conséquent un système de diagnostic ne peut pas employer
les différents modèles de chaque défaut pour modéliser l'effet combiné des défauts. Par
ailleurs, l'énumération des différentes combinaisons de défauts multiples et l’analyse de leurs
effets devient combinatoire pour des procédés de grande échelle.

I.5.7 Estimation d’erreur de classification

Une condition pratique importante pour un système de diagnostic, est l’établissement des,
niveaux de confiance. Ceci pourrait être considérablement facilité si le système de diagnostic
pouvait fournir une évaluation a priori sur l'erreur de classification qui peut se produire.
De telles mesures d'erreur seraient utiles pour projeter des niveaux de confiance sur les
décisions produites par le système donnant à l'utilisateur une meilleure sensation à propos de
la fiabilité des recommandations proposées par le système de diagnostic.

I.5.8 Facilité d’explication

En dehors de la capacité d'identifier la source d’un défaut de fonctionnement, un système de
diagnostic devrait également fournir des explications sur la façon dont le défaut a commencé
et a évolué vers la situation actuelle. C'est un facteur très important lorsque l’on conçoit des
systèmes interactifs d'aide à la décision en ligne. Ceci exige la capacité de raisonner au sujet
de la cause et des effets dans un procédé. Un système de diagnostic doit justifier ses
recommandations de sorte que l'opérateur puisse en conséquence évaluer et agir en utilisant
son expérience. On voudrait que le système de diagnostic justifie non seulement pourquoi
certaines hypothèses ont été proposées, mais expliquer également pourquoi certaines autres
hypothèses n'ont pas été proposées.

[Link] des mesures dans un système de diagnostic

Généralement , le processus de diagnostic peut être défini comme un processus décisionnel.
La décision finale est constituée de l’ensemble des défauts possibles. La Figure 1.2 , montre
les diverses transformations que les données du procédé subissent pendant le processus de
diagnostic.

11
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

Figure I.2- Transformations de mesures dans un système de diagnostic

L'espace des mesures est défini par les sorties du système, ces sorties peuvent avoir subi des
prétraitements (élimination de données aberrantes, filtrage,…). Ce sont les entrées du système
de diagnostic.

L'espace des caractéristiques est l’espace obtenu en fonction des mesures en utilisant une
connaissance a priori du problème. Ici, les mesures sont analysées et combinées à l'aide d’une
connaissance a priori du procédé pour extraire les caractéristiques utiles concernant le
comportement du procédé, pour faciliter le diagnostic. Il y a deux manières de développer
l'espace des caractéristiques à partir de l'espace des mesures, à savoir, la sélection et
l'extraction des caractéristiques. L'extraction de caractéristiques est un procédé qui facilite,
par l'utilisation de la connaissance antérieure du problème, une transformation de l'espace de
mesure dans un espace de dimension inférieure.

Le passage de l'espace des caractéristiques à l'espace des décisions est habituellement fait
pour satisfaire une certaine fonction objective (telle que réduire au minimum la non
détection). Cette transformation est réalisée en employant une fonction discriminante, dans
certains cas en employant des fonctions simples de seuil.

L'espace des décisions est l’espace correspondant au nombre de variables de décision

obtenu par des transformations appropriées de l'espace des caractéristiques.

L'espace des classes correspond au nombre de classes de défaut. En classant les défauts, nous
indiquons catégoriquement à quelle classe un défaut appartient. L'espace des classes est ainsi
l'interprétation finale du système de diagnostic fourni à l'utilisateur.

12
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

I.7. Classification des méthodes de diagnostic

Compte tenu de l’importance des enjeux en termes de productivité (arrêt inutile des
installations), de sécurité (anomalie non détectée) ou de qualité de production (mesure
incorrecte d’une grandeur à contrôler), de nombreuses approches ont été utilisées pour
apporter une contribution à la solution du problème de diagnostic. On distingue cependant,
parmi les différentes méthodes, deux familles principales ; celles qui utilisent un modèle du
système à surveiller, et utilisent généralement les techniques de l’automatique ; et celles pour
qui seules les données acquises sur le procédé considéré, qui se limitent à son observation
passée et présente, permettent de caractériser son mode de fonctionnement. Parmi les
méthodes basées sur l’utilisation d’un modèle du procédé, deux sous-familles importantes
peuvent également être dégagées ; celles utilisant un modèle qualitatif et celles recourant à
un modèle quantitatif (dit aussi modèle numérique). Pour certaines méthodes, cette
distinction est cependant sans fondement, car elles utilisent de manière mixte les deux types
de modèles en essayant de tirer profit des avantages et s’affranchir des limitations de chacun
d’entre eux. Toutes les méthodes nécessitent une phase d’apprentissage qui permet de
connaître une référence de fonctionnement normal ou dégradé. Cette référence permet de
générer en temps réel les symptômes lors de la phase de diagnostic. Les sections suivantes
présentent les grandes familles de méthodes de diagnostic représentées par la Figure 1.3.

I.3- Les méthodes de diagnostic

13
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

I.7.1 Approche à base de modèle qualitatif

C’est une connaissance s’exprime avec des symboles et peut être traitée, en simulant la façon
dont un humain raisonne logiquement, par des techniques informatiques relevant de
l’intelligence artificielle. Pour construire des représentations qui permettent une meilleure
compréhension des phénomènes physiques mis en jeux dans un procédé, l’intelligence
artificielle s’intéresse au milieu des années 1980 au raisonnement qualitatif. Le
développement de techniques de représentation qualitative du monde pour en comprendre
schématiquement le fonctionnement est une idée très naturelle et par conséquent très
ancienne. On trouve dans la proche qualitatif (Modèles causaux, Abstractions hiérarchiques).

I.7.2 Approche à base de modèle quantitatif

L’utilisation de modèle mathématique pour le diagnostic est très largement répandue. Les
méthodes de surveillance à base de modèle (Model-based- FDI : Fault Detection and
Isolation) utilisent un modèle et reposent sur la génération et l’étude d’un signal particulier
appelé "indicateur de défaut" ou "résidu".

Un résidu est un signal qui reflète la cohérence des données mesurées vis-à-vis d’un modèle
comportemental du système. Les modèles utilisés peuvent être de nature et de complexité
différentes. Ils peuvent être : à temps continu ou à temps discret, linéaires ou non linéaires,
représentant le bon fonctionnement ou tenant compte des défaillances. Les sorties du modèle
se représentent classiquement sous la forme de valeurs numériques.

La localisation qui suit la détection peut être effectuée à l’aide d’une table de signature ou
matrice d’incidence. Les colonnes de cette table sont représentatives des différents défauts et
les lignes des différents résidus. Un 1 représente le symptôme qui indique que le résidu est
sensible au défaut, un 0 le contraire. La signature d’un défaut (colonne de la matrice) définit
l’état des symptômes lorsque ce défaut affecte le système (l’état est une grandeur logique ou
symbolique).

En surveillance, par contre, il est parfois nécessaire de compléter le modèle afin de

caractériser le comportement défaillant du système. Suivant la classe de modèle considérée,
les techniques de génération de résidus seront différentes. Suivant le niveau de connaissance
sur les défaillances, les résultats de la procédure de surveillance seront plus ou moins précis et
complets.

14
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

I.7.3 Approche à base d’extraction de caractéristiques qualitatives

Ces méthodes ont l’avantage de ne pas nécessiter la connaissance d’un modèle mathématique
ou structurel du procédé, contrairement aux méthodes à base de modèles. Seulement, la
disponibilité de grandes quantités de données historiques du procédé est nécessaire. Il y a
diverses façons de traiter ces données et de les présenter en tant que connaissance a priori
pour le diagnostic. Ceci est connu comme l’extraction de caractéristiques. Parmi les méthodes
d’extraction de caractéristiques qualitatives, se trouvent les méthodes à base de règles
(systèmes experts), et les méthodes de représentation et analyse qualitatives des tendances.

I.7.4 Approche à base d’extraction de caractéristiques quantitatives

Dans la section précédente, on a vu la méthode basée sur l’extraction de caractéristiques
qualitatives. Dans cette section, on va évoquer la méthode qui extrait les caractéristiques
quantitatives.
On peut analyser les signaux issus de l’instrumentation et proposer pour chacun de ces
signaux un modèle.

La méthode statistique
Comme : tests d’hypothèses, classification probabilistes, Analyse en composante principale &
Moindres Carrés Partiels, cette méthodes est en fait une référence du comportement d’un
signal, qui peut être représentative de l’état normal de l’installation ou d’un défaut particulier.

Les méthodes non statistiques

Telle que la classification floue, Réseaux de neurones, Méthodes de signature, Méthodes
fréquentielles, temps-fréquence et temps-échelle.
La méthode non statistique est une autre méthode de représentation pour un signal qui peut
très bien être utilisée comme référence pour le diagnostic, est sa représentation par un spectre
de fréquences, une signature graphique, ou encore l’extraction d’attributs (ou
caractéristiques), à partir des données suivie de leur classification, qui permettent, pour des
procédés dont on veut analyser plusieurs signaux liés entre eux, de discriminer divers modes
de fonctionnement, dont le mode normal et certains modes de défaillance.

15
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

Figure I.4- Les différentes méthodes de diagnostic

I.8 Les différents types de maintenance

Les différents types de maintenance utilisés dans les applications industrielles dépendent
principalement du niveau de criticité du système considéré. Plus la criticité est importante,
plus les moyens et donc le coût associé seront élevés. On distingue :
•La maintenance préventive. Celle-ci consiste à changer les systèmes à date fixe
prédéterminée. Le coût associé est important.
•La maintenance corrective : L’évaluation du système est réalisée fréquemment grâce à
une instrumentation supplémentaire. Notons que plus cette évaluation est fréquente plus
les coûts d‘immobilisation du système seront réduits. C‘est dans le cadre de la maintenance
corrective que l‘utilisation de méthodes de diagnostic devient fondamentale. En effet, plus
le diagnostic sera précis, plus l‘intervention sera courte.
•La maintenance prédictive : L'intervention est effectuée avant l'apparition des défauts,
actuellement, pour des raisons financières, les industriels s’orientent de plus en plus vers la
maintenance prédictive.

16
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

Nous allons maintenant aborder le diagnostic de défauts des machines électriques. Les
méthodes que nous allons présenter sont basées sur l’approche heuristique ce modèle
repose sur une reconnaissance de signatures déjà observées.

I.9 Le diagnostic appliqué aux machines électriques

Pour la détection d'un dysfonctionnement de machine électrique, il est nécessaire de
posséder :
♦ une bonne connaissance du comportement de la machine, de son état et de son
fonctionnement.
♦ une bonne connaissance des défaillances et de leurs conséquences sur le fonctionnement de
la machine.

I.10 Les grandeurs analysées par les méthodes de diagnostic :

Les méthodes de diagnostic heuristiques permettent maintenant au monde industriel de
développer des systèmes de surveillance performants et efficaces. Elles peuvent impliquer
plusieurs domaines scientifiques et techniques. Elles reposent sur l‘analyse de grandeurs de
type :
• électrique: analyse des grandeurs électriques aux bornes de la machine (tensions,
courants, …);
• mécanique: analyse de la vitesse, du couple ou même des vibrations (accéléromètres,
…);
• thermique: analyse de la température de la machine (sonde de température, mesure
infra rouge);
• électromagnétique: analyse des flux dans l‘entrefer (spires de mesure de flux), des flux
axiaux, des flux de fuite (bobines externes) ou des champs électriques externes (mesure de
décharges partielles);
• chimique: analyse des huiles et des gaz.

Dans notre mémoire nous avons optés pour la détection des défauts sur l’approche "signal"
des grandeurs mesurables. L’analyse fréquentielle des signaux nous permettra de mettre en
lumière les nouvelles fréquences liées à l’apparition des défauts. Trois grandeurs seront
analysées à savoir le couple moteur, le courant statorique et le flux magnétique. L’analyse de
ces trois signaux sera faite sur deux types de machine : la machine à réluctance variable et la
machine synchro_réluctante.

17
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

II. Différents types de défauts dans les machines électriques

Les défauts qui apparaissent dans la machine électrique ont des causes variées. Ils peuvent
êtres dus au vieillissement des constituants du moteur ou bien aux conditions d'utilisation ou,
tout simplement à un défaut de fabrication dont l'effet serait imperceptible lors de la mise en
service.

II.1 Au niveau de la machine

Donc dans cette partie, nous allons exposer les principaux types de défauts couramment
rencontrés dans les machines électriques conventionnelles, leurs origines.

♦ Les défauts de bobinage statorique [2] [3]

Ceux-ci peuvent être soit une ouverture, soit un court circuit d'une ou de plusieurs
phases du bobinage statorique. Dans le cas d‘une spire court-circuitée, la machine peut
continuer à fonctionner, mais si le courant de court-circuit est important il peut entraîner une
surchauffe de l'enroulement statorique et ainsi amplifier le défaut en produisant des courts-
circuits entre deux phases ou entre une phase et le neutre de la machine. Il est donc
très important de détecter un tel défaut très rapidement avant la destruction complète de la
machine. Comme montre la figure il dans les enroulements statoriques, il peut y avoir cinq
cas de défauts à savoir : court-circuit enter spires, court-circuit à la tête de la bobine,
ouverture de la phase, court-circuit phase à phase, court-circuit phase- terre.

Figure I.5- Les courts-circuits dans les enroulements statoriques

18
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

♦ Les défauts de bobinage rotorique [2] [4]

Pour les machines synchrones, il s‘agit essentiellement de courts-circuits au bobinage du
rotor. La conséquence d‘un court-circuit au rotor est un échauffement local qui peut
déformer le rotor et provoquer des vibrations. Les vibrations du rotor peuvent influencer
le fonctionnement de l'alternateur, en l'empêchant d'atteindre le facteur de puissance et la
puissance active désirés. Pour les machines asynchrones, on peut également rencontrer
des courts circuits. Pourtant, la principale source de défaillance semble être des fissures
ou des ruptures dans les barres ou les anneaux. Ces ruptures peuvent provoquer des
vibrations de la machine. Les barres cassées peuvent alors toucher le stator et endommager
ainsi l'isolation statorique. Les étincelles qui se produisent lors de la dégradation de la
cage rotorique peuvent également être critiques pour l’équipement.

♦Les excentricités
Les excentricités peuvent être statiques, dynamiques ou mixtes. Il s’agit d’un mauvais
positionnement du rotor par rapport au stator qui peut avoir pour origine un mésalignement
du à un pliement de l‘axe ou un défaut de roulement. La conséquence de ce type de défauts
est le déséquilibre de la force électromagnétique appliquée au rotor ce qui augmente
l'excentricité, la conséquence finale pouvant être le frottement entre le stator et le rotor.

Figure I.6- Différents types d’excentricités, noté × le centre géométrique du rotor

♣Centre du rotation du rotor

19
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

En résume dans la figure 1.7 les différents défauts pouvant intervenir dans une machine
électrique. cette étude statistique établit par IEEE [2].

Figure I.7- Répartition des défauts dans une machine électrique suivant leur localisation

II.2 Les défauts au niveau de l’onduleur [5]

On peut distincte les défauts selon leurs origine en deux parties : au niveau de la machine
(statorique comme les court circuit ou rotorique les excentricités) et au niveau de l’onduleur ;
Les interrupteurs de l’onduleur doivent supporter de grandes variations de température (cycles
thermiques). Pour les transistors bipolaires type IGBTs, la disparité entre le coefficient de
dilatation thermique et les coefficients du fil en aluminium et de la puce de silicone est la
cause d’efforts de cisaillement sur les interfaces de liaison qui peuvent mener à la rupture du
fil. Afin de prévoir la durée vie d'un IGBT, le nombre de cycles thermiques qu’il supporte lors
du fonctionnement doit être correctement déterminé. Indépendamment des commutateurs
d'inverseur, le condensateur C est une pièce importante de l'onduleur et devrait être pris en
considération en estimant la fiabilité de l'onduleur. Afin d'améliorer la fiabilité de l’onduleur,
les IGBTs et les diodes correspondantes ainsi que les commandes de porte (et les circuits
additionnels de commande) peuvent être fusionnés dans un module unique de puissance. De
tels modules de puissance sont souvent désignés sous le nom de modules intelligents de
puissance (IPMs). Les circuits additionnels de commande peuvent fournir des fonctions de
protection contre la surchauffe, le court-circuit, ….

20
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

II.2.1 Les défauts au niveau de l’onduleur

a. Un Onduleur Arrêté b. Une Phase court circuitée

c. Un Court circuit équilibré d. Une Phase ouverte

II.3 Les origines des défauts [2]

sont diverses ; nous citerons à titre indicatif des origines type:

♦ mécaniques: mauvaise fabrication, vibrations de la machine, forces électromagnétiques
déséquilibrées, force centrifuge, fluctuations de la charge;
♦électriques: dégradations d'isolement, décharges partielles, étincelles;
♦thermiques: pertes dans le cuivre, manque de refroidissement général ou localisé;
♦environnementales: humidité, ozone, poussière.

Ces défauts pouvant induire :

♦un déséquilibre des courants de lignes et des tensions;
♦ Une augmentation des oscillations du couple;
♦ Une diminution du couple moyen;
♦ Une augmentation des pertes et donc une réduction de l'efficacité énergétique;
♦Un échauffement excessif et donc un vieillissement accéléré.

21
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

III. Analyse des défauts par la mise en œuvre de techniques de diagnostic appropriées[6]
III.1 Analyse fréquentielle
Le diagnostic et la surveillance des systèmes industriels nécessitent parfois l'interprétation des
signaux mesurés tant dans le domaine temporel que fréquentiel. Dans ce second domaine,
c'est l'aspect périodique des signaux qui est privilégié. Un des outils mathématiques adapté à
l'analyse des signaux périodiques est la Transformée de Fourier.

a. Transformée de Fourier
Une alternative à la représentation temporelle du signal consiste à représenter son contenu
fréquentiel. Celui-ci est obtenu en calculant sa Transformée de Fourier (FT :Fourier
Transform), définie pour un signal x(t) d’énergie finie par :

L (ℜ ) → L (ℜ )
2 2

x(t ) → FT (v) = ∫ x(t ) e − i 2πvtdt = x, ei 2πvvt

(I.1)
ℜ

Cette transformée étant inversible on a :

(I.2)
x(t ) = ∫ FT x (v) e − i 2πvtdv
ℜ

La FT mesure donc un indice de similarité entre le signal et une collection d’ondes de

(exponentielles complexes) de longueur infinie, FTx(v) représentant la contribution de l’onde
fréquence ν au signal x( t) . La FT est en fait une transformée en ondes du signal dans laquelle
toute information temporelle disparaît.

b. Transformée de Fourier à temps discret

Par définition, la transformée de Fourier à temps discret (DTFT : Discrete Time Fourier
Transform) d’une suite { x [k] } est la fonction de ν , périodique de période 1 :

+∞
X (v ) = ∑ x[ K ] e − i 2πkv
k = −∞
(I.3)

C’est une fonction à fréquence continue, c’est-à-dire ν. Il est d’usage de la représenter sur les
intervalles [−1/ 2,1/ 2] ou [0,1], du fait de sa périodicité.
On déduit la formule inverse :

22
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

1/ 2
X [k ] = ∫ X (v) ei 2πkv dv
−1 / 2
(I.4)

2 2
La fonction ⏐X (v) ⏐est désignée par le terme de spectre. ⏐X (v) ⏐ s’interprète comme la
distribution de l’énergie le long de l’axe des fréquences.

+∞ 1/ 2

∑ x[k ] = ∫ X (v) dv
2 2
(I.5)
k = −∞ −1 / 2

c. Transformée de Fourier discrète

Le calcul sur ordinateur de la DTFT à partir de la valeur des échantillons de x [k] impose une
charge de calcul infinie puisque la série comporte un nombre infini de termes et que la
fréquence ν varie continûment sur l’intervalle [0,1]. C’est la raison pour laquelle la notion de
transformée de Fourier discrète (DFT : Discrete Fourier Transform) a été introduite. Son
calcul est limité à un nombre fini de valeurs de k et à un nombre fini de valeurs de ν . Son
intérêt pratique est très largement dû à la découverte d’une méthode de calcul rapide connue
sous le nom de transformée de Fourier rapide (FFT : Fast Fourier Transform).
Considérons la suite finie {x [k]}ou k=0 :N-1, le plus souvent cette suite provient de
l’échantillonnage d’un signal à temps continu. Partant de la définition, la DTFT s’écrit:

N =1 (I.6)
X (v) = ∑ x[k ] e − i 2πkv
k =0

où ν ∈[0,1] est la fréquence normalisée. En se limitant à un nombre fini de L de valeurs de la

fréquence, à savoir ν = n/L avec n =[ 0 : L −1] , on aboutit à la notion de transformée de
Fourier discrète. L et N ne jouent pas le même rôle.
L est le nombre de points de calcul de la DTFT et N est le nombre de points de la suite
temporelle. Comme nous le verrons par la suite, L influe sur la précision du tracé du spectre,
tandis que N est lié à la résolution. Toutefois prendre L = N ne restreint en rien la généralité
de la notion introduite. Cela aboutit à la définition de la transformée de Fourier discrète pour
la suite {X [n]} :

23
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

N −1
X [n] = ∑ x[k ] Φ Nnk , n = 0 : N − 1 (I.7)
k =0

Avec

Φ N = e − i 2πnk / N
nk (I.8)

La formule inverse est :

N −1 N −1
1 i 2πnk / N = 1
x[k ] =
N
∑ X [ n ] e
N
∑ X [n] Φ−Nnk , n = 0 : N − 1 (I.9)
n=0 n=0

d. Transformée de Fourier rapide

p
La transformée de Fourier discrète (DFT) d’un jeu de N = 2 échantillons.

x[k], k = 0,..., N −1 , est un ensemble de N coefficients :

N −1
x[k ] = ∑ x[k ] ei 2πnk / N , n = 0 : N − 1 (I.10)
n=0

La DFT est donc (à un facteur de normalisation près) une transformation orthogonale associée
à une matrice N×N .Le principe est le suivant. Notons :

x pair [k ] = x[2k ], ximpair [k ] = x[2k + 1], k = 0 : N / 2 − 1

(I.11)

24
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

On a :

X [n] = X pair [n] + e − 2iπn / N X impair [n], (I.12)

Le célèbre algorithme de transformée de Fourier rapide (FFT : Fast Fourier Transform) de

calcul de DFT est fondé sur l’utilisation récursive de ces formules.

X [n + N / 2] = X pair [n] + e − 2iπn / N X impair [n]. (I.13)

Lorsque N n’est pas une puissance de deux, des calculs analogues en utilisant une
décomposition en facteurs premiers conduisent au même ordre de grandeur O( N log2N ) du
nombre de multiplications et additions. Beaucoup de travaux ont contribué à l’optimisation
de cet algorithme, parmi lesquels on notera le programme FFTW (Fastest Fourier Transform
in the West) disponible sur le Net.

25
CHAPITRE I LE DIAGNOSTIC : ETAT DE L’ART

Conclusion

Dans ce chapitre, nous avons présenté un état de l’art sur la thématique de diagnostic. vu
que chaque famille de machines possède des caractéristiques différentes, les méthodes de
diagnostic qui leur sont associées sont également différentes. De même, la puissance et
le coût des machines à surveiller jouent un rôle important dans la méthode de diagnostic
choisie. Quatre catégories de méthodes ont été présentées : les méthodes à base de modèles
quantitatifs, les méthodes à base de modèles qualitatifs, les méthodes à base d’historique de
données qualitatif, et les méthodes à base d’historique de données quantitatif. Le choix d’une
de ces méthodes dépend essentiellement des connaissances disponibles sur le procédé, sans
oublier les considérations techniques et économiques.
Nous avons présenté, également, les différents défauts qui apparaissent dans les machines
électriques tournantes. Nous avons fais un survol sur les différentes méthodes de diagnostic et
nous avons opté sur l’approche signal via l’analyse spectrale. La variation du spectre est donc
un moyen correct de détecter et même parfois de localiser un défaut. Les défauts donnent
naissance à des raies spectrales bien particulières.

Dans le second chapitre nous allons abordé la structure de la machine à réluctance variable
orientée vers la modélisation et la simulation des défauts.

26
 Références du Chapitre I

Références du Chapitre I

[1] S. TALEB, «Génération des symptômes flous par fft et swt application à la détection de
vibrations», Thèse Doctorat, L’institut national polytechnique de Grenoble, 14 Février 2007.

[2] BUI. Viet Phuong, «Diagnostic des machines électriques par l’analyse du champ
magnétique application a l’identification de défauts rotoriques d’un alternateur à vide», Thèse
Doctorat, L’INP Grenoble, 4 Octobre 2007.

[3] A. LEBAROUD, «Modélisation et techniques orientées vers le diagnostic de la machine

asynchrone associée à des sources variables», Thèse Doctorat, Univ de Constantine, 23
Septembre 2006.

[4] M. BLODT, «Condition Monitoring of Mechanical Faults in Variable Speed Induction

Motor Drives», Thèse Doctorat, L’INP Grenoble ,14 Septembre 2006.

[5] O. WALLMARK, «Control of Permanent-Magnet Synchronous Machinesin Automotive

Applications», Thèse Doctorat, Goteborg , Sweden 2006.

[6] G. Heit, «Cours traitement du signale_ transformées discrète», Cours, Estaca, 17 Octobre
2007.

[7] [Link], «Analyse spectrale grandeur spectrales», Cours, L’INP Grenoble ,2008-2009.

27
CHAPITRE : II

CHOIX ET ETUDE THÉORIQUE

D’UNE STRUCTURE

A RÉLUCTANCE VARIABLE

28
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

INTRODUCTION

Les travaux sur la modélisation et le diagnostic des machines électriques ont souvent porté sur
les machines «conventionnelles» surtout asynchrones. C’est pourquoi nous allons nous
intéresser ici à des machines moins classiques type machine à réluctance variable (MRV) et
machine synchrone à réluctance variable (MSRV) qui couvrent actuellement beaucoup
d’applications ceci grâce à la simplicité de leurs structure ainsi qu’à leurs performances. Ces
atouts leur ont permis de faire un bon chemin dans le domaine des entraînements à vitesse
variable.

Dans ce chapitre nous allons présenter une étude théorique et comparative des deux types de
machines à savoir la MRV et la MSRV ; ce qui permettra leur modélisation ainsi que leur
implantation sous environnements FEMM ,flux2d,et Matlab/Simulink.

Nous entamons ce chapitre par un historique sur les machines à réluctance variable. Nous
étalerons sur les différentes structures et types, actuellement, utilisées pour les machines
synchro-réluctantes. Nous passerons en revue sur le principe de fonctionnement et les
équations classiques du modèle de chaque machine. Nous arborons ensuite sur les avantages
et les inconvénients de chaque machines.

29
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

II.1. LA MACHINE A RELUCTANCE VARIABLE

II.1.1. Historique de la MRV [1]

Le principe de réluctance variable a été également choisi pour la génération de courants
hautes fréquences (de quelques kHz à 20 kHz) dans l'alimentation de fours à induction et dans
l'émission des ondes radio. Les premiers alternateurs "à fer tournant" apparurent dans les
années 1888 ; une structure remarquable, de type homopolaire, fut conçue et commercialisée à
cette époque par la société OERLIKON , sa puissance était de 70 chevaux. Au XX siècle,
d'autres alternateurs à réluctance variable, homopolaires, ou hétéropolaires (structures de
Lorentz- Schmitt et Guy) furent réalisés et utilisés.
Il semble que les premières structures à double saillance, telles que nous les connaissons
aujourd'hui, soient apparues dans les années 1920. Un article de 1927 sur les applications de
l'électricité dans les bateaux de guerre décrit un moteur de type 6/4 pareil à celui que nous
étudions actuellement. Il était utilisé comme transmetteur de mouvement entre deux points
éloignés.

II.1.2. Les différents types de MRV [2]

Généralement, les MRV possèdent au moins trois phases afin de résoudre les problèmes de
démarrage et de réduire les fluctuations du couple mais au prix d’un convertisseur plus
volumineux !

a) Machine à stator lisse « sans balais »

La machine à réluctance commutée la plus simple possède un stator lisse et un rotor à 2p
saillances; ces saillances peuvent être dues à une anisotropie du matériau (rotor laminé).
Les enroulements statoriques constituent q phases ; ils sont non imbriqués, ont une ouverture
approximativement égale à π/q et sont alimentés successivement pendant 1/q période
(alimentation dite à 120° en triphasé). La forme d’onde des courants doit être adaptée à celle
de la perméance pour assurer un couple optimal constant.

30
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

[Link].1- Machine à stator lisse et rotor saillant commutée à 120°

[Link].2- Machine à stator lisse et rotor laminé commutée à 120°

b) Machine polyphasée à galettes monophasées (multistack)

En associant sur un même arbre q machines identiques monophasées à plots, on réalise une
machine polyphasée. Le stator est en général commun et les dentures rotoriques doivent être
décalées d’un angle électrique (2π/q) l’une par rapport à l’autre ; la configuration inverse, à
rotor commun et stators décalés, est aussi envisageable. Le stator et les rotors possèdent le
même nombre de dents. La figure (2.3) représente le schéma de principe d’une telle machine
multigalette ou multistack triphasée. La nécessité de découpler les galettes rend ces machines
assez encombrantes et elles ne sont réellement intéressantes que si elles sont munies d’une
excitation homopolaire par des aimants logés dans l’arbre rotorique.

31
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

[Link].3- Machine polyphasée à galettes

c) Machines polyphasées à plots

Le stator porte P plots ; PP le nombre de paires de plots de chacune des q phases et coïncide
avec le nombre de paires de pôles dans le cas simple et courant ou les plots ont des polarité
alternées Nord et Sud. Le nombre P est donc égal à :

P=[Link].q (II.1)

Ces plots peuvent être dentés, comme dans les machines monophasées. Quelle que soit la
configuration, on constate, en alimentant successivement chaque phase par un courant
constant, que le rotor peut prendre un nombre de positions différentes égal à :

S = [Link] (II.2)

S est souvent appelé nombre de pas par tour. Cette dénomination est d’ailleurs ambiguë, le
tour désignant la séquence d’alimentation et non le déplacement du rotor, celui-ci pouvant
avoir à effectuer plusieurs tours pour prendre toutes les positions possibles.

(a) machine à 6 et 4 dents (b) machine à plots dentés

[Link].4 - Machines triphasées bipolaires à plots

32
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

[Link].5 - Machines dissymétriques polyphasées

[Link].6 - Structures de machines à réluctance variable à double saillance

II.1.3 Choix de la structure [3] [4]

Le chois d’une machine dépend essentiellement de l’application et le cahier de charge si en
prend le problème des véhicules électrique on trouve que la conception de la motorisation
d'un V.E. doit elle tenir compte de nombreuses contraintes techniques, économiques et
opérationnelles du marché de l'automobile telles que : rendement élevé, couple élevé à vitesse
quasi-nulle, possibilité de récupération d'énergie au freinage, puissance spécifique élevée,
réduction du nombre de semi-conducteurs de puissance,……
Nous avons proposé d'expérimenter pour cette application un moteur à réluctance variable à
double saillance (MRVDS) répondant au cahier des charges et parmi les différentes structures
de MRV existantes, nous avons choisi d'utiliser une machine à réluctance variable à double
saillance (MRVDS) à grosses dents (Switched Reluctance Machines pour les Anglo-Saxons).

33
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

II.1.3. Le principe de fonctionnement d’une MRV [5] [6]

La machine à réluctance variable qui fait l’objet de notre étude, possède une double saillance
(MRVDS) avec Nr dents au rotor et un stator bobiné à Ns dents. Cette machine à la
particularité de ne pas comporter de circuit d’excitation (MRV pure). Les deux pièces polaires
sont réalisées par un empilage de tôles en matériau ferromagnétique.

[Link].7- Machine à réluctance variable à double saillance MRVDS 6/4.

Le principe de fonctionnement d'une machine à réluctance variable peut être décrit à partir
d'une structure monophasée élémentaire identique à celle représentée ci-dessous (Fig. II.8) :

[Link].8- Structure élémentaire d’une MRV monophasée.

Cette structure possède deux positions extrêmes remarquables :

34
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

♦ une position d'opposition pour laquelle le circuit magnétique présente une réluctance
maximale ou inductance minimale ;
♦ une position de conjonction où la réluctance est minimale donc d’inductance maximale.

[Link].9 - Positions extrêmes de la MRV élémentaire

Si le rotor se trouve dans une position intermédiaire entre l'opposition et la conjonction,

l'injection d'un courant dans l'enroulement d'excitation modifie l'état du système qui tend alors
à présenter une réluctance minimale (règle de flux maximal) : on observe ainsi une rotation du
rotor vers la position de conjonction. Lorsque le courant est supprimé, si l'énergie cinétique
emmagasinée durant la rotation est suffisante pour l'amener jusqu'à l'opposition, alors on peut
répéter le cycle et obtenir un mouvement continu de rotation.
Ce principe est généralisable à la machine triphasée de notre application ([Link].10)

[Link].10 - Séquences d’alimentation des phases pour produire un mouvement de

rotation

35
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

L'alimentation de la phase 1 aligne AA' avec 11'.

L'alimentation de la phase 2 aligne BB' avec 22'.
L'alimentation de la phase 3 aligne AA' avec 3'3.

II.1.4. Principes de base des MRV [7]

Cette partie rappelle les principes de base de la conversion électromécanique d’énergie
appliqués aux machines à réluctance variable. Considérons le circuit magnétique simple de la
figure 2.11 constitué d'une culasse ferromagnétique fixe et d’une pièce ferromagnétique
mobile dont l'axe est repéré par l’angle θ. Le passage d'un courant I dans les n spires de
l'enroulement crée une force magnétomotrice (fmm) définie par: F = N.I.

[Link].11- Caractéristiques Flux-At d’une MRV monophasée non excitée

Le flux traversant les N spires de l'enroulement varie en fonction de la réluctance magnétique

des différentes portions du circuit, à savoir la réluctance ℜ qui est la somme des réluctances
de la culasse stator, de la pièce mobile rotor et de l'entrefer qui est généralement dominante.
La réluctance peut être représentée par un autre paramètre dit perméance définie par :

1
℘= (II.3)
ℜ

et en relation avec le flux totalisé :

(II.4)
Ψ = L ⋅ Ι = N 2 ⋅ Ι ⋅℘

L : étant l'inductance propre de l'enroulement.

Pour un courant constant I, la variation de perméance entre les 2 positions extrêmes
(Figure.2.11). se traduira par une variation de flux donc d’énergie magnétique : pour θ = 0
(modulo π), le flux est maximal ; il est minimal pour θ = ± π/2 (modulo π).

36
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

La perméance possède donc une périodicité de 2θ ; en se limitant au premier harmonique, elle

s’exprime par :
℘ = ℘0 + ℘1cos(2θ ) (II.5)

De manière similaire, la relation entre l'inductance de l'enroulement et l'angle θ est :

L = Ν 2℘0 + Ν 2℘1cos(2θ ) = L 0 + L1COS (2θ ) (II.6)

Pour des variations du courant I entre 0 et If et du flux ψ entre 0 et ψf , nous pouvons

déterminer l'énergie magnétique Wm ainsi que la coénergie magnétique Wc (figure II.12).

[Link].12 - Energie et coénergie magnétiques.

L’énergie magnétique est définie par:

Wm = ∫ i.∂ψ (II.7)

Et la coénergie par :
(II.8)
Wc = ∫ψ .∂i

A partir de la caractéristique flux-At, on peut calculer le couple développé par une machine
donnée en suivant la méthode des travaux virtuels qui considère soit une variation de l’énergie
magnétique soit de la coénergie ; nous pouvons ainsi exprimer le couple instantané par :

ce =
∂Wm
∂θm
= − Νr
∂Wm
( )
∂θe
(II.9)
ψ = cte
.

37
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

En régime linéaire, le flux est lié aux ampères tours (ni) par la perméance ℘(θ)

Donc :
ψ (θ ) = ℘(θ ).ni (II.10)

Wc = ∫℘(θ ).[Link] = ℘(θ ).ni

1
(II.11)
2

Ce qui permet d'exprimer le couple par

Nr 2 d℘
ce = .n i . (II.12)
2 dθ

∂Wc ⎛ ∂Wc ⎞
ce = = Nr. ⎜ ⎟ (II.13)
∂θm ⎝ ∂θe ⎠ni =cte

Ou bien par rapport à l'inductance de la phase L(θ)

Nr 2 d (n2℘) Nr 2 dL
ce = .i . = .i . (II.14)
2 dθ 2 dθ

De la même manière, si on impose le flux on obtient à partir de l'expression de l'énergie

magnétique le couple :

Nr 2 2 dℜ (II.15)
ce = − . n .ψ
2 dθ
Où ℜ(θ) est la réluctance.
Lorsque la machine est saturée l'inductance de la phase ne varie plus seulement en fonction de
la position mais aussi du courant. Il faut alors calculer le couple à partir de l'énergie convertie
lors d'un cycle d'alimentation.

38
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

La puissance électromagnétique convertie Pem s'exprime alors par :

Ρem = f .W (II.16)

[Link].13 - Energie convertie lors d'un cycle d'alimentation.

La fréquence électrique d'alimentation f est liée à la pulsation de rotation Ω par :

Ω
f = Nr. (II.17)
2.π
Ainsi, le couple électromagnétique Cem d’un moteur à q phases s’exprime par :

[Link] W
Cem = = q. N r . (II.18)
Ω 2.π

Equations électriques :
Lorsque l’on applique une tension V aux bornes d’une phase de la machine, on peut décrire
son fonctionnement par les équations suivantes :

dψ (θ , i ) dφ (θ , i ) (II.19)
V = r.i + = ri + n.
dt dt
⎡ di ∂φ (θ , i ) dθ ∂φ (θ , i ) ⎤
V = r.i + n.⎢ . + .
⎣ dt ∂i dt ∂θ ⎥⎦ (II.20)

dθ
ω= = Ω.Nr
(II.21)
dt

39
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

En régime linéaire (Ψ = L(θ). I ), on peut réécrire la dernière équation sous la forme :

⎡ di dL ⎤
V = r.i + n.⎢ L. + i.ω . ⎥ (II.22)
⎣ dt dθ ⎦

Le terme dû à la variation d’inductance en fonction de la position en i.ω.(dL /dθ ) est

classiquement dénommé force contre-électromotrice.

1.1.5- Le mode d’alimentation

Théoriquement, il existe 2 modes fondamentaux d'alimentation en moteur :

a) Alimentation en courant unidirectionnel

Type haché synchronisé ([Link].14) avec le mouvement du rotor (capteur position) ; on
obtient l’équivalent d’une machine synchrone.
Ce courant en créneaux est équivalent à la superposition d’un courants continu DC (excitation
série) plus un courant alternatif AC dont la pulsation ω est reliée à la pulsation de rotation du
rotor Ω par :

(II.23)
ω = Nr.Ω .

[Link].14 - Alimentation en courant unidirectionnel

40
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

La figure II.15 montre la forme d’onde du courant, l’inductance, le flux et le couple en

fonction de l’angle électrique :

10
Courant [A]

0
0 50 100 150 200 250 300 350 400
Linductance [H]

0.5
0 50 100 150 200 250 300 350 400
10
Flux [Wb]

0
0 1 2 3 4 5 6 7 8
20
Couple [Nm]

-20
0 50 100 150 200 250 300 350 400
teta(°)

[Link].15- Courant, inductance, flux et couple pour une alimentation avec un courant
unidirectionnel..

b) Alimentation en courant alternatif

En remarquant que le sens du courant n'intervient pas dans le fonctionnement précédent (sens
de rotation et signe du couple ne dépendent que des valeurs de θ pour lesquelles le bobinage
est alimenté ou non), on devrait obtenir un fonctionnement identique au précédent avec une
alimentation purement alternative ([Link].16) de pulsation:

ω’= (Nr/2).Ω. (II.24)

41
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

[Link].16 - Alimentation en courant bidirectionnel

Comme le cas précédant la figure II.17 montre la forme d’onde du courant, l’inductance, le
flux et le couple en fonction de l’angle électrique pour une alimentation en courant alternatif.

5
Courant [A]

-5
0 50 100 150 200 250 300 350 400
Linductance [H]

0.5
0 50 100 150 200 250 300 350 400
5
Flux [Wb]

-5
-4 -3 -2 -1 0 1 2 3 4
5
Couple [Nm]

-5
0 50 100 150 200 250 300 350 400
teta(°)
[Link].17 - Courant, inductance, flux et couple pour une alimentation avec un courant
bidirectionnel.

42
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

Photos du prototype de la machine à réluctance variable (fig. II.18): [8]

Fig. II.18

Avantages et inconvénients d’une MRV [5] [9]

Avantages Inconvénients

* Fabrication simple * Mauvais facteur de puissance

* Machine robuste * Electronique de commande assez compliquée
* Possibilité de fonctionner sur une large (nécessité d’un capteur de position)

plage de vitesse * Présence non négligeable de bruit et de

* Couple massique élevé vibrations (couple pulsatoire)

* Coût attrayant, multiples applications, .. * Coût du convertisseur.

43
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

II.2. LA MACHINE SYNCHRONE A RELUCTANCE VARIABLE

II.2.1. Constitution [10] [11]

La machine synchro-réluctante est structurellement une machine synchrone à pôles saillants
dépourvue d’excitation. Son stator est identique à celui des machines à courant alternatif
ordinaires. Le couple électromagnétique est constitué exclusivement du couple de saillance.
La conversion d’énergie dans la machine s’effectue par variation des inductances propres et
mutuelles de ses enroulements de phases due à la rotation d’un rotor magnétiquement
dissymétrique entre les axes direct (de réluctance minimale) et en quadrature (de réluctance
maximale). Il s’agit bien d’une machine à réluctance variable et plus précisément d’une
machine à réluctance pure polyphasée à stator lisse avec une alimentation sinusoïdale.
Comme la machine tourne au synchronisme avec son alimentation et comme elle fonctionne
par variation de réluctance, nous avons choisi l’appellation « machine synchro-réluctante »
(Synchronous Reluctance Motor).

Le rotor d’une machine synchrone à réluctance variable ne contient ni aimants, ni bobinage

d’excitation ; le couple est créé grâce à l’effet de réluctance. Il est construit de manière à ce
que le rapport entre l’inductance dans l’axe direct et l’axe en quadrature (Ld /Lq ) soit le plus
important possible. La plage de fonctionnement à puissance constante est directement liée à
ce rapport. Il en est de même pour le facteur de puissance : plus ce rapport est élevé, plus le
facteur de puissance l’est. L’obtention d’un rapport Ld/Lq élevé induit des contraintes au
niveau de la fabrication qui se répercutent négativement sur le coût.

[Link].17- Structure d’une MSRV

44
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

II.2.2. Caractéristiques
La machine synchro-réluctante convient donc aux applications à forte puissance et vitesse
élevée, domaine couvert actuellement par la machine asynchrone. Mais cette dernière présente
des pertes Joule et des pertes fer au rotor en régime permanent. La machine synchro-
réluctante est donc sur ce point une concurrente de la machine asynchrone. De plus, le rotor,
pouvant être massif, est robuste et joue le rôle d’amortisseur. En outre, du fait de l’absence
d’excitation, son courant de court-circuit est largement plus faible que celui d’une machine
synchrone. Cela diminue le coût du système de protection et la taille du convertisseur
d’alimentation.
La machine synchro-réluctante couvre une large gamme de puissance de 750 [W] à 100 [kW]
et de vitesse: de 3000 [tours/min] à 48000 [tours/min].

II.2.3. Historique [1]

Dans les années 1930, les "moteurs synchrones à pôles saillants non excités" (synchrones à
réluctance variable) ont commencé à être étudiés. On avait besoin, pour certaines applications,
d'une vitesse précise et constante mais aussi d'un démarrage autonome. Ainsi, le moteur
synchrone (à champ tournant) à pôles saillants mais dépourvu d'excitation, à condition qu'il
soit muni d'une cage d'écureuil de démarrage, pouvait satisfaire à ces exigences. Ces moteurs
furent d'un emploi assez restreint car leur facteur de puissance et leur rendement restaient
faibles, la structure électromagnétique (stator à pôles lisses et présence de la cage) ne
permettaient pas un rapport de saillance suffisant (grande variation d'inductance). Dans les
années 1960, en Angleterre, on s'intéressa de nouveau à ces moteurs ; ce fut probablement le
point de départ des travaux des équipes de plusieurs universités anglaises, notamment celles
de Leeds LAWRENSON (Professeur Peter J. ...) et Nottingham. Ensuite, sont apparus les
moteurs synchrones à réluctance variable à rotor segmenté, à barrières de flux et à rotor
axialement laminé : ces techniques avaient pour but d'augmenter le rapport de saillance et, par
là, le facteur de puissance et les performances. À JARRET la même époque, naissait en
France, avec les frères , un nouvel engouement pour les moteurs à réluctance variable Vernier
à grand nombre de dents et fort couple massique : ces moteurs semblaient particulièrement
adaptés à la réalisation d'entraînements directs à basse vitesse, par exemple, pour des roues de
véhicules électriques... Quant au moteur à réluctance variable à double saillance auto-
commuté, le terme équivalent anglo-saxon : "Switched Reluctance Motor" semble être apparu
en 1969 ; c'est aujourd'hui le terme employé dans la littérature scientifique internationale pour
qualifier ces machines. Plusieurs laboratoires universitaires français, dont le LÉSiR, ont

45
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

travaillé récemment dans le domaine des machines à réluctance variable à double saillance à
alimentation électronique. En 1994, plusieurs industriels en commercialisent, ce sont
Allenwest Ltd. au Royaume Uni, Sicme-Motori en Italie. Les applications envisagées
aujourd'hui sont les entraînements industriels à vitesse variable, et celles potentielles sont
certains accessoires automobiles et électroménagers, les démarreurs alternateurs de turbines
(aéronautiques)...

II.2.4. Les différents types de machines synchrones à réluctance variable [10] [12]
Dans ce paragraphe, nous donnons une description des différentes structures de rotor qui ont
été développées jusqu’à présent pour ce type de machine. L’optimisation de la structure du
rotor pour augmenter le rapport de saillance (Ld/Lq) et améliorer les performances de la
machine présente un des principaux axes de recherche dans ce domaine.
On distingue trois principales structures de rotors de machines synchro-réluctantes : massive,
avec barrières de flux et axialement laminée.

a) Machine synchro-réluctante massive :

Le rotor massif est constitué d’un bloc de matériau ferromagnétique. En général, il s’agit d’un
cylindre sur lequel ont été faites des découpes pour créer la saillance du rotor. Cette structure
est la plus simple à fabriquer est donc aussi la plus économique. Elle est également la plus
robuste car constituée d’un seul bloc non assemblé. Elle est ainsi pourvue de la meilleure
tenue mécanique et au moins de ce point de vue elle est idéale pour un fonctionnement à haute
vitesse. Par ailleurs, la section de passage du flux d’axe d est assez large pour que le rotor soit
peu saturé. Pour réduire les pertes aérodynamiques, le rotor peut être caréné par du matériau
non magnétique ou peut être à trous comme le montre la figure 1.18. Mais son principal
défaut est la faiblesse relative du rapport de saillance (Ld/Lq).

46
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

[Link].18 - Machines synchro-réluctantes massives.

b) Machine synchro-réluctante avec rotor de barrières de flux :

Les barrières de flux sont un des dispositifs qui permettent d’augmenter le rapport de saillance
(Ld/Lq) (jusqu’à environ 13). Le rapport de saillance est augmenté surtout par la réduction de
l’inductance Lq , c’est à dire par l’augmentation de la réluctance du chemin du flux d’axe q
dans le rotor. Le rotor devient alors un assemblage de segments ferromagnétiques et non
magnétiques. On règle le rapport de saillance en jouant sur les largeurs relatives des segments.
Les barrières de flux peuvent être constituées d’air (trous dans le rotor). Un dimensionnement
soigneux est alors nécessaire pour assurer la solidité du rotor et la réduction du flux d’axe q.

[Link].19 - Machines synchro-réluctantes avec barrières de flux rotoriques.

47
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

Pour préserver l’équilibre mécanique du rotor, les barrières de flux peuvent aussi être en acier
non magnétique. Un technique d’assemblage robuste (soudage par explosion) est alors
indispensable pour permettre au rotor de fonctionner en haute vitesse.

c) Machine synchro-réluctante avec rotor axialement laminé :

Un deuxième dispositif permettant d’augmenter le rapport de saillance est le laminage axial
du rotor. Le rotor est alors constitué d’une succession de feuilles ferromagnétiques et non
magnétiques. Le rotor se comporte alors comme un matériau homogène anisotrope. C’est
cette anisotropie du matériau du rotor qui assure la dissymétrie entre l’axe direct et l’axe en
quadrature. On règle le rapport de saillance en jouant sur les épaisseurs relatives des feuilles.
On peut atteindre ainsi un rapport de saillance avoisinant 20 en deux pôles et 10 en quatre
pôles. Notons que cette structure est généralement utilisée pour les machines synchro-
réluctantes à quatre pôles. Les feuilles assemblées sont alors en forme d’hyperboles comme
le montre la figure. Du fait du nombre élevé d’éléments assemblés et de la relativement faible
tenue mécanique qui en découle, cette structure est limitée aux faibles vitesses (inférieures à
5000 tr/min) et faibles puissances (inférieures à 2kW) mais elle possède des performances
plus élevées que les deux premières.

[Link].20 - Machines synchro-réluctantes avec rotor axialement laminé.

48
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

¾ Assistance par aimants

Dans le but de réduire davantage le flux d’axe q, on peut utiliser des aimants permanents
logés dans les barrières de flux. Contrairement à une excitation de type ordinaire, le flux créé
par les aimants s’oppose au flux en quadrature. La majeure partie du couple reste due à la
saillance du rotor. Mais du fait de la polarité des aimants, la machine perd sa réversibilité en
sens de rotation. Bien entendu, l’utilisation des aimants augmente le coût de la machine.

[Link].21-a- Assistance par aimants d’une machine synchro-réluctante. Configurations à

deux et à quatre pôles.

¾ Utilisation de supraconducteurs comme isolants magnétiques

Le supraconducteur à haute température critique est de loin le meilleur dispositif pour réduire
le flux d’axe q. Il est utilisé en tant qu’isolant magnétique. Il peut être utilisé pour caréner un
rotor massif ou à la place des barrières de flux comme montré sur la figure. L’introduction des
barrières de flux en matériau supraconducteur présente l’avantage de réduire
considérablement Lq en préservant une valeur élevée de Ld. L’isolation magnétique entraîne
en conséquence une augmentation considérable du rapport de saillance et par conséquent une
amélioration importante du facteur de puissance. Contrairement à l’assistance par aimants,
l’utilisation de supraconducteurs permet de conserver la réversibilité de la machine par
rapport au sens de rotation. Mais elle a deux inconvénients majeurs : le coût du matériau
supraconducteur et la complexité de la structure et de la mise en œuvre.

49
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

[Link].21. b- Utilisation de supraconducteurs comme écran contre le flux d’axe q.

II.2.5. Choix de la structure étudier [11]

Nous avons choisi, pour notre étude diagnostic et simulation de défauts, le rotor axialement
laminé car le rapport de saillance (Ld/Lq) qui est un paramètre déterminant des performances
intrinsèques de la machine atteint des ratios supérieurs à 8, ce qui fait que ce type de rotor est
utilisé dans une « pile-volant » de parc éolien.

50
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

[Link].22 - Le rotor axialement laminé du modèle de MSRV choisie.

II.2.6. Principe de fonctionnement d’une machine synchrone à réluctance variable

[10] [12]

Le stator de la machine synchro-réluctante est bobiné de la même façon que celui des
machines synchrones ordinaires. Il s’agit d’un bobinage triphasé avec p paires de pôles
alimenté par un système triphasé équilibré de courants de pulsation ω. Le bobinage triphasé
crée alors une force magnétomotrice (f.m.m) tournante avec une vitesse angulaire de ω/P.
Son rotor saillant présente une « dissymétrie » entre l’axe direct et l’axe en quadrature. Il se
positionne par rapport à la f.m.m tournante de manière à ce que la réluctance traversée par le
flux d’induction magnétique dans l’entrefer soit la plus petite que lui permet la charge qu’il
entraîne. En tournant, la force magnétomotrice entraîne ainsi le rotor à la même vitesse ω/P.
L’angle θ0 entre le maximum de la fmm et l’axe d du rotor est appelé « angle de charge ».

51
 CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

[Link].23 - Machine diphasée équivalente.

En régime linéaire, l’énergie magnétique de la machine s’exprime en fonction des

inductances, de la valeur efficace I2 du courant diphasé et de l’angle de charge par :

W = I 22 (L d cos2 θ 0 + L d sin 2 θ 0 ) =
1 2
2
( ) 1 2
(
I 2 L d − L q cos 2θ 0 + I 2 L d + L q
2
) (II.25)

Nous pouvons constater que l’énergie possède une partie non convertible indépendante de θ0,
dite énergie magnétisante, et une partie convertible en cos(2θ0) dont l’amplitude est
proportionnelle à la différence (Ld-Lq) . Nous obtenons le couple en dérivant l’énergie par

Rapport à l’angle de charge :

ce =
∂W
⎛θ ⎞
( )
= −Ρ I 22 L d − L q sin 2θ 0 (II.26)
∂⎜⎜ 0 ⎟⎟
⎝ Ρ ⎠
Soit :

( ) (
ce = p L d − L q i d i q = 2 p L d − L q I d I d) (II.27)

52
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

Où Id et Iq sont les valeurs efficaces des courants direct et en quadrature id et iq, liés par :

id (II.28)
Id =
2

iq
Iq = (II.29)
2

Ce couple est maximum pour -θ0 = 45 , c’est à dire quand Id=Iq=I2/√2 . Dans ce cas, le
couple s’écrit :

⎛ Ld ⎞
⎜ ⎟
( ) 2L ⎜ Lq ⎟ (II.30)
ce = p L d − L q I 2 = Ρ
2 I 2 d⎜ Ld ⎟
⎟
⎜
⎜ L q ⎟⎠
⎝

Le point de fonctionnement à couple maximum est particulièrement intéressant. En ce point,

le facteur de puissance peut approximativement s’exprimer par :

Ld
−1
1 Lq
cos ϕ ce max ≈
2 2 (II.31
⎛ Ld ⎞
⎜L ⎟ +1
⎝ q⎠

Nous avons représenté sur la figure II.24 , les variations du facteur de puissance en fonction
du rapport de saillance Ld/Lq .
Nous pouvons observer sur cette figure II. que le facteur de puissance commence à devenir
intéressant pour des rapports de saillance supérieur à 6. La prise en compte des pertes
électriques (pertes Joule et pertes dans le fer) augmente légèrement ce facteur de puissance.

53
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

facteur de puissance en fonction du rapport de saillance

0.9

0.8

0.7

0.6
cos(phimax)

0.5

0.4

0.3

0.2

0.1

0
0 2 4 6 8 10 12 14
le rapport(Ld/Lq)

Figure II.24 - Facteur de puissance en fonction du rapport de saillance Ld /Lq

La production de couple vient donc directement de la dissymétrie entre l’axe direct et l’axe en
quadrature, autrement dit de la saillance du rotor ainsi le facteur de puissance dépend du
rapport Ld/Lq ,donc en essayant de préserver une valeur suffisante de Ld .

II.2.7. Modèle mathématique d’une machine synchrone à réluctance [13] [12] [15]
a) Hypothèses
Pour pouvoir développer le modèle électrique équivalent de la machine, il est nécessaire de
faire certaines hypothèses. Ces hypothèses classiques sont les suivantes :
♦On considère une distribution sinusoïdale des forces magnétomotrices d’entrefer créées par
les bobinages statoriques et rotoriques.
♦On néglige les phénomènes de saturation et d’hystérésis magnétique.
♦On néglige les pertes dans le fer de la machine.
♦On néglige l’effet de peau (en particulier dans les barres du rotor).
♦On néglige l’effet de la température sur la valeur des résistances.

54
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

On donne sur la figure II.25, une représentation symbolique de la machine bipolaire

équivalente. Le stator est composé de trois bobinages déphasés de 2p/3 dans l’espace. La cage
du rotor peut être modélisée par deux enroulements en quadrature, l’un placé suivant l’axe d
(axe de faible entrefer) et l’autre suivant l’axe q. Cette représentation simplifiée du rotor ne
permet pas d’accéder à la connaissance du courant circulant effectivement dans chaque barre
mais conduit à une traduction assez fidèle de l’influence des barres rotoriques sur le
comportement de la machine. L’angle θ représente la position mécanique du rotor et p le
nombre de paires de pôles.

Fig. II.25 - Représentation symbolique de la machine

b) Equations dans un repère lié au stator

Nous allons rappeler les équations principales du modèle de la machine non saturée dans le
repère du stator.
Les expressions des flux à travers les bobinages statoriques et rotoriques sont les suivantes :

55
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

⎡ψ ⎤
⎢ sa ⎥ ⎡ L sa (θe) L sasb (θe) L sasc (θe) L sard (θe) L sarq (θe )⎤ ⎡i sa ⎤
⎢ ⎥
⎢ψ sb ⎥ ⎢ L sbsa (θe ) Lsb(θe ) L sbsc (θe ) L sbrd (θe ) L sbrq (θe )⎥⎥ ⎢i sb ⎥
⎢ ⎥ ⎢
⎢ sc ⎥ = ⎢ L scsa (θe )
ψ L scsb (θe ) L sc (θe ) L scrd (θe ) L scrq (θe )⎥.⎢i sc ⎥
⎥ ⎢ L rdsa (θe ) ⎥⎢ ⎥
⎢
ψ L rdsb (θe ) L rdsc (θe ) 0 ⎥ ⎢i rd ⎥
⎢ rd ⎥ ⎢
L rd
⎢ ⎥
⎢ψ ⎥ ⎢⎣ L rqsa (θe ) L rqsb (θe ) L rqsc (θe ) 0 L rq ⎥⎦ ⎢i rq ⎥
⎣ rq ⎦ ⎣ ⎦

… (II.32)

Les inductances qui apparaissent dans l’équation (II.7) dépendent pour la plupart de la
position du rotor. avec l’hypothèse du premier harmonique d’espace, les expressions des
inductances propres et mutuelles du stator s’écrivent:

Lsa (θe ) = L 0 + L 2 cos(2θe ) ⎛ ⎛ 2π ⎞⎞

L sasb (θe ) = M 0 + M 2 cos⎜⎜ 2⎜θe + ⎟ ⎟⎟
⎝ ⎝ 3 ⎠⎠
⎛ ⎛ 2π ⎞⎞
L sb (θe ) = L 0 + L 2 cos⎜⎜ 2⎜θe − ⎟ ⎟⎟ L sbsc (θe ) = M 0 + M 2 cos(2θe )
⎝ ⎝ 3 ⎠⎠
⎛ ⎛ 2π ⎞⎞
⎛ ⎛ 2π ⎞ ⎞ L scsa (θe ) = M 0 + M 2 cos⎜⎜ 2⎜θe − ⎟ ⎟⎟
L sc (θe ) = L 0 + L 2 cos⎜⎜ 2⎜θe + ⎟ ⎟⎟ ⎝ ⎝ 3 ⎠⎠
⎝ ⎝ 3 ⎠⎠
… (II.33)

Avec les valeurs des coefficients de la décomposition en série de Fourier des inductances
données dans le tableau II.1 [12] ,on représente les inductance et mutuelles statorique sur la
figure II.26.

L0 L2 M0 M2
0.270 0.050 -0.095 0.125

56
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

0.32

0.31
inductances statoriques Lsabc(teta)en H

0.3

0.29

0.28

0.27

0.26

0.25

0.24

0.23

0.22
0 20 40 60 80 100 120 140 160 180
position du rotor en degre

0.1

0.05
mutuelles statoriques Lsabc(teta)en H

-0.05

-0.1

-0.15

-0.2

-0.25
0 20 40 60 80 100 120 140 160 180
position du rotor en degre

Figure II.26 - Les inductances propres et mutuelles du stator.

57
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

Les expressions des inductances mutuelles entre les enroulements statoriques et les
enroulements rotoriques sont données par :

L sard (θe ) = M srd cos(2θe ) L sarq (θe ) = M srq sin (θe )

⎛ 2π ⎞ ⎛ 2π ⎞
L sbrd (θe ) = M srd cos⎜θe − ⎟ L sbrq (θe ) = M srq sin ⎜θe − ⎟
⎝ 3 ⎠ ⎝ 3 ⎠
⎛ 2π ⎞ ⎛ 2π ⎞
L scrd (θe ) = M srd cos⎜θe + ⎟ L scrq (θe ) = M srq sin ⎜θe + ⎟
⎝ 3 ⎠ ⎝ 3 ⎠

… (II.34)

Il existe une relation entre les coefficients L 0 et M0 et entre les coefficients L2 et M2 ainsi
apparaissant dans l’expression des inductances propres et mutuelles du stator :

M0 1 M2
=− et =1
L0 2 L2 (II.35)

Les équations générales des tensions s’obtiennent en écrivant la loi de Faraday pour chacun
des enroulements en considérant la chute de tension ohmique. Comme les enroulements
rotoriques équivalents à la cage sont en court-circuit, la tension appliquée est nulle :

d ψ sa
V sa = R s . i sa +
dt d ψ rd
d ψ sb 0 = R rd . i rd +
V sb = R s . i sb + dt
dt (II.36)
d ψ rq
d ψ sc 0 = R rq . i rq +
V sc = R s . i sc + dt
dt

L’expression du couple électromagnétique est obtenue par dérivation de la coénergie :

1 t ∂[L]
cem =
2
[i] ⋅ ∂θe ⋅ [i ] (II.37)

Ce modèle, valable dans le cadre des hypothèses définies précédemment, peut s’écrire plus
simplement dans un repère commun lié au rotor.

58
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

c) Modèle de la machine dans le repère (d-q) lié au rotor

Pour simplifier les équations du modèle, on se place dans un repère lié au rotor (axes d-q).
Les grandeurs statoriques sont ramenées dans le repère du rotor en appliquant la
transformation de Park. Cette transformation mathématique est rappelée ci-dessous pour les
courants. Elle est parfaitement définie quel que soit θe. Elle est orthogonale et possède une
matrice inverse.

La matrice de Park normée est :

⎡ 2π 2π ⎤
⎢cosθ cos(θ −
3
) cos(θ − )⎥
3
⎢ 2π 2π ⎥
[Ρ] = 2 ⎢ sin θ sin(θ − ) sin(θ − ) ⎥
3⎢ 3 3 ⎥ (II.38)
⎢ 1 1 1 ⎥
⎢⎣ 2 2 2 ⎥⎦

Avec :

[ Ρ]− 1 = [ Ρ]t (II.39)

Toutes les grandeurs électriques triphasées sont projetées sur le repère lié au rotor par la
transformation de Park, nous écrivons alors dans le cas général:

⎡ x a⎤ ⎡xd ⎤ (II.40)
⎢ ⎥ t⎢ ⎥
⎢ x b ⎥ = [ Ρ] ⎢ x q ⎥
⎢⎣ x c ⎥⎦ ⎢⎣ x h ⎥⎦

⎡ 2π 2π ⎤
⎢ cos(θe ) cos(θe − ) cos(θe + )
⎡i sd ⎤ 3 3 ⎥ ⎡i sa ⎤
⎢ ⎥ 2⎢ 2π 2π ⎥ ⎢ ⎥
⎢i sq ⎥ = ⎢− sin(θe) − sin(θe − ) − sin(θe + )⎥ i (II.41)
⎢ ⎥ 3⎢ 3 3 ⎥ ⎢ sb ⎥
⎢ ⎥
⎢⎣i s 0 ⎥⎦ ⎢ 1 2 1 2 1 2 ⎥ ⎣i sc ⎦
⎢⎣ ⎥⎦

59
 CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

Devient alors: le système d’équations

[v dqh] = R s[Ρ] [ ] [
t d Ψ dqh d [ Ρ]]
t
[ ]
[ Ρ]
t t
Ι dqh + [ Ρ]
dt
+
dt
Ψ dqh
(II.42)

Ou [ Xdqh] désigne tout vecteur de grandeurs exprimées dans le référentiel lie au rotor. d,q,h
les indices d, q et h désignent respectivement les composantes directe, en quadrature et
homopolaire.

La multiplication des deux membres de (II.6) par [P] nous donne:

] [ ] [ ]
d [ Ρ]
t
[
V dqh = Rs Ι dqh +
d Ψ dqh
dt
+ pΩ
dθ
Ψ dqh [ ] (II.43)

Avec

⎡0 − 1 0 ⎤
d [ Ρ]
t
[Ρ] = ⎢⎢1 0 0⎥⎥
(II.44)
dθ
⎢⎣0 0 0⎥⎦
Ce qui donne:

⎡ Ld 0 0⎤ ⎡0 − Lq 0⎤
[ ] [ ] ⎢
v dqh = Rs Ι dqh + ⎢ 0 Lq
d
dt
[ ]
0 ⎥⎥ Ι dqh + pΩ ⎢⎢ Ld 0 [ ]
0⎥⎥ Ι dqh (II.45)
⎢⎣ 0 0 Lh ⎥⎦ ⎢⎣ 0 0 0⎥⎦

Ou :

⎡ Ld 0 0 ⎤
(II.46)
⎢0 0 ⎥⎥ = [Ρ ][L ][ Ρ]
t
⎢ Lq
⎢⎣ 0 0 Lh ⎥⎦

60
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

Avec:

⎧ 1
⎪ Ld = L0 + Lf − M 0 + 2 L 2 + M 2
⎪
⎪ Lq = L0 + Lf − 1
⎨ M 0 − L2 − M 2 (II.47)
⎪ 2
⎪ Lh = L 0 − M 0
⎪
⎩

Le neutre de la machine étant isole, ce qui implique naturellement ih=0,on peut écrire:

⎡vd ⎤ ⎡ Rs − pΩLq ⎤ ⎡id ⎤ ⎡ Ld 0 ⎤ d ⎡id ⎤

⎢ ⎥=⎢ + (II.48)
⎢⎣vq ⎥⎦ ⎣ pΩLd Rs ⎥⎦ ⎢⎣ iq ⎥⎦ ⎢⎣ 0 Lq ⎥⎦ dt ⎢⎣ iq ⎥⎦

Ou encore, sous la forme d’équations d’état:

⎡ − Rs pΩLq ⎤ ⎡ 1 ⎤
0 ⎥ vd
d ⎡id ⎤ ⎢ Ld Ld ⎡id ⎤ + ⎢ Ld
⎥ ⎡ ⎤
=⎢ − Rs ⎥ ⎢iq ⎥ ⎢ 1 ⎥ ⎢ vq ⎥
(II.49)
dt ⎢⎣iq ⎥⎦ ⎢− pΩLd ⎥⎣ ⎦ ⎢ 0 ⎥⎣ ⎦
⎣⎢ Lq Lq ⎦⎥ ⎣⎢ Lq ⎦⎥

Equations mécaniques
Le couple électromagnétique est dérivé, selon le principe des travaux virtuels, de l’expression
de la coénergie et en régime linéaire de fonctionnement:

t ⎡ ∂L ⎤
p [ Iabc] ⎢ ⎥[Iabc]
1
cem =
2 ⎣ ∂θ ⎦ (II.50)

61
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

Le remplacement du vecteur courant par sa transformée de Park nous donne:

⎧ t ⎡ ∂L ⎤
⎪cem = 2 [ Idqh] [ P] ⎢ ∂θ ⎥[Idqh]
1 t
⎪ ⎣ ⎦
⎪ ⎡ 0 − Lq 0⎤ (II.51)
⎪ 1 t⎢
⎨= p [ Idqh] ⎢ Ld 0 0⎥⎥[Idqh]
⎪ 2 ⎢⎣ 0 0 0⎥⎦
⎪
⎪= p ( Ld − Lq )idiq
⎪
⎩

On associe à l’équation du couple électromagnétique, la relation fondamentale de la

dynamique des systèmes en rotation :

dΩ
J + fΩ = cem − cr (II.52)
dt

Où :
Ω : représente la vitesse de rotation du rotor,
J : représente le moment d’inertie du système,
f :le coefficient de frottements visqueux.
Cem : est le couple électromagnétique
Cr : le couple résistant.

Le modèle électrique équivalent de la machine dans le repère commun du rotor est représenté
figureII.26. Dans ce repère, on obtient un modèle relativement simple où la machine est
représentée par deux circuits électriques couplés, l’un suivant l’axe d et l’autre suivant l’axe
q. Dans la suite, nous ne considérerons pas les équations relatives à la composante
homopolaire (machine couplée en étoile sans neutre). Les équations des tensions statoriques
aux bornes des bobinages équivalents d’axe d et d’axe q s’écrivent :

62
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

dψsd
U sd = R s . i sd + − ωe.ψsq
dt
(II.53)
dψsq
U sq = R s . i sq + + ωe.ψsd
dt

[Link].26 - Modèle équivalent dans le repère (dq) .

Les flux totalisés dans les enroulements statoriques sont liés aux courants par les relations
suivantes :

ψsd = Ld .isd + Md .ird (II.54)

ψsq = [Link] + [Link]

63
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

De la même façon, on peut écrire les équations des tensions rotoriques (les enroulements
équivalents sont en court-circuit) :

dψrd
0 = Rrd .ird +
dt (II.55)
dψrq
0 = [Link] +
dt

Les flux totalisés dans les enroulements rotoriques s’écrivent :

ψrd = Lrd .ird + Md .isd (II.56)

ψrq = [Link] + [Link]

On montre que l’expression du couple électromagnétique développé par la machine peut

s’écrire de la façon suivante :

cem = Ρ.(ψsd .isq − ψ[Link] ) (II.57)

64
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

Le modèle électromécanique résultant

La MSRV sous la forme d’équations d’état:

⎡ − Rs pLq ⎤ ⎡ 1
id + iqΩ ⎤
⎡ ⎤
id ⎢ Ld Ld ⎥ ⎢ 0 ⎥
⎢ ⎥ ⎢ − Rs ⎥ ⎢ Ld
d ⎢iq ⎥ ⎢ iq −
pLd
idΩ ⎥ ⎢ 0 1 ⎥ ⎡vd ⎤
=⎢ ⎥
⎥+⎢
(II.58)
dt ⎢Ω ⎥ ⎢
Lq Lq Lq ⎥ ⎢⎣ vq ⎥⎦
⎢ ⎥ p ( Ld − Lq )idiq f ⎥
θ ⎢ − Ω − Tr ⎥ ⎢ 0 0 ⎥
⎣ ⎦ J J ⎢ ⎥
⎢ Ω ⎥ ⎣ 0 0 ⎦
⎣ p ⎦

Note : Le modèle électrique que nous avons rappelé est valable dans le cadre des hypothèses
définies précédemment (modélisation au sens du premier harmonique) sans tenir compte de la
saturation magnétique des matériaux.
Le modèle au sens du premier harmonique d’espace est suffisant pour représenter le
fonctionnement d’une machine à réluctance lorsque les enroulements statoriques sont couplés
en étoile à neutre isolé.

65
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

Prototype de la machine étudiée : [10]

Les avantages et inconvénients d’une MSR [5] [14]

Avantages Inconvénients

*Rotor passif permettant un *Nécessité d’un rapport de saillance élevé, donc

fonctionnement à vitesse élevée, fabrication délicate en grande série,

*Rendement relativement meilleur par % à *Facteur de puissance relativement faible,

celui des machines asynchrones. *Electronique de commande assez compliquée
(nécessité d’un capteur de position).

66
CHAPITRE II Choix et étude théorique d’une structure à réluctance variable

Conclusion

Ce chapitre nous a permis de mettre en évidence le principe de conversion d’énergie dans les
machines à réluctance variable (dans la famille des machines à réluctance variable la MSRV
bien sur) cette machine présente un certain nombre d’avantages faisant d’elle une alternative
viable pour les entraînements électriques.

Les différentes structures proposées cherchent à accentuer la dissymétrie entre l’axe direct et
l’axe en quadrature, qui est à l’origine même du couple de réluctance. L’amélioration de
performance revient alors à minimiser l’inductance en quadrature tout en essayant de garder
une inductance directe la plus élevée possible ;

Et comme la MRV a fait l’objet d’un regain d’intérêt notoire dans le domaine des véhicules
électrique et la production d’énergie. ces dernières années la plupart des recherches ont été
basés sur l’optimisation de la structure «le rapport de saillance » ceci afin d’obtenir avec ce
type de machine, des performances comparables à celles de la machine asynchrone.
Nous avons opté en fin de compte sur la MRV 6/4 et la MSRV avec un rotor axialement
laminé. Ces modèles de machines ont été orientés vers la simulation et le diagnostic des
défauts.

67
 Références du Chapitre II

Références du Chapitre II

[1] Bernard. MULTON, «Historiques des machine électromagnétiques et plus

particulièrement des machines à réluctance variable», Article, Ecole Normale Supérieure de
Cachan , 17 Mai 1994 .

[2] Alain. Mailfert, Françaois. Michel SARGOS, « Machines à réluctance variable (MRV)
Principes des MRV. Machines à commutation», D 3680, L’institut national polytechnique de
Lorraine.

[3] Jean Claude. Mouchoux, «Etude et réalisation de l’alimentation électronique d’un moteur
à réluctance variable pour vheicule électrique expérimentation du moteur», Mémoire
d’ingénieur CNAM, 18 Octobre. 1994.

[4] Bernard. Multon, «Conception et alimentation électronique des machines à réluctance

variable à double saillance», Habilitation, INS Cachan, 17 Mai 1994.

[5] Yacine. AMARA, «Contribution à la conception et à la commande des machines

synchrones à double excitation Application au véhicule hybride », Thèse Doctorat, Univ de
PARIS, 21 December 2001.

[6] Boudouda. MALIK, «Commande de la machine à réluctance variable en vue d’une

application alterno démarreur intégré», Thèse Magister en Electrotechnique, Univ de Batna,
17 December 2006.

[7] B. Multon, « Moteur pas à pas », Note de cours Agrégation génie électrique ,ENS Cachan
Antenne de Bretagne.

[8] Matthieu. LEROY, «Etude et mise au point de moto ventilateurs a hautes performances
pour laéronautique », Thèse Doctorat, Univ de Toulouse, 15 Novembre2006.

68
 Références du Chapitre II

[9] Serhoud. Hicham, «Contribution à l’étude de la machine synchrone à réluctance variable»,

Thèse Magister en Electrotechnique, Univ de Biskra, 01 Juillet 2009.

[10]Tsarafidy. RAMINOSOA, «Optimisation des performances des machine synchro-

réluctantes par réseaux de permeances », Thèse Doctorat, INPL Dea Protee, 05 Octobre
2006.

[11] Ned. Mohan , «Flywell batery with a synchronious reluctance motor/ generator», Article,
Univ de Minnesota, September 2004.

[12] Thierry. Lubin , «Modélisation et commande de la machine synchrone à réluctance

variable. Prise en compte de la saturation magnétique», Thèse Doctorat, Univ de Henri
Poincaré, Nancy-I, 18 Avril 2003.

[13] Mohand Ouramdane. HAMITI, «Réduction des ondulations de couple d’une machine
synchrone à réluctance variable Approches par la structure et par la commande», Thèse
Doctorat, Univ de Henri Poincaré, Nancy-I, 15 Juin2009.

[14] T. J. Miller, senior Member IEEE, Alain Hutton, Caulum Cossar, and David A. Staton
«Design of synchronous reluctance motor driver», IEEE Transactions on Industrial
Applications, VOL.27, NO.4, Juillet / Octobre 1991.

[15] David A. Staton, Wen L. Soong, Member IEEE, Timothy J.E. Miller, senior Member
IEEE, «Unified theory of torque production in switched reluctance and synchronous
reluctance motors», IEEE Transactions on Industrial Applications, VOL.31, NO.2, Mars/Avril
1995.

69
CHAPITRE : III

ANALYSE DE DEFAUTS
PAR
ELEMENTS FINIS

(SOUS FEMM ET FLUX 2D)

70
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

INTRODUCTION

L’étude de tout système physique nécessite une modélisation. Celle-ci nous permet de simuler
le comportement de ce système face à différentes sollicitations et d’appréhender ainsi les
mécanismes régissant son fonctionnement.
L’objectif de toute démarche est d’affiner le modèle en analysant avec précision les
phénomènes électriques et magnétiques. On distingue souvent deux méthodes: analytique et
par éléments finis (parfois une troisième approche consiste à coupler les deux !).

La modélisation permet de s’affranchir de ces difficultés si l'on possède des modèles de

machine électrique dédiés au diagnostic, c'est-à-dire des modèles de fonctionnement sains
mais aussi des modèles de fonctionnement en présence de défaut. En effet, les méthodes de
modélisation jouent actuellement un rôle important dans la capitalisation des connaissances
de l'état de la machine, du procédé et de la caractérisation des défaillances. Il existe à
ce jour plusieurs méthodes de modélisation adéquates comme la méthode des éléments finis,
la méthode des réseaux de perméances ou des méthodes analytiques. Grâce à ces
approches, il est possible de décrire le comportement de la machine en présence d'un défaut et
de permettre ainsi de caractériser son influence.

71
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.1 Choix et méthode d’étude [1] [2]

La méthode analytique permet d’exprimer les principes physiques mis en jeu lors de la
conversion électromécanique, sous forme de relations analytiques. Les résultats obtenus sont
approchés par la simplification des hypothèses de travail, nécessaire pour le calcul, ce qui
pénalise limite la précision. De plus, pour choisir ces hypothèses, il est souvent nécessaire de
connaître à l’avance les interactions entre les différents constituants.

La résolution de l’équation (III.1) par des méthodes analytiques est impossible du fait de la
géométrie complexe des machines et des propriétés non linéaires des matériaux magnétiques
les constituant. Seules les méthodes numériques peuvent être utilisées. Le domaine d’étude est
alors discrétisé en un nombre finis de points.

a- Les équations de Maxwell

r r r ∂D
∇×Η = j + (III.1)
∂t
r
r r ∂B
∇× E = − (III.2)
∂t
r r (III.3)
∇.B = 0

r r
∇.D = ρ (III.4)

r r r r (III.5)
B = µ 0 ( H + M ) + Br

r r r (III.6)
B = B(H )

r r r (III.7)
J = σE + J s

Les relations (III.1), (III.2), (III.3), (III.4) représentent les équations de Maxwell.
Les relations (III.5), (III.6) représentent les lois de comportement magnétique.

72
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

Les courants de déplacement sont négligeables devant les courants de conduction, la loi
d’Ampère (III.1) s’écrit alors :

r r r
∇×Η = j (III.8)

La loi de Gauss s’écrit (dans le cadre des systèmes étudiés, on peut considérer que les charges
volumiques sont négligées)

r r
∇.D = 0 (III.9)

Le potentiel vecteur magnétique A est défini par :

r r r
B = ∇. A
(III.10)

Pour garantir l’unicité de la solution nous utiliserons la jauge de Coulomb, soit

r r
∇. A = 0 (III.11)

La substitution de (III.10) dans (III.2) donne :

r
r r ∂A
∇ × (E + ) = 0 (III.12)
∂t

Que l’on peut écrire en introduisant le potentiel scalaire φ :

r
r ∂A r (III.13)
E+ = −∇φ
∂t
r
r ∂A r (III.14)
E=− − ∇φ
∂t

La substitution de par son expression (III.14) dans l’équation (III.7) donne :

r
r ∂A r r (III.15)
J = −σ − σ∇φ + J s
∂t

73
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

A partir de (III.5), (III.8), (III.10), (III.15), l’équation aux dérivées partielles qui régit le
problème s’écrit :

r
r 1 r r r rr ∂A r r r r 1 r
∇ × ( ∇ × Α ) = J s − σ∇φ − σ + ∇ × M ( A) + ∇ × ( B r ) (III.16)
µ0 ∂t µ0

De (III.9) on peut déduire :

r r (III.17)
∇.J = 0

L’introduction de (III.15) dans (III.17) et en tenant compte de (III.11) donne :

r r (III.18)
∇.(σ∇φ ) = 0

La résolution des équations (III.16) et (III.18) en bidimensionnel. Dans ces conditions, le

système est contenu dans le plan OXY , et on suppose qu’il est invariant dans la direction OZ.
La forme de la solution recherchée vérifie l’équation

r r ∂ ( x, y , t )
∇. A = Az =0 (III.19)
∂Z

Le fait que J ait une seule composante selon la direction OZ nous permet de choisir doncφ =0.
Ainsi la relation (III.16) s’écrit

∂ ⎛⎜ 1 ∂A ⎞⎟ ∂ ⎛⎜ 1 ∂A ⎞⎟ ∂A ⎛⎜ ∂ M y ∂ M x ⎞⎟ ⎛⎜ ⎛⎜ ∂ 1 ⎞⎞ ⎛ ∂ ⎛ 1
⎟−⎜ ⎜
⎞⎞
⎟⎟
+ = − J s +σ − − − ( B ) ⎟ ( B )
∂x ⎜⎝ µ 0 ∂x ⎟⎠ ∂y ⎜⎝ µ 0 ∂x ⎟⎠ ∂t ⎜⎝ ∂x ∂y ⎟⎠ ⎜⎝ ⎜⎝ ∂x µ 0 r y ⎟⎠ ⎟⎠ ⎜⎝ ∂y ⎜⎝ µ 0 r x ⎟⎠ ⎟⎠

……..(III.20)

La modélisation analytique sera l’étude du chapitre IV.

74
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.2. Les méthodes numériques [3]

Nous distinguons deux techniques principales de résolution des équations des champs
électromagnétiques :

a) Par Différences finies (D.F) : le maillage est, souvent, un quadrillage rectangulaire sur les
nœuds pour lesquels est effectuée la discrétisation de l’équation différentielle associée à la
décomposition en série de Taylor du potentiel scalaire.

b) Par Eléments Finis (E.F) : autorisent l’emploi d’éléments de toute forme à l’intérieur
desquels le potentiel est approché par un polynôme ; la résolution se ramène à la minimisation
d’une fonctionnelle liée à l’énergie emmagasinée dans les éléments.
La méthode des éléments finis permet de représenter finement les différentes grandeurs dans
le système. Cependant, les résultats obtenus restent « gourmands » en ressources
informatiques, et ne peuvent donc pas être directement exploités dans les structures actuelles
de commande ou de diagnostic. La détermination des champs électromagnétiques est d’un
grand intérêt pour la conception des machines électriques. L’analyse de ces phénomènes
électromagnétiques est basée sur la résolution des équations de Maxwell dont la complexité
est telle qu’elle est exclue dans la plupart des cas du formalisme analytique.

III.3. Le principe du calcul par éléments finis

Le principe du calcul par éléments finis est de diviser la structure de la machine étudiée en
nombre important d'éléments de dimensions finies puis de résoudre les équations de Maxwell
sur chacun de ces éléments. Les conditions aux limites de chaque élément sont fixées par les
éléments voisins. La combinaison de tous ces éléments de calcul permet alors de connaître
l'état magnétique de la structure complète et donc d'effectuer des calculs de flux, de force, et
d'inductance...L’intérêt cette méthode de calcul est qu'elle permet, comme tout calcul
théorique, d'avoir accès aux trois phases de manière indépendante. La méthode des éléments
finis discrétise une formulation intégrale de l’équation aux dérivées partielles pour conduire à
un système d’équations algébriques qui fournit une solution approchée du problème étudié.
Le domaine d’étude est décomposé en un nombre fini d’éléments polygonaux qui forment le
maillage. La valeur du potentiel vecteur est déterminée sur tous les sommets des polygones
(les sommets sont appelés les nœuds du maillage). En employant des fonctions d’interpolation
appropriées, la solution en tout point du domaine sera déterminée en fonction des valeurs aux
sommets de l’élément. Pour transformer un système d’équations aux dérivées partielles par

75
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

une formulation intégrale, les processus les plus souvent utilisés sont la méthode des résidus
pondérés et la méthode variationnelle.

Nous étudions dans ce chapitre la modélisation des deux machines en régime statique. Cette
étude permet d’avoir les caractéristiques magnétostatique de chaque machine. pour y parvenir
nous avons mis en place le code de calcul le logiciel FEMM, simple à utiliser , moins
adaptable mais plus rapide à mettre en œuvre, permettant de calculer les différentes grandeurs
pour plusieurs positions sans pour autant être obligé à chaque fois de redéfinir la géométrie et
le maillage ;

Une seconde étude à été réalisée en régime dynamiques, à l’aide du logiciel FLUX 2d, avec
lequel deux modélisation on été réalisées l’une pour le cas sain, la seconde pour le cas de
défaut.

III.3.1. Le logiciel FEMM est une suite de programmes modulaires utilisés pour résoudre les
problèmes magnétostatique et électrostatique en deux dimensions (axisymétrique ou plan) ;
le femm peut divisé en trois parties essentielles : [4]

o Preprocessor ([Link]) : C'est un programme de CAO pour la définition de

géométrie du problème à résoudre, et définissant les conditions aux limites. Les
fichiers DXF d'Autocad peuvent être importés pour faciliter des géométries.
o Solver([Link]) : Le solver pend en compte les fichiers de données décrivant le
problème et résout les équations de Maxwell pour obtenir les valeurs de l'induction
magnétique dans domaine de résolution.
o Postprocessor ([Link]) : C'est un programme graphique qui affiche les
résultants sous forme de champ, de conteur et densité. Le programme permet aussi à
l'utilisateur d'inspecter le champ aux points arbitraires, et d'évaluer plusieurs
différentes intégrales et tracer diverses quantités d'intérêt le long des contours définis
par l'utilisateur. Deux programmes supplémentaires sont aussi appelés pour exécuter
des tâches spécialisées. Ceux-ci sont:

1. [Link] : Triangle découpe la région avec en grand nombre de triangles, c'est une
partie vitale du procédé d'éléments finis.
2. [Link] : Ce programme est utilisé pour afficher les diverses solutions en 2D.

76
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

Pour traiter des problèmes magnétodynamiques on doit intégrer un sous programme par
langage LUA .
Lua est un langage de programmation étendu, conçu pour des procédures générales de
programmation avec des facilités de description des données, Lua est prévu pour être employé
comme un langage puissant et facile à concevoir, pour différents programmes. Lua est
implémenté comme une bibliothèque en C.
Le script LUA est une partie de programme directement interprétée par FEMM, contenant des
fonctions propres au logiciel FEMM.
L’utilisation de la méthode numérique de calcul par éléments finis nous apporte plusieurs
avantages : la prise en compte de la géométrie de la machine, de la saturation des matériaux
magnétiques, ainsi que l'effet de peau dans les barres rotoriques. D'autre part, d'un point de
vue pratique, il est très facile de faire varier les conditions de fonctionnement de la machine
(fréquence et l’alimentation), en changeant simplement les données du programme de calcul
d’où une grande souplesse d’utilisation .Ainsi ce code de calculs permet de calculer les
différentes grandeurs pour plusieurs positions sans pour autant être obligé à chaque fois de
redéfinir la géométrie et le maillage. La modélisation par élément finis permet d’étudier avec
précision l'effet des défauts sans autant détruire la machine.

III.3.2. Le logiciel FLUX 2d [6] [7]

Le logiciel permet de calculer et de visualiser les grandeurs utiles à l’ingénieur, pour des
dispositifs bidimensionnels ou à symétrie de révolution comportant des matériaux à
caractéristiques linéaires ou non, isotropes ou non. C’est un logiciel complet ayant l’avantage
de permettre le couplage avec les équations de circuits ainsi que l’ajout d’une région
surfacique particulière dite “ bande de roulement ” pour l’étude des machines tournantes avec
différentes positions du rotor, sans avoir à modifier la géométrie et le maillage.

La résolution d’un problème fait appel à des modules spécialisés:

¾ Module pré-processeur
Permet, à travers plusieurs modules, de définir la géométrie du dispositif à étudier (PREFLU),
de choisir et/ou de construire une banque de matériaux (CLSMAT), d’affecter les propriétés
physiques aux différentes régions géométriques prédéfinies (PROPHY) et de définir le

77
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

schéma et/ou les données du circuit électrique (CIRFLU). Il permet également un maillage
automatique d’une géométrie 2D prédéfinie (PREFLU).

¾ Module processeur
Constitué principalement d’un module de résolution 2D (SOLVER_2D) des différents
modèles usuels de l’électromagnétisme et des problèmes thermiques.

¾ Module post-processeur (PostPro_2D)

Permet, entre autres, de tracer les équipotentielles ou les lignes de flux, le maillage, la
géométrie et les courbes 2D ou 1D selon un chemin prédéfini. Il permet aussi de calculer des
grandeurs globales telles que le couple ou la force appliqués à un contour fermé, les
inductions, les flux, les inductances, etc.

III.4. Modélisation par élément finis

III.4.1. Régime statique par la mise en œuvre du logiciel FEMM

III.4.1.1 Définition et maillage de la structure d’une MRV 6/4

Avant de modéliser la machine avec défauts en 2D, nous mènerons une étude préliminaire sur
le fonctionnement de la machine (sain puis avec défauts).On peut exploiter la symétrie
centrale géométrique de la machine étudiée et l’antisymétrie de son alimentation par rapport
à l’axe central pour réduire la région à mailler, et par suite le temps de calcul ; cela consiste à
travailler sur un seul pôle et à et à utiliser des conditions aux limites anti-périodique.

78
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.[Link] Dimensions géométriques de la MRV [8]

Les principales dimensions du prototype dans la tableau III.1

Intitulé Notation valeur

Longueur active L 40 [mm]

Rayon extérieur Rext 60 [mm]
Rayon statorique Rs 29.5 [mm]
Rayon rotorique Rr 29 [mm]
Hauteur des dents statoriques hs 21.5 [mm]
Hauteur des dents rotoriques hr 11 [mm]
Epaisseur d’entrefer e 0.5 [mm]
Epaisseur du culasse statorique ec 9 [mm]
Arc polaire statorique Bs 30°
Arc polaire rotorique Br 30°
Rayon de l’arbre Rr 9 [mm]
Nombre de spire par pole Ns 120

III.[Link]. Matériau utilisé

Le matériau utilisé pour construire les tôles du stator et rotor c’est le fer.

Caractéristiques
La tension V=800 volt.
Epaisseur = 0.5 mm.
Perte a 1.5 T = 8W/Kg.
Le fer choisi afin de construire le prototype le moins coûteux ;

79
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

Fig. III.1 - Courbe de magnétisation de fer .

Maintenant pour le fil de bobinages statoriques :

Diamètre Facteur de remplissage Maximum efficace de la densité de N de spire / pole enroulement
du cuivre
du fil courant
1.07 0.4 10 A/mm^2 120
mm

L’arbre et le stator ont des dimensions standards comme celle d’une machine asynchrone.

III.[Link]. La machine modélisée par le FEMM

Stator

Les encoches
statorique

rotor

Dents statoriques

80
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

Fig. III.2 - MRV 6/4 (2-D Planaire ,longeure utile 40 mm)

III.[Link]. Le maillage
Le maillage complet comporte ainsi 7688 nœuds et 15302 éléments .

81
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.[Link]. Exploitation des résultats

a) Les lignes iso-potentiels vecteurs

b) La densité du flux magnétique

c) Le vecteur d’induction

82
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.[Link] Caractéristiques électromagnétiques

a) Caractéristique magnétique du flux ϕ(Ni,θ) :

La caractéristique magnétique du flux en fonction des ampères tours et la position du rotor
obtenue est représentée par figure III.3 ,pour des positions variant de 0 à 45° par pas de 5°.

flux=f(teta;i)
0.4

téta =45° CONJONCTION +45°

0.35

0.3 +35°

+30°
0.25
flux [Wb]

+25°
0.2
+20°
0.15
+15°
0°
0.1

0.05
téta=0° OPPOSITION

0
0 2 4 6 8 10 12 14 16 18 20
i [A]

Fig. III.3 Le flux magnétique

b) Caractéristique magnétique de l’inductance L(θ,Ni) :

inductance=f(teta;i)
0.055
2A
0.05 5A

0.045

0.04
8A
0.035
inductance [H]

0.03 11 A

0.025 14 A

0.02 17 A
20 A
0.015

0.01

0.005
0 5 10 15 20 25 30 35 40 45 50
teta [degres]

Fig. III.4 L’inductance

83
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

c) Caractéristique magnétique du couple ce(θ,Ni)

couple=f(teta;ni)
9
i=20 A
8
i=17 A
7

6 i=14 A
couple costa [Nm]

i=11 A
4

3
i=8 A
2
i=5 A
1
i=2 A
0

-1
0 5 10 15 20 25 30 35 40 45
teta [degres]

Fig. III.5- Le couple statique par phase

84
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.4.1.2. Définition et maillage de la structure d’une MSRV

III.[Link]. Dimensions géométriques de la MSRV [9]

a) Les dimensions du stator b) Les dimensions du rotor à barrières flux

[Link].6- Structure et paramètres géométriques de la machine synchro-réluctante
étudiée.

Le tableau (III-2) présente les paramètres associatif aux prototypes.

Intitulé Notation Valeur

Longueur utile L 200 [mm]

Nombre de pôles P 2
Diamètre extérieur du stator Rs 450 [mm]
Diamètre du rotor Rr 282.12 [mm]
Epaisseur de segments ferromagnétiques epfer 13.48 [mm]

Epaisseur des barrières de flux ecfer 21.61 [mm]

Nombre de barrières de flux Nlam 6

Epaisseur d’entrefer e 2 [mm]
Nombre d’encoches Ns 36
Nombre d’encoches par pôle et par phase m 6
Nombre de spires par encoche et par phase na 3
Hauteur d’une dent hd 24.09 [mm]
Hauteur de l’épanouissement d’une dent Embase rectangulaire : b=0.5 2.5 [mm]
Embase trapézoïdale : f=2

Largeur d’une dent ld 10.27 [mm]

Ouverture d’encoche αep 2°
Pas dentaire τd 10°

85
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.[Link]. Matériau utilisé

Matériau utilisé Acier AISI 1010.
Coefficients de pertes fer dans le modèle de Bertotti :
Hystérésis : KH = 15,45.10^-3 [SI].
Supplémentaires : KE=3.2.10^-3 [SI].
Résistivité : 17,6.10^8 [Ωm].
Epaisseur d’une tôle : 0,66.10^-3 [m].
Masse volumique : 7,85.10^3 [kg/m^3].

[Link].7- Caractéristique B=f(H) de l’acier 1010.

86
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.[Link] La machine modélisée par le FEMM

Statotr steel 1010

Rotor axialement
laminé

Materiau
amagnétique Al

36 Encoches en ciuvre

III.[Link] Le maillage
Le maillage complet comporte 3842 nœuds et 7610 éléments .

steel

Air Air Air

Air Air Air Air
Air Air
Air
Air
Air
coil
coil [+A:11]
[-A:11] coil
coil [+A:11]
[-A:11] coil
coil steelsteelsteelsteel
steel
steel [+A:11]
[-A:11] steel coil
coil Air [+A:11]
[-A:11] Air Al Al Al Al Al Al coil
coil [+A:11]
[-A:11] coil
coil [+A:11]
[-A:11] Air
Air
Air Air
Air Air
Air Air
Air Air Air Air Air

87
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.[Link] Exploitation des résultats

a) Les lignes iso-potentiels vecteurs

b) La densité du flux magnétique

88
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

c) Le vecteur d’induction

III.[Link]. Caractéristiques électromagnétiques

a) Caractéristique magnétique du flux Ψ(Ni,θ)

flux=f(teta;i)
1.8

1.6

1.4

1.2

1
flux [Wb]

0.8

0.6

0.4

0.2

0
0 50 100 150 200 250 300 350 400 450
i [A]

[Link].8- le flux total magnétique

89
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

c) Caractéristique magnétique du couple ce(θ,Ni)

couple=f(teta;ni)
1200

1000

800
couple [Nm]

600

400

200

0
0 10 20 30 40 50 60 70 80 90
teta [degres]

[Link].9- Le couple statique par phase

La figure III.9 représente le couple en fonction de l’angle de charge avec l’amplitude Im du

courant triphasé en paramètre. Nous pouvons remarquer que la limite de stabilité se décale
légèrement sensiblement vers les angles de charge plus élevés à mesure que l’amplitude du
courant Im du courant triphasé augmente. Elle se situe à 65° pour un courant égal au courant
nominal.
On assiste à une forte désaturation de la machine quand l’angle de charge varie de 0 à 90°.
De ce fait, la différence (Ld-Lq) est croissante quand l’angle de charge varie de 0 à 90. Le
maximum du produit (Ld-Lq).sin(2θο), se trouve ainsi repoussé loin au delà de 45°.

90
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.4.2. REGIME DYNAMIQUE PAR LA MISE EN ŒUVRE DU LOGICIEL FLUX 2D

L’étude du comportement des systèmes électromagnétiques dépend des conditions de leur

fonctionnement. En effet, plusieurs formulations sont possibles tel que: l’électrostatique,
l’électrodynamique, le magnétostatique, le magnétodynamique, etc. Vu que les moteurs à
réluctance variable comporte une dynamique que nous ne pouvons négliger lors de l'étude,
alors notre choix s'est porté sur une étude en magnéto-évolutif ; ceci permet de suivre le
comportement transitoire et temporel de la machine.

III.4.2.1. LA MRV 6/4

III.[Link]. La géométrie

Stator
B- C+

B+ C-

A+
A+

L’arbre Rotor

III.[Link]. Le maillage

91
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.[Link]. Le maillage dans l’entrefer

Une seule couche

Le maillage dans l’entrefer doit être bien réglé, Car on doit étudier la variation des grandeurs
électromagnétiques et mécaniques en fonction de la position du rotor afin d’avoir les
caractéristiques de la machine modélisée , on fait tourner d’un angle quelconque le rotor,
cependant au fur et mesure que l’angle de déplacement augmente, la distorsion des éléments
de la région de l’entrefer (la bande de roulement) augmente aussi, ce qui provoque des
difficultés d’ordre numérique donc le maillage de la région de l’entrefer doit être régulier (la
méthode des éléments finis donne de bons résultats avec des éléments réguliers).
l'utilisation de la bande de roulement, une fonction du logiciel FLUX 2d, nous a permet de
considérer la rotation du rotor en étude magnéto-évolutif sans pour autant effectuer un
nouveau maillage de la machine à chaque position du rotor.

III.[Link]. Les conditions aux limites

_
_

dirichlet

Bounds conditions
User's dirichlet
Dirichlet
Anticy clic
ref of anticyclic condition
Cy clic
ref of cy clic condition
Periodic
ref of periodic condition
Translation
ref of translation condition
Floating

anticuclique
bande
roulement

92
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.[Link]. Le circuit d’alimentation

La partie alimentation de la machine est représentée par un circuit électrique représentant le

bobinage statorique et son alimentation ( Figure. III.10).

MRV_SAINE

Switch1
A+ Switch3
Switch2

Switch4

B- Switch6
B+
Switch5

Switch7
Switch9
C- C+ Switch8

Figure. III.4 - L’onduleur d’alimentation

Fig. III.10 - Topologie de l’onduleur par flux 2d

La loi de command envisagée pour piloter ce moteur est représentée sur la figure III.11

93
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

MRV_SAINE

99,999 CURVE va
Circuit / Voltage
Time
BOBA ;

-100
CURVE C2D_22
Circuit / Current
Time
BOBA ;

0,005 0,01

[Link].11 – La tension et le courant de la MRV

III.[Link]. Analyse des résultats

La figure III.12 présente les formes d’ondes caractéristiques du prototype, soit les tension de
phase Va,Vb,Vc ,les courants de phase ia,ib,ic ,ainsi que le couple dynamique. Concernant
de la vitesse de rotation est imposée au début constante de valeur 2500 [tr /mn].

Isovalues Results
Quantity : Equi flux Weber

Time (s.) : 92,5E-6 Pos (deg): -21,611

Line / Value
1 / -271,48286E-6
2 / -242,90571E-6
3 / -214,32857E-6
4 / -185,75143E-6
5 / -157,17428E-6
6 / -128,59714E-6
7 / -100,0191E-6
8 / -71,44285E-6
9 / -42,86571E-6
10 / -14,28857E-6
11 / 14,28857E-6
12 / 42,86571E-6
13 / 71,44285E-6
14 / 100,0191E-6
15 / 128,59714E-6
16 / 157,17428E-6
17 / 185,75143E-6
18 / 214,32857E-6
19 / 242,90571E-6
20 / 271,48286E-6

[Link].12 - La MRV à t=9.25e-005 [s] pour position initiale = 23°.

94
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

MRV_SAINE MRV_SAINE

CURVE C2D_4 30
CURVE C2D_7
Circuit / Voltage Circuit / Current
Time Time
BOBA ; BOBA ;
99,999

20

CURVE C2D_5 CURVE C2D_8

0 Circuit / Voltage Circuit / Current
Time Time
BOBB; BOBB;
10

-100

CURVE C2D_7
Circuit / Voltage CURVE C2D_9
Time Circuit / Current
BOBC; 0 Time
BOBC;
0,005 0,01
0,005 0,01

a1. Les tensions va, vb, vc [V] a2. Les courant ia, ib, ic [A]

MRV_SAINE MRV_SAINE

2,5
(E3) rpm N.m
3

2,5 2

1 CURVE C2D_4
2,5 CURVE C2D_6
Mechanics / Angular velocity Mechanics / Axis torque
Time Time

2,499

-1

(E-3) s. (E-3) s.

5 10
5 10

a3. La vitesse [tr/mn] a4. Le couple moteur [Nm]

MRV_SAINE MRV_SAINE
(E-3) Tesla
(E-6) Tesla 200

100
200
CURVE C2D_2
0 Flux density / Normal component
airgap
Time (s.) : 92,5E-6

-100
150

CURVE C2D_14 -200

mm
Flux density / Magnitude
Time 0 100 200 300 400
100 Point(0,0)
(E-3) Tesla
40

30
50
20
CURVE C2D_3
Flux density / Tangent component
(E-3) s. 10 airgap
Time (s.) : 92,5E-6

0
5 10
-10
mm
0 100 200 300 400

a5. La densité du flux magnétique[A/m] a.6 La composante normale et tangentielle

95
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

MRV_SAINE MRV_SAINE
Volt
-149,5
Ampere
-149,75
CURVE vs
Circuit / Voltage 25
-150
Time
-150,25 V; 20

-150,5
(E-3) s.
5 10 15
CURVE C2D_29
Circuit / Current
Ampere Time
30 Switch1 ;

10
20
CURVE is
10 Circuit / Current 5
Time
0
V;
0
(E-3) s.
-10
(E-3) s.
5 10 5 10

a7. La tension et le courant dans la source V a8. Le courant dans l’interrupteur1

Fig. III.13 – Les caractéristiques de la MRV 6/4.

96
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.[Link]. Modélisation de défauts

La simulation des différents défauts de court circuit (détaille voir le chapitre I), tel que
monophasé, court circuit triphasé, circuit ouvert, les excentricités ont étés réalisés on utilisant
l’option transite start up du flux 2d.

a) A vide et a vitesse constante

1. Un court circuit monophasé dans la phase a de 25 %

La simulation est réalisée par diminution la résistance et de l’inductance de la phase a.

1.1. Les grandeurs globales

_
_
M R V_CO URT C20T ES T

A mpere

30

20
CU R VE C2 D _ 1 4
Circ uit / Cur re nt
T im e
BO BA ;

10

0
(E-3) s .

20 30 40 50

Fig.b1 – Le courant dans la phase court circuité

_
_ MRV_COURTC20TEST

CURVE ia
Circuit / Current
30 Time
BOBA ;

20 CURVE ib
Circuit / Current
Time
BOBB ;
10

CURVE ic
Circuit / Current
0 Time
BOBC ;
0,02 0,03 0,04 0,05

Fig.b2 – Les courants ia,ib et ic

97
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

_
_
MRV_COURTC20TEST

N.m

CURVE C2D_19
Mechanics / Axis torque
Time
3
Augmentation du
couple de la phase
a court circiuté
2

(E-3) s.

20 30 40 50

Fig.b3 – Le couple moteur

1.2. Les grandeurs locales

MRV_VOURTC20TEST MRV_VOURTC20TEST
(E-3) Tesla
(E-6) Tesla
500

150
CURVE C2D_22
0 Flux density / Normal component
airgap
Time (s.) : 13,875E-3

-500

mm
CURVE C2D_21
Flux density / Magnitude 0 100 200 300 400
Time
100 Point(0,0) (E-3) Tesla
100

50

0
CURVE C2D_23
Flux density / Tangent component
-49,999 airgap
Time (s.) : 13,875E-3

(E-3) s. -100

-150
20 30 40 50 mm

0 100 200 300 400

Fig.b4 – Le flux magnétique Fig.b5– La composante normale et tangentielle

1.3. Interprétation des résultats

La figure b2. présente les courants statoriques dans le cas du court-circuit monophasé de 25
%. On constate que l'amplitude du courant dans la phase détériorée augmente par rapport aux
deux autres phases. Le courant de la phase saine présente une valeur égale à 29.16[A] cette
valeur augmente à 33.69[A] l’ors d’un court circuit de 25 % donc un court circuit monophasé
de 25 % fait une augmentation de 4.53 % .

98
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

Toutefois, les courants qui parcourent ces phases sont augmentés par rapport au cas du moteur
sain ,les courants ib = ic= 29.27 [A] en cas sain augmente à 10 %. Notons que même si le
défaut concerne uniquement une phase statorique, toutes les phases sont affectées par le court-
circuit.
Plus important est le nombre de spires en court-circuit, plus forte est l’augmentation des
courants.

Pour l’évolution du couple figure b3. La valeur du couple sain égale à 3.26 [Nm],en cas de
court circuit de 25% la valeur augmente dans la phase infecté à 5.42 [Nm] c’est à dire
augmentation de 66.25% et pour les deux autres phases la valeur augmente à 4.78 [Nm] c’est
à dire 46.625 %.

99
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

2. Un court circuit triphasé de 25 %

2.1. Les grandeurs globales

MRV_CC03PHASE

CURVE ia
Circuit / Current
30 Time
BOBA ;

20 CURVE ib
Circuit / Current
Time
BOBB ;
10

CURVE ic
Circuit / Current
0 Time
BOBC ;
0,02 0,03 0,04 0,05

Fig.c1 – Les courants ia, ib et ic

MRV_CC03PHASE MRV_CC03PHASE

N.m N.m

5 5

4 4

CURVE couple CURVE C2D_12

Mechanics / Axis torque Mechanics / Magnetic torque
Time Time
3 3

2 2

(E-3) s. (E-3) s.

20 30 40 50
20 30 40 50

Fig.c2– Le couple moteur Fig.c3– Le couple électromagnétique

100
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

2.2. Les grandeurs locales

MRV_CC03PHASE MRV_CC03PHASE
(E-3) Tesla
(E-6) Tesla
500

150
CURVE B
0 Flux density / Normal component
airgap
Time (s.) : 13,875E-3

-500

mm
CURVE C2D_14 0 100 200 300 400
Flux density / Magnitude
Time
100 Point(0,0) (E-3) Tesla
100

50

0
CURVE TANGENTE
Flux density / Tangent component
-49,999 airgap
Time (s.) : 13,875E-3

-100
(E-3) s.
-150
mm
20 30 40 50
0 100 200 300 400

Fig.c4– Le flux magnétique Fig.c5– La composante normale

et tangentielle

2.3. Interprétation des résultats

La figure III.c1 présente les courants statoriques dans le cas du court-circuit triphasé de 25 %.
Par apport à la Fig. III.7.b) la valeur des courant ia,ib et ic augmente à 34.85[A] en
pourcentage 19.5%.
La valeur du couple moteur voire une augmentation de 68% qui corresponde à 5.48 [Nm].

101
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

3. Un circuit ouvert dans la phase a

Pour simuler ce de défaut il faut augmenter la résistance de la phase a dans notre cas
Ra=10*Rb .

3.1. Les grandeurs globales statoriques le courant ia de la phase infecté, les courants des
phases ia, ib et ic et les tensions des phases va,vb,vc.

MRV_OPENCIRCIUT MRV_OPENCIRCIUT

Ampere
15 CURVE C2D_12
Circuit / Current
30 Time
BOBA ;

10 20

CURVE ia CURVE C2D_13

Circuit / Current
Circuit / Current Time
BOBB;
Time
BOBA ; 10
5

CURVE C2D_14
Circuit / Current
0 Time
BOBC;
0
(E-3) s.
0,02 0,03 0,04 0,05

20 30 40 50

Fig.d1 – Le courant de la phase A ouverte Fig.d2- Les trois courants de phases

MRV_OPENCIRCIUT

N.m

3
CURVE C2D_18
Mechanics / Axis torque
Time

(E-3) s.

20 30 40 50

Fig.d3– Le couple moteur

102
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

3.2. Les grandeurs locales la composante normale, tangentielle et le flux magnétique

MRV_OPENCIRCIUT MRV_OPENCIRCIUT
(E-3) Tesla
(E-6) Tesla
150
500

CURVE C2D_21
0 Flux density / Normal component
airgap
Time (s.) : 13,875E-3

-500
100
mm

0 100 200 300 400

CURVE C2D_23
Flux density / Magnitude
Time
(E-3) Tesla Point(0,0)
100

50

0 50
CURVE C2D_22
Flux density / Tangent component
-49,999 airgap
Time (s.) : 13,875E-3
-100
(E-3) s.
-150
mm
20 30 40 50
0 100 200 300 400

Fig.d4 - La composante normale et tangentielle Fig.d5 - Le flux magnétique

3.3. Interprétation des résultats

Quand une phase ouverte dans une machine à réluctance variable aucune excitation n'est
disponible à cette phase défectueuse et il n'y a aucune contribution du couple dans cette phase
la figure d3. présente le couple moteur en cas ou la phase A ouvert le couple de la phase A
est nul et pour les deux phases saines le couple augmente à 4.86 [Nm] , par apport un cas sain
le couple fait une augmentation de 49 %.
Le courant dans la phase A diminue à 9.20 [A] en pourcentage 68.449 %, pour une résistance
égale à 10 fois la résistance d’une phase saine, figure d1.
Et les courants dans les phases saines augmente à 33 [A] en pourcentage 13.16 % voire la
figure d2.
La phase peut ouvrir en raison des défauts de fabrication ou grillage d'un lien faible dans
l'enroulement.
La simulation d’un défaut de circuit ouvert peut être conduite en débranchant les signaux de
la gâchette de phases ouverts à l'instant de défaut.

103
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

4. Les défauts au niveau de l’onduleur

Tous les défauts de la machine à réluctance variable au niveaux du stator étant maintenant
simulé, il nous reste de simulé les défauts au niveau de l’alimentation électrique, la machine
est alimentée à travers un convertisseur statique.
L’onduleur qui alimente la machine présente à la figure III.10. ce dernier constitue 09
interrupteurs ,06 diodes, et une source de tension fournie une tension continue égale à 150
[V].
Les différents défauts de l’onduleur (onduleur Arrêté, court circuit monophasé , court circuit
équilibré, circuit ouvert) cités dans le chapitre 1 dans ce chapitre on va présenté le cas ou
l’interrupteur IGTB de la phase A court circuité, les diodes peut supporté une tension
inverse plus grande que les IGTB, donc afin de prévoir la vie d'un IGBT, le nombre de cycles
thermiques que l'IGBT supporte lors du fonctionnement doit être correctement déterminé.

4.1. Un IGBT court circuit

S1 court circiuté

[Link].14 - La topologie de l’onduleur qui alimente la géométrie complète de la mrv.

104
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

4.1.1. Les grandeurs globales statoriques

MRV_INVERTERFAUT

CURVE ia
50 Circuit / Current
Time
BOBA ;
40

30 CURVE ib
Circuit / Current
Time
20
BOBB ;

10
CURVE ic
Circuit / Current
0 Time
BOBC ;
0,02 0,03 0,04 0,05

Fig.e1- Les courants de phases

MRV_INVERTERFAUT

N.m
6

4
CURVE couple
Mechanics / Axis torque
Time

(E-3) s.

20 30 40 50

Fig.e2- Le couple total

105
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

4.1.2. Les grandeurs locales

MRV_INVERTERFAUT MRV_INVERTERFAUT

(E-3) Tesla (E-3) Tesla

100

500

49,999

CURVE tangent e
CURVE C2D_5 0 Flux density / Tangen t compon ent
0 Flux density / Normal component airgap
Time (s.) : 13,875E-3
airgap
Time (s.) : 13,875E-3

-50

-500 -100

mm
mm
0 100 200 300 400

0 100 200 300 400

Fig.e3- La composante normale Fig.e4 - La composante tangentielle

4.1.3. Interprétation des résultats

La figure e1. montre la l’allure des courants statoriques des trois phases superposés pour
mieux voire l’influence de défaut, on vois clairement une grande augmentation égale à 4 fois
l’augmentation due au court circuit triphasé dans le courant de la phase ou l’interrupteur est
court circiuté, et les courants des phases saines ib et ic augmente à 33 [A]

Pour la grandeur du couple fig.e2 non seulement une valeur augmente à 100% qui correspond
à 6 [Nm] mais aussi une déformation du forme par apport le cas sain on peut dire une
modulation causé par un court circuit dans une phase de l’onduleur.

106
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

4.2. Une phase de l’onduleur court circuit

IGBT CC DIODE CC

[Link].15 - Un onduleur avec une phase court circuit

4.2.1. Les grandeurs globales

MRV_INVERTERFAULT02

Ampere
15

10

CURVE ia
Circuit / Current
Time
BOBA ;
5

0 (E-3) s.

20 30 40 50

Fig.f1- Le courant dans la phase a

107
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

MRV_INVERTERFAULT02

CURVE ia
30 Circuit / Current
Time
BOBA ;
20

10
CURVE ib
Circuit / Current
Time
BOBB ;
0

-10
CURVE ic
Circuit / Current
Time
BOBC ;
0,02 0,03 0,04 0,05

Fig. f2 - Les courants de trois phases

MRV_INVERTERFAULT02

N.m

CURVE couple
Mechanics / Axis torque
Time
2

(E-3) s.

20 30 40 50

Fig.f3 - Le couple moteur

108
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

4.2.2. Les grandeurs locales

MRV_INVERTERFAULT02 MRV_INVERTERFAULT02
Tesla
A/m
0,5

100 CURVE normale

0 Flux density / Normal component
airgap
Time (s.) : 13,875E-3

-0,5

mm
CURVE flux
0 100 200 300 400
Field / Magnitude
50 Time Tesla
Point(0,0) 0,1

49,999E-3
CURVE tangente
0 Flux density / Tangent component
airgap
Time (s.) : 13,875E-3

-0,05
(E-3) s.
-0,1
mm
20 30 40 50
0 100 200 300 400

Fig.f4- Le flux magnétique Fig. f 5- la composante normale et tangentielle

4.2.3. Interprétation des résultats

Le cas ou une phase de l’onduleur court circuité [Link].15 le courant de la phase infecté est
devenu nul donc cette phase cesse à produira le couple totale ,l’allure du couple fig.f3 traduira
bien cette observation le couple total est constitue de deux phases b et c, la phase a est
disparue.

Le scénario de défaut dépend de la topologie de l’onduleur, il existe une variété d’onduleur

pour les machine à réluctance variable. [13]

109
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

b) Démarrage à vide avec J=0.0013 et F=0.0183

Cette paragraphe nous permis de voire l’influence de défauts sur la vitesse ,on prend le cas du
court circuit monophasé.

1. Le cas sain

MRV_NOLOAD

N.m

CURVE couple
Mechanics / Axis torque
Time

-5 (E-3) s.

50 100

Fig.1b - Le couple moteur

MRV_NOLOAD

(E3) rpm

1,5

CURVE vitesse
Mechanics / Angular velocity
Time

0,5

(E-3) s.

50 100

Fig.2b - La vitesse

110
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

MRV_NOLOAD

CURVE IA
Circuit / Current
40
Time
BOBA ;

30

CURVE IB
Circuit / Current
20 Time
BOBB ;

10

CURVE IC
Circuit / Current
0
Time
BOBC ;
0,05 99,999E-3

Fig.3b - Les courants de phases

2. Un court circuit monophasé dans la phase a de 25 %

__
MRV_NOLOADCC25

N.m
Phase a court circuitée
10

CURVE COUPLE
Mechanics / Axis torque
Time

-5 (E-3) s.

50 100

Fig.4b – Le couple moteur en cas du court circuit

111
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

MRV_NOLOADCC25

(E3) rpm

1,5

CURVE vit
Mechanics / Angular velocity
Time

0,5

(E-3) s.

25 50 75

Fig.5b – La vitesse en cas du court circuit monophasé

1800

1600

1400

1200 vitesse court ciruit 25%

la vitesse [tr/mn]

vitesse saine
1000

800

600

400

200

0
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
le temps [s]

Fig.6b – La vitesse en cas du court circuit monophasé et saine

112
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.4.2.2. La MSRV avec un rotor Axialement Laminé

III.[Link]. La géométrie

Barres en fer

Stator

C-
B+

A+
A-

Barres en Al

III.[Link]. Le maillage

113
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.[Link]. Le maillage dans l’entrefer

III.[Link]. Les conditions aux limites

Bounds conditions
User's dirichlet
Dirichlet
Anticyclic
ref of anticyclic condition
Cyclic
ref of cyclic condition
Periodic
ref of periodic condition
Translation
ref of translation condition
Floating

114
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.[Link]. Le circuit d’alimentation

MSRV_AXIAL

Switch1 Switch2 Switch3

PA
Switch4 Switch5 PB
Switch6 MC
MA

Fig. III.16- Topologie de l’onduleur par flux 2d

La commande appliqué sur les interrupteurs de l’onduleur est de type décalé on commande
S1 S2 S3 en complémentaires et on décale S6 S4 S5 qui sont commandés aussi en
complémentaires , ce type de commande permet de diminuer le taux d’harmonique .

S1 S2 S3

S6 S4 S5

115
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.[Link]. Analyse des résultats

Nous allons présenter les caractéristiques électromagnétiques globales et locales du prototype
simulé :

Isovalues Results
Quantity : Equi flux Weber

Time (s.) : 666E-6 Pos (deg): -9,999E-3

Line / Value
1 / -860,71162E-6
2 / -770,1104E-6
3 / -679,50919E-6
4 / -588,90792E-6
5 / -498,3067E-6
6 / -407,70549E-6
7 / -317,10428E-6
8 / -226,50305E-6
9 / -135,90184E-6
10 / -45,30061E-6
11 / 45,30061E-6
12 / 135,90184E-6
13 / 226,50305E-6
14 / 317,10428E-6
15 / 407,70549E-6
16 / 498,3067E-6
17 / 588,90792E-6
18 / 679,50919E-6
19 / 770,1104E-6
20 / 860,71162E-6

Figure III.17 - La MSRV en conjonction à t= 0.000666[s]

MSRV_AXIAL MSRV_AXIAL
Volt
50
CURVE IA CURVE VA
Circuit / Current 0 Circuit / Voltage
Time Time
20 PA ;
PA ; -49,999
s.

0,01 0,02 0,03 0,04 0,05 0,06

Volt
10
100
CURVE IB 49,999 CURVE VB
Circuit / Current
0 Circuit / Voltage
0 Time Time
PB ; -50 PB ;
-100 s.

-10 0,01 0,02 0,03 0,04 0,05 0,06

Volt
100
CURVE IC
-20
Circuit / Current 49,999 CURVE VC
Time Circuit / Voltage
0
MC ; Time
0,01 0,02 0,03 0,04 0,05 0,06
-50 MC ;
s.

0,01 0,02 0,03 0,04 0,05 0,06

Fig.g1- Les courants de phases ia, ib, ic Fig.g2- Les tension va, vb, vc

116
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

MSRV_AXIAL MSRV_AXIAL
Ampere
Ampere
20 20

10 CURVE isource
Circuit / Current
Time 15
0 VSOURCE ;

-10
(E-3) s. 10
CURVE ISWITCH
10 20 30 40 50 60 Circuit / Current
Time
Volt Switch1 ;
-149,5 5

-149,75

CURVE vsource 0
-150 Circuit / Voltage
Time
VSOURCE ; (E-3) s.
-150,25

-150,5
(E-3) s. 10 20 30 40 50 60

10 20 30 40 50 60

Fig.g3 - Le courant et la tension Fig.g4 - Le courant l’interrupteur 1

dans la source V

MSRV_AXIAL

N.m

CURVE COUPLE
Mechanics / Axis torque
0 Time

-5 (E-3) s.

10 20 30 40 50 60

Fig.g5 - Le couple moteur

117
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

MSRV_AXIAL

(E-3) Tesla

400

300

CURVE FLUX
Flux density / Magnitude
Time
200 Point(0,0)

100

(E-3) s.

10 20 30 40 50 60

Fig.g6 - Le flux magnétique

MSRV_AXIAL MSRV_AXIAL

(E-3) Tesla (E-3) Tesla

30
10

20

10 4,999

0 CURVE NORMALE CURVE TANGENTE

Flux density / Normal component Flux density / Tangent component
airgap airgap
0 Time (s.) : 333E-6
Time (s.) : 333E-6
-9,999

-19,999 -5

-30

mm -10 mm

0 100 200 300 400 500 600 699,999

0 100 200 300 400 500 600 699,999

Fig.g7- la composante normale Fig.g8 - composante tangentielle

Fig. III.18 – Les caractéristiques de la MSRV.

118
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.[Link]. Modélisation de défauts pour la MRSV

a) A vide et a vitesse constante

1. Un court circuit monophasé dans la phase a de 25 %

1.1. Les grandeurs globales

AXIAL_CCMONOPH25

Ampere

10

0
CURVE ia
Circuit / Current
Time
PA ;

-10

(E-3) s.

100 150 200 250

Fig.h1- Le courant dans la phase court circuit

AXIAL_CCMONOPH25

N.m

2,5

CURVE couple
Mechanics / Axis torque
Time

-2,5 (E-3) s.

100 150 200 250

Fig.h2- Le couple moteur

119
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

1.2. Les grandeurs locales

MSRV_AXIAL

(E-3) Tesla

400

300

CURVE FLUX
Flux density / Magnitude
Time
200 Point(0,0)

100

(E-3) s.

10 20 30 40 50 60

Fig.h3- Le flux magnétique

AXIAL_CCMONOPH25 AXIAL_CCMONOPH25
Tesla
A/m 199,999E-3

0,1
40 CURVE C2D_7
Flux density / Normal component
0 airgap
Time (s.) : 66,599999E-3

-0,1
35

-0,2 mm
CURVE flux
0 100 200 300 400 500 600 699,999
30 Field/ Magnitude
Time Tesla
Point(0,0)
25
0,05

CURVE tangente
Flux density / Tangent component
20 0 airgap
Time (s.) : 66,599999E-3

(E-3) s.
-49,999E-3
mm
100 150 200 250
0 100 200 300 400 500 600 699,999

Fig. h4- le flux magnétique Fig. h5 - la composante normale et

tangentielle

1.3. Interprétation des résultats

Le courant ia augmente de 14 [A] à 15 [A] l’ors un défaut de 25 % si on compare l’influence
d’un court circuit monophasé de même dégrée d’une msrv et mrv on vois d’après les calcules
statistique précédente le court circuit augmente la valeur du couple de la mrv 3 fois celle de la

120
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

msrv et pour le couple de la msrv on vois une modulation et augmentation de valeur aller à 6
[Nm], le couple de la msrv par contre la mrv est ondulé donc le défaut apparaît sur le long du
couple comme la machine asynchrone ou synchrone, dans la mrv le couple est constitue de
portions le défauts claire sur la portion infecté.

121
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

2. Un court circuit triphasé de 25 %

2.1. Les grandeurs globales

AXIAL_CCTRIPHASE25

CURVE C2D_2
Circuit / Current
Time
PA ;
10

0
CURVE ib
Circuit / Current
Time
PB;

-10

CURVE ic
Circuit / Current
Time
MC;
0,1 0,15 0,2 0,25

Fig. i1- Les courants dans les trois phases

AXIAL_CCTRIPHASE25

N.m

2,5

CURVE couple
Mechanics / Axis torque
Time
0

-2,5

(E-3) s.

100 150 200 250

Fig.i2 - Le couple moteur

122
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

2.2. Les grandeurs locales

AXIAL_CCTRIPHASE25 AXIAL_CCTRIPHASE25
(E-3) Tesla
(E-3) Tesla 200

400 100
CURVE normale
Flux density / Normal component
0 airgap
Time (s.) : 66,599999E-3
350
-100

-200 mm
300 CURVE flux
Flux density / Magnitude 0 100 200 300 400 500 600 699,999
Time
Point(0,0) (E-3) Tesla

250
50

CURVE tangente
200 Flux density / Tangent component
0 airgap
Time (s.) : 66,599999E-3

(E-3) s.
-50
mm
100 150 200 250
0 100 200 300 400 500 600 699,999

Fig.i3 - Le flux magnétique Fig.i4 - La composante normale et

tangentielle

2.3. Interprétation des résultats

Les courant de trois phases fait une augmentation à 15 [A] par apport le cas sain et le couple
plus il augmente est modulé, l’influence du court circuit dans une machine à réluctance
variable est plus fort qu’une machine synchrone à réluctance variable.

123
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

3. Un circuit ouvert dans la phase a

Toujours avec la même procédure appliqué sur la MRV on augmente la résistance de la phase
a pour empêcher la circulation du courant dans cette phase un défaut de circuit ouvert.

3.1. Les grandeurs globales

__
AXIAL_OPENCIRCUIT

Ampere

10

5
CURVE ia
Circuit / Current
Time
0 PA ;

-5

(E-3) s.

100 150 200 250

Fig.j1- Le courant de la phase a

AXIAL_OPENCIRCUIT

CURVE ia
Circuit / Current
Time
PA ;
10

CURVE ib
Circuit / Current
0
Time
PB ;

-10
CURVE ic
Circuit / Current
Time
MC ;
0,1 0,15 0,2 0,25

Fig.j2- Les courant des phases ia,ib,ic

124
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

AXIAL_OPENCIRCUIT

N.m
4

CURVE couple
1
Mechanics / Axis torque
Time

-1

(E-3) s.

100 150 200 250

Fig.j3 - Le couple moteur

3.2. Les grandeurs locales

AXIAL_OPENCIRCUIT AXIAL_OPENCIRCUIT
Tesla
(E-3) Tesla 199,999E-3
300
0,1
CURVE normale
Flux density / Normal component
0 airgap
Time (s.) : 66,599999E-3
250
-0,1

-0,2 mm
200 CURVE flux
Flux density / Magnitude 0 100 200 300 400 500 600 699,999
Time
Point(0,0) Tesla

150
0,05

CURVE tangente
Flux density / Tangent component
100 0 airgap
Time (s.) : 66,599999E-3

(E-3) s.
-49,999E-3
mm
100 150 200 250
0 100 200 300 400 500 600 699,999

Fig.j4- Le flux magnétique Fig.j5- La composante normale et

tangentielle

3.3. Interprétation des résultats

Le courant de la phase a ouvert diminue par contre les courants des phases b et c augmente et
le couple comme présenté dans la figure.j3 apparaît une déformation de forme et une
augmentation de valeur aller à 4.35 [Nm].

125
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

4. Les défauts au niveau de l’onduleur

Dans cette section on présente les principaux défauts au niveau de l’onduleur et puis on les
commentes.

4.1. Un IGBT circuit ouvert

IGBT CO

[Link].19 – Onduleur avec un interrupteur circuit ouvert

4.1.1. Les grandeurs globales

AXIAL_IGBTCO

CURVE IA
Circuit / Current
10 Time
PA ;

0
CURVE IB
Circuit / Current
Time
-5 PB ;

-10
CURVE IC
Circuit / Current
-15
Time
MC ;
0,1 0,15 0,2 0,25

Fig.k1- Les courants de phases

126
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

AXIAL_IGBTCO

N.m

2,5

0
CURVE couple
Mechanics / Axis torque
Time

-2,5

(E-3) s.

100 150 200 250

Fig.k2 – Le couple moteur

4.1.2. Les grandeurs locales

AXIAL_IGBTCO AXIAL_IGBTCO
(E-3) Tesla
(E-3) Tesla
200
350
100
CURVENORMAL
Flux density / Normal component
300 0 AIRGAP
Time(s.) :66,599999E-3

-100

250
-200 mm
CURVE C2D_7
Flux density / Magnitude
Time 0 100 200 300 400 500 600 699,999
200 Point(0,0)
(E-3) Tesla

150 50

CURVETANGENTE
Flux density / Tangent component
100 0 AIRGAP
Time(s.) : 66,599999E-3
(E-3) s.

100 150 200 250 -50

mm
0 100 200 300 400 500 600 699,999

Fig.k3 – Le flux magnétique Fig.k4 – La composante normale et

tangentielle

127
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

4.2. Un IGBT court circuit

IGBT CC

[Link].20 – Onduleur avec un interrupteur court circuit

4.2.1. Les grandeurs globales

__
AXIAL_IGBT CC

Ampere

-10

-20
CURVE IA
Circuit / Current
Time
-30
PA ;

-40

(E-3) s.

100 150 200 250

Fig.l1 – Le courant dans la phase a

128
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

_
_ AXIAL_IGBT CC

CURVE IA
Circuit / Current
Time
25 PA ;

0
CURVE IB
Circuit / Current
Time
PB ;
-25

CURVE IC
Circuit / Current
Time
MC ;
0,1 0,15 0,2 0,25

Fig.l2 – Les courants dans les trois phases

__
AXIAL_IGBT CC

N.m

CURVE C2D_1
-10 Mechanics / Axis torq
Time

-20

(E-3) s.

100 150 200 250

Fig.l3 – Le couple moteur

129
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

4.2.2. Les grandeurs locales

AXIAL_INVERTER AXIAL_INVERTER
(E-3) Tesla
(E-3) Tesla 200
500
100
CURVE NOR
Flux density / Normal component
0 airgap
400 Time (s.) : 66,599999E-3

-100

300 -200 mm
CURVE F
Flux density / Magnitude
Time 0 100 200 300 400 500 600 699,999
Point(0,0)
(E-3) Tesla
200

50

100 CURVE TANGE

Flux density / Tangent component
0 airgap
Time (s.) : 66,599999E-3
(E-3) s.

100 150 200 250 -50

mm
0 100 200 300 400 500 600 699,999

Fig.l4 – Le flux magnétique Fig.l5 – La composante normale et

tangentielle

4.3. Interprétation des résultats

On commence par le cas d’un onduleur avec un interrupteur ouvert [Link].19 , ce cas est
semblable au cas d’une phase ouverte on vois la même allure du couple sauf l’augmentation
de la valeur du couple est moins importante par apport un circuit ouvert et ce défaut peut être
simulé en arrêtant des signaux de porte à cette phase à l'instant de défaut, le courant dans la
phase a correspondante est diminue à 2 [A] l’interrupteur ouvert arrête le courant d’excitation
à cette phase les courant ib et ic diminues de 28 % par apport un cas sain.

Dans le cas ou un interrupteur de l’onduleur est court circuit [Link].20 ,le courant dans la
phase infecté est diminue par contre le courant dans la phase b aller à 42.5 [A], la phase a
cesse à produire le couple, la valeur du couple total augmente à 7.75 [Nm] avec une
déformation de ca forme par apport le cas sain dont la forme est ondulée.

130
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.5. Les défauts mécaniques [10] [11] [12]

Ils existe diffèrent types de balourd mécanique ou le mauvais positionnement du rotor par
rapport au stator dans cette étude on prend le cas de l’excentricité dynamique.

Nous allons modéliser avec le logiciel éléments finis FEMM la machine synchrone à
réluctance variable et la machine à réluctance variable présentant une excentricité dynamique
(figure III.23 et III.24 )

[Link].22 - Le positionnement du centre rotorique « Or »

Pour réaliser ce défaut il faut décrire le déplacement du rotor par des coefficients a et b

g (θ ) = g − a cos(θ ) − b sin(θ ) (III.21)

0

a = g 0 γ d cosθ r (III.22)

b = g 0 γ d sni θ r

g 0 = Rs − Rr (III.23)

Et les dégrées statique et dynamiques respectivement

γ d +γ s <1 (III.24)

131
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

La procédure de réalisation est suffit de positionner le rotor par apport le stator et mesure le
couple statique sur phase à chaque position si on veut voire le comportement dynamique de ce
défauts on peut construire des blocs simulink et avec les méthodes d’interpolation kreigeage
ou bi cubic spline on peut voire le comportement dynamique dans notre étude on ce limite a
voire l’influence sur le comportement statique du couple.

couple=f(teta;ni) couple=f(teta;ni)
9 16

8 14

7 12

6
10

couple costa [Nm]

couple costa [Nm]

5
8
4
6
3
4
2
2
1
0
0

-1 -2
0 5 10 15 20 25 30 35 40 45 0 5 10 15 20 25 30 35 40 45
teta [degres] teta [degres]

(a) (b)

[Link].23 - (a) Le couple statique de la MRV dans la cas sain et (b) Le couple statique en cas
d’une excentricité dynamique de 25%.

couple=f(teta;ni) couple=f(teta;ni)
1200 900

800
1000
700

800 600
couple costa [Nm]
couple [Nm]

500

600
400

300
400

200

200
100

0
0 0 5 10 15 20 25 30 35 40 45
0 10 20 30 40 50 60 70 80 90 teta [degres]
teta [degres]

(a) (b)

[Link].24 - (a) Le couple statique de la MSRV dans la cas sain et (b) Le couple statique en
cas d’une excentricité dynamique de 25%.

132
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.6. L’analyse spectrale

On vas montrer les spectres du couple et courant statorique dans le cas de différents défauts
simuler par flux 2d pour par la mrv les spectre de la msrv donnés dans l’annexe.

III.6.I. LA MRV
III.6.1.1. L’analyse spectrale du couple
Pour voire le comportement fréquentielle pour chaque défaut il faut une analyse minutieuse
des résultats :

MRV_STATRTSAINE MRV_VOURTC20TEST
N.m N.m
5

4
4
CURVE COUP CURVE COUP
Mechanics / Axis torque Mechanics / Axis torque
3 Time Time
3
Full cycle / Normal Full cycle / Normal

2 2
(E-3) s. (E-3) s.
15 20 25 30 35 40 15 20 25 30 35 40

(E-3) N.m (E-3) N.m

400 500

400
300
SPECTRUM Spectr_45 300
SPECTRUM Spectr_29
From COUP From COUP
Fundamental 35,713 Fundamental 35,713
200 200

100
100

1 2 3 4 5 6 7 1 2 3 4 5 6 7

Fig.n1- Couple sain Fig.n2- Couple avec cc monophasé

_
_
_
MRV_CC03PHASE MRV_OPENCIRCIU
N.m N.m T
5
4
4
CURVE CC
Mechanics / Axis torque
Time
3 CURVE couple
Mechanics / Axis torque
3 Time
Full cycle / Normal
2 Full cycle / Normal
2
(E-3) s. 1
(E-3) s.
15 20 25 30 35
15 20 25 30 35
(E-3) N.m
550
N.m
1,5
500

SPECTRUMSpectr_7
From CC
450 Fundamental 37,036 1 SPECTRUM Spectr_17
From couple
Fundamental 37,036

400
0,5

1 2 3 4 5 6 7

1 2 3 4 5 6 7

Fig.n3 - Couple en cas cc triphasé Fig.n4 - Couple circuit ouvert

133
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

MRV_INVERTERFAULT02 MRV_INVERTERFAUT
N.m N.m
6

4
5

3 CURVE couple CURVE couple

Mechanics / Axis torque
4 Mechanics / Axis torque
Time Time
2
Full cycle / Normal 3 Full cycle / Normal

1
2
(E-3) s. (E-3) s.

15 20 25 30 35 15 20 25 30 35

N.m N.m
1

1,5
0,75

SPECTRUMSpectr_1 SPECTRUMSpectr_3
1 From couple From couple
0,5 Fundamental 37,036
Fundamental 37,036

0,5 0,25

1 2 3 4 5 6 7
1 2 3 4 5 6 7

Fig.n5- Couple en cas d’un interrupteur cc Fig.n6- Couple en cas ou une phase
de l’onduleur cc

On vois qu’un court circuit monophasé de 25% fait une amplification de l’harmonique 3 et
4 ;et pour un court circuit triphasé 25% on vois l’apparition d’une nouvelle harmonique
l’harmonique 2 et amplification de l’harmonique 3 ;

Au niveau de l’onduleur l’amplification de l’harmonique 4 trois fois par apport un cas sain et
l’apparition de l’harmonique 2 indique qu’un interrupteur est court circuit et l’amplification
de l’harmonique 5 trois fois et l’harmonique 6 plus de deux fois par apport un cas sain
indique d’une phase complète de l’onduleur est court circuit.

Le cas d’un circuit ouvert est semblable au cas ou une phase de l’onduleur court circuit.

134
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

III.6.1.2. L’analyse spectrale du Courant statorique ia

MRV_STATRTSAINE MRV_VOURTC20TEST
Ampere Ampere
30
30

20 CURVE ibsain 20 CURVE isain

Circuit / Current Circuit / Current
Time Time
BOBB ; BOBA ;
Full cycle / Normal Full cycle / Normal
10 10

0 (E-3) s. 0 (E-3) s.
15 20 25 30 35 15 20 25 30 35

Ampere Ampere
10

10
7,5
SPECTRUM Spectr_43 SPECTRUM Spectr_23
From ibsain From isain
Fundamental 37,036 Fundamental 37,036
5 5

2,5

1 2 3 4 5 6 7 1 2 3 4 5 6 7

Fig.O1- Le courant sain Fig.O2-Le courant avec cc monophasé

MRV_CC03PHASE MRV_INVERTERFAULT02
Ampere Ampere
15
30

CURVE isain 10
20 Circuit / Current
CURVE isain
Circuit / Current
Time Time
BOBA ; BOBA ;
Full cycle / Normal Full cycle / Normal
10 5

0 (E-3) s. (E-3) s.
0
15 20 25 30 35 15 20 25 30 35
Ampere Ampere
3

10
2,5
SPECTRUMSpectr_9 SPECTRUMSpectr_3
From isain
Fundamental 37,036 2 From isain
5 Fundamental 37,036

1,5

1
1 2 3 4 5 6 7
1 2 3 4 5 6 7

Fig.O3- Le courant avec cc triphasé Fig.O4- Le courant avec circuit ouvert

MRV_INVERTERFAUT MRV_OPENCIRCIUT
Ampere Ampere
15
50

40 10
CURVE isain CURVE isain
Circuit / Current Circuit / Current
Time Time
30 BOBA ; BOBA ;
Full cycle / Normal Full cycle / Normal
5

20

(E-3) s. 0 (E-3) s.

15 20 25 30 35 15 20 25 30 35

Ampere Ampere
5
10
4

7,5 SPECTRUMSpectr_5 3 SPECTRUM Spectr_19

From isain From isain
Fundamental 37,036
Fundamental 37,036
2
5
1

2,5
1 2 3 4 5 6 7 1 2 3 4 5 6 7

Fig.O5- Le courant en cas ou Fig.O6- Le courant en cas ou

un interrupteur de l’onduleur cc une phase de l’onduleur cc

135
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

Au niveau des courant statoriques on va analysé seulement le courant de la phase a le cas d’un
court circuit monophasé de 25% Fig.O2 fait apparaître une nouvelle harmonique
l’harmonique 2 et amplification de l’harmonique 5 quartes fois le cas sain Fig.O1;
Le cas d’un court circuit triphasé Fig.O3 de 25% par apport un court circuit monophasé ne
donne pas de nouvelles harmoniques juste une amplification de valeurs biens sure dans ce cas
au niveau de trois courant statoriques.
Le cas d’un circuit ouvert Fig.O4 la disparition de l’harmonique 2 et 3 et la diminution de
l’harmonique 4 8 fois l’apparition de l’harmonique 6 indique qu’il y a un circuit ouvert à la
phase a.
Au niveau de l’onduleur un interrupteur court circuit Fig.O5 traduit par l’amplification de
l’harmonique 5et 4 ;Si une phase de l’onduleur court circuit Fig.O6 la diminution de
l’harmonique 4 8 fois le cas sain et l’apparition de l’harmonique 6 un bon indicateur de
défaut.

136
Chapitre III ANALYSE DE DEFAUTS PAR ELEMENTS FINIS

Conclusion

Nous avons traité dans ce chapitre la modélisation et la simulation des machines MRV et
MSRV par la méthode des éléments finis à l’aide de deux logiciel le femm et le flux 2d. La
modélisation a été réalisée en deux régimes, à savoir, le régime statique et le régime
dynamique. La simulation est effectuée dans le cas sain et dans le cas de défaut.

Ce chapitre comporte trois parties essentielles : le régime statique où nous avons utilisés le
femm pour voir les caractéristiques statiques de chaque machine. C’est grâce à ce logiciel
que nous avons simulés le défaut d’excentricité dynamique à diffèrent degré. La second partie
de ce chapitre porte sur le régime dynamique, celui-ci a été réalisé par le logiciel flux 2d .
Dans cette partie les deux machines sont modélisées et couplées avec un onduleur de tension
dont les interrupteurs sont commandés par une commande décalés. Nous avons appliqué
différents types de défauts soit au niveau du stator tel que le court circuit , circuit ouvert.
Soit au niveau du circuit d’alimentation « l’onduleur » tel que interrupteur court circuité,
circuit ouvert et une phase court circuitée.

L’analyse fréquentielle réalisée par la FFT a été appliquée sur le courant statorique dans le
cas de défaut et sur le couple moteur de la machine à vide et à vitesse constante.

Le dernier chapitre est consacre à la modélisation analytique et à la simulation de la machine

à réluctance variable en régime dynamique.

137
 Références du Chapitre III

Références du Chapitre III

[1] Zaatar. MAKNI, «Contribution au Développement d’un Outil d’Analyse Multiphysique

pour la Conception et l’Optimisation d’Actionneurs Electromagnétiques», Thèse Doctorat,
Univ Paris-SUD XI , d’Orsay, 12 Décembre 2006.

[2] Redoune. Maaziz, «Développement et implantation d’un modèle prenant en compte

les aspects scalaires vectoriels et dynamiques du modèle d’hystérésis», Thèse Magister, Univ
de Batna, 17 Décembre 2006.

[3] Malek. Bouharkate, «Etude et l’évolution des courants rotoriques d’une machine
asynchrone à cage en regime dynamique», Thèse Doctorat, Univ Batna, 5 Février 2006.

[4]D. MEEKER, « Finite Elements Method Magnetics », User’s Manual, v3.4, Feb.
2005. [Link] [Link].

[5] Ghislain. Remy, «Commande optimisée d’un actionneur lineaire synchrone pour un axe
de positionnement rapide», Thèse Doctorat, Univ de Paris, L’institut des arts et mitiers, 07
Décembre 2007.

[6] Flux-2D (version 7.6) commercialisée par Cedrat, Zirst, 38 240 Meylan, France.

[7] Guillaume. Lacomb, «Définition et realisation d’une nouvelle génération de logiciel pour
la conseption des moteur au futur», Thèse Doctorat, Univ de Paris, L’INP Gronoble, 29
Novembre 2007.

138
 Références du Chapitre III
[8] J. Nascimento, L. Rolim, P. Heidrich et al, «Design and simulation aspects of a switched
reluctance drive», Proc. of 4º Brazilian Power Electronics Conference (COBEP96) ”, pp. 79-
83,São Paulo, 1996.

[9] Tsarafidy . Raminosoa « Optimisation des performances des machines synchro-reluctantes

par reseaux de permeances », Thèse Doctorat, L’INPL de L’orraine,5 Octobre 2006.

[10] H. Torkamanand, [Link], «Magneto static field regarding analysis regarding the effect
of eccentricity in switched reluctance motor», Shahid Beheshti University, Tehran Iran,
Vol.8,163180,2009.

[11] Juan R. Briso-Montiano, Roland. Karrelmeyer and Elmar. Dilger Robert. Bosch GmbH.
Stuttgart, «Simulation of Faults by means of Finite Element Analysis in a Switched
Reluctance Motor», Conference 2005.

[12] Martin. BLODT, «Condition Monitoring of Mechanical Faults in Variable Speed

Induction Motor Drives», Thèse Doctorat, L’INPL de Toulouse,14 Septembre 2006.

[13] Iqbal. Husain, M.N. Anwa, « Fault Analysis of Switched Reluctance Motor Drives »
Article, [Link] Electrical Engineering, University of Akron.

139
CHAPITRE : IV

MODÉLISATION, COMMANDE
&
SIMULATION DE DÉFAUTS
PAR
MATLAB / SIMULINK

140
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

INTRODUCTION

Nous avons déjà évoqué certains outils numériques pour la résolution des équations de
l’électromagnétisme ; leur principal intérêt est qu’ils nécessitent de faire peu hypothèses
Simplificatrices. mais d’une part les coûts calculatoires restent prohibitifs par rapport aux,
objectifs du modèle (contrainte de temps) et d’autres part ce type d’outil demande un
investissement pour leur prise en main qu’ils réserve a des utilisateurs « confirmés ».

Apres avoir interprété les résultats numériques ce chapitre décrit analytiquement la simulation
des défauts en régime linéaire et non linéaire de la machine à réluctance variable, à l’aide du
logiciel MATLAB/SIMULINK : une simulation dynamique saine et puis avec défauts a été
réalisée avec deux stratégies de commande en hystérésis et en tension.

Les logiciels d’utilisation plus générale pour le calcul matriciel (MATALB-Simulink,

ACSL ,MATRIXx….) leur avantage est d’être beaucoup plus ouvert par contre leur utilisation
nécessite généralement de trouver manuellement les équations différentielles régissant le
système. Pour des systèmes complexes cette recherche peut être laborieuse.

IV.1. Présentation de l'environnement MATLAB/SIMULINK™

MATLAB est un logiciel destiné principalement au calcul scientifique, à la modélisation et à
la simulation. Le noyau de calcul est associé à l'environnement de SMULINK permettant une
modélisation basée sur des schémas-blocs. Des bibliothèques spécialisées sont disponibles
(les "Toolboxes") pour la plupart des domaines scientifiques nécessitant des moyens de calcul
importants : automatique, traitement de signal, mathématiques appliquées,
télécommunications, etc. Des modules temps réel, développés autour des produits d’SPACE
(Real Time Interface, TRACE, COCKPIT) sont également proposés et rendent
l'environnement de travail particulièrement attractif : l'utilisateur dispose en effet de tous les
outils nécessaires à l'élaboration d'une application sur procédé réel, de la simulation à la
supervision.

141
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.2. Modèle de simulation

Pour simuler la machine a réluctance variable avec sa commande à hysteresis et en tension, le
logiciel Matlab/Simulink connue utilise par de nombreuses communautés nous assemble un
bon compromis.
La simulation de la machine est alors fondée sur L’intégration numérique d’une
représentation d’état continue de la machine qui nécessite :

[Link] trouver manuellement l’ensemble des équations différentielles indépendantes régissant

le volution de la machine et de son alimentation ;

[Link] déterminer les valeurs numériques des paramètres du modèle dans notre cas particulier il
s’agit des inductances .

Dans les sections suivantes nous présentons les méthodes et les simulations réalisées sous
MATLAB versions 6.5 et 7.

IV.3. Les équations électromagnétiques [1]

Le modèle de la MRV est élaboré sur la base d’un certain nombre d’hypothèses
simplificatrices :
♦On néglige les effets des mutuelles et de la saturation (l’inductance de la machine
indépendante du courant est idéalisée).
♦Les interrupteurs du convertisseur sont parfaits.
♦La tension est parfaitement continue.
♦ On suppose que la résistance ne dépend pas de la température.

La tension instantanée appliquée sur l'enroulement d'une phase du moteur MRV à une relation
avec le flux par la loi de Faraday:
dψ
VS = Ri + (IV-1)
dt
VS: est la tension appliquée.
i: est le courant dans la phase.
R: est la résistance de l'enroulement de phase.
Ψ: est le flux.

142
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

À cause de la structure à double saillance du moteur MRV et l'effet de la saturation

magnétique, le flux dans une phase MRV varie en fonction de la position du rotor θ et le
courant de la phase. L'équation (II-1) peut être développée comme suit:

dψ di dψ dθ
V = Ri + + (IV-2)
di dt dθ dt
∂Ψ ∂Ψ ∂θ
est défini comme l'inductance instantanée L(θ,I) et le terme est le FCEM
∂i ∂θ ∂t
instantané.

IV.4. Production du couple

La MRV peut être décrite par une fonction convexe qui dépend seulement de la position du
rotor θ et le courant dans le neme phase. Cette fonction est la co-énergie W (ψ,θ). Dans une
manière semblable, la fonction de l'énergie W (ψ, θ), où les variables sont les flux de n phase
et la position de rotor, aussi permettent de décrire la MRV. Tenant compte des vecteurs ψ et I,
la fonction de co-énergie et énergie vérifie l'inégalité suivante:

W ( I, θ ) + W ( Ψ , θ ) ≥ Ψ t I (IV-3)

Grâce à son double saillance, la MRV peut avoir une variation de l'énergie magnétique et par
conséquence la production du couple. La dérivée partielle de la fonction énergie par rapport à
la position du rotor donne le couple de la machine :
∂W
Ce ( Ψ1 ,..., Ψ n , θ ) = ( Ψ1,..., Ψ n , θ ) (IV-4)
∂θ

Appliquant cette relation à une MRV 6/4, on a :

∂W
Ce ( Ψ1 , Ψ 2 , Ψ 3 , θ ) = ( Ψ1, Ψ 2 , Ψ 3 , θ ) (IV-5)
∂θ
Quand on excite une phase, le couple apparaît afin que le rotor évolue dans la direction des
augmentations de l'inductance. Par conséquent, le couple sera dans la direction de la position
alignée la plus proche.

143
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.5. Modèle Linéaire de la MRV [1] [3]

Les caractéristiques électromagnétiques utilisées pour la simulation de la MRV 6/4 sont
obtenus par éléments finis.
L'équation électrique de chaque phase est donnée par:

dψi (θ,ii )
Vi = Rii + , i = {1, 2,3} (IV-6)
dt

Le flux dans chaque phase est donné par l'équation linéaire:

Ψ i (θ,ii ) = L(θ)ii (IV-7)

L'énergie totale associée avec les trois phases est donnée par :

1 3 (IV-8)
Wtotal = ∑ L(θ + (n − i − 1)θs )ii2 .
2 i =1

Et le couple total du moteur est donné par :

1 3 L(θ + (n − i − 1)θs ) 2
Ce = ∑
2 i =1 dθ
ii (IV-9)

Les équations mécaniques sont :

dω (IV-10)
J = Ce − Cr − Fω .
dt

dθ
=ω (IV-11)
dt

Cr : représente le couple de charge,

F : Le coefficient du frottement de la machine.
J : Moment d’inertie.

La figure (IV.3) représente le diagramme général de simulation utilisé pour le modèle linéaire
de la MRV alors que la figure (IV.4) montre le contenu du bloc de la phase 1. Il contient trois
autres blocs, chacun est associé à une fonction MATLAB spécifique. Les paramètres de la
machine sont donnés en annexe C.

144
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

Comme ce fut le cas en numérique avec flux 2d, il est nécessaire de choisir une position
initiale du rotor, l’angle choisi doit être supérieure à 15° car si l’angle au départ inférieur a
15° le couple de charge sera nul et par conséquence le rotor n'évolue pas. Quand le couple de
la charge n'est pas nul, le rotor ce déplace pour établir une vitesse où Cr = Ce .
Noter que la position initiale choisit selon le profil de l’inductance.

IV.5.1. La stratégie de commande par hysteresis [1] [3]

La stratégie de commande utilisée pour le régime linéaire est la commande en courant par
hystérésis. Cette technique est valable pour une plage de fonctionnement à couple constant et
à basses vitesses, donc la MRV est alimentée par des créneaux de courant de hauteur
réglable en fonction du couple demandé. Ces créneaux peuvent êtres simplement de forme
rectangulaire appliqués durant la croissance de l'inductance si l'on veut un fonctionnement en
mode moteur et durant la phase décroissante de l'inductance si l'on veut un fonctionnement en
mode générateur. Pour mettre en œuvre le système de régulation des courants, plusieurs
solutions sont envisageables.

Cette méthode de contrôle des courants, la plus utilisée dans l’industrie actuellement, consiste
à commander les interrupteurs de l’onduleur de telle sorte que les courants de phase statorique

du moteur ne puissent évoluer en dehors d’une bande de largeur ∆i encadrant leur courbe de
référence.

Ce contrôle se fait donc par une comparaison permanente entre les courants réels et leurs
références. La différence entre les deux valeurs permet à la logique de commande de
l’onduleur d’imposer une commutation des interrupteurs à chaque fois que cette différence

∆i
dépasse, en valeur absolue, .Les résultats montrés ont été obtenus pour un courant de
2
référence Iref= 8[A], avec fonctionne sans charge (à vide). θon =0° et θoff ajusté à la valeur
de l'angle 38°, pour diminuer les oscillations de vitesse il est nécessaire de produire plus de
couple.

145
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

[Link].1 - Le principe de commande en courant par hysteresis

Forme de tension et courant.

0 .0 6

0 .0 5

0 .0 4

0 .0 3

0 .0 2

0 .0 1

0
0 .0 2 0 .0 4 0 .0 6 0 .0 8 0 .1 0 .1 2 0 .1 4

1 5 0

1 0 0

5 0

-5 0

-1 0 0

-1 5 0
0 .0 2 0 .0 4 0 .0 6 0 .0 8 0 .1 0 .1 2 0 .1 4

0
0 .0 2 0 .0 4 0 .0 6 0 .0 8 0 .1 0 .1 2 0 .1 4

[Link].2 - Le principe de commande en courant par hysteresis

Forme de tension ,courant et de l’inductance par Simulink.

146
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.5.2. Schémas bloc de simulation

t
Clock
To Workspace4
iabc

torque
Iref
ind1
Step
current1
T1
teta1
tension torque
total torque
Phase1 0

torque2 vitesse [rd/s]

Iref 1 1 0
1/J
ind2 s s Te
Sum
Sum1 speed 1/J teta teta1
current2 To Workspace
teta2
tension v b L
F
Phase2 L,VA,I To Workspace1
f
all i
torque3
Iref To Workspace2
ind3 L, VA, I
vabc
w
current3
teta To Workspace3
vc

Phase3
L, VA, I
0

ind2
0

ind3

Fig. IV.3_ Diagramme de simulation de MRV en regime lineaire.

R 3
current1
R i1
To Workspace2
va 1
1
f lux i i
v s
Iref Sum voltage torq 1
teta Sum1 a
ind teta torque
switch1 t1
inductance torqu
To Workspace3

MATLAB
Function
modulo pi/2

flux

4
0 tension

teta s V f
2
2 To Workspace1 To Workspace5 ind1
teta1 Sum2

Fig. IV.4_ Diagramme de simulation de MRV par phase en régime linéaire.

147
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.5.3. Simulation du régime linéaire sain

6
le courant [A]

0
0.2 0.22 0.24 0.26 0.28 0.3 0.32
le temps[s]

Fig.a1 - Le courant de la phase a

150

100

50
la tension va [V]

-50

-100

-150
0.2 0.22 0.24 0.26 0.28 0.3 0.32
le temps [s]

Fig.a2 - La tension de la phase a

148
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

3.5

2.5

2
le couple [Nm]

1.5

0.5

-0.5 couple negatif

-1
0.2 0.22 0.24 0.26 0.28 0.3 0.32
le temps [s]

Fig.a3 - Le couple de la phase a

3.5

2.5
le couple total [Nm]

1.5

0.5

0
0.2 0.22 0.24 0.26 0.28 0.3 0.32
le temps [s]

Fig.a4 - Le couple électromagnétique total

149
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

160

140

120

100
la vitesse [rd/s]

80

60

40

20

-20
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
le temps [s]

Fig.a5 - La vitesse

0.06

0.05

0.04
l'inductance [H]

0.03

0.02

0.01

0
0.2 0.22 0.24 0.26 0.28 0.3 0.32
le temps [s]

Fig.a6 - L’inductance de la phase a

150
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.5.4. Simulation de défauts

IV.5.4.1. Le cas d’un court circuit monophasé 25% dans la phase C à t= 0.167 [s]

140

120
Les fluctuations

100
la vitesse [rd/s]

80
avant cc apres le court circuit
60

40

20
L'INSTANT DE COURT CIRCUIT
0

-20
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s]

zoom

138

137

136

135

134

0.2 0.25 0.3 0.35 0.4 0.45 0.5

Fig.b1 - La vitesse

151
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

4.5

le couple total electromagnétique [Nm]

La phase court circité
3.5

2.5

1.5

0.5

0
0.12 0.14 0.16 0.18 0.2 0.22 0.24
le temps [s]

Fig.b2 - Le couple total

5
couple de la phase a
couple de la phase b
4 couple de la phase c
les couples de phase [Nm]

-1
0.12 0.14 0.16 0.18 0.2 0.22 0.24
le temps [s]

Fig.b3 - Le couple dans les phases a , b et c

152
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

12

10
le courant de la phase c ic[A]

0
0.12 0.14 0.16 0.18 0.2 0.22 0.24
le temps [s]

Fig.b4 - Le courant dans la phases c

12
ic
ib
10 ia
les courant ia ib et ic [A]

0
0.12 0.14 0.16 0.18 0.2 0.22 0.24
LE TEMPS [S]

Fig.b5 - Le courant dans les phases a , b et c

153
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.5.4.2. Le cas d’un court circuit triphasé 25% à t= 0. 3269 [s]

160

140

120

100
zoom
la vitesse [rd/s]

80
145

140
60
135

130
40
0.34 0.36 0.38 0.4 0.42 0.44 0.46

20
L'instant de défaut =0.3269
0

-20
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
le temps [s]

Fig.c1 - La vitesse

8
le couple éléctromagnétique total [Nm]

0
0.3 0.31 0.32 0.33 0.34 0.35 0.36 0.37
le temps [s]

Fig.c2 - Le couple total

154
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

15
ia
ib
ic

les courants ia,ib,ic [A]

10

0
0.3 0.31 0.32 0.33 0.34 0.35 0.36 0.37
le temps [s]

Fig.c3 - Les courants de phases

9
couple phase c
8 couple phase b
couple phase a
7
les couples phases a,b,c [Nm]

-1
0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s]

Fig.c4 - Le couple de phases a,b,c

155
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.5.4.3. Le cas d’un circuit ouvert

160

140

120
circuit ouvert phase a
avant défaut
100 circuit ouvert phase b et a
la vitesse [rd/s]

80

60

40

20

-20
0 0.5 1 1.5
le temps [s]

Fig.d1- La vitesse

3.5
phase a ouvert phase a et b ouvert
3

2.5
le couple total [Nm]

1.5

0.5

0
0.86 0.87 0.88 0.89 0.9 0.91 0.92 0.93 0.94 0.95 0.96
le temps [s]

Fig.d2- Le couple électromagnétique total

156
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

couple phase c
3.5 couple phase b
couple phase a
3
le couple de phase Ca,Cb,Cc [Nm]
2.5

1.5

0.5

-0.5

-1
0.86 0.87 0.88 0.89 0.9 0.91 0.92 0.93 0.94 0.95 0.96
le temps [s]

Fig.d3 - Le couple de phases

ia
9 ib
phase a ouverte phase a et b ouvertes ic
8

7
les courants ia,ib,ic[A]

0
0.86 0.87 0.88 0.89 0.9 0.91 0.92 0.93 0.94 0.95 0.96
le temps [s]

Fig.d4 - Les courants de phases

157
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.5.5. Commentaire

L’application d’un défaut cc monophasé de 25 % induit une augmentation « pic » égale à 17

fois la valeur saine à l’instant de défaut dans la phase infectée les autres phase reste saine et
après l’instant de défaut la phase infecté cesse à produire le couple on voie dans fig. b2 la
valeur du couple diminue cette diminution dépend du court circuit . Pour le courant fig.b4 sa
valeur augmente à l’instant de défaut et comme nous avons commandé la machine en courant
sa valeur reste sans changement après l’instant de défaut. La vitesse fig.b1 présente des
fluctuation et une diminution légère.

Le cas d’un cc triphasé fait diminue la vitesse 10% les couples de phases diminues de 3[Nm]
à 2.75[Nm] et le courant présente des modulations.

Le cas d’un circuit ouvert la vitesse fig.d1 diminue, le courant fig.d4 présente une
déformation de la forme , nous avons simuler le cas ou une phase ouverte et à un certain
temps deux phases ouverts pour mieux voir l’influence de défaut, le couple moteur fig.d2 et
les couples de phase fig.d3 diminues ce défaut au niveau du couple présente une modulation
dans la phase infecté et diminution de valeur.

158
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.6. Régime non linéaire [1] [4]

Le caractère non–linéaire de la caractéristique magnétique B = f(H) est maintenant considéré.

La meilleure solution qui permet d’aboutir à une bonne simulation d'une machine à
réluctance variable réside dans l’utilisation d’une méthodologie qui permet de prendre en
compte la non–linéarité de sa caractéristique magnétique en minimisant le temps de
simulation. La procédure que nous avons adoptée avec MATLAB–SIMULINK consiste
d'éviter toutes les dérivées partielles qui constituent des sources d'erreurs. La technique utilise
un tableau « look–up table » qui approxime des relations flux/courant/position et courant
/couple /position, comme montré dans la figure (IV.5).

Dans le chapitre III, notre machine MRV 6/4 a été modélisée en magnétostatique par une
analyse à base d’éléments finis avec logiciel FEMM qui nous a fournie ses différentes
données magnétiques pour la réalisation de tableau « look–up-table ». Les données qui ont été
fournies sous forme d'une matrice, sont obtenir par la méthode de Kreigeage [5],

cette dernière est mieux adaptée à la modélisation de la MRV permet d’obtenir un réseau de
courbe flux (I,teta) suffit à caractériser la MRV.

cou
ple 5
éle 8
ctr 4
om
ag
néti 3 6
qu
courant(A)

e(N 2
.m) 4
1

0
2
-1
8
0
6 50
0.8
40
4 30 0.6 50
2 20 40
10 0.4 30
0 0
courant(A) position(°) 20
0.2
10
0

(a) (b)
Figure. IV.5 - Model non linéaire de la MRV: (a) caractéristique flux/téta/courant et (b) la
caractéristique table courant/téta/couple.

159
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

Le couple généré par la MRV est en fonction du flux magnétique donc le couple d’une phase
n’est pas constant pour un courant constant pour différents position du rotor.
Une commande en courant optimisé a permis d'atteindre un taux d'ondulation du couple,
relativement faible, dans la plage à couple maximal constant. Cependant, le couple reste
intrinsèquement très pulsatoire dans la plage des hautes vitesses ce qui peut conduire à une
fatigue du réducteur mécanique.

a. La méthode d’interpolation de Kriegeage

Le Krigeage est la méthode optimale, au sens statistique du terme, d’estimation. On peut

l’utiliser autant pour l’interpolation que l’extrapolation. Le Krigeage porte le
spatiales en deux dimensions. Le Krigeage porte le nom de son précurseur, l’ingénieur minier
sud-africain D.G. Krige. Dans les années 50, Krige développé une série de méthodes
statistiques empiriques afin de déterminer la distribution spatiale de minerais à partir d’un
ensemble de forages. nous présenterons une description du variogramme, le cœur du Krigeage

b. Le principe du Krigeage
Soit une fonction de deux variables f (x,y) à approximer . On peut l’approximer par:

f ( x, y ) = a1 + a 2.x + a3 y + ∑ N α j.g .( ( x − xi) 2 + ( y − yi) 2 (IV-12)

j =1

Avec α1 , α2, α3 et a1 , a2 ,……an solutions du système (IV.13). Ces coefficients sont

appelées coefficients de Krigeage. La fonction f(x,y) passe par les N points de coordonnées
f(xi ,yi ),xi ,yi ).

160
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

⎡ k11 k12 ... k1i ... k1N 1 x1 y1 ⎤ ⎡ α1 ⎤ ⎡ f ( x1, y1) ⎤

⎢ k 21 k 22 ... k 2i ... k 2 N 1 x2 y 2 ⎥⎥ ⎢⎢α 2 ⎥⎥ ⎢⎢ f ( x 2, y 2)⎥⎥
⎢
⎢ : : : : : : : : : ⎥⎢ : ⎥ ⎢ : ⎥
⎢ ⎥⎢ ⎥ ⎢ ⎥
⎢ ki1 ki 2 ... kii ... kiN 1 xi yi ⎥ ⎢ : ⎥ ⎢ : ⎥
⎢ : : : : : : : : : ⎥ ⎢ αi ⎥ = ⎢ : ⎥ (IV-13)
⎢ ⎥⎢ ⎥ ⎢ ⎥
⎢kN1 kN 2 ... kNi ... kNN 1 xN yN ⎥ ⎢αN ⎥ ⎢ f ( xi, iy ) ⎥
⎢ 1 1 ... 1 ... 1 0 0 0 ⎥ ⎢ α1 ⎥ ⎢ 0 ⎥
⎢ ⎥⎢ ⎥ ⎢ ⎥
⎢ x1 x2 ... xi ... xN 0 0 0 ⎥ ⎢α 2 ⎥ ⎢ 0 ⎥
⎢ x1 ⎥ ⎢
0 ⎦ ⎣α 3 ⎦ ⎣ ⎥ ⎢ ⎥
⎣ x2 ... xi ... xN 0 0 0 ⎦

Avec f(xi ,yi )les valeurs mesurées ou calculées de la fonction à interpoler aux points (xi ,yi ).
Les coefficients ki,j sont donnés par la relation suivante:

⎛ 2 2⎞
kij = g.⎜⎜ ( xi − xj ) + ( yi − yj ) ⎟⎟ = g (h) (IV-14)
⎝ ⎠

La fonction g(h) doit être choisie en fonction du problème a résoudre. La fonction g(h)
suivante à été utilisée pour une MRV à une dent par plot

g (h) = h 2 ln(h) (IV-15)

En résume, pour réaliser l’interpolation par Krigeage d’une fonction de deux variables f(x,y),
il est nécessaire de mesurer ou de calculer pour N points la fonction f(xi ,yi).Ensuite les
système (IV.13) est résolu et permet de calculer les coefficients a1 ,a2 ,a2 et ,α1 , α2…..αn
avec la fonction g(h). A partir des coefficients calcules et des points de calcul ou de mesure de
la fonction f(xi ,yi ),il est alors possible avec la fonction de Krigeagede déterminer f(xi ,yi )pour
tout point (xi ,yi )de l’intervalle d’interpolation.

161
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.6.1. Simulation [6] [1]

Comparée avec le schéma de SIMULINK du modèle linéaire dans la figure (IV.3), la

structure de SIMULINK du modèle non–linéaire de la MRV n'a pas changé. D'autre part, les
blocs des phase 1, phase 2 et phase 3 deviennent différents. La figure (IV.12) par exemple,
montre le nouveau bloc de la phase 1. Elle montre deux tableaux « look–up-table »:
flux/position/courant (tableau du courant) et courant/position/couple (tableau du couple).
Ces tableaux ont été obtenus on utilisant le bloc « look–up-table » disponible dans la
catégorie non–linéaire de SIMULINK.

(a) (b)
Figure. IV.6- Caractéristique statique du couple : (a) Données sous forme un tableau
(b) données sous forme graphique.

(a) (b)
Figure. IV.7- Caractéristique statique du flux : (a) Données sous forme un tableau (b) données sous
forme graphique.

162
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

Les résultats de simulation relatifs au modèle non linéaire ont été obtenus à partir de la
commande en tension et la commande en courant par hystérésis.

IV.6.2. La technique de la commande en tension

IV.6.2.1. Topologie de l’onduleur [7]

La topologie la plus utilisée pour le convertisseur est une structure d’onduleur de tension en
demi pont en H asymétrique ce qui se nomme également configuration " classique " par
quelques auteurs. Le schéma structurel d’un tel convertisseur statique est illustré par la figure
( IV.8). Les avantages et les inconvénients de cette topologie sont résumés comme suit :

¾ Il peut être employé avec des machines de n de phase.

¾ Ce convertisseur fournit un degré élevé de flexibilité de commande, parce que le courant

la commande en peut être totalement indépendante pour chaque phase, permettant
l'opération avec recouvrement désiré de phase.

¾ Elle fournit également la fiabilité supplémentaire à la commande de SR système, en plus

de la fiabilité plus élevée inhérente déjà fournie par le SR machine, parce qu'elle ne
permet pas un court-circuit de la barre omnibus courant continu par les commutateurs
principaux (tirer-à travers), l’enroulement de phase est connecté en série entre les
commutateurs principaux.
Cependant, cette configuration n'est pas la plupart économique, parce qu'il exige deux
commutateurs de puissance pour chaque phase. Il y a d'autres topologies qui utilisent moins
de deux commutateurs par phase, mais eux généralement exhibent quelques limitations par
rapport à la topologie " classique ". Voire Figure (IV.8) ci dessus.

163
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

(a) (b)

[Link].8 – La topologie de l’onduleur : (a) Les trois phases de l’onduleur et (b) Une phase
d’onduleur.

IV.6.2.2. Le type des interrupteurs constituent l’onduleur

Chaque phase porte deux interrupteurs de type IGBT, ce composant, qui combine les
caractéristiques suivantes
• La grande simplicité de commande,
• Les faibles pertes par conduction,
• A permis de nombreux progrès dans les applications de l’électronique de puissance, aussi
bien en ce qui concerne la fiabilité que sur l’aspect économique.
•Les transistors IGBT ont permis d’envisager des développements jusqu’alors non viables en
particulier dans la vitesse variable ainsi que dans les applications des machines électriques et
des convertisseurs de puissance qui nous accompagnent chaque jour et partout, sans que nous
en soyons particulièrement conscients : automobiles, trains, métros, bus, avions, bateaux,
ascenseurs, électroménager, télévision, domotique, etc.

164
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

[Link].9 - Interrupteur IGBT

IV.6.2.3. Le type de la commande appliqué [8] [9]

Il existe différents types de la commande plein onde, à savoir la commutation dure (hard
chopping), et la commutation douce (soft chopping), Figure IV.10.

a_ Commande douce b_ Commande dure

c_ Commande pleine onde

Figure IV.10 - Les différents types de commutation.

165
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

La commande choisit pour simulink c’est commutation dure.

Selon le type de la commande choisit dans l’étude les séquences de son fonctionnement sont
résumées dans le tableau suivant :

Cas Q1 Q2 D1 D2 U
1 conducteur conducteur 0 0 +V
2 bloqué conducteur 0 1 0
3 conducteur bloqué 1 0 0
4 bloqué bloqué 1 1 -V

Tableau IV.1 La loi de commande des interrupteurs en commutation dure.

IV.6.2.4. La commande en tension

Supposant une forme d'inductance idéale, les courbes de simulation dans les figures
(e1,e2,e3,e4) illustrent le fonctionnement de la MRV quand elle est commandée par tension.
Le contrôle a lieu en appliquant une tension à une bobine de phase à un angle θon jusqu'à
l'angle θoff. Après cela, la tension appliquée est inversée jusqu'à un certain angle d'extinction
θq qui permet le retour du flux magnétique vers zéro.

Pour appliquer la tension V à une phase, les deux IGBT Q1 et Q2 doivent être ON. Au
contraire, pour appliquer la tension -V et pour assurer la continuité de courant, les deux
diodes D1 et D2 sont conductrices.
La relation de la tension de la phase:
dψ
V = Ri + (IV-16)
dt

∂Ψ ∂Ψ
dΨ = dθ + di
∂θ ∂i (IV-17)

Avec ψ = L(θ)I, on a
dΨ dL di
Et =i ω + L(θ) (IV-18)
dt dθ dt

166
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

di dL
V = Ri + L(θ) + iω (IV-19)
dt dθ

dL
Dans l'équation (VI-18), le terme iω est la tension induite « FEM » et qui sera grande pour
dθ
les grandes vitesses. Pour augmenter la croissance du courant et éviter une grande FCEM,
l'angle θon doit être choisi choisir quand l’inductance et la FEM sont minimums.

Les figures (e1,e2,e3,e4) montrent les résultats pour θon = 0° et θoff = 38° avec la machine
fonctionnant à vide. La différence considérable entre les résultats du modèle linéaire et du
modèle non–linéaire, réside dans le signal du couple. Dans le modèle linéaire, quand l'effet
de la saturation n'est pas pris en compte, le couple de la machine présente de grandes valeurs,
comme montré dans les figure ci dessus. D'autre part, l'amplitude du couple avec le modèle
non linéaire est moins grand dû à la saturation tel que montré dans la figure (e2).

IV.6.3. Le schéma bloc de simulation de MRV plus L’onduleur

tem p
Cl ock
T o Workspace4

control l eur
i abc

U0 U 0 Iref i abc
im puls ion A
VA i
current1
va
teta1
im puls ion B total torque
torque1
vb Iref te w

puls e a total torque w

im puls ion C 1 pos ition
vc 0.1346
posi ti on
Phase1
speed vitesse [rd/s]
converti sseur VB
current2 1 1
1/J 9.031
teta2 s s
Sum
torque2 Sum 1 1/J teta1
teta
Iref

positionpuls e b
F
P hase2
f
VC
torque3

teta 3
c urrent3
Iref Abs

|u|
pos itionpulse c

Phase3
Gai n2

-K - pos ition

Capteur de posi ti on

Fig. IV.11.- Diagramme de simulation de MRV et l’onduleur

167
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

1
current1
va R i1
To Workspace2
R To Workspace3

1 va 1
VA v s
teta Sum1 a
i
switch1 torq 2
table current good
teta torque1
table torque costa3

MATLAB torque signe

Function te1
angle compatibility1
0 To Workspace1
teta s MATLAB
2 Function
teta1 Sum2 modulo pi/2

flux

3
Iref

Relay

switchc

position pulse A 3
4 Position pulse a

-C- Switch logic

Fig. IV.12 - Diagramme d’une phase de MRV

ground
1 impulsion A
impulsion A
g
C

Out1 + va
v 2 Scope1
-
m
E

U1
m
k

Out2
a

convertisseur 1

impulsion B
ground
2 impulsion A
+
v 1
-
Out1 vb
+ U0
v 3
-
U
Out2

i
convertisseur 2 +
-
Scope
Scope2

impulsion C
ground
3 impulsion A

Out1 vc
+
v 4
-
U
Out2

convertisseur 3

Fig. IV.13- Diagramme de simulation de L’onduleur

168
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

U1
2

source to signal VA

g
C
k

m
IGBT2 +
DIODE v

m
-

E
a
2
1
Out1
3
impulsion A
Out2

m
g DIODE
C

a
IGBT1
m
E

ground
1

Fig. IV.14- Diagramme de simulation d’une phase de l’onduleur

IV.6.4. Simulation du régime non linéaire sain

250

200

150
la vitesse [rd/s]

100

50

-50
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s]

Fig.e1 - La vitesse

169
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

le couple total [Nm]

4

0
0.25 0.255 0.26 0.265 0.27 0.275 0.28 0.285 0.29 0.295 0.3
le temps [s]

Fig.e2 – Le couple total

18
ia
16 ib
ic
14
les courants de phases [A]

12

10

0
0.25 0.255 0.26 0.265 0.27 0.275 0.28 0.285 0.29 0.295 0.3
le temps [s]

Fig.e3 – Les courants de phases

170
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

couple phase a
7 couple phase b
couple phase c
6

les couples de phases [Nm] 5

-1
0.25 0.255 0.26 0.265 0.27 0.275 0.28 0.285 0.29 0.295 0.3
le temps [s]

Fig.e4 – Le couple de phases

150

va
vb
100 vc
les tensions de phases [v]

50

-50

-100

-150
0.25 0.255 0.26 0.265 0.27 0.275 0.28 0.285 0.29 0.295 0.3
le temps [s]

Fig.e5 – Les tensions de phases va vb vc

171
 CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.6.5. Simulation avec défauts

IV.6.5.1. Le cas d’un court circuit monophasé 25% dans la phase a à t= 0.2675 [s]

230
250
225

200 220 sain

cc
215
150
la vitesse [rd/s]

210

100 205
0.38 0.4 0.42 0.44 0.46 0.48

50

0
t fault =0.2675 (s)

-50
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s]

(a) (b)
Fig.f1 - La vitesse (a) En cas de court circuit 25% monophasé et - (b) Superposition de la
vitesse de cas sain et court circuit

zoom
25 25
ia avant le cc apres le cc
ib la phase court circuité
ic
20 20
les courants ia ib ic [A]
les courants de phases [A]

15 15

10 10

5 5

0 0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.18 0.2 0.22 0.24 0.26 0.28 0.3 0.32
le temps [s] le temps [s]

Fig.f2 – Les courants de phases

172
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

9
couple phase a
8 couple phase b
couple phase c
7

les couples de phases [Nm]

6

-1
0.18 0.2 0.22 0.24 0.26 0.28 0.3 0.32
le temps [s]

Fig.f3 – Le couple de phases

zoom
9 9
avant le court circuit apres le cc
8 8

7 7
le couple total [Nm]

6 6

5 5

4 4

3 3

2 2

1 1

0 0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.24 0.25 0.26 0.27 0.28 0.29 0.3 0.31 0.32 0.33
le temps [s]

Fig.f4 – Le couple total

173
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.6.5.2. Le cas d’un court circuit triphasé 25%

Le défaut et appliqué à t = 0[s]

300

250

200
la vitesse [rd/s]

150

100

50

-50
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s]

Fig.g1 - La vitesse

6
le couple total [Nm]

0
0.29 0.3 0.31 0.32 0.33 0.34
le temps [s]

Fig.g2 – Le couple total

174
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

20
ia
18 ib
ic
16
les courants de phases [A]

14

12

10

0
0.29 0.3 0.31 0.32 0.33 0.34
le temps [s]

Fig.g3 – Les courants

couple phase a
8 couple phase b
couple phase c
7

6
le coupl des trois phases [Nm]

-1
0.29 0.3 0.31 0.32 0.33 0.34
le temps [s]

Fig.g4 – Le couple de phases

175
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.6.5.3. Le cas d’un circuit ouvert

250
cas sain avant le co la phase a ouverte a et b ouvertes

200

la vitesse [rd/s] 150

100

50

-50
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s]

Fig.h1 – La vitesse

6
le couple total [Nm]

0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s]

Fig.h2 – Le couple total

176
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

ia
20 ib
ic
18

16

14

les courants ia ib ic [A]

12

10

0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s]

Fig.h3 – Les courants de phases

IV.6.5.4. Le cas d’un IGTB d’une phase tombe en panne

Le défaut est apiqué à t=0 (s)

U1
2

IG B T
m
k

diode
a

2
1
Ou
3 t1
i m p u l si o n A
O u t2
m
k
g
C

DIO DE
a

IG B T 1 IG B T
m
E

g ro u n d
1
+
-

so u rce to si g n a l
v

S co p e 2

[Link].15 - Un interrupteur d’une phase de l’onduleur en panne

177
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

9 250

8
200
7

6
150
le couple total [Nm]

la vitesse [rd/s]
5

4 100

3
50
2

1
0
0

-1 -50
0.38 0.39 0.4 0.41 0.42 0.43 0.44 0.45 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s] le temps [s]

Fig.i1 - Le couple total Fig.i2 – La vitesse

-21 couple phase a

x 10
14 9 couple phase b
couple phase c
12 8

7
le couple de la phase a infecté [Nm]

10
les couples des phases [Nm]

6
8
5

6
4

4 3

2
2

1
0
0
-2
0.38 0.39 0.4 0.41 0.42 0.43 0.44 0.45 -1
le temps [s] 0.38 0.39 0.4 0.41 0.42 0.43 0.44 0.45
le temps [s]

Fig.i3 - Le couple de la phase a Fig.i4 – Les couples de phases

-20
x 10 ia
3 18
ib
ic
16
2.5
14
les courants de phases ia ib ic [A]

2 12
le courant ia [A]

10

1.5
8

6
1

0.5
2

0
0 0.38 0.39 0.4 0.41 0.42 0.43 0.44 0.45
0.38 0.39 0.4 0.41 0.42 0.43 0.44 0.45 le temps [s]
le temps [s]

Fig.i5 - Le courant de la phase a Fig.i6 – Les courants de phases ia,ib,ic

178
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

Commentaire

Le cas d’un court circuit monophasé ou triphasé on va commencer par le courant statorique
dans la phase infecté car ces grandeurs fournis des informations accessibles a l’extérieure de
la machine.
En cas de court circuit si on compare le courant de la phase infecté avec les deux autres
phases saines on voie bien que ia augmente de 16[A] à 20[A] donc une amplification de 25%
pour un court circuit de 25% Fig.f2.
Le couple total présente une augmentation juste dans la phase infecté par contre une machine
asynchrone qui présente l’ors d’un court circuit une perturbation sur le long du couple.

La vitesse par contre le régime linéaire en cas de court circuit augmente le même résultat est
obtenu par flux 2d. la vitesse augmente de 212 [rd/s] à 226 [rd/s] fig.1f - (b)

Le cas d’un circuit ouvert ce type de défaut engendre une diminution et fluctuation au niveau
de la vitesse et modulation pour les courant et le couple.

179
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.6.6. Simulation sain et en défauts en cas de commande en courant par hystéresis

Le principe de la commande est le même que celui en un régime linéaire.

t
Clock
Iref time

Iref current1
iabc

vitesse

teta1 torque total torque vit

torque
0
Phase1
vitesse [rd/s]
1 1 0
Iref current2 1/J
s s
Sum
Sum1 speed 1/J teta1
teta

teta2torque2

F
Phase2
f

Abs teta

|u|
Iref current3

teta3torque3

Phase3

[Link].16 - Diagramme de simulation de MRV

1
cu rre
i n t1
R1
i
R

v a1 1
1
Ire f1 s
S um vol ta g e i1
teta1 S um 1 a torq1 2
swi tch 1 teta1 torqu e
ta bl e cu rre n t co sta ta bl e to rq u e costa
si g ne torq ue te 1
te
M AT L AB
Fu n cti o n
a ng l e co m p a ti b i l i ty

0 M AT L AB
Fu n cti o n
teta s
2 Su m 2 m o d u l o p i /2
teta1

v
f
T o Wo rkspa ce 2
T o Wo rksp ace 3

[Link].17 - Diagramme de simulation de MRV par phase

180
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.6.7. Le cas sain

3.5

le couple total éléctromagnétique [Nm]

2.5

1.5

1
0.025 0.03 0.035 0.04 0.045 0.05 0.055
le temps [s]

Fig.j1 – Le couple total

160

140

120

100
la vitesse [rd/s]

80

60

40

20

-20
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
le temps [s]

Fig.j2 – La vitesse

181
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

couple phase c
3.5 couple phase b
couple phase a

le couple des phases a,b,c [Nm]

2.5

1.5

0.5

0
0.025 0.03 0.035 0.04 0.045 0.05 0.055
le temps [s]

Fig.j3 – Les couples de phase

12 ia
ib
ic
10
les courants ia,ib,ic [A]

0
0.025 0.03 0.035 0.04 0.045 0.05 0.055
le temps [s]

Fig.j4 – Les courants de la phases

182
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.6.8. Simulation de défauts

Pour appliqué un défaut de court circuit en régime non linéaire il faut suivre les étapes
suivantes : on doit simulé la MRV par logiciel femm on diminue le nombre de spire (un court
circuit) on fait construire les tables du courant et couple par la méthode d’interpolation
kreigeage ou bicubic spline.

1
current1
i
R
To Workspace1
R

va 1
1
v s
Iref Sum voltage
teta Sum1 a i1
switch1 torq 2
table current good1
teta1 torque
table torque costa1 te1
torque To Workspace4
MATLAB
Function
angle compatibility

0 MATLAB
Function
teta s
2 Sum2 modulo pi/2
teta1

v f
To Workspace2 To Workspace3

[Link].18 – schéma bloc du régime non linéaire avec un défaut de court circuit

183
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.6.8.1. Un court circuit monophasé 25 % dans la phase a à t=0 [s]

court circuit monophasé 25%

180

160

140

120 VITESSE
la vitesse [rd/s]

100

80

60
court circuit monophasé 25%

la vitesse [rd/s]
40 170
168
20 166
164
162
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.3 0.35 0.4 0.45 0.5
le temps [s] le temps [s]

Fig.k1 – La vitesse

7
le couple total electromagnétique [Nm]

3
8
La phase a court circuité
2
6
1 Les deux phses saines

4
0
2
-1
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s] 0
0.39 0.4 0.41 0.42 0.43 0.44 0.45 0.46 0.47 0.48

Fig.k2 – Le couple total

184
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

150
va
vb
100 v3

50
les tensions [v]

-50

-100

-150
0.44 0.445 0.45 0.455 0.46 0.465 0.47
le temps [s]

Fig.k3 – Les tensions

12 ia
ib
ic
10
les courants de phases [A]

0
0.39 0.4 0.41 0.42 0.43 0.44 0.45 0.46 0.47 0.48
le temps [s]

Fig.k4 – Les courants de phases

185
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

12
Des pics dans la phase infectée
10

le courant ia [A] 8

0
0.39 0.4 0.41 0.42 0.43 0.44 0.45 0.46 0.47 0.48
le temps [s]

Fig.k5 – Le courant de la phase a

Le 5
co
upl
e 4
[N
m] 3

0
0.452 0.454 0.456 0.458 0.46 0.462 0.464 0.466 0.468
le temps [s]

7.2
2.2

2
6.8

1.8 6.6
le courant ia [A]

le
cou 6.4
rant
1.6 ia
[A]
6.2

1.4 6

5.8
1.2

5.6

1
5.4
0.4585 0.459 0.4595 0.46 0.4605 0.457 0.4575 0.458 0.4585 0.459
le temps [s] le temps [s]

Fig.k6 – Le couple de la phase a

186
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

couple de la phase a
8 couple de la phase b
couple de la phase c
7

les couples de phases [Nm]

5

-1
0.39 0.4 0.41 0.42 0.43 0.44 0.45 0.46 0.47 0.48
le temps [s]

Fig.k7 – Les couples de phases

IV.6.8.2. Le cas d’un court circuit triphasé 25% à t= 0 [s]

250

200
la vitesse [rd/s]

150

100

50

0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s]

Zoom
203
202.5
la vitesse [rd/s]

202

201.5
201

200.5

200
0.3 0.35 0.4 0.45

Fig.L1 – La vitesse

187
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

6
le couple total [Nm]

0
0.43 0.435 0.44 0.445 0.45 0.455 0.46
le temps [s]
le couple total [Nm]

2
0.448 0.4485 0.449 0.4495 0.45 0.4505
le temps [s]

Fig.L2 – Le couple total

couple phase a
8 couple phase b
couple phase c
7

6
les couples de phases [Nm]

-1
0.43 0.435 0.44 0.445 0.45 0.455 0.46
le temps [s]

Fig.L3 – Les couples de phases

188
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

12
ia
ib
ic
10
les courant de phases ia ib ic [A]

0
0.43 0.435 0.44 0.445 0.45 0.455 0.46
le temps [s]

Fig.L4– Les courants de phases

IV.6.8.3. Le cas d’un circuit ouvert

140
avant le circuit ouvert une phse ouverte deux phse ouvertes a et b
120

100
la vitesse [rd/s]

80

60

40

20

0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s] zoom

Fig.M1 – La vitesse

189
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

cas sain phase a ouverte phase a et b ouverte

3.5

2.5
le couple total [Nm]

1.5

0.5

0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s]

Fig.M2 – Le couple total

3.5

2.5
couple phase c
les couple de phase [Nm]

2 couple phase a
couple phase b

1.5

0.5

-0.5
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s]

Fig.M3 – Les couples de phases

190
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

ia
12 ib
ic

10
les courants de phases [A]

0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s]

Fig.M4 – Les courants de phases

IV.6.9. Commentaire

Le cas d’un régime non linéaire commande en courant pour commenter ces défaut en reste
toujours avec la même théorie appliqué en cas d’une commande en tension sauf que pour ce
cas les défauts présentent des fluctuations et modulations beaucoup plus le cas d’une
commande en tension l’observation importante est l’apparition des pics au niveau du courant
et couple la figure de couple total f3 et k2 explique cette observation et la figure des courant
k6 et f2 ;Pour la vitesse en cas du régime permanent sans défaut sa valeur par commande en
courant égale 150 [rd/s] et en cas d’une commande en tension atteint 250 [rd/s] le cas d’un cc
de 25% fait une augmentation de 150[rd/s] à 170 [rd/s].

191
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

IV.7. Synthèse et Analyse spectrale de défauts par la FFT

Apres avoir simulé les différents défauts en régime linéaire et non linéaire et par deux types
de commande on va analysé le spectre du couple, courant et vitesse en cas d’un court circuit
monophasé en régime non linéaire,

IV.7. 1. Régime non linéaire

IV.7. 1.1. Commande en courant le cas d’un court circuit monophasé 25%

torque sain
8 torque short circuit 80
spectre couple sain
7 70 spectre couple court ciruit

6 60

50
5
le couple total [Nm]

le spectre du couple
40
4
30
3
20
2
10
1
0
0
-10

-1
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 -20
0 50 100 150 200 250 300 350 400 450 500
le temps [s]
la frequence [Hz]

(a) (b)
Fig.N1 – (a) Le couple sain et avec défaut et (b) son spectre

180 90
spectre sain de la vitese
160 spectre de la vitesse en cc
80

140
70
le spectre de la vitesse

120
vitesse saine
la vitesse [rd/s]

60
vitesse court ciruit
100

80 50

60 40

40
30

20
20
0 0 50 100 150 200 250 300 350 400 450 500
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 la frequence [Hz]
le temps [s]

(a) (b)

Fig.N2 – (a) La vitesse saine et avec défaut et (b) son spectre

192
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

12 ia sain 80
ia cc monophasé courant sain
70 courant court circuit
10
60

le spectre du courant
8
50
le courant ia [A]

40
6

30

4
20

2 10

0
0 0 50 100 150 200 250 300 350 400 450 500
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 la frequence [Hz]
le temps [s]

(a) (b)
Fig.N3 – (a) Le courant sain et avec défaut et (b) son spectre

IV.7. 1.2. Commande en tension le cas d’un court circuit monophasé 25%

couple sain
9 couple en court ciruit 70
couple sain
8
60 couple avec cc

7
50
6
le couple total [Nm]

le spectre du couple

40
5

4 30

3
20

2
10
1
0
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
le temps [s] -10
0 50 100 150 200 250 300 350 400 450 500
la frequence [Hz]

(a) (b)

Fig.O1 – (a) Le couple sain et avec défaut et (b) son spectre

193
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

25 80
ia sain courant sain
ia cc 70 courant cc

20 60
le courant de la phse court circuitée

50

le spectre du courant ia
15 40

30

10
20

10
5
0

-10
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 -20
le temps [s] 0 50 100 150 200 250 300 350 400 450 500
la frequence [Hz]

(a) (b)
Fig.O2 – (a) Le courant sain et avec défaut et (b) son spectre

250 100
spectre vitesse saine
spectre du vitesse en court cruit
90
200

80
le spectre de la vitesse

vitesse saine
la vitesse [rd/s]

150 vitesse cc 70

60
100

50

50
40

0 30
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0 50 100 150 200 250 300 350 400 450 500
le temps [s] la frequence [Hz]

(a) (b)

Fig.O3 – (a) La vitesse saine et avec défaut et (b) son spectre

194
CHAPITRE IV MODELISATION, COMMANDE ET SIMULATION DE DEFAUTS PAR MATLAB/SIMULINK

Conclusion
Dans ce chapitre nous avons traité la modélisation et la simulation des machines MRV et
MSRV en régime dynamique dans le cas sain et dans le cas de défauts. Nous avons envisagé
trois régimes à savoir, le régime linéaire commandé en courant, le régime non linéaire
commandé en courant et le régime non linéaire commandé en tension.
Une comparaison des résultats obtenus a été faite en régime linéaire et non linéaire de même
type de commande. De ce faite nous avons établi l’influence de la caractéristique non linéaire
du matériau sur le régime de fonctionnement et sur le comportement de la machine, que ce
soit dans le cas sain ou dans le cas de défauts.
La comparaison du régime non linéaire commandé en courant par celui commandé en tension
montre l’influence de la commande sur le comportement non linéaire.
Les différences entre la méthode analytique et la méthode numérique se manifestent
particulièrement dans le cas du régime non linéaire commandé en tension avec les mêmes
défauts appliqués dans le chapitre précédent à savoir, court circuit monophasé de 25% , court
circuit triphasé de 25 %, circuit ouvert et un interrupteur de l’onduleur en panne.
Nous achevons ce chapitre par une analyse spectrale du courant et du couple en régime non
linéaire.

195
 Références du Chapitre IV

Références du Chapitre IV

[1] F. Soares and P.J. Costa Branco, «Simulation of a 6/4 Switched Reluctance Motor
Based on Matlab/Simulink Environment», Article, Institue Supérieure de Mecatrónica.

[2] F. Soares, B. Costa, «Simulation a 6/4 of motor based on matlab/simulink environment,

IEEE [Link] and Electronic Systems,vol.37, no.3,Jul.2001, pp.989-1009.

[3] Semsudin. Masic, Jasmin. corda, Senad. Smaka, «Computation ofStatic, Steady-state and
Dynamic Characteristics of the Switched Reluctance Motor», ISSN 00051144 ATKAAF
43(3-4),109-117(2002).

[4] Jií. Skalický, «Mathematical Model Switched Relectance Motor»,Article, Dept. of Power
Electrical and Electronic Engineering, FEEC, BUT.

[5] Xiang Dang. Xue, [Link], Member, IEEE, and [Link], «Simulation of Switched
Reluctance Motor Drives Using Two Dimensional BicubicSpline», IEEE Transaction Energie
Conversion, VOL.17,NO.4,Decembre 2002 471.

[6] Alexey. Matveev, «Development of Methods, Algorithms and Software for Optimal
Design of Switched Reluctance Drives», Technische Universiteit Eindhoven, ISBN 90-386-
1783-6, 2006.

[7] K. De Brabandere, J. Driesen, and R. Belmans, «The control of switched reluctance drives
and their use for flywheel energy storage»,

[8] Bernard. Multon, «Possibilités du moteur à réluctance variable à double saillance pour la
motorisation de véhicules électriques. Bilan des essais d'un prototype de 27 kW», C-VELEC
95, Grenoble 1-2 février 1995.

[9] Hans Kuß, Torsten Wichert and Bernard Szymanski, «Design of a high speed Switched
Reluctance Motor for spindle drive»,Article,CPE 2007.

196
 Références du Chapitre IV

[10] [Link], [Link], [Link], [Link], [Link], «Variable-speed

Switched Reluctance Motors», IEE proceedings vol. 128, Pt B, no. 4, pp. 253265, July 1980.

[11] P.O. Rasmussen and G. Andersen, L. Helle, J.K. Pedersen and F. Blaabjerg, «Fully
automatic characterezation system for switched reluctance motors», Institute of Energy
Technology Pontoppidanstraede 101.

197
CONCLUSION
GÉNÉRALE

198
 CONCLUSION GÉNÉRALE

Conclusion Générale et Perspectives

Bien que le principe de fonctionnement des moteurs à réluctance variable (MRV) soit connu
depuis longtemps, leur développement a été freiné pour diverses raisons. Mais aujourd’hui,
les puissants outils de calcul informatiques permettant leur dimensionnement ainsi que les
avancées dans les domaines de l’électronique de puissance et de commande ont contribué à
leur essor avec de multiples applications dans la traction électrique, la conversion des EnR, ....

Le travail présenté dans cette thèse touche au domaine de la maintenance prédictive des
machines électriques qui est devenu un concept stratégique dans tout processus industriel de
production. Notre contribution au diagnostic de défauts a porté sur des machines à réluctance
variable, contrairement à la plupart des études qui se sont intéressées exclusivement à la
machine asynchrone.

Pour aborder notre étude nous avons commencé par un état de l’art sur les principales
techniques de diagnostic et les différents défauts qui peuvent apparaitre dans la machine puis
nous avons modélisé une MSR et une MRV selon deux approches, analytique (sous
Matlab/Simulink) et numérique par éléments finis (sous FEMM et Flux2D), ce qui a permis
de représenter différents défauts pouvant survenir lors d'un fonctionnement normal.

Ces modèles et les simulations qui ont suivi nous ont permis de mettre en évidence l'effet des
défauts sur le fonctionnement de la machine, les signatures associées à chaque défaut et de
suivre les éventuelles évolutions des pannes électriques qui peuvent être générées suite à ces
anomalies.

C’est ainsi que les courts-circuits au bobinage statorique augmentent les amplitudes des
courants dans les phases affectées ; par contre, les défauts d’un circuit ouvert, au contraire,
induisent une diminution du courant de la phase affectée. Pour les excentricités, nous avons
simulé l’excentricité dynamique par logiciel FEMM couplé avec programme LUA ; ce dernier
permet de visualiser des grandeurs inaccessibles dans d'autres types de modèles, évolutions
magnétiques dans la machine (induction, champs, flux ...). Cette analyse a permis de connaitre
le comportement du moteur suite à une anomalie.

199
 CONCLUSION GÉNÉRALE

Par la suite, différentes techniques de diagnostic (avec et sans modèle) ont été étudiées et pour
certaines testées. Ceci a permis de définir les domaines d'utilisation et de voir les limites
rencontrées pour chaque méthode.

En conclusion, ce travail nous a permis de nous familiariser avec les techniques de diagnostic
et de les appliquer à des machines autres que la machine asynchrone afin de détecter les
défauts qui peuvent survenir et de contribuer ainsi au diagnostic précoce et à la maintenance
prédictive des machines réluctantes.

200
ANNEXES

201
 ANNEXE

Annexe A

Le maillage de l’entrefer
La résolution d'un problème par la méthode des éléments finis nécessite le découpage du
domaine d'étude en éléments finis tout en respectant les interfaces de la géométrie.

Car on doit étudier la variation des grandeurs électromagnétiques et mécaniques en fonction

de la position du rotor afin d’avoir les caractéristiques de la machine modélisée , on fait
tourner d’un angle quelconque le rotor, cependant au fur et mesure que l’angle de
déplacement augmente, la distorsion des éléments de la région de l’entrefer (la bande de
roulement) augmente aussi, ce qui provoque des difficultés d’ordre numérique donc le
maillage de la région de l’entrefer doit être régulier (la méthode des éléments finis donne
de bons résultats avec des éléments réguliers).
Il existe différents méthodes pour régler le maillage dans la région d’entrefer :
♦La bande de roulement (fig1)
♦ Le couplage d’interface(fig.2)
♦ Méthode du double entrefer (fig.3)
♦ La ligne de glissement (fig.4)
Sont quelques méthodes rencontres dans une modélisation par éléments finis.

Fig.1 fig.2

Fig.3 Fig.4

202
 ANNEXE
Maillage d'une machine tournante
Pour calculer correctement le couple reluctant d'une machine tournante, son entrefer doit être
maillé avec des éléments de très bonne qualité. En effet, ce résultat dépend de la composante
tangentielle de l'induction dans cet entrefer, composante qui est généralement très faible par
rapport à la composante normale.

Maillage d'une région bande de roulement

La bande de roulement ne doit comprendre qu'une seule couche d'éléments triangulaires,
qui doivent être aussi équilatéraux que possible, et ne peut donc être maillée qu’avec le
mailleur automatique. Le nombre de nœuds sur les contours interne et externe de cette région
peut être, différent, mais la répartition des nœuds tout autour de la bande de roulement doit
respecter les symétries du problème.

Maillage de la bande de roulement

203
 ANNEXE

Annexe B
Les caractéristique de la mrv et la msrv par flux 2d. Les figures (II.8), (II.9), (II.10) illustrent
respectivement pour les deux positions extrêmes: les lignes de flux, la densité de flux, la carte
d’induction.

1. La machine à réluctance variable

Isov alues Result s

Quantity : Equi flux Weber
Isovalues Results
Time (s.) : 1,48E-3 Pos (deg): -799,998E-3 Quantity : Equi flux Weber
Line / Value
1 / -2,83724E-3 Time (s.) : 4,532499E-3 Pos (deg): 44,987
2 / -2,53859E-3 Line / Value
3 / -2,23993E-3 1 / -1,79559E-3
4 / -1,94127E-3 2 / -1,60658E-3
5 / -1,64261E-3 3 / -1,41757E-3
6 / -1,34396E-3 4 / -1,22856E-3
7 / -1,0453E-3 5 / -1,03955E-3
8 / -746,64311E-6 6 / -850,54274E-6
9 / -447,98586E-6 7 / -661,53321E-6
10 / -149,32861E-6 8 / -472,52374E-6
11 / 149,32861E-6
9 / -283,51424E-6
12 / 447,98586E-6
10 / -94,50475E-6
13 / 746,64311E-6
14 / 1,0453E-3 11 / 94,50475E-6
15 / 1,34396E-3 12 / 283,51424E-6
16 / 1,64261E-3 13 / 472,52374E-6
17 / 1,94127E-3 14 / 661,53321E-6
18 / 2,23993E-3 15 / 850,54274E-6
19 / 2,53859E-3 16 / 1,03955E-3
20 / 2,83724E-3 17 / 1,22856E-3
18 / 1,41757E-3
19 / 1,60658E-3
20 / 1,79559E-3

Ligne iso-valeurs pour les deux positions extrêmes

Color Shade Results

Qu an tity : | F lu x d en sity | T esla

T ime (s. ) : 4 , 5 3 2 4 9 9 E -3 P o s (d eg ): 4 4 , 9 8 7
S cale / Co l o r
1 , 3 5 1 5 5 E -6 / 9 0 , 2 1 0 1 E -3
9 0 , 2 1 0 1 E -3 / 1 8 0 , 4 2 0 6 5 E -3
Color Shade Results 1 8 0 , 4 2 0 6 5 E -3 / 2 7 0 , 6 3 0 3 E -3
Qu an tit y : | F lu x d en si ty | T esl a
2 7 0 , 6 3 0 3 E -3 / 3 6 0 , 8 3 9 9 6 E -3
3 6 0 , 8 3 9 9 6 E -3 / 4 5 1 , 0 4 9 5 7 E -3
Time (s. ) : 1 , 4 8 E-3 P o s (d eg ): -7 9 9 , 9 9 8 E -3
4 5 1 , 0 4 9 5 7 E -3 / 5 4 1 , 2 5 9 2 3 E -3
S cale / Co l o r
5 4 1 , 2 5 9 2 3 E -3 / 6 3 1 , 4 6 8 8 9 E -3
1 , 9 0 1 2 1 E -6 / 9 9 , 0 6 8 2 5 E -3
6 3 1 , 4 6 8 8 9 E -3 / 7 2 1 , 6 7 8 5 E -3
9 9 , 0 6 8 2 5 E -3 / 1 9 8 , 1 3 4 5 7 E -3
7 2 1 , 6 7 8 5 E -3 / 8 1 1 , 8 8 8 1 6 E -3
1 9 8 , 1 3 4 5 7 E-3 / 2 9 7 , 2 0 0 9 2 E -3
8 1 1 , 8 8 8 1 6 E -3 / 9 0 2 , 0 9 7 8 2 E -3
2 9 7 , 2 0 0 9 2 E-3 / 3 9 6 , 2 6 7 2 7 E -3
9 0 2 , 0 9 7 8 2 E -3 / 9 9 2 , 3 0 7 4 8 E -3
3 9 6 , 2 6 7 2 7 E-3 / 4 9 5 , 3 3 3 5 8 E -3 9 9 2 , 3 0 7 4 8 E -3 / 1 , 0 8 2 5 2
4 9 5 , 3 3 3 5 8 E-3 / 5 9 4 , 3 9 9 9 9 E -3 1 ,0 8 2 5 2 / 1 ,1 7 2 7 3
5 9 4 , 3 9 9 9 9 E-3 / 6 9 3 , 4 6 6 3 1 E -3 1 ,1 7 2 7 3 / 1 ,2 6 2 9 4
6 9 3 , 4 6 6 3 1 E-3 / 7 9 2 , 5 3 2 6 2 E -3 1 ,2 6 2 9 4 / 1 ,3 5 3 1 5
7 9 2 , 5 3 2 6 2 E-3 / 8 9 1 , 5 9 8 1 E -3 1 ,3 5 3 1 5 / 1 ,4 4 3 3 6
8 9 1 , 5 9 8 1 E -3 / 9 9 0 , 6 6 5 3 8 E -3
9 9 0 , 6 6 5 3 8 E-3 / 1 , 0 8 9 7 3
1 ,0 8 9 7 3 / 1 ,1 8 8 8
1 ,1 8 8 8 / 1 , 2 8 7 8 6
1 ,2 8 7 8 6 / 1 ,3 8 6 9 3
1 ,3 8 6 9 3 / 1 ,4 8 5 1
1 ,4 8 5 1 / 1 , 5 8 5 0 6

Dégradé de l’induction pour les deux positions extrêmes

Region vectors results

Region vectors results Quantity : Flux density Tesla
Quantity : Flux density Tesla
Time (s.) : 4,532499E-3 Pos (deg): 44,987
Time (s.) : 1,48E-3 Pos (deg): -799,998E-3 Colour scale
Colour scale 17,58776E-6 / 84,29772E-3
22,23604E-6 / 90,78197E-3 84,29772E-3 / 168,57785E-3
90,78197E-3 / 181,5417E-3 168,57785E-3 / 252,85799E-3
181,5417E-3 / 272,30143E-3 252,85799E-3 / 337,13812E-3
272,30143E-3 / 363,06116E-3 337,13812E-3 / 421,41825E-3
363,06116E-3 / 453,8209E-3 421,41825E-3 / 505,69839E-3
505,69839E-3 / 589,97852E-3
453,8209E-3 / 544,58063E-3
589,97852E-3 / 674,25865E-3
544,58063E-3 / 635,34036E-3
674,25865E-3 / 758,53879E-3
635,34036E-3 / 726,10009E-3 758,53879E-3 / 842,81892E-3
726,10009E-3 / 816,85982E-3 842,81892E-3 / 927,09905E-3
816,85982E-3 / 907,61955E-3 927,09905E-3 / 1,01138
907,61955E-3 / 998,37929E-3 1,01138 / 1,09566
998,37929E-3 / 1,08914 1,09566 / 1,17994
1,08914 / 1,1799 1,17994 / 1,26422
1,1799 / 1,27066 1,26422 / 1,3485
1,27066 / 1,36142
1,36142 / 1,45218

La densité de flux pour les deux positions extrêmes

Conjonction à t=0.00148 (s) Opposition à t= 0.0045325 (s)

204
 ANNEXE
2. La machine synchrone à réluctance variable

Isovalues Results
Quantity : Equi flux Weber

Time (s.) : 666E-6 Pos (deg): -9,999E-3 Isovalues Results

Line / Value Quantity : Equi flux Weber
1 / -860,71162E-6
2 / -770,1104E-6 Time (s.) : 6,66E-3 Pos (deg): 89,9
3 / -679,50919E-6 Line / Value
4 / -588,90792E-6 1 / -3,53901E-3
2 / -3,16649E-3
5 / -498,3067E-6
3 / -2,79396E-3
6 / -407,70549E-6 4 / -2,42143E-3
7 / -317,10428E-6 5 / -2,0489E-3
8 / -226,50305E-6 6 / -1,67638E-3
9 / -135,90184E-6 7 / -1,30385E-3
10 / -45,30061E-6 8 / -931,31949E-6
11 / 45,30061E-6 9 / -558,79168E-6
12 / 135,90184E-6 10 / -186,26389E-6
13 / 226,50305E-6 11 / 186,26389E-6
12 / 558,79168E-6
14 / 317,10428E-6 13 / 931,31949E-6
15 / 407,70549E-6 14 / 1,30385E-3
16 / 498,3067E-6 15 / 1,67638E-3
17 / 588,90792E-6 16 / 2,0489E-3
18 / 679,50919E-6 17 / 2,42143E-3
19 / 770,1104E-6 18 / 2,79396E-3
20 / 860,71162E-6 19 / 3,16649E-3
20 / 3,53901E-3

Ligne iso-valeurs pour les deux positions extrêmes

Color Shade Results

Qu an ti t y : | F l u x d en si t y | T esl a

T ime (s. ) : 6 6 6 E -6 P o s (d eg ): -9 ,9 9 9 E -3
S cal e / Co l o r
1 6 , 2 2 4 1 E -6 / 5 8 ,8 6 5 2 3 E -3
5 8 , 8 6 5 2 3 E -3 / 1 1 7 , 7 1 4 2 4 E -3
Color Shade Results
Qu an t it y : | F l u x d en si ty | T esl a
1 1 7 , 7 1 4 2 4 E -3 / 1 7 6 , 5 6 3 2 5 E -3
1 7 6 , 5 6 3 2 5 E -3 / 2 3 5 , 4 1 2 2 4 E -3
2 3 5 , 4 1 2 2 4 E -3 / 2 9 4 , 2 6 1 2 5 E -3 T ime (s. ) : 6 ,6 6 E -3 P o s (d eg ): 8 9 ,9
2 9 4 , 2 6 1 2 5 E -3 / 3 5 3 , 1 1 0 2 5 E -3 S cale / Co l o r
3 5 3 , 1 1 0 2 5 E -3 / 4 1 1 , 9 5 9 2 6 E -3 1 2 ,4 6 8 7 6 E -6 / 4 5 0 ,0 0 9 1 1 E -3
4 1 1 , 9 5 9 2 6 E -3 / 4 7 0 , 8 0 8 3 E -3
4 5 0 ,0 0 9 1 1 E -3 / 9 0 0 ,0 0 5 7 6 E -3
4 7 0 , 8 0 8 3 E -3 / 5 2 9 , 6 5 7 2 4 E -3
9 0 0 ,0 0 5 7 6 E -3 / 1 ,3 5
5 2 9 , 6 5 7 2 4 E -3 / 5 8 8 , 5 0 6 3 4 E -3
5 8 8 , 5 0 6 3 4 E -3 / 6 4 7 , 3 5 5 3 2 E -3 1 ,3 5 / 1 ,7 9 9 1
6 4 7 , 3 5 5 3 2 E -3 / 7 0 6 , 2 0 4 3 E -3 1 ,7 9 9 1 / 2 , 2 4 9 1
7 0 6 , 2 0 4 3 E -3 / 7 6 5 , 0 5 3 2 7 E -3 2 ,2 4 9 1 / 2 , 6 9 9 9 9
7 6 5 , 0 5 3 2 7 E -3 / 8 2 3 , 9 0 2 3 7 E -3 2 ,6 9 9 9 9 / 3 , 1 4 9 9 9
8 2 3 , 9 0 2 3 7 E -3 / 8 8 2 , 7 5 1 3 5 E -3
3 ,1 4 9 9 9 / 3 , 5 9 9 9 9
8 8 2 , 7 5 1 3 5 E -3 / 9 4 1 , 6 0 0 3 2 E -3
3 ,5 9 9 9 9 / 4 , 0 4 9 9 8
4 ,0 4 9 9 8 / 4 , 4 9 9 9 8
4 ,4 9 9 9 8 / 4 , 9 4 9 9 8
4 ,9 4 9 9 8 / 5 , 3 9 9 9 7
5 ,3 9 9 9 7 / 5 , 8 4 9 9 7
5 ,8 4 9 9 7 / 6 , 2 9 9 9 7
6 ,2 9 9 9 7 / 6 , 7 4 9 9 6
6 ,7 4 9 9 6 / 7 , 1 9 9 9 6

Dégradé de l’induction pour les deux positions extrêmes

Region vectors results

Quantity : Flux density Tesla

Time (s.) : 666E-6 Pos (deg): -9,999E-3 Region vectors results

Colour scale Quantity : Flux density Tesla
86,15456E-6 / 11,57525E-3
11,57525E-3 / 23,06435E-3 Time (s.) : 6,66E-3 Pos (deg): 89,9
23,06435E-3 / 34,55345E-3 Colour scale
34,55345E-3 / 46,04255E-3
46,04255E-3 / 57,53164E-3 222,20323E-6 / 50,62255E-3
57,53164E-3 / 69,02074E-3 50,62255E-3 / 101,0229E-3
69,02074E-3 / 80,50984E-3 101,0229E-3 / 151,42324E-3
80,50984E-3 / 91,99894E-3 151,42324E-3 / 201,82359E-3
91,99894E-3 / 103,48804E-3 201,82359E-3 / 252,22394E-3
103,48804E-3 / 114,97713E-3 252,22394E-3 / 302,62428E-3
114,97713E-3 / 126,46623E-3 302,62428E-3 / 353,02463E-3
126,46623E-3 / 137,95533E-3 353,02463E-3 / 403,42498E-3
137,95533E-3 / 149,44443E-3
403,42498E-3 / 453,82532E-3
149,44443E-3 / 160,93353E-3
160,93353E-3 / 172,42262E-3 453,82532E-3 / 504,22567E-3
172,42262E-3 / 183,91172E-3 504,22567E-3 / 554,62601E-3
554,62601E-3 / 605,02636E-3
605,02636E-3 / 655,42671E-3
655,42671E-3 / 705,82705E-3
705,82705E-3 / 756,2274E-3
756,2274E-3 / 806,62775E-3

La densité de flux pour les deux positions extrêmes

Conjonction à t=0.000666 (s) Opposition à t= 0.00666 (s)

205
 ANNEXE

Annexe C
L’analyse spectrale
LA MSRV
1. L’analyse spectrale du couple

AXIAL_START AXIAL_CCMONOPH25
N.m N.m
5
5

2,5
CURVE C2D_1 2,5
Mechanics / Axis torque CURVE C2D_13
Time Mechanics / Axis torque
Time
Full cycle / Normal
0 Full cycle / Normal
0

-2,5 (E-3) s.
-2,5 (E-3) s.
75 100 125
75 100 125
N.m
N.m

SPECTRUMSpectr_1 1
From C2D_1
Fundamental 11,904 SPECTRUMSpectr_9
From C2D_13
0,5 Fundamental 11,904
0,5

5 10 15 20 25 30

5 10 15 20 25 30

Fig.c1 Fig.c2

AXIAL_CCTRIPHASE25 AXIAL_OPENCIRCUIT
N.m N.m
4
5
3

2,5 2
CURVE cc CURVE cc
Mechanics / Axis torque Mechanics / Axis torque
Time 1 Time
0 Full cycle / Normal Full cycle / Normal
0

-1
-2,5
(E-3) s. (E-3) s.

75 100 125 75 100 125

N.m N.m
1,5 1,5

1 1
SPECTRUMSpectr_3 SPECTRUMSpectr_3
From cc From cc
Fundamental 11,904 Fundamental 11,904

0,5
0,5

5 10 15 20 25 30 5 10 15 20 25 30

Fig.c3 Fig. c4

AXIAL_INVERTER AXIAL_INVERTER02
N.m N.m
5
2,5
0
CURVE cc CURVE cc
Mechanics / Axis torque Mechanics / Axis torque
0 Time -5 Time

Full cycle / Normal Full cycle / Normal

-10

-2,5
-15
(E-3) s. (E-3) s.
75 100 125 75 100 125

N.m N.m
2 5

1,5 4

3
SPECTRUMSpectr_3 SPECTRUMSpectr_3
1 From cc From cc
Fundamental 11,904 Fundamental 11,904
2

0,5
1

5 10 15 20 25 30 5 10 15 20 25 30

Fig.c5 Fig.c6

206
 ANNEXE

AXIAL_INVERTER03
N.m

CURVE cc
Mechanics / Axis torque
-10 Time

Full cycle / Normal

-20
(E-3) s.
75 100 125

N.m

4
SPECTRUMSpectr_3
From cc
3 Fundamental 11,904

5 10 15 20 25 30

Fig.c7

Par contre la mrv un défaut de 25% pour une machine synchrone à reluctance variable en cas
de défaut monophasé et triphasé Fig.2c et Fig.3c respectivement ne fait apparaisses aucunes
nouvelle fréquence mais il amplifie les fréquences du cas sain ce critère est infiable pour le
diagnostic. Le cas d’un circuit ouvert Fig.4c amplification de l’harmonique 7 est disparition
de l’harmonique 5 et 6.
Le cas d’un circuit ouvert et Fig.c4 et semblable à celui d’un interrupteur court circuit.

207
 ANNEXE
2. L’analyse spectrale du courant statorique ia

AXIAL_START AXIAL_CCMONOPH25
Ampere Ampere

10 10

CURVE i CURVE i
0 Circuit / Current 0 Circuit / Current
Time Time
PA ; PA ;
Full cycle / Normal Full cycle / Normal

-10 -10

(E-3) s. (E-3) s.
75 100 125 75 100 125

Ampere Ampere

7,5 7,5

5 SPECTRUMSpectr_7 5 SPECTRUM Spectr_15

From i From i
Fundamental 11,904 Fundamental 11,904

2,5 2,5

5 10 15 20 25 30 5 10 15 20 25 30

Fig.1d Fig.2d

AXIAL_CCTRIPHASE25 AXIAL_OPENCIRCUIT
Ampere
Ampere
10
10

CURVE i
0 Circuit / Current
Time 5 CURVE i
PA ; Circuit / Current
Full cycle / Normal Time
PA ;
-10 Full cycle / Normal
0

(E-3) s.
-5
75 100 125 (E-3) s.
Ampere
75 100 125

7,5
Ampere

4
5 SPECTRUM Spectr_5
From i
Fundamental 11,904
3
2,5 SPECTRUMSpectr_5
From i
2 Fundamental 11,904

5 10 15 20 25 30
1

5 10 15 20 25 30

Fig.3d Fig.4d

AXIAL_INVERTER03
Ampere

-10

-20 CURVE i
Circuit / Current
Time
PA ;
-30 Full cycle / Normal

-40

(E-3) s.
75 100 125

Ampere
15

10
SPECTRUMSpectr_5
From i
Fundamental 11,904

5 10 15 20 25 30

Fig.5d phase court circuit

208
 ANNEXE

ANNEXE D

La méthode d’interpolation bi-cubic spline

Fig. Flow chart of computation of the interpolation coefficients in bicubic Spline function.

⎡ gj1(i ) ⎤
⎢ gj 2(i )⎥
ψ (θ , i ) = ( gk1(θ ), gk 2(θ ), gk 3(θ ), gk 4(θ )) A⎢ ⎥ ……. 1
⎢ gj 3(i ) ⎥
⎢ ⎥
⎣ gj 4(i )⎦

⎧ gk1 = 1
⎪ gk 2 = (θ − θk ) ⎧ gj1 = 1
⎪ ⎪ gj 2 = (i − ij )
⎨ ⎪
⎪ gk 3 = (θ − θk ) ⎨
2 …….2
⎪ gj 3 = (i − ij )
2

⎪kg 4 = (θ − θk )3
⎩ ⎪kj 4 = (i − ij )3
⎩

209
 ANNEXE

⎛ akj11 ak12 ak13 ak14 ⎞

⎜ ⎟
⎜ akj 21 ak 22 ak 23 ak 24 ⎟
A=⎜ ………….3
akj31 ak 32 ak 33 ak 34 ⎟
⎜ ⎟
⎜ akj 41 ak 44 ⎟⎠
⎝ ak 42 ak 43

⎛ ψkj qkj Ψk , j + 1 qk , j + 1 ⎞ ……………4

⎜ ⎟
⎜ pkj rkj pk , j + 1 rk , j + 1 ⎟
C =⎜
Ψk + 1, j qk + 1, j ψk + 1, j + 1 qk + 1, j + 1⎟
⎜ ⎟
⎜ pk + 1, j rk + 1, j pk + 1, j + 1 rk + 1, j + 1 ⎟
⎝ ⎠

⎛ gk1(θk ) gk 2(θk ) gk 3(θk ) gk 3(θk ) ⎞

⎜ ⎟
⎜ gk '1(θk ) g ' k 2(θk ) g ' k 3(θk ) g ' k 4(θk ) ⎟ …………..5
V (θk ) = ⎜
gk1(θk + 1) gk 2(θk + 1) gk 3(θk + 1) gk 4(θk + 1) ⎟
⎜ ⎟
⎜ g ' k1(θk + 1) g ' k 2(θk + 1) g ' k 3(θk + 1) g ' k 4(θk + 1) ⎟
⎝ ⎠

=
⎛1 0 0 ⎞ 0
⎜ ⎟ ………..6
⎜0 1 0 ⎟ 0
⎜ 1 ∆θk ∆θ 2
∆θ 3 ⎟
⎜ k k⎟
⎜0 1 2∆θk 3∆ θ 2 ⎟
⎝ k⎠

C = V (θk ) A[V (ij )]t ……………7

…………….8
A = [V (θk )]− 1 C [V (ij )]− t

1 1 1 1 1 1 1 1
pk − 1, j + 2[ + ] pkj + pk + 1, j = qk , j − 1 + 2[ + ]qkj + qk , j + 1 =
∆θk − 1 ∆θk − 1 ∆θk ∆θk ∆ij − 1 ∆ij − 1 ∆ji ∆ij
⎡ψkj −ψkj − 1, j ψk + 1, j − ψkj ⎤ ⎡ψkj −ψkj − 1, j ψk + 1, j −ψkj ⎤
3⎢ + ⎥ 3⎢ +
⎣ ∆θk − 1
2
∆θk
2
⎦ ⎥
⎣ ∆i j − 1 ∆ij
2 2
⎦
k = 2,......N − 1 ….9
j = 2,......M − 1 ….10
210
 ANNEXE

1 1 1 1
rk − 1, j + 2[ + ]rkj + rkj + 1, j =
∆θk − 1 ∆θk − 1 ∆θk ∆θk
⎡ qkj − qkj − 1, j qk + 1, j − qkj ⎤ ………11
3⎢ + ⎥
⎣ ∆θk − 1 ∆θk ⎦
2 2

k = 1,......, N − 1

1 1 1 1
rk , j − 1 + 2[ + ]rkj + rk , j + 1 =
∆ ij − 1 ∆ij − 1 ∆ij ∆ij ………..12
⎡ pkj − pkj − 1, j pk + 1, j − pkj ⎤
3⎢ + ⎥
⎣ ∆i j − 1 ∆i j
2 2
⎦
j = 2,......, Μ − 1

∆i j = i j + 1 − i j ………………..13
∆θk = θk + 1 − θk

211
 ANNEXE

ANNEXE E

Le programme de la fft

SIGNAL %%%%%%%%Data sous forme un vecteur obtenir par matlab

N=512;
N2=N/2;
x1=fft(SIGNAL ,N);
fm=1000; %%%%fm=1/dt avec dt=do/n ou n est le nombre d'échontillons;

px1=x1.*conj(x1)/N;
f=fm*(0:N2)/N;

figure (1)
plot(f,px1(1:N2+1));

figure (2)
plot(f,abs(x1(1:257)));

figure (3)
plot(f,20*log10(abs(x1(1:257))),'g');

212
 ‫ﻣﻠﺨﺺ‬
‫اﻝﻐﺮض ﻣﻦ هﺪﻩ اﻻاﻃﺮوﺣﺔ هﻮ ﻧﻤﺪﺟﺔ و ﺗﻤﺜﻴﻞ ﻣﺨﺘﻠﻒ اﻝﺤﻮادث اﻝﻜﻬﺮﺑﺎءیﺔ و اﻝﻤﻴﻜﺎﻧﻴﻜﻴﺔ اﻝﻤﻤﻜﻦ ﺣﺼﻮﻝﻬﺎ ﻋﻠﻰ اﻝﻤﺎآﻨﻴﺎت‬
‫دات اﻝﻤﻤﺎﻧﻌﺔ اﻝﻤﺘﻐﻴﺮة هﺪﻩ اﻻﺣﻴﺮة ﻋﺮﻓﺖ اهﺘﻤﺎم واﺳﻊ ﻓﻲ اﻝﻔﺘﺮة اﻻﺥﻴﺮة و اﺟﺘﺎﺣﺖ ﻣﺠﺎﻻت ﻣﻬﻤﺔ آﺎﻝﻈﺎﻗﺔ و اﻝﻨﻘﻞ و‬
‫ﻋﻠﻰ ﻏﺮار اﻻﻝﺔ اﻝﻼﻣﺘﺰاﻣﻨﺔ اﻝﺘﻲ ﺗﻌﺮف اﺳﺘﻌﻤﺎل واﺳﻊ ﻓﻲ اﻝﺼﻨﺎﻋﺔ و ﺗﻄﻮر آﺒﻴﺮ ﻓﻲ اﻻﺑﺤﺎث ﻣﺜﻞ اﻝﺘﺤﻠﻴﻞ اﻝﻄﻴﻔﻲ‬
‫ﺑﻤﺨﺘﻠﻒ ﺗﻘﻨﻴﺎﺗﻪ وﻝﻬﺪا اﺥﺘﺮﻧﺎ ﻝﺘﻤﺜﻴﻞ اﻝﺤﻮادث اﻝﻤﺪروﺳﺔ و اﻝﻤﻤﺜﻠﺔ ﺑﺎﻝﻄﺮیﻘﺘﻴﻦ ﺗﺎﻋﺪدیﺔ ﺑﻮاﺳﻄﺔ اﻝﻌﻨﺎﺹﺮ اﻝﻼﻣﺘﻨﺎهﻴﺔ‬
. ‫ﺳﺮیﻊ‬FFT ‫ ﺗﺤﻮیﻞ ﺳﻬﻞ اﻝﺘﻄﺒﻴﻖ‬MATLAB/SIMULINK ‫و اﻝﺘﺤﻠﻴﻠﻴﺔ ﻋﻦ ﻃﺮیﻖ‬
‫اﻝﺘﺸﺨﻴﺺ یﻬﺪف ﻻیﺠﺎد اﻝﺪاء و اﻝﺪواء ﻝﺤﻤﺎیﺔ اﻻﻝﺔ و اﻝﻤﻨﺘﻮج ﺥﺎﺹﺔ ﻝﻬﺪا ﻓﻬﺪﻩ اﻝﺪراﺳﺔ هﻲ اﻧﻄﻼﻗﺔ ﻧﺤﻮ ﺗﺤﻘﻴﻖ‬
‫ﺗﺸﺨﻴﺺ ﻣﻦ درﺟﺔ ﻋﻠﻴﺎ اﻝﺘﻨﺒﺎ ﺑﺎﻣﻜﺎﻧﻴﺔ ﺣﺼﻮل ﺥﻠﻞ ﻣﺎ و هﺪا یﻌﺘﻤﺪ ﻋﻠﻰ اﺥﺪ ﺑﻌﻴﻦ‬. ‫اﻻﻋﺘﺒﺎر‬
‫ﺿﺮوف ﻋﻤﻞ اﻻﻝﺔ‬

Abstract
The main purpose of this thesis “dissertation” is to modeling and simulating electrical and
mechanical fault of switched reluctance machine and synchronous reluctance machine,
therefore we used analytical method with MATLAB and numerical method with flux 2d.
These two kinds of machine are from step by step machines family ; this type of machine
knows a wide use in recent years at important fields : energy, transport, as soon as
asynchronous machine witch knows a large utility in industry and a large development in
searching field particular spectral analysis through its various tools and technics, we have
applied fast fourier transformation « FFT» in order to the influence of faults and to achieve
the ability of prediction to confirm about the fault and to determine its kind and therefore we
save the production and the machine as much as possible and for this purpose we focused on
current and moment field.
 Résumé

Ce mémoire de magistère porte sur la modélisation et la simulation des défauts électriques

de la machine à réluctance variable (MRV) et de la machine synchrone à réluctance variable
(MSRV). La modélisation est réalisée selon deux approches : l’une est numérique avec le
logiciel flux 2d, la seconde est analytique avec le logiciel MATLAB. La modélisation de la
MRV et de la MSRV par la méthode des éléments finis nous a permit de simuler différents
types de défauts tel que l’excentricité, circuits ouverts, court-circuits dans les bobinages et
défaut d’interrupteurs de convertisseur. L’analyse spectrale par FFT des signatures des
défauts simulés a permis de faire la lumière entre les raies saines de celles de défauts.