Diesel Turbo PDF
Diesel Turbo PDF
Polytechnique de Lorraine
Centre de Recherche
en Automatique de Nancy
Ecole
doctorale IAEM Lorraine
D
epartement de Formation Doctorale en Automatique
Diagnostic `
a base de mod`
ele :
application `
a un moteur diesel
suraliment
e`
a injection directe
`
THESE
presentee et soutenue publiquement le 20 decembre 2006
pour lobtention du
Zahi SABEH
Composition du jury
President :
D. MAQUIN
Professeur `a lINPL
Rapporteurs :
A. OUSTALOUP
G. GISSINGER
Examinateurs :
J. RAGOT
F. KRATZ
P. DUPRAZ
Professeur `a lINPL
Professeur `a lENSI de Bourges
DELPHI DIESEL SYSTEMS
Remerciements
Cette thse a t ralise au sein de lentreprise Delphi Diesel Systems Blois, en
collaboration avec le Centre de Recherche en Automatique de Nancy (CRAN) et le
Laboratoire de Vision et Robotique (LVR) de Bourges.
Ce travail naurait pu tre men bien sans le soutien de nombreuses personnes que je tiens
remercier et je mexcuse davance auprs de celles que jaurais omis de citer.
Je tiens remercier tout particulirement Monsieur Jos RAGOT, Professeur lInstitut
National Polytechnique de Lorraine et responsable du thme Diagnostic au CRAN,
davoir dirig cette thse avec ses grandes qualits tant sur le plan humain que scientifique.
Ses conseils, sa confiance et son intrt pour mes travaux ont t le carburant de mon
avancement durant les quatre annes de thse.
Jexprime galement ma reconnaissance la direction de Delphi Diesel Systems Blois
reprsente par Monsieur Philipe BERCHER, responsable du centre technique, et Monsieur
Jean MARTIN responsable de lquipe EMS (Engine Management Systems) qui ont permis la
ralisation de cette thse au sein de lentreprise.
Mes sincres remerciements vont Monsieur Frdric KRATZ, Professeur lEcole
Nationale Suprieure dIngnieurs de Bourges et responsable du thme Diagnostic au
LVR, pour son soutien et son enthousiasme vis--vis de la thmatique de mes travaux. Je tiens
souligner ses efforts considrables ainsi que ceux de Monsieur Yves PARMENTIER,
animateur du Ple Capteurs et Automatismes Bourges, qui ont permis le financement de
cette thse grce aux Fonds Europens (FEDER) et aux Fonds de Restructuration de la
Dfense (FRED) dans la rgion Centre.
Je remercie Monsieur Pierre DUPRAZ, responsable de lquipe EMS Gnrique, pour son
encadrement et son soutien tout au long de ce travail ainsi que pour les longues discussions
techniques enrichissantes et ses remarques et corrections du mmoire.
Je tiens exprimer toute ma gratitude Monsieur Grard GISSINGER, Professeur
lUniversit de Haute-Alsace, et Monsieur Alain OUSTALOUP, Professeur lUniversit
de Bordeaux I, pour avoir accept de rapporter sur mon mmoire. Leurs questions et
apprciations mont permis dacqurir un meilleur recul sur le travail effectu. Laccueil
quils mont rserv dans leurs bureaux respectifs, Mulhouse et Bordeaux, ma donn la
chance de dcouvrir deux grandes personnalits scientifiques dans le domaine
automatique/automobile.
Un grand merci galement Monsieur Didier MAQUIN, Professeur lInstitut National
Polytechnique de Lorraine, pour avoir prsid le jury de cette thse et pour ses quelques
corrections minutieuses qui ont servi durant ldition de la version dfinitive de ce mmoire.
Je tiens citer aussi quelques camarades avec qui jai eu de longues discussions et que je
remercie pour tout ce quils ont partag avec moi durant mes annes de thse Blois :
Christophe GAUTHIER, Arnaud GIRARD, Nicolas VACHER, Olivier CAGE, Guillaume
GRATON, Nicolas DUVERNOIS.
Enfin, tous mes remerciements vont mes parents et toute ma famille a qui je ddie ce travail
car ils nont cess de me soutenir et de mencourager tout au long du chemin.
Introduction
7
7
9
9
10
11
11
11
11
13
14
15
17
19
21
21
22
23
24
25
26
27
27
27
28
31
33
35
35
36
37
37
38
38
38
39
41
42
43
44
45
45
46
47
50
50
51
52
53
55
55
56
58
58
62
63
64
65
67
69
71
72
73
74
75
75
77
77
80
81
81
84
85
86
87
89
91
94
95
95
96
96
98
101
103
105
106
106
107
107
108
109
111
112
113
117
119
119
120
121
121
122
124
128
132
133
134
135
135
136
137
138
141
142
145
146
147
148
149
149
150
151
157
157
Conclusion et perspectives
161
Annexes
165
Rfrences
175
Introduction
Introduction
Introduction
Au cours des dernires dcennies, des exigences de plus en plus rigoureuses vis--vis des
missions polluantes des moteurs de vhicules ont pouss les constructeurs dautomobiles
mettre en place des solutions de plus en plus complexes. Par consquent, les risques de
dysfonctionnement de ces systmes ont augment, ce qui a ncessit lintroduction des
systmes de diagnostic embarqu (OBD).
Les systmes OBD, introduits en dbut des annes 90 travers lEurope, le Japon et les Etats
Unis sont obligatoires sur les nouveaux vhicules depuis 2000 [1]. Ces systmes visent
dtecter nimporte quel dfaut capable daugmenter le niveau des missions polluantes audel des seuils rglementaires. La dtection de ce type de dfauts doit donc tre effectue tout
au long de la dure de vie du vhicule et non seulement lorsquil est neuf.
Du fait de ses bonnes performances conomiques et de sa durabilit, le moteur diesel est
dsormais la norme pour les poids lourds en Europe et est de plus en plus rpandu pour les
voitures particulires et les utilitaires lgers. Dans certains pays europens, le march des
voitures particulires se partage dj galit entre moteur diesel et moteur essence. Les
moteurs diesel modernes intgrent : technologie 4 soupapes, turbocompresseur gomtrie
variable et injection directe. La technologie rampe commune , dcline dans une grande
varit de systmes, reprsente de plus en plus ltat de lart [2].
Le dveloppement des moteurs diesel est essentiellement concentr sur la rduction des
missions lchappement, notamment des oxydes dazote (NOx) et des particules. Les
technologies de rduction utilises actuellement sont le recyclage des gaz dchappement
(EGR), pour la rduction de NOx, et, si ncessaire, les filtres particules. Pour lavenir, des
processus de combustion diffrents (par exemple : lallumage par compression charge
homogne HCCI) sont en cours de dveloppement pour les moteurs diesel afin de rduire
encore plus les missions de NOx et de particules.
La complexit des moteurs diesel modernes contrls par llectronique, avec un nombre
croissant de capteurs et dactionneurs, ncessite forcment des systmes de diagnostic
amliors et automatiss bord des vhicules utilisant ces moteurs. La surveillance
approprie et la dtection prcoce des dfauts permettent de rduire la frquence des
oprations de maintenance des vhicules et daider leur rparation lorsque cela est
ncessaire. Ainsi, la fiabilit et la disponibilit des vhicules peuvent tre amliores
considrablement.
Les systmes de diagnostic embarqus actuels sont principalement bass sur les tests de
plausibilit qui consistent vrifier lamplitude de quelques signaux mesurs par rapport
des valeurs limites considres valides. Ces mthodes souffrent aujourdhui de nombreuses
faiblesses cause de leur grande tolrance aux dfauts, ce qui rend ces techniques de moins
en moins adaptes aux exigences futures de diagnostic de plus en plus rigoureux.
Cest pourquoi les techniques de diagnostic base de modle, dveloppes fortement au cours
des dernires annes, semblent trs prometteuses pour lamlioration de la qualit de dtection
des dfauts voire la possibilit de les localiser.
Introduction
Comme son nom lindique, lapproche de diagnostic base de modle repose sur lutilisation
dun modle (physique, mathmatique) du systme tudi (le moteur diesel dans notre cas).
Une comparaison entre les sorties du modle et les signaux mesurs de celles-ci, provenant du
systme rel, peut tre ensuite effectue afin de dtecter et de localiser dventuels
dysfonctionnements de ce systme.
Une des parties importantes du moteur diesel surveiller est la boucle des gaz de celui-ci.
Cette dernire comporte plusieurs capteurs et actionneurs responsables du contrle de lair
dadmission et des gaz dchappement du moteur. Une fuite dans la boucle des gaz ou un
dfaut sur lun de ses capteurs ou actionneurs conduit souvent une augmentation des
missions polluantes du moteur et probablement une baisse des performances ou une
augmentation de la consommation.
Lobjectif de ce travail de thse consiste dans un premier temps modliser la boucle des gaz
dun moteur diesel suraliment injection directe. La suite de ltude vise appliquer une
approche de diagnostic utilisant le modle dvelopp prcdemment dans le but de dtecter et
de localiser trois types principaux de dfauts : ceux des capteurs, des actionneurs et les fuites
de la boucle des gaz.
Aprs cette brve introduction, le droulement de ce mmoire sorganisera de la manire
suivante :
Le chapitre I rappelle les notions de base du fonctionnement des moteurs combustion
interne et notamment des moteurs diesel.
Le chapitre II introduit les diffrents termes et principes relatifs au diagnostic base de
modle. Laccent sera mis, travers ce chapitre, sur les techniques de dtection et disolation
des dfauts en utilisant lapproche des rsidus structurs. Des travaux de diagnostic base de
modles dans le domaine automobile et portant en particulier sur les moteurs combustion
interne seront galement cits.
Le chapitre III voque les aspects matriels et exprimentaux qui ont aid la ralisation de
ce travail de thse. Nous prsentons dans cette partie le vhicule de recherche utilis ainsi que
les diffrents quipements instrumentaux qui ont permis de modliser le moteur du vhicule et
ensuite de valider lapproche de diagnostic utilisant ce modle en temps rel.
Le chapitre IV est entirement consacr la modlisation du moteur diesel et plus
particulirement de la boucle des gaz de celui-ci en utilisant une approche valeur moyenne
MVEM (Mean Value Engin Model). Le systme a t dcompos en plusieurs composants qui
ont t modliss par la suite en se basant soit sur des principes de physique et de
thermodynamique, soit sur des modles de comportement du type boite noire . Des
rsultats de validation de ces modles seront galement prsents dans ce chapitre.
Introduction
Chapitre I
Moteur diesel
suraliment
injection directe
1.2
1.3
Moteur diesel........................................................................................................................................... 9
1.4
Cycle du moteur...................................................................................................................................... 9
1.5
1.6
Rendement volumique.......................................................................................................................... 11
1.7
1.8
1.9
Emissions polluantes............................................................................................................................. 11
1.10
1.11
1.12
1.13
Le but de ce chapitre est de rappeler brivement la faon dont fonctionnent les moteurs
combustion interne et plus particulirement les moteurs diesel suraliments injection directe
(dont le moteur utilis dans cette tude). Les principes thoriques prsents dans ce chapitre et
des informations dtailles sur les moteurs combustion interne peuvent tre trouvs dans [3].
Laccent sera mis sur certains phnomnes et procds fondamentaux sur lesquels nous nous
baserons ultrieurement (cf. Chapitre IV) afin dtablir un modle valeur moyenne du
moteur diesel tudi.
1.1
Les types d'entranement les plus couramment employs dans les vhicules sont les moteurs
combustion interne. Ils dveloppent leur puissance en convertissant en chaleur l'nergie
chimique contenue dans le carburant, puis en transformant cette chaleur en travail mcanique.
La conversion de l'nergie chimique en chaleur s'effectue par combustion et transformation de
l'nergie calorifique en travail mcanique par l'action de celle-ci sur un fluide moteur dont la
pression augmente tout d'abord et dont la dtente successive produit un travail.
Les fluides moteurs envisageables sont les liquides qui, par vaporisation, accroissent la
pression de travail ou les gaz dont la pression de travail peut tre augmente par compression.
La combustion des carburants, gnralement composs d'hydrocarbures, requiert de loxygne
normalement introduit avec lair dadmission. Si la combustion seffectue dans la chambre de
travail, il sagit dune combustion interne. Dans ce cas, les gaz de combustion servent
directement de fluide moteur.
Une production continue de travail mcanique ne peut s'obtenir que par un processus cyclique
(moteur pistons) ou par un processus continu d'absorption de chaleur, de dtente (production
de travail) et de retour du fluide moteur son tat initial (cycle de fonctionnement).
Si le fluide moteur est modifi lors de l'absorption de chaleur, par exemple par l'utilisation
d'une partie de ses composants comme agents d'oxydation, il ne peut retrouver son tat initial
que par remplacement. Ce processus est appel cycle ouvert et se caractrise par un
remplacement cyclique des gaz (limination des gaz de combustion et introduction d'une
charge frache). C'est pourquoi la combustion interne ncessite toujours un cycle ouvert.
1.2
Dans tous les moteurs alternatifs, de l'air ou un mlange air-carburant est priodiquement
comprim dans le cylindre de travail. Soit le mlange est allum, soit le carburant est inject
dans l'air chaud comprim et la combustion du carburant est amorce. Le travail utile gnr
est converti par un mcanisme d'embiellage en un couple disponible l'extrmit de laxe du
moteur appel vilebrequin.
Un moteur est gnralement constitu dun nombre de pistons (4 pistons dans notre cas) qui
montent et descendent dans leurs cylindres respectifs (figure 1.1). La bielle relie le piston au
vilebrequin pour permettre le mouvement rotatif. Le vilebrequin est un arbre commun aux
pistons qui permet de transmettre un mouvement rotatif, liss par un volant dinertie, la
boite de vitesse et aux roues motrices du vhicule.
Injecteur
ou bougie
Arbre cames
Arbre cames
Soupape
Came
Conduit
Piston
PMH
Cylindre
PMB
Bielle
Volant dinertie
Vilebrequin
1.3
Moteur diesel
Conu par Rudolf Diesel au dbut des annes 1900, le moteur diesel est un moteur
combustion interne dont l'allumage n'est pas command mais spontan par phnomne
d'autoallumage (auto-inflammation). Pendant le temps de compression (figure 1.3), lair est
comprim une pression comprise entre 30 et 55 bar (moteurs atmosphriques) ou entre 80 et
110 bar (moteurs suraliments par turbocompresseur) et simultanment chauff une
temprature comprise entre 700 et 900 C. Cette temprature suffit pour provoquer lautoinflammation du carburant inject peu avant la fin de la compression au voisinage du point
mort haut du piston.
1.4
Cycle du moteur
Dans le cadre de notre tude nous utilisons un moteur diesel quatre temps, savoir quun
cycle moteur (correspondant une rotation de 2 x 360) comporte les quatre phases suivantes
(cf. figure 1.3) :
ADMISSION :
- Le piston descend et la soupape d'admission s'ouvre.
- L'air frais s'engouffre dans le cylindre.
COMPRESSION :
- La soupape d'admission se ferme et le piston remonte.
- L'air est comprim et s'chauffe fortement.
- Lorsque le piston atteint les 9/10me de sa course environ,
une quantit dfinie de carburant est injecte.
DETENTE :
- Les soupapes sont fermes.
- Le mlange d'air et de carburant contenu dans le cylindre est sous pression.
- La temprature est suffisamment leve pour que ce mlange brle en
crant une surpression qui pousse le piston vers le bas.
ECHAPPEMENT :
- La soupape d'chappement s'ouvre.
- Le piston remonte chassant les gaz brls vers la sortie d'chappement.
Dans le cas dun moteur quatre cylindres (soit 4 pistons), chacun des pistons effectue les 4
phases dcrites prcdemment et est dcal de 180 (angle vilebrequin) par rapport au piston
avoisinant. Par exemple (cf. figure 1.1), pendant que le premier piston est en fin de phase
d'admission, le deuxime est en fin de phase de compression (dbut de dtente), le troisime
en fin de phase d'chappement et le quatrime en fin de phase d'explosion (dbut
dchappement).
1.5
Les cylindres du moteur sont habituellement relis de chaque ct, admission et chappement,
avec un groupe de tubes appels collecteurs. Il peut y avoir, dans un moteur, un ou plusieurs
collecteurs dadmission et dchappement (cf. figure 1.2). Au maximum, un moteur peut avoir
un collecteur de chaque type par cylindre. Par la suite nous allons fusionner, du point de vue
volumique, les collecteurs multiples du moteur tudi pour ne considrer quun seul collecteur
dadmission et un seul collecteur dchappement.
10
1.6
Rendement volumique
Comme nous lavons vu dans le paragraphe prcdent, un cylindre fonctionne comme une
pompe volumique. Ainsi aprs chaque cycle, le cylindre aspire idalement une masse de gaz
quivalente la variation de volume, cause par le dplacement du piston entre le PMH et le
PMB, multiplie par la densit des gaz dans le collecteur dadmission. Plusieurs facteurs
comme les restrictions de flux formes par les soupapes dadmission, le croisement
douverture (cf. figure 1.7) des soupapes dadmission et dchappement ainsi que les gaz
rsiduels (une partie des gaz brls restant dans les cylindres) contribuent la rduction de
cette masse thorique aspire [4]. Le rendement volumique du moteur, dont la valeur est
normalement suprieure 0.8, permet dexprimer la masse relle des gaz aspirs par le
cylindre selon la formule suivante :
m cyl = vol V d
o :
mcyl
est la masse relle aspire, vol est le rendement volumique, est la densit des gaz dans
1.7
Pertes de pompage
En gnral, dans les moteurs combustion interne la pression dadmission est infrieure
celle dchappement. Le moteur, considr comme une pompe volumique, permet de
comprimer les gaz en les transfrant du collecteur dadmission vers le collecteur
dchappement. Ce transfert cote bien videmment de lnergie mcanique rcupre
directement du vilebrequin du moteur en question.
1.8
Rendement de la combustion
1.9
Emissions polluantes
Selon les conditions de la combustion, le moteur diesel peut rejeter lchappement des gaz
ayant diffrentes compositions chimiques. Idalement, lorsque la combustion diesel est
complte, savoir que tout loxygne et le carburant ncessaires une raction chimique de
combustion quilibre sont utiliss, les seuls rejets seraient de la vapeur deau (H2O) et du
dioxyde de carbone (CO2). En pratique, la combustion nest jamais totale ce qui multiplie le
nombre de composants gazeux ou solides rejets par le moteur diesel dans latmosphre.
Nous mentionnons ci-dessous les polluants dont le rejet est rglement par les autorits
gouvernementales en Europe :
- les hydrocarbures imbrls (HC),
- les oxydes dazote (NOx),
11
Figure 1.4 Valeurs limites dmission en Europe pour voitures moteur diesel.
12
13
Dans la pratique, c'est le calculateur de contrle moteur (ECU) qui joue le rle de calculateur
d'injection, et utilise ce titre l'ensemble des informations issues des capteurs du contrle
moteur : position de la pdale d'acclration, pression dans le collecteur d'admission,
temprature et dbit de l'air l'admission, rgime moteur et position angulaire de l'arbre
cames et de nombreux autres paramtres.
14
15
La deuxime faon de recycler les gaz dchappement passe par une voie externe. Elle
consiste relier le collecteur dchappement celui dadmission en utilisant un tube
mtallique dot dune vanne et souvent dun changeur de refroidissement. La vanne,
contrle par un actionneur, permet de rgler le dbit des gaz recycls. Quant
lchangeur de refroidissement, il utilise le liquide de refroidissement du moteur afin
de diminuer la temprature des gaz recycls ce qui permet de rduire davantage la
formation des NOx dans les cylindres. La figure 1.8 montre un diagramme synoptique
dun systme de recyclage EGR au sein du moteur. En effet, ce concept est utilis dans
le moteur tudi dans le cadre de cette thse.
Figure 1.8 Diagramme synoptique dun systme de recyclage EGR par voie externe.
16
17
18
Chapitre II
19
Prambule.............................................................................................................................................. 21
2.2
2.3
Diagnostic base de modles ............................................................................................................... 22
2.3.1
Terminologie et dfinitions................................................................................................................ 23
2.3.2
Principe de la redondance analytique................................................................................................. 24
2.3.3
Mthode des rsidus........................................................................................................................... 25
2.3.4
Modles des dfauts ........................................................................................................................... 26
2.3.5
Isolation des dfauts........................................................................................................................... 27
2.3.5.1
Les rsidus directionnels .......................................................................................................... 27
2.3.5.2
Les rsidus structurs ............................................................................................................... 27
2.4
20
2.1 Prambule
Notre socit dpend fortement dun nombre important de processus technologiques assez
complexes (processus industriels, nergtiques, informatiques). Outre les exigences de
performances, les exigences de scurit des individus, de lenvironnement et des processus
eux-mmes doivent tre galement satisfaites. Pour ce faire, il est important que tous les
composants dun processus fonctionnent correctement suivant leur domaine nominal. Dans le
cas contraire, il convient de dtecter de faon prcoce tout dysfonctionnement du systme de
faon pouvoir ragir de manire approprie.
Un dfaut par dfinition est un phnomne qui a pour effet de changer le comportement dun
composant dun processus de telle sorte quil ne remplisse plus sa fonction nominale attendue
[11]. Dans les processus complexes et automatiss, les effets dun dfaut peuvent rapidement
se propager conduisant alors la dgradation des performances de ces processus ou, mme
pire, des dfaillances catastrophiques. Cest pourquoi, les dfauts de fonctionnement des
processus doivent tre rapidement dtects et localiss, et des dcisions doivent tre prises
afin dviter la dgradation des performances ou la dfaillance de ces processus.
21
La structure gnrale dune application technique de diagnostic est illustre sur la figure 2.1
o le systme de diagnostic reoit toutes les informations disponibles (appeles aussi
observations) du processus. Ce dernier est souvent rgul par un contrleur qui utilise les
signaux de capteurs et envoie ses commandes aux actionneurs du processus. Ainsi, lensemble
des signaux de sortie de capteurs et de contrleur (ou commandes dactionneurs) constitue les
variables connues (observations) du processus.
Contrleur
Processus
Mesures
Commandes
Dfauts
Perturbations
Systme de
Diagnostic
Diagnostic
des dfauts
Figure 2.1 Structure gnrale dune application technique de diagnostic.
Comme nous pouvons le constater sur cette figure, les perturbations, considres aussi comme
entre inconnue diffrente des dfauts, influencent galement le processus surveill. Lune
des difficults majeures auxquelles le systme de diagnostic est confront consiste
discriminer linfluence des perturbations de celle des dfauts du processus. De plus, le
systme tant boucl au moyen de la loi de commande, certains dfauts ont tendance tre
masqus ou attnus par cette loi de commande et leur dtection peut tre alors rendue
difficile.
Dans le cadre de ce travail de thse, nous nous intressons au diagnostic dun systme
technique (moteur de vhicule) quip dun certain nombre de capteurs et dactionneurs. Les
observations sont typiquement les signaux de sortie des capteurs, tandis que la connaissance
du comportement normal du systme est drive des entres de celui-ci en combinaison avec
son modle mathmatique do lexpression : diagnostic base de modles.
22
couramment utiliss. Au dbut des annes 70, les premiers rapports de recherche sur le
diagnostic base de modles ont t publis notamment dans le domaine arospatial. Puis, la
recherche dans ce domaine na cess de sintensifier durant les annes 80 et 90. Aujourdhui,
le diagnostic base de modles constitue un vaste thme de recherche avec beaucoup de
questions ouvertes et de nombreux ouvrages de rfrence dont [12] [13] [14] et [15].
Les mthodes de diagnostic base de modles ont t dveloppes pour diffrentes familles
de modles parmi lesquelles nous pouvons citer :
- les modles du domaine de lIntelligence Artificielle qui sont souvent bass sur la logique
[16],
- les modles des systmes dynamiques vnements discrets [17][18].
- les modles utilisant des variables continues temps discret ou continu. Ces modles sont
habituellement exprims travers dquations diffrentielles, de rcurrence ou statiques. Il
faut noter que les stratgies de diagnostic base de modles ont des performances
directement lies la qualit des modles, cest--dire leur aptitude dcrire
correctement le fonctionnement du systme surveiller. Cela dit, un modle mme
imprcis peut tre utilis condition dassortir le diagnostic dun intervalle de confiance.
Dans le cadre de cette thse, seule cette dernire famille de modles est utilise.
24
25
Aprs la gnration des rsidus, il est ncessaire de dcider si leur amplitude est suffisamment
importante pour rvler une anomalie de fonctionnement du systme. Il sagit ici de ltape
dvaluation des rsidus dont la valeur de sortie est de nature binaire (0 ou 1), car elle infirme
ou confirme la prsence dun dfaut dans le systme.
Une solution simple consiste comparer le niveau du rsidu un seuil prdfini (qui prend en
compte par exemple le niveau de bruit des capteurs et/ou une estimation de lintervalle de
confiance du modle de systme utilis) au-del duquel la prsence dun dfaut est
confirme. Dautres solutions plus complexes bases sur des approches statistiques peuvent
tre trouves par exemple dans [13].
Dans le cadre de cette thse, un systme de diagnostic base de modles sera propos afin de
dtecter les dfauts dun moteur de vhicule. Ce systme utilisera un gnrateur de rsidus
entirement bas sur les notions thoriques prsentes prcdemment dans ce paragraphe. La
solution dvaluation des rsidus par seuil sera galement retenue pour notre application.
y = Cx + Du
Les dfauts additifs des capteurs et des actionneurs peuvent tre modliss en utilisant
respectivement les signaux additifs fc et fa de la faon suivante :
x& = Ax + Bu + Bf a
y = Cx + Du + f c
y = Cx + Du
Dans la suite de cette thse, nous ne considrons que les modles de dfauts additifs pour
tudier linfluence des dfauts de capteurs et dactionneurs du moteur. De plus, dautres types
26
de dfauts, comme la fuite des gaz dans le moteur par exemple, peuvent tre tudis en
utilisant des modles additifs. Lexemple suivant montre brivement la faon dont les fuites
peuvent tre modlises :
Considrons un systme dadmission dair dans un moteur de vhicule. La dynamique de la
pression dans le collecteur dadmission du moteur peut tre dcrite par lquation suivante :
p& =
RT
m& thr m& cyl
V
o p est la pression dans le collecteur, m& thr est le dbit massique de lair passant par le volet
dadmission, m& cyl est le dbit massique de lair entrant dans les cylindres,
R, T, V sont des constantes.
Afin de modliser, par exemple, une fuite dans le collecteur dadmission, un dbit dair
additionnel f sortant du collecteur peut tre intgr dans lquation prcdente :
p& =
RT
m& thr m& cyl f
V
27
Figure 2.3 Exemples de structures dinfluence utilises pour lisolation des dfauts
Ainsi, aprs la dtection dun dfaut, lisolation de ce dernier peut idalement tre accomplie
en identifiant la colonne de structure dinfluence associe qui correspond aux valeurs des
rsidus obtenues (aprs ltape dvaluation). Par exemple, si seuls les rsidus r1 et r3 ont
ragi un dfaut du systme, lisolation selon la deuxime structure dinfluence de la figure
2.3 peut conclure que le dfaut f2 est prsent.
La condition ncessaire pour que tous les dfauts inclus dans une structure dinfluence soient
isolables est que les colonnes de cette structure soient parfaitement distinctes. La structure
dinfluence satisfaisant cette condition est alors appele structure isolante. Les structures 2 et
3 de la figure 2.3 sont donc isolantes tandis que la structure 1 ne lest pas puisque les dfauts
f1 et f2 ont la mme signature.
Enfin, cette approche disolation base sur les rsidus structurs sera retenue pour notre
application de diagnostic du moteur diesel dveloppe notamment dans le Chapitre VI de ce
mmoire. Dautres informations et dveloppements intressants autour des mthodes
structures disolation des dfauts et de leurs applications peuvent tre trouves dans [23] et
[24].
28
nombreuses tudes ont t menes [25] [26] et qui prouvent le gain considrable que peut
offrir le diagnostic base de modles aux moteurs de vhicules.
Dans le domaine des moteurs essence nous pouvons citer [27] [28] [29] et [30] pour le
diagnostic des dfauts du circuit dadmission en utilisant lapproche des tests dhypothses
structurs. Des travaux bass sur la mme approche ont t mens sur les moteurs diesel pour
le diagnostic de la boucle des gaz (dfauts de capteurs et fuites) [31] et du systme dinjection
directe rampe commune [32]. Une mthode base sur des modles de rseaux de neurones et
qui utilise lapproche des rsidus structurs est prsente dans [33] pour le diagnostic des
systmes dinjection et de suralimentation dun moteur diesel.
Dans le domaine du post-traitement des missions polluantes, par exemple, une approche de
diagnostic dun pot catalytique dun vhicule est prsente dans [34]. Cette approche est
base sur des tests statistiques de lerreur de prdiction du rapport air/carburant dans le pot.
Dans le domaine de la suspension et du comportement dynamique des vhicules, des
mthodes de diagnostic base de modles semi-physiques utilisant des rseaux de neurones et
une approche de rsidus structurs sont proposes dans [35] et [36]. Une approche base sur
lestimation paramtrique est propos dans [42] pour dtecter des dysfonctionnements
damortisseurs. Des travaux intressants de modlisation de systmes de suspension peuvent
tre trouvs galement dans [46].
Enfin, des mthodes destimation et de modlisation dans le domaine dadhrence de
vhicules, qui peuvent tre trs utiles pour le diagnostic, sont proposes dans [43] [44] et [45].
29
30
Installation exprimentale
Chapitre III
Installation exprimentale
31
Installation exprimentale
Motivations............................................................................................................................................ 33
3.2
3.3
Calculateur ............................................................................................................................................ 35
3.4
Module de drivation............................................................................................................................ 36
3.5
PC portable ........................................................................................................................................... 37
3.6
Capteurs additionnels........................................................................................................................... 37
3.6.1
Pression.............................................................................................................................................. 38
3.6.2
Temprature ....................................................................................................................................... 38
3.6.3
Vitesse angulaire................................................................................................................................ 38
32
Installation exprimentale
Le moteur que nous avons modlis dans le cadre de cette thse est le DaimlerChrysler
OM646. Cest un moteur diesel suraliment quatre cylindres et quip dun systme
dinjection directe rampe commune fourni par Delphi Diesel Systems.
Equipant la Mercedes 220 E-class, vhicule de recherche que nous avons utilis (voir figure
3.1), ce moteur a une cylindre de 2.2 litres et dispose dun systme dadmission papillon.
Il est galement quip dun circuit de recyclage des gaz dchappement EGR et de
turbocompresseur gomtrie variable.
Dans le cadre de cette thse, beaucoup de travail a t consacr la mise en place dune
installation exprimentale dans le vhicule de recherche ainsi quau choix des capteurs
additionnels et la recherche de solutions diffrents problmes techniques.
3.1 Motivations
Plusieurs capteurs additionnels ont t installs dans le moteur tudi et plus prcisment au
niveau de la boucle des gaz de ce dernier. Ces modifications ont t ncessaires pour tudier
lallure de certaines grandeurs physiques (pression, temprature) dans cette boucle et pour
aider la construction et la validation du modle de cette boucle.
33
Installation exprimentale
Ainsi, la ncessit dune centrale dacquisition capable de recueillir la fois les signaux des
capteurs additionnels et ceux dorigine nous a sembl vidente.
Par ailleurs, le calculateur lectronique est responsable de la gestion du moteur. Il reoit les
informations des capteurs dorigine et effectue, entre autres, le contrle des actionneurs de la
boucle des gaz.
Par consquent, il tait ncessaire de connatre les signaux de commande de ces actionneurs
parmi dautres car ils constituent, comme nous le verrons par la suite, les entres du modle
de la boucle des gaz.
Enfin, le but final de ce travail est de proposer une solution de diagnostic de la boucle des gaz
en temps rel.
Il tait donc important de pouvoir appliquer, tester et valuer directement les algorithmes de
diagnostic mis au point en utilisant le vhicule de recherche et sans en modifier le
fonctionnement dorigine.
Compte tenu de toutes les rflexions mentionnes ci-dessus, nous avons dcid dembarquer
dans le vhicule de recherche linstallation exprimentale, illustre sur la figure 3.2, dont nous
allons dtailler les diffrentes fonctions.
Ethernet
CAN
PC Portable
Module de
drivation
XPC Target Box
Habitacle
Calculateur (ECU)
Moteur
Capteurs
additionnels
Capteurs
dorigine
Actionneurs
Par la suite nous parlerons uniquement des lments rajouts la configuration dorigine du
vhicule de recherche dans cette installation exprimentale.
34
Installation exprimentale
3.3 Calculateur
Appel souvent ECU (Electronic Control Unit), le calculateur lectronique est responsable du
contrle du moteur. Ainsi, ce calculateur est lorigine de toute action mene pour dfinir le
point de fonctionnement du moteur (ouverture dinjecteur, positionnement des ailettes de la
turbine, ouverture du volet dadmission).
Cest pourquoi, il est primordial de pouvoir accder aux diffrents paramtres des contrleurs
implments de faon logicielle dans le calculateur pour pouvoir les modifier selon les
besoins. Pour ce faire, le calculateur standard est remplac par un calculateur modifi. Ce
dernier est quip de mmoire modifiable (RAM), accessible par lutilisateur, permettant ainsi
la modification des paramtres du programme de contrle. Cette modification peut tre
effectue depuis un PC matre grce un logiciel de mise au point et une liaison
CAN (Control Area Network) avec le calculateur.
Plusieurs modifications de paramtres ont t effectues dans le calculateur durant les travaux
de cette thse dont nous pouvons citer quelques exemples :
- La modification du positionnement statique ou dynamique (chelon, squence binaire
pseudo alatoire) des actionneurs de la boucle des gaz en vue lidentification des
paramtres de son modle.
- La modification des paramtres de conversion des signaux lectriques provenant des
capteurs de la boucle des gaz (pression, temprature, dbit) pour lanalyse de leurs
dfauts en ligne.
35
Installation exprimentale
1
2
CAM
Q10
V
D1
D2
D3
L1
L2
L3
Dfinition
Vecteur de temps du calculateur
Vecteur de temps de la XPC Target Box
Pression atmosphrique
Pression aprs filtre air
Pression aprs compresseur
Pression dans le collecteur dadmission
Pression dans le collecteur dchappement
Pression aprs turbine
Dbit dair traversant le compresseur
Temprature aprs filtre air
Temprature aprs compresseur
Temprature aprs changeur de refroidissement dair
Temprature dans le collecteur dadmission
Temprature dans le collecteur dchappement
Temprature aprs turbine
Temprature aprs changeur de refroidissement dEGR
Temprature de liquide de refroidissement
Temprature aprs conduit de retour EGR
Vitesse angulaire du turbocompresseur
Source
ECU
XPC TB
XPC TB/ECU (Capteur dorigine)
XPC TB/ECU (Capteur dorigine)
XPC TB (Capteur additionnel)
XPC TB/ECU (Capteur dorigine)
XPC TB/ECU (Capteur dorigine)
XPC TB (Capteur additionnel)
XPC TB/ECU (Capteur dorigine)
XPC TB/ECU (Capteur dorigine)
XPC TB (Capteur additionnel)
XPC TB/ECU (Capteur dorigine)
XPC TB (Capteur additionnel)
XPC TB/ECU (Capteur dorigine)
XPC TB (Capteur additionnel)
XPC TB (Capteur additionnel)
XPC TB/ECU (Capteur dorigine)
XPC TB (Capteur additionnel)
XPC TB (Capteur additionnel)
Rgime du moteur
Tableau 3.1 Liste des variables connues issues de linstallation exprimentale de la boucle des gaz.
36
Installation exprimentale
3.5 PC portable
Trs pratique pour les applications embarques, le PC portable nous a permis de superviser le
fonctionnement des deux plateformes temps rel : le calculateur de contrle du moteur et la
XPC Target Box durant les diffrentes phases du projet.
Un logiciel de supervision et de mise au point (VISU), dvelopp par Delphi Diesel Systems,
permet au PC portable de communiquer avec le calculateur via un bus CAN.
Grce cette communication, certains paramtres du calculateur peuvent tre modifis
comme nous lavons expliqu prcdemment. De plus, il est possible deffectuer et de
conserver sur le disque dur du PC portable des enregistrements de variables connues dans
lenvironnement du calculateur (telles que le dbit dinjection Q10 ou le rgime du moteur
2).
De faon similaire, la boite outils XPC Target dans lenvironnement Matlab/simulink
permet au PC portable de communiquer avec la cible XPC Target Box grce une liaison
Ethernet.
Ainsi, il est possible de superviser le fonctionnement en temps rel et de transfrer les
applications excuter vers la cible ou les donnes de mesures effectues par celle-ci vers le
PC portable.
Ces donnes peuvent tre stockes sous forme de fichiers sur le disque dur du PC portable
pour tre traites ultrieurement.
Deux jeux de donnes sont produits pour chaque enregistrement du fait que les variables
mesures (dfinies dans le tableau 3.1) peuvent provenir de deux sources diffrentes : le
calculateur (ECU) et la XPC Target Box. Par consquent, une synchronisation temporelle
entre les deux jeux de donnes savre ncessaire.
Pour ce faire, un signal de rfrence doit tre connect aux deux sources. Le signal CAM
dcrit dans le tableau 3.1 nous a sembl trs intressant pour la synchronisation des donnes
laide des deux vecteurs de temps : tECU et tXPC.
De plus, la priode de ce signal binaire correspond un cycle complet du moteur ce qui
permet le calcul des valeurs moyennes de toutes les variables utilises dans le modle.
37
Installation exprimentale
3.6.1 Pression
Nous avons utilis des capteurs de pression identiques ceux dorigine pour les points de
mesure additionnels. Deux modles diffrents de ces capteurs ont t utiliss en fonction des
plages de pression mesurer ainsi qu leur emplacement dans la boucle des gaz :
-
le premier modle caractris par une plage de mesure de [0 3] bar ; il supporte des
tempratures jusqu 140C et il est ddi ladmission pour mesurer par exemple la
pression aprs compresseur P2.
le deuxime modle est plutt utilis pour la partie dchappement car il peut supporter
jusqu 700C et peut mesurer des pressions atteignant 5 bar.
Les deux modles utiliss possdent une bande passante de 2 kHz environ, ce qui couvre
largement les frquences dues aux pompages du moteur (150 Hz 4500 tours/min).
3.6.2 Temprature
Les capteurs qui ont t utiliss pour les mesures de temprature additionnelles sont des
thermocouples du type K ayant un diamtre de 0.5mm. Ces capteurs peuvent mesurer des
tempratures situes dans une plage de [-20 1150] C avec une prcision de 0.3%.
Quant au comportement dynamique de ces capteurs, il correspond assez bien un systme du
premier ordre avec une constante de temps denviron 1.4 s.
38
Modlisation du systme
Chapitre IV
Modlisation du systme
39
Modlisation du systme
4.2
Dcomposition du systme ................................................................................................................... 42
4.2.1
Filtre air........................................................................................................................................... 43
4.2.1.1
Modle de dbit ........................................................................................................................ 44
4.2.1.2
Modle de pression................................................................................................................... 45
4.2.2
Compresseur ...................................................................................................................................... 45
4.2.2.1
Modle de dbit ........................................................................................................................ 46
4.2.2.2
Modle de temprature ............................................................................................................. 47
4.2.2.3
Modle de puissance................................................................................................................. 50
4.2.2.4
Modle de pression................................................................................................................... 50
4.2.3
Echangeur de refroidissement de lair comprim............................................................................... 51
4.2.3.1
Modle de dbit ........................................................................................................................ 52
4.2.3.2
Modle de temprature ............................................................................................................. 53
4.2.3.3
Modle de pression................................................................................................................... 55
4.2.4
Volet dadmission .............................................................................................................................. 55
4.2.4.1
Modle de dbit ........................................................................................................................ 56
4.2.5
Collecteur dadmission ...................................................................................................................... 58
4.2.5.1
Modle de pression................................................................................................................... 58
4.2.5.2
Modle de temprature ............................................................................................................. 62
4.2.5.3
Modle de densit..................................................................................................................... 63
4.2.6
Chambre de combustion..................................................................................................................... 64
4.2.6.1
Modles de dbit ...................................................................................................................... 65
4.2.6.2
Modle de temprature ............................................................................................................. 67
4.2.6.3
Modle du rapport Air/Carburant () ....................................................................................... 69
4.2.7
Collecteur dchappement.................................................................................................................. 71
4.2.7.1
Modle de pression................................................................................................................... 72
4.2.7.2
Modle de temprature ............................................................................................................. 73
4.2.8
Systme de recyclage des gaz dchappement (EGR)........................................................................ 74
4.2.8.1
Conduit de retour...................................................................................................................... 75
4.2.8.1.1 Modle de temprature ........................................................................................................ 75
4.2.8.2
Electrovanne............................................................................................................................. 77
4.2.8.2.1 Modle de dbit ................................................................................................................... 77
4.2.8.2.2 Modle de temprature ........................................................................................................ 80
4.2.8.3
Echangeur de refroidissement .................................................................................................. 81
4.2.8.3.1 Modle de temprature ........................................................................................................ 81
4.2.9
Turbine............................................................................................................................................... 84
4.2.9.1
Modle de dbit ........................................................................................................................ 85
4.2.9.2
Modle de temprature ............................................................................................................. 86
4.2.9.2.1 Sparation des rendements................................................................................................... 87
4.2.9.2.2 Rendement arodynamique.................................................................................................. 89
4.2.9.2.3 Rendement d lchange de chaleur ................................................................................. 91
4.2.9.3
Modle de pression................................................................................................................... 94
4.2.9.4
Modle de puissance................................................................................................................. 95
4.2.9.5
Modle de vitesse ..................................................................................................................... 95
4.2.10
Systme dchappement................................................................................................................ 96
4.2.10.1
Modle de dbit ........................................................................................................................ 96
4.2.11
Actionneurs lectriques ................................................................................................................. 98
40
Modlisation du systme
Le but de cette thse est dtablir un diagnostic de la boucle des gaz dans le moteur en
utilisant un modle de bon fonctionnement de celle-ci. Nous dcrivons donc dans ce chapitre
la modlisation de la boucle des gaz du moteur diesel DaimlerChrysler OM646.
Le modle utilis dans cette tude est un modle valeur moyenne, ce qui veut dire que les
variables et les signaux traits sont moyenns sur un cycle moteur complet.
Lobjectif principal durant la construction de ce modle a t de le rendre le plus simple
possible. Un modle destin tre utilis dans une application automobile en temps rel
doit combiner prcision et simplicit et tre fiable sans pour autant tre complexe.
Figure 4.1 Vue globale de la boucle des gaz du moteur diesel DaimlerChrysler OM646.
41
Modlisation du systme
1
n
i =1
xi xi
Ermax = max
1 i n
xi
xi xi
xi
o :
x
x,
n
42
Modlisation du systme
Ermoy , Ermax sont respectivement lerreur relative moyenne et lerreur relative maximum entre
Figure 4.2 Diagramme schmatique de la boucle des gaz du moteur diesel DaimlerChrysler OM646.
43
Modlisation du systme
C T (t )
P0 (t ) P1 (t ) = C1 Q0 (t ) 14 1
P0 (t )
avec P0 (t ) P1 (t )
(4.1)
P0 (t )
C14 T1 (t )
[P0 (t ) P1 (t )]
(4.2)
C1
o :
C1 est le coefficient de dcharge prenant en compte des notions physiques lies la restriction
telles que la section effective de passage et les frictions avec les parois,
C14 est la constante du gaz parfait de lair dadmission,
P0
T1
Q0
La figure 4.3 montre les rsultats de validation de la relation entre le carr du dbit volumique
de lair et la perte de charge travers le filtre dcrite par la formule (4.1).
Mesure
Modle
2000
1800
1600
P0-P1 (Pa)
1400
1200
1000
800
600
400
200
0
0.002
0.004
0.006
0.008
0.01
0.012
0.014
T0 2
(Q0 C14
P0 )
44
Modlisation du systme
Etant donn que lhypothse dun gaz incompressible est retenue ici, nous considrons quil
ny a pas de diffrence entre la temprature en entre et celle en sortie du filtre air.
C14 T1 (t )
C0
[Q (t ) Q (t )] dt
0
(4.3)
t0
o :
C 0 est le volume du filtre,
C14
P1
T1
Q0
Q1
A partir des quations 4.2 et 4.3, la figure 4.4 indique les variables dentres/sorties du
modle global du filtre air, ce dernier faisant ensuite lobjet dune simulation.
Q1 (Kg/s)
T1 (K)
P1 (Pa)
P0 (Pa)
4.2.2 Compresseur
Dans lensemble turbocompresseur, lnergie des gaz dchappement du moteur rcupre
grce la turbine et transmise par larbre turbocompresseur est utilise pour entraner le
compresseur. Cette opration a pour but de comprimer lair dadmission et de suralimenter le
moteur diesel.
45
Modlisation du systme
Le modle du compresseur comporte quatre parties : trois parties statiques ddies au dbit
dair, la temprature en sortie du compresseur et la puissance absorbe par celui-ci. Quant
la quatrime partie, elle est dynamique et a pour but de dcrire la pression en sortie du
compresseur.
P1 (t ) Tref
Pref T1 (t )
(4.4)
avec :
P (t )
1cor (t ) = 1 (t )
(4.5)
Tref
T1 (t )
o :
P1 , P2
T1
P2 (t )
P2 (t )
P2 (t )
P2 (t )
Q1cor (t ) = f Q1cor
, 1cor (t ) = a
b
c
d
+
+
+
P (t )
P (t )
P1 (t )
P1 (t )
1
1
(4.6)
46
Modlisation du systme
o :
a, b, c, d, e, f sont les constantes des polynmes utiliss,
P1 , P2
Les courbes pointilles traces sur la figure 4.4 illustrent les rsultats de modlisation de dbit
corrig Q1cor en fonction du rapport de pressions
P2
P1
compresseur. La superposition des rsultats du modle aux mesures nous permet de qualifier
la prcision obtenue comme tant trs satisfaisante (Ermoy = 1.8%, Ermax = 3.6%).
0.25
Nominal
Limitation
de dbit
0.2
N1
N2
N3
N4
N5
N6
N7
0.15
0.1
0.05
Instabilit
0
1.5
2.5
Rapport de pressions
Figure 4.5 Domaine de fonctionnement du compresseur dfini par des points de mesures,
chaque rgime de rotation 1=Ni (avec N1< N2 < Ni) est reprsent par un
groupe de points, les courbes pointilles reprsentent les rsultats de modlisation.
47
Modlisation du systme
T2isent (t ) P2 (t )
=
T1 (t )
P1 (t )
C7
(4.7)
o :
T2isent est la temprature thorique en sortie du compresseur dans le cas isentropique,
T1
C7 =
1
, = c p / cv
volume constant.
Du fait de lirrversibilit travers le compresseur (due en particulier aux pertes thermiques
ou frottement), la compression, en ralit, nest pas vraiment isentropique.
Ainsi, le rendement isentropique du compresseur 1 est introduit pour dterminer la relation
entre la monte thorique de la temprature due la compression et celle obtenue en pratique:
1 (t ) =
T2isent (t ) T1 (t )
, 0 < 1 (t ) 1
T2 (t ) T1 (t )
(4.8)
o :
T2 est la temprature relle de lair observe en sortie du compresseur.
T2isent est la temprature thorique en sortie du compresseur dans le cas isentropique,
T1
En remplaant T2isent dans (4.8) par sa valeur de (4.7) nous obtenons la relation modlisant la
temprature T2 :
T (t ) P (t )
T2 (t ) = T1 (t ) + 1 2
1 (t ) P1 (t )
C7
(4.9)
o :
1 est le rendement isentropique du compresseur.
Llment principal dans ce modle de temprature propos est bien entendu le rendement
isentropique qui doit tre identifi, de faon similaire au dbit Q1cor , dans le domaine de
fonctionnement nominal du compresseur.
Pour cela, nous proposons la fonction f1 qui relie le rendement isentropique du compresseur
au dbit corrig de ce dernier ainsi qu son rgime de rotation corrig.
Nous utilisons pour cette fonction un modle polynomial dont nous avons identifi les
paramtres par optimisation par rapports aux valeurs de rendement calcules.
Ces calculs sont naturellement bass sur des mesures de pressions et de tempratures rcoltes
tout en sinspirant de lquation (4.9). La formule suivante dcrit le modle du rendement
voqu ci-dessus :
48
Modlisation du systme
1cor
(t ) + c )(d 3 (t ) + e 2 (t ) +
1cor
1cor
1cor
(t ))
(4.10)
o :
a, b, c, d, e sont les constantes des polynmes utiliss,
La figure 4.6 illustre les rsultats de modlisation (en trait pointill) du rendement
isentropique du compresseur par rapport aux valeurs calcules issues de mesures.
Enfin, afin de valider le modle de temprature propos dans lquation (4.9), nous traons
sur la figure 4.7 les points de mesures de temprature en superposition avec les valeurs
calcules laide du modle. Nous constatons galement une trs bonne concordance entre les
rsultats du modle et les valeurs mesures.
0.8
N1
N2
N3
N4
N5
N6
0.7
Rendement isentropique
0.6
0.5
0.4
0.3
0.2
0.1
0.02
0.04
0.06
0.08
0.1
0.12
Dbit corrig (Kg/s)
0.14
0.16
0.18
0.2
Figure 4.6 Rendement isentropique du compresseur illustr par des points de mesures,
chaque rgime de rotation 1=Ni (avec N1< N2 < Ni) est reprsent par un
groupe de points, les courbes pointilles reprsentent les rsultats de modlisation.
49
Modlisation du systme
500
Mesure
Modle
480
460
440
420
400
380
360
340
320
300
0.02
0.04
0.06
0.08
Dbit corrig (Kg/s)
0.1
0.12
0.14
(4.11)
o :
W1 est la puissance absorbe par le compresseur,
C5
Q1
T1 ,T2
50
Modlisation du systme
Ainsi, nous considrons que la diffrence entre le dbit traversant le compresseur et le dbit
traversant lchangeur de refroidissement de lair permet daugmenter ou de rduire la
pression forme aprs le compresseur.
Cette pression peut tre donc calcule en utilisant la loi du gaz parfait comme suit :
P2 (t ) P2 (t 0 ) =
C14T2 (t )
C16
[Q (t ) Q (t )] dt
1
(4.12)
t0
o :
P2 est la pression en sortie du compresseur,
C14
C16
Q1
Q2
T2
La figure 4.8 visualise le bloc fonctionnel liant les variables dentres et de sorties du modle
global du compresseur.
Q1 (Kg/s)
T1 (K)
P2 (Pa)
P1 (Pa)
Q2 (Kg/s)
T2 (K)
w1 (rad/s)
W1 (W)
51
Modlisation du systme
Par consquent, nous pouvons proposer pour lchangeur les mmes modles de dbit et de
pression utiliss pour le filtre air.
Quant la temprature en sortie de lchangeur, elle sera calcule grce un modle
spcifique.
C T (t )
P2 (t ) P3 (t ) = C2 Q2 (t ) 14 2
P2 (t )
Q2 (t ) =
P2 (t )
C14T2 (t )
[P2 (t ) P3 (t )]
C2
avec P2 (t ) P3 (t )
(4.13a)
avec P2 (t ) P3 (t )
(4.13b)
o :
Q 2 est le dbit traversant lchangeur de refroidissement de lair comprim,
C2
P3
La figure 4.9 montre les rsultats de validation de la relation (4.13a) entre le carr du dbit
volumique de lair et la perte de charge travers lchangeur de refroidissement.
12000
Mesure
Modle
10000
P2-P 3 (Pa)
8000
6000
4000
2000
0.01
T2 2
(Q2 C14
P2 )
Figure 4.9 Validation de la perte de charge travers lchangeur de refroidissement de lair comprim,
Ermoy = 5.4%, Ermax = 8.6%.
52
Modlisation du systme
Pour un gaz parfait, la capacit thermique est seulement fonction de la temprature ce qui
permet de reformuler lexpression prcdente comme suit :
2 (t ) =
T2 (t ) T3 (t )
T2 (t ) T1 (t )
(4.14)
o :
2 est lefficacit de lchangeur.
T1 ,T2
T3
(4.15)
Lefficacit de lchangeur est donc llment cl dans le modle de temprature propos cidessus.
Afin de dterminer lefficacit nous proposons un modle polynomial dpendant du dbit
dair traversant lchangeur dune part, et de la vitesse du flux refroidissant (dans ce cas la
vitesse du vhicule assimile la vitesse du vent) dautre part.
Ces deux variables ont t identifies comme tant les grandeurs dinfluence principales qui
agissent sur lefficacit de lchangeur de refroidissement de lair comprim. La structure
suivante sest rvle adapte la modlisation :
2 (t ) = f 2 (Q2 (t ), V (t )) = (aQ2 (t ) + b )(V (t ) + c )
(4.16)
o :
a, b, c sont les constantes des polynmes utiliss,
53
Modlisation du systme
La figure 4.10 illustre les rsultats de modlisation de lefficacit 2 par rapport aux valeurs
calcules de celle-ci laide des mesures et de lquation (4.15). Except pour le cas
correspondant 90 Km/h, les rsultats sont satisfaisants.
1
70 Km/h
90 Km/h
110 Km/h
130 Km/h
Efficacit de l'changeur
0.95
0.9
0.85
0.8
0.75
0.04
0.05
0.06
0.07
0.08
0.09
0.1
Dbit de l'changeur Q2 (Kg/s)
0.11
0.12
0.13
Les rsultats de validation du modle de temprature propos par lquation (4.15) sont
illustrs sur la figure 4.11. On constate une bonne correspondance entre la temprature fournie
par le modle et celle mesure.
350
Mesure
Modle
340
330
320
310
300
290
0.04
0.05
0.06
0.07
0.08
0.09
0.1
Dbit de l'changeur Q2 (Kg/s)
0.11
0.12
0.13
Figure 4.11 Validation du modle de temprature des gaz sortant de lchangeur de refroidissement de
lair comprim, Ermoy = 3.4%, Ermax = 5.2%.
54
Modlisation du systme
C14 T3 (t )
C17
[Q
(t ) Q3 (t )] dt
(4.17)
t0
o :
P3 est la pression aprs lchangeur,
C17
C14
Q2
Q3
dadmission.
T3 est la temprature de lair comprim en sortie de lchangeur.
La figure 4.12 rappelle quelles sont les variables dentres/sorties du modle global de
lchangeur de refroidissement de lair comprim.
Q3 (Kg/s)
P3 (Pa)
P2 (Pa)
T2 (K)
T3 (K)
T1 (K)
Q2 (Kg/s)
V (Km/h)
55
Modlisation du systme
C 4 +1
2
P3 (t )
2 P4 (t ) C 4 P4 (t ) C 4
Athr (t )
C14T3 (t ) C7 P3 (t )
P3 (t )
Q3 (t ) =
+
C
1
2
4
2 2 C 4 1 2 C 4 1
A (t ) P3 (t )
thr
C4 + 1
C14T3 (t ) C7 C4 + 1
C4
Si
P4 (t ) 2 C 4 1
>
P3 (t ) C4 + 1
(4.18)
Sinon
o :
Q3 est le dbit dair traversant le volet dadmission,
Athr
est le produit de la section effective de passage du volet Athr eff par son coefficient de
dcharge Cthr ,
T3
P3 , P4
C4
C14
L1
Nous pouvons distinguer deux cas dcrits par les expressions (4.18) :
C4
P4
P3
2 C 4 1
critique
0.53 ,
C4 + 1
cela signifie que le flux dair na pas atteint la vitesse du son. Dans ce cas, la premire
expression est applique. Ainsi, le dbit du volet dpend la fois de langle douverture et
du rapport des pressions travers le volet.
P4
0.53 ),
P3
56
Modlisation du systme
(4.19)
12
x 10
Mesure
Modle
10
Athr
-2
10
15
20
25
30
35
Angle d'ouverture normalis L1 (%)
40
45
50
Les rsultats de validation du modle de dbit corrig du volet dadmission sont illustrs sur
la figure 4.14 o lon a reprsent ce dbit en fonction du rapport de pression et de langle
douverture. On constate une bonne cohrence entre les valeurs mesures et celles fournies
par le modle.
Mesure
Modle
-4
x 10
Q3
C14T3
P3
50
40
0
30
1
0.9
20
0.8
0.7
10
0.6
0.5
Rapport de pressions P 4/P3
0.4
Position d'ouverture
normalise L1 (%)
57
Modlisation du systme
La figure 4.15 montre un synoptique des entres et sorties du modle global du volet
dadmission.
P3 (Pa)
P4 (Pa)
Q3 (Kg/s)
L1 (%)
T3 (K)
m egr (t )
(4.20)
m air (t ) + m egr (t )
o :
mair
Modlisation du systme
t
t
.
m air (t ) m air (t 0 ) = m air dt = Q3 (t ) Q 4 (t ) 1 egr (t ) dt
t0
t0
t
t
.
=
=
m
(
t
)
m
(
t
)
m
dt
Q6 (t ) Q 4 (t ) egr (t ) dt
egr
egr
egr 0
t0
t0
)]
(4.21)
o :
est le dbit dair sortant de lchangeur de refroidissement et traversant le volet
Q3
dadmission,
Q 4 est le dbit moyen total des gaz entrant dans les cylindres,
est le dbit des gaz traversant le circuit dEGR,
Q6
En cas dinversion de sens de circulation des gaz travers le circuit dEGR, le calcul des
masses peut tre tabli selon les quations suivantes :
t
t
.
m air (t ) m air (t 0 ) = m air dt = Q3 (t ) (Q 4 (t ) + Q6 (t )) 1 egr (t ) dt
t0
t0
t
t
.
t0
t0
)]
(4.22)
Afin de poursuivre le dveloppement du modle de pression nous considrons dabord que les
deux composantes gazeuses sont parfaitement mlanges et ont le mme volume et la mme
temprature dans le collecteur dadmission.
Ainsi, en drivant la loi du gaz parfait par rapport au temps nous pouvons crire :
.
P 4 = P egr + P air
(4.23)
o :
T3 est la temprature de lair comprim refroidi (aprs changeur),
T8
C14
C15
C18
P4
Dans un premier temps nous considrons que les gaz traversant le circuit dEGR circulent
dans le sens direct savoir : du collecteur dchappement vers le collecteur dadmission.
En utilisant les expressions (4.21) dans lquation (4.23) nous obtenons :
59
Modlisation du systme
+
+
P 4 (t ) = 6
C18
C18
C18
C18
.
Q4 (t ) 1 egr (t ) R4 (t )T4 (t )
C18
m (t )C14 T 3 (t )
+ air
=
C18
(4.24)
.
+
+
C18
C18
C18
C18
C18
.
o :
P4 est la pression dans le collecteur dadmission,
T4
R4
est la constante du gaz parfait du mlange gazeux dans le collecteur dadmission dfinie
(4.25)
Pour calculer les variations de temprature des gaz dchappement, dans le collecteur
dadmission, nous utilisons la dfinition de lnergie interne dun gaz parfait ayant une
capacit thermique constante :
U egr (t ) = megr (t )C10T8 (t )
(4.26)
o :
T8 est la temprature des gaz en sortie de lchangeur du circuit EGR,
est la capacit thermique volume constant des gaz dchappement (suppose constante),
C10
En drivant lquation (4.26) avec la prise en compte de la deuxime quation des expressions
(4.21) nous pouvons crire :
.
U egr (t ) = m egr (t )C10T8 (t ) + megr (t )C10 T 8 (t ) = Q6 (t )C10T8 (t ) Q4 (t ) egr (t )Cv 4 (t )T4 (t ) + megr (t )C10 T 8 (t )
(4.27)
o :
Cv4
o :
C 6 est la capacit thermique, suppose constante, volume constant de lair.
60
(4.28)
Modlisation du systme
Nous pouvons galement exprimer, par analogie, la drive de lnergie interne de lair
comme suit :
.
(4.29)
o :
U air
Par ailleurs, le principe de la conservation dnergie pour un systme ouvert nous permet
dcrire les expressions suivantes :
.
.
.
.
.
U egr = H in egr H out egr q egr = Q6 (t )C 9T8 (t ) Q4 (t ) egr (t )Cp 4 (t )T4 (t ) q egr (t )
.
.
.
.
.
U air = H in H out q = Q (t )C T (t ) Q (t ) 1 (t ) Cp (t )T (t ) q (t )
air
air
egr
3
5 3
4
4
4
air
air
(4.30)
o :
.
q egr
reprsente lchange de chaleur entre les gaz dchappement et les parois du collecteur
C 5 m air (t ) + C 9 m egr (t )
m air (t ) + m egr (t )
= C 5 1 egr (t ) + C 9 egr (t )
(4.31)
o :
.
q air reprsente lchange de chaleur entre lair et les parois du collecteur dadmission (ce
1
Q6 (t )T8 (t )C15 egr (t )Q4 (t )T4 (t )R4 (t ) q egr (t )
T 8 (t ) =
megr (t )C10
.
.
air
6
61
(4.32)
Modlisation du systme
P 4 (t ) =
) ]
1
Q6 (t )T8 (t )C15 C 8 + Q3 (t )T3 (t )C14 C 4 Q 4 (t )T4 (t )R 4 (t ) egr (t )C 8 + 1 egr (t ) C 4 f . (t )
q4
C18
(4.33)
avec :
f . (t ) =
q4
C15 .
C .
q egr (t ) + 14 q air (t )
C10
C6
o :
est une fonction exprimant lchange global de chaleur avec les parois du collecteur
f.
q4
dadmission. Cette fonction est positive lorsquil y a une perte de chaleur et vice-versa.
Ainsi, nous pouvons calculer, par intgration, la pression dans le collecteur dadmission :
P4 (t ) P4 (t 0 ) =
1
C18
Q (t )T (t )C
6
15 C 8
t0
(4.34)
o :
4 (t ) = egr (t )C8 + (1 egr (t ))C 4 est le rapport des capacits thermiques spcifiques des gaz
1
C18
Q (t )T (t )C
3
t0
14 C 4
(Q 4 (t ) + Q6 (t ))T4 (t )R 4 (t ) 4 f . (t ) dt
q4
(4.35)
Dans le cadre de cette tude nous allons ngliger lchange de chaleur avec les parois et donc
considrer f . nul.
q4
P4 (t )C18
R4 (t )[mair (t ) + megr (t )]
(4.36)
o :
est le volume du collecteur dadmission,
C18
P4
R4
mair
62
Modlisation du systme
m air (t ) + megr (t )
(4.37)
C18
o :
C18
mair
Enfin, la figure 4.16 rappelle les diffrentes variables intervenant dans le modle global du
collecteur dadmission.
Q4 (Kg/s)
Gamma4
Q3 (Kg/s)
R4
R5
Rho4 (Kg/m3)
T3 (K)
tegr
Gamma5
T4 (K)
Q6 (Kg/s)
P4 (Pa)
T8 (K)
63
Modlisation du systme
n cyl 2 (t )
(4.38)
2 60
o :
ncyl est le nombre de cylindres (4 dans notre cas),
2 (t ) est le rgime moteur exprim en tours/minute.
Par exemple pour un rgime de 3000 tr/min le moteur effectue 100 phases par seconde.
Ainsi, pour convertir par exemple un dbit quelconque Q (Kg/coup) exprim dans le domaine
vnementiel en un dbit Q (Kg/s) exprim dans le domaine temporel, il suffit dappliquer la
formule de conversion suivante :
Q(t ) = Q e (t )
(4.39)
64
Modlisation du systme
Par la suite, nous allons dvelopper des modles statiques dans le but de dterminer le dbit
total moyen des gaz entrant dans les cylindres, le dbit et la temprature des gaz sortant de
ceux-l ainsi que le rapport air/carburant prsent dans la chambre de combustion.
)(
4 (t ) = f 4 ( 4 (t ), 2 (t )) = a 43 (t ) + b 42 (t ) + c 4 (t ) + d 22 (t ) + e 2 (t ) + f
(4.40)
o :
a, b, c, d, e, f sont les constantes des polynmes utiliss.
Connaissant le rendement volumique du moteur, nous pouvons dcrire le dbit moyen total
des gaz entrant dans les cylindres par lquation suivante :
Q 4 (t ) = 4 (t ) e (t ) 4 (t )
C19
4
(4.41)
o :
C19
e est la frquence des phases du moteur exprim en coup/seconde (voir quation 4.38).
La figure 4.17 illustre le rendement volumique du moteur utilis, dtermin laide des
mesures, ainsi quune surface reprsentant le modle polynomial identifi de ce rendement.
65
Modlisation du systme
Modle
Mesure
1.4
Rendement volumtrique
1.3
1.2
1.1
1
0.9
0.8
0.7
5000
4000
3
3000
2.5
2
2000
1.5
1000
Rgime moteur (tr/min)
1
0
0.5
Densit (Kg/m3)
Le modle obtenu du rendement volumique est ainsi utilis dans lquation (4.41) afin de
valider le modle de dbit entrant dans les cylindres. Les rsultats de validation pour les
points de fonctionnement disponibles sont ainsi illustrs sur la figure 4.18. Le dbit
Q4 calcul grce au modle en fonction du rgime du moteur et de la densit des gaz est tout
fait en cohrence avec les mesures de ce dbit.
Durant la phase dchappement, les gaz librs par la combustion sont expulss en dehors des
cylindres vers le collecteur dchappement. Se basant sur le principe de la conservation de la
masse, le dbit total des gaz sortant des cylindres est gal la somme du dbit des gaz entrant
et le dbit du carburant inject dans ceux-l.
Ainsi, la formule (4.42) dcrit le modle du dbit total des gaz sortant des cylindres.
Q5 (t ) = Q 4 (t ) + Q10 (t )
(4.42)
o :
Q 4 , Q5
Q10
sont respectivement les dbits des gaz entrant et sortant des cylindres,
66
Modlisation du systme
Mesure
Modle
0.18
0.16
Dbit Q4 (Kg/s)
0.14
0.12
0.1
0.08
0.06
0.04
0.02
0
5000
4000
3
2.5
3000
2
1.5
2000
1
Rgime moteur (tr/min)
1000
0.5
Densit (Kg/m3)
Figure 4.18 Validation du modle de dbit des gaz entrant dans les cylindres,
Ermoy = 2.9%, Ermax = 4.6%.
Etant donn que le dbit du carburant inject fait partie des lments de combustion, il est
souvent exprim dans le domaine vnementiel moteur. Afin dobtenir lquivalent temporel
de ce dbit et dtablir le calcul dcrit par lquation (4.42), nous utilisons la formule de
conversion (4.39) comme suit :
Q10 (t ) = e (t )
Qinj
(4.43)
10 6
o :
Q inj
T4
T12
67
Modlisation du systme
C9
C3
5 est une sorte de rendement dcrivant la portion dnergie du carburant qui contribue
)(
(4.45)
o :
a, b, c, d, e, f sont les constantes des polynmes utiliss.
Sachant 5 , nous pouvons enfin dcrire le modle de temprature des gaz dchappement en
considrant lquation (4.44) :
T12 (t ) = T4 (t ) +
C 3 Q10 (t ) 5 (t )
C 9 Q5 (t )
(4.46)
La figure 4.19 illustre le rendement 5 du moteur utilis, dtermin laide des mesures, ainsi
quune surface reprsentant le modle polynomial identifi de ce rendement. La prcision
moyenne obtenue peut tre explique par le manque de points de mesure dans certaines zones
du fonctionnement.
Modle
Mesure
0.6
Rendement l'chappement
0.5
0.4
0.3
0.2
0.1
0
-0.1
-0.2
4000
3500
5
3000
4
3
2500
2000
Rgime moteur (tr/min)
-3
x 10
1
1500
68
Modlisation du systme
Le modle obtenu du rendement 5 est utilis dans lquation (4.46) afin de valider le modle
de temprature des gaz dchappement.
Ne disposant pas de moyen de mesure directe du taux dEGR instantan, la validation du
modle de temprature a t effectue avec un taux dEGR nul.
Les rsultats de validation de ce modle pour les points de fonctionnement disponibles sont
ainsi illustrs sur la figure 4.20. Malgr la prcision moyenne obtenue par le modle de
rendement 5 les erreurs de modlisation de la temprature des gaz dchappement restent
acceptables.
Mesure
Modle
1000
900
800
700
600
500
400
4000
3500
5
3000
4
3
2500
2000
Rgime moteur (tr/min)
-3
x 10
1
1500
Figure 4.20 Validation du modle de temprature des gaz sortant de la chambre de combustion,
Ermoy = 8.7%, Ermax = 12.9%.
(4.47)
o :
a , b reprsentent respectivement le nombre datomes de carbone et dhydrogne prsents
69
Modlisation du systme
Pour prendre en compte une combustion complte, lquation (4.47) doit tre quilibre
travers les portions de carburant et dair mises en jeu. Lorsque cet quilibre entre la masse
dair et celle du carburant est atteint, le rapport Air/Carburant est dit stchiomtrique. Ce
dernier est driv de lquation (4.47) de la faon suivante :
mO2 + 3.773m N 2
1 +
a
4
A
=
b
F s
mC + m H
a
)
(4.48)
o :
mO2 , mN 2 , mC , mH sont les masses des diffrentes molcules,
A
F s
Air
dsigne le rapport :
stchiomtrique.
Fuel
A
En situation normale, la valeur de est autour de 14.7, cest--dire que pour quilibrer la
F s
raction de combustion la masse dair doit tre 14.7 fois suprieure celle du carburant.
Dans le domaine de la combustion on parle plus souvent de rapport Air/Carburant relatif
celui stchiomtrique et qui est dfini comme suit :
.
m air
= .
A
m fuel
F s
.
(4.49)
egr
1 egr
(4.50)
.
En remplaant m air et m fuel dans lquation (4.49) puis dans (4.50) par leurs valeurs
correspondantes (cf. figure 4.2) nous obtenons le modle final du rapport Air/Carburant relatif
instantan dans les cylindres :
70
Modlisation du systme
Q 4 (t ) [ 1 egr (t )
total (t ) =
Q10 (t )C 20
1 egr (t )
egr (t )
(4.51)
o :
Q4 est le dbit moyen total des gaz entrant dans les cylindres,
Q10
La figure 4.21 synthtise les entres et sorties du modle global de la chambre de combustion.
w2 (rpm)
Q4 (Kg/s)
Qinj (mg/coup)
Q5 (Kg/s)
Rho4 (Kg/m3)
T12 (K)
tegr
Lambda
T4 (K)
71
Modlisation du systme
m ig (t )
(4.52)
m exh (t ) + m ig (t )
o :
m exh
.
m exh (t ) m exh (t 0 ) = m exh dt = Q5 (t ) (Q6 (t ) + Q7 (t )) 1 ig (t ) dt
t0
t0
t
t
.
m ig (t ) m ig (t 0 ) = m ig dt = (Q6 (t ) + Q7 (t )) ig (t ) dt
t0
t0
)]
(4.53)
o :
Q 5 est le dbit des gaz dchappement expulss par la chambre de combustion,
Q6
Q7
ig est le taux instantan des gaz dadmission prsents dans le collecteur dchappement.
En cas dinversion de sens de circulation des gaz travers le circuit dEGR, le calcul des
masses peut tre accompli en utilisant le systme dquations suivant :
t
t
.
m exh (t ) m exh (t 0 ) = m exh dt = Q5 (t ) Q 7 (t ) 1 ig (t ) dt
t0
t0
t
t
.
m ig (t ) m ig (t 0 ) = m ig dt = Q6 (t ) Q7 (t ) ig (t ) dt
t0
t0
)]
(4.54)
En suivant la dmarche utilise pour le collecteur dadmission, nous obtenons les quations
suivantes pour la pression dans le collecteur dchappement :
72
Modlisation du systme
1
P5 (t ) P5 (t0 ) =
C21
Q (t )T
5
12
q5
t0
(4.55)
avec :
R 4 m ig (t ) + C15 m exh (t )
m exh (t ) + m ig (t )
):
= R 4 ig (t ) + C15 1 ig (t )
C 21
T12
Q5
Q6
Q7
f.
est une fonction exprimant lchange global de chaleur avec les parois du collecteur
q5
dchappement. Cette fonction est positive lorsquil y a une perte de chaleur et vice-versa.
En cas dinversion concernant les gaz du circuit dEGR, lquation suivante remplace
lquation (4.55) pour le calcul de la pression dans ce collecteur comme suit :
P5 (t ) P5 (t 0 ) =
1
C 21
Q (t )T
5
12
q5
t0
(4.56)
o :
R 4 est la constante du gaz parfait du mlange gazeux dans le collecteur dadmission,
T11
est la temprature des gaz prleve entre le conduit de retour EGR et le collecteur
dchappement.
Dans le cadre de cette tude nous allons ngliger galement lchange de chaleur avec les
parois et donc considrer f . nul.
q5
Modlisation du systme
T5 (t ) =
P5 (t )C 21
R5 (t )[ m exh (t ) + mig (t )]
(4.57)
o :
P5 est la pression des gaz dans le collecteur dchappement,
R5
C 21
m exh
T12 (K)
P5 (Pa)
Q5 (Kg/s)
Q7 (Kg/s)
T5 (K)
Q6 (Kg/s)
T11 (K)
R5
R4
Gamma5
Gamma4
74
Modlisation du systme
Cette rduction de temprature se traduit naturellement par une rduction dmission des
oxydes dazote (Nox) et par une augmentation dmission des particules.
Sous certaines conditions, telles quune coupure brutale du carburant inject dans la chambre
de combustion, une circulation des gaz dans le sens inverse de celui voqu ci-dessus peut se
produire dans le circuit durant une courte priode.
Notre but ici est de dterminer le dbit ainsi que la temprature des gaz reconduits quel que
soit leur sens de circulation laide de modles statiques. Pour ce faire, nous dcomposons le
circuit dEGR en trois parties selon le sens direct des gaz comme suit :
- conduit de retour
- lectrovanne
- changeur de refroidissement
(4.58)
o :
6 est lefficacit de refroidissement du conduit,
T9
T5
T10
Afin de dterminer lefficacit nous proposons un modle polynomial simple utilisant le dbit
des gaz travers le conduit :
6 (t ) = f6 (Q6 (t )) = aQ6 3 (t ) + bQ6 2 (t ) + cQ6 (t ) + d
(4.59)
o :
a, b, c, d
Q6
75
Modlisation du systme
La figure 4.23 illustre les rsultats de modlisation de lefficacit 6 par rapport aux valeurs
calcules de celle-ci laide des mesures et de lquation (4.58). Nous constatons une bonne
cohrence entre le modle et les mesures.
0.8
Mesure
Modle
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0.002
0.004
0.006
0.008
0.01
0.012
0.014
Dbit des gaz traversant le conduit de retour Q6 (Kg/s)
0.016
0.018
Les rsultats de validation du modle de temprature propos dans lquation (4.58) sont
illustrs sur la figure 4.24. Une prcision plutt satisfaisante est obtenue.
Mesure
Modle
750
700
650
600
550
500
450
400
800
750
0.02
700
0.015
650
0.01
600
0.005
550
Figure 4.24 Validation du modle de temprature des gaz sortant du conduit de retour dEGR,
Ermoy = 5.7%, Ermax = 8.6%.
76
Modlisation du systme
Dans le cas o le sens de circulation des gaz est invers dans le circuit, la mme structure
dquations peut tre utilise simplement en prenant soin de rattribuer les variables dentre
et de sortie. Ainsi, lquation (4.58) peut tre reformule dans ce cas comme suit :
T11 (t ) = T10 (t ) 6 (t )[T10 (t ) T9 (t )]
(4.60)
o :
T11
est la temprature des gaz prleve entre le conduit de retour dEGR et le collecteur
dchappement, et qui sera utilise comme temprature finale des gaz provenant du circuit
dEGR.
T10 reprsente dans ce cas la temprature des gaz en sortie de llectrovanne dEGR et en
entre du conduit de retour.
La figure 4.25 illustre le diagramme synoptique du modle ddi au conduit de retour dEGR
dans les deux cas.
Q6 (Kg/s)
Q6 (Kg/s)
T5 (K)
T10 (K)
T10 (K)
T9 (K)
T11 (K)
T9 (K)
(A)
(B)
4.2.8.2 Electrovanne
Il sagit de lactionneur principal du circuit dEGR qui permet de rgler le dbit souhait des
gaz dchappement recycls vers ladmission. Cet actionneur est quip dun capteur de
position permettant la surveillance et le contrle de sa leve par le calculateur lectronique.
Ainsi, nous avons utilis ce signal de position en plus des mesures de pression et dbit pour
construire les modles que nous allons prsenter par la suite.
77
Modlisation du systme
5 +1
2
P5 (t )
2 P4 (t ) 5 P4 (t ) 5
Aegr (t )
P (t )
R5T10 (t ) 5 1 P5 (t )
5
Q6 (t ) =
5 +1
2
2 2 5 1 2 5 1
A (t ) P5 (t )
egr
1
1
+
+
(
)
R
T
t
5
5
5
5 10
Si
P4 2 5 1
>
P5 5 + 1
(4.61)
Sinon
o :
Q6 est le dbit des gaz traversant llectrovanne dEGR,
Aegr est le produit de la section effective de passage de llectrovanne Aegreff par son coefficient
de dcharge Cegr,
P4 est la pression des gaz dans le collecteur dadmission,
P5
T10
R5
5 est le rapport des capacits thermiques spcifiques des gaz dans le collecteur
dchappement,
L2 est la position douverture normalise de llectrovanne.
Comme nous lavons propos pour le volet dadmission, un modle polynomial peut tre
galement utilis pour dcrire la fonction Aegr comme suit :
Aegr (t ) = f Aegr (L2 (t )) = aL2 3 (t ) + bL2 2 (t ) + cL2 (t ) + d
(4.62)
o :
a, b, c, d sont les constantes du polynme utilis.
La figure 4.26 montre les rsultats de modlisation de cette fonction par rapport aux valeurs
mesures de Aegr dduites des quations (4.61). Le modle obtenu pour la fonction Aegr est
utilis par la suite dans le modle de dbit corrig driv des quations (4.61). Les rsultats de
validation de ce modle sont illustrs sur la figure 4.27.Les rsultats obtenus nous semblent
tout fait satisfaisants.
78
Modlisation du systme
-4
x 10
Aegr
Mesure
Modle
10
20
30
40
50
60
70
Position d'ouverture normalise L2 (%)
80
90
100
Mesure
Modle
-4
x 10
1
0.6
Q6
C15T10
P5
0.8
0.4
100
0.2
80
0
60
1
0.9
40
0.8
0.7
0.6
20
0.5
0.4
Position d'ouverture
normalise L2 (%)
Figure 4.27 Validation du modle de dbit corrig des gaz traversant llectrovanne dEGR,
Ermoy = 3.2%, Ermax = 6.9%.
79
Modlisation du systme
En cas dinversion de sens de circulation des gaz dans le circuit, la pression dans le collecteur
dadmission est suprieure celle dans le collecteur dchappement. Ainsi, les quations
(4.61) doivent tre rcrites en respectant ce cas de figure comme suit :
4 +1
2
P4 (t )
2 P5 (t ) 4 P5 (t ) 4
Aegr (t )
R 4 T6 (t ) 4 1 P4 (t )
P4 (t )
4 +1
2
Q6 (t ) =
2 2 4 1 2 4 1
A (t ) P4 (t )
egr
4 1 4 + 1
4 + 1
(
)
R
T
t
4
6
A = f (L ) = A
Aegr
2
egr eff (L 2 ) C egr (L 2 )
egr
Si
P5 2 4 1
>
P4 4 + 1
(4.63)
Sinon
o :
T6 est la temprature des gaz en sortie de lchangeur de refroidissement de lEGR et en
entre de llectrovanne (selon le sens inverse des gaz),
R 4 est la constante du gaz parfait du mlange gazeux dans le collecteur dadmission,
4 est le rapport des capacits thermiques spcifiques des gaz dans le collecteur dadmission.
5 1
5
(4.64)
o :
T10 , T6
EGR,
P5 , P4
5 est le rapport des capacits thermiques spcifiques des gaz dans le collecteur
dchappement.
En cas dinversion de sens de circulation des gaz dans le circuit dEGR, lquation (4.64)
devient :
P (t )
T10 (t ) = T6 5
P4 (t )
4 1
4
(4.65)
80
Modlisation du systme
o :
T6 , T10
EGR,
P4 , P5
4 est le rapport des capacits thermiques spcifiques des gaz dans le collecteur dadmission.
La figure 4.28 illustre le diagramme synoptique du modle ddi llectrovanne dEGR dans
les deux cas.
P5 (Pa)
P5 (Pa)
Q6 (Kg/s)
T10 (K)
L2 (%)
L2 (%)
P4 (Pa)
P4 (Pa)
R5
Q6 (Kg/s)
T6 (K)
R4
T6 (K)
Gamma5
T10 (K)
Gamma4
(A)
(B)
(4.66)
o :
81
Modlisation du systme
est la temprature des gaz prleve entre lchangeur et llectrovanne du circuit dEGR,
T9
Afin de dterminer lefficacit 7 , nous proposons un modle polynomial simple qui utilise le
dbit des gaz traversant lchangeur :
7 (t ) = f7 (Q6 (t )) = aQ6 2 (t ) + bQ6 (t ) + c
(4.67)
o :
a, b, c sont les constantes du polynme utilis.
Q6
La figure 4.29 illustre les rsultats de modlisation de lefficacit 7 par rapport aux valeurs
calcules de celle-ci laide des mesures et de lquation (4.66). Nous constatons une bonne
correspondance entre le modle et les valeurs mesures.
1.15
Mesure
Modle
1.1
1.05
Efficacit de l'changeur
0.95
0.9
0.85
0.8
0.75
0.7
0.65
0.002
0.004
0.006
0.008
0.01
0.012
Dbit des gaz traversant l'changeur Q6 (Kg/s)
0.014
0.016
0.018
Les rsultats de validation du modle de temprature propos dans lquation (4.66) sont
galement illustrs sur la figure 4.30. Ils montrent une bonne concordance entre le modle et
les mesures.
82
Modlisation du systme
Dans le cas dun flux de gaz invers dans le circuit dEGR, lquation (4.66) est remplace
par celle correspondant ce cas de figure comme suit :
T6 (t ) = T4 (t ) 7 (t )[T4 (t ) T9 (t )]
(4.68)
o :
T4 est la temprature des gaz dans le collecteur dadmission et dans ce cas en entre de
lchangeur de refroidissement du circuit dEGR.
Mesure
Modle
480
460
440
420
400
380
360
340
800
700
0.02
0.015
600
0.01
500
Temprature aprs vanne EGR T6 (K)
0.005
400
83
Modlisation du systme
Q6 (Kg/s)
Q6 (Kg/s)
T6 (K)
T8 (K)
T4 (K)
T9 (K)
T6 (K)
T9 (K)
(A)
(B)
P4 (Pa)
Q6 (Kg/s)
T4 (K)
P5 (Pa)
T5 (K)
R5
T8 (K)
Gamma5
R4
Gamma4
L2 (%)
T11 (K)
T9 (K)
4.2.9 Turbine
Lnergie cinmatique des gaz dchappement est utilise par la turbine afin dentraner le
compresseur. En gnral, la turbine fonctionne de faon similaire au compresseur mais dans le
sens oppos.
Ainsi :
un phnomne de dtente est constat,
un abaissement de temprature est caus,
une puissance utile est produite et une rotation de larbre turbocompresseur est obtenue.
Le principe dextraction dnergie cinmatique dans la turbine consiste restreindre le flux
des gaz traversant cette dernire grce ses ailettes. Ainsi, du fait de la diffrence de pression
84
Modlisation du systme
cre, une acclration du flux des gaz est produite entrainant par la mme occasion
lacclration de lensemble ailettes-arbre turbocompresseur.
Dans le cadre de cette tude une turbine dite gomtrie variable est utilise. Elle offre la
possibilit de faire varier langle dincidence de ses ailettes par rapport au flux des gaz afin de
faire varier le rendement de la turbine et contrler ainsi la vitesse de rotation du
turbocompresseur.
En rsum, afin de dcrire le fonctionnement de cette partie du systme nous devons
modliser les grandeurs suivantes :
le dbit des gaz traversant la turbine,
la temprature et la pression des gaz en sortie de la turbine,
la puissance fournie par la turbine et sa vitesse de rotation.
P5 (t ) Tref
Pref T5 (t )
(4.69)
o :
P5 , T5 sont respectivement la pression et la temprature des gaz dans le collecteur
dchappement,
Q7 cor est le dbit corrig des gaz traversant la turbine,
Pref , Tref sont respectivement la pression et la temprature atmosphrique de rfrence.
Un modle polynomial est ici propos pour dterminer le dbit corrig des gaz traversant la
turbine :
3
2
P5 (t )
P5 (t )
P5 (t )
P (t )
+ b
+ c 5 + d L3 2 (t ) + eL3 (t ) + f
, L3 (t ) = a
Q7 cor (t ) = f Q7 cor
P7 (t )
P7 (t )
P7 (t )
P7 (t )
o :
P7 est la pression en sortie de la turbine,
L3
85
(4.70)
Modlisation du systme
Les rsultats de validation du modle dcrit par lquation (4.70) sont illustrs sur la figure
4.33. Nous constatons une trs bonne corrlation entre les valeurs mesures et celles estimes
par le modle.
Mesure
Modle
0.1
0.09
0.08
0.07
0.06
0.05
0.04
2.5
80
60
40
1.5
20
Rapport de pressions P5 / P7
Figure 4.33 Validation du modle de dbit corrig des gaz traversant la turbine,
Ermoy = 3.2%, Ermax = 4.9%.
C11
(4.71)
o :
P7 est la pression en sortie de la turbine,
P5 , T5
C8 1
C8
86
Modlisation du systme
A cause des diffrentes pertes, la turbine possde galement son propre rendement dcrivant
la relation entre labaissement de la temprature thorique et celui obtenu en pratique comme
le montre la formule suivante :
3 (t ) =
T5 (t ) T7 (t )
T5 (t ) T7isent (t )
(4.72)
o :
T7 , T7isent
sortie de la turbine,
T5 est la temprature des gaz dans le collecteur dchappement.
En remplaant T7isent par sa valeur de (4.71) et en isolant le terme T7 dans lquation (4.72),
nous obtenons le modle de la temprature en sortie de la turbine :
P (t ) C11
T7 (t ) = T5 (t ) 3 (t )T5 (t )1 5
P7 (t )
(4.73)
87
Modlisation du systme
3
L3 = 0%
L3 = 11%
L3 = 23%
L3 = 35%
2.5
Rendement de la turbine
94 krpm
89 krpm
2
76 krpm
107 krpm
1.5
108 krpm
94 krpm
114 krpm
106 krpm
114 krpm
144 krpm
1
90 krpm
115 krpm
0.5
1.5
2.5
Rapport de pressions P5 / P7
En revanche, leffet de lchange de chaleur devient ngligeable dans les zones o les rapports
de pressions et les rgimes de rotation de la turbine sont plus levs. Ainsi, les valeurs de
rendement obtenues convergent dans ce cas vers les valeurs fournies par le fabricant du
turbocompresseur et souvent releves dans les zones de fonctionnement mentionnes
prcdemment.
Pour des raisons de clart, la dpendance du rendement de la turbine au rgime de rotation de
celle-ci ( rapport de pressions
P5
P7
la figure 4.34. Cependant, le rendement de la turbine dpend bel et bien de trois lments : le
rapport de pressions travers la turbine, son rgime de rotation et enfin langle douverture de
ses ailettes. De plus, leffet dchange de chaleur augmente ltendue dchelle du rendement
de faon considrable (jusqu 300% pour des rapports de pressions proches de 1.15).
Tous ces aspects voqus ci-dessus, rajouts au manque de donnes exprimentales (du fait de
la difficult de ralisation de ce type dessais sur vhicule) rendent inenvisageable la
modlisation par rgression non-linaire directe des donnes. Cest pourquoi nous proposons
ici une mthode de modlisation base sur la sparation du rendement arodynamique de la
turbine de celui rajout par leffet dchange de chaleur voqu prcdemment. Ainsi, le
rendement arodynamique sera modlis en utilisant une approche physique base sur le
rapport entre la vitesse des ailettes et celle des gaz traversant la turbine.
88
Modlisation du systme
(4.74)
Cette quation reflte le fait que la temprature des gaz sortant de la turbine est gale celle
obtenue en utilisant uniquement son rendement arodynamique ( T7 aero ), moins la perte de
temprature lie lchange de chaleur ( T7 ).
En isolant le terme 3 dans lquation (4.73) tout en y remplaant T7 par sa valeur de
lquation (4.74) nous obtenons :
3 (t ) =
T7 aero (t )
T5 (t )
P (t )
1 5
P7 (t )
C11
T7 (t )
+
T5 (t )
P (t )
1 5
P7 (t )
C11
= 3aero (t ) + 3 (t )
(4.75)
Modlisation du systme
C u (t ) =
2C 9 Tref
(4.76)
P (t ) C11
1 5
P7 (t )
avec :
~1cor (t ) = 1 (t )
Tref
(4.77)
T5 (t )
o :
~1cor est le rgime de rotation corrig de la turbine,
Tref ,T5 sont respectivement la temprature atmosphrique de rfrence et celle des gaz dans le
C9
C 23
Le rendement arodynamique de la turbine peut ainsi tre dfini comme une fonction
quadratique de Cu :
2
C u (t )
(
)
C
t
u
3aero (t ) = 3 max (t )2
C uopt (t ) C uopt (t )
(4.78)
o :
3 max , Cuopt sont respectivement le rendement arodynamique maximum et le rapport optimal
de vitesse des ailettes. Pour ces deux grandeurs nous proposons galement deux modles
polynomiaux du type bote grise dcrits par les formules suivantes :
C uopt (t ) = f C
uopt
o :
L3 est langle douverture normalis des ailettes,
90
(4.79)
Modlisation du systme
a, b, c, d, e, f, g, h sont les coefficients des polynmes utiliss et qui sont identifies en utilisant
0.9
1.3
0.85
1.2
0.8
1.1
Rapport optimal de vitesse des ailettes
0.75
0.7
0.65
0.6
0.55
0.9
0.8
0.7
0.6
0.5
0.5
0.45
0.4
20
40
60
80
Angle d'ouverture normalis L3 (%)
0.4
100
20
40
60
80
Angle d'ouverture normalis L3 (%)
100
T5 (t )
P (t )
1 5
P7 (t )
(4.80)
C11
91
Modlisation du systme
(4.81)
o :
8 est lefficacit de refroidissement de la structure mtallique du turbocompresseur,
T2
(4.82)
o :
a, b, c sont les constantes du modle dtermines grce aux donnes exprimentales,
Q7 cor
La figure 4.36 visualise quelques donnes exprimentales de lefficacit 8 ainsi que des
valeurs de celle-ci issues de son modle donn par lquation (4.82).
92
Modlisation du systme
0.18
Mesure
Modle
0.16
0.14
0.12
0.1
0.08
0.06
0.04
0.02
0.045
0.05
0.055
0.06
0.065
0.07
Dbit corrig des gaz Q7cor (Kg/s)
0.075
0.08
0.085
Enfin, les figures 4.37 et 4.38 montrent respectivement quelques rsultats de validation des
modles de rendement global donn par lquation (4.75) et de temprature des gaz en sortie
de la turbine donn par lquation (4.73).
3
L3 = 11%
L3 = 23%
L3 = 35%
Rendement de la turbine
2.5
1.5
0.5
1.1
1.2
1.3
1.4
1.5
1.6
Rapport de pressions P5 / P7
1.7
1.8
1.9
93
Modlisation du systme
900
L3 = 11 %
L3 = 23 %
850
L3 = 35 %
800
750
700
650
600
550
0.7
0.8
0.9
1.8
1.6
1.1
1.4
1.2
1.2
1.3
1
1.4
Rapport de pressions P5 / P7
x 10
C15 T7 (t )
C 22
[Q
(t ) Q9 (t )] dt
(4.83)
t0
o :
C 22
C15
T7
Q7 , Q9
dchappement.
94
Modlisation du systme
(4.84)
o :
C 9 est la capacit thermique des gaz dchappement pression constante,
est le dbit massique des gaz traversant la turbine,
Q7
T5 , T7
C12 1 (t 0 )
[W
(t ) W1 (t )] dt
(4.85)
t0
o :
C12
W1 , W 2
turbine.
La figure 4.39 est un diagramme synoptique du modle global de la turbine rappelant les
diffrentes variables dentres et de sorties.
95
Modlisation du systme
W1 (W)
w1 (rad/s)
T2 (K)
Q7 (Kg/s)
L3 (%)
T5 (K)
P7 (Pa)
P5 (Pa)
T7 (K)
Q9 (Kg/s)
C 8 +1
2
P7 (t )
2 P0 (t ) C8 P0 (t ) C8
C13
C15T7 (t ) C11 P7 (t )
P7 (t )
Q9 (t ) =
C 8 +1
2
P7 (t )
2 2 C8 1 2 C8 1
C
13 C T (t ) C C + 1
C8 + 1
11 8
15
7
96
C8
Si
P0 2 C8 1
>
P7 C8 + 1
(4.86)
Sinon
Modlisation du systme
o :
C13
C11 =
C8 1
.
C8
La figure 4.40 illustre la validation de la relation linaire (dont la pente est gale C13 )
dcrite par lquation (4.86) entre le terme contenant la racine carre et le dbit corrig :
Q9
C15T7
P7
.
-4
x 10
Q9
C15 T7
P7
0.1
0.2
0.3
f(P0/P7)
0.4
0.5
0.6
2
f (P0 / P7 ) =
C11
C8 +1
2
P0 (t ) C8 P0 (t ) C8
P (t )
P7 (t )
Enfin, la figure 4.41 illustre le modle global du systme dchappement par un diagramme
synoptique des entres/sorties.
97
Modlisation du systme
T7 (K)
P7 (Pa)
Q9 (Kg/s)
P0 (Pa)
L1 (t ) L1 (t 0 ) = D1 (t )1 e C24
t t
L 2 (t ) L 2 (t 0 ) = D 2 (t )1 e C25
t t
L3 (t ) L3 (t 0 ) = D3 (t )1 e C26
(4.87)
o :
D1 , D 2 , D3
98
Modlisation du systme
D1 (%)
L1(%)
D2 (%)
L2 (%)
D3 (%)
L3 (%)
99
Modlisation du systme
100
Chapitre V
101
5.2
5.3
Simulation............................................................................................................................................ 106
5.3.1
Simulation pas fixe........................................................................................................................ 106
5.3.2
Choix des intgrateurs temps discret............................................................................................. 107
5.4
Recalage............................................................................................................................................... 107
5.4.1
Recalage statique ............................................................................................................................. 108
5.4.2
Recalage dynamique ........................................................................................................................ 109
5.5
Validation ............................................................................................................................................ 111
5.5.1
Cycle de conduite extra-urbain (Extra Urban Driving Cycle) ......................................................... 112
5.5.2
Rsultats exprimentaux .................................................................................................................. 113
102
Nous prsenterons dans ce chapitre limplmentation du modle complet de la boucle des gaz
du moteur OM646 dont nous avons modlis les composants dans le chapitre prcdent.
Nous voquerons ensuite trois aspects primordiaux pour tout modle dvelopp en vue de
lapplication en temps rel :
La simulation qui est ncessaire pour vrifier la cohrence des relations entres/sorties de
chaque sous-systme, le comportement global des modles (temps de rponse, stabilit,
sens de variation, sensibilit) et donc la possibilit dimplmentation des modles en
temps rel.
Le recalage des constantes du modle qui permet, grce des mesures, la ressemblance
des comportements statique et dynamique du modle ceux du systme physique.
La validation du modle complet qui, aprs la simulation et le recalage des paramtres,
permet dvaluer la prcision obtenue compte tenu des erreurs de modlisation en
comparant les rsultats de simulation aux donnes de mesure et de donner les limites du
modle.
Inlet
air(Kg/s)
flow
Qcomp
Ambient
Pamb (Pa)pressure
Pf ilt (Pa)
After filter pressure
Vehicle
speed
Vwind (m/s)
Engine
speed
Neng (rpm)
Tinter (K)
After intercooler temperature
Fuel
flow
Qinj (mg/st)
Intake manifold
pressure
Pintake
(Pa)
Coolant
Tcool (K) temperature
Throttle setpoint
Throttle
Legr (%) actuator lift
Exhaust manifold temperature
Texhaust (K)
EGR setpoint
EGR
actuator lift
Lv gt (%)
VGT setpoint
VGT
actuator lift
Lthr (%)
Exhaust manifold
pressure
Pexhaust
(Pa)
Electric actuators
Figure 5.1 Vue globale du modle complet de la boucle des gaz du moteur diesel OM646
avec ses actionneurs lectriques dans lenvironnement Matlab/simulink.
103
Figure 5.2 Architecture du modle complet de la boucle des gaz implment dans lenvironnement Matlab/simulink.
104
Dfinition
Temprature dair aprs filtre air
Pression atmosphrique
Vitesse du vhicule
Rgime du moteur
Nature
Entre 1
Entre 2
Entre 3
Entre 4
Q10
T9
L1
L2
L3
Q1
P1
T3
P4
T5
P5
Entre 5
Entre 6
Entre 7
Entre 8
Entre 9
Sortie 1
Sortie 2
Sortie 3
Sortie 4
Sortie 5
Sortie 6
105
Variable
D1
D2
D3
L1
L2
L3
Dfinition
Consigne dangle douverture normalis du volet dadmission
Consigne dangle douverture normalis de llectrovanne EGR
Consigne dangle douverture normalis des aubes de la turbine
Angle douverture normalis du volet dadmission
Angle douverture normalis de llectrovanne EGR
Angle douverture normalis des aubes de la turbine
Nature
Entre 1
Entre 2
Entre 3
Sortie 1
Sortie 2
Sortie 3
5.3 Simulation
En gnral, lorsque la simulation dun modle est envisage, il est prfrable dutiliser un
intgrateur pas de calcul variable. Ainsi, la progression de la simulation peut sadapter la
dynamique des variables simules ce qui peut tre plus conomique quun pas fixe de petite
amplitude compatible avec lensemble des dynamiques.
Cependant, lorsque le modle est utiliser par la suite dans une application automobile en
temps rel, il doit tre simul dans ce cas, pour des raisons de simplicit, avec un intgrateur
pas fixe.
Dans les paragraphes suivants, nous exposerons les raisons pour lesquelles nous avons utilis
ce type dintgrateurs et de la faon dont les simulations ont t faites.
106
maximum du pas fixe pour un modle non-linaire comme celui de la boucle des gaz
dvelopp ici. Le seul moyen de dterminer cette largeur de pas maximum est de tester et
simuler le modle jusquaux limites de stabilit.
Plusieurs essais de simulation ont t effectus durant la phase de rglage des constantes du
modle (que nous dvelopperons dans la section 5.4 de ce chapitre) en utilisant des donnes
de mesures exprimentales. Ces essais ont permis de dfinir la largeur maximum du pas fixe
de simulation (4 ms) au-del de laquelle la simulation donnait des rsultats conduisant des
instabilits.
1
s
h(z + 1)
dans le domaine discret.
2(z 1)
Cette mthode semble intressante pour des dynamiques dites raides (stiff) o le rapport
de constantes de temps entre deux sous-systmes du modle est grande (facteur de 100).
En revanche, la mthode trapzodale cre des boucles algbriques dans le modle qui
ncessitent une puissance de calcul importante.
Cest la raison pour laquelle nous avons plutt opt pour la mthode dEuler rcursive o
lintgrateur
1
s
h
z 1
Aprs avoir discrtis tous les tats continus du modle, notamment les pressions dans les
volumes de la boucle des gaz ainsi que la vitesse du turbocompresseur, nous avons utilis tout
simplement lintgrateur discret (pas dtat continu) de Simulink.
5.4 Recalage
Le but de cette phase de dveloppement du modle est de faire correspondre au mieux le
comportement du modle simul celui du systme rel, savoir la boucle des gaz du
moteur.
Pour ce faire, nous avons utilis des donnes exprimentales provenant de diffrents types
dessais effectus laide du vhicule de recherche prsent au cours du chapitre III.
Une partie de ces donnes, correspondant aux variables dentre du modle, a t injecte
dans le simulateur dj prsent pour tablir la simulation.
107
Quant la deuxime partie des donnes, elle a servi valuer lerreur entre les valeurs des
sorties simules et celles issues du systme rel en vue du recalage du modle.
Description
Coefficient de dcharge du filtre air
Nature
Scalaire
Matrice 64x64
f1
Matrice 32x32
C2
Scalaire
f 2
Matrice 16x16
f Athr
Vecteur 1x16
f 4
Matrice 32x32
f5
Matrice 32x32
f6
Vecteur 1x16
f Aegr
Vecteur 1x16
f 7
Vecteur 1x16
f Q7 cor
Matrice 32x32
fC
uopt
Vecteur 1x10
f3 max
Vecteur 1x10
f8
Vecteur 1x16
C13
Scalaire
f Q1cor
Tableau 5.3 Paramtres statiques du modle de la boucle des gaz du moteur diesel.
108
Description
Volume constitu par le filtre air et son raccord au compresseur
Volume situ entre le compresseur et son changeur de refroidissement
Volume situ entre lchangeur de refroidissement et le volet dadmission
Volume du collecteur dadmission
Volume du collecteur dchappement
Volume situ entre la turbine et le systme d'chappement
Moment dinertie de larbre turbocompresseur
Constante de temps de la rponse du volet dadmission
Constante de temps de la rponse de la vanne EGR
Constante de temps de la rponse des ailettes de la turbine
Tableau 5.4 Paramtres dynamiques du modle de la boucle des gaz du moteur OM646.
Pour identifier les constantes de temps des actionneurs lectriques, nous avons appliqu des
chelons de consigne de position chacun de ces actionneurs.
Diffrents points de fonctionnement et amplitudes dchelons ont t utiliss durant les tests
des actionneurs. Les donnes exprimentales dentres/sorties rcoltes ont permis
didentifier les valeurs des constantes de temps des trois actionneurs tests.
109
Ensuite, tant donn que le dbit dair dadmission affecte pratiquement toutes les parties
mcaniques et arodynamiques du modle de la boucle des gaz, nous avons dcid de
procder des essais du type rponse un chelon de dbit afin de rvler les diffrentes
dynamiques contenues dans cette boucle.
Pour raliser des chelons de dbit dair les plus francs possibles nous avons appliqu des
chelons de consigne de position au volet dadmission compte tenu de la rapidit daction de
ce dernier par rapport aux autres possibilits de sollicitation (chelons de rgime, dbit de
carburant ou de position dailettes de la turbine).
Enfin, nous avons utilis la boite outils Simulink Parameter Estimation de Mathworks afin
dajuster automatiquement les valeurs des paramtres dynamiques du modle complet de la
boucle des gaz.
Cette boite outils permet simplement de recaler la rponse dun modle de simulation (celui
de la boucle des gaz dans notre cas) celle du vrai systme physique en utilisant diffrents
algorithmes doptimisation (gradient, moindres carrs non-linaire).
Notons que des valeurs typiques des paramtres dynamiques du modle (provenant du
constructeur) ainsi que des donnes exprimentales dentres/sorties du modle ont t
injectes initialement dans la boite outils avant de procder au recalage.
Diffrents points de fonctionnement du moteur (rgime, dbit de carburant) ont t utiliss
pour rcolter les donnes exprimentales ncessaires au recalage dynamique.
Les valeurs des paramtres dynamiques obtenues ont t pratiquement les mmes quel que
soit le point de fonctionnement utilis.
Ce rsultat peut tre expliqu par le fait que les paramtres recherchs, de par leur nature
(volumes, moment dinertie), ne dpendent que de larchitecture matrielle de la boucle des
gaz dans le moteur.
La figure 5.1 illustre le recalage de la rponse du modle de simulation de la boucle des gaz
travers un extrait de quelques signaux mesurs mettant en valeur les proprits dynamiques du
systme.
Pour des raisons lies au fonctionnement (risque de calage moteur), les chelons de position
de volet dadmission applicables ne permettaient pas dobtenir des variations de dbit dair
suprieures 15 g/s.
Par consquent, les variations de pressions P1 , P3 et P7 obtenues ntaient pas suffisamment
importantes pour identifier et rajuster les paramtres dynamiques viss ( C 0 , C17 , C 22 ).
Ces constantes ont donc t fixes leurs valeurs typiques durant le recalage dynamique
automatique du modle.
110
L1 (%), Q1 (g/s)
100
L1
80
60
40
Q1
20
P4 (Pa)
x 10
1.2
Mesure
Modle
1
0.8
P5 (Pa)
x 10
1.6
Mesure
Modle
1.4
1.2
w1 (Krpm)
90
Mesure
Modle
80
70
60
20
25
30
35
40
45
50
Temps (s)
Figure 5.1 Extrait des signaux mesurs ayant servi la validation du recalage dynamique
du modle de la boucle des gaz du moteur OM646.
5.5 Validation
Pour valider les diffrents sous-modles de la boucle des gaz du moteur nous avons compar
leurs valeurs de sorties simules avec celles mesures.
Ainsi, nous avons prsent dans le chapitre IV des rsultats de validation pour les sousmodles statiques sous forme graphique au fur et mesure du dveloppement de ceux-ci.
Nous avons montr galement (figure 5.1) un extrait des signaux disponibles illustrant les
rsultats de recalage des sous-modles dynamiques du systme et qui reprsentent une sorte
de validation de ceux-l.
Nous tenons rappeler que les donnes exprimentales ayant servi au recalage statique ou
dynamique du modle sont diffrentes de celles utilises pour la validation de celui-l.
Afin de combiner validation statique et dynamique et dvaluer les prcisions des sorties
globales du modle (dfinies dans les tableaux 5.1, 5.2 et utilises par la suite pour le
diagnostic du systme), nous avons effectu diffrents types dessai de circulation avec le
vhicule de recherche : en ville, sur autoroute et mixte.
Nous avons choisi de prsenter par la suite un exemple de validation du modle complet
effectue durant une partie dun cycle de conduite standard (EUDC).
111
En ce qui nous concerne, nous nous sommes focaliss sur le dernier cycle du test NEDC
savoir le cycle extra-urbain EUDC (Extra Urban Driving Cycle) pour sa bonne couverture en
termes de points de fonctionnement du systme (vitesse du vhicule, rgime du moteur,
pressions dadmission et dchappement, dbits dair et de carburant).
De plus, ce type de tests (NEDC) tait souvent effectu dans les laboratoires de lentreprise
pour optimiser la mise au point du contrle du moteur en fonction des rsultats dmissions.
Cela nous permettait donc de nous procurer des donnes exprimentales reprsentatives de
ces tests sans avoir les effectuer spcialement pour le dveloppement de notre tude.
112
V (Km/h)
150
100
50
w2 (rpm)
3000
2000
Q10 (mg/coup)
1000
40
20
0
D1 (%)
100
50
0
D2 (%)
100
50
0
D3 (%)
20
10
0
1050
1060
1070
1080
1090
1100
Temps (s)
1110
1120
1130
1140
1150
Figure 5.3 Signaux dentre des modles extraits des donnes exprimentales dun cycle EUDC.
113
L1 (%)
60
40
20
0
100
Mesure
Modle
80
L2 (%)
60
40
20
0
20
Mesure
Modle
L3 (%)
15
10
5
0
1050
1060
1070
1080
1090
1100
Temps (s)
1110
1120
1130
1140
1150
Figure 5.4 Validation du modle dactionneurs lectriques de la boucle des gaz du moteur
par rapport des donnes exprimentales dun cycle EUDC.
0.07
Mesure
Modle
Q1 (Kg/s)
0.06
0.05
0.04
0.03
0.02
P1 (Pa)
1.015
x 10
Mesure
Modle
1.01
1.005
304
Mesure
Modle
T3 (K)
302
300
298
296
1050
1060
1070
1080
1090
1100
Temps (s)
1110
1120
1130
1140
114
1150
x 10
Mesure
Modle
P4 (Pa)
1.8
1.6
1.4
1.2
1
1000
Mesure
Modle
T5 (K)
800
600
400
200
x 10
Mesure
Modle
P5 (Pa)
1.8
1.6
1.4
1.2
1
1050
1060
1070
1080
1090
1100
Temps (s)
1110
1120
1130
1140
1150
Nous avons constat sur les figures 5.4, 5.5a, 5.5b, une trs bonne concordance entre les
valeurs de sorties des modles simuls avec Matlab/Simulink et celles issues des donnes
exprimentales obtenues durant un cycle EUDC.
En revanche, nous avons observ un dcalage significatif entre les signaux de temprature
T3 , T5 issus du modle de simulation et ceux mesurs, notamment lors des transitions
importantes (suite des variations importantes de dbit dinjection par exemple).
Ce dcalage peut tre expliqu par le fait que les modles de temprature utiliss sont bass
directement sur les phnomnes thermodynamiques (compression, dtente, apport
nergtique), tandis que les mesures proviennent de sondes de temprature souffrant de
rponse frquentielle plutt mdiocre.
Mais comme nous allons le voir dans le chapitre VI ddi limplmentation du diagnostic,
une compensation frquentielle sous forme dun filtre passe bas sera applique aux sorties du
modle afin dliminer le dcalage en question.
A noter que cette compensation a t applique avant dtablir le calcul derreurs de
modlisation, concernant les variables de temprature, prsent dans le tableau 5.4.
Enfin, le modle complet de la boucle des gaz a t implment dans sa version temps rel et
embarqu dans le vhicule de recherche grce la plateforme XPC Target Box et les boites
outils Real Time Workshop et XPC Target de Mathworks.
115
Ermoy (%)
Ermax (% )
Q1
P1
T3
P4
T5
P5
L1
L2
L3
3.2
0.5
2
1.35
1.74
2.33
1.1
1.4
1.2
4.85
1.1
2.75
2.25
2.5
2.5
1.5
1.8
1.3
116
Chapitre VI
Dveloppement et application
du diagnostic
117
118
Dans ce chapitre, le modle de la boucle des gaz du moteur diesel dvelopp prcdemment
sera utilis dans un systme de diagnostic base de modle. Le but est de dterminer les
diffrents types de dfauts de la boucle des gaz que nous sommes capables de dtecter et
disoler en appliquant la mthode de diagnostic que nous allons dvelopper.
La validation du systme de diagnostic obtenu portera sur trois types de dfauts principaux au
sein de la boucle des gaz du moteur diesel :
- dfauts de capteurs
- dfauts dactionneurs
- fuites dans la boucle des gaz
Des rsultats exprimentaux seront illustrs pour chaque type de dfauts grce aux tests
effectus en temps rel avec un vhicule de recherche embarquant le systme de diagnostic
propos.
119
En revanche, la pression dadmission doit rester infrieure une certaine limite (dpendant du
point de fonctionnement du moteur) au-del de laquelle la consommation de carburant risque
daugmenter. En effet, llvation de la pression dadmission est ncessairement assure par
une rotation rapide du turbocompresseur. Cette rotation implique une restriction importante
des gaz dchappement du moteur par les aubes de la turbine ce qui a pour effet dlever la
pression dchappement. Ainsi, les pertes de pompage du moteur sont augmentes leur tour
ce qui entrane lexcs de consommation de carburant.
De son ct, le recyclage des gaz dchappement (EGR) joue un rle prcurseur dans la
rduction des oxydes dazote. Ce recyclage permet, selon le dbit dEGR appliqu, de
remplacer une partie de lair dadmission par des gaz dchappement. En revanche, tout excs
de dbit dEGR non surveill risque dappauvrir les gaz dadmission en oxygne ce qui
conduit la gnration de fumes lissue dune combustion incomplte. Ainsi, une
rgulation rigoureuse du dbit dEGR est un lment trs important dans le contrle de la
boucle des gaz et plus gnralement dans le contrle du moteur diesel.
A ce jour, aucun moteur diesel moderne nest quip de capteur permettant de mesurer
directement le dbit dEGR afin den tablir la rgulation.
Souvent dans ce domaine, une approche simple de contrle est applique. Celle-ci se base sur
le fait qu un point de fonctionnement donn du moteur, il existe une relation unique entre le
dbit dair, le dbit dEGR et la pression dans le collecteur dadmission. Cela signifie que
pour chaque point de fonctionnement du moteur, une seule configuration dactionneurs
(responsables du contrle de la boucle des gaz) peut raliser le triplet : {dbit dair, dbit
dEGR, pression dans le collecteur dadmission} dsir.
Traditionnellement, la tche du contrle de la pression dadmission est ralise laide du
turbocompresseur en rgulant louverture des aubes de la turbine tout en ayant le retour de
pression mesure P4. La tche du contrle du couple : {dbit dair, dbit dEGR}, quant elle,
seffectue en rgulant uniquement louverture de llectrovanne dEGR tout en ayant le retour
de dbit dair mesur Q1.
120
une capacit accrue du systme dEGR (dbits suprieurs ralisables) sans en changer le
dimensionnement.
une meilleure gestion de la boucle des gaz car il est possible datteindre les consignes de
rgulation tout en vitant au maximum la saturation des actionneurs.
une diminution de la pression dchappement (gain en consommation de carburant) est
dsormais possible car diminuer la pression dadmission avec le volet permet de raliser
les taux dEGR ncessaires avec des pressions dchappement infrieures.
Moteur
Calculateur
D1
Actionneur
Electrique
Intelligent
L1
Papillon
Q1 ref
D2
Gnrateur
de consignes
Rgulateur
de dbit
P4 ref
Actionneur
Electrique
Intelligent
L2
Vanne EGR
Rgulateur
de pression
D3
Actionneur
Electrique
Intelligent
L3
Turbocompresseur
P4
2 , Q10 , P5
Figure 6.1 Diagramme synoptique de le contrle de la boucle des gaz du moteur diesel.
121
dcrire par la suite chaque type de dfauts en illustrant ses consquences sur le
fonctionnement du systme laide de donnes exprimentales.
122
0.044
0.042
0.04
Q1 (Kg/s)
0.038
0.036
fin du dfaut
0.034
0.032
0.03
0.028
dbut du dfaut
(biais de -5 g/s)
0.026
1064
1065
1066
1067
1068
Temps (s)
1069
1070
1071
Figure 6.2 Rsultats exprimentaux dun test de dfaut sur le capteur de dbit dair.
x 10
1.5
1.45
1.4
P4 (Pa)
1.35
1.3
Dbut du dfaut
(biais de 5 KPa)
1.25
Fin du dfaut
1.2
1.15
1.1
1092
1094
1096
1098
Temps (s)
1100
1102
1104
Figure 6.3 Rsultats exprimentaux dun test de dfaut sur le capteur de pression dadmission.
123
x 10
Rponse nominale
Rponse errone
304
Rponse nominale
Rponse errone
1.8
303
1.7
302
1.6
Fin du dfaut
Dbut du dfaut
(biais de 3 K)
T3 (K)
P5 (Pa)
301
1.5
Dbut du dfaut
(biais de -5 KPa)
300
Fin du dfaut
1.4
299
1.3
298
1.2
297
1075
1080
1085
1090 1095
Temps (s)
1100
1105
1110
Les dfauts de capteurs de la boucle des gaz impactent gnralement les missions polluantes,
la consommation du carburant et lagrment de conduite du vhicule.
Un biais sur le dbit dair mesur, par exemple, modifie invitablement le rapport de la
combustion. Ainsi, en fonction du sens de la modification une augmentation des missions
des NOx ou de la gnration des fumes peut tre produite.
Quant aux impacts dune pression dadmission biaise, ils se manifestent par une
augmentation de la consommation lorsque cette pression est sous-estime par le capteur (
cause des pertes de pompage). Tandis quune pression surestime conduit des performances
dynamiques mdiocres du vhicule, accompagnes par une gnration excessive de fumes
par manque dair durant les acclrations soutenues.
124
Comme illustr sur la figure 6.1, le contrle de la boucle des gaz du moteur (par la rgulation
du dbit et de la pression dadmission) revient gnrer trois consignes de positionnement :
D1, D2, D3.
Dans le cas du bon fonctionnement, les positions des actionneurs correspondent leurs
consignes respectives (avec des temps de rponse de 80 ms environ) ce qui permet de
configurer les lments de rglage de la boucle des gaz (volet, vanne, aubes de la turbine).
Cette configuration permet au systme datteindre les objectifs de dbit et de pression
dadmission prdfinis (Q1ref et P4ref) et peut tre connue instantanment grce aux capteurs
intgrs aux actionneurs : L1, L2, L3.
Chaque actionneur intelligent dispose dun systme dautodiagnostic intgr capable de
communiquer directement au calculateur (via le rseau CAN) des informations sur ltat de
fonctionnement de lactionneur. Cependant, nous avons dcid dutiliser notre mthode de
diagnostic que nous dvelopperons ultrieurement dans ce chapitre pour la dtection des
dfauts dactionneurs. En effet, cette mthode devrait pouvoir sappliquer des actionneurs
dpourvus dintelligence (donc moins coteux) utiliss dans dautres applications.
Les dfauts dactionneurs lectriques souvent rencontrs peuvent tre classs dans deux
catgories principales :
- dfauts dorigine mcanique : actionneur drgl, bloqu
- dfauts dorigine lectrique : capteur de positionnement biais, caractristiques modifies
du moteur lectrique, court-circuit, circuit ouvert
En ce qui concerne les effets des dfauts mentionns ci-dessus, ils peuvent tre analyss selon
deux points de vue diffrents :
Point de vue de lactionneur :
En gnral, quelle que soit lorigine du dfaut (sauf en cas de blocage), la boucle
dasservissement interne de lactionneur ragit instantanment afin de le ramener au point
de fonctionnement correspondant la consigne Dx (ce qui revient occulter le dfaut en
question). En principe, aucun problme particulier nest pos par cette raction au sein de
lactionneur. Nanmoins, le fait de modifier le point de fonctionnement de celui-l (
cause dun dfaut) change invitablement la dfinition de ses positions de butes ou en
dautres termes son domaine utile de fonctionnement. En revanche, le dfaut de blocage,
lui, est plus gnant car il provoque un chauffement anormal du moteur lectrique de
lactionneur.
Point de vue de la boucle des gaz :
Comme nous lavons vu dans le point prcdent, les dfauts de lactionneur (sauf le
blocage) sont compenss par la boucle dasservissement de ce dernier de telle faon quils
deviennent invisibles par lactionneur.
Du point de vue de la boucle des gaz, cette compensation de dfauts est perue travers le
changement de point fonctionnement rel de lactionneur dfaillant (section de passage
125
Dans le but de montrer exprimentalement les deux cas prcdents, nous avons dcid de
crer des dfauts dactionneurs en temps rel sur le vhicule de recherche. Notre mthode
consistait modifier lectriquement le signal dun capteur de position intgr lun des
actionneurs lectriques durant le fonctionnement du moteur.
La figure 6.5 illustre le cas dun dfaut cr au niveau de lactionneur du turbocompresseur o
la compensation du dfaut (par la rgulation de la pression dadmission) ne conduit pas la
saturation de lactionneur.
En revanche, la figure 6.6 montre le cas contraire o la compensation de la pression
dadmission est impossible, suite au dfaut de lactionneur, cause de la saturation.
126
1.55
x 10
Pression relle
Pression nominale
1.5
P4 (Pa)
1.45
1.4
1.35
1.3
1.25
25
Position L3
L3 , D3 (%)
20
Consigne D3
15
Fin du dfaut
10
5
0
Dbut du dfaut
(biais de -10%)
-5
1065
1070
1075
Temps (s)
1080
1085
Figure 6.5 Rsultats exprimentaux dun test de dfaut dactionneur sans saturation.
1.55
x 10
Pression relle
Pression nominale
1.5
P4 (Pa)
1.45
1.4
1.35
1.3
1.25
25
L3 , D3 (%)
20
15
Position L3
Dbut du dfaut
(biais de 10%)
Consigne D3
10
Fin du dfaut
5
0
-5
1065
1070
1075
Temps (s)
1080
1085
Figure 6.6 Rsultats exprimentaux dun test de dfaut dactionneur avec saturation.
127
128
129
Nominal
Rel
Q1 (Kg/s)
0.036
0.034
0.032
0.03
0.028
x 10
Nominal
Rel
1.35
P 4 (Pa)
Dbut du dfaut
1.3
1.25
1.9
Nominal
Rel
Lambda
1.8
1.7
1.6
1.5
Nominal
Rel
L1 (%)
30
25
20
85
L2 (%)
Nominal
Rel
80
75
Dbut du dfaut
Nominal
Rel
L3 (%)
10
0
1064
1064.5
1065
1065.5
1066
Temps (s)
1066.5
1067
1067.5
Figure 6.7 Rsultats exprimentaux dun test de fuite dans le circuit dadmission.
130
1068
Q1 (Kg/s)
0.035
Nominal
Rel
0.03
0.025
P4 (Pa)
1.4
x 10
Nominal
Rel
1.3
1.2
1.1
Lambda
Nominal
Rel
2.5
2
1.5
5
x 10
P5 (Pa)
1.4
Nominal
Rel
1.3
1.2
Dbut du dfaut
45
Nominal
Rel
40
L1 (%)
35
30
25
20
15
100
Nominal
Rel
L2 (%)
90
80
70
L3 (%)
10
Nominal
Rel
Dbut du dfaut
0
1094
1096
1098
1100
1102
Temps (s)
1104
1106
1108
Figure 6.8 Rsultats exprimentaux dun test de fuite dans le circuit dchappement.
131
1110
132
Dfauts
Dfauts
Bruits
Dfauts
Capteurs
Actionneurs
Moteur Diesel
U
Gnrateur
des rsidus
Modle de bon
fonctionnement
Dtecteur
des
symptmes
Logique
de dcision
Diagnostic
des dfauts
133
Comme nous lavons dit au dbut de ce chapitre, notre mthode de diagnostic ncessite au
pralable la recherche des redondances analytiques dans le modle global du systme. Pour ce
faire, nous devons vrifier la solubilit de ce modle c'est--dire examiner si les quations
permettent de dterminer toutes les grandeurs inconnues (variables intermdiaires non
mesures) en nutilisant que des signaux connus (entres/sorties connues ou mesures).
6.3.1.1 quations du modle de bon fonctionnement
Pour faciliter le suivi des diffrentes tapes de solution des quations du modle, nous allons
rcapituler ici ces dernires (les modles complets ont t tablis au cours du chapitre 4) sous
forme de fonctions reliant les variables du modle entre elles comme suit :
Q0 = f (P0 , T1 , P1 )
(6.1)
P1 = f (T1 , Q0 , Q1 )
(6.2)
Q1 = f (P1 , P2 , 1 , T1 )
(6.3)
T2 = f (T1 , Q1 , 1 , P2 , P1 )
(6.4)
P2 = f (T2 , Q1 , Q 2 )
(6.5)
W1 = f (Q1 , T1 , T2 )
(6.6)
Q 2 = f (P2 , T2 , P3 )
(6.7)
T3 = f (T1 , T2 , Q 2 , V )
(6.8)
P3 = f (Q 2 , Q3 , T3 )
(6.9)
Q3 = f (L1 , P3 , T3 , P4 )
(6.10)
(6.11)
(
)
megr = f (Q4 , Q6 , egr )
R 4 = f ( egr )
4 = f ( egr )
(6.12a)
m air = f Q3 , Q 4 , egr
(6.12b)
(6.13)
(6.14)
P4 = f (Q3 , T3 , Q4 , T4 , Q6 , T8 )
(6.15)
T4 = f P4 , m air , m egr
(6.16)
(
)
4 = f (m air , m egr )
(6.17)
Q 4 = f ( 4 , 2 )
(6.18)
Q5 = f (Q 4 , Q10 )
(6.19)
(6.20)
ig = f (mig , m exh )
(6.21)
m exh = f Q5 , Q6 , Q 7 , ig
(6.22)
(
)
5 = f ( 4 , ig )
(6.23)
R5 = f R 4 , ig
(6.24)
P5 = f (Q5 , T5 , Q6 , Q7 , T12 , R5 , 5 )
(6.25)
T5 = f P5 , m exh , mig
(6.26)
134
(6.27)
Q6 = f (L2 , P4 , P5 , T10 )
(6.28)
T6 = f (P4 , P5 , T10 )
(6.29)
T8 = f (T6 , Q6 , T9 )
(6.30)
Q7 = f (P5 , T5 , P7 , L3 )
(6.31)
T7 = f (1 , T2 , L3 , P5 , T5 , P7 )
(6.32)
P7 = f (T7 , Q7 , Q9 )
(6.33)
W 2 = f (T5 , Q7 , T7 )
(6.34)
1 = f (W1 , W 2 )
(6.35)
Q9 = f (P0 , P7 , T7 )
(6.36)
L1 = f (D1 )
(6.37)
L2 = f (D2 )
(6.38)
L3 = f (D 3 )
(6.39)
Cette approximation a pour consquence dliminer les expressions (6.13), (6.14), (6.23)
et (6.24) de lensemble dquations rsoudre.
A partir des quations (6.3), (6.5), (6.4), (6.7), (6.9), (6.10) nous obtenons :
(6.41)
{P2, 1, T2, Q2, P3, Q3} = (P1, T1, Q1, T3, P4, L1)
A partir des quations (6.15), (6.16), (6.28), (6.27), (6.29), (6.30), (6.12a), (6.12b), (6.11)
nous obtenons :
{egr, mair, megr, Q4, T4, T8, Q6, T10, T6} = (Q3, T3, P4, T5, P5, T9, L2)
(6.42)
Compte tenu des solutions (6.41) nous obtenons :
{egr, mair, megr, Q4, T4, T8, Q6, T10, T6} = (P1, T1, Q1, L1, T3, P4, T5, P5, T9, L2)
(6.43)
A partir de lquation (6.17) nous obtenons :
(6.44)
4 = (mair, megr)
4 = (P1, T1, Q1, L1, T3, P4, T5, P5, T9, L2)
Q5 = (Q4, Q10)
Q5 = (P1, T1, Q1, L1, T3, P4, T5, P5, T9, L2, Q10)
T12 = (P1, T1, Q1, L1, T3, P4, T5, P5, T9, L2, Q10, 2)
(6.52)
Q7 = (P1, T1, Q1, L1, T3, P4, T5, P5, T9, L2, Q10)
P7 = (P1, T1, Q1, L1, T3, P4, T5, P5, T9, L2, Q10, L3)
T7 = (P1, T1, Q1, L1, T3, P4, T5, P5, T9, L2, Q10, L3)
W2 = (P1, T1, Q1, L1, T3, P4, T5, P5, T9, L2, Q10, L3)
Q9 = (P0, P1, T1, Q1, L1, T3, P4, T5, P5, T9, L2, Q10, L3)
Ainsi, nous pouvons constater que toutes les grandeurs inconnues du systme peuvent tre
dtermines laide des quations du modle de rfrence et des grandeurs mesures ou
connues de celui-ci.
137
138
r1
r2
r3
r4
r5
r6
P0
T1
x
x
x
x
x
x
Q10
T9
L1
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
L2
L3
Q1
P1
T3
P4
T5
P5
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
D1
r7
r8
r9
D2
D3
L1
L2
L3
x
x
x
x
Tableau 6.2 Rsidus lmentaires des actionneurs lectriques de la boucle des gaz
et leurs dpendances aux variables dentre/sortie du modle.
Normalement, les rsidus obtenus ( base de relations de consistance) sont gnrs de telle
faon quils soient nuls lorsque le systme est sans dfaut. En pratique, les rsidus sont
lgrement diffrents de zro pour un systme en bon fonctionnement ( cause de bruits,
dincertitudes de modle) et scartent significativement de zro lorsquun dfaut survient.
En fonction du dfaut prsent, une ou plusieurs variables du systme peuvent tre influences.
Par consquent, les rsidus utilisant ces variables risquent de ragir (diverger de zro) en
fonction de leurs sensibilits celles-ci. Ainsi, certains rsidus peuvent tre sensibles un ou
plusieurs dfauts et insensibles dautres en fonction de leurs dpendances aux variables
impactes par le dfaut prsent.
Tous les aspects prcdents mettent en valeur ce que lon appelle lisolation des dfauts qui
est une notion trs importante dans le diagnostic des systmes. Lisolation des dfauts de la
boucle des gaz sera traite plus loin dans ce chapitre travers la structuration des rsidus et la
mise en place dune logique de dcision capable didentifier les dfauts selon leurs origines.
Cependant, afin daccrotre les possibilits disolation des dfauts, nous devons essayer
dobtenir des rsidus nutilisant pas (ou peu) de variables communes.
Dans le cadre de cette tude, le systme de diagnostic propos nest valable que lorsque les
donnes dentres des modles sont supposes sans dfauts.
En effet, notre systme de diagnostic utilise un modle de processus qui ne dcrit que le
fonctionnement de la boucle des gaz du moteur.
139
Par consquent, des dfauts sur les signaux dentre tels que le rgime du moteur 2 ou le
dbit de carburant inject Q10 ne peuvent en aucun cas tre dtects par notre systme de
diagnostic. Car ce dernier ne dispose pas de modles dcrivant le fonctionnement des
systmes dont sont issues les variables dentres en question (modle dinjection, modle de
couple).
En rsum, les dpendances des rsidus obtenus prcdemment doivent tre dsormais
examines uniquement par rapport aux signaux propres notre systme (zones surlignes en
gris dans les tableaux 6.1 et 6.2).
En ce qui concerne les actionneurs lectriques, le tableau 6.2 montre que les rsidus obtenus
sont parfaitement indpendants ce qui semble avantageux pour lisolation des dfauts de ces
actionneurs.
En revanche, le tableau 6.1 ddi la boucle des gaz montre certains rsidus (comme r1 et r2)
qui utilisent des signaux communs (Q1 et P1), et dautres (comme r5 et r6) utilisant les mmes
signaux.
Cette situation est contraignante du point de vue de lisolation des dfauts. Mais dautres
rsidus additionnels peuvent tre gnrs partir des rsidus lmentaires afin dliminer
linfluence de certains signaux et de rendre les rsidus plus indpendants.
Pour ce faire, nous considrons les six relations de consistance correspondant aux rsidus
r1r6. Nous supposons inconnues les variables de sortie : Q1, P1, T3, P4, T5, P5 dans les six
quations prcdentes. Ainsi, nous avons six quations six inconnues que nous pouvons
rsoudre analytiquement.
La solution de ces quations permettra dcrire les six variables de sortie en fonction des
entres du modle comme suit :
Q1 = 1(P0, T1, V, 2, Q10, T9, L1, L2, L3)
(6.72)
P1 = 2(P0, T1, V, 2, Q10, T9, L1, L2, L3)
(6.73)
T3 = 3(P0, T1, V, 2, Q10, T9, L1, L2, L3)
(6.74)
P4 = 4(P0, T1, V, 2, Q10, T9, L1, L2, L3)
(6.75)
T5 = 5(P0, T1, V, 2, Q10, T9, L1, L2, L3)
(6.76)
P5 = 6(P0, T1, V, 2, Q10, T9, L1, L2, L3)
(6.77)
En sachant que toutes les variables utilises dans les expressions prcdentes peuvent tre
connues (entres connues, sorties mesures), les quations prcdentes peuvent enfin servir
de nouvelles relations de consistance.
Six rsidus additionnels ceux obtenus auparavant peuvent donc tre gnrs en utilisant les
nouvelles relations de consistance :
(6.78)
(6.79)
140
(6.80)
(6.81)
(6.82)
(6.83)
Le tableau 6.3 rcapitule les nouveaux rsidus additionnels obtenus et leurs dpendances aux
signaux dentres/sorties du modle.
r10
r11
r12
r13
r14
r15
P0
T1
Q10
T9
L1
L2
L3
Q1
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
P1
T3
P4
T5
P5
x
x
x
x
x
Lavantage majeur des rsidus additionnels obtenus prcdemment rside dans le fait quils
soient dcoupls les uns des autres vis--vis des variables de sortie (Q1, P1, T3) utilises et par
consquent lisolation des dfauts qui affectent les capteurs correspondants sera facilite.
La prsence commune des variables L1, L2 et L3 dans la formulation des rsidus additionnels ne
pose pas de problme particulier car les dfauts capables dinfluencer ces variables seront
couverts par les rsidus lmentaires associs (cf. tableau 6.2).
6.3.2.2 Effet du temps de rponse des capteurs
La gnration des rsidus dans un systme de diagnostic base de modle repose sur deux
lments essentiels : un modle du systme et des capteurs associs.
La qualit des informations provenant de ces deux lments ainsi que leur cohrence
dterminent directement la qualit des rsidus obtenus et les performances en termes de
dtection du systme de diagnostic utilisant ces rsidus.
Comme nous lavons dj indiqu (cf. Chapitre V, section 5.5.2), un dcalage important est
constat entre les sorties du modle et les mesures des tempratures T3 et T5 cause du temps
de rponse de leurs capteurs. Ce dcalage impacte les rsidus utilisant les signaux de
temprature mentionns prcdemment suivant la faon dont ces derniers interviennent dans
le calcul des rsidus. Dans le cas o la gnration des rsidus est base, par exemple, sur la
diffrence entre les valeurs mesures des signaux T3 ou T5 et celles issues du modle (le cas de
r12 et r14) : lamplitude de ces rsidus devient importante uniquement durant les fortes
variations de temprature dans un systme sans dfauts.
Une solution consiste utiliser dans la formule du rsidu un filtre passe bas permettant de
recaler le terme reprsentant la sortie du modle celui reprsentant la mesure.
141
200
r14 (K)
100
0
-100
-200
-300
-400
900
Mesure
Modle
800
T5 (K)
700
600
500
400
300
1050
1060
1070
1080
1090
1100
Temps (s)
1110
1120
1130
1140
1150
142
Selon lquation (6.63), le rsidu r1 est gnr partir des quations (6.2) et (6.40). Lquation
(6.40) peut tre exprime dans le domaine discret de la faon suivante (cf. Annexe C,
quation C.1) :
Q0 (n ) = 1 P0 (n ) P1 (n )
avec : 1 (n ) =
P0 (n )
(6.85)
C14 T1 (n ) C1
Par ailleurs, lquation (6.2) peut tre crite dans le domaine discret de la manire suivante
(cf. Annexe C, quation C.2) :
P1 (n ) P1 (n 1) = 2 [Q0 (n ) Q1 (n )]
avec : 2 (n ) =
C14 T1 (n )
Te
C0
(6.86)
(6.87)
(r1 = 0.92) .
Un dfaut de biais de -5 % sur le dbit dair (Q1 = 0.0645 Kg / s ) conduit une amplitude de
rsidu : r1 = 146.8 . Tandis que pour le mme point de fonctionnement, un dfaut de biais de 5 % sur la pression (P1 = 95475 Pa ) fait passer lamplitude du rsidu : r1 = 6774 .
Lexemple prcdent montre clairement que le rsidu r1 est beaucoup plus sensible un dfaut
sur P1 qu un dfaut sur Q1. Nous avons test de la mme manire dautres points de
fonctionnement qui ont donn des rsultats similaires celui de lexemple prcdent.
Ltude de la sensibilit des rsidus peut conduire en liminer, en pratique, des variables peu
influentes. Par consquent, les possibilits de dtection et disolation des dfauts peuvent tre
diffrentes en pratique de celles bases sur les tableaux de dpendances thoriques (tableaux
6.1, 6.2 et 6.3).
Dans notre cas, nous allons nous intresser particulirement la sensibilit des rsidus :
r7r15 aux variables de sortie des modles (boucle des gaz et actionneurs lectriques)
143
indiques par les zones surlignes en gris dans les tableaux 6.2 et 6.3. Dans les cas simples,
lvaluation de cette sensibilit peut-tre faite de faon analytique en calculant les drives
partielles des rsidus vis--vis des variables et ceci autour dun point de fonctionnement. La
complexit des modles nous a incit privilgier une approche exprimentale.
Notre mthode de dtermination de la sensibilit des rsidus consiste ajouter un biais
prdfini chacune des variables de sortie et calculer ensuite, par simulation, les amplitudes
de rsidus correspondantes.
Nous avons appliqu cette mthode en utilisant diffrents points de fonctionnement et niveaux
de biais qui ont conduit pratiquement aux mmes conclusions. Le tableau 6.4 montre un
exemple de rsultats obtenus pour un biais de 10% appliqu aux valeurs nominales des
variables suivantes :
L1 = 25%, L2 = 80%, L3 = 8%, Q1 = 0.04 Kg/s, P1 = 101 kPa, T3 = 299 K, P4 = 140 kPa, T5 = 780 K, P5 = 145
kPa.
r7
r8
r9
r10
r11
r12
r13
r14
r15
Q1
P1
T3
P4
T5
P5
L1
L2
L3
0
0
0
0.0038
0
0
0
0
0
0
0
0
0
10000
0
0
0
0
0
0
0
0
0
2.5
0
0
0
0
0
0
0
0
0
12500
0
0
0
0
0
0
0
0
0
54
0
0
0
0
0
0
0
0
0
-16000
2
0
0
0.0015
0
0.2
2500
-12
2000
0
8
0
0
0
0
0
0
0
0
0
2
0
0
0
-2500
-2
-3000
Tableau 6.4 Amplitudes des rsidus obtenues par simulation dun biais de 10%
sur les variables de sortie des modles utilises.
Le tableau 6.5 rsume les rsultats dtude de la sensibilit des rsidus aux variables de sortie
des modles. Dans ce tableau, les croix utilises signifient que le rsidu en question est
sensible aux variables indiques dans les colonnes. Nous avons galement renomm les
rsidus en utilisant des indices rappelant leurs origines.
Q1
rL1 (r7)
rL2 (r8)
rL3 (r9)
rQ1 (r10)
rP1 (r11)
rT3 (r12)
rP4 (r13)
rT5 (r14)
rP5 (r15)
P1
T3
P4
T5
P5
L1
L2
L3
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
Tableau 6.5 Sensibilit des rsidus aux variables de sortie des modles.
144
On peut utilement comparer les occurrences indiques dans ce tableau et celles issues des
quations (6.78) (6.83) montrant ainsi la diffrence entre sensibilit structurelle (occurrence)
et sensibilit numrique.
Rsidu
1
Grandeur
estime
(modle)
100
0.00037692
z-0.99962
Filtre passe bas
1
Seuils de dtection
Normalisation
Dtermination
du signe
Symptme
2
Grandeur
mesure
(capteur)
145
La figure 6.12 illustre les rsultats dapplication des diffrentes phases de dtection de
symptmes en temps rel (sur le vhicule de recherche) dun dfaut de capteur de dbit dair
(biais de -5 g/s).
0.06
Fin
du dfaut
Dbut
du dfaut
0.04
0.02
Rsidu
rQ1 (Kg/s)
0.01
-0.01
Rsidu filtr normalis
rQ1 (%)
10
0
-10
Seuils
de dtection
-20
Symptme
S Q1
1
0
-1
1055
1060
1065
1070
1075
1080
1085
1090
1095
1100
Temps (s)
Il existe plusieurs approches permettant dtablir lisolation des dfauts. Dans le cadre de
cette tude, nous utilisons une mthode base sur la structuration des rsidus. Cette mthode
consiste avoir un groupe de rsidus dans lequel chaque rsidu est sensible un sousensemble de dfauts. Ainsi, en dterminant quel rsidu est sensible quel dfaut, lisolation
des dfauts est effectue. Ce rsultat peut tre dcrit par une structure ou matrice dinfluence
(ou dincidence) o habituellement les dfauts surveills Fj : j = 1,2n sont placs dans
des colonnes dune matrice tandis que les diffrents rsidus ri : i = 1,2m utiliss
constituent les lignes de celle-ci.
Lorsque la valeur 1 est mise dans la case (i, j) de la structure dinfluence, cela signifie que
le rsidu ri est sensible au dfaut Fj. linverse, lorsque la case en question contient la valeur
0 , cela signifie que le rsidu ri ne ragit pas au dfaut Fj.
146
Le tableau 6.6 montre un exemple de structure dinfluence ddie aux dfauts des actionneurs
lectriques de la boucle des gaz et base essentiellement sur ltude de la sensibilit rsume
dans le tableau 6.5.
Une structure dinfluence est dite isolante , cest--dire quelle permet disoler tous les
dfauts traits, lorsque les colonnes de celle-ci sont parfaitement distinctes. Dans ce cas,
chacun des dfauts fait ragir un sous-ensemble diffrent de rsidus dfinissant une sorte de
signature propre chaque dfaut. Ainsi, une logique de dcision peut tre mise en place de
faon identifier cette signature et par consquent dterminer lorigine du dfaut dtect.
Lexemple donn par le tableau 6.6 reprsente une structure dinfluence parfaitement isolante
car chacun des dfauts constituant les colonnes de cette matrice nimpacte quun seul rsidu
la fois (matrice dinfluence diagonale).
rL1
rL2
rL3
FL1
FL2
FL3
1
0
0
0
1
0
0
0
1
Tableau 6.6 Structure dinfluence des dfauts dactionneurs lectriques de la boucle des gaz.
de la boucle des gaz du moteur diesel dfinis dans le tableau 6.8. Les valeurs contenues dans
cette table seront justifies travers des rsultats exprimentaux dapplication que nous
prsenterons dans le paragraphe suivant.
Remarque :
Pour certains dfauts, notamment ceux des capteurs ou des actionneurs, le signe du symptme
ddi dpend directement du sens de drive du composant dfaillant en question. Par
consquent, les deux symboles : + et - peuvent figurer simultanment dans la case
correspondant ce cas de figure dans la table des symptmes (obtenus selon la figure 6.12).
SL1
SL2
SL3
SQ1
SP1
ST3
SP4
ST5
SP5
F0
FQ1
FP1
FT3
FP4
FT5
FP5
FL1
FL2
FL3
Fint
Fexh
0
0
0
0
0
0
0
0
0
0
0
0
+0
0
0
0
0
0
0
0
0
+0
0
0
0
0
0
0
0
0
+0
0
0
0
0
0
0
0
0
+0
0
0
0
0
0
0
0
0
+0
0
0
0
0
0
0
0
0
+-
+0
0
+0
++-+
+-
0
+0
0
0
0
0
0
0
0
0
+-0
0
-+
+-+
0
0
0
+
0
0
+
-
0
0
0
0
+
-
Tableau 6.7 Table exprimentale des symptmes de dfauts de la boucle des gaz.
Notation
F0
FQ1
FP1
FT3
FP4
FT5
FP5
FL1
FL2
FL3
Fint
Fexh
Dfinition
Fonctionnement nominal (0 dfaut)
Dfaut de capteur de dbit dair dadmission
Dfaut de capteur de pression aprs filtre air
Dfaut de capteur de temprature dair comprim refroidi
Dfaut de capteur de pression dans le collecteur dadmission
Dfaut de capteur de temprature des gaz dans le collecteur dchappement
Dfaut de capteur de pression dans le collecteur dchappement
Dfaut dactionneur lectrique du volet dadmission
Dfaut dactionneur lectrique de la vanne de recyclage des gaz EGR
Dfaut dactionneur lectrique des aubes du turbocompresseur
Fuite de lair comprim dans le circuit dadmission (avant collecteur)
Fuite des gaz dans le circuit dchappement (avant turbine)
Tableau 6.8 Liste rcapitulative des dfauts traits dans le cadre du diagnostic
de la boucle des gaz du moteur diesel OM646.
148
Afin de tester la raction de ces indicateurs de pannes, diffrents types de dfauts ont t
provoqus en temps rel soit par modification du programme de contrle dans le calculateur
lectronique (changement des paramtres de conversion des signaux des capteurs), soit par
modification physique de la boucle des gaz du moteur (cration dorifice afin de provoquer
des fuites dans les circuits dadmission et dchappement).
Les tests de dfauts mentionns prcdemment, dont nous allons prsenter les rsultats par la
suite, ont t effectus durant un cycle de circulation extra-urbaine de rfrence (EUDC).
Pour chaque cas de figure (avec ou sans dfauts) nous allons donc tracer les rsidus filtrs et
normaliss avec leurs seuils de dtection (en pointill) ainsi que les symptmes
correspondants obtenus. Ainsi, nous justifierons le contenu de la table des symptmes que
nous avons dresse par anticipation dans le paragraphe prcdent (tableau 6.7).
Rsidu
rL1
Seuil de dtection
1
rL2
rL3
rQ1
rP1
0.5
rT3
rP4
2.5
rT5
2.5
rP5
2.5
149
symptmes 6.7 nous dduisons les points suivants concernant linfluence des dfauts de
capteurs sur les symptmes :
- Parmi les six symptmes ddis la boucle des gaz (SQ1, SP1, ST3, SP4, ST5, SP5), seul un
symptme est impact la fois par le dfaut de capteur correspondant. Par exemple :
lorsquun dfaut de capteur de pression FP4 est prsent, seul le symptme SP4 est influenc.
- Le signe du symptme impact par le dfaut de capteur prsent dpend directement du
signe de ce dfaut. Par exemple : lorsquil y a surestimation de la grandeur mesure par le
capteur, le signe du symptme est positif et vice-versa.
- Les symptmes associs aux actionneurs : SL1, SL2, SL3 ne sont pas influencs par les dfauts
de capteurs de la boucle des gaz
En effet, les dfauts de capteurs de la boucle des gaz du moteur sont dtectables condition
que leurs amplitudes soient suprieures aux seuils de dtection correspondants.
Il est galement possible de dterminer le capteur dfaillant grce la signature spcifique de
son dfaut. Cette signature consiste avoir un seul symptme (positif ou ngatif) de dfaut de
capteur en association avec des valeurs nulles pour les autres symptmes comme lindique la
table des symptmes 6.7.
Dfaut
FQ1
FP1
FT3
FP4
FT5
FP5
Dbut
1065 s
1075 s
1085 s
1095 s
1105 s
1115 s
Fin
1070 s
1080 s
1090 s
1100 s
1110 s
1120 s
Amplitude
-8 %
2%
10 %
4%
-5 %
-4 %
Tableau 6.10 Caractristiques des dfauts de capteurs utiliss durant les tests exprimentaux.
150
Dbut
1070 s
1100 s
1130 s
Fin
1080 s
1110 s
1140 s
Amplitude
-7 %
7%
-7 %
Tableau 6.11 Caractristiques des dfauts dactionneurs utiliss durant le test exprimental.
6.3.4.2.4 Fuites
Afin de mettre en vidence linfluence des fuites sur les symptmes de dfauts, nous illustrons
sur les figures 6.16 et 6.17 les rsultats de deux tests de fuite effectus respectivement dans
les circuits dadmission et dchappement. Les caractristiques de ces dfauts sont dfinies
dans le tableau 6.12.
Compte tenu des rsultats de tests de fuite et de la table des symptmes 6.7 nous pouvons
mettre les remarques suivantes :
- Lorsquune fuite est prsente dans la boucle des gaz du moteur, plusieurs symptmes de
dfauts de capteurs peuvent tre affects simultanment.
- Les symptmes de dfauts dactionneur ne sont pas influencs par la prsence des fuites.
- Certaines combinaisons de symptmes caractrisent clairement le cas dune fuite. Par
exemple, les symptmes ngatifs de dfaut de capteur de pressions dadmission et
dchappement SP4, SP5 ainsi que le symptme positif de dfaut de capteur de temprature
dchappement ST5 sont lis la sollicitation croissante du turbocompresseur pour
compenser les fuites.
- Dautres symptmes affects simultanment avec ceux voqus ci-dessus offrent une aide
importante la discrimination entre une fuite dans le circuit dadmission et celle dans le
circuit dchappement. Nous pouvons citer par exemple le symptme de dfaut de capteur
qui est positif en cas de fuite dadmission et ngatif en cas de fuite dchappement.
Test 1
Test 2
Dfaut
Fint
Fexh
Dbut
1070 s
1070 s
Tableau 6.12 Caractristiques des fuites introduites durant les tests exprimentaux.
151
0.5
rL1 (%)
rP1 (%)
rQ1 (%)
0.5
0
0
-0.5
-5
-0.5
-1
rL2 (%)
0.5
rP4 (%)
rT3 (%)
1
0
-1
0
-0.5
-2
-5
-1
rL3 (%)
rP5 (%)
rT5 (%)
0.5
0
0
-0.5
1080
1100
1120
Temps (s)
-5
1060
1140
-1
1060
1140
0.5
0.5
0.5
0
-0.5
-1
-1.5
1.5
1.5
1.5
0.5
0.5
0.5
S L2
-1
-1.5
0
-0.5
-1
-1.5
1.5
1.5
1.5
0.5
0.5
0.5
SL3
-1
-1.5
SP5
-1
0
-0.5
-1
-1
-1
-1.5
1060
-1.5
1060
1100
1120
Temps (s)
1140
1100
1120
Temps (s)
1140
0
-0.5
-1.5
1060
1080
1080
0
-0.5
-1.5
-0.5
1140
-1
1100
1120
Temps (s)
-0.5
-1.5
-0.5
1080
1.5
SL1
1.5
S P4
S T3
1100
1120
Temps (s)
-0.5
ST5
1080
1.5
SP1
S Q1
-5
1060
1080
1100
1120
Temps (s)
1140
152
10
-10
0.5
rL1 (%)
rP1 (%)
rQ1 (%)
-0.5
-1
-20
-2
-1
20
rP4 (%)
rL2 (%)
0.5
10
rT3 (%)
0
-0.5
-10
-5
-1
rL3 (%)
rP5 (%)
rT5 (%)
0.5
0
0
-0.5
1080
1100
1120
Temps (s)
-5
1060
1140
-1
1060
1140
0.5
0.5
0.5
0
-0.5
-1
-1.5
1.5
1.5
1.5
0.5
0.5
0.5
S L2
-1
-1.5
0
-0.5
-1
-1.5
1.5
1.5
1.5
0.5
0.5
0.5
SL3
-1
-1.5
SP5
-1
0
-0.5
-1
-1
-1
-1.5
1060
-1.5
1060
1100
1120
Temps (s)
1140
1100
1120
Temps (s)
1140
0
-0.5
-1.5
1060
1080
1080
0
-0.5
-1.5
-0.5
1140
-1
1100
1120
Temps (s)
-0.5
-1.5
-0.5
1080
1.5
SL1
1.5
S P4
S T3
1100
1120
Temps (s)
-0.5
ST5
1080
1.5
SP1
S Q1
-5
1060
1080
1100
1120
Temps (s)
1140
Figure 6.14 Rsultats exprimentaux des rsidus et des symptmes correspondants obtenus
dans le cas de dfauts de capteurs de la boucle des gaz du moteur diesel.
153
-20
0
-5
10
-2
rL2 (%)
rP4 (%)
-10
-10
-5
-20
-10
20
10
10
10
1080
1100
1120
Temps (s)
-10
-20
1060
1140
rL3 (%)
-4
1080
1100
1120
Temps (s)
0
-5
1060
1140
0.5
0.5
0.5
0
-0.5
-1
-1.5
1.5
1.5
1.5
0.5
0.5
0.5
S L2
-1
-1.5
0
-0.5
-1
-1.5
1.5
1.5
1.5
0.5
0.5
0.5
SL3
-1
-1.5
SP5
-1
0
-0.5
-1
-1
-1
-1.5
1060
-1.5
1060
1100
1120
Temps (s)
1140
1100
1120
Temps (s)
1140
0
-0.5
-1.5
1060
1080
1080
0
-0.5
-1.5
-0.5
1140
-1
1100
1120
Temps (s)
-0.5
-1.5
-0.5
1080
1.5
SL1
1.5
SP1
1.5
S P4
S T3
-0.5
-0.5
ST5
-40
rP5 (%)
rT3 (%)
rT5 (%)
rP1 (%)
rQ1 (%)
-10
1060
S Q1
0.5
rL1 (%)
20
1080
1100
1120
Temps (s)
1140
Figure 6.15 Rsultats exprimentaux des rsidus et des symptmes correspondants obtenus
dans le cas de dfauts dactionneurs de la boucle des gaz du moteur diesel.
154
20
0.5
15
rP1 (%)
rQ1 (%)
rL1 (%)
0.5
10
0
-0.5
0
-5
-0.5
-1
0
-1
0.5
rL2 (%)
rP4 (%)
rT3 (%)
-5
0
-0.5
-2
-3
-10
10
-1
0.5
rL3 (%)
rP5 (%)
rT5 (%)
-2
-0.5
-4
1065
1070 1075
Temps (s)
1080
-6
1060
1085
1080
-1
1060
1085
0.5
0.5
0.5
0
-0.5
-1
-1.5
1.5
1.5
1.5
0.5
0.5
0.5
S L2
-1
-1.5
0
-0.5
1.5
1.5
1.5
0.5
0.5
0.5
SL3
-1
-1.5
SP5
-1
-1.5
0
-0.5
-1
-1
-1
-1.5
1060
-1.5
1060
1070 1075
Temps (s)
1080
1085
1070 1075
Temps (s)
1080
1085
0
-0.5
-1.5
1060
1065
1065
-1
1085
-0.5
-1.5
-0.5
1080
-1
1070 1075
Temps (s)
-0.5
-1.5
-0.5
1065
1.5
SL1
1.5
S P4
S T3
1070 1075
Temps (s)
-0.5
ST5
1065
1.5
SP1
S Q1
-5
1060
1065
1070 1075
Temps (s)
1080
1085
Figure 6.16 Rsultats exprimentaux des rsidus et des symptmes correspondants obtenus
dans le cas dune fuite dans le circuit dadmission du moteur diesel.
155
10
0.5
rL1 (%)
0.5
rP1 (%)
rQ1 (%)
5
0
-0.5
-5
rP4 (%)
rT3 (%)
-5
-0.5
-1
10
0.5
rL2 (%)
-10
0
-5
-10
0
-5
10
0.5
rL3 (%)
rP5 (%)
5
rT5 (%)
-1
-15
10
0
-5
1070 1075
Temps (s)
1080
-15
1060
1085
1080
-1
1060
1085
0.5
0.5
0.5
0
-0.5
-1
-1.5
1.5
1.5
1.5
0.5
0.5
0.5
S L2
-1
-1.5
0
-0.5
1.5
1.5
1.5
0.5
0.5
0.5
SL3
-1
-1.5
SP5
-1
-1.5
0
-0.5
-1
-1
-1
-1.5
1060
-1.5
1060
1070 1075
Temps (s)
1080
1085
1070 1075
Temps (s)
1080
1085
0
-0.5
-1.5
1060
1065
1065
-1
1085
-0.5
-1.5
-0.5
1080
-1
1070 1075
Temps (s)
-0.5
-1.5
-0.5
1065
1.5
SL1
1.5
S P4
S T3
1070 1075
Temps (s)
-0.5
ST5
1065
1.5
SP1
S Q1
1065
0
-0.5
-10
-10
0
-0.5
-10
1060
1065
1070 1075
Temps (s)
1080
1085
Figure 6.16 Rsultats exprimentaux des rsidus et des symptmes correspondants obtenus
dans le cas dune fuite dans le circuit dchappement du moteur diesel.
156
6.3.4.3 Synthse
Le but de cette partie est de dfinir, partir des rsultats exprimentaux et de la table des
symptmes prsents prcdemment, les rgles de dtection et disolation appliquer en
temps rel afin dtablir le diagnostic du systme.
Ces rgles seront exprimes sous forme de combinaisons logiques de conditions utilisant les
neuf symptmes de dfauts prsents dans le tableau 6.7.
Ltablissement des rgles logiques de dtection et disolation des dfauts est bas
essentiellement sur lanalyse de leurs signatures reprsentes par les diffrentes colonnes de la
table exprimentale des symptmes. Le tableau 6.13 dfinit les diffrentes rgles constituant
la logique de dcision utilise dans le diagnostic de la boucle des gaz du moteur diesel.
Remarque :
Labsence du symptme de dfaut de capteur T3 des conditions de dtection des fuites est lie
la faible sensibilit de ce symptme ce type de dfauts compare celle des autres
symptmes utiliss (cf. figure 6.16). Ainsi, lintgration du symptme en question dans les
conditions de dtection des fuites dadmission par exemple empcherait la dtection de petites
fuites de faon prcoce.
Diagnostic
F0
FQ1
FP1
FT3
FP4
FT5
FP5
FL1
FL2
FL3
Fint
Fexh
Rgle utilise
(SL1 = 0)(SL2 = 0)(SL3 = 0)(SQ1 = 0)(SP1 = 0)(ST3 = 0)(SP4 = 0)(ST5 = 0)(SP5 = 0)
(SL1 = 0)(SL2 = 0)(SL3 = 0)(SQ1 0)(SP1 = 0)(ST3 = 0)(SP4 = 0)(ST5 = 0)(SP5 = 0)
(SL1 = 0)(SL2 = 0)(SL3 = 0)(SQ1 = 0)(SP1 0)(ST3 = 0)(SP4 = 0)(ST5 = 0)(SP5 = 0)
(SL1 = 0)(SL2 = 0)(SL3 = 0)(SQ1 = 0)(SP1 = 0)(ST3 0)(SP4 = 0)(ST5 = 0)(SP5 = 0)
(SL1 = 0)(SL2 = 0)(SL3 = 0)(SQ1 = 0)(SP1 = 0)(ST3 = 0)(SP4 0)(ST5 = 0)(SP5 = 0)
(SL1 = 0)(SL2 = 0)(SL3 = 0)(SQ1 = 0)(SP1 = 0)(ST3 = 0)(SP4 = 0)(ST5 0)(SP5 = 0)
(SL1 = 0)(SL2 = 0)(SL3 = 0)(SQ1 = 0)(SP1 = 0)(ST3 = 0)(SP4 = 0)(ST5 = 0)(SP5 0)
SL1 0
SL2 0
SL3 0
(SL1 = 0)(SL2 = 0)(SL3 = 0)(SQ1 > 0)(SP4 < 0)(ST5 > 0)(SP5 < 0)
(SL1 = 0)(SL2 = 0)(SL3 = 0)(SQ1 < 0)(SP4 < 0)(ST5 > 0)(SP5 < 0)
Tableau 6.13 Rgles logiques de diagnostic de la boucle des gaz du moteur diesel OM646.
157
Le point de fonctionnement correspondant du vhicule est dcrit par les valeurs suivantes :
2 = 2500 tr/min, Q10 = 25 mg/coup, V = 75 Km/h.
Il est important de noter que les dfauts dtects sont de faibles amplitudes (10% pour le plus
important dentre eux). Nanmoins, ces dfauts peuvent avoir des consquences sensibles sur
le niveau dmissions polluantes du moteur (augmentation des Nox ou des particules) do
limportance de leur dtection. Bien videmment, pour des dfauts damplitude plus
importante, par exemple de lordre de 25%, la dtection est quasi parfaite mme en prsence
158
de bruit. Au vu de ces rsultats, nous concluons de faon trs positive quant la capacit de
lapproche propose dtecter et isoler correctement les dfauts tudis.
Q1 = 0.04 Kg/s
FQ1 = 5 %
FQ1 = 8 %
FQ1 = 10 %
Bruit 1 %
62 / 0
89 / 0
100 / 0
Bruit 3 %
37 / 8
76 / 12
82 / 10
Bruit 5 %
23 / 30
65 / 28
70 / 23
Tableau 6.14 Taux de dtection / fausse dtection dun dfaut de capteur FQ1.
L2 = 80%
FL2 = 5 %
FL2 = 8 %
FL2 = 10 %
Bruit 1 %
85 / 0
93 / 0
100 / 0
Bruit 3 %
72 / 13
81 / 9
88 / 14
Bruit 5 %
59 / 48
70 / 24
66 / 29
Tableau 6.15 Taux de dtection / fausse dtection dun dfaut dactionneur FL2.
Fint : 2mm
Fint : 3mm
Fint : 4mm
Bruit 1 %
Bruit 3 %
Bruit 5 %
53 / 0
73 / 0
93 / 0
39 / 0
59 / 7
74 / 10
24 / 8
40 / 20
50 / 23
Tableau 6.16 Taux de dtection / fausse dtection dune fuite dans le circuit dadmission Fint.
159
160
Conclusion et Perspectives
Conclusion et
Perspectives
161
Conclusion et Perspectives
Le travail men dans le cadre de cette thse, au sein de Delphi Diesel Systems quipementier
automobile, a t dintroduire et de promouvoir lapproche de diagnostic base de modles
travers une application un moteur diesel suraliment injection directe.
Dans un contexte de dveloppement incessamment croissant de stratgies de contrle et de
gestion des moteurs diesel grce llectronique embarque, le sujet de cette thse a t initi
pour rpondre deux besoins principaux :
- construire un modle de simulation complet de la boucle des gaz du moteur diesel dans le
but de comprendre et danalyser les effets de diffrents types de dfauts dans cette boucle,
- proposer et exprimenter en temps rel une nouvelle approche de diagnostic plus robuste
et plus subtile en termes de dtection et de discrimination des dfauts que les approches
actuelles plutt basiques voire insuffisantes.
Afin datteindre le premier objectif de modlisation, un recourt vers un modle du moteur
valeur moyenne a t motiv par la simplicit de cette famille de modles et par sa forte
prsence dans la littrature.
Une contribution importante de ce travail de thse a t de construire un modle de simulation
valeur moyenne dun moteur diesel destin particulirement la boucle des gaz de celui-ci.
Le recalage paramtrique et la validation du modle obtenu ont t galement effectus avec
succs laide dun vhicule de recherche quip dun moteur diesel moderne. La boucle des
gaz de ce moteur est dote dun volet dadmission, dun turbocompresseur gomtrie
variable et dun systme de recyclage des gaz dchappement. Des capteurs dorigine
fournissent un certain nombre de mesures de dbit, de pression et de temprature en plus des
positions instantans des actionneurs de la boucle.
Grce au modle obtenu dans la premire phase, des simulations ont pu tre effectues dans le
but de constater et danalyser linfluence des diffrents types de dfauts, de capteurs et
dactionneurs ainsi que des fuites, sur le systme.
Une deuxime contribution de cette thse a t donc de proposer une approche de diagnostic
base de modles utilisant lapproche des rsidus structurs. Dabord, cette mthode a t
teste par simulation avec le modle de la boucle des gaz du moteur.
Au vu des rsultats de simulation encourageants, la mthode de diagnostic propose a t
ensuite implmente et teste en temps rel sur le vhicule de recherche grce aux techniques
de prototypage rapide. Les rsultas exprimentaux de lapplication ont montr la capacit de
la mthode dtecter et isoler les trois types de dfauts indiqus prcdemment avec des
performances et robustesse vis--vis du bruit satisfaisantes.
Comme tout travail de recherche industrielle, cette thse na fait que marquer quelques pas en
avant notamment sur le plan de lapplication et du transfert de la thorie vers le domaine
technologique et technique. Plusieurs ides de dveloppement futur pourraient faire suite ce
travail :
162
Conclusion et Perspectives
Sur le plan de la modlisation, une piste damlioration du modle pourrait par exemple
enrichir la partie chambre de combustion du modle actuel. Ceci permettrait de
prendre en compte les effets des nouvelles techniques dinjection (multi-injections) sur la
temprature des gaz dchappement.
Sur le plan du diagnostic, une rflexion sur la prise en compte du phnomne de
vieillissement des composants du systme dans la conception de la mthode du diagnostic
pourraient tre fortement utile.
Sur le plan de lapplication, des optimisations et des simplifications diverses (de modle,
de calcul de rsidus) pourraient tre apportes dans le but de sadapter aux contraintes
de performances et de cot des calculateurs lectroniques utiliss dans lautomobile.
163
Conclusion et Perspectives
164
Annexes
Annexes
165
Annexes
Annexe A : Abrviations
CAN
ECU
EGR
EUDC
HCCI
MVEM
NEDC
NOx
OBD
PMB
PMH
RAM
VGT
166
Annexes
Annexe B : Notations
Notation
Nature
Unit
Description
T1
Variable
T2
Variable
T3
Variable
T4
Variable
T5
Variable
T6
Variable
T7
Variable
T8
Variable
T9
Variable
T10
Variable
T11
Variable
T12
Variable
P0
Variable
Pa
Pression atmosphrique
P1
Variable
Pa
P2
Variable
Pa
P3
Variable
Pa
P4
Variable
Pa
P5
Variable
Pa
P7
Variable
Pa
Q0
Variable
Kg/s
Q1
Variable
Kg/s
Q2
Variable
Kg/s
Q3
Variable
Kg/s
Q4
Variable
Kg/s
Q5
Variable
Kg/s
Q6
Variable
Kg/s
Q7
Variable
Kg/s
Q9
Variable
Kg/s
Q10
Variable
Kg/s
Qinj
Variable
mg/coup
Variable
Km/h
Variable
rad/s
167
Annexes
Variable
rpm
Variable
coup/s
Variable
Variable
Variable
Variable
Variable
Pourcentage de lnergie apporte par le carburant et permettant de rchauffer les gaz dans la
chambre de combustion
Variable
Variable
Variable
W1
Variable
watt
W2
Variable
watt
D1
Variable
D2
Variable
D3
Variable
L1
Variable
L2
Variable
L3
Variable
mair
Variable
Kg
megr
Variable
Kg
egr
Variable
Variable
R4
Variable
J/(Kg.K)
cv 4
Variable
J/(Kg.K)
c p4
Variable
J/(Kg.K)
Variable
Kg/m3
mexh
Variable
Kg
mig
Variable
Kg
ig
Variable
Variable
Rapport des capacits thermiques spcifiques des gaz dans le collecteur dchappement
R5
Variable
J/(Kg.K)
cv5
Variable
J/(Kg.K)
Capacit thermique spcifique volume constant des gaz dans le collecteur dchappement
c p5
Variable
J/(Kg.K)
Capacit thermique spcifique pression constante des gaz dans le collecteur dchappement
C0
constante
m3
C1
constante
Pa2.s2/Kg2.K
Rgime moteur
Cadence des phases du cycle moteur
168
Annexes
C2
constante
Pa2.s2/Kg2.K
C3
constante
J/Kg
C4
constante
C5
constante
J/(Kg.K)
C6
constante
J/(Kg.K)
C7
constante
C8
constante
C9
constante
J/(Kg.K)
C10
constante
J/(Kg.K)
C11
constante
C12
constante
Kg.m2
C13
constante
Kg.J/ Pa2.s2
C14
constante
J/(Kg.K)
C15
constante
J/(Kg.K)
C16
constante
m3
C17
constante
m3
C18
constante
m3
C19
constante
m3
Cylindre du moteur
C20
constante
C21
constante
m3
C22
constante
m3
C23
constante
C24
constante
C25
constante
C26
constante
Pref
constante
Pa
Tref
constante
= c p / cv
cp
cv
C 4 1
C4
169
Annexes
Filtre air
Q0 (t ) =
P0 (t )
C14 T1 (t )
[P0 (t ) P1 (t )]
C14 T1 (t )
C0
P1 (t ) P1 (t 0 ) =
(C.1)
C1
t
[Q (t ) Q (t )] dt
0
(C.2)
t0
Compresseur
P (t ) Tref
P (t )
Q1 (t ) = f Q1cor 2 , 1 (t ) 1
Pref T1 (t )
P1 (t )
T2 (t ) = T1 (t ) +
(C.3)
P (t ) C7
T1 (t )
2 1
f 1 (Q1 (t ), 1 (t )) P1 (t )
P2 (t ) P2 (t 0 ) =
C14T2 (t )
C16
(C.4)
[Q (t ) Q (t )] dt
1
(C.5)
t0
W1 (t ) = C 5 Q1 (t )[T2 (t ) T1 (t )]
(C.6)
P2 (t )
C14T2 (t )
[P2 (t ) P3 (t )]
(C.7)
C2
T3 (t ) = T2 (t ) f 2 (Q2 (t ), V (t ))[T2 (t ) T1 (t )]
P3 (t ) P3 (t 0 ) =
C14T3 (t )
C17
(C.8)
[Q (t ) Q (t )] dt
2
(C.9)
t0
170
Annexes
Volet dadmission
2
C4 +1
P3 (t )
2 P4 (t ) C4 P4 (t ) C4
Athr (t )
P (t )
C14T3 (t ) C 7 P3 (t )
2
C
1
+
Q3 (t ) =
4
P3 (t )
2 2 C4 1 2 C4 1
Athr (t )
C 7 C 4 + 1
C 4 + 1
(
)
C
T
t
14
3
C4
Si
P4 2 C4 1
>
P3 C 4 + 1
(C.10)
Sinon
Collecteur dadmission
egr (t ) =
m egr (t )
(C.11)
m air (t ) + m egr (t )
t
t
.
m air (t ) m air (t 0 ) = m air dt = Q3 (t ) Q4 (t ) 1 egr (t ) dt
t0
t0
t
t
.
(
)
(
)
=
=
m
t
m
t
m
dt
Q6 (t ) Q4 (t ) egr (t ) dt
egr
egr
egr 0
t0
t0
)]
(C.12)
Ou :
t
t
.
m air (t ) m air (t 0 ) = m air dt = Q3 (t ) (Q 4 (t ) + Q6 (t )) 1 egr (t ) dt
t0
t0
t
t
.
t0
t0
R 4 (t ) =
)]
m air (t ) + m egr (t )
1
C18
(C.14)
[Q (t )T (t )C
6
15C8
(C.13)
(C.15)
t0
Ou :
P4 (t ) P4 (t 0 ) =
T4 (t ) =
1
C18
[Q (t )T (t )C
3
14 C 4
(Q 4 (t ) + Q6 (t ))T4 (t )R 4 (t ) 4 ] dt
t0
P4 (t )C18
R4 (t )[mair (t ) + megr (t )]
(C.16)
171
Annexes
4 (t ) =
m air (t ) + m egr (t )
(C.17)
C18
Chambre de combustion
2 (t )
Q4 (t ) = f 4 ( 4 (t ), 2 (t ))
30
4 (t )
C19
4
(C.18)
Q5 (t ) = Q4 (t ) + Q10 (t )
T12 (t ) = T4 (t ) +
(C.19)
C 3 Q10 (t ) f5 (Q10 (t ), 2 (t ))
(C.20)
C 9 Q5 (t )
Collecteur dchappement
ig (t ) =
m ig (t )
(C.21)
m exh (t ) + m ig (t )
t
t
.
m exh (t ) m exh (t 0 ) = m exh dt = Q5 (t ) (Q6 (t ) + Q7 (t )) 1 ig (t ) dt
t0
t0
t
t
.
m ig (t ) m ig (t 0 ) = m ig dt = (Q6 (t ) + Q7 (t )) ig (t ) dt
t0
t0
)]
(C.22)
Ou :
t
t
.
m exh (t ) m exh (t 0 ) = m exh dt = Q5 (t ) Q7 (t ) 1 ig (t ) dt
t0
t0
t
t
.
m ig (t ) m ig (t 0 ) = m ig dt = Q6 (t ) Q 7 (t ) ig (t ) dt
t0
t0
R5 (t ) =
R 4 m ig (t ) + C15 m exh (t )
m exh (t ) + m ig (t )
)]
= R 4 ig (t ) + C15 1 ig (t )
5 (t ) = ig (t ) 4 + (1 ig (t ))C8
(C.24)
[Q (t )T
1
P5 (t ) P5 (t 0 ) =
C 21
12
t0
Ou :
P5 (t ) P5 (t 0 ) =
1
C 21
[Q (t )T
5
(C.23)
12
t0
172
(C.25)
Annexes
T5 (t ) =
P5 (t )C 21
R5 (t )[m exh (t ) + mig (t )]
(C.26)
Ou :
(C.27)
5 +1
2
P5 (t )
2 P4 (t ) 5 P4 (t ) 5
Aegr (t )
P (t )
R5T10 (t ) 5 1 P5 (t )
5 +1
2
Q6 (t ) =
2 2 5 1 2 5 1
A (t ) P5 (t )
egr
5 1 5 + 1
5 + 1
(
)
R
T
t
5
10
A = f
Aegr (L 2 ) = Aegr eff (L 2 ) C egr (L 2 )
egr
Si
P4 2 5 1
>
P5 5 + 1
(C.28)
Sinon
Ou :
4 +1
2
P4 (t )
2 P5 (t ) 4 P5 (t ) 4
Aegr (t )
R 4 T6 (t ) 4 1 P4 (t )
P4 (t )
4 +1
2
Q6 (t ) =
2 2 4 1 2 4 1
A (t ) P4 (t )
egr
4 1 4 + 1
4 + 1
(
)
R
T
t
4
6
A = f (L ) = A
Aegr
2
egr eff (L 2 ) C egr (L 2 )
egr
P (t )
T6 (t ) = T10 4
P5 (t )
Si
P5 2 4 1
>
P4 4 + 1
Sinon
5 1
5
(C.29)
Ou :
P (t )
T10 (t ) = T6 5
P4 (t )
4 1
4
T8 (t ) = T6 (t ) f 7 (Q6 (t ))[T6 (t ) T9 (t )]
Ou
(C.30)
T6 (t ) = T4 (t ) f 7 (Q6 (t ))[T4 (t ) T9 (t )]
Annexes
Turbine
P (t )
P (t ) Tref
Q7 (t ) = f Q7 cor 5 , L3 (t ) 5
P7 (t )
Pref T5 (t )
(C.31)
P (t ) C11
P5 (t )
T7 (t ) = T5 (t ) f 3
, 1 (t ), T5 (t ), T2 (t ), L3 (t )T5 (t )1 5
P7 (t )
P7 (t )
(C.32)
P7 (t ) P7 (t0 ) =
C15T7 (t )
C22
[Q (t ) Q (t )] dt
7
(C.33)
t0
W 2 (t ) = C 9 Q7 (t )[T5 (t ) T7 (t )]
1 (t ) 1 (t0 ) =
C121 (t0 )
(C.34)
[W (t ) W (t )] dt
2
(C.35)
t0
Systme dchappement
C8 +1
2
P7 (t )
2 P0 (t ) C8 P0 (t ) C8
C13
P (t )
C15 T7 (t ) C11 P7 (t )
C
+
2
Q9 (t ) =
8 1
P7 (t )
2 2 C8 1 2 C8 1
C13
C15T7 (t ) C11 C 8 + 1
C8 + 1
C8
Si
P0 2 C8 1
>
P7 C 8 + 1
(C.36)
Sinon
Actionneurs lectriques
t t
L1 (t ) L1 (t 0 ) = D1 (t )1 e C24
t t
L 2 (t ) L 2 (t 0 ) = D 2 (t )1 e C25
t t
L3 (t ) L3 (t 0 ) = D3 (t )1 e C26
(C.37)
174
Rfrences
Rfrences
175
Rfrences
[1]
M. Nyberg, T. Stutte, and V. Wilhelmi. Model based diagnosis of the air path of an
automotive diesel engine. IFAC Workshop : Advances in Automotive Control,
Karlsruhe, Germany, March 2001.
[2]
[3]
[4]
[5]
[6]
[7]
[8]
[9]
F. Swartling. Gas flow observer for Diesel Engines with EGR. Masters thesis,
Linkping University, Sweeden, 2005.
176
Rfrences
177
Rfrences
[29] M. Nyberg, A. Perkovic. Model Based Diagnosis of Leaks in the Air-Intake System of
an SI-Engine. Volume 980514 of SAE Paper, 1998.
[30] M. Nyberg. SI-Engine Air-Intake System Diagnosis by Automatic FDI-Design. IFAC
Workshop Advances in Automotive Control. Columbus, Ohio, pp. 225230, 1998.
[31] M. Nyberg, T. Stutte, V. Wilhelmi. Model Based Diagnosis of the Air Path of an
Automotive Diesel Engine. IFAC Workshop: Advances in Automotive Control.
Karlsruhe, Germany, pp. 629634, 2001.
[32] Z. Sabeh, J. Ragot, F. Kratz. Structured hypothesis tests based diagnosis : application
to common rail diesel injection system. Advances in Vehicle Control and Safety,
Genova, Italy, 2004.
[33] R. Isermann, A. Schwarte, F. Kimmich. Model-based fault detection of a Diesel engine
with turbo charger - a case study. 1st IFAC Symposium on Advances in Automotive
Control, University of Salerno, Italien, 2004.
[34] J C. Peyton Jones, K R. Muske, J W. Howse Model-Based Fault Detection for ThreeWay Automotive Catalyst Systems. 1st IFAC Symposium on Advances in Automotive
Control, University of Salerno, Italien, 2004.
[35] D. Fischer, H.-P. Schner, R. Isermann, Model-based fault detection for an active
vehicle suspension. IFAC Symposium on Advances in Automotice Control, Salerno,
Italiten, 2004.
[36] R. Isermann, D. Fischer, M. Brner, J. Schmitt, Fault detection for lateral and vertical
vehicle dynamics. 3rd IFAC Symposium on Mechatronic Systems, Sidney, Australien,
2004.
[37] http://www.mathworks.fr/products/xpctargetbox/
[38] http://www.mathworks.fr/products/rtw/
[39] http://www.mathworks.fr/products/xpctarget/
[40] http://www.micro-epsilon.com/en/Sensors/turboSPEED-DZ135/
[41] G. Zwingelstein. Diagnostic des dfaillances : Thorie et pratique pour les systmes
industriels. Trait des nouvelles technologies Srie diagnostic et maintenance. Hermes.
1995.
[42] J. Ragot, D. Maquin, O. Adrot, A. Boatas, D. Jampi Dtection de dysfonctionnement
d'un systme amortisseur de vhicule automobile. 5me Congrs International
Pluridisciplinaire Qualit et Sret de Fonctionnement, Qualita 2003, Nancy, France,
18 - 20 mars 2003.
178
Rfrences
179