Brevet FR2916243A3 — publié le 20081121

FR 2916243 A3 20081121 FR 0703493 A 20070515 La présente invention se situe dans le domaine des moteurs rotatifs à combustion interne (Diesel), à injection directe ou indirecte de carburant, et à N nombre de cylindres. Elle se rapporte plus précisément à un procédé de recalage de l'avance à l'injection en fonction de l'indice de cétane du carburant d'un moteur à combustion interne, sur un nombre n de cycles de déplacement d'un piston dans un cylindre. Elle est plus particulièrement applicable aux moteurs avec des combustions Diesel à fort taux de gaz d'échappement recirculés (EGR) et à faible 10 température de combustion. Les normes antipollution deviennent de plus en plus sévères, aussi bien aux Etats-Unis qu'en Europe. Il est important pour les constructeurs de véhicules automobiles de bien maîtriser le pilotage de l'injection de carburant dans les moteurs pour limiter 15 les émissions polluantes. L'un des paramètres important sur ces systèmes d'injection est la nature du carburant. Notamment pour les moteurs Diesel, l'indice de cétane du gasoil est à prendre en compte. Cet indice est différent aux Etats-Unis et en Europe. 20 Il conditionne les propriétés d'auto-inflammation du carburant (autrement dit le délai d'auto-inflammation), ainsi que le phasage de la combustion dans le cycle thermodynamique, par rapport à une consigne d'avance à l'injection donnée. Une combustion mal phasée entraînera une élévation des émissions polluantes. Il est donc très important de prendre en compte l'indice de cétane du 25 carburant. Actuellement, cet indice est pris en compte lors de la mise au point du moteur, en utilisant un carburant de référence. Les délais d'auto-inflammation qui en dérivent sont pris en compte (au même titre que les caractéristiques des injecteurs) dans la calibration de l'avance à injection. Ainsi, on calibre un phasage 30 de l'injection qui assurera un bon phasage de la combustion par rapport au carburant de référence utilisé. Mais si le véhicule est utilisé avec un autre carburant (avec un indice de cétane différent), le moteur ne fonctionne plus avec ses réglages optimums. 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 l'é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 de permettre le recalage à l'avance à l'injection de manière à ce que le moteur fonctionne selon des réglages optimums, quel que soit l'indice de cétane du carburant utilisé. L'invention propose donc un procédé de recalage à l'avance à l'injection en fonction de l'indice de cétane du carburant d'un moteur à combustion interne, sur un nombre n de cycles de déplacement d'un piston dans un cylindre. Ce procédé est remarquable par le fait qu'il comprend les étapes suivantes pour différents angles successifs que prend le vilebrequin du piston dans chaque cycle : a) calculer, sur banc moteur, lors de la mise au point dudit moteur, les valeurs de délai d'auto-inflammation de z carburants de référence présentant chacun un indice de cétane "I1, I2, .,IX,..., I," , x étant compris entre 1 et z ; b) en déduire les valeurs d'avance à l'injection correspondantes "SOI con" , c) calculer, pour un carburant en train d'être consommé et présentant un indice de cétane I,t, les valeurs de délai d'auto-inflammation correspondantes ; d) en déduire les valeurs d'avance à l'injection correspondantes "SOI base" e) affecter à chaque "SOI base", la "SOI co," correspondante pour un même indice de cétane, et l'enregistrer comme nouvelle valeur d'avance à injection. Selon d'autres caractéristiques avantageuses et non limitatives de ce 30 35 procédé : - on calcule les valeurs de délai d'auto-inflammation à partir d'une valeur de début de combustion, qui est elle-même déterminée à partir de l'évaluation d'une variable " Co" à un angle "0" du vilebrequin dans la chambre de combustion du cylindre, "Co" correspondant à la dérivée angulaire de la quantité d'énergie "Q" à un angle "0", qui est calculée à partir de l'expression suivante : dQ V p dT dV d8 y-1 T d9 + p dB "p" étant la pression dans le cylindre à l'angle "0", "V" étant le volume de la chambre de combustion à l'ange "0", "T" étant la température dans la chambre de combustion à l'angle "0", 10 "y" étant le rapport de chaleur spécifiques " cp/c,, ". Avantageusement, ce procédé consiste à : comparer chaque "Co" à une valeur seuil "Si" au-delà de laquelle la combustion pourrait avoir commencé ; - comparer "1 Co - Co_oo 1" à une valeur "S2", "Co_oo" étant une 15 variable précédemment évaluée pour l'angle "O-00", "S2" étant une valeur limite au-delà de laquelle "Co" est considérée comme une valeur non suffisamment stable pour identifier une combustion ; - s'il s'avère, à l'issue des deux précédentes étapes de comparaison, que le seuil "Si" a été dépassé et que le seuil "S2" n'a pas été dépassé pour un Co, la 20 valeur "0" est enregistrée comme étant la valeur angulaire de début de combustion pour le cycle considéré. Selon une variante de réalisation, ce procédé consiste à calculer la çdQ Jdû amplitude maximale AM, cette valeur CA10 étant enregistrée comme étant la valeur 25 angulaire de début de combustion pour le cycle considéré. Préférentiellement, ladite amplitude maximale AM est mesurée en établissant la différence entre : - les valeurs maximale et minimale de Q : ou 30 - la valeur maximale de Q et zéro ; ou - zéro ou la valeur minimale de Q, et l'angle de fin d'intégration de Q si celui-ci correspond à la valeur maximale de Q. valeur d'angle CA10 qui est égale à l'angle où Q = atteint 10% de son D'autres caractéristiques apparaîtront à la lecture de la description détaillée qui va suivre d'un mode de réalisation préférentiel. Cette description sera faite en référence aux dessins annexés dans lesquels : - la figure 1 représente schématiquement une vue en coupe d'un 5 cylindre d'un moteur à combustion interne en fonctionnement ; - la figure 2 est une représentation graphique, pour un carburant présentant un indice de cétane donné, des correspondances entre les valeurs de base d'avance à l'injection et les valeurs corrigées. En référence à la figure 1, un dispositif auquel s'applique l'invention 10 est illustré. Il s'agit d'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 "0" de rotation du vilebrequin 8, un capteur de pression 4 placé dans la chambre de combustion 3, et des moyens de traitement 7 15 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 oeuvre un procédé selon l'invention. Lors du fonctionnement du moteur, une estimation de l'indice de cétane est faite. Elle s'effectue avantageusement dans des conditions de 20 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 25 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"). Une détermination possible du SOC est basée sur le calcul et 30 l'analyse d'une variable C, image de la combustion. A partir de la valeur de pression po régnant dans la chambre de combustion, et de la valeur de la position angulaire 0 du vilebrequin, qui ont été mesurées lors d'une première étape, une variable Co est en permanence calculée lors d'une étape ultérieure. 35 Co 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 O...DTD: Cette évolution angulaire de la quantité d'énergie est avantageusement la dérivée d8 Cette valeur Co est calculée en permanence sur la base du dégagement d'énergie dans le cylindre, pour un angle 0 évoluant entre -180 et 5 +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 Co 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 Co n'a pas eu lieu, la combustion est 10 estimée comme n'ayant pas commencé. Selon une première possibilité, Co se calcule à partir de : C = F d9 (1) où "Q" est l'énergie dans le cylindre, "0" est l'angle vilebrequin, et 15 "F" est une fonction de "0" 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. Une fonction "triangle" possède de telles propriétés, et peut être définie de la façon suivante : 20 F = Fo + K.8 (2) où "K" et "Fe" sont des constantes. En pratique, le calcul du dégagement d'énergie ~~ donne une information bruitée. L'effet d'une telle fonction "F" est d'atténuer le signal de dégagement 25 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). Cette opération augmente donc le rapport signal sur bruit et diminue les contraintes sur le filtrage de la pression à mettre en place. Selon une deuxième possibilité, Co se calcule à partir de l'expression 30 simplifiée du dégagement d'énergie volumique dans le cylindre suivante : C=V dB (3) où V représente le volume de la chambre de combustion 3 pour un angle O. Comme le volume varie cycliquement avec 0, la multiplication par V 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. En pratique, le calcul de d e 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 : dQ 1 pdm pdT pdV y dV V + )+ p (4) dO yù1 \m dB T dB V dB yù1 dO avec : - Q : quantité d'énergie présente dans les gaz de la chambre de combustion 3 ; - 0 : angle vilebrequin ; - y : rapport des chaleurs spécifiques (cp et c,,) à 0 ; - V : volume de la chambre de combustion à 0 ; - p : pression dans la chambre de combustion à 0 ; - m : masse de gaz enfermés dans le cylindre à 0 ; - T : température dans la chambre de combustion à 0. Il est à noter que l'expression (4) pourrait être écrite en remplaçant 0 25 par le temps t. Il est toutefois préféré l'utilisation de l'angle 0, 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 30 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, seules la pression p et l'angle 0 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 CùVd9yù1 Tde +Vde (5) Le rapport des chaleurs spécifiques y est considéré constant avant la combustion ainsi que dans les premiers instants de la combustion. On pourra par 10 exemple l'estimer é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 O. Pour estimer le volume V et sa dérivée dV , une dt9 table de valeurs donnant le volume V pour différentes valeurs de 0 est utilisée. La température T et sa dérivée d sont calculées sur la base desdits 15 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) avec : To, po, Vo : température, pression et volume initiaux (i.e. à 25 l'injection) To, po,V9 : température, pression et volume à l'angle 0 La température est donc estimée pour chaque valeur de l'angle 0, à 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. 30 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. T Donc à deux instants ou à deux angles différents 20 poVo = poVo (7) To To On en déduit : To=ToP9V0 (8) Po Vo Cette méthode peut être préférée car elle nécessite peu de calculs, étant donné que T n'est calculé qu'une seule fois par cycle et qu'il suffit de Po Vo multiplier ce résultat par p9V9 à chaque nouvelle valeur de l'angle O. Selon un deuxième mode de calcul de T, on suppose aussi que le gaz 5 est parfait (pV = mrT, où mr est contant), mais on en déduit cette fois : 1dT ldp 1dV ±- (9) T dO pdO V dO Cette expression pourrait être réintroduite dans l'expression de C, mais cela ferait réapparaître la dérivée de la pression dé , qui est une source de bruit pour le calcul. Un autre moyen pour exprimer T ~~ 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 0 et O-AO, où A0 est constant et représente le pas du calcul. L'expression (9) devient alors : 1 To ù To-Ao ù 1 Po ù Po-oe + 1 Vo ù Vo-oo (10) To 09 po 09 Vo 09 L'expression suivante de la température T à l'angle 0 est alors déduite, en fonction de la température à l'angle 0-DO : To To-o 1 (11) Po-Ao + VO-DO -1 ^ Po Vo La température est donc estimée de manière récursive pour chaque valeur de l'angle 0, à partir d'une température initiale correspondant à la température 20 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 0 est donnée par la formule : 25 Co = 1 Po To ûTo-oo + Po Ve-oo (12) ù yû1\ TG 09 ) Vo 09 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 V0). 10 15 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 0 une valeur de volume. Dans un système embarqué il est cependant intéressant de diminuer 5 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 Vp donnée par l'expression : 10 Vp=Vo+KpO2, où Vo est constant et est le volume au PMH (volume mort), Kp un coefficient à ajuster de telle façon que la courbe de Vp soit très proche à la courbe de V au voisinage du PMH. Vp peut ainsi être partout substitué à v dans les calculs exposés 15 précédemment : pour le calcul des températures To et To_oo, et dans le terme po Ve ù Vo_Ao Vo 00 Du temps de calcul est ainsi gagné et des ressources sont économisées. Une fois que l'étape de calcul de Co a été réalisée, une détection du 20 SOC est misè en oeuvre, au moyen de deux étapes de comparaison de Co avec des valeurs seuil respectives "SI" et "S2". Deux tests sont ainsi successivement effectués. Le premier test observe si la variable Co dépasse une valeur S1 au-delà de laquelle la combustion pourrait avoir commencé. En effet, la pente prise par 25 Q en 0 est alors suffisamment grande pour estimer que, a priori, la combustion est détectée. Le deuxième test compare Co et C0_oe en observant si la valeur absolue de la différence entre Co et Co_oo est inférieure à une valeur seuil S2, S2 étant une valeur limite au-delà de laquelle on considère que la combustion ne peut 30 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 Co à sa suivante. Le deuxième test permet 35 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 "Si" 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 AO est plus importante. Il est à remarquer qu'on peut appliquer le principe précédemment expliqué au cas multi-injection avec au moins 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, et appliquer des seuils différents aux deux détections. Une fois la détection du SOC effectuée, l'angle vilebrequin O correspondant est enregistré. 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. L'indice de cétane est alors estimé. En premier lieu, une étape de moyennage des valeurs de ces n angles mémorisées est mise ne oeuvre, donnant une valeur O,,,oy. En second lieu, le délai 8 d'inflammation ou d'auto-inflammation est ensuite calculé. Ce délai S 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 8. S =SOC -SOI (13) 10 15 Une autre méthode de calcul de la valeur de début de combustion peut être utilisée, toujours à partir de la variable Co = d B Cette méthode consiste à calculer la valeur d'angle CA10 qui est égale à l'angle où Q = d8 atteint 10% de son amplitude maximale AM, cette valeur CA10 étant enregistrée comme étant la valeur angulaire de début de combustion pour le cycle considéré. Avantageusement, cette amplitude maximale AM est mesurée en établissant la différence entre : -les valeurs maximale et minimale de Q : ou - la valeur maximale de Q et zéro ; ou - zéro ou la valeur minimale de Q, et l'angle de fin d'intégration de Q si celui-ci correspond à la valeur maximale de Q. Grâce à cette valeur de délai d'auto-inflammation, le carburant qui est en train d'être consommé par le moteur est classé parmi un certain nombre de catégories de carburants : - Carburant donnant un délai Dl correspondant à un indice de cétane I1 ; 20 cétane I2 ;Carburant donnant un délai D2 correspondant à un indice de - Carburant donnant un délai D3 correspondant à un indice de - Carburant donnant un délai Dx correspondant un indice de cétane - Carburant donnant un délai Dz correspondant à un indice de Pour chacune de ces catégories, une courbe de correction de l'angle d'avance à l'injection (SOI con) aura été identifiée lors de phase de mise au point du 30 moteur sur banc moteur. Cela signifie que l'on aura préalablement déterminé, lors de la mise au point du moteur sur banc, les valeurs d'avance à l'injection pour z carburants présentant autant d'indice de cétane différents. L'hypothèse est faite ici qu'il y a une relation linéaire entre SOI corrigé (SOI con) et SOI de base (SOI base) pour différentes valeurs de D. Il est 510 également pris comme hypothèse (vérifiée expérimentalement dans le cas de moteur Diesel à fort taux d'EGR et faible température de combustion) que cette relation est linéaire quelque soit la charge et le régime du moteur. Ce SOI corrigé sera celui qui sera optimum en terme de polluants et de bruit émis par le moteur. La cartographie de la figure 2 montre un exemple d'identification de cette cartographie sur banc moteur pour un carburant donné. Sur cette figure, la droite dénommée "linéaire (série 1)" est celle qui résulte de la détermination des SOI con sur banc moteur, tandis que la courbe "série 1" est celle qui résulte de la détermination des SOI base en conditions réelles. La sortie de cette cartographie fournie la nouvelle valeur de SOI qui sera optimum en terme d'émissions polluante et de bruit du moteur par rapport au carburant en train d'être consommé. Elle sera appliquée par le moteur tant que les conditions d'activation de la stratégie ne seront pas à nouveau remplies. 15 REVENDICATIONS 1. Procédé de recalage de l'avance à l'injection en fonction de l'indice de cétane du carburant d'un moteur à combustion interne, sur un nombre n de cycles de déplacement d'un piston dans un cylindre, caractérisé par le fait qu'il comprend les étapes suivantes pour différents angles successifs que prend le vilebrequin du piston dans chaque cycle : a) calculer, sur banc moteur, lors de la mise au point dudit moteur, les valeurs de délai d'auto-inflammation de z carburants de référence présentant chacun un indice de cétane "I1, I2, ...,IX,..., I," , x étant compris entre 1 et z ; b) en déduire les valeurs d'avance à l'injection correspondantes "SOI eOfC" . c) calculer, pour un carburant en train d'être consommé et présentant un indice de cétane IX, les valeurs de délai d'auto-inflammation correspondantes ; d) en déduire les valeurs d'avance à l'injection correspondantes 15 "SOI base" e) affecter à chaque "SOI base", la "SOI cari correspondante pour un même indice de cétane, et l'enregistrer comme nouvelle valeur d'avance à injection. 2. Procédé selon la revendication 1, caractérisé par le fait que l'on calcule les valeurs de délai d'auto-inflammation à partir d'une valeur de début de 20 combustion, qui est elle-même déterminée à partir de l'évaluation d'une variable "Ce" à un angle "0" du vilebrequin dans la chambre de combustion du cylindre, "Co" correspondant à la dérivée angulaire de la quantité d'énergie "Q" à un angle "0", qui est calculée à partir de l'expression suivante : dQ V 'pdT' dV de y-1T de) de 25 - "p" étant la pression dans le cylindre à l'angle "0", - "V" étant le volume de la chambre de combustion à l'ange "0", - "T" étant la température dans la chambre de combustion à l'angle "0", - "y" étant le rapport de chaleur spécifiques " cp/c, ". 30 3. Procédé selon la revendication 2, caractérisé par le fait qu'il consiste à :- comparer chaque "Co" à une valeur seuil "Sl" au-delà de laquelle la combustion pourrait avoir commencé ; - comparer "1 Co - Co_oe 1" à une valeur "S2", "C0_oe" étant une variable précédemment évaluée pour l'angle "O-A0", "S2" étant une valeur limite au-delà de laquelle "Co" est considérée comme une valeur non suffisamment stable pour identifier une combustion ; - 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 Co, la valeur "0" est enregistrée comme étant la valeur angulaire de début de combustion pour le cycle considéré. 4. Procédé selon la revendication 2, caractérisé par le fait qu'il consiste à calculer la valeur d'angle CA10 qui est égale à l'angle où dQ de enregistrée comme étant la valeur angulaire de début de combustion pour le cycle 15 considéré. 5. Procédé selon la revendication 4, caractérisé par le fait que ladite amplitude maximale AM est mesurée en établissant la différence entre : -les valeurs maximale et minimale de Q : ou 20 - la valeur maximale de Q et zéro ; ou - zéro ou la valeur minimale de Q, et l'angle de fin d'intégration de Q si celui-ci correspond à la valeur maximale de Q. Q= atteint 10% de son amplitude maximale AM, cette valeur CA10 étant