La masa de la Tierra es M y su radio es R. En el espacio hay una varilla de masa despreciable de longitud 2l cuyo centro dista d del centro de la Tierra. En los extremos de la varilla hay dos masas m puntuales iguales. El sistema puede oscilar libremente en el plano de la figura alrededor del centro de masas de la varilla.

Las masas puntuales distan del centro de la Tierra r₁ y r₂ respectivamente

${\displaystyle \begin{array}{l}{r}_1^2={\left(R+l·\cos \theta \right)}^2+{\left(l·\text{sin}\theta \right)}^2=R^2+l^2+2lR\cos \theta \\ r_2^2={\left(R-l·\cos \theta \right)}^2+{\left(l·\text{sin}\theta \right)}^2=R^2+l^2-2lR\cos \theta \end{array}}$

Las fuerzas que ejerce la Tierra sobre cada una de las masas puntuales tienen por módulo

${\displaystyle F_1=G\frac{Mm}{r_1^2}{F}_2=G\frac{Mm}{r_2^2}}$

la dirección es a lo largo de la recta que une el centro de la Tierra y cada una de las masas puntuales, y de sentido hacia el centro de la Tierra.

El momento de las fuerzas $\overrightarrow{F_1}$ y $\overrightarrow{F_2}$ respecto del centro de la varilla es

M=F₁·l·sin(θ-α₁)-F₂·l·sin(θ+α₂)

Aplicando el teorema del seno

${\displaystyle \frac{r_1}{\text{sin(180-}\theta \text{)}}=\frac{d}{\text{sin}(\theta -{\alpha }_1)}\frac{r_2}{\text{sin}\theta }=\frac{d}{\text{sin}(180-\theta +{\alpha }_2)}}$

La ecuación del movimiento de rotación alrededor del centro de la varilla es

${\displaystyle I\frac{d^2\theta }{d{t}^2}=G\frac{Mm}{r_1^2}l\frac{d\text{sin}\theta }{r_1}-\frac{Mm}{r_2^2}l\frac{d\text{sin}\theta }{r_2}}$

El momento de inercia de las dos masas puntuales, respecto del eje de rotación perpendicular a la varilla y que pasa por su c.m. es  I=2ml²

${\displaystyle \frac{d^2\theta }{d{t}^2}=\frac{GM·d·\text{sin}\theta }{2l}\left(\frac{1}{r_1^3}-\frac{1}{r_2^3}\right)}$

Se resuelve la ecuación diferencial del movimiento por procedimientos numéricos con las siguientes condiciones iniciales: en el instante t=0, θ=π/6, dθ/dt=0

Convertimos la ecuación diferencial en otra más manejable tomando como unidad el radio de la Tierra, R=6.37·10⁶ m

${\displaystyle \frac{d^2\theta }{d{t}^2}=\frac{GM}{R^3}\frac{\left(d/R\right)\sin \theta }{2\left(l/R\right)}\left(\frac{1}{{\left(r_1/R\right)}^3}-\frac{1}{{\left(r_2/R\right)}^3}\right)}$

Definimos la función pendulo_marea del sistema de dos ecuaciones diferenciales de primer orden, donde x(1) representará los sucesivos valores de la variable θ y x(2) representará a la derivada dθ/dt.

function y=pendulo_marea(t,x,p)
    %x(1) posición , x(2) velocidad
    r=p(1); lon=p(2); %parámetros
    r1=1.0+lon^2/r^2-2*lon*cos(x(1))/r;
    r2=1.0+lon^2/r^2+2*lon*cos(x(1))/r;
    y=[x(2); -(1.543e-6*sin(x(1))/(2*lon*r^2))*(1.0/r1^1.5-1.0/r2^1.5)];
end

Resolvemos la ecuación diferencial mediante ode45 de MATLAB

x0=[pi/6 0]; %valores iniciales [posición, velocidad]
tspan=[0 20000]; 
p=[2.0 0.5]; %parámetros [distancia al centro, longitud del péndulo]
%no pasamos nada [] en el parámetro options de ode45
[t,x]=ode45(@pendulo_marea,tspan,x0,[],p);
plot(t,x(:,1),'r')
grid on
xlabel('t')
ylabel('ángulo');
title('Péndulo marea')

El periodo es aproximadamente 8287 s

Aproximaciones

Supongamos que la longitud de la varilla 2l es pequeña en comparación con la distancia d entre el centro de la Tierra y el c.m. de la varilla.

${\displaystyle \frac{1}{r_1^3}=\frac{1}{d^3}{\left(1+\frac{l^2}{d^2}+2\frac{l}{d}\cos \theta \right)}^{-3/2}\approx \frac{1}{d^3}\left(1-3\frac{l}{d}\cos \theta \right)}$

despreciando los términos en l²/d² en adelante. Del mismo modo

${\displaystyle \frac{1}{r_2^3}=\frac{1}{d^3}{\left(1+\frac{l^2}{d^2}-2\frac{l}{d}\cos \theta \right)}^{-3/2}\approx \frac{1}{d^3}\left(1+3\frac{l}{d}\cos \theta \right)}$

La ecuación del movimiento se puede aproximar a

${\displaystyle \frac{d^2\theta }{d{t}^2}=-\frac{3GM·\text{sin}\theta ·\cos \theta }{d^3}}$

o bien

${\displaystyle \frac{d^2(2\theta )}{d{t}^2}+\frac{3GM·}{d^3}\text{sin(2}\theta )=0}$

Esta es una ecuación similar a la del péndulo simple, en vez del ángulo θ, que hace el hilo con la vertical aparece el ángulo 2θ.

Suponiendo que el péndulo oscila con una pequeña amplitud angular, aproximamos sin(2θ)≈(2θ). Obtenemos la ecuación diferencial de un MAS

${\displaystyle \frac{d^2(2\theta )}{d{t}^2}+\frac{3GM·}{d^3}\text{(2}\theta )=0}$

de periodo

${\displaystyle P=\frac{2\pi }{\omega }=2\pi \sqrt{\frac{d^3}{3GM}}}$

El periodo es independiente de la longitud l del péndulo, esto indica que las fuerzas de marea se manifiestan también en sistemas cuyo tamaño es pequeño.

Si el péndulo se coloca en la superficie de la Tierra d=6370 km, M=5.98·10²⁴ kg, G=6.67·10^-11 Nm²/kg²

El periodo P=2920.2 s=48.7 min

Referencias

Grǿn Ǿ. A tidal force pendulum. Am. J. Phys. 51 (5) May 1983, pp. 429-431