Equivalent Mechanical Models for Sloshing

Resumen Técnico: Modelos Mecánicos Equivalentes para el Sloshing (Oscilación de Fluidos)

Planteamiento del Problema
El sloshing de propelente sigue siendo un desafío crítico para la estabilidad y el control de las naves espaciales, particularmente en dos regímenes distintos: entornos de microgravedad que requieren una alta estabilidad de apuntamiento (por ejemplo, misiones interplanetarias) y condiciones de alta aceleración ($g$) encontradas durante el lanzamiento, la reentrada y el aterrizaje. Si bien los modelos basados en péndulos y de masa-resorte-amortiguador son herramientas estándar para que los analistas de control evalúen las perturbaciones inducidas por el sloshing, es necesaria una formulación matemática rigurosa que conecte estas representaciones, considerando específicamente los efectos de las fuerzas longitudinales nominales en regímenes de alto-$g$, para un diseño de guía, navegación y control (GNC) robusto.

Metodología
El artículo presenta una derivación rigurosa de modelos mecánicos equivalentes para el sloshing, progresando desde la dinámica no lineal hacia representaciones linealizadas aptas para el análisis de control.

1. Formulación No Lineal:

   • El estudio comienza con un péndulo esférico simple de masa $m_P$ y longitud $\ell$ acoplado a un cuerpo rígido $B$ de masa $m_B$.
   • Utilizando un enfoque Newtoniano, los autores derivan las ecuaciones de movimiento acopladas para el sistema de 8 grados de libertad (traslación/rotación del cuerpo rígido + ángulos del péndulo).
   • La derivación maneja explícitamente la fuerza de tensión $N$ y la elimina para reducir las variables del sistema, resultando en una forma matricial compacta $\mathcal{M}\ddot{\mathbf{x}} = \mathbf{u}$.
   • La formulación se extiende para acomodar múltiples péndulos, representando múltiples tanques o modos de sloshing, utilizando una matriz de masa con estructura de bloques para explotar la dispersión (sparsity) para una solución numérica eficiente.
   • Se incorpora un modelo de amortiguamiento viscoso mediante la adición de términos proporcionales a las velocidades angulares, los cuales afectan la dinámica del péndulo sin alterar directamente las ecuaciones del cuerpo rígido.

2. Linealización y Análisis Modal:

   • Reconociendo que los sistemas de control a bordo típicamente estiman el centro de masa asumiendo una distribución de propelente rígida, los autores reformulan la dinámica con relación a un centro de masa global nominal ($\bar{G}$) en lugar del centro de masa del cuerpo rígido.
   • El sistema se linealiza alrededor de un equilibrio caracterizado por una actitud constante alineada con el marco inercial y una fuerza longitudinal nominal prescrita $\bar{F}_z$.
   • Se analizan dos casos: fuerzas definidas en el marco inercial y fuerzas definidas en el marco del cuerpo fijo.
   • A través de transformaciones de coordenadas y manipulación algebraica, la dinámica linealizada se traslada a una "forma modal". Esto implica definir nuevas coordenadas modales ($\eta_\theta, \eta_\phi$) que desacoplan la dinámica del péndulo en dos osciladores planos independientes.
   • El sistema resultante se expresa como un sistema lineal MIMO de segundo orden ($M\delta\ddot{\mathbf{x}} + K\delta\mathbf{x} = \delta\mathbf{u}_{ext}$), identificando explícitamente la matriz de participación modal y los términos de rigidez inducidos por la fuerza nominal.

3. Equivalencia con Modelos de Masa-Resorte-Amortiguador:

   • El artículo deriva las ecuaciones linealizadas para un sistema de masa-resorte-amortiguador donde la masa de sloshing está restringida a moverse lateralmente mediante resortes de rigidez $k$.
   • Al establecer que la rigidez del resorte $k = m_P \omega_0^2$ (donde $\omega_0$ es la frecuencia natural derivada del modelo de péndulo), los autores demuestran que la dinámica linealizada de masa-resorte-amortiguador es matemáticamente idéntica a la dinámica linealizada del péndulo.
   • El análisis destaca que la presencia de una fuerza longitudinal nominal $\bar{F}_z$ es esencial para que el péndulo exhiba un comportamiento oscilatorio análogo al de un resorte; en ausencia de esta fuerza ($\bar{F}_z=0$), el modelo de péndulo requiere una modificación (por ejemplo, añadir un tercer grado de libertad modal) para representar el sloshing en 0-g.

Resultados Clave

• Equivalencia Matemática: El estudio establece formalmente la equivalencia exacta entre el modelo de péndulo linealizado y el modelo clásico de masa-resorte-amortiguador bajo condiciones de alto-$g$, aclarando cómo la fuerza nominal actúa como el mecanismo de restauración en la formulación del péndulo.
• Representación Modal: La derivación produce una representación modal limpia donde los modos de sloshing se desacoplan de la dinámica del cuerpo rígido en la matriz de rigidez, con el acoplamiento apareciendo solo en la matriz de masa a través de la matriz de participación modal.
• Validación:
  • Validación No Lineal: El modelo no lineal fue validado contra una implementación de MATLAB Simscape Multibody (cubo rígido con cuatro péndulos). Las simulaciones en el dominio del tiempo mostraron una excelente concordancia, con errores máximos inferiores al 0.002% de la amplitud de la respuesta, con discrepancias atribuidas únicamente a las diferencias en los esquemas de integración numérica.
  • Validación Lineal: Los diagramas de Bode que comparan el modelo linealizado analítico con la planta linealizada de Simscape demostraron una excelente concordancia en todo el rango de frecuencias para entradas de fuerza tanto en el marco inercial como en el del cuerpo, con menores discrepancias de baja frecuencia en el caso del marco inercial atribuidas a los límites de precisión numérica en el proceso de linealización de Simscape.

Significancia y Reivindicaciones
El artículo afirma proporcionar una base matemática rigurosa para modelar el sloshing de propelente en entornos de alta aceleración. Su contribución principal es la derivación explícita de las ecuaciones de movimiento no lineales para un sistema de múltiples péndulos y la subsiguiente linealización que tiene en cuenta las fuerzas longitudinales nominales. Al probar la equivalencia matemática entre los modelos de péndulo y de masa-resorte-amortiguador, el trabajo valida el uso de este último para el diseño de control, fundamentándolo en la dinámica de péndulo más fundamental. Los autores enfatizan que los modelos propuestos están bien planteados para aplicaciones realistas de naves espaciales donde los ángulos de sloshing permanecen dentro de límites razonables, y la validación confirma la corrección de las formulaciones derivadas tanto para la simulación en el dominio del tiempo como para el análisis en el dominio de la frecuencia.