Espira en un campo magnético no homogéneo

Autoinducción nula (R≠ 0, L=0)

Si la espira de anchura a, se ha introducido x en la región en la que existe campo magnético que apunta perpendicularmente al plano del dibujo y hacia dentro, el flujo a través de la espira será

Φ= B · S =B(ax)

Aplicando la ley de Faraday se obtiene la fem

V ε = dΦ dt =Bav

De acuerdo a la ley de Lenz el sentido de la corriente inducida i es antihorario, ya que el flujo aumenta.

Para una espira de resistencia R, la ecuación del circuito es Vε =iR

i= Bav R

La fuerza sobre el lado derecho de la espira es

F =i· u ^ t × B ·lF=iBa= B 2 a 2 R v

de sentido contrario a la velocidad v de la espira.

La ecuación del movimiento de la espira será

m dv dt = B 2 a 2 R v   

y la solución de la ecuación diferencial, con las condiciones iniciales t=0, v=v0, es

v 0 v dv v = B 2 a 2 Rm 0 t dt v= v 0 exp( t τ ),τ= Rm B 2 a 2

La velocidad v de la espira disminuye exponencialmente con el tiempo.

A partir de la expresión de la velocidad en función del tiempo obtenemos la posición x del móvil que parte del origen en el instante inicial t=0, x=0

x= 0 t vdt= v 0 τ( 1exp( t τ ) )

La intensidad en función del tiempo es

i= Ba R v 0 exp( t τ )

La intensidad i que circula por la espira disminuye exponencialmente con el tiempo.

Alternativamente, expresamos la velocidad v en función de la posición x

dv dx dx dt = B 2 a 2 mR v dv dx = B 2 a 2 mR v 0 v dv = 1 τ 0 x dx v= v 0 x τ

La velocidad disminuye linealmente con x

Estas ecuaciones se aplican mientras que el lado derecho de la espira esté en el seno del campo magnético. Cuando ambos lados estén dentro de la región x≥ 0, la intensidad de la corriente inducida será cero y la fuerza que ejerce el campo magnético sobre la espira será nula y por tanto, la velocidad de la espira será constante.

Estudio energético

Comprobamos que la energía cinética inicial de la espira Ek, se disipa en la resistencia de la espira. La ley de Joule afirma que la energía por unidad de tiempo (potencia) disipada en la resistencia es i2R.

En cualquier instante, la suma de la energía cinética Ek de la espira y de la energía disipada en la resistencia ER es igual a la energía cinética inicial de la espira.

E k = 1 2 m v 2 = 1 2 m v 0 2 exp(2t/τ) E R = 0 t i 2 R·dt= 1 2 m v 0 2 ( 1exp(2t/τ) ) E k + E R = 1 2 m v 0 2

Ejemplo

Una espira cuadrada de lado a entra en una región de anchura b>a donde hay un campo magnético constante B uniforme apuntando hacia dentro de la pantalla. La velocidad inicial de la espira es v0.

Vamos a calcular la velocidad final v2 de la espira cuando sale de dicha región y el tiempo T empleado en atravesarla.

Se trata de una variante de la situación estudiada en la página titulada Una espira que se mueve a través de un campo magnético uniforme

Representamos la velocidad v de la espira en función de la posición x del lado derecho de la espira

El tiempo total T=t1+t2+t3 es

T=τln( v 0 τ v 0 τa )+ ba v 0 τa τ+τln( v 0 τa v 0 τ2a )= τ{ ba v 0 τa +ln( v 0 τ v 0 τ2a ) }

Resistencia nula (L≠ 0, R=0)

Consideremos que la espira está hecha de un material superconductor de modo que R≈0. En este caso, no podemos ignorar la autoinducción L, que produce una fem

V L =L di dt

La ecuación del circuito (suma de fems igual a intensidad por resistencia) se escribe ahora VL+Vε =0.

L di dt +Ba dx dt =0i= Ba L x

que satisface la condición inicial x=0, i=0.

La fuerza que ejerce el campo magnético sobre el lado derecho de la espira será de nuevo F=-iaB, por lo que la ecuación del movimiento se escribirá

m d 2 x d t 2 + B 2 a 2 L x=0

Esta es la ecuación diferencial de un Movimiento Armónico Simple, x=A·sin(ω0 t+φ ).

A partir de las condiciones iniciales para t=0, x=0, dx/dt=v0 calculamos la amplitud A y la fase inicial φ .

x= v 0 ω 0 sin ω 0 t ω 0 = Ba mL v= v 0 cos ω 0 t

Siempre que la amplitud de la oscilación v00 no sea mayor que la longitud de la espira, de modo que la espira no esté completamente introducida en la región x≥ 0, la espira describirá un MAS con frecuencia angular ω0 y periodo P=2π/ω0.

La intensidad de la corriente vale

i= Ba L x= v 0 m L sin( ω 0 t)

La máxima intensidad no depende de la anchura a de la espira, ni del campo magnético B.

Estudio energético

La energía cinética de la espira más la energía acumulada en la autoinducción en forma de campo magnético debe de ser constante e igual a la energía cinética inicial.

E k = 1 2 m v 2 = 1 2 m v 0 2 cos 2 ( ω 0 t) E L = 1 2 L i 2 = 1 2 B 2 a 2 L 2 x 2 = 1 2 m v 0 2 sin 2 ( ω 0 t) E k + E L = 1 2 m v 0 2

Actividades

En el programa interactivo que viene a continuación, se estudia el comportamiento de una espira que viaja hacia una región (x≥ 0) en la que existe un campo magnético uniforme perpendicular al plano de la espira: hacia adentro, negativo, en color azul claro o hacia fuera, positivo, en color rosa.

Se introduce

Tenemos dos opciones

Se pulsa el botón titulado Nuevo para situar la espira en la posición de partida, y el botón   ►   para poner en movimiento la espira.

En el caso de que la espira tenga autoinducción no nula, el lado izquierdo de la espira estará siempre en la región x<0, en la que no hay campo magnético. Si se introduce una velocidad v0 excesiva de la espira el programa no prosigue y nos invita a modificar el valor de la velocidad inicial de la espira. Para ello, la longitud de la espira ha de ser mayor que la amplitud v00 de las oscilaciones de la espira.

Ejemplo 1.

La constante de tiempo τ vale

τ= Rm B 2 a 2 = 0.5· 10 3 ·0.1 0.1 2 · 0.25 2 =0.08

La espira se introduce en la región en la que existe campo magnético.

En el instante t=0.2 s la velocidad de la espira es

v=2.0·exp(0.2/0.08)=0.16m/s

La posición del lado de la espira es

x=2·0.08( 1exp(0.2/0.08) )=0.15m

La intensidad de la corriente en este instante vale

i= Bav R = 0.1·0.25·0.16 0.5· 10 3 =8.21A

Como el flujo aumenta, la intensidad (en sentido antihorario) se opone al aumento de flujo

Ejemplo 2.

La frecuencia angular ω0 del MAS vale

ω 0 = Ba mL = 0.1·0.25 0.1·8· 10 5 =8.84rad/s

El semiperiodo, tiempo que tarda en entrar y salir la espira es P/2=π/ ω0=0.36 s

La espira penetra hasta una distancia

x m = v 0 ω 0 = 2.0 8.84 =0.23m

La intensidad máxima en este instante es

i m = Ba L x m = 0.1·0.25 8· 10 5 0.23=70.7A



Resistencia y autoinducción no nulas (L≠ 0, R≠ 0)

En este apartado se estudian las oscilaciones amortiguadas, críticas y sobreamortiguadas de la espira, siempre que que cumpla la condición de que solamente su lado superior esté en el interior del campo magnético x≥0.

Consideremos una región en la que existe un campo magnético uniforme B perpendicular al plano de la espira, en todos los puntos x≤0, y está libre de campo para x>0. Una espira de dimensiones a y h, se mueve hacia la región donde el campo es nulo x>0. El flujo del campo magnético a través de la espira disminuye y se generará una corriente inducido en el sentido contrario a las agujas del reloj tal como se incia en la figura.

El flujo total a través de la espira es la suma del flujo producido por el campo magnético externo, Ba(h-x) y el flujo propio, Li. Donde L es la autoinducción d ela espira

Φ=Ba(hx)+Li

La fem inducida es

V ε = dΦ dt =Ba dx dt L di dt

Despejando la intensidad i en la ecuación del movimiento e introduciéndola en la ecuación del circuito, llegamos a la siguiente ecuación diferencial de segundo orden

d 2 v d t 2 + R L dv dt + B 2 a 2 mL v=g R L

Que describe las oscilaciones amortiguadas

d 2 v d t 2 +2γ dv dt + ω 0 2 v=2γg ω 0 2 = B 2 a 2 mL γ= R 2L

La solución particular de la ecuación diferencial es una constante C, cuyo valor es

C= 2γg ω 0 2 = mgR B 2 a 2

La ecuación diferencial homogénea con las condiciones iniciales:

t=0,v=0, dv dt =0

tiene tres posibles soluciones:

Oscilaciones amortiguadas

Si la resistencia R no es muy grande de modo que γ<ω0

v=exp(γt)( Acos(ωt)+Bsin(ωt) )+ 2γg ω 0 2 , ω 2 = ω 0 2 γ 2 dv dt =γexp(γt)( Acos(ωt)+Bsin(ωt) )+ωexp(γt)( Asin(ωt)+Bcos(ωt) )

donde ω es la frecuencia de las oscilaciones amortiguadas. Las constantes A y B se determinan a partir de las condiciones iniciales

t=0{ 0=A+ 2γg ω 0 2 g=γA+ωB v=g{ 2γ ω 0 2 +exp(γt)( 1 ω ( 1 2 γ 2 ω 0 2 )sin(ωt) 2γ ω 0 2 cos(ωt) ) } dv dt =gexp(γt)( γ ω sin(ωt)+cos(ωt) )

Despejamos la intensidad i de la ecuación del movimiento

i= 1 Ba ( mgm dv dt ) i= mg Ba ( 1exp(γt)( γ ω sin(ωt)+cos(ωt) ) )

Para obtener la posición x del lado inferior de la espira integramos dx/dt=v, sabiendo que para t=0, x=0, parte del origen

x=g{ 2γ ω 0 2 t+ 1 ω ( 1 2 γ 2 ω 0 2 ) 0 t exp(γt) sin(ωt)·dt 2γ ω 0 2 0 t exp(γt) cos(ωt)·dt }

Integrando dos veces por partes,

exp(γt)cos(ωt)·dt= exp(γt) ω 0 2 { γcos(ωt)+ωsin(ωt) } exp(γt)sin(ωt)·dt= exp(γt) ω 0 2 { γsin(ωt)+ωcos(ωt) }

>> syms g w t;
>> int(exp(-g*t)*sin(w*t),t)
ans =-(exp(-g*t)*(g*sin(t*w) + w*cos(t*w)))/(g^2 + w^2)
>> int(exp(-g*t)*cos(w*t),t)
ans =-(exp(-g*t)*(g*cos(t*w) - w*sin(t*w)))/(g^2 + w^2)

obtenemos la posición x, del lado inferior de la espira en función del tiempo

x= g ω 0 2 { 2γ·t+exp(γt)( ( 4 γ 2 ω 0 2 1 )cos(ωt) γ ω ( 3 4 γ 2 ω 0 2 )sin(ωt) ) 4 γ 2 ω 0 2 +1 }

Después de un tiempo muy grande, t→∞, la velocidad v y la intensidad i tienden a

v = 2γg ω 0 2 = mgR B 2 a 2 i = mg Ba

Ejemplo

Fijamos los siguientes parámetros

Cambiamos los valores de los parámetros

Conocidos estos, calculamos el campo magnético B= ω 0 a mL y la resistencia de la espira R=2

En el script definimos las funciones que nos calculan: la velocidad v, intensidad de la corriente inducida en la espira, i y posición x del lado inferior de la espira rectangular, en función del tiempo t

m=0.1; %masa
a=0.25; %ancho de la espira rectangular
L=8e-5; %autoinducción

w0=10; %frecuencia propia
gamma=2; %constante de amortiguamiento
R=2*L*gamma; %resistencia
B=w0*sqrt(m*L)/a; %campo magnético

w=sqrt(w0^2-gamma^2); %frecuencia de la oscilación amortiguada
%velocidad
v=@(t) 9.8*(2*gamma/w0^2+exp(-gamma*t).*((1-2*gamma^2/w0^2)*sin(w*t)/w-
2*gamma*cos(w*t)/w0^2));
%intensidad
i=@(t) m*9.8*(1-exp(-gamma*t).*(gamma*sin(w*t)/w+cos(w*t)))/(B*a);
%posición del lado inferior de la espira
x=@(t) 9.8*(2*gamma*t+exp(-gamma*t).*((4*gamma^2/w0^2-1)*cos(w*t)+
gamma*(4*gamma^2/w0^2-3)*sin(w*t)/w)-4*gamma^2/w0^2+1)/w0^2;
fplot(v,[0,5*2*pi/w0]) %velocidad
fplot(i,[0,5*2*pi/w0]) %intensidad
fplot(x,[0,3*2*pi/w0]) %posición
grid on
xlabel('t')
ylabel('v')
title('Amortiguadas, velocidad')

Modificamos la sentencia fplot para representar la intensidad i

Modificamos la sentencia fplot para representar la posición x del lado inferior de la espira rectangular

Superconductor

La resistencia de una espira hecha de un material superconductor es nula, R=0, γ=0, ω=ω0

v= g ω 0 sin( ω 0 t)

Como v=dx/dt, integramos esta ecuación diferencial, sabiendo que en el instante t=0, x=0

x= g ω 0 2 ( 1cos( ω 0 t) )

Obtenemos esta misma ecuación, suponiendo que en el caso de un superconductor, el flujo a través de la espira es constante, -Bax+Li=constante. En el instante t=0, x=0, i=0, por lo que i=Bax/L. Sustituyendo en la ecuación del movimiento

m dv dt =mg Ba L xBa d 2 x d t 2 + ω 0 2 x=g

Cuya solución es (particular más homogénea)

x=Acos( ω 0 t)+Bsin( ω 0 t)+ g ω 0 2 v= dx dt = ω 0 ( Asin( ω 0 t)+Bcos( ω 0 t) ) t=0{ 0=A+ g ω 0 2 0= ω 0 B x= g ω 0 2 ( 1cos( ω 0 t) )

El lado inferior de la espira oscila entre x=0 y

x máx = 2g ω 0 2

con frecuencia angular ω0. La intensidad varía de i=0 cuando x=0 a

i máx = 2mg Ba

Energías

En el instante t=0, la posición del lado inferior de la espira es x=0, parte del reposo, v=0, y la intensidad inducida i=0, por lo que la energía inicial es cero. La energía de la espira en el instante t es la suma de

La suma de las tres contribuciones es, Ek+Ep+EL=0

Ejemplo

Un posible experimento para mostrar las oscilaciones de la espira superconductora R=0, sería el siguiente:

m=1.4e-3; %masa
a=0.01; %ancho de la espira rectangular
B=0.05; %campo magnético
L=1e-8; %autoinducción

w0=sqrt(B^2*a^2/(m*L)); %frecuencia angular
x=@(t) 9.8*(1-cos(w0*t))/w0^2;
fplot(x,[0,5*2*pi/w0]) %cinco periodos
grid on
xlabel('t (s)')
ylabel('x (m)')
title('Espira superconductora, posición')

El periodo de las oscilaciones sería, 2π/ω0=0.047 s. Su lado inferior oscilaría entre las posiciones x=0, y xmáx=0.0011 m. El lado inferior estaría en la región x>0, donde no hay campo magnético y el lado superior estaría en la región x<0, donde hay campo magnético B. La máxima intensidad que circularía por la espira, imáx=54.9 A

Oscilaciones críticas

Cuando la resistencia aumenta, puede ocurrir que γ=ω0.

La solución de la ecuación diferencial es

v=exp(γt)( At+B )+ 2g γ dv dt =exp(γt)( γAt+AγB )

Donde A y B se determinan a partir de las condiciones iniciales

t=0{ 0=B+ 2g γ g=AγB v= 2g γ exp(γt){ gt+ 2g γ } dv dt =g( γt+1 )exp(γt)

Despejamos la intensidad i de la ecuación del movimiento

i= 1 Ba ( mgm dv dt ) i= mg Ba { 1( γt+1 )exp(γt) }

Para obtener la posición x del lado inferior de la la espira integramos dx/dt=v, sabiendo que para t=0, x=0, parte del origen

x= 2g γ tg 0 t texp(γt)·dt 2g γ 0 t exp(γt)·dt

Integrando por partes,

0 t texp(γt)·dt = 1 γ ( t+ 1 γ )exp(γt)

obtenemos la posición x, del lado inferior de la espira en función del tiempo.

x= 2g γ t+ g γ texp(γt)+ 3g γ 2 exp(γt) 3g γ 2

Ejemplo

Cambiamos los valores de los parámetros

Conocidos estos, calculamos el campo magnético B y la resistencia de la espira R

En el script definimos las funciones que nos calculan: la velocidad v, intensidad de la corriente inducida en la espira, i y posición x del lado inferior de la espira rectangular, en función del tiempo t

m=0.1; %masa
a=0.25; %ancho de la espira rectangular
L=8e-5; %autoinducción

w0=10; %frecuencia propia
gamma=w0; %constante de amortiguamiento
R=2*L*gamma; %resistencia
B=w0*sqrt(m*L)/a; %campo magnético

%velocidad
v=@(t) 9.8*(2/gamma-(t+2/gamma).*exp(-gamma*t));
%intensidad
i=@(t) m*9.8*(1-(gamma*t+1).*exp(-gamma*t))/(B*a);
%posición
x=@(t) 9.8*(2*t/gamma+t.*exp(-gamma*t)/gamma+3*(exp(-gamma*t)-1)/gamma^2);
fplot(v,[0,3*2*pi/w0]) %velocidad
grid on
xlabel('t (s)')
ylabel('v (m/s)')
title('Críticas, velocidad')

Representamos la velocidad de la espira en función del tiempo

La velocidad alcanza rápidamente el valor límite constante, v=g/γ. Modificamos la sentencia fplot para representar la intensidad i o la posición x del lado inferior de la espira

Oscilaciones sobreamortiguadas

Cuando la resistencia es grande, puede ocurrir que γ>ω0.

La solución de la ecuación diferencial es

v=exp(γt)( Acosh(ωt)+Bsinh(ωt) )+ 2γg ω 0 2 , ω 2 = γ 2 ω 0 2 dv dt =γexp(γt)( Acosh(ωt)+Bsinh(ωt) )+ωexp(γt)( Asinh(ωt)+Bcosh(ωt) )

Donde A y B se determinan a partir de las condiciones iniciales

t=0{ 0=A+ 2γg ω 0 2 g=γA+ωB v=g{ 2γ ω 0 2 +exp(γt)( 1 ω ( 1 2 γ 2 ω 0 2 )sinh(ωt) 2γ ω 0 2 cosh(ωt) ) } dv dt =gexp(γt)( γ ω sinh(ωt)+cosh(ωt) )

Despejamos la intensidad de la ecuación del movimiento

i= 1 Ba ( mgm dv dt ) i= mg Ba ( 1exp(γt)( γ ω sinh(ωt)+cosh(ωt) ) )

Para obtener la posición x del lado inferior de la la espira integramos dx/dt=v, sabiendo que para t=0, x=0, parte del origen

x=g{ 2γ ω 0 2 t+ 1 ω ( 1 2 γ 2 ω 0 2 ) 0 t exp(γt) sinh(ωt)·dt 2γ ω 0 2 0 t exp(γt) cosh(ωt)·dt }

Integrando por partes

exp(γt)cosh(ωt)·dt= exp(γt) ω 0 2 { γcosh(ωt)+ωsinh(ωt) } exp(γt)sin(ωt)·dt= exp(γt) ω 0 2 { γsinh(ωt)+ωcosh(ωt) }

obtenemos la posición x, del lado inferior de la espira en función del tiempo.

x= g ω 0 2 { 2γ·t+exp(γt)( ( 4 γ 2 ω 0 2 1 )cosh(ωt) γ ω ( 3 4 γ 2 ω 0 2 )sinh(ωt) ) 4 γ 2 ω 0 2 +1 }

Ejemplo

Cambiamos los valores de los parámetros

Conocidos estos, calculamos el campo magnético B y la resistencia de la espira R

En el script definimos las funciones que nos calculan: la velocidad v, intensidad de la corriente inducida en la espira, i y posición x del lado inferior de la espira rectangular, en función del tiempo t

m=0.1; %masa
a=0.25; %ancho de la espira rectangular
L=8e-5; %autoinducción

w0=10; %frecuencia propia
gamma=12; %constante de amortiguamiento
R=2*L*gamma; %resistencia
B=w0*sqrt(m*L)/a; %campo magnético

w=sqrt(gamma^2-w0^2); 
%velocidad
v=@(t) 9.8*(2*gamma/w0^2+exp(-gamma*t).*((1-2*gamma^2/w0^2)*sinh(w*t)/w-
2*gamma*cosh(w*t)/w0^2));
%intensidad
i=@(t) m*9.8*(1-exp(-gamma*t).*(gamma*sinh(w*t)/w+cosh(w*t)))/(B*a);
%posición
x=@(t) 9.8*(2*gamma*t+exp(-gamma*t).*((4*gamma^2/w0^2-1)*cosh(w*t)+
gamma*(4*gamma^2/w0^2-3)*sinh(w*t)/w)-4*gamma^2/w0^2+1)/w0^2;
fplot(v,[0,3*2*pi/w0]) %velocidad
grid on
xlabel('t (s)')
ylabel('v (m/s)')
title('Sobreamortiguadas, velocidad')

Representamos la velocidad de la espira en función del tiempo

La velocidad alcanza rápidamente el valor límite constante, v=2γg/ω02. Modificamos la sentencia fplot para representar la intensidad i o la posición x del lado inferior de la espira

Estudio energético

Las energías que inteviene son:

La energía inicial, t=0, x=0, v=0, i=0, es nula, Ei=0

La energía final es la suma

E f = 1 2 m v 2 + 1 2 L i 2 mgx

La diferencia, es decir, la energía final, es la que se disipa en la resistencia

Comprobamos que la energía Ef en el instante t, es igual (en valor absoluto) a la energía disipada en la resistencia. Añadimos las siguientes líneas a cualquiera de los tres scripts (amortiguadas, críticas o sobreamortiguadas)

...%amortiguadas
t=0.5; %instante
Ef=m*v(t)^2/2+L*i(t)^2/2-m*9.8*x(t); %energía final
f=@(t) i(t).^2*R;
resis=integral(f,0,t); %disipada en la resistencia
disp([Ef,resis])
   -0.2388    0.2388

Actividades

Se introduce

El programa interactivo ha fijado los valores de los siguientes parámetros:

Conocidos estos, determinamos el campo magnético B= ω 0 a mL y la resistencia de la espira R=2

Se suguiere al lector probar las siguintes situaciones:

La espira se mueva bajo la acción de su peso y de la fuerza hacia arriba, que ejerce el campo magnético sobre la corriente inducida en el lado superior de la espira. Los vectores señalan la intensidad, vector u ^ t (color rojo), vector campo magnético B , (color azul) y el vector fuerza F m =i u ^ t × B a , (color negro). La fuerza magnética es proporcional a la intensidad de la corriente inducida i en la espira

Cuando la espira sale del campo magnético, la fuerza magnética desaparece y la espira se mueve bajo la acción de su propio peso (movimiento uniformemente acelerado). Supondremos que la intensidad decrece rápidamente debido a la resistencia del circuito

Los puntos de color rojo sobre la espira representan portadores de carga positivos y señalan el sentido de la corriente inducida en la espira. Estos puntos se mueven con una velocidad proporcional a la intensidad de la corriente. En la parte derecha, se representa la intensidad (en valor absoluto) de la corriente inducida en la espira.

En la parte superior derecha, se proporcionan los datos del tiempo t en s, la posición x en cm del lado inferior de la espira y su velocidad v en cm/s

Referencias

Saslow W. M., Electromechanical implications of Faraday’s law: A problem collection. Am. J. Phys. 55 (11) November 1987, pp. 986-993.

Robert H Romer. The motion of a superconducting loop in an inhomogeneous magnetic field: the harmonic oscillator equation in an unfamiliar setting. Eur. J. Phys. 11 (1990), pp. 103-106