FR 3040156 A1 20170224 FR 1659974 A 20161014 « Procédé de commande d'actionneur de frein automatisé et véhicule ou frein de véhicule mettant en oeuvre ce procédé » L'invention se rapporte à un procédé de commande d'un actionneur de 5 frein de véhicule, mis en oeuvre par un calculateur numérique embarqué au sein dudit véhicule. Ce procédé de commande comprend une itération des étapes suivantes : à partir d'une valeur de l'effort de serrage recherché, détermination d'une valeur de consigne y*(t0) visée pour une grandeur de sortie ; application 10 d'une grandeur de commande de valeur u(t0) à un instant tO ; mesure d'une valeur constatée y(t0) de la grandeur de sortie, produite au sein de l'actionneur sous l'effet de ladite application, et tant qu'une condition d'ajustement est vraie, calcul d'une nouvelle valeur u(t1) pour la grandeur de commande, à un nouvel instant t1, de façon à réaliser une commande 15 sans modélisation de l'actionneur par correction Proportionnelle-Intégrale de type « Intelligent », dite « commande sans modèle ». Elle concerne en outre un tel procédé appliqué à un frein de stationnement automatisé. L'invention concerne aussi un véhicule ou frein de véhicule mettant en 20 oeuvre ce procédé. Etat de la technique Pour faire commander le fonctionnement d'un frein motorisé dans un véhicule, par exemple par un calculateur embarqué, il est connu d'utiliser 25 une méthode consistant à utiliser une modélisation mathématique du comportement du frein et de l'actionneur en fonction des paramètres de commandes qu'il reçoit, par exemple tension et intensité électriques. À partir de la modélisation de ce comportement, on crée une loi de contrôle qui sera utilisée pour commander l'actionneur en lui envoyant les 30 paramètres électriques adéquats. Cette loi doit cependant elle-même être ajustée et paramétrée pour que le comportement réel, en situation, réponde aux attentes du concepteur, sans pour autant nécessiter de trop grandes performances de la part des 3040156 - 2 - composants, par exemple en puissance moteur et en résistance mécanique des pièces, en raison des coûts et des contraintes associées à une telle exigence en matière de composants. Actuellement, l'estimation de l'effort de serrage appliqué est basée sur un rétrocalcul utilisant des équations physiques qui décrivent la structure physique du frein et de l'actionneur, par exemple le taux de démultiplication des différents composants mis en jeu. Un but de l'invention est de pallier tout ou partie des inconvénients de l'art antérieur, et en particulier d'améliorer la fiabilité et la précision du comportement obtenu en situation réelle, en particulier en conditions sévères et pour des variations environnementales importantes. Exposé de l'invention Ce but est atteint avec, selon un premier aspect de l'invention, un procédé de commande d'un actionneur d'un frein de véhicule, mis en oeuvre par un calculateur numérique embarqué au sein dudit véhicule pour produire une application par ledit frein d'un effort de serrage déterminé, comprenant une itération des étapes suivantes : - à partir d'une valeur de l'effort de serrage recherché, détermination d'une valeur de consigne y*(t0) visée pour une grandeur de sortie, application d'une grandeur de commande de valeur u(t0) à un instant tO en entrée de l'actionneur, - mesure d'une valeur constatée y(t0) de la grandeur de sortie, produite au sein de l'actionneur sous l'effet de ladite application, et tant qu'une condition d'ajustement est vraie, calcul d'une nouvelle valeur u(n) pour la grandeur de commande, à un nouvel instant t1, et application de cette nouvelle valeur en entrée de l'actionneur, laquelle valeur dépend de : - la valeur u(t0) de la grandeur de commande, la dérivée ST(tO) de la grandeur de sortie, - le taux de variation 3;*(t1) de la valeur de consigne, entre tO et 3040156 3- - un facteur proportionnel à l'erreur de suivi de la grandeur de sortie, - un facteur proportionnel à l'intégrale de l'erreur de suivi de la grandeur de sortie, ladite nouvelle valeur u(t1) pouvant s'écrire sous la forme suivante : ri G (t1) 3;* (t1) u(t1) = a + Kpe(t1) + f a e(t)dt ro avec G(t1) = (t0) - a u(t0) 10 où : - a, Kp, et Ki sont des coefficients prédéterminés, - e(t) = y(t) - y*(t) est l'erreur de suivi ; de façon à réaliser une commande sans modélisation de l'actionneur par correction Proportionnelle-Intégrale de type « Intelligent », dite 15 « commande sans modèle ». L'invention propose ainsi un procédé de commande d'un actionneur de frein motorisé, en particulier mais non obligatoirement pour commander le fonctionnement d'un frein de stationnement automatisé. La condition d'ajustement peut être variée. Elle est par exemple définie 20 par rapport à l'erreur de suivi. De préférence, le procédé se poursuit tant que la valeur absolue de l'erreur de suivi diminue. Il s'agit d'un modèle que l'on pourrait qualifié d'empirique, au sens où il ne nécessite pas, ou moins, de modéliser le comportement du frein de façon mathématiquement précise. 25 Cette méthode utilise des données déterminées et mises au point par une forme d'apprentissage sur un prototype, avec enregistrement des résultats constaté en association avec les paramètres appliqués pour obtenir ces résultats. Il n'est plus nécessaire de connaître toutes les équations de fonctionnement, ou moins nécessaire. 30 Ce type de commande sans modèle mathématique est possible en particulier car il s'applique à un système au fonctionnement linéaire par 3040156 4- exemple qui ne comporte pas de came à changement brusque ou de composants fonctionnant en tout ou rien. Cette méthode peut être mise en oeuvre sans capteur "préféré", en "open loop". Elle peut aussi être mise en oeuvre avec choix d'un capteur préféré, 5 par exemple avec un capteur à effet Hall sur la vitesse, qui est lié à la position de l'actionneur. Le mode de commande du moteur électrique se fait de préférence en modulation de largeur d'impulsion ("Pulse Width Modulation"), en particulier pour un moteur à courant continu (DC), par exemple avec une commande 10 sériée en courant. Il peut se faire aussi en commande vectorielle, par exemple pour un moteur sans balais ("brushless"). Cette méthode permet une meilleure robustesse et fiabilité des résultats obtenus, en particulier en situations réelles, et une meilleure précision du niveau d'effort appliqué, avec moins d'influence des variations 15 environnementales. Des modes de réalisation variés de l'invention sont prévus, intégrant selon l'ensemble de leurs combinaisons possibles les différentes caractéristiques optionnelles exposées ici. Dans un mode de réalisation, la consigne y* est déterminée à partir d'une 20 liste de correspondance comprenant une pluralité de valeurs d'intensité associées chacune à une valeur d'un effort de serrage produite par ladite intensité au sein dudit actionneur. Avantageusement, le procédé de commande comprend, à partir d'une valeur initiale d'intensité choisie en fonction d'un effort de serrage à 25 appliquer, une itération des étapes suivantes : comparaison de ladite valeur initiale avec une base de données formant la liste de correspondances et contenant une série de valeurs d'intensités associées à un effort de serrage obtenu lors d'une phase préalable d'apprentissage, pour calculer une valeur modifiée 30 d'intensité selon un coefficient d'ajustement, envoi d'une commande d'alimentation du moteur (par exemple par réglage de la tension) basée sur cette valeur modifiée d'intensité, selon au moins une méthode déterminée de gestion du moteur 3040156 - 5 - notamment en modulation de largeur d'impulsion ou en commande vectorielle, mesure d'au moins un paramètre de résultat au sein de l'actionneur, estimation d'un gain constaté pour ladite valeur modifiée d'intensité, 5 en fonction dudit gain constaté, comparaison avec ladite base de données pour calculer une valeur modifiée du coefficient d'ajustement. De préférence, le procédé de commande une mesure d'un effort de serrage obtenu après application de la grandeur de commande u à partir de 10 la consigne y* et une mise à jour de l'effort de serrage associé à ladite consigne y* dans la liste de correspondance. Selon une particularité, la consigne y* est une consigne d'intensité de courant à obtenir dans l'actionneur de frein. Selon une possibilité, la valeur constatée y est obtenue par mesure de 15 l'intensité de courant telle que générée dans l'actionneur de frein. Cette mesure est de préférence directe et/ou sans grandeur intermédiaire. De préférence, la dérivée y de la valeur constatée y est déterminée sous la forme d'une estimée y(t), où : 3! y(t) = - 7-73- (2T(t ny(r)ch- t-T 20 où T est un paramètre prédéterminé. Selon un mode de réalisation, l'estimée est calculée par discrétisation de son intégrale, par exemple par la méthode des trapèzes. Toute autre méthode peut être utilisée, notamment la méthode des rectangles. Selon un autre aspect de l'invention, il est proposé un véhicule ou frein 25 de véhicule, comprenant au moins un calculateur programmé pour commander un ou plusieurs actionneurs selon un procédé selon le premier aspect de l'invention ou l'un de ses perfectionnements. Description des figures 30 D'autres particularités et avantages de l'invention apparaîtront à la lecture de la description détaillée de mises en oeuvre et de modes de 3040156 -6 réalisation nullement limitatifs, au regard de figures annexées sur lesquelles : - la figure 1 est un diagramme illustrant un mode de réalisation d'un procédé 100 selon l'invention, 5 - les figures 2 à 6 sont des planches de présentation qui précisent l'application de la méthode à un actionneur de frein en général, - les figures 7 à 12 sont des planches de présentation qui précisent l'application de la méthode à un actionneur de frein de stationnement automatisé. 10 Description de l'invention Ces modes de réalisation n'étant nullement limitatifs, on pourra notamment réaliser des variantes de l'invention ne comprenant qu'une sélection de caractéristiques décrites par la suite, telles que décrites ou 15 généralisées, isolées des autres caractéristiques décrites, si cette sélection de caractéristiques est suffisante pour conférer un avantage technique ou pour différencier l'invention par rapport à l'état de la technique. Le procédé 100 selon l'invention commande un bloc actionneur 1, aussi nommé « System » sur la figure, d'un frein de véhicule. 20 Selon un mode de réalisation d'une mise en oeuvre du procédé 100, le procédé 100 est mis en oeuvre par un calculateur numérique embarqué au sein du véhicule, notamment pour produire une application par le frein d'un effort de serrage déterminé. Le procédé 100 comprend une itération des étapes suivantes : 25 à partir d'une valeur de l'effort de serrage recherché, une détermination 102 d'une intensité y*(t0), réprésentée par la lettre I* sur la figure 1, visée pour une grandeur de sortie, une application 104 d'une tension de commande de valeur u(tO) à un instant tO, 30 une mesure 106 d'une valeur constatée y(t0), représentée par la lettre I sur la figure 1, de la grandeur de sortie, produite au sein de l'actionneur 1 sous l'effet de ladite application, et 3040156 7- tant qu'une condition d'ajustement est vraie, calcul d'une nouvelle valeur u(t1) pour la grandeur de commande, à un nouvel instant t1, laquelle valeur dépend de : - la valeur u(t0) de la grandeur de commande, 5 - la dérivée jr(tO) de la grandeur de sortie, - le taux de variation 3;*(t1) de la valeur de consigne, entre tO et t1, - un facteur proportionnel à l'erreur de suivi de la grandeur de sortie, 10 - un facteur proportionnel à l'intégrale de l'erreur de suivi de la grandeur de sortie. Dans l'exemple représenté sur la figure 1, à l'index k, c'est-à-dire un instant t1, le calcul de la nouvelle valeur est déterminé par la mise en oeuvre de trois blocs fonctionnels, respectivement un bloc correcteur 2 de 15 type intégrateur proportionnel intelligent, un bloc générateur de tension 3 et un bloc estimateur 4, en plus de l'actionneur 1. Le bloc actionneur 1, aussi nommé « System » sur la figure, comporte, selon la représentation de la figure 1, une entrée U (non représentée) qui est une tension, reliée en entrée à une tension notée « tension [V] » sur la 20 figure, et présente une sortie I, qui est une intensité, reliée à une sortie notée « intensity [A] » sur la figure. Le bloc correcteur 2, aussi nommé « i-PI » sur la figure, comporte trois entrées, respectivement notée I*, G(k), et I, ainsi qu'une sortie notée u(k). L'entrée I* est reçue de l'étape 102 de détermination de consigne. L'entrée 25 G(k) est reliée à une sortie G(k) du bloc estimateur 3. La sortie u(k) du bloc correcteur 2 est égale à une sommation : - d'un signal provenant d'un correcteur de type proportionnel intégral sur une commande d'entrée I*, qui peut s'écrire : t1 Kpe(t1)+ KL e(t)dt, to 30 avec Kr, et Ki qui sont des coefficients prédéterminés de gain du correcteur de type proportionnel intégral, et e(t) = I(k) -r(k) qui représente l'erreur de suivi ; et 3040156 - 8 - d'un signal provenant de la multiplication par un gain égal à l'inverse d'un coefficient a, d'un signal d'écart entre celui provenant d'un taux de variation de la consigne et le signal reçu sur l'entrée G(k), qui peut s'écrire : 5 G (k) (i*) (k) a a avec (i*)(k) -1* (k-1) ou Ts est une durée sur laquelle le taux de Ts variation est estimée. En pratique, la durée Ts, qui est égale à la durée d'un intervalle d'échantillonnage, est de l'ordre la microseconde. 10 Ainsi, la sortie u(k) du bloc correcteur 2 peut s'écrire : u(k) = - G (t1) + (y. *)(t1) + Kpe (t1) + Kif e (t)dt a a Le bloc générateur de tension 3, aussi nommé « PWM » sur la figure, comporte une entrée notée u(k), et deux sorties, respectivement notée u(k) et U. La sortie u(k) est directement reliée à l'entrée u(k). Dans l'exemple 15 représenté, la batterie électrique de stockage, notée « battery », utilisée pour générer la tension U est de 12 Volts. Aussi, la sortie U est égale à 12 fois la valeur de l'entrée u(k). Le bloc estimateur 4, aussi nommé « Estimation of G » sur la figure, comporte deux entrées, notées respectivement u(k) et I. L'entrée u(k) est 20 reliée à la sortie u(k) du bloc générateur de tension 3. L'entrée I est reliée à la sortie I du bloc actionneur 1. La sortie G(k) du bloc estimateur 4 est égale à la différence entre : un signal reçu d'un bloc de différentiation numérique, noté « numerical differentiation » sur la figure 1, qui peut s'écrire i(t0), 25 et - une multiplication par le gain a d'un signal décalé du signal d'entrée u(k), d'un registre d'une unité temporelle. Ainsi, la sortie G(k) du bloc estimateur 4 peut s'écrire : G(k) = (k) - a u(k - 1) 30 Ainsi, la valeur 1(k) du bloc estimateur 4 peut être estimée par : 3040156 9 t --D- ! /(t) = -T3 f (2T(t - r) - T)y(T)d-c t-T Dans l'exemple représenté, cette intégrale est typiquement calculée sur une fenêtre temporelle T de l'ordre de quelques dizaine de fois la durée de la période d'échantillonnage. 5 Toujours à tire d'exemple, l'intégrale peut être calculée par la méthode numérique connue sous le nom de méthode des trapèzes. L'utilisation de cette commande sans modèle est justifiée ci-après. Elle repose sur la possibilité de pouvoir écrire sur un court laps de temps, 10 l'équation « phénoménologique » : y(v) = G + au où : - y est l'ordre de dérivation de y, en général égal à 1 ou 2, choisi par l'opérateur (ce n'est pas l'ordre de dérivation maximum de y qui 15 reste inconnu), - le paramètre constant a est fixés par l'opérateur pour que les valeurs numériques \au et y(v) aient le même ordre de grandeur, - G qui contient toutes les informations structurelles sur le laps de temps choisi. 20 La valeur numérique de G, pour un index k peut être estimée à partir de l'équation suivante : [G(k)le = [y(v) (k)Je- au(k - 1) où [G(k)je et ly(v) (k)le sont estiméès à l'index k, c'est-à-dire à l'instant kTs où Ts est la durée de la période d'échantillonnage. 25 La valeur numérique de Ly(v) (k)Je peut être estimée. Celle-ci peut être obtenue par une mesure de valeur de la sortie y. Ces mesures peuvent être affectée par du bruit. Il est ainsi préférable d'utiliser une différentiation numérique, qui comprend un filtrage du bruit et une dérivation. 30 3040156 - 10 - La figure 2 illustre le principe général de la mise en oeuvre d'une commande sans modèle qui se base sur l'équation donnée par le modèle phénoménologie : 3,0') = G + au 5 Sur la figure 2, la consigne est notée y*, la commande à l'index k est notée u(k) et la sortie, sous l'effet de la commande u(k), d'un système, noté « unmodeled plant » dont le modèle n'est pas connu, est notée y. L'erreur de poursuite y* - y est notée e. Toujours sur la figure 2, la commande est déterminée, par un bloc, noté 10 « Model-free control », à partir du signal y* et du signal y. Plus précisément, ledit bloc comprend quatre sous-bloc : un premier sous-bloc, noté « i-PI », recevant les signaux y* et y ainsi qu'un signal [G(k)]e, et générant un signal u(k), à partir de l'équation : 15 G(k) (y*)(k) u(k) - + Kpe(k) + f e(k) a a un deuxième sous-bloc, noté « Delay », recevant un signal u(k) et fournissant un signal u(k - 1), un troisième sous-bloc, noté « Numerical Differentiation », recevant 20 un signal y et générant un signal [y(v) (k)]e, à partir de l'équation : J (2T(t - r) - T)y(r)dr dx = T3 25 t-T signal un quatrième sous bloc, noté « Estimation of G », recevant le u(k - 1) et générant le signal [G(k)],, à partir de l'équation : [G(k)]e = [y(v) (k)]e- au(k - 1) /CO = (2T(t - r) - T)y(r)dr t-T Différentiation numérique : la figure 3 permet de justifier l'équation mise en oeuvre le troisième sous-bloc de la figure 2. 30 Une linéarisation à l'instant t du signal y permet d'écrire : t dy 3! 3040156 dy y(t) = ao + ait, avec ai =t (t) Une transformation de Laplace de cette équation donne : ao, al S S2 Y(s) = + -, d ou sYm = ao - En dérivant l'équation de droite par rapport à la variable s, on obtient : 5 Y(s) + s d(s) = - al, d'où s-2Y(s) + s-1Y(s) = -s-4a1 ds s2 D'où : t 3! a1 = -T3 (2T(t - t) - T)y(r)d-c T3 t-T D'où : T dy 3! dt (t) = - (T - 2t)y(t)d-c o Où l'intervalle [0, T] est une fenêtre de temps assez petite qui est déplacée pour estimer la dérivée à chaque instant. La méthode des trapèzes peut être utilisée pour calculer, en pratique, cette intégrale. 15 Cette méthode pour déterminer la dérivée a pour avantage de pouvoir être utilisée pour des déterminations de dérivées d'ordre plus élevé. La figure 4 est une illustration de l'intérêt de l'utilisation de cette méthode par rapport à une méthode d'estimation de la dérivée plus 20 classique, telle que celle mettant en oeuvre un filtre passe bas du premier ordre suivi d'une détermination de la dérivée selon la méthode d'Euler. Le graphe de gauche présente trois courbes d'intensité entre -1,5 A et 2 A, en fonction du temps sur l'intervalle allant de 0 à 8 secondes. La graphe de droite est un zoom sur l'intervalle allant de 5 à 7 secondes et -0,2 A à 25 1,6 A dwgraphe de gauche. On observe que la courbe illustrant la mise en oeuvre de la différentiation numérique est en moyenne plus proche de la courbe parfaite de la dérivée attendue, c'est à dire d'un cosinus, que la courbe mettant en oeuvre un filtre 10 3040156 - 12 - passe-bas du premier ordre suivi d'un calcul de la dérivée selon la méthode d'Euler. Cette courbe correspond donc bien à la dérivée attendue, et il n'est pas nécessaire de placer un filtre avant la dérivation. La figure 5 illustre les dispositifs mis en oeuvre pour produire les courbes de la figure 4. Plus précisément, on observe la génération d'un signal bruité à dériver, généré à partir d'un signal de sinusoïde auquel un signal de bruit blanc est ajouté. 10 On observe en sortie les trois signaux représentés sur la figure 4, c'est-à- dire un signal de cosinus (qui est la dérivée attendue d'un signal de sinusoïde), le signal selon la méthode d'Euler et le signal mettant en oeuvre une différentiation numérique. Le signal selon la méthode d'Euler met en oeuvre un filtre du premier 15 ordre, générant un signal filtré à partir du signal bruité, le signal filtré étant ensuite dérivé par application de la transformée en z de la dériation discrète, c'est-à-dire qe1. Tsz La dérivation selon la méthode présentée plus haut utilise une période T égale à 30 fois le temps Ts, c'est-à-dire 30 fois la durée de l'intervalle 20 d'échantillonnage. Choix de la fenêtre de différentiation numérique : la figure 6 illustre l'influence du choix de la valeur de la période T sur la détermination de la dérivée par la méthode de différentiation numérique. 25 Le graphe de gauche de la figure 6 représente une dérivée attendue, la courbe cosinus et 3 courbes obtenues pour différentes valeurs de T, respectivement 20 fois Ts, 30 fois Ts et 40 fois Ts. Le graphe de droite représente les distances des trois courbes obtenues à la dérivée attendue. 30 On observe que plus la valeur de T est grande, moins le bruit a d'influence sur la détermination de la dérivée. Toutefois, plus la valeur de T est grande, plus le délai est long avant d'obtenir la valeur attendue. 3040156 - 13 - En pratique, un compromis doit être trouvé pour déterminer la valeur de T. Application à un frein de stationnement automatique (« APB ») 5 Les figures 7 et 8 étudient la faisabilité de la mise en oeuvre d'une commande sans modèle pour un moteur DC. La figure 7 illustre le fait qu'un système dont la fonction de transfert présente un zéro dans le plan positif est un système qui n'a pas de minimum de phase. 10 La figure 8 utilise la modélisation d'un moteur DC selon le logiciel Simulink pour en déterminer la fonction de transfert et vérifier qu'elle ne présente pas de zéro dans le plan positif, et donc que le moteur DC est un système à minimum de phase. Il est donc possible de mettre en oeuvre une commande sans modèle pour le moteur DC. 15 Simulation : la figure 9 a pour objet d'illustrer le comportement du modèle sans commande par rapport à une variation au cours du temps d'une résistance d'un système comportant un moteur. Le graphe de gauche illustre trois courbes, respectivement une intensité 20 de référence, une intensité mesurée avec une résistance erronée et une intensité mesurée avec une résistance exacte, au cours du temps de 0,5 secondes à 4 secondes et de -5 A à 35 A. On observe que les trois courbes sont très proches. Ceci est confirmé par le graphe de droite, qui est un zoom du graphe de 25 gauche entre 0,95 secondes et 1,15 secondes et entre de 0 A à 40 A. On observe que même si la résistance varie de 0.35 ohms à 1 ohms, le comportement du processus de commande sans modèle est performant. En conséquence, le processus de commande sans modèle se comporte de façon satisfaisante même si le système évolue. 30 La figure 10 illustre l'influence de la détermination de calcul utilisé dans la différentiation numérique sur les performances de la méthode, pour une courbe relative à un frein de stationnement électrique. 3040156 - 14 - La figure 11 comporte deux graphes qui illustrent les erreurs obtenues entre intensité estimée et intensité de référence pour plusieurs ensemble de paramètres a, Kp et Ki. 5 La figure 12 illustre un mode de réalisation d'un système mettant en oeuvre un procédé de commande selon l'invention. Le système illustré sur la figure 12 comporte un moteur DC d'un actionneur de frein, pouvant recevoir une tension comprise entre 0 et 16 10 volts en provenance d'un circuit électronique. Le circuit électronique est capable de générer une tension régulée, à partir d'une batterie de 16 V, et d'une entrée PWM provenant d'un module de calcul électronique, programmé par exemple en Matlab ou Simulink. Le module de calcul électronique est programmé pour mettre en oeuvre, 15 lorsqu'il est exécuté par un ordinateur ou une unité de traitement, une partie du procédé de commande selon l'invention, c'est-à-dire tant qu'une condition d'ajustement est vraie, calcul d'une nouvelle valeur u(t1) pour la grandeur de commande, à un nouvel instant t1, et application de cette nouvelle valeur en entrée de l'actionneur, laquelle valeur dépend de : 20 la valeur u(t0) de la grandeur de commande, - la dérivée 3r(t0) de la grandeur de sortie, laquelle est déterminée à partir d'une valeur de courant I mesurée par un module capteur d'intensité au sein de l'actionneur, le taux de variation (y.*)(t1) de la valeur de consigne, entre tO et t1, 25 - un facteur proportionnel à l'erreur de suivi de la grandeur de sortie, - un facteur proportionnel à l'intégrale de l'erreur de suivi de la grandeur de sortie, ladite nouvelle valeur u(t1) pouvant s'écrire sous la forme suivante : G(t1) (y*)(t1) u(t1) = + Kpe(t1) + f e(t)dt a a to 30 avec où : G(t1) = 31(t0)- a u(t0) 3040156 - 15 - a, Kp, et Ki sont des coefficients prédéterminés, de façon à réaliser une commande sans modélisation de l'actionneur par correction Proportionnelle-Intégrale de type Intelligent >>, dite « commande sans modèle ». 5 Bien sûr, l'invention n'est pas limitée aux exemples qui viennent d'être décrits et de nombreux aménagements peuvent être apportés à ces exemples sans sortir du cadre de l'invention. De plus, les différentes caractéristiques, formes, variantes et modes de réalisation de l'invention 10 peuvent être associés les uns avec les autres selon diverses combinaisons dans la mesure où ils ne sont pas incompatibles ou exclusifs les uns des autres. REVENDICATIONS1. Procédé de commande d'un actionneur (1) d'un frein de véhicule, mis en oeuvre par un calculateur numérique embarqué au sein dudit véhicule pour produire une application par ledit frein d'un effort de serrage déterminé, comprenant une itération des étapes suivantes : à partir d'une valeur de l'effort de serrage recherché, détermination (102) d'une valeur de consigne y*(t0) (I*) visée pour une grandeur de sortie, application (104) d'une grandeur de commande de valeur u(tO) à un instant tO en entrée de l'actionneur (1), mesure d'une valeur constatée y(t0) (106) de la grandeur de sortie, produite au sein de l'actionneur sous l'effet de ladite application, et tant qu'une condition d'ajustement est vraie, calcul d'une nouvelle valeur u(t1) pour la grandeur de commande, à un nouvel instant t1, et application de cette nouvelle valeur en entrée de l'actionneur (1), laquelle valeur dépend de : - la valeur u(t0) de la grandeur de commande, - la dérivée ST(tO) de la grandeur de sortie, - le taux de variation 3;*(t1) de la valeur de consigne, entre tO et t1, - un facteur proportionnel à l'erreur de suivi de la grandeur de sortie, - un facteur proportionnel à l'intégrale de l'erreur de suivi de la grandeur de sortie, ladite nouvelle valeur u(t1) pouvant s'écrire sous la forme suivante : ti G (t1) 3;* (t1) u(t1) = + Kpe(t1) + e(t)dt a a to avec G(t1) = ji(t0) - a u(t0) où : 3040156 - 17 - - a, Kp, et Ki sont des coefficients prédéterminés, - e(t) = y(t) - y*(t) est l'erreur de suivi ; de façon à réaliser une commande sans modélisation de l'actionneur par correction Proportionnelle-Intégrale de type « Intelligent », dite 5 « commande sans modèle ». 2. Procédé selon la revendication 1, dans lequel la consigne y* est déterminée à partir d'une liste de correspondance comprenant une pluralité de valeurs d'intensité associées chacune à une valeur 10 d'effort de serrage produite par ladite intensité au sein dudit actionneur. 3. Procédé selon la revendication 2, caractérisé en ce qu'il comprend, à partir d'une valeur initiale d'intensité choisie en fonction d'un effort 15 de serrage à appliquer, une itération des étapes suivantes : comparaison de ladite valeur initiale avec une base de données formant la liste de correspondances et contenant une série de valeurs d'intensités associées à un effort de serrage obtenu lors d'une phase préalable d'apprentissage, pour calculer une valeur 20 modifiée d'intensité selon un coefficient d'ajustement, envoi d'une commande d'alimentation du moteur basée sur cette valeur modifiée d'intensité, selon au moins une méthode déterminée de gestion du moteur notamment en modulation de largeur d'impulsion ou en commande vectorielle, 25 mesure d'au moins un paramètre de résultat au sein de l'actionneur, estimation d'un gain constaté pour ladite valeur modifiée d'intensité, en fonction dudit gain constaté, comparaison avec ladite base de 30 données pour calculer une valeur modifiée du coefficient d'ajustement 3040156 - 18 - 4. Procédé selon la revendication 2 ou 3, comprenant en outre une mesure d'un effort de serrage obtenu après application de la grandeur de commande u à partir de la consigne y* et une mise à jour de l'effort de serrage associé à ladite consigne y* dans la liste de 5 correspondance. 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel la consigne y* est une consigne d'intensité de courant à obtenir dans l'actionneur de frein. 10 6. Procédé selon la revendication 5, dans lequel la valeur constatée y est obtenue par mesure de l'intensité de courant telle que générée dans l'actionneur de frein. 15 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel la dérivée )*/ de la valeur constatée y est déterminée sous la forme d'une estimée y(t), où : y(t) = -T3 (2T(t - - nyWch- t-T où T est un paramètre prédéterminé. 20 8. Procédé selon la revendication 7, dans lequel l'estimée est calculée par discrétisation de son intégrale, par exemple par la méthode des trapèzes. 25 9. Véhicule ou frein de véhicule, comprenant au moins un calculateur programmé pour commander un ou plusieurs actionneurs selon un procédé selon l'une quelconque des revendications précédentes.