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

Résumé Technique : Réduction de l'Ordre de Modèle d'une Poutre Coulissante à l'aide d'une Base Globale

Énoncé du Problème
Les structures télescopiques, omniprésentes dans les machines lourdes telles que les grues et les chariots élévateurs, reposent sur des poutres coulissantes dont les longueurs effectives et les conditions aux limites (CL) varient de manière continue. La modélisation de ces systèmes présente un défi fondamental : les techniques classiques de réduction de l'ordre de modèle (MOR) supposent des topologies fixes, calculant les formes modales une seule fois pour une configuration statique. Dans le cas des poutres coulissantes, à mesure que la structure s'étend ou se rétracte, la longueur effective et les CL évoluent, rendant invalides les coordonnées modales fixes. La mise à jour de la base modale pendant la simulation pour capturer ces changements nécessite généralement de recalculer les coordonnées, ce qui annule les gains d'efficacité computationnelle de la réduction. Bien que la réduction de l'ordre de modèle paramétrique (PMOR) offre des solutions, les approches locales existantes (interpolant entre des valeurs de paramètres discrètes) entraînent souvent des modèles réduits qui changent avec les paramètres, compliquant l'interprétation physique des coordonnées. De plus, les approches de PMOR globales, qui projettent le système sur un sous-espace fixe, ont reçu une attention limitée dans le contexte de la dynamique multi-corps (MBD) à contraintes pour les systèmes coulissants.

Méthodologie
Les auteurs proposent une approche de PMOR globale basée sur des instantanés (snapshots) intégrée dans un formalisme multi-corps à contraintes. La méthodologie se déroule en deux étapes principales :

1. Construction d'une Base de Réduction Globale :

   • Instantanés (Snapshots) : Le système est discrétisé à l'aide d'éléments finis d'Euler-Bernoulli (EB) plans. Les instantanés sont définis comme des matrices modales ($\Psi_i$) contenant un ensemble de modes propres ($\phi$) calculés pour des positions de coulisse ($s$) spécifiques.
   • Paramétrage : Deux paramètres régissent la génération des instantanés : $n_{snap}$ (nombre de positions de coulisse distinctes échantillonnées) et $r_{snap}$ (nombre de modes propres conservés par instantané). La CL de glissement est imposée algébriquement à des nœuds spécifiques lors du calcul de ces modes propres.
   • Compression : Tous les instantanés sont empilés pour former une matrice « pool » ($\Psi_{all}$). Une décomposition en valeurs singulières (SVD) de l'analyse en composantes principales (POD) est appliquée à ce pool pour extraire une base de réduction globale compacte ($\Psi_{glob}$) composée des vecteurs singuliers gauches les plus significatifs.

2. Intégration dans le Formalisme Multi-corps :

   • Imposition des Contraintes : La poutre coulissante est modélisée à l'aide d'un cadre de Lagrange de Type 1 avec des contraintes algébriques. Le vecteur de contrainte $g$ impose un déplacement nul à l'extrémité fixe et un déplacement transversal nul à la position du curseur mobile. La position du curseur est dépendante du temps (rhéonomique) et holonome.
   • Réduction : Les coordonnées complètes du système $q$ sont exprimées par $q(t) = \Psi_{glob} p(t)$, où $p$ représente les coordonnées modales réduites. Cette projection est appliquée aux équations différentielles algébriques (DAE).
   • Implémentation : Les matrices de masse ($M_{red}$) et de rigidité ($K_{red}$) réduites sont constantes. Cependant, le Jacobien de contrainte réduit ($C_{q,red}$) et le vecteur de force ($f_{red}$) sont recalculés à chaque pas de temps pour tenir compte du mouvement du curseur et des charges externes. Le Jacobien de contrainte est réduit via $C_{q,red} = \Psi_{glob}^T C_q$.

Contributions Clés

• Base Globale en MBD à Contraintes : Ce travail démontre l'intégration d'une base de réduction globale dans un formalisme multi-corps à contraintes pour les poutres coulissantes, une combinaison qui n'a pas été systématiquement évaluée dans la littérature antérieure.
• Formulation de Contrainte Algébrique : Les auteurs formulent les contraintes de glissement algébriques directement à l'aide de fonctions de forme de Hermite dans le cadre EB, garantissant une continuité $C^1$ lorsque le curseur se déplace entre les éléments.
• Évaluation Systématique : Une étude exhaustive examine l'impact de la quantité d'instantanés ($n_{snap}$) et du nombre de modes propres par instantané ($r_{snap}$) sur la précision, particulièrement sous un budget fixe de modes globaux ($r_{glob}$).
• Étalonnage (Benchmarking) : La méthode proposée est validée par rapport à une solution de référence à haute fidélité utilisant la formulation des coordonnées nodales absolues (ANCF) avec une implémentation de joint coulissant, comparant à la fois les réponses cinématiques et la conservation de l'énergie.

Résultats

• Précision : Dans un cas de test exigeant impliquant une poutre hautement flexible avec un chargement dépendant du temps et un mouvement rapide du curseur, la base de réduction globale (utilisant $n_{snap}=5, r_{snap}=6, r_{glob}=20$) a maintenu une erreur de déplacement RMS normalisée inférieure à 2 % par rapport à la solution EB non réduite. Plus précisément, l'erreur était de 1,82 % pour le déplacement de la poutre et de 1,57 % pour le déplacement de l'extrémité dans le cas de test « évolué ».
• Conservation de l'Énergie : L'évolution de l'énergie du modèle EB réduit correspond au benchmark ANCF, confirmant que l'imposition de la contrainte algébrique n'introduit ni dissipation ni injection d'énergie artificielle.
• Efficacité Computationnelle : La réduction de l'espace de coordonnées de 153 degrés de liberté à 20 a entraîné une réduction du temps de calcul d'environ 90 % par rapport à l'approche EB non réduite. Comparée au benchmark ANCF, l'approche réduite est environ 50 % plus rapide.
• Sensibilité aux Paramètres : L'étude a révélé un compromis dans la sélection des instantanés. Sous un budget fixe de modes globaux, un équilibre entre le nombre d'instantanés et les modes par instantané est critique. Trop peu d'instantanés produisent une base mal informée, tandis que trop d'instantanés (avec peu de modes chacun) peuvent conduire à une troncature inefficace et des erreurs plus élevées. Un produit optimal de $n_{snap} \times r_{snap} \approx 26$ a été identifié pour ce cas de test spécifique.

Signification et Revendications
L'article affirme que la base de réduction globale proposée répond avec succès au défi de l'évolution de la topologie dans les poutres coulissantes sans sacrifier la signification physique des coordonnées modales. En projetant le système sur un sous-espace fixe dérivé d'un pool d'instantanés, la méthode permet une simulation efficace de systèmes multi-corps flexibles avec des frontières mobiles. Les auteurs affirment que la méthode est adaptée à la simulation en temps réel des structures télescopiques, offrant une accélération significative tout en maintenant une précision acceptable (erreurs < 2 %). Ce travail établit une base pour l'application de la PMOR globale au sein de la dynamique multi-corps à contraintes, comblant le fossé entre les techniques de réduction classiques et la complexité des systèmes mécaniques dépendants des paramètres. Il est noté que les travaux futurs étendront cette approche aux formulations de référentiel flottant pour gérer les grands mouvements de corps rigides.