^{1}

^{*}

^{1}

^{1}

The paper describes a simulated experiment that focuses on the numeric computation of magnetic loss in the laminated core of a single-phase power transformer. The students’ laboratory work is part of the library of experiments of the Electrical Machines virtual laboratory and makes use of the two-dimensional open-access electromagnetic field analysis software Finite Element Method Magnetics. The idea of the simulated exercise is to demonstrate how the magnetic loss caused by time-varying excitations affects the magnetic permeability,
*μ*, of the laminated core and the terminal quantities of the energizing winding. A parametric analysis employing different values for the electrical conductivity and maximum hysteresis-induced angle of the laminated material yields five different field problems with increasing magnetic loss. Electric circuits characterized by the (
*I-V*) operating point and reflected impedance of the energizing winding provide the information required to compute the changes in real power Δ
*P*, reactive power Δ
*Q* and magnetically stored energy Δ
*W*
_{m} between successive problems characterized by increasing magnetic loss. The concept of reflected impedance helps to explain the physical meaning of the changes in power dissipation and energy storage in the laminated core.

Most power transformers employ thin laminations in their ferromagnetic cores in the attempt to reduce eddy current effects. The material used for transformer cores is grain-oriented silicon-iron. The laminated material is hot-rolled then cold worked twice, followed by an anneal to improve the grain orientation and increasing the permeability along the rolling direction [

In two-dimensional (2-D) finite-element analysis, transformers are modeled as planar structures with longitudinal symmetry. In the analysis, both driving currents and induced currents flow in the longitudinal or z-direction. A detailed discussion on the numeric modeling of magnetic devices with longitudinal symmetry has been recently published and appears in [

Finite Element Method Magnetics (FEMM) is a suite of programs for solving low-frequency electromagnetic problems on two-dimensional planar and axisymmetric domains. The program currently addresses linear and non-linear magnetostatic problems, linear and non-linear time-harmonic magnetic problems, linear electrostatic problems, and steady-state heat flow problems. Its main advantage is the low computational cost and its ease of use. Despite these advantages, the software does not address three-dimensional planar problems, cannot be used to analyze transient AC problems, and can only operate as a voltage-driven program under batch-run analysis mode [

his discussion about the use of the package as an educational tool, K.B. Baltzis [

The proposed CAD experiment focuses on the computation of the no-load magnetic loss in the laminated core of a single-phase shell-type power transformer using the non-linear finite element time-harmonic technique. This technique is an enhancement of modern electromagnetic field analysis programs that approximately models the effects of saturation and hysteresis on the fundamental component of the B-field, for different levels of the sinusoidal excitation represented by the H-field.

The main objectives of the experiment are:

· Inspection of computed values of the complex permeability, µ, at different locations of the ferromagnetic core;

· Determination of terminal quantities of the energizing winding, viz: terminal voltage, inductive impedance, active power and reactive power;

· Determination of the change in real power and reactive power that occurs in successive field problems with increasing magnetic loss.

The design of the test transformer including the relationship between the core’s geometry and rated power, choice of windings and assembly of the laminated core is based on the design rules presented in [

Micro and macro-eddy currents present in the magnetic core of an electric equipment produce ohmic loss and their own magnetic field—known as reaction field—that opposes the flux change originated by alternating excitations [

Design parameter | High-voltage winding | Low-voltage winding |
---|---|---|

N˚ of turns Apparent power (VA) RMS voltage (V) RMS current (A) Peak current (A) Type of conductor | 1100 58.0 220.0 0.26364 0.37284 28 AWG | 150 58.0 30.0 1.93333 2.73415 24 AWG |

In a magnetically linear time-harmonic solution, the permeability µ of the core is complex-valued and has the same value in the whole volume occupied by the laminated core. In non-linear problems, the value of the permeability varies from point to point in the region of the laminated core. The effect of increasing magnetic losses on terminal quantities like the active power P and reactive power Q can be readily inspected using the concept of reflected impedance of the energizing winding introduced in section 5.3.

In periodic field problems formulated with the H- and B-fields expressed in terms of phasor vector quantities, two distinct situations should be considered: (i) the effect of magnetic loss is negligible and, as illustrated in _{h} of the magnetic material. This effect is known as hysteresis-induced lag, and is illustrated in the diagram of

The vector fields H and B that produce the approximate elliptical hysteresis loop are calculated from two terminal quantities of the energizing winding, viz: the driving current I and the flux linkage λ, respectively. The value of the flux linkage is furnished by the finite-element solver together with the other winding’s circuit properties like total voltage drop and reflected impedance. Initially, let us consider the no-load operation of the test transformer. In the experiment, the 1100-turn high-voltage winding is kept open-circuited, and the 150-turn low-voltage winding is supplied by a 60 Hz sinusoidal current I. The magnitude of the no-load driving current is usually very small, and ranges from 1% to 5% of the winding’s rated current. Let H and B represent the phasor vectors computed from the driving current I and computed flux linkage λ with the exciting winding, respectively. According to Ampère’s circuital law, the magnetic field strength H is directly proportional to the driving current I and given by

H = ( N l m ) I , (1)

where N is the number of turns of the exciting winding and l_{m} is the mean value of the closed path that encircles the magnetomotive force NI. The magnetic induction B is directly proportional to the flux linkage λ and given by

B = 1 ( N A ) λ , (2)

where A is the cross-sectional area of the shell-type transformer’s central limb.

Important bulk magnetic properties of the material employed in the ferromagnetic core of machines and transformers can be obtained from the hysteresis loop. During a cyclic magnetization process, magnetic energy is converted into thermal energy and the amount of loss can be computed from the B-H plot or hysteresis loop. The area enclosed by the hysteresis loop gives the density of thermal energy, W_{h}, dissipated on completion of a magnetization cycle and is computed by

W h = ∮ H d B . (3)

If the operating frequency is f, the power loss, P_{h}, in units of watt is computed by

P h = ( ∮ H d B ) V f , (4)

where V is the volume of magnetic material. When post-processing the field solution in the FEMM CAD system, it is necessary to select the region occupied by the laminated core and choose the task “total losses” from the drop list of numeric integrations.

Electromagnetic field simulators contain a built-in library to facilitate the access to materials data files. The library for problems of magnetism contains several material property files of materials commonly used in the industry of electric equipment separated in file directories for conductors, permanent magnets and soft magnetic materials. The data file of a soft magnetic material may contain all necessary information to fully characterize the ferromagnetic core of a power transformer, i.e., 1^{st} quadrant B-H curve, loss curves, initial relative permeability, electric conductivity of laminations, lamination thickness and lamination fill factor. The physical properties and special attributes that characterize the M-36 steel laminated core of the test transformer are summarized in the block scheme of

To investigate the effect of magnetic loss on the various electromagnetic quantities of the test transformer, a parametric analysis has been carried out. The two parametric quantities are the electric conductivity σ of the laminations and the specified maximum hysteresis-induced angle f_{hmax}. Five different configurations of the laminated core have been considered, and the parameters that characterize each field problem are highlighted in the block scheme of _{hmax} and the electric conductivity σ of the steel laminations are specified as zero. The other four field problems concern a lossy laminated core and are distinguished by the value specified for the maximum hysteresis angle f_{hmax}. In these four problems, the electric conductivity of the steel laminations is specified as σ = 2 × 10^{6} S/m.

In the simulated work, the five time-harmonic problems are distinguished by different operating points (I, V) consisting of driving current I and estimated terminal voltage V. In a current-driven finite-element program, the simulation of no-load operating conditions in transformers and electric machines involves difficulties, and these are mainly related to the choice of the driving current. When the classical open-circuit test is carried out in the laboratory, the transformer operates with sinusoidal applied voltage but nonsinusoidal current. Whereas the winding rated voltage is well known as part of the transformer’s ratings, the magnetizing current must be measured or estimated, either analytically or numerically using, e.g., the method described in [

Five different sequences of problems driven by sinusoidal currents with increasing magnitudes have been used to identify the best fit match between the driving current—taken as the independent variable—and the rated voltage. In each sequence of problems, the magnitude of the sinusoidal driving current varies from 1.0% to 5.0% of the low-voltage winding’s rated current. Each sequence is related to one of the core configurations indicated in the block scheme of

In the graph of _{rated} = 30 V. Each characteristic representing the increase in terminal voltage is associated to a different parametric analysis. To

simulate a permeable but lossless magnetic core, both the maximum hysteresis angle f_{hmax} and the electric conductivity σ of the steel laminations are specified as zero. The dashed characteristic represents the results of the lossless core configuration. The four characteristics that represent a lossy laminated core are distinguished by the value specified to the maximum hysteresis angle: (i) for the characteristic marked with “*”, f_{hmax} = 0˚; (ii) for the characteristic marked with “·”, f_{hmax} = 10˚; (iii) for the characteristic marked with “▲”, f_{hmax} = 20˚; and (iv) for the characteristic marked with “¨”, f_{hmax} = 30˚.

Observation of the graph presented in _{rated} = 30 V, rms) differs for each sequence of problems. For the lossless core configuration (σ = 0; f_{hmax} = 0˚), the best fit match to the rated voltage occurs when the driving current represents 4.1% of the winding’s rated current. For each sequence of simulations, the value of the required driving current that yields a terminal voltage close to the rated value is indicated in the 4^{th} and 5^{th} columns of

Results of the parametric analysis are presented in the following. Magnitudes and phase angles of the computed flux linkages are indicated in the 6^{th} and 7^{th} columns of _{1} to λ_{5}, are basically distinguished by their lag phase angle with respect to the reference vector I. The

Type of analysis | Parametric quantities | Driving current | Flux linkage λ | |||
---|---|---|---|---|---|---|

Conductivity (S/m) | Max. hyst. angle (degree) | % of I_{rated} | Peak value (mA) | |λ|, (peak) (milliweber) | λ (angle) (degree) | |

Lossless | 0 | 0 | 4.10 | 112.10 | 112.92 | −8.3 × 10^{−7 } |

Hysteretic | 2 × 10^{6} | 0 | 4.20 | 114.83 | 112.29 | −11.7 |

Hysteretic | 2 × 10^{6} | 10 | 4.30 | 118.94 | 111.40 | −20.6 |

Hysteretic | 2 × 10^{6} | 20 | 4.50 | 123.04 | 113.37 | −29.3 |

Hysteretic | 2 × 10^{6} | 30 | 4.60 | 125.77 | 113.05 | −37.7 |

Type of analysis | Parametric quantities | Driving current | Estimate to rated voltage | ||||
---|---|---|---|---|---|---|---|

Conductivity (S/m) | Max. hyst. angle (degree) | % of I_{rated} | Peak value (mA) | |V|, (peak) (V) | Angle (degree) | Error in |V| (%) | |

Lossless | 0 | 0 | 4.10 | 112.10 | 42.57 | 89.9 | 0.33 |

Hysteretic | 2 × 10^{6} | 0 | 4.20 | 114.83 | 42.35 | 78.3 | 0.19 |

Hysteretic | 2 × 10^{6} | 10 | 4.30 | 118.94 | 42.02 | 69.4 | 0.97 |

Hysteretic | 2 × 10^{6} | 20 | 4.50 | 123.04 | 42.77 | 60.7 | 0.79 |

Hysteretic | 2 × 10^{6} | 30 | 4.60 | 125.77 | 42.65 | 52.3 | 0.51 |

radius of the arc that appears in the illustration—in the same length unit used to plot the vectors—is numerically equal to the mean value of the magnitude of flux linkages, ||λ_{i}||_{(mean) }= 112.61 milliweber. The diagram clearly shows that, higher magnetic losses – originated by increasing ohmic loss and reaction fields in the laminated core -, lead to an increase in the time phase difference between the driving current and resulting winding’s flux linkage.

Additional results of the parametric analysis are presented in _{rated}_{ }= 42.43 V, are presented in the 6^{th} and 7^{th} columns of ^{th} column shows that, in the five field solutions, the error in the magnitude of the estimates to the rated voltage is very small and below 1.0%. The graphical representation of driving currents and terminal voltages is shown in the vector diagram of _{1} to V_{5}, are basically distinguished by their leading phase angle with respect to the reference vector I. The radius of the arc that appears in the illustration—in the same length unit used to plot the vectors—is numerically equal to the mean value of the magnitude of voltage estimates, ||V_{i}||, mean = 42.47 V.

Eddy currents and hysteresis make the flux linkage λ lag the terminal current I,

and make the inductance of the energizing winding complex-valued [

L = λ / I , (5)

where λ is the flux linkage with the winding and I is the terminal current. The complex inductance L can be decomposed as

L = L r − j L i , (6)

where L_{r} and L_{i }are the real and imaginary components of L, respectively. The implied inductive impedance, Z, is given by

Z = j ω L , (7)

where ω is the angular frequency of excitation. Substitution of the right-hand side of (6) into (7) leads to

Z = ω L i + j ω L r . (8)

The electric circuit that represents the energizing winding of the test transformer at no-load operation is shown in _{S} denotes the specified sinusoidal driving current, and V_{S} denotes the computed terminal voltage. According to the terminology used by Stoll [

· the ωL_{i} term is the frequency-dependent circuit’s resistance. This circuit element is known as the winding’s “reflected resistance” because it represents the ohmic loss set up by eddy currents in the magnetic core and drawn from the energy supply;

· the ωL_{r} term is the frequency-dependent circuit’s inductive reactance. This circuit element is known as the winding’s “modified reactance” because its calculation considers the reduction in the flux-carrying capacity of the core caused by the reaction field created by the eddy currents present in the laminated core.

For the five different configurations of the laminated core, the values of the driving currents, I_{S}, are those indicated in the 5^{th} line of _{S}, are indicated in the 6^{th} and 7^{th} columns of _{i} and ωL_{r} are indicated in the 3^{rd} and 4^{th} columns of ^{th} and 7^{th} columns of

The illustration shown in the bar graph of _{i} and a reduction in the modified inductive reactance ωL_{r}. This change in reflected impedance implies: 1) an increase in the active power P drawn from the source supply; 2) and a reduction in the flow of reactive power Q in the energizing winding. Let I denote the “rms” value of the prescribed current I_{S} in the circuit of

P = ( ω L i ) I 2 (9)

and the reactive power, Q, is computed by

Q = ( ω L r ) I 2 . (10)

Let us now consider the changes in successive configurations, from the fourth to the fifth one, for example. In the fourth core configuration one has

P 4 = ( 169.93 ) ︸ ω L i , 4 ( 0.12304 / 2 ) 2 = 1.286 W (11)

and

Q 4 = ( 303.24 ) ︸ ω L r , 4 ( 0.12304 / 2 ) 2 = 2.295 var . (12)

In the fifth core configuration, one obtains P_{5} = 1.642 W and Q_{5} = 2.121 var. These figures show that, the alteration in core configuration yields an increase of 27.7% in real power P and a reduction of 7.58% in reactive power Q. According to data presented in the 8^{th} column of _{m4} = 3.04332 mJ and W_{m5} = 2.81253 mJ, respectively. As expected, both the magnetic stored energy W_{m} and the reactive power Q undergo the same percent reduction of 7.58%. The percent variations in real power, reactive power and magnetically stored energy that occur in the pairs of successive core configurations are summarized in

Type of analysis | Inductive impedance | Complex power | ||||
---|---|---|---|---|---|---|

ωL_{i} (W) | ωL_{r}_{ } (W) | |Z| (W) | Real, P (W) | Reactive, Q (var) | Stored energy (mJ) | |

Lossless | 0.47 | 379.76 | 379.76 | 2.93 × 10^{−3} | 2.38611 | 3.16403 |

Hysteretic, f_{hmax}=0˚ | 74.70 | 361.12 | 368.76 | 4.92 × 10^{−1} | 2.38084 | 3.15701 |

Hysteretic, f_{hmax}=10˚ | 125.63 | 334.58 | 357.39 | 8.68 × 10^{−1} | 2.31238 | 3.06618 |

Hysteretic, f_{hmax}=20˚ | 169.93 | 303.24 | 347.60 | 1.28621 | 2.29520 | 3.04332 |

Hysteretic, f_{hmax}=30˚ | 207.59 | 268.20 | 339.16 | 1.64185 | 2.12122 | 2.81253 |

Initial Configuration | Final Configuration | Increase in real power ΔP (%) | Reduction in reactive power ΔQ (%) | Reduction in stored energy ΔW_{m} (%) |
---|---|---|---|---|

Lossless | Hysteretic, f_{hmax} = 0˚ | 168.0 | 0.22 | 0.22 |

Hysteretic, f_{hmax} = 0˚ | Hysteretic, f_{hmax} = 10˚ | 76.0 | 2.88 | 2.88 |

Hysteretic, f_{hmax} = 10˚ | Hysteretic, f_{hmax} = 20˚ | 48.0 | 0.74 | 0.74 |

Hysteretic, f_{hmax} = 20˚ | Hysteretic, f_{hmax} = 30˚ | 27.7 | 7.58 | 7.58 |

Electromagnetic CAD systems can be effectively used as complimentary tools in the teaching of subjects around electromagnetics. Electrical machines is considered one of the most abstract and conceptually difficult subjects of the electrical engineering undergraduate curriculum, and this subject suffers from lack of student interest. The benefits of field simulations such as visualization of electromagnetic field distributions and enhanced capabilities to inspect phenomena of high complexity like magnetic saturation and distribution of eddy currents are determinant to overcome the difficulties inherent to the teaching of this subject.

The electrical machines virtual laboratory contains a series of experiments designed to teach complex electromagnetic phenomena. The experiments have been designed following recommendations of experts in engineering education and the authors’ previous experience. In the laboratory of simulations, it is given emphasis on demonstrating how to set-up field problems for solution, and how to examine the numeric results critically.

In the experiment described in the paper, the students get acquainted with the finite element time-harmonic technique. This technique is very compact, computationally cheap, and can be used to advantage whenever the excitations are sinusoidal, and only the steady-state solution is needed. At the pre-processing stage, students are faced with the difficulties inherent to the choice of driving currents in current-driven finite element programs. The most important post-processing task consists of building and analyzing the electric circuits that represent the sequence of field problems with increasing magnetic loss. Certainly, the experience gained in this laboratory meeting will be relevant to the student’s future work in industry or academia.

The authors acknowledge the financial support, in the form of scholarships, from the Brazilian Federal Agency for Postgraduate Studies (CAPES). The authors give thanks to David Meeker (dmeeker@ieee.org) for the use of the finite element CAD system. The authors also give thanks to Dr. Benedito Antonio Luciano (benedito@dee.ufcg.edu.br) for the useful comments on the manuscript.

The authors declare no conflicts of interest regarding the publication of this paper.

Nogueira, A.F.L., Weinert, R.L. and Maldonado, L.J.A.S. (2021) Finite Element CAD Experiments on the Effect of Magnetic Loss in Power Transformers with Laminated Cores. Journal of Electromagnetic Analysis and Applications, 13, 1-19. https://doi.org/10.4236/jemaa.2021.131001

In time-harmonic problems, the complex permeability can be defined in different ways. In the following, two different definitions are briefly discussed.

A1. The Complex Permeability Is Independent of the FrequencyLet H and B represent the phasor vectors computed from the prescribed driving current I and computed flux linkage λ with the exciting winding, respectively. In one of the definitions of complex permeability, it is assumed that magnetic hysteresis creates a “fixed” time phase difference, f_{h}, between vectors B and H that is independent of the frequency of excitation. This situation is illustrated in _{h}, known as the hysteresis angle of the magnetic material. The complex-valued magnetic permeability, µ_{h}, is then defined by

μ h = ‖ B ‖ ‖ H ‖ e − j ϕ h = μ r μ 0 e − j ϕ h . (A.1)

where μ_{r} is the magnetostatic relative permeability, and μ_{0} is the permeability of free space. In this approximate model, the hysteresis loop becomes an ellipse with its major axis making an angle of f_{h} radians with the H-axis [

In the electromagnetic field simulator FEMM, the time phase difference, f_{h}, is part of the definition of a linear time-harmonic problem. According to the illustration of FigureA1, it is firstly necessary to specify the magnitude of the relative magnetic permeability in the x- and y-directions. This is followed by the specification of the hysteresis-induced angles in the x- and y-directions. The default hysteresis lag angles are f_{hx} = f_{hy} = 0˚. In this example, the magnitude of the relative permeabilities for the M-36 steel is the same in both directions, i.e., µ_{x} = µ_{y} = 1616.

To make the matter concrete, let us consider the no-load operation of the test transformer with a 60 Hz driving current I_{s} = 0.123 mA supplying the low-voltage winding, and a set of numeric simulations distinguished by the user-defined hysteresis angle f_{h}. The driving current represents 4.5% of the winding’s rated

current and the transformer operates at rated terminal voltage. In this simplified analysis, it assumed that the hysteresis angles defined as part of the linear material properties are f_{hx} = f_{hy} = f_{h}. In the sequence of simulations, the hysteresis lag angle varies from 0˚ to 30˚, in steps of 5˚. In the problem definition, the special attributes of the M-36 steel laminations are: (i) electric conductivity of 2.0 × 10^{6} S/m; (ii) thickness of 0.635 mm; (iii) and fill factor of 0.98. In the analysis, the figure of prime interest is the computed complex permeability. Results of the simulations are presented in TableA1 together with FigureA2 and FigureA3.

Values of the computed complex permeabilities are presented in columns 4 and 5 of TableA1, and their graphical representation is shown in FigureA2. In the plot, the idealized, user-defined relative permeabilities μ h = 1616 e − j ϕ h are marked with an asterisk, whereas the perceived, computed relative permeabilities are marked with a square.

It is important to note that, for each linear time-harmonic problem, this approximate model produces a unique value for the computed complex permeability,

Specified hysteresis angle (degree) | User-defined permeability | Computed permeability | Error in |µ| (%) | ||
---|---|---|---|---|---|

Magnitude (p.u.) | Angle (degree) | Magnitude (p.u.) | Lag Angle (degree) | ||

0˚ | 1616 | 0˚ | 1580.76 | −2.9 | 2.2 |

−5˚ | 1616 | 5˚ | 1573.78 | −7.9 | 2.6 |

−10˚ | 1616 | 10˚ | 1566.95 | −12.9 | 3.0 |

−15˚ | 1616 | 15˚ | 1560.34 | −17.8 | 3.4 |

−20˚ | 1616 | 20˚ | 1553.98 | −22.7 | 3.8 |

−25˚ | 1616 | 25˚ | 1547.91 | −27.6 | 4.2 |

−30˚ | 1616 | 30˚ | 1542.17 | −32.5 | 4.6 |

µ_{h}, in the whole volume occupied by the laminated core. For the field problem where f_{hx} = f_{hy} = 20˚, the time-harmonic solution gives µ_{x} = µ_{y} = 1433.61-j599.68 p.u. at any point situated in the region of the laminated core. Observation of the results shows small discrepancies between computed and user-defined complex permeabilities. According to data presented in the 6^{th} column of TableA1, the percent error in computed magnitude varies from 2.2% to 4.6%. The averaged absolute error in computed phase angle is 2.7˚.

The approximate elliptical hysteresis loop for the field problem where f_{hx} = f_{hy} = 20˚ is exhibited in FigureA3. The vector fields H and B are calculated using (1) and (2). For a terminal driving current, I, given by

I = 123.037 ∠ 0 ∘ mA, (A.2)

and a computed flux linkage, λ, given by

λ = 30.220 ∠ − 22.7 ∘ mWb, (A.3)

the computed magnetic induction B is

B = 0.48 ∠ − 22.7 ∘ T, (A.4)

and magnetic field strength H is

H = 100 ∠ 0 ∘ A / m . (A.5)

A2. The Complex Permeability Is Frequency-DependentIn time-harmonic problems involving thin laminations, the effect of eddy currents and hysteresis can be encapsulated in the effective, frequency-dependent permeability defined by

μ eff = μ r e − j ϕ h 2 tanh [ − j ϕ h 2 j ω σ μ r μ 0 d 2 ] j ω σ μ r μ 0 d 2 , (A.6)

where f_{h} represents a constant phase lag between the H- and B-fields, σ is the electric conductivity of the laminations, d is the thickness of individual laminations, and ω is the angular frequency of excitation in rad/s. The term “ωσμ_{r}μ_{0}” present on the RHS of (A.6) reveals that, the calculation of the effective permeability considers the attenuation of the H-field below the surface of laminations due to the action of induced eddy currents.

The magnetization characteristic B = B(H) for the M-36 steel expressed as a 1^{st} quadrant curve is shown in FigureB1. The increments along the horizontal axis are 100.0 A/m between ticks, and the increments along the vertical axis are 0.2 T. Values of the magnetic induction B extend up to 1.50 T.