Oscilaciones

Osciladores acoplados

Dos acoplados

Tres acoplados

Cadena monoatómica

Modos normales
de vibración

Cadena diatómica

Varilla que pende 
 de dos muelles

Péndulo de Wilberforce

El péndulo doble

Péndulo-muelle

De las oscilaciones
a las ondas

Péndulos no acoplados
de distinta longitud

Modos normales de vibración

Cuando los péndulos son iguales

Actividades

Referencias

En esta página, se estudia un sistema de dos osciladores acoplados, el péndulo doble. Como en los ejemplos, de las páginas previas “Dos osciladores acoplados” y “El péndulo de Wilberforce” vamos a resolver las ecuaciones del movimiento, a calcular las frecuencias de los modos normales de oscilación, y las condiciones iniciales que hacen que el sistema describa un modo normal de oscilación.

El péndulo doble como se muestra en la figura, está formado por dos péndulos simples de longitudes l1 y l2, de los que cuelgan partículas de masas m1 y m2.  En un instante determinado t, los hilos inextensibles forman ángulos θ1 y θ2 con la vertical.

Ecuaciones del movimiento

La energía del sistema es la suma de la energía potencial y de la energía cinética de las dos partículas.

Situamos el nivel cero de energía potencial en el punto de suspensión del primer péndulo. La energía potencial es

E_p=-m₁·gl₁cos θ₁-m₂g(l₁cos θ₁ +l₂cos θ₂)

El módulo de la velocidad del primer péndulo es v₁=l₁·dθ₁/dt (velocidad angular por el radio de la circunferencia que describe), sus componentes son

El módulo de la velocidad de la segunda partícula en un sistema de referencia que se mueve con la velocidad v₁ de la primera partícula es v₂=l₂·dθ₂/dt. Sus componentes son

La velocidad de la segunda partícula respecto al sistema de referencia inercial situado en O es la suma vectorial de ambas velocidades

Calculamos los módulos de las velocidades de las dos partículas. La energía cinética del sistema es

Si nos restringimos a pequeños valores de los ángulos θ₁ y θ₂, las ecuaciones del movimiento se hacen mucho más simples. Desarrollamos en serie de Taylor cosθ y tomamos los dos primeros términos no nulos del desarrollo en serie

La energía potencial con esta aproximación, se expresa

La energía cinética con esta aproximación, se expresa

Las ecuaciones del movimiento de Lagrange nos llevan al sistema de dos ecuaciones diferenciales de segundo orden. La lagrangiana L=E_k-E_p con el símbolo

Las ecuaciones del movimiento son

Calculamos las derivadas segundas respecto del tiempo de estas dos ecuaciones diferenciales

Entre estas cuatro ecuaciones, eliminamos θ₂, obteniendo la ecuación diferencial de cuarto orden en θ₁.

Suponiendo una solución de la forma

θ₁=Asen(ωt)+Bcos(ωt) o bien, θ₁=Asen(ωt+φ)

Calculamos la derivada segunda y la derivada cuarta de θ₁ y las introducimos en la ecuación diferencial de cuarto orden, obteniendo la siguiente ecuación bicuadrada en ω.

Las dos raíces reales de esta ecuación son las frecuencias ω₁ y ω₂ de los modos normales de vibración

La forma general del ángulo θ₁ en función del tiempo t es una combinación de los dos modos normales de vibración

θ₁(t)=A₁sen(ω₁t)+ B₁cos(ω₁t)+ C₁sen(ω₂t)+ D₁cos(ω₂t)

Lo mismo para θ₂

θ₂(t)=A₂sen(ω₁t)+ B₂cos(ω₁t)+ C₂sen(ω₂t)+ D₂cos(ω₂t)

Las velocidades angulares son

Relacionamos los coeficientes A₁ y A₂, B₁ y B₂, C₁ y C₂, D₁ y D₂, calculando la derivada segunda de θ₂, la derivada segunda de θ₁ e introduciéndolos en una de las dos ecuaciones diferenciales del movimiento

Tenemos cuatro relaciones

Las condiciones iniciales determinan los valores de los cuatro coeficientes A₁, B₁, C₁ y D₁ y por tanto, de los coeficientes A₂, B₂, C₂ y D₂, a través de las relaciones establecidas anteriormente.

Desplazamos el primer péndulo un ángulo θ₁₀ con respecto a la vertical, y el segundo péndulo un ángulo θ₂₀ a continuación los soltamos. La velocidad inicial es  dθ₁/dt=0, dθ₂/dt=0 en el instante t=0.

Tenemos que resolver el sistema de cuatro ecuaciones junto a las relaciones establecidas anteriormente.

θ₁₀=B₁+D₁
0=ω₁·A₁+ ω₂·C₁
θ₂₀=B₂+D₂
0=ω₁·A₂+ ω₂·C₂

La solución de este sistema es

A₁=C₁=A₂=C₂=0

La posición angular en función del tiempo t de cada una de las partículas es

Modos normales de vibración

• El modo 1 de frecuencia angular ω₁, se obtiene cuando D₂=0, la relación entre los ángulos iniciales de desviación de los dos péndulos es

B₁=θ₁₀, B₂=θ₂₀

La posición angular en función del tiempo t de cada una de las partículas es

• El modo 2 de frecuencia angular ω₂, se obtiene cuando B₂=0, la relación entre los ángulos iniciales de desviación de los dos péndulos es

D₂=θ₂₀, y D₁=θ₁₀

La posición angular en función del tiempo t de cada una de las partículas es

Cuando los péndulos son iguales

Cuando m₁=m₂ y l₁=l₂ las frecuencias ω₁ y ω₂ de los modos normales de vibración tienen una expresión mucho más simple

Los coeficientes B₂, B₁, D₂ y D₁ valen

Es importante analizar el caso de que θ₂₀=0. Se desplaza el primer péndulo un ángulo θ₁₀ respecto de la vertical y se suelta

 Las ecuaciones del movimiento de cada unos de las partículas se escriben

En la figura, se representa en color rojo θ₁(t) y en color azul la amplitud modulada

En la figura, se representa en color rojo θ₂(t) y en color azul la amplitud modulada

• El modo 1 de frecuencia angular ω₁, se obtiene cuando D₂=0, la relación entre los ángulos iniciales de desviación de los dos péndulos es

B₂=θ₂₀, y B₁=θ₁₀

La posición angular en función del tiempo t de cada una de las partículas es

 

• El modo 2 de frecuencia angular ω₂, se obtiene cuando B₂=0, la relación entre los ángulos iniciales de desviación de los dos péndulos es

D₂=θ₂₀, y D₁=θ₁₀

La posición angular en función del tiempo t de cada una de las partículas es

Actividades

Se introduce

• La masa m₂ de la partícula del péndulo inferior en kg, en el control de edición titulado Masa 2.

• La masa de la partícula del péndulo superior se ha fijado en m₁=1 kg

• La longitud l₂ del péndulo inferior en m, en el control de edición titulado Longitud 2.

• La longitud del péndulo superior se ha fijado en l₁=1 m

Se pulsa el botón titulado Inicio

• El ángulo de desviación inicial θ₁₀ del péndulo superior en grados, actuando en la barra de desplazamiento titulada Angulo 1.

• El ángulo de desviación inicial θ₂₀ del péndulo inferior en grados, actuando en la barra de desplazamiento titulada Angulo 2.

Se pulsa el botón titulado Empieza

Se observa el movimiento de los dos péndulos, y la representación gráfica del ángulo de desviación de cada péndulo θ₁ y θ₂ en función del tiempo t.

En la parte superior derecha del applet, se muestra el valor de la energía total que permanece constante durante el movimiento de los péndulos.