FR 2681426 A1 19930319 FR 9111274 A 19910912 PROCEDE ET DISPOSITIF DE MESURE DU COUPLE D'UN MOTEUR THERMIQUE A COMBUSTION INTERNE TENANT COMPTE, NOTAMMENT, DE LA RECIRCULATION DES GAZ D'ECHAPPEMENT ET DES GAZ BRULES RESIDUELS ET DE L'EXCES DE COMBURANT L'invention se rapporte à un procédé et un dispositif de mesure de la valeur du couple gaz moyen par combustion d'un moteur thermique à combustion interne et, plus précisément, à des procédés et dispositifs de correction de la valeur de ce couple gaz moyen par combustion en fonction de certains paramètres de fonctionnement tels que le taux de recirculation des gaz d'échappement, le taux des gaz brûlés résiduels et l'excès de comburant (air). Sur un moteur thermique à combustion interne, le couple produit est la résultante de nombreux paramètres et reflète l'adéquation des réglages Il constitue la grandeur de sortie d'un système dynamique complexe Sa mesure est relativement simple et courante sur un banc d'essais, mais le coût de cette mesure est relativement élevé En outre, les bancs de mesure usuels ne délivrent que les valeurs moyennes sur plusieurs cycles du couple du moteur en fonctionnement stabilisé De telles mesures de valeurs moyennes sur un grand nombre de combustions sont insuffisantes à bien des égards, par exemple pour optimiser certains réglages du moteur ou pour en diagnostiquer certains défauts de fonctionnement Parmi ces défauts, on notera les ratés de combustion dont la détection et l'évaluation statistique sont prévues par les nouvelles réglementations internationales. Pour atteindre ces derniers objectifs, l'analyse quantitative du couple gaz moyen par combustion produit par les combustions du mélange gazeux dans les différents cylindres des moteurs thermiques s'impose. Une telle analyse n'a jusqu'à présent été pratiquée qu'en laboratoire ou sur des moteurs de très grande puissance et elle a généralement fait appel à des relevés de pression dans les chambres de combustion. Cette technique présente l'inconvénient majeur de ne pouvoir être mise en oeuvre que sur des moteurs conçus (ou spécialement modifiés) pour permettre une implantation de capteurs de pression En outre, elle n'est, à l'évidence, pas immédiatement applicable aux moteurs montés en série sur véhicule et, par ailleurs, elle le demeurera aussi longtemps que la fiabilité, le prix, la durée de vie et la commodité de mise en oeuvre des capteurs de pression ne seront pas conformes aux impératifs économiques de l'industrie automobile. Dans une demande internationale de brevet PCT, déposée par Motorola Inc et publiée le 20 juin 1990 sous le No WO 90/07051, est décrit un système de gestion électronique du fonctionnement d'un moteur thermique à combustion interne Ce système est fondé sur le postulat que la valeur instantanée de la période de défilement des dents d'une couronne de mesure, solidaire du volant d'inertie du moteur, devant un capteur fixe, correspond à la mesure de la puissance instantanée successivement produite dans chacun des cylindres du moteur Le signal ainsi produit par le capteur est traité en conséquence. Ce traitement consiste: ( 1) à mesurer la période instantanée di de défilement des dents de la couronne devant le capteur, ( 2) à respectivement multiplier les périodes di des dents afférentes à chacun des cylindres par des facteurs de pondération Pl déterminés appartenant à une suite correspondant à un critère de fonctionnement particulier du moteur tel que le cliquetis ou la puissance (la procédure de détermination de ces facteurs n'est pas décrite) ( 3) à faire la somme des résultats di Pl obtenus pour chaque cylindre, ( 4) à comparer cette somme à une valeur particulière prise comme référence, ( 5) à déduire, le cas échéant, du résultat de cette comparaison la présence de cliquetis dans un cylindre ou le déséquilibre de puissance d'un cylindre par rapport aux autres et ( 6) à modifier en conséquence l'alimentation du cylindre en mélange air-carburant. Le but du système ainsi décrit est de supprimer tout cliquetis dans les cylindres du moteur et/ou d'équilibrer les puissances instantanées fournies par chacun d'eux De ce fait, la valeur absolue des différentes grandeurs intermédiaires obtenues n'est ni recherchée ni trouvée. Dans la demande de brevet français déposée ce jour par la demanderesse et intitulée:"PROCEDE ET DISPOSITIF DE MESURE DU COUPLE D'UN MOTEUR THERMIQUE A COMBUSTION INTERNE", il a été décrit un dispositif pour produire une valeur Cg, représentative du couple gaz moyen engendré par chaque combustion du mélange gazeux dans le ou les cylindre(s) d'un moteur thermique à combustion interne Ce dispositif comprend: des repères de mesure disposés sur une couronne solidaire du volant d'inertie du moteur ou solidaire du vilebrequin; des moyens pour définir une référence d'indexation par tour de volant ou par tour d'arbre à cames pour les repères; un capteur de défilement des repères, monté fixe au voisinage de la couronne; des moyens pour produire une valeur primaire di représentative de la durée de défilement de chacun des repères devant le capteur; des moyens pour traiter cette valeur primaire di et ainsi produire deux valeurs secondaires respectivement représentatives de la vitesse angulaire moyenne nm sur des repères au cours d'une période des combustions dans le moteur et, pour la même période de combustion, de la projection Ecos 4, sur la ligne de référence de phase, des groupes de repères respectivement afférents aux combustions dans les différents cylindres du moteur, de la composante alternative E de la vitesse angulaire instantanée ni des repères, à la fréquence des combustions dans le moteur; des moyens pour combiner ces deux valeurs secondaires, suivant une relation: Cg a lu Ecos 4 + b 2 m et ainsi produire la valeur recherchée, les termes a et b étant des constantes déterminées expérimentalement. Dans le cas particulier d'un moteur à essence quatre cylindres à quatre temps qui fonctionne suivant le cycle de Beau de Rochas, chacun des quatre temps du cycle (admission, compression, combustion- détente, échappement) se produit au cours d'un demi-tour particulier du volant d'inertie solidaire de vilebrequin du moteur L'énergie cinétique acquise par le système concerné (vilebrequin, volant et masses alternatives) est le résultat des différents couples instantanés positifs et négatifs, exercés sur lui pendant chacun des différents temps du cycle de fonctionnement du moteur. Le couple Cg, calculé par demi-tour dans le cas particulier ci- dessus, est engendré pendant les phases de compression et de combustion- détente des mélanges gazeux respectivement emprisonnés dans deux cylindres à combustions contiguës Le moteur étant équipé de quatre cylindres dans le cas envisagé, les deux autres cylindres sont en phases admission et échappement et on peut démontrer, en première approximation, qu'ils n'interviennent pas sur la valeur de Cg du fait de la projection cosi Les phases de compression étant moins sujettes à des variations, à des dispersions que les phases de combustion-détente, on peut dire pour des questions de clarté de présentation, que la valeur de Cg calculée sur un demi-tour est celle relative au cylindre qui était en combustion On peut dire aussi que le couple Cg est une estimation du couple gaz moyen sur les deux phases, compression et combustion-détente, du même cylindre emprisonnant donc le même mélange gazeux. Quel que soit le moteur, Cg est, d'une part, normalement positif et, d'autre part, équilibré par la somme des couples moyens résistants provenant des variations des inerties rotatives du moteur, des différents frottements internes au moteur, des dépressions d'admission et des surpressions d'échappement des autres cylindres et, bien entendu, du couple utile appliqué à la charge, c'est-à-dire au véhicule. Le théorème de l'énergie cinétique enseigne que la variation élémentaire de l'énergie cinétique d'un système égale le travail élémentaire des couples auxquels il est soumis En appliquant ce théorème au système principalement concerné, à savoir le vilebrequin, le volant d'inertie et les masses alternatives impliquées et en supposant que l'ensemble vilebrequin- volant est indéformable, en faisant une analyse harmonique du résultat obtenu, la charge disposée en aval du volant d'inertie étant considérée comme totalement découplée pour la fréquence considérée, on démontre que la composante alternative Ca du couple gaz à la fréquence des combustions dans le moteur, est liée à la vitesse angulaire moyenne am du volant d'inertie et à la projection, sur une ligne de référence de phase, de la composante alternative E de la vitesse angulaire instantanée ai de ce volant, à la fréquence des combustions du moteur, par la relation linéaire suivante: Ca = al m Ecos 9 + bl flm dans laquelle al et bl sont deux termes constants, le premier proportionnel à l'inertie rotative du moteur et le second, fonction du moment d'inertie des masses alternatives et, dans certains cas de mesure, des irrégularités de position des repères de la couronne de mesure. Par ailleurs, une étude expérimentale au banc d'essais de plusieurs moteurs thermiques à combustion interne a montré que le rapport Cg/Ca de la valeur moyenne Cg du couple gaz et de la composante alternative Ca du même couple gaz, varie peu pour les différents régimes et pour les différentes charges de fonctionnement du moteur concerné et ce d'autant mieux que la richesse (r) du mélange gazeux est constante ou, supérieure ou égale à l'unité. En conséquence, l'expérimentation et le calcul pour un système vilebrequin-volant indéformable, ont permis d'établir la relation Cg = -a nm Ecos 4 + b 2 Qm qui est à l'origine du procédé et du dispositif de mesure pour les bas régimes du moteur tels que décrits et revendiqués dans la demande de brevet précitée. Dans le cas des hauts régimes de fonctionnement du moteur, cette relation de base doit toutefois subir une correction En effet, à haut régime, les liaisons élastiques de torsion qui existent tout particulièrement celles entre le vilebrequin et le volant d'inertie du moteur, présentent une fonction de transfert fréquentiel qui modifie le vecteur E représentatif de la composante alternative fondamentale de la vitesse angulaire instantanée du volant, à la fréquence fi des combustions dans le moteur Comme cette liaison comporte un terme d'amortissement négligeable, sa fonction de transfert peut s'écrire Z = 1 + (fi/fr)2, le terme fi étant la fréquence de combustion proportionnelle à la fréquence fv de rotation du volant (fi = fv 2 k/t) avec k le nombre de cylindres et t le nombre de temps du cycle de fonctionnement du moteur et fr, la fréquence de résonance mécanique en torsion de la liaison vilebrequin/volant En pratique, cette fréquence fr est égale à quatre à six fois la fréquence de rotation maximale en Hertz du moteur, de sorte que le terme z = (fi/fr)2 relatif à la fréquence de combustion fi sera souvent faible et qu'un simple facteur de correction ( 1-z) pourra être appliqué à la constante "a" de la relation de base visée ci-dessus, pour que la valeur Cg calculée soit corrigée de l'effet "torsion vilebrequin/volant". Dans la demande de brevet français précitée, il n'a pas été tenu compte de l'influence, pour le calcul de Cg de la recirculation des gaz d'échappement (RGE), de la présence des gaz brûlés résiduels (GBR) et de l'excès de comburant (air) dans les cylindres. Aussi, le but de la présente invention est de compléter le procédé et le dispositif de la demande de brevet français précitée pour calculer la valeur Cg du couple gaz moyen par combustion en tenant compte du taux de recirculation des gaz d'échappement (RGE) et du taux des gaz brûlés résiduels (GBR). L'invention concerne donc, dans un système pour produire une valeur Cg représentative du couple gaz moyen par combustion, un dispositif pour corriger la valeur Cg en fonction du taux total (q) de recirculation des gaz d'échappement et des gaz brûlés résiduels dans les cylindres du moteur, caractérisé en ce qu'il comprend des moyens de calcul pour produire un facteur de correction Q 2 = ( 1-p 2 q), dans lequel P 2 est une constante proche de 0,5 à déterminer expérimentalement et q, la valeur dudit taux total, puis appliquer ce facteur Q 2 à la valeur Cg, le cas échéant préalablement corrigée en fonction de la richesse r du mélange gazeux, pour fournir une valeur de mesure corrigée (Cg)r Ou (Cg)rq- L'invention concerne également un dispositif pour corriger la valeur Cg en cas d'excès de comburant dans les cylindres, caractérisé en ce qu'il comprend: des premiers moyens de calcul pour produire une valeur intermédiaire Ca=Cg/u, avec u description d'une forme de réalisation particulière de l'invention, faite ci-après à titre d'exemple non limitatif, en référence aux dessins annexés dans lesquels la figure 1 est un schéma des différents éléments constitutifs du dispositif de mesure du couple gaz moyen par combustion selon la demande de brevet français précitée, et les figures 2, 3 et 4 sont des schémas d'étages effectuant des corrections selon l'invention. Pour permettre la compréhension de la présente invention, il est nécessaire de décrire, en relation avec la figure 1, le dispositif ou système de calcul de la valeur Cg représentative du couple gaz moyen par combustion correspondant à la demande de brevet français précitée. Sur cette figure 1, est représenté le circuit 10 de mesure du couple gaz moyen Cg, produit par chaque combustion du mélange gazeux dans un moteur thermique à quatre temps (t= 4) et à quatre cylindres (k= 4), équipé d'une couronne de mesure dentée 12, solidaire du volant d'inertie, par exemple la couronne de type standard montée sur les moteurs à allumage électronique. A titre d'exemple, la couronne 12 comporte sur sa périphérie trente-huit dents identiques, telles que celles référencées 14 et 16, qui sont réparties en deux séries de dix-neuf dents séparées par deux dents de référence, telles que celles référencées 18 et 20, qui sont diamétralement opposées et qui ont une largeur triple de celle des autres dents En fait, la couronne comporte: n = 2 x 19 + 2 x 3 = 44 repères équidistants appelés dents réelles ou virtuelles de même module dont deux sont diamétralement opposées pour servir de référence ou d'origine d'indexation permettant de numéroter les dents et, notamment de déterminer la dent do qui sera définie ci-après. A la couronne 12 est associé un capteur fixe 22, par exemple à réluctance variable, adapté à délivrer un signal alternatif 24 de fréquence proportionnelle à la vitesse de défilement des dents de la couronne, c'est-à-dire proportionnelle à la vitesse instantanée du volant. La position angulaire de ce capteur 22 par rapport aux dents d'indexation 18 au moment o le piston est au point mort haut est connue ou relevée, ce qui permet de déterminer la dent do comme étant celle qui défile devant le capteur pendant le passage du piston du cylindre concerné à son point mort haut de combustion. Le signal, délivré par le capteur 22, est appliqué à l'entrée d'un circuit de mise en forme 26, adapté à délivrer des signaux 28 à fronts raides, de durée égale à la période instantanée di des signaux incidents 24 l'indice i variant de O à 21 au fur et à mesure que les dents défilent devant le capteur Chaque période di du signal ainsi produit correspond à la durée de passage d'une dent, soit un plein et un creux, devant le capteur 22 Pour ce qui concerne les signaux incidents produits par les dents d'indexation 18, le circuit de mise en forme 26 les transforme de la même façon en un signal à flancs raides, de durée sensiblement triple de celle des signaux afférents aux autres dents ( 14, 16). Les signaux 28 sont appliqués à un étage 30 de mesure et de calcul des périodes instantanées di de défilement des dents réelles et virtuelles de la couronne de mesure 12 devant le capteur 22. A cet effet, à titre d'exemple de réalisation de l'étage 30, un compteur reçoit de façon continue des impulsions chronométriques de fréquence élevée ( 1 à 20 M Hz, par exemple), produites par une horloge 22, par exemple, à quartz, et il délivre en sortie, des valeurs numériques représentatives des nombres il d'impulsions d'horloge comptées, donc de la "date" (modulo débordement du compteur) de chaque transition creux-plein des signaux à fronts raides, produits par l'étage de mise en forme 26 Ce compteur transmet alternativement la "date" ti 1 qu'il contient à un premier registre et la date suivante ti à un second registre Ces deux registres sont suivis par un étage de calcul adapté à calculer la durée di = (ti-ti) de la dent d'indice i et à en adresser le résultat à une mémoire-tampon 34 De la sorte, aucune impulsion d'horloge n'est perdue et l'erreur de mesure des durées di est minimisée. Pour ce qui est du traitement de chaque signal 28 produit par l'une des dents longues d'indexation 18, 20 dans le cadre de l'exemple décrit sa valeur sera divisée par trois et le résultat, affecté aux trois dents virtuelles correspondantes avec, le cas échéant, prise en compte des restes De la sorte, l'étage 30 de mesure et de calcul des périodes di, adresse à la mémoire-tampon 34, des séries successives de nc=n/2 (soit 22 dans le cas de l'exemple décrit) valeurs numériques di, respectivement associées à n/2 coefficients de pondération contenus dans une mémoire permanente 41 et définis par le rang i de la dent réelle ou virtuelle concernée. Dans le cas du moteur à quatre temps et quatre cylindres de l'exemple, à chaque demi-tour du volant, les valeurs numériques contenues dans la mémoire-tampon 34 sont remises à jour Les valeurs di, ainsi mémorisées en 34, à la fin d'un demi-tour de combustion-détente (dans l'exemple décrit), sont appliquées à un étage de calcul 36, adapté à calculer la durée T 4 (l'indice 4 correspond aux quatre combustions par cycle) de chaque demi-tour du volant, en accord avec la relation T 4 = E di o La dernière valeur numérique T 4 calculée dans l'étage 36 est appliquée à une mémoire-tampon 38 et elle y demeure jusqu'à ce qu'elle soit remplacée par la suivante La vitesse angulaire moyenne du volant d'inertie de la phase combustion-détente qui vient de s'achever, est nm = v/T 4- Les valeurs di mémorisées en 34, sont également appliquées à un autre étage de calcul 40, adapté à réaliser, sous forme numérique, une détection angulaire synchrone d'une projection notée D 4 cos 4 '4 du terme alternatif fondamental de di, rapportée à une position de référence de phase constituée par la première dent d'indice i = O soit do de la série de n/2 = 22 dents affectée à un cylindre donné Pour ce faire, l'étage 40 est adapté à calculer: D 4 cos 4 '4 = ( 1/11) Z di cos(ilr/11) o et à l'appliquer à une mémoire-tampon 42. Les n/2 facteurs de pondération (soit 22 dans l'exemple décrit) de l'expression visée ci-dessus, sont constitués par les onze valeurs successivement prises deux fois, mais avec des signes opposés, par le terme (cos i v/11) lorsque le rang i de la dent concernée varie de O à 10 puis de Il à 21 Ces valeurs sont stockées dans la mémoire permanente 41, associée à l'étage de calcul 40. Le terme Se étant le pas angulaire nominal des dents de la couronne de mesure 12 (dans l'exemple décrit ici 68 = v/22), la vitesse angulaire instantanée ni du volant est ni = Se/di La valeur projetée sur l'axe de référence de la composante alternative fondamentale E 4 de ni est: E 4 cos 14 = ( 1/11) E fl ni cos(iv/11) o En remplaçant ai par sa valeur 68/di et en tenant compte du fait que di varie peu autour de sa valeur moyenne dm = T 4/22 dans l'exemple, on démontre que: E 4 cos'4 =-(am /$e) D 4 cos 4 '4 et que: Cg = -a nm E 4 cos 44 + b am Si les dents (Sei) ne sont pas régulièrement disposées sur la couronne, on démontre qu'il apparaitra un terme proportionnel à am dans E 4 cosg 4 mais ce terme n'est pas gênant car son influence dans Cg sera proportionnelle à 2 m et pourra être compensée en modifiant la valeur du coefficient b. De la sorte, si on décompose en série de Fourier les deux termes di et ni et si on fait sous forme numérique une détection synchrone d'une projection de leur composante alternative fondamentale, on obtiendra Dcose' et Ecos 4 qui sont deux variables dépendantes Dans l'exemple décrit d'un moteur à quatre temps et à quatre cylindres, soit deux combustions par tour de volant et donc quatre combustions par cycle de fonctionnement du moteur, ce qui est rappelé par l'indice 4 utilisé, la composante alternative de valeur E 4 cose 4 de ai est reliée à la composante alternative de valeur D 4 cos 4 '4 de di par la relation -2 E 4 cos P 4 = _m ( 22/ff) D 4 cos I'4 Comme l'étude expérimentale visée plus haut l'a montré pour les bas régimes du moteur, 2 3 2 Cg a= am E'4 + b m = a nm ( 22/z) D'4 + b am en posant E'4 = E 4 cos 4 et D'4 = D 4 cos I'4 Comme nm = v/T 4 on a: Cg = A D'4/T 4 + B/T 4 ce qui est l'expression cherchée avec 21 21 T 4 = E di et D'4 = ( 1/11) E di cos(iv/11) o o les termes A et B étant des constantes propres soit au moteur concerné, soit à tous ceux de son type. En pratique, ces constantes A et B sont déterminées par une mesure du couple du moteur en deux points particuliers ou éventuellement, par une analyse statistique (méthode des moindres carrés) d'un grand nombre de mesures réalisées en laboratoire sur un moteur de référence, conforme à la série des moteurs auxquels est destiné le modèle de couplemètre selon l'invention. En conséquence de ce qui précède, dans le cadre de l'exemple décrit, les grandeurs T 4 et D 14, respectivement élaborées dans les étages de calcul 36 et , sont, à la fin de chaque demi-tour du volant d'inertie du moteur, (synchronisé par au moins une dent d'indexation par cycle moteur), appliquées à un étage de calcul 44 qui élabore Cg, à partir de ces deux grandeurs d'entrée T 4 et D'4 et des constantes A et B mémorisées à cet effet, le terme A étant préalablement corrigé avant usage pour être applicable aux hauts régimes du moteur. En effet, comme cela a été indiqué plus haut, la relation ci-dessus Cg = A D'4/T 4 + B/T 4 n'est valable en toute rigueur que pour les bas régimes du moteur Pour les hauts régimes, une correction du terme A est nécessaire Pour ce faire, dans le cas d'une couronne 12 solidaire du volant d'inertie, un étage de calcul complémentaire 48 est prévu qui élabore un terme corrigé Ac = A( 1-z) avec z = (fi/fr)2, à partir de la valeur de A stockée dans une mémoire permanente 46, de la valeur T 4 fournie par la mémoire-tampon 38 (fi = 1/T 4 = fréquence demi-tour, dans l'exemple décrit) et d'un terme fr (fréquence de résonance en torsion de la liaison vilebrequin/volant) également stocké dans la mémoire 46 Le terme corrigé Ac est appliqué au calculateur 44 qui reçoit T 4 de la mémoire-tampon 38, B de la mémoire permanente 50 et D'4 de la mémoire- tampon 42 et qui calcule le couple gaz moyen sur un demi-tour, selon la relation, valable à tous les régimes du moteur Cg = Ac D'4/T 4 + B/T 4 Dans le but d'obtenir une estimation plus précise de Cg, les coefficients A et B peuvent être identifiés par zone de fonctionnement moteur; ils peuvent aussi être rendus dépendants des réglages moteur, par exemple la richesse r, le taux de recirculation des gaz d'échappement (RGE). Mais une autre façon de faire est de maintenir A et B constants, ce qui permet d'obtenir une première valeur brute de Cg puis, soit de faire intervenir une cartographie de correction (de type additif), soit d'appliquer des termes correctifs de façon à affiner la valeur de Cg Quelques exemples pratiques de correction seront décrits ci-après. Dans le cas d'un moteur fonctionnant en mélange pauvre (richesse carburant/air, r de la richesse r du mélange lorsque r Selon la présente invention, en plus de cet étage optionnel 52, un autre étage de calcul additionnel 56 (figure 2) est prévu pour corriger le terme Cg en fonction du taux total q de recirculation des gaz d'échappement (RGE) et de présence des gaz brûlés résiduels (GBR) dans les cylindres du moteur Ce taux q est toujours petit par rapport à l'unité et sa valeur sera soit mesurée au moyen d'un capteur particulier prévu à cet effet (non représenté), soit prélevée dans une mémoire permanente 58 contenant une cartographie déterminée en laboratoire en fonction des variables d'entrée concernées et notamment du régime du moteur déduit de T 4 et de la valeur Cg calculée En fonction des valeurs T 4 et Cg, appliquées en adresses à la mémoire 58, la valeur q à retenir est entrée dans l'étage de calcul 56 qui élabore un facteur de correction Q 2 = (l-p 2 q), avec P 2 proche de 0,5, puis on affecte Cg ou (Cg)r pour produire un couple (Cg)q corrigé en fonction de q ou (Cg)r,q corrigé en fonction de r et de q. Une solution autre que celle décrite dans la demande de brevet précitée de la demanderesse est possible pour effectuer les corrections afférentes au fonctionnement d'un moteur alimenté avec excès de comburant (air) et donc fonctionnant en gaz pauvre (r En effet, pour les moteurs à essence fonctionnant en mélange pauvre à richesse non constante ou pour les moteurs Diesel fonctionnant toujours avec excès d'air, le rapport entre Cg, le couple gaz moyen produit par chaque combustion, et Ca, le couple gaz alternatif à la fréquence des combustions, n'est plus constant. Toutefois, la connaissance de la façon dont il évolue permet d'effectuer la correction nécessaire (correction excès d'air ou correction de remplissage) pour passer de la valeur de Ca, qui va être effectivement mesurée, à la valeur Cg corrigée qui est recherchée. Tout d'abord, en plein gaz pour tous moteurs, le total des taux (RGE +GBR) est minimal et par ailleurs, dans le cas des moteurs à essence r> 1 et dans celui des moteurs Diesel, r est maximale Donc, en plein gaz, les corrections de richesse r du mélange et de gaz brûlés (RGE + GBR) sont minimales Si on les néglige, on retrouve le cas Cg/Ca = constante de la demande de brevet précitée et on peut écrire u (Ca)p g = (cg)p*g. quel que soit N, avec les indices p g signifiant "plein gaz" et le terme N, le nombre de tours/minute du moteur Cette relation permet, comme on l'a vu plus haut, d'identifier les termes A et B de l'expression définissant Cg 3 2 c'est-à-dire Cg = Ac D'4/T 4 + B/T 4 (étage 44) par les mesures de deux points "plein gaz" pris à deux régimes différents Quant au coefficient u, il doit être déterminé à partir de mesures de pressions de chambres de combustion sur un moteur de chaque type: u = (Cg/Ca)p g, le terme u est constant au premier ordre, quel que soit N. En plein gaz et avec r>l, le remplissage R des cylindres est proportionnel à Ca (ou Cg) au premier ordre, ce qui explique (voir ci-après) que le rapport u = (Cg/Ca)p g. soit sensiblement constant quel que soit N en plein gaz. Hors le cas du plein gaz, en faisant fonctionner le moteur en charge partielle à plusieurs régimes, on a constaté expérimentalement que (Cg)p = Ca k R R, avec kp R, un terme constant au premier ordre et R, le remplissage du cylindre en phase compression Le terme R représente le poids total du mélange gazeux dans le cylindre (carburant, comburant, gaz brûlés résiduels ou recyclés) Il peut être estimé à partir d'une cartographie mise en mémoire et/ou de mesures faites par les capteurs adéquats déjà installés sur le moteur Le terme k R est, à tous régimes N et à toutes charges, une constante au premier ordre que l'on peut identifier,en particulier, en plein gaz (p g) par la relation k R = (l-u) (Ca)p g /Rp g. Quant au calcul de Ca, il reste inchangé en charge partielle: Ca = Cg/u = (Ac D'k/T? + B/Tk) u La figure 3 représente le schéma fonctionnel des circuits complémentaires d'un dispositif de mesure du couple gaz moyen par combustion, assurant la correction de remplissage en cas de fonctionnement en mélange pauvre à richesse non constante Selon la figure 3, les circuits 52 à 58 des figures 1 et 2 sont remplacés par des circuits 60 à 66 Le circuit de calcul 60 reçoit de l'étage 44, un signal Cg non corrigé: Cg = A D'k/Tk + B/Tk et d'une mémoire 62, le coefficient u Le circuit 60 calcule un terme intermédiaire Ca = Cg/u et le résultat est appliqué à un autre étage de calcul 64 Une mémoire 66 contient le terme k R visé ci-dessus et, en outre, la cartographie des remplissages R, exprimés à partir de grandeurs disponibles, calculées ou mesurées par des capteurs adéquats La mémoire 66 applique les valeurs k R et R au circuit de calcul 64 qui produit un signal de mesure du couple gaz moyen par combustion d'un moteur fonctionnant en mélange pauvre: 1 k 3 2 (Cg)p = Ca k R R (Ac D'k/Tk + B/Tk) k R R Dans le cas d'un moteur Diesel, le remplissage est constant au premier ordre et l'expression de Cg en fonction de Ca se ramène à (Cg)D = Ca KR, le terme KR étant une constante au premier ordre mais, pour plus de précision, KR peut être cartographiée De la sorte, pour un moteur Diesel, le couple gaz moyen par combustion est: (Cg)D = (Ac D'k/Tk + B/Tk) KR. u La figure 4 représente le schéma des modifications qui découlent de ce qui précède Selon la fig 4, les circuits 64-66 de la fig 3 sont remplacés par des circuits 68-70 Le circuit 68 est une mémoire qui contient la constante KR ou sa cartographie et le circuit 70, un étage de calcul qui produit le signal de mesure (Cg)D = Ca KR, applicable aux moteurs Diesel. Les différents étages de calcul visés ci-dessus sont, au fur et à mesure de leur intervention dans le processus décrit, constitués par un microprocesseur opérant sous la dépendance d'un logiciel d'application L développé à cet effet et contenu dans la mémoire permanente 50 Une autre manière de faire est, bien entendu, de réaliser au moins une partie des opérations prévues, au moyen d'un calculateur électronique analogique dont les différents étages assurent chacun une fonction particulière. L'ensemble des circuits représentés peut être réalisé sous forme de circuit intégré ou de circuit hybride, spécifique du problème à résoudre. Dans l'exemple décrit ci-dessus, le moteur concerné est un moteur à essence à quatre temps et à quatre cylindres et le nombre de repères réels ou virtuels de la couronne de mesure est 44 L'invention n'est, bien entendu, pas limitée à cet exemple particulier Elle concerne en effet tous les moteurs thermiques à combustion interne quel qu'en soit le type: à quatre temps, à deux temps, à essence, à gazole (diesel) ou à tout autre carburant et quel que soit le nombre de cylindres ou le nombre de repères de la couronne de mesure. Si l'on désigne par t, le nombre de temps par cycle de fonctionnement du moteur, par k, le nombre de ses cylindres et par n, le nombre de repères réels ou virtuels de la couronne de mesure, le nombre de repères qui correspond à chacune des périodes angulaires (vt/k) des combustions dans le moteur et qui sera pris en compte pour la détection angulaire synchrone de la vitesse angulaire instantanée à la fréquence des combustions est nc = n t/( 2 k), ce qui donne combustion par combustion et donc cylindre par cylindre: (nc) Tk = E di. o Dans ces conditions, le terme 21 D 4 cosf'4 = ( 1/11) E di cos (i v/11) =D 4 o visé ci-dessus, s'écrit sous la forme généralisée: (nc-1) Dkcos P'k = ( 2/nc) E di cos (i 2 v/nc) = D'k o L'invention a été décrite en relation avec un exemple particulier de réalisation qui fait appel aux variables intermédiaires di, T 4 et D'4 et dans lequel les calculs de di, T 4 et D'4 étaient effectués de manière numérique mais le calcul de Cg et de ses corrections peut être réalisé par des circuits analogiques classiques. Ainsi, le signal primaire 24 dont la fréquence est modulée par les variations de vitesse n, contient le vecteur E à la fréquence fi Pour extraire le vecteur E et plus précisément sa projection Ecose, il faut utiliser un convertisseur fréquence-tension suivi d'un détecteur synchrone angulaire à la fréquence fi avec une phase convenable puis calculer: Cg = -a Sm Ecos 4 + bfm L'invention a été décrite en référence à des dispositifs de correction mais elle concerne aussi des procédés de correction de la valeur Cg en fonction de certaines caractéristiques de fonctionnement du moteur. Les opérations de ces procédés découlent de manière évidente de la description des dispositifs de correction faite en relation avec les figures 2, 3 et 4 et ne seront donc pas décrites de manière séparée. REVENDICATIONS 1 Système de mesure de la valeur numérique Cg, représentative du couple gaz moyen par combustion du mélange gazeux dans les k cylindres d'un moteur à combustion interne comportant t temps par cycle de fonctionnement, ledit système comprenant: un nombre N de repères de mesure ( 14-16) qui sont disposés sur une couronne ( 12) solidaire du volant d'inertie du moteur ou solidaire du vilebrequin; des moyens ( 18- 20) pour définir au moins une référence d'indexation des repères ( 14-16); un capteur ( 22) de défilement des repères, monté fixe au voisinage de la couronne ( 12); des moyens de calcul ( 26-30-32) pour produire une valeur numérique primaire di représentative de la durée instantanée de défilement devant le capteur ( 22) de chacun des repères ( 14-16) réels ou virtuels; des moyens de traitement numérique de la valeur numérique primaire di comprenant, notamment, une unité de calcul ( 36-40-44), un logiciel d'application (L), des mémoires permanentes ( 41-46-50) et des mémoires temporaires ( 34-38-42) l'ensemble de ces moyens de calcul et de traitement numérique étant organisé pour constituer en temps opportun: des moyens de calcul ( 36) pour élaborer à partir de la valeur numérique primaire di, une première valeur numérique secondaire Tk représentative de la durée totale de défilement devant le capteur ( 22) de chaque série de nc repères ( 14-16), afférente à la période angulaire d'une combustion dans le moteur, (nc-1) avec nc = (n t/2 k) et Tk = E di o des moyens de calcul ( 40) pour élaborer, à partir de la valeur numérique primaire di et de facteurs de pondération cos(i 2 r/nc) stockés dans une mémoire permanente ( 41), une seconde valeur numérique secondaire Dkcos V'k représentative de la projection, sur une ligne de référence de phase des repères ( 14-16) correspondant à l'origine des périodes angulaires des combustions, de l'amplitude de la composante alternative des durées instantanées di de défilement devant le capteur ( 22) des repères de mesure ( 14-16), à la fréquence des combustions dans le moteur, avec (nc-l) Dkcosi'k = ( 2/nc) E di cos(i 2 v/nc). o des moyens de calcul ( 44) pour élaborer la valeur numérique recherchée Cg, défini par la relation Cg = A (Dkcos 4 'k)/Tk + B/Tk =A D'k/Tk + B/Tk à partir des deux valeurs numériques secondaires Dkco 54 'ket Tk et de deux constantes A et B expérimentalement déterminées et stockées dans les mémoires permanentes ( 46-50), caractérisé en ce que ledit système comprend, en outre, un dispositif de correction de la valeur numérique Cg en fonction du taux total q de recirculation des gaz d'échappement (RGE) et des gaz brûlés résiduels (GBR), ledit dispositif de correction comprenant des moyens de calcul ( 56) pour produire un facteur de correction Q 2 = ( 1-p 2-q), dans lequel P 2 est une constante proche de 0,50 à déterminer expérimentalement et q, la valeur dudit taux total, puis appliquer ce facteur Q 2 à la valeur Cg, le cas échéant préalablement corrigée en fonction de la richesse r du mélange gazeux, pour fournir une valeur numérique corrigé (Cg)q Ou (Cg)rq. 2.Système selon la revendication 1, caractérisé en ce que dans ledit dispositif de correction, les valeurs du taux total q, afférentes à un type donné de moteurs, sont mesurées ou cartographiées notamment en fonction du couple Cg produit et du régime moteur déduit de Tk puis stockées dans une mémoire permanente ( 58) associée aux moyens de calcul ( 56). 3 Système selon le préambule de la revendication 1, caractérisé en ce que ledit système comprend, en outre, un dispositif de correction de la valeur numérique Cg dans le cas d'un moteur fonctionnant avec excès de comburant, c'est-à-dire dans le cas d'un mélange gazeux pauvre, r constante mémorisée ( 62), déterminée expérimenta- lement; des moyens de calcul ( 64) sont prévus pour produire la valeur recherchée (Cg)p = Ca-k R R, avec k R une constante mémorisée ( 66) déterminée expérimentalement et R, le remplissage des cylindres, déterminé expérimentalement en fonction de grandeurs déjà calculées ou mesurées, et formant une cartographie mémorisée ( 66). 4 Système selon le préambule de la revendication 1, caractérisé en ce que ledit système comprend, en outre, un dispositif de correction de la valeur numérique Cg dans le cas d'un moteur Diesel pour obtenir une valeur corrigée (Cg)D, ledit dispositif de correction comprenant: des moyens de calcul ( 60) sont prévus pour produire une valeur intermédiaire Ca = Cg/u, avec u constante mémorisée ( 62), déterminée expérimenta- lement; des moyens de calcul ( 70) sont prévus pour produire la valeur recherchée (Cg)D = Ca-KR, avec KR, une valeur déterminée expérimentalement, mémorisée ( 68) comme une constante ou mieux comme une cartographie, en fonction de grandeurs déjà calculées ou mesurées. 5 Procédé de correction de la valeur numérique Cg, représentative du couple gaz moyen par combustion et obtenue à l'aide du système du préambule de la revendication 1, en fonction du taux total q de recirculation des gaz d'échappement et des gaz brûlés résiduels dans les cylindres du moteur, caractérisé en ce qu'il comprend les opérations suivantes: (el) produire un facteur de correction Q 2 = (l-p 2 q), dans lequel P 2 est une constante proche de 0,50 à déterminer expérimentalement et q, la valeur dudit taux total, puis (e 2) appliquer ce facteur Q 2 à la valeur de mesure Cg, le cas échéant préalablement corrigé en fonction de la richesse r du mélange gazeux quand r ou (Cg)rq- 6 Procédé selon la revendication 5, caractérisé en ce que les valeurs du taux total q, afférentes à un type donné de moteurs, résultent d'une opération de mesure ou de cartographie notamment en fonction du couple (Cg) produit et du régime du moteur déduit de Tk. 7 Procédé de correction de la valeur numérique Cg, représentative du couple gaz moyen par combustion et obtenue à l'aide du système tel que défini dans le préambule de la revendication 1, caractérisé en ce que, -1 pour produire une valeur de mesure du couple gaz moyen par combustion (Cg)p afférent à un moteur fonctionnant avec excès de comburant c'est-à-dire d'un mélange gazeux pauvre (r remplissage des cylindres, déterminé expérimenta- lement en fonction de grandeurs déjà calculées ou mesurées, et formant une cartographie mémorisée ( 66). 8 Procédé de correction de la valeur numérique Cg, représentative du couple gaz moyen par combustion et obtenue à l'aide du système tel que défini dans le préambule de la revendication 1, caractérisé en ce que, pour produire une valeur de mesure du couple gaz moyen par combustion d'un moteur Diesel, (Cg)D, il comprend, en outre, les opérations suivantes: (i) calcul ( 60) pour produire une valeur intermédiaire Ca = Cg/u, avec u fonction de grandeurs déjà calculées ou mesurées.