4) Demostración primera ley de Kepler
Por medio de la primera ley de Newton y la ley de gravitación universal, podemos
establecer la primera ley de Kepler. La primera ley de Kepler describe las órbitas planetarias como
elípticas, con el Sol ubicado en uno de los focos de la elipse.
Para demostrar la primera ley de Kepler, también conocida como: ley de las órbitas elípticas, se
utiliza la dinámica de movimiento bajo la influencia de la fuerza gravitatoria descrita por la ley de
gravitación universal de Newton. Esta ley implica una fuerza central que actúa sobre el planeta y
dirige su movimiento hacia el Sol.
Segunda ley de Newton F=m∗a G = Constante de gravitación universal
GMm
Ley de la Gravitación: F= 2 M= Masa del sol
r
m = Masa de un planeta R = distancia entre el sol y el planeta
Para nuestra demostración tenemos que tomar en cuenta que al aplicar la ley de gravitación
universal para calcular trayectorias orbitales de planetas alrededor del Sol, consideramos la
dirección de la fuerza gravitatoria como una atracción hacia el Sol, conforme a la convención de
que la fuerza gravitatoria es siempre atractiva entre masas positivas por tanto la consideramos
negativa.
−GMm
F= 2
r
La relación entre la fuerza gravitatoria y la posición vectorial de nuestro planeta puede describirse
mediante la ley de gravitación universal, multiplicada por el vector posición r de nuestro planeta.
Antes de esto, es necesario transformar r en un vector unitario, denotado como u.
rr r
u= → →u=
√r r 2 r
Entonces al relacionar un vector unitario por medio de la ley de gravitación obtenemos:
−GMm
2
∗r
−GMm r −GMm
F= 2
U → F= → F= 3
r
r r r
Para avanzar en los cálculos, es importante obtener una expresión que relacione las fuerzas
gravitatorias con las leyes del movimiento, permitiéndonos comprender cómo varían la velocidad y
la posición del planeta. Esto implica empezar a relacionar la segunda ley de Newton con la ley de
gravitación universal.
GMm GMm GM
ma= 2
→ ma= 2 → a= 2
r r r
d d
También es importante enmarcar que: (r )=v , 2
( r )=a , podemos notar por medio de la
dt dt
siguiente grafica que a es paralelo a r:
Es importante mencionar que un planeta se mueve en un plano, para hallar dicho plano
necesitamos un vector perpendicular a r, entonces podemos notar que para obtener un vector
perpendicular a r se debe realizar el producto cruz entre dos vectores, para este caso entre r x v,
en consecuencia, a que a es paralela a r esto debido a:
d d
(r x v )→ ¿
dt dt
d
( r x v )=r ' x v +r x v ' → r ' x v +r x v '
dt
r x v=h
De esta se infiere a que el vector r es perpendicular a h para todos los valores t así que el planeta
que esta orbitando alrededor del sol siempre se ubica en el plano que pasa por el origen
perpendicular a h.
El vector h es el momento cinético especifico el cual se refiere a una magnitud física que es
utilizada en la mecánica celeste para caracterizar el movimiento de un cuerpo.
Podemos reescribir h en base a datos anteriores como:
h=r x v → r x r ' → ru x ( ru ) → ru x ( r u' + r ´ u )
'
→ r 2 ( u u' ) +r r ' ( u x u )=r 2 (u x u' )
Entonces ahora lo que necesitaríamos es un vector que nos represente la semidistancia a partir del
foco, gráficamente este vector lo podemos obtener por medio del producto cruz entre a x h
−GM
u x ( r u x u ) ( por regla del producto cruz se cancelan r )
2 ' 2
a x h= 2
r
¿−GM u x ( u x u ' )
−GM [ ( u ∙ u ) u−( u ∙u ) u ] Por regla del producto cruz se puede llegar a la simplificación
' '
'
u ∙u =0
De lo anterior se obtiene que:
'
a × h=GM u
O bien si
'
a=v =¿ Aceleracion es igual a la derivada de velocidad
Entonces se puede decir que (v × h)' =a ×h . Y como h es una constante entonces:
( v ' × h )=a× h
De este modo también obtenemos que (v × h)' =(v¿ ¿' ×h)=a× h=GM u ' ¿
Ahora al integrar esa ecuación obtenemos:
'
a × h=GMu
' '
v ×h=GM u
v × h=GMu+ c
Donde c es un vector constante.
En este punto es conveniente elegir los ejes de coordenadas de tal manera que el vector
de base estándar k apunte en la dirección del vector h. Entonces, el planeta se mueve
en el plano xy. Puesto que tanto v X h como u son perpendiculares a h, la ecuación
v × h=GMu+ c
indica que c se tiende en el plano xy. Esto significa que se puede elegir los ejes “x” y “y” de
tal forma que el vector i se tienda en la dirección de c, como se muestra en la figura:
Si θ es el ángulo entre c y r, entonces (r, θ ) son las coordenadas polares del planeta.
De la ecuación se obtiene que
r ∙ ( v ×h ) =r ∙ ( GMu+c ) =GMr ∙ u+r ∙ c
Para entender la ecuación v × h=GMu+ c, y la expresión r ( v × h )=r ∙GMu+ r ∙ c , comencemos
analizando cada parte paso a paso.
Dado que:
1. a es la aceleración.
2. v es la velocidad.
3. a=v′, donde el apóstrofe denota la derivada temporal.
4. h=r×v es el momento angular específico.
5. GM es el producto de la constante gravitacional y la masa del cuerpo central.
6. c es un vector constante.
Primero, verifiquemos la ecuación v × h=GMu+ c :
Para hacerlo, podemos usar la definición de h y la relación entre aceleración, velocidad y posición
en el contexto de una órbita kepleriana.
Paso 1: Producto Cruz de velocidad y momento Angular
v × h=v ×(r × v )
Usamos la identidad del producto cruz:
v × ( r × v )=r ( v ∙ v )−v ( r ∙ v )
2
Dado que v ∙ v=v y r ∙ v es la proyección de la velocidad en la dirección del radio vector, notamos
que:
r ∙ v=r v r
Así que:
2
v × ( r × v )=r v −v (r v r )
Paso 2: Ecuación de la Velocidad Cruzada con el Momento Angular
La ecuación dad v × h=GMu+ c sugiere que producto se puede descomponer en dos
componentes: una componente que depende de la gravitación GMU y una constante adicional C.
Si consideramos una órbita kepleriana, podemos usar la relación de la aceleración centrípeta para
deducir esta forma. Recordemos que:
−GM
a= 2
u
r
Como a=v '
' −GM
v= 2
u
r
Podemos relacionar esta aceleración con la conservación del momento angular específico en la
órbita.
Paso 3: Comprobar la Relación
r ∙( v × h)=r ∙(GMu+c)
Usamos la propiedad distributiva del producto punto:
r ∙ ( v ×h ) =r ∙ ( GMu ) +r ∙ c
Paso 4: Simplificación
r
La expresión r ∙ (GMu ) se simplifica porque u =
r
r
r ∙ (GMu )=GMr =GMr
r
Por lo tanto:
r ∙ ( v ×h ) =GMr+ r ∙ c
Esto verifica que:
r ∙ ( v ×h ) =r ∙ GM u+r ∙ c
r ∙ ( v ×h ) =GMru ∙u+|r||c|cosθ=GMr +rc cosθ
Así, la deducción confirma la relación inicial, mostrando cómo la combinación de las propiedades
del momento angular, la aceleración gravitacional y las identidades vectoriales llevan a la expresión
deseada.
Donde c = |c|, entonces
2
eh
r ∙(v ×h) c
r= =
GM + c cosθ 1+ e cosθ
Donde h = |h|, así:
2
eh e h2
(GM ) c
r= =
1+e cosθ 1+e cosθ
2
h
Al escribir d= se obtiene la ecuación
c
ed
r=
1+e cosθ
Primera ley de Kepler
