Model Order Reduction of a Sliding Beam using a Global Basis: Formulation and Evaluation

Resumen Técnico: Reducción de Orden de Modelo de una Viga Deslizante mediante una Base Global

Planteamiento del Problema
Las estructuras telescópicas, ubicuas en maquinaria pesada como grúas y manipuladores de carga, dependen de vigas deslizantes donde las longitudes efectivas y las condiciones de contorno (BC) cambian continuamente. El modelado de estos sistemas presenta un desafío fundamental: las técnicas clásicas de Reducción de Orden de Modelo (MOR) asumen topologías fijas, computando las formas modales una sola vez para una configuración estática. En las vigas deslizantes, a medida que la estructura se extiende o se retrae, la longitud efectiva y las BC evolucionan, invalidando las coordenadas modales fijas. La actualización de la base modal durante la simulación para capturar estos cambios requiere típicamente la recomputación de las coordenadas, lo que anula las ganancias de eficiencia computacional obtenidas por la reducción. Si bien la PMOR paramétrica ofrece soluciones, los enfoques locales existentes (interpolación entre valores de parámetros discretos) suelen dar lugar a modelos reducidos que cambian con los parámetros, complicando la interpretación física de las coordenadas. Además, los enfoques de PMOR global, que proyectan el sistema sobre un subespacio fijo, han recibido poca atención dentro del contexto de la dinámica de sistemas multibosque de restricciones (MBD) para sistemas deslizantes.

Metodología
Los autores proponen un enfoque de PMOR global basado en instantáneas (snapshots) integrado dentro de un formalismo de multibosque de restricciones. La metodología procede en dos etapas principales:

1. Construcción de una Base de Reducción Global:

   • Instantáneas (Snapshots): El sistema se discretiza utilizando elementos finitos de Euler–Bernoulli (EB) planos. Las instantáneas se definen como matrices modales ($\Psi_i$) que contienen un conjunto de modos propios ($\phi$) computados para posiciones de deslizamiento ($s$) específicas.
   • Parametrización: Dos parámetros gobiernan la generación de las instantáneas: $n_{snap}$ (número de posiciones de deslizamiento distintas muestreadas) y $r_{snap}$ (número de modos propios retenidos por instantánea). La condición de contorno de deslizamiento se impone algebraicamente en nodos específicos durante la computación de estos modos propios.
   • Compresión: Todas las instantáneas se apilan para formar una matriz de "reserva" (pool) ($\Psi_{all}$). Se aplica la Descomposición de Órdenes Propios (POD) mediante la Descomposición en Valores Singulares (SVD) a esta reserva para extraer una base de reducción global compacta ($\Psi_{glob}$) que consiste en los vectores singulares izquierdos más significativos.

2. Integración en el Formalismo de Multibosque:

   • Imposición de Restricciones: La viga deslizante se modela utilizando un marco Lagrangiano de Tipo 1 con restricciones algebraicas. El vector de restricción $g$ impone desplazamiento cero en el extremo fijo y desplazamiento transversal cero en la posición del deslizador móvil. La posición del deslizador es dependiente del tiempo (reonómica) y holonómica.
   • Reducción: Las coordenadas completas del sistema $q$ se expresan como $q(t) = \Psi_{glob} p(t)$, donde $p$ son las coordenadas modales reducidas. Esta proyección se aplica a las Ecuaciones Diferenciales Algebraicas (DAE).
   • Implementación: Las matrices de masa ($M_{red}$) y rigidez ($K_{red}$) reducidas son constantes. Sin embargo, la matriz Jacobiana de restricción reducida ($C_{q,red}$) y el vector de fuerza ($f_{red}$) se recomputan en cada paso de tiempo para contabilizar el movimiento del deslizador y las cargas externas. La matriz Jacobiana de la restricción se reduce mediante $C_{q,red} = \Psi_{glob}^T C_q$.

Contribuciones Clave

• Base Global en MBD de Restricciones: El trabajo demuestra la integración de una base de reducción global en un formalismo de multibosque de restricciones para vigas deslizantes, una combinación que no ha sido evaluada sistemáticamente en la literatura previa.
• Formulación de Restricción Algebraica: Los autores formulan restricciones de deslizamiento algebraicas directamente utilizando funciones de forma de Hermite dentro del marco EB, asegurando la continuidad $C^1$ a medida que el deslizador se mueve entre elementos.
• Evaluación Sistemática: Un estudio exhaustivo investiga el impacto de la cantidad de instantáneas ($n_{snap}$) y el conteo de modos propios por instantánea ($r_{snap}$) sobre la precisión, particularmente bajo un presupuesto fijo de modos globales ($r_{glob}$).
• Benchmarking: El método propuesto se valida contra una solución de referencia de alta fidelidad utilizando la Formulación de Coordenadas Nodales Absolutas (ANCF) con una implementación de junta deslizante, comparando tanto las respuestas cinemáticas como la conservación de la energía.

Resultados

• Precisión: En un caso de prueba desafiante que involucra una viga altamente flexible con carga dependiente del tiempo y movimiento rápido del deslizador, la base de reducción global (usando $n_{snap}=5, r_{snap}=6, r_{glob}=20$) mantuvo un error de raíz cuadrática media (RMS) de desplazamiento inferior al 2% en comparación con la solución EB no reducida. Específicamente, el error fue de 1.82% para el desplazamiento de la viga y 1.57% para el desplazamiento de la punta en el caso de prueba "evolucionado".
• Conservación de la Energía: La evolución de la energía del modelo EB reducido coincidió con el benchmark de ANCF, confirmando que la imposición de la restricción algebraica no introduce ni disipación ni inyección de energía artificial.
• Eficiencia Computacional: La reducción del espacio de coordenadas de 153 grados de libertad a 20 resultó en una reducción del tiempo de computación de aproximadamente el 90% en comparación con el enfoque EB no reducido. Comparado con el benchmark de ANCF, el enfoque reducido fue aproximadamente un 50% más rápido.
• Sensibilidad de Parámetros: El estudio reveló un compromiso (trade-off) en la selección de instantáneas. Bajo un presupuesto fijo de modos globales, el equilibrio entre el número de instantáneas y los modos por instantánea es crítico. Demasiadas pocas instantáneas producen una base mal informada, mientras que demasiadas instantáneas (con pocos modos cada una) pueden conducir a una truncación ineficiente y mayores errores. Se identificó un producto óptimo de $n_{snap} \times r_{snap} \approx 26$ para el caso de prueba específico.

Significancia y Reivindicaciones
El artículo sostiene que la base de reducción global propuesta aborda con éxito el desafío de la evolución de la topología en vigas deslizantes sin sacrificar el significado físico de las coordenadas modales. Al proyectar el sistema sobre un subespacio fijo derivado de una reserva de instantáneas, el método permite la simulación eficiente de sistemas multibosque flexibles con fronteras móviles. Los autores afirman que el método es adecuado para la simulación en tiempo real de estructuras telescópicas, ofreciendo una aceleración significativa manteniendo una precisión aceptable (errores < 2%). El trabajo establece una base para aplicar la PMOR global dentro de la dinámica de multibosque de restricciones, cerrando la brecha entre las técnicas de reducción clásicas y las complejidades de los sistemas mecánicos dependientes de parámetros. Se señala como trabajo futuro la extensión de este enfoque a formulaciones de marco de referencia flotante para manejar grandes movimientos de cuerpo rígido.