Desviación hacia el este de un cuerpo que cae

Supongamos un planeta de masa M y radio R, que tiene un movimiento de rotación con velocidad angular ω . Un observador situado en la superficie del planeta ve como cae un cuerpo de masa m desde una altura h por encima del observador.

Para un observador inercial no ligado al movimiento de rotación del planeta, el cuerpo se lanza desde una distancia r=R+h del centro del planeta con una velocidad inicial v1=ω ·r.

Desde el punto de vista de este observador, el cuerpo se mueve bajo la única influencia de la fuerza de atracción del planeta. Describirá por tanto, una órbita elíptica en uno de cuyos focos estará el centro del planeta. Los pasos de la presente discusión serán los siguientes:

Ecuación de la trayectoria elíptica

La ecuación de una elipse en coordenadas polares es

r= d 1+εcosθ

Los valores del parámetro d y de la excentricidad ε se calculan a partir de la energía E y del momento angular L de la partícula

E= 1 2 m v 1 2 GMm r 1 L=m v 1 r 1

Ejemplo: Consideremos el planeta Tierra con los siguientes datos

Supongamos que el cuerpo se deja caer desde una altura h=0.1·R=637.8 km o bien, desde una distancia r1=7.02·106 m,

Los parámetros d y ε de la trayectoria se obtienen mediante las fórmulas

d= L 2 GM m 2 =32122.8m ε= 1+ 2 L 2 E G 2 M 2 m 3 =0.9954

Conocidos los parámetros d y ε de la ecuación de la trayectoria, se obtiene r2 y r1

Cuando θ=0º, r2=16098 m
Cuando θ=180º, r1=7.02·106 m

Fuerza central y conservativa

No es necesario utilizar estas dos fórmulas para calcular la ecuación de la trayectoria, existe otro camino alternativo

Conocido r1 y v1 se calcula r2 y v2 aplicando la constancia del momento angular y de la energía

m r 1 · v 1 =m r 2 · v 2 1 2 m v 1 2 GMm r 1 = 1 2 m v 2 2 GMm r 2

Como en la ecuación de la trayectoria, el valor de r1 se obtiene para θ =180º, r 1 = d 1ε
y el valor de r2 se obtiene para θ =0º, r 2 = d 1+ε

De estas dos ecuaciones despejamos d y ε .

Intersección con la superficie del planeta

Calculamos la intersección de una elipse y una circunferencia de radio R.

Poniendo en la ecuación de la trayectoria r=R, despejamos el ángulo θ

cosθ= dR εR

Con los datos anteriores cos θi=-0.9995, lo que da un ángulo de intersección de θi =178.3º.

Tiempo que tarda en chocar con la superficie

Para calcular el tiempo que tarda el cuerpo desde que se deja caer hasta que choca con la superficie del planeta en el punto P, emplearemos la ley de las áreas.

En coordenadas polares el momento angular se expresa

L=m r 2 dθ dt

Como hemos visto el área barrida por el radio vector entre el instante t y el instante t+dt es un triángulo diferencial de área r2· /2.

El área barrida por el radio vector en el tiempo t es

A= 0 θ r 2 dθ 2 = L 2m 0 t dt

Si calculamos el área A sombreada en color azul claro, obtenemos el tiempo t.

El área A sombreada es la suma del área de un triángulo y el área de la porción de elipse de la figura.

El área del triángulo es

A 1 = R·cos(π θ i )·Rsin(π θ i ) 2 = 1 4 R 2 sin( 2(π θ i ) )

El área de la porción de elipse A2 se puede calcular sumando las áreas de los elementos infinitesimales ydx comprendidos entre x1=–a y x2= –R·cos(π -θi)+c. Donde a es el semieje mayor de la elipse a=(r1+r2)/2, y c es la semidistancia focal c=ε ·a.

La ecuación de la elipse en coordenadas rectangulares es

x 2 a 2 + y 2 b 2 =1 A 2 = x 1 x 2 ydx = x 1 x 2 1 b 1 x 2 a 2 dx = ab 2 ( arcsin x a + x a 1 x 2 a 2 ) x 1 x 2

como x1=-a, la expresión se  reduce a

A 2 = ab 2 ( π 2 +arcsin x 2 a + x 2 a 1 x 2 2 a 2 )

con x2= –R·cos(π -θi)+ε a

Si x2=+a obtenemos la mitad del área de la elipse πab/2

Ejemplo:

Siguiendo con los datos anteriores tenemos que

Solamente nos queda despejar el tiempo de la ecuación

A= L 2m t , por lo que t=391.6 s

Determinación de su desviación respecto de la dirección radial por el observador no inercial, o en rotación con el planeta.

El cuerpo después de un tiempo t choca con la superficie del planeta en P. La posición angular del punto P es πi. Mientras tanto el observador se ha desplazado al punto O cuya posición angular es ω·t.

La distancia a lo largo de la superficie del planeta entre O y P es la longitud del arco

s=R·(πi-ω·t)

Con los datos que disponemos s=11837 m

Comparación con la desviación obtenida aplicando la fórmula de la aceleración de Coriolis

Vamos a comparar la desviación hacia el este de un cuerpo que se deja caer desde una altura h en el Ecuador mediante el procedimiento explicado en esta página, con la desviación obtenida aplicando la fórmula de la aceleración de Coriolis.

  1. La fuerza de atracción es central y conservativa. La trayectoria que sigue el cuerpo en su caída es elíptica

  2. Sea h=0.01·R=63780 m. cerca de la superficie de la Tierra

    Los datos iniciales necesarios para determinar la trayectoria elíptica son:

    r1=R+h=6.44·106 m
    v1=ω ·r1=468.5 m/s

    De la constancia del momento angular y de la energía obtenemos r2 y v2. Solamente nos interesa r2=11436 m.

    El parámetro d=a(1- e2)=22831.2 m

    El ángulo θi=179.52º

    El área total A=A1+ A2=1.715·1011+2.30·109=1.74·1011

    El tiempo t que tarda en llegar el móvil a la superficie de la Tierra t=2A/(r1·v1)=115.21 s

    La desviación del cuerpo respecto del observador en la superficie terrestre.

    s=R·(πi-ω·t)=355.43 m

  3. Caída de un cuerpo descrita por un observador en rotación (no inercial).

  4. Tiempo que tarda en llegar a la superficie de la Tierra cayendo desde una altura de h=63780 m.

    Con h=gt2/2, tomando g=9.8, se obtiene t=114.1 s

    Aplicando la fórmula de la desviación para la latitud λ =0º.

    s= 1 3 gωcosλ· t 3 = 1 3 9.8· 2π 24·60·60 · t 3 =352.8m

    Obtenemos un valor muy parecido, a pesar de que g es más pequeño que 9.8 debido a la aceleración centrífuga y a que disminuye con la altura.

Cálculo con MATLAB

R=6.37e6; %radio de la Tierra
M=5.98e24; %masa de la Tierra
G=6.67e-11; %constante G
w=2*pi/(24*3600); %velocidad angular de rotación

r1=(1+0.01)*R; %un cuerpo se deja caer desde una altura de 0.01*R
v1=w*r1;

L=v1*r1; %momento angular
E=v1^2/2-G*M/r1; %energía
ex=sqrt(1+2*L^2*E/(G*M)^2); %excentricidad
d=L^2/(G*M);
a=d/(1-ex^2);  %semieje mayor de la elipse
b=a*sqrt(1-ex^2); %semieje menor
angulo=acos((d-R)/(R*ex));
area1=R^2*sin(2*(pi-angulo))/4;
x2=-R*cos(pi-angulo)+ex*a;
area2=a*b*(pi/2+asin(x2/a)+x2*sqrt(1-(x2/a)^2)/a)/2;
t=2*(area1+area2)/L; %tiempo
alcance=R*(pi-angulo-w*t); %alcance
fprintf('Alcance (m) es %3.1f y tiempo (s) de vuelo %4.1f\n'
,alcance, t)
Alcance (m) es 354.3 y tiempo (s) de vuelo 115.0

Actividades

Se pulsa el botón titulado Nuevo.

El objeto que está a una altura h por encima del observador no inercial situado en la superficie del planeta, empieza a caer. Su trayectoria para un observador inercial es una porción de una elipse, al mismo tiempo que el objeto cae el observador no inercial describe un movimiento circular. El observador inercial mide el desplazamiento de ambos durante el tiempo de caída del cuerpo.

El observador no inercial situado en la superficie del planeta mide el desplazamiento relativo del cuerpo, la longitud del arco de circunferencia entre la posición que ocupa el observador (un punto de color rojo) y el punto de impacto del cuerpo sobre la superficie del planeta.

La desviación hacia el este del cuerpo que cae es pequeña para los planetas con velocidad de rotación muy baja como Venus y es muy acusada, para planetas con elevada velocidad angular de rotación como Júpiter.