]> Oscilaciones del péndulo-muelle.

Oscilaciones del péndulo-muelle.

En esta página, se estudia un sistema oscilante simple, que consiste en una partícula de masa m en el extremo inferior de un muelle elástico de constante k, fijado por el extremo superior. El muelle, inicialmente en equilibrio en posición vertical, se desplaza a una posición (x0, y0) y se suelta.

Observamos la combinación de dos modos de oscilación: la del péndulo simple y la del muelle elástico, ambos modos tienen su frecuencia característica y están acoplados de forma no lineal.

Como consecuencia, es difícil mantener oscilando hacia arriba y hacia abajo, el sistema formado por la partícula y el muelle elástico, el acoplamiento hace que pronto la partícula se desvíe hacia un lado y hacia el otro (como un péndulo). Se observa que estas oscilaciones son muy pronunciadas con una elección adecuada de la masa de la partícula, de la constante elástica o de la longitud del muelle sin deformar.

Este sistema no lineal de dos grados de libertad, muestra el papel importante que juegan las inestabilidades en la transferencia de energía desde un grado de libertad del sistema hacia el otro.

Ecuaciones del movimiento

Colocamos muelle elástico de constante k y longitud L0 sin deformar en posición vertical, lo fijamos por el extremo superior y colocamos una partícula de masa m en su extremo libre inferior.

El muelle elástico se estira una longitud ye=mg/k, la longitud del muelle es Le=L0+ye.

Situamos el origen del sistema de coordenadas en la posición de la partícula en equilibrio. Desplazamos la partícula a la posición (x0, y0) y la soltamos. Cuando la partícula está en la posición (x, y), la longitud del muelle deformado es L, las fuerzas sobre la partícula son:

Aplicamos la segunda ley de Newton

m d 2 x d t 2 =Tsinθ m d 2 y d t 2 =Tcosθmg

Sustituyendo T=k(L-L0), mg=k(Le-L0), y teniendo en cuenta que senθ=x/L, y cosθ=(Le-y)/L, el sistema de ecuaciones diferenciales se escribe en función de la posición (x, y) de la partícula.

d 2 x d t 2 = k m (L L 0 ) x L d 2 y d t 2 = k m (L L 0 ) ( L e y) L k m ( L e L 0 )

La longitud del muelle deformado es

L= ( L e y) 2 + x 2

Es importante hacer notar, que la no linealidad de las ecuaciones diferenciales procede de la geometría no del muelle cuyo comportamiento se supone lineal.

Se resuelve el sistema de ecuaciones diferenciales por procedimientos numéricos con las siguientes condiciones iniciales: en el instante t=0, la posición inicial es x0, y0, y parte del reposo.

Aproximaciones

Cuando los desplazamientos x e y de la posición de equilibrio son pequeños en comparación con la longitud del muelle L0, el sistema de ecuaciones diferenciales se puede expresar de forma más simple.

En la primera ecuación diferencial, el término

(L L 0 ) x L =x L 0 L x L 0 L = 1 x 2 L 0 2 + y 2 L 0 2 2 L e y L 0 2 + L e 2 L 0 2 1 2 L e y L 0 2 + L e 2 L 0 2 = L 0 L e ( 1 2y L e ) 1/2 L 0 L e + L 0 L e 2 y

La primera ecuación diferencial, se escribe

d 2 x d t 2 = k m ( 1 L 0 L e )x+ k m L 0 L e 2 xy d 2 x d t 2 + g L e x= k m L 0 L e 2 xy

En la segunda ecuación diferencial, el término

(L L 0 ) ( L e y) L ( L e L 0 )=y L 0 L ( L e y)+ L 0 =y L 0 ( L e y) ( L e y) 2 + x 2 + L 0 = y L 0 ( 1+ x 2 ( L e y) 2 ) 1/2 + L 0 y L 0 ( 1 1 2 x 2 L e 2 )+ L 0 =y+ 1 2 L 0 L e 2 x 2

La segunda ecuación diferencial, se escribe

d 2 y d t 2 + k m y= 1 2 k m L 0 L e 2 x 2

Definimos ωx la frecuencia de oscilación a largo del eje X, ωy la frecuencia de oscilación a lo largo del eje Y, y λ la constante de proporcionalidad del acoplamiento como

ω x 2 = g L e ω y 2 = k m λ= k m L 0 L e 2

El sistema de dos ecuaciones diferenciales de segundo orden acopladas queda

d 2 x d t 2 + ω x 2 x=λxy d 2 y d t 2 + ω y 2 y= λ 2 x 2

Resonancia autoparamétrica

Supongamos que comenzamos a hacer oscilar la partícula muy próxima al eje Y, con x<<y. La solución de la segunda ecuación diferencial (se desprecia el término λx2) es, aproximadamente

y=Acos(ωy·t)

La primera ecuación diferencial se convierte en

d 2 x d t 2 +( ω x 2 λAcos( ω y t) )x=0

Esta ecuación diferencial se denomina del tipo Mathieu. Cuando

ω y ω x = 2 n n=1, 2..

la solución es inestable incluso para pequeños valores del parámetro A (amplitud de la oscilación a lo largo del eje Y). Cuando menor sea n, más pronunciada es la inestabilidad. Así, para n=1, la relación entre las frecuencias angulares es

ωy=2ωx.

Los valores de la constante elástica k, de la longitud L0 del muelle sin deformar, y de la masa m de la partícula que hacen que se cumpla esta relación son

ω y 2 =4 ω x 2 k m =4 g L e L e = L 0 + mg k L 0 = 3mg k

Energía del sistema

En la situación inicial, t=0,  la partícula unida al muelle elástica se libera (v=0) en la posición (x0, y0). Situamos el nivel cero de energía potencial Ep=0, en la posición de la partícula cuando se encuentra en equilibrio, tal como se muestra en la figura.

La energía inicial es

E=mg y 0 + 1 2 k ( x 0 + ( L e y 0 ) 2 L 0 ) 2

En el instante t la partícula se encuentra en la posición (x, y) y lleva una velocidad v=vxi+vyj

La energía del sistema es

E= 1 2 m( v x 2 + v y 2 )+mgy+ 1 2 k ( x 2 + ( L e y) 2 L 0 ) 2

Actividades

Se introduce

Se pulsa el botón titulado Inicio

La partícula se sitúa en la posición de equilibrio a lo largo del eje vertical, el muelle se deforma ye=mg/k.

Se pulsa el botón titulado Empieza

Se observa el movimiento de la partícula, y su trayectoria

Nota: En el programa interactivo, el origen se sitúa en la posición del extremo libre del muelle sin deformar, L0=1 m por debajo del punto de sujeción del muelle, mientras que en los fundamentos físicos, se ha situado en la posición de equilibrio de la partícula, ye=mg/k por debajo del origen, o Le=L0+mg/k por debajo del punto de sujeción.

En el programa interactivo resuelve el sistema de dos ecuaciones diferencial de segundo orden por el método de Runge-Kutta. Este procedimiento numérico produce, en general, buenos resultados, que podemos confirmar observando que la energía se mantiene constante.

Ejemplo:

El valor de la constante k del muelle que hace que ωy=2ωx es k=8.82 N/m

En la posición de equilibrio el muelle se deforma ye=mg/k=1/3 m. La longitud del muelle en la posición de equilibrio es Le=1+1/3=4/3 m

Si partimos de una posición inicial muy próxima al eje vertical Y, por ejemplo, x0=0.01 m e y0=0.0, observamos como crecen las oscilaciones hacia uno y otro lado del eje vertical, del tipo péndulo, y disminuyen las oscilaciones tipo muelle. Luego, sucede el proceso inverso.

En la parte superior izquierda del applet, se muestra la energía total de la partícula y la posición de la partícula (x, y) de la partícula en cada instante t.

En la parte derecha del applet, se representan los tres tipos de energías mediante un diagrama de barras:

La suma de las tres clases de energía se mantiene constante.

LibresApplet 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

Rusbridge M.G., Motion of the sprung pendulum. Am. J. Phys. 48 (2) February 1980, pp. 146-151.

Olsson M. G. Why does a mass on a spring sometimes misbehave?. Am. J. Phys. 44 (12) December 1976, pp. 1211-1212.