Stability analysis of admittance control using asymmetric stiffness matrix

Análisis de Estabilidad del Control de Admitancia utilizando una Matriz de Rigidez Asimétrica

1. Planteamiento del Problema

En tareas de automatización que implican un contacto rico con el entorno (como el ensamblaje de precisión), el ajuste de la rigidez del sistema de control es un factor crítico. Aunque el uso de una matriz de rigidez asimétrica en el control de admitancia puede extender el rango de inducción de movimiento y mejorar la capacidad del robot para manejar restricciones mecánicas variables, esto presenta un desafío fundamental: la inestabilidad.

Los métodos tradicionales de análisis de estabilidad, basados en la pasividad, fallan cuando la matriz de rigidez es asimétrica. La parte asimétrica genera un campo de fuerza no conservativo (un campo de "curl") que puede excitar oscilaciones espirales, impidiendo que se garantice la estabilidad mediante los enfoques estándar. Hasta este estudio, no existía una demostración teórica rigurosa sobre cómo asegurar la estabilidad en estos sistemas asimétricos.

2. Metodología

Los autores proponen un enfoque basado en el análisis de la estabilidad de convergencia utilizando el método del lugar geométrico de las raíces (root locus) y el análisis de autovalores, en lugar de depender exclusivamente de la pasividad.

• Modelo: Se utiliza un modelo de admitancia planar (2D) con un bucle de control de posición interno estable. La ecuación de admitancia es $F_{ext} = M\ddot{x} + D\dot{x} + Kx$, donde $K$ es una matriz de rigidez asimétrica ($K = K_s + K_a$) y $D$ es una matriz de amortiguamiento diagonal.
• Descomposición: La matriz de rigidez asimétrica se descompone en una parte simétrica ($K_s$), que crea un campo de fuerza convergente, y una parte asimétrica ($K_a$), que genera un campo de fuerza rotacional.
• Análisis de Autovalores: Se deriva la matriz de estado del sistema y se analizan sus autovalores ($\lambda$). Se demuestra que la presencia de componentes imaginarios en los autovalores de la matriz de rigidez es la causa de las oscilaciones espirales.
• Derivación de Condiciones:
  • Se establece una condición necesaria para que los autovalores sean reales (evitando oscilaciones espirales): $|k_{am}| < \sqrt{k_{sm}^2 + \frac{1}{4}(k_{xm} - k_{ym})^2}$.
  • Se deriva una condición para el coeficiente de amortiguamiento ($d$) que garantiza que todos los autovalores del sistema tengan parte real negativa, asegurando la estabilidad asintótica. La fórmula clave derivada es $d > \sqrt{\frac{k_{am}^2}{k_x}}$ (para el caso simplificado $k_x = k_y$).

3. Contribuciones Clave

1. Identificación de la causa de oscilación: Se clarifica que las oscilaciones espirales en el control de admitancia con matrices asimétricas son excitadas por la parte imaginaria de los autovalores de la matriz de rigidez.
2. Condiciones de Estabilidad: Se demuestra que, aunque la pasividad no se cumple globalmente, la estabilidad de convergencia puede garantizarse si se cumplen ciertas condiciones en los parámetros de la matriz de rigidez y el amortiguamiento.
3. Criterio de Amortiguamiento Crítico: Se establece que, bajo ciertas condiciones de la matriz de rigidez, es posible alcanzar un amortiguamiento crítico (sin oscilación), algo que antes se consideraba imposible con matrices asimétricas.
4. Validación Experimental: Se presenta un método para diseñar parámetros que evitan la inestabilidad, extendiendo el rango de operación segura más allá de lo permitido por las matrices simétricas.

4. Resultados

• Simulaciones:
  • Se compararon respuestas temporales con diferentes relaciones de amortiguamiento ($\zeta$). Cuando la condición de autovalores reales no se cumplía, el sistema mostraba oscilaciones persistentes incluso con alto amortiguamiento.
  • Cuando se cumplían las condiciones derivadas (autovalores reales y amortiguamiento suficiente), el sistema convergía sin oscilaciones.
  • Se validó que si el amortiguamiento $d$ es menor que el umbral calculado, la energía del sistema aumenta localmente (violando la pasividad local), pero si $d$ supera el umbral, la energía total converge a cero.
  • Los mapas de estabilidad generados mostraron una coincidencia casi perfecta entre la región teórica derivada y la región estable observada en la simulación.
• Experimentos:
  • Se realizó una prueba con un brazo robótico que interactuaba con una placa metálica.
  • Se comparó una matriz de rigidez simétrica frente a una asimétrica (triangular).
  • Resultado: Con la matriz asimétrica, el robot logró un seguimiento de trayectoria más preciso (desviación de 6 mm vs. 8 mm en la simétrica) y mantuvo una fuerza de contacto más alta y estable, demostrando la efectividad del método propuesto en un escenario real.

5. Significado e Impacto

Este trabajo es significativo porque rompe la barrera teórica que limitaba el uso de matrices de rigidez asimétricas en el control de robots debido a preocupaciones de estabilidad.

• Ampliación de Capacidades: Permite a los ingenieros diseñar controladores de admitancia con rangos de movimiento inducido más amplios y flexibles, aprovechando la asimetría para tareas complejas de ensamblaje.
• Seguridad Garantizada: Proporciona un marco matemático riguroso para configurar los parámetros de control, asegurando que el sistema no se vuelva inestable a pesar de la falta de pasividad estricta.
• Futuro: El estudio sienta las bases para extender estas teorías a sistemas de 3 o más grados de libertad y a estrategias de control de rigidez variable, utilizando conceptos como "tanques de energía" para generalizar el análisis de estabilidad.

En resumen, el artículo demuestra que, mediante un diseño cuidadoso de los parámetros basado en el análisis de autovalores y el lugar de las raíces, es posible utilizar matrices de rigidez asimétricas de manera segura y efectiva, mejorando el rendimiento en tareas de contacto rico.