Brevet FR3137132A1 — publié le 20231229

Le procédé de réglage du dispositif d’asservissement (MD_ABP) de l’actionneur comprend une phase d’identification d’un modèle de l’actionneur (BP) sous la forme d’une première fonction de transfert (G(s)) du second ordre et une phase de calcul d’un correcteur (C(s)). Conformément à l’invention, la phase d’identification comprend les étapes de a) commander l’actionneur en boucle fermée, avec un simple correcteur proportionnel, par un échelon et en appliquant un signal vibratoire anti-frottement, le système en boucle fermée étant modélisé sous la forme d’une deuxième fonction de transfert du second ordre ; b) déterminer des paramètres de la deuxième fonction de transfert à partir de mesures sur un signal de position sortant du système en boucle fermée en réponse à la commande appliquée ; et c) identifier la première fonction de transfert en calculant ses paramètres à partir des paramètres déterminés à l’étape b). Fig.1 PROCÉDÉ DE RÉGLAGE D’UN ASSERVISSEMENT EN POSITION D’UN ACTIONNEUR, TEL QU’UN ACTIONNEUR DE VÉHICULE AUTOMOBILE L’invention concerne de manière générale le domaine des asservissements de position des actionneurs. Plus particulièrement, l’invention se rapporte à un procédé mis en œuvre par ordinateur de réglage d’un asservissement en position d’un actionneur, tel qu’un actionneur dans un véhicule automobile. Les actionneurs commandés en position sont très utilisés dans les véhicules de type thermique pour le pilotage de leur groupe motopropulseur. Un réglage précis de la commande de ces actionneurs contribue à une optimisation du fonctionnement du groupe motopropulseur, notamment pour ce qui est du rendement de celui-ci, de l’agrément de conduite et du respect des normes antipollution. A titre d'exemples d'actionneurs, on peut citer des actionneurs présents dans la boucle d’admission d’air d'un moteur thermique de véhicule, comme ceux dans la vanne de recirculation des gaz d'échappement, dite vanne EGR, le boîtier de papillon des gaz, la soupape de décharge d’un turbocompresseur et autres. Les actionneurs utilisés dans les véhicules sont fabriqués en grandes quantités et ont fréquemment des caractéristiques relativement dispersées. Par ailleurs, ils sont assujettis à des frottements secs importants, de l’encrassement, des variations de température et des efforts aérauliques rendant leur asservissement délicat et complexe. Dans l’état de l’art, il est connu d’utiliser une méthode dite du « régleur », par essai et erreur, pour régler le correcteur de l’asservissement de commande en position jusqu’à obtenir la réponse souhaitée du système. Cette méthode nécessite un grand nombre d’essais physiques jusqu’à obtenir la performance voulue, ainsi qu’un savoir-faire important et de l’expérience pratique. Cependant, elle ne garantit pas la robustesse des performances obtenues. Il est nécessaire de réaliser des essais supplémentaires de validation avec des systèmes dispersés et des conditions extérieures variables, comme l’altitude et la température, pour s’assurer de la robustesse des performances. Les actionneurs sont habituellement pilotés par un calculateur, avec des techniques de commande issues de l’automatique avancée, en exploitant des algorithmes de commande modernes du type retour d’état, polynomiale, adaptative ou autres. La mise en œuvre des algorithmes de commande avancés susmentionnés requiert la plupart du temps un modèle mathématique dynamique du système à commander. Une première approche utilisée par l’homme du métier est de construire un modèle dynamique dit de « connaissance » en appliquant directement les lois de la physique. Cette approche nécessite une connaissance parfaite des éléments constitutifs du système à commander. De plus, elle demande de nombreux essais physiques pour adapter les paramètres du modèle, de façon à obtenir une représentation fidèle de la réponse du système réel. Une autre approche moins coûteuse est d’identifier un modèle mathématique dit de « comportement » à partir de la connaissance uniquement des entrées et des sorties du système à commander. L’identification d’un modèle dynamique pour les actionneurs de la boucle d’admission d’air présente des complexités inhérentes comme une réponse en boucle ouverte divergente, des frottements secs prépondérants et un retard pur. Ainsi, en boucle ouverte, pour une entrée constante, le comportement divergent de l’actionneur se traduit par une sortie qui augmente continuellement jusqu’à atteindre des butées mécaniques du système. Les frottements secs empêchent l’algorithme d’identification du modèle mathématique de converger vers les bons paramètres de celui-ci. L’identification du retard pur est essentielle dans la synthèse des asservissements de position, en permettant d’augmenter considérablement les performances obtenues, ainsi que leur robustesse. Par le document EP3341604A1, il est connu un système d’asservissement d’un actionneur à une consigne de position, dans lequel est intégré un compensateur de retard pur. L’invention vise à apporter une solution à la problématique exposée ci-dessus de l’état de la technique notamment pour les actionneurs commandés en position dans la boucle d’admission d’air d’un véhicule thermique, en fournissant un procédé mis en œuvre par ordinateur de réglage d’un asservissement de position d’un actionneur, grâce à une identification aisée et automatisable d’un modèle dynamique comportemental de l’actionneur du type fonction de transfert du second ordre avec un retard pur. Selon un premier aspect, l’invention concerne un procédé mis en œuvre par ordinateur de réglage d’un dispositif d’asservissement en position d’un actionneur comprenant une phase d’identification d’un modèle mathématique de l’actionneur sous la forme d’une première fonction de transfert du second ordre avec un retard pur d et une phase de calcul d’un correcteur du dispositif d’asservissement en position en utilisant le modèle mathématique identifié. Conformément à l’invention, la phase d’identification comprend les étapes de a) commander l’actionneur dans un système en boucle fermée avec une consigne en échelon et en appliquant un signal vibratoire anti-frottement, le système en boucle fermée comprenant un simple correcteur proportionnel et étant modélisé mathématiquement sous la forme d’une deuxième fonction de transfert du second ordre avec le retard pur d ; b) déterminer des paramètres de gain statique K, d’amortissement П, de pulsation naturelle Ω n et de retard pur d de la deuxième fonction de transfert à partir de mesures sur un signal de position sortant du système en boucle fermée délivré par celui-ci en réponse à la commande avec la consigne en échelon et le signal vibratoire anti-frottement appliquée à l’étape a) ; c) identifier la première fonction de transfert avec le retard pur d mesuré à l’étape b) et des paramètres de gain statique k, d’amortissement ξ, et de pulsation naturelle ω n calculés algébriquement à partir des paramètres de gain statique K, d’amortissement П et de pulsation naturelle Ω n déterminés à l’étape b). Selon une caractéristique particulière, la modélisation mathématique du système en boucle fermée est obtenue à l’aide d’une approximation en série de Taylor du second ordre. Selon une autre caractéristique particulière, la phase d’identification comprend une étape préalable de réglage du signal vibratoire anti-frottement comprenant une mesure en boucle ouverte d’une hystérésis de l’actionneur due à des frottements secs. Selon encore une autre caractéristique particulière, l’étape préalable de réglage comprend un réglage du signal vibratoire anti-frottement avec une amplitude au moins égale à une largeur de l’hystérésis et une fréquence haute ne perturbant pas le fonctionnement de l’actionneur. L’invention concerne aussi un calculateur comportant une mémoire stockant des instructions de programme pour la mise en œuvre du procédé tel que décrit brièvement ci-dessus. Selon une forme de réalisation particulière, le calculateur est un calculateur de contrôle moteur d'un véhicule automobile. L’invention concerne aussi un ensemble comprenant un moteur thermique et un calculateur de commande, le moteur thermique intégrant au moins un actionneur commandé en position dans sa boucle d’admission d’air et le calculateur de commande étant un calculateur comme indiqué ci-dessus. L’invention concerne aussi un véhicule comprenant un ensemble comme indiqué ci-dessus. D’autres avantages et caractéristiques de la présente invention apparaîtront plus clairement à la lecture de la description détaillée ci-dessous d’un mode de réalisation particulier de l’invention, en référence aux dessins annexés, dans lesquels : La est un bloc-diagramme de principe relatif à la mise en œuvre du procédé selon l’invention dans un moteur thermique de véhicule. La est un logigramme général de différentes étapes comprises dans un mode de réalisation particulier du procédé selon l’invention. La est un diagramme de principe relatif à la mesure d’une hystérésis due à des frottements secs dans un actionneur commandé en position. La est une courbe montrant une hystérésis due à des frottements secs dans un actionneur de position. La est un bloc-diagramme d’un système en boucle fermée utilisé dans le procédé de l’invention pour l’identification d’un modèle mathématique d’un actionneur commandé en position. La montre des équations mathématiques mises en application dans le procédé de l’invention. La montre des formes d’onde d’une commande du système en boucle fermée de la par une consigne en échelon, avec un signal vibratoire anti-frottement, et d’un signal de position sortant délivré en réponse par le système en boucle fermée. La montre un agrandissement de la forme d’onde du signal de position sortant de la et des mesures faites sur ce signal de position sortant pour l’identification d’un modèle mathématique de l’actionneur. En référence à la , dans la forme de réalisation particulière décrite ici, le procédé selon l’invention est mis en œuvre dans un calculateur de commande CCM d’un véhicule, tel qu’un calculateur de contrôle moteur, ayant à charge la commande d’un moteur thermique MT d’un véhicule. Comme montré schématiquement à la , le moteur thermique MT comporte une boucle d’admission d’air ayant notamment un boîtier de papillon des gaz BP, un turbocompresseur TC et une vanne EGR ayant des actionneurs commandés en position par des systèmes d’asservissement numériques. Un système logiciel embarqué SLE est contenu dans une mémoire MEM du calculateur CCM et a pour fonction la commande générale du moteur thermique MT. Les communications de données pour la commande du moteur thermique MT, entre le système logiciel embarqué SLE et différents actionneurs et capteurs du moteur thermique MT, sont réalisées typiquement à travers un réseau de communication de données (non représenté) du véhicule, par exemple un bus de type CAN. Le système logiciel embarqué SLE comprend différents modules logiciels dont un module logiciel MD_ABP, représenté à titre d’exemple à la , implémentant un dispositif d’asservissement numérique dédié à la commande en position de l’actionneur du boîtier de papillon des gaz BP. Conformément à l’invention, le système logiciel embarqué SLE comprend également un module logiciel MD_R coopérant avec le dispositif d’asservissement numérique MD_ABP. Le module logiciel MD_R est chargé de la mise en œuvre du procédé selon l’invention par l’exécution d’instructions de code de programme par un processeur (non représenté) du calculateur CCM. Le dispositif d’asservissement numérique MD_ABP est configurable dans plusieurs modes, sous la commande du module logiciel MD_R chargé de la mise en œuvre du procédé de l’invention, à savoir une configuration fonctionnelle et des configurations de réglage. A la , le dispositif d’asservissement numérique MD_ABP est montré schématiquement dans une configuration fonctionnelle, en boucle fermée. Dans cette configuration fonctionnelle, le modèle mathématique, dynamique et comportemental, de l’actionneur de position du boîtier de papillon des gaz BP a été identifié précédemment, comme cela sera décrit plus bas, au moyen du procédé de l’invention. La connaissance du modèle de l’actionneur du boîtier BP, représenté par la fonction de transfert G(s) avec s étant l’opérateur de Laplace, a permis de définir un correcteur C(s) typiquement de type dit « PID » (pour « Proportionnel-Intégral-Dérivé »). La définition du correcteur C(s), une fois identifiée la fonction de transfert G(s), est réalisée par des méthodes de calcul connues de l’homme du métier, qui ne seront pas décrites ici. Typiquement, une fonction de précommande P_Cde est intégrée également dans le dispositif d’asservissement numérique MD_ABP. De plus, une fonction anti-frottement (non représentée) pourra être intégrée dans le dispositif d’asservissement numérique MD_ABP, sous la forme d’un signal vibratoire de fréquence et amplitude prédéfinies. Comme visible à la , le dispositif d’asservissement numérique MD_ABP reçoit en entrée une consigne de position, désignée R(s) dans l’espace de Laplace, et délivre en sortie une commande de position, désignée U(s) et U respectivement dans l’espace de Laplace et l’espace temporel. En contre-réaction, le dispositif d’asservissement numérique MD_ABP reçoit également en entrée une information de retour de position Y délivrée par le boîtier de papillon des gaz BP, désignée Y(s) dans l’espace de Laplace, et représentant la position effective du papillon des gaz dans le boîtier BP suite à la consigne de position R(s). Comme visible à la , la fonction de transfert G(s) de l’actionneur du boîtier BP, à identifier par le procédé de l’invention, est une fonction de transfert du second ordre avec un retard pur, dans laquelle les paramètres k, ξ, ω n et d représentent respectivement le gain statique, l’amortissement, la pulsation naturelle et le retard pur de l’actionneur susmentionné en boucle ouverte. En référence aussi maintenant aux Figs.2 à 8, il est décrit en détail un mode de réalisation particulier du procédé de l’invention, dans le cadre traité ici du réglage du dispositif d’asservissement numérique MD_ABP assurant la commande en position de l’actionneur du boîtier de papillon des gaz BP. Comme visible à la , le procédé selon l’invention comprend essentiellement trois grandes étapes S1 à S3 dont l’exécution est gérée par le module logiciel MD_R. Le processus de réglage du dispositif d’asservissement numérique MD_ABP est lancé typiquement par un opérateur humain, chargé de la mise au point du dispositif, qui interagit avec le module logiciel MD_R via une interface homme-machine telle qu’un ordinateur ORD, ou un banc d’essai, hébergeant une application logicielle dédiée. Les étapes S1 à S3 correspondent respectivement à A) une première étape de test en boucle ouverte de l’actionneur du boîtier BP, avec le dispositif d’asservissement numérique MD_ABP dans une première configuration de réglage, pour la définition d’un signal vibratoire anti-frottement utilisé dans le procédé de l’invention, B) une deuxième étape de test en boucle fermé de l’actionneur du boîtier BP, avec le dispositif d’asservissement numérique MD_ABP dans une deuxième configuration de réglage, pour identifier les paramètres k, ξ, ω n et d à partir de mesures sur une réponse à un échelon avec signal anti-frottement et de calculs algébriques, et C) une troisième étape de définition de la configuration fonctionnelle du dispositif d’asservissement numérique MD_ABP, connaissant les paramètres k, ξ, ω n et d. En référence plus particulièrement aux Figs.3 et 4, à l’étape S1, l’actionneur du boîtier BP est commandé directement en boucle ouverte, sous le contrôle du processus de traitement du module logiciel MD_R, par une consigne en double pente Rdp délivrée par le dispositif d’asservissement numérique MD_ABP dans la première configuration de réglage susmentionnée. La réponse de position effective Ydp de l’actionneur du boîtier BP, pour la consigne en double pente Rdp émise, est lue par le processus du module logiciel MD_R. La consigne en double pente Rdp et la réponse de position effective Ydp de l’actionneur du boîtier BP sont analysées par le processus du module logiciel MD_R, de façon à identifier une relation, sous la forme d’une hystérésis, représentative des frottements secs s’opposant au mouvement fluide de l’actionneur du boîtier BP. La largeur de cette hystérésis permet de dimensionner l’amplitude que doit avoir le signal vibratoire anti-frottement pour vaincre les frottements secs. La fréquence du signal vibratoire anti-frottement à appliquer est extraite typiquement d’une base de connaissance, en fonction du type d’actionneur, cette fréquence devra être suffisamment haute pour ne pas perturber le fonctionnement de l’actionneur. Un exemple de courbe d’hystérésis obtenue, Hdp, est montrée à la . Les amplitudes en pourcents (%) de la consigne en double pente Rdp et de la réponse de position effective Ydp sont indiquées en abscisses et en ordonnées, respectivement, de la courbe d’hystérésis Hdp. La largeur LH de l’hystérésis est ici égale à environ 20% de l’amplitude de la consigne Rdp. Comme visible à la , l’étape S2 comprend essentiellement trois sous-étapes S20 à S22. A la sous-étape S20, le processus du module logiciel MD_R commande la mise du dispositif d’asservissement numérique MD_ABP dans la deuxième configuration de réglage susmentionnée, désignée R_ABP à la . Dans cette configuration de réglage, l’actionneur du boîtier BP est commandé en boucle fermée avec un correcteur proportionnel k p et le signal vibratoire anti-frottement identifié à l’étape S1 et désigné ici SAF. La commande en boucle fermée avec le correcteur proportionnel k p vise à apporter une solution au comportement divergent en boucle ouverte de l’actionneur du boîtier BP, en forçant celui-ci à converger. Le signal vibratoire anti-frottement SAF est ajouté dans la commande de position Ue(s) appliquée à l’actionneur du boîtier BP. Le signal vibratoire anti-frottement SAF s’oppose à l’effet des frottements secs, en évitant un arrêt de l’actionneur à une position qui ne correspondrait pas à la commande, ce qui introduirait une erreur dans l’identification du gain statique de la fonction de transfert G(s). Dans cette configuration de réglage de la , le système en boucle fermée, formé de l’actionneur du boîtier BP et du dispositif R_ABP, a pour fonction de transfert la fonction H(s) représentée par l’équation E1 à la . Dans la fonction de transfert H(s) =Ye(s)/Re(s), Re(s) et Ye(s) sont respectivement l’entrée et la sortie du système bouclée, à savoir, la consigne de position et la réponse en position de l’actionneur à cette consigne, et les paramètres k, ξ, ω n et d sont, comme susmentionné, le gain statique, l’amortissement, la pulsation naturelle et le retard pur à identifier de l’actionneur en boucle ouverte. La fonction de transfert H(s) reste une fonction du second ordre grâce à l’utilisation d’un correcteur proportionnel k p . Cependant, cette fonction de transfert H(s) n’est pas rationnelle du fait de la présence du terme e -ds dans le dénominateur de la fonction H(s). Pour pouvoir réaliser des calculs algébriques, il est nécessaire d’éliminer le terme e -ds dans le dénominateur de la fonction H(s) en le remplaçant par une fonction rationnelle. Le choix est fait dans l’invention d’utiliser l’approximation en série de Taylor du second ordre, donnée par l’équation E2 à la , pour remplacer le terme e -ds . Des simulations comparatives, réalisées par l’entité inventive et basées sur le diagramme de Bode (module et phase), entre la fonction de transfert originale et la fonction de transfert approximée, ont montrées que l’erreur due à cette approximation est négligeable dans la bande passante du système bouclée. Le choix d’approximer le terme e -ds par la série de Taylor du second ordre permet de conserver le second ordre pour la fonction de transfert H(s) approximée, ce qui ne serait pas le cas avec l’approximation de Padé du premier ordre communément utilisé pour le terme e -ds . La fonction de transfert H(s) obtenue avec l’approximation en série de Taylor est donnée par l’équation E3 à la . Cette équation E3 est du second ordre et comprend des paramètres de gain statique K, d’amortissement П, de pulsation naturelle Ω n et de retard pur d. Par identification terme à terme des équations E1 et E3 de la fonction de transfert H(s), il est obtenu les équations E4 à la qui lient les paramètres K, П et Ω n aux paramètres k, ξ et ω n à identifier. Par inversion des équations E4, il est obtenu les égalités E5 à la qui donnent les paramètres à identifier k, ξ et ω n en fonction des paramètres K, П, Ω n , d et k p . Les sous-étapes S21 et S22 sont respectivement une sous-étape d’obtention par mesure et calcul des paramètres K, П, Ω n et d et une sous-étape de calcul des paramètres à identifier k, ξ et ω n en fonction des valeurs déterminées de K, П, Ω n , d et du gain k p connu. En référence à la et aux Figs.7 et 8, les paramètres K, П, Ω n et d sont obtenus à la sous-étape S21 à partir de mesures effectuées sur un signal de position sortant Ye (position de l’actionneur) délivré par le système bouclé configuré à la sous-étape S20 (R_ABP, ) en réponse à une consigne en échelon Re appliquée en entrée du système bouclé, avec le signal vibratoire anti-frottement SAF présent dans la commande Ue de l’actionneur. A la sont montrés des exemples de formes d’onde de la consigne en échelon Re, du signal de position sortant Ye, ainsi que de la commande Ue, SAF, et de la commande moyennée correspondante Ue_moy. La position P en pourcents (%) est indiquée dans l’axe des ordonnées de ces formes d’onde. Comme visible à la , sur le signal de position sortant Ye sont mesurés le retard d, un temps de pic tp, un dépassement de pic Dp et un delta sortant Ds correspondant à un delta entrant De de la consigne en échelon Re. Le temps de pic tp est le temps de montée du signal de position sortant Ye jusqu’à atteindre le dépassement de pic Dp mesuré par rapport au delta sortant Ds. Le retard d à identifier est donc mesuré directement sur le signal de position sortant Ye. Les paramètres K, П et Ω n sont calculés avec les égalités connues E6) montrées à la , à partir des valeurs de tp, Dp, Ds et De. A la sous-étape S22, les paramètres k, ξ et ω n de la fonction de transfert G(s) sont calculés avec les égalités E5 montrées à la , à partir des valeurs de K, П, Ω n et d obtenues à la sous-étape S21. A l’issue de la sous-étape S22, les différents paramètres de la fonction de transfert G(s) sont identifiés et la modélisation mathématique de l’actionneur du boîtier BP est donc achevée. A l’étape S3, le processus du module logiciel MD_R calcule, avec la fonction de transfert G(s) identifiée à l’étape S2, un correcteur C(s) adapté pour le dispositif d’asservissement numérique MD_ABP dans sa configuration fonctionnelle montrée à la . Les méthodes de calcul des correcteurs sont connues de l’homme du métier et ne sont donc pas détaillées ici. Le procédé de l’invention trouve une application privilégiée, mais non exclusive, dans la commande en position des actionneurs des moteurs thermique de type essence et Diesel. La présente invention permet d’améliorer la robustesse et les performances des asservissements de position des actionneurs, ainsi que leurs temps de mise en œuvre et de validation. Le procédé de l’invention est totalement automatisable et le processus de réglage qui en résulte peut-être mis en pratique aisément par une personne peu expérimentée. L’invention ne se limite pas au mode de réalisation particulier qui a été décrit ici à titre d’exemple. L’homme du métier, selon les applications de l’invention, pourra apporter différentes modifications et variantes entrant dans le champ de protection de l’invention. Procédé mis en œuvre par ordinateur de réglage d’un dispositif d’asservissement en position (MD_ABP) d’un actionneur (BP) comprenant une phase d’identification (S1, S2) d’un modèle mathématique dudit actionneur (BP) sous la forme d’une première fonction de transfert (G(s)) du second ordre avec un retard pur d et une phase de calcul (S3) d’un correcteur (C(s)) dudit dispositif d’asservissement en position (MD_ABP) en utilisant ledit modèle mathématique (G(s)), caractérisé en ce que ladite phase d’identification (S1, S2) comprend les étapes de a) commander (S21) ledit actionneur (BP) dans un système en boucle fermée (R_ABP, BP) avec une consigne en échelon (Re) et en appliquant un signal vibratoire anti-frottement (SAF), ledit système en boucle fermée (R_ABP, BP) comprenant un simple correcteur proportionnel (k p ) et étant modélisé mathématiquement sous la forme d’une deuxième fonction de transfert du second ordre avec ledit retard pur d (H(s)) ; b) déterminer (S22, E6) des paramètres de gain statique K, d’amortissement П, de pulsation naturelle Ω n et de retard pur d de ladite deuxième fonction de transfert (H(s)) à partir de mesures (De, Ds, d, tp, Dp) sur un signal de position sortant (Ye) dudit système en boucle fermée (R_ABP, BP) délivré par celui-ci en réponse à la commande avec ladite consigne en échelon (Re) et ledit signal vibratoire anti-frottement (SAF) appliquée à l’étape a) ; et c) identifier (S23, E5) ladite première fonction de transfert (G(s)) avec ledit retard pur d mesuré à l’étape b) et des paramètres de gain statique k, d’amortissement ξ, et de pulsation naturelle ω n calculés algébriquement (E5) à partir des paramètres de gain statique K, d’amortissement П et de pulsation naturelle Ω n déterminés à l’étape b). Procédé selon la revendication 1, caractérisé en ce que la modélisation mathématique dudit système en boucle fermée (R_ABP, BP) est obtenue à l’aide d’une approximation en série de Taylor du second ordre (E2). Procédé selon la revendication 1 ou 2, caractérisé en ce que ladite phase d’identification comprend une étape préalable de réglage (S1) dudit signal vibratoire anti-frottement (SAF) comprenant une mesure en boucle ouverte (Rdp, Ydp) d’une hystérésis (Hdp) dudit actionneur (BP) due à des frottements secs. Procédé selon la revendication 3, caractérisé en ce que ladite étape préalable de réglage (S1) comprend un réglage dudit signal vibratoire anti-frottement (SAF) avec une amplitude au moins égale à une largeur (LH) de ladite hystérésis (Hdp) et une fréquence haute ne perturbant pas le fonctionnement dudit actionneur (BP). Calculateur (CCM) comportant une mémoire (MEM) stockant des instructions de programme (MD_R) pour la mise en œuvre du procédé selon l’une quelconque des revendications 1 à 4. Calculateur selon la revendication 5, caractérisé en ce que ledit calculateur est un calculateur de contrôle moteur (CCM) d'un véhicule. Ensemble comprenant un moteur thermique (MT) et un calculateur de commande (CCM), ledit moteur thermique (MT) intégrant au moins un actionneur (BP, EGR, TC) commandé en position dans sa boucle d’admission d’air, caractérisé en ce que ledit calculateur de commande est un calculateur (CCM) selon la revendication 5 ou 6. Véhicule caractérisé en ce qu’il comprend un ensemble (MT, CCM) selon la revendication 7.