]> El oscilador de Fermi
Siguiente

El oscilador de Fermi

Consideremos el siguiente sistema oscilante: una partícula de masa m, se deja caer desde una altura h, sobre una plataforma de masa M unida a un muelle elástico de constante k. La partícula cae y choca elásticamente con la plataforma, el bloque asciende y desciende, la plataforma oscila y se vuelven a encontrar al cabo de un cierto tiempo, chocan y así, sucesivamente.

Se aplica la segunda ley de Newton para estudiar el movimiento del sistema formado por el muelle y la plataforma, y las ecuaciones del movimiento de caída libre para estudiar el movimiento del bloque. En el momento del choque, se aplica el principio de conservación del momento lineal y de la energía.

El sistema es conservativo, la energía del sistema formado por el bloque, el muelle y la plataforma se mantiene constante.

Ecuaciones del movimiento

Estudiamos las distintas etapas del movimiento del bloque y de la plataforma, para ello establecemos el origen en la posición del extremo del muelle sin deformar. Vamos a disponer de dos relojes uno que mide el tiempo parcial t de cada una de las etapas del movimiento y otro que mide el tiempo total tt, desde el momento en el que se libera el bloque a una altura h sobre el origen.

Situación inicial

En la figura, se muestra la situación inicial. El bloque de masa m está a una altura h sobre el origen. La plataforma de masa M descansa sobre el muelle elástico de constante k. La plataforma desciende una altura ye, hasta que se equilibran el peso Mg y la fuerza que ejerce el muelle k·ye.

ye=Mg/k

La energía inicial del sistema formado por el bloque, el muelle y la plataforma es

E 0 =mghMg y e + 1 2 k y e 2

Movimiento del bloque

Se libera el bloque, las ecuaciones del movimiento son:

x=h 1 2 g t 2 v x =gt

El bloque permanece en reposo, en la posición de equilibrio, tal como se muestra en la figura

y=-ye
vy
=0

El choque ente el bloque y la plataforma se produce en la posición x0=-ye, en el instante t1.

t 1 = 2(h+ y e ) g

La velocidad del bloque es

vx=-g·t1

Choque elástico

 

Supondremos que el bloque y la plataforma forman un sistema aislado en el momento del choque. El momento lineal se conserva. Si además el choque es elástico, la energía cinética antes y después del choque es la misma.

m v x +M v y =m V x +M V y 1 2 m v x 2 + 1 2 M v y 2 = 1 2 m V x 2 + 1 2 M V y 2

donde vx y Vx son las velocidades del bloque antes y después del choque, y vy y Vy son las velocidades de la plataforma antes y después del choque.

Escribimos las ecuaciones de la forma alternativa

m( v x V x )=M( v y V y ) m( v x 2 V x 2 )=M( v y 2 V y 2 )

La diferencia de los cuadrados de dos números es igual al producto de su suma por su diferencia

m( v x V x )=M( v y V y ) m( v x V x )( v x + V x )=M( v y V y )( v y + V y )

Nos queda un sistema de dos ecuaciones más fácil de resolver

m( v x V x )=M( v y V y ) ( v x + V x )=( v y + V y )

Despejamos las velocidades del bloque Vx y de la plataforma Vy después del choque

V x = 2M v y (Mm) v x m+M V y = 2m v x +(Mm) v y m+M

Movimiento del bloque después del choque

En el instante t=0, (tiempo total tt=t1) el bloque parte de la posición x0=-ye, con velocidad inicial v0x=Vx

x= x 0 + v 0x t+ 1 2 g t 2 v x = v 0x gt

Movimiento de la plataforma después del choque.

En el instante t=0, (tiempo total tt=t1) la plataforma parte de la posición y0=-ye, con velocidad inicial v0y=Vy

Las fuerzas sobre la plataforma son el peso y la fuerza que ejerce el muelle. La segunda ley de Newton se escribe

Ma=-ky-Mg

En forma de ecuación diferencial

d 2 y d t 2 + ω 2 y=g ω 2 = k M

La solución de esta ecuación diferencial es la suma de la solución de la ecuación diferencial homogénea

y1=Asin(ωt)+Bcos(ωt

donde A y B son constantes que se determinan a partir de las condiciones iniciales

y de la solución particular

y2=C

Introduciendo la solución particular en la ecuación diferencial determinamos el valor de la constante C.

ω2C=-g

La solución de la ecuación diferencial completa es y=y1+y2

y= g ω 2 +Asin(ωt)+Bcos(ωt) v y = dy dt =Aωcos(ωt)Bωsin(ωt)

Las condiciones iniciales son las siguientes: la plataforma parte en el instante t=0 (tiempo total tt=t1) de la posición y0, con velocidad v0y. Resolvemos el sistema de dos ecuaciones con dos incógnitas para determinar A y B.

y= g ω 2 + v 0y ω sin(ωt)+( y 0 + g ω 2 )cos(ωt) v y = v 0y cos(ωt)( y 0 ω+ g ω )sin(ωt)

Segundo choque

El próximo choque se produce en el instante t2, (tiempo total tt=t1+t2) cuando coinciden la posición x del bloque y la posición y de la plataforma. Para ello, hay que resolver por procedimientos numéricos la ecuación trascendente 

x 0 + v 0x t+ 1 2 g t 2 = g ω 2 + v 0y ω sin(ωt)+( y 0 + g ω 2 )cos(ωt)

Se obtiene la raíz t2 de la ecuación y se calcula la posición xc=yc de encuentro en dicho instante.

x c = y c = x 0 + v 0x t 2 + 1 2 g t 2 2

La velocidad del bloque antes del choque es

v x = v 0x g t 2

La velocidad de la plataforma antes del choque es

v y = v 0y cos(ω t 2 )( y 0 ω+ g ω )sin(ω t 2 )

Se calculan las velocidades Vx y Vy del bloque y de la plataforma después del choque. Estas son las velocidades v0x=Vx y v0y=Vy iniciales de ambos cuerpos en el instante t=0 (tiempo total tt=t1+t2) cuando se encuentran en la posición x0=xc, y0=yc calculada en el instante t2.

El bloque asciende y la plataforma oscila, vuelven a chocar y así, sucesivamente…

La energía del sistema formado por el bloque, el muelle y la plataforma en cualquier instante tt es

E= 1 2 m v x 2 +mgh+ 1 2 M v y 2 +Mgy+ 1 2 k y 2

La energía total permanece constante

Sensibilidad a un pequeño cambio de las condiciones iniciales

Unos de los aspectos más interesantes de los sistemas que estudiamos en esta parte del capítulo "Introducción al régimen caótico" es la de examinar la respuesta del sistema ante un pequeño cambio de las condiciones iniciales.

Cómo evoluciona el sistema si cambiamos ligeramente la altura h desde la que se deja caer el bloque. Veamos algunos ejemplos

Ejemplo 1:

En la figura, vemos la representación gráfica x-t del bloque para dos alturas iniciales distintas

 

Las posiciones del bloque en el instante tt=10 s son muy parecidas

Ejemplo 2:

En la figura, vemos la representación gráfica x-t del bloque para dos alturas iniciales distintas

 

Las posiciones del bloque en el instante tt=10 s difieren bastante

Actividades

Se introduce

Se pulsa el botón titulado Inicio

El programa verifica si la máxima deformación inicial que experimenta el muelle xm después de chocar por primera vez el bloque y la plataforma es superior a 0.5 m, en cuyo caso se invita al usuario a disminuir la masa del bloque o a aumentar el valor de la constante elástica.

Se pulsa el botón titulado Empieza

Se observa el movimiento del bloque y de la plataforma en la parte izquierda del applet, y la representación gráfica de la posición x del bloque en función del tiempo total tt .

El programa calcula en cada instante el tanto por ciento de error relativo en la energía o el cociente

| E E 0 E 0 |·100

donde E es la energía del sistema en cualquier instante tt, y E0 es la energía inicial del sistema.

Este valor se proporciona en caracteres de color rojo en la parte inferior izquierda del applet. Su valor debe ser siempre cero, o un valor muy pequeño lo que indica que la energía del sistema permanece constante y el programa realiza los cálculos correctamente. Si el procedimiento numérico empleado, es incapaz de resolver dos choques consecutivos muy próximos, el error se incrementa, y el programa debe de interrumpirse manualmente ya que no realiza los cálculos correctamente.

Podemos comparar la posición del bloque en función del tiempo para dos alturas iniciales distintas del siguiente modo (véase las dos figuras más arriba):

Cuando se completa una gráfica x-t (en color rojo o azul) en el intervalo deseado 0-10, 10-20, 20-30 s…, se copia la imagen del applet en el portapapeles y se pega en el programa de dibujo Paint que viene con el sistema operativo Windows.

Se cambia ligeramente la altura inicial h y se completa la gráfica x-t (en color azul o rojo), en el mismo intervalo. Se copia la imagen del applet en el portapapeles y se pega en el programa de dibujo Paint, en la misma posición que el applet anterior.

Se debe desactivar antes de realizar estas operaciones la opción del menú Imagen/Dibujar figuras opacas del programa Paint.

CaoticoApplet aparecerá en un explorador compatible con JDK 1.1.

Referencias

Lopac V., Dananic V. Energy conservation and chaos in the gravitationally driven Fermi oscillator. Am. J. Phys. 66 (10) October 1998, pp. 892-902Siguiente