/

Simulado 3; Movimiento Unidimensional Equivalente para el

Potencial Kepleriano
Una función de energía potencial utilizada para describir la interacción gravitacional entre
dos masa
k
V=- ;
r
donde r es la separación entre las dos masas.
A) Genere un diagrama de energía del sistema para diferentes energías, y clasifique los
diferentes tipos de movimiento para dichas energías
B) Genere una animación que muestre un diagrama de energía del sistema y la forma
cómo evoluciona r para el movimiento ligado (acotado)
C) Halle la energía necesaria para que el movimiento sea circular uniforme

Solución A)
1  l2
m r2 +
k
La energía de este sistema es E = 2 2 μ r2
- r ; el potencial efectivo es

Lz2 k
Vef (r) = - ;
2 m r2 r
Adimensionando el problema:
1 1
Vef (r) = - ;
2 r2 r

1 1 1
In[1]:= Vef[r_] := - ; energia[r_][Vr_] := Vr2 + Vef[r]
2 r2 r 2
1 1
In[2]:= E1 → energia [2], E2 → energia[1][1], E3 → energia[2] 
2 2
1
Out[2]= E1 → 2, E2 → 0, E3 → - 
4
2 SIM 3 MUE PARA [Link]

1 1
In[3]:= PlotVef[r], energia [2], energia[1][1], energia[2] ,
2 2
{r, 0, 6}, AxesLabel → {"r(?)", "Vef (?)"}, AxesOrigin → {0, 0},
PlotStyle → {Automatic, Dashed, Dashed, Dashed},
PlotLabel → "DIAGRAMA DE ENERGIA",
PlotRange → {{0, 4}, {-0.5, 2.5}}
DIAGRAMA DE ENERGIA
Vef (?)
2.5

2.0

1.5
Out[3]=

1.0

0.5

r(?)
1 2 3 4

-0.5

TIPOS DE MOVIMIENTO:
Condiciones Iniciales (r0 = 1/2, Vr0 = -2 ); E = 2 > 0; Hay un solo punto de retorno y el movimiento no
es ligado (No acotado).
Condiciones Iniciales (r0 = 1, Vr0 = 1 ); E = 0 ; Hay un solo punto de retorno y el movimiento no es
ligado (No acotado)
Condiciones Iniciales (r0 = 2, Vr0 = 1/2 ); E = -1/4 < 0; Hay dos puntos de retorno y el movimiento es
ligado ( acotado)

Solución B)
La fuerza se define como F = - d Vef ;
dr

La posición en función del tiempo la obtenemos integrando la segunda ley

d2 r d2 r d Vef
= F[r]; =- ;
dt dt dr
In[4]:= ecuMovi = r''[t] ⩵ -D[Vef[r[t]], r[t]]
1 1
Out[4]= r′′ [t] ⩵ -
r[t]3 r[t]2

Solución numerica para E = -1/4, esto es para cuando r[0] = 2, r [0] = 1/2
 SIM 3 MUE PARA [Link] 3

In[5]:= sol =
1 1 1
NDSolver′′ [t] ⩵ - , r[0] ⩵ 2, r'[0] ⩵ , r, {t, 0, 100}
r[t]3 r[t]2 2

Out[5]= r → InterpolatingFunction Domain: {{0., 100.}} 

Output: scalar

In[6]:= rr[t_] := r[t] /. sol[[1]];

Puntos de retorno
-1 1
In[7]:= FindRoot ⩵ Vef[r], {r, 0.5}, FindRoot- ⩵ Vef[r], {r, 4}
4 4
Out[7]= {{r → 0.585786}, {r → 3.41421}}

In[8]:= r1 = 0.585786437626905`; r2 = 3.414213562373095`;

1
In[11]:= diagrama = PlotVef[r], energia[2] ,
2
{r, 0, 6}, AxesLabel → {"r(?)", "V (?)"},
AxesOrigin → {0, 0}, PlotStyle → {Automatic, Dashed},
Ticks → {{{r1, "r1 "}, {r2, "r2 "}}, Automatic},
PlotLabel → "DIAGRAMA DE ENERGIA"
DIAGRAMA DE ENERGIA
V (?)

r(?)
r1 r2

-0.1

Out[11]= -0.2

-0.3

-0.4

-0.5

In[12]:= MUE[t_] := Show[{diagrama, Graphics[{Red, PointSize[0.03],

Point[{rr[t], 0}], Dashed, Line[{{rr[t], 0}, {rr[t], Vef[rr[t]]}}]}]}];
4 SIM 3 MUE PARA [Link]

In[13]:= Manipulate[MUE[t], {t, 0, 50, 0.1}]

DIAGRAMA DE ENERGIA
V (?)

r(?)
r1 r2

-0.1
Out[13]=

-0.2

-0.3

-0.4

-0.5

Movimiento Visto desde un Marco de Referencia Inercial (MRI)

Recuerde que eI Movimiento Unidimensional Equivalente, es el que percibe un obser-
vador que se encuentra en un sistema de referencia no inercial, que está rotando con

respecto a un MRI con velocidad angular θ. Para determinar el movimiento visto desde un
MRI integremos la ecuación para la coordenada θ
 1
P θ = r 2 θ = L z ; Lz = r vr = 2 * = 1
2
In[14]:= ecuPθ = θ '[t] ⩵ 1  rr[t]2 ;

In[15]:= solucionθ = NDSolve[{ecuPθ, θ[0] ⩵ 0}, θ, {t, 0, 100}]

Out[15]= θ → InterpolatingFunction Domain: {{0., 100.}} 

Output: scalar

In[16]:= θθ[t_] := θ[t] /. solucionθ[[1]]

In[17]:= x[t_] := rr[t] Cos[θθ[t]]; y[t_] := rr[t] Sin[θθ[t]];

 SIM 3 MUE PARA [Link] 5

In[18]:= trayectoria = ParametricPlot[Evaluate[{x[t], y[t]}],

{t, 0, 80}, AxesLabel → {"X", "Y"}, PlotLabel → "Trayectoria de μ"]
Trayectoria de μ
Y

2.5

2.0

1.5

Out[18]=
1.0

0.5

X
-0.5 0.5 1.0 1.5 2.0 2.5

-0.5

In[19]:= orbita[t_] :=
Show[{trayectoria, Graphics[{Red, PointSize[0.03], Point[{x[t], y[t]}]}]}];
6 SIM 3 MUE PARA [Link]

In[20]:= Manipulate[orbita[t], {t, 0, 80, 0.1}]

Trayectoria de μ
Y

2.5

2.0

1.5
Out[20]=

1.0

0.5

X
-0.5 0.5 1.0 1.5 2.0 2.5

-0.5

Solución C) Movimiento Circular

La condición para que se dé este movimiento es que la energía coincida con el mínimo del potencial
efectivo
ⅆVef
E = Vef [R]; donde =0 (1)
ⅆr R

In[21]:= ecu3 = D[Vef[r], r] ⩵ 0 // Simplify

-1 + r
Out[21]= ⩵0
r

In[22]:= Solve[ecu3, r]
Out[22]= {{r → 1}}

In[23]:= rmin = 1; e = Vef[rmin]

1
Out[23]= -
2
 SIM 3 MUE PARA [Link] 7

In[24]:= diagrama = Plot{Vef[r], e}, {r, 0, 6}, AxesLabel → {"r(?)", "V (?)"},
AxesOrigin → {0, 0}, PlotStyle → {Automatic, Dashed},
Ticks → {{{1, "R"}}, Automatic}, PlotLabel → "DIAGRAMA DE ENERGÍA"
DIAGRAMA DE ENERGÍA
V (?)

0.2

r(?)
R
Out[24]=

-0.2

-0.4

El movimiento unidimensional equivalente, para estas condiciones iniciales, nos muestra que la solu-

ción vista por un observador que rota con velocidad angular θ, para la coordenada radial, es r[t] = rmin ,
esto es, que esta coordenada se mantiene constante a medida que transcurre el tiempo. Integremos
ahora la ecuación para la coordenada θ

2
  1
Pθ = rmin θ = 1; θ = = 1; θ[t] = θ0 + t
2
rmin
De modo que la coordenada angular cambia linealmente con el tiempo, lo que caracteriza al MCU
SOLUCIÓN NUMÉRICA DEL MCU
Partimos del lagrangiano del sistema
1 2
2 + r 2 θ k
L= μ r + (2)
2 r
Adimensionado
1 1
In[26]:= L= r '[t]2 + r[t]2 θ '[t]2  + ;
2 r[t]

In[27]:= ecur = D[D[L, r '[t]], t] - D[L, r[t]] ⩵ 0 // Simplify

1
Out[27]= + r′′ [t] ⩵ r[t] θ′ [t]2
r[t]2

In[28]:= ecuθ = D[D[L, θ '[t]], t] - D[L, θ[t]] ⩵ 0 // Simplify

Out[28]= r[t] (2 r′ [t] θ′ [t] + r[t] θ′′ [t]) ⩵ 0

Note que para que se dé el MCU, la energía se iguala al Vef y esto corresponde a las condiciones
iniciales r[0] = 1 Vr[0] = 0;

In[29]:= sist = {ecur, ecuθ}; cond = {r[0] ⩵ 1, r '[0] ⩵ 0, θ[0] ⩵ 0, θ '[0] ⩵ 1};
8 SIM 3 MUE PARA [Link]

In[30]:= solNMCU = NDSolve[sist ⋃ cond, {r, θ}, {t, 0, 50}]

Out[30]= r → InterpolatingFunction Domain: {{0., 50.}} ,

Output: scalar

θ → InterpolatingFunction Domain: {{0., 50.}} 

Output: scalar

In[31]:= rc[t_] := r[t] /. solNMCU[[1]];

In[32]:= θc[t_] := θ[t] /. solNMCU[[1]]

In[33]:= xc[t_] := rc[t] Cos[θc[t]]; yc[t_] := rc[t] Sin[θc[t]];

In[34]:= ParametricPlot[Evaluate[{xc[t], yc[t]}], {t, 0, 80},

AxesLabel → {"X", "Y"}, PlotLabel → "Trayectoria de μ"]
Trayectoria de μ
Y

1.0

0.5

Out[34]=
X
-1.0 -0.5 0.5 1.0

-0.5

-1.0