FR 3090913 A1 20200626 FR 1874139 A 20181224 Description Titre de l'invention : METHODE DE CARACTERISATION D'UN SIMULATEUR DE MOUVEMENTS, PROCEDE DE SUR- VEILLANCE ET DISPOSITIF CORRESPONDANT Domaine technique 100011 La présente invention concerne de manière générale le domaine des systèmes électro- mécaniques/mécatroniques et plus particulièrement des appareils de test dans ce domaine. Elle concerne plus spécifiquement une méthode de caractérisation d'un simulateur de mouvements, un procédé de surveillance et un dispositif correspondant permettant d'effectuer une surveillance du fonctionnement dudit simulateur de mouvements ainsi que de la prévention et des diagnostics. Arrière-plan technologique [0002] Les simulateurs de mouvements sont des équipements destinés à tester des composants inertiels, en particulier les gyromètres et les centrales inertielles. Ils sont le plus souvent constitués d'un ou de plusieurs axes rotatifs inclus dans une géométrie de type cardan. Il existe également. des simulateurs comprenant des axes en translation. Les actionneurs de ces équipements sont généralement des moteurs électriques, parfois hydrauliques. Les vitesses et/ou les positions angulaires de chaque axe sont mesurées par des capteurs ad-hoc, par exemple les positions angulaires peuvent être mesurées par des codeurs. [0003] Ces simulateurs sont largement utilisés dans des procédés de production dans le cadre de tests systématiques. Dans ce cadre, la fiabilité du simulateur de mouvements se révèle de la plus haute importance, attendu qu'une défaillance totale du moyen de test est susceptible de bloquer complétement une chaire de production, ce qui est très coûteux pour l'industriel ou est susceptible de donner de fausses informations sur les produits testés. [0004] Afin d'éviter une défaillance inattendue du simulateur de mouvements, des opérations de maintenance préventive sont réalisées périodiquement. Néanmoins, il manque à l'heure actuelle un dispositif automatique permettant de diagnostiquer les axes des simulateurs de mouvements, et en particulier il manque un moyen d'estimer précisément un certain nombre de paramètres physiques de l'équipement, en vue de détecter une évolution anormale desdits paramètres. Il est certain que la détection de l'évolution anormale des paramètres physiques du simulateur avant qu'elle ne devienne critique, permettrait d'entreprendre unc mesure de réparation anticipée de l'équipement, avant que le dysfonctionnement naissant empêche de se servir dudit simulateur de mouvements. Un tel dispositif permettrait ainsi une très significative 2 réduction des risques de défaillance du moyen de test. Exposé de l'invention [0005] La présente invention vise à pallier à ce manque et consiste notamment à équiper le simulateur de mouvements d'un moyen de diagnostic et d'estimation de ses paramètres physiques. [0006] A cette fin, l'invention concerne une méthode de caractérisation d'un simulateur de mouvements dans laquelle on estime des paramètres liés au fonctionnement dudit simulateur de mouvements, le simulateur de mouvements comportant un ensemble d'axes, chaque axe étant entraîné en rotation ou translation par un organe effecteur recevant un signal de commande U(t) et comportant au moins un capteur de mesure de la rotation ou de la translation dudit axe permettant de mesurer au moins la position t ) angulaire ou À ( t) linéaire dudit axe, chaque organe effecteur d'axe étant dans unc boude de correction fermée reliée par une première entrée, au dit au moins un capteur de mesure et produisant, en sortie, le signal de commande U(t), ladite boucle de correction comportant unc seconde entrée de signal destinée à recevoir un signal de consigne pour commander ledit axe, [0007] dans laquelle méthode, les paramètres estimés sont au moins un paramètre sé- lectionné parmi : [0008] - un paramètre (ri a) dépendant du moment d'inertie sur l'axe et dont l'effet sur le fonctionnement est fonction d'une variable d'accélération angulaire (t) ou linéaire Â (t ) dudit axe, [0009] - un paramètre (17 v) dépendant du coefficient de frottements visqueux sur l'axe et dont l'effet sur le fonctionnement est fonction d'une variable de vitesse angulaire é ( t ) ou linéaire À ( t ) dudit axe, [0010] - un paramètre (17 s) dépendant du coefficient de frottements secs sur l'axe et dont l'effet sur le fonctionnement est fonction d'une variable de vitesse angulaire e (,) ou linéaire À ( t ) dudit axe, [0011] - au moins un paramètre () ) (r7 g fis) ) (t) ) dépendant du coefficient de balourd c (couples de gravité) et dont l'effet sur le fonctionnement est fonction d'une variable de position angulaire et ( t ) ou linéaire À ( t ) dudit axe, et 100121 dans laquelle méthode, on effectue une identification dynamique de modèle inverse (1MM) par la méthode de la variable instrumentale dans laquelle on détermine une matrice instrumentale Z, ladite identification dynamique de modèle inverse (1DIM) s'effectuant sous une excitation de l'axe par un signal d'excitation envoyé dans la seconde entrée de signal, 100131 dans laquelle méthode, on acquiert : 3 [0014] - le signal de commande U(t), [0015] - des mesures dudit au moins un capteur de mesure et on calcule des fonctions des mesures acquises, [0016] on filtre le signal de commande U(t) acquis, une/des mesures acquises et des fonctions des mesures acquises afin d'obtenir des données correspondantes acquises filtrées. [0017] Selon l'invention, le signal d'excitation est produit par un générateur de consigne calculant une consigne de position e (t) ou À Jo, une consigne de vitesse é Jr) ou À À t) et une consigne d'accélération é À t ) ou À À t) qui correspondent au signal d'excitation, et pendant l'excitation on acquiert en outre : 100181 - la consigne de vitesse () (t) ou À ,(t) et la consigne d'accélération e (t) ou À c.( t ) du générateur de consigne, et [0019] on filtre en outre les consignes de vitesse et d'accélération acquises afin d'obtenir des données correspondantes acquises filtrées, [0020] on calcule la matrice instrumentale Z en utilisant des variables, lesdites variables étant : [0021] - la consigne d'accélération acquise filtrée ef t r ou À t ) pour ce qui c concerne le paramètre (r) (1) dépendant du moment d'inertie sur l'axe, [0022] - la consigne de vitesse acquise filtrée ( t ) ou À cf( t ) pour ce qui concerne c le paramètre (17 ,.') dépendant du coefficient de frottements visqueux sur l'axe, [0023] - une des fonctions des mesures acquises filtrées pour ce qui concerne le paramètre ) dépendant du coefficient de frottements secs sur l'axe, [0024] - au moins une des fonctions des mesures acquises filtrées pour ce qui concerne ledit au moins un paramètre (Ti ) c s ) (11 ) dépendant du coefficient de balourd, et [0025] on estime les paramètres par une méthode de calcul mettant en oeuvre la matrice ins- trumentale Z. [0026] D'autres caractéristiques non limitatives et avantageuses de la méthode, prises indivi- duellement ou selon toutes les combinaisons techniquement possibles, sont les suivantes : [0027] - la boucle de correction fermée est en rétroaction, [0028] - le paramètre dépendant du coefficient de balourd est omis de la méthode dans le cas où l'axe entraîné en rotation est vertical ou dans le cas où l'axe est entraîné en translation, [0029] - les fonctions des mesures acquises filtrées concernant ledit au moins un paramètre ( ) 4_) (17 ) (n .) dépendant du coefficient de balourd sont prises nulles dans le cas où l'axe en rotation est vertical ou dans le cas où l'axe est en translation, [0030] - chaque exécution de la méthode s'effectue pour un axe donné du simulateur de mouvements, [0031] - on exécute la méthode pour chaque axe devant être caractérisé parmi les axes du si- mulateur de mouvements, [0032] - les frottements secs correspondent à un modèle dit de Coulomb, [0033] - le simulateur de mouvements comporte un générateur de consigne, [0034] - le signal de consigne et le signal d'excitation sont des signaux de position angulaire ou linéaire, [0035] - le signal de consigne est un signal de vitesse angulaire ou linéaire, [0036] - le capteur de mesure de la rotation de l'axe permet de mesurer la position angulaire 0(t) dudit axe, [0037] - le capteur de mesure de la translation de l'axe permet de mesurer la position linéaire M') dudit axe, [0038] - le capteur de mesure de la rotation de l'axe permet de mesurer la position 0(0 et la vitesse é (t ) angulaires dudit axe, [0039] - le capteur de mesure de la translation de l'axe permet de mesurer la position À0) et la vitesse À ( t ) linéaires dudit axe, 100401 - le capteur de mesure de la rotation de l'axe permet en outre de mesurer l'accélération angulaire 6 (t) ou linéaire À (t) 100411 - pour le coefficient de balourd, on estime deux paramètres en quadrature, un cor- respondant à un sinus (ri Lis c ) et l'autre au cosinus (1) ) de la position angulaire de [0042] - le capteur de mesure de l'axe permet de mesurer la position angulaire 0(t) ou linéaire À ( t ) dudit axe, et pour le coefficient de balourd, on estime deux paramètres en quadrature, un correspondant à un sinus (q s) et l'autre au cosinus (ti )(le la position angulaire 0(t), et en ce que les fonctions des mesures dudit. capteur de mesure sont choisies parmi au moins : [0043] - un calcul de la vitesse angulaire Q ( t ) ou linéaire À ( t) à partir de la position angulaire 0(t) ou linéaire À ( t ) mesurée, [0044] - un calcul de l'accélération angulaire Q (t ) ou linéaire À (t ) à partir de la position angulaire 0(t) ou linéaire À ( t ) mesurée, [0045] - un calcul de détermination du signe de la vitesse angulaire ou linéaire, cette dernière étant obtenue par un calcul de la vitesse à partir de la position angulaire 0(t) ou linéaire À ( t ) mesurée, [0046] - un calcul de cosinus de la position angulaire 0(t) mesurée, [0047] - un calcul de sinus de la position angulaire 0(t) mesurée, [0048] - le capteur de mesure de l'axe permet de mesurer en outre la vitesse angulaire ( t ) ou linéaire À ( t ) dudit axe, et pour le coefficient de balourd, on estime deux paramètres en quadrature, un correspondant à un sinus (17 Ys) et l'autre au cosinus ( c) de la position angulaire 00), et les fonctions des mesures dudit capteur de g mesure sont choisies parmi au moins : [0049] - un calcul de l'accélération angulaire 0 (t ) ou linéaire A (t ) à partir de la vitesse angulaire é ( t ) ou linéaire À ( t ) mesurée, [0050] - un calcul de détermination du signe de la vitesse angulaire ou linéaire à partir de la vitesse angulaire é ( t ) ou linéaire À ( t ) mesurée, [0051] - la vitesse angulaire é ( t ) ou linéaire À ( t ) est calculée par la dérivée première de la position angulaire 0(t) ou linéaire ii ( t ) mesurée, [0052] - l'accélération angulaire 0 (t ) ou linéaire A ( t ) est calcifiée par la dérivée seconde de la position angulaire 0(t) ou linéaire À ( t ) mesurée, 100531 - l'accélération angulaire a(,) ou linéaire A ( t ) est calculée par la dérivée première de la vitesse angulaire é ( t ) ou linéaire À ( t ) mesurée, [0054] - clans le cas où le paramètre (1-1 s) dépendant du coefficient de frottements secs sur l'axe doit être estimé, alors, pendant l'excitation on met en oeuvre un signal d'excitation qui ne correspond qu'à un seul sens de rotation de l'axe, [0055] - dans le cas où le paramètre (17 s) dépendant du coefficient de frottements secs sur l'axe doit être estimé, alors, pendant l'excitation on met en oeuvre un signal d'excitation qui ne correspond qu'à un seul sens de rotation de l'axe, et on assigne d'office une valeur constante à la variable correspondant à la fonction des mesures acquises filtrées de calcul de détermination du signe de la vitesse de l'axe, [0056] - les mesures sont échantillonnées et numérisées, les calculs des fonctions des mesures sont numériques et le filtrage est un filtrage numérique par un filtre numérique et pour calculer la matrice instrumentale Z on détermine un vecteur ligne j ( t ) regroupant les variables, ledit vecteur ligne jz ( t ) étant [0057] lz(t)- [Oc (t) écf(t) signt,(é(o) sef(t) cei(t)1 dans le cas d'un axe non vertical entraîné en rotation, ou [0058] Ut) = rocin) sem) signr(ew) dans le cas d'un axe vertical entraîné en rotation, ou [0059] (t) [À ci( t ) À c,(t ) sign t(A (t)) dans le cas d'un axe entraîné en translation, où 19cf f (t)À (0 correspondent à la consigne 6 J. (t À cf ( ) correspondent à la consigne de vitesse acquise filtrée, sugn/ (e (t))sign (À (t)) correspondent au filtrage du calcul de détermination du signe de la vitesse angulaire ou linéaire respectivement, ce r(t ) correspond au filtrage du calcul de cosinus de la position angulaire et se r(t) correspond au filtrage du calcul de sinus de la position angulaire, la matrice Z étant : [0060] pour t = 1 N, où N est le nombre d'échantillons, /Z() 100611 - la méthode de calcul mettant en oeuvre la matrice instrumentale Z utilisée pour estimer les paramètres correspond à la résolution de z I z y, où la matrice X est la concaténation des vecteurs-lignes bk- ( t )l, [0062] X avec 1y(N) [0063] ix (t) [91. t) e f,(11 sign t,(ë (t)) ) ce 1,(t)] dans le cas d'un axe non vertical entraîné en rotation, ou [0064] = [eif(t) 6 f(t) sign 1,(a (t)) dans le cas d'un axe vertical entraîné en rotation, ou [0065] ix (t) = [À r(t)À r(t) sign t(À (t)) dans le cas d'un axe entraîné en translation, et [0066] où 0"f( t )À"(t) correspondent au filtrage de l'accélération. 9 f ( At) cor- respondent au filtrage de la vitesse, sign f(0 )sigrz f(À (t)) cor- respondent au filtrage du calcul de détermination du signe de la vitesse angulaire ou linéaire respectivement, cet.(t ) correspond au filtrage du calcul de cosinus de la position angulaire et se f(t) correspond au filtrage du calcul de sinus de la position angulaire, et [0067] 7 ,(1) pour t = 1 -N , où -N est le nombre d'échantillons, Ul'(2) U (N) [0068] - clans le cas où le capteur de mesure de la rotation de l'axe permet de mesurer la position angulaire 0(t) dudit axe, alors dans le vecteur ligne jz(t) et dans le vecteur ligne j ( t ) , é « ( t ) correspond à la consigne d'accélération acquise filtrée, éc ( t ) correspond à la consigne de vitesse acquise filtrée, eir(t) correspond au filtrage du résultat du calcul de l'accélération de l'axe ë (t) à partir de la position angulaire 0(t) mesurée, e f(t) correspond au filtrage du résultat du calcul de la vitesse angulaire de l'axe 9 ( t ) à partir de la position angulaire 0(t) mesurée, sign (0(0) correspond au filtrage du calcul de détermination du signe de la vitesse angulaire de l'axe, cette dernière é ( t ) étant obtenue par un calcul de la vitesse angulaire de l'axe à partir de la position angulaire 0(t) mesurée, c ef' ( t ) correspond au filtrage du calcul de cosinus de la position angulaire 0(t) mesurée et S e ,.(t) correspond au filtrage du calcul de sinus de la position angulaire 0(0 mesurée, [0069] - dans le cas où le capteur de mesure de la rotation de l'axe permet de mesurer la vitesse angulaire é ( t ) dudit axe, alors dans le vecteur ligne 17 ( t ) et dans le vecteur ligne j y ( t ) , é cf (T ) correspond à la consigne d'accélération acquise filtrée_- (t) correspond à la consigne de vitesse acquise filtrée, 0.i.(t) correspond au filtrage du résultat du calcul de l'accélération de l'axe O (t) à partir de la vitesse angulaire é ( t ) mesurée, 6 f(t) correspond au filtrage de la vitesse angulaire e ) mesurée, sign (o (t ) ) correspond au filtrage du calcul de détermination du signe de la vitesse angulaire mesurée é (t ),ce j.( t ) correspond au filtrage du calcul de cosinus de la position angulaire O(t), cette dernière HO-) étant obtenue par un calcul de la position angulaire de l'axe à partir de la vitesse angulaire ( t ) mesurée, et s e f(t) correspond au filtrage du calcul de sinus de la position angulaire 0(t), cette dernière 0(0 étant obtenue par un calcul de la position angulaire de l'axe à partir de la vitesse angulaire e t ) mesurée, [0070] - les paramètres s'appliquent à un modèle inverse d'un axe de simulateur de mouvements, ledit modèle inverse étant : 8 [0071] jr fl, f . f Sin (e + go ) dans le K Ll M e + r 0 + k.- s z g n ( 0 ) + J K K a m - a ni -11- ei ni - if cas d'un axe non vertical entraîné en rotation, ou [0072] f - f . fs -- dans le cas d'un axe vertical K'a ni 0 + K a inv 0 +K s 1 g n (0) = U CI in entraîné en rotation, ou [0073] I f fssign (A ) = udans le cas d'un axe entraîné K- À + K JimÀ + a In a in en translation, [0074] où 6 est la position angulaire de l'axe, 0 la vitesse angulaire de l'axe, ) la vitesse linéaire de l'axe, 0 l'accélération angulaire de l'axe et À l'accélération linéaire de l'axe, K {à,,, est un scalaire, est le moment d'inertie, f est un coefficient de frottements visqueux, f co- s est un coefficient de frottements secs, f est un co- g efficient de balourd, 9 est un angle en relation avec la position du centre de gravité par rapport à la verticale et les paramètres sont au moins un de : [0075] qa-Kon, [0076] f. ri Kam [0077] f' Kani [0078] g - KtI Tri [0079] rig e et ri '1',, vérifiant 17 A 2 q 2 et rÏ gs atan2 g'-, [0080] - K a in est un scalaire égal à la constante de couple de l'organe effecteur, notamment un moteur, entraînant l'axe en rotation par le gain de la chaine de commande de l'organe effecteur, notamment le gain de l'amplificateur de puissance actionnant ledit moteur, [0081] - dans le cas où l'axe entraîné en rotation est vertical, on assigne la valeur 0 à /7 soit n = n = n = o, 0 il c g s [0082] - dans le cas où l'axe est entraîné en translation, on assigne la valeur 0 à n , soit q = Cl C n a = 0, ,i ri s , [0083] - le filtrage est effectué par un filtre passe bande ou passe bas, 9 [0084] - le filtrage est analogique, [0085] - les calculs des fonctions des mesures sont analogiques, [0086] - le filtrage est un filtrage numérique, [0087] - le filtre numérique a une fonction de transfert correspondant à la discrétisation de la fonction de transfert continue 1 avec 7 une constante de temps de ( T S + 1 ) 2 réglage, [0088] - le filtre numérique est choisi parmi les filtres stables passe bas, [0089] - le simulateur de mouvements comporte un mode maintenance et la méthode est mise en oeuvre lorsque le simulateur de mouvements est dans le mode maintenance, [0090] - le simulateur de mouvements comporte une interface utilisateur avec un dispositif d'affichage et on affiche sur le dispositif d'affichage les paramètres estimés obtenus lors de l'exécution de la méthode, [0091] - le simulateur de mouvements comporte une interface utilisateur avec un dispositif d'affichage et une mémoire mémorisant les paramètres estimés obtenus lors d'exécutions successives au cours du temps de la méthode et on affiche sur le dispositif d'affichage les évolutions des paramètres estimés obtenus, [0092] - le simulateur de mouvements comporte un des moyens de calculs statistiques et de tests comparatifs sur les paramètres estimés. [0093] L'invention concerne également un procédé de surveillance du fonctionnement d'un simulateur de mouvements comportant un ensemble d'axes, chaque axe étant entraîné en rotation ou translation par un organe effecteur recevant un signal de commande U(t) et comportant au moins un capteur de mesure de la rotation ou de la translation dudit axe permettant de mesurer au moins la position 9 ( t ) angulaire ou À ( t ) linéaire dudit axe, chaque organe effecteur d'axe étant dans une boucle de correction fermée reliée par une première entrée, au dit au moins un capteur de mesure et produisant, en sortie, le signal de commande U(t), ladite boucle de correction comportant une seconde entrée de signal destinée à recevoir un signal de consigne pour commander ledit axe, dans lequel procédé, on exécute d'une manière répétée au cours du temps la méthode de caractérisation de l'invention permettant d'obtenir des estimations des paramètres liés au fonctionnement dudit simulateur de mouvements et on compare entre eux les paramètres estimés obtenus lors des exécutions successives de la méthode et/ou on compare les paramètres estimés obtenus lors d'une exécution de la méthode avec des valeurs nominales de paramètres. 100941 D'autres caractéristiques non limitatives et avantageuses du procédé, prises indivi- duellement ou selon toutes les combinaisons techniquement possibles, sont les suivantes : 100951 - les valeurs nominales de paramètres sont des valeurs des estimations des paramètres obtenues par la méthode de l'invention sur un simulateur de mouvements ayant un fonctionnement nominal, de préférence le même simulateur de mouvements, [0096] - le simulateur de mouvements comporte des moyens pour l'exécution de la méthode et du procédé et on mémorise dans une mémoire du simulateur de mouvements les paramètres estimés obtenus lors des exécutions successives de la méthode, [0097] - on raccorde au simulateur de mouvements un équipement distinct comportant des moyens pour l'exécution de la méthode et du procédé et on mémorise dans une mémoire de l'équipement distinct les paramètres estimés obtenus lors des exécutions successives de la méthode, [0098] - on effectue des calculs statistiques sur les paramètres estimés mémorisés obtenus lors des exécutions successives de la méthode, [0099] - les calculs statistiques sont choisis parmi des moyennes, des régressions, des ajustements de courbes... [0100] - les exécutions successives de la méthode sont séparées entre elles d'au moins une heure, [0101] - les exécutions successives de la méthode sont séparées entre elles d'au moins une journée, [0102] - les exécutions successives de la méthode sont séparées entre elles d'au moins une semaine, [0103] - à chaque exécution de la méthode de caractérisation de l'invention, on estime le même/les mêmes paramètres d'au moins un même des axes et on compare pour ledit ou pour chacun des paramètres et pour ledit ou pour chaque axe, les paramètres estimés obtenus lors des exécutions successives de la méthode, [0104] - pour un même axe, on compare entre elles les estimations du paramètre n obtenues lors des exécutions successives de la méthode, [0105] - pour un même axe, on compare entre elles les estimations du paramètre riv obtenues lors des exécutions successives de la méthode, [0106] - pour un même axe, on compare entre elles les estimations du paramètre 17 obtenues lors des exécutions successives de la méthode, [0107] - pour un même axe, on compare entre elles les estimations du paramètre n g obtenues lors des exécutions successives de la méthode, 101081 - pour un même axe, on compare entre elles les estimations du paramètre I7 g c obtenues lors des exécutions successives de la méthode, [0109] - pour un même axe, on compare entre elles les estimations du paramètre ri ys obtenues lors des exécutions successives de la méthode, [0110] - pour un même axe, on compare entre elles des mêmes combinaisons linéaires de 11 deux ou plusieurs des estimations des paramètres 17 a v, obtenues lors des exécutions successives de la méthode, 10111] - lors de la comparaison, on calcule une différence entre paramètres estimés de deux exécutions de la méthode et on compare par rapport à un seuil ou à une fourchette de valeurs, ladite différence, 10112] - lors de la comparaison, on calcule une différence entre paramètres estimés de deux exécutions successives de la méthode et on compare par rapport à un seuil ou à une fourchette de valeurs, la valeur absolue de ladite différence, 10113] - le seuil ou la fourchette de valeurs sont déterminés à partir de paramètres estimés sur un simulateur de mouvements de référence, 10114] - le seuil ou la fourchette de valeurs sont déterminés à partir de paramètres estimés sur le même simulateur de mouvements alors que ce dernier est connu pour être dans un état de fonctionnement adéquat, de préférence optimal, 10115] - lorsque la comparaison par rapport au seuil ou à la fourchette de valeurs est signi- ficative d'une anomalie pour un paramètre estimé, alors un signalement est effectué, 10116] - la comparaison est effectuée dans le simulateur de mouvements, 10117] - la comparaison est effectuée en dehors du simulateur de mouvements, dans un équipement informatique distinct ayant récupéré les paramètres estimés mémorisés. [0118] L'invention concerne aussi un simulateur de mouvements comportant un ensemble d'axes, chaque axe étant entraîné en rotation ou translation par un organe effecteur recevant un signal de commande U(t) et comportant au moins un capteur de mesure de la rotation ou de la translation dudit axe permettant de mesurer au moins la position e (t) angulaire ou À ( t) linéaire dudit axe, chaque organe effecteur d'axe étant dans une boucle de correction fermée reliée par une première entrée, au dit au moins un capteur de mesure et produisant, en sortie, le signal de commande U(t), ladite boucle de correction comportant une seconde entrée de signal destinée à recevoir un signal de consigne pour commander ledit axe, ledit simulateur de mouvements comportant un moyen de contrôle configure pour exécuter la méthode de caractérisation et/ou le procédé de de surveillance de l'invention. Le moyen de contrôle peut être un circuit électronique analogique et numérique et, de préférence, à programmation informatique. [0119] L'invention concerne enfin un système de gestion d'un ensemble de simulateurs de mouvements comportant chacun un ensemble d'axes, chaque axe étant entraîné en rotation ou translation par un organe effecteur recevant un signal de commande U(t) et comportant au moins un capteur de mesure de la rotation ou de la translation dudit axe permettant de mesurer au moins la position e (t) angulaire ou À ( t) linéaire dudit axe, chaque organe effecteur d'axe étant dans une boucle de correction fermée gnon g g gs 12 reliée par une première entrée, au dit au moins un capteur de mesure et produisant, en sortie, le signal de commande U(11, ladite boucle de correction comportant une seconde entrée de signal destinée à recevoir un signal de consigne pour commander ledit axe, dans lequel système, les simulateurs de mouvements comportent des moyens de contrôle et sont. connectés à un réseau informatique relié à un serveur distant, les simulateurs de mouvements et le serveur pouvant échanger des données, ledit serveur distant comportant des moyens informatiques, et les moyens de contrôle des simulateurs dc mouvements et les moyens informatiques du serveur sont configurés pour exécuter : [0120] - la méthode de caractérisation de l'invention permettant d'obtenir des estimations des paramètres liés au fonctionnement desdits simulateurs de mouvements, ou [0121] - le procédé de surveillance de [0122] et les paramètres estimés sont stockés dans le serveur. [0123] Dans diverses variantes du système de gestion : [0124] - la méthode est exécutée dans le simulateur de mouvements, [0125] - la méthode est exécutée dans le serveur, [0126] - la méthode est exécutée en partie dans le simulateur de mouvements et en partie dans le serveur, [0127] - le procédé est exécuté dans le simulateur de mouvements, [0128] - le procédé est exécuté dans le serveur, [0129] - la procédé est exécuté en partie dans le simulateur de mouvements et en partie dans le serveur. Brève description du dessin [0130] [fig.1] de l'état de l'art représente le schéma dc principe de l'identification par la méthode IDIM-IV faisant en particulier utilisation d'un modèle dynamique direct du système. Description détaillée d'un exemple de réalisation [0 131] La description qui va suivre en regard des dessins annexés, donnés à titre d'exemples non limitatifs, fera bien comprendre en quoi consiste l'invention et comment elle peut être réalisée. [0132] Afin de pouvoir présenter l'invention dans tous ses détails, la description qui suit va commencer par présenter les principes généraux d'un simulateur de mouvements puis une solution existante permettant de caractériser un simulateur de mouvements et, enfin, la solution proposée permettant une optimisation de la caractérisation d'un simulateur de mouvements. Pour ces explications on considérera des axes entraînés en rotation dans le simulateur de mouvements mais l'invention peut aussi s'appliquer à un axe entraîné en translation rectiligne (ce qui entraîne une simplification du fait de 13 l'absence d'effet du balourd) voire mixte si on peut séparer les effets entre la rotation et la translation rectiligne. [0133] Dans ce document, pour ce qui concerne les axes entraînés en rotation, le symbole 0 et similaires est utilisé. Pour ce qui concerne les axes entraînés en translation, le symbole 7v et similaires est utilisé. Concernant les similaires : les symboles 0 et X, sans point sur le dessus, se rapportent à des positions qu'elles soient calculées ou mesurées. Les symboles é et À (donc avec un seul point sur le dessus) se rapportent à des vitesses, qu'elles soient calculées ou mesurées. Les symboles é et À (donc avec deux points sur le dessus) se rapportent à des accélérations, qu'elles soient calculées ou mesurées. [0134] On peut mettre en oeuvre les principes suivants dans le cadre de la modélisation d'un axe d'un simulateur de mouvements. Tout d'abord, un simulateur de mouvements peut être considéré comme un robot d'une forme particulière, et le formalisme mathématique pour sa modélisation résulte des équations d'Euler-Lagrange. Cependant, il est toujours possible de caractériser un axe séparément des autres, les paramètres résultant des couplages inter-axes n'étant alors pas pris en compte. En outre, les conditions d'essai font qu'un axe entraîné en rotation peut être testé lorsqu'il est vertical ou non vertical ce qui fait que la modélisation peut être plus ou moins simplifiée (on verra que le calcul est simplifié pour un axe vertical du fait de l'absence d'effet du balourd). [0135] Soit 6 ( t ) le signal de commande envoyé sur l'amplificateur de l'actionneur de l'axe, actionneur généralement électrique. La plupart du temps, l'assemblage amplificateur, électronique de puissance et moteur est établi de sorte que le couple moteur ni (0 est pratiquement linéaire par rapport à ladite commande [I(f), les dynamiques électriques propres à cet actionneur étant ici négligées. On peut alors écrire : [0136] Fm = Kamti(t) [0137] Où K (in, est un scalaire, produit de la constante de couple du moteur par le gain de l'amplificateur de puissance. Considérons par ailleurs la position angulaire e (t) de l'axe, ici en radians. La simple application de la seconde loi de Newton (théorème fondamental de la dynamique) amène aux équations différentielles suivantes pour un axe vertical : [0138] Je (t) = - fve(t) - fssign(O) + J ,U(t) [0139] Et pour un axe horizontal : [0140] J (t) = - fve (0 - fssign(0)- f,' sin (9 + ço ) IC '29t) [0141] J : étant le moment d'inertie (k g m 2) [0142] f: étant un coefficient de frottements visqueux ( N. in.$) 14 [0143] f s : étant un coefficient de frottements secs (N. ) (les frottements secs étant ici représentés par le modèle dit de Coulomb). 101441 f : étant un coefficient de balourd (couples de gravité) ( N.m ) résultant de ce g que le centre de gravité de l'axe ne se situe pas sur l'axe de rotation. [0145] (1) : un angle en relation avec la position du centre de gravité par rapport à la verticale (r a d). [0146] Les équations précédentes peuvent être réécrites pour un axe vertical : [0147] v If; ( ) L' K v kd- Sign = n tTI a a m m [0148] Et pour un axe horizontal : [0149]f sin (6 + ri) K e +f1' e + Kfs sign(9 ) g UW UW U rt7 Ci 117 [0150] Ces équations constituent ce que l'on appelle le modèle inverse d'un axe de si- mulateur de mouvements. [0151] Posons : [0152] n - K - a [0153] f ', . m [0154] f 's KR 111 [0155] g K Li )71 [0156] Il est clair que pour unc inertie invariante, la comparaison des paramètres n g obtenus par estimation, dans une situation normale et dans une situation où l'un des composants se serait dégradé, permet d'obtenir des informations afin de localiser l'origine d'un défaut. [0157] Ainsi, une différence significative de la valeur du paramètre ri il obtenue par es- timation entre un cas normal et un cas pathologique à inertie identique, révèle une variation de la constante de couple du moteur, ou du gain de l'amplificateur, ce qui permet d'inférer qu'une anomalie est survenue au niveau du variateur et/ou de l'amplificateur du moteur, ou bien au niveau du moteur lui-même. [0158] De même, une variation significative des coefficients liés aux frottements visqueux ou secs rl s permet de déduire une anomalie au niveau des roulements associés à cet axe. [0159] Afin qu'un diagnostic fiable puisse être établi, il est nécessaire que le test puisse être toujours effectué à inertie invariante, en général sans charge sur le plateau du simulateur. D'autre part il est nécessaire également d'avoir à disposition un algorithme d'estimation des paramètres en question, dénommé algorithme d'identification, qui soit précis. [0160] On connait des algorithmes d'estimation des paramètres d'un simulateur par identi- fication mais l'identification des paramètres T7 a na;,, ns, n a présente une difficulté majeure : Cette opération nécessite qu'un axe du simulateur soit stimulé afin d'enregistrer son comportement. Les axes des simulateurs ayant par construction un caractère intégrateur, la stimulation doit être effectuée en boucle fermée, c'est-à-dire que durant cette phase Faxe est déjà contrôlé par un correcteur. Or il est bien connu que d'une façon générale, l'identification des systèmes en boucle fermée est un domaine bien plus difficile que l'identification des systèmes en boucle ouverte. En particulier, l'application des algorithmes d'identification en boucle ouverte appliqués à des données acquises lors d'essais effectués en boucle fermée produisent des estimés biaisés, dès lors que le système est entaché de perturbations aléatoires, ce qui est systématiquement le cas sur un simulateur de mouvements. 101611 Aussi, la mise au point d'une technique fiable produisant des paramètres estimés non biaisés et précis est un enjeu crucial pour la pertinence d'un diagnostic d'un simulateur de mouvements. [0162] Le développement d'algorithmes spécifiques à l'identification des systèmes électro- mécaniques a fait l'objet de nombreuses recherches. En particulier, on peut mentionner la méthode d'identification dynamique du modèle inverse en temps continu, connue dans la littérature sous le sigle IDIM (Inverse Dynamic Identification Model), qui peut être utilisée, associée à une technique d'estimation par la méthode des moindres carrés explicites (IDIM-LS). Cette méthode existe depuis les années 1980, mais, comme toutes les méthodes incluant une minimisation basée sur l'algorithme des moindres carrés explicites aboutissent généralement à des paramètres estimés biaisés en cas de perturbations aléatoires, elle n'est pas satisfaisante. [0163] Aussi plusieurs techniques d'identification des systèmes électromécaniques ont 6t6 développées beaucoup plus récemment afin de pallier aux inconvénients de la méthode IDIM-LS. Parmi celles-ci figure une méthode publiée en 2014 et basée sur le principe de la variable instrumentale (IDIM-IV), et décrite en particulier dans la publication « A generic Instrumental variable approach for industrial robot identification » dont les auteurs sont A. Janet et ai, parue dans IEEE Transactions on control systems technology (vol. 22, pp.132-145). [0164] Cette méthode va être décrite dans la suite de ce document. Néanmoins, l'un de ses inconvénients est qu'elle nécessite la mise en oeuvre non seulement d'un modèle 16 inverse, en l'occurrence du simulateur de mouvements, niais également d'un modèle direct du système muni de son correcteur, c'est-à-dire du modèle complet de la boucle fermée. [0165] On peut déjà annoncer à ce stade de la description, compte tenu des spécificités d'un simulateur de mouvements et en particulier de la possibilité de caractériser chaque axe séparément, qu'on propose pour l'invention une version dérivée de la méthode IDIMIV et qui est adaptée aux simulateurs de mouvements. Cette version dérivée et adaptée apporte des simplifications drastiques dans l'algorithme, en faisant l'économie d'un modèle de la boucle fermée, et ne nécessite pas de connaissance à priori du correcteur utilisé. Ces deux éléments de simplifications favorisent grandement la mise en oeuvre industrielle de l'algorithme, par la réduction des coûts de calcul qu'il engendre. [0166] Concernant la méthode IDIM-LS en elle-même, reprenons le modèle inverse de l'axe horizontal, celui de l'axe vertical pouvant être considéré comme un cas dégénéré avec 0 : [0167] n 8 + n ,,0 + ri sign(0)+ ti sin (e + =U [0168] Qui, du fait d'une propriété trigonométrique élémentaire, peut encore s'écrire : [0169] ri ('0 (1+ ,6 (t)+ r f g n(B( t )) + c sin ( 6 ( t ) ) + scos (G ( t )) U(t) [0170] Avec [0171] =,7 2 g gs [0172] 1.7 yo = a tarz2( L' [0173] La méthode IDIM-LS, comme toutes les méthodes d'identification, consiste d'abord à exciter l'axe, ici en boucle fermée, et d'acquérir les entrées et les sorties, en l'occurrence les signaux U(t) et e (t). [0174] On procède ensuite à des calculs par application de fonctions aux signaux acquis U(t) et e (t) et de filtrages des résultats de ces calculs, de façon à obtenir les données filtrées suivantes dans le cas où les calculs et filtrages sont effectués en commun [0175] [0176] [0177] [0178] U (rs+1)2U(t) e ro-) (Ts1±',e(t) érw_ (Ty-±n, s e(t) erm s2 ( (rs+1)- 17 [0179] S ig t )s 1+ )2s ig (à) (t)) [0180]1 se /- (t) (rs +1) 2 sin (e (t) [0181] cOf(t) - r s 1 ) 2cos to(t)) [0182] On note que ces fonctions s'appliquent à cc cas où un capteur de mesure de la rotation permet seulement de mesurer la position angulaire e t ) de l'axe et il est nécessaire de calculer la vitesse et l'accélération angulaires à partir de la position angulaire mesurée. Si le capteur de mesure de la rotation permet de mesurer en plus la vitesse angulaire e t ) de l'axe alors l'accélération angulaire pourra être calculée par une fonction adéquate à partir de la vitesse angulaire. On peut prévoir le cas de capteurs de mesure de la rotation permettant de mesurer la position angulaire O t ) et la vitesse angulaire Û ( t ) et, alors, on ne calculera que l'accélération angulaire. Dans ce dernier cas, on peut envisager de calculer la vitesse à partir de la position pour faire des comparaisons et éventuellement des corrections entre les calculs et mesures. 101831 Le filtrage est naturellement réalisé en temps discret par discrétisation des filtres continus. Ici est une constante de temps qui constitue un paramètre de réglage. [0184] Il est important de noter que la fonction de transfert de filtrage présentée n'est donnée qu'à titre d'exemple, l'emploi d'autres fonctions de transfert est possible, en fait il peut s'agir d'un filtre linéaire stable, passe bas, d'ordre éventuellement plus élevé et de gain statique éventuellement différent de 1. En outre les fonctions appliquées sont données en référence à l'exemple présenté où le capteur de mesure de la rotation de l'axe permet de mesurer la position angulaire O t ) dudit axe. Il en résulte que la vitesse et l'accélération de l'axe sont calculées par application d'une fonction dérivée première et seconde respectivement. On comprend que si le capteur de mesure de la rotation de l'axe permet de mesurer la vitesse angulaire e t ) . voire l'accélération angulaire 9 (t ), en complément de la position et/ou l'accélération, les fonctions (dont une absence de fonction complète du fait qu'on mesure directement la variable concernée et qu'on a pas besoin de la calculer, il suffit seulement de la filtrer) appliquées aux mesures pourront être différentes pour les calculs de la position et/ou la vitesse et/ou l'accélération de Faxe. Ainsi, si le capteur de mesure de la rotation de l'axe permet de mesurer la position et la vitesse, seule l'accélération devra être calculée par application d'une fonction sur une des deux mesures de position ou de vitesse. [0185] On doit noter que si on peut appliquer des fonctions donnant la dérivée d'une mesure pour obtenir la variable qui sera filtrée, en revanche, l'intégration d'une mesure n'est 18 pas envisagée dans les fonctions applicables. [0186] On peut écrire : [0187] f(t)+ (t)± r ssign f(6) (t)) + acs0 f(t)+ tLisct9 f(t) - [0188] Définissons le vecteur ligne : [0189] / y ( ) = [0 f(t) e f(t) sign/(f) (0) se f(t) Fe f(t)] [0190] dans le cas d un axe horizontal et 10191] 1 ( = [a r(t) e r(t) sign rte (t))1 [0192] dans le, cas d'un axe vertical ou en translation. [0193] On peut alors constituer une matrice X par concaténation des vecteurs-lignes (t ) pour t = 1 N, où N le nombre d'échantillons, telle que : [0194] Ix(l) Ix(2) 1\.(N) [0195] De même, on peut constituer un vecteur Y tel que : [0196] Ut(1) Ur(2) Uf(N) [0197] Poson ri le vecteur des paramètres : [0198] n n n , ngc ri iris [0199] pour ut axe horizontal, et [0200] q n 10201] pour un axe ertical ou en translation. 10202] On alors : 10203] Y = Xri P [0204] 19 p (1) est un vecteur de perturbations aléatoires (p (0). (2) N 102051 L'estimation de n par la méthode des moindres carrés explicites, ce qui correspond à l'algorithme IDIM-LS, fournit le vecteur estimé iï de '7 par la relation classique : [0206] 15. = (xTx)-Ix.TY [0207] Comme expliqué plus haut, il est bien connu que la méthode des moindres carrés aboutit à un vecteur estimé biaisé à partir du moment où les composantes de la matrice X sont corrélées avec le vecteur de perturbations aléatoires, cc qui est le cas dans un simulateur de mouvements. [0208] Maintenant, si on s'intéresse à la méthode de la variable instrumentale (algorithme IDIM-IV), on utilise dans le calcul de i' (estimation de (7 par une méthode de calcul utilisant une matrice Z dont les composantes seraient fortement corrélées avec les composantes de X si ces composantes n'étaient pas perturbées et dont, en revanche, les composantes de Z sont fortement décorrélées des composantes du vecteur P. [0209] L'estimée par la méthode de la variable instrumentale est donnée par la relation : [0210] n =. (ZTX)-1-ZTY 1021 1 Dans le document « A generic Instrumental variable approach for industrial robot identification » déjà cité, la matrice instrumentale Z est constituée à partir des signaux internes à un modèle de la boucle fermée du système (incluant un modèle direct du système du correcteur et générateur de trajectoire associé). La mise en oeuvre de ce modèle de la boucle fermée tel que l'on peut le voir dans la figure 4 de ce document et qui correspond à la figure 1 de la présente demande, alourdit la charge de calcul de l'algorithme d'estimation, ce qui peut constituer un frein à son implantation industrielle. 102121 A noter que dans ces explications, on a considéré tous les paramètres, c'est-à-dire ceux en relation avec le moment d'inertie et les frottements visqueux et les frottements secs et le balourd. Or, dans certaines configurations de simulateurs de mouvements, certains de ces paramètres peuvent être considérés d'effets négligeables et/ou comme non évolutifs. On a déjà vu que dans le cas d'un axe vertical en rotation, le balourd n'intervenait pas. Il peut en être de même pour ce qui concerne les frottements en cas d'utilisation de paliers sans contact et dont on considère qu'ils resteront sans contact, non évolutifs. Dans de tels cas, les éléments correspondants peuvent être éliminés ou être assignés à une valeur par défaut constante et ne pas nécessiter des calculs les concernant par exemple on peut mettre en oeuvre un palier magnétique sans contact. [0213] Attendu que l'identification d'un simulateur de mouvements se fait de façon séparée pour chaque axe, et que l'on manipule de ce fait des modèles découplés, on se retrouve dans un cas favorable, permettant des simplifications de l'algorithme IDIM-IV. Cc sont ces simplifications qui font l'objet de la présente invention. [0214] Ces simplifications portent sur la composition du vecteur instrumental et de la matrice Z associée. On a vu que dans la matrice X, chaque ligne s'écrit : [0215] ity( t)= [gf(t) 6(t) Sign (É) ('t)) Se f(t) Ce f(t)] [0216] Il apparaît tout (l'abord que les mouvements destinés à la caractérisation d'un axe de simulateur de mouvements peuvent être choisis de telle sorte qu'ils ne présentent pas de changement de sens. Cela implique alors que s ig n (() (1-.)) peut être considéré comme indépendant de la perturbation aléatoire p(0. [0217] Par ailleurs, le niveau de perturbation sur e (t) est généralement extrêmement faible au regard de l'amplitude de e (t) avec un rapport signal sur bruit pouvant être supérieur à 60 dB (écart type). En conséquence, les termes cet( ), se 1.( t) sont extrêmement peu affectés par les perturbations aléatoires. [0218] Dans le vecteur /A- ( t ) restent alors les ternies 8 R f / et O. r(r) qui eux sont mani- festement corrélés avec les perturbations. Cependant, il est à noter que les systèmes de commande des simulateurs de mouvements sont pourvus de générateurs de trajectoire calculant en temps récl non seulement la consigne de position a c(t ), mais également une consigne de vitesse 0 JO, et également une consigne d'accélération 0 e(t ). Il est évident que ces consignes ne sont pas du tout corrélées avec les perturbations aléatoires affectant le système. En revanche, les signaux 9 c(t) et a coi sont fortement corrélés avec 0 ( t ) et 0 (t) respectivement et ils peuvent donc être considérés pour établir la matrice Z. [0219] A partir de e '.(t) et Û ,,(t) il est aisé de calculer les signaux filtrésf(t)et e cr(t) tels que : 102201 e( t ) - 1 e .(t) c f (rs +1) [0221] e( t ) _ 1 (rs +1)20 c(t) c f 102221 Les signaux Ût,(t) et 0 f(i) .ont totalement décorrélés d'avec les perturbations. [0223] On peut alors constituer une matrice instrumentale Z simplifiée par concaténation de vecteurs-lignes 1z ( t ) tels que : [0224] /z(t) = [O (t)ef(t) sign .(t9 (t)) se f(t) ce f(t)] 21 [0225] où les composantes sont fortement corrélées à X (s'il n'était pas contaminé par les perturbations aléatoires) et sont fortement décorrélées des perturbations aléatoires et où les termes ce (t) , se t( t) sont fortement décorrélés des perturbations aléatoires. 102261 Et la matrice Z simplifiée est alors constituée suivant : 102271 Z lz(N) [0228] On peut alors estimer les paramètres par la méthode de la variable instrumentale pré- cédemment vue et en utilisant la matrice Z simplifiée, avec : [0229] = (Z1)(YEZ1Y [0230] On voit que l'avantage de la méthode proposée ici est de permettre d'obtenir à bon compte le vecteur instrumental jz( t ) en faisant l'économie du modèle de la boucle fermée de la méthode IDIM-IV ce qui, encore une fois, constitue une réduction de coût significative en termes de charge de calcul. [0231] Avec la méthode proposée, il est possible d'obtenir un vecteur nonn biaisé avec une précision d'estimation inférieure au %, voire inférieure au dixième de %. La grande précision sur les paramètres estimés et le caractère non biaisé de ces estimées permet l'implémentation dc fonctionnalités de prévention et/ou dc diagnostic utilisant ces résultats d'identification par un suivi de l'évolution des paramètres estimés au cours du temps. [0232] En outre, avec la méthode proposée, on traite chacun des axes indépendamment des autres, ce qui apporte aussi une simplification : le signal d'excitation est envoyé dans la boucle du correcteur de l'axe étudié et les mesures et calculs ne concernent que cet axe pour obtenir les paramètres estimés de cet axe. [0233] Dans un mode de réalisation particulier, le simulateur de mouvements comporte un moyen de contrôle de son fonctionnement demandant à l'utilisateur ou même contraignant l'utilisateur à lancer périodiquement, par exemple tous les mois, ou toutes les semaines, une caractérisation du simulateur de mouvements dans le but d'obtenir le vecteur estimé i3 des paramètres. Les paramètres estimés sont ensuite stockés dans une base de données. Cette acquisition périodique permet alors de se livrer à des statistiques afin d'observer des tendances concernant l'évolution desdits paramètres estimés. [0234] Le moyen de contrôle du simulateur de mouvements, qui peut être du type microor- dinateur, peut comporter le générateur de consigne et être configure pour calculer les paramètres estimés. Dans ce cas, le moyen de contrôle a directement accès à la consigne de vitesse ( ,(t) ou À «(t) et à la consigne d'accélération Û ou t) et peut utiliser directement ces données pour les calculs. Par contre, si le générateur de consigne est indépendant du moyen de contrôle du simulateur alors cc dernier peut ne recevoir que la consigne correspondant au signal d'excitation et, dans cc cas, le moyen de contrôle devra calculer à partir de la consigne reçue par le simulateur de mouvements la consigne de vitesse et la consigne d'accélération. [0235] Lc simulateur de mouvements peut comporter un mode maintenance dans lequel les paramètres estimés du simulateur de mouvements et l'évolution de ces paramètres peut être affichés sur un dispositif d'affichage dudit simulateur de mouvements et/ou, dans une version évoluée, une alerte et/ou un diagnostic automatique produit en cas d'anomalie des paramètres estimés. [0236] De préférence, le simulateur de mouvements qui est sous le contrôle d'un moyen de contrôle à programmation informatique/logiciel pour son fonctionnement, comporte des moyens de comptage temporel et/ou de nombre d'utilisations qui permettent de déclencher périodiquement (en fonction du temps et/ou du nombre d'utilisations depuis la précédente caractérisation) ou de demander le déclenchement par l'utilisateur de la méthode dc caractérisation du simulateur de mouvements. Dans une variante, le simulateur de mouvements est connecté sur un réseau informatique, par exemple ETHERNET et/ou INTERNET, et le déclenchement est commandé à distance depuis un serveur distant ayant un moyen de suivi du fonctionnement du simulateur dc mouvements. [0237] On peut prévoir divers modes de fonctionnement en fonction de la présence ou non d'une connexion du simulateur de mouvements sur réseau informatique. Par exemple en mode local, le moyen de contrôle du simulateur de mouvements permet d'exécuter localement la méthode et de calculer en local les paramètres estimés et de faire des comparaisons pour alerte et/ou diagnostic, les résultats, en particulier les paramètres estimés, pouvant éventuellement être transmis à un serveur distant directement si une connexion est disponible ou sinon indirectement, par exemple par déchargement des paramètres estimés mémorisés localement dans une clef USB ou un ordinateur portable d'un technicien d'entretien. En mode distant, un serveur distant exécute la méthode et pour cela reçoit les mesures et reçoit (si le générateur de consigne du simulateur de mouvements est utilisé) ou transmet (si le serveur peut comporter ou simuler un générateur de consigne) le signal d'excitation et, avec les résultats des paramètres estimés calculés par le serveur, le serveur peut faire des comparaisons pour alerte et/ou diagnostic. Dans d'autres variantes de modes mixtes, ces tâches sont partagées entre le simulateur de mouvements et le serveur distant. [0238] Ainsi, grâce à une telle connexion sur réseau informatique, une exploitation sta- 23 tistique des paramètres estimés recueillis et mémorisés peut être faite à distance par des équipes de maintenance ou par le constructeur, en vue de détecter le plus tôt possible une anomalie de fonctionnement susceptible d'amener une interruption du fonctionnement du simulateur. Dans des variantes, l'exploitation statistique peut être faite sur site, dans le simulateur de mouvements lui-même ou dans un ordinateur, possiblement portable, relié au simulateur de mouvements. [0239] Un dispositif d'alerte de l'utilisateur et du constructeur peut donc être prévu en cas d'anomalie sur les paramètres estimés. [0240] Dans un mode de mise en oeuvre avantageux, plusieurs simulateurs de mouvements sont sur des réseaux informatiques et envoient lors des caractérisations les mesures et signaux d'excitation et/ou directement les résultats de calcul des paramètres estimés dans une hase de donnée commune sur un serveur distant où la méthode peut être (re)exécutée, possiblement améliorée/affinée pour éventuellement, en retour, mise à jour des programmes de la méthode de caractérisation dans les simulateurs de mouvements, des comparaisons et calculs statistiques peuvent être effectués dans le serveur pour suivi d'une flotte de simulateur de mouvements. Un fabriquant de simulateurs de mouvements pourra donc détecter les problèmes potentiels et les résoudre en amont et sur les nouveaux produits qu'il fabriquera. En outre, la gestion du service après-vente et des contrats d'entretien en sera améliorée. [0241] Parmi les anomalies détectables et alertes, et à titre d'exemple, une tendance haussière du paramètre A peut laisser supposer qu'un défaut progressif au niveau de " l'actionneur ou de son amplificateur est en train de survenir. Dc même, une tendance haussière des paramètres A A peut laisser présager de l'installation d'un défaut y' '" s au niveau (les roulements ou des joints tournants (lu simulateur de mouvements. [0242] Lorsqu'une tendance haussière d'un des paramètres est détectée, un message (l'alerte est envoyé à l'utilisateur du simulateur (le mouvements. [0243] De même une valeur de A ou n significativement différente de la valeur I s nominale du paramètre en question, permet de diagnostiquer une anomalie immédiate sur l'un des organes du simulateur de mouvements et (l'envoyer à l'utilisateur un message d'alerte. [0244] Ces diverses alertes permettent alors à l'utilisateur d'engager une procédure de maintenance préventive et/ou corrective de l'équipement, avant que son fonctionnement ne soit totalement interrompu, ce qui d'un point de vue économique est très appréciable. [0245] Dans un mode de réalisation avancé, les valeurs estimés des paramètres, , obtenues périodiquement, sont envoyées via un réseau informatique vers un serveur d'équipes 24 de maintenance ou du constructeur du simulateur de mouvements afin de leur permettre une exploitation statistique des données récoltées, et de détecter une anomalie susceptible de mener à une interruption du simulateur, de préférence sans recourir à l'utilisateur de l'équipement. [Revendication I]