The analytical solution of 2D electromagnetic wave equation for eddy currents in the cylindrical solid rotor structures
Lale T. Ergene^{1} , Yasemin D. Ertuğrul^{2}
^{1, 2}Istanbul Technical University, Istanbul, Turkey
^{1}Corresponding author
Journal of Vibroengineering, Vol. 19, Issue 5, 2017, p. 39113925.
https://doi.org/10.21595/jve.2017.18257
Received 16 February 2017; received in revised form 11 April 2017; accepted 30 April 2017; published 15 August 2017
JVE Conferences
The paper presents the closedform solution of two dimensional (2D) electromagnetic wave equation for eddy current problems in cylindrical structures. The magnetic field calculation is a complex issue for electrical machines especially. The paper provides an analytical solution for solid cylindrical structures with linear magnetization. This approach is extended to the electrical machines with the solid rotors. Multilayered cylindrical geometry is used for the solution. Two dimensional parabolic partial differential equation is solved for each layer. According to the analytical solution, eddy current and eddy current losses are calculated in the solid region and compared to the results of the finite element model.
Keywords: electromagnetic wave equation, diffusion equation, parabolic partial differential equations, eddy current, closed form solution, solid rotor, boundary conditions.
1. Introduction
The magnetic field calculation is a very important issue for the electrical machine designers. In the literature, there are different approaches to obtain the magnetic field by using analytical methods, equivalent circuit model, reluctance mesh method or numerical methods. Solving the governing partial differential equations for different boundary conditions analytically is not very easy problem. Because of that the numerical methods such as Finite Element Analysis (FEA) is preferred over the analytical methods for the electrical machine most of time [1]. However, closed form solution of magnetic field provides a fast and approximated solution of the field analysis for the beginning of the design. According to magnetic field calculation with analytical method, sizing equations of motors can be refined easily.
Magnetic field analyses are achieved for different type of motors with different methodologies. The induction machines are analyzed with the analytical methods in the literature. For example, Poynting’s theorem is used for the calculating output power and joule losses [2]. In another study, reluctance network method is used for to determine eddy current in the stator core [3]. Equivalent circuit method is used with new approach in canned solid rotor induction machine [4]. The eddy current calculation is realized in solid conducting region using transverse alternating magnetic field [5]. In some of applications, finite element method and analytical method are combined to decrease the analytical method complexity [6]. The motor with permanent magnet rotor has analytical solution using closed form solution of Maxwell equations [7, 8]. In this study, there is no current source and stator is solid conductor with eddy current. Double Fourier series are applied on the motor with slotless permanent magnet rotor with current layer, which correspond to time and space harmonics [9].
In this paper, the induction motor is modelled as three concentric cylindrical layers representing the different regions showing different material properties. The governing two dimensional partial differential equations are solved for the source, conducting, and nonconducting regions including air gap. The MATLAB^{®} is used for the analytical calculations. Accuracy of the method is verified with the finite element simulations. This method is applied in rotating cylinder excited by a stationary coil. Then the importance of this method is emphasized for different motor applications in future works [10].
2. Electromagnetic vibration
The magnetic forces between stator and rotor because of the attraction and the repulsion cause the magnetic vibration. These magnetic forces can be in low and high frequency ranges. The low frequency range includes slip frequency and once, twice and higher multiples of the fundamental frequency. The high frequency range includes frequencies related to the product of the number of the rotor slots and the rotational speed of the rotor. These frequencies may be up to kHz for small machines while they may be as low as a few hundred Hz for large machines with a high number of poles. The magnetic vibrations are dominated by the flux flowing between air gap producing rotating force waves. Other combinations of stator and rotor space harmonics and saturation harmonics can also lead to important magnetic vibrations.
The magnetic stored energy in the air gap can be obtained from the total flux density, which is the sum of the rotor and stator flux densities. The total stored energy in the air gap changes while rotor rotates, and a torque occurs. The total flux density in the air gap with respect to time and space is [11]:
${\mathbf{B}}_{rotor}\left(\theta ,t\right)$and ${\mathbf{B}}_{stator}\left(\theta ,t\right)$ are the form of:
The source of electromagnetic vibration can be traced to the harmonics of the magnetic flux density distribution. Moreover, the flux density analysis can also give information about the torque profile. In a symmetrical stator winding, the MMF harmonic components are phasebelt harmonics with order of ($2q\pm 1$), stator slot harmonics with order of (${N}_{s}\pm p$), and rotor slot harmonics with order of (${N}_{r}\pm p$). In the case, the stator doesn’t have the slots so stator harmonic components will not be considered. The MMF distribution in the air gap can be given as Eq. (3) [11]:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+B\mathrm{c}\mathrm{o}\mathrm{s}\left(\left(2q1\right)p\theta +\omega t\right)+C\mathrm{c}\mathrm{o}\mathrm{s}\left(\left(2q+1\right)p\theta \omega t\right)+D\mathrm{c}\mathrm{o}\mathrm{s}\left(\left({N}_{s}p\right)\theta +\omega t\right)+\mathrm{E}\mathrm{c}\mathrm{o}\mathrm{s}\left(\left({N}_{s}+p\right)\theta \omega t\right)+F\mathrm{c}\mathrm{o}\mathrm{s}\left(\left({N}_{r}p\right)\theta \left[\left(1s\right){N}_{r}p\right]\frac{\omega}{p}t\right)+G\mathrm{c}\mathrm{o}\mathrm{s}\left(\left({N}_{r}+p\right)\theta \left[\left(1s\right){N}_{r}+p\right]\frac{\omega}{p}t\right)+\dots .$
Magnetic flux can cross the interface between air and infinitely permeable material only in a direction perpendicular to the surface of the material. A magnetic force is occurred in the direction of flux. The magnitude of the force is given by in units of pressure or force per unit area:
where $B$ is the magnetic flux density in Wb/m^{2} and ${\mu}_{0}$ is the permeability of free space (4$\pi $×10^{7} H/m). The analysis of the $B$ field is very critical to understand the harmonic components of electromagnetic forces since the electromagnetic forces are often found to be the major noise and vibration source in many cases. Because of this reason, the analytical approach is focused on the determining the magnetic vector potentials and flux density for the solid rotor structures. Eddy currents and resulting extra losses because of the solid structure will be affecting the system efficiency. The analytical calculation of eddy current losses is also given in the paper.
3. Physical description
The crosssectional model of the solid rotor induction motor is shown in Fig. 1. The rotor region is an eddy current region. Moment is created by the rotor magnetic field which arises from eddy current in the rotor. The region lying between ${r}_{2}<r<{r}_{3}$ in Fig. 1 is the stationary region called the stator (Region 3). The region lying between ${r}_{1}<r<{r}_{2}$ is a bounded free space called air gap between stator and rotor (Region 2). The eddy current region is defined where $0<r<{r}_{1}$ (Region 1) and called ‘rotor’. Beyond the stator where ${r}_{3}<r<\infty $ is the free region (Region 4). Stator magnetic field passes through rotor over the air gap. The current density is given between air gap and stator yoke (current sheet).
Fig. 1. Two dimensional crosssectional of a solid rotor motor model
Calculations are made in cylindrical coordinate system. The following is a summary of some of main limitations:
• The problem is considered as two dimensional and linear.
• End effects are neglected.
• The permeability (${\mu}_{r1}$, ${\mu}_{r3}$) and conductivity ($\sigma $) of materials are time invariant and independent of the field.
• Saturation and hysteresis effects are neglected.
• The stator is laminated and induced eddy current are neglected.
• The stator current source is represented by a current sheet on the inner bore of the stator.
• The rotor is solid and the permeability and conductivity are isotropic and single valued.
• The field solution is obtained in cylindrical polar coordinates.
• The vector potential has only a ‘$z$’ directed component equal to ${A}_{z}$.
• The flux density has two components ${B}_{r}$ and ${B}_{\theta}$.
4. Closed form solution
4.1. Formulation
Derivation of the two dimensional timeharmonic eddy current is supplied using diffusion equations in the cylindrical polar coordinates. Maxwell equations are written for every region of the solid rotor motor:
Ampere’ Law is given in Eq. (5) and Eq. (6) where $\mathbf{H}$ is the magnetic field intensity, $\mathbf{B}$ is the magnetic flux density, $\mathbf{A}$ is the magnetic vector potential, ${\mathbf{J}}_{e}$ and ${\mathbf{J}}_{s}$ are the eddy current density and source current density, and ${\mu}_{0}$ and ${\mu}_{r}$ represents permeability of space and relative permeability of the material. According to Faraday’s Law, relationship of the electric field and the magnetic flux density are given in Eq. (7). The relationship of the electric field and the magnetic flux density are written in Eq. (7) and Eq. (8) by using $\mathbf{B}=\nabla \times \mathbf{A}$ and $\mathbf{J}=\sigma \mathbf{E}$ equations, where $\mathbf{E}$ is the electric field:
Using vector identity rule, curl of the magnetic field becomes in Eq. (9):
The magnetic vector potential is assumed to be $\mathbf{A}=\left\mathbf{A}\right{e}^{j\omega t}$. The defined problem is in twodimensional and magnetic flux density has two components which are ${\mathbf{B}}_{r}$ and ${\mathbf{B}}_{\theta}$ in polar coordinates and magnetic vector potential has only $z$ component according to relationship with $\mathbf{B}$ and $\mathbf{A}$. The magnetic vector potential is written using trigonometric equivalent in Eq. (10) and then obtained derivative using trigonometric identities in Eq. (11):
In polar coordinates, $\mathbf{A}$ has only $z$ – component and curl of magnetic field in $z$ direction is represented in Eq. (12). The curl of magnetic vector potential is equal to the difference of the current densities which are eddy current and source current density. In polar coordinates:
The motor regions have different magnetic and electrical properties which are given Table 1. These properties of materials require solving partial differential equations in Eq. (13).
Table 1. Properties of motor regions
Regions

Properties

Region 1 (rotor)

${J}_{s}=$ 0, ${\mu}_{r1}$, $\sigma $

Region 2 (air gap)

${\mu}_{0}$

Region 3 (stator)

${J}_{s}=$ 0,${\mu}_{r3}$, $\sigma =$ 0

Region 4 (free space)

${\mu}_{0}$

Region 1 is the eddy current region where ${\mathbf{J}}_{s}=$ 0, characterized by two dimensional diffusion equations in polar coordinates. It is represented in Eq. (14):
The diffusion equation which is parabolic partial differential equation solved with separation of variables method. Magnetic vector potential can be separate two different variables and then diffusion equations is refreshed which is given Eq. (15) and Eq. (16):
In closed form, the equation can be written as Eq. (17):
Eq. (17) can be simplified with using coefficients which is ${n}^{2}$ and ${k}^{2}$ to transform into a second order ordinary differential equation given in Eq. (18):
Eq. (18) can be expressed as a parametric Bessel equation which is written in general form in Eq. (19) [12]. ‘$j{k}^{2}$’ term is simplified as ${\beta}^{2}$ to simulate parametric modified Bessel equation:
The parametric modified Bessel equation solution has two components which are given in Eq. (20):
${I}_{n}$ and ${K}_{n}$ are the modified Bessel’s functions and $n$ is the order of equation and represented as $ber$ and $bei$ functions:
First part of the Eq. (15) is arranged in Eq. (23) with modified Bessel’s functions:
Using the separation of variables method, second part of the separated variable is given in Eq. (24):
This is ordinary differential equation and solved as Eq. (25):
The general solution of Eq. (25) is given as summation of cosine and sinus functions in Eq. (26) which is arranged and simplified:
The solution of magnetic vector potential is related with $r$ and$\theta $ by use of the separation of variables method shown in Eq. (27):
The magnetic vector potential equation is given in Eq. (27). This equation is represented by trigonometric equations comes from modified Bessel functions.
4.2. Calculation for regions
The closed form of the diffusion equation is solved with the separation of variables method. The found magnetic vector potential equation can be used for every regions and boundaries. Eq. (27) can be applied every region of the model given in Fig. 1. According to regions magnetic and electrical properties, coefficients of the magnetic vector potential are simplified.
Region 1: ($0<r<{r}_{1}$) In rotor region, the field at $r=$ 0 is finite and has only cosine components, we have ${D}_{n}=$ 0 and ${F}_{n}=$ 0. The rotor has electrical conductivity and it is called the eddy current region. The magnetic vector potential of rotor regions is given in Eq. (28):
Region 2: (${r}_{1}<r<{r}_{2}$) The air gap region is between ${r}_{1}$ and ${r}_{2}$ and magnetic vector potential is in Eq. (29):
Region 3: (${r}_{2}<r<{r}_{3}$) The stator region has no electrical conductivity and it is eddy current free region and represented magnetic vector potential in Eq. (30):
Region 4: (${r}_{3}<r<\infty $) The last air region begins from the stator outer diameter to infinite. The magnetic vector potential is in below as in Eq. (31):
Normal component of the magnetic flux density (${B}_{r}$), and tangential components of the magnetic field (${H}_{\theta}$) are equal to each other at boundaries. Therefore, the boundary conditions:
The magnetic flux densities and the field intensities for each region are obtained from magnetic vector potential equations and given in Eqs. (3744).
Region 4 equations:
Region 3 equations:
Region 2 equations:
Region 1 equations:
The boundary conditions are considered to find coefficients of the magnetic flux densities. The normal component of magnetic flux densities (${B}_{r}$), and the tangential components of magnetic field intensities (${H}_{\theta}$) are equal to each other at boundaries related to each region.
At the boundary of $r={r}_{3}$, the boundary conditions of ${B}_{r3}={B}_{r4}$, ${H}_{\theta 3}={H}_{\theta 4}$ should be satisfied. The related equations are given in Eq. (45) to find out ${D}_{3n}$ and ${C}_{3n}$ coefficients in terms of ${D}_{4n}$. The revised equations are given in Eqs. (4647):
$\sum _{1}^{n}\frac{n}{{\mu}_{0}{\mu}_{r3}}\left({{C}_{3n}{r}_{3}^{n1}D}_{3n}{r}_{3}^{n1}\right)\mathrm{cos}\left(n\theta \right)=\sum _{1}^{n}\left(\frac{n}{{\mu}_{0}}{D}_{4n}{r}_{3}^{n1}\right)\mathrm{cos}\left(n\theta \right),$
${D}_{3n}={D}_{4n}\left(\frac{1+{\mu}_{r3}}{2}\right),{C}_{3n}={D}_{4n}\left(\frac{1{\mu}_{r3}}{2}\right){r}_{3}^{2n},$
At the boundary of $r={r}_{1}$, the boundary conditions of ${B}_{r1}={B}_{r2}$, ${H}_{\theta 1}={H}_{\theta 2}$ should be satisfied. The related equations are given in Eq. (48) to find out ${D}_{2n}$ and ${C}_{2n}$ coefficients in terms of ${C}_{1n}$. The revised equations are given in Eqs. (5152):
$\sum _{1}^{n}\frac{nk{C}_{1n}}{{\mu}_{0}{\mu}_{r1}}\left[Be{r}_{n}^{\text{'}}\left(k{r}_{1}\right)+jBe{i}_{n}^{\text{'}}\left(k{r}_{1}\right)\right]\mathrm{cos}\left(n\theta \right)$
$=\sum _{1}^{n}\frac{n}{{\mu}_{0}}\left({{C}_{2n}{r}_{1}^{n1}D}_{2n}{r}_{1}^{n1}\right)\mathrm{cos}\left(n\theta \right),$
${C}_{2n}=\frac{{C}_{1n}}{2}\left(\frac{P}{{r}_{1}}+\frac{kQ}{{\mu}_{r1}}\right){r}_{1}^{n+1},{D}_{2n}=\frac{{C}_{1n}}{2}\left(\frac{P}{{r}_{1}}\frac{kQ}{{\mu}_{r1}}\right){r}_{1}^{n+1},$
where:
At the boundary of $r={r}_{2}$, the boundary conditions of ${B}_{r2}={B}_{r3}$, ${H}_{\theta 3}{H}_{\theta 2}={J}_{s}$ should be satisfied. The related equations are given in Eq. (54) to find out ${D}_{4n}$ and ${C}_{1n}$ coefficients. $E$ and $F$ expressions given in Eq. (53) are used to simplify the equations:
With the boundary conditions, all of the coefficients are obtained which are involved in the equations of magnetic vector potential, magnetic field intensity and also magnetic flux density. Coefficients are given in Eqs. (5960):
In these equations, some reductions carried out to arrange equations which are ‘$U$’ and ‘$V$’ which are given as follows:
Since all the integration constants have been determined, the respective vector potentials, the radial components of the fluxdensities and the tangential components of the magnetic fields in the different regions are obtained as follows.
Region 1:
Region 2:
Region 3:
Region 4:
The eddy current density in eddy current region is given in the Eq. (75), by using the Eq. (63):
After deriving the current density equation, the power loss in eddy current region can be obtained by using Eq. (76) in polar coordinates:
Substituting $Je$ from Eq. (75) and its conjugate in Eq. (76) for power loss, one obtains after some algebraic manipulation:
5. Finite element model
The finite element model of the solid rotor structure is constructed to verify and compare the analytical model. The magnetic material used in both rotor (${\mu}_{r1}$) and the stator (${\mu}_{r3}$)_{}is given in Fig. 2. The flux density through the air gap is calculated by the finite element results of solid rotor induction motor. Calculations are obtained for different operating conditions of the motor. The rotor has the conductivity also. The motor has 4 kW 380 V 1459 min^{1}, 50 Hz power ratings.
Fig. 2. Magnetization curve of the stator and rotor material
The current density of the solid rotor motor is created by three phase stator currents. The current source of the solid rotor motor is given as the first harmonic of the signal in the solid rotor machine:
The line current density injected in the FEM model to make the model similar to the analytical. The injected line current in the motor and the motor mesh are shown in Fig. 3. The rotor conductivity σ is chosen as 3×10^{6} S/m in both analytical and numerical calculations. ${\mu}_{r1}$ and_{}${\mu}_{r3}$ are defined in the linear region as 750. The air region permeability is ${\mu}_{0}=$ 4$\pi $10^{7} H/m. The length of the rotor is 145 mm. The diameters are as follows: ${r}_{1}=$ 58.25 mm, ${r}_{2}=$ 58.5 mm, ${r}_{3}=$ 85 mm. The frequency is 50 Hz. The output torques at the rated value for the both analytical and numerical models are coherent and the relative error is less than 5 %.
The phase currents for one cycle from the FEM are given in Fig. 4.
Fig. 3. The quarter model of the line current density and its mesh
a)
b)
Fig. 4. Stator phase currents of solid rotor induction motor.
6. Results
Slip speed of an induction motor is defined as the difference between the synchronous speed and the actual rotor speed. The solid rotor induction machine is analyzed at steady state for different slip values. The equiflux lines are given in Fig. 5. The skin effect and phase delay because of the eddy currents can be observed with increasing slip value.
Fig. 5. Equiflux lines for different slip values
The flux density variation in the air gap for different slip values are given in Fig. 6. The four pole distribution is obvious for each slip value.
Fig. 6. Flux density in the air gap for different slip values
In the analytical model of the solid rotor machine, first harmonic of the air gap flux density is calculated. Fig. 7 shows normal component of the flux density and its first harmonic. Flux density and its harmonics are calculated by the finite element method. The finite element result of the solid rotor induction motor is taken as the reference for analytical calculations. Air gap flux densities are compared with each other in Fig. 8. Calculated first harmonic is compared with the analytical result of the air gap flux density. The results are coherent and the relative error is less than 10 %.
Fig. 7. Air gap flux density and its first harmonic with FEM
Fig. 8. Comparison of the analytical and numerical results of the flux density around the air gap
The finite element and analytical results of eddy current density and eddy current loss for the solid rotor are given in Table 2. They are coherent.
Table 2. Comparison of the eddy current density and the power loss in the rotor at the rated operation condition (${s}_{N}=$ 2.73 %)
Method

${J}_{e}$ (pu)

$Pe$ (pu)

FEM model

1

1

Analytical model

1.08

1.12

7. Conclusions
Solid rotor cylindrical structure is analyzed using closed form solution of magnetic field in the case study of a solid rotor induction motor. The analysis is based on the solution of the modified Bessel functions of the parabolic wave equation for an eddy current medium. The obtained analytical expressions have been compared to the finite elements yielding almost the same results for the cylindrical same geometry. Closedform solution of magnetic fields in induction motors provides a fast and easy method of field analysis, which can help in sizing studies in design and evaluation of performance and in the determination of equivalent circuit parameters. Although a rigorous representation of the geometry is not within the scope of closedform studies, a quick and easy solution for the eddy current region is obtained. In the paper, a method of evaluating magnetic fields in a cylindrical solid structures are presented and formulae for eddy currents and power losses are derived.
References
 Ahfock T., Hewitt A. Modeling of the solid rotor induction motor. 15th Australasian Universities Power Engineering Conference (AUPEC 2005), Hobart, Australia, 2005. [Search CrossRef]
 Markoviç M., Perriard Y. An analytical solution for the torque and power of a solidrotor induction motor. Proceedings of IEEE International Electric Machines and Drives Conferences (IEMDC 2011), 2011, p. 10531057. [Search CrossRef]
 Mirzayee M., Mehrjerdi H., Tsurkerman I. Analysis of a high speed solid rotor induction motor using coupled analytical method and reluctance networks. IEEE/ACES International Conference on Wireless Communications and Applied Computational Electromagnetics, Honolulu, Hawaii, 2005. [Publisher]
 Ergene L. T., Salon S. J. Determining the equivalent circuit parameters of canned solid rotor induction motors. IEEE Transactions on Magnetics, Vol. 41, Issue 7, 2005, p. 22812286. [Search CrossRef]
 Perry M., Jones T. Eddy current induction in a solid conducting cylinder with a transverse magnetic field. IEEE Transactions on Magnetics, Vol. 14, Issue 4, 1978, p. 227232. [Search CrossRef]
 Theodoropoulos I., Tatis K., Kladas A., Tegopoulos J. Solid rotor induction machine optimization by using analytical and finite element techniques. Prace Instytutu Elektrotechniki, Vol. 216, 2003, p. 716. [Search CrossRef]
 Markovic M., Perriard Y. An analytical determination of Eddycurrent losses in a configuration with a rotating permanent magnet. IEEE Transactions on Magnetics, Vol. 43, Issue 8, 2007, p. 33803386. [Search CrossRef]
 Markovic M., Perriard Y. Analytical solution for Eddy current losses in a slotless permanent magnet motor: the case of current sheet excitation. IEEE Transactions on Magnetics, Vol. 44, Issue 3, 2008, p. 386393. [Search CrossRef]
 Markovic M., Pfister P.D., Perriard Y. Analytical solution for rotor Eddycurrent losses in a slotless permanentmagnet motor: the case of current layer excitation. International Conference on Electrical Machines and Systems ICEMS, Tokyo, Japan, 2009. [Search CrossRef]
 Chari M. K. V., Bedrosian, G. Electromagnetic field analysis of Eddy currents effects in rotating electrical apparatus and machinery. IEEE Transactions on Magnetics, Vol. 18, Issue 6, 1982, p. 17131715. [Search CrossRef]
 Ergene L. T., Donmezer Y. Improved design for minimizing torque ripples at BLDC motors used in washers. International Review of Electrical Engineering, Vol. 5, Issue 5, 2010, p. 20222032. [Search CrossRef]
 McLachlan N.W. Bessel Functions for Engineers. Clarendon Press, Oxford, 1954. [Search CrossRef]