Ecuación de la trayectoria http://www.sc.ehu.es/sbweb/fisica3/celeste/ecuaci...

Inicio Dinámica celeste Fuerza central y conservativa

Ecuación de la trayectoria
Las fuerza de interacción gravitatoria es central y conservativa. Por tanto, la energía y el
momento angular se mantienen constantes en todos los puntos de la trayectoria

Expresando el vector posición, el vector velocidad, la energía y el momento angular de la

partícula en coordenadas polares, llegamos a la ecuación de la trayectoria, una cónica:
elipse, parábola o hipérbola

La energía y el momento angular en coordenadas

polares
La posición del punto P es

x=r·cosθ
y=r·sinθ

Expresamos la velocidad de la partícula en coordenadas

polares

1 of 9 7/16/20, 7:16 PM
Ecuación de la trayectoria http://www.sc.ehu.es/sbweb/fisica3/celeste/ecuaci...

→
→ dr dr dr̂
v = = r̂ +r
dt dt dt
Calculamos las componentes rectangulares de los vectores unitarios r̂, θ̂

r̂ = î cos θ + ĵ sin θ
θ̂ = − î sin θ + ĵ cos θ
vemos que

dr̂ dθ dθ dθ
= − sin θ î + cos θ ĵ = θ̂
dt dt dt dt
Las componentes del vector velocidad en coordenadas polares son, por tanto

→ dr dθ
v = r̂ + r θ̂
dt dt
La expresión de la energía en coordenadas polares es

1
E= 2
mv2 − GMm
r
2 2
E= 1
2
m( dr
dt
) + 1
2
mr (
2 dθ
dt
) − GMm
r

Donde -GMm/r es la energía potencial correspondiente a la fuerza conservativa F=GMm/r2.

Expresamos el momento angular en coordenadas polares

→ →
L = r × m v = m r × ( r̂ + r θ̂) = mr2 k̂
→ → dr dθ dθ
dt dt dt
Despejamos dθ /dt en la expresión del momento angular y la introducimos en la expresión de
la energía. Tenemos dos ecuaciones
2
E= 1
2
m( dr
dt
) + L2
2mr2
− GMm
r

L = mr2 dθ
dt

Movimiento en la dirección radial

Analizamos la primera ecuación de la energía. Decimos que la partícula se mueve en una
región unidimensional r>0 bajo un potencial efectivo, cuyo aspecto se muestra en la figura

L2 GMm
Ep (r) = −
2mr2 r
Representamos una función de la forma f(r)=a/r2-b/r con a=1 y b=5. Calculamos los puntos de

2 of 9 7/16/20, 7:16 PM
Ecuación de la trayectoria http://www.sc.ehu.es/sbweb/fisica3/celeste/ecuaci...

intersección con la recta horizontal e=-2

a=1;
b=5;
e=-2;
f=@(r) a./r.^2-b./r;
fplot(f,[0.2,3],'color','r')
line([0,3],[e,e],'color','k')

%raíces
r1=(-b+sqrt(b^2+4*a*e))/(2*e);
r2=(-b-sqrt(b^2+4*a*e))/(2*e);
line([r1,r1],[0,e], 'lineStyle','--')
line([r2,r2],[0,e],'lineStyle','--')
grid on
xlabel('r')
ylabel('V(r)')
title('Energía potencial')

El mínimo se encuentra en la posición rm=2a/b

>> 2*a/b
ans = 0.4000

La energía E de la partícula puede ser positiva o negativa. El valor de la energía total no puede
ser menor que el mínimo de la energía potencial efectiva.

El valor mínimo de obtiene derivando Ep(r) respecto a r e igualando a cero

3 of 9 7/16/20, 7:16 PM
Ecuación de la trayectoria http://www.sc.ehu.es/sbweb/fisica3/celeste/ecuaci...

L2 GMm L2
− + =0 rc =
mr3 r2 GMm2
El cuerpo describe una trayectoria circular de radio rc con velocidad vc. El momento angular
L=mrcvc.

GM
v2c =
rc
La misma expresión, que la que se deduciría a partir de la dinámica del movimiento circular
uniforme. La energía vale

GMm
E=−
2rc
Cuando la energía total E es mayor que el mínimo y menor que cero, como vemos en la
gráfica, la recta horizontal que señala la energía total corta a curva que describe la energía
potencial Ep(r), en dos puntos r1 y r2, soluciones de la ecuación de segundo grado

L2 GMm
E= −
2mr2 r
Llamando u=1/r

L2 2
2m
u − GMmu − E = 0
1±√1+ 2L2 E
G2 M 2 m3 1±ε
u1,2 = L2
= d
GMm2
d d
r1 = 1+ε
r2 = 1−ε

Donde ε es la excentricidad y d un factor de escala, que obtendremos más abajo al deducir la

ecuación de la trayectoria.

Si la energía E<0 ó ε<1 hay dos raices positivas. La partícula se mueve entre r1 y r2
Si E>0 ó ε>1 solamente hay una raíz positiva, la partícula se mueve entre r1 e ∞

Fuerza central y conservativa

Utilizamos las propiedades de la fuerza de atracción: central y conservativa, para obtener la
ecuación de la trayectoria

En la primera ecuación de la energía sustituímos dr/dt por

dr dr dθ dr L
= =
dt dθ dt dθ mr2
Despejamos dθ/dr para obtener la ecuación de la trayectoria

4 of 9 7/16/20, 7:16 PM
Ecuación de la trayectoria http://www.sc.ehu.es/sbweb/fisica3/celeste/ecuaci...

L
dθ mr2
=
dr
√ m2 (E − L2
+ r )
GMm
2mr2

Para integrar se hace el cambio u=1/r

L −du
θ − θ0 = ∫
√2m √E − L 2
2m
u + GMmu

Donde θ0 es la constante de integración que se determina a partir de las condiciones iniciales.

Tenemos una integral del tipo

−du −du
∫ =∫
√−au2 + bu + c
√a [ ac +
2
b2
4a2
− (u − b
2a
) ]

L2
a= 2m
b = GMm c=E

Hacemos un nuevo cambio de variable

= √ + 2 cos ξ
b c b2
u−
2a a 4a

Con este cambio la integral es inmediata

θ − θ0 = L 1
∫ dξ
√2m √a

θ − θ0 = ξ
Deshacemos los cambios

1
r − b
2a
= √ ac + b2
4a2
cos (θ − θ0 )
1
r = GMm2
L2
(1 + √1 + 2L2 E
G2 M 2 m3
cos (θ − θ0 ))

La ecuación de la trayectoria es una cónica (elipse, parábola o hipérbola) dependiendo del

valor de la excentricidad ε .

ε = √1 +
d 2L2 E L2
r= d=
1 + ε cos(θ − θ0 ) m3 G2 M 2 GMm2

Clase de cónica Descripción geométrica Descripción física

5 of 9 7/16/20, 7:16 PM
Ecuación de la trayectoria http://www.sc.ehu.es/sbweb/fisica3/celeste/ecuaci...

Elipse ε<1 E<0

Parábola ε=1 E=0
Hipérbola ε>1 E>0

Las secciones cónicas (elipse, hipérbola y parábola) es el resultado de la intersección de una

superficie cónica y un plano

Dibujamos las cónicas de excentricidades (0, 0.3,0.5,0.7,0.8,0.85) elipses, (1) parábola y

(1.1,1.3,1.5,2.5) hipérbolas

hold on
for e=[0,0.3,0.5,0.7,0.8,0.85,0.9,1,1.1,1.3,1.5,2.5]
f=@(x) 1./(1+e*cos(x));
ang=pi;
if e>1 %hipérbolas
ang=acos(-1/e)-pi/180;
elseif e==1 %parábola
ang=5*pi/6;
end
fplot(@(x) f(x).*cos(x), @(x) f(x).*sin(x) ,[-ang,ang],
'displayName',num2str(e))
end
legend('-DynamicLegend','Location','northwest')
plot(0,0,'ko','markersize',3,'markeredgecolor','k','markerfacecolor','k')
hold off
axis equal
axis off

6 of 9 7/16/20, 7:16 PM
Ecuación de la trayectoria http://www.sc.ehu.es/sbweb/fisica3/celeste/ecuaci...

Para ver bien las trayectorias es necesario utilizar los botones Zoom In y Pan de la ventana
gráfica

Ecuación del movimiento

En este apartado vamos a obtener la ecuación de la trayectoria aplicando la segunda ley de
Newton.

Fuerzas atractivas
En el primer apartado hemos obtenido la expresión del vector velocidad en coordenadas
polares. La expresión del vector aceleración es:

→ →
d2r
+ (r ddt2θ + ) θ̂
dv dr dr̂ 2 dr dθ dθ̂
a = dt
= dt2
r̂ + dt dt dt dt
+ r dθ
dt dt
=
2
d2r
+ (r ddt2θ + ) θ̂ ) r̂ =
dr dθ 2 dr dθ
dt2
r̂ + dt dt
θ̂ dt dt
− r( dθ
dt
2
( ddt2r − r( dθ ) ) r̂ + (r ddt2θ + 2 dr ) θ̂
2 2 dθ
dt dt dt

Sobre la partícula de masa m actúa una fuerza atractiva inversamente proporcional al

cuadrado de la distancia r al centro de fuerzas. Las ecuaciones del movimiento son

2
m ( ddt2r − r( dθ ) ) = −G Mm
2

dt r2

m (r ddt2θ + 2 dθ )
2 dr
dt dt
=0

La segunda ecuación nos indica que el momento angular L (expresado en coordenadas polares
es constante)

(mr2 ) = 0 L = mr2
d dθ dθ
= cte
dt dt dt
Teniendo en cuenta la constancia del momento angular L, escribimos la primera ecuación de la
forma

d2r L2 M
− = −G
dt2 m 2 r3 r2
Para calcular la ecuación de la trayectoria r=r(θ), eliminamos el tiempo t del siguiente modo
dr dr dθ L dr
dt
= dθ dt
= mr2 dθ
d2r
dt2
= d
( L dr )
dt mr2 dθ
= d
( L dr ) dθ
dθ mr2 dθ dt
= L d
( L dr )
mr2 dθ mr2 dθ
2
d2r
( ddθr2 − 2r ( dθ ) )
2
L2 dr
dt2
= m2 r4

7 of 9 7/16/20, 7:16 PM
Ecuación de la trayectoria http://www.sc.ehu.es/sbweb/fisica3/celeste/ecuaci...

Obtenemos una ecuación en términos de r y de sus derivadas respecto del ángulo θ

2
d2r 2 dr GMm2 2
− ( ) −r=− r
dθ2 r dθ L2
Hacemos el cambio de variable u=1/r
dr
dθ
= − u12 du
dθ
2
d2r
dθ
= 2
u3
( du
dθ
) − 1 d2u
u2 dθ2

la ecuación del movimiento se convierte en

d2u GMm2
+ u =
dθ2 L2
La solución de la ecuación diferencial es

GMm2
u = A sin θ + B cos θ +
L2
En forma equivalente, escribiendo
2 2
L L
A GMm 2 = ε sin θ0 B GMm 2 = ε cos θ0

GMm2
u= L2
(ε cos(θ − θ0 ) + 1)

Calculamos la constante ε a partir de la conservación de la energía.

2 2
E= 1
2
m {( dr
dt
) + r2 ( dθ
dt
) }− GMm
r

2
E= 1 L2
2 mr4
{( dθ
dr
) + r2 } − GMm
r

2
E= 1 L2
2 m
u4 { u14 ( du
dθ
) + 1
u2
} − GMmu
2
E= 1 L2
2 m
{( du
dθ
) + u2 } − GMmu

Introduciendo la expresión de u y su derivada du/dθ, obtenemos

ε = √1 +
2L2
E
G2 M 2 m3

La ecuación de la trayectoria es r=1/u

8 of 9 7/16/20, 7:16 PM
Ecuación de la trayectoria http://www.sc.ehu.es/sbweb/fisica3/celeste/ecuaci...

1 GMm2
= (ε cos(θ − θ0 ) + 1)
r L2
Fuerzas repulsivas
Sobre la partícula actúa una fuerza repulsiva inversamente proporcional al cuadrado de la
distancia r al centro de fuerzas. Por ejemplo, un núcleo fijo de carga Q que repele a una
partícula alfa de masa m y carga q. Las ecuaciones del movimiento son

2
m ( ddt2r − r( dθ ) )=
2 k Qq
dt r2
k= 4πε0

m (r ddt2θ + 2 dθ )
2 dr
dt dt
=0

Siguiendo los mismos pasos que para la fuerza atractiva, la solución es

km
u = A sin θ + B cos θ −
L2
De forma equivalente

km
u= (ε cos(θ − θ0 ) − 1)
L2
Calculamos la constante ε a partir de la conservación de la energía.

2 2
E= 1
2
m {( dr
dt
) + r2 ( dθ
dt
) }+ k
r

2
E= 1 L2
2 m
{( du
dθ
) + u2 } + ku

Introduciendo la expresión de u y su derivada du/dθ, obtenemos

ε = √1 +
2L2
E
k2 m

La ecuación de la trayectoria es r=1/u

1 km
= 2 (ε cos(θ − θ0 ) − 1)
r L

Energías Renovables   Curso Interactivo de Física en Internet

Angel Franco García, Copyright © 2016

9 of 9 7/16/20, 7:16 PM