Anterior

Caída de un imán

En las páginas tituladas “Demostración de la ley de Faraday”, hemos estudiado la corriente inducida producida en una espira, cuando un imán se mueve a lo largo de su eje con velocidad constante. La corriente inducida no afecta al movimiento del imán.

En esta página, vemos a estudiar un ejemplo algo más complicado. El imán se sitúa a cierta altura, se libera y cae bajo la acción de la gravedad hacia la espira a lo largo de su eje. Se originan corrientes inducidas en la espira que van a modificar el movimiento del imán, como veremos a continuación.

En la página "Movimiento de un imán en un tubo metálico vertical" no se han considerado los efectos de la autoinducción de la espira. En esta página, la espira va a tener una autoinducción no nula y esto hace cambiar la descripción del movimiento del imán.

Fuerzas sobre el imán

El campo magnético producido por una espira de radio a por la que circula una corriente eléctrica de intensidad I, en un punto z de su eje es

B z = μ 0 I 2 a 2 ( z 2 + a 2 ) 3/2

Consideramos el imán como un dipolo de momento μ= μk

La energía potencial de un dipolo de momento magnético μ en un campo magnético B que tiene la dirección del eje Z es el producto escalar

U=-μ·B=-μ·Bz

Como B es variable a lo largo del eje de la espira, el dipolo magnético experimenta una fuerza

F z = U z = 3μ μ 0 I a 2 2 z ( z 2 + a 2 ) 5/2

Si la corriente I en la espira es negativa (en el sentido de las agujas del reloj) la fuerza es repulsiva (las corrientes se repelen si tienen sentido contrario y se atraen, si tienen el mismo sentido).

Aplicamos la segunda ley de Newton al movimiento del imán

m d 2 z d t 2 =mg 3μ μ 0 I a 2 2 z ( z 2 + a 2 ) 5/2

Ecuación del circuito (espira)

Si consideramos que el imán es un dipolo magnético de momento magnético μ=iπa2, el campo magnético producido por el imán tiene la siguientes componentes. En la figura, se muestran las líneas de campo magnético producido por el imán

B y = μ 0 ·μ 4π r 3 ( 3yz r 2 ) B z = μ 0 ·μ 4π r 3 ( 3 z 2 r 2 1 )

El flujo del campo producido por el imán a través de una espira de radio a es.

Φ= S BdS = S ( B y j ^ + B z k ^ )dS k ^ = S B z ·dS

Dado que el plano de la espira es perpendicular al eje Z, el flujo de la componente Y del campo es nulo.

El elemento diferencial de superficie dS, es el área de un anillo de radio y y de espesor dy, su valor es dS=2πy·dy

Φ= 0 a μ 0 ·μ 4π r 3 ( 3 z 2 r 2 1 )·2πy·dy · Φ= π μ 0 ·μ 4π ( 3 z 2 0 a 2ydy ( y 2 + z 2 ) 5 0 a 2ydy ( y 2 + z 2 ) 3 )= μ 0 μ· a 2 2 ( a 2 + z 2 ) 3/2

Aplicando la ley de Faraday

V ε = dΦ dt = dΦ dz dz dt = 3 μ 0 μ a 2 2 z ( a 2 + z 2 ) 5/2 dz dt

La espira tiene una resistencia R y una autoinducción L. La espira es equivalente al circuito de la figura, cuya ecuación es

V ε L dI dt =IR

Recuérdese que la autoinducción actúa como una batería que se opone a los cambios en la corriente I de la espira

L dI dt +IR= 3 μ 0 μ a 2 2 z ( a 2 + z 2 ) 5/2 dz dt

Balance energético

Las energías del sistema formado por la espira y el imán son:

Se deberá cumplir que

d dt ( mgz+ 1 2 m v 2 + 1 2 L I 2 )= I 2 R

Si la espira es un superconductor R=0, la suma de las tres clases de energía permanece constante.

Solución de las ecuaciones del movimiento

Para determinar el movimiento del imán tenemos que resolver un sistema de dos ecuaciones diferenciales acopladas, con las condiciones iniciales siguientes: en el instante t=0, la posición del imán es z0, su velocidad inicial es dz/dt=v0 y la intensidad inicial en la espira es I0.

m d 2 z d t 2 =mg 3μ μ 0 I a 2 2 z ( z 2 + a 2 ) 5/2 L dI dt +IR= 3 μ 0 μ a 2 2 z ( a 2 + z 2 ) 5/2 dz dt

Escribimos las ecuaciones diferenciales en forma adimensional, definiendo las nuevas variables x, τ e i.

z=xat=τ a g I= mg a 2 μ 0 μ i

El sistema de dos ecuaciones diferenciales en términos de las variables adimensionales x, τ e i, se escribe.

d 2 x d τ 2 =1 3 2 x ( x 2 +1) 5/2 i di dτ + R L a g i= (μ μ 0 ) 2 mg a 3 L 3 2 x ( x 2 +1 ) 5/2 dx dτ

Definimos los parámetros α y β

d 2 x d τ 2 =1 3 2 x ( x 2 +1 ) 5/2 i di dτ =αi+ 3β 2 x ( x 2 +1 ) 5/2 dx dτ α= R L a g β= (μ μ 0 ) 2 mg a 3 L

El programa interactivo que viene a continuación, resuelve el sistema de dos ecuaciones diferenciales por el procedimiento de Runge-Kutta, con las siguientes condiciones iniciales, en el instante τ=0, la posición inicial del imán es x0, y parte del reposo, su velocidad inicial dx/dτ=0, la intensidad inicial en la espira i0=0.

En función de las variables adimensionales x, τ e i, los distintos tipos de energía se escriben:

mgz=-(mga)x

1 2 m ( dz dt ) 2 =(mga) 1 2 ( dx dτ ) 2

1 2 L I 2 = 1 2 L ( mg a 2 μ 0 μ ) 2 i 2 =(mga) 1 2 i 2 β

Actividades

Se pulsa el botón titulado Inicio

Se introduce

Se pulsa el botón titulado Empieza

Observamos el movimiento del imán, las fuerzas sobre el mismo, que en términos de las variables adimensionales son:

3 2 x ( x 2 +1 ) 5/2 i

Positiva si la corriente i es en el sentido de las agujas del reloj y negativa, en caso contrario.

Los puntos rojos en la espira, representan cargas positivas cuyo movimiento nos sugiere el sentido de la corriente inducida.

En la parte derecha del applet, se representa:

En el extremo derecho del applet, un diagrama de barras representa el balance energético.

La energía inicial del imán es x0

La energía en un instante τ es la suma de la energía

Cuando el imán cae, la energía potencial gravitatoria, se convierte por una parte en energía cinética y por la otra, en energía del campo magnético producido por la espira. Si la espira tiene resistencia, una parte de la energía se disipa en forma de calor por efecto Joule, la energía total (suma de las tres clases de energías) va disminuyendo con el tiempo.

Cuando la espira es superconductora, R=0, y α=0, el imán puede levitar sobre la espira, se mueve hacia abajo y hacia arriba, la energía se conserva, y se puede medir el periodo de las oscilaciones. Si la resistencia es pequeña, se observan oscilaciones amortiguadas, cuya amplitud va disminuyendo con el tiempo.

En la siguiente tabla, se sugieren algunos valores de los parámetros α y β, y la posición inicial x0

Alfa Beta Posición inicial
1 35 3.0
0 25 3.0
0 35 3.0
0 18 0.0
SolenoideApplet aparecerá en un explorador compatible JDK 1.1

Arrastre la flecha horizontal de color azul con el puntero del ratón

Referencias

Belcher J. W. The falling magnet. web.mit.edu/8.02t/www/802TEAL3D/visualizations/resources/fallingMagnetEqs.pdf

Anterior