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

Sintesi Tecnica: Riduzione dell'Ordine del Modello di una Trave Scorrevole mediante una Base Globale

Enunciato del Problema
Le strutture telescopiche, onnipresenti in macchinari pesanti come gru e telehandler, si affidano a travi scorrevoli le cui lunghezze efficaci e condizioni al contorno (BC) cambiano continuamente. La modellazione di questi sistemi presenta una sfida fondamentale: le tecniche classiche di Riduzione dell'Ordine del Modello (MOR) assumono topologie fisse, calcolando le forme modali una sola volta per una configurazione statica. Nelle travi scorrevoli, man mano che la struttura si estende o si ritrae, la lunghezza efficace e le BC evolvono, rendendo invalide le coordinate modali fisse. L'aggiornamento della base modale durante la simulazione per catturare tali cambiamenti richiede tipicamente il ricalcolo delle coordinate, il che nega i vantaggi in termini di efficienza computazionale derivanti dalla riduzione. Sebbene la PMOR (Riduzione Parametrica dell'Ordine del Modello) offra soluzioni, gli approcci locali esistenti (che interpolano tra valori di parametri discreti) spesso producono modelli ridotti che cambiano con i parametri, complicando l'interpretazione fisica delle coordinate. Inoltre, gli approcci di PMOR globale, che proiettano il sistema su un sottospazio fisso, hanno ricevuto limitata attenzione nel contesto della dinamica multibody (MBD) vincolata per sistemi scorrevoli.

Metodologia
Gli autori propongono un approccio PMOR globale basato su snapshot integrato all'interno di un formalismo multibody a vincoli. La metodologia procede in due fasi principali:

1. Costruzione di una Base di Riduzione Globale:

   • Snapshot: Il sistema è discretizzato utilizzando elementi finiti planari di Euler–Bernoulli (EB). Gli snapshot sono definiti come matrici modali ($\Psi_i$) contenenti un insieme di autovettori ($\phi$) calcolati per specifiche posizioni dello slider ($s$).
   • Parametrizzazione: Due parametri governano la generazione degli snapshot: $n_{snap}$ (numero di posizioni distinte dello slider campionate) e $r_{snap}$ (numero di autovettori mantenuti per ogni snapshot). La condizione al contorno (BC) di scorrimento è imposta algebricamente in nodi specifici durante il calcolo di questi autovettori.
   • Compressione: Tutti gli snapshot sono impilati per formare una matrice "pool" ($\Psi_{all}$). La Decomposizione Ortogonale Propria (POD) tramite Decomposizione ai Valori Singolari (SVD) viene applicata a questo pool per estrarre una compatta base di riduzione globale ($\Psi_{glob}$) composta dai vettori singolari sinistro più significativi.

2. Integrazione nel Formalismo Multibody:

   • Applicazione dei Vincoli: La trave scorrevole è modellata utilizzando un framework Lagrangiano di Tipo 1 con vincoli algebrici. Il vettore di vincolo $g$ impone spostamento nullo all'estremità fissa e spostamento trasversale nullo alla posizione dello slider mobile. La posizione dello slider è tempo-dipendente (reonica) e olofona.
   • Riduzione: Le coordinate complete del sistema $q$ sono espresse come $q(t) = \Psi_{glob} p(t)$, dove $p$ sono le coordinate modali ridotte. Questa proiezione è applicata alle Equazioni Differenziali Algebriche (DAE).
   • Implementazione: Le matrici di massa ($M_{red}$) e rigidezza ($K_{red}$) ridotte sono costanti. Tuttavia, la matrice Jacobiana del vincolo ridotta ($C_{q,red}$) e il vettore di forza ($f_{red}$) sono ricalcolati ad ogni passo temporale per tenere conto dello slider in movimento e dei carichi esterni. La Jacobiana del vincolo è ridotta tramite $C_{q,red} = \Psi_{glob}^T C_q$.

Contributi Chiave

• Base Globale nella MBD Vincolata: Il lavoro dimostra l'integrazione di una base di riduzione globale in un formalismo multibody vincolato per travi scorrevoli, una combinazione non valutata sistematicamente nella letteratura precedente.
• Formulazione del Vincolo Algebrico: Gli autori formulano i vincoli di scorrimento algebrici direttamente utilizzando funzioni di forma di Hermite all'interno del framework EB, garantendo la continuità $C^1$ mentre lo slider si muove tra gli elementi.
• Valutazione Sistematica: Uno studio esaustivo investiga l'impatto della quantità di snapshot ($n_{snap}$) e del numero di autovettori per snapshot ($r_{snap}$) sull'accuratezza, in particolare sotto un budget fisso di modi globali ($r_{glob}$).
• Benchmarking: Il metodo proposto è validato rispetto a una soluzione di riferimento ad alta fedeltà utilizzando la Formulazione delle Coordinate Nodali Assolute (ANCF) con un'implementazione di giunto scorrevole, confrontando sia le risposte cinematiche che la conservazione dell'energia.

Risultati

• Accuratezza: In un caso di test impegnativo che coinvolge una trave altamente flessibile con carico tempo-dipendente e rapido movimento dello slider, la base di riduzione globale (utilizzando $n_{snap}=5, r_{snap}=6, r_{glob}=20$) ha mantenuto un errore di spostamento RMS normalizzato inferiore al 2% rispetto alla soluzione EB non ridotta. Nello specifico, l'errore è stato dell'1,82% per lo spostamento della trave e dell'1,57% per lo spostamento della punta nel caso di test "evoluto".
• Conservazione dell'Energia: L'evoluzione dell'energia del modello EB ridotto ha corrisponduto al benchmark ANCF, confermando che l'applicazione del vincolo algebrico non introduce né dissipazione né iniezione di energia artificiale.
• Efficienza Computazionale: La riduzione dello spazio delle coordinate da 153 gradi di libertà a 20 ha comportato una riduzione del tempo di calcolo di circa il 90% rispetto all'approccio EB non ridotto. Rispetto al benchmark ANCF, l'approccio ridotto è stato circa il 50% più veloce.
• Sensibilità ai Parametri: Lo studio ha rivelato un compromesso nella selezione degli snapshot. Sotto un budget fisso di modi globali, l'equilibrio tra il numero di snapshot e i modi per snapshot è critico. Troppi pochi snapshot producono una base non informata, mentre troppi snapshot (con pochi modi ciascuno) possono portare a una troncamento inefficiente e a errori più elevati. È stato identificato un prodotto ottimale $n_{snap} \times r_{snap} \approx 26$ per il caso di test specifico.

Significatività e Rivendicazioni
L'articolo sostiene che la base di riduzione globale proposta affronta con successo la sfida della topologia evolutiva nelle travi scorrevoli senza sacrificare il significato fisico delle coordinate modali. Proiettando il sistema su un sottospazio fisso derivato da un pool di snapshot, il metodo consente la simulazione efficiente di sistemi multibody flessibili con confini mobili. Gli autori affermano che il metodo è adatto alla simulazione in tempo reale di strutture telescopiche, offrendo un incremento significativo della velocità pur mantenendo un'accuratezza accettabile (errori < 2%). Il lavoro stabilisce una base per l'applicazione della PMOR globale all'interno della dinamica multibody vincolata, colmando il divario tra le tecniche di riduzione classiche e la complessità dei sistemi meccanici dipendenti dai parametri. Si nota come lavoro futuro l'estensione di questo approccio alle formulazioni di frame di riferimento flottante per gestire grandi movimenti di corpo rigido.