Stability analysis of admittance control using asymmetric stiffness matrix

1. Problématique

Dans les tâches d'automatisation impliquant de nombreux contacts (comme l'assemblage de précision type « pêne dans trou »), le contrôle de la compliance est crucial. Les méthodes de contrôle d'impédance et d'admittance permettent d'ajuster la rigidité du robot pour gérer les erreurs de positionnement.

• Limitation actuelle : Les approches basées sur la passivité sont couramment utilisées pour garantir la stabilité, mais elles imposent que la matrice de rigidité soit symétrique.
• Opportunité manquée : L'utilisation d'une matrice de rigidité asymétrique peut étendre la plage de mouvement induite et améliorer les performances de l'assemblage. Cependant, l'asymétrie introduit un champ de force non conservatif (champ de rotation ou « curl ») qui peut exciter des oscillations spirales et violer la condition de passivité, rendant la stabilité difficile à prouver.
• Objectif : Déterminer les conditions de stabilité pour un contrôle d'admittance utilisant une matrice de rigidité asymétrique, sans se limiter à l'analyse de passivité.

2. Méthodologie

Les auteurs se concentrent sur le contrôle d'admittance d'un bras robotique en 2D, en supposant une boucle de contrôle de position interne stable.

• Modélisation :
  • Le modèle d'admittance est défini par $F_{ext} = M\ddot{x} + D\dot{x} + Kx$.
  • La matrice de rigidité $K$ est décomposée en une partie symétrique $K_s$ (champ de force convergent) et une partie asymétrique $K_a$ (champ de force rotatif).
  • La matrice d'amortissement $D$ est supposée diagonale.
• Analyse de stabilité par Lieu des Racines (Root Locus) :
  • Au lieu de s'appuyer uniquement sur la passivité, les auteurs analysent les valeurs propres (eigenvalues) de la matrice du système d'état.
  • Ils dérivent les conditions mathématiques pour lesquelles les valeurs propres du système restent réelles et négatives (assurant la convergence sans oscillation).
• Déduction des conditions critiques :
  • Ils identifient que l'apparition de composantes imaginaires dans les valeurs propres de la matrice de rigidité est la cause des oscillations spirales.
  • Ils établissent une relation entre le coefficient d'amortissement visqueux ($d$) et les paramètres de rigidité ($k_x, k_y, k_s, k_a$) pour garantir que toutes les valeurs propres du système soient réelles et négatives.

3. Contributions Clés

1. Identification de la cause des oscillations spirales : Les auteurs démontrent que les oscillations spirales sont excitées par la partie asymétrique de la matrice de rigidité lorsque les valeurs propres de cette matrice ne sont pas réelles.
2. Condition d'absence d'oscillation : Ils prouvent que la condition pour éviter les oscillations spirales est équivalente à l'exigence que les valeurs propres de la matrice de rigidité soient des nombres réels.
3. Dérivation de la limite de stabilité : En utilisant la méthode du lieu des racines, ils dérivent une inégalité précise pour le coefficient d'amortissement minimal nécessaire ($d$) en fonction de la rigidité asymétrique ($k_a$) et des rigidités diagonales ($k_x, k_y$).
   • La condition clé est : $d > \sqrt{\frac{k_a^2}{k_x}}$ (dans le cas simplifié $k_x = k_y$).
   • Ils montrent également que si la condition sur les valeurs propres réelles est satisfaite, un amortissement suffisant peut garantir la convergence même si la passivité locale n'est pas respectée à chaque instant (notamment lors des phases de faible vitesse).
4. Validation par simulation et expérience : La théorie est validée par des simulations numériques et des expériences réelles sur un bras robotique.

4. Résultats

• Simulations :
  • Lorsque la condition de valeurs propres réelles est violée (ex: $k_a$ trop élevé par rapport à $k_x$), le système oscille même avec un amortissement élevé, et le point critique d'amortissement n'est jamais atteint.
  • Lorsque la condition est respectée, un amortissement suffisant permet d'éliminer les oscillations et d'atteindre un amortissement critique ($\zeta = 1$).
  • Les simulations montrent que l'énergie totale peut augmenter localement lors des phases de faible vitesse (violation de la passivité locale), mais que la convergence globale est garantie si les paramètres respectent les limites dérivées.
• Expériences :
  • Un robot a été programmé pour suivre une surface métallique.
  • Comparaison entre une matrice symétrique et une matrice asymétrique (triangulaire).
  • Avec la matrice asymétrique (paramétrée selon les limites de stabilité), l'erreur de positionnement finale était réduite (6 mm contre 8 mm pour la matrice symétrique) et la force de contact était mieux maintenue, prouvant l'efficacité de l'approche pour des tâches de contact précis.

5. Signification et Perspectives

• Impact : Ce travail élargit considérablement la plage de conception des contrôleurs d'admittance. Il permet d'utiliser des matrices de rigidité asymétriques pour améliorer les performances des tâches d'assemblage, là où les méthodes traditionnelles basées sur la passivité imposaient des contraintes trop restrictives.
• Innovation théorique : Il déplace le paradigme de la garantie de stabilité de la « passivité stricte » vers une analyse basée sur les valeurs propres et le lieu des racines, offrant une solution pour les systèmes non passifs mais stables.
• Travaux futurs : Les auteurs prévoient d'étendre cette théorie aux systèmes à 3 degrés de liberté (ou plus) et d'intégrer le concept de « réservoirs d'énergie » (energy tanks) pour généraliser l'analyse de stabilité aux systèmes à rigidité variable.

En résumé, l'article fournit un cadre mathématique rigoureux pour exploiter l'asymétrie de la rigidité en robotique, en définissant précisément les limites de stabilité qui permettent d'éviter les oscillations indésirables tout en améliorant la précision des tâches de contact.