FR 3062884 A1 20180817 FR 1751206 A 20170215 1 PROCEDE DE CONTROLE PERMETTANT LE CALCUL D'UNE SENSIBILITE D'UN MODELE DE COMPORTEMENT D'UN ACTIONNEUR DE MOTEUR THERMIQUE DE VEHICULE AUTOMOBILE [0001] La présente invention porte sur un procédé de contrôle permettant le calcul d'une sensibilité d'un modèle de comportement d'un actionneur de moteur thermique de véhicule automobile. [0002] Du fait des dispersions de fabrication, de l'usure, et de l'encrassement des actionneurs du moteur thermique, leur comportement physique peut différer des modèles de comportement intégrés dans le calculateur moteur. Ce décalage des modèles d'actionneurs, notamment ceux relatifs à la branche d'admission et à la branche d'injection, amène à des dérives de richesse, donc à une surconsommation ou à une augmentation des émissions polluantes, et à des impacts sur l'agrément de conduite ressenti par le conducteur. [0003] Le calculateur moteur doit donc, tout au long de la vie du véhicule, exploiter les informations des différents capteurs présents sur le système pour recaler les modèles de comportement et les adapter au mieux au fonctionnement réel des actionneurs. Pour cela, le calculateur moteur doit pouvoir établir, sur le point de fonctionnement courant, le lien entre les grandeurs mesurées par les capteurs et les paramètres des modèles de comportement à recaler. Ce lien est donné par les sensibilités, c'est-à-dire pour un modèle, par le rapport entre la variation d'une sortie du modèle et la variation d'une de ses entrées (ou paramètre interne). [0004] Les sensibilités peuvent être déterminées en fonction de leur complexité par une première méthode de dérivée mathématique des équations d'état du système, une deuxième méthode de calcul d'une variation locale autour du point courant, ou une troisième méthode d'identification des sensibilités du moteur lors de la phase de mise au point. [0005] En utilisant la première méthode, on peut déterminer par exemple la sensibilité du modèle de comportement de la vanne papillon liant un écart de débit (débit mesuré par le capteur de débit moins le débit estimé par le modèle de comportement de la vanne papillon) à la variation de la section de la vanne papillon. On part de l'équation connue qui modélise le comportement de la vanne papillon: 3062884 2 PO Qpapillon = Section papillonx ,x f (-P ) -N/ TO PO avec: - Qpapillon étant le débit de la vanne papillon; - f(-)étant une modélisation connue de la perte de charge au niveau de la vanne PO 5 papillon; - PO étant la pression en amont du volet de la vanne papillon; - P étant la pression en aval du volet de la vanne papillon; - TO étant la température en amont du volet de la vanne papillon; - Section papillon étant la section de passage découverte par l'ouverture du volet de la 10 vanne papillon. [0006] En dérivant cette équation par rapport à la section de la vanne papillon, on obtient la sensibilité suivante : PO = x Cfe(-) eection papillon .NITO PO [0007] Cette équation permet, pour chaque point de fonctionnement, de déterminer le 15 lien entre une variation de section de la vanne papillon M'ection papillon et la variation de débit Qpapillon résultante. [0008] La deuxième méthode peut être appliquée par exemple dans le cas du calcul de la sensibilité de la modélisation du débit d'un cylindre à la position du déphaseur d'arbre à cames d'admission. Cette deuxième méthode consiste à comparer, sur le point courant, la 20 valeur de débit du cylindre Qcylindre issue de la modélisation, et celle de ce même modèle calculée sur le même point courant mais en rajoutant un biais sur l'entrée pour laquelle on souhaite connaitre la sensibilité, noté Qcylindre biais Vi7TA dans l'équation ci-dessous. Dans le cas de l'exemple, on introduit un biais biais WTA sur la position du déphaseur d'arbre à cames d'admission. On obtient la sensibilité de l'écart de débit du cylindre 25 &cylindre à la variation de la position du déphaseur d'arbre à cames d'admission &/VTA par la formule suivante : (Qcylindre Qcylindre - Qcylindre biais VVTA SVVTA -biais VVTA 3062884 3 [0009] La troisième méthode consiste, au moment de la phase de mise au point du moteur thermique, à réaliser des essais spécifiques sur un moteur thermique placé au banc d'essai. Les essais consistent à réaliser des variations des différents paramètres sur chacun des points de fonctionnement sur lesquels se déroule la mise au point, et d'en 5 déduire la sensibilité réelle du moteur thermique aux différents paramètres. Ces sensibilités doivent être ensuite mémorisées par point de fonctionnement (ou zone de fonctionnement) dans le calculateur moteur. [0010] Toutefois, l'application de la première méthode peut s'avérer relativement lourde et difficile à mettre en oeuvre. En effet, le calcul de la dérivée mathématique de la fonction 10 de transfert du modèle pour lequel on cherche la sensibilité n'est pas toujours aisé. Ainsi, la dérivée mathématique de l'équation du débit du cylindre par rapport à la position du déphaseur d'arbre à cames d'admission est relativement complexe puisque cette position intervient à de multiples endroits dans la modélisation du débit du cylindre. Dans d'autres cas, la détermination de la dérivée mathématique est même impossible. C'est le cas par 15 exemple si la fonction de transfert n'est pas modélisée par une équation dérivable. [0011] L'application de la deuxième méthode implique généralement une duplication du modèle dont on souhaite connaître la sensibilité. On a, en effet, besoin une première fois du modèle pour déterminer la valeur courante de la sortie et une seconde fois pour obtenir la valeur courante de la sortie modifiée par le biais sur l'entrée pour chaque sensibilité à 20 calculer. Suivant la taille du modèle dont on souhaite identifier les sensibilités, cela peut s'avérer coûteux en taille de mémoire occupée dans le calculateur mais également alourdir la charge du processeur par la multiplication des calculs. [0012] L'application de la troisième méthode est également relativement lourde puisqu'elle nécessite des essais spécifiques au banc moteur thermique en phase de mise 25 au point. De plus, la troisième méthode nécessite un système complexe de mémorisation des sensibilités des modèles d'actionneurs dans le calculateur moteur. [0013] L'invention vise à remédier efficacement à ces inconvénients en proposant un procédé de contrôle permettant le calcul d'une sensibilité d'au moins un modèle de comportement d'un actionneur de moteur thermique d'un véhicule automobile ayant au 30 moins une sortie déterminée en fonction d'au moins deux entrées, ledit modèle étant implémenté dans un calculateur moteur, ledit procédé comportant: - une étape de calcul de la sortie dudit modèle sans biais et de stockage de ladite valeur de sortie comme valeur de référence, 3062884 4 - une étape de calcul de la sortie dudit modèle avec un biais appliqué sur une première entrée de manière à déterminer une sensibilité de ladite sortie du modèle à la variation de ladite première entrée, au moins une fonction du calculateur moteur utilisant ledit modèle étant alimentée par ladite valeur de référence de ladite sortie déterminée pendant l'étape 5 précédente, - une étape de calcul de la sortie dudit modèle avec un biais appliqué sur une deuxième entrée de manière à déterminer la sensibilité de ladite sortie dudit modèle à la variation de ladite deuxième entrée. [0014] L'invention permet ainsi, du fait du calcul de la sensibilité pendant l'exploitation du 10 modèle par le calculateur moteur, d'éviter la duplication du ou des modèles dont on souhaite déterminer la sensibilité. L'invention consomme donc peu de ressources mémoires et a donc peu d'impact sur la charge du calculateur moteur. [0015] Selon une mise en oeuvre, le procédé comporte de plus une étape de commande dudit actionneur de moteur thermique en tenant compte desdites sensibilités dudit modèle 15 préalablement déterminées. [0016] Selon une mise en oeuvre, entre les deux étapes de calcul de la sortie dudit modèle avec un biais appliqué sur une de ses entrées, ledit procédé comporte: - une étape de mise à jour de ladite valeur de référence de ladite sortie dudit modèle sans biais en retirant le biais sur ladite première entrée, 20 - la ou les fonctions dudit calculateur moteur utilisant ledit modèle étant alimentées par ladite valeur de référence de la sortie mise à jour lors de la mise en oeuvre de l'étape de calcul de la sortie du modèle avec le biais appliqué sur la deuxième entrée. [0017] Selon une mise en oeuvre, chaque étape de calcul est mise en oeuvre lors d'un appel de tâche distinct correspondant du modèle de comportement de l'actionneur de 25 moteur thermique. [0018] Selon une mise en oeuvre, dans le cas où le modèle de l'actionneur de moteur thermique nécessite une durée de convergence, les étapes de calcul sont mises en oeuvre sur un laps de temps défini pour ledit modèle. [0019] Selon une mise en oeuvre, ledit laps de temps est variable en fonction d'un point 30 de fonctionnement courant, par exemple, en fonction d'un régime moteur thermique. 3062884 5 [0020] Selon une mise en oeuvre, des valeurs des biais appliqués sur les entrées du modèle sont fixes. [0021] Selon une mise en oeuvre, des valeurs des biais appliqués sur les entrées du modèle sont variables dynamiquement. 5 [0022] Selon une mise en oeuvre, une limite minimale des biais est telle qu'une introduction desdits biais génère une réponse dudit modèle. [0023] L'invention a également pour objet un calculateur moteur pour véhicule automobile comportant une mémoire stockant des instructions logicielles pour la mise en oeuvre du procédé de contrôle d'une sensibilité d'au moins un modèle de comportement 10 d'un actionneur de moteur thermique du véhicule automobile tel que précédemment défini. [0024] 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'à titre illustratif mais nullement limitatif de l'invention. [0025] La figure 1 est une représentation graphique illustrant une sensibilité d'un débit 15 d'un cylindre de moteur thermique à une variation d'un déphaseur d'arbre à cames d'admission; [0026] La figure 2 est une représentation schématique d'un calculateur moteur comportant une mémoire stockant des instructions logicielles pour la mise en oeuvre du procédé selon la présente invention; 20 [0027] La figure 3 un diagramme fonctionnel des différentes étapes du procédé selon la présente invention. [0028] L'invention est un procédé de contrôle permettant le calcul, sur un point de fonctionnement courant, d'une sensibilité d'une sortie Y1, Y2 d'un modèle de comportement H(x) d'un actionneur de moteur thermique d'un véhicule automobile à une 25 de ses entrées X1, X2, X3, X4, tel que cela est montré sur la figure 1. Le procédé est basé sur l'analyse de l'impact sur la sortie Y1, Y2 lorsque l'on crée une variation locale d'une entrée X1, X2, X3, X4 sur le modèle H(x) considéré. [0029] Par exemple, on peut appliquer ce procédé dans le cas du calcul de la sensibilité de la modélisation du débit d'un cylindre à la position du déphaseur d'arbre à cames 30 d'admission. Cela consiste à comparer sur le point courant, la valeur de débit du cylindre 3062884 6 Qcylindre issue de la modélisation et celle de ce même modèle calculée sur le même point courant mais en rajoutant un biais sur l'entrée pour laquelle on souhaite connaitre la sensibilité, notée Qcylindre biais WTA. [0030] Dans le cas de l'exemple, on introduit un biais sur la position du déphaseur 5 d'arbre à cames d'admission biais VV7'A. Par la formule suivante, on obtient la sensibilité de l'écart de débit du cylindre 6Qcylindre à la variation de la position du déphaseur d'arbre à cames d'admission 8WTA: (Qcylindre Qcylindre- Qcylindre biais WTA SVVTA -biais WTA [0031] Un exemple de représentation graphique de la sensibilité du débit du cylindre 10 Ev Qcyl à une variation de la position du déphaseur d'arbre à cames d'admission Ev VVTA est représenté par la courbe C sur la figure 1. Dans cet exemple, la valeur de débit du cylindre Qcylindre est déjà présente et disponible dans le calculateur moteur 1. Elle n'est pas spécifique au calcul de la sensibilité puisqu'elle est nécessaire à l'estimation de la masse d'air enfermée dans le cylindre qui va servir notamment à gérer la commande 15 d'ouverture de l'injecteur. Par contre, la valeur de débit du cylindre impactée par un biais sur le déphaseur d'arbre à cames d'admission Qcylindre biais WTA est spécifique au besoin de calcul de la sensibilité. On note que sur la figure 1, VVTA_pc correspond à la valeur de la position du déphaseur d'arbre à cames d'admission au point de fonctionnement courant, et que VVTA biais correspond à la valeur de la position du 20 déphaseur d'arbre à cames d'admission au point de fonctionnement courant additionnée du biais. [0032] L'invention propose un procédé pour déterminer la sortie du modèle impactée par un biais sur une des entrées sans avoir à dupliquer le calcul du modèle H(x). Cela minimise les impacts sur la charge du calculateur moteur 1 et l'occupation de l'espace de 25 la mémoire 2 du calculateur moteur 1. [0033] Le procédé selon l'invention consiste à utiliser le modèle H(x) déjà présent dans le calculateur moteur 1 et d'appliquer ponctuellement sur un pas de calcul un biais sur une des entrées X1, X2, X3, X4. L'analyse de la sortie Y1, Y2 du modèle H(x) avec le biais comparé à cette même sortie Y1, Y2 sans le biais (retenue du pas de calcul précédent) 30 permet de déterminer la sensibilité du modèle H(x) à cette entrée X1, X2, X3, X4. 3062884 7 [0034] Pour le modèle donné H(x) comportant les quatre entrées X1, X2, X3, X4 et les deux sorties Y1, Y2, dès qu'un besoin d'estimer les sensibilités est émis, le procédé détermine la sensibilité de la sortie Y1 aux variations des entrées X1 puis X2 à travers le déroulement de plusieurs étapes comme cela est représenté sur la figure 3. 5 [0035] Lors d'une première étape 100, le calculateur moteur 1 calcule et stocke les valeurs de référence des sorties Y1, Y2 du modèle H(x), c'est-à-dire la valeur sur le point de fonctionnement courant des sorties Y1, Y2. Le calculateur moteur 1 les stocke respectivement dans les variables Yrefl et Y ref2. Cette étape 100 est réalisée sur l'appel de tâche courant, noté n, du modèle H(x) au moment où le besoin d'estimer les sensibilités 10 est émis. A la fin de cette étape 100, le calculateur moteur 1 fige à leur valeur courante les entrées Xl, X2, X3, X4 du modèle H(x). [0036] Une deuxième étape 101 consiste à calculer la sortie Y1 du modèle H(x) avec un biais appliqué sur la première entrée X1 lors d'un deuxième appel de tâche du modèle H(x), c'est-à-dire l'appel de tâche n+1. Ainsi, par exemple, si le modèle H(x) est appelé et 15 mis à jour toutes les 10ms, la première étape 100 est réalisée au temps zéro et la deuxième étape 101 est réalisée au temps t=0+10ms. [0037] Lors de cette étape 101, on applique un biais additif, noté biais Xl, sur l'entrée X1 du modèle H(x). Les autres entrées X2, X3, X4 du modèle H(x) ne sont pas modifiées et sont conservées figées comme demandé lors de la première étape 100. La sortie Y1 du 20 modèle H(x) pour cette configuration est notée Y biaisl. Durant cette étape 101, les autres fonctions Fonc du calculateur moteur 1 qui consomment les sorties Y1, Y2 du modèle H(x) sont alimentées avec les valeurs Yrefl et Y ref2 calculées lors de la première étape 100. Cela permet d'éviter de perturber le système avec l'utilisation de sorties modifiées par l'application du biais. 25 [0038] A partir de Y biais1 (qui est la sortie du modèle H(x) avec application d'un biais additif biais X1 sur l'entrée X1), de Y ref1 (sortie du modèle H(x) sans biais), et de la valeur du biais biais X1 appliquée sur l'entrée X1, on calcule d'après la formule suivante, la sensibilité de la sortie Y1 à la variation de l'entrée X1 du modèle H(x) : 6171 Y refl-Y biaisl -biais X1 3062884 8 [0039] Une troisième étape 102, réalisée sur l'appel de tâche suivant celui de la deuxième étape 101, c'est-à-dire sur l'appel de tâche n+2, consiste à mettre à jour le stockage des valeurs de référence des sorties du modèle Y1, Y2. A cet effet, le calculateur moteur 1 retire le biais additif biais X/ sur l'entrée X1, libère le figeage des entrées X2, X3, 5 X4, et remplace Yrefl et Yref2 par les sorties courantes Y1, Y2 du modèle H(x). Cette étape 102 permet de prendre en compte le fait que le système peut avoir évolué en terme de point de fonctionnement depuis la première étape 100 et qu'il est nécessaire de mettre à jour les nouvelles valeurs de Y1, Y2, X1, X2, X3 et X4. A la fin de cette étape 102, le calculateur moteur 1 fige à leur valeur courante les entrées X1, X2, X3, X4 du modèle 10 H(x). [0040] Au cours d'une quatrième étape 103 se déroulant lors de l'appel de tâche du modèle H(x) suivant celui de la troisième étape 102, c'est-à-dire l'appel de tâche n+3, le calculateur moteur 1 applique uniquement un biais additif biais X2 sur l'entrée X2 du modèle H(x). Les autres entrées X1, X3, X4 du modèle H(x) ne sont pas modifiées et sont 15 conservées figées. La sortie du modèle H(x) est noté Y biais2. Durant cette étape 103, les autres fonctions Fonc du calculateur moteur 1 qui consomment les sorties Y1, Y2 du modèle H(x) sont alimentées avec les valeurs de Yrefl et Yref2 calculées lors de la troisième étape 102. Cela permet d'éviter de perturber le système avec l'utilisation de sortie modifiée par l'application du biais. 20 [0041] A partir de Y biais2 (correspondant à la sortie du modèle H(x) avec application d'un biais additif (biais X2) sur l'entrée X2), de Yref2 (sortie du modèle H(x) sans biais), et de la valeur du biais biais X2 appliqué sur l'entrée X2, le calculateur moteur 1 calcule d'après la formule suivante, la sensibilité de la sortie Y1 à la variation de l'entrée X2 du modèle H(x) : 25 6Y1 Y refl-Y biais2 5X2 -biais X 2 [0042] Au terme des quatre étapes 100, 101, 102, 103 précédentes, le procédé fournit les sensibilités de la sortie Y1 du modèle H(x) aux variations des entrées X1 et X2 sur le point courant. Le calculateur moteur 1 retire ensuite les biais additifs sur les différentes entrées (biais X1 et biais X2), libère le figeage des entrées X1, X2, X3, X4, et alimente les 30 autres fonctions Fonc par les sorties courantes Y1, Y2 du modèle H(x). L'actionneur du moteur est ensuite commandé dans une cinquième étape 104 en tenant compte des sensibilités du modèle H(x) ainsi déterminées. 3062884 9 [0043] Le processus pourra être relancé depuis la première étape 100 si un nouveau besoin d'estimer les sensibilités est émis par le calculateur moteur 1. [0044] En fonction des besoins en termes de calcul des sensibilités, on peut rallonger le procédé en rajoutant des étapes pour calculer, par exemple, la sensibilité de la sortie Y1 5 du modèle H(x) aux variations de l'entrée X3 ou bien la sensibilité de la sortie Y2 du modèle H(x) aux variations de l'entrée X1. Dans ce cas, on rajoute alternativement, une étape de mise à jour, identique à la troisième étape 102, puis une étape d'introduction d'un biais sur une des entrées X1, X2, X3, X4, comme la deuxième étape 101 ou la quatrième étape 103. Le processus de calcul de l'ensemble des sensibilités est alors plus long 10 puisque chaque étape doit être faite successivement sur les appels de tâche du modèle considéré. [0045] Une variante du procédé pourra être de considérer l'étape 102 de mise à jour des sorties Y1, Y2 du modèle H(x) comme optionnelle. Le calculateur moteur 1 conserve alors la valeur de référence des sorties calculées lors de la première étape 100 et enchaîne à la 15 suite sur chaque appel de tâche l'ajout et le retrait d'un biais sur une des entrées X1, X2, X3, X4. Cette variante réduit le temps total d'exécution du procédé, notamment lorsque l'on doit estimer un grand nombre de sensibilité d'un même modèle, mais rajoute un risque vis-à-vis des autres procédés du calculateur moteur 1 consommant les sorties Y1, Y2 du modèle considéré. En effet, les sorties Y1, Y2 n'étant plus mises à jour aussi 20 régulièrement, il y a un risque, notamment lors d'un transitoire de charge, que la valeur figée ne soit plus suffisamment représentative de ce qui se passe réellement sur le moteur thermique. [0046] Par exemple, si le modèle considéré H(x) est le modèle de débit du cylindre dont une des sorties Y1, Y2 est l'estimation de la masse d'air enfermée dans le cylindre 25 permettant de déterminer le temps d'ouverture de l'injecteur, le calculateur moteur 1 risque de ne plus commander la bonne quantité de carburant et de ne plus respecter la consigne de richesse si cette sortie Y1, Y2 du modèle H(x) n'est pas réactualisée alors qu'un transitoire de charge survient. Ceci risque d'avoir un impact sur la consommation et les émissions polluantes. Cette variante est donc à privilégier dans le cas où le procédé peut 30 être exécuté sur des points stabilisés. Par point stabilisé, on entend un point de fonctionnement présentant des entrées et des sorties du modèle considéré H(x) comme relativement invariantes dans un délai donné. En cas d'interruption du critère de stabilité 3062884 10 du point, le principe de figeage des variables et d'application des biais doit être interrompu instantanément pour permettre de retrouver un calcul classique du modèle considéré H(x). [0047] Par contre, cette variante peut être nécessaire s'il y a un besoin impératif d'avoir les sensibilités du modèle H(x) calculées sur un même point de fonctionnement, par 5 exemple pour alimenter une stratégie d'adaptation répartissant une erreur mesurée sur le moteur thermique à travers plusieurs sources en utilisant les sensibilités pour définir cette répartition. Il est alors important de ne pas réaliser l'étape 102 de mise à jour pour s'assurer que les calculs de chaque sensibilité sont effectués avec exactement les mêmes entrées X1, X2, X3, X4. 10 [0048] En fonction du modèle H(x) considéré pour le calcul des sensibilités, il est possible que ce dernier nécessite un temps de convergence, c'est-à-dire que lorsque l'on modifie une de ses entrées X1, X2, X3, X4, il faut un certain temps pour que sa sortie Y1, Y2 ait le temps d'évoluer vers sa valeur finale prenant en compte l'impact de l'entrée X1, X2, X3, X4 modifiée. C'est notamment le cas de modèle comportant des éléments de 15 filtrage ou de bouclage de sorte qu'une sortie Y1, Y2 du modèle impacte directement une des entrées X1, X2, X3, X4 par exemple. Pour ce type de modèle, il est possible d'adopter une autre variante de la stratégie proposée. [0049] Cette variante reprend le procédé et les étapes 100, 101, 102, 103, 104 présentées précédemment, mais au lieu d'enchaîner les étapes 100, 101, 102, 103 à 20 chaque pas de calcul du modèle H(x), on les enchaîne avec un laps de temps défini pour le modèle H(x). Ce laps de temps peut être variable en fonction du point de fonctionnement courant, par exemple, en fonction du régime moteur thermique. Ainsi, par exemple, si l'on considère un modèle H(x) appelé toutes les 10ms mais qui nécessite un temps de convergence moyen de 150ms, la première étape 100 commençant au temps 25 zéro, la deuxième étape 101 où l'on introduit un biais sur une des entrées X1, X2, X3, X4 ne devra pas être exécutée au temps t=0+10ms mais au temps t=0+150ms, et la troisième étape 102 où l'on met à jour les sorties Y1, Y2 au temps t=0+150ms+150ms. Les deux variantes présentées précédemment pourront bien entendu être combinées entre elles. [0050] On précise qu'en fonction des besoins d'utilisation des sensibilités, le procédé 30 selon l'invention peut appliquer des valeurs de biais fixes ou bien variables dynamiquement. D'une manière générale, pour la plupart des cas, le biais peut être choisi fixe et de valeur relativement faible. La limite minimale du biais est telle que l'introduction du biais doit générer une réponse du modèle H(x). REVENDICATIONS: 1. Procédé de contrôle permettant le calcul d'une sensibilité d'au moins un modèle de comportement (H(x)) d'un actionneur de moteur thermique d'un véhicule automobile ayant au moins une sortie (Y1, Y2) déterminée en fonction d'au moins deux entrées (X1, X2, X3, X4), ledit modèle (H(x)) étant implémenté dans un calculateur moteur (1), ledit procédé comportant: - une étape (100) de calcul de la sortie (Y1) dudit modèle (H(x)) sans biais et de stockage de ladite valeur de sortie (Y1) comme valeur de référence, - une étape (101) de calcul de la sortie (Y1) dudit modèle (H(x)) avec un biais appliqué sur une première entrée (X1) de manière à déterminer une sensibilité de ladite sortie du modèle (Y1) à la variation de ladite première entrée (X1), au moins une fonction (Fonc) du calculateur moteur (1) utilisant ledit modèle (H(x)) étant alimentée par ladite valeur de référence de ladite sortie (Y1) déterminée pendant l'étape précédente (100), - une étape (103) de calcul de la sortie (Y1) dudit modèle (H(x)) avec un biais appliqué sur une deuxième entrée (X2) de manière à déterminer la sensibilité de ladite sortie (Y1) dudit modèle (H(x)) à la variation de ladite deuxième entrée (X2). 2. Procédé selon la revendication 1, caractérisé en ce qu'il comporte encore une étape (104) de commande dudit actionneur de moteur thermique en tenant compte desdites sensibilités dudit modèle (H(x)) préalablement déterminées. 3. Procédé selon la revendication 1 ou la revendication 2, caractérisé en ce que, entre les deux étapes (101, 103) de calcul de la sortie dudit modèle (H(x)) avec un biais appliqué sur une de ses entrées, ledit procédé comporte: - une étape (102) de mise à jour de ladite valeur de référence de ladite sortie (Y1) dudit modèle (H(x)) sans biais en retirant le biais sur ladite première entrée (X1), - la ou les fonctions (Fonc) dudit calculateur moteur (1) utilisant ledit modèle (H(x)) étant alimentées par ladite valeur de référence de la sortie (Y1) mise à jour lors de la mise en oeuvre de l'étape (103) de calcul de la sortie (Y1) du modèle (H(x)) avec le biais appliqué sur la deuxième entrée (X2). 3062884 12 4. Procédé selon l'une des revendications 1 à 3, caractérisé en ce que chaque étape de calcul (100-103) est mise en oeuvre lors d'un appel de tâche distinct correspondant du modèle de comportement (H(x)) de l'actionneur de moteur thermique. 5. Procédé selon l'une des revendications 1 à 3, caractérisé en ce que dans le cas 5 où ledit modèle de comportement (H(x)) de l'actionneur de moteur thermique nécessite une durée de convergence, les étapes de calcul (100-103) sont mises en oeuvre sur un laps de temps défini pour ledit modèle (H(x)). 6. Procédé selon la revendication 5, caractérisé en ce que ledit laps de temps est 10 variable en fonction d'un point de fonctionnement courant, par exemple, en fonction d'un régime moteur thermique. 7. Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que des valeurs des biais appliqués sur les entrées du modèle (H(x)) sont fixes. 15 8. Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que des valeurs des biais appliqués sur les entrées du modèle (H(x)) sont variables dynamiquement. 9. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce qu'une 20 limite minimale des biais est telle qu'une introduction des biais génère une réponse dudit modèle (H(x)). 10. Calculateur moteur (1) pour véhicule automobile comportant une mémoire (2) stockant des instructions logicielles pour la mise en oeuvre du procédé de contrôle d'une 25 sensibilité d'au moins un modèle de comportement d'un actionneur de moteur thermique du véhicule automobile tel que défini selon l'une quelconque des revendications précédentes.