El péndulo de Foucault

A la izquierda, una fotografía del péndulo de Foucault para una demostración en una Feria de la Ciencia. El péndulo de Foucault del Museo de las Artes y las Ciencias de Valencia es mucho más largo

En la figura (más abajo), se muestra el ángulo girado por el plano de oscilación del "péndulo" durante el periodo de una oscilación. El péndulo parte de A y regresa a B, para iniciar una nueva oscilación. El ángulo girado es Δθ =Ω·P. Siendo P=p el periodo de una oscilación

El ángulo girado por el plano de oscilación del péndulo en una hora, es el producto de Δθ por el número de oscilaciones que da el péndulo en una hora.

Δθ·60·60/P=Ω·60·60=15º a la hora

Teniendo en cuenta que la velocidad angular de rotación Ω de la Tierra es de 360º en 24 h.

Para un lugar de latitud λ, el ángulo girado por el plano de oscilación del péndulo en una hora vale 15º·sinλ. La razón estriba en que el vector velocidad angular de rotación Ω forman un ángulo 90º-λ con la dirección perpendicular al plano local, tal como se ve en la figura. Recuérdese que la aceleración de Coriolis responsable de este fenómeno es el producto vectorial 2 Ω × v

Sabiendo que la latitud de Paris es de aproximadamente 49º, el plano de oscilación del péndulo de Foucault gira a razón de 11.3º cada hora.

Para describir el péndulo de Foucault, se pueden seguir los mismos pasos que para el péndulo esférico, salvo que el ahora está situado en un Sistema de Referencia Local, no inercial.

Sin embargo, el péndulo de Foucault es muy largo y el ángulo de desviación es pequeño por lo que es aplicable la aproximación lineal. En esta aproximación, el péndulo se comporta de forma análoga a una masa unida a un muelle elástico. Este sistema es más sencillo y ya hemos estudiado en una página anterior, el comportamiento de una masa unida a un muelle elástico situado sobre una plataforma que gira con velocidad angular constante

Ecuaciones del movimiento

En la página titulada Movimiento relativo en el Sistema de Referencia Local estudiamos el movimiento de una partícula en el Sistema de Referencia Local, llegando a las siguientes ecuaciones:

m d 2 x d t 2 2mΩ dy dt cosαm Ω 2 ( xcosα+(R+z)sinα )cosα+ V x =0 m d 2 y d t 2 +2mΩ( dx dt cosα+ dz dt sinα )m Ω 2 y+ V y =0 m d 2 z d t 2 2mΩ dy dt sinαm Ω 2 ( xcosα+(R+z)sinα )sinα+ V z =0

Consideremos una partícula de masa m unida a uno de los extremos de un muelle elástico de constante k mientras que el otro está sujeto al origen del Sistema de Referencia Local. El muelle está en el plano XY, su energía potencial será

V(x,y)= 1 2 k( x 2 + y 2 )= 1 2 m ω 2 ( x 2 + y 2 )

La partícula se mueve en el plano XY, con z=0 y dz/dt=0, las ecuaciones del movimiento son ahora

d 2 x d t 2 2Ω dy dt cosα+( ω 2 Ω 2 cosα )x= 1 2 Ω 2 Rsin(2α) d 2 y d t 2 +2Ω dx dt cosα+( ω 2 Ω 2 )y=0

Eliminamos el término constante en la primera ecuación, con un cambio de variable x'=x+c. Resolvemos el sistema de dos ecuaciones diferenciales con las siguientes condiciones iniciales

t=0{ dx dt = x ˙ 0 dy dt = y ˙ 0 x= x 0 y= y 0

w=1; %muelle elástico
W=w/10; %velocidad angular de rotación
alfa=pi/6; %posición del observador
x0=[1,0,0,0]; %condiciones iniciales
tspan=[0,50];
% x(1)=x, x(2)=dx/dt,x(3)=y, x(4)=dy/dt
fg=@(t,x)[x(2); 2*W*x(4)*cos(alfa)+(W^2*cos(alfa)-w^2)*x(1); x(4);
 -2*W*x(2)*cos(alfa)+(W^2-w^2)*x(3)];
[t,x]=ode45(fg,tspan,x0);
 
plot(x(:,1),x(:,3))
grid on
axis equal
xlabel('x')
ylabel('y')
title('Muelle en el plano local')

Aproximaciones

Si despreciamos los términos en Ω2, el sistema de ecuaciones diferenciales se puede resolver de forma análoga al caso del muelle sobre una plataforma en rotación

d 2 x d t 2 2Ω dy dt cosα+ ω 2 x=0 d 2 y d t 2 +2Ω dx dt cosα+ ω 2 y=0

Utilizando la notación compleja ξ=x+iy, el sistema de dos ecuaciones diferenciales se transforma en una única ecuación diferencial en ξ

d 2 ξ d t 2 +2Ωi dξ dt cosα+ ω 2 ξ=0

Las raíces de la ecuación característica son

s 2 +2Ωicosα·s+ ω 2 =0{ s 1 =i( Ωcosα+ ω 2 + Ω 2 cos 2 α ) s 2 =i( Ωcosα ω 2 + Ω 2 cos 2 α )

Escribimos las raíces s1=i(-a+b) y s2=-i(a+b)

La solución de la ecuación diferencial lineal homogénea es

ξ=Cexp( i(a+b)t )+Dexp( i(a+b)t )=( Cexp( ibt )+Dexp( ibt ) )exp( iat )

Un número complejo se expresa en forma polar, ξ=x+iy=r·exp(iθ)

Identidicamos θ=-at o bien, θ=-Ωcosα·t. El plano de oscilación del péndulo de Foucault gira con velocidad angular constante Ω·cosα en el sentido de las agujas del reloj, tal como puede verse en la simulación al final de la página

Los coeficientes C y D de la solución de la ecuación diferencial, son números complejos cuya parte real e imaginaria se determinan a partir de las condiciones iniciales:

t=0{ x= x 0 dx dt = x ˙ 0 y= y 0 dy dt = y ˙ 0

El resultado es

C= 1 2 ( ( 1+ a b ) x 0 + y ˙ 0 b )+i( ( 1+ a b ) y 0 x ˙ 0 b ) D= 1 2 ( ( 1 a b ) x 0 y ˙ 0 b )+i( ( 1 a b ) y 0 + x ˙ 0 b )

Las expresiones de x e y en función del tiempo son:

x= C r cos(a+b)t+ D r cos(a+b)t C i sin(a+b)t+ D i sin(a+b)t y= C i cos(a+b)t+ D i cos(a+b)t+ C r sin(a+b)t D r sin(a+b)t

Donde Cr es la parte real del coeficiente C, Ci es la parte imaginaria y lo mismo cabe decir del coeficiente D

alfa=pi/6; %posición
w=1;  %sqrt(k/m)
W=w/10;  %velocidad rotación
a=W*cos(alfa);
b=sqrt(w^2+(W*cos(alfa))^2);
x0=1;   %posición inicial
y0=0;
vx0=0; %velocidad inicial
vy0=0;

%coeficientes
C=((1+a/b)*x0+vy0/b+1i*((1+a/b)*y0-vx0/b))/2;
D=((1-a/b)*x0-vy0/b+1i*((1-a/b)*y0+vx0/b))/2;
x=@(t) real(C)*cos((b-a)*t)+real(D)*cos((a+b)*t)-
imag(C)*sin((b-a)*t)+imag(D)*sin((a+b)*t);
y=@(t) imag(C)*cos((b-a)*t)+imag(D)*cos((b+a)*t)+
real(C)*sin((b-a)*t)-real(D)*sin((b+a)*t);
fplot(x,y,[0,50])
grid on
axis equal
xlabel('x')
ylabel('y')
title('Muelle en el plano local')

Comparando las dos figuras, vemos que son similares. En la primera hemos incluido términos proporcionales a Ω2 y en la segunda los hemos despreciado. En el caso del péndulo de Foaucault de longitud 67 m, la frecuencia ω de oscilación del péndulo y la velocidad angular de rotación Ω de la Tierra son, respectivamente.

ω= g/l = 9.8/67 =0.382rad/s Ω= 2π 24·60·60 =7 .27·10 5 rad/s

Lo que justifica la eliminación de los términos proporcionales a Ω2 en el sistema de dos ecuaciones diferenciales

En este programa, la trayectoria de la partícula unida a un muelle elástico en el Sistema de Referencia Local, se calcula resolviendo el sistema de dos ecuaciones diferenciales (incluyendo los términos proporcionales a Ω2) por el procedimiento de Runge-Kutta.

Actividades

Especificamos los parámetros del sistema

Especificamos las condiciones iniciales

Si la trayectoria aparece grande o se ve pequeña, se puede reducir o aumentar la escala para ajustarla a la ventana gráfica, en el control titulado Escala