Caída de un 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)

En la página, Campo magnético producido por una imán hemos supuesto que un imán se comporta como un dipolo magnético de momento μ. Las componentes del campo son

B ρ = μ 0 μ 4π 3ρz ( ρ 2 + z 2 ) 5 2 B z = μ 0 μ 4π 2 z 2 ρ 2 ( ρ 2 + z 2 ) 5 2

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

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

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

Φ= 0 b μ 0 ·μ 4π ( 3 z 2 ( ρ 2 + z 2 ) 5 1 ( ρ 2 + z 2 ) 3 )2πρ·dρ Φ= μ 0 μ 4 ( 3 z 2 0 b 2ρ·dρ ( ρ 2 + z 2 ) 5 0 b 2ρ·dρ ( ρ 2 + z 2 ) 3 )= μ 0 μ· a 2 2 ( a 2 + z 2 ) 3

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 de la corriente 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, 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:

function a = caida_iman_1(t,x,k) 
    a=zeros(3,1);
    a(1)=x(2);    
%k(1) es alfa y k(2) es beta
%x(1) es posición, x(2) es velocidad y x(3) es intensidad
    a(2)=-1-1.5*x(1)*x(3)/(1+x(1)^2)^(5/2);
	a(3)=-k(1)*x(3)+1.5*k(2)*x(1)*x(2)/(1+x(1)^2)^(5/2);
end
x0=zeros(1,3);
x0(1)=4; %valor inicial de x
x0(2)=0; %valor inicial de v;
x0(3)=0; %intensidad inicial en la espira
tf=10; %tiempo final;
k=[0 35]; %valor de los parámetros alfa y beta
tspan=[0 tf]; 
[t,x]=ode45(@caida_iman_1,tspan,x0, [],k);

plot(t,x(:,1))  %posición en función del tiempo
xlabel('t')
ylabel('x');
grid on
title('Movimiento de un imán')

Para representar la intensidad modificamos en el script la sentencia plot

....
plot(t,x(:,3),'r')  %intensidad en función del tiempo
ylabel('i');
...

Actividades

Se introduce

Se pulsa el botón titulado Nuevo

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

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

En la parte derecha, se representa:

En el extremo derecho, 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, 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

Referencias

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