Oscilaciones

Osciladores (II)

Oscilador no lineal (I)

Oscilador no lineal (II)

El péndulo

Oscilaciones
anarmónicas

Molécula diatómica

El oscilador Morse

Oscilaciones de un
 cilindro que rueda

Movimiento de
un pistón

Movimiento del cilindro con un imán en su interior

Actividades

Referencias

Consideremos un cilindro con un imán en su interior, situado sobre un plano inclinado en una región en la que existe un campo magnético uniforme. Supondremos que en todo momento el coeficiente de rozamiento estático μe es suficientemente grande para que se cumpla la condición que el cilindro rueda sin deslizar.

Movimiento del cilindro sobre el plano inclinado

En la página “Movimiento de rodar sobre un plano inclinado” hemos estudiado el movimiento de un cilindro baja rodando (deslizando y sin deslizar) por un plano inclinado.

Las fuerzas que actúan sobre el cilindro son: el peso mg la reacción del plano inclinado N la fuerza de rozamiento Fr en el punto de contacto entre la rueda y el plano.

Las ecuaciones del movimiento serán

• Movimiento de traslación del c.m.

mg·senq -F_r=ma_c

• Movimiento de rotación alrededor de un eje que pasa por el c.m.

F_rR=I_ca

• Relación entre el movimiento de traslación y rotación (rueda sin deslizar)

a_c=a R

Para un cilindro de momento de inercia I_c=mR²/2, la aceleración del centro de masa es

Movimiento del cilindro con un imán en su interior

Consideremos ahora un cilindro situado sobre un plano inclinado con un imán de momento magnético μ en su interior. En la región en la que se mueve el cilindro existe un campo magnético  uniforme B, que apunta verticalmente hacia arriba, producido por un para de bobinas de Helmholtz.

Supongamos que en el instante inicial, el cilindro se encuentra en el origen, el campo magnético B y el momento magnético del imán μ forman un ángulo π/2. El cilindro rueda sobre el plano inclinado, y en el instante t el ángulo que forma el momento magnético con la dirección de B es  φ+π/2. El campo magnético B ejerce un momento sobre el imán, dado por el producto vectorial M=μ×B

• El módulo del producto vectorial es M=μB·sen(φ+π/2)= μB·cosφ
• La dirección a lo largo del eje de rotación
• El sentido se determina aplicando la regla del sacacorchos

El momento M no es constante, sino que su módulo y su sentido cambian a medida que el cilindro rueda a lo largo del plano inclinado tal como se muestra en la figura

Las ecuaciones del movimiento se escriben ahora

• Movimiento de traslación del c.m.

mg·senq -F_r=ma_c

• Movimiento de rotación alrededor de un eje que pasa por el c.m.

F_rR-μ·B·cosφ=I_ca

• Relación entre el movimiento de traslación y rotación (rueda sin deslizar)

a_c=a R

Despejamos la aceleración del centro de masas a_c, o mejor, el producto ma_c

Cuando el cilindro rueda sin deslizar, su c.m. se desplaza una distancia x a lo largo del plano inclinado. El ángulo girado por el momento magnético del imán μ respecto de su dirección inicial es φ=x/R, siendo R el radio del cilindro, tal como se muestra en la figura.

La aceleración del c.m. a_c no es constante sino que es función de la posición x. La ecuación del movimiento no se puede integrar fácilmente.

Curva de la energía potencial

El centro de masas del cilindro se mueve como una partícula cuya masa es igual a la masa m del cilindro, bajo la acción de la fuerza F=ma_c.

La fuerza F es una función de la posición x del c.m. del cilindro. Se trata por tanto, de una fuerza conservativa. El trabajo de la fuerza conservativa F, es igual a la diferencia entre la energía potencial inicial y la final.

Tomando E_p(0)=0.

La energía potencial la podemos escribir en términos de las variables adimensionales z=x/R, y A=μB/(mgRsenθ), y V(z)=E_p(x)/(2mgR/3), quedando

V(z)=-z+Asen(z)

Calculamos los extremos de esta función

Son máximos cuando la derivada segunda es negativa y son mínimos cuando la derivada segunda es positiva. Llamando f=arcos(1/A).

• Los máximos se producen para z_n=f+2nπ
• Los mínimos para z_n= -f+2nπ

En los extremos (máximos y mínimos) la fuerza F neta que actúa sobre el c.m. del cilindro es cero. Estas son las posiciones de equilibrio: estable en los mínimos e inestable en los máximos.

La representación gráfica de V(z) en función de la variable adimensional z, nos permite describir cualitativamente el movimiento del c.m. del cilindro para una energía total ε dada.

Como se aprecia en la figura, para la energía total ε hay dos regiones en las que se pueden mover el cilindro, aquellas en las que su energía cinética es positiva o bien aquellas en las que la energía total ε es mayor que la energía potencial V(z). La condición ε≥V(z), se cumple para los puntos z en el segmento z₁≤z≤z₂ y en la semirrecta z≥z₃ (señaladas con trazo grueso en la figura). Siendo z₁, z₂ y z₃ las raíces de la ecuación trascendente ε=V(z).

La posición inicial del c.m. del cilindro determina la región en la que se moverá posteriormente. Imaginemos que el cilindro se encuentra en la posición z≥z₃, señalada en la figura. Su energía potencial es el segmento CD, la energía total es el segmento CA, luego la energía cinética (de traslación y de rotación) es el segmento AB. El punto de intersección z₃ entre la recta V(z)= ε,  y la curva de energía potencial marca el punto de retorno, es decir, aquél en el que la velocidad del c..m. del cilindro es nula.

Cuando A=1, no hay máximos ni mínimos locales, los puntos z=0, ±2π, ±4π, son de inflexión ya que la primera y segunda derivadas son nulas en dichos puntos. dV/dz=0, y d2V/dz2=0.

Ecuación del movimiento del cilindro

La aceleración del c.m. a_c de cilindro no es constante, sino que depende de la posición x del cilindro. Escribimos la ecuación diferencial del movimiento en función del parámetro adimensional z=x/R.

Se resuelve la ecuación diferencial mediante el procedimiento numérico de Runge-Kutta con las condiciones iniciales z=z₀ y dz/dt=0 en el instante t=0, siendo z₀ un punto de retorno o una de las raíces de la ecuación trascendente ε=V(z).

El comportamiento del sistema viene determinado por el parámetro A, ya que el parámetro C actúa de factor de escala. Fijado el valor de C, lo que equivale a fijar el ángulo θ del plano inclinado para un valor dado del radio R del cilindro, se puede cambiar el parámetro A modificando el valor del campo magnético B, es decir, modificando la intensidad de la corriente que produce dicho campo.

1. Cuando A>1 se pueden dar dos comportamientos

• El cilindro oscila entre dos posiciones.
• El cilindro baja rodando sin deslizar a lo largo del plano inclinado.

3. Cuando A=1, la partícula permanece en equilibrio en z=0, 2π, -2π, .. si la energía total es ε=0,-2π, 2π,...
   Para otros valores de la energía, el cilindro baja rodando sin deslizar a lo largo del plano inclinado.

3. Para A<1 el cilindro siempre baja rodando a lo largo del plano inclinado.

Pequeñas oscilaciones alrededor de una posición de equilibrio

Si desarrollamos en serie la función f(z)=C(1-Acos(z)), alrededor de la posición de equilibrio (mínimo de la energía potencial) z0=-arccos(1/A).

La ecuación del movimiento se transforma en la ecuación diferencial de un M.A.S.

La frecuencia angular ω del M.A.S. es

Actividades

Se introduce

1. El valor del parámetro A, en el control de edición titulado Parámetro
2. El valor de la energía total ε, en el control de edición titulado Energía

Se pulsa el botón titulado Empieza

El programa interactivo calcula las raíces de la ecuación trascendente ε=V(z), situando al cilindro aleatoriamente es en una de las posiciones de retorno con velocidad inicial nula. A continuación, resuelve la ecuación diferencial del movimiento por el procedimiento de Runge-Kutta, para calcular la posición del cilindro en función del tiempo.

El parámetro C se ha fijado en el valor C=1, el radio del cilindro se ha fijado en R=1, lo que equivale a fijar el ángulo del plano inclinado en el valor θ=8.8º.

Cuando se introduce el valor de A, se representa gráficamente de V(z) en función de z, la energía total ε se representa mediante una línea horizontal, y un segmento vertical de color azul en la posición del c.m. del cilindro indica su energía cinética (de traslación y de rotación).

En la parte inferior del applet se observa el movimiento del cilindro sobre el plano inclinado. Con origen en su c.m. se han dibujado.

• El campo magnético (una flecha de color azul que apunta hacia arriba).
• El momento magnético del imán (una flecha de color rojo, que gira con el cilindro).
• La fuerza F neta que actúa sobre el c.m. del cilindro, (una flecha de longitud variable de color negro).

Fijarse que la fuerza F es nula en los máximos y en los mínimos de V(z), es positiva cuando la pendiente de V(z) es negativa y es negativa cuando la pendiente es positiva.

Ejemplos

Para A=4 observamos una sucesión de máximos y mínimos locales en las posiciones siguientes:

• Máximos: -4.97, 1.32, 7.6…

• Mínimos: -7.60, -1.32, 4.97

Introducimos una energía ε=2.2. Las raíces de la ecuación trascendente ε=V(z) son

z₁=-2.95, z₂=0.88, z₃=1.74

El cilindro tiene dos posibles movimientos

• Oscila entre las posiciones z₁≤z≤z₂
• Baja rodando sin deslizar desde la posición z₃ en adelante

Al pulsar el botón titulado Empieza, se selecciona aleatoriamente una de las tres posiciones de partida.

Cuando la energía total es ε=-1.5 próxima al mínimo local de V(z) el cilindro oscila con un periodo del 3.3 s entre las posiciones z₁=-2.05 y z₂=-0.53. Este periodo es un poco mayor que el que se obtiene mediante la fórmula aproximada.

Los valores del tiempo t, la posición z y la velocidad v del c.m. del cilindro se muestra en la esquina superior izquierda del applet. Cuando el cilindro sale de una de las posiciones extremas se empieza a contar el tiempo. Se puede medir el periodo, utilizando los botones titulados Pausa/Continua y Paso, deteniendo el movimiento del cilindro cuando su c.m. regrese a la posición de partida.

Para A=1. Si introducimos una energía ε=0, el cilindro permanece en esta posición de equilibrio.

Para A=0.3. Si introducimos una energía ε=0, el cilindro baja rodando a lo largo del plano inclinado desde la posición inicial z₀=0. Este tipo de movimiento se produce para cualquier valor de la energía ε.