FR 3002048 A1 20140815 FR 1351268 A 20130214 Procédé de détection améliorée de collision d'un robot avec son environnement, système et produit programme d'ordinateur mettant en oeuvre le procédé L'invention concerne le domaine des robots et des systèmes de contrôle associés. Plus précisément, l'invention a pour objet un procédé de détection des collisions d'un robot avec son environnement. L'invention s'applique avantageusement à tout type de robot manipulateur, en particulier des robots manipulateurs sériels mais également 10 des robots manipulateur à structure arborescente ou à structure cinématique fermée. Ces dernières années, de nombreuses architectures de robots légers ont été proposées. Cette nouvelle génération de robot a pour but avoué 15 d'interagir avec des opérateurs humains lors de la réalisation d'une tâche. Le robot et l'opérateur partagent donc le même espace de travail, et ceci implique de nouvelles problématiques portant sur la sécurité de l'opérateur. En effet, des mouvements à grande dynamique du robot peuvent notamment limiter la zone d'accès rapprochée du robot. 20 Une problématique générale concernant les robots manipulateurs porte sur la détection de collisions entre le robot et son environnement. En effet, dans le but d'améliorer la sûreté de fonctionnement du robot, il est important de pouvoir détecter rapidement des collisions entre le robot et son environnement afin de minimiser les dommages éventuels en appliquant des 25 stratégies post-impact adaptées. Les algorithmes connus de détection de collision sont le plus souvent basés sur la comparaison de mesures avec un modèle, permettant de créer un signal appelé « résidu » et qui constitue une image de la collision. La 30 modélisation mathématique du système n'étant jamais parfaitement représentative du comportement réel du robot, le résidu est entaché d'erreurs, et les stratégies de détection doivent imposer l'utilisation de marges de sécurité (matérialisées dans la pratique par des seuils) afin d'être robustes vis-à-vis de ces erreurs. Cela évite l'occurrence de fausses alertes. Mais du fait de ces marges conservatives, le robot perd sa sensibilité aux collisions. La présente invention vise ainsi à améliorer les méthodes de détection de collision basées sur la génération et l'évaluation d'un résidu en supprimant ou limitant l'impact des erreurs de modélisation du système. Comme cela sera développé par la suite, les méthodes de détection de collision entre un robot et son environnement sont généralement conçus selon une approche en deux étapes : une première étape de génération du résidu, terme qui contient l'information associée au phénomène de collision, et une seconde étape d'évaluation de ce résidu qui consiste notamment à prendre une décision sur l'occurrence d'une collision en fonction de la valeur du résidu. Concernant la première étape de génération du résidu, plusieurs approches sont inspirées du domaine scientifique de l'automatique. On peut notamment citer le document [1] qui exploite un filtre estimateur de Kalman, dont la grandeur estimée est utilisée comme résidu. Le document [2] décrit une approche plus spécifique qui est issue de la détection de faute. L'article [3] décrit un procédé d'observation de diagnostic utilisé sur un bras de robot afin de générer le résidu. Dans le domaine de la robotique, la stratégie de génération de résidu la plus utilisée consiste à comparer une mesure des couples articulaires T avec leur estimation î déduite à partir du modèle dynamique inverse du robot. Cette technique a été améliorée en utilisant l'équation du moment généralisé d'un robot dans les articles [4] et [5], ce qui évite le calcul en ligne de l'accélération et réduit l'influence des bruits de mesures. Une autre stratégie décrite dans [6], basée sur le filtrage du couple, permet d'éviter le calcul de l'accélération. Toutes ces méthodes sont basées sur un modèle du robot, et ceci implique donc une sensibilité aux erreurs de modélisation. Concernant la seconde étape de détection de collision, l'évaluation du 5 résidu permet de prendre en compte les incertitudes du modèle. Elle consiste à vérifier que le résidu est en dessous d'un certain seuil. Si le résidu dépasse ce seuil, l'algorithme détecte la collision. Ce seuil permet d'établir une marge afin de rendre la méthode insensible aux erreurs de modélisation et d'éviter les fausses alertes. La technique la plus simple, décrite par exemple dans le 10 document [5], consiste à comparer le résidu avec un seuil constant ou statique, qui représente donc l'erreur de modélisation maximum. Cette méthode réduit considérablement la sensibilité du robot à une collision car un seuil statique ne permet pas de discriminer les variations du résidu liées aux erreurs de modélisation. 15 Des stratégies plus évoluées utilisent des seuils dynamiques, dont le niveau est adapté en ligne. On peut citer l'article [3] dans lequel un seuil adaptatif est généré à l'aide de techniques fondées sur la logique floue. Cette solution est assez complexe à mettre en place, car elle nécessite un recueil de données afin de développer les lois de logique floue. Dans les articles [6] 20 et [7], plusieurs seuils dynamiques sont proposés, prenant en compte les incertitudes paramétriques. Ces solutions sont uniquement basées sur le modèle dynamique inverse du robot qui, même si les paramètres du modèle sont parfaitement estimés, ne constitue pas une modélisation suffisante. En effet, il ne prend pas en compte certains phénomènes physiques propres aux 25 robots manipulateurs, en particulier les robots sériels légers, tels que les phénomènes de flexibilités. L'invention propose un procédé de détection de collisions qui prend en 30 compte les erreurs de modélisation du résidu et permet de résoudre les limitations des solutions connues exposées ci-dessus. Les erreurs de modélisation, qui comprennent à la fois des incertitudes liées au paramètres du modèle du robot mais aussi les incertitudes non paramétriques liées à des facteurs extérieurs au modèle choisi, sont identifiées et filtrées afin de générer un résidu filtré. Un seuil dynamique est déterminé en ligne pour adapter le critère de détection aux incertitudes liées aux erreurs de modélisation. L'application du procédé selon l'invention permet de conserver une bonne sensibilité aux collisions tout en diminuant les phénomènes de fausse alarme liés aux erreurs de modélisation du modèle dynamique inverse utilisé 10 pour générer le résidu. Les performances du procédé selon l'invention sont notamment améliorées pour des robots soumis à des dynamiques fortes en vitesse ou accélération. 15 L'invention a pour objet un procédé de détection de collision entre un robot composé d'une pluralité de corps reliés entre eux par au moins une articulation et son environnement ledit procédé comprenant les étapes suivantes : - Générer, à partir d'un modèle dynamique du robot, un signal 20 représentatif des collisions entre le robot et son environnement, ledit signal étant appelé résidu r et comprenant autant de composantes que d'articulations du robot, ledit procédé étant caractérisé en ce qu'il comporte en outre les étapes suivantes : 25 - Réaliser un filtrage passe-haut adaptatif du résidu r pour rendre le résidu indépendant d'incertitudes paramétriques ou non paramétriques liées à des phénomènes basses-fréquences, - Déterminer un seuil adaptatif T composé d'au moins un premier terme Tot, dynamique, égal à une information sur les incertitudes 30 paramétriques entre ledit modèle et le comportement réel du robot, lesdites incertitudes paramétriques étant liées à une première variable ei caractéristique de l'état dudit robot parmi les variables suivantes : la position, la vitesse ou l'accélération d'une articulation du robot, ou une fonction linéaire ou non linéaire de l'une de ces variables ou d'une combinaison de ces variables, - Comparer le résidu filtré au seuil adaptatif T pour en déduire l'existence ou non d'une collision. Selon un aspect particulier de l'invention, ledit seuil adaptatif T est composé de la somme de plusieurs termes Tot, Tot, Toi, dynamiques, chacun égaux à une information sur les incertitudes paramétriques entre ledit modèle et le comportement réel du robot, lesdites incertitudes paramétriques étant liées à une variable différente, caractéristique de l'état dudit robot parmi les variables suivantes : la position, la vitesse ou l'accélération d'un point fixe du robot, ou une fonction linéaire ou non linéaire de l'une de ces variables ou d'une combinaison de ces variables. Selon un autre aspect particulier de l'invention, ledit seuil adaptatif T comporte en outre un terme Tstatic statique configuré de sorte à être supérieur à un niveau de bruit de mesure. Selon un autre aspect particulier de l'invention, l'étape de filtrage est réalisée au moyen d'un algorithme des moindres carrés récursifs. Cette étape de filtrage peut comporter les sous-étapes suivantes : - Estimer, de façon récursive, les coefficients de la fonction de transfert du filtre passe-haut G0-1, - Filtrer le résidu r avec le filtre passe-haut estimé. Selon un autre aspect particulier de l'invention, une étape 25 supplémentaire de filtrage temporel est appliquée au résidu filtré pour chacune de ses composantes. L'étape supplémentaire de filtrage temporel peut être une étape de calcul de moyenne quadratique. Selon un autre aspect particulier de l'invention, l'étape de détermination d'au moins un premier terme Tot, Tot, Toi, dynamique du seuil 30 adaptatif T est réalisée au moyen d'un algorithme des moindres carrés récursifs. L'étape de détermination d'au moins un premier terme Tot, Tot, Toi, dynamique du seuil adaptatif T peut comporter les sous-étapes suivantes : - Estimer, de façon récursive, les coefficients de la fonction de transfert A modélisant les incertitudes paramétriques liées à ladite première variable ei caractéristique de l'état dudit robot, - Filtrer ladite première variable ei caractéristique de l'état dudit robot avec un filtre de fonction de transfert Ai estimée à l'étape précédente pour obtenir un terme TA-1, Tot, Toi, dynamique du seuil adaptatif T. Selon un autre aspect particulier de l'invention, une étape supplémentaire de filtrage temporel est appliquée audit terme TA-1, Tot, Toi, dynamique du seuil adaptatif T pour chacune de ses composantes. L'étape supplémentaire de filtrage temporel peut être une étape de calcul de moyenne quadratique. Selon un autre aspect particulier de l'invention, l'étape de comparaison du résidu filtré au seuil adaptatif T pour en déduire l'existence ou non d'une collision comprend les sous-étapes suivantes : - Comparer, pour chaque composante articulaire, le résidu filtré au seuil adaptatif T, - Conclure à l'existence d'une collision si, pour au moins K composantes, le résidu filtré est supérieur au seuil adaptatif T, K étant un nombre entier prédéterminé strictement positif et inférieur ou égal au nombre d'articulations du robot. Selon un autre aspect particulier de l'invention, l'étape de génération 25 du résidu est constituée des sous-étapes suivantes : - Déterminer, à partir d'une information sur l'état du robot et par le biais d'un modèle dynamique, une estimation des couples articulaires du robot, - Effectuer une mesure de l'état du robot, par exemple une mesure 30 des couples articulaires, - Calculer le résidu comme la différence entre l'estimation et la mesure de l'état du robot. Selon un autre aspect particulier de l'invention, les incertitudes paramétriques liées à l'accélération articulaire du robot sont des incertitudes 5 sur la matrice d'inertie du robot. Selon un autre aspect particulier de l'invention, les incertitudes paramétriques liées à la vitesse articulaire du robot sont des incertitudes sur la matrice des termes centrifuges et de Coriolis du robot et/ou sur les frottements visqueux. 10 Selon un autre aspect particulier de l'invention, une fonction non linéaire est la fonction signe ou la fonction exponentielle ou la fonction valeur absolue. Selon un autre aspect particulier de l'invention, les incertitudes paramétriques liées au signe de la vitesse articulaire du robot sont des 15 incertitudes sur les frottements secs. L'invention a également pour objet un programme d'ordinateur comportant des instructions pour l'exécution du procédé de détection de collision selon l'invention, lorsque le programme est exécuté par un processeur. 20 L'invention a également pour objet un support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution du procédé de détection de collision selon l'invention, lorsque le programme est exécuté par un processeur. L'invention a également pour objet un système de contrôle d'un robot 25 comprenant un organe de commande pour la manipulation du robot, une interface pour l'échange d'informations sur l'état du robot et un module de détection de collisions adapté pour mettre en oeuvre le procédé selon l'invention. D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit en relation aux dessins annexés qui représentent : - La figure 1, un organigramme illustrant les étapes de mise en oeuvre du procédé selon l'invention, - La figure 2, un schéma représentant plus en détail l'étape d'évaluation du résidu, - La figure 3, un diagramme illustrant l'amélioration de la fiabilité de la décision de détection de collision par application du procédé selon l'invention, - La figure 4, un schéma d'un système de contrôle d'un robot manipulateur selon l'invention. La figure 1 représente, sur un organigramme, les différentes étapes 15 de mise en oeuvre du procédé de détection de collisions selon l'invention. Le procédé selon l'invention reçoit en entrée, une information E sur l'état du robot, à tout instant ou sur un intervalle de temps donné. L'information E sur l'état du robot est constituée de variables d'état du robot. Ces variables peuvent par exemple prendre la forme du triplet lq4ell 20 comprenant la position q, la vitesse et l'accélération q pour chaque articulation du robot. Comme cela sera explicité plus en détail par la suite, l'information E sur l'état du robot peut également être constituée d'une seule des trois variables d'état précitées ou de deux de ces variables ou encore de toute fonction linéaire ou non linéaire d'une ou de plusieurs de ces trois 25 variables. En particulier mais pas uniquement, les fonctions non linéaires possibles sont la fonction signe, la fonction exponentielle, la fonction valeur absolue ou toute combinaison d'une ou plusieurs de ces fonctions. En outre l'information E sur l'état du robot comprend également une mesure -c des couples articulaires du robot ou une information permettant de 30 déduire cette mesure, par exemple une mesure du courant moteur sur chaque articulation du robot. Une première partie 100 du procédé selon l'invention consiste à générer un résidu r qui contient l'information sur les évènements associés à une collision du robot avec son environnement. Cette première partie 100 comprend une première étape 101 de mesure r des couples articulaires du robot, c'est-à-dire des couples mesurés au niveau de chaque articulation du robot. Comme expliqué ci-dessus, cette étape peut, par exemple, être réalisée à partir de mesures des courants moteurs sur chaque articulation. La première partie 100 du procédé comprend également une 10 deuxième étape 102 de détermination d'un modèle du comportement du robot conduisant à une estimation î des couples articulaires du robot. Différents modèles sont possibles et connus de l'Homme du métier. Un exemple particulier d'un modèle dynamique d'un robot sériel rigide est défini par les relations suivantes : 15 je = Â(qM + R (q, 0 + Tc (1) avec : I1(q,) = e(q,0 + Q(q) + ff (2) 20 q, ,g sont les vecteurs des estimations ou mesures respectives des positions, vitesses et accélérations articulaires. --È est un vecteur des estimées des couples articulaires. r, est un vecteur des couples de collision. Les forces appliquées sur l'outil terminal du robot peuvent être reflétées au niveau moteur par exemple par le 25 biais de la relation : Tc = j(q)TFc où F' est le torseur des forces extérieures appliquées au robot et J(q) est la matrice Jacobienne du robot. Â(q) est une estimée de la matrice d'inertie du robot. e(q, ) est une estimée de la matrice des termes centrifuges et de Coriolis. 30 Q(q) est un vecteur des estimées des couples de gravité. ff est un vecteur des estimées des couples de frottement. Les frottements peuvent, par exemple, être modélisés par un modèle de Coulomb selon la relation suivante : Tf = + Fssign(i') (2') où Fv est le coefficient des frottements visqueux et Fs le coefficient des frottements secs, sign() désigne la fonction signe. Le terme représente une estimée des frottements visqueux. Le terme Fssign(i') représente une estimée des frottements secs. L'inconnue dans l'équation (1) est le vecteur des couples de collision, 10 l'étape 102 du procédé selon l'invention consiste donc à estimer le vecteur des couples articulaires selon la relation (3) : = Â(qM +11(q,) (3) Le modèle dynamique du robot présenté ci-dessus est donné à titre 15 d'exemple illustratif mais ne saurait être interprété comme limitatif de la portée de l'invention. En effet, d'autres modèles peuvent être utilisés, par exemple, un modèle similaire à celui de l'équation (1) dans lequel l'un des deux termes A(q)ei ou H(q,) est ignoré. De façon générale une estimée d'un couple articulaire du robot est composée d'au moins un terme 20 quelconque parmi les termes suivants : un terme dépendant de l'accélération articulaire, un terme dépendant de la vitesse articulaire, un terme dépendant de la position articulaire, un terme dépendant d'une fonction linéaire ou non linéaire de l'un quelconque des trois termes précédents ou d'une combinaison de ces termes, par exemple la fonction signe, la fonction 25 exponentielle ou la fonction valeur absolue. La première partie 100 du procédé selon l'invention comporte en outre une troisième étape 103 de génération du résidu r par le calcul de la différence entre la mesure des couples articulaires produite par la première 30 étape 101 et l'estimée des couples articulaires produite par la deuxième étape 102 : r = - . Toute méthode connue de l'Homme du métier, alternative à celle décrite ci-dessus pour la première partie 100 du procédé selon l'invention, peut être utilisée en remplacement dans la mesure où elle permet d'obtenir un résidu, c'est-à-dire un signal représentatif d'une collision entre le robot et son environnement. Par exemple, le document [5] décrit une méthode de génération du résidu basée sur la comparaison de quantités de mouvement et non des couples articulaires. Cette méthode peut être utilisée en remplacement de la première partie 100 du procédé selon l'invention. A partir des relations (1) et (3), le résidu r peut être exprimé comme la somme des couples de collision -c, et d'au moins un terme correspondant aux erreurs de modélisation sur les modèles utilisés pour estimer les couples articulaires. Les erreurs de modélisation comprennent d'une part des incertitudes paramétriques et d'autre part des incertitudes non paramétriques. Les incertitudes paramétriques sont les erreurs bornées pouvant affecter les paramètres d'un modèle de structure connue, par exemple un modèle représenté par les relations (1) et (2) précitées. Ces erreurs sont dues à l'identification expérimentale ou à la variation intrinsèque des paramètres du modèle en fonction des conditions de fonctionnement, par exemple la variation de la charge utile, la température, le vieillissement. Les incertitudes non paramétriques regroupent les erreurs qui ne sont pas liées au modèle lui-même mais à d'autres phénomènes tels que les 25 flexibilités du robot ou le modèle de frottement utilisé. Dans l'exemple donné ci-dessus où l'estimation des couples articulaires est obtenue à l'aide de la relation (3), le résidu r peut être exprimé à l'aide des relations suivantes : 30 r = T - Î = Tc + AA(q)q + AH(q,) (4) Avec : DA(q) = A(q) - Â(q) AI-1(q,) = H(q,) - rl(q,) AA(q) et AH(q,) sont les erreurs entre les valeurs réelles et les valeurs estimées des termes employés pour modéliser le comportement du robot, dans l'exemple choisi il s'agit de la matrice d'inertie A du robot et du 5 vecteur H définie par la relation (2). Ces erreurs ont un impact direct sur la valeur du résidu qui ne peut plus être considéré comme strictement égal au vecteur des couples de collisions (pour chaque articulation). En conséquence, même en l'absence de collision, la valeur du résidu n'est pas toujours nulle ce qui pose le 10 problème du critère à appliquer pour évaluer le résidu et produire une décision sur l'occurrence d'une détection. Un but de l'invention est de réaliser un filtrage du résidu et une élaboration d'un seuil de détection dynamique qui permet de s'affranchir des 15 fluctuations liées aux erreurs entre le modèle dynamique du robot et son comportement réel. Ces traitements sont regroupés dans une seconde partie 200 du procédé selon l'invention qui comporte une étape 104 de filtrage du résidu r pour obtenir un résidu filtré rf, une étape 105 d'élaboration en ligne d'un seuil dynamique T de détection en fonction du résidu r et de l'état du 20 robot E et enfin une étape 106 de comparaison du résidu filtré rf au seuil de détection T afin d'élaborer une décision D sur la présence ou l'absence de collisions. Avant de décrire les étapes mises en oeuvre dans la seconde partie 25 200 du procédé selon l'invention, on introduit la modélisation suivante du résidu afin de mieux comprendre le principe à la base de l'invention. Pour cela, on se place dans le domaine des transformées en Z. La relation (5) donne la modélisation du résidu qui sert de base à l'invention. rm(z) = T c(z) + EiN=1 Si(z)ei(z) + G0(z)b(z) (5) Le vecteur du modèle du résidu rm(z) se compose de trois éléments. La première composante est le vecteur des couples de collision Tc(z) qui est 5 le signal que l'on cherche à détecter. La deuxième composante _Si(z)ei(z) correspond aux erreurs de modélisation, en particulier dues aux incertitudes paramétriques. Les vecteurs ei(z) représentent les informations sur l'état du robot, à savoir les variables d'état du robot tel que l'accélération, la vitesse, la position ou toute fonction linéaire ou non linéaire 10 de l'une de ces variables ou d'une combinaison de ces variables. La fonction Si(z) est la fonction de transfert des erreurs de modélisation associées à l'information ei(z). Autrement dit, dans l'exemple décrit ci-dessus où le modèle dynamique du robot est conforme à la relation (1), on voit que les incertitudes 15 paramétriques sont fonctions de trois termes distincts, l'accélération articulaire la vitesse articulaire et le signe de la vitesse articulaire sign(cj). Le modèle issu de la relation (5) s'écrit alors, pour cet exemple particulier : rm(z) = Tc (Z) + (z) (z) + S 2(z) e 2(z) + S 3 (z) e3 (z) + G 0 (z) b (z) (6) 20 Les termes ei(z), e2(z), e3(z) correspondent respectivement aux transformées en Z de l'accélération articulaire de la vitesse articulaire et du signe de la vitesse articulaire. La fonction de transfert Si(z) vise à modéliser les erreurs de 25 modélisation liées aux termes d'accélération dans la relation (4), à savoir le terme AA(q)ei . Ce terme est particulièrement significatif lorsque le robot est soumis à des trajectoires à fortes dynamiques. La fonction de transfert S2(z) vise à modéliser les incertitudes paramétriques liées aux termes de vitesse dans la relation (4), autrement dit 30 les composantes du terme AH (q , qui dépendent de la vitesse articulaire. Ce terme prend en compte les non linéarités dépendant de la vitesse, comme les phénomènes liés aux frottements ainsi que des incertitudes paramétriques sur les vecteurs de Coriolis et centrifuges. En référence aux relations (2) et (2'), on peut voir que les termes dépendant de la vitesse, dans l'expression du modèle dynamique du robot, sont dépendants des estimées utilisées pour la matrice des termes centrifuges et de Coriolis C(q,) (relation (2)) et pour les couples de frottement xf (relation (2')). Enfin la fonction de transfert S3(z) vise à modéliser les erreurs de modélisation liées aux termes correspondant au signe de la vitesse dans la relation (4) permettant ainsi de tenir compte des incertitudes sur les 1 0 frottements secs (voir relation (2')). La troisième composante du modèle du résidu rm(z), à savoir le terme Go(z)b(z) représente un vecteur de bruits blancs b(z) filtrés, qui permet de tenir compte du bruit de mesure et de dynamiques non modélisées comme 15 les flexibilités qui peuvent être présentes dans les robots sériels légers. De telles flexibilités sont caractérisées par des modes résonnants en basses fréquences. Dans le modèle proposé, selon l'invention, par les relations (5) et (6), le bruit blanc b(z) impactant le résidu est filtré par un filtre passe-bas Go(z). Le terme Go(z)b(z) équivalent à cette opération de filtrage correspond 20 en fait à un bruit dont le spectre fréquentiel est majoritairement situé dans les basses-fréquences. En effet, les bruits de mesure qui impactent le modèle dynamique d'évolution du robot sont généralement liés à des phénomènes basse-fréquence tels que des effets gravitationnels ou des flexibilités basses fréquences. Ces erreurs entrent dans la catégorie des incertitudes non 25 paramétriques évoquées précédemment. Au contraire, les phénomènes de collision se traduisent généralement par des fréquences couvrant un large spectre et plus le contact, lors de la collision, est caractérisé par une interaction de grande raideur avec l'environnement, plus le signal résultant comporte des composantes haute-fréquences. Un objectif du procédé selon 30 l'invention consiste donc à distinguer les composantes hautes fréquences, caractéristiques d'une collision, des composantes basses fréquences, caractéristiques de bruits de mesure. La figure 2 représente, plus en détail, les traitements mis en oeuvre 5 dans la seconde partie 200 du procédé. Dans la première étape 104 de la seconde partie 200 du procédé, le résidu r, obtenu à l'issu de la première partie 100 du procédé, est filtré par un filtre passe-haut Go-1(z) afin de rendre le résidu indépendant des incertitudes non paramétriques liées aux phénomènes basse-fréquence comme explicité 1 0 ci-dessus. La fonction de transfert du filtre passe-haut est, à titre illustratif, représentée par le terme Go-1(z) qui est l'inverse du filtre passe-bas Go(z) utilisé dans le modèle du résidu selon les relations (5) et (6). En pratique les coefficients du filtre passe-haut utilisé sont déterminés à l'aide d'un algorithme adaptatif, basé, par exemple, sur le modèle connu des moindres 15 carrés récursifs ou « recursive least squares » (RLS) en anglais. La fonction réalisée par un tel algorithme consiste à calculer les coefficients de la fonction de transfert visée, ici le filtre passe-haut, en tenant compte des mesures passées et du jeu de coefficients calculé à l'instant précédent. Le filtrage adaptatif du résidu est représenté à la figure 2 par la boite 20 201 qui correspond à l'apprentissage des coefficients du filtre et à la boite 202 qui correspond au filtrage lui-même réalisé sur le résidu r afin de produire en sortie un résidu filtré rf. Le résidu filtré obtenu est un vecteur dont les composantes sont les résidus filtrés respectifs correspondant à chaque articulation du robot, étant entendu que le robot comporte au moins une 25 articulation. Optionnellement, un module de filtrage 203 supplémentaire peut être ajouté. Ce module opère un filtrage temporel de chaque composante du vecteur des résidus filtré rf pour éliminer le bruit pouvant résulter d'une dérivation numérique. Le module de filtrage 203 peut, par exemple, être 30 réalisé par un calcul de moyenne quadratique ou « root mean square » en anglais. Alternativement, le module de filtrage 203 peut être remplacé par la fonction valeur absolue. On décrit à présent plus en détail un exemple de réalisation du filtrage 5 adaptatif du résidu r par le biais d'une méthode des moindres carrés récursifs (RLS). On rappelle tout d'abord les équations bien connues de mise en oeuvre de l'algorithme RLS. Soit ek le vecteur des paramètres à estimer, dans le cas présent il 10 s'agit du vecteur des coefficients du filtre G01 qui est appris en temps réel. Un paramètre de l'algorithme RLS est l'ordre des numérateurs et des dénominateurs de la fonction de transfert à estimer, ici G01(z) . Plus l'ordre de la fonction de transfert est grand, plus son pouvoir descriptif est grand et meilleure est la précision de représentation des résultats. En revanche, cela 15 se fait au détriment du temps de calcul qui peut devenir prohibitif. Le choix de l'ordre est donc un compromis entre précision du modèle et temps de calcul. Soit (1)k le vecteur des mesures d'entrées, dans le cas présent il s'agit des mesures successives passées du résidu r. 20 Les équations de l'algorithme RLS peuvent s'écrire de la façon suivante : Fk_iok-i(g_iFk-i) (7) Fk = 7lrk-1 Ek = Yk ek-11j k-1(8) ek = ek-1+ FkOk-1£k (9) 25 Fk est la matrice de covariance de la grandeur d'entrée de l'algorithme (ici le résidu r). À est le paramètre de facteur d'oubli. Yk est la grandeur d'entrée de l'algorithme, égale au résidu rk dans le cas présent. 30 Ek est l'erreur de prédiction qui correspond au critère à minimiser. Dans le cas présent, on cherche à minimiser la différence entre la valeur courante rk du résidu, produite en entrée de l'algorithme, et le résultat du filtrage, par la fonction de transfert aux coefficients estimés ek-1, des valeurs passées du résidu. L'indice k représente l'indice de la valeur courante d'une grandeur, les indices k-i représentent les indices des valeurs passées d'une grandeur. Pour cela, on pose Y k = rk Ok = [rk-i---rk-nGol T Le résidu filtré est calculé à l'aide de la relation suivante : rfilt = rk - goOk, où go désigne le vecteur transposé des coefficients estimés de la fonction de transfert du filtre G0-1 et (1)k le vecteur des valeurs passées du résidu rk. Les équations (7),(8),(9) données ci-dessus le sont à titre illustratif pour expliciter une mise en oeuvre possible de l'estimation des coefficients du filtre G0-1. Ces équations ne doivent pas être interprétées comme limitatives, étant entendu que l'Homme du métier sait implémenter tout autre algorithme adaptatif ou toute variante d'implémentation de l'algorithme adaptatif RLS pour arriver au même résultat visé par la présente invention. Dans la deuxième étape 105 de la seconde partie 200 du procédé, un seuil de détection de collision dynamique T est déterminé en ligne, c'est-à-dire au fur et à mesure du déplacement du robot. Le seuil dynamique T est calculé afin de prendre en compte les différentes erreurs de modélisation affectant le modèle dynamique du robot, 25 telles que décrites ci-dessus. Dans l'exemple de la figure 2, le modèle utilisé est du type de celui de la relation (1), autrement dit il est constitué essentiellement de composantes dépendantes de l'accélération articulaire (représentée par la grandeur ei(z)), de la vitesse articulaire (représentée par la grandeur e2(z)) et du signe de la vitesse articulaire (représenté par la grandeur e3(z)). En reprenant la relation (6) et en multipliant chaque membre par la fonction de transfert du filtre passe-haut G01(z), on obtient un modèle du résidu filtré : rf (z) = G01 (z)x, (z) + (z)ei (z) + 3,2 (Z)e2 (z) + 3,3 (Z)e3 (z) + b(z) (7) Avec Ai(z) = G0-1(z)Si(z) Selon le modèle de la relation (7), le résidu filtré est donc composé de trois éléments: - le vecteur des couples de collision filtré, qui représente le signal à détecter ; - les erreurs de modélisation dépendant des signaux ei(z) du système ; - un vecteur de bruits blancs. Un objectif de la deuxième étape 105 de la seconde partie 200 du procédé selon l'invention consiste à estimer en ligne les fonctions de transfert Ai(z), pour chaque variable d'entrée ei(z), et à filtrer cette variable 20 par un filtre reproduisant cette fonction de transfert. Pour cela, une méthode possible consiste à utiliser, tout comme pour le filtrage du résidu, un algorithme adaptatif du type moindres carrés récursifs (RLS). En reprenant les équations (7),(8) et (9) données ci-dessus, on applique 25 de la même façon que pour le filtrage du résidu, l'algorithme RLS avec les paramètres suivants pour chacune des variables d'entrée ei(z). La grandeur d'entrée de l'algorithme Yk est égale au résidu filtré rfk. ek est le vecteur des coefficients de la fonction de transfert Ai(z) à 30 estimer. (1)k est une matrice qui contient, pour chaque articulation du robot, un vecteur de la composante courante et des composantes passées (sur un horizon donné dépendant de l'ordre de la fonction de transfert choisi) de la variable d'entrée ei(z), par exemple l'accélération articulaire ou la vitesse articulaire. (1)k contient également les valeurs passées du seuil dynamique Toi. L'erreur de prédiction Ek que l'on cherche à minimiser est l'erreur entre la valeur courante du résidu filtré et le résultat du filtrage de la variable d'entrée ei par l'estimée de la fonction de transfert Ai(z). La valeur courante du seuil dynamique TAi(k) est calculée à l'aide de la 10 relation suivante Toi(k) = OT kOk, avec ek l'estimée des coefficients du filtre Ai(z) calculée à l'aide de la relation (9). L'algorithme adaptatif est donc appliqué pour chaque variable d'entrée ei(z) et consiste comme expliqué ci-dessus en une première partie 15 211,212,213 qui correspond à l'apprentissage des coefficients de la fonction de transfert du filtre Ai(z) et à une seconde partie 221,222,223 qui correspond au filtrage lui-même réalisé sur la variable d'entrée ei(z) afin d'obtenir une composante Tpi du seuil dynamique. En sortie de l'algorithme adaptatif, la composante Tpi est une grandeur 20 vectorielle qui contient les composantes associées aux différentes articulations du robot. Comme déjà explicité ci-dessus, optionnellement, un module de filtrage temporel supplémentaire 231,232,233 peut être ajouté pour éliminer des bruits résiduels dans les composantes Toi. Ces modules de filtrage 25 231,232,233 peuvent, par exemple, être réalisés par un calcul de moyenne quadratique ou « root mean square » en anglais. Alternativement, ils peuvent être remplacés par la fonction valeur absolue. In fine, le seuil dynamique T global est obtenu en sommant les 30 composantes Tpi obtenues pour chaque variable d'entrée ei(z) et en ajoutant une composante statique Tstatic dont la valeur est réglée en fonction du niveau de bruit. Dans l'exemple de la figure 2 qui porte sur trois variables d'entrée, le seuil dynamique global obtenu par application de l'invention est donné par la 5 relation suivante : T = Tstatic + Tpl + TA2 TA3 La dernière étape 106 du procédé selon l'invention consiste à comparer le résidu filtré rf en sortie de l'étape 104 avec le seuil dynamique global T et en déduire l'existence d'une collision si le résidu filtré dépasse ledit seuil. 10 Cette comparaison peut être effectuée pour chaque composante du seuil T correspondant à chaque articulation du robot. Une collision peut être détectée si le résidu filtré dépasse le seuil dynamique T pour au moins une articulation. D'autres variantes de réalisation de l'étape 106 finale de détection sont possibles comme le déclenchement d'une alerte de détection 15 seulement si un nombre K prédéterminé parmi les N articulations du robot considérés satisfait au critère de comparaison entre le résidu filtré et le seuil dynamique T. Toute autre variante de réalisation permettant d'élaborer une décision sur l'occurrence d'une collision en fonction de la comparaison du vecteur des résidus filtrés et du vecteur des seuils dynamiques T est 20 envisageable et sera comprise par l'Homme du métier comme faisant partie à part entière de l'invention. La figure 3 illustre, sur un diagramme représentant l'évolution du résidu filtré rf obtenu par application du procédé selon l'invention en fonction du 25 temps, l'apport de l'invention par rapport aux solutions de l'art antérieur. Sur la figure 3 est représentée l'évolution 301 du résidu filtré en fonction du temps. Pendant une première période de temps 320, le résidu filtré subit une première variation d'amplitude liée aux erreurs de modélisation sur le modèle dynamique du robot. Pendant une seconde période de temps 321, le 30 résidu filtré subit une seconde variation d'amplitude liée à une collision avec un objet. Dans le cas des solutions de l'art antérieur qui utilisent un seuil de détection statique, on remarque qu'il n'est jamais possible, avec un même seuil, à la fois de détecter la collision tout en évitant de déclencher une alarme pour des phénomènes liés aux erreurs de modélisation. Si le seuil 310 est trop élevé, il ne permet pas de détecter la collision 321, si le seuil 311 est trop bas, il engendre une fausse détection sur l'évènement 320 lié aux erreurs de modélisation. Sur la figure 3 est également représenté le seuil dynamique 312 obtenu par application du procédé selon l'invention et qui permet à la fois de correctement détecter le phénomène de collision 321 tout en évitant une fausse détection sur le premier évènement 320 source de variation d'amplitude du résidu. La figure 4 représente un schéma d'un système de contrôle 400 selon 15 l'invention d'un robot manipulateur 500 comprenant une ou plusieurs articulations A,B. Le système 400 selon l'invention comporte un organe de commande 401 de la trajectoire du bras du robot manipulateur 500, une interface physique 402 d'entrée/sortie entre le système 400 et le robot 500 et un module 403 de 20 détection des collisions. L'organe de commande 401 reçoit en entrée un ordre de commande haut-niveau 410 pour piloter le robot 500 et transmet en sortie 411, via l'interface 402, des informations au robot 500 pour le manipuler, par exemple un couple de référence. En retour, une mesure 412 des positions et des 25 couples articulaires du robot 500 est fournie au système 400. Cette mesure 412 peut être directement prise égale au couple de référence 411. En effet, dans le cas de moteurs à courants continus, un couple de référence est envoyé aux amplificateurs de puissance pilotant les moteurs du robot et asservis en courant par une boucle de rétroaction interne indépendante du 30 contrôleur principal. Cet asservissement en courant étant très rapide à l'échelle de la dynamique du robot et du temps d'échantillonnage principal, on peut en général supposer que le couple effectivement appliqué au robot est égal au couple de référence. Grâce à cette hypothèse il est donc possible d'utiliser le couple de référence directement comme entrée à l'algorithme de détection, sans avoir besoin de mesurer explicitement le couple articulaire. Alternativement, si le couple de référence 411 n'est pas disponible, par exemple parce que l'organe de commande 401 et le module de détection de collisions 403 sont implémentés dans deux équipements distincts, une mesure 412 des couples articulaires du robot est nécessaire. Le module de détection de collisions 403 reçoit la mesure 412 des 10 couples articulaires du robot et établi une décision 413 sur la présence ou l'absence de collision. Cette décision 413 est par exemple une décision binaire qui est ensuite fournie à une interface 414 pour être exploitée. Le module de détection des collisions 403 peut être implémenté à partir d'éléments matériel et/ou logiciel. Il peut notamment être mis en oeuvre 15 en tant que programme d'ordinateur comportant des instructions pour son exécution. Le programme d'ordinateur peut être enregistré sur un support d'enregistrement lisible par un processeur. Références 20 [1] V. Venkatasubramanian, R. Rengaswamy, K. Yin, and S. N. Kavuri, "A review of process fault detection and diagnosis: Part i: Quantitative modelbased methods," Computers Chemical Engineering, pp. 293-311,2003. [2] P. Frank and X. Ding, "Survey of robust residual generation and 25 evaluation methods in observer-based fault detection systems," Journal of process control, vol. 7, no. 6, pp. 403-424,1997. [3] H. Sneider and P. Frank, "Observer-based supervision and fault detection in robots using nonlinear and fuzzy logic residual evaluation," Control Systems Technology, IEEE Transactions on, vol. 4, no. 3, pp.274-282, 1996. 30 [4] A. De Luca and R. Mattone, "Actuator failure detection and isolation using generalized momenta," in Robotics and Automation, 2003. Proceedings. ICRA'03. IEEE International Conference on, vol. 1. IEEE, 2003, pp. 634-639. [5] S. Haddadin, A. Albu-Schaffer, A. De Luca, and G. Hirzinger, "Collision detection and reaction: A contribution to safe physical humanrobot 5 interaction," in Intelligent Robots and Systems, 2008. IROS 2008.IEEE/RSJ International Conference on. IEEE, 2008, pp. 3356-3363. [6] W. Dixon, I. Walker, D. Dawson, and J. Hartranft, "Fault detection for robot manipulators with parametric uncertainty: a prediction-errorbased approach," Robotics and Automation, IEEE Transactions on, vol. 16, no. 6, 10 pp. 689-699, 2000. [7] A. De Luca and R. Mattone, "An adapt-and-detect actuator fdi scheme for robot manipulators," in Proceedings. ICRA'04. 2004 IEEE International Conference on, vol. 5. IEEE, 2004, pp. 4975-4980. 15 REVENDICATIONS1. Procédé de détection de collision entre un robot composé d'une pluralité de corps reliés entre eux par au moins une articulation et son environnement ledit procédé comprenant les étapes suivantes : - Générer (100), à partir d'un modèle dynamique du robot, un signal représentatif des collisions entre le robot et son environnement, ledit signal étant appelé résidu r et comprenant autant de composantes que d'articulations du robot, ledit procédé étant caractérisé en ce qu'il comporte en outre les étapes suivantes : - Réaliser (104) un filtrage passe-haut adaptatif du résidu r pour rendre le résidu indépendant d'incertitudes paramétriques ou non paramétriques liées à des phénomènes basses-fréquences, - Déterminer (105) un seuil adaptatif T composé d'au moins un premier terme Tot, dynamique, égal à une information sur les incertitudes paramétriques entre ledit modèle et le comportement réel du robot, lesdites incertitudes paramétriques étant liées à une première variable ei caractéristique de l'état dudit robot parmi les variables suivantes : la position, la vitesse ou l'accélération d'une articulation du robot, ou une fonction linéaire ou non linéaire de l'une de ces variables ou d'une combinaison de ces variables, - Comparer (106) le résidu filtré au seuil adaptatif T pour en déduire l'existence ou non d'une collision. 2. Procédé de détection de collision selon la revendication 1 dans lequel ledit seuil adaptatif T est composé de la somme de plusieurs termes Tot, Tot, Toi, dynamiques, chacun égaux à une information sur les incertitudes paramétriques entre ledit modèle et le comportement réel du robot, lesdites incertitudes paramétriques étant liées à une variable différente, caractéristique de l'état dudit robot parmi les variables suivantes : la 3002 04 8 25 position, la vitesse ou l'accélération d'un point fixe du robot, ou une fonction linéaire ou non linéaire de l'une de ces variables ou d'une combinaison de ces variables, 5 3. Procédé de détection de collision selon l'une des revendications précédentes dans lequel ledit seuil adaptatif T comporte en outre un terme Tstatic statique configuré de sorte à être supérieur à un niveau de bruit de mesure. 10 4. Procédé de détection de collision selon l'une des revendications précédentes dans lequel l'étape de filtrage (104) est réalisée au moyen d'un algorithme des moindres carrés récursifs. 5. Procédé de détection de collision selon la revendication 4 dans lequel 15 l'étape de filtrage (104) comporte les sous-étapes suivantes : - Estimer (201), de façon récursive, les coefficients de la fonction de transfert du filtre passe-haut G0-1, - Filtrer (202) le résidu r avec le filtre passe-haut estimé. 20 6. Procédé de détection de collision selon l'une des revendications 4 ou 5 dans lequel une étape supplémentaire de filtrage temporel (203) est appliquée au résidu filtré pour chacune de ses composantes. 7. Procédé de détection de collision selon la revendication 6 dans lequel 25 l'étape supplémentaire de filtrage temporel (203) est une étape de calcul de moyenne quadratique (RMS). 8. Procédé de détection de collision selon l'une des revendications précédentes dans lequel l'étape de détermination (105) d'au moins un 30 premier terme Tot, Tot, Toi, dynamique du seuil adaptatif T est réalisée au moyen d'un algorithme des moindres carrés récursifs. 9. Procédé de détection de collision selon la revendication 8 dans lequel l'étape de détermination (105) d'au moins un premier terme Tot, T42, T435 dynamique du seuil adaptatif T comporte les sous-étapes suivantes : - Estimer (211,212,213), de façon récursive, les coefficients de la fonction de transfert A modélisant les incertitudes paramétriques liées à ladite première variable ei caractéristique de l'état dudit robot, - Filtrer (221,222,223) ladite première variable ei caractéristique de l'état dudit robot avec un filtre de fonction de transfert A estimée à l'étape précédente (211,212,213) pour obtenir un terme Tm, Tm, Toi, dynamique du seuil adaptatif T. 10. Procédé de détection de collision selon l'une des revendications 8 ou 9 dans lequel une étape supplémentaire de filtrage temporel (231,232,233) est appliquée audit terme T3,1, Tm, Toi, dynamique du seuil adaptatif T pour chacune de ses composantes. 11. Procédé de détection de collision selon la revendication 10 dans lequel l'étape supplémentaire de filtrage temporel (231,232,233) est une étape de calcul de moyenne quadratique (RMS). 12. Procédé de détection de collision selon l'une des revendications précédentes dans lequel l'étape (106) de comparaison du résidu filtré au seuil adaptatif T pour en déduire l'existence ou non d'une collision comprend les sous-étapes suivantes : - Comparer, pour chaque composante articulaire, le résidu filtré au seuil adaptatif T, - Conclure à l'existence d'une collision si, pour au moins K composantes, le résidu filtré est supérieur au seuil adaptatif T, K étant un nombre entier prédéterminé strictement positif et inférieur ou égal au nombre d'articulations du robot. 13. Procédé de détection de collision selon l'une des revendications précédentes dans lequel l'étape (100) de génération du résidu est constituée des sous-étapes suivantes : - Déterminer (101), à partir d'une information sur l'état du robot et par le biais d'un modèle dynamique, une estimation des couples articulaires du robot, - Effectuer une mesure (102) de l'état du robot, par exemple une mesure des couples articulaires' - Calculer (103) le résidu comme la différence entre l'estimation et la mesure de l'état du robot. 14. Procédé de détection de collision selon l'une des revendications précédentes dans lequel les incertitudes paramétriques liées à l'accélération articulaire du robot sont des incertitudes sur la matrice d'inertie du robot. 15. Procédé de détection de collision selon l'une des revendications précédentes dans lequel les incertitudes paramétriques liées à la vitesse articulaire du robot sont des incertitudes sur la matrice des termes centrifuges et de Coriolis du robot et/ou sur les frottements visqueux. 16. Procédé de détection de collision selon l'une des revendications précédentes dans lequel une fonction non linéaire est la fonction signe ou la fonction exponentielle ou la fonction valeur absolue. 17. Procédé de détection de collision selon la revendication 16 dans lequel les incertitudes paramétriques liées au signe de la vitesse articulaire du robot sont des incertitudes sur les frottements secs. 18. Programme d'ordinateur comportant des instructions pour l'exécution du procédé de détection de collision selon l'une quelconque desrevendications 1 à 17, lorsque le programme est exécuté par un processeur. 19.Support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution du procédé de détection de collision selon l'une quelconque des revendications 1 à 17, lorsque le programme est exécuté par un processeur. 20.Système (400) de contrôle d'un robot (500) comprenant un organe de commande (401) pour la manipulation du robot (500), une interface (402) pour l'échange d'informations sur l'état du robot (500) et un module (403) de détection de collisions adapté pour mettre en oeuvre le procédé selon l'une des revendications 1 à 17.15