Real-time Dynamics of Soft Manipulators with Cross-section Inflation: Application to the Octopus Muscular Hydrostat

Resumen Técnico: Dinámica en Tiempo Real de Manipuladores Blandos con Inflación de la Sección Transversal

Planteamiento del Problema
Los brazos robóticos blandos, inspirados en los hidrostatos musculares biológicos como el del pulpo, poseen una alta destreza e inteligencia incorporada, lo que les permite realizar tareas complejas en entornos no estructurados. Sin embargo, modelar su dinámica presenta desafíos significativos. Si bien la teoría clásica de la varilla de Cosserat modela eficazmente la flexión, torsión, cizallamiento y estiramiento de vigas esbeltas, esta se basa en la suposición de secciones transversales rígidas. Esta suposición no logra capturar las deformaciones en el plano, específicamente la inflación o contracción de la sección transversal, que son críticas para tareas que involucran tracción lateral activa (por ejemplo, el alargamiento del brazo de un pulpo mediante la contracción muscular transversal). Los modelos existentes de orden reducido suelen sacrificar precisión en favor de la velocidad computacional, mientras que los métodos de elementos finitos (FEM) de 3D completos son computacionalmente prohibitivos para el control en tiempo real. Además, las extensiones previas de la teoría de Cosserat para incluir la deformación de la sección transversal se han visto limitadas a menudo a contextos estáticos o han dependido de leyes constitutivas simplificadas que se desvían de la realidad física.

Metodología
Los autores proponen una teoría de la varilla de Cosserat extendida que incorpora la inflación de la sección transversal como una variable dinámica. La metodología se estructura de la siguiente manera:

1. Cinemática y Variables de Deformación: El modelo introduce una variable de deformación adicional, $\rho(s)$, que representa la relación de inflación de la sección transversal de una geometría axisimétrica. Esta variable modifica el vector de posición de los puntos materiales, permitiendo que la sección transversal se expanda o contraiga de manera homogénea. El tensor de gradiente de deformación se deriva para incluir esta inflación, lo que conduce a un tensor de deformación de Green-Lagrangian modificado.
2. Ley Constitutiva: Para describir con precisión el comportamiento viscoelástico de los materiales blandos, la ley constitutiva integra el modelo hiperelástico de Saint-Venant Kirchhoff (SVK) con el modelo viscoso de Kelvin-Voigt. Esta combinación da cuenta tanto del almacenamiento de energía elástica como de la disipación viscosa. Los autores derivan expresiones explícitas para la densidad de energía de deformación y la tasa de disipación, incorporando la nueva variable de inflación y asegurando la compatibilidad con las restricciones de incompresibilidad (relación de Poisson $\approx 0.5$).
3. Formulación de la Dinámica: Las ecuaciones de movimiento se derivan de la conservación del momento lineal, del momento angular y del momento material. El balance de momento material produce una ecuación específica que gobierna la evolución de la relación de inflación $\rho$, acoplada con las ecuaciones estándar de equilibrio de fuerzas y momentos. Las cargas externas, incluyendo la gravedad, la flotabilidad, la resistencia hidrodinámica y la masa añadida, se integran en el modelo. Las cargas de actuación interna se categorizan en transversales (que afectan el tamaño de la sección transversal), longitudinales/oblicuas (que afectan la tensión y la flexión) y cargas puntuales.
4. Solución Numérica (GVS): Para permitir la simulación en tiempo real, los autores adaptan el enfoque de Deformación de Variable Geométrica (GVS, por sus siglas en inglés). Los campos de deformación continuos (incluyendo la nueva variable de inflación) se parametrizan utilizando funciones de base polinómicas (polinomios de Legendre y Hermite) mapeadas a un conjunto finito de coordenadas generalizadas. Esto reduce las ecuaciones diferenciales parciales (PDE) a un sistema de ecuaciones diferenciales ordinarias (ODE), permitiendo una computación recursiva eficiente de la cinemática, los Jacobianos y la dinámica.

Contribuciones Clave

• Teoría de Cosserat Extendida: La derivación de una teoría de la varilla de Cosserat dinámica que incluye explícitamente una variable de deformación de inflación en el plano, abordando las limitaciones geométricas de la teoría clásica respecto a la deformación de la sección transversal.
• Ley Constitutiva Viscoelástica: El desarrollo de una ley constitutiva que combina la hiperelasticidad SVK con el amortiguamiento de Kelvin-Voigt, proporcionando una descripción físicamente fundamentada del comportamiento del material blando bajo grandes deformaciones.
• Método Numérico de Orden Reducido: La adaptación del marco GVS a la teoría extendida, permitiendo soluciones rápidas y precisas para la dinámica de varillas blandas con secciones transversales variables.
• Modelado Integral de Cargas: Un marco unificado que da cuenta de las cargas de actuación transversales, longitudinales y oblicuas, así como de las interacciones ambientales complejas (hidrodinámica).

Resultos y Aplicaciones
El artículo valida el modelo a través de tres aplicaciones de simulación:

1. Ajuste de Rigidez Axial: El modelo demuestra la capacidad de ajustar la rigidez axial mediante la modulación de la actuación transversal. A diferencia del modelo clásico, que no puede resistir la compresión axial sin pandeo, el modelo extendido muestra que la contracción transversal genera una tensión interna, permitiendo que el manipulador resista fuerzas de compresión externa y logre un comportamiento rígido. Los valores de rigidez simulados coinciden con las soluciones analíticas.
2. Alcance del Brazo de Pulpo: El modelo simula el movimiento de alcance de un brazo de pulpo, caracterizado por una onda de co-contracción de los músculos longitudinales y transversales. La simulación reproduce el fenómeno biológico del alargamiento del brazo impulsado por la contracción muscular transversal. Las comparaciones con el modelo de Cosserat clásico y un modelo extendido sin actuación transversal muestran que el modelo completo es necesario para capturar el alargamiento correcto y la preservación del volumen (cambio de volumen mínimo $\approx -2.99\%$ frente a una pérdida significativa en los modelos clásicos).
3. Recogida del Brazo de Pulpo: El modelo simula un movimiento de recogida que implica acortamiento y flexión hacia atrás. Incorpora un músculo oblicuo para prevenir la autocolisión durante la flexión de alta curvatura. La simulación destaca la necesidad de la actuación oblicua para evitar la autointersección irreal, una característica que no capturan los modelos que carecen de estos modos de actuación específicos.

Significación y Reivindicaciones
El artículo afirma que el modelo de varilla de Cosserat extendido propuesto cierra la brecha entre la precisión geométrica requerida para movimientos biológicos complejos y la eficiencia computacional necesaria para el control en tiempo real. Al capturar la inflación en el plano, el modelo ofrece una representación más realista de los hidrostatos musculares, particularmente en lo que respecta a la preservación del volumen y la modulación de la rigidez. Los autores sostienen que este marco permite el estudio de movimientos complejos de robótica blanda, como los comportos de alcance y recogida de los pulpos, que anteriormente eran difíciles de modelar con precisión debido a las limitaciones de las suposiciones de sección transversal rígida. El trabajo proporciona una base para la investigación futura en el control y optimización de manipuladores blandos, aprovechando la naturaleza basada en la energía de la ley constitutiva para estrategias de control de bucle cerrado. Los autores reconocen que el modelo sigue siendo una simplificación, ya que asume secciones transversales axisimétricas y no impone plenamente la invarianza de volumen global, sugiriendo que trabajos futuros podrían involucrar leyes constitutivas más generales y restricciones de volumen estrictas.