Equivalent Mechanical Models for Sloshing

Sintesi Tecnica: Modelli Meccanici Equivalenti per lo Sloshing

Definizione del Problema
Lo sloshing del propellente rimane una sfida critica per la stabilità e il controllo dei veicoli spaziali, in particolare in due regimi distinti: ambienti di microgravità che richiedono un'elevata stabilità di puntamento (ad es., missioni interplanetarie) e condizioni di accelerazione elevata ($g$) incontrate durante il lancio, il rientro e l'atterraggio. Sebbene i modelli basati su pendolo e massa-molla-smorzatore siano strumenti standard per gli analisti di controllo per valutare le perturbazioni indotte dallo sloshing, è necessaria una rigorosa formulazione matematica che colleghi tali rappresentazioni, tenendo conto specificamente degli effetti delle forze longitudinali nominali nei regimi ad alta-$g$, per una progettazione robusta della guida, navigazione e controllo (GNC).

Metodologia
Il documento presenta una rigorosa derivazione di modelli meccanici equivalenti per lo sloshing, progredendo dalla dinamica non lineare alle rappresentazioni linearizzate adatte all'analisi del controllo.

1. Formulazione Non Lineare:

   • Lo studio inizia con un singolo pendolo sferico di massa $m_P$ e lunghezza $\ell$ attaccato a un corpo rigido $B$ di massa $m_B$.
   • Utilizzando un approccio Newtoniano, gli autori derivano le equazioni del moto accoppiate per il sistema a 8 gradi di libertà (traslazione/rotazione del corpo rigido + angoli del pendolo).
   • La derivazione gestisce esplicitamente la forza di tensione $N$ ed elimina tale variabile per ridurre le variabili del sistema, risultando in una forma matriciale compatta $\mathcal{M}\ddot{\mathbf{x}} = \mathbf{u}$.
   • La formulazione è estesa per accomodare più pendoli, rappresentando più serbatoi o modi di sloshing, utilizzando una matrice di massa a struttura a blocchi per sfruttare la sparsità per una soluzione numerica efficiente.
   • Un modello di smorzamento viscoso viene incorporato aggiungendo termini proporzionali alle velocità angolari, che influenzano la dinamica del pendolo senza alterare direttamente le equazioni del corpo rigido.

2. Linearizzazione e Analisi Modale:

   • Riconoscendo che i sistemi di controllo di bordo tipicamente stimano il centro di massa assumendo una distribuzione del propellente rigida, gli autori riformulano la dinamica rispetto a un centro di massa globale nominale ($\bar{G}$) piuttosto che rispetto al centro di massa del corpo rigido.
   • Il sistema viene linearizzato attorno a un equilibrio caratterizzato da un'assetto costante allineato con il sistema di riferimento inerziale e una forza longitudinale nominale prescritta $\bar{F}_z$.
   • Vengono analizzati due casi: forze definite nel sistema di riferimento inerziale e forze definite nel sistema di riferimento del corpo fisso.
   • Attraverso trasformazioni di coordinate e manipolazioni algebriche, la dinamica linearizzata viene ricondotta a una "forma modale". Ciò comporta la definizione di nuove coordinate modali ($\eta_\theta, \eta_\phi$) che disaccoppiano la dinamica del pendolo in due oscillatori planari indipendenti.
   • Il sistema risultante è espresso come un sistema MIMO lineare del secondo ordine ($M\delta\ddot{\mathbf{x}} + K\delta\mathbf{x} = \delta\mathbf{u}_{ext}$), identificando esplicitamente la matrice di partecipazione modale e i termini di rigidezza indotti dalla forza nominale.

3. Equivalenza con i Modelli Massa-Molla-Smorzatore:

   • Il documento deriva le equazioni linearizzate per un sistema massa-molla-smorzatore in cui la massa soggetta a sloshing è vincolata a muoversi lateralmente tramite molle di rigidezza $k$.
   • Impostando la rigidezza della molla $k = m_P \omega_0^2$ (dove $\omega_0$ è la frequenza naturale derivata dal modello del pendolo), gli autori dimostrano che la dinamica linearizzata massa-molla-smorzatore è matematicamente identica alla dinamica linearizzata del pendolo.
   • L'analisi evidenzia che la presenza di una forza longitudinale nominale $\bar{F}_z$ è essenziale affinché il pendolo esibisca un comportamento oscillatorio analogo a una molla; in assenza di questa forza ($\bar{F}_z=0$), il modello del pendolo richiede una modifica (ad es., l'aggiunta di un terzo grado di libertà modale) per rappresentare lo sloshing in assenza di gravità (0-g).

Risultati Chiave

• Equivalenza Matematica: Lo studio stabilisce formalmente l'esatta equivalenza tra il modello del pendolo linearizzato e il classico modello massa-molla-smorzatore sotto condizioni di alta-$g$, chiarendo come la forza nominale agisca come meccanismo di ripristino nella formulazione del pendolo.
• Rappresentazione Modale: La derivazione produce una pulita rappresentazione modale dove i modi di sloshing sono disaccoppiati dalla dinamica del corpo rigido nella matrice di rigidezza, con l'accoppiamento che appare solo nella matrice di massa tramite la matrice di partecipazione modale.
• Validazione:
  • Validazione Non Lineare: Il modello non lineare è stato validato contro un'implementazione MATLAB Simscape Multibody (hub rigido con quattro pendoli). Le simulazioni nel dominio del tempo hanno mostrato un eccellente accordo, con errori massimi inferiori allo 0,002% dell'ampiezza della risposta, con discrepanze attribuite esclusivamente alle differenze negli schemi di integrazione numerica.
  • Validazione Lineare: I diagrammi di Bode che confrontano il modello analitico linearizzato con il plant linearizzato di Simscape hanno dimostrato un eccellente accordo nell'intero intervallo di frequenza sia per gli input di forza nel sistema inerziale che in quello del corpo fisso. Lievi discrepanze a bassa frequenza nel caso del sistema di riferimento inerziale sono state attribuite ai limiti di accuratezza numerica nel processo di linearizzazione di Simscape.

Significatività e Rivendicazioni
Il documento rivendica di fornire una base matematica rigorosa per la modellazione dello sloshing del propellente in ambienti ad alta accelerazione. Il suo contributo primario è la derivazione esplicita delle equazioni del moto non lineari per un sistema multi-pendolo e la successiva linearizzazione che tiene conto delle forze longitudinali nominali. Dimostrando l'equivalenza matematica tra il modello del pendolo e il modello massa-molla-smorzatore, il lavoro valida l'uso di quest'ultimo per la progettazione del controllo, fondandolo però sulla dinamica più fondamentale del pendolo. Gli autori sottolineano che i modelli proposti sono ben definiti per applicazioni spaziali realistiche in cui gli angoli di sloshing rimangono entro limiti ragionevoli, e la validazione conferma la correttezza delle formulazioni derivate sia per la simulazione nel dominio del tempo che per l'analisi nel dominio della frequenza.