Problema dels dos cossos

determinar el moviment de dos partícules puntuals que interactuen només entre elles

En mecànica clàssica, el problema dels dos cossos té per objectiu determinar el moviment de dues partícules puntuals que només interactuen entre si.

Dos cossos orbitant al voltant del seu centre de masses en òrbites el·líptiques.
Dos cossos amb una petita diferència de massa orbitant al voltant del seu centre de massa, les mides dibuixats són similars als del sistema Plutó-Caront.

Les lleis de Newton permeten reduir el problema dels dos cossos a un problema d'un cos equivalent, és a dir, resoldre el moviment d'una partícula sotmesa a un camp gravitatori conservatiu i que per tant deriva d'un potencial extern. Atès que el problema es pot resoldre de forma exacta, el problema dels dos cossos corresponent també es pot resoldre amb exactitud. Per contra, el problema dels tres cossos (i, més generalment, el problema dels n cossos amb ) no pot es pot resoldre analíticament de forma exacta, excepte en casos especials.

Descripció del problema

modifica
 
Coordenades de Jacobi en el problema dels dos cossos; les coordenades són   i   on  .[1]

Siguin   i   les posicions de dos cossos, i   i   les seves masses.

La segona llei de Newton determina que

  (equació 1)
  (equació 2)

on   és la força que experimenta la massa 1 a causa de la seva interacció amb la massa 2,   és la força que experimenta la massa 2 a causa de la seva interacció amb la massa 1, i   representen la derivada segona respecte el temps de la posició dels cossos 1 i 2 respectivament.

La suma de les dues equacions produeix una equació que descriu el moviment del centre de masses (baricentre) del sistema.[2] La resta de la segona equació de la primera produeix una equació que descriu la variació amb el temps del vector r = x1 − x₂ entre les dues masses.[2] Les solucions d'aquestes dues equacions independents de problemes d'un cos poden combinar-se per a obtenir les trajectòries x1(t) i x(t).

Moviment de les dues masses

modifica

Moviment del centre de masses (primer problema d'un cos)

modifica

La suma de les dues equacions (1) i (2) produeix:[2]

 

on s'ha emprat la tercera llei de Newton, F₁₂ = −F21, i on

 

  és la posició del centre de masses del sistema.

L'equació resultant

 

mostra com la velocitat V = dR/dt del centre de masses és constant, d'on s'extreu que la quantitat de moviment m1 v1 + mv₂ també es constant (conservació del moment). Per tant, la posició R(t) del centre de masses pot ser determinada en qualsevol instant de temps a partir de les posicions i velocitats inicials.

Moviment del vector de desplaçament (segon problema d'un cos)

modifica

Restant l'equació (2) de l'equació (1), s'obté:[2]

 

on s'ha emprat   segons la tercera llei de Newton.

El vector de posició relativa entre les dues masses,  , és

 

L'equació pot reescriure's com:

 

on   és la massa reduïda i s'expressa com

 

Moviment de les dues masses

modifica

Les equacions d'  i   permeten obtenir finalment les equacions de moviment de cada un dels dos cossos del sistema:[2]

 
 

Ambdues equacions es verifiquen de forma senzilla mitjançant la substitució en elles d'  i  .

Propietats del moviment

modifica

El moviment de dos cossos és pla

modifica

El moviment de dos cossos que interactuen entre si sempre està en un pla.[3]

El moment angular del sistema,   s'expressa com

 

on   és el moment lineal i   és la massa reduïda.   denota un producte vectorial.

La variació del moment angular amb el temps és igual al moment de força  

 

on

 

La força entre les dues partícules està en la mateixa línia que les uneix i és paral·lela al radi vector,  . Per tant, utilitzant la propietat del producte vectorial que estableix que el producte vectorial de dos vectors que apunten a la mateixa direcció és nul,  ,

 

El moment angular   és constant (es conserva). Per tant, el vector de desplaçament   i el vector de velocitat   sempre es troben en un pla perpendicular al vector constant  .

Llei de les àrees

modifica

Per demostrar la llei de les àrees o equació de Kepler de forma senzilla, és útil canviar a les coordenades polars, ja que el moviment està en un pla i, per a molts problemes físics, la força   només és una funció del radi   (és una força central).

En moure's durant un instant de temps el vector de posició   descriu una àrea elemental

 

de manera que la velocitat areolar o àrea escombrada pel vector de posició en la unitat de temps és:

 .

El mòdul del moment angular

 

on  .

Per tant, es pot expressar la velocitat areolar en funció del moment angular

 

amb   és la "constant de les àrees".

Aquesta llei de les àrees o segona llei de Kepler va ser enunciada empíricament per primera vegada el 1609 per Johannes Kepler i explica el moviment dels planetes al voltant del Sol. Aquest fenomen és una propietat general del moviment de les forces centrals. Per tant, és per tant més general que les forces de la gravitació inversament proporcionals al quadrat de la distància.

El moviment d'un planeta en el pla de la seva òrbita es compon de dos moviments: un l'angle que gira el radi vector i l'altre el seu acostament o allunyament del primari, és a dir la variació del mòdul del radi vector amb el temps. La llei de les àrees determina que, un cos gira més de pressa quan és a prop i lent quan està lluny i ho fa quantitativament, com per poder establir l'angle de gir, encara que és difícil. Per obtenir l'angle de gir E amb el temps cal expressar aquesta fórmula d'una altra manera:

 

Aquesta fórmula s'anomena equació de Kepler, on M és l'anomalia mitjana, e és l'excentricitat i E l'anomalia excèntrica.

L'òrbita

modifica

Tot objecte en l'univers atrau a un altre objecte al llarg de la línia que uneix el centre dels objectes, proporcional a les masses de cada objecte, i inversament proporcional al quadrat de la distància entre ells.[4]

En la segona llei de Newton l'acceleració a és de la forma

 

En coordenades polars, i considerant que l'òrbita està en el pla OXY, la velocitat i l'acceleració s'expressen

 
 

Descomposant l'acceleració per components, i donat que només existeix força en la component radial:

 
 

Substituint   i   en la segona equació, s'obté

 

Reordenant i agrupant les variables, s'obté

 

La integració d'aquesta equació és:

 

on   és la constant d'integració, essent   el moment angular específic (per unitat de massa).[1]

A continuació, s'introdueix el següent canvi de variable:

 

I s'obtenen les noves expressions per la velocitat i l'acceleració, respectivament:

 
 

En l'expressió de l'acceleració, s'ha tingut en compte que  . L'equació de moviment en  ,

 

es reescriu com

 

La llei de la gravitació universal indica que la força per unitat de massa és

 

on G és la constant de gravitació universal i M és la massa de l'estrella. Per tant,

 

Aquesta equació diferencial té la solució general

 

on e i θ0 són constants arbitràries d'integració. Desfent el canvi de variable, i considerant  ,

 

Aquesta és l'equació d'una cònica amb excentricitat e i origen en un focus. Per tant, la primera llei de Kepler és un resultat directe de la llei de la gravitació i de la segona llei de Newton del moviment.   rep el nom d'anomalia veritable i és l'angle que forma el radi vector amb el periastre. L'excentricitat i l'anomalia veritable són elements orbitals que defineixen una òrbita.


Referències

modifica
  1. 1,0 1,1 Betounes, David. Springer. Differential Equations, 2001, p. 58; Figura 2.15. ISBN 0-387-95140-7. 
  2. 2,0 2,1 2,2 2,3 2,4 Robert Taylor, John. «Two-Body Central Force Problems». A: University Science Books. Classical Mechanics (en anglès), 2005. ISBN 189138922X. 
  3. Montenbruck, Oliver; Gill, Eberhard. «Introductory Astrodynamics». A: Springer Science & Business Media. Satellite Orbits: Models, Methods, and Applications (en anglès), 2000. 
  4. Newton, Isaac. Philosophiæ Naturalis Principia Mathematica (en llatí), 1686/1687.  Traducció d'Andrew Motte (anglès)

Vegeu també

modifica

Bibliografia

modifica
  • Landau, Lev; Lifshitz, Evgeny. Pergamon Press. Mechanics (en anglès). 3a edició, 1976. ISBN 0-08-029141-4. 
  • Goldstein H, Herbert. Addison-Wesley. Classical Mechanics (en anglès). 2a edició, 1980. ISBN 0-201-02918-9.