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
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
Derivando con respecto del tiempo
Aplicamos la ley de Faraday al secundario
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
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
El coeficiente de inducción mutua y de autoinducción son
Derivando con respecto del tiempo
Aplicamos la ley de Faraday y la ley de Ohm al primario (véase el apartado titulado 'Ecuación del circuito RL'
El campo eléctrico total, mueve las cargas en el conductor produciendo una intensidad i1. La primera integral es el producto de la intensidad por la resistencia, i1R1.
es el campo eléctrico producido por la diferencia de potencial Us de la batería, la segunda integral, es Us
Para N1 espiras
El resultado es
Sustituyendo en la primera ecuación, el primer término
El comportamiento del transformador está descrito por dos ecuaciones, una para el circuito primario y otra, para el secundario
Gráficas
Los datos del transformador son los siguientes:
- Primario: autoinducción L1=3.69 H, resistencia, R1=4.5 Ω, número de espiras N1=900
- Secundario: autoinducción L2=0.41 H, resistencia, R2=1.4 Ω, número de espiras N2=300
- El primario está conectado a una fem alterna de 60 Hz, ω=2π·60 rad/s, de amplitud Us=50.2 V
- Coeficiente de inducción mutua, M=1.176 H
Las ecuaciones de los circuitos en forma compleja son
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
Multiplicamos la primera por i1, la segunda por i2 y sumamos
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
El primer término es positivo, el segundo será positivo o nulo
Se define el coeficiente de acoplamiento como el cociente
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
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
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
Introducimos en la ecuación diferencial
Calculamos la amplitud I0 y el desfase φ
Sustituimos a1, a2, a3 y a4 por sus valores
En un transformador se cumple aproximadamente que M2≈L1L2
Cuando la frecuencia angular ω se hace grande I0 tiende hacia el valor
Gráficas
Los datos del transformador son los siguientes:
- Primario: autoinducción L1=3.69 H, resistencia, R1=4.5 Ω
- Secundario: autoinducción L2=0.41 H, resistencia, R2=1.4+20 Ω
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