Física de Campos.

Notas de Clase Interacción Gravitacional y Campo Gravitacional 2015-02 Presencial.

ANALISIS DE TRAYECTORIAS.

Partamos de la ley gravitacional de Newton:

−𝒢𝑚𝑚′
⃗ 𝑔 (𝒓
𝑭 ⃗)= 𝒊̂ ′
|𝒓 ⃗ ′ |2 (𝒓⃗−𝒓⃗ )
⃗ −𝒓

Ahora si la Fuente de Campo Gravitacional está ubicada en el centro del marco de

referencia o sistema coordenado tendríamos:

⃗ ′ = 0𝒊̂𝑥 + 0𝒊̂𝑦 + 0𝒊̂𝑧 = 𝓸

𝒓 ⃗⃗

−𝒢𝑚𝑚′ −𝒢𝑚𝑚′
⃗ 𝑔 (𝒓
𝑭 ⃗)= ⃗ 𝑔 (𝒓
𝒊̂ , 𝑭 ⃗)= 𝒓̂
⃗ |2 𝒓⃗
|𝒓 𝑟2

Aplicando la Segunda Ley de Newton:

−𝐷
⃗ 𝑁𝑒𝑡𝑎 (𝒓
𝑭 ⃗ ) = 𝑚𝒂 ⃗ 𝑔 (𝒓
⃗, 𝑭 ⃗)= 𝒊̂𝒓⃗ , 𝐷 = 𝒢𝑚𝑚′
|𝒓
⃗| 2

De la cual se deriva la Ecuación Diferencial de Movimiento:

−𝐷 ⃗
𝑑𝒗 −𝐷 𝑑2𝒓⃗ −𝐷
⃗ =
𝒂 𝒊̂ ⟹ = 𝒊̂ ⟹ = 𝒊̂
⃗ |2 𝒓⃗
𝑚|𝒓 𝑑𝑡 𝑚|𝒓⃗ |2 𝒓⃗ 𝑑𝑡 2 𝑚|𝒓⃗ |2 𝒓⃗

𝑑2𝒓⃗ −𝒢𝑚′
⃗⃗ (𝒓
=𝒈 ⃗ ); 𝒈
⃗⃗ (𝒓
⃗)= 𝒊̂
𝑑𝑡 2 ⃗ |2 𝒓⃗
|𝒓

Cuya solución es la ecuación de la elipse sino que se generaliza a la de cualquier

sección cónica:

2
⃗ | ⁄𝐷𝑚)
(|𝑳
𝑟=
1 + 𝑒 cos(𝜃)
La fuerza gravitacional es una fuerza central, por tanto es paralela al vector posición
pero de sentido contrario. Ahora bien aplicando segunda ley de Newton que tendría el
mismo comportamiento, podremos afirmar que el torque neto es nulo dado que la
fuerza y el vector posición son paralelos, tenemos:

⃗
𝑑𝒗 𝑑(𝑚𝒗⃗) ⃗
𝑑𝒑
⃗𝑴
⃗⃗ 𝑁𝑒𝑡𝑜 = ⃗𝒓 × ⃗𝑭𝑁𝑒𝑡𝑎 = ⃗𝒓 × 𝑚𝒂
⃗ = ⃗𝒓 × 𝑚 = ⃗𝒓 × = ⃗𝒓 ×
𝑑𝑡 𝑑𝑡 𝑑𝑡
Profesor: Julián Fernando Ruiz Roa.
Septiembre 04 2015
1|7
Universidad de Antioquia.
Instituto de Física. Facultad de Ingenierías.
⃗
𝑑𝑳 ⃗ ×𝒑
𝑑(𝒓 ⃗) ⃗
𝑑𝒑 ⃗
𝑑𝒓 ⃗
𝑑𝒑
⃗𝑴
⃗⃗ 𝑁𝑒𝑡𝑜 = = ⃗ ×
=𝒓 + ⃗ =𝒓
×𝒑 ⃗ × ⃗ × (𝑚𝒗
+𝒗 ⃗)
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡
⃗
𝑑𝑳 ⃗
𝑑𝒑 ⃗
𝑑𝒑
⃗⃗⃗ 𝑁𝑒𝑡𝑜 =
𝑴 ⃗ ×
=𝒓 + 𝑚(𝒗 ⃗)=𝒓
⃗ ×𝒗 ⃗ ×
𝑑𝑡 𝑑𝑡 𝑑𝑡
Sabemos que el torque es nulo, decimos que el momento angular se conserva:

⃗
𝑑𝑳 ⃗
𝑑𝒑
⃗𝑴
⃗⃗ 𝑁𝑒𝑡𝑜 = 𝓸
⃗⃗ ⟹ ⃗ ×
=𝒓 ⃗⃗ ⟹ ⃗𝑳: 𝑐𝑡𝑒
=𝓸
𝑑𝑡 𝑑𝑡
Lo cual implica que la trayectoria hace que el móvil se mantenga constante en un
plano. Y explica porque las áreas barridas en tiempos iguales sean las mismas.

⃗⃗
𝑑𝑨 1 ⃗
𝑑𝒗 1 ⃗𝑳
⃗ × )= ( )
= (𝒓
𝑑𝑡 2 𝑑𝑡 2 𝑚

La ecuación general de las secciones cónicas es:

𝑒𝑑
𝑟=
1 + 𝑒 cos(𝜃)

Ahora geométricamente y físicamente tenemos:

⃗⃗
𝑑𝑨 1 𝑑𝜃
= 𝑟 2 𝜔𝒊̂𝑧 ; 𝜔 =
𝑑𝑡 2 𝑑𝑡
𝑑𝜃 𝑑𝜃
⃗𝑳 = 𝑚𝑟 2 ⃗ |2
𝒊̂𝑧 = 𝑚|𝒓 𝒊̂
𝑑𝑡 𝑑𝑡 𝑧
Al combinar la acción de la fuerza y la conservación del momentum angular tenemos:

−𝐷 𝑑𝜃
⃗ =
𝒂 𝒊̂ ⃗
𝒓 , ⃗𝑳 = 𝑚|𝒓
⃗ | 2
𝒊̂
⃗ |2
𝑚|𝒓 𝑑𝑡 𝑧

−𝐷 𝑑𝜃 −𝐷𝑚|𝒓 ⃗ |2 𝑑𝜃 𝑑𝜃
⃗ × ⃗𝑳 =
𝒂 𝒊̂ × ⃗
𝑚|𝒓 | 2
𝒊̂ = (𝒊̂ × 𝒊̂ ) = 𝐷 (𝒊̂ × 𝒊̂𝒓⃗ )
⃗ |2 𝒓⃗
𝑚|𝒓 𝑑𝑡 𝑧 ⃗ |2 𝑑𝑡 𝒓⃗
𝑚|𝒓 𝑧
𝑑𝑡 𝑧

𝑑𝜃 𝑑𝒊̂𝒓⃗ 𝑑𝒊̂𝒓⃗ ⃗
𝑑𝒗 𝑑𝒊̂𝒓⃗
⃗𝒂 × ⃗𝑳 = 𝐷 =𝐷 ⟹ × ⃗𝑳 = 𝐷
𝑑𝑡 𝑑𝜃 𝑑𝑡 𝑑𝑡 𝑑𝑡

Profesor: Julián Fernando Ruiz Roa.

Septiembre 02 2015
2|7
 Física de Campos.
Notas de Clase Interacción Gravitacional y Campo Gravitacional 2015-02 Presencial.
⃗
𝑑𝒗 𝑑 𝑑𝒊̂𝒓⃗
⟹ × ⃗𝑳 = (𝒗⃗ × ⃗𝑳) = 𝐷
𝑑𝑡 𝑑𝑡 𝑑𝑡
Al integrar a ambos lados de la igualdad tenemos:

⃗ ) = 𝐷(𝒊̂𝒓⃗ + 𝒆
⃗ ×𝑳
(𝒗 ⃗)

⃗ es un vector constante y cuya magnitud es la excentricidad.

Donde 𝒆

Ahora al hacer una proyección con el vector posición tenemos:

⃗ × ⃗𝑳) = 𝒓
⃗ ∙ (𝒗
𝒓 ⃗ ∙ [𝐷(𝒊̂𝒓⃗ + 𝒆
⃗ )] ⟹ (𝒓 ⃗ ) ∙ ⃗𝑳 = 𝐷𝒓
⃗ ×𝒗 ⃗ ∙ (𝒊̂𝒓⃗ + 𝒆
⃗)

𝑚 1 1
⟹ (𝒓
⃗ ×𝒗 ⃗ = (𝒓
⃗ )∙𝑳 ⃗ × (𝑚𝒗 ⃗ = (𝒓
⃗ )) ∙ 𝑳 ⃗ = 𝐷𝒓
⃗ )∙𝑳
⃗ ×𝒑 ⃗ ∙ (𝒊̂𝒓⃗ + 𝒆
⃗)
𝑚 𝑚 𝑚
2
1 1 2 ⃗|
|𝑳
⃗ ⃗ ⃗ (𝒊̂
⃗ ∙ 𝒓⃗ + 𝒆
⟹ 𝑳 ∙ 𝑳 = |𝑳| = 𝐷𝒓 )
⃗ ⟹𝒓 (𝒊̂
⃗ ∙ 𝒓⃗ + 𝒆)
⃗ =
𝑚 𝑚 𝑚𝐷
2
⃗|
|𝑳
⃗ ∙ (𝒊̂𝒓⃗ + 𝒆
⟹𝒓 ⃗)=𝒓
⃗ ∙ 𝒊̂𝒓⃗ + 𝒓
⃗ ∙𝒆
⃗ = 𝑟 + 𝑟𝑒 cos(𝜃) =
𝑚𝐷
2 2
⃗|
|𝑳 ⃗ | ⁄𝐷𝑚)
(|𝑳
⟹ 𝑟(1 + 𝑒 cos(𝜃)) = ⟹𝑟=
𝑚𝐷 1 + 𝑒 cos(𝜃)
Pero cuál posible trayectoria opta por tomar, la respuesta está en la energía,
recordemos que la energía potencial gravitacional es:

−𝒢𝑚𝑚′
𝑈𝑔 (𝒓
⃗)= [J]
|𝒓 ⃗ ′|
⃗ −𝒓

Ahora si la Fuente de Campo Gravitacional está ubicada en el centro del marco de

referencia o sistema coordenado tendríamos:

⃗ ′ = 0𝒊̂𝑥 + 0𝒊̂𝑦 + 0𝒊̂𝑧 = 𝓸

𝒓 ⃗⃗

−𝒢𝑚𝑚′ −𝒢𝑚𝑚′
𝑈𝑔 (𝒓
⃗)= , 𝑈𝑔 (𝒓
⃗)=
|𝒓
⃗| 𝑟
−𝐷 −𝐷
𝑈𝑔 (𝒓
⃗)= , 𝑈𝑔 (𝒓
⃗)= 𝐷 = 𝒢𝑚𝑚′
|𝒓
⃗| 𝑟
Profesor: Julián Fernando Ruiz Roa.
Septiembre 04 2015
3|7
Universidad de Antioquia.
Instituto de Física. Facultad de Ingenierías.
Al remplazar la ecuación de la órbita en la definición de energía potencial
gravitacional tenemos:

2
⃗ | ⁄𝐷𝑚)
(|𝑳 1 1 𝐷𝑚[1 + 𝑒 cos(𝜃)]
𝑟= ⟹ = = 2
1 + 𝑒 cos(𝜃) 𝑟 |𝒓⃗| ⃗|
|𝑳

−𝐷 𝐷𝑚[1 + 𝑒 cos(𝜃)]
𝑈𝑔 (𝒓
⃗)= ⟹ 𝑈𝑔 = −𝐷 { 2 }
|𝒓
⃗| ⃗
|𝑳|

𝐷2 𝑚
⟹ 𝑈𝑔 (𝜃) = − 2
[1 + 𝑒 cos(𝜃)]
⃗|
|𝑳
Ahora calculemos la energía cinética:

1
𝐾𝑇𝑅 = 𝑚|𝒗 ⃗ × ⃗𝑳) = 𝐷(𝒊̂𝒓⃗ + 𝒆
⃗ |2 , (𝒗 ⃗)
2
De la ecuación de la excentricidad tenemos:

⃗ × ⃗𝑳) ∙ (𝒗
(𝒗 ⃗ × ⃗𝑳) = [𝐷(𝒊̂𝒓⃗ + 𝒆
⃗ )] ∙ [𝐷(𝒊̂𝒓⃗ + 𝒆
⃗ )]
2
⃗ × ⃗𝑳| = 𝐷2 (𝒊̂𝒓⃗ + 𝒆
⟹ |𝒗 ⃗ )2
2 2 2
⃗ | = [|𝒗
⃗ ×𝑳
⟹ |𝒗 ⃗ |(𝒊̂𝜃 × 𝒊̂𝑧 )] = |𝒗
⃗ ||𝑳 ⃗ | = 𝑣 2 𝐿2
⃗ |2 |𝑳

2 𝐷2 (𝒊̂𝒓⃗ + 𝒆
⃗ )2
⟹ ⃗ |2
|𝒗 ⃗| =𝐷
|𝑳 2 |𝒊̂
⃗
𝒓 ⃗
+𝒆 |2 ⟹ ⃗ |2
|𝒗 = 2
⃗|
|𝑳

𝐷2 𝐷2
⃗ |2 =
⟹ |𝒗 2
[(𝒊̂𝒓⃗ + 𝒆
⃗ ) ∙ (𝒊̂𝒓⃗ + 𝒆
⃗ )] = 2
(𝒊̂𝒓⃗ ∙ 𝒊̂𝒓⃗ + 2𝒊̂𝒓⃗ ∙ 𝒆
⃗ +𝒆 ⃗)
⃗ ∙𝒆
⃗|
|𝑳 ⃗|
|𝑳

𝐷2
⟹ ⃗ |2
|𝒗 = 2
(1 + 2𝑒 cos(𝜃) + 𝑒 2 )
⃗|
|𝑳

𝐷2 𝑚
𝐾𝑇𝑅 (𝜃) = 2
(1 + 2𝑒 cos(𝜃) + 𝑒 2 )
⃗|
2|𝑳

Profesor: Julián Fernando Ruiz Roa.

Septiembre 02 2015
4|7
 Física de Campos.
Notas de Clase Interacción Gravitacional y Campo Gravitacional 2015-02 Presencial.
Como la energía mecánica se conserva tenemos que:

𝐸(𝜃) = 𝐾𝑇𝑅 (𝜃) + 𝑈𝑔 (𝜃)

𝐷2 𝑚 2)
𝐷2 𝑚
⟹ 𝐸(𝜃) = 2
(1 + 2𝑒 cos(𝜃) + 𝑒 − 2
[1 + 𝑒 cos(𝜃)]
⃗|
2|𝑳 ⃗|
|𝑳

𝐷2 𝑚 1 𝑒2
⟹ 𝐸(𝜃) = 2 [2 + 𝑒 cos(𝜃) + 2 − 1 − 𝑒 cos(𝜃)]
⃗
|𝑳|

𝐷2 𝑚
⟹𝐸= 2
(𝑒 2 − 1) ∶ 𝑐𝑡𝑒
⃗|
2|𝑳
Como puede darse cuenta la relación entre la excentricidad y la energía es una
parábola desplazada.

𝑒 = 0 Circunferencia 𝐸 negativa y mínima

𝑒<1 Elipse 𝐸 negativa
𝑒=1 Parabola 𝐸 es cero
𝑒>1 Hiperbola 𝐸 positiva

Por tanto el momentum angular determina la excentricidad y la excentricidad

determina la órbita.

Ahora demostraremos la tercera ley de Kepler, partiendo de su segunda ley:

⃗⃗
𝑑𝑨 1 ⃗
𝑑𝒗 ⃗
1 𝑳
⃗ × )= ( )
= (𝒓
𝑑𝑡 2 𝑑𝑡 2 𝑚

La integramos con respecto al tiempo, durante un periodo completo:

⃗⃗ 𝑇
𝑨 𝑇⃗𝑳 ⃗𝑨⃗ ⃗𝑳 ⃗𝑳
∫ 𝑑𝑨⃗⃗ = ∫ ( ) 𝑑𝑡 ⟹ 𝑨
⃗⃗ | = ( ) 𝑡| ⟹ (𝑨
⃗⃗ − 𝓸
⃗
⃗ ) = ( ) (𝑇 − 0)
⃗⃗
𝓸 0 2𝑚 ⃗⃗
𝓸 2𝑚 0
2𝑚

⃗𝑳𝑇
⃗⃗ =
⟹𝑨
2𝑚
Luego:

⃗ |𝑇
|𝑳 2|𝑨⃗⃗ |
⃗⃗ | ⟹ 𝑇 =
= |𝑨
2𝑚 ⃗ |⁄ 𝑚 )
(|𝑳
Profesor: Julián Fernando Ruiz Roa.
Septiembre 04 2015
5|7
Universidad de Antioquia.
Instituto de Física. Facultad de Ingenierías.

Un caso particular, sistema solar circular como un cuerpo rígido, período 𝑇 constante
para todos los planetas.

Otro caso particular es que el sistema solar circular se comportara como un vórtice
𝑇 ∝ 𝑎2

Ahora tomemos la ecuación de trayectoria:

2
⃗ | ⁄𝐷𝑚)
(|𝑳
𝑟=
1 + 𝑒 cos(𝜃)
La generalizada:

𝑒𝑑
𝑟=
1 + 𝑒 cos(𝜃)
Que expresada en función del semieje mayor es:

𝑎(1 − 𝑒 2 )
𝑟=
1 + 𝑒 cos(𝜃)

Al igualar esta última con la primera:

2
⃗ | ⁄𝐷𝑚)
(|𝑳 𝑎(1 − 𝑒 2 ) 2
= ⃗ | ⁄𝐷𝑚) = 𝑎(1 − 𝑒 2 )
⟹ (|𝑳
1 + 𝑒 cos(𝜃) 1 + 𝑒 cos(𝜃)
2 2
⃗|
|𝑳 ⃗|
𝑚 |𝑳 ⃗|
|𝑳 𝐷𝑎
⟹ = ( ) = 𝑎(1 − 𝑒 2 ) ⟹ = √ (1 − 𝑒 2 )
𝐷𝑚 𝐷 𝑚 𝑚 𝑚

Luego el periodo será:

2|𝑨⃗⃗ | ⃗⃗ |
2|𝑨
⟹𝑇= ⟹𝑇=
⃗ |⁄ 𝑚 )
(|𝑳 √(𝐷𝑎⁄𝑚)(1 − 𝑒 2 )

Área de la elipse se define como:

⃗⃗ | = 𝐴 = 𝜋𝑎𝑏 = 𝜋𝑎2 √1 − 𝑒 2
|𝑨

Profesor: Julián Fernando Ruiz Roa.

Septiembre 02 2015
6|7
 Física de Campos.
Notas de Clase Interacción Gravitacional y Campo Gravitacional 2015-02 Presencial.
Al remplazar en el periodo tenemos:

⃗⃗ |
2|𝑨 2𝜋𝑎2 √1 − 𝑒 2
⟹𝑇= =
√(𝐷𝑎⁄𝑚)(1 − 𝑒 2 ) √(𝐷𝑎⁄𝑚)√1 − 𝑒 2

2𝜋𝑎2 2𝜋
⟹𝑇= ⟹𝑇= 𝑎3⁄2
√(𝐷⁄𝑚)𝑎1⁄2 √(𝐷⁄𝑚)

4𝜋 2 3
2 2
4𝜋 2 3
⟹𝑇 = 𝑎 ⟹𝑇 = 𝑎
(𝐷⁄𝑚) 𝒢𝑚′

Profesor: Julián Fernando Ruiz Roa.

Septiembre 04 2015
7|7