FR 3027961 A1 20160506 FR 1460609 A 20141104 PROCEDE D'ESTIMATION D'UNE POSITION D'UN PAPILLON D'ARRIVEE DES GAZ POUR LE CONTROLE D'UN MOTEUR A COMBUSTION INTERNE [0001] La présente invention porte sur un procédé d'estimation d'une position d'un papillon d'arrivée des gaz pour le contrôle d'un moteur à combustion interne. [0002] De façon connue en soi, pour assurer une bonne combustion dans les moteurs à essence à injection directe et indirecte, on commence à injecter le carburant durant la phase d'échappement du cylindre considéré. Pour connaître la quantité de carburant à injecter sur ce cylindre, il faut déterminer dès le début de cette phase d'échappement quelle sera la quantité d'air présente dans la chambre de combustion délimitée par le cylindre au moment de la fermeture de la soupape admission, c'est-à-dire la quantité d'air qui va participer à la combustion. On réalise donc pour chaque cylindre, durant la phase d'échappement et durant la phase d'admission, une prédiction de la masse d'air présente dans la chambre de combustion associée au cylindre au moment de la fermeture de la soupape admission. [0003] Cette prédiction de la masse d'air peut être réalisée à partir d'un modèle standard d'estimation de la masse d'air dans le cylindre qu'on alimente avec une pression du répartiteur admission prédite au moment de la fermeture de la soupape admission. [0004] Deux systèmes connus consistent à prédire la pression qui va régner dans le répartiteur admission à l'horizon temporel donné. Cette pression prédite va ensuite alimenter un modèle de débit de cylindre qui va fournir un débit prédit à l'horizon donné. [0005] Un des systèmes de prédiction de pression se base sur l'analyse du gradient de la pression mesurée. On extrapole alors la valeur de pression qui règnera dans le répartiteur admission à l'horizon de la fermeture soupape admission en considérant ce gradient constant et en partant de la pression mesurée courante. [0006] L'autre système se base sur l'utilisation d'un modèle du répartiteur d'admission modélisant les débits entrants (débits issus du papillon et du canister principalement) et sortant (débit dans le cylindre) courants. Ce modèle permet d'estimer la pression courante dans le répartiteur par l'intégration de l'écart entre ces débits. Pour la prédiction de pression, on reprend ce principe de calcul en intégrant l'écart entre ces débits à un horizon correspondant au moment de la fermeture de la soupape admission. [0007] Les deux systèmes connus présentent la même limitation. Ils se basent principalement sur l'analyse des mesures courantes (de pression dans le premier système, et des positions actionneurs pour la modélisation des débits pour le deuxième système) et les extrapolent jusqu'à l'horizon demandé. Or, en tenant compte uniquement de ces mesures, ces systèmes prédisent la pression en considérant les phénomènes figés. L'évolution est donc supposée identique jusqu'à l'horizon visé, ce qui n'est pas forcément le cas tout au long d'un transitoire. [0008] En effet, comme on peut le voir sur l'exemple de la figure 1, à l'instant de calcul t, le gradient de pression Grad _P est calculé à partir de la mesure de pression au point courant et au point précédent. Autrement dit, le gradient d'évolution est mesuré entre deux pas de calcul. La pression prédite P pred est obtenue par extrapolation à l'horizon désiré du gradient calculé Grad P. Sur cet exemple, on voit que l'évolution physique de la pression amène, au terme de l'horizon H, à une pression P_pred très différente de la pression cible P c mesurable à l'instant cible t c. Les deux systèmes connus ne sont donc pas très performants en termes de qualité de prédiction pour les horizons longs H et pour les variables présentant des évolutions non linéaires. [0009] L'invention vise à remédier efficacement à cet inconvénient en proposant un procédé d'estimation d'une position d'un papillon d'arrivée des gaz pour le contrôle d'un moteur à combustion interne, caractérisé en ce qu'il comporte les étapes suivantes: - une étape de mesure d'une position du papillon d'arrivée des gaz, et - une étape de prédiction, sur un horizon temporel de calcul donné, de la position du papillon d'arrivée des gaz en fonction de la position du papillon mesurée et d'un gradient de position entre une position courante et une position précédente, - ladite étape de prédiction de la position du papillon d'arrivée des gaz étant adaptée notamment en fonction de l'horizon de calcul donné et/ou d'informations de consigne et/ou de paramètres physiques du papillon d'arrivée des gaz. [0010] Le procédé selon l'invention qui utilise des informations de consigne, des mesures, et les paramètres physiques du comportement du papillon permet d'obtenir un signal de la position du papillon prédit plus précis et réactif (notamment en début de transitoire) que la méthode de prédiction selon l'état de l'art. [0011] Selon une mise en oeuvre, ladite étape de calcul du gradient de position est constituée par le choix entre: - la mise en oeuvre d'une première stratégie de calcul basée sur un calcul du gradient de la consigne de position du papillon d'arrivée des gaz, - la mise en oeuvre d'une deuxième stratégie de calcul basée sur un calcul du gradient de la mesure de la position du papillon d'arrivée des gaz, ou - la mise en oeuvre d'une combinaison de la première et de la deuxième stratégies de calcul précitées. [0012] De préférence, le choix entre les différentes mises en oeuvre de stratégies de calcul ou de combinaison des stratégies de calcul est effectué en comparant une valeur de l'horizon temporel de calcul donné par rapport à une dynamique connue de déplacement dudit papillon d'arrivée des gaz. [0013] Selon une mise en oeuvre, ledit procédé comporte une étape de saturation du gradient de la position du papillon basée sur des paramètres physiques du papillon d'arrivée des gaz. [0014] Selon une mise en oeuvre, ledit procédé comporte en outre une étape de saturation de la prédiction de position du papillon d'arrivée des gaz. [0015] Selon une mise en oeuvre, l'étape de saturation de la position du papillon d'arrivé des gaz prédite est réalisée de telle façon qu'en phase d'ouverture du papillon d'arrivée des gaz, la position prédite ne peut pas être supérieure à une consigne du papillon, et en phase de fermeture papillon, la position prédite ne peut pas être inférieure à la consigne du papillon. [0016] Selon une mise en oeuvre, ledit procédé comporte une étape de filtrage de la position du papillon prédite. [0017] Selon une mise en oeuvre, ledit procédé comporte une étape d'application d'un filtre du premier ordre. Cela permet de lisser les variations brusques du gradient de pression. [0018] L'invention sera mieux comprise à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent. Ces figures ne sont données qu'a titre illustratif mais nullement limitatif de l'invention. [0019] La figure 1, déjà décrite, est un graphique illustrant la prédiction de pression dans le répartiteur d'admission suivant un procédé selon l'état de la technique au cours d'une évolution transitoire de pression; [0020] La figure 2 est un schéma fonctionnel du procédé selon l'invention illustrant la réalisation d'un bilan des débits entrants/sortants du répartiteur d'admission en vue d'estimer le débit de gaz dans le cylindre; [0021] La figure 3 est un schéma fonctionnel du procédé selon l'invention illustrant le calcul de la prédiction de la position du papillon d'arrivée des gaz; [0022] Les figures 4a et 4b sont des graphiques illustrant la prédiction de la position du papillon respectivement dans le cas d'un transitoire long et d'un transitoire court; [0023] Les figures 5a et 5b sont des graphiques illustrant les deux étapes permettant d'estimer la pression à l'intérieur du répartiteur d'admission ainsi que le débit correspondant des gaz dans le cylindre. [0024] La figure 6 est un schéma fonctionnel illustrant l'utilisation du débit dans le cylindre pour déterminer les consignes de commande des injecteurs. [0025] Les éléments identiques, similaires, ou analogues conservent la même référence d'une figure à l'autre. [0026] La figure 2 illustre de manière schématique la modélisation du bilan du répartiteur d'admission dans le cadre de son application au calcul d'un débit de gaz prédit dans le cylindre Dcyl_pred. [0027] Cette modélisation est basée notamment sur l'utilisation d'un module prédicteur de débits entrants 10 correspondant ici au débit papillon prédit Dpap pred à l'horizon temporel H considéré et d'un module prédicteur de débits sortants 11 correspondant au débit dans le cylindre prédit Dcyl pred à l'horizon de calcul H considéré. [0028] Le module prédicteur de débit entrant 10 est alimenté par un module 12 de prédiction de la position du papillon d'arrivée des gaz Ppap_pred. Comme cela est illustré sur la figure 3, la prédiction de la position angulaire du papillon Ppap pred à un horizon donné H se base essentiellement sur l'observation du gradient de position Grad pos_pap entre la position courante et celle du pas précédent. Ce gradient Grad pos pap pourra le cas échéant être saturé par le module 124 comme expliqué ci-après. La formule relative à l'évolution de la prédiction de la position du papillon Ppap pred est : Ppap_pred(n) = P pap mes(n) + Grad pos_pap* H [0029] Plus précisément, le calcul du gradient de position Grad pos pap est décomposé en deux stratégies de calcul. La première stratégie de calcul mise en oeuvre par le module 121 est basée sur le calcul du gradient de la consigne de position du papillon Grad cons pap consistant à calculer l'écart entre deux positions successives de la position du papillon de consigne Ppap cons. La deuxième stratégie de calcul mise en oeuvre par le module 122 est basée sur le calcul du gradient de la mesure de position du papillon Grad mes_pap consistant à calculer l'écart entre deux positions successives de la position du papillon mesurée Ppap mes. Le choix du gradient de la position du papillon Grad pos_pap est effectué par le module 123 en fonction de la valeur de l'horizon temporel H. [0030] La première stratégie de calcul est ainsi sélectionnée pour les horizons temporels H relativement longs, à faible régime par exemple. En effet, on considère que dans ce cas, l'actionneur aura eu le temps d'atteindre sa position de consigne et donc la position du papillon prédite Ppap_pred sera plus proche de la consigne que de la mesure actuelle. Par contre, pour les horizons temporels courts H à fort régime par exemple, la deuxième stratégie de calcul est sélectionnée. En effet, dans ce cas, la position prédite Ppap pred sera plus proche de la position mesurée actuelle que de la consigne. [0031] Il est à noter qu'il est également possible de choisir la mise en oeuvre d'une combinaison des deux stratégies de calcul précitées. Par exemple, on peut sélectionner la première stratégie de calcul pour le régime minimum (au ralenti, à 700tr/min) et la deuxième stratégie de calcul pour le régime maximal (6000tr/min par exemple), et pour un régime médian (3000tr/min), il est possible de choisir une combinaison de type 50% de la valeur de la première stratégie ajoutée à 50% de la valeur de la deuxième stratégie. [0032] En outre, le choix de la stratégie de calcul ne dépend pas uniquement du régime. Ainsi, dans la pratique, on met en oeuvre le procédé sur trois horizons différents. A un même régime, pour chacun des trois horizons, on ne va pas nécessairement choisir la même stratégie de calcul ou combinaison des deux stratégies. Le choix est effectué en comparant la valeur de l'horizon H exprimé en secondes (prise en compte du régime et de la valeur angulaire de l'horizon) par rapport à la dynamique connue de déplacement du papillon d'arrivée des gaz. [0033] Sur les figures 4a et 4b, on a représenté par la courbe Cl l'évolution de la consigne de la position du papillon en fonction du temps. Par ailleurs, la courbe 02 représente l'évolution de la mesure de la position du papillon. Le point Ppap Gcons correspond à la position prédite à l'horizon H donné à partir du gradient de consigne Grad cons pap. Le point Ppap Gmes correspond à la position prédite à l'horizon H donné à partir du gradient Grad mes pap mesuré entre la position mesuré à l'instant t de calcul et la valeur précédente. [0034] Dans l'exemple d'un horizon H long illustré à la figure 4a, à l'instant de calcul t, il ressort que l'utilisation du gradient de consigne Grad cons_pap (et la saturation par la consigne du papillon) permet de réaliser une prédiction de position Ppap Gcons plus proche de la position cible Ppap c que l'on aura à l'instant visé t c que le calcul utilisant le gradient mesuré Grad mes_pap. [0035] Dans l'exemple d'un horizon H court illustré à la figure 4b, à l'instant de calcul t, on voit que l'utilisation du gradient mesuré Grad mes_pap permet de réaliser une prédiction de position plus proche de la position Ppap c que l'on aura à l'instant visé t c que le calcul utilisant le gradient de consigne Grad cons pap. [0036] Une fois le gradient adapté sélectionné, le module 124 assure une saturation du gradient de la position du papillon Grad_pos pap. Cela permet de sécuriser le gradient de la position du papillon par une limitation d'évolution. Ces saturations sont basées sur des paramètres physiques du papillon d'arrivée des gaz, tels que ses vitesses maximales de déplacement en ouverture et en fermeture. [0037] La prédiction brute de la position du papillon Ppap br est ensuite calculée à partir de la position du papillon mesuré Ppap mes, du gradient de position saturé Grad pos_pap sat et de l'horizon de prédiction H. Les modules 125 et 126 assurent la combinaison de ces données suivant la formule précitée. [0038] La saturation de la prédiction brute de la position du papillon Ppap br est ensuite réalisée par le module 127 en fonction de deux cas pour obtenir Ppap br sat. Dans le premier cas correspondant à une phase d'ouverture du papillon, la position prédite Ppap_pred ne peut pas être supérieure à la consigne de position du papillon. Dans le second cas correspondant à une phase de fermeture du papillon, la position prédite Ppap_pred ne peut pas être inférieure à la consigne de position du papillon. [0039] On finalise ensuite le calcul par un filtrage de la prédiction brute de position du papillon Ppap br sat. A cet effet, un module 128 pourra par exemple appliquer un filtre du premier ordre pour lisser les variations brusques du gradient de pression. [0040] Ce principe de calcul utilisant les informations de consigne Ppap cons, de mesure Pap mes, et les paramètres physiques de comportement du papillon, permet d'obtenir un signal de position du papillon prédit beaucoup plus réactif (notamment en début de transitoire) et plus précis que les techniques de l'état de l'art. [0041] Par ailleurs, comme on peut le voir sur la figure 2, le module 11 prédicteur de débits sortants est alimenté par un module 13 prédicteur de position angulaire des 10 déphaseurs d'arbre à came Pdeph adm_pred, Pdeph ech pred et un module 14 prédicteur de régime moteur Wmth pred. [0042] Plus précisément, le module 13 estime la position angulaire du déphaseur d'arbre à cames à l'admission Pdeph adm pred ainsi que la position angulaire du déphaseur d'arbre à cames à l'échappement Pdeph ech pred. Ces estimations sont effectuées en 15 fonction de mesures des positions angulaires des déphaseurs d'arbre à cames à l'admission Mes deph adm et à l'échappement Mes deph ech et de consignes des positions angulaires des déphaseurs d'arbre à cames à l'admission Cons deph adm et à l'échappement Cons deph ech, ainsi que de l'horizon de prédiction H. [0043] Par ailleurs, le module 14 estime le régime moteur Wmth pred en fonction d'une 20 mesure de régime moteur Wmth mes et de l'horizon de prédiction H. [0044] La réalisation d'un bilan des débits entrants et sortants du répartiteur admission permet de déterminer la pression P_pred qui règne dans le répartiteur d'admission. Le principe de base de ce bilan est qu'en régime stabilisé, la loi physique de conservation des masses appliquée au répartiteur admission fait que le bilan entre ses débits entrants et 25 sortants est nul et que donc la pression est stable. [0045] A cet effet, le module 15 calcule l'écart Delta E/S entre ces débits Dpap_pred et Dcyl pred. Le module 16 transforme alors cet écart Delta E/S en gradient de pression Grad P pred en tenant compte de la température Temp. [0046] Ce gradient Grad P pred permet, à partir du pas de calcul Pcalc de la fonction, 30 et de la valeur précédente de pression prédite Ppred(n-1), de calculer, via les modules 17, 18, une nouvelle valeur de pression prédite Ppred(n). La formule d'évolution est la suivante: P_pred(n) = P_pred(n-1) + Grad P pred * Pcalc [0047] Cette valeur de pression prédite P_pred(n) est ensuite renvoyée vers les modules prédicteurs de débits 10 et 11. Les débits sont ainsi mis à jour et le calcul continue jusqu'à trouver une équilibre, c'est-à-dire lorsque la différence entre les débits entrants et sortants devient nul. [0048] Comme cela est illustré sur la figure 5a, au pas d'initialisation, les estimations de débit papillon Dpap pred init et de débit dans le cylindre Dcyl pred init réalisés à partir de la pression prédite initiale P pred init sont assez différents. La différence des débits Delta E/S est intégrée pour calculer une évolution de la pression prédite P_pred. Cette nouvelle pression P pred est renvoyée vers les modules prédicteurs de débits 10, 11. Il est à noter que sur les figures 5a et 5b, la droite C3 correspond à la caractéristique de débit du cylindre (la droite de remplissage), tandis que la courbe C4 correspond à la caractéristique de débit du papillon à pression constante. [0049] Comme cela est illustré sur la figure 5b, en intégrant progressivement l'écart entre les débits prédits papillon Dpap_pred et cylindre Dcyl_pred, la pression prédite converge jusqu'à une valeur stable P_pred conv, déterminant l'équilibre des débits prédits. A cet instant, on obtient une estimation de la prédiction du débit dans le cylindre Dcyl pred à l'horizon considéré H. [0050] La prédiction de débit cylindre Dcyl pred au moment de la fermeture soupape admission du cylindre considéré est calculé durant toute la phase d'échappement et poursuivit pendant la phase d'admission. [0051] Comme cela est illustré sur la figure 6, le débit cylindre prédit Dcyl pred est transformé, par le module 21, en masse d'air prédite Mair_pred. A cet effet, le module 21 tient compte du régime moteur Wmth mes, ainsi que de la concentration en air dans le répartiteur d'admission Cair. [0052] Cette masse d'air prédite Mair pred est ensuite convertie en masse de carburant de consigne, ce qui sert à calculer et à commander le temps d'ouverture des injecteurs. Le module 22 établit ainsi des consignes de commande Cons comm en fonction de la masse d'air prédite Mair pred et également de la richesse de consigne Rich cons et du facteur d'enrichissement Fenr. [0053] Ainsi, le débit cylindre prédit Dcyl_pred étant calculé à partir du modèle de débit cylindre alimenté par des paramètres principaux prédits (pressions répartiteur admission prédites P pred, positions des déphaseurs arbres à cames prédits et régime moteur prédit Wmth_pred), l'invention permet d'améliorer la finesse de prédiction de débit cylindre Dcyl pred. [0054] En outre, la justesse de la prédiction de pression P pred est sensiblement améliorée par l'utilisation des modèles de prédiction des positions des actionneurs (papillon, déphaseurs arbres à cames admission et échappement) et du régime Wmth à l'horizon donné. En effet, l'horizon de prédiction H est directement intégré dans les modèles de position des actionneurs qui servent à réaliser le bilan du répartiteur d'admission. [0055] Ces modèles s'avèrent robustes en termes de prédiction, dans la mesure où ils peuvent prédire plus précisément le comportement des actionneurs en se basant sur: - la connaissance de la physique des actionneurs, de la prise en compte de leurs paramètres physiques (vitesse d'ouverture par exemple), - l'analyse des gradients d'évolution courants (gradient de la position du papillon mesurée par exemple), ainsi que - l'analyse de la commande qui leur est appliquée (gradient de la position du papillon de consigne par exemple). [0056] On obtient ainsi une prédiction de pression plus juste qu'une simple extrapolation à l'horizon donné d'un gradient de pression considéré comme constant. REVENDICATIONS: 1. Procédé d'estimation d'une position d'un papillon d'arrivée des gaz pour le contrôle d'un moteur à combustion interne, caractérisé en ce qu'il comporte les étapes suivantes: - une étape de mesure d'une position du papillon d'arrivée des gaz (Ppap mes), et - une étape de prédiction, sur un horizon temporel de calcul donné (H), de la position du papillon d'arrivée des gaz (Ppap pred) en fonction de la position du papillon mesurée (Ppap mes) et d'un gradient de position (Grad pos_pap) entre une position courante et une position précédente, - ladite étape de prédiction de la position du papillon d'arrivée des gaz (Ppap pred) étant adaptée notamment en fonction de l'horizon de calcul donné et/ou d'informations de consigne et/ou de paramètres physiques du papillon d'arrivée des gaz. 2. Procédé selon la revendication 1, caractérisé en ce que ladite étape de calcul du gradient de position (Grad pos_pap) est constituée par le choix entre: - la mise en oeuvre d'une première stratégie de calcul basée sur un calcul du gradient de la consigne de position du papillon d'arrivée des gaz (Grad cons pap), - la mise en oeuvre d'une deuxième stratégie de calcul basée sur un calcul du gradient de la mesure de la position du papillon d'arrivée des gaz (Grad mes_pap), ou - la mise en oeuvre d'une combinaison de la première et de la deuxième stratégies de calcul précitées. 3. Procédé selon la revendication 2, caractérisé en ce que le choix entre les différentes mises en oeuvre de stratégies de calcul ou de combinaison des stratégies de calcul est effectué en comparant une valeur de l'horizon temporel de calcul donné (H) par rapport à une dynamique connue de déplacement dudit papillon d'arrivée des gaz. 4. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce qu'il comporte une étape de saturation du gradient de la position du papillon (Grad_pos pap) basée sur des paramètres physiques du papillon d'arrivée des gaz. 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce qu'il comporte en outre une étape de saturation de la prédiction de position du papillon d'arrivée des gaz (Ppap_pred). 6. Procédé selon la revendication 5, caractérisé en ce que l'étape de saturation de la position du papillon d'arrivé des gaz prédite (Ppap_pred) est réalisée de telle façon qu'en phase d'ouverture du papillon d'arrivée des gaz, la position prédite (Ppap pred) ne peut pas être supérieure à une consigne du papillon, et en phase de fermeture papillon, la position prédite (Ppap pred) ne peut pas être inférieure à la consigne du papillon. 7. Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce qu'il comporte une étape de filtrage de la position du papillon prédite (Ppap_pred). 8. Procédé selon la revendication 7, caractérisé en ce qu'il comporte une étape d'application d'un filtre du premier ordre.