El péndulo de Foucault

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

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

Ecuaciones del movimiento

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
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;
t=0:0.01:50;
x=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=imag(C)*cos((b-a)*t)+imag(D)*cos((b+a)*t)+real(C)*sin((b-a)*t)
-real(D)*sin((b+a)*t);
plot(x,y)
grid on
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