FR 2935154 A1 20100226 FR 0804641 A 20080819 La présente invention concerne le domaine du contrôle moteur et plus particulièrement le contrôle de combustion des moteurs à allumage par compression, tels que les moteurs diesels ou les moteur CAI (� Controlled Auto Ignition �). Le fonctionnement du moteur diesel est basé sur l'auto inflammation d'un mélange d'air, de gaz brûlés et de carburant. Le cycle du moteur se décompose en plusieurs phases (figure 1) : - Lors de l'admission (ADM), la soupape d'admission (SAOM) laisse entrer dans la chambre (CHB) le mélange d'air et de gaz brûlés. L'air est prélevé dans l'environnement extérieur du moteur. Les gaz brûlés sont prélevés dans la tubulure d'échappement (ECH) et ramené vers la tubulure d'admission (re-circulation de gaz brûlés EGR). Ce mélange gazeux rempli la chambre de combustion (CHB) et se mélange avec les gaz brûlés restés dans la chambre depuis la précédente combustion (EGR interne). - La soupape d'admission se ferme (IVC : intake valve closing). Le piston (PIS) comprime les gaz. - L'injecteur de carburant (INJ) injecte une masse précise de carburant. Après un bref délai d'auto inflammation, le mélange air, gaz brûlés, carburant s'enflamme créant ainsi une surpression qui repousse le piston. - Une fois le piston redescendu, la soupape d'échappement (SECH) s'ouvre, le mélange gazeux est alors évacué par la tubulure d'échappement. Après la fermeture de la soupape d'échappement, il reste une partie des gaz dans le cylindre (EGR interne). Les gaz évacués par la tubulure d'échappement sont divisés en deux. Une partie re-circule vers l'admission (EGR) tandis que le reste est évacué hors du moteur (via l'échappement). Le but du contrôle moteur est de garantir au conducteur le couple qu'il demande tout en minimisant le bruit et les émissions de polluants. II faut donc régler aussi finement que possible le pilotage des quantités des différents gaz et du carburant. État de la technique Pour réaliser un contrôle de combustion d'un moteur à allumage par compression, on connaît des méthodes permettant de déterminer le milieu de combustion à l'aide de capteurs montés sur le moteur. Le plus précis étant d'utiliser un capteur de pression dans la chambre de combustion. Une telle méthode est décrite par exemple dans le document suivant : J. Bengtsson, P. Strandh, R. Johansson, P. Tunestal and B. Johansson, "Control of Homogeneous Charge Compression Ignition (HCCI) Engine Dynamics", Proceeding of the 2004 American Control Conference, Boston, June 30 - July 2, 2004 Cependant l'utilisation de tels capteurs sur des véhicules de série n'est pas envisageable du fait du coût très important de tels capteurs. De plus, ces capteurs sont généralement sujets à des dérives relativement rapides. On connaît également des méthodes dans lesquelles les quantités et les phasages sont optimisés sur chaque point de fonctionnement statique (régime et couple) de manière à dégager une stratégie idéale à chaque point. On utilise donc un calibrage sur banc moteur pour dégager les valeurs optimales des deux jeux de données principaux : les masses d'air Mai,. et de gaz brûlés Me nécessaire dans la chambre de combustion, qu'on représentera parXair = (Mai,., Mgb) . La masse de carburant M1. ainsi que l'angle vilebrequin O auquel le carburant est injecté, qu'on représentera parXfue, = (Mf,9f) . Cependant ces stratégies se révèlent insuffisantes dans les phases transitoires. En effet, lors des phases de transition d'un point de fonctionnement à un autre (changement de vitesse du véhicule ou du profil de la route), le contrôle moteur supervise les différents actionneurs présents dans le moteur pour garantir le couple désiré tout en minimisant le bruit, les émissions de polluants et la consommation. Cela se traduit donc par le passage des valeurs des paramètres du point initial aux valeurs des paramètres du point final: X initial X final air au- x initial x final fuel Jijel Or, il existe dans le moteur deux échelles de temps. La plus rapide (50 Hz) correspond à l'ensemble du phénomène de combustion (1 cycle moteur). A cette échelle, on est capable de changer la stratégie de l'injection ()(fuel) pour piloter la combustion. C'est la boucle de fuel. (Cf. (b)). La plus lente (1 Hz) correspond à la dynamique des gaz dans les tubulures du moteur (admission, échappement, re-circulation de gaz brûlés). On ne peut changer plus vite la stratégie de cette boucle d'air ()(air ). (Cf. (a)) Avec les méthodes actuelles, les variables pilotées (Xair, Xji,el) n'arrivent donc pas en même temps à leurs valeurs de consigne à cause de cette différence de dynamique. Les objectifs en production de couple, consommation, polluants, bruit sont donc respectés dans (a) (b) les phases statiques (les deux boucles dynamiques sont stabilisées à leurs valeurs de référence), en revanche, si on ne prend pas de précaution dans les phases transitoires, une partie des paramètres atteignant presque instantanément sa valeur de consigne finale alors que l'autre partie est encore aux valeurs de consigne initiale, le moteur produit alors plus de polluants ou de bruit, et peut même dans certain cas s'éteindre. De plus, sans capteurs de pressions cylindre, les méthodes connues ne permettent pas de contrôler le phasage de la combustion pendant les phases transitoires. Or comme l'illustre les figures 2 et 3, ceci est insuffisant pour assurer le fonctionnement du moteur en transitoire. L'objet de l'invention concerne un procédé pour assurer le contrôle de la combustion d'un moteur à allumage par compression, notamment en phase transitoire, tout en s'affranchissant des problèmes de l'art antérieur. Le procédé y parvient, d'une part en contrôlant les deux boucles dynamiques de manière séparée, et d'autre part en corrigeant la valeur de référence de l'angle d'injection via un contrôle de l'angle CA,,. Le procédé selon l'invention Ainsi, l'invention concerne un procédé de contrôle de combustion d'un moteur à allumage par compression, dans lequel : - on détermine des valeurs de consigne de paramètres physiques liés à l'admission de comburant gazeux dans une chambre de combustion, ainsi qu'une valeur de consigne (9;,,zj )reg. d'un angle vilebrequin auquel un carburant doit être injecté dans la chambre de combustion, lesdites valeurs de consigne étant déterminées de façon à optimiser la combustion, - un système de contrôle moteur pilote des actionneurs de façon à ce que les valeurs desdits paramètres physiques soient égales auxdites valeurs de consigne. Le procédé comporte les étapes suivantes : - on corrige la valeur de consigne (9;,j )reg. avant que les paramètres physiques n'atteignent leurs valeurs de consigne, en calculant une correction d9,,f à appliquer à la valeur de consigne (9;,y)re(, de façon à ce qu'un angle vilebrequin CA,, , auquel y pour cent du carburant est consommé lors de la combustion, soit égale à une valeur de référence de cet angle pour une combustion optimisée ; le système de contrôle moteur pilote une injection de carburant dans la chambre de combustion lorsque l'angle vilebrequin est égal à la valeur corrigée de consigne (Oiw. ),.ef afin de conserver la combustion optimale. Selon un mode de réalisation, on détermine la correction d91,y. en tenant compte des différences entre des valeurs réelles des paramètres physiques et les valeurs de consignes des paramètres physiques, et en contrôlant l'angle vilebrequin CA, au moyen d'une modélisation de la combustion comportant un premier modèle pour modéliser un phénomène d'auto inflammation, et un second modèle pour modéliser un dégagement d'énergie en fonction de l'angle vilebrequin CA,.. On peut alors déterminer la correction d8;,u. en appliquant les étapes suivantes : - on détermine les valeurs réelles des paramètres physiques ; - on calcule les différences entre les valeurs réelles et les valeurs de consignes ; - on calcule des premiers coefficients de linéarisation en linéarisant au premier ordre le premier modèle ; - on calcule des seconds coefficients de linéarisation en linéarisant au premier ordre le second modèle, et en considérant que l'angle vilebrequin CAS, est égale à sa valeur de référence ; et on calcule la correction de,,,j au moyen d'une combinaison linéaire des différences, les coefficients de cette combinaison linéaire étant définis à partir des premiers et seconds coefficients de linéarisation. Selon l'invention, les paramètres physiques peuvent être choisis parmi au moins les paramètres suivants au moment de la fermeture soupape : pression dans la chambre de combustion (Pwc ), température dans la chambre de combustion (T,,,c ), rapport (X,,r ) entre une masse de gaz brûlés et une masse de gaz totale dans la chambre de combustion, et masse ) de gaz totale dans le cylindre. D'autres caractéristiques et avantages du procédé selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après. Présentation succincte des figures - la figure 1 montre les différentes phases d'un cycle de combustion d'un moteur à allumage par compression. - la figure 2 illustre une chronologie de combustion, en fonction de l'angle vilebrequin selon trois situations de contrôle de combustion : contrôle optimal (réalisé en phase stabilisée), contrôle actuel en phase transitoire sans contrôle du CA,, et contrôle souhaité en phase transitoire avec contrôle du CA,,. - la figure 3 illustre les trois courbes de dégagement d'énergie Q en fonction de l'angle vilebrequin 0 pour les trois situations décrites en figure 2. la figure 4 illustre un schéma de calcul de la correction d9inf de l'angle d'injection du carburant. Description détaillée du procédé Le procédé selon l'invention permet de contrôler le déroulement de la combustion d'un moteur à allumage par compression, en phase statique comme en phase transitoire. II comporte un contrôle séparé et indépendant de la boucle d'air (boucle lente) et de la boucle de fuel (boucle rapide), puis une adaptation de la dynamique de la boucle de fuel de manière à être cohérent avec la boucle d'air. Le procédé permet ainsi une adaptation de Xfue, pour conserver les caractéristiques de la combustion demandée (à travers la demande de couple du conducteur). On limite ainsi l'impact sur les émissions de polluants et de bruit tout en garantissant au conducteur le couple demandé. Selon ce procédé, le contrôle de la combustion d'un moteur à allumage par compression s'effectue en quatre étapes : 1- Détermination de valeurs de consiqne pour différents paramètres physiques. Lors des phases de transition d'un point de fonctionnement à un autre (changement de vitesse du véhicule ou du profil de la route), le contrôle moteur supervise les différents actionneurs présents dans le moteur pour garantir le couple désiré tout en minimisant le bruit, les émissions de polluants et la consommation. Cela se traduit donc par le passage des valeurs des paramètres Xair et Xfuel d'un point initial aux valeurs des paramètres d'un point final: X initial x final air air initial y final x fuel fuel Les valeurs finales sont définies de façon à optimiser la combustion, c'est-à-dire à brûler le maximum de carburant de façon à minimiser les émissions de polluants et la consommation tout en minimisant le bruit. Ces valeurs finales optimisant la combustion sont (a) (b) appelées valeurs de consigne. Le contrôle moteur est chargé de faire respecter ses valeurs consignes. Les paramètres physiques importants régulés par la boucle d'air qu'il faut connaître sont la pression, la température et la composition chimique des gaz dans la chambre. Idéalement, ces paramètres atteignent instantanément leur valeur consigne. En réalité, la lenteur de la boucle d'air fait qu'il existe une erreur sur ces paramètres Xajr entre leur valeur de consigne et leur valeur réelle pendant toute la phase de transition. En conséquence, les paramètres thermodynamiques (masse, pression, température et taux de gaz brûlés) de la charge gazeuse aspirée dans le cylindre sont différents de leur valeur de consigne. On adapte le contrôle de la boucle de fuel aux erreurs des paramètres suivants: P : La pression dans la chambre de combustion. Elle est fonction de l'angle vilebrequin 9 . T : La température dans la chambre de combustion. Elle est fonction de l'angle vilebrequin O . X : Le rapport entre la masse de gaz brûlés et la masse de gaz totale dans la chambre de combustion (paramètre entre 0 et 1). Ce taux n'évolue pas pendant tout le cycle et est donc égal au taux de gaz brûlés au moment où on ferme la soupape d'admission. M : La masse totale de gaz (air + gaz brûlés) enfermée dans le cylindre. La masse étant conservée, ce paramètre est constant. On distingue dans notre analyse la valeur de ces paramètres au moment de la fermeture soupape (IVC): - Pivc : La pression dans la chambre de combustion au moment de la fermeture soupape. - T,,,.c : La température dans la chambre de combustion au moment de la fermeture soupape. - X,,,c : Le rapport entre la masse de gaz brûlé et la masse de gaz totale dans la chambre de combustion au moment de la fermeture soupape. - M,,,c : La masse de gaz totale (air + gaz brûlés) dans le cylindre. Les valeurs de ces quatre paramètres sont déterminées en continue. Pour ce faire, nous supposons que la composition (X,vc) et la pression (P,vc) dans le cylindre à l'ivc sont les même que celles dans le collecteur d'admission où l'on dispose de mesures (par des Rmivc où R est la constante des gaz parfait (R=287) et M,vc est la masse aspirée par le cylindre qui est mesurée par un débitmètre. Pour ces quatre paramètres physiques liés à l'admission de comburant gazeux dans la chambre de combustion du moteur, les valeurs de consigne sont respectivement notées : Per.,T.e/ etXre. et M1ej Ces valeurs de consignes sont obtenues à partir d'une cartographie de consigne établie sur banc moteur : les valeurs de consigne de ces paramètres sont données par le point optimal cartographié au banc moteur (valeur que ces paramètres doivent atteindre). Ces valeurs de consigne sont déterminées de façon à optimiser la combustion. La valeur référence de ces paramètres, donnée par le point optimal cartographié au banc moteur (valeur que ces paramètres doivent atteindre) : La pression référence dans la chambre de combustion au moment de la fermeture soupape (obtenu pour le point optimal de référence). Trek La température référence dans la chambre de combustion au moment de la fermeture soupape (obtenu pour le point optimal de référence). La valeur référence du rapport entre la masse de gaz brûlé et la masse de gaz totale dans la chambre de combustion au moment de la fermeture soupape (obtenu pour le point optimal de référence). ù Mref. : La valeur référence de la masse de gaz enfermée dans le cylindre (obtenu pour le point optimal de référence). Ces paramètres sont liés par la relation des gaz parfait (PV=MRT) mais pour des raisons de simplicité, on n'explicite pas directement cette relation. Cela n'affecte en rien la méthode proposée. Selon l'invention, le paramètre de la stratégie de fuel important pour adapter la boucle de fuel est l'angle vilebrequin, noté , pour lequel du carburant est injecté. On note B1ei sa ug valeur de consigne. Cette valeur est également donnée par le point optimal cartographié au banc moteur. Elle correspond aux valeurs de consigne PYet,Tret , XYet et Mrer capteurs ou des estimateurs). Nous estimons T. par la loi des gaz parfait I. = PnrVwc30 2- Contrôle de la boucle d'air (boucle lente). Une fois déterminées les valeurs de consigne Pref,TYef , Xref et MYef , un système de contrôle moteur pilote des actionneurs de façon à ce que les valeurs des paramètres physiques P,vc, Tivc et X,vc et M,vc soient égales à ces valeurs de consigne (Fret ,Tref , X,.ef , Mref ). Idéalement, les quatre paramètres P,vc, Tivc et X,vc et M,,,c atteignent instantanément leurs valeurs de consigne Pef,TYef , Xref et Mref. En réalité, la lenteur de la boucle d'air fait qu'il existe une erreur sur ces paramètres entre la consigne et leur valeur réelle pendant toute la phase de transition. On adapte donc la boucle de fuel aux erreurs des ces paramètres. Pour ce faire on adapte la valeur de consigne (8;,,zj ),.ef de l'angle d'injection. 3- Correction de la valeur de consiqne de l'angle d'injection (B;nj )ref . Si le contrôle de la boucle d'air était parfait, les quatre paramètres P,,,c,T,,,c , M,,,c et X1,,c atteindraient leurs valeurs références Pref,Tref., MPef et Xref instantanément. En réalité, on ne peut rendre la boucle d'air aussi rapide que désiré. En phase transitoire, les paramètres P1 ,T,,,c , Mnr et X,,,c sont donc différents de leur valeur référence. Le contenu du cylindre au moment de la fermeture soupape est donc différent du contenu référence pour lequel a été cartographié la stratégie d'injection. Il faut donc prendre en compte les erreurs de ces paramètres au moment de la fermeture soupape (dP = P,vc - Pre f ), (dT = T,vc - Tref ), (dM = M,vc - Mref) et (dX = X,vc - Xref) pour modifier l'angle d'injection de fuel de manière à conserver une combustion aussi proche que la combustion référence. Selon l'invention, on contrôle le phasage de la combustion, en contrôlant l'angle CA, , qui correspond à l'angle vilebrequin auquel y% du carburant a été consommé par la combustion. Dans toute la suite, y est donc un réel quelconque entre 0 et 100. (y=0 correspond à un contrôle de début de combustion). Pour cela, on cherche un nouvel angle d'injection corrigé + pour que l'angle CA, soit à sa valeur de référence (dcA = CA û (CAr.) e~. = 0). On cherche donc dB;,y . tel que (voir figures 2 et 3 pour les situations): si il n'y a pas d'erreur, c'est-à-dire si tous les paramètres ont atteint leur valeur de référence ((dP, dT, dM, dX) = (0,0,0,0) ), on est exactement dans la situation du point de fonctionnement référence, on a donc dO,,, = 0 (situation OO ). Si les paramètres n'ont pas atteint leur valeur de référence ((dP,dT,dM,dX) ≠ (0,0,0,0) ), la durée du délai d'auto inflammation et la vitesse de la combustion ne seront pas identique à celles de la combustion référence. On a donc un déphasage de tout le processus de combustion et le CA n'est pas atteint à sa valeur référence (situation OO ). Pour contrebalancer les erreurs (dP,dT,dM,dX) ≠ (0,0,0,0) , on introduit donc une correction angulaire en] ≠ 0 sur l'angle d'injection pour avoir le même phasage CAy . La correction dO;,,i ≠ 0 va aussi entraîner un décalage de l'angle de début de combustion. L'angle de début de combustion devient alors (Osoireg + dOso, (situation OO ). La figure 2 illustre une chronologie de combustion selon trois situations. Pour chaque situation l'axe horizontal représente l'angle vilebrequin O. Sur ces axes sont repérés : la valeur de consigne (0;,j )rer de l'angle d'injection, la valeur de consigne (Osoe )rer de l'angle de combustion, l'angle CAy, la valeur effective (Osoc) de l'angle de combustion ; la valeur effective de l'angle d'injection. On note également le délai d'auto inflammation (DAI) qui est compris entre (°j,y )ref et (Osoc)ref , et la combustion se produit à partir de (Osoc)reg La figure 3 illustre les trois courbes de dégagement d'énergie Q en fonction de l'angle vilebrequin O pour les trois situations décrites précédemment (figure 2). On réalise ensuite une modélisation du système de combustion. Selon l'invention, cette modélisation comporte deux modèles : l'un permet de modéliser le phénomène d'auto inflammation, le second permet de modéliser le dégagement d'énergie. Selon un exemple particulier de réalisation, on peut utiliser les modèles décrits ci-après.30 Modèle d'auto inflammation Pour ce faire on se base dans un premier temps sur le modèle de � l'intégrale de Knock �. Ce modèle est décrit dans le document suivant : K. Swan, M.Shahbakhti and C.R. Koch, "Predicting Start of Combustion Using a modified Knock Integral Method for an HCCI Engine", in Proc. Of SAE Conference, 2006. Il est important de noter que la méthode de synthèse du contrôleur d'angle d'injection proposée est applicable à tout modèle de délai d'auto inflammation présentant la même forme intégrale que celui de Swan et al. Selon ce modèle, le début de la combustion ne se fait pas immédiatement après l'injection de fuel. Il y a un délai d'auto inflammation qui se modélise sous la forme de l'intégrale de Knock. Celui-ci permet, à partir des valeurs des paramètres P,T,X et , de déterminer l'angle de début de combustion 0,0, : f g(P(0), T (0), X, Ne) d O= 1 avec : i - g(P,T,X,Ne) = 6N ci + C X .P exp TA e 2 - A, C, , C2, n et TA sont des paramètres physiques fixes à calibrer. - Ne est le régime moteur. Classiquement, celui-ci est considéré constant. (Cette hypothèse étant justifié par le fait que les variations de régimes sont beaucoup plus lentes que les variations des autres paramètres considérés). 9 est l'angle du vilebrequin. Puis, on considère qu'avant le début de la combustion, le mélange de gaz est en compression adiabatique. On peut donc ramener facilement la connaissance de P(9), T(9) à celle de Pwc, I. et du volume de la chambre V(0), celui-ci étant parfaitement connu. De plus, le taux de gaz brûlés n'évolue pas tout au long de la phase de compression sans combustion X(9) = . On peut donc ramener l'intégrale de Knock à l'intégrale suivante : f ' f (Prvc Tivc Xrvc Ne99) de =1 10 Avec f une fonction entièrement connue est définie à l'annexe 3. Modèle de dégagement d'énergie Pour contrôler l'ensemble du déroulement de la combustion, via l'angle CAy, et non seulement le début de la combustion, il est nécessaire de considérer un modèle de dégagement d'énergie. On peut utiliser le modèle de dégagement d'énergie Diesel de Chmela (F.G. Chmela and G.C. Orthaber, "Rate of Heat Release Prediction for Direct Injection Diesel Engines Based on Purely Mixing Controlled Combustion", in Proc. Of SAE Conference, 1999). Ce modèle donne le dégagement d'énergie lié à la combustion Diesel sous la forme d'une équation différentielle: d B Cmod e Mù Q exp Craie 3 v k(e) QLHV / ~Vcy! (e) avec : Mt. est la masse de fuel totale injectée. QLHV est l'énergie massique disponible dans le carburant. Q est l'énergie libérée par la combustion - V (9) est le volume du cylindre. - k(9) est la densité d'énergie cinétique turbulente présente dans le cylindre. Elle satisfait une certaine équation différentielle présentée en annexe 4. - Cmode et Craie sont deux constantes. - Ne le régime moteur. On apporte une modification à ce modèle, afin de prendre en compte la présence de gaz brûlés re-circulés dans la chambre de combustion. Pour cela, on ajoute un terme en fonction du taux de gaz brûlé X,,,c présent avant la combustion sous la forme: exp Craie h(X1) (1) Vcy / \1k(e,Bi,y.,m ,mlvc'Ne) A'LmodeeX P QLHV Cra,e (e) h(Xwc)dO= zoo dx =ûln(1û Y ) (1û x) 100 (2) Estimation de la correction de la valeur de consiqne de l'anqle d'infection (9;,�rer Le calcul de la correction se fait en deux étapes identiques appliquées aux deux modèles (auto inflammation et dégagement d'énergie). On linéarise les modèles en Pitic,TIvc , Mrvc , X wc , O,,, et 95oe aux alentours de leurs valeurs références en introduisant des écarts dP, dT, dM, dX , d 9;äi et d8soe . De cette manière on obtient deux équations permettant d'obtenir d9,.,,i en fonction des erreurs de la boucle d'air dP, dT, dM, dX . Première étape: En linéarisant l'intégrale de Knock au premier ordre, on peut simplement lier la correction dO1, à la variation d'angle de début de combustion desoe (on fait de ce fait l'hypothèse de petites perturbations suivante : et,. - Deuxième étape: On considère de même l'équation (2) comme une équation implicite en (M,f,c , X,,.c , et 95oe ). On la linéarise autour de (,Tre1 , Mref , Xret , (9;ä1)ref et (eSOe )ref) ce qui permet d'introduire les écarts (dM, dX , de,y. et dOroc) (on fait de ce fait l'hypothèse de petites perturbations suivante : - . J ( ref,Tef,Xref,Ne,e) = Vrvc V(9) V(0) = . Pre 5 ANNEXE 1 Expression des coefficients de linéarisation (aP>ar,ax) On a: dOSOe=a;d9;+apdP+a,dT+axdX (3) Avec: (per,Tef,X,.el.,Ne,8) d8 ). / J (Prel ,Tref,Xref,Ne,(esoc)ref ) n 1 Pef / (Pref,Tre(,Xref,Ne,(esoc)reO je",,,f ÔT ))éf ~f (PreJ,Tel.,XreJ,Ne,O) dO f(Pref , T ef , )(rd , Ne , (esoc )ref ) TA (~8., ),... Tivc Î(Per,Tref,XYef,Ne,9).v' '(9) de J (Pref,Tef,Xref,Ne, eso,)ref ) 9,),,., a.Î (Pei Tef,X,.ef,Ne,O) dO J (P ef , Tief , Xref , Ne , (esoc )r CZ Aä.),,f(Pf,Tef,Xref ,Ne,0).ln(v(9)) n+ TA v'-T(6) dO dy(T,Xref) te,,; ),,., Tm_ dX rer f (Pref,Tef ,Xref ,Ne,(ein; /ref) ainj = (P'ef , Tref X ref , Ne , (esoc )ref ) aP a,. ax ) ( Cl + C2 X ref J (PrefTref,Xref,Ne,(esoc)ref ) Nroc = ùexr ("rate Vvcy(((ûsoc)ref ) h�Xrvc)1n(1 _ y ) 100 V k((Bsoc )ref, )ref' Mf ' (MI VC)ref ) ANNEXE 2 Expression des coefficients de linéarisation (/'inj'/3M,PX'/3soc) Ona: +/3MdM+/XdX+flacdOsa, =0 (4) Jk(6', (9inj )ref' Mf' (MIVC)ref ) VVV Cm"' hX ( ) ),..r :VI 1 de J QLHV a8^l C ç no, / 2 exr C ~k(e'(B^j )ref~~f(MlvC)tef)\ e,ss mode h(XIVC )C d9 rate rate 6NeQLIIV 3 IV r(e) i ,jk(û' Br^j ref,Mf ,(MIVC)ref) aexp Cro,e '~ _ ("mode hc~.L, v,,(6) I dB QLHV h(Xlvc)~eä l aMlvc '(MIVC)ref ) \ik(a'(anJ)ref'Mf '(MIVC)ref ) C �X, mode VC ~("ra,e ex~ ("ra,e 3 QLHV 2(Mlvc+Mf) ,., VV~,i(B) 3V,,i(B) /~ ("mode .).r,ah(Xlve)eXP ("~k(8,(Bnj)ref,Mf,(M1VC)ref Px d9 QLHV �-�"r aXlvc ra,e \1Vy,;(a) V ko(Mf,N) /~ 6Nr~(O a,,,) "(B,t9nj ,Mf ,M/vc) C,rb Mlvc+Mf e avec Rnj aexr Cra,e d9 et ANNEXE 3 Intégrale De Knock Le début de la combustion ne se fait pas immédiatement après l'injection de fuel. Il y a un délai d'auto inflammation qui se modélise sous la forme de l'intégrale de Knock. Celle-ci relie de manière implicite l'angle d'injection et l'angle de début de combustion : 4),,,, g(P(0),T(0),X,Ne) dO =1 Avec g(P,T,X,Ne) = .P .exp Ne C, C,X - 9 est l'angle du vilebrequin. - A, C, , C2 ,n et TA sont des paramètres fixes et connus du modèle. ù Ne le régime moteur. Exemple de jeu de paramètre: Paramètre A C, C2 n TA Unité (hPa)-" s' 1 1 1 °K Valeur 0.22 1 10.5 1.13 1732 En considérant que le mélange gazeux est parfait et est soumis à une compression adiabatique, on a, avant le début de la combustion : - P(0).V (0)r = cste P(B) = P,,,c .vr (0) - T(0).V(0)r-' = cste = T(0) = T,vc.vr-1(0) Avec: - v(0) = V(9) : fonction du volume entièrement connue. - v,vc : le volume de la chambre à la fermeture des soupapes. ù V(0) : le volume de la chambre en fonction de l'angle du vilebrequin. T " 4 _ Tl y(X) : paramètre de compression adiabatique. Il dépend notamment de la composition chimique X . X : La composition n'évolue pas dans le cylindre. Celle-ci est donc égale à la composition au moment de la fermeture soupape admission: X = X,vc On obtient donc au final l'équation générale suivante : Î(P,vc,Tivc,Xwc,Ne,0) d9 =1 f avec 1 A { ( T / J (Pive T1 , X1 vc, Ne, 9) = PIVC~~ Ur(xly ).n (8) PXp ù .v' 7(Xnc~ (0) 6Ne C, +CZX,vc \ Tivc \ / ANNEXE 4 Obtention de la densité d'énergie cinétique turbulente En combustion diesel à injection directe, il est courant de considérer que l'énergie cinétique turbulente est majoritairement due à l'injection de carburant (de l'ordre de 95%). La densité d'énergie cinétique turbulente k a donc une évolution gouvernée par l'équation différentielle suivante: dk Cd,.ss k + C,,,rb dL:cin de 6 Ne M r+ M;vc dB dEc,n = CD (Ni, ).Dcarb3 L dB avec: CD(Ne) une fonction qui dépend du régime (donc considérée comme constante). C,,,,.b est une constante. Dcarb est le débit de carburant instantané. Il dépend de la masse injectée (qui détermine la durée et le profil du débit de carburant). De plus, si on décale l'injection dans le temps, on ne fait que décaler ce profil. Le débit de carburant ne dépend donc que de Mf et de l'intervalle de temps entre l'instant présent (angle vilebrequin B) et B-B..` l'injection (angle vilebrequin Oh,/ ). Ce qui donne: Dcarb = Dcarb M f, nf Ne ~ Exemple de jeu de paramètre: Paramètre Cmrb Cdiss Co(Ne) Unité [-] s' mzkg-z Valeur 1 200 9e 6 /Ne(rpm) L'équation différentielle du premier ordre peut s'intégrer avec la méthode de variation de la constante: __ k0 (M,,Ne) k(6) Mf +Mwc \3 c,,,, Z ù e,,, avec ko (M f , 6, , Ne) = Ctärb f CD (Ne )Dcarb M f , ,, e Ne dz Ne i L'intégrale ci dessus a un intégrande non nul uniquement pendant l'injection (en dehors, le débit de carburant est nul). Si on appelle AO,,,; la durée de l'injection, on a donc On s'intéresse à l'expression de k pendant la combustion (car k intervient uniquement dans l'expression du dégagement d'énergie). Or on a fait l'hypothèse que l'injection était finie lorsque la combustion commence. Donc VO>Bsoc>+49D,;, ko(Mf,0,,Ne)ùCn,rbCD(Ne)•~'l 'Dcarb Mf, ; eNe dz Ne ) / 3 c Bq Z z =CIurbCD(Ne).f Dcarb Mf,ù e Ne dZ=este Ne / 10 Donc, pendant la combustion, ko dépend de la manière dont le carburant est injecté (profil de Dcnrb) mais ne dépend pas directement de 8,,; ni de O . Donc: k 8 = k o (M f , Ne) e N O M f + M1,, VO > +Ako(M f,O,Brä;,Ne) =CturbCD(Ne)• Dcarb Mf ' e Ne dZ Ne 21 (B-Br., ) B-0,,,, ) avec ko (M f , Ne) = Ca,rbCD (Ne z e^'e dz 33 15 ANNEXE 5 Dégagement d'énergie Pendant la combustion, le dégagement d'énergie se fait selon la loi: dQ _ d0 Cmode .Jk(0) Crate QLHV / Q \ exp h(Xwc ) Exemple de jeu de paramètre Paramètre Cmod e QLHV Crase Unité kJ.kg-'deg-' kJ.kg-' s Valeur 350 42789 1,5e-3 Si on appelle x la fraction de masse de fuel consommée. On a la relation classique suivante, qui lie l'énergie libérée Q , l'énergie chimique totale contenue dans le carburant MfQLHV et x : Q = xMfQLHV Au final en reportant la deuxième équation dans la première, on a la relation: dx = Cùde (1û x)exp QLHVVk(0) rase h(X,vc)d0 C 3V,I(e) a dx = Cmode ex C. ~k(0) p rate (1û x) QLHV vv,yl (0) Le CAy, est l'angle vilebrequin auquel y% du carburant a été consommé. Il correspond donc à x = 100 . On peut donc intégrer séparément les deux membres de la dernière équation 15 pour obtenir une relation implicite entre le CAti, et le début de combustion Bso~ : h(Xlvc)dd zoo dx = ûln(1û Y ) = h(Xrvc (1û x) 100 Cùde expC ~k(0) v ÇAr rase de QLHV V V i! (e) REVENDICATIONS1. Procédé de contrôle de combustion d'un moteur à allumage par compression, dans lequel : - on détermine des valeurs de consigne de paramètres physiques liés à l'admission de comburant gazeux dans une chambre de combustion, ainsi qu'une valeur de consigne (8;,,i),.ef d'un angle vilebrequin auquel un carburant doit être injecté dans la chambre de combustion, lesdites valeurs de consigne étant déterminées de façon à optimiser la combustion, - un système de contrôle moteur pilote des actionneurs de façon à ce que les valeurs desdits paramètres physiques soient égales auxdites valeurs de consigne, caractérisé en ce que le procédé comporte les étapes suivantes : - on corrige la valeur de consigne (O,,,f. )ref. avant que lesdits paramètres physiques n'atteignent leurs valeurs de consigne, en calculant une correction d9;,,/ à appliquer à la dite valeur de consigne (O,. )ref , de façon à ce qu'un angle vilebrequin CA). , auquel y pour cent du carburant est consommé lors de la combustion, soit égale à une valeur de référence de cet angle pour une combustion optimisée ; ledit système de contrôle moteur pilote une injection de carburant dans la chambre de combustion lorsque l'angle vilebrequin est égal à ladite valeur corrigée de consigne (8;,y )7ef. afin de conserver la combustion optimale. 2. Procédé selon la revendication 1, dans lequel on détermine ladite correction dO,,i en tenant compte des différences entre des valeurs réelles desdits paramètres physiques et lesdites valeurs de consignes desdits paramètres physiques, et en contrôlant l'angle vilebrequin CA), au moyen d'une modélisation de la combustion comportant un premier modèle pour modéliser un phénomène d'auto inflammation, et un second modèle pour modéliser un dégagement d'énergie en fonction de l'angle vilebrequin CAti 3. Procédé selon la revendication 2, dans lequel on détermine ladite correction dO1,v en appliquant les étapes suivantes : - on détermine lesdites valeurs réelles desdits paramètres physiques ; - on calcule lesdites différences entre lesdites valeurs réelles et lesdites valeurs de consignes ;on calcule des premiers coefficients de linéarisation en linéarisant au premier ordre ledit premier modèle ; on calcule des seconds coefficients de linéarisation en linéarisant au premier ordre ledit second modèle, et en considérant que l'angle vilebrequin CA, est égale à ladite valeur de référence ; et - on calcule la correction d0,,,i au moyen d'une combinaison linéaire desdites différences, les coefficients de ladite combinaison linéaire étant définis à partir des premiers et seconds coefficients de linéarisation. 4. Procédé selon l'une des revendications précédentes, dans lequel lesdits paramètres physiques sont choisis parmi au moins les paramètres suivants au moment de la fermeture soupape : pression dans la chambre de combustion (Pnr ), température dans la chambre de combustion (Tlvc ), rapport (X,vc) entre une masse de gaz brûlés et une masse de gaz totale dans la chambre de combustion, et masse (M1,c) de gaz totale dans le cylindre.