Decoupling Torque and Stiffness: A Unified Modeling and Control Framework for Antagonistic Artificial Muscles

Resumen Técnico

1. Planteamiento del Problema

Los músculos artificiales antagonistas (como músculos neumáticos PAM, actuadores electrohidráulicos HASEL y actuadores de elastómero dieléctrico DEA) tienen el potencial de desacoplar el torque y la rigidez de una articulación, imitando el comportamiento biológico donde el "sesgo" (bias) genera torque neto y la "co-contracción" ajusta la rigidez mecánica.

Sin embargo, en la práctica, mantener este control independiente durante transitorios de contacto (impactos repentinos con el entorno) es extremadamente difícil debido a:

• La dependencia no lineal de la fuerza del actuador respecto a la deformación y la entrada.
• Límites de velocidad y efectos viscoelásticos que son más pronunciados en los primeros milisegundos tras el contacto.
• La dificultad de distinguir si el comportamiento observado se debe a la programación de la impedancia en el bucle externo o a la modulación real de la rigidez a nivel del actuador (planta).

El objetivo es lograr un control robusto que permita rastrear comandos de torque y rigidez de forma independiente, incluso durante impactos rápidos y cambios abruptos en la dinámica de contacto.

2. Metodología y Marco Unificado

Los autores proponen un marco unificado en tiempo real aplicable a múltiples familias de músculos artificiales. Este marco consta de tres componentes principales:

A. Modelado Unificado de Baja Orden

• Ley de Fuerza Base Separable (Padé [2/1]): Se utiliza un modelo de fuerza activa $F(\alpha, x) = \alpha F_{base}(x)$, donde $\alpha$ es la activación normalizada y $x$ es la deformación interna. La función $F_{base}$ se modela mediante una fracción racional de Padé [2/1] con cuatro coeficientes.
  • Este modelo captura el comportamiento no lineal (incluyendo el "rollover" de fuerza en HASELs) con baja complejidad computacional.
  • Se adapta a diferentes actuadores mediante mapeos de activación específicos (ej. presión para PAMs, voltaje al cuadrado para HASELs/DEAs).
• Envoltura Dinámica de Dos Estados: Para capturar la banda de frecuencia limitada y la viscoelasticidad interna, se añade un modelo dinámico mínimo:
  1. Un estado de "impulso" ($a$) que filtra la entrada física (voltaje/presión).
  2. Una rama de Kelvin-Voigt en serie que introduce una deformación interna ($x$) y amortiguamiento, separando la deformación externa ($\epsilon$) de la interna.

B. Control en Cascada con Coordenadas Desacopladas
Se implementa un controlador en cascada que opera en coordenadas de co-contracción ($c_\alpha$) y sesgo ($b_\alpha$):

• Bucle Externo (Opcional): Ley de impedancia que genera un torque deseado ($T_{des}$) basado en errores de posición.
• Bucle Interno (Solver de Comando): Resuelve inversamente para encontrar las entradas de los actuadores que logren simultáneamente el torque deseado ($T_{des}$) y la rigidez mecánica deseada ($K_{mech,des}$).
  • Utiliza un solucionador de Gauss-Newton amortiguado (Levenberg-Marquardt) con un paso de predicción de un solo paso (1 ms).
  • Incorpora retroalimentación PI en las coordenadas desacopladas para corregir errores de torque y rigidez.

C. Política de Interacción Adaptativa a la Profundidad
Se propone una política bio-inspirada que ajusta la rigidez en función de la profundidad de penetración ($\delta$) en el entorno, sin necesidad de etiquetas explícitas de la superficie:

• La rigidez aumenta linealmente desde un valor bajo ($K_{low}$) hasta uno alto ($K_{high}$) a medida que aumenta la penetración.
• Esto permite un compromiso automático: baja rigidez para absorber impactos (reducir choque) en contactos duros, y alta rigidez para mantener la estabilidad en superficies blandas.

3. Contribuciones Clave

1. Modelo Unificado: Un modelo de fuerza Padé [2/1] con una envoltura dinámica de dos estados que es válido para PAMs, HASELs y DEAs, validado experimentalmente en hardware HASEL.
2. Control Desacoplado en Tiempo Real: Un controlador en cascada que logra un rastreo independiente de torque y rigidez en menos de 1 ms por ciclo de control, demostrando robustez ante transitorios de contacto.
3. Evaluación de Interacción: Un protocolo de contacto simulado que incluye:
   • Una prueba de "escalón de rigidez con torque fijo" para validar la capacidad de cambiar la rigidez sin perturbar el torque.
   • Un estudio de auditoría de acoplamiento que descompone los efectos de la programación de impedancia virtual frente a la modulación de rigidez mecánica real.

4. Resultados

• Precisión del Modelo: El modelo identificado en HASELs alcanzó un error cuadrático medio (RMSE) de 0.78 N (aprox. 10% en el rango de operación) y un $R^2$ de 0.83, capturando tanto tendencias cuasi-estáticas como respuestas transitorias.
• Desacoplamiento: El controlador de lazo cerrado mantuvo el desacoplamiento de torque y rigidez incluso bajo perturbaciones externas. En comparación con un control de lazo abierto, el error de torque se redujo en un factor de 29x ante perturbaciones de 1.0 N·m.
• Prueba de Escalón de Rigidez: El sistema pudo realizar transiciones rápidas de rigidez manteniendo el torque constante, reduciendo significativamente los "golpes" de torque y la carga acumulada durante el cambio.
• Política Adaptativa vs. Fija:
  • En entornos rígidos, la política de rigidez fija baja (Fixed-Low) y la adaptativa redujeron el pico de impacto (shock) y la carga acumulada en comparación con la rigidez fija alta.
  • En entornos blandos, la rigidez fija alta y la adaptativa mostraron mayor estabilidad (menor oscilación) que la rigidez fija baja.
  • La política adaptativa logró un rendimiento intermedio óptimo, ofreciendo baja carga en impactos duros y alta estabilidad en superficies blandas, sin necesidad de identificar la rigidez del entorno explícitamente.
• Análisis de Atribución: Los experimentos de ablación (2x2) demostraron que tanto la programación de la ganancia virtual como la modulación mecánica contribuyen al rendimiento, siendo la combinación de ambas la más efectiva en condiciones nominales.

5. Significado e Impacto

Este trabajo proporciona una base fundamental para el control de robots musculoesqueléticos que utilizan músculos artificiales antagonistas. Al lograr un desacoplamiento robusto de torque y rigidez en tiempo real, incluso durante contactos dinámicos, se habilita el desarrollo de robots más seguros y adaptables para entornos no estructurados.

La propuesta demuestra que es posible implementar comportamientos de impedancia adaptativa complejos (como la mitigación de impactos y la estabilidad en superficies variables) utilizando un modelo unificado de baja complejidad y un controlador eficiente, sentando las bases para la próxima generación de robots blandos interactivos.

Limitaciones y Trabajo Futuro:
Actualmente, los resultados de contacto se basan en simulaciones con entornos lineales simplificados y la identificación se realizó solo en HASELs. Los siguientes pasos incluyen la replicación en hardware real (PAMs, DEAs), la extensión a sistemas multi-DOF y la comparación con políticas basadas en aprendizaje automático.