Orbite dei corpi celesti e tipi di traiettorie

Come stabilito dalla prima legge di Keplero le traiettorie dei pianeti del Sistema Solare sono ellittiche. Tale proprietà vale per qualsiasi corpo celeste? Più in generale, come si determina la traiettoria di un corpo celeste qualsiasi e quali grandezze dobbiamo analizzare per individuare il tipo di traiettoria?

In questa lunga e densissima lezione, dedicata ai soli studenti universitari, vedremo che le traiettorie dei corpi celesti sono coniche e che analizzandone l'energia è possibile stabilire se esse siano ellittiche, paraboliche o iperboliche.

Come determinare la traiettoria di un corpo celeste

Abbiamo già verificato che i principi della Dinamica di Newton ci permettono di dare una giustificazione teorica per le leggi di Keplero, ma lo abbiamo sempre fatto nell'ipotesi di una traiettoria circolare perché tale presupposto semplifica di molto la trattazione.

Ora vogliamo dimostrare che sotto l'effetto della forza di gravitazione universale, che è una forza centrale dipendente dall'inverso del quadrato della distanza, l'orbita di un pianeta è un'ellisse; più in generale, che l'orbita di un corpo celeste attorno a un altro è una conica.

Premettiamo brevemente che una conica (ellisse, iperbole, parabola) è descritta in matematica come il luogo geometrico dei punti (x,y) tali per cui il rapporto tra la distanza da un punto, detto fuoco, e la distanza da una retta, detta direttrice, è pari a una costante positiva ε chiamata eccentricità.

In particolare se si sceglie di usare un sistema di riferimento in coordinate polari, allora le coordinate sono (r,θ).

Traiettoria corpi celesti

Con r indichiamo la distanza del punto rispetto al fuoco e con θ l'angolo che r forma con la direzione perpendicolare alla direttrice. L'equazione di una conica in coordinate polari è della forma

(1)/(r) = (1)/(ε d)-(1)/(d) cos(θ)

Si distinguono tre casi:

- se ε < 1 si ottiene un'ellisse;

- se ε = 1 si ha una parabola;

- se ε > 1 si ottiene un'iperbole.

Andrebbe considerato anche il caso limite ε = 0 in cui si ottiene una circonferenza, ma essa non è altro che un caso particolare dell'ellisse (asse maggiore e asse minore di pari lunghezza).

Premesso ciò, consideriamo due corpi che interagiscono gravitazionalmente; entrambi saranno soggetti a una forza, che dal secondo principio della Dinamica potremo esprimere nella forma

F = m a_(m) ; F = -M a_(M)

Il segno meno della seconda forza si deve al principio di azione-reazione: i due corpi sono soggetti a forze uguali e contrarie. L'accelerazione cui è soggetto il corpo di massa m rispetto a quello di massa M è:

a = a_(m)-a_(M) = F ((1)/(m)+(1)/(M)) = F ((M+m)/(mM))

Di conseguenza, la forza agente su entrambi i corpi è data da:

F = ((mM)/(M+m)) a = μ a

dove la grandezza μ

μ: = (mM)/(M+m)

viene detta massa ridotta del sistema. Quindi il moto relativo di due punti che interagiscono tra di loro con una certa forza F è uguale a quello di un punto di massa pari alla massa ridotta e sottoposto alla medesima forza F.

Visto che stiamo lavorando con la forza gravitazionale, possiamo scrivere:

F = -G (mM)/(r^(2)) u_(r)

dove u_r è il versore lungo la direzione di r. Per una forza centrale è possibile scrivere l'accelerazione in coordinate polari secondo la cosiddetta formula di Binet (che non dimostriamo):

a = -(L^(2))/(m^(2)r^(2)) [ (dr^(2))/(d θ^(2)) ((1)/(r))+(1)/(r) ] u_(r)

dove L indica come sempre il momento angolare.

Tornando allora alla formula della forza in funzione della massa ridotta, vale a dire F = μa, abbiamo:

-G (mM)/(r^(2)) u_(r) = -μ (L^(2))/(μ^(2)r^(2)) [ (dr^(2))/(d θ^(2)) ((1)/(r))+(1)/(r) ] u_(r)

Riordinando l'equazione e considerando solo i moduli dei vettori, si ottiene la seguente equazione differenziale:

GmM = (L^(2))/(μ) [ (dr^(2))/(d θ^(2)) ((1)/(r))+(1)/(r) ]

o, in forma canonica

(dr^(2))/(dθ^(2)) ((1)/(r))+(1)/(r) = (G μ m M)/(L^(2))

La soluzione ci dà la posizione r = r(θ) in funzione dell'angolo θ, ossia la traiettoria descritta dal corpo m rispetto a M. Ecco come si presenta la soluzione:

(1)/(r) = (G μ mM)/(L^(2))+A cos(θ)

dove A è una costante che può essere determinata mediante le condizioni iniziali.

Come potete vedere, la funzione che descrive la traiettoria ha la stessa identica struttura dell'equazione generale delle coniche che abbiamo scritto all'inizio. La forza di gravitazione ci ha portato a dire che il corpo di massa m rispetto a quello di massa M descrive una traiettoria data da una conica. 

Tipo di orbita, momento angolare ed energia

Fin qui abbiamo compreso che le traiettorie dei pianeti e dei satelliti sono coniche, ma ciò non significa che siano necessariamente ellittiche: potrebbe anche essere paraboliche o iperboliche.

Confrontiamo l'equazione generale di una conica in coordinate polari con l'equazione ottenuta in precedenza:

 (1)/(r) = (1)/(ε d)-(1)/(d) cos(θ) ; (1)/(r) = (G μ mM)/(L^(2))+A cos(θ)

Dal confronto emerge che:

L^(2) = G μ m M ε d

Il momento angolare del corpo orbitante dipende dai parametri ε,d dell'orbita; una volta fissati questi parametri, il momento angolare diventa una costante, come già sapevamo.

Ora proviamo invece a ricavare l'energia meccanica in funzione dei parametri dell'orbita, dove nell'energia cinetica usiamo la massa ridotta μ del sistema.

E = (1)/(2) μ v^(2)-G (mM)/(r)

Visto che la forza gravitazionale è una forza centrale, il moto avviene sicuramente su un piano (come abbiamo visto nella lezione sulle forze centrali). Ciò ci permette di riscrivere la velocità in funzione delle coordinate polari.

v = √(((dr)/(dt))^(2)+r^(2) ((d θ)/(dt))^(2))

Da cui possiamo ricavarne il quadrato

v^(2) = ((dr)/(dt))^(2)+r^(2) ((d θ)/(dt))^(2)

Dunque, tornando all'energia, abbiamo

E = (1)/(2) μ [((dr)/(dt))^(2)+r^(2) ((d θ)/(dt))^(2)]-G (mM)/(r)

ossia

E = (1)/(2) μ ((dr)/(dt))^(2)+(1)/(2) μ r^(2) ((d θ)/(dt))^(2)-G (mM)/(r)

Ora dobbiamo trasformare questa formula facendo alcuni passaggi. Innanzitutto ricordiamoci che il momento angolare può essere riscritto nella seguente forma

L = μ rv = μ r^(2) ω = μ r^(2) (d θ)/(dt) (•)

Poi possiamo ricavare il termine (dr)/(dt) differenziando l'equazione generale delle coniche in coordinate polari:

(d)/(dt) ((1)/(r)) = -(1)/(r^(2)) (dr)/(dt) = (sin(θ))/(d) (d θ)/(dt)

da cui

(dr)/(dt) = -(r^(2) sin(θ))/(d) (d θ)/(dt)

Eleviamo l'ultima equazione al quadrato:

((dr)/(dt))^(2) = (r^(4) sin^(2)(θ))/(d^(2)) ((d θ)/(dt))^(2)

Dalla formula del momento angolare (•) ricaviamo:

 L = μ r^(2) (d θ)/(dt) ; L^(2) = μ^(2) r^(4) ((d θ)/(dt))^(2) ; r^(4) ((d θ)/(dt))^(2) = (L^(2))/(μ^(2))

Di conseguenza:

((dr)/(dt))^(2) = (sin^(2) (θ) : L^(2))/(d^(2) μ^(2))

A questo punto siamo pronti per poter riscrivere il primo termine dell'equazione dell'energia:

(1)/(2) μ ((dr)/(dt))^(2) = (1)/(2) μ (sin^(2) (θ) : L^(2))/(d^(2) μ^(2)) = (1)/(2) (sin^(2)(θ) L^(2))/(d^(2) μ)

Trasformiamo anche il secondo termine dell'equazione dell'energia, facendo comparire il momento angolare:

(1)/(2) μ r^(2) ((d θ)/(dt))^(2) = (1)/(2) (L^(2))/(μ r^(2))

Non ci resta che trasformare il terzo temine dell'energia e per farlo rimaneggiamo l'equazione che abbiamo trovato all'inizio per il quadrato del momento angolare.

L^(2) = G μ m M ε d → GmM = (L^(2))/(μ ε d)

In questo modo il terzo termine dell'equazione dell'energia diventa:

-G (mM)/(r) = -(L^(2))/(μ ε d r)

Ora possiamo riscrivere l'intera equazione dell'energia:

E = (1)/(2) (sin^(2) (θ) L^(2))/(d^(2) μ)+(1)/(2) (L^(2))/(μ r^(2))-(L^(2))/(μ ε d r)

Trasformando i termini (1)/(r), (1)/(r^2) mediante l'equazione generale delle coniche, e utilizzando l'equazione del quadrato del momento angolare, è possibile riscrivere l'energia in funzione dei parametri dell'orbita ε,d e in una forma più compatta.

E = G (mM)/(2d ε) (ε^(2)-1)

È stato un conto lungo e impegnativo, ma siamo arrivati a un risultato che ci permette di capire una cosa importante. Quando l'energia del sistema è negativa, allora l'eccentricità deve essere minore di 1 (ε < 1) e ciò significa che il corpo di massa m descrive un orbita ellittica. Quando l'energia è nulla, allora ε = 1 e l'orbita è parabolica e infine, quando l'energia è positiva, allora ε > 1 e l'orbita è iperbolica.

È dunque l'energia del sistema a determinare il tipo di traiettoria seguita, ad esempio da un corpo celeste attorno al Sole. Il sistema Terra-Sole è legato è questo significa che l'energia potenziale è in valore assoluto maggiore di quella cinetica; ciò fa sì che l'energia totale sia negativa, pertanto l'orbita della Terra intorno al Sole è ellittica. Il fatto che l'energia sia negativa significa che la Terra non ha abbastanza energia cinetica (e quindi velocità) per poter sfuggire al campo gravitazionale solare (per nostra fortuna!).

Per un corpo in orbita parabolica invece l'energia cinetica e quella potenziale si equivalgono perfettamente, cosicché l'energia totale è nulla. Se consideriamo una cometa, per esempio, essa curverebbe la propria traiettoria in prossimità del Sole per poi allontanarsi descrivendo una parabola.

Un'orbita iperbolica è invece descritta da corpi la cui energia cinetica è maggiore in valore assoluto di quella potenziale; la cometa sarebbe così lanciata a una velocità tale per cui il campo gravitazionale solare non sarebbe in grado di trattenerla e più la velocità è alta, più l'orbita è un'iperbole prossima a una retta, cioè un'iperbole con un alto valore di eccentricità.


Siamo in dirittura d'arrivo! Nella lezione successiva - la penultima del capitolo dedicato alla Gravitazione Universale, tratteremo il teorema di Gauss per il campo gravitazionale. Nel caso voleste curiosare un po' vi rimandiamo al tool di Astronomia online dedicato alle orbite. ;)

Buona Fisica a tutti!

Alessandro Catania (Alex)

Lezione precedente.....Lezione successiva


Tags: energia e traiettorie di pianeti e satelliti - come determinare il tipo di orbita di un pianeta o di un satellite.

Ultima modifica: