Soluzione definiva (!) dell'orbita attorno a Papalterra (di Fabrizio e Umberto) ***

Questo articolo mostra la soluzione definitiva del problema dell'orbita minore possibile di un satellite artificiale, attorno alla Papalterra, una Terra con il centro di gravità spostato rispetto al centro geometrico. Tutte le informazioni e un paio di soluzioni si trovano QUI e QUI. In un successivo articolo (QUI) proponevo un metodo generale che portasse alla soluzione.

Fabrizio e Umberto lo hanno affrontato di petto e sono giunti alla conclusione. Nessuna delle soluzioni precedenti era la migliore... Devo dire che hanno fatto un lavoro encomiabile e decisamente laborioso dal punto di vista dei passaggi algebrici. Non hanno desistito e hanno trovato lo stesso risultato! Congratuliamoci vivamente con loro, una risorsa fondamentale di questo Circolo.

Il problema è stato impostato così:

Si è considerato un punto generico O (che stia sull'asse congiungente baricentro e centro di Papalterra) come centro dell'ellissi.

Si scrive , perciò l'equazione dell'ellisse generica che abbia centro in O e fuoco in G, e poi quella del cerchio di raggio R con centro in C. Si ponga come variabile fondamentale:

c = OC

La condizione fondamentale è che cerchio ed ellissi siano tangenti (ossia non si intersechino e si tocchino in solo due punti). In tal modo si ottiene l'equazione a = f(c) e si cerca il valore di c che dia il minimo semiasse (minimo semiasse vuol dire minimo periodo di rivoluzione). Un perfetto studio di funzione (da costruire) e di applicazione della derivata di un qualcosa non molto semplice.

La figura base è quella che segue:

La distanza OC, come già detto, è indicata con c, mentre d è la distanza (nota) tra il centro del pianeta e il centro di gravità.

Equazione circonferenza con centro in C:

$(x-c)^{2}+y^{2}=R^{2}$

equazione canonica dell'ellisse con centro negli assi:

ponendo c+d=e, si ha:

$b^{2}=a^{2}-e^{2}$ , (attenzione: questa e non ha niente a che vedere con l'eccentricità!)

$\frac{x^{2}}{a^{2}}+\frac{y^{2}}{a^{2}-e^{2}}=1$

queste sono tutte le possibili ellissi con fuoco in G.

Ricordiamo che un orbita ellittica è possibile se il centro di gravità G coincide con il fuoco dell'ellisse. E questo è stato adeguatamente tenuto in conto.

Limitiamoci al semipiano positivo delle y (possiamo farlo per questioni di simmetria); affinché la semi-ellisse contenga la semi-circonferenza, deve essere:

$(1-\frac{x^{2}}{a^{2}})(a^{2}-e^{2})\geq R^{2}-(x-c)^{2}$

$a^{2}-e^{2}-x^{2}+x^{2}\cdot \frac{e^{2}}{a^{2}}\geq R^{2}-x^{2}-c^{2}+2xc$

$x^{2}(\frac{e^{2}}{a^{2}})-x^{2}+a^{2}-e^{2}-R^{2}+x^{2}+c^{2}-2xc\geq 0$

$x^{2}(\frac{e^{2}}{a^{2}})+a^{2}-e^{2}-R^{2}+c^{2}-2xc\geq 0$

$x^{2}\cdot \frac{e^{2}}{a^{2}} -2xc + a^{2}-e^{2}-R^{2}-c^{2}\geq 0$

questa è una parabola rivolta verso l'alto (il coefficiente di $x^{2}$ è positivo) ; studiamo il discriminante:

$\Delta =c^{2}-\frac{e^{2}}{a^{2}}\cdot (a^{2}-e^{2}-R^{2}-c^{2})$ ; ponendo $\Delta =0$ (in tal modo essendo la parabola rivolta verso l'alto e l'intersezione con l'asse delle x unica, essa risulta sempre >=0 )e sostituendo a e=c+d; otteniamo una espressione per a (semiasse maggiore) in funzione di c;

$c^{2}-\frac{e^{2}}{a^{2}}\cdot (a^{2}-e^{2}-R^{2}-c^{2})=0$

moltiplichiamo per $a^{2}$ :

$a^{2}c^{2}- e^{2}\cdot (a^{2}-e^{2}-R^{2}+c^{2})=0$

$a^{2}c^{2}= e^{2}\cdot (a^{2}-e^{2}-R^{2}+c^{2})=0$

$a^{2}c^{2}-e^{2}a^{2}= e^{2}(-e^{2}-R^{2}+c^{2})=0$

$a^{2}(c^{2}-e^{2})= e^{2}(-e^{2}-R^{2}+c^{2})$

$a^{2}=\frac{ e^{2}(-e^{2}-R^{2}+c^{2})}{(c^{2}-e^{2})}$

$a^{2}=\frac{ (c+d)^{2}(-(c+d)^{2}-R^{2}+c^{2})}{(c^{2}-e^{2})}$

$a^{2}=\frac{ (c+d)^{2}(-c^{2}-d^{2}-2dc-R^{2}+c^{2})}{(c^{2}-e^{2})}=\frac{ (c+d)^{2}(-d^{2}-2dc-R^{2})}{(c^{2}-e^{2})}$

$c^{2}-e^{2}=c^{2}-(c+d)^{2}=c^{2}-c^{2}-2cd-d^{2}=-d(2c+d)$

$a^{2}=\frac{ (c+d)^{2}(-d^{2}-2dc-R^{2})}{-d(2c+d)}$

$a=(c+d)\cdot \sqrt{\frac{ (-d^{2}-2dc-R^{2})}{-d(2c+d)}}$

Calcolandola per un valore utilizzato in passato (c = 0), si ottiene il risultato voluto, che però non è il valore relativo al minimo valore di a, come si dimostrerà tra poco.

Forse, alcuni passaggi precedenti non sono chiarissimi a tutti (Umberto è decisamente un matematico "puro" e rigoroso!), e mi permetto di inserire i passaggi che ho svolto personalmente, che danno, ovviamente lo stesso risultato:

Consideriamo come origine il centro dell’ellissi generica.

La figura è quella precedente di Umberto

x²/a²+ y²/b² = 1

sappiamo anche però che:

OG = c + d

Da cui:

b² = a² – (c + d)²

Ponendo e = c + d

Si ha:

b² = a² – e²

da cui:

x²/a²+ y²/(a² – e²) = 1 …. (1)

La circonferenza di centro C (OC = c) ha equazione:

(x – c)² + y² = R² …. (2)

Facendo sistema tra (1) e (2) trovo le intersezioni tra ellissi e circonferenza

(x – c)² + y² = R²

x²/a²+ y²/(a² – e²) = 1

Ricavo y² dalla prima e lo sostituisco nella seconda:

y² = R² – (x-c)²

x²/a²+ (R² – (x-c)² )/(a² – e²) = 1

questa è l’equazione che ci deve dare le ascisse delle intersezioni, se esistono, tra ellissi e cerchio.

Svolgiamola per bene:

x²(a² – e²) + a²(R² – (x-c)²) = a²(a² – e²)

x²a² – x² e² + a²R² – a²x² - a²c² + 2a²xc - a⁴ + a²e² = 0

semplificando x²a²

si ottiene:

- e²x² + 2 a²xc + a²R² - a²c²- a⁴ + a²e² = 0

Cambiamo di segno:

e²x² - 2 a²xc - a²R² + a²c²+ a⁴ - a²e² = 0

A noi interessa solo il caso in cui le due curve siano tangenti, ossia che la soluzione sia unica. Questo si ottiene ponendo uguale a zero il discriminante di un’equazione di secondo grado ridotta:

a⁴c² – e²(- a²R² + a²c²+ a⁴ - a²e²) = 0

a⁴c² + e²a²R² – e²a²c² –a⁴e²+ a²e⁴ = 0

a²(a²c² + e²R² – e²c² – a²e² + e⁴) = 0

a non può essere uguale a zero, quindi:

a²c² + e²R² – e²c² – a²e² + e⁴ = 0

e è funzione di c e d, d è noto, R è noto. La funzione precedente vede a come funzione della sola c.

a²(c² – e²) = e²c² – e²R² –e⁴

La funzione va esplicitata in a:

a² = (e²c² – e²R² – e⁴)/( c² – e²)

a = e ((c² –R² – e²)/(c² – e²))^1/2

Prendiamo solo il segno più dato che a deve essere positiva

Poniamo, finalmente, (c + d) al posto di e

a = (c + d)((c² – R² – (c + d)²)/(c² – (c + d)²)^1/2

a = (c + d)((c² – R² – c² – d² – 2dc)/(c² – c² – d² – 2dc))^1/2

a = (c + d)((R² + d² + 2dc)/(d² + 2dc))^1/2

Questa sarebbe la funzione di cui si deve trovare il minimo…

Ovviamente, per

c = 0

vale a = (R² + d²)^1/2

Prendiamo la soluzione riportata nei commenti da Fabrizio:

a = (c + d) (1 + R²/(d(2c + d)))^1/2

e la sviluppiamo

a = (c + d) ((d(2c + d) + R2)/(d(2c +d))^1/2

a = (c + d) ((2dc + d² + R²)/(2dc + d²))^1/2

Che è perfettamente il risultato precedente…

La prima parte è completata: abbiamo la funzione di cui dobbiamo trovare il minimo. Essa è proprio del tipo:

a = f(c) (tutti gli altri parametri sono noti)

$a(c)=(c+d)\cdot \sqrt{\frac{ (d^{2}+2dc+R^{2})}{d(2c+d)}}$

$a(0)=\sqrt{d^{2}+R^{2}}$

Per facilitare i calcoli poniamo c = x e la faccenda diventa più... "normale".

Seguiamo, adesso, il procedimento di Umberto che porta allo stesso risultato di Fabrizio:

$a(x)=(x+d)\cdot \sqrt{\frac{ (d^{2}+2dx+R^{2})}{d(2x+d)}}$

Deve essere, ovviamente: $2x+d\neq 0$ ovvero $x\neq -d/2$

il numeratore è positivo per $d^{2}+2dx+R^{2}\geq 0$

ovvero $2dx\geq -R^{2}- d^{2}$

$x\geq \frac{-R^{2}- d^{2}}{2d}$ ,

il denominatore è positivo per $x>-d/2$ , il segno del rapporto sarà negativo quando i segni sono diversi, quindi in $\frac{-R^{2}- d^{2}}{2d}\leq x<- \frac{d}{2}$ la funzione non è definita

scriviamo la funzione come Fabrizio , dividendo per d(2x+d):

$a(x)=(x+d)\cdot \sqrt{1+\frac{ R^{2}}{d(2x+d)}}$

$Da(x)=\sqrt{1+\frac{R^{2}}{d(2x+d))}}-(x+d) \frac{R^{2}}{d(2x+d)^{2}\sqrt{1+\frac{R^{2}}{d(2x+d))}}}$

(derivata di un prodotto)

poniamo Da(x)=0

$\sqrt{1+\frac{R^{2}}{d(2x+d))}}-(x+d) \frac{R^{2}}{d(2x+d)^{2}\sqrt{1+\frac{R^{2}}{d(2x+d))}}}=0$

$d(2x+d)^{2}\sqrt{1+\frac{R^{2}}{d(2x+d))}})\sqrt{1+\frac{R^{2}}{d(2x+d))}}-(x+d) R^{2}=0$

${d(2x+d)^{2}+(2x+d)R^{2}}-(x+d) R^{2}=0$

${d(4x^{2}+d^{2}+4xd)+(2x+d)R^{2}}-(x+d) R^{2}=0$

$4dx^{2}+d^{3}+4xd^{2}+2xR^{2}+dR^{2}-xR^{2}-R^{2}d =0$

che porta all'equazione di 2° grado:

$4dx^{2}+(R^{2}+4d^{2})x+d^{3}=0$

Calcoliamo prima il discriminante:

$\Delta =R^{4}+16d^{4}+8R^{2}d^{2}-16d^{4}=R^{4}+8R^{2}d^{2}=R^{2}(R^{2}+8d^{2})$

poi le radici:

$x_{1,2}=\frac{-(R^{2}+4d^{2})\pm R\sqrt{R^{2}+8d^{2}}}{8d}$

Osserviamo che per x = c=0, la derivata non si annulla, quindi $a(0)=\sqrt{d^{2}+R^{2}}$ non è il minimo cercato.

La prima radice , è:

$\frac{-R^{2}-4d^{2}- R\sqrt{R^{2}+8d^{2}}}{8d}<-\frac{d}{2}$ proviamo a vedere

$-R^{2}-4d^{2}- R \cdot \sqrt{R^{2}+8d^{2}}<-4d^{2}$

$R^{2}+4d^{2}+ R \cdot \sqrt{R^{2}+8d^{2}}>4d^{2}$

che è vero, quindi la scartiamo

consideriamo la seconda soluzione e vediamo se

$\frac{-R^{2}-4d^{2}+ R\sqrt{R^{2}+8d^{2}}}{8d}>-\frac{d}{2}$

$-R^{2}-4d^{2}+ R \cdot \sqrt{R^{2}+8d^{2}}>-4d^{2}$

$R \cdot \sqrt{R^{2}+8d^{2}}>R^{2}$

che è vero .

Quindi, questa è la soluzione da prendere, ed è un minimo per la funzione a(x)

$a(x)=(x+d)\cdot \sqrt{1+\frac{ R^{2}}{d(2x+d)}}$

calcoliamo (x₂+d) e (2x₂+d)

$x_{2}+d=\frac{-R^{2}-4d^{2}+ R\sqrt{R^{2}+8d^{2}}}{8d} +d=\frac{-R^{2}+4d^{2}+ R\sqrt{R^{2}+8d^{2}}}{8d}$

$2x_{2}+d=2\frac{-R^{2}-4d^{2}+ R\sqrt{R^{2}+8d^{2}}}{8d} +d=\frac{-R^{2}+ R\sqrt{R^{2}+8d^{2}}}{4d}$ $a(x_{2})=\frac{-R^{2}+4d^{2}+ R\sqrt{R^{2}+8d^{2}}}{8d} \sqrt{1+\frac{R^{2}}{2d\frac{-R^{2}+ R\sqrt{R^{2}+8d^{2}}}{4d} }}$

$a(x_{2})=\frac{-R^{2}+4d^{2}+ R\sqrt{R^{2}+8d^{2}}}{8d} \sqrt{1+\frac{R^{2}}{\frac{-R^{2}+ R\sqrt{R^{2}+8d^{2}}}{2} }}$

$a(x_{2})=\frac{-R^{2}+4d^{2}+ R\sqrt{R^{2}+8d^{2}}}{8d} \sqrt{1+\frac{2R^{2}}{\frac{-R^{2}+ R\sqrt{R^{2}+8d^{2}}}{1} }}$

$a(x_{2})=\frac{-R^{2}+4d^{2}+ R\sqrt{R^{2}+8d^{2}}}{8d} \sqrt{\frac{R^{2}+R\sqrt[]{R^{2}+8d^{2}}}{-R^{2}+\sqrt{R^{2}+8d^{2}}}}$

Chi ha voglia può divertirsi a disegnarla e a trovare la posizione del perigeo, ma pensiamo si sia già lavorato abbastanza...

L	M	M	G	V	S	D
	1	2	3	4	5	6
7	8	9	10	11	12	13
14	15	16	17	18	19	20
21	22	23	24	25	26	27
28	29	30

Soluzione definiva (!) dell'orbita attorno a Papalterra (di Fabrizio e Umberto) ***

Soluzione definiva (!) dell'orbita attorno a Papalterra (di Fabrizio e Umberto) ***

1 commento

Lascia un commento Cancel Reply

Soluzione definiva (!) dell'orbita attorno a Papalterra (di Fabrizio e Umberto) ***

Articoli correlati

Perché il problema dei tre corpi è caotico ? **

Soluzione (una delle tante possibili) al problema del numero primo ****

Facciamo sul serio ... ****

La spirale aurea è un nostra creazione... accontentiamoci dei numeri irrazionali **

1 commento

Lascia un commento Cancel Reply