El transformador (II)

En esta apartado, vamos a proporcionar una descripción más detallada del funcionamiento del transformador

En la figura, se muestra un transformador con N1 espiras en el primario y N2 en el secundario. φ1 es el flujo a través de una espira producido por la corriente i1 en el primario y φ2 es el flujo a través de una espira producido por la corriente i2 en el secundario.

Secundario

El flujo a través de una espira del secundario es la diferencia φ=φ2-φ1. El flujo total en el secundario es

N 2 φ= N 2 φ 2 N 2 φ 1

El coeficiente de inducción mutua M es el cociente entre el flujo a través del secundario (producido por el primario) y la corriente en el primario. El coeficiente de autoinducción L2 es el cociente del flujo propio a través de un circuito de la derecha y la intensidad que circula por el mismo

M= N 2 φ 1 i 1 , L 2 = N 2 φ 2 i 2 N 2 φ= L 2 i 2 M i 1

Derivando con respecto del tiempo

N 2 dφ dt = L 2 d i 2 dt M d i 1 dt

Aplicamos la ley de Faraday al secundario

N 2 dφ dt = A B E i · dl = N 2 E i · dl = i 2 ( R 2 + R L )

R2 es la resistencia del secundario y RL es la resistencia de un circuito externo (dispositivo) conectado a los extremos A y B del secundario

La ecuación del secundario es

i 2 ( R 2 + R L )= L 2 d i 2 dt M d i 1 dt M d i 1 dt = L 2 d i 2 dt + i 2 ( R 2 + R L )

Primario

Repetimos los pasos para obtener la ecuación del primario, intercambiando los papeles de primario y secundario

El flujo a través de una espira del primario es la diferencia φ=φ1-φ2. El flujo total en el primario es

N 1 φ= N 1 φ 1 N 1 φ 2

El coeficiente de inducción mutua y de autoinducción son

M= N 1 φ 2 i 2 , L 1 = N 1 φ 1 i 1 N 1 φ= L 1 i 1 M i 2

Derivando con respecto del tiempo

N 1 dφ dt = L 1 d i 1 dt M d i 2 dt

Aplicamos la ley de Faraday y la ley de Ohm al primario (véase el apartado titulado 'Ecuación del circuito RL'

E = E s + E i A B E · dl = A B E s · dl + A B E i · dl

El resultado es

i 1 R 1 = U s N 1 dφ dt

Sustituyendo en la primera ecuación, el primer término

N 1 dφ dt = U s i 1 R 1 L 1 d i 1 dt M d i 2 dt = U s i 1 R 1 U s = i 1 R 1 + L 1 d i 1 dt M d i 2 dt

El comportamiento del transformador está descrito por dos ecuaciones, una para el circuito primario y otra, para el secundario

{ i 1 R 1 + L 1 d i 1 dt M d i 2 dt = U s i 2 ( R 2 + R L )+ L 2 d i 2 dt M d i 1 dt =0

Gráficas

Los datos del transformador son los siguientes:

Las ecuaciones de los circuitos en forma compleja son

{ I 1 R 1 +iω L 1 I 1 iωM I 2 = U s I 2 ( R 2 + R L )+iω L 2 I 2 iωM I 1 =0 I 1 = R 2 + R L +iω L 2 ( R 1 +iω L 1 )( R 2 + R L +iω L 2 )+ ω 2 M 2 U s I 2 = iωM ( R 1 +iω L 1 )( R 2 + R L +iω L 2 )+ ω 2 M 2 U s

I1 e I2 son números complejos que contienen la información acerca de la amplitud y la fase

L1=3.69; %primario
N1=900;
R1=4.5;
L2=0.41; %secundario
N2=300;
R2=1.4;
M=1.176; %inducción mutua
Us=50.2; %fuente de corriente alterna
w=2*pi*60;

i1=@(x)  Us*(1i*w*L2+R2+x)./((R1+1i*w*L1)*(1i*w*L2+R2+x)+w^2*M^2);
i2=@(x) 1i*Us*w*M./((R1+1i*w*L1)*(1i*w*L2+R2+x)+w^2*M^2);
I1=@(x) abs(i1(x));
I2=@(x) abs(i2(x));
figure %intensidades I1 e I2
hold on
fplot(I1,[0,100])
fplot(I2,[0,100])
hold off
grid on
xlabel('R_L')
ylabel('I_1, I_2')
legend('I_1', 'I_2', 'location', 'best')
title('Intensidades en el primario y secundario')

figure %diferencia de potencial en el secundario
V=@(x) abs(i2(x).*x);
fplot(V,[0,100])
grid on
xlabel('R_L')
ylabel('V')
title('ddp en el secundario')

figure %energía por unidad de tiempo
ang1=@(x) angle(i1(x));
producida=@(x) Us*I1(x).*cos(ang1(x)); 
consumida=@(x) (I2(x).^2).*x;
hold on
yyaxis left
fplot(producida,[0,100])
fplot(consumida,[0,100])
ylabel('Potencia')
yyaxis right
rendimiento=@(x) consumida(x)./producida(x);
fplot(rendimiento,[0,100])
ylabel('Rendimiento')
hold off
grid on
xlabel('R_L')
title('Potencia, producida, consumida y redimiento')

Representamos las amplitudes de las intensidades I1 e I2 en función de la resistencia RLdel circuito (dispositivo) conectado al secundario

Representamos la diferencia de potencial en el secundario I2RL en función de la resistencia RL

Representmos la potencia (energía por unidad de tiempo) producida en el circuito primario, la consumida en el circuito secundario y el redimiento del transformador

Representamos el cociente I2/I1, V1/V2 y la recta N1/N2 en función de RL.

L1=3.69; %primario
N1=900;
R1=4.5;
L2=0.41; %secundario
N2=300;
R2=1.4;
M=1.176; %inducción mutua
Us=50.2; %fuente de corriente alterna
w=2*pi*60;

i1=@(x)  Us*(1i*w*L2+R2+x)./((R1+1i*w*L1)*(1i*w*L2+R2+x)+w^2*M^2);
i2=@(x) 1i*Us*w*M./((R1+1i*w*L1)*(1i*w*L2+R2+x)+w^2*M^2);
I_r=@(x) abs(i2(x)./i1(x));
hold on
fplot(I_r,[0,100])
V_r=@(x) Us./abs(i2(x).*x);
fplot(V_r,[0,100])
line([0,100],[N1/N2, N1/N2],'lineStyle','--')
grid on
ylim([0,5])
xlabel('R_L')
ylabel('Cocientes')
title('Cocientes, I2/I1, V1/V2, N1/N2')

Coeficiente de acoplamiento

Supongamos que el circuito primario y el secundario no tienen resistencias. Las ecuaciones que los describen son

{ L 1 d i 1 dt M d i 2 dt = U s L 2 d i 2 dt M d i 1 dt =0

Multiplicamos la primera por i1, la segunda por i2 y sumamos

L 1 i 1 d i 1 dt + L 2 i 2 d i 2 dt M i 1 d i 2 dt M i 2 d i 1 dt = U s i 1 d dt ( 1 2 L 1 i 1 2 + 1 2 L 2 i 2 2 M i 1 i 2 )= U s i 1

Entre paréntesis, los dos primeros términos son las energías almacenadas por las dos bobinas Li2/2, el tercer término corresponde al acoplamiento. La cantidad entre parétesis es positiva o nula

1 2 L 1 i 1 2 + 1 2 L 2 i 2 2 M i 1 i 2 0 1 2 ( L 1 i 1 L 2 i 2 ) 2 + i 1 i 2 ( L 1 L 2 M )0

El primer término es positivo, el segundo será positivo o nulo

L 1 L 2 M0

Se define el coeficiente de acoplamiento como el cociente

k= M L 1 L 2 ,0k1

k=1, corresponde a un acoplamiento ideal y k=0, a una falta de acoplamiento entre los dos circuitos

Se conecta el primario a una fuente alterna

Para simplificar las ecuaciones, designamos R2 a la resistencia total del secundario, la suma de la resistencia externa RL y de la propia de la bobina que hasta ahora hemos denominado R2

El sistema de dos ecuaciones es

i 1 R 1 + L 1 d i 1 dt M d i 2 dt = V 0 sin( ωt ) i 2 R 2 + L 2 d i 2 dt M d i 1 dt =0

Derivamos la primera respecto del tiempo y utilizamos la segunda y su derivada respecto de tiempo, para eliminar las derivadas de i1 en la primera

R 1 d i 1 dt + L 1 d 2 i 1 d t 2 M d 2 i 2 d t 2 = V 0 ωcos( ωt ) { M d i 1 dt = i 2 R 2 + L 2 d i 2 dt M d 2 i 1 d t 2 = R 2 d i 2 dt + L 2 d 2 i 2 d t 2 ( L 1 L 2 M 2 ) d 2 i 2 d t 2 +( R 1 L 2 + R 2 L 1 ) d i 2 dt + R 1 R 2 i 2 =M V 0 ωcos( ωt ) a 1 d 2 i 2 d t 2 + a 2 d i 2 dt + a 3 i 2 = a 4 cos( ωt ),{ a 1 = L 1 L 2 M 2 a 2 = R 1 L 2 + R 2 L 1 a 3 = R 1 R 2 a 4 =M V 0 ω

La solución de la ecuación diferencial homogénea describe el estado transitorio y la solución particular, el estado estacionario, que es el que nos interesa

i 2 (t)=Asin( ωt )+Bcos( ωt ) d i 2 dt =ω( Acos( ωt )Bsin( ωt ) ) d 2 i 2 d t 2 = ω 2 ( Asin( ωt )+Bcos( ωt ) )

Introducimos en la ecuación diferencial

a 1 ω 2 ( Asin( ωt )+Bcos( ωt ) )+ a 2 ω( Acos( ωt )Bsin( ωt ) )+ a 3 ( Asin( ωt )+Bcos( ωt ) )= a 4 cos( ωt ) ( A a 1 ω 2 B a 2 ω+ a 3 A )sin( ωt )+( a 1 ω 2 B+ a 2 ωA+ a 3 B a 4 )cos( ωt )=0 { A a 1 ω 2 B a 2 ω+ a 3 A=0 a 1 ω 2 B+ a 2 ωA+ a 3 B= a 4 A= a 4 a 2 ω ( a 2 ω ) 2 + ( a 3 a 1 ω 2 ) 2 ,B= a 4 ( a 3 a 1 ω 2 ) ( a 2 ω ) 2 + ( a 3 a 1 ω 2 ) 2

Calculamos la amplitud I0 y el desfase φ

i 2 (t)= I 0 sin( ωt+φ )= I 0 sin( ωt )cosφ+ I 0 cos( ωt )sinφ, { A= I 0 cosφ B= I 0 sinφ ,{ I 0 = A 2 + B 2 tanφ= B A I 0 = a 4 a 1 2 ω 4 +( a 2 2 2 a 3 a 1 ) ω 2 + a 3 2 ( a 2 ω ) 2 + ( a 3 a 1 ω 2 ) 2 ,tanφ= a 3 a 1 ω 2 a 2 ω

Sustituimos a1, a2, a3 y a4 por sus valores

I 0 = M V 0 ω ( L 1 L 2 M 2 ) 2 ω 4 +( ( R 1 L 2 + R 2 L 1 ) 2 2 R 1 R 2 ( L 1 L 2 M 2 ) ) ω 2 + R 1 2 R 2 2 ( ( R 1 L 2 + R 2 L 1 )ω ) 2 + ( R 1 R 2 ( L 1 L 2 M 2 ) ω 2 ) 2 = M V 0 ω ( L 1 L 2 M 2 ) 2 ω 4 +( R 1 2 L 2 2 + R 2 2 L 1 2 +2 R 1 R 2 M 2 ) ω 2 + R 1 2 R 2 2 tanφ= R 1 R 2 ( L 1 L 2 M 2 ) ω 2 ( R 1 L 2 + R 2 L 1 )ω

En un transformador se cumple aproximadamente que M2L1L2

I 0 V 0 ω L 1 L 2 ( R 1 L 2 + R 2 L 1 ) 2 ω 2 + R 1 2 R 2 2 tanφ R 1 R 2 ( R 1 L 2 + R 2 L 1 )ω

Cuando la frecuencia angular ω se hace grande I0 tiende hacia el valor

lim ω I 0 = V 0 L 1 L 2 R 1 L 2 + R 2 L 1

Gráficas

Los datos del transformador son los siguientes:

Representamos la amplitud de la intensidad en el secundario, I0 en función de la frecuencia angular ω de la fuente alterna conectada al primario

L1=3.69; %primario
R1=4.5;
L2=0.41; %secundario
R2=1.4+20; %resistencia total en el secundario
V0=50.2; %fuente de corriente alterna

I2=@(x) V0*x*sqrt(L1*L2)./sqrt((R1*L2+R2*L1)^2*x.^2+R1^2*R2^2);
fplot(I2,[0,10])
If= V0*sqrt(L1*L2)/(R1*L2+R2*L1);
line([0,10],[If,If],'lineStyle','--')
grid on
ylim([0,1])
xlabel('\omega')
ylabel('I_2')
title('Intensidad en el secundario')

La amplitud tiende hacia un valor constante cuando se incrementa la frecuencia ω. La amplitud final disminuye con la resistencia R2

Representamos la fase de la intensidad en el secundario, φ en función de la frecuencia angular ω de la fuente alterna conectada al primario

L1=3.69; %primario
R1=4.5;
L2=0.41; %secundario
R2=1.4+20;  %resistencia total en el secundario
V0=50.2; %fuente de corriente alterna

phi=@(x) atan(R1*R2./((R1*L2+R2*L1)*x));
fplot(phi,[eps,10])
grid on
set(gca,'YTick',0:pi/4:pi)
set(gca,'YTickLabel',{'0','\pi/4','\pi/2','3\pi/4','\pi'})
xlabel('\omega')
ylabel('\phi')
title('Desfase en el secundario')

Referencias

Fabio S El Hage, Fabio P B Vieira, Rodrigo Carareto. Deriving and validating an electrodynamic model of a voltage transformer from scratch. Eur. J. Phys. 42 (2021) 035207

J. Pierrus. Solved Problems in Classical Electromagnetism. Analytical and numerical solutions with comments. Oxford University Press (2018). Question 6.23, pp. 327-328.

B. Evstatiev, Mutually coupled inductors. Lecture notes in Theory of electrical engineering