9. Corpo rigido - Gruppo di fisica dei sistemi complessi di Bologna
by user
Comments
Transcript
9. Corpo rigido - Gruppo di fisica dei sistemi complessi di Bologna
9. Corpo rigido Il corpo rigido è il più significativo tra i sistemi di punti vincolati, poiché descrive molti sistemi macroscopici e presenta alcuni moti integrabili. La cinematica del corpo rigido con un punto fisso è quella del gruppo delle rotazioni; in assenza di vincoli esterni, alle rotazioni si aggiungono le traslazioni, che rendono elicoidale il moto istantaneo. La dinamica è descritta dalle equazioni per la quantità di moto ed il suo momento, dette equazioni cardinali o di Eulero, oppure dalle equazioni di Lagrange e di Hamilton. Le equazioni del moto si scrivono nel sistema di riferimento solidale, perché la distribuzione di massa del corpo, i cui momenti primo e secondo sono il centro di massa e la matrice d’inerzia, non varia. Il momento della quantità di moto e l’energia cinetica dipendono dalla matrice d’inerzia, che è diagonale nel sistema di assi principali, determinati dalle simmetrie del corpo. I moti di precessione di un corpo non soggetto a forze (moto per inerzia) e di un corpo con un asse di simmetria soggetto alla forza peso (moto giroscopico) vengono analizzati attraverso gli integrali primi. Il moto per inerzia viene anche illustrato tramite la elegante costruzione geometrica introdotta da Poinsot. 9.1. CINEMATICA Il corpo rigido è un insieme di punti materiali le cui distanze sono fisse. I vincoli che realizzano questa condizione si dicono interni, o di rigidità. Se non vi sono vincoli esterni, il corpo rigido si dice libero ed ha sei gradi di libertà. Infatti, nota la posizione di tre punti del corpo non allineati, che fissano un sistema di assi cartesiani, la posizione di qualunque altro punto risulta univocamente determinata. Le nove coordinate cartesiane dei tre punti sono vincolate da tre equazioni, che esprimono la costanza delle reciproche c 88-08- 9820 9.1. Cinematica del corpo rigido 195 distanze. Se i punti che formano il corpo rigido sono N , le loro coordinate cartesiane sono 3N e le relazioni che esprimono l’invariabilità delle reciproche distanze sono N2 = N (N2−1) . Queste non sono tra loro indipendenti se N ≥ 5; se P1 , P2 , P3 fissano il riferimento, solo le distanze di ogni successivo punto Pk con k > 3 da P1 , P2 , P3 sono indipendenti e fissano la posizione di Pk , invece le distanze tra Pk e Pj con k, j ≥ 4 risultano già determinate. Quindi i vincoli indipendenti sono ||ri − rj || = cij , ( i = 1, i = 2, i = 3, j = 2, . . . , N j = 3, . . . , N j = 4, . . . , N (9.1.1) ed il loro numero è esattamente 3N − 6. Per determinare la configurazione spaziale del corpo rigido è conveniente scegliere un riferimento cartesiano S ′ solidale con il corpo (ad esempio l’origine O in P1 e gli assi x, y nel piano P1 , P2 , P3 ). I sei parametri necessari per individuare S ′ rispetto a un riferimento fisso S sono le tre coordinate (x0 , y0 , z0 ) di O′ e i tre angoli di Eulero (φ, θ, ψ). I vincoli esterni più comuni sono un punto fisso e un asse fisso; nel primo caso i gradi di libertà diventano tre ed è conveniente scegliere l’origine di S ′ nel punto fisso sicché i tre parametri che individuano S ′ rispetto a S sono gli angoli di Eulero. Nel caso di un asse fisso si ha un sol grado di libertà poiché la posizione del corpo è individuata dall’angolo di rotazione rispetto all’asse. Si considerano anche corpi rigidi che si muovono in modo tale che una loro sezione con un piano fisso si mantenga nel piano; il moto del corpo è determinato da quello della sua sezione cioè di una figura piana nel piano. I gradi di libertà sono tre: le coordinate (x0 , y0 ) dell’origine O′ del sistema solidale S ′ rispetto al sistema fisso S e l’angolo tra gli assi x′ e x dei due sistemi. Se un punto della figura piana è obbligato a muoversi su una curva assegnata, i gradi di libertà si riducono a due. Vincoli esterni di tipo anolonomo sono quelli imposti dalla condizione che un corpo rotoli senza strisciare su di una superficie. Campi di velocità e accelerazione Per un corpo rigido libero lo stato di moto più generale è costituito dalla sovrapposizione di una traslazione e di una rotazione. Essendo S ′ rigidamente connesso al corpo rigido, le posizioni, rispetto ad S ′ , di tutti i suoi punti sono costanti; la velocità dei punti del corpo rigido è quindi semplicemente la velocità di trascinamento v = v0 + ω × (r − r0 ) (9.1.2) Notiamo che v è un campo vettoriale somma di un campo uniforme v0 e di uno solenoidale ω × (r − r0 ), le cui linee di forza sono circonferenze. Chiamiamo asse di rotazione la retta passante per l’origine O′ di S ′ e parallela ad ω. Se si cambia l’origine del sistema solidale S ′ il termine di traslazione cambia, e cambia l’asse di rotazione pur restando parallelo al precedente perché la velocità angolare ω rimane inalterata. 196 c 88-08- 9820 9. Corpo rigido CCCCCCCCC ω CCCCCCCCC CCCCCCCCC CCCCCCCCC CCCCCCCCC CCCCCCCCC CCCCCCCCC v CCCCCCCCC CCCCCCCCC CCCCCCCCC CCCCCCCCC CCCCCCCCC v CCCCCCCCC CCCCCCCCC CCCCCCCCC v CCCCCCCCC CCCCCCCCC CCCCCCCCC CCCCCCCCC CCCCCCCCC CCCCCCCCC ω M vM 0 Figura 9.1.1. r !!! !!! !!! !!! Asse di Mozzi. Se O1 è la nuova l’origine di S ′ e r1 = P − O1 e v1 = v0 + ω × (r1 − r0 ) la sua velocità si ha v = v0 + ω × (r − r1 + r1 − r0 ) = v1 + ω × (r − r1 ) (9.1.3) Il punto O1 può esser scelto in modo che la velocità di traslazione v1 sia parallela all’asse di rotazione. Teorema di Mozzi. Ad ogni istante esiste un asse di rotazione, detto asse di Mozzi, i cui punti hanno velocità parallela all’asse stesso. Rispetto ad un riferimento con l’origine sull’asse di Mozzi il moto istantaneo del corpo rigido è elicoidale. Siano i punti di tale asse dati da rM = r0 + r⊥ + λω, con r⊥ ortogonale ad ω, e sia vM = v0 + ω × r⊥ la velocità di questi punti. La velocità di un punto rispetto ad un sistema S ′ con l’origine su un punto M dell’asse di Mozzi è v = vM + ω × (r − rM ). Per determinare rM imponiamo che vM sia parallelo ad ω: ω × vM = ω × v0 + ω × (ω × r⊥ ) = ω × v0 − ω 2 r⊥ = 0 Ne segue che (9.1.4) ω ω ω × v0 (9.1.5) , v = · v 0 M ω2 ω ω Quindi vM è semplicemente la proiezione di v0 sull’asse di rotazione. Indicando con v⊥ la componente di v0 perpendicolare all’asse di rotazione da (9.1.5) segue che il suo modulo è dato da v⊥ = ωr⊥ Il moto è di tipo elicoidale perché gli archi di traiettoria descritti in intervalli di tempo infinitesimi corrispondono a quelli di un’elica cilindrica, dato che la traslazione avviene lungo l’asse di rotazione, vedi figura 9.1.1. Nel moto di una figura piana sul piano l’asse di Mozzi interseca il piano in un punto C, detto centro istantaneo di rotazione, la cui velocità istantanea è nulla. Rispetto al punto C, il cui vettore posizione è rC = r0 + r⊥ , ha un moto di pura rotazione. Infatti la velocità r⊥ = c 88-08- 9820 9.1. Cinematica del corpo rigido 197 di C è perpendicolare al piano, in quanto parallela all’asse di rotazione, e quindi non può che esser nulla, perché altrimenti la figura uscirebbe dal piano. La velocità rispetto ad un sistema con l’origine in C si scrive quindi v = ω × (r − rC ) (9.1.6) Il punto C cambia da istante a istante descrivendo una curva nel sistema fisso S che si chiama rulletta, e un’altra curva nel sistema solidale S ′ , detta base. Per individuare C basta conoscere le traiettorie di due punti del corpo; infatti r − rC è ortogonale a v e la normale alla traiettoria di ogni punto è una retta che contiene C. Nel caso di un corpo rigido piano, la condizione che questo rotoli senza strisciare su una curva fissa γ è espressa dall’annullarsi della velocità istantanea del punto di contatto C, vedi figura 9.1.2. Questo punto coincide con il centro di istantanea rotazione e la condizione di vincolo si scrive vC = v0 + ω × (rC − r0 ) = 0 (9.1.7) CCCCCCCCC CCCCCCCCC CCCCCCCCC CCCCCCCCC O CCCCCCCCC PCCCCCCCCC 0 CCCCCCCCC CCCCCCCCC C CCCCCCCCC Figura 9.1.2. Corpo che rotola senza strisciare. Come esempio di moto di una figura piana mostriamo nella figura 9.1.3 un’asta AB di lunghezza 2ℓ i cui estremi sono vincolati a muoversi su due guide ortogonali, lungo cui scegliamo gli assi x, y di S: la base è un cerchio di raggio ℓ e centro nel punto medio dell’asta, la rulletta un cerchio di raggio 2ℓ e centro in O. Infatti detto θ l’angolo tra l’asta e l’asse x le coordinate di C nel sistema fisso sono x = 2ℓ cos θ, y = 2ℓ sin θ, mentre nel sistema solidale si ha x′ = ℓ 1 + cos(2θ) , y ′ = ℓ sin(2θ). L’accelerazione del corpo rigido è data dall’accelerazione di trascinamento a = a0 + dω × (r − r0 ) + ω × [ω × (r − r0 )] dt (9.1.8) e, nel caso di pura rotazione uniforme, si riduce alla sola accelerazione centripeta data dall’ultimo termine in (9.1.8) 198 c 88-08- 9820 9. Corpo rigido y y’ CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB B CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB D CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB O CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB A CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB base CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC BBBBBBBBBBBBBBBBB CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC rulletta CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC Figura 9.1.3. x x’ Centro istantaneo di rotazione C per un’asta, rulletta e base. 9.2. SISTEMI DI FORZE Per un corpo rigido, costituito da un sistema di N punti con vincoli di rigidità ed eventuali vincoli esterni, le equazioni del moto per i singoli punti sono date da (6.1.3). Le equazioni del moto per la quantità di moto P ed il suo momento L date da (6.1.8) determinano completamente il moto del corpo rigido e sono note come equazioni cardinali. Per il corpo (e) (e) rigido libero F vinc = 0 e Ω vinc = 0 e le equazioni del moto diventano dP = F(e) dt (9.2.1) dL = Ω (e) dt Insieme con le (8.2.17) e M drG /dt = P costituiscono un sistema di 12 equazioni differenziali del primo ordine nelle coordinate del baricentro, negli angoli di Eulero, nella velocità del baricentro e in quella angolare, poiché L = Jω dove J è un operatore lineare, detto tensore di inerzia la cui rappresentazione nel sistema di assi solidali è una matrice costante, vedi paragrafo 9.3. Quando il corpo è vincolato il numero di equazioni indipendenti si riduce. Se il corpo rigido ha un punto fisso O si sceglie in O l’origine di S ed S ′ ; la reazione (e) vincolare è applicata in O per cui Ω vinc =0e dL = Ω (e) dt (9.2.2) insieme con (8.2.17) costituisce un sistema di 6 equazioni differenziali del primo ordine negli angoli di Eulero e nelle componenti della velocità angolare. c 88-08- 9820 9.2. Sistemi di forze 199 Se il corpo rigido ha un asse fisso z, le reazioni vincolari sono ortogonali all’asse e il (e) momento delle reazioni lungo l’asse è nullo Ωzvinc = 0. L’equazione del moto è dLz = Ω(e) z dt (9.2.3) dove Lz = J φ̇ con φ angolo di rotazione e J momento d’inerzia rispetto all’asse z. Sistemi di forze equivalenti Nelle equazioni del moto compaiono solo la somma vettoriale F e il momento risultante Ω di un sistema di forze applicato al corpo rigido. Due sistemi di forze con gli stessi vettori F ed Ω producono lo stesso moto e si dicono equivalenti. Un sistema di forze con F = Ω = 0 non altera il moto e la sua applicazione o rimozione conduce ad un sistema di forze equivalente a quello dato. Il risultato è intuitivo: il sistema più semplice con F = Ω = 0 è quello di due forze uguali e opposte (F = 0) e con stessa linea di azione (Ω = 0). F4 C F1 F2 F3 F Figura 9.2.1. Forze equivalenti. L’effetto di queste due forze sarebbe quello di allontanare (o avvicinare) i punti P1 e P2 cui sono applicate, con velocità diretta lungo la congiungente: ciò è impedito dalla condizione di rigidità |ri − rj |2 = cij , derivando la quale si ha (vi − vj ) · (ri − rj ) = 0. Di qui seguono tutte le regole relative alla composizione delle forze di un corpo rigido. Come applicazione determiniamo il centro di un sistema di forze parallele, vedi figura 9.2.1. Siano F1 , F2 , . . . , FN le forze parallele applicate a P1 , P2 , . . . , PN con Fk = nFk . Vogliamo determinare il luogo dei punti r rispetto al quale il momento totale si annulla Ω= N X k=1 (rk − r) × nFk = 0, −→ r= N X F k rk k=1 N X k=1 + αn Fk (9.2.4) 200 c 88-08- 9820 9. Corpo rigido y A CCCCCCCCCC x CCCCCCCCCC CCCCCCCCCC z CCCCCCCCCC CCCCCCCCCC CCCCCCCCCC CCCCCCCCCC CCCCCCCCCC C CCCCCCCCCC !!! CCCCCCCCCC !!! CCCCCCCCCC !!! !!! CCCCCCCCCC !!! !!! B CCCCCCCCCC CCCCCCCCCC !!! Figura 9.2.2. F Centro di forze parallele (centrifughe) per un’asta rotante. Per ogni direzione n fissata si ha una retta parallela ad n rispetto ai punti della quale Ω = 0; tutte le rette si incontrano in un punto rC , definito da (9.2.4) per α = 0, che risulta indipendente dall’orientamento delle forze e vien detto centro delle forze parallele. Nel caso della forza peso Fk = mk g tale punto coincide con il centro di massa rG definito da (4.2.1). Poiché rispetto a G il momento è nullo, il sistema delle forze peso è equivalente a una sola forza di modulo F = M g applicata al baricentro G. Calcolo di centri di forza Molti corpi rigidi si possono considerare costituiti da una distribuzione continua di massa con densità assegnata; in questo caso la somma in (9.2.4) o (4.2.1) è sostituita da un integrale di volume (di superficie o di linea se il corpo è una lamina o un filo). Come esempio consideriamo il campo di forze centrifughe su un’asta AB rotante con velocità angolare ω attorno all’asse y passante per l’estremo A, vedi figura 9.2.2. Detta ξ la distanza di un punto P dell’asta da A se θ è l’angolo tra l’asta e l’asse di rotazione, la distanza del punto dall’asse è x = ξ sin θ e la forza centrifuga dF agente su un elemento di lunghezza dξ è dato da dF = ω 2 xdm = ω 2 sin θ ξρ(ξ)dξ e quindi il centro del sistema di forze è dato da Z Z ℓ ℓ ξ 2 ρ(ξ)dξ ξdF ξC = Z0 ℓ 0 dF = Z0 ℓ (9.2.5) ξρ(ξ)dξ 0 dove ρ(ξ) è la densità lineare dell’asta e ℓ la sua lunghezza. Se l’asta è omogenea ξC = 32 ℓ, la forza totale applicata al centro è F = 12 mω 2 ℓ sin θ ed il suo momento Ωz = F ξC cos θ = 1 2 2 6 mω ℓ sin(2θ) . Il lavoro corrispondente è dW = Ωz dθ = −dV /dθ ed il potenziale delle 1 forze centrifughe V = 12 mω 2 ℓ2 (1 − 2 sin2 θ). Notiamo quindi che la lagrangiana dell’asta nel sistema fisso è L = 16 mℓ2 (θ̇2 + ω 2 sin2 θ), se usiamo (9.3.16), (9.4.1) e (8.2.19) con ψ̇ = 0, φ̇ = ω e può essere letta come L = T ′ − V , dove T ′ è l’energia cinetica nel sistema rotante e V il potenziale delle forze centrifughe sopra calcolato. c 88-08- 9820 9.2. Sistemi di forze 201 y BBBBBBBBBBBBBB BBBBBBBBBBBBBB BBBBBBBBBBBBBB BBBBBBBBBBBBBB BBBBBBBBBBBBBB BBBBBBBBBBBBBB BBBBBBBBBBBBBB BBBBBBBBBBBBBB G BBBBBBBBBBBBBB BBBBBBBBBBBBBB BBBBBBBBBBBBBB BBBBBBBBBBBBBB BBBBBBBBBBBBBB BBBBBBBBBBBBBB Figura 9.2.3. x Centro di massa di un disco con un foro. Per il calcolo del centro di massa osserviamo che se il corpo ha un piano di simmetria il centro di massa appartiene al piano. Infatti se x, y è tale piano per ogni punto (x, y, z) si ha un punto speculare (x, y, −z) con la stessa massa e quindi zG = 0 poiché nella somma si hanno contributi a due a due uguali ed opposti. Se un corpo è formato da sottosistemi, di cui conosciamo il centro massa G1 , G2 , . . . e la massa M1 , M2 , . . ., allora il centro di massa totale G è quello del sistema di punti G1 , G2 , . . .. Tale procedimento è utile se il corpo è omogeneo ed ha una cavità: si considera il centro di massa G1 del corpo senza cavità che ha massa M1 e poi il centro di massa G2 della cavità attribuendogli massa negativa M2 dove M = M1 + M2 è la massa del corpo assegnato. Come esempio calcoliamo il centro di massa di un disco omogeneo di raggio R con un foro circolare di raggio r < R/2 il cui centro dista r dal centro del disco, vedi figura 9.2.3. Per simmetria il centro di massa sta sull’asse y passante per i due centri. Scelta l’origine nel centro di massa del disco pieno che ha massa M1 = ρπR2 il centro di massa del foro è in (0, r) e la sua massa è M2 = −ρπr2 . Il centro di massa del disco con il foro è quindi xG = 0, yG = r3 M2 r =− 2 M1 + M2 R − r2 (9.2.6) Equilibrio Un corpo rigido è in equilibrio se vk = 0 e se rk = rk c soddisfano le equazioni F(e) + F vinc Ω (e) + Ω (e) =0 vinc (e) (9.2.7) =0 Si può provare che rk = rk c , vk = 0 sono un’orbita del sistema, cioè scelto rk (0) = rk c , vk (0) = 0 si ha rk (t) = rk c , vk (t) = 0. Infatti dalle equazioni (6.1.8) segue che P(t) = P(0) = 0 e che L(t) = L(0) = 0. Quindi vG = 0 e da (4.2.4) segue che è nullo il momento della quantità di moto L′ rispetto al centro di massa G ove scegliamo l’origine 202 c 88-08- 9820 9. Corpo rigido di S’. Si annulla quindi la velocità angolare poiché L′ = Jω = 0 implica ω = 0 essendo J un operatore lineare non singolare, detto tensore di inerzia, vedi paragrafo 9.3. Come esercizio ricaviamo le condizioni per l’equilibrio usando il principio dei lavori virtuali. Per il corpo rigido libero gli spostamenti virtuali coincidono con gli spostamenti reali infinitesimi δrk = vk δt = δr0 + δα × (r − r0 ) (9.2.8) dove δα = nδα essendo n il versore dell’asse di rotazione e δα l’angolo infinitesimo di rotazione. Il lavoro virtuale è δW = N X k=1 Fk · δrk = δr0 · F + δα · Ω (9.2.9) e poiché i vettori δr0 e δα = nδα sono arbitrari, da δW = 0 segue F = 0, Ω = 0. Se il corpo ha un punto fisso nel quale scegliamo l’origine O del riferimento fisso e di quello solidale, δW = δα · Ω e l’annullarsi del lavoro virtuale implica Ω = 0. Se c’è un asse fisso diretto come n, solo δα è arbitrario l’annullarsi di δW implica n · Ω = 0. Questi risultati sono (e) (e) sono in accordo con (9.2.7) perché se il corpo è libero risulta F vinc = Ω vinc = 0, se vinc (e) vinc (e) ha un punto fisso Ω = 0, se ha un asse fisso, n · Ω = 0. 9.3. TENSORE D’INERZIA Il problema del moto di un corpo rigido va affrontato mediante le equazioni cardinali oppure le equazioni di Lagrange. Occorre quindi esprimere il momento della quantità di moto nel primo caso, l’energia cinetica nel secondo, come funzione delle coordinate che individuano il sistema solidale con il corpo e delle loro derivate prime. Se il corpo rigido è libero scegliamo nel centro di massa G l’origine del sistema solidale, che risulta individuato dalle coordinate di rG e dagli angoli di Eulero. Se la forma differenziale F · rG + Ω · dα, dove esprimiamo dα = ωdt come differenziale negli angoli di Eulero tramite la (8.2.17), risulta esatta, la si scrive come −dV dove V = V (xG , yG , zG , φ, θ, ψ) è il potenziale. La energia cinetica ed il momento della quantità di moto si separano nella componente traslazionale ed in quella rotazionale avendo scelto l’origine nel centro di massa. Infatti detti T ′ e L′ la energia cinetica ed il momento della quantità di moto rispetto al sistema con l’origine in rG e gli assi paralleli al sistema fisso, che rappresentano il contributo rotazionale, valgono le relazioni (4.2.4). D’ora in poi intenteremo che T ed L si riferiscano a un corpo rigido con un punto fisso poiché le stesse considerazioni si applicano a T ′ ed L′ nel caso di un corpo rigido libero. c 88-08- 9820 9.3. Tensore di inerzia 203 Momento d’inerzia Se l’asse di rotazione è fisso ω = nω con n costante, valgono le relazioni elementari T = 1 J(n)ω 2 2 L · n = J(n)ω (9.3.1) dove J(n) è il momento d’inerzia rispetto all’asse J(n) = N X ms ρ2s (9.3.2) s=1 e ρs è la distanza del punto Ps rispetto dall’asse di rotazione. Infatti da vs = kω×rs k = ωρs seguono (9.3.1) e (9.3.2). La proiezione di L sull’asse di rotazione è data da N N L·ω 1X 2T 1X L·n= ms ω · rs × vs = ms ω × rs · vs = = = J(n)ω ω ω s=1 ω s=1 ω (9.3.3) Se l’asse di rotazione non è costante, n e J(n) variano nel tempo. È allora conveniente esprimere T ed L mediante le componenti di ω assi del sistema solidale S ′ , rispetto ai quali la distribuzione delle masse del corpo rigido è costante nel tempo. Matrice di inerzia La velocità di un punto di un corpo rigido rotante con velocità angolare ω è vs = ω × rs = e′x (ωy′ zs′ − ωz′ ys′ ) + e′y (ωz′ x′s − ωx′ zs′ ) + e′z (ωx′ ys′ − ωy′ x′s ) (9.3.4) e l’energia cinetica è una forma quadratica in ωx′ , ωy′ , ωz′ T = 1 2 2 2 [Jxx ωx′ + Jyy ωy′ + Jzz ωz′ + 2Jxy ωx′ ωy′ + 2Jxz ωx′ ωz′ + 2Jyz ωy′ ωz′ ] 2 I coefficienti della forma quadratica formano Jxx J J = xy Jxz la matrice d’inerzia J. Jxy Jxz Jyy Jyz Jyz Jzz (9.3.5) (9.3.6) La matrice J è simmetrica definita positiva: infatti T è positiva e si annulla se e soltanto se ω = 0. Le componenti della matrice di inerzia si calcolano partendo dalla definizione di energia cinetica, ovvero sostituendo (9.3.4) nella definizione di T N N 1X 1X T = ms vs · vs = ms [(ωy′ zs′ − ωz′ ys′ )2 + (ωz′ x′s − ωx′ zs′ )2 + (ωx′ ys′ − ωy′ x′s )2 ] (9.3.7) 2 s=1 2 s=1 204 c 88-08- 9820 9. Corpo rigido e sono date da Jxx = N X 2 2 ms (ys′ + zs′ ) Jxy = − s=1 N X Jyy = Jzz = s=1 N X 2 2 ms (x′s + zs′ ) 2 Jyz = − 2 ms (x′s + ys′ ) Jxz = − s=1 N X s=1 N X s=1 N X ms x′s ys′ ms ys′ zs′ (9.3.8) ms x′s zs′ s=1 Gli elementi diagonali Jxx , Jyy , Jzz sono i momenti d’inerzia rispetto agli assi coordinati. Il momento d’inerzia rispetto a un asse arbitrario individuato dal versore n si ottiene immediatamente ponendo ω = ωn e confrontando (9.3.5) con (9.3.1). 2 2 2 J(n) = Jxx n′x + Jyy n′y + Jzz n′z + 2Jxy n′x n′y + 2Jxz n′x n′z + 2Jyz n′y n′z (9.3.9) L’espressione per L′x si ricava osservando che N N N X X X ∂T ∂vs ′ ′ rs × ms vs = L′x ms vs · ms vs · ex × rs = ex · ′ = ′ = ∂ωx ∂ωx s=1 s=1 s=1 (9.3.10) e analoghe espressioni si hanno per L′y e L′z . Ne segue che le componenti di L sul sistema solidale sono date da ′ ′ ωx Lx L′y = J ωy′ (9.3.11) ′ ′ ωz Lz La matrice di inerzia trasforma le componenti della velocità angolare in quelle del momento della quantità di moto. Indicando con J il corrispondente operatore lineare, detto anche tensore di inerzia, scriveremo L = Jω, T = 1 1 ω · Jω = ω · L 2 2 (9.3.12) Assi principali d’inerzia L’asse x si dice principale d’inerzia se Jxy e Jxz sono nulli; il momento di inerzia rispetto all’asse x che indichiamo con Jx si dice momento principale di inerzia. Se l’asse x è principale di inerzia, la matrice d’inerzia ha la struttura a blocchi Jx 0 J = 0 0 Jyy Jyz 0 Jyz Jzz (9.3.13) c 88-08- 9820 9.3. Tensore di inerzia 205 Analogamente y e z sono assi principali d’inerzia se Jxy = Jyz = 0 e se Jxz = Jyz = 0 rispettivamente; se due assi sono principali di inerzia anche il terzo asse lo è e la matrice J è diagonale. Gli assi principali di inerzia sono gli autovettori della matrice J ed i momenti principali di inerzia Jx , Jy , Jz i corrispondenti autovalori. Gli autovettori di J sono ortonormali poiché J è simmetrica e quindi la trasformazione agli assi principali si ottiene con una rotazione. Se il corpo ha particolari simmetrie gli assi principali si individuano senza bisogno di diagonalizzare la matrice. Ogni asse ortogonale ad un piano di simmetria è principale d’inerzia; un asse di simmetria e gli assi ad esso ortogonali sono principali di inerzia. Infatti se xy è un piano di simmetria, per ogni punto di massa m e coordinate (x, y, z) esiste un punto (x, y, −z) speculare rispetto al piano con la stessa massa; nella somma che definisce Jxz per ogni termine mxz c‘è un termine −mxz di segno opposto e quindi Jxz = 0. Nello stesso modo si prova che Jyz = 0, mentre Jxy 6= 0. Un asse di simmetria è l’intersezione di piani di simmetria. Se z è un asse di simmetria i piani xz e yz sono di simmetria: quindi tutti gli assi sono principali d’inerzia e Jx = Jy . I solidi di rotazione omogenei o con densità che dipende solo dalla distanza dall’asse hanno questa simmetria e si dicono a struttura giroscopica. Figura 9.3.1. Corpo con un piano di simmetria (lato sinistro), con un asse di simmetria (lato destro). Metodi di calcolo Le componenti della matrice d’inerzia sono state qui definite come somme sui punti del corpo rigido. Spesso conviene considerare il corpo come un solido con una distribuzione di massa continua anziché discreta. Allora le somme vengono sostituite da integrali di volume oppure di superficie o di linea se il corpo è una lamina o un filo. Le masse dei punti diventano elementi di massa espressi dal prodotto della densità per l’elemento infinitesimo volume o di superficie o di linea.. 206 c 88-08- 9820 9. Corpo rigido Per semplificare il calcolo della matrice di inerzia è conveniente scegliere inizialmente il sistema di assi principali con l’origine nel centro di massa: infatti se il corpo ha un piano od un asse di simmetria il centro di massa appartiene al piano o all’ asse ed è semplice determinare gli assi principali ed i loro momenti d’inerzia JxG , JyG , JzG . Per un corpo generico senza simmetrie la matrice d’inerzia rispetto ad un qualunque sistema si esprime attraverso questi momenti e le coordinate del centro di massa Jxx = JxG + M (yG2 + zG2 ) Jxy = −M xG yG Jyy = JyG + M (x2G + zG2 ) Jxz = −M xG zG Jzz = JzG + M (x2G + yG2 ) Jyz = −M yG zG (9.3.14) che seguono dalle equazioni di trasformazione delle coordinate xs = xG + x′s notando che P ′ s ms xs = 0 e da relazioni analoghe per y e z. Teorema di Steiner. Il momento di inerzia J(n) rispetto ad un asse n è uguale al momento di inerzia J G (n) rispetto ad un asse parallelo passante per il baricentro, sommato con il prodotto tra la massa totale ed il quadrato della distanza ℓ tra i due assi J(n) = J G (n) + M ℓ2 (9.3.15) Il risultato, evidente se n è uno degli assi coordinati. In generale si prova sostituendo (9.3.14) in (9.3.9) che dà J(n) = J G (n) + M kn × rG k2 . Per un corpo omogeneo con una cavità è conveniente calcolare il momento di inerzia del corpo senza cavità, quello della cavità con la stessa densità di massa presa con segno negativo e sommare i due risultati, seguendo lo stesso procedimento usato per il calcolo del centro di massa. A titolo di esempio calcoliamo i momenti di un’asta omogenea lunga ℓ di massa m rispetto ad un sistema con l’origine nel centro di massa e l’asse z orientato lungo l’asta. Gli assi sono principali di inerzia ed essendo la densità lineare dell’asta ρ = m/ℓ Jx = Jy = Z ℓ/2 z 2 ρdz = −ℓ/2 mℓ2 , 12 Jz = 0 (9.3.16) Rispetto ad un sistema con l’origine in un estremo dell’asta, i momenti si ottengono col teorema di Steiner e valgono Jx = Jy = mℓ2 /3. Per una lamina omogenea quadrata di lato ℓ e massa m consideriamo il riferimento cartesiano con origine nel centro di massa ed assi x, y paralleli ai lati. Per la simmetria gli assi coordinati sono principali di inerzia e si ha Jx = Jy = mℓ2 /12, Jz = Jx + Jy . Per una sfera i tre momenti di inerzia rispetto agli assi cartesiani con origine nel centro sono uguali Jx = Jy = Jz . Se la densità di massa è ρ = 3m(4πR3 )−1 si ha usando coordinate polari Jx + Jy + Jz = 8π Z 0 R ρ r4 dr = 8π 5 6 ρR = mR2 5 5 (9.3.17) c 88-08- 9820 9.4. Moto per inerzia 207 e quindi Jx = 52 mR2 . Per un cerchio omogeneo di raggio R e massa m i momenti di inerzia rispetto ad un sistema di assi baricentrico con l’asse z perpendicolare al piano del cerchio sono m Jx = πR2 Z 0 2π dφ Z R r2 sin2 φ rdr = 0 mR2 , 4 Jy = Jx , Jz = Jx + Jy = 2Jx (9.3.18) 9.4. MOTO PER INERZIA Scegliendo come assi del riferimento solidale con il corpo rigido quelli principali di inerzia, l’energia cinetica e il momento della quantità di moto risultano quindi espressi da T = 1 2 2 2 (Jx ωx′ + Jy ωy′ + Jz ωz′ ) 2 L′x = Jx ωx′ L′y = Jy ωy′ (9.4.1) L′z = Jz ωz′ Le equazioni del moto (9.2.2) proiettate sugli assi del sistema solidale, si scrivono utilizzando, per il vettore L, la trasformazione (8.3.3) tra le derivate temporali, che vale per un vettore generico ′ dL dL = +ω×L (9.4.2) dt dt Equazioni di Eulero Le equazioni del moto (9.4.2) proiettate sul sistema di assi solidale sono note come equazioni di Eulero. Jx ω̇x′ + ωy′ ωz′ (Jz − Jy ) = Ω′x Jy ω̇y′ + ωx′ ωz′ (Jx − Jz ) = Ω′y (9.4.3) Jz ω̇z′ + ωx′ ωy′ (Jy − Jx ) = Ω′z Insieme con le equazioni (8.2.17) costituiscono un sistema di sei equazioni differenziali del primo ordine nelle variabili ωx′ , ωy′ , ωz′ , φ, θ, ψ. Analoghe equazioni si ottengono per il corpo rigido non vincolato; infatti, scelta l’origine del sistema solidale nel baricentro G, si ha dv M G = F(e) dt dL′ dt ′ + ω × L′ = ΩG(e) (9.4.4) dove L′ è il momento angolare nel sistema del centro di massa, vedi (4.2.6). Per corpo rigido con un punto fisso, libero o soggetto a forze con momento Ω(e) nullo, le G ′ equazioni di Eulero (9.4.3), per le componenti ωi della velocità angolare, si disaccoppiano 208 c 88-08- 9820 9. Corpo rigido da quelle per gli angoli (8.2.18) e si integrano per quadrature. Ciò è possibile perché l’energia cinetica 1 2 2 2 E = (Jx ωx′ + Jy ωy′ + Jz ωz′ ) (9.4.5) 2 ed il momento della quantità di moto L sono integrali primi del moto. La conservazione di L2 si scrive 2 2 2 L2 = Jx2 ωx′ + Jy2 ωy′ + Jz2 ωz′ , (9.4.6) La curva che ω descrive sulla intersezione dei due elissoidi (9.4.5) e (9.4.6) solidali con il corpo si parametrizza esprimendo ωy′ e ωz′ in funzione di ωx′ cioè 2 2 2 2 Jy (Jz − Jy )ωy′ = 2EJz − L2 − Jx (Jz − Jx )ωx′ Jz (Jy − Jz )ωz′ = 2EJy − L2 − Jx (Jy − Jx )ωx′ (9.4.7) La prima equazione (9.4.3) con Ω′x = 0 si integra per separazione di variabili e la soluzione, esprimibile attraverso le funzioni ellittiche, risulta periodica in t. In modo analogo si trovano le soluzioni per ωy (t) e ωz (t). Gli angoli di Eulero si ottengono tramite una sola quadratura (per l’angolo φ), se scegliamo l’asse z del sistema fisso diretto come L; in questo caso le componenti di L sul sistema solidale sono L′x = L sin θ sin ψ = Jx ωx′ L′y = L sin θ cos ψ = Jy ωy′ (9.4.8) L′z = L cos θ = Jz ωz′ poiché θ è l’ angolo polare e π/2 − ψ l’ angolo azimultale dell’asse z rispetto al sistema solidale. Da (9.4.8) si ottengono θ e ψ cos θ = Jz ′ ω , L z tan ψ = Jx ωx′ Jy ωy′ (9.4.9) e l’equazione per la φ si ha invertendo (8.2.17) ωx′ sin ψ + ωy′ cos ψ dφ = dt sin θ (9.4.10) Moltiplicando in (9.4.10) il numeratore ed il denominatore per L2 sin θ ed usando (9.4.8) si ottiene 2 2 Jx ω ′x + Jy ωy′ dφ =L 2 ′2 (9.4.11) dt Jx ω x + Jy2 ωy′ 2 Stabilità degli assi principali Tra le soluzioni delle equazioni di Eulero per il moto libero ci sono quelle di equilibrio ωx′ = ω, ωy′ = ωz′ = 0 e le altre due ottenute per permutazione ciclica degli indici che c 88-08- 9820 9.4. Moto per inerzia 209 corrispondono a rotazioni con velocità angolare costante attorno agli assi principali. Un problema che si pone è quello della stabilità lineare degli assi. Consideriamo allora una soluzione vicina a quella di equilibrio in cui le velocità angolari sono ωx′ + ω, ωy′ , ωz′ dove ωx′ , ωy′ , ωz′ sono cosı̀ piccole da poter trascurare i termini quadratici. Le equazioni di Eulero diventano Jx ω̇x′ = 0 Jy ω̇y′ + ωz′ ω(Jx − Jz ) = 0 (9.4.12) Jz ω̇z′ + ωy′ ω(Jy − Jx ) = 0 La soluzione è periodica e l’orbita nel piano ωx′ , ωy′ è una ellisse se (Jx − Jy )(Jx − Jz ) > 0. Se invece (Jx − Jy )(Jx − Jz ) < 0 l’orbita è una iperbole e la soluzione è esponenziale in t. Infatti derivando la terza equazione (9.4.12) e sostituendo nella seconda si ha ω2 (Jx − Jz )(Jx − Jy ) = 0 Jy ω¨y′ + ω˙z′ (Jx − Jz ) = Jy ω¨y′ + ωy′ Jz (9.4.13) Permutando ciclicamente gli indici si ottengono le condizioni di stabilità per le rotazioni attorno agli assi y e z: se Jx < Jy < Jz le rotazioni attorno agli assi x′ , z ′ sono stabili mentre la rotazione attorno all’asse y ′ è instabile. z’ ω m r y’ x’ Figura 9.4.1. Elissoide di inerzia. È possibile determinare le traiettorie che il vettore ω traccia su un elissoide solidale con il corpo, detto elissoide di inerzia mediante una costruzione geometrica basata sugli invarianti del moto. L’elissoide di inerzia è definito da 2 2 2 F (x′ , y ′ , z ′ ) = Jx x′ + Jy y ′ + Jz z ′ = 1 (9.4.14) La conservazione dell’energia (9.4.5) implica che il vettore ω r= √ 2E (9.4.15) 210 c 88-08- 9820 9. Corpo rigido ′ ωx′ √ωy √ωz′ , 2E , 2E ) nel sistema solidale appartiene all’elissoide, di componenti (x′ , y ′ , z ′ ) = ( √2E vedi figura 9.4.1. Il vettore r individua il punto in cui l’asse di rotazione istantanea interseca l’elissoide e ad ogni istante appartiene anche ad un secondo elissoide, che si ottiene dalla conservazione del momento angolare (9.4.6) 2EJx2 ′ 2 2EJy2 ′ 2 2EJz2 ′ 2 x + y + z =1 L2 L2 L2 (9.4.16) Pertanto la traiettoria dell’asse di rotazione è definita dalla intersezione dei due elissoidi. Vicino alle intersezioni con gli assi x′ e z ′ le intersezioni sono curve chiuse la cui proiezione sui piani ortogonali sono ellissi, mentre nel caso dell’asse y ′ si ottiene una curva la cui proiezione è una iperbole se Jx < Jy < Jz . Ad esempio la proiezione delle orbite prossime ad x′ ottenuta eliminando x′ tra (9.4.14) e (9.4.16) 2 2 (Jx − Jy )Jy y ′ + (Jx − Jz )Jz z ′ = Jx − L2 /(2E) (9.4.17) è una ellisse se (Jx − Jy )(Jx − Jz ) > 0. Costruzione geometrica di Poinsot Dalla soluzione delle equazioni di Eulero sopra delineata non emergono in modo chiaro le caratteristiche del moto per inerzia nel sistema fisso data la complessità della soluzione analitica. Esiste tuttavia una costruzione di natura puramente geometrica, da cui risulta che il corpo rigido è animato da un moto di precessione. Teorema di Poinsot. Nel moto per inerzia di un corpo rigido con un punto fisso l’elissoide di inerzia rotola senza strisciare √ su di un piano fisso la cui normale è diretta come L e la cui distanza dall’origine è d = 2E/L. L’asse di rotazione descrive un cono nel sistema fisso ed un cono nel sistema solidale; l’asse di rotazione è ad ogni istante la semiretta comune ai due coni, che rotolano senza strisciare l’uno sull’altro. Il versore m normale all’ellissoide in r si calcola tenendo conto che L′i ∂F Ji ωi′ ′ √ √ = 2J x = 2 = 2 i i ∂x′i 2E 2E ed è dato da m= L grad F = ||grad F || L (9.4.18) (9.4.19) Il vettore L è costante rispetto al sistema fisso, non in quello solidale. Il modulo L è invece costante in entrambi i sistemi poiché è invariante per rotazione. Si riconosce quindi che la proiezione di r sull’asse avente la direzione di m è costante. c 88-08- 9820 9.5. Moti giroscopici ω 211 z CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC m CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC L CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC z’ CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC P’ P CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC r CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC O y’ x’ Figura 9.4.2. Piano invariante per l’elissoide di inerzia. ω·L = r·m= √ L 2E √ 2E L (9.4.20) Tale prodotto scalare rappresenta la distanza dall’origine del piano tangente all’ellissoide. Durante il moto del corpo rigido l’elissoide d’inerzia si muove in modo solidale con questo poiché i suoi assi coincidono con gli assi principali d’inerzia; inoltre, se il corpo rigido non è soggetto a forze T ed L sono costanti e quindi anche r · m è costante. Perciò si ha un piano fisso, perpendicolare ad L, con cui elissoide è a contatto e che ha dal centro di questo una distanza fissa, vedi figura 9.4.2. Su questo piano, detto invariabile, l’elissoide rotola senza strisciare poiché la velocità del punto di contatto vale v = ω × r = 0. L’asse di rotazione istantaneo è la retta che congiunge il centro dell’elissoide con il suo punto di contatto col piano invariabile. Se il corpo ha un asse di simmetria z ′ , allora l’elissoide è di rotazione attorno a z ′ e il punto di contatto durante il moto descrive una circonferenza sia sul piano invariabile sia sull’elissoide stesso. Perciò la velocità angolare ω nel sistema fisso si sposta descrivendo un cono rotondo (la sezione normale all’asse è una circonferenza); anche nel sistema solidale ω descrive un cono rotondo con l’asse in z ′ , vedi figura 9.4.3. Durante il moto i due coni rotolano senza strisciare l’uno sull’altro. La loro intersezione con la sfera unità sono due cerchi che rotolano senza strisciare ed il punto di contatto si muove di moto uniforme su un cerchio, vedi figura 9.6.1. Il moto è una precessione regolare perché l’asse di rotazione ruota attorno ad una direzione fissa in ciascuno dei due sistemi di riferimento. 212 c 88-08- 9820 9. Corpo rigido L ω z’ Figura 9.4.3. Coni di Poinsot. 9.5. MOTI GIROSCOPICI Il moto per inerzia è particolarmente semplice quando il corpo possiede un asse di simmetria z ′ cioè se ha struttura giroscopica. I coni di Poinsot sono rotondi e le loro intersezioni con la sfera unitaria sono circonferenze che rotolano senza strisciare l’una sull’altra. L’asse di simmetria z ′ ha θ, π/2 − φ come angoli polare ed azimultale rispetto al sistema S mentre ψ è l’angolo di rotazione del corpo attorno a z ′ . La traiettoria descritta dall’asse z ′ sulla sfera unitaria è anch’essa una circonferenza il cui centro è sull’asse avente la direzione di L. La soluzione analitica è semplice in questo caso: infatti dato che Jx = Jy le equazioni di Eulero diventano ′ ′ ′ Jx ω̇x + ωy ωz (Jz − Jx ) = 0 Jx ω̇y′ + ωz′ ωx′ (Jx − Jz ) = 0 (9.5.1) Jz ω̇z′ = 0 L’ultima equazione implica ωz′ = ωz′ (0) e, posto ω = ωz′ (0) Jx − Jz Jx (9.5.2) le prime due equazioni diventano ω̇x′ = ω ωy′ ω̇y′ = −ω ωx′ (9.5.3) La soluzione ωx′ = A sin(ωt + α) e ωy′ = A cos(ωt + α) può essere riscritta nella forma seguente ′ ′ ωx 0 ωx (0) ωy′ = Rz (ωt + α) A = Rz (ωt) ωy′ (0) (9.5.4) ωz′ ωz′ (0) ωz′ (0) e mostra che il vettore ω ruota attorno a z ′ con velocità angolare ω nel sistema solidale. c 88-08- 9820 9.5. Moti giroscopici 213 Per determinare il moto nel sistema fisso, confrontiamo la soluzione ottenuta (9.5.4) con l’espressione (8.2.17) delle componenti ωi′ della velocità angolare in funzione degli angoli di Eulero. Dal confronto non emerge una soluzione unica; l’arbitrarietà risiede nella scelta degli assi del sistema fisso rispetto al vettore costante L. La soluzione più semplice si ottiene scegliendo ψ = ωt + α da cui segue θ̇ = 0, φ̇ sin θ = A, φ̇ cos θ + ψ̇ = ωz′ (0) (9.5.5) In questo caso l’asse z ′ ha inclinazione costante rispetto a z, il corpo ruota con velocità angolare costante ψ̇ attorno all’asse z ′ , che a sua volta ruota con velocità angolare costante φ̇ attorno a z. Con la scelta fatta l’asse z non è qualsiasi ma ha la stessa direzione di L. Per verificarlo è sufficiente considerare la espressione che hanno le componenti di L nel sistema solidale quando l’asse z è diretto lungo L. Da (9.4.8) tenendo conto della condizione di simmetria Jx = Jy segue che ′ ωx 0 ωy′ = Rz (ψ) sin θ L/Jx (9.5.6) ′ ωy cos θ L/Jz Confrontando (9.5.6) e (8.2.17) emerge che θ̇ = 0, φ̇ = L/Jx e φ̇ cos θ + ψ̇ = cos θ L/Jz . Il moto nel riferimento fisso è una precessione come si vede scrivendo ω = φ̇ez + ψ̇e′z dove l’angolo θ tra i due assi ez , e′z è costante e le rispettive componenti, detto L′z = L cos θ sono 1 L 1 ′ − , φ̇ = (9.5.7) ψ̇ = ω = Lz Jz Jx Jx Moti nel campo di gravità Pendolo fisico. È cosı̀ chiamato un corpo rigido con un asse z fisso. Detta ℓ la distanza del baricentro dall’asse e φ l’angolo che il piano passante per l’asse z, che supponiamo orizzontale, e il baricentro forma con il piano verticale yz, vedi figura 9.5.1, la lagrangiana si scrive 1 L = J φ̇2 + M gℓ cos φ (9.5.8) 2 dove J è il momento d’inerzia rispetto all’asse di rotazione. L’equazione del moto J φ̈ + M gℓ sin φ = 0 y z Figura 9.5.1. CCCCCCCCCCCCC CCCCCCCCCCCCC O CCCCCCCCCCCCC CCCCCCCCCCCCC x CCCCCCCCCCCCC CCCCCCCCCCCCC φ G CCCCCCCCCCCCC CCCCCCCCCCCCC CCCCCCCCCCCCC CCCCCCCCCCCCC CCCCCCCCCCCCC CCCCCCCCCCCCC Pendolo fisico: sezione con il piano verticale passante per il baricentro. (9.5.9) 214 c 88-08- 9820 9. Corpo rigido è quella di un pendolo semplice di lunghezza J/M ℓ, che prende il nome di lunghezza ridotta. La equazione (9.5.9) si ottiene anche a partire da (9.2.3) dove le proiezioni sull’asse z del momento della quantità di moto e della forza peso sono dati da Lz = Jz ωz = Jz φ̇ e da Ωz = −M gℓ sin φ. Dal confronto con (9.5.9) segue Ωz = −dV /dφ, relazione identica a (9.2.9), che in questo caso diventa dW = Ωz dφ = −dV . Se l’asse di rotazione non è orizzontale indichiamo con x′ , z ′ gli assi ottenuti ruotando di α rispetto all’asse x gli assi y, z indicati in nella figura 9.5.1. Nel sistema x′ , z ′ , ove z ′ è l’asse di rotazione del pendolo, l’accelerazione di gravità è g = −g cos α j′ + g sin α k′ e la posizione del centro di massa rG = sin φ i − cos φ j′ . Il potenziale della forza peso è V = −M g · rG = −M gℓ cos φ cos α da cui si ricava Ω′z = −dV /dφ e quindi la equazione del moto è (9.5.9) dove g viene moltiplicato per cos α. Trottola di Lagrange. È costituita da corpo rigido con un asse di simmetria z ′ e un punto fisso appartenente a z ′ , vedi figura 9.5.2. La sua lagrangiana è data da 1 Jx 2 Jx 2 2 Jz 1 2 2 2 Jx (ωx′ + ωy′ ) + Jz ωz′ − V = θ̇ + φ̇ sin θ + (ψ̇ + φ̇ cos θ)2 − M gℓ cos θ 2 2 2 2 2 (9.5.10) dove ℓ è la distanza tra il punto fisso e il baricentro e per ωx′ , ωy′ , ωz′ si sono usate le (8.2.17). L= z’ z G ϕ θ y ψ x L Figura 9.5.2. x’ Giroscopio pesante o trottola di Lagrange. Poiché L non dipende da ψ, φ e t i tre integrali primi del moto pψ , pφ , H = T + V riconducono la soluzione a quadrature. Essendo φ e ψ gli angoli di rotazione attorno agli assi z e z ′ , l’invarianza della lagrangiana per rotazione attorno ad essi implica (teorema c 88-08- 9820 9.5. Moti giroscopici 215 di Nöther) che che le proiezioni del momento della quantità di moto L su questi assi sono costanti L′z = pψ = Jz (ψ̇ + φ̇ cos θ) Lz = pφ = Jx φ̇ sin2 θ + Jz cos θ(ψ̇ + φ̇ cos θ) = Jx φ̇ sin2 θ + pψ cos θ (9.5.11) Dalla seconda delle (9.5.11) si ricava φ̇ pφ − pψ cos θ Jx sin2 θ (9.5.12) p2ψ 1 1 (pφ − pψ cos θ)2 + + M gℓ cos θ Jx θ̇2 + 2 2 2Jz Jx sin2 θ (9.5.13) φ̇ = e infine per l’energia si trova E= Introdotta la variabile u = cos θ l’equazione del moto si scrive du dt 2 p2ψ (pφ − upψ )2 2 2 (1 − u ) E − − M gℓu − = P (u) ≡ Jx 2Jz Jx2 (9.5.14) Osserviamo che P (±1) < 0 e che per u → ∞ si ha P (u) ∝ u3 ; quindi P (u) nell’intevallo [−1, 1] ha due zeri distinti oppure uno zero doppio. Detti u± gli zeri e θ± i valori corrispondenti per θ, il moto di θ(t) risulta periodico tra questi due valori. Descriviamo il moto attraverso l’intersezione dell’asse di simmetria z ′ con la sfera unitaria, ove θ e π/2−φ sono le coordinate polari. La traiettoria è compresa tra due cerchi paralleli θ = θ− , θ+ ed il moto in φ è determinato dall’equazione (9.5.12). Perciò se il numeratore non si annulla il moto in φ è monotono e la traiettoria è quella indicata nella figura 9.5.3, mentre se si annulla φ non è più monotona e si ha inversione del moto; nel caso di transizione φ(t) è ancora monotona ma la traiettoria presenta delle cuspidi. L’oscillazione dell’asse z ′ rispetto all’asse z, che è assente nel moto per inerzia, viene detta nutazione. z y x Figura 9.5.3. Traiettorie tracciate sulla sfera dall’asse di simmetria del giroscopio. 216 c 88-08- 9820 9. Corpo rigido 9.6. FORMULAZIONE HAMILTONIANA Prendiamo in esame due formulazioni hamiltoniane per il moto del corpo rigido basate sugli angoli di Eulero e sugli angoli di Deprit. Nel primo caso si ottiene una hamiltoniana esplicitamente integrabile se si sceglie l’asse z diretto come il momento della quantità moto L, nel secondo invece la scelta del sistema fisso può farsi in modo indipendente da L. Hamiltoniano con angoli di Eulero Calcoliamo i momenti coniugati pk degli angoli di Eulero qk esprimendoli attraverso le componenti del momento della quantità di moto L′i nel sistema solidale. Dalla lagrangiana 2 2 2 L = 12 (Jx ωx′ + Jy ωy′ + Jz ωz′ ), tenendo conto che L′k = Jk ωk′ , si ottiene 3 3 X X ∂ω ′ ∂L ∂ωi′ pk = Ji ωi′ L′i i = = ∂ q̇k ∂ q̇ ∂ q̇k k i=1 i=1 (9.6.1) e usando (8.2.17) si ricava pθ = L′x cos ψ − L′y sin ψ pφ = L′x sin θ sin ψ + L′y sin θ cos ψ + L′z cos θ (9.6.2) pψ = L′z Le componenti di L in entrambi i riferimenti si esprimono attraverso gli angoli di Eulero ed i loro momenti coniugati pθ Lx Ly = Rz (−φ)Rx (−θ) pφ − pψ cos θ sin θ Lz pψ (9.6.4) L’hamiltoniana del sistema assume la forma pθ L′x L′y = Rz (ψ) pφ − pψ cos θ , sin θ L′z pψ 2 2 2 2 L′y L′x pφ − pψ cos θ L′z 1 H= + + = sin ψ + pθ cos ψ + 2Jx 2Jy 2Jz 2Jx sin θ 2 p2ψ pφ − pψ cos θ 1 cos ψ − pθ sin ψ + 2Jy sin θ 2Jz (9.6.5) ed il momento della quantità di moto totale è dato da L2 = (pφ − pψ cos θ)2 + p2ψ + p2θ 2 sin θ (9.6.6) c 88-08- 9820 9.6. Formulazione hamiltoniana 217 Corpo con asse di simmetria Se il corpo ha struttura giroscopica, detto z ′ l’asse di simmetria per cui Jx = Jy l’hamiltoniana diventa p2ψ 1 p2ψ 1 (pφ − pψ cos θ)2 L2 1 2 H= = + − (9.6.7) + pθ + 2Jx 2Jz 2Jx 2 Jz Jx sin2 θ Oltre ad H gli integrali primi sono L, pφ pψ . L’angolo θ è periodico e φ è monotono se il cono descritto da z ′ contiene l’asse z, vedi figura 9.6.1; se invece l’asse z è esterno al cono, θ e φ sono periodici con lo stesso periodo. Nel primo caso infatti φ aumenta di 2π mentre θ ritorna al valore iniziale, quando z ′ descrive la circonferenza sulla sfera, nel secondo caso entrambi ritornano al loro valore iniziale. La soluzione θ(t) è una funzione periodica di t che oscilla tra due valori θ− e θ+ , corrispondenti ai due paralleli, che delimitano la traiettoria dell’asse di simmetria. Infatti posto u = cos θ, da (9.6.7) segue u soddisfa l’equazione u̇ = P (u) dove P (u) è un polinomio di grado 2 i cui zeri u± = cos θ± sono i due punti di inversione del moto. L’angolo φ che soddisfa la equazione (9.5.12), è somma di una funzione lineare in t e di una funzione periodica con lo stesso periodo di θ(t). La scelta dell’asse z lungo L rende θ costante e φ lineare in t, eliminando i problemi derivanti dalla parametrizzazione della circonferenza fuori asse. ω z z L z’ z’ ω L Figura 9.6.1. !!!! !!!! !!!! Orbita dell’asse z ′ sulla sfera unitaria: L×ez =0 (lato sinistro), L×ez 6=0 (lato destro). 218 c 88-08- 9820 9. Corpo rigido Allineamento con il momento angolare La scelta dell’asse z lungo L consente una semplice soluzione analitica anche per il moto di un corpo senza simmetria Jx < Jy < Jz . Infatti in questo caso si ha pφ = Lz = L, pψ = L′z = L cos θ e sostituendo pψ = pφ cos θ in (9.6.6) si trova L2 = p2φ + p2θ che implica pθ = 0. Allo stesso risultato si perviene notando che pθ è la proiezione di L lungo la linea dei nodi, ad esso ortogonale. L’hamiltoniana (9.6.5) assume quindi la forma seguente p2φ sin2 θ H= 2 sin2 ψ cos2 ψ + Jx Jy p2ψ p2φ − p2ψ + = 2Jz 2 sin2 ψ cos2 ψ + Jx Jy p2ψ + 2Jz (9.6.8) Siccome pφ è un integrale primo, l’hamiltoniana per ψ è simile a quella di un pendolo e se ne può tracciare il ritratto di fase nel piano ψ, pψ . I punti critici sono pψ = 0, ψ = 0, π/2, π dove il primo ed il terzo sono iperbolici e corrispondono a rotazioni attorno all’asse y ′ , mentre il secondo è ellittico e corrisponde alla rotazione attorno all’asse x′ . Per verificarlo è sufficiente valutare il gradiente di H ∂H = pψ ∂ψ 1 sin2 ψ cos2 ψ − − Jz Jx Jy , ∂H sin(2ψ) 2 = (pφ − p2ψ ) ∂ψ 2 1 1 − Jx Jy (9.6.9) e l’hessiano di H che è diagonale nei punti critici ove vale ∂ 2 H/∂p2ψ = (Jz−1 − Jy−1 ) e ∂ 2 H/∂ψ 2 = p2φ (Jx−1 − Jy−1 ) se ψ = 0, π, mentre è dato da ∂ 2 H/∂p2ψ = (Jz−1 − Jx−1 ) e ∂ 2 H/∂ψ 2 = p2φ (Jy−1 − Jx−1 ) se ψ = π/2. Nel primo caso H ha una sella nel secondo un massimo. Notiamo inoltre che pψ = pφ cos θ = 0 implica θ = π/2 e quindi l’asse z ′ appartiene al piano x, y. Siccome φ è l’unico angolo che varia, la rotazione avviene attorno all’asse z; se ψ = 0 si ha z = y ′ poiché l’asse dei nodi OL coincide con x′ , se invece ψ = π/2, l’asse dei nodi è ortogonale a x′ e quindi z = x′ . Se x fosse l’asse con il momento di inerzia massimo anziché minimo, H avrebbe un minimo anziché un massimo. Quando il corpo ha un asse z ′ di simmetria l’hamiltoniana diventa p2ψ p2φ + H= 2Jx 2 1 1 − Jz Jx (9.6.10) e coincide con (9.6.7). Le frequenze ψ̇ di rotazione attorno all’asse di simmetria, φ̇ di precessione attorno all’asse z, risultano in accordo con (9.5.7), tenuto conto che pφ = L, pψ = L′z . L ’angolo θ e le frequenze φ̇, ψ̇ si esprimono attraverso l’energia E ed il momento angolare totale L ricordando che pψ = L cos θ Jz − 2EL−2 Jz Jx , cos θ = Jz − Jx 2 dove 2EL−2 varia tra Jz−1 e Jx−1 . L φ̇ = . Jx ψ̇ = L cos θ 1 1 − Jz Jx (9.6.11) c 88-08- 9820 9.6. Formulazione hamiltoniana 219 Angoli di Deprit Gli angoli di Deprit α, β, γ che consentono una scelta unica degli assi del sistema fisso, qualunque sia la direzione di L sono definiti nella figura 9.6.2 z L z’ γ BBBBBBBBBBBBBBBBy’ BBBBBBBBBBBBBBBB BBBBBBBBBBBBBBBBQ BBBBBBBBBBBBBBBB CCCCCCCCCCC BBBBBBBBBBBBBBBB O CCCCCCCCCCC BBBBBBBBBBBBBBBB y CCCCCCCCCCC BBBBBBBBBBBBBBBB CCCCCCCCCCC BBBBBBBBBBBBBBBB CCCCCCCCCCC BBBBBBBBBBBBBBBB αBBB CCCCCCCCCCC BBBBBBBBBBBBBBBB β γ CCCCCCCCCCC BBBBBBBBBBBBBBBB BBB CCCCCCCCCCC BBBBBBBBBBBBBBBB BBB CCCCCCCCCCC BBBBBBBBBBBBBBBB BBB CCCCCCCCCCC BBBBBBBBBBBBBBBB CCCCCCCCCCC BBBBBBBBBBBBBBBB M CCCCCCCCCCC BBBBBBBBBBBBBBBB x CCCCCCCCCCC x’ BBBBBBBBBBBBBBBB N L CCCCCCCCCCC BBBBBBBBBBBBBBBB BBBBBBBBBBBy’ BBBBBBBBBBB O O BBBBBBBBBBB CCCCCCCCCCCCCC y BBBBBBBBBBB CCCCCCCCCCCCCC BBBBBBBBBBB CCCCCCCCCCCCCC BBBBBBBBBBB CCCCCCCCCCCCCC CCC φ BBBBBBBBBBB CCCCCCCCCCCCCC CCC ψ BBBBBBBBBBB CCCCCCCCCCCCCC x BBBBBBBBBBB CCCCCCCCCCCCCC !!!! !!!!BBBBBBBBBBB CCCCCCCCCCCCCC M L !!!! BBBBBBBBBBB CCCCCCCCCCCCCC !!!! x’ L N !!!! !!!! BBBBBBBBBBB Figura 9.6.2. Angoli di Deprit e relazione con angoli di Eulero. Consideriamo la retta OM perpendicolare all’asse z e ad L e la retta ON perpendicolare a z ′ e L ed indichiamo con eM , eN i versori rispettivi. Gli angoli α, γ, β sono quelli formati dalle rette x, OM ; ON, x′ ; OM, ON rispettivamente. Indicando con eM = ez × L , ||ez × L|| eN = L × e′z , ||L × e′z || (9.6.12) i versori degli assi OM e ON si ha cos α = ex · eM , cos β = eM · eN cos γ = e′x · eN , (9.6.13) 220 c 88-08- 9820 9. Corpo rigido Usando l’identità vettoriale (a ×b)· (c ×d) = (a · c)(b · d)−(a · d)(b · c) possiamo esprimere gli angoli di Deprit attraverso gli angoli di Eulero ed i loro momenti coniugati. cos α = − q cos β = p cos γ = q Ly L2x + L2y = −p Ly L2 − L2z Lz L′z − L2 cos θ q L2 − L2z L2 − L′z 2 (9.6.14) L′y L2 − L′z 2 Si noti che Lz = pφ , L′z = pψ , e che L′y , Ly sono espresse attraverso gli angoli di Eulero ed i loro momenti coniugati da (9.6.4). Per avere le trasformazioni inverse osserviamo che l’asse x e le rette OM, OL sono su uno stesso piano e quindi φ − α e l’angolo tra OM e OL; anche le rette OL, ON e l’asse x′ sono su uno stesso piano e ψ − γ è l’angolo tra OL e ON , vedi figura 9.6.2. Essendo il versore dell’asse OL dato da eL = ez × e′z | sin θ|−1 si ha s r ′ 2 Lx L′z 2 Lz Lz cos θ = − cos β 1 − 2 1 − 2 L2 L L ′ Lz − Lz cos θ p cos(φ − α) = eL · eM = (9.6.15) | sin θ| L2 − L2z Lz − L′z cos θ q cos(ψ − γ) = eL · eN = | sin θ| L2 − L′z 2 dove la prima equazione si ottiene invertendo la seconda equazione in (9.6.14). Per esprimere H in funzione degli angoli α, β, γ e dei loro momenti coniugati, Lz , L, L′z , si noti che il piano formato da L e z ′ interseca il piano x′ y ′ lungo la retta OQ, che forma un angolo γ 2 con l’asse y ′ e π/2 − γ con l’asse x′ . La proiezione di L lungo OQ ha modulo (L2 − L′z )1/2 , per cui q q L2 − L′z 2 sin γ, L′x = L′y = − L2 − L′z 2 cos γ (9.6.16) e l’hamiltoniana si scrive 2 L′ L2 − L′z H= z + 2Jz 2 2 sin2 γ cos2 γ + Jx Jy (9.6.17) in accordo con (9.6.8). I momenti coniugati di α, β, γ sono Lz , L, L′z in quanto componenti del momento angolare lungo i rispettivi assi di rotazione, secondo il teorema di Nöther, enunciato nella versione hamiltoniana nel paragrafo 15.6. La rappresentazione delle rotazioni tra riferimenti, definita dagli angoli di Deprit o di Eulero, fornisce una carta della della varietà a 3 dimensioni che definisce le configurazioni del corpo rigido. Questa varietà va rappresentata su almeno due carte compatibili, date dagli angoli di Deprit (α, γ, β) e (α′ , γ ′ , β ′ ) riferiti a due diversi sistemi fissi.