16. Transformer#

_images/Slide2101.PNG

In the last lecture, the ideal transformer is studied. This lecture will introduce the non-ideal part of the transformer: leakages and loss, which will be reflected as leakage inductance and resistances.

After taking the lecture, we should be able to understand the practical construction of a transformer and apply the equivalent circuit model to solve transformer problems.

16.1. Practical transformers and applications#

_images/Slide823.PNG

We first start with various applications of transformers and how a transformer is made in practice.

_images/Slide924.PNG

Transformers are all around us and come in many different shapes and sizes. In power grid we have distribution power transformer which operates at 400 V and 50 Hz in the Netherlands, as shown on the left. They are often mounted on poles, or sometimes in electrical rooms. You may have encountered them in your local community.

For the power transportation system, the power transformers used are much larger and have higher voltage ratings. What you see on the right hand side is a typical high voltage transformer operates in substations.

_images/Slide1023.PNG

Top right is the ideal model of the transformer. Bottom left shows the structure of a real power transformer. The brown poles on the transformer are called bushings, which are hollow electrical insulators that allow the electrical conductors to pass safely through the grounded tank of the transformer. There are two sets of phase bushings: the larger ones are for high voltage windings, and the tiny ones are for low voltage windings. The bushings with concentric rings on the far side are surge arresters, which are able to absorb electrical transients to protect the transformer. The concentric rings on top of them are called corona rings, which is there to prevent corona discharge.

The three phase primary and secondary windings are wounded around magnetic cores made from silicon steel laminations, which is shown on the bottom right.

The magnetic core and windings are submerged in transformer oils enclosed by the tank. The transformer oil is used to insulate, suppress corona discharge and arcing, and to serve as a coolant. The conservator tank of transformer provides adequate space to the expansion of transformer oil when it is heated up by losses or ambient.

The radiators and fans installed on the tank are used to dissipate heat generated by winding and core losses inside the transformer.

_images/Slide1124.PNG

Here it shows the manufacturing process of a power transformer. The first step is to attach copper conductors to the coil former. Depending on the type of transformers, the coil may have a helical shape or a disk shape.

_images/Slide1224.PNG

The copper conductors used are not necessarily to have round shapes. For example, flat copper condutors are used to construct the coil on this slide.

_images/Slide1325.PNG

Then the silicon steel laminations are cut or punched into different widths to make the magnetic core.

_images/Slide1425.PNG

The punched laminations are then stacked together and fixed by mechanical structures to form the magnetic core.

_images/Slide1526.PNG

Then the coils will be mounted to the magnetic core to form the primary and secondary three phase windings, and extra insulation will be added so that the conductors are well insulated from the supporting structure and the magnetic core.

_images/Slide1626.PNG

Transformers are also extensively used in power electronics. These transformers usually operate at high frequencies: from several kHz to MHz. High frequency is used because to convert the same power, a higher frequency helps reduce the size of transformer. For high frequency transformers, pre-made standard ferrite cores are used.

There are some examples as shown on the right hand side of the slide.

_images/Slide1725.PNG

Transformer is also used for contactless power transfer. For these transformers, “air core” is used instead of ferromagnetic core, so the magnetic flux flows in air and is not regulated by the magnetic core. As a consequence of that, the coupling factor between the two coils is very low. To efficiently transfer power through the two loosely coupled coils, compensation capacitors are usually added to create a resonance power conversion. This will be further elaborated in the module 1 of the course lab.

16.2. Mutual inductances in transformers#

_images/Slide1825.PNG

Let us now deal with the non-ideal behaviours of the transformer.

_images/Slide1923.PNG

In practice, not all flux can be short-circuited by the magnetic core: there is always part of the flux flows to surroundings. Therefore, there would be part of the flux generated from one coil does not couple with the secondary coil.

_images/Slide2022.PNG

As we can see here, the flux lines can be divided into three paths:

  • \(\phi_{11}\): those are linked with coil 1 through surrounding air, but do not go to coil 2 through the magnetic core;

  • \(\phi_{12}\): those link both coil 1 and coil 2 through the magnetic core;

  • \(\phi_{22}\): those are linked with coil 2 through surrounding air, but do not go to coil 1 through the magnetic core.

For the three magnetic flux paths, we can define their permeance respectively: \(\mathcal{P}_{11}\), \(\mathcal{P}_{12}\) and \(\mathcal{P}_{22}\). If the total MMF given by the coils is \(Ni\), we have

\[\begin{split} \begin{align} \phi_{11} &= Ni\mathcal{P}_{11}, \\ \phi_{12} &= Ni\mathcal{P}_{12}, \\ \phi_{22} &= Ni\mathcal{P}_{22}. \end{align} \end{split}\]
_images/Slide2126.PNG

Here we used the magnetic principles we learnt in the last lecture, as shown in this flow chart. Let us continue to apply this flow chart to derive the inductances of the coils based on the permeance.

_images/Slide2223.PNG

To study how to quantify the coupling between the two coils, let us do a thought experiment. Let current flow in coil 1 to be \(i_1\), and keep the second coil open, so \(i_2 = 0\). The total flux linked with coil 1 is composed of the flux links coil 2 through the magnetic core \(\phi_{21}\) and the part only links with coil 1 itself \(\phi_{11}\).

\[ \phi_1 = \phi_{11} + \phi_{21}. \]

From the magnetic circuit principle we know

\[ \phi_1 = \mathcal{P}_1N_1i_1 = \mathcal{P}_{11}N_1i_1 + \mathcal{P}_{21}N_1i_1. \]

The permeance \(\mathcal{P}_{21}\) is corresponding to the main path of the flux, and \(\mathcal{P}_{11}\) is corresponding to the leakage path of coil 1.

_images/Slide2322.PNG

Then the flux linkages of the two coils are

\[\begin{split} \begin{align} \lambda_1 &= N_1^2i_1\mathcal{P}_1, \\ \lambda_{21} &= N_2N_1i_1\mathcal{P}_{21}, \end{align} \end{split}\]

where \(\lambda_{21}\) represents the flux linkage linked by coil 2 caused by current in coil 1. The coil 2 itself does not generate flux since it is open at the moment. \(\mathcal{P}_{21}\) is the permeance of the magnetic circuit path from coil 1 to coil 2.

Then based on Faraday’s law, we are able to derive the voltage and the inductance of coil 1, which are

\[\begin{split} \begin{align} v_1 &= {\mathrm{d} \lambda_1 \over \mathrm{d} t} \\ &= N_1 {\mathrm{d} \phi_1 \over \mathrm{d} t} = N_1 \frac{\mathrm{d}}{\mathrm{d} t}(\phi_{11}+\phi_{21}) \\ &= N_1^2(\mathcal{P}_{11} + \mathcal{P}_{21})\frac{\mathrm{d}i_1}{\mathrm{d} t} \\ &= N_1^2\mathcal{P}_{1}\frac{\mathrm{d}i_1}{\mathrm{d} t} = L_1 \frac{\mathrm{d}i_1}{\mathrm{d} t} \end{align} \end{split}\]
\[ L_1 = N_1^2\mathcal{P}_1. \]

\(L_1\) is called self inductnace coil 1 since it corresponds to all flux linked with coil 1 itself.

_images/Slide2420.PNG

Similarly, according to Faraday’s law, there is also voltage induced in coil 2.

\[\begin{split} \begin{align} v_2 &= \frac{\mathrm{d}\lambda_2}{\mathrm{d} t} \\ &= N_2\frac{\mathrm{d}\phi_{21}}{\mathrm{d} t} = N_2 \mathcal{P}_{21} N_1 \frac{\mathrm{d} i_1}{\mathrm{d} t} \\ &= M_{21} {\mathrm{d} i_1 \over \mathrm{d} t}. \end{align} \end{split}\]

Here we define \(M_{21} = N_1N_2\mathcal{P}_{21}\) as the mutual inductance from coil 1 to coil 2, since it represents the capability of coil 1 induces voltage in coil 2.

_images/Slide2517.PNG

Similarly, by setting coil 1 open, and inject a current \(i_2\) in coil 2, following the current and voltage definitions given in the figure, the voltages in coil 1 and coil 2 can be solved as

\[\begin{split} \newcommand{\diff}{\mathrm{d}} \begin{align} v_1 &= \frac{\diff \lambda_1}{\diff t} = -N_1 \frac{\diff \phi_{12}}{\diff t} \\ &= -N_1N_2\mathcal{\mathcal{P}_{12}}\frac{\diff i_2}{\diff t}. \end{align} \end{split}\]
\[\begin{split} \newcommand{\diff}{\mathrm{d}} \begin{align} v_2 &= \frac{\diff \lambda_2}{\diff t} = -N_2 (\mathcal{P}_{22}+\mathcal{P}_{12})\frac{\diff i_2}{\diff t} \\ &= -N_2^2\mathcal{P_{2}}\frac{\diff i_2}{\diff t} = -L_2 \frac{\diff i_2}{\diff t}. \end{align} \end{split}\]

Here we define \(M_{12} = N_1N_2\mathcal{P}_{12}\) as the mutual inductace from coil 2 to coil 1. \(L_2\) is the self inductance of coil 2.

Apparently the flux path from coil 1 to coil 2 are exactly the same as that from coil 2 to coil 1, so the two mutual inductance is the same

\[ M_{12} = M_{21} = M. \]
_images/Slide2615.PNG

In the above two thought experiments, one of the coils is kept open. Now let’s check a more general case, where both \(i_1\) and \(i_2\) are non zero.

Following the notation given in the figure on the slide, the flux linkages of the two coils should include the contributions from coil itself and the other coil

\[\begin{split} \lambda_1 &= N_1(\phi_{1} - \phi_{12}), \\ \lambda_2 &= -N_2(\phi_{2} - \phi_{21}). \end{split}\]

Then the voltage equations can be written as

\[\begin{split} \newcommand{\diff}{\mathrm{d}} \begin{align} v_1 &= {\diff \lambda_1 \over \diff t} \\ &= N_1 \left( {\diff \phi_1 \over \diff t} - {\diff \phi_{12} \over \diff t}\right) \\ &= L_1 \frac{\diff i_1}{\diff t} - M\frac{\diff i_2}{\diff t}. \end{align} \end{split}\]
\[\begin{split} \newcommand{\diff}{\mathrm{d}} \begin{align} v_2 &= {\diff \lambda_2 \over \diff t} \\ &= -N_2 \left( {\diff \phi_2 \over \diff t} - {\diff \phi_{21} \over \diff t}\right) \\ &= -L_2 \frac{\diff i_1}{\diff t} + M\frac{\diff i_2}{\diff t}. \end{align} \end{split}\]

or as the matrix form shown on the slide.

Attention

There are many different definitions for the current and voltage directions and many textbooks will choose a different secondary voltage definition such that the signs in the second row of the matrix reverses. Here we choose to follow the definition shown on the slide so that is it consistent with the standard definitions of an ideal transformer.

_images/Slide2715.PNG

The two voltage equations or their matrix form do not give an intuitive feeling about how well the two coils are coupled. Hereby we introduce the concept of coupling factor.

The product of the two self inductances expressed by permeance is shown on the slide. Keeping in mind that the magnetic circuit remains the same, so \(\mathcal{P}_{12} = \mathcal{P}_{21}\). By taking \(\mathcal{P}_{12}\) out of the parentheses, we have

\[\begin{split} \begin{align} L_1L_2 &= M^2\left (1+\frac{\mathcal{P}_{11}}{\mathcal{P}_{12}} \right) \left ( 1+ \frac{\mathcal{P}_{22}}{\mathcal{P}_{12}} \right) \\ &= M^2 \frac{1}{k^2}. \end{align} \end{split}\]

Here the coefficient of coupling, or coupling coefficient \(k\) is defined as \(1/\sqrt{(1+\mathcal{P}_{11}/\mathcal{P}_{12})(1+\mathcal{P}_{22}/\mathcal{P}_{12})}\). Therefore we have \(M = k\sqrt{L_1L_2}\).

_images/Slide2814.PNG

\(k = 0\) means there is no common flux between the two coils and the mutual inductance is zero. \(k = 1\) means the coupling between the two coils is perfect and the leakage flux \(\phi_{11} = \phi_{22} = 0\). The allowable range for the coefficient of coupling is \(0 \leq k \leq 1\). In practical systems it is impossible to obtain \(k=1\). However, in certain cases, such as power transformers, where the intention is to maximize the coupling, it is possible to achieve a coupling coefficient in excess of 0.99. In contactless power transfer applications, however, \(k\) can be as low as several percents.

_images/Slide2911.PNG

Using the voltage equations and mutual inductance is not always straightforward for circuit analysis, since we have to deal with two equations. It is handy for analysis of coupled circuits if we introduce the leakage and magnetising inductances.

_images/Slide3011.PNG

We know the self inductance of coil 1 can be written as the sum of two parts: the one corresponds to the flux which only links coil 1 itself \(N_1^2\mathcal{P}_{11}\), and the one corresponds to the flux couples with both coils \(N_2^2\mathcal{P}_{21}\).

THe former one represents the flux which is leaked to the surroundings, so it is defined as the leakage inductance \(L_{L1}\), and the second one is associated with the coupled flux which magnetise the magnetic core, so it is defined as the magnetising inductance \(L_m\).

This way we are able to divide the coil 1 inductance into two inductances: magnetising inductance and leakage inductance.

Similar approach can be applied to coil 2.

_images/Slide3115.PNG

In practice, it is difficult to measure the permeance \(\mathcal{P}_{11}\) or \(\mathcal{P}_{21}\), so it is difficult to obtain the value of \(L_m\) and \(L_{L1}\) from permeances. However, recalling the definition of the mutual inductance and coupling coefficient, we are able to obtain the relative relationship between the inductances.

For the primary side (coil 1), the magnetising inductance is

\[\begin{split} \begin{align} L_{mp} &= N_1^2\mathcal{P}_{21} = \frac{N_1}{N_2}N_1N_2\mathcal{P}_{21} \\ &= \frac{N_1}{N_2} M = \frac{N_1}{N_2}k\sqrt{L_1L2}. \end{align} \end{split}\]

Using the same approach, if magnetising inductance is placed on the secondary side, we have the secondary magnetising inductance as

\[\begin{split} \begin{align} L_{ms} &= N_2^2\mathcal{P}_{21} = \frac{N_2}{N_1}N_1N_2\mathcal{P}_{21} \\ &= \frac{N_2}{N_1} M = \frac{N_2}{N_1}k\sqrt{L_1L2}. \end{align} \end{split}\]

You may notice that

\[ \frac{L_{mp}}{L_{ms}} = \frac{N_1^2}{N_2^2}, \]

which is corresponding to the impedance transfer ratio when the magnetising inductance is referred from the secondary side to the primary side.

_images/Slide3211.PNG

Based on the defined magnetising inductance and the leakage inductances, we will introduce the equivalent circuit of the non ideal transformer.

_images/Slide3311.PNG

Now we can approximate a non ideal transformer by adding the magnetising inductance and the leakage inductances to an ideal transformer, as shown on the slide.

If we assume the coupling coefficient \(k\) is high enough, since \(L=N^2\mathcal{P}\), the transformer ratio of the ideal transformer can be approximated as the square root of the self inductance ratio

\[ a \approx \frac{N_1}{N_2} = \sqrt{\frac{L_1}{L_2}}. \]
_images/Slide347.PNG

Then the magnetising inductance is calculated as

\[\begin{split} \begin{align} L_{mp} &= \frac{N_1}{N_2}M = \frac{N_1}{N_2}k\sqrt{L_1L2} \\ &= \sqrt{\frac{L_1}{L_2}}k\sqrt{L_1L2} = kL_1. \end{align} \end{split}\]

The primary leakage inductance is therefore

\[ L_{L1} = L_1 - L_{mp} = (1-k)L_1. \]

Likewise, the secondary side inductance is solve

\[ L_{L1} = L_2 - L_{ms} = (1-k)L_2. \]

16.3. Example#

_images/Slide356.PNG

Then let us use this example to practice problem solving with equivalent circuit. Please try to work it out yourself and click the block below for the detailed answer. The Python code below also shows the calculation procedure.

Hide code cell source
%reset -f

import numpy as np
from sympy import *
import cmath
from IPython.display import display, Markdown, Math, Latex

k  = 0.9
L1 = 400.0e-3
L2 = 100.0e-3
Vp = 200.0
R1 = 1.0
R2 = 100.0
f  = 50.0
omega = 2*np.pi*f

M  = k*np.sqrt(L1*L2)

a  = sqrt(L1/L2)

# first use the coupled inductors method

i_1, i_2 = symbols('i_1, i_2')
eq1 = Vp-i_1*R1-1j*omega*(i_1*L1-i_2*M)
eq2 = 1j*omega*M*i_1 - 1j*omega*L2*i_2 - R2*i_2

sols = solve([eq1,eq2], [i_1,i_2])

I_1  = sols[i_1].evalf()
I_2  = sols[i_2].evalf()

print('From the coupled inductors method:')
display(Math('$\mathbf{{I_1}}={:.2f}\\angle{:.2f}^\circ$'.format(abs(I_1), cmath.phase(I_1)/np.pi*180)))
display(Math('$\mathbf{{I_2}}={:.2f}\\angle{:.2f}^\circ$'.format(abs(I_2), cmath.phase(I_2)/np.pi*180)))

# then use the equivalent circuit method

Lm  = k*L1
LL1 = (1-k)*L1
LL2 = (1-k)*L2

Z2  = 1j*omega*LL2+R2

# refer secondary impedance to primary side
Z2_prime = Z2*a**2

Zm  = 1j*omega*Lm
Z1  = R1 + 1j*omega*LL1

# total impedance
Ztot = Z1 + 1/(1/Zm+1/Z2_prime)
I1   = Vp/Ztot

# voltage across Lm
Vm   = (1/(1/Zm+1/Z2_prime))/Ztot*Vp
# referred to secondary side
V2   = Vm/a
I2  = V2/Z2

print('From the equivalent circuit method:')
display(Math('$\mathbf{{I_1}}={:.2f}\\angle{:.2f}^\circ$'.format(abs(I1), cmath.phase(I1)/np.pi*180)))
display(Math('$\mathbf{{I_2}}={:.2f}\\angle{:.2f}^\circ$'.format(abs(I2), cmath.phase(I2)/np.pi*180)))
From the coupled inductors method:
\[\displaystyle \mathbf{I_1}=1.66\angle-75.51^\circ\]
\[\displaystyle \mathbf{I_2}=0.90\angle-2.95^\circ\]
From the equivalent circuit method:
\[\displaystyle \mathbf{I_1}=1.66\angle-75.51^\circ\]
\[\displaystyle \mathbf{I_2}=0.90\angle-2.95^\circ\]

In the code above, we solved the problem using two approaches: 1) analyse it as coupled inductor and solve the voltage equations, 2) use the equivalent circuit to simplify the circuit.

Apparently both methods give use the same results. The equivalent circuit is easier for hand calculation since it does not require linear equation solving.

16.4. Transformer parameter measurements#

_images/Slide394.PNG

The circuit elements in the equivalent circuit are not directly accessible, however, it is possible to obtain their values indirectly by doing some tests.

First an open circuit test is carried out. We supply the primary side with rate voltage \(V_{1o}\), then measure the input current \(I_{1o}\) and secondary side open circuit voltage \(V_{2o}\). The open circuit reactance \( X_{1o}\) can be obtained from \(V_{1o}\) and \(I_{1o}\). From them the following equations are obtained

\[ a = \frac{V_{1o}}{V_{2o}}. \]
\[ \omega(L_{L1} + L_m) = X_{1o}. \]

Then a short circuit test is carried out. We short circuit the secondary side, and supply primary side with a rated current \(I_{1s}\). Then primary side voltage \(V_{Is}\) is measured. The short circuit reactance \(X_{1s}\) can be obtained from \(V_{1s}\) and \(I_{1s}\). If we assume \(k\) is high enough, we can reasonably neglecting current going through \(L_m\), so we have

\[ \omega(L_{L1} + a^2L_{L2}) = X_{1s}. \]

There are three equations but four unknows, we need to make further assumptions to make it solvable. One approach is to assume \(\mathcal{P}_{11} = \mathcal{P}_{22}\), i.e., the leakage paths have the same permeance on both sides, so we have

\[ L_{L1} = a^2 L_{L2}. \]

Then we are able to solve all four parameters in the equivalent circuits from measurement results.

The method described above is a simplified version based on the standard approach. In practice, the measurements are more complicated since we have to deal with more non-ideal factors including winding resistances and core losses etc.