# Transformer Magnetization Losses Using a Nonfiltered Voltage-Source Inverter

Authors: Rickard Ekström; Senad Apelfröjd; Mats Leijon
Journal: Advances in Power Electronics
 (2013)

---

## Abstract

Results from the magnetization of an 80 kVA power
transformer, using a directly coupled nonfiltered three-phase
voltage-source inverter (VSI), are presented. The major benefits
of this topology are reduction in switching filter size as well as
filter losses. Drawbacks include higher stress on the transformer
windings and higher transformer magnetization losses. In this
paper, the total magnetization losses are presented for different
voltage levels. The transformer has been magnetized with the
rated frequency of 50 Hz at various voltage levels. The saturation
characteristics as well as the magnetizing resistance are derived
as functions of the voltage. These are used as inputs for the
simulations. The magnetization losses have been experimentally
measured and simulated for three different DC levels. Results
from the simulations show good agreement with the experimental
results. As expected, the pulsed voltage waveforms generate
larger magnetization losses than the corresponding 50 Hz case. The losses are strongly dependent on the DC level.


---

## Body

## 1. Introduction

Many grid connected voltage source inverters (VSIs) have an associated LCL filter to eliminate the switching frequencies from the output [1–3]. If a step-up transformer is placed after the filter, its leakage inductance can reduce the size of the second inductor leg or completely replace it [4]. Alternatively, the transformer can be put directly after the inverter, with the LCL filter on the transformer secondary side. This will reduce the ampacity of the filter inductors, reducing both size and losses in the filter. The drawback of this is an increased voltage transient stress on the transformer windings that has to be taken into account in the transformer design. However, the total efficiency of the system may be improved, and this topology could be considered for grid connection of renewable energy sources. Stray capacitance between windings, which normally is neglected at steady state, may play a more central role. Semiconductor switching into magnetic components is not unusual in power electronics, for example, the flyback converter. However, these circuits are based on ferrite cores, optimized for the switching frequencies in the system. Little research is found in the area of continuous pulse modulation into laminated steel-core transformers and its associated losses.This paper presents a VSI with bipolar sinusoidal pulse-width modulation (SPWM) directly connected to a three-phase power transformer and the measured total magnetization losses. To get a reference case, the transformer is magnetized with 50 Hz at different voltage levels. The transformer saturation curve and magnetizing resistance are derived as functions of the induced voltage. A simulation model of the system is developed and uses the measured data at 50 Hz to set the internal parameters. Experiments with variations in DC voltage level and modulation index are conducted. The simulation model is verified, and the magnetization losses are evaluated for different amplitude modulations.


## 2. Theory

Starting from Faraday’s law, the Magnetic field intensity,B¯, in the transformer core is derived from the induced voltage V¯ind according to
(1)B(t)=1NA∫0t‍Vind(t)dt+Br(0),
where N is the number of winding turns, A is the average cross-sectional area of the magnetic core, and Br(0) is the remanescent flux of the magnetic core. B is coupled to the magnetic field strength H¯ according to
(2)H¯=μrμ0B¯,
where μ0=4π·10-7Vs/Am is the magnetic permeability of free space, and μr is the relative permeability constant, a nonlinear material-dependent parameter. If the magnetizing current I¯m is known, H¯ is calculated by
(3)H(t)=I¯mNℜ,
where ℜ is the magnetic reluctance of the core.The core losses of any magnetic component are distributed between eddy current losses, excess eddy current losses, and hysteresis losses. All of these have nonlinear behaviour and depend on the magnetic material properties and the operating temperature, as well as the waveform and frequency of the induced voltage. There are two established ways of calculating the core losses, either from the theoretical derivations of each loss separately [5, 6] or empirically by the Steinmetz equation [7] as follows:
(4)P=kfαBβ,
where k is a material constant and f is the frequency of the magnetic flux. The exponents α and β are nonlinear functions of the other parameters.Steinmetz equation is preferable, since the needed parameters usually are provided by the manufacturers data sheets. However, it is limited to sinusoidal waveforms and specific operating conditions. Also, having nonsinusoidal waveforms may cause large deviations in the loss calculations [8–10]. Splitting the waveform into its Fourier components and calculating the losses at the individual frequencies does not give very accurate results due to the complex interdependence of the losses. Also, it is difficult to accurately obtain all the Steinmetz parameters for a broad range of frequencies. Modifications of the Steinmetz equation are found in, for example, [11]. In [12, 13], attempts have been made to generalize the core loss calculations for an arbitrary current shape.The classical eddy current losses may be described by [14]
(5)Peddy=A8πρ(dB(t)dt)2,
while hysteresis losses are calculated by integrating the B-H-curve over one cycle as
(6)Phys=f∫tt+T‍B(t)H(t)dt.The high-frequency model of a single-phase transformer is shown for the no load case in Figure1. Rp and Rs are the winding resistances on the primary and secondary sides, respectively. Only the leakage inductance of the primary side, Lp, is included, since the secondary leakage flux is so small at no load. Publications on high-frequency modelling of power transformers mainly discuss the pulse response, as in, for example, [15].Figure 1
Transformer high-frequency equivalent circuit.The internal stray capacitances of the windings,Cp and Cs, as well as the stray capacitance between the primary and secondary winding, Cm, become very important in transient analysis. Detailed calculations on transformer stray capacitance is found in [16, 17]. The stray capacitance will be in the order of pF and is neglected throughout this work since it will not significantly alter the power dissipation.The magnetizing resistance,Rm, represents the total core losses, and the magnetizing inductance, Lm, represents the magnetic flux in the core. Assuming sinusoidal input voltage, both Rm and Lm, for rated voltage and frequency of the transformer, can be found in the manufacturers data sheets. However, these results are not valid at voltages and frequencies deviating from the transformer ratings.One common way of inverter control signal generation is by the sinusoidal pulse-width modulation (SPWM), where a sinusoidal reference signal is compared with a carrier waveform. Assuming symmetry, the SPWM output voltage of the inverter can be expressed by its Fourier components according to(7)VSPWM(t)=∑h=1N‍Vhsin(hωt),
where h is the order of the harmonic. Assuming ideal switching with no delay time between switching states, the three-phase rms voltage of the fundamental frequency is calculated as (remembering the neutral is floating)
(8)V1=VDC2ma
and may be considered independent of the switching frequency, fs. Equation (8) is only valid in the linear range (0-1) of the amplitude modulation index, ma. In SPWM, the significant output harmonics are predominantly grouped around multiples of fs [18].The rms voltage of the unfilteredVSPWM can be calculated as
(9)VSPWM,rms=∑h=1N‍Vh,rms,
where rms values of the harmonics can be found in [18].


## 3. Simulations

The Simulink model of the system is illustrated in Figure2. Each block, from left to right, is described in the following. In the model, an ideal DC voltage source is used, and the DC power drawn by the circuit is measured. The three-phase inverter uses six IGBT blocks in parallel with six free-wheeling diodes. Device parameters come from the data sheets for the devices used in the experimental set-up. The inverter is operated with SPWM control, comparing a perfect sine wave to a sawtooth carrier wave at fs=6 kHz.Figure 2
Simulink model block diagram.To limit the inrush currents at start-up and allow the system to settle, damping resistors are initially connected in series with the inverter and the transformer. When the system has settled, in approximately 100 ms, a breaker is closed to bypass the resistors in the circuit.The three-phase transformer is modelled as three linear transformers in accordance with the equivalent model shown in Figure1, but neglecting the stray capacitance of the system. The transformer parameters from Table 1 are used as input together with the saturation curve interpolated in Figure 4. A detailed description of the saturation model used is found in [19, 20]. By estimation of VSPWM,rms from (9), a value of Rm is taken from the interpolated curve in Figure 3. Rm is kept constant for each full simulation together with ma and the DC-link voltage.Table 1
Transformer characteristics.
Power rating
80 kVA
Voltage ratio
345 V Y/1 kV Y
R
p
0.0030535 pu
R
s
0.0096 pu
L
p
0.19858 pu
L
s
0.57303 pu
L
m (50 Hz, 1 kV)
500 pu
R
m (50 Hz, 1 kV)
89 puFigure 3
Magnetizing reluctance versusVrms.Figure 4
Saturation curve used in simulations.The model does not include any temperature dependence, since the power loss is considered low for the no load tests. Also, the model assumes uniformly distributed flux in the core.Initially, the transformer at no load is magnetized with the 50 Hz grid voltage at different voltage levels. The total losses are shown in Table2. From this, a value of Rm is obtained, as shown in Figure 3. A curve fit is done by piece-wise shape preserving cubic interpolation (pchip). This curve is used as simulation input, using the derived rms value of the inverter voltage from (9). Rm varies widely over the measured range due to the magnetic properties and the transformer design. The interpolated curve is considered valid within the range of measured data points, that is, for Vrms∈[01.15]. The value of the magnetizing resistance represents for both hysteresis losses and eddy current losses.Table 2
Magnetization losses at normal 50 Hz magnetization.
V
s
[
p
u
]
Experimental
P
loss
[
p
u
]
Simulated
P
loss
[
p
u
]
0.130
0.00010
0.000088
0.425
0.00040
0.00042
0.700
0.0011
0.0011
1.00
0.0037
0.0037
1.15
0.0098
0.0098To get the saturation curve of the transformer, no load phase currents and voltages were measured at various voltage levels. The induced voltage is proportional to the induced flux, and they coincide at 1 pu. The results are used to plot the upper part of the saturation curve, shown in Figure4. Symmetry is assumed for the lower quadrant, and the saturation curve is extrapolated by the use of an arctan function.


## 4. Experimental Set-Up

A single-line diagram of the experimental set-up is depicted in Figure5. The rectified grid voltage is used as DC source in the experiment. By changing tap on the tap transformer, different DC levels are obtained. The DC-capacitor bank is rated at 19 mF. Additional capacitors with small equivalent series resistance (ESR) are also used to reduce the high-frequency ripple in the DC voltage. The three-phase inverter consists of six 400GB126D IGBTs with 2SC0108T2Ax-17 Concept Driver Boards, as shown in Figure 6(a). Bipolar SPWM control is used with a sawtooth wave as the carrier waveform. The switching frequency, fs, is set to 6 kHz, and the delay time between switching transitions is set to 2.5 μs. The inverter is connected to the low side of the YY 345/1 kV 80 kVA three-phase transformer shown in Figure 6(b). The neutrals of the transformer are left floating to avoid the flow of zero sequences. The transformer ratings at rated voltage and frequency are displayed in Table 1. Base power and base voltages used for per-unit calculations are taken from this. The three-phase voltages into the magnetized transformer are measured by high-frequency differential voltage probes SI-9002, with an accuracy of ±1%. To calculate the total power losses dissipated in the transformer for no load. The three-phase currents are measured with current clamps Fluke i310s, with an accuracy of ±2%. Measurements are sampled at 400 kHz, and the power dissipation per phase is calculated. Transformer temperature was not measured during the experiment.Figure 5
Single-line diagram of the experimental set-up. Magnetization losses are measured on the transformer to the right.(a) Inverter, (b) three-phase transformer.
(a)
(b)


## 5. Results and Discussion

Transformer data at rated voltage and 50 Hz frequency are given by the manufacturer data sheet. The transformer is magnetized in no load at voltages lower than rated for 50 Hz. Measured and simulated power losses are displayed in Table2. Simulations are performed for the same set-up, using the saturation curve in Figure 4 and the value of Rm in Figure 3. The winding parameters given for rated voltage are assumed constant and independent of frequency for the entire voltage range at no load. There is very good agreement between simulations and experiments, and the magnetization losses are slightly less than 1% of the transformer rating at rated voltage.The saturation curve shows the classical trendline of an arctan function and will not affect magnetization losses significantly, as the winding resistance is small. For more accurate simulations, more measured values ofRm are required to give a better curve fit, especially at lower and higher voltage levels, where the curve is expected to stagnate. The approximated value of Rm puts a limit on the simulation range, which is kept within (0.2–0.8). The estimation of Rm is of course not exact, as it represents frequency-dependent eddy current losses and hysteresis losses. However, it will give a fair approximation of the system losses that is simple to implement.Next, the model is used to simulate the magnetization losses with the nonfiltered inverter.VSPWM,rms is estimated for the different scenarios and used to get an average value of Rm. The true value of VSPWM,rms will deviate slightly from the estimated, as the voltage will be partially filtered by the transformer. The measured magnetizing current and induced phase voltage for one phase are displayed in Figure 4. In Figure 7(a), 50 Hz magnetization is shown. In Figure 7(b) the same fundamental voltage is induced using the SPWM-controlled inverter instead. Excluding switching frequency harmonics, the harmonic content of the two currents is very similar. However, the inverter current is almost twice as large as the 50 Hz current. The pulsed current does not induce magnetic flux in the magnetic core as efficiently as the 50 Hz current. Thus it has to be higher for the same flux, and the magnetic losses will be higher.(a) Magnetizing current for 50 Hz, (b) magnetizing current for SPWM.
(a)
(b)The DC voltage is set at three different levels, and the total power loss is measured, while the amplitude modulation is varied. In Figure8(a), the DC level is set to 145 V, in Figure 8(b) is set to 335 V, and in Figure 8(c) is set to 550 V, respectively. There is rather good agreement between simulations and experiments. The limitation of the model is the estimation of Rm by the use of VSPWM,rms. For VDC=145V, the simulated power is slightly higher than the experimental. This is also true for VDC=550V. Here, the “knee” seen in the power curve is sharper in the experiments than the simulations. The possible explanation to this is the approximation of the Rm curve and using VSPWM,rms to represent the voltage input. A squared interpolation of the Rm curve could possibly give more accurate simulation results.Comparison of magnetization losses between simulations and experiments, at different values ofVDC.
(a)
V
DC
=
145 V
(b)
V
DC
=
335 V
(c)
V
DC
=
550 VFor the lower DC voltages, the power loss is approximately proportional toma, whereas it starts to increase with the square of ma for the highest voltage. This is reasonable, since the losses increase with the square of the rms voltage. Also, the magnetic core will be in its linear region for lower voltages and magnetic fluxes, whereas it will go into partial saturation, as the induced voltage is increased. In Figure 8(c), it can be seen that the losses at an induced voltage of 0.8 pu are less than half the losses at 1 pu. This indicates that a significant improvement in the magnetization losses can be made by shifting the operating point for the transformer, which is equivalent to having a larger core at nominal voltage. When the transformer is SPWM-magnetized at its rated voltage, (VDC=550V, ma=1), the magnetization losses are ten times higher than the 50 Hz magnetization at rated voltage given in Table 2.In Figure9, the overlap of the experimental results are plotted from the data in Figure 8. The same fundamental voltage output is obtained for different DC levels by varying the amplitude modulation. The corresponding losses at 50 Hz-magnetization are also plotted as a reference. Nonfiltered inverter magnetization always dissipate higher losses in the transformer compared to the sine wave. This makes perfect sense, since the transformer is optimized for this frequency. Also, both eddy current losses and hysteresis losses increase with frequency. In this experiment, the lower DC voltages are not sufficient to give the transformer full magnetization, as can be seen in Figure 9. However, there is a general trend that the losses increase for higher DC voltages with the same Vs. This is in agreement with that VSPWM,rms increase for an increased VDC while maintaining the product ma·VDC constant. At rated voltage, least magnetization losses are obtained for the lowest possible VDC, that is, when ma=1.Figure 9
Magnetization losses for the same transformer voltage with differentVDC.


## 6. Conclusions

In this paper the losses of a three-phase power transformer, magnetized by a nonfiltered inverter, have been evaluated. A simulation model has been developed based on the measured parameters of the system. The total core losses are lumped into an equivalent magnetizing resistance based on the rms value of the induced voltage, and this method shows good resemblance with the experimental measurements. For the studied cases, the magnetization losses are always higher using direct VSI magnetization compared to normal 50 Hz voltages. For VSI magnetization, the losses increase rapidly close to the transformer nominal voltage. This is because the switched voltage waveform pushes the transformer to operate more in its saturated region. At rated voltage, the magnetization losses are approximately ten times higher for SPWM than for normal 50 Hz. The results show that by shifting the transformer operating point, the magnetization losses can be significantly reduced.To minimize the magnetization losses for a given secondary voltage, the amplitude modulation should be kept as high as possible. The increased transformer losses have to be evaluated in conjunction with decreased LCL filter losses to evaluate the total system efficiency in this set-up compared with more conventional topologies.



---

*Source: 261959-2013-05-08.xml*