Mathematical Problems in Engineering
Volume 2010, Article ID 976930, 24 pages
http://dx.doi.org/10.1155/2010/976930
Research Article

## A Perturbative Method for Calculating the Impedance of Coils on Laminated Ferromagnetic Cores

Program of Postgraduates and Research in Electrical Engineering, Instituto Tecnológico de Morelia, Avenida Tecnológico 1500 Col. Lomas de Santiaguito, 58120 Morelia, MICH, Mexico

Received 24 February 2010; Revised 11 August 2010; Accepted 6 September 2010

Copyright © 2010 Serguei Maximov et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

A new set of formulas for calculating the self and mutual impedances of coils on straight and closed laminated ferromagnetic cores of circular cross-section has been derived. The obtained formulas generalize the well-known formulas for impedances of coils on homogeneous ferromagnetic cores, for the case of laminated cores, and improve the previously known formulas for laminated cores. The obtained formulas are fully consistent with Maxwell's equations and, therefore, offer an excellent accuracy. The perturbation theory and the average field technique are used to solve Maxwell's equations inside and outside the core. The solution inside the core can also be used in the analysis of thermal effects occurring inside the laminated core.

#### 1. Introduction

Analytic calculation of the self and mutual impedances for coils on ferromagnetic cores is not a new topic [17]. A rigorous analytic solution of Maxwell’s equations with the relevant boundary conditions is the natural basis for deriving accurate formulas for the coil impedances. In this approach, eddy currents induced into a laminated core and the frequency dependence of the impedance due to skin effect are properly represented. The mathematical complexity of the approach is related with the geometry and structure of the ferromagnetic core, so, rigorous analysis of a realistic case may become complicated. Nevertheless, the problem may be solved to any required degree of accuracy by using different approximate techniques and modern computational methods.

In the paper, an analytic study of electromagnetic field in a transformer with an infinitely long laminated ferromagnetic core of circular cross-section is presented. The formula for self and mutual impedance for filamentary turns is derived. Due to the laminated structure of the ferromagnetic core, some of the parameters involved in the Maxwell’s equations, such as magnetic permeability, permittivity, and electrical conductance, are discontinuous functions. As a result, the electromagnetic field becomes discontinuous, it loses the axial symmetry, and the use of numerous boundary conditions on the surface between each lamination and dielectric is required. The problem becomes too complicated.

In this paper, a new approach to the impedance calculation for transformers, which considers the laminated core as a whole, is presented. In order to simplify the problem, preserving a good accuracy, the use of the average electromagnetic field is proposed. In this approach, the problem becomes similar to the case of a uniform but anisotropic ferromagnetic core, where the magnetic permeability, permittivity, and electrical conductance are slightly different in the and directions. In turn, it means that there exists a small parameter which can be taken, for example, as the relation between the difference of the magnetic permeability in the and directions and the magnetic permeability of the uniform ferromagnetic core. Then, Maxwell’s equations can be solved by employing the perturbation theory [8, 9]. In this approach, every function involved in Maxwell’s equations is considered as a power series in the small parameter. The leading term of this series is the solution of the exactly solvable problem with the axial symmetry. The higher-order terms are the corrections that complement the leading term to the full solution of Maxwell’s equations. An approximate solution of the full problem is obtained by truncating the series, keeping only the first two terms of the power series. The higher orders become successively less important, and, therefore, they can be neglected.

#### 2. Equations for the Average Electromagnetic Field

In Figure 1, a ferromagnetic laminated core with filamentary turns is presented. The ferromagnetic core is taken to be infinitely long; each lamination is of permeability and conductivity . The laminations are separated by a dielectric of permittivity . The medium around the core is of permittivity and permeability . The core radius is . An energizing coil of filamentary turns of radius , placed at the coordinates , is considered. The ferromagnetic lamination width is , and the distance between two near laminations is .

Figure 1: Filamentary turns on an infinite ferromagnetic core.

The permeability, permittivity, and electrical conductance can be modeled by introducing the function presented in Figure 2, such that inside and outside the ferromagnetic lamination. Once the function is defined, the electrical conductance in the whole space can be represented in the following form: where is the Heaviside step function [10], such that inside and outside the core (the medium 1 is nonconducting).

Figure 2: Function .

The permittivity in terms of the function can be modeled as follows: and the permeability takes the form:

The mutual impedance between the energizing turn and a second filamentary turn of radius , as well as the self impedance of a nonfilamentary turn of small cross-section, follows directly from the solution for the electric field intensity. We assume that the electromagnetic field is quasistationary, and all capacitive effects may be taken into account separately.

As it has been mentioned before, the problem has not cylindrical symmetry. Moreover, due to unhomogeneities of the ferromagnetic core, the structure of Maxwell’s equations becomes complicated. In order to avoid such difficulty, we consider the electromagnetic field averaged with respect to the direction over . The procedure of transforming the exact electromagnetic fields to the average fields comprises having on hand, for example, a bump function [9], that is, compactly supported smooth (in the sense of ) function, such that(a) for ;(b) for and ;(c).

Let be one of the components of the electromagnetic field. Then, the averaged component is obtained as the following integral:

As a result, the averaged function has derivatives of all orders

The averaged Maxwell’s equations in the quasistationary approximation together with the continuity equation take the form where the averaged current density in the energizing coil is modeled by using the Dirac delta-function [10] as follows: where . The energizing filamentary turn carries a sinusoidal current .

The system of (2.6) is not closed due to the presence of the correlation functions , , and . In the medium 1, the parameters , , are constants. Then, these correlation functions take the following simple form: , , . In the medium 2, the correlation functions are to be estimated.

Let us approximate the bump function by the following discontinuous function:

Then, we obtain where is the average electric field intensity inside the ferromagnetic lamination, which is given by the equation

Similarly, the average electric field intensity in the dielectric medium between the laminations is

Applying the boundary condition for the average components and , we can approximately write

On the other hand, the full average electric field intensity has the following form:

Equations (2.12) and (2.13) form a system of equations which may be elementary solved with respect to the average . Substituting into (2.9), we finally obtain approximately where . The rest of the correlation functions can be calculated similarly by using the appropriate boundary conditions between a lamination and the dielectric. As a result, we obtain, that after averaging, the medium 2 becomes anisotropic, and the permeability, permittivity, and electrical conductance of the laminated ferromagnetic core acquire tensor properties, that is, , , and , where

Finally, for both mediums, the tensors , , and can be represented in the forms where (a)for medium 1, , , and ;(b)for medium 2,

is the identity matrix, and the tensor has the following form:

Then, after averaging, the system of Maxwell’s equations (2.6) takes the following form:

#### 3. The Zero-Order Perturbation Theory

The perturbation theory is an asymptotic method of solving differential equations [8, 9]. This method is commonly used in quantum mechanics [11, 12]. In order to develop the perturbation theory for the system of (2.19), the relation is taken as a small parameter. Then, the average electromagnetic field components should be expanded in a power series in as follows: where and are the leading terms of the series (3.1), and are the -order terms of the perturbation theory. In this paper, we were limited only to the first order of the perturbation theory.

Substituting the expansions (3.1) into the system of (2.19), we obtain for the zero-order perturbation theory in the frequency domain

Equations (3.2)–(3.4) are similar to the Maxwell’s equations for the case of the homogeneous ferromagnetic core with the core permeability , conductivity and permittivity . This problem carries axial symmetry, which makes the use of the cylindrical coordinates appropriate. Maxwell’s equations for coils on homogeneous ferromagnetic cores are solved in [1]. The complete solution of the system of (3.2)–(3.4) for the electric field intensity can be represented in the form where is the Fourier-transform of the electric field intensity with respect to the -coordinate

The electric field intensities outside and inside the core, that is, and , are given by the following equations: where , and , , and are the well-known Bessel functions [13]. The difference between the solutions (3.5), (3.7), and (3.8) and the results obtained in [1] is that (3.5), (3.7), and (3.8) contain the averaged core parameters and , which are smaller than the ferromagnetic parameters and .

The solutions (3.5), (3.7), and (3.8) are going to be used in the solution of equations for the first-order perturbation theory.

#### 4. The First-Order Perturbation

In the first-order perturbation theory, we obtain from the system of (2.19) in the frequency domain

Equations (4.1) are equivalent to Maxwell’s equations for a polarized medium with the electric current density , magnetization , and the polarization density , given by the following equations: where is the unit vector codirectional with the axis. These vector fields are calculated based on the solutions for the zero-order perturbation theory, and, therefore, in the first-order perturbation theory, they are considered as given functions. The contribution of the vectors and to the fields and can be taken into account in a future research.

From the definition of the coefficients and given above it follows, that in the case , these two constants are zero, that is, . As a result, the vector fields and become zeros, and therefore, the fields and are caused only by the magnetization . Further, only this case, in which the permittivities and are very close or equal, will be analyzed.

Substituting into (4.1), we obtain where is the vacuum permittivity. From (4.3) it follows that . Then, excluding the field from (4.3), we can obtain the equations for the electric field intensity for mediums 1 and 2, respectively,

Equation (4.4) is to be solved with the evident boundary condition at the core surface

The calculation of mutual impedance between the inducing coil and the filamentary turn, which is the aim of this paper, requires of the angular component of the intensity , integrated over the turn of the radius [1], that is, where is the mutual impedance in the zero-order approximation, which can be represented in the following form [1]: and is the first-order perturbation of the whole mutual impedance. Let us denote the component of the electric field intensity , averaged over all the directions , as , where

Then, the first-order perturbation term of the full mutual impedance can be calculated as follows:

Let us apply the averaging procedure (4.8) to (4.4). The advantage of the use of the averaging (4.8) is that the equations for acquire axial symmetry, that is, this averaged component loses the dependence on the angle . Then, after averaging procedure (4.8) and (4.4) can be written in the cylindrical coordinates as follows:for medium 1, for medium 2, where the right-hand side of (4.11) can be easily calculated as follows:

The magnetic field in the zero-order perturbation theory in virtue of (3.3) can be obtained from the solution for medium 2 (3.8) in the form

Then, (4.10) and (4.11) become where the Fourier transform of the field with respect to the -coordinate is taken. Equations (4.14) and (4.15) are to be complemented with the averaged boundary conditions at the core surface

The solution of (4.14), taking into account the limit requirement as , takes the following form:

Equation (4.15) can be solved by means of the method of variation of constants. The solution of (4.15), after applying the condition as , is the following: where

The integral in the right-hand side of (4.18) can be calculated giving as a result the Meijer -function [13]

The constants , are obtained from the boundary conditions (4.16). Substituting (4.17) and (4.18) into (4.16), we get and

Substituting constant into (4.17) and then, the inverse Fourier transform of the field into (4.9), we obtain as a result the correction for the mutual impedance

Finally, the whole mutual impedance of coils on laminated ferromagnetic core is obtained by substituting (4.7) and (4.23) into (4.6).

Note that both terms and take into account the laminated structure of the core.

Formulas (4.6), (4.7), and (4.23) generalize the results of [1] in the sense that the result for the mutual impedance obtained in [1] is a particular case of the formulas obtained in the present paper.

The self impedance of a coil of turns of small cross section can be obtained from (4.6), (4.7), and (4.23) by taking and successively to coincide with the point on the surface of each turn of the coil and then summing the results for each turn of the coil. The self impedance of a single turn is a particular case which can be obtained by substituting into the formula for self impedance of a coil.

#### 5. Special Cases

##### 5.1. Closed Cores

An approximate formula for the mutual impedance for turns on a closed toroidal core can be obtained by placing the finite-length core between two infinite plates of perfect magnetic materials, that is, infinite permeability and conductance [1], with respective boundary conditions. In order to provide a reasonable accuracy to this approximation, the core is to be sufficiently long (), and the turns should be taken close to the core surface. In practice, these assumptions are equivalent to the case of an infinitely long core with an infinite number of energizing turns distributed periodically along the core surface [1]. The formulas (4.6), (4.7), and (4.23) are still applicable in this case, provided that the number of turns becomes infinite, is the closed core length, and the Fourier transform is to be replaced by a Fourier series, due to the periodical structure of the problem. Formally, the transition from the Fourier transform to the Fourier series in (4.7) and (4.23) can be fulfilled by substituting [1]

As a result, the formula for the complete mutual impedance takes the form where is the inductance in air [1] and and and .

In Figures 3 and 4, contours of constant are presented for different frequencies and as a function of and . For illustrative purposes, the following parameters values have been taken:  S/m, , ,  m,  m,  m, μm, and μm. Figure 3 represents contours of constant at the frequency 1 MHz. Figure 4 illustrates contours of constant inductive part of the mutual impedance at  rad/s. The mutual inductance increases from the dark parts to the light regions. Figures 3 and 4, respectively, illustrate the effects of different ranges of flux penetration into the core. The behavior of the resistive part of the mutual impedance is somewhat similar, and therefore, the respective figures have been omitted.

Figure 3: Contours of constant at 1 MHz.
Figure 4: Contours of constant at 60 Hz.

In Figures 5, 6, 7, and 8, the behavior of the correction term at different frequencies is presented. At 60 Hz, the resistive part of the correction impedance brings a negative contribution to the complete mutual impedance near the energizing filamentary turn, whereas at 1 MHz it is positive.

Figure 5: Contours of constant at 1 MHz.
Figure 6: Contours of constant at 1 MHz.
Figure 7: Contours of constant at 60 Hz.
Figure 8: Contours of constant at 60 Hz.

From Figures 3 and 4 and [1], we can see that the behavior of the mutual impedance of turns on homogeneous and laminated ferromagnetic cores is somewhat similar. In order to compare the mutual impedance of coils on a laminated ferromagnetic core and on a homogeneous ferromagnetic core , and to estimate the maximum contribution of the laminated core to the mutual impedance, let us plot the relative increment of the mutual impedance and resistance. The core parameters are the same as those in Figures 38.

In Figure 9, the curve versus the frequency is shown at  m and  m. A maximum difference (up to 66%) under the conditions mentioned above occurs at the limit .

Figure 9: curve versus frequency.

In Figure 10, the curve at  m and  m is shown. The maximum contribution (17.5%) of the laminated structure of the core, for the given conditions, occurs at approximately 16.5 kHz.

Figure 10: curve versus frequency.
##### 5.2. Coils on an Infinite Core

The mutual and self impedance for coils on an infinite core can be obtained under the same assumptions that are given in [1]. We consider two coils of the radius and , widths and , thicknesses and , and the turn numbers and . The distance between two coils is . Then, the mutual impedance caused by the laminated ferromagnetic core is the following: where where is the modified Struve function [13]

In the case of the self impedance calculations of a coil on an infinite ferromagnetic laminated core we should put: , , and .

##### 5.3. Coils on a Closed Core

For the case of coils on a closed core, we can fulfill the same considerations that have been done in the previous section. Then, the formula for the mutual impedance for two coils on a ferromagnetic laminated toroidal core takes the form where

In Figures 11 and 12, the theoretical curves and experimental points for the resistive (Figure 11) and inductive (Figure 12) components of self impedance in the frequency range 102 Hz to 1 MHz for test coil on an outer limb of the 25 kVA core are presented. The experimental data and the core and coil parameters were taken from [3]:  Ω−1, μm, μm, , ,  m,  m, and the test core has 12 turns of heavy-gauge wire, and so forth. The inductance in air is calculated following the recommendations given in [3]. Both theoretical curves show a good agreement between measurements and calculations, which confirms once again the idea expressed in [3] about the validity of treating the laminated core as solid.

Figure 11: Resistive component of self impedance (test coil on an outer limb of the 25 kVA core). Dashed line corresponds to the solid core and the continuous line represents the laminated core.
Figure 12: Inductive component of self impedance (test coil on an outer limb of the 25 kVA core). Dashed line corresponds to the solid core and the continuous line represents the laminated core.

Also, in [3] it has been expressed that in spite of a good accuracy of obtained formulas, the prediction of transient behavior of transformer under general conditions requires some more accurate considerations. In particular, the consideration of the differential mutual impedance is required. Any substantial errors obtained in the calculated quantities and for the solid core (dotted lines in Figures 13 and 14) may result in significant errors in transient oscillations. Therefore, it is important to keep these errors small.

Figure 13: Resistive component of differential impedance between two test coils (coils on outer limb of the 25 kVA core at separation of 0.2 m). Dotted line corresponds to the solid core, the continuous line represents the laminated core with one parameter calculated and the dashed line is for two calculated parameters and .
Figure 14: Inductive component of differential impedance between two test coils (coils on outer limb of the 25 kVA core at separation of 0.2 m). Dotted line corresponds to the solid core, the continuous line represents the laminated core with one parameter calculated and the dashed line is for two calculated parameters and .

To improve these results, it was proposed in [3] to enhance the core resistivity. This proposal seems reasonable because this is that really occurs inside the core. Indeed, the averaged conductivity, that is, presented in (5.9), is smaller than the conductivity of lamination, that is, , and therefore, the core resistivity increases. Nevertheless, on the contrary, a simple numerical analysis shows that the reduction of the averaged conductivity results in a larger error between the theoretical and experimental results. To corroborate this, it is enough to substitute in (5.6) and plot the differential resistivity and inductance curves. Therefore, a simple change of the core resistivity cannot improve the agreement between the theoretical and experimental results.

On the other hand, taking into account the first-order perturbative term, that is, assuming , improves the theoretical results only a little. A reasonable explanation of this fact can lie in how the correlations , , and have been calculated. The correct calculation of these correlations implies having on hand the solutions for and before solving Maxwell’s equations, which is impossible. Then, there are two possible ways to deal with these terms. Both methods consider the parameters , , and as unknown constants. In the first one, these parameters should be recalculated after solving the exact Maxwell’s equations. The second way is to calculate these constants by comparing the theoretical and the experimental results. However, the first way is not feasible in this particular case because the technique proposed in this paper has the purpose of avoiding the necessity of calculating the exact fields and . Then, the second program should be carried out.

In Figures 13 and 14, the theoretical and experimental behavior of the differential resistivity and inductance, respectively, is presented. The dotted line shows the theoretical behavior of the differential impedance in the case of the solid core [3]. The continuous line represents the theoretical differential impedance for the laminated core with the parameter calculated to obtain a good agreement with the experimental differential inductance. In this case, to obtain the parameter , one experimental point as a reference is enough. However, to obtain the best result, the least squares method can be used. The obtained result is  H, and the theoretical result obtained with this parameter value shows an excellent agreement with the experimental data in the case of the differential inductance, and a better approximation for the differential resistance.

In turn, in this method, the rest of parameters, that is, and cannot be taken as zero, and therefore, the impedance terms corresponding to the vectors and , cannot be neglected. Figures 13 and 14 show the differential resistivity and inductance calculated (the dashed line), taking into account the terms and in the first-order perturbation theory (the corresponding contribution of these terms to the impedance is calculated similarly as the contribution corresponding to the term ). The parameters  H,  F/m, and Ω−1 are calculated to obtain a better correlation with the experimental data by using the least squares method. We can see that, involving the vectors and into the differential impedance calculation, we obtain a better approximation of the theoretical differential resistivity to the measured data, especially in the high frequency range.

In [3], it was proposed to adjust the theoretical results with measurements by introducing an additional empiric term in the formula for the differential impedance. On the contrary, the formulas obtained in this paper are theoretical and completely consistent with Maxwell’s equations. Therefore, it seems more logical and consistent to compare our results with the experimental data and the respective theoretical results for the solid core [3], which has been carried out in this section.

#### 6. Conclusions

New formulas for calculating the self and mutual impedances of coils on laminated ferromagnetic cores of circular cross-section have been derived. These formulas generalize the previously established formulas for the homogeneous ferromagnetic cores [1], for the case of the laminated cores. In the limit and , the obtained formulas are transformed into the well-known formulas for the self and mutual impedances for coils on homogeneous ferromagnetic cores published in [1].

The obtained result for the frequency-dependent self and mutual impedances for coils on laminated core is compared with the experimental and theoretical results for coils on a practical transformer core, published in [3]. The presented numerical results show a very good correspondence of the new formulas to the measured data. These new formulas do not require introducing the additional empiric terms in order to obtain a better correspondence to the experimental points; all the new terms, involved in the impedance formula, are fully consistent with Maxwell’s equations. Moreover, the accuracy, reached in the new formulas, is not the limit; it can be improved by a more thin analysis of the correlation terms , , and . Thus, all the stated above can be considered as a method for further theoretical research.

As the proposed method implies the strict solving of Maxwell’s equations, it can be used not only for impedance calculation, but also for the electromagnetic fields analysis inside the core. For example, substituting (4.22) into (4.18), we can obtain the averaged electric field intensity inside the core, which, in turn, can be used in the analysis of the laminated core heating. In turn, it calls for numerical evaluation of different special functions such as the Meijer -function [13]. However, with available computers and established numerical techniques, this is no longer a problem [1]. The programming of the Meijer -function is not a great obstacle; this function can be calculated with any required degree of accuracy, which is not particularly time consuming.

Further theoretical work could be done to establish the exact formulas for the toroidal cores and cores of different geometries. Also, it is of interest the theoretical analysis of the thermal effects occurring inside the laminated core during the transformer operation. The proposed methods, taken as a principle for solving this problem, make possible such type of analysis.

#### Acknowledgments

This paper was supported by the Dirección General de Educación Superior Tecnológica and the Consejo Nacional de Ciencia y Tecnología-Sistema Nacional de Investigadores of México.

#### References

1. D. J. Wilcox, M. Conlon, and W. G. Hurley, “Calculation of self and mutual impedances for coils on ferromagnetic cores,” IEE Proceedings A, vol. 135, no. 7, pp. 470–476, 1988.
2. J. Avila-Rosales and F. L. Alvarado, “Nonlinear frequency dependent transformer model for electromagnetic transient studies in power systems,” IEEE Transactions on Power Apparatus and Systems, vol. 101, no. 11, pp. 4281–4288, 1982.
3. D. J. Wilcox, W. G. Hurley, and M. Conlon, “Calculation of self and mutual impedances between sections of transformer windings,” IEE Proceedings C, vol. 136, no. 5, pp. 308–314, 1989.
4. H. Akçay and D. G. Ece, “Modeling of hysteresis and power losses in transformer laminations,” IEEE Transactions on Power Delivery, vol. 18, no. 2, pp. 487–492, 2003.
5. S. Maximov and V. H. Coria, Métodos analíticos en el cálculo de impedancias de transformadores con diferentes núcleos, M.S. thesis, Programa de Graduados en Investigación en Ingeniería Eléctrica, México, Mexico, 2008.
6. J. H. Harlow, Electric Power Transformer Engineering, IEEE, CRC Press LLC, Boca Raton, Fla, USA, 2004.
7. L. L. Grigsby, Electric Power Engineering Handbook, IEEE, CRC Press LLC, Boca Raton, Fla, USA, 1998.
8. M. Reed and B. Saĭmon, Methods of Modern Mathematical Physics, Academic Press, San Diego, Calif, USA, 1978.
9. J. Mathews and R. L. Walker, Mathematical Methods of Physics, Addison-Wesley, Redwood City, Calif, USA, 2nd edition, 1969.
10. V. S. Vladimirov, Generalized Functions in Mathematical Physics, Mir, Moscow, Russia, 1979.
11. L. I. Schiff, Quantum Mechanics, McGraw-Hill, New York, NY, USA, 2nd edition, 1965.
12. L. D. Landau, M. D. Lifshitz, and L. P. Pitaevskii, Electrodynamics of Continuos Media, U.S.S.R. Academy of Sciences, Moscow, Russia, 2nd edition, 1984.
13. L. C. Andrews, Special Functions for Engineers and Applied Mathematicians, Macmillan, New York, NY, USA, 1985.