Analisis de Trayectorias
Analisis de Trayectorias
Analisis de Trayectorias
ANALISIS DE TRAYECTORIAS.
−𝒢𝑚𝑚′
⃗ 𝑔 (𝒓
𝑭 ⃗)= 𝒊̂ ′
|𝒓 ⃗ ′ |2 (𝒓⃗−𝒓⃗ )
⃗ −𝒓
−𝒢𝑚𝑚′ −𝒢𝑚𝑚′
⃗ 𝑔 (𝒓
𝑭 ⃗)= ⃗ 𝑔 (𝒓
𝒊̂ , 𝑭 ⃗)= 𝒓̂
⃗ |2 𝒓⃗
|𝒓 𝑟2
−𝐷
⃗ 𝑁𝑒𝑡𝑎 (𝒓
𝑭 ⃗ ) = 𝑚𝒂 ⃗ 𝑔 (𝒓
⃗, 𝑭 ⃗)= 𝒊̂𝒓⃗ , 𝐷 = 𝒢𝑚𝑚′
|𝒓
⃗| 2
−𝐷 ⃗
𝑑𝒗 −𝐷 𝑑2𝒓⃗ −𝐷
⃗ =
𝒂 𝒊̂ ⟹ = 𝒊̂ ⟹ = 𝒊̂
⃗ |2 𝒓⃗
𝑚|𝒓 𝑑𝑡 𝑚|𝒓⃗ |2 𝒓⃗ 𝑑𝑡 2 𝑚|𝒓⃗ |2 𝒓⃗
𝑑2𝒓⃗ −𝒢𝑚′
⃗⃗ (𝒓
=𝒈 ⃗ ); 𝒈
⃗⃗ (𝒓
⃗)= 𝒊̂
𝑑𝑡 2 ⃗ |2 𝒓⃗
|𝒓
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 𝑚
𝑒𝑑
𝑟=
1 + 𝑒 cos(𝜃)
⃗⃗
𝑑𝑨 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 𝑑𝑡 𝒓⃗
𝑚|𝒓 𝑧
𝑑𝑡 𝑧
𝑑𝜃 𝑑𝒊̂𝒓⃗ 𝑑𝒊̂𝒓⃗ ⃗
𝑑𝒗 𝑑𝒊̂𝒓⃗
⃗𝒂 × ⃗𝑳 = 𝐷 =𝐷 ⟹ × ⃗𝑳 = 𝐷
𝑑𝑡 𝑑𝜃 𝑑𝑡 𝑑𝑡 𝑑𝑡
⃗ ) = 𝐷(𝒊̂𝒓⃗ + 𝒆
⃗ ×𝑳
(𝒗 ⃗)
⃗ × ⃗𝑳) = 𝒓
⃗ ∙ (𝒗
𝒓 ⃗ ∙ [𝐷(𝒊̂𝒓⃗ + 𝒆
⃗ )] ⟹ (𝒓 ⃗ ) ∙ ⃗𝑳 = 𝐷𝒓
⃗ ×𝒗 ⃗ ∙ (𝒊̂𝒓⃗ + 𝒆
⃗)
𝑚 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]
|𝒓 ⃗ ′|
⃗ −𝒓
−𝒢𝑚𝑚′ −𝒢𝑚𝑚′
𝑈𝑔 (𝒓
⃗)= , 𝑈𝑔 (𝒓
⃗)=
|𝒓
⃗| 𝑟
−𝐷 −𝐷
𝑈𝑔 (𝒓
⃗)= , 𝑈𝑔 (𝒓
⃗)= 𝐷 = 𝒢𝑚𝑚′
|𝒓
⃗| 𝑟
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|𝑳
𝐷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.
⃗⃗
𝑑𝑨 1 ⃗
𝑑𝒗 ⃗
1 𝑳
⃗ × )= ( )
= (𝒓
𝑑𝑡 2 𝑑𝑡 2 𝑚
⃗⃗ 𝑇
𝑨 𝑇⃗𝑳 ⃗𝑨⃗ ⃗𝑳 ⃗𝑳
∫ 𝑑𝑨⃗⃗ = ∫ ( ) 𝑑𝑡 ⟹ 𝑨
⃗⃗ | = ( ) 𝑡| ⟹ (𝑨
⃗⃗ − 𝓸
⃗
⃗ ) = ( ) (𝑇 − 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
2
⃗ | ⁄𝐷𝑚)
(|𝑳
𝑟=
1 + 𝑒 cos(𝜃)
La generalizada:
𝑒𝑑
𝑟=
1 + 𝑒 cos(𝜃)
Que expresada en función del semieje mayor es:
𝑎(1 − 𝑒 2 )
𝑟=
1 + 𝑒 cos(𝜃)
2
⃗ | ⁄𝐷𝑚)
(|𝑳 𝑎(1 − 𝑒 2 ) 2
= ⃗ | ⁄𝐷𝑚) = 𝑎(1 − 𝑒 2 )
⟹ (|𝑳
1 + 𝑒 cos(𝜃) 1 + 𝑒 cos(𝜃)
2 2
⃗|
|𝑳 ⃗|
𝑚 |𝑳 ⃗|
|𝑳 𝐷𝑎
⟹ = ( ) = 𝑎(1 − 𝑒 2 ) ⟹ = √ (1 − 𝑒 2 )
𝐷𝑚 𝐷 𝑚 𝑚 𝑚
2|𝑨⃗⃗ | ⃗⃗ |
2|𝑨
⟹𝑇= ⟹𝑇=
⃗ |⁄ 𝑚 )
(|𝑳 √(𝐷𝑎⁄𝑚)(1 − 𝑒 2 )
⃗⃗ | = 𝐴 = 𝜋𝑎𝑏 = 𝜋𝑎2 √1 − 𝑒 2
|𝑨
⃗⃗ |
2|𝑨 2𝜋𝑎2 √1 − 𝑒 2
⟹𝑇= =
√(𝐷𝑎⁄𝑚)(1 − 𝑒 2 ) √(𝐷𝑎⁄𝑚)√1 − 𝑒 2
2𝜋𝑎2 2𝜋
⟹𝑇= ⟹𝑇= 𝑎3⁄2
√(𝐷⁄𝑚)𝑎1⁄2 √(𝐷⁄𝑚)
4𝜋 2 3
2 2
4𝜋 2 3
⟹𝑇 = 𝑎 ⟹𝑇 = 𝑎
(𝐷⁄𝑚) 𝒢𝑚′