Attitude dynamics: angoli di Eulero e quaternioni senza confusione

Attitude Dynamics è l’esame in cui la maggior parte degli studenti di ingegneria aerospaziale scopre che “capire l’orientamento di un corpo rigido” è molto più sottile di quanto sembri sulla carta. Il programma mette insieme tre modi diversi di descrivere una rotazione — angoli di Eulero, matrici, quaternioni — e li chiede quasi sempre insieme, nello stesso esercizio, perché ognuno ha un vantaggio e un difetto che l’altro copre. Aggiungo a questo le equazioni di Eulero per il corpo rigido, che spiegano perché un satellite (o una racchetta da tennis) può iniziare a ruotare in modo scomposto anche senza nessuna forza esterna. Ti porto attraverso il programma con la logica che uso a lezione con i miei studenti di ingegneria aerospaziale, ti mostro dove si perdono i punti, e risolviamo insieme tre esercizi con i conti fatti fino in fondo.

Cosa chiede davvero l’esame

Angoli di Eulero e sequenze di rotazione. Un orientamento nello spazio si può descrivere con tre angoli successivi attorno ad assi che cambiano ad ogni rotazione (per questo si parla di sequenze, tipicamente 3-2-1 o “yaw-pitch-roll” in aerospaziale: prima imbardata attorno a $z$, poi beccheggio attorno al nuovo $y$, poi rollio attorno al nuovo $x$). La matrice di rotazione totale è il prodotto ordinato delle tre matrici elementari:

L’ordine conta: le rotazioni non commutano, quindi $R_x{R}_y{R}_z\ne R_z{R}_y{R}_x$. Il problema di questa rappresentazione è il gimbal lock: quando l’angolo centrale (il beccheggio $\theta$ nella sequenza 3-2-1) raggiunge $\pm 90°$, il primo e il terzo asse di rotazione si allineano, e il sistema perde un grado di libertà — non riesci più a distinguere imbardata e rollio, solo la loro somma o differenza. Non è un bug del tuo conto, è un limite geometrico intrinseco della rappresentazione a tre angoli.

Quaternioni. Un quaternione $q=q_0+q_{1}i+q_{2}j+q_{3}k$ (o in notazione compatta $q=(q_0,\vec{q})$ con $q_0$ scalare e $\vec{q}$ vettoriale) rappresenta una rotazione senza soffrire di gimbal lock, a costo di essere meno intuitivo. Tre proprietà che l’esame chiede sempre:

• Norma unitaria: un quaternione di rotazione ha $\Vert q\Vert =\sqrt{q_0^2+q_1^2+q_2^2+q_3^2}=1$. Se il tuo risultato non ha norma 1 (a meno di errori di arrotondamento), hai sbagliato qualcosa.
• Composizione di due rotazioni tramite prodotto di Hamilton, $q=q_1\otimes q_2$ — non commutativo, esattamente come le matrici.
• Cinematica, cioè come varia $q$ nel tempo se il corpo ruota con velocità angolare $\vec{\omega }$:

dove $\omega$ qui è il quaternione puro $(0,\vec{\omega })$. Questa equazione è quella che si integra numericamente a bordo di un satellite reale per tenere aggiornato l’assetto stimato.

Equazioni di Eulero per il corpo rigido. Nel sistema di riferimento principale d’inerzia (assi allineati con gli autovettori del tensore d’inerzia, momenti principali $I_1,I_2,I_3$), la dinamica rotazionale senza coppie esterne è:

Da qui nasce il risultato più bello (e più chiesto) del corso: la stabilità della rotazione attorno agli assi principali. Un corpo rigido asimmetrico ($I_1\ne I_2\ne I_3$) ruota in modo stabile attorno all’asse di inerzia massima e a quello di inerzia minima, ma è instabile attorno all’asse intermedio. È il cosiddetto effetto racchetta da tennis (o teorema di stabilità dell’asse intermedio): lancia in aria una racchetta facendola ruotare attorno al manico (asse intermedio) e vedrai che dopo un giro o due comincia a ruotare in modo scomposto, cosa che non succede se la fai ruotare attorno agli altri due assi.

Il corso chiude quasi sempre con cenni su gravity-gradient stabilization (un satellite allungato tende ad allinearsi con l’asse di inerzia minima verso la Terra, per via del gradiente di gravità) e ruote di momento/reazione come attuatori per il controllo attivo d’assetto.

Dove si perdono i punti

1. Sbagliare l’ordine della sequenza di Eulero. $R_x(\varphi )R_y(\theta )R_z(\psi )$ non è $R_z(\psi )R_y(\theta )R_x(\varphi )$. Prima di scrivere il prodotto, scrivi a parole la sequenza fisica (“prima ruoto attorno a $z$, poi attorno al nuovo $y$…”) e traducila in matrici nello stesso ordine.

2. Non riconoscere il gimbal lock quando l’esercizio te lo mette davanti. Se $\theta =90°$ e il sistema di equazioni diventa singolare (o le variabili si combinano solo come $\varphi +\psi$ o $\varphi -\psi$), il professore sta verificando se sai riconoscere il fenomeno, non solo calcolarlo numericamente.

3. Confondere il prodotto di Hamilton con un prodotto scalare o vettoriale qualunque. Il prodotto di quaternioni ha una formula precisa (parte scalare meno il prodotto scalare delle parti vettoriali, più termini incrociati col prodotto vettoriale): non è commutativo e non si “intuisce”, si applica la formula.

4. Dimenticare che le equazioni di Eulero sono scritte nel sistema principale d’inerzia. Se il tensore d’inerzia dato nell’esercizio non è diagonale, va prima diagonalizzato (autovalori = momenti principali, autovettori = assi principali) prima di applicare le equazioni nella forma vista sopra.

5. Dimostrare la stabilità dell’asse intermedio “a memoria” invece che linearizzando. Il professore vuole vedere la linearizzazione delle equazioni di Eulero attorno alla rotazione nominale e il segno del coefficiente che ne risulta, non la frase fatta “l’asse intermedio è instabile per il teorema della racchetta da tennis” senza conti.

Il metodo che uso

1. Fissa sempre la convenzione (sequenza di Eulero, quaternione scalare-primo o scalare-ultimo, sistema principale o no) prima di scrivere la prima equazione: metà degli errori nascono da una convenzione cambiata a metà esercizio.
2. Passa dalla descrizione fisica alla matematica, non il contrario: descrivi a parole la sequenza di rotazioni o la configurazione del corpo, poi scrivi le equazioni.
3. Verifica sempre la norma di un quaternione di rotazione dopo ogni composizione: è un controllo di trenta secondi che intercetta la maggior parte degli errori di conto.
4. Per la stabilità, linearizza sempre attorno alla rotazione nominale e guarda il segno del coefficiente risultante: positivo (esponenziale, instabile), negativo (oscillatorio, stabile).

Tre esercizi risolti

Esercizio 1 — Matrice di rotazione 3-2-1 e gimbal lock

Sequenza $\psi =30°$ (imbardata), $\theta =20°$ (beccheggio), $\varphi =10°$ (rollio). La matrice totale $R=R_x(\varphi )R_y(\theta )R_z(\psi )$ vale (arrotondata a quattro cifre):

Applicata al versore $\hat{x}=(1,0,0{)}^T$ dà $(0,8138,-0,4410,0,3785{)}^T$: la nuova direzione dell’asse $x$ del corpo, vista dal sistema di riferimento inerziale.

Ora ripeti con $\theta =90°$ (stessi $\psi ,\varphi$): il gimbal lock.

Nota la prima riga: $(0,0,-1)$, indipendente sia da $\varphi$ che da $\psi$. E gli elementi rimanenti dipendono solo dalla combinazione $\varphi -\psi$ (o $\varphi +\psi$ a seconda della sequenza), non dai due angoli separatamente: hai perso un grado di libertà. È esattamente il motivo per cui i sistemi di controllo d’assetto reali usano i quaternioni per la propagazione, tenendo gli angoli di Eulero solo come rappresentazione “umana” per l’output.

Esercizio 2 — Composizione e cinematica di quaternioni

Componi una rotazione di $90°$ attorno a $z$ con una di $90°$ attorno a $x$, applicando prima $q_2$ e poi $q_1$, cioè $q=q_1\otimes q_2$. I due quaternioni sono

Con la formula del prodotto di Hamilton il risultato è

Verifica di norma:

Corretto. Convertendo in angolo-asse ($q_0=\cos (\Theta /2)$): $\Theta =2\arccos (0,5)=120°$, attorno all’asse $(1,1,1)/\sqrt{3}$ — la diagonale del cubo. Due rotazioni di $90°$ attorno ad assi ortogonali si compongono in un’unica rotazione di $120°$ attorno alla diagonale: un risultato che con le matrici richiederebbe più passaggi per essere visto altrettanto chiaramente.

Cinematica: un corpo con quaternione corrente $q=(1,0,0,0)$ (orientamento di riferimento) e velocità angolare $\vec{\omega }=(0,1,0,2,-0,05)$ rad/s ha

In quell’istante, $\dot{q}$ è semplicemente metà di $\vec{\omega }$ (con parte scalare nulla) perché $q$ è l’identità; un passo di integrazione numerica userebbe questo $\dot{q}$ per aggiornare $q$ e poi rinormalizzarlo a norma 1, dato che l’integrazione numerica introduce piccoli errori che allontanano $q$ dalla sfera unitaria.

Esercizio 3 — Equazioni di Eulero e stabilità dell’asse intermedio

Corpo rigido con momenti principali $I_1=10$, $I_2=20$, $I_3=25$ kg·m² ($I_1<I_2<I_3$: asse 1 minore, asse 2 intermedio, asse 3 maggiore). Verifica la stabilità di una rotazione nominale $\Omega =5$ rad/s attorno all’asse intermedio (asse 2), e attorno all’asse minore (asse 1).

Rotazione attorno all’asse 2 (intermedio). Poni ${\omega }_2=\Omega$ costante, ${\omega }_1,{\omega }_3$ piccole perturbazioni. Linearizzando le equazioni di Eulero per ${\omega }_1$ e ${\omega }_3$ e derivando una seconda volta si ottiene

Con i numeri dati:

Coefficiente positivo $\Rightarrow$ soluzione esponenziale ${\omega }_1(t)\propto e^{\sqrt{5}t}$: la perturbazione cresce, con tempo caratteristico $1/\sqrt{5}\approx 0,45$ s. Instabile, come previsto dal teorema della racchetta da tennis.

Rotazione attorno all’asse 1 (minore). Stesso procedimento con ${\omega }_1=\Omega$ costante, perturbazioni su ${\omega }_2,{\omega }_3$:

Coefficiente negativo $\Rightarrow$ soluzione oscillatoria, ${\omega }_2(t)\propto \cos (\sqrt{7,5}t)$, con pulsazione $\approx 2,74$ rad/s (periodo di “wobble” di circa 2,3 secondi) che resta limitata nel tempo: stabile. Lo stesso conto, ripetuto per l’asse 3 (maggiore), dà lo stesso segno negativo: entrambi gli assi estremi sono stabili, solo quello intermedio non lo è.

Come lavoriamo nelle lezioni

Con Attitude Dynamics parto sempre dal capire quale rappresentazione l’esercizio ti sta chiedendo e perché — angoli di Eulero per l’intuizione fisica, quaternioni quando serve robustezza numerica o quando l’esercizio parla esplicitamente di assetto vicino al gimbal lock, equazioni di Eulero quando si parla di dinamica libera o stabilità. Lavoriamo sui temi d’esame degli anni precedenti del tuo corso specifico, perché il modo in cui ogni professore presenta convenzioni (ordine delle sequenze, quaternione scalare-primo o scalare-ultimo) cambia da un corso all’altro, ed è lì che si annidano gli errori di distrazione più costosi.

Se stai preparando anche Meccanica Orbitale o Meccanica Razionale nello stesso semestre, guarda la pagina delle ripetizioni universitarie con tutte le simulazioni per corso, oppure quella di fisica per il resto del programma di base. E se un esercizio di quaternioni o di stabilità non ti torna, prenota una diagnosi gratuita: mandami il testo e vediamo insieme dov’è il punto in cui il modello fisico ti sfugge.