]> El péndulo de Wilberforce
Anterior

El péndulo de Wilberforce

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 kθ a constante en las oscilaciones torsionales. Sea x el desplazamiento vertical del muelle de la posición de equilibrio, y θ  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 εxθ /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

E= 1 2 m ( dx dt ) 2 + 1 2 I ( dθ dt ) 2 + 1 2 k x x 2 + 1 2 k θ θ 2 + 1 2 εxθ

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 θ ˙ = dθ dt x ˙ = dx dt se escribe

L= 1 2 m x ˙ 2 + 1 2 I θ ˙ 2 1 2 k x x 2 1 2 k θ θ 2 1 2 εxθ d dt L x ˙ L x =0 d dt L θ ˙ L θ =0

Las ecuaciones del movimiento son:

d 2 x d t 2 + k x m x+ ε 2m θ=0 d 2 θ d t 2 + k θ I θ+ ε 2I x=0

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

ω x 2 = k x m ω θ 2 = k θ I

Eliminamos x en el sistema de dos ecuaciones diferenciales

d 4 θ d t 4 + ω θ 2 ( d 2 θ d t 2 )+ ε 2I ( d 2 x d t 2 )= d 4 θ d t 4 + ω θ 2 ( d 2 θ d t 2 )+ ε 2I ( ω x 2 x ε 2m θ )=0 d 4 θ d t 4 +( ω θ 2 + ω x 2 )( d 2 θ d t 2 )+( ω θ 2 · ω x 2 ε 2 4mI )θ=0

Suponiendo una solución de la forma

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

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

ω 4 ω 2 ( ω x 2 + ω θ 2 )+ ω x 2 · ω θ 2 ε 2 4mI =0

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

ω 1 2 = 1 2 ( ω x 2 + ω θ 2 + ( ω x 2 + ω θ 2 ) 2 4( ω x 2 · ω θ 2 ε 2 4mI ) )= 1 2 ( ω x 2 + ω θ 2 + ( ω x 2 ω θ 2 ) 2 + ε 2 mI ) ω 2 2 = 1 2 ( ω x 2 + ω θ 2 ( ω x 2 ω θ 2 ) 2 + ε 2 mI )

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)=Asin(ω1t)+ Bcos(ω1t)+ Csin(ω2t)+ Dcos(ω2t)

La velocidad angular de rotación es

dθ dt = ω 1 ( Acos( ω 1 t)Bsin( ω 1 t) )+ ω 2 ( Ccos( ω 2 t)Dsin( ω 2 t) )

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.

d 2 θ d t 2 + k θ I θ+ εx 2I =0 x(t)= 2I ε { ( ω 1 2 ω θ 2 )( Asin( ω 1 t)+Bcos( ω 1 t) )+ ( ω 2 2 ω θ 2 )( Csin( ω 2 t)+Dcos( ω 2 t) )

La velocidad del c.m. del cilindro es

dx dt = 2I ε { ( ω 1 2 ω θ 2 ) ω 1 ( Acos( ω 1 t)Bsin( ω 1 t) )+ ( ω 2 2 ω θ 2 ) ω 2 ( Ccos( ω 2 t)Dsin( ω 2 t) )

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

x 0 = 2I ε { ( ω 1 2 ω θ 2 )B+( ω 2 2 ω θ 2 )D } 0= 2I ε { ( ω 1 2 ω θ 2 ) ω 1 A+( ω 2 2 ω θ 2 ) ω 2 C } θ 0 =B+D 0= ω 1 A+ ω 2 C

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

A=C=0

B= ε x 0 2I ( ω 2 2 ω θ 2 ) θ 0 ( ω 2 2 ω 1 2 ) D= ε x 0 2I ( ω 1 2 ω θ 2 ) θ 0 ( ω 2 2 ω 1 2 )

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

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

dθ dt = ω 1 Bsin( ω 1 t) ω 2 Dsin( ω 2 t) x(t)= 2I ε { ( ω 1 2 ω θ 2 )Bcos( ω 1 t)+( ω 2 2 ω θ 2 )Dcos( ω 2 t) } dx dt = 2I ε { ( ω 1 2 ω θ 2 ) ω 1 Bsin( ω 1 t)+( ω 2 2 ω θ 2 ) ω 2 Dsin( ω 2 t) }

Modos normales de vibración

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

x 0 = 2I ε ( ω 1 2 ω θ 2 ) θ 0

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

θ(t)= θ 0 cos( ω 1 t) x(t)= 2I ε ( ω 1 2 ω θ 2 ) θ 0 cos( ω 1 t)

x 0 = 2I ε ( ω 2 2 ω θ 2 ) θ 0

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

θ(t)= θ 0 cos( ω 2 t) x(t)= 2 ε ( ω 1 2 ω θ 2 ) θ 0 cos( ω 2 t)

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 ωx y la frecuencia angular de las oscilaciones torsionales ωθ sean iguales ω =ωxθ

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

ω 1 2 = ω 2 + ε 2 4mI ω 2 2 = ω 2 ε 2 4mI   

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

θ(t)= θ 0 2 ( cos( ω 1 t)+cos( ω 2 t) ) x(t)= θ 0 2 I m ( cos( ω 1 t)cos( ω 2 t) )

x 0 =+ I m θ 0

x 0 = I m θ 0

Ejemplo:

 Frecuencia angular del movimiento de traslación es

ω x = 16.98 0.43 =6.284rad/s

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 1.28·10-4 kg·m2 más el momento de inercia de las dos esferas consideradas como masas puntuales de 3.4·10-3 kg cada una.

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

La frecuencia angular ωθ 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.

ω θ = 5.74· 10 3 I

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 ωθ=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 ωx, y en la parte superior ωθ.

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

16.98 0.43 = 5.74· 10 3 1.28· 10 4 +2·3.4· 10 3 d 2 d=0.051m

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

La frecuencias angulares ωx=ωθ=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ías correspondientes a las oscilaciones torsionales

Energía de acoplamiento 1 2 εxθ

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

Se pulsa el botón titulado Empieza

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

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

stokesApplet aparecerá en un explorador compatible con JDK 1.1.

Arrastre 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.

Anterior