El péndulo de Wilberforce

prev.gif (1231 bytes)home.gif (1232 bytes)next.gif (1211 bytes)

Oscilaciones

Osciladores acoplados
Varilla que pende 
 de dos muelles
marca.gif (847 bytes)Péndulo de Wilberforce
El péndulo doble
Péndulo-muelle
De las oscilaciones
a las ondas
Péndulos no acoplados
de distinta longitud
Ecuaciones del movimiento

Modos normales de vibración

Cuando las frecuencias son iguales

Balance energético

java.gif (886 bytes)Actividades

Referencias

 
El péndulo de Wilberforce es un conocido dispositivo para verificar el principio de la conservación de la energía. Se puede además, mostrar las oscilaciones acopladas de los modos longitudinales y torsionales de un cuerpo de forma cilíndrica que cuelga de un muelle en forma de hélice.

En esta página, vamos a resolver las ecuaciones del movimiento: de traslación y de rotación alrededor del eje vertical, a calcular las frecuencias de los modos normales de oscilación, y a determinar las condiciones iniciales que hacen que el sistema describa un modo normal de oscilación.

 

 

Ecuaciones del movimiento

Sea kx la constante elástica del muelle en las oscilaciones longitudinales y kq a constante en las oscilaciones torsionales. Sea x el desplazamiento vertical del muelle de la posición de equilibrio, y q  el ángulo de rotación alrededor del eje vertical.

El acoplamiento entre los dos modos de oscilación está descrito por una función lineal de la forma εxq /2, donde ε se denomina constante de acoplamiento.

La energía del sistema es la suma de la energía cinética de traslación del bloque, de rotación alrededor del eje vertical, la energía potencial del muelle cuando se deforma una longitud x, y cuando gira un ángulo θ, y la energía de acoplamiento

Las ecuaciones del movimiento de Lagrange nos llevan al sistema de dos ecuaciones diferenciales de segundo orden. La lagrangiana L=Ek-Ep, con los símbolos se escribe

Las ecuaciones del movimiento son:

En ausencia del término de acoplamiento ε=0 las ecuaciones diferenciales describen dos Movimientos Armónicos Simples de frecuencias angulares

Eliminamos x en el sistema de dos ecuaciones diferenciales

Suponiendo una solución de la forma

θ=A·sen(ωt)+ B·cos(ωt)

e insertándola en la ecuación diferencial de cuarto orden en θ, obtenemos la ecuación bicuadrada

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

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

θ(t)=Asen(ω1t)+ Bcos(ω1t)+ Csen(ω2t)+ Dcos(ω2t)

La velocidad angular de rotación es

Para calcular la posición x(t), se introduce θ(t) y su derivada segunda en la ecuación diferencial de segundo orden que describe la rotación del cilindro.

La velocidad del c.m. del cilindro es

Las condiciones iniciales determinan los valores de los coeficientes A, B, C, D.

Desplazamos el cilindro x0 y lo giramos un ángulo θ0 y a continuación lo soltamos. La velocidad inicial es  dx/dt=0, dθ/dt=0 en el instante t=0.

Tenemos que resolver el sistema de cuatro ecuaciones

Despejamos las incógnitas A, B, C, y D

A=C=0

Las expresiones de la posición x y θ y velocidad dx/dt y dθ/dt en función del tiempo son

θ(t)= Bcos(ω1t)+ Dsen(ω2t)

 

Modos normales de vibración

  • El modo 1 de frecuencia angular ω1, se obtiene cuando D=0,

La relación entre la posición inicial x0 del cilindro y el ángulo inicial de giro θ0 es

La posición x del cilindro y el ángulo girado θ en función del tiempo t son

θ(t)= θ0cos(ω1t)

  • El modo 2 de frecuencia angular ω2, se obtiene cuando B=0,

La posición x del cilindro y el ángulo girado θ en función del tiempo t son

θ(t)= θ0cos(ω2t)

 

Cuando las frecuencias son iguales

Si modificamos el momento de inercia I, cambiando la distancia d de las esferas al eje de rotación podemos conseguir que la frecuencia angular de las oscilaciones longitudinales wx y la frecuencia angular de las oscilaciones torsionales wq sean iguales w =wx=wq

Las frecuencias ω1 y ω2 de los modos normales de vibración tienen una expresión mucho más simple

Es importante analizar el caso de que x0=0. Se gira el cilindro un ángulo inicial θ0 y luego, se suelta.

Los coeficientes B=D=θ0/2. Las ecuaciones del movimiento son

  • La relación entre la posición inicial x0 del cilindro y el ángulo inicial de giro θ0 es en el modo normal 1 de frecuencia angular w1 es

  • La relación entre la posición inicial x0 del cilindro y el ángulo inicial de giro θ0 es en el modo normal 2 de frecuencia angular w2 es

 

Ejemplo:

  • Constante de acoplamiento ε=0.03 kg·m/s2
  • Masa del cilindro y las esferas m=0.43 kg

  • Constante elástica kx =16.98 N/m

 Frecuencia angular del movimiento de traslación es

  • Constante elástica kq=5.74·10-3 Nm

  • Momento de inercia I es variable

El momento de inercia I respecto del eje de rotación se puede variar, moviendo dos pequeñas esferas de igual masa m situadas a la misma distancia d del eje de rotación.

El momento de inercia del sistema, será igual al momento de inercia del cuerpo cilíndrico más el momento de inercia de las dos esferas consideradas como masas puntuales.

I=1.28·10-4+2·3.4·10-3·d2 kg m2

La frecuencia angular wq de las oscilaciones torsionales se puede cambiar modificando el momento de inercia I es decir, la distancia d de las pequeñas esferas, al eje de rotación.

Para d=3.5 cm=0.035 m, el momento de inercia I=1.3633·10-4 kg m2, la frecuencia angular de las oscilaciones torsionales es wq=6.489 rad/s

Las frecuencias de los modos de los modos normales de oscilación 1 y 2 valen

ω1=6.569 rad/s
ω
2
=6.200 rad/s

Si la posición inicial del cilindro x0=-0.04, el ángulo que hemos de girar para observar el primer modo normal de oscilación es

θ0=-4.20 rad=-240.7º

para observar el segundo modo normal de oscilación es

θ0=1.20 rad=68.8º

Con el puntero del ratón podemos modificar la distancia d o el momento de inercia I hasta lograr que ambas frecuencias sean iguales. Los valores de dichas frecuencias aparecen en el applet: en la parte inferior la frecuencia wx, y en la parte superior wq.

El valor de d que hace que ambas frecuencias sean iguales es

El momento de inercia vale I=1.45·10-4 kg m2

La frecuencias angulares wx=wq=6.284 rad/s

Las frecuencias de los modos normales de oscilación 1 y 2 valen

ω1=6.433 rad/s
ω
2
=6.131 rad/s

Si la posición inicial del cilindro x0=-0.04, el ángulo que hemos de girar para observar el primer modo normal de oscilación es

θ0=-2.18 rad=-124.6º

para observar el segundo modo normal de oscilación es

θ0=2.18 rad=124.6º

Es importante examinar el caso de que x0=0 y θ0≠0, por ejemplo, θ0=0.70 rad=40º

θ(t)=20(cos(ω1t)+ cos(ω2t))   grados
x
(t)=0.64(cos(ω1t)-cos(ω2t))  cm

 

Balance energético

Se puede observar los cambios energéticos en el diagrama en forma de tarta a la derecha del applet. En distintos colores se representan las energías potenciales elásticas y las energías cinéticas correspondientes a ambos modos de oscilación.

Energías correspondientes a las oscilaciones longitudinales

  • Energía potencial
  • Energía cinética

Energías correspondientes a las oscilaciones torsionales

  • Energía potencial
  • Energía cinética

Energía de acoplamiento

En la parte superior derecha del applet, se proporciona el valor numérico de la energía total que se mantiene constante durante el movimiento del péndulo.

 

Actividades

Se pulsa el botón titulado Inicio

Se arrastra con el puntero del ratón el círculo de color rojo situado a la izquierda del cilindro que cuelga del muelle, hacia el centro para disminuir el momento de inercia, hacia fuera para aumentarlo.

Se introduce

  • La posición inicial vertical x0 del cilindro en cm, en el control de edición titulado Posición (cm)
  • El ángulo de giro inicial θ0 en grados, en el control de edición titulado Angulo de giro

Se pulsa el botón titulado Empieza

Se observa el movimiento del péndulo, o la representación gráfica de:

  • El desplazamiento x en función del tiempo, activando la casilla traslación y pulsando el botón titulado Gráfica
  • El desplazamiento angular q en función del tiempo, activando la casilla rotación y pulsando el botón titulado Gráfica

Para comenzar una nueva "experiencia", se pulsa el botón titulado Inicio

stokesApplet aparecerá en un explorador compatible con JDK 1.1.

Arrastrar con el puntero del ratón el pequeño círculo de color rojo.

 

Referencias

Berg R. H, Marshall T. Wilberforce pendulum oscillations and normal modes. Am. J. Phys. 59 (1) January 1991, pp. 32-38.