Brevet FR2904660A1 — publié le 20080208

--> L'invention se situe dans le domaine des moteurs à combustion interne (Diesel ou Essence) à injection directe ou indirecte de carburant à N nombre de cylindres. L'invention concerne plus précisément la détermination d'un début de combustion d'un carburant injecté dans un moteur à combustion interne, sur un nombre n de cycles (n étant un entier supérieur ou égal à 1 ) de déplacement d'un piston dans un cylindre. L'invention concerne tout particulièrement l'estimation de l'indice de cétane du carburant injecté, trouvé à partir de la détermination du début de combustion. Les normes antipollution devenant de plus en plus sévères, il est important pour les constructeurs automobiles de bien maîtriser le pilotage de l'injection de carburant dans leurs moteurs pour limiter les émissions polluantes. L'un des paramètres importants sur ces systèmes d'injection est la nature du carburant. Notamment pour les moteurs Diesel ou Essence, l'indice de cétane du carburant est à prendre en compte. En effet, le cétane influence notamment le délai d'inflammation du carburant. Or la maîtrise de ce dernier délai est important car il conditionne le phasage de la combustion dans le cycle thermodynamique par rapport à une consigne d'avance à l'injection donnée, et une combustion mal phasée entraîne une élévation des émissions polluantes. La connaissance de l'indice de cétane du carburant s'avère donc très utile. Le document WO 2005/119034 divulgue une estimation de l'indice de cétane du carburant sur la base du délai d'auto-inflammation trouvé entre un temps de référence et le moment d'auto-inflammation. Le moment d'auto-inflammation est déterminé à partir de l'étude de évolution temporelle de la quantité d'énergie existant dans les gaz présents dans la chambre de combustion au cours d'un cycle, la quantité d'énergie étant -->déduite, à un moment déterminé du fonctionnement du cylindre, du volume de la chambre de combustion et de la pression mesurée dans le cylindre. S'il est détecté que ladite quantité d'énergie augmente soudainement d'une valeur supérieure à une valeur seuil, la combustion est alors considérée comme détectée, et le moment associé à cette augmentation d'énergie est estimé être ledit moment d'auto-inflammation. L'indice de cétane est alors calculé. Un objet de l'invention est d'augmenter le rapport signal sur bruit dans le calcul de ladite évolution de quantité d'énergie, et ainsi de fiabiliser encore plus la détection du moment d'inflammation. Un autre objet de l'invention est de diminuer les risques d'erreurs quant à la valeur du moment d'inflammation. Pour atteindre ces objets, l'invention propose, selon un premier aspect, un procédé de détermination d'un début de combustion d'un carburant lors de son injection dans un moteur à combustion interne, sur un nombre n de cycles (n étant un entier supérieur ou égal à 1 ) de déplacement d'un piston dans un cylindre, comprenant les étapes suivantes pour différents angles successifs que prend le vilebrequin du piston dans chaque cycle : (a) évaluer une variable C θ à un angle θ du vilebrequin dans la chambre de combustion du cylindre, C θ étant fonction d'une évolution angulaire de la quantité d'énergie dans les gaz présents dans la chambre de combustion à un angle θ ; (b) comparer chaque C θ à une valeur seuil S1 au-delà de laquelle la combustion pourrait avoir commencé ; (c) comparer | C θ - C Θ -A& | à une valeur S2 , C â- Δ θ étant une variable précédemment évaluée pour l'angle 0-ΔΘ , S2 étant une valeur limite au-delà de laquelle C θ est considérée comme une valeur non suffisamment stable pour identifier une combustion ; (d) s'il s'avère, à l'issue des deux précédentes étapes de comparaison, que le seuil S1 a été dépassé et que le seuil S2 n'a pas -->été dépassé pour un C θ , la valeur θ est enregistrée comme étant la valeur angulaire de début de combustion pour le cycle considéré. D'autres caractéristiques optionnelles de ce procédé de détermination d'un début de combustion sont : — ladite évolution angulaire de la quantité d'énergie à un angle θ correspond à la dérivée angulaire de la quantité d'énergie Q à un angle θ , qui est calculée à partir de l'expression suivante : dQ _ _ P_ dT_ dV_ dθ ~ y - 1 T dθ + P dθ p étant la pression dans le cylindre à l'angle θ , V étant le volume de la chambre de combustion à l'angle θ , T étant la température dans la chambre de combustion à l'angle θ , étant le rapport de chaleur spécifiques c p /c v . — C θ est calculée à partir de : une mesure de l'angle de vilebrequin θ , prises dans le cylindre en fonctionnement ; une mesure de la pression à l'angle de vilebrequin θ , prise dans le cylindre en fonctionnement ; une sélection du volume de la chambre de combustion à l'angle θ dans une liste connue de valeurs volumiques classées par angle vilebrequin ou un calcul du volume selon une loi mathématique ; un calcul de la température de la chambre de combustion à l'angle θ sur la base desdits volume et pression ; - S1 et ou S2 est une valeur fixe ou une valeur déterminée à partir d'une cartographie fonction de conditions d'injection du moteur ; - le moteur comprend une deuxième injection, de sorte qu'il y ait ainsi une injection pilote et une injection principale. Le procédé comprend en outre la mise en œuvre des étapes (a) et (b) pour la deuxième injection, et une étape finale d'enregistrement de la valeur θ comme étant la valeur angulaire de début de combustion pour le cycle considéré si le seuil S1 a été dépassé ; -->— alternativement, le moteur comprend une deuxième injection, de sorte qu'il y ait ainsi une injection pilote et une injection principale, le procédé comprenant alors la mise en œuvre des étapes (a) à (d) pour la deuxième injection ; Selon un deuxième aspect, l'invention propose un procédé d'estimation de l'indice de cétane d'un carburant injecté dans un moteur à combustion interne, sur un nombre n de cycles (n étant un entier supérieur ou égal à 1 ), comprenant les étapes suivantes : - déterminer un début de combustion d'un carburant injecté pour chacun des cycles, conformément audit procédé de détermination d'un début de combustion; - calculer l'angle moyen de début de combustion θ moy à partir des θ de début de combustion relevés pour l'ensemble des n cycles ; - calculer le délai de combustion entre un temps de référence déterminé et le temps début de combustion trouvé à partir de θ moy ; - déduire l'indice de cétane à partir du délai de combustion, sur la base d'un modèle déterminé ou d'une liste de correspondance entre des délais de combustion et des indices de cétane. D'autres caractéristiques optionnelles de ce procédé d'estimation d'indice de cétane sont: - l'indice de cétane est déduit à partir d'un ajustement de sa valeur afin de minimiser un écart entre ledit délai de combustion calculé et un délai de combustion estimé, le délai de combustion estimé étant calculé à partir de la loi d'Arrhenius prenant comme variables des pression, température, volume et masse de gaz mesurés dans le cylindre à chaque angle vilebrequin ; - alternativement, l'indice de cétane est déduit en faisant correspondre ledit délai de combustion calculé avec une expression empirique du délai de combustion, basée sur la loi d'Arrhenius, conduisant à l'expression suivante : -->B. r \ D C CN = - \ "pmh J Tpmh - où CN, T pmh et Ppm h sont respectivement l'indice de cétane, la température et la pression au point mort haut sans combustion, les coefficients A, B, C et D sont des constantes prédéterminées, et δ est ledit délai de combustion, avec : r = 7* rv?-' pnilt oclm p pmh = p adm rv? où rv=V adm /V P mh, où T adm et P a d m sont la température et la pression dans la chambre de combustion à la fermeture de la soupape d'admission, et la valeur de (γ ) est identifiée à chaque cycle pendant la phase de compression, c'est-à-dire entre la fermeture de la soupape d'admission et l'injection, de la manière suivante : = _ln(p -lnfo 7 W^Ϋ ) où les angles de fin ( fin ) et de début ( deb ) sont pris assez éloignés, soit respectivement juste avant l'injection et juste après la fermeture de la soupape d'admission. D'autres caractéristiques, buts et avantages de l'invention apparaîtront mieux à la lecture de la description détaillée qui suit, qui est non limitative, illustrée par les figures suivantes : La figure 1 représente schématiquement une vue en coupe d'un cylindre d'un moteur à combustion en fonctionnement. La figure 2 est un organigramme représentant différentes étapes d'un procédé de détermination d'un moment d'inflammation selon l'invention. La figure 3A représente une fonction triangle. La figure 3B représente un graphe comparatif entre une dérivée angulaire d'une quantité d'énergie dans le cylindre pour différents angles vilebrequin et la même dérivée multipliée par la fonction triangle de la figure 3A. Tpmh --> Les figures 4A et 4B sont deux graphes représentant chacun une dérivée angulaire d'une quantité d'énergie dans le cylindre pour différents angles vilebrequin, prise seule (figure 4A) ou multipliée par l'inverse du volume de la chambre de combustion (figure 4B). La figure 5A est un graphe représentant l'évolution du volume de la chambre de combustion au cours d'un cycle, en comparaison avec un modèle parabolique. Les figures 5B et 5C sont deux graphes représentant chacun une dérivée angulaire d'une quantité d'énergie dans le cylindre pour différents angles vilebrequin, multipliée par l'inverse du volume réel de la chambre de combustion (figure 5A) ou multipliée par l'inverse du volume selon un modèle parabolique de la chambre de combustion (figure 5B). La figure 6 représente différentes étapes menant au calcul de Ce. La figure 7 est un schéma fonctionnel d'un détecteur de début de combustion selon l'invention. La figure 8 est un schéma représentant deux détecteurs de début de combustion selon l'invention, mis en œuvre respectivement dans un moteur à double injection. La figure 9A est un graphe donnant des résultats d'estimation de moments de début de combustion en fonction du nombre de cycles piston. La figure 9B est un graphe donnant des résultats d'estimation délais d'auto-inflammataion en fonction de différentes avances à l'injection. En référence à la figure 1 , un dispositif selon l'invention y est illustré. Il s'applique à un moteur à combustion interne Diesel ou Essence comportant un piston mobile 5, une bielle 6, une chambre de combustion 3, un dispositif d'admission de gaz 2, un dispositif d'échappement de gaz 1 , un moyen de mesure apte à mesurer différents angles θ de rotation du vilebrequin 8, un capteur de pression 4 placé dans la chambre de combustion 3, et des moyens de traitement 7 aptes à recevoir les signaux provenant dudit capteur de pression 4 et du moyen de mesure 8, à synchroniser ces signaux pour pouvoir -->les corréler, et à traiter ces signaux pour mettre en œuvre un procédé selon l'invention. Lors du fonctionnement du moteur, une estimation de l'indice de cétane selon l'invention est faite. Elle s'effectue avantageusement dans des conditions de fonctionnement fixées à l'avance lors de la mise au point du moteur. L'indice de cétane du carburant est estimé, comme on le verra plus loin, grâce à l'estimation du délai d'inflammation ou d'auto-inflammation du carburant. Le délai d'inflammation ou d'auto inflammation est le temps qui s'écoule entre un temps de référence, e.g. le moment où l'on injecte le carburant, et l'instant où la combustion démarre. Il est donc avant tout essentiel de déterminer le moment du début de la combustion (encore noté SOC ). En référence à la figure 2, sont représentées différentes étapes permettant de réaliser cette détermination du SOC, basée sur le calcul et l'analyse d'une variable C, image de la combustion. A partir de la valeur de pression p θ régnant dans la chambre de combustion, et de la valeur de la position angulaire θ du vilebrequin, qui ont été mesurées lors de l'étape 10, une variable C β est en permanence calculée lors de l'étape 20. C θ est fonction d'une évolution angulaire de la quantité d'énergie Q dans les gaz présents dans la chambre de combustion 3 à un angle θ. Cette évolution angulaire de la quantité d'énergie est avantageusement la . , . . dQ dérivée -^ . dθ Cette valeur C β est calculée en permanence sur la base du dégagement d'énergie dans le cylindre, pour un angle θ évoluant entre -180° et +180° (i.e. dans un cycle), ou dans un intervalle plus réduit. Le principe est d'identifier un changement significatif dans la valeur de C θ pour diagnostiquer le début de la combustion. La valeur de l'angle vilebrequin -->est alors enregistrée, et l'instant de début de combustion est ainsi déduit. Tant que ce changement significatif dans la valeur de C 9 n'a pas eu lieu, la combustion est estimée comme n'ayant pas commencé. Selon un premier mode de réalisation, C θ se calcule à partir de : C = -β (1 ) F dθ où Q est l'énergie dans le cylindre, θ est l'angle vilebrequin, et F est une fonction de θ telle que F est maximale au point mort bas (encore noté PMB), minimale au point mort haut (encore noté PMH), monotone entre les deux, et toujours strictement positive. En référence à la figure 3A, une fonction "triangle" possède de telles propriétés, et peut être définie de la façon suivante : F = F 0 + K\Θ\ (2) où K et F 0 sont des constantes. En pratique, le calcul du dégagement d'énergie — donne une dθ 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 PMB), par rapport aux zones où il est susceptible de contenir l'information qui intéresse (du coté du PMH), comme on le voit sur la figure 4B. Cette opération augmente donc le rapport signal sur bruit et diminue les contraintes sur le filtrage de la pression à mettre en place. Selon un deuxième mode de réalisation, C θ se calcule à partir de l'expression simplifiée du dégagement d'énergie volumique dans le cylindre suivante : C = ^ (3) V dθ où V représente le volume de la chambre de combustion 3 pour un angle θ. -->Comme le volume varie cycliquement avec θ, la multiplication par — a pour avantage de diminuer l'amplitude du bruit sur le dégagement d'énergie là où il est typiquement le plus fort (vers les PMB). Ce calcul amplifie aussi le signal là ou il est utile pour la détection de SOC (typiquement rencontré au voisinage du PMH). Cette opération augmente donc le rapport signal sur bruit et diminue les contraintes sur le filtrage de la pression à mettre en place (voir figure 2). En pratique, le calcul de — peut s'avérer complexe et peut donc nécessiter de la 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 réaliser les simplifications exposées ci-après. On part de l'expression, bien connue de l'homme du métier, du dégagement d'énergie apparent — du gaz présent dans la chambre : dθ dθ γ - l [m dθ T dθ V dθ J γ - V dθ avec : - Q : quantité d'énergie présente dans les gaz de la chambre de combustion 3 ; - θ : angle vilebrequin ; - γ : rapport des chaleurs spécifiques (c p et c v ) à θ ; - V : volume de la chambre de combustion à θ ; - p : pression dans la chambre de combustion à θ ; - m : masse de gaz enfermés dans le cylindre à θ ; - T : température dans la chambre de combustion à θ. Il est à noter que l'expression (4) pourrait être écrite en remplaçant θ par le temps t. Dans le cadre de l'invention, il sera toutefois préféré l'utilisation de l'angle θ, car le système devient alors indépendant de la vitesse de rotation du moteur. -->Il est à noter par ailleurs que l'expression (4) prise comme point de départ contient la dérivée angulaire de la température et non pas la dérivée de la pression. Ceci va permettre de minimiser le bruit, car la dérivée de la pression p mesurée introduit un bruit important dans le résultat. Avantageusement selon l'invention, seules la pression p et l'angle Θ sont des valeurs mesurées par le dispositif selon la figure 1. Les autres variables de l'expression (4) sont négligées ou estimées, comme nous allons le voir. Concernant la masse m, l'hypothèse est prise de la considérer comme constante (la masse injectée est donc négligée). L'expression de C devient alors : 1 dQ 1 ( p dTλ p dV V dθ γ -W T dθ V dθ Le rapport des chaleurs spécifiques γ est considéré constant avant la combustion ainsi que dans les premiers instants de la combustion. On pourra par exemple l'estimer égal à 1 ,4. Pour une géométrie donnée d'un moteur, il y a une relation immuable dV entre le volume V et l'angle θ. Pour estimer le volume V et sa dérivée — , une dθ table de valeurs donnant le volume V pour différentes valeurs de θ est utilisée. La température T et sa dérivée — sont calculées sur la base desdits dθ volume V et pression p. Selon un premier mode de calcul de T, on suppose que le gaz est parfait. On peut alors écrire à tout moment : pV - mr = constante (6) Donc à deux instants ou à deux angles différents : C = T dθ + ^- V dθ (5) --> n P avec : To, po. V 0 : température, pression et volume initiaux (i.e. à l'injection) Te, Pβ, θ : température, pression et volume à l'angle θ La température est donc estimée pour chaque valeur de l'angle θ, à 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 pourra 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. Cette méthode peut être préférée car elle nécessite peu de calculs, étant T donné que — — n'est calculé qu'une seule fois par cycle et qu'il suffit de P V o multiplier ce résultat par p β \f θ à chaque nouvelle valeur de l'angle θ. Selon un deuxième mode de calcul de T, on suppose aussi que le gaz est parfait (pV = mrT, où mr est contant), mais on en déduit cette fois : \ dT 1 dp 1 dV , n . = — — + (9) T dθ p dθ V dθ Cette expression pourrait être réintroduite dans l'expression de C, mais cela ferait réapparaître la dérivée de la pression — , qui est une source de bruit pour le calcul. Un autre moyen pour exprimer sera donc préféré. Pour cela, l'expression (9) est discrétisée, et les valeurs prises par V, p, T sont considérées aux angles θ et 0-ΔΘ, où Δθ est constant et représente le pas du calcul. L'expression (9) devient alors : Pβ'e _ Λ o r\ T θ T 0 On en déduit : T θ =T 0 ^ i- (8) --> 1 T a -AΘ 1 PΘ-PΘ -Aθ + - 1 V 0 - Θ-AΘ T θ Aθ p θ Aθ V θ Aθ L'expression suivante de la température T à l'angle θ est alors déduite, en fonction de la température à l'angle 0-ΔΘ : (10) τ a =τ 1 Θ-AΘ Pθ-AΘ + v , (11) -1 Θ ,-AΘ ' PΘ ' Θ La température est donc estimée de manière récursive pour chaque valeur de l'angle θ, à partir d'une température initiale correspondant à la température dans le cylindre à l'instant de la fermeture de la soupape d'admission. Elle peut être prise constante et égale à 300K, ou bien dépendante d'une température d'air admis mesurée par des moyens appropriés. Pour terminer, l'expression de C est également discrétisée, et la valeur de C pour chaque valeur de θ est donnée par la formule : C=- 1 Pθ T θ - T 6 6 θ-AΘ γ-\ Aθ + PeV t t θ-AΘ V a Aθ (12) En outre, l'expression (9), (11 ) ou (12) peut être encore plus simplifiée, en termes de calcul, en simplifiant l'expression de V (ou W e ). 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 valeurs, associant à chaque valeur de θ une valeur de volume. Dans un système embarqué il est cependant intéressant de diminuer la taille des tables de valeurs car elles consomment des ressources de mémoire. La simplification exposée ci-après va permettre pour le moins de diminuer la taille de tables de valeurs en mémoire, et pour le plus de s'en affranchir. Plutôt que d'utiliser le volume vrai dans les calculs exposés précédemment, on lui substituera une loi parabolique V p donnée par l'expression : V p = V 0 + K p θ 2 , où V 0 est constant et est le volume au PMH (volume mort), -->K p un coefficient à ajuster de telle façon que la courbe de V p soit très proche à la courbe de V au voisinage du PMH (voir figure 5A pour illustration). L'influence de cette simplification sur C est illustrée en comparant les figures 5B et 5C, dans lesquelles on voit l'effet négligeable visible de l'approximation parabolique. V p peut ainsi être partout substitué à V dans les calculs exposés précédemment : pour le calcul des températures T θ et T θ-Δθ , et dans le terme Pθ * θ ~ "θ-AΘ V θ Aθ Du temps de calcul est ainsi gagné et des ressources sont économisées. En synthèse, compte tenu des simplifications énoncées précédemment, la chaîne de calcul menant à la variable C β peut se schématiser par la figure 6. En référence à la figure 2, une fois que l'étape de calcul 20 de C e a été réalisée, une détection du SOC est mise en œuvre, au moyen de deux étapes 30 et 40 de comparaison de C θ avec des valeurs seuil respectives S1 et S2 . Deux tests 30-40 sont ainsi successivement effectués. Le premier test 30 observe si la variable C θ dépasse une valeur S1 audelà de laquelle la combustion pourrait avoir commencé. En effet, la pente prise par Q en θ est alors suffisamment grande pour estimer que, a priori, la combustion est détectée. Le deuxième test 40 compare C θ et Cρ. Δθ en observant si la valeur absolue de la différence entre C θ et C Ô- ΔΘ est inférieure à une valeur seuil S2, S2 étant une valeur limite au-delà de laquelle on considère que la combustion ne peut démarrer. Si ces deux conditions sont vérifiées, le SOC est détecté. Le deuxième test a notamment pour objet d'éviter les fausses détections dues au bruit de mesure, ce bruit ayant pour caractéristique une forte variation en amplitude d'une valeur de C θ à sa suivante. Le deuxième test permet donc d'éviter les dépassements accidentels de S1. En outre, étant donné que ce -->deuxième test offre une garantie supplémentaire contre des erreurs de détection du SOC, il est permis d'ajuster plus finement le premier seuil S1 pour améliorer la détection. Les seuils S1 et S2 peuvent être fixes ou issus d'une cartographie en fonction par exemple de paramètres d'injection du moteur, tels que la charge moteur (ou la quantité de carburant injectée). Ces valeurs seuil peuvent être basées par exemple sur des expérimentations préalables. Ainsi, sachant que, lors de fortes injections, l'amplitude de C est importante, il est judicieux d'augmenter S1 pour diminuer le risque de fausses détections dues au bruit de mesure. Dans ce cas là, il est aussi utile d'augmenter S2 car l'amplitude de variation sur un Δθ est plus importante. Le dispositif logique de détection du SOC est illustré sur la figure 7, recevant en entrée 29 la variable C β, et comprenant un comparateur 31 permettant la comparaison de C θ avec S1 , un moyen 41 agencés pour lire la variable C Θ- Δ Θ préalablement calculée, des éléments de calculs 42-43 aptes à réaliser l'opération de soustraction en valeur absolue entre C θ et C Θ -Δ Θ , un comparateur 44 réalisant la comparaison de \C Θ - C Θ - Δ Θ| avec S2, et l'élément logique AND 35 fournissant 1 en sortie 51 si les deux conditions sont remplies (i.e. SOC détecté) et 0 en sortie 51 si zéro ou une des deux conditions ne sont pas remplies (i.e. SOC non détecté). Selon une configuration particulière de l'invention, le moteur est à double injection (une injection pilote et une injection principale), dont la première (la pilote) est différente en quantité injectée de la seconde (principale). Il est à remarquer qu'on peut appliquer le principe précédemment expliqué au cas multi-injection avec plus de deux injections, en appliquant à chaque injection (ex. : pilote, principale et post) ledit algorithme. Pour une double injection, on peut vouloir détecter séparément le SOC de chacune des deux injections. Pour cela il suffit de dupliquer ladite détection (selon la figure 7), et appliquer des seuils différents aux deux détections (voir figure 8 pour illustration). -->Optionnellement, on peut modifier le comparateur 31 pour l'injection principale en lui faisant appliquer un =. Ceci a pour avantage d'éviter que l'algorithme de détection de la combustion principale puisse détecter la combustion de l'injection pilote. Dans ce cas, le seuil S2 de l'injection principale peut alors être fixé égal au seuil S2 de la pilote. En variante, on peut aussi appliquer deux seuils à l'injection pilote et un seuil à l'injection principale (le SOC de cette dernière étant plus simple à détecter car plus facile à discriminer du bruit de par son amplitude plus élevée). En référence à la figure 2, une fois la détection (40) du SOC effectuée, l'angle vilebrequin θ correspondant est enregistré en (50). Pour un nombre n de cycles de vilebrequin (n étant un entier supérieur ou égal à 1 ), n angles de SOC peuvent ainsi être mémorisés. En référence à la figure 10, l'indice de cétane est alors estimé. En premier lieu, une étape de moyennage (60) des valeurs de ces n angles mémorisées est mise ne œuvre, donnant une valeur 0 moy . En second lieu, le délai δ d'inflammation ou d'auto-inflammation est ensuite calculé (70). Ce délai δ est pris égal à la différence dudit angle moyen de SOC avec une valeur de référence, qui peut être la valeur de l'avance à l'injection (encore notée SOI), l'avance à l'injection étant la durée en angle (ou en temps) entre le point mort haut (PMH) et l'angle du début de l'injection. Grâce à la connaissance du début de combustion (SOC) et du phasage de l'injection (SOI) on est alors capable de reconstruire le délai d'auto inflammation δ. δ = SOC - SOI (13) Les figures 9A et 9B montrent les résultats d'estimation de SOC et de durées d'auto-inflammation δ pour un moteur Diesel déterminé. Le SOC est donné en fonction du nombre de cycles (figure 9A), et le délai d'autoinflammation δ en fonction du de l'avance à l'injection (figure 9B). Le SOC et le délai δ ont été estimés pour différentes valeurs d'avances à l'injection (SOI), -->exprimées en degré vilebrequin, et variant entre -15 et +15, avec un pas de 5. La durée d'auto-inflammation est quant à elle exprimée à la fois en degré vilebrequin et en millisecondes. Le PMH correspondant à 360°, on peut constater en premier lieu que, comme dit précédemment, l'auto-inflammation est déclenchée à proximité du PMH. On peut aussi constater que le délai d'injection évolue entre 1 ,18 ms et 2,15 ms selon la valeur de l'avance à l'injection, les valeurs les plus faibles étant ici trouvées pour les faibles valeurs de l'avance à l'injection. En référence à la figure 10, l'indice de cétane est alors estimé (80), en incluant la valeur du délai d'auto-inflammation δ dans un modèle donnant une estimation de l'indice de cétane en fonction du délai d'auto-inflammation. Selon un premier mode de calcul de l'indice de cétane, celui-ci est estimé sur plusieurs cycles n en régime stabilisé, en ajustant sa valeur afin de minimiser l'écart entre ledit délai de combustion calculé δ et un délai de combustion estimé δe. Le délai de combustion estimé δe est calculé à partir de la loi d'Arrhenius suivante (comme on le verra plus loin) : Ea τ=P 'n .Aexp — (i 4) où 7 représente un temps caractéristique qui sera intégré au cours du cycle moteur afin de vérifier qu'il atteint un critère et donc que la combustion est déclenchée ; P et T sont respectivement la pression et la température des gaz dans la chambre de combustion à un angle θ ; les coefficients n et A sont des coefficients à calibrer ; les termes Ea et r représentent l'énergie d'activation du carburant et la constante des gaz parfaits. -->Il existe, dans la littérature une multitude de valeurs données aux paramètres n, A, Ea et r de la loi d'Arrhenius. Ces paramètres ont été typiquement déterminés à partir d'expériences de type tube à choc. Nous donnons ci-dessous des valeurs assez classiques utilisées pour du gazole conventionnel (Stringer, Clark et Clark - The spontaneous ignition of hydrocarbon fuels in a following System, Proc. Instn. Mech. Engrs vol 184 pt 3J, 1969) : n A Ea/r 0.757 0.0405 5473 Pour calculer T, les variables sont donc la pression, la température, le volume et la masse de gaz qui peuvent tous être mesurés dans le cylindre à chaque angle vilebrequin. Cependant, on peut faire l'économie des moyens de détection de la température (T), en retrouvant la valeur de cette dernière à partir de l'acquisition de la pression dans le cylindre (P) et de la connaissance du volume (V), grâce à la loi des gaz parfaits suivante. P.V T = (15) m.r où V, m et r sont respectivement le volume de la chambre de combustion 3, la masse de gaz enfermée dans la chambre 3 et la constante des gaz parfaits. Le volume (V) est par exemple déterminé à partir d'une cartographie ou calculé avec une loi classique de type : V cyl = f{θ), comme vu précédemment. La masse (m) enfermée dans le cylindre est déterminée par exemple à partir d'une cartographie de remplissage (fonction du régime et de la pression collecteur) ou du débit mesuré par le débitmètre auquel on ajoute l'EGR (acronyme anglo-saxon de Exhaust Gas Recirculation , pour gaz d'échappement recirculé à l'admission, afin de favoriser la diminution des polluants comme les NOx) présent dans le collecteur d'admission. Les valeurs de pression, température et volume sont enregistrés à chaque degré vilebrequin (ou sur un pas d'échantillonnage différent). -->Le délai d'auto-infammation r est donc déterminé à chaque instant (i.e. à chaque angle vilebrequin θ), au pas d'échantillonnage Δθ près. La combustion est estimée être déclenchée lorsque l'intégrale de l'inverse de T atteint un critère (ici 1 ). "? 1 ) ~.dt = l (16) t T où ti n i et t S0Ce sont respectivement l'instant de début d'injection et de début de combustion. Cette équation peut alors être transcrit en variable angulaire : „ γ i f -.—dθ (17) θ m, τ (O où ω est la vitesse angulaire du vilebrequin, et 0ini et θ S0Ce sont respectivement l'angle de phasage de l'injection et l'angle de début de combustion. Connaissant 0 ini , l'angle θ soc est ainsi obtenu à partir de r. Le délai d'auto inflammation estimé δ e est alors calculé: δ e = SOC e -SOI (18) où SOC e représente θ soc ou sa correspondance temporelle SOI est le temps de début d'injection, exprimée en angle vilebrequin ou en temps (ms). Une étape de comparaison entre le délai d'auto inflammation calculé δ et le délai d'auto inflammation estimé δ e est alors mise en œuvre. S'il existe une différence significative entre les deux délais, ledit paramètre d'énergie d'activation (Ea) de la loi d'Arrhenius est alors estimé erroné. Or le paramètre (Ea) est directement fonction du taux de cétane du carburant (CN) : Fh Ea = — (19) CN -->Sur les cycles suivants la valeur du taux de cétane (CN) est alors ajustée et l'opération du calcul de l'estimation du délai (δ e ) est à nouveau mise en œuvre afin de minimiser l'écart avec le délai d'auto inflammation calculé (δ). ε ≈ δ-δ. (20) Si l'écart (e) est positif le coefficient (CN) doit être augmenté, sinon diminué. Afin de déterminer l'évolution du (CN), une méthode de dichotomie ou autre de type gradient peut être utilisée. L'ajustement des paramètres des algorithmes est fait en fonction du degré de précision et de vitesse de convergence souhaité. Selon un deuxième mode de calcul de l'indice de cétane, celui-ci est estimé à partir de la formule empirique du délai d'auto-inflammation suivante, déduite à partir de la littérature de Hardenberg and Hase (SAE 790493) : δ = A. exp B D\ Ppmh J CN.T , p ,mh (21) où CN, T et P sont respectivement l'indice de cétane, la température et la pression au PMH sans combustion. Les coefficients A, B, C et D doivent être identifiés. Ils sont typiquement identifiés sur des essais au banc moteur avec des carburants dont l'indice de cétane est connu. De la dernière formule, une loi inverse permet d'obtenir l'indice de cétane pour un délai donné : B, CN-- C "pmh J (22) pmh .In A La pression et la température sont calculées à partir des lois de compression adiabatique suivantes : T pmh =T adm .rv-> { ) -->P p ,„„ = P π(tm -rv r (24) où rv = (Vadm pmh), où T a d m et P a d m sont la température et la pression dans la chambre de combustion à la fermeture de la soupape d'admission 2. La valeur de (γ ) est identifiée à chaque cycle pendant la phase de compression, c'est-à-dire entre la fermeture de la soupape d'admission 2 et la première injection. ^_mfej in(^) ^ )Mv deb ) où les angles (fin) et (deb) sont pris assez éloignés, soit respectivement juste avant la première injection (ou avant le PMH) et juste après la fermeture de la soupape d'admission 2. L'identification du γ permet de prendre en compte de meilleur façon les pertes aux parois. Le calcul de l'indice de cétane (CN) n'intervient qu'une fois par cycle après avoir identifié le coefficient (γ ) et que le délai de combustion δ a été calculé comme précédemment décrit. Cette approche est peu coûteuse en terme de calcul. Selon un troisième mode d'estimation de l'indice de cétane, une liste de correspondance entre des délais de combustion et des indices de cétane est établis à l'avance, sur la base par exemple d'expérimentations. Ainsi, connaissant la valeur du délai d'auto-inflammation δ, on peut lui faire correspondre un indice de cétane. Cet indice de cétane peut alors être pris en compte lors de la mise au point du moteur. Ceci est effectué pour un carburant de référence. L'invention permet en outre d'adapter l'injection selon des modifications trouvées dans les propriétés d'injection, telles que par exemple la nature du carburant (qui a nécessairement un indice de cétane différent), pour que le moteur fonctionne avec des réglages optimums. -->Les délais d'auto-inflammation qui dérivent de cet indice de cétane peuvent également être pris en compte (au même titre que les caractéristiques des injecteurs) dans la calibration de l'avance à l'injection. Un phasage de l'injection peut ainsi être calibré de sorte à assurer un bon phasage de la combustion par rapport au carburant de référence utilisé. --> REVENDICATIONS 1. Procédé de détermination d'un début de combustion d'un carburant lors de son injection dans un moteur à combustion interne, sur un nombre n de cycles (n étant un entier supérieur ou égal à 1 ) de déplacement d'un piston dans un cylindre, comprenant les étapes suivantes pour différents angles successifs que prend le vilebrequin du piston dans chaque cycle : (a) évaluer une variable C θ à un angle θ du vilebrequin dans la chambre de combustion du cylindre, C θ étant fonction d'une évolution angulaire de la quantité d'énergie dans les gaz présents dans la chambre de combustion à un angle θ ; (b) comparer chaque C θ à une valeur seuil S1 au-delà de laquelle la combustion pourrait avoir commencé ; (c) comparer | C θ - C Θ- ΔΘ | à une valeur S2 , C Θ- Δ Θ étant une variable précédemment évaluée pour l'angle 0-ΔΘ , S2 étant une valeur limite au-delà de laquelle C θ est considérée comme une valeur non suffisamment stable pour identifier une combustion ; (d) s'il s'avère, à l'issue des deux précédentes étapes de comparaison, que le seuil S1 a été dépassé et que le seuil S2 n'a pas été dépassé pour un C Θ , la valeur θ est enregistrée comme étant la valeur angulaire de début de combustion pour le cycle considéré. 2. Procédé selon la revendication précédente, caractérisée en ce que ladite évolution angulaire de la quantité d'énergie à un angle θ correspond à la dérivée angulaire de la quantité d'énergie Q à un angle θ , qui est calculée à partir de l'expression suivante : dQ _ _ d dV_ dθ ~ γ - 1 T dθ + P dθ p étant la pression dans le cylindre à l'angle θ , --> V étant le volume de la chambre de combustion à l'angle θ , T étant la température dans la chambre de combustion à l'angle θ , γ étant le rapport de chaleur spécifiques c p /c v . 3. Procédé selon l'une des revendications précédentes, caractérisé en ce que C θ est calculée à partir de : - une mesure de l'angle de vilebrequin θ , prises dans le cylindre en fonctionnement ; - une mesure de la pression à l'angle de vilebrequin θ , prise dans le cylindre en fonctionnement ; - une sélection du volume de la chambre de combustion à l'angle θ dans une liste connue de valeurs volumiques classées par angle vilebrequin, ou un calcul du volume selon une loi mathématique ; - un calcul de la température de la chambre de combustion à l'angle θ sur la base desdits volume et pression. 4. Procédé selon l'une des revendications précédentes, caractérisé en ce que S1 et/ou S2 est une valeur fixe ou une valeur déterminée à partir d'une cartographie fonction de conditions d'injection du moteur. 5. Procédé selon l'une des revendications précédentes, caractérisé en ce que le moteur comprend une deuxième injection, de sorte qu'il y ait ainsi une injection pilote et une injection principale, et en ce qu'il comprend en outre la mise en œuvre des étapes (a) et (b) pour la deuxième injection, et une étape finale d'enregistrement de la valeur θ comme étant la valeur angulaire de début de combustion pour le cycle considéré si le seuil S1 a été dépassé. 6. Procédé selon l'une des revendications 1 à 4, caractérisé en ce que le moteur comprend une deuxième injection, de sorte qu'il y ait ainsi une injection -->pilote et une injection principale, et en ce qu'il comprend en outre la mise en œuvre des étapes (a) à (d) pour la deuxième injection. 7. Procédé d'estimation de l'indice de cétane d'un carburant injecté dans un moteur à combustion interne, sur un nombre n de cycles (n étant un entier supérieur ou égal à 1 ), comprenant les étapes suivantes : - déterminer un début de combustion d'un carburant injecté pour chacun des cycles, conformément au procédé selon l'une des revendications précédentes ; - calculer l'angle moyen de début de combustion 0 moy à partir des θ de début de combustion relevés pour l'ensemble des n cycles ; - calculer le délai de combustion entre un temps de référence déterminé et le temps début de combustion trouvé à partir de θ moy ; - déduire l'indice de cétane à partir du délai de combustion, sur la base d'un modèle déterminé ou d'une liste de correspondance entre des délais de combustion et des indices de cétane. 8. Procédé d'estimation de l'indice de cétane selon la revendication précédente, caractérisé en ce que l'indice de cétane est déduit à partir d'un ajustement de sa valeur afin de minimiser un écart entre ledit délai de combustion calculé et un délai de combustion estimé, le délai de combustion estimé étant calculé à partir de la loi d'Arrhenius prenant comme variables des pression, température, volume et masse de gaz mesurés dans le cylindre à chaque angle vilebrequin. 9. Procédé d'estimation de l'indice de cétane selon la revendication 7, caractérisé en ce que l'indice de cétane est déduit en faisant correspondre ledit délai de combustion calculé avec une expression empirique du délai de combustion, basée sur la loi d'Arrhenius, conduisant à l'expression suivante : -->B. ' C Λ ° CN V P ' 1 ' J T pml,-to \A où CN, T pmh et Pp mh sont respectivement l'indice de cétane, la température et la pression au point mort haut sans combustion, les coefficients A, B, C et D sont des constantes prédéterminées, et δ est ledit délai de combustion, en ce que : pmh = T adm -rv 7-1 pmh P adm .rv' où rv = Vadm/Vpmh, où T ad m et P adm sont la température et la pression dans la chambre de combustion à la fermeture de la soupape d'admission, et en ce que la valeur de ( r ) est identifiée à chaque cycle pendant la phase de compression, c'est-à-dire entre la fermeture de la soupape d'admission et l'injection, de la manière suivante : r = - \n(v βll )-\n{V deb ) où les angles de fin ( fin ) et de début ( deb ) sont pris assez éloignés, soit respectivement juste avant l'injection et juste après la fermeture de la soupape d'admission.