Brevet FR2905411A1 — publié le 20080307

DOMAINE TECHNIQUE GENERAL L'invention concerne un procédé de détermination d'un début d'une combustion dans une chambre à combustion d'au moins un cylindre d'un moteur à combustion interne, le cylindre comportant un piston mobile entre un point mort haut (PMH) et un point mort bas (PMB), le piston étant en outre relié à un vilebrequin monté à rotation dans le moteur. L'invention concerne également un véhicule mettant en oeuvre le 10 procédé. ETAT DE L'ART Les normes européennes antipollution deviennent de plus en plus sévères sur les moteurs d'automobiles. II est important pour les constructeurs automobiles de réduire les émissions polluantes de leurs 15 moteurs, sans négliger leur consommation en carburant et leurs performances. Or la phase prépondérante du cycle quatre temps sur ces émissions est la phase de combustion. II est donc important de bien connaître cette phase et ses caractéristiques. 20 L'une de ces caractéristiques est l'instant où débute la combustion, c'est-à-dire l'instant où les réactions chimiques irréversibles commencent, correspondant à la transformation de l'énergie chimique en énergie mécanique. Il est un paramètre d'ordre un sur le rendement de cycle thermodynamique et donc sur le couple. Il agit aussi sur le rendement de 25 combustion. En pratique, plutôt que de parler d'instant de début de combustion, on parle d'angle de début de combustion. En effet, plutôt que de prendre le temps comme repère, on préfère prendre l'angle du vilebrequin. Le zéro du repère est par exemple pris à la position où le piston est au point mort 30 haut (PMH). Sur un moteur fonctionnant correctement, le début de combustion se situe au voisinage du PMH, donc ainsi au voisinage de 0 2905411 2 vilebrequin. C'est ce repère et cette origine que nous prendrons dans la suite de ce document. Aujourd'hui, le début de combustion est rarement estimé en ligne sur un moteur implanté dans un véhicule. Lorsque cette détection est effectuée, 5 elle l'est grâce à des méthodes numériques basées sur des principes thermodynamiques ou de traitement numérique du signal, difficiles à mettre en oeuvre. WO 02/48522 Al divulgue un procédé de détermination de début de combustion ( Start Of Combustion ou SOC selon la terminologie anglo- saxonne généralement utilisé par l'homme du métier) qui utilise la dérivée de la pression de combustion pour détecter le SOC, par comparaison à un seuil. C'est une méthode simple, mais peu robuste. L'angle de début de combustion (SOC) est difficilement identifiable à l'aide de cette méthode. En effet, la dérivée de la pression comporte un pic représentatif de la brusque élévation de pression, trahissant le début de la combustion. II est alors possible par l'application d'une valeur seuil de détecter un instant représentatif du commencement de la phase de combustion. Cependant, l'application d'une valeur de seuil est rendue difficile par 20 plusieurs phénomènes. Tout d'abord, l'opération de dérivation du signal de pression entraîne l'amplification du bruit de mesure. Ce bruit est présent sur le signal pendant toute la phase de compression et contraint le choix de la valeur de seuil. En outre, l'élévation rapide de la pression due à la simple compression 25 des gaz enfermés dans la chambre de combustion engendre une valeur de dérivée non négligeable, et impose des contraintes supplémentaires sur le choix de la valeur de seuil. Enfin, cette élévation normale de la pression pendant la phase de compression dépend du point de fonctionnement du moteur (pression et 30 température des gaz admis, température des parois de la chambre, masse admise, etc.). De ce fait, lors de combustions peu vives, le pic de la dérivée 2905411 3 dû à la combustion peut être inférieur au pic de la dérivée dû à la compression. Un autre procédé possible consiste à étudier le comportement de la dérivée seconde de la pression cylindre pendant la phase de combustion. 5 Cependant, la double dérivation engendre un niveau de bruit sur le signal très important. US 5 359 883 divulgue un autre procédé, s'appuyant quant à lui sur le calcul du rapport des fonctions logarithmes de pression et de volume. Le calcul de ce rapport permet de détecter un changement de comportement 10 du coefficient polytropique y. Ce coefficient est à peu près constant pendant la phase de compression du cycle thermodynamique du moteur. Puis, pendant la phase de combustion, le comportement du coefficient polytropique change, ce qui se traduit par un changement de pente du tracé du cycle. Cette rupture de pente peut être détectée par exemple en 15 appliquant un seuil sur la grandeur : d /log p\ dt , logV~ Cependant, l'opération engendre aussi un bruit de mesure important. US 2003/0145829 divulgue un procédé selon lequel on propose de détecter le commencement de la phase de combustion en étudiant le profil de la pression de cylindre, auquel il est appliqué une transformation en ondelettes. Ce procédé est particulièrement lourd à mettre en oeuvre et nécessite un temps de calcul élevé, ce qui rend très difficile son utilisation à bord d'un véhicule. D'autres méthodes détectent le début de combustion, non pas à l'aide d'un capteur de pression mais grâce à un accéléromètre, comme par exemple dans WO 02/25237. Ces méthodes ne sont pas robustes ni précises, et nécessitent un traitement de signal complexe. En effet, un accéléromètre placé sur le moteur capte toute sorte de signaux parasites (les ouvertures et fermetures de soupapes par exemple). 2905411 4 PRESENTATION DE L'INVENTION L'invention propose de pallier au moins un des inconvénients précités. A cet effet, on propose selon l'invention un procédé de détermination d'un début d'une combustion dans une chambre à combustion d'au moins 5 un cylindre d'un moteur à combustion interne, le cylindre comportant un piston mobile entre un point mort haut (PMH) et un point mort bas (PMB), le piston étant en outre relié à un vilebrequin monté à rotation dans le moteur, caractérisé en ce qu'il comporte les étapes consistant, grâce à un boîtier électronique, à : 10 a) calculer au moins deux valeurs d'une variable (Ce , Co_Ao ) caractéristique de la combustion dans la chambre (3) de combustion, la variable caractéristique étant notamment une fonction de l'énergie dégagée dans la chambre, une première valeur (CB_4e) étant calculée pour un premier angle (9-ie) de rotation du vilebrequin ou un premier instant (t-nt), 15 une deuxième valeur (Co) étant calculée pour un deuxième angle de rotation (9) du vilebrequin ou un deuxième instant (t); b) déterminer si la deuxième valeur (Co) de la variable caractéristique dépasse un premier seuil (Si) ; et c) déterminer si la différence entre la deuxième valeur (C9) de la 20 variable caractéristique et la première valeur (CB_OB) de la variable caractéristique est inférieure à un deuxième seuil (S2), le début de la combustion étant déterminé si les conditions des étapes b) et c) sont vérifiées. L'invention est avantageusement complétée par les caractéristiques 25 suivantes, prises seules ou en une quelconque combinaison techniquement possible : - le deuxième angle (9) ou le deuxième instant (t) est proche du premier angle (9-LA) ou du premier instant (t-nt) ; - le calcul de la variable caractéristique C est effectué par : 30 C= 1 dQ F dB 2905411 5 où Q est l'énergie dans la chambre, 0 est l'angle du vilebrequin, et F est une fonction de 0 telle que F est maximale au point mort bas, minimale au point mort haut, monotone entre les deux, et toujours 5 strictement positive ; - la fonction F est donnée par : F(0) = Fo +K.101 où K et Fo sont des constantes ; - le calcul de la variable caractéristique C est effectué par : 10 C=1dQ V d9 où V représente le volume de la chambre pour une valeur 0 de l'angle, Q est l'énergie dans la chambre, et 0 est l'angle du vilebrequin ; 15 -ona: 1 V + + Y p d6 Yù1 m deTdOV d0/ yù1 d9 Avec : Q la quantité d'énergie présente dans les gaz de la chambre de combustion ; 20 0 l'angle du vilebrequin ; y le rapport des chaleurs spécifiques (Cp et ci') ; V le volume de la chambre de combustion ; p la pression dans la chambre de combustion ; m la masse de gaz enfermés dans la chambre; et 25 T la température dans la chambre de combustion ; -on a: dQ V p dT " dV d9 yù1 T de+pd8 Avec : 2905411 6 Q la quantité d'énergie présente dans les gaz de la chambre de combustion ; 0 l'angle du vilebrequin ; y le rapport des chaleurs spécifiques (Cp et CV) ; 5 V le volume de la chambre de combustion ; p la pression dans la chambre de combustion ; et T la température dans la chambre de combustion ; - Te = To PeVe PoVo avec To, po,Vo les température, pression et volume initiaux ; - Te = Te-oe 1 P0-oe + Ve-oe _1 \ PO Vo ~ 10 - C ù 1 po Te ùTo,o + Po Ve ù Ve-oe 9 y-1,Te AO ) Vo 08 - le volume est tiré d'une table de valeurs donnant le volume V et sa dérivée temporelle pour chaque valeur d'angle du vilebrequin ; - le calcul du volume vrai V est remplacé par une approximation Vp 15 suivant une loi parabolique donnée par l'expression : Vp(0) = Vo + KpO2 où: Vo est constant et est le volume de la chambre au point mort haut (PMH), Kp un coefficient à ajuster de telle façon que la courbe de Vp soit au 20 plus près de la courbe de V au voisinage du point mort haut ; - le premier seuil (Si) et le deuxième seuil (S2) sont fixes ou issus d'une cartographie fonction de la charge du moteur ou de la quantité de carburant injectée ; - le procédé comporte, dans le cas d'un moteur à au moins deux 25 injections de quantités injectées différentes dans la chambre, une étape consistant à détecter séparément le début de combustion de chacune des injections ; 2905411 7 - la détection de chaque début de combustion s'effectue en multipliant les détections et en appliquant des seuils différents pour chaque détection ; - la détection de chaque début de combustion s'effectue en multipliant les détections, et en déterminant si la différence entre la deuxième valeur 5 (Co) de la variable caractéristique et la première valeur (CB_ào) de la variable caractéristique pour l'injection de plus grande quantité est supérieure à un deuxième seuil (S2) ; - le deuxième seuil (S2) est identique pour chaque détection ; et - la détection de chaque début de combustion s'effectue en appliquant 10 deux seuils pour l'injection de plus faible quantité, et un seuil unique à l'injection de plus grande quantité. L'invention concerne également un véhicule mettant en oeuvre le procédé. L'invention présente de nombreux avantages. 15 Elle permet l'estimation en temps réel, sur un moteur essence ou Diesel, de début de combustion (SOC) et notamment de l'angle de début de combustion. L'invention a une meilleure robustesse et une meilleure précision de détection. 20 L'invention est simple à mettre en oeuvre. Une détermination précise du début de combustion présente elle aussi des avantages. En effet, la détection de début de combustion peut servir à des stratégies d'estimation ou de régulation du bruit, d'optimisation du couple 25 par régulation du placement de la combustion dans le cycle thermodynamique d'un moteur, de recalage de la caractéristique des injecteurs pour ajuster les quantités injectées. Cette détection peut entrer dans la composition de système d'aide à la mise au point des moteurs sur bancs, mais il sera surtout très utile dans les systèmes de contrôle moteur 30 embarqués. Concernant les stratégies d'estimation ou de régulation du bruit, l'invention est particulièrement utile et avantageuse, car elle résout le 2905411 8 problème du fenêtrage, surtout lors d'une combustion homogène (moteurs dits HCCI). En effet, les systèmes de mesure du bruit de combustion basés sur l'observation de la dérivée de la pression dans la chambre (divulgué par exemple dans FR 2 857 410) sont contraints d'observer la pression dans 5 une fenêtre temporelle ouverte au moment de l'injection de carburant, afin de ne pas confondre le gradient de pression dû à la combustion avec le gradient de pression naturel dû à la compression du gaz dans la chambre. Ils ne sont pas utilisables sur des moteurs HCCI car dans ce type de moteur, l'injection est très précoce dans le cycle thermodynamique et ne 10 peut pas servir à déclencher l'ouverture de cette fenêtre d'observation. PRESENTATION DES FIGURES D'autres caractéristiques, buts et avantages de l'invention ressortiront de la description qui suit, qui est purement illustrative et non limitative, et qui doit être lue en regard des dessins annexés sur lesquels : 15 - la figure 1 représente schématiquement un moteur d'un véhicule mettant en oeuvre l'invention; la figure 2 représente schématiquement les étapes principales d'un procédé selon l'invention pour calculer le début de la combustion, et par conséquent l'état de la combustion dans une chambre de combustion ; 20 la figure 3 représente une fonction triangle ; - la figure 4 représente schématiquement l'effet d'une fonction selon la figure 3 dans le calcul d'une variable caractéristique selon l'invention ; - la figure 5 représente schématiquement l'effet d'une fonction en IN dans le calcul de la variable caractéristique selon l'invention ; 25 - la figure 6 représente un exemple d'approximation du volume d'une chambre de combustion ; - les figures 7 et 8 représentent schématiquement l'effet d'une fonction selon la figure 6 dans le calcul d'une variable caractéristique selon l'invention par rapport à un calcul effectué avec un vrai volume ; 30 - la figure 9 représente schématiquement les étapes principales d'un procédé selon l'invention pour calculer le début de la combustion avec des approximations avantageuses ; 2905411 9 - la figure 10 représente schématiquement un montage logique mettant en oeuvre l'invention ; et - la figure 11 représente schématiquement une mise en oeuvre de l'invention avec au moins deux injections de quantités différentes. 5 Sur l'ensemble des figures, les éléments similaires portent des références numériques identiques. DESCRIPTION DETAILLEE La figure 1 représente schématiquement un moteur d'un véhicule mettant en oeuvre l'invention. Le moteur à combustion interne est du type 10 Diesel ou essence. Le moteur comporte au moins un cylindre et un vilebrequin 8 monté à rotation dans le moteur. Le cylindre comporte une chambre de combustion 3 et un piston 5. La chambre 3 de combustion est comprise dans la partie haute du 15 cylindre, dans un espace compris entre le piston 5, un dispositif 2 d'admission de gaz et un dispositif 1 d'échappement de gaz. Le piston 5 est relié au vilebrequin 8 par une bielle 6, et est mobile entre un point mort haut (PMH) et un point mort bas (PMB). Le moteur comporte en outre un détecteur 9 de l'angle du vilebrequin 20 8 et un capteur 4 de la pression dans la chambre 3 de combustion. Le véhicule comporte de plus un boîtier électronique 7 de commande du système. Le boîtier 7 comporte un calculateur 71, un différenciateur 74, deux comparateurs 72 et 73 et une porte logique ET (AND) 75. 25 La figure 2 représente schématiquement les étapes principales d'un exemple d'un procédé selon l'invention pour calculer le début de la combustion, et par conséquent l'état de la combustion dans une chambre de combustion. Le principe d'un procédé est l'analyse d'une variable calculée C, 30 image de la combustion, au moyen de critères de discrimination. A partir de la valeur de pression p régnant dans la chambre 3 de combustion, et de la valeur 0 de la position du vilebrequin 8, on calcule en 2905411 10 permanence des valeurs d'une variable C (étape a , C qui est le plus souvent basée sur le dégagement d'énergie dans la chambre . d8 Ces valeurs sont analysées en permanence, pour un angle 0 compris entre -180 et +180 , ou dans un intervalle plus réduit (étapes b) et c . 5 Tant que certaines conditions ne sont pas remplies, on estime que la combustion n'a pas commencé. Dès qu'elles le sont, on estime que la combustion vient de débuter, et on repère la valeur de l'angle vilebrequin. On connaît donc l'instant de début de combustion. Ainsi, un procédé selon l'invention comporte les étapes consistant à, 10 grâce au boîtier électronique 7, à : a) calculer au moins deux valeurs Co et CB_49 de la variable C caractéristique de la combustion dans la chambre 3 de combustion, une première valeur CB_49 étant calculée pour un premier angle e-t e de rotation du vilebrequin 8 ou un premier instant t-At, une deuxième valeur Co étant 15 calculée pour un deuxième angle de rotation 0 du vilebrequin 8 ou un deuxième instant t; b) déterminer si la deuxième valeur Co de la variable C caractéristique dépasse un premier seuil SI ; et c) déterminer si la différence entre la deuxième valeur Co de la 20 variable C caractéristique et la première valeur CB_49 de la variable C caractéristique est inférieure à un deuxième seuil S2, le début de la combustion étant déterminé si les conditions des étapes b) et c) sont vérifiées. Bien entendu, le deuxième angle 0 ou le deuxième instant t est proche 25 du premier angle 0-A8 ou du premier instant t-Ot (par exemple, le deuxième angle est séparé du premier angle de 1 ou moins). PREMIER MODE DE MISE EN OEUVRE Le calcul de la variable caractéristique C peut être effectué par l'expression : 2905411 11 1 dQ Fde où Q est l'énergie dans la chambre, 0 est l'angle du vilebrequin, et F est une fonction de 0 telle que F est maximale au point mort bas, 5 minimale au point mort haut, monotone entre les deux, et toujours strictement positive. Une fonction F "triangle" possède de telles propriétés, et peut être réalisée très simplement. II suffit de faire : F=Fo+K.10 10 où K et Fo sont des constantes, comme cela est représenté à la figure 3. dQ En pratique, le calcul du dégagement d'énergie d9 donne une information bruitée. L'effet d'une telle fonction F est d'atténuer le signal de dégagement d'énergie dans les zones où il est le plus bruité (du coté du 15 PMB), par rapport aux zones où il est susceptible de contenir l'information qui nous intéresse (du coté du PMH), comme on le voit sur la figure 4. Cette opération augmente donc le rapport signal sur bruit et diminue les contraintes sur le filtrage de la pression à mettre en place. DEUXIEME MODE DE MISE EN OEUVRE 20 Le calcul de la variable caractéristique C peut également être effectué par: 1 dQ C= Vde où V représente le volume de la chambre pour une valeur 0 de l'angle, 25 Q est l'énergie dans la chambre, et 0 est l'angle du vilebrequin. 1 Comme le volume varie cycliquement avec 0, la multiplication par V a pour avantage de diminuer l'amplitude du bruit sur le dégagement 2905411 12 d'énergie, là où il est le plus fort (vers les points morts bas). Ce calcul amplifie aussi le signal là ou il est utile pour la détection du SOC (au point mort haut). Cette opération augmente le rapport signal sur bruit et diminue les contraintes sur le filtrage de la pression à mettre en place, comme le 5 montre la figure 5. dQ En pratique, le calcul de do est complexe et nécessite trop de puissance de calcul pour être réalisé tel quel dans un système électronique embarqué de contrôle moteur. Mais la division par le volume est intéressante car elle permet de 10 réaliser des simplifications. Dans un premier temps, on peut prendre : dQ 1 r p dm p dT p dV y dV V --+ + dO yù1 m dO TdO V dO yù1pdO Avec : Q la quantité d'énergie présente dans les gaz de la chambre de 15 combustion ; 0 l'angle du vilebrequin ; y le rapport des chaleurs spécifiques (cp et C ) ; V le volume de la chambre de combustion ; p la pression dans la chambre de combustion ; 20 m la masse de gaz enfermés dans la chambre; et T la température dans la chambre de combustion. Notons que l'expression pourrait être écrite en remplaçant 0 par le temps t. on a alors : d0=wdt. 25 Cependant il est plus intéressant d'utiliser l'angle 0 car le calcul devient alors indépendant de la vitesse de rotation w du moteur. dQ De plus, l'expression do prise comme point de départ est celle faisant intervenir la dérivée angulaire de la température et non pas la dérivée de la 2905411 13 pression. Ceci encore une fois par souci de minimisation du bruit, car la dérivée de la pression introduit un bruit important dans le résultat. En toute rigueur, V, p, T, y, m, sont des grandeurs variables en fonction de l'angle 0, et seuls la pression p et l'angle 0 sont des valeurs 5 mesurées par le moteur. Il faut donc définir ou estimer V, T, y, m. Concernant la masse m, nous faisons ici l'hypothèse qu'elle est constante (masse injectée négligée). On peut alors considérer que dQ V p dT dV dO y -1 T de + p de 10 Avec : Q la quantité d'énergie présente dans les gaz de la chambre de combustion ; 0 l'angle du vilebrequin ; y le rapport des chaleurs spécifiques (cp et ci' ) ; 15 V le volume de la chambre de combustion ; p la pression dans la chambre de combustion ; et T la température dans la chambre de combustion. Concernant le rapport des chaleurs spécifiques y, nous considérons qu'il est constant avant la combustion, et qu'il le reste encore dans les 20 premiers instants de la combustion. Nous le prenons par exemple égal à 1,4 Pour une géométrie donnée d'un moteur, il y a une relation immuable entre le volume V et l'angle 0. Calculer en permanence ce volume V d'après 0 n'est pas 25 envisageable dans un système de contrôle moteur embarqué car l'expression de V est trop complexe. dV A la place, pour estimer le volume V et sa dérivée d9 , il est possible d'utiliser une table de valeurs donnant le volume V pour chaque valeur de 0. 2905411 14 dT Pour estimer la température T et sa dérivée de , il existe plusieurs méthodes. Une première méthode préférée utilise la loi des gaz parfaits pour écrire : 5 p V = mr = constante. T Donc à deux instants ou deux angles différents : peVe po Vo TG To On en déduit : 10 To = To avec : To, po, Vo : température, pression et volume initiaux To , pe , Vo : température, pression et volume à l'angle e 15 La température est donc estimée pour chaque valeur de l'angle e , à partir de conditions initiales de pression, température et volume en début de cycle, de préférence juste après la fermeture de la soupape d'admission. On peut prendre la température initiale constante et égale à 300K, ou bien la faire dépendre de la température d'air admis, si on la mesure. 20 Cette méthode est préférée sur la suivante, car elle nécessite moins de calculs: T n'est calculé qu'une seule fois par cycle et il suffit de multiplier po Vo ce résultat par peVe à chaque nouvelle valeur de l'angle O. On expose ici une deuxième méthode pour estimer la température T dT 25 et sa dérivée de . Toujours en supposant le gaz parfait (pV = mrT, où mr est contant), on peut écrire : peVe po Vo 2905411 15 1dT 1 dp 1dV T dO p dO V dO On pourrait réintroduire cette expression dans l'expression de C, mais dp cela ferait réapparaître la dérivée de la pression d9 , qui est une source de bruit pour le calcul. Il est donc préférable de trouver un autre moyen pour l dT 5 exprimer T d8 . Pour cela, on discrétise l'expression, et on considère les valeurs prises par V, p, T aux angles 0 et ùA0 , où 3,0 est constant et représente le pas de calcul. L'expression devient : 1 To ù Te-oo = 1 Po ù Pe-AO + 1 Vo ù Ve-oe Te 09 po M VB M 10 On en tire une expression de la température T à l'angle en fonction de celle à l'angle 0ùA0 : Te = Te-oe Pe-oe + Ve-oe _1 Po Ve ) La température est donc estimée de manière récursive pour chaque valeur de l'angle e, à partir d'une température initiale correspondant à la 15 température dans le cylindre à l'instant de la fermeture de la soupape d'admission. On peut la prendre constante et égale à 300K, ou bien la faire dépendre de la température d'air admis, si on la mesure. Pour terminer, l'expression de C est également discrétisée, et la valeur de C pour chaque valeur de 0 est donnée par la formule : 1 po Te ùTe-oe + Po vo -1-é-oe yù1\TB A0 ) Vo Il est possible d'aller encore plus loin dans la simplification du calcul de C, en optant pour une simplification de la détermination du volume Ve . 1 Co = 20 2905411 16 En effet, comme il est dit plus haut, la loi de volume V, qui est immuable pour une géométrie de moteur donnée, peut être déterminée par une table de valeur, associant à chaque valeur de 9 une valeur de volume. Dans un système embarqué, il est cependant intéressant de diminuer 5 la taille des tables de valeurs car ceci consomme des ressources de mémoire. La simplification exposée ici va permettre de s'affranchir de tables de valeurs en mémoire. Plutôt que d'utiliser le volume vrai dans les calculs exposés précédemment, on lui substituera une loi parabolique Vp donnée par 10 l'expression. Vp =Vo +Kp62 où : Vo est constant et est le volume au point mort haut (volume mort), Kp un coefficient à ajuster de telle façon que la courbe de Vp soit au plus près de la courbe de V au voisinage du PMH. Un tel exemple de courbe est représenté à la 6. 15 Vp sera partout substitué à v dans les calculs exposés précédemment : pour le calcul des températures Te et T0-ee , et dans le PO VO Ve-oe terme VB 09 Bien que Vp soit très éloigné de V en dehors de la zone du PMH, on peut démontrer que l'influence sur le résultat est minime. En pratique 20 l'erreur est très peu perceptible, comme on le voit sur les figures 7 et 8. La figure 9 montre une synthèse des principales étapes d'un procédé selon l'invention. A partir de l'angle du vilebrequin 8 donné par le détecteur 9, le boîtier 7 peut calculer le volume de la chambre 3, ou de son substitut parabolique, 25 à au moins deux valeurs de l'angle du vilebrequin (ou à deux instants). La valeur du volume, et la valeur de la pression dans la chambre de combustion donnée par le capteur 4 permettent au boîtier 7 de mesure la température dans la chambre de combustion, à au moins deux valeurs de l'angle du vilebrequin (ou à deux instants). 2905411 17 Le volume, la pression et la température permettent au boîtier 7 de calculer la valeur de C, à au moins deux valeurs de l'angle du vilebrequin (ou à deux instants). La figure 10 représente schématiquement un montage logique 5 permettant de mettre en oeuvre un procédé selon l'invention. Le boîtier 7 comporte ainsi un calculateur 71 calculant au moins deux valeurs Co et Ce_oe de la variable C caractéristique. La deuxième valeur Co de la variable caractéristique est envoyée à un comparateur 72 qui détermine si la variable caractéristique dépasse un 10 premier seuil SI. Un différenciateur 74 calcule la différence entre la deuxième valeur Co de la variable caractéristique et la première valeur CB_48 de la variable caractéristique. La différence est envoyée à un comparateur 73 qui détermine si la 15 différence est inférieure à un deuxième seuil S2. La porte logique ET (AND en anglais) 75 détermine si les deux conditions des comparateurs 72 et 73 sont vérifiées. Si c'est le cas, le début de la combustion est déterminé. Le différenciateur 74 calcule de préférence la valeur absolue de la 20 différence. Les deux seuils ont pour effet d'éviter les fausses détections dues au bruit de mesure. Ce bruit a pour caractéristique une forte variation en amplitude d'une valeur de Co à sa suivante. Les seuils SI et S2 peuvent être fixes ou issus d'une cartographie en 25 fonction de la charge moteur (ou quantité de carburant injectée) par exemple. Lors de fortes injections, on sait que l'amplitude de C va être importante, il est alors judicieux d'augmenter S1 pour diminuer le risque de fausses détections dû au bruit de mesure. Dans ce cas là, il est aussi 30 nécessaire d'augmenter S2 car l'amplitude de variation sur un AO est plus importante. 2905411 18 La figure 11 montre que dans le cas d'un moteur à au moins deux injections de quantités injectées différentes dans la chambre (par exemple une injection pilote de plus faible quantité et une injection principale de plus forte quantité), un procédé selon l'invention comporte une étape consistant 5 à détecterséparément le début de combustion de chacune des injections. Une première possibilité est alors d'effectuer la détection de chaque début de combustion en multipliant les détections et en appliquant des seuils différents chaque détection. Une autre possibilité est d'effectuer la détection de chaque début de 10 combustion en multipliant les détections, et en déterminant si la différence entre la deuxième valeur Ca de la variable caractéristique et la première valeur Ce_4B de la variable caractéristique pour l'injection de plus grande quantité est supérieure à un deuxième seuil S2. Ceci a pour avantage d'éviter que la détection de la combustion principale ne détecte la 15 combustion de l'injection pilote. Dans ce cas, le deuxième seuil S2 est identique pour chaque détection. Une autre possibilité est que la détection de chaque début de combustion s'effectue en appliquant deux seuils pour l'injection de plus 20 faible quantité, et un seuil unique à l'injection de plus grande quantité. • REVENDICATIONS 1. Procédé de détermination d'un début d'une combustion dans une chambre (3) à combustion d'au moins un cylindre d'un moteur à combustion interne, le cylindre comportant un piston (5) mobile entre un point mort haut (PMH) et un point mort bas (PMB), le piston (5) étant en outre relié à un vilebrequin (8) monté à rotation dans le moteur, caractérisé en ce qu'il comporte les étapes consistant, grâce à un boîtier électronique (7), à : a) calculer au moins deux valeurs d'une variable (Ce , Ce-oe ) caractéristique de la combustion dans la chambre (3) de combustion, la variable caractéristique étant notamment une fonction de l'énergie dégagée dans la chambre (3), une première valeur (CB_49) étant calculée pour un premier angle (9-A8) de rotation du vilebrequin (8) ou un premier instant (t-15 At), une deuxième valeur (Co) étant calculée pour un deuxième angle de rotation (0) du vilebrequin (8) ou un deuxième instant (t); b) déterminer si la deuxième valeur (Co) de la variable caractéristique dépasse un premier seuil (SI) ; et c) déterminer si la différence entre la deuxième valeur (C9) de la 20 variable caractéristique et la première valeur (CB_49) de la variable caractéristique est inférieure à un deuxième seuil (S2), le début de la combustion étant déterminé si les conditions des étapes b) et c) sont vérifiées. 25 2. Procédé selon la revendication 1, dans lequel le deuxième angle (0) ou le deuxième instant (t) est proche du premier angle (O-A8) ou du premier instant (t-At). 3. Procédé selon l'une des revendications 1 ou 2, dans lequel le calcul 30 de la variable caractéristique C est effectué par : 2905411 20 1 dQ F dO où Q est l'énergie dans la chambre, 0 est l'angle du vilebrequin, et F est une fonction de e telle que F est maximale au point mort bas, 5 minimale au point mort haut, monotone entre les deux, et toujours strictement positive. 4. Procédé selon la revendication 3, dans lequel la fonction F est donnée par: 10 F(9)=Fo+K.9l où K et Fo sont des constantes. 5. Procédé selon l'une des revendications 1 ou 2, dans lequel le calcul de la variable caractéristique C est effectué par : C=1dQ Vd9 où V représente le volume de la chambre pour une valeur 0 de l'angle, Q est l'énergie dans la chambre, et 0 est l'angle du vilebrequin. 6. Procédé selon la revendication 5, dans lequel on a : dQ - 1 ( p dm p dT p dV y dV de V ù+ù- + yù1 md9 Td9 Vd9 yù1p d9 Avec : Q la quantité d'énergie présente dans les gaz de la chambre de 25 combustion ; 0 l'angle du vilebrequin ; y le rapport des chaleurs spécifiques (et ci') ; V le volume de la chambre de combustion ; p la pression dans la chambre de combustion ; 2905411 21 m la masse de gaz enfermés dans la chambre; et T la température dans la chambre de combustion. 7. Procédé selon la revendication 5, dans lequel on a : dQù V pdT dV de yù1. T de ,+p dB Avec : Q la quantité d'énergie présente dans les gaz de la chambre de combustion ; 0 l'angle du vilebrequin ; 10 y le rapport des chaleurs spécifiques (cp et C ) ; V le volume de la chambre de combustion ; p la pression dans la chambre de combustion ; et T la température dans la chambre de combustion. 15 8. Procédé selon l'une des revendications 5 à 7, dans lequel Te =To po Vo avec To, po,Vo les température, pression et volume initiaux. 9. Procédé selon l'une des revendications 5 à 7, dans lequel i 1 Po-Ae + VB-4B -1 Po Vo 10. Procédé selon les revendications 7 et 9, dans lequel i l 1 Pe Te ùT ce 0-40 + PO VO ù Ve -oe yù1 \ Te 09 ) Vo A0 11. Procédé selon l'une des revendications 5 à 10, dans lequel le volume est tiré d'une table de valeurs donnant le volume V et sa dérivée temporelle pour chaque valeur d'angle du vilebrequin. 5 PO VO 20 Te ùTe Ao 2905411 22 12. Procédé selon l'une des revendications 5 à 11, dans lequel le calcul du volume vrai V est remplacé par une approximation Vp suivant une loi parabolique donnée par l'expression : 5 Vp(0) = Vo +KpO2 où : Vo est constant et est le volume de la chambre au point mort haut (PMH), Kp un coefficient à ajuster de telle façon que la courbe de Vp soit au plus près de la courbe de V au voisinage du point mort haut. 10 13. Procédé selon l'une des revendications 1 à 12, dans lequel le premier seuil (Si) et le deuxième seuil (S2) sont fixes ou issus d'une cartographie fonction de la charge du moteur ou de la quantité de carburant injectée. 14. Procédé selon l'une des revendications 1 à 13, comportant, dans le cas d'un moteur à au moins deux injections de quantités injectées différentes dans la chambre, une étape consistant à détecter séparément le début de combustion de chacune des injections. 15. Procédé selon la revendication 14, dans lequel la détection de chaque début de combustion s'effectue en multipliant les détections et en appliquant des seuils différents pour chaque détection. 25 16. Procédé selon la revendication 14, dans lequel la détection de chaque début de combustion s'effectue en multipliant les détections, et en déterminant si la différence entre la deuxième valeur (Co) de la variable caractéristique et la première valeur (Ce_oe) de la variable caractéristique pour l'injection de plus grande quantité est supérieure à un deuxième seuil 30 (S2). 15 20 2905411 23 17. Procédé selon la revendication 16, dans lequel le deuxième seuil (S2) est identique pour chaque détection. 18. Procédé selon la revendication 14, dans lequel la détection de 5 chaque début de combustion s'effectue en appliquant deux seuils pour l'injection de plus faible quantité, et un seuil unique à l'injection de plus grande quantité. 19. Véhicule comportant un indicateur d'un début d'une combustion 10 dans une chambre (3) à combustion, le véhicule comportant un moteur à combustion interne, le moteur comportant au moins un cylindre et un vilebrequin (8) monté à rotation dans le moteur, le cylindre comportant une chambre de combustion (3) et un piston 15 (5), le piston (5) étant relié au vilebrequin (8) et étant mobile entre un point mort haut (PMH) et un point mort bas (PMB), caractérisé en ce qu'il comporte un boîtier électronique (7) comportant: a) un calculateur (71) d'au moins deux valeurs d'une variable 20 (Co , Co-4e) caractéristique de la combustion dans la chambre (3) de combustion, la variable caractéristique étant notamment une fonction de l'énergie dégagée dans la chambre (3), une première valeur (CB_oe) étant calculée pour un premier angle (9-ie) de rotation du vilebrequin (8) ou un premier instant (t-nt), une deuxième valeur (Co) étant calculée pour un 25 deuxième angle de rotation (9) du vilebrequin (8) ou un deuxième instant (t); b) un comparateur (72) pour déterminer si la deuxième valeur (C9) de la variable caractéristique dépasse un premier seuil (SI) ; c) un différenciateur (74) pour calculer la différence entre la 30 deuxième valeur (Co) de la variable caractéristique et la première valeur (CB_oe ) de la variable caractéristique ; et • 2905411 5 24 d) un comparateur (73) pour déterminer si la différence est inférieure à un deuxième seuil (S2), le début de la combustion étant déterminé si les conditions déterminées par les comparateurs (72, 73) sont vérifiées. 20. Véhicule selon la revendication 19, dans lequel le moteur comporte en outre : - un détecteur (9) de l'angle du vilebrequin (8) ; et - un capteur (4) de la pression dans la chambre (3) de combustion. 10