En general, un péndulo doble o doble péndulo es un sistema compuesto por dos péndulos, con el segundo colgando del extremo del primero. En el caso más simple, se trata de dos péndulos simples, con el inferior colgando de la masa pendular del superior.

Normalmente se sobreentiende que nos referimos a un péndulo doble plano, con dos péndulos planos coplanarios. Este sistema físico posee dos grados de libertad y exhibe un rico comportamiento dinámico. Su movimiento está gobernado por dos ecuaciones diferenciales ordinarias acopladas. Por encima de cierta energía, su movimiento es caótico.

Análisis del movimiento del péndulo doble plano

Cinemática

En la cinemática solo estamos interesados en encontrar las expresiones de la posición, la velocidad, la aceleración y en términos de las variables que especifican el estado del péndulo doble, sin interesarnos por las fuerzas actuantes. Nos serviremos de las siguientes coordenadas:

  • x,y = posición horizontal y vertical de la masa de un péndulo
  • θ = ángulo de un péndulo respecto a la vertical (0 = vertical hacia abajo, antihorario es positivo)
  • l = longitud de la varilla (constante)

Asociaremos al péndulo superior el subíndice 1, y al inferior el subíndice 2. Pondremos el origen de coordenadas en el punto de pivote del péndulo superior. El sentido de las ordenadas crecientes se toma hacia arriba.

A partir de consideraciones trigonométricas escribimos las expresiones de las posiciones x1, y1, x2, y2 en términos de los ángulos θ1, θ2:

x 1 = l 1 sen θ 1 {\displaystyle x_{1}=l_{1}\operatorname {sen} \theta _{1}\,}
y 1 = l 1 cos θ 1 {\displaystyle y_{1}=-l_{1}\cos \theta _{1}\,}
x 2 = x 1 l 2 sen θ 2 {\displaystyle x_{2}=x_{1} l_{2}\operatorname {sen} \theta _{2}\,}
y 2 = y 1 l 2 cos θ 2 {\displaystyle y_{2}=y_{1}-l_{2}\cos \theta _{2}\,}

Derivando con respecto al tiempo obtenemos:

x ˙ 1 = θ ˙ 1 l 1 cos θ 1 {\displaystyle {\dot {x}}_{1}={\dot {\theta }}_{1}l_{1}\cos \theta _{1}}
y ˙ 1 = θ ˙ 1 l 1 sen θ 1 {\displaystyle {\dot {y}}_{1}={\dot {\theta }}_{1}l_{1}\operatorname {sen} \theta _{1}}
x ˙ 2 = x ˙ 1 θ ˙ 2 l 2 cos θ 2 {\displaystyle {\dot {x}}_{2}={\dot {x}}_{1} {\dot {\theta }}_{2}l_{2}\cos \theta _{2}}
y ˙ 2 = y ˙ 1 θ ˙ 2 l 2 sen θ 2 {\displaystyle {\dot {y}}_{2}={\dot {y}}_{1} {\dot {\theta }}_{2}l_{2}\operatorname {sen} \theta _{2}}

Y derivando una segunda vez:

x ¨ 1 = θ ˙ 1 2 l 1 sen θ 1 θ ¨ 1 l 1 cos θ 1 {\displaystyle {\ddot {x}}_{1}=-{\dot {\theta }}_{1}^{2}l_{1}\operatorname {sen} \theta _{1} {\ddot {\theta }}_{1}l_{1}\cos \theta _{1}}
y ¨ 1 = θ ˙ 1 2 l 1 cos θ 1 θ ¨ 1 l 1 sen θ 1 {\displaystyle {\ddot {y}}_{1}={\dot {\theta }}_{1}^{2}l_{1}\cos \theta _{1} {\ddot {\theta }}_{1}l_{1}\operatorname {sen} \theta _{1}}
x ¨ 2 = x ¨ 1 θ ˙ 2 2 l 2 sen θ 2 θ ¨ 2 l 2 cos θ 2 {\displaystyle {\ddot {x}}_{2}={\ddot {x}}_{1}-{\dot {\theta }}_{2}^{2}l_{2}\operatorname {sen} \theta _{2} {\ddot {\theta }}_{2}l_{2}\cos \theta _{2}}
y ¨ 2 = y ¨ 1 θ ˙ 2 2 l 2 cos θ 2 θ ¨ 2 l 2 sen θ 2 {\displaystyle {\ddot {y}}_{2}={\ddot {y}}_{1} {\dot {\theta }}_{2}^{2}l_{2}\cos \theta _{2} {\ddot {\theta }}_{2}l_{2}\operatorname {sen} \theta _{2}}

Fuerzas

Definimos las variables:

Usaremos la ley de Newton F = m a {\displaystyle F=ma} , escribiendo por separado las ecuaciones de las componentes verticales y horizontales de las fuerzas.

Sobre la masa m 1 {\displaystyle m_{1}} actúan la tensión en la parte superior de la varilla T 1 {\displaystyle T_{1}} , la tensión en la parte inferior de la varilla T 2 {\displaystyle T_{2}} , y la gravedad -m1g:

m 1 x ¨ 1 = T 1 sen θ 1 T 2 sen θ 2 {\displaystyle m_{1}{\ddot {x}}_{1}=-T_{1}\operatorname {sen} \theta _{1} T_{2}\operatorname {sen} \theta _{2}}
m 1 y ¨ 1 = T 1 cos θ 1 T 2 cos θ 2 m 1 g {\displaystyle m_{1}{\ddot {y}}_{1}=T_{1}\cos \theta _{1}-T_{2}\cos \theta _{2}-m_{1}g}

Sobre la masa m 2 {\displaystyle m_{2}} , actúan la tensión T 2 {\displaystyle T_{2}} y la gravedad –m2g:

m 2 x ¨ 2 = T 2 sen θ 2 {\displaystyle m_{2}{\ddot {x}}_{2}=-T_{2}\operatorname {sen} \theta _{2}}
m 2 y ¨ 2 = T 2 cos θ 2 m 2 g {\displaystyle m_{2}{\ddot {y}}_{2}=T_{2}\cos \theta _{2}-m_{2}g}

Ecuaciones de movimiento

A partir de las ecuaciones anteriores, tras realizar numerosas operaciones algebraicas con la finalidad de encontrar las expresiones de θ 1 ¨ {\displaystyle {\ddot {\theta _{1}}}} , θ 2 ¨ {\displaystyle {\ddot {\theta _{2}}}} en términos de θ 1 {\displaystyle \theta _{1}\,} , θ 1 ˙ {\displaystyle {\dot {\theta _{1}}}\,} , θ 2 {\displaystyle \theta _{2}\,} , θ 2 ˙ {\displaystyle {\dot {\theta _{2}}}\,} , llegaríamos a las ecuaciones de movimiento para el péndulo doble:

θ ¨ 1 = g ( 2 m 1 m 2 ) sen θ 1 m 2 g sen ( θ 1 2 θ 2 ) 2 sen ( θ 1 θ 2 ) m 2 ( θ ˙ 2 2 l 2 θ ˙ 1 2 l 1 cos ( θ 1 θ 2 ) ) l 1 ( 2 m 1 m 2 m 2 cos ( 2 θ 1 2 θ 2 ) ) {\displaystyle {\ddot {\theta }}_{1}={\frac {-g(2m_{1} m_{2})\operatorname {sen} \theta _{1}-m_{2}g\operatorname {sen}(\theta _{1}-2\theta _{2})-2\operatorname {sen}(\theta _{1}-\theta _{2})m_{2}({\dot {\theta }}_{2}^{2}l_{2} {\dot {\theta }}_{1}^{2}l_{1}\cos(\theta _{1}-\theta _{2}))}{l_{1}(2m_{1} m_{2}-m_{2}\cos(2\theta _{1}-2\theta _{2}))}}}
θ ¨ 2 = 2 sen ( θ 1 θ 2 ) ( θ ˙ 1 2 l 1 ( m 1 m 2 ) g ( m 1 m 2 ) cos θ 1 θ ˙ 2 2 l 2 m 2 cos ( θ 1 θ 2 ) ) l 2 ( 2 m 1 m 2 m 2 cos ( 2 θ 1 2 θ 2 ) ) {\displaystyle {\ddot {\theta }}_{2}={\frac {2\operatorname {sen}(\theta _{1}-\theta _{2})({\dot {\theta }}_{1}^{2}l_{1}(m_{1} m_{2}) g(m_{1} m_{2})\cos \theta _{1} {\dot {\theta }}_{2}^{2}l_{2}m_{2}\cos(\theta _{1}-\theta _{2}))}{l_{2}(2m_{1} m_{2}-m_{2}\cos(2\theta _{1}-2\theta _{2}))}}}

Energía

La energía cinética viene expresada por:

T = 1 2 m 1 ( x ˙ 1 2 y ˙ 1 2 ) 1 2 m 2 ( x ˙ 2 2 y ˙ 2 2 ) = 1 2 m 1 l 1 2 θ ˙ 1 2 1 2 m 2 [ l 1 2 θ ˙ 1 2 l 2 2 θ ˙ 2 2 2 l 1 l 2 θ ˙ 1 θ ˙ 2 cos ( θ 1 θ 2 ) ] {\displaystyle T={\frac {1}{2}}m_{1}({\dot {x}}_{1}^{2} {\dot {y}}_{1}^{2}) {\frac {1}{2}}m_{2}({\dot {x}}_{2}^{2} {\dot {y}}_{2}^{2})={\frac {1}{2}}m_{1}l_{1}^{2}{\dot {\theta }}_{1}^{2} {\frac {1}{2}}m_{2}[l_{1}^{2}{\dot {\theta }}_{1}^{2} l_{2}^{2}{\dot {\theta }}_{2}^{2} 2l_{1}l_{2}{\dot {\theta }}_{1}{\dot {\theta }}_{2}\cos(\theta _{1}-\theta _{2})]}

La energía potencial:

V = m 1 g y 1 m 2 g y 2 = ( m 1 m 2 ) g l 1 cos θ 1 m 2 g l 2 cos θ 2 {\displaystyle V=m_{1}gy_{1} m_{2}gy_{2}=-(m_{1} m_{2})gl_{1}\cos \theta _{1}-m_{2}gl_{2}\cos \theta _{2}\,} .

Por tanto, el movimiento se regirá por la lagrangiana

L = T V = 1 2 ( m 1 m 2 ) l 1 2 θ ˙ 1 2 1 2 m 2 l 2 2 θ ˙ 2 2 m 2 l 1 l 2 θ ˙ 1 θ ˙ 2 cos ( θ 1 θ 2 ) ( m 1 m 2 ) g l 1 cos θ 1 m 2 g l 2 cos θ 2 {\displaystyle {\mathcal {L}}=T-V={\frac {1}{2}}(m_{1} m_{2})l_{1}^{2}{\dot {\theta }}_{1}^{2} {\frac {1}{2}}m_{2}l_{2}^{2}{\dot {\theta }}_{2}^{2} m_{2}l_{1}l_{2}{\dot {\theta }}_{1}{\dot {\theta }}_{2}\cos(\theta _{1}-\theta _{2}) (m_{1} m_{2})gl_{1}\cos \theta _{1} m_{2}gl_{2}\cos \theta _{2}}

Ecuaciones de movimiento de Lagrange

Usando las ecuaciones de Lagrange en este caso particular son:

d d t ( L θ ˙ 1 ) L θ 1 = 0 , d d t ( L θ ˙ 2 ) L θ 2 = 0 {\displaystyle {\frac {d}{dt}}\left({\frac {\partial {\mathcal {L}}}{\partial {\dot {\theta }}_{1}}}\right)-{\frac {\partial {\mathcal {L}}}{\partial \theta _{1}}}=0,\qquad {\frac {d}{dt}}\left({\frac {\partial {\mathcal {L}}}{\partial {\dot {\theta }}_{2}}}\right)-{\frac {\partial {\mathcal {L}}}{\partial \theta _{2}}}=0}

Calculando explícitamente las derivadas de la expresión anterior se llega a:

{ ( m 1 m 2 ) l 1 2 θ ¨ 1 m 2 θ ¨ 2 l 1 l 2 cos ( θ 1 θ 2 ) m 2 θ ˙ 2 l 1 l 2 ( θ ˙ 1 θ ˙ 2 ) sen ( θ 1 θ 2 ) m 2 θ ˙ 1 θ ˙ 2 l 1 l 2 sen ( θ 1 θ 2 ) ( m 1 m 2 ) g l 1 sen θ 1 = 0 m 2 l 2 2 θ ¨ 2 m 2 θ ¨ 1 l 1 l 2 cos ( θ 1 θ 2 ) m 2 θ ˙ 1 l 1 l 2 ( θ ˙ 1 θ ˙ 2 ) sen ( θ 1 θ 2 ) m 2 θ ˙ 1 θ ˙ 2 l 1 l 2 sen ( θ 1 θ 2 ) m 2 g l 2 sen θ 2 = 0 {\displaystyle {\begin{cases}(m_{1} m_{2})l_{1}^{2}{\ddot {\theta }}_{1} m_{2}{\ddot {\theta }}_{2}l_{1}l_{2}\cos(\theta _{1}-\theta _{2})-m_{2}{\dot {\theta }}_{2}l_{1}l_{2}({\dot {\theta }}_{1}-{\dot {\theta }}_{2})\operatorname {sen}(\theta _{1}-\theta _{2}) m_{2}{\dot {\theta }}_{1}{\dot {\theta }}_{2}l_{1}l_{2}\operatorname {sen}(\theta _{1}-\theta _{2}) (m_{1} m_{2})gl_{1}\operatorname {sen} \theta _{1}=0\\m_{2}l_{2}^{2}{\ddot {\theta }}_{2} m_{2}{\ddot {\theta }}_{1}l_{1}l_{2}\cos(\theta _{1}-\theta _{2})-m_{2}{\dot {\theta }}_{1}l_{1}l_{2}({\dot {\theta }}_{1}-{\dot {\theta }}_{2})\operatorname {sen}(\theta _{1}-\theta _{2})-m_{2}{\dot {\theta }}_{1}{\dot {\theta }}_{2}l_{1}l_{2}\operatorname {sen}(\theta _{1}-\theta _{2}) m_{2}gl_{2}\operatorname {sen} \theta _{2}=0\end{cases}}}

Simplificando obtenemos:

{ ( m 1 m 2 ) l 1 2 θ ¨ 1 m 2 θ ¨ 2 l 1 l 2 cos ( θ 1 θ 2 ) m 2 θ ˙ 2 2 l 1 l 2 sen ( θ 1 θ 2 ) ( m 1 m 2 ) g l 1 sen θ 1 = 0 m 2 l 2 2 θ ¨ 2 m 2 θ ¨ 1 l 1 l 2 cos ( θ 1 θ 2 ) m 2 θ ˙ 1 2 l 1 l 2 sen ( θ 1 θ 2 ) m 2 g l 2 sen θ 2 = 0 {\displaystyle {\begin{cases}(m_{1} m_{2})l_{1}^{2}{\ddot {\theta }}_{1} m_{2}{\ddot {\theta }}_{2}l_{1}l_{2}\cos(\theta _{1}-\theta _{2}) m_{2}{\dot {\theta }}_{2}^{2}l_{1}l_{2}\operatorname {sen}(\theta _{1}-\theta _{2}) (m_{1} m_{2})gl_{1}\operatorname {sen} \theta _{1}=0\\m_{2}l_{2}^{2}{\ddot {\theta }}_{2} m_{2}{\ddot {\theta }}_{1}l_{1}l_{2}\cos(\theta _{1}-\theta _{2})-m_{2}{\dot {\theta }}_{1}^{2}l_{1}l_{2}\operatorname {sen}(\theta _{1}-\theta _{2}) m_{2}gl_{2}\operatorname {sen} \theta _{2}=0\end{cases}}}

Estas son las ecuaciones de Lagrange para un péndulo doble en el que hemos escogido como coordenadas generalizadas las polares y en el que hay dos ligaduras( l 1 {\displaystyle l_{1}} y l 2 {\displaystyle l_{2}} constantes).

Véase también

  • Péndulo
  • Péndulo balístico
  • Péndulo cicloidal
  • Péndulo cónico
  • Péndulo de Foucault
  • Péndulo de Newton
  • Péndulo de Pohl
  • Péndulo de torsión
  • Péndulo esférico
  • Péndulo físico
  • Péndulo simple
  • Péndulo simple equivalente

Bibliografía

  • Marion, Jerry B. (1996). Dinámica clásica de las partículas y sistemas. Barcelona: Ed. Reverté. ISBN 84-291-4094-8

Enlaces externos

  • Weisstein, Eric W. «Double Pendulum». ScienceWorld (en inglés). Wolfram Research. 
  • Animación que muestra movimiento del péndulo doble y reparto de energía entre uno y otro péndulo Archivado el 22 de enero de 2008 en Wayback Machine.
  • Curso Interactivo de Física en Internet. Ángel Franco García.

Péndulo doble Físicadcyt

Opiniones de doble pendulo

El Péndulo Doble A PDF Péndulo Velocidad

El experimento del péndulo doble ¿Es posible predecir el azar? NeoTeo

Pendulo Doble PDF Física Mecanica clasica