Péndulo doble


Introducción

Un péndulo doble es uno de los sistemas físicos más simples que exhibe caos. Consiste en un primer brazo de péndulo colgando de un pivote fijo, con un segundo brazo de péndulo colgando del extremo del primero. Bajo la gravedad, este arreglo produce un movimiento extraordinariamente sensible a sus condiciones iniciales — dos péndulos dobles liberados desde posiciones que difieren incluso en una fracción de grado trazarán, en cuestión de segundos, trayectorias completamente distintas. Este comportamiento no es resultado de aleatoriedad ni de ruido; es una consecuencia determinista de las ecuaciones no lineales que rigen el sistema. El péndulo doble se sitúa en la encrucijada de la mecánica clásica, la teoría de oscilaciones y la teoría del caos, lo que lo convierte en uno de los sistemas más instructivos y visualmente impactantes de toda la física.


La física explicada

Un péndulo simple se balancea de un lado a otro en un patrón regular y repetitivo. Para ángulos pequeños el movimiento es casi sinusoidal, y el período depende solo de la longitud del brazo y de la intensidad de la gravedad — no de la masa ni de la amplitud. Por eso los péndulos se usaron en relojes durante siglos. El péndulo simple es un oscilador bien comportado y predecible.

El péndulo doble rompe esa predictibilidad. Como el segundo brazo está unido a un punto en movimiento en lugar de a un pivote fijo, los dos brazos interactúan entre sí a través de sus fuerzas mutuas. El movimiento del brazo superior cambia la gravedad efectiva que experimenta el inferior, y el balanceo del brazo inferior tira del superior a través de la tensión en la varilla de unión. Estas interacciones son no lineales — las fuerzas dependen de los ángulos de una manera que no puede simplificarse en una relación proporcional limpia — y es esta no linealidad la que produce el caos.

El marco correcto para analizar el péndulo doble es la mecánica lagrangiana, una reformulación de la mecánica newtoniana que trabaja en términos de energía en lugar de fuerzas. En vez de resolver todas las tensiones y reacciones en los puntos pivote, escribimos la energía cinética total y la energía potencial total del sistema, y luego aplicamos un procedimiento matemático llamado ecuaciones de Euler-Lagrange para derivar las ecuaciones de movimiento. El resultado es un par de ecuaciones diferenciales acopladas, no lineales, de segundo orden — una para cada ángulo. Estas ecuaciones no tienen solución analítica cerrada; para encontrar cómo cambian los ángulos en el tiempo, deben integrarse numéricamente.

El caos en el péndulo doble significa que pequeñas diferencias en las condiciones iniciales crecen exponencialmente con el tiempo, una propiedad cuantificada por el exponente de Lyapunov. Si dos péndulos dobles parten con ángulos que difieren en una cantidad pequeña, la diferencia entre sus trayectorias crece aproximadamente como e elevado al exponente de Lyapunov por el tiempo. Esta divergencia exponencial es la firma matemática del caos, y pone un límite absoluto a qué tan lejos en el futuro puede predecirse el movimiento en la práctica. A pesar de esto, el sistema conserva la energía perfectamente — el caos aquí no significa aleatoriedad en la energía, solo en la trayectoria.

Para oscilaciones pequeñas alrededor de la posición de equilibrio inferior, los términos no lineales son despreciables y el sistema se comporta como dos osciladores lineales acoplados. En este régimen hay dos modos normales: uno en el que ambos brazos se balancean en la misma dirección a la misma frecuencia, y otro en el que los brazos se balancean en direcciones opuestas a una frecuencia más alta. A medida que la amplitud crece, los términos no lineales entran en juego y el movimiento transita de la oscilación regular al caos.


Ecuaciones clave

Ángulo del brazo superior desde la verticalθ₁ — medido en radianes desde la vertical inferior
Ángulo del brazo inferior desde la verticalθ₂ — medido en radianes desde la vertical inferior
Energía cinética totalKE = ½·(m₁ + m₂)·L₁²·θ₁'² + ½·m₂·L₂²·θ₂'² + m₂·L₁·L₂·θ₁'·θ₂'·cos(θ₁ − θ₂)
Energía potencial totalPE = −(m₁ + m₂)·g·L₁·cos(θ₁) − m₂·g·L₂·cos(θ₂)
Ecuación de movimiento del brazo superior (de Euler-Lagrange)(m₁ + m₂)·L₁·θ₁'' + m₂·L₂·θ₂''·cos(θ₁ − θ₂) + m₂·L₂·θ₂'²·sin(θ₁ − θ₂) + (m₁ + m₂)·g·sin(θ₁) = 0
Ecuación de movimiento del brazo inferior (de Euler-Lagrange)m₂·L₂·θ₂'' + m₂·L₁·θ₁''·cos(θ₁ − θ₂) − m₂·L₁·θ₁'²·sin(θ₁ − θ₂) + m₂·g·sin(θ₂) = 0
Conservación de la energía mecánica totalE = KE + PE = constante (sin fricción ni resistencia del aire)
Divergencia exponencial de trayectorias cercanas (caos)|δ(t)| ≈ |δ(0)|·exp(λ·t), donde λ es el exponente de Lyapunov

Variables clave

SímboloNombreUnidadSignificado
θ₁Ángulo del brazo superiorradÁngulo del brazo superior medido desde la vertical inferior
θ₂Ángulo del brazo inferiorradÁngulo del brazo inferior medido desde la vertical inferior
θ₁'Velocidad angular superiorrad/sTasa de cambio del ángulo del brazo superior
θ₂'Velocidad angular inferiorrad/sTasa de cambio del ángulo del brazo inferior
θ₁''Aceleración angular superiorrad/s²Tasa de cambio de la velocidad angular superior
θ₂''Aceleración angular inferiorrad/s²Tasa de cambio de la velocidad angular inferior
m₁Masa de la lenteja superiorkgMasa de la lenteja al final del brazo superior
m₂Masa de la lenteja inferiorkgMasa de la lenteja al final del brazo inferior
L₁Longitud del brazo superiormLongitud del brazo superior
L₂Longitud del brazo inferiormLongitud del brazo inferior
gAceleración gravitatoriam/s²9,81 m/s² en la superficie de la Tierra
KEEnergía cinéticaJEnergía cinética total de ambos brazos y lentejas
PEEnergía potencialJEnergía potencial gravitatoria total de ambas lentejas
EEnergía totalJEnergía mecánica total; se conserva a lo largo del movimiento
λExponente de Lyapunov1/sMide la tasa de divergencia de las trayectorias cercanas

Ejemplos del mundo real


Cómo funciona la simulación

La simulación muestra dos brazos de péndulo conectados, con masas en las puntas. Cuatro deslizadores te permiten fijar la longitud de cada brazo y los ángulos iniciales de cada uno (medidos en grados desde la vertical inferior). Pulsa Iniciar para liberar el sistema desde el reposo en los ángulos elegidos. La simulación integra las ecuaciones de Euler-Lagrange completas y no lineales en cada paso de tiempo usando un integrador numérico estable, así que la trayectoria es precisa incluso en el régimen caótico de gran amplitud.

Un trazo de la posición de la lenteja inferior se dibuja en el lienzo, revelando los hermosos patrones que el péndulo doble puede generar. Pulsa Reiniciar y vuelve a lanzar después de cambiar uno de los ángulos por solo un grado para ver de forma directa cómo dos trayectorias casi idénticas divergen rápidamente — la firma del caos. La energía total se mantiene constante a lo largo de toda la integración, una verificación de que el método numérico es estable y de que el sistema es realmente conservativo.


Lecturas adicionales