*Article* **Magnetic Field Analytical Solution for Non-homogeneous Permeability in Retaining Sleeve of a High-Speed Permanent-Magnet Machine**

#### **Gabriel A. Mendonça \*, Thales A. C. Maia and Braz J. Cardoso Filho**

Departamento de Engenharia Elétrica (DEE), Universidade Federal de Minas Gerais (UFMG), Belo Horizonte, MG 31270-901, Brazil<EMAIL_ADDRESS>(T.A.C.M<EMAIL_ADDRESS>(B.J.C.F.) **\*** Correspondence<EMAIL_ADDRESS>Tel.: +55-31-98823-6263

Received: 9 October 2018; Accepted: 7 November 2018; Published: 10 November 2018

**Abstract:** This work presents a novel solution for magnetic field calculation in two-dimensional problems in which one region is defined with space-varying magnetic parameter. The proposed solution extends the well-established Maxwell–Fourier method for calculating magnetic fields in surface-mounted cylindrical high-speed permanent-magnet machines. This contribution is effective to evaluate more realistic magnetic parameters, where measurements of a high-speed permanent-magnet generator prototype indicate saturation in the retaining sleeve due to pole-to-pole leakage flux. The saturation profile is a function of mechanical angle and can be modeled with the aid of a space-varying relative permeability, expressed in terms of a Fourier series. As an example, the presented solution has been applied to a surface-mounted PM machine at no-load condition. Magnetic field calculations show that a simple saturation profile, with low order space-varying permeability in the retaining sleeve significantly affects the magnetic flux density distribution in the air-gap. The analytical solution is confronted with finite-element method, which confirms validity of the proposed methodology.

**Keywords:** analytical model; high-speed; sleeve; non-homogeneous permeability; permanent-magnet

#### **1. Introduction**

An increasing interest in high-speed (HS) machines are noticeable throughout several applications, ranging from distributed energy resources to medical equipment [1,2]. These machines commonly use permanent-magnet (PM) due to high torque per unit volume and low rotor losses [3].

For HS applications, surface-mounted PM machines usually require banding with a sleeve for retaining the PMs against centrifugal forces [4]. One important disadvantage of introducing such a component is an increased effective air-gap, where the added reluctance decreases PM excitation performance [5]. Furthermore, the addition of a conducting retaining sleeve (RS), which works as damper winding or shielding cylinder, reduces PM eddy-current losses in HS operating range [6]. Therefore, an accurate modeling of the RS is crucial for the evaluation of the PM machine overall electromagnetic behavior [7].

The modeling method plays an important role in HS machines electromagnet performance analysis. Numerical methods, such as finite-elements, provide accurate results for magnetic field, but with great computational cost. Alternatively, (semi-)analytical modeling based on the formal solution of Maxwell's equations provides precise calculation of magnetic field distribution in different machine regions. Compared with finite-element method (FEM), this approach is based on strong and limiting assumptions for both geometrical and physical parameters, however, with lower computation time and greater insight into the problem. (Semi-)analytical techniques are considered powerful tools, being vastly used in HS machine early design stage [1–3,6,8,9].

The available literature on (semi-)analytical technique for PM machine evaluation is extensive, with published papers ranging over twenty-five years, to the best of the authors' knowledge [10–12]. Moreover, there is still great effort towards methods with better accuracy and less restricting assumptions [4,13–22] to list a few recently published works. Dubas et al. [23] realized an overview on the existing (semi-)analytical models in Maxwell–Fourier methods (vis., multi-layer models, eigenvalues model, subdomain technique and hybrid models) with the effect of local/global saturation. Further details, advantages and disadvantages of these techniques can be found in [23–25]. Ramakrishnan et al. [13] presented a comprehensive comparison of analytical methods, where a subdomain technique exceeded the others in terms of accuracy with reasonable calculation time.

Furthermore, extending the limitations discussed in [25], recent developments include more realistic geometric structures, such as tooth-tips [14–16,26–28] or rotor eccentricity [27,29]. Chebak et al. [30] and Rahdeh et al. [31] defined more realistic magnetic parameters, but in slotless topologies. Qian et al. [29] and Ortega et al. [32] evaluate several geometries and physical imperfections expected from manufacturing processes. Spranger et al. [33] and Dubas et al. [23,24] have recently developed new techniques to account for finite soft-magnetic permeabilities, respectively: (i) in the multi-layer model using the Cauchy's product theorem [33], and (ii) in the subdomain technique by applying the superposition principle in both directions [23,24]. As discussed by Hannon et al. [34], both methodologies are very effective since they enable the magnetic field calculation in the ferromagnetic material of slotted geometries. The subdomain technique has been improved to consider soft-magnetic material permeability in radial-flux electrical machines, such as switched reluctance machines [21] and spoke-type PM machines [4]. In addition, this technique has been extended to the thermal modeling for the steady-state temperature distribution in rotating electrical machines [22]. According to Pfister et al. [35], the first work on the semi-analytical subdomain technique was published by [36] with the state-of-the-art in the introduction. Moreover, the method based on the Cauchy's product theorem has been improved to account for saturation in stator tooth [19] and a configurable model, adaptable to many machine geometries [20].

However, among the presented references, a small number of authors give a proper realistic evaluation for RS [3,6,9,14,18,28]. Furthermore, only a few of these works evaluate all electromagnetic properties of the retaining sleeve, where relative permeability is different from the air [14,18] and eddy-currents are considered [3,6,9,18,28].

In this paper, a novel formulation is presented to extend a two-dimensional analytical solution from Laplace's and Poisson's equation in polar coordinates. The proposed technique allows the definition of periodic space-varying magnetic parameter, which is defined through Fourier series. A similar approach was used by Sprangers et al. [33], where permeability variation in slotted regions was evaluated directly into the field solution. Such method, based on the convolution theorem, allows the calculation of magnetic field in the tooth region, where interactive solutions can be applied to evaluation of magnetic saturation [19]. The technique presented in this paper, based on well-established Maxwell–Fourier method solutions, evaluates the space-varying permeability by reexamining the boundary conditions. Therefore, the proposed analysis provides an intuitive and fast solution to problems involving magnetic saturation of an annular region and can be extended to other analytical techniques, e.g., the subdomain method for evaluating slotting effects.

Finally, validation of the proposed model is carried out based on a two-pole surface-mounted permanent-magnet generator prototype [37], where experimental data indicated saturation in the retaining sleeve. The geometry and magnetic parameters of the PM generator are approximated, thus providing simplified formulation. The derived model is, then, evaluated with different saturation characteristics and the results are compared with those obtained from FEM.

#### **2. Geometry and Mathematical Formulation**

Throughout the literature, great efforts are invested in developing and improving analytical models focused on the early design stages, where optimization routines play an important role. This work is based on the performance of a HS machine prototype and justifies the reevaluation of common approximations. Thus, the model presented is implemented to include a space-varying permeability to evaluate magnetic saturation effects of the retaining sleeve. Moreover, several other machine aspects not related to the present method are simplified, such as machine geometry and some magnetic parameters, as discussed in this section.

The proposed method is derived according to well-established analytical solution of Maxwell's equations. The saturation phenomenon in the RS is mainly due to pole-to-pole leakage flux, thus changing the RS relative permeability in the tangential direction. In a first approximation, this behavior is considered as independent on loading conditions and the machine is studied evaluating only PM magnetization. Furthermore, stator iron core is assumed to present linear and very large relative permeability, assumed to be infinite. The PM do not demagnetize. All electrical conductivities are assumed to be null. Finally, the magnetic field solutions are derived considering a slotless geometry, as illustrated in Figure 1.

**Figure 1.** A quarter cross section of PM synchronous machine for analytical calculation.

These approximations are necessary to simplify the analytical model and must be used with care. The objective is to derive simpler mathematical formulation, allowing better insight of the proposed method. However, to fully investigate PM machine performance, further research must be done to evaluate more realistic geometry and electromagnetic parameters, such as slotting effects, conducting regions and loading influence.

The remainder of assumptions are necessary to evaluate the two-dimensional analytical problem using Fourier series representation. Magnetic field end effects are neglected and all parameters are periodic in respect to the tangential coordinate. In the RS region, saturation occurs from pole-to-pole leakage flux, where the resultant relative permeability is a function of rotor mechanical angle and varies periodically from pole-to-pole. Therefore, it can be characterized as periodic in the circumferential direction with fundamental space period equals to 2*p*, where *p* is the number of pole pairs. The relative permeability distribution is defined in terms of a Fourier series as

$$
\mu\_r\left(\boldsymbol{\uprho}\right) = \sum\_{k=0,1,2,3,\dots}^{\infty} \hat{\mu}\_{r,k} \cos(2pk\boldsymbol{\uprho}\_r),\tag{1}
$$

where *μr* is the space-varying relative permeability.

The Fourier coefficients *μ*ˆ*r*,*k* can be calculated accordingly for the evaluated phenomenon. For example, to estimate saturation in an annular region, an iterative method could be used to account for the nonlinear magnetic material.

#### *The Poisson's Equation Solution*

Considering previous approximations, the basis for the analytical solution for the air-gap magnetic field is set. Using the magnetic vector potential, and using the Coulomb gauge, the general governing equation for the defined regions is

$$
\nabla^2 \mathbf{A} = -\nabla \times \mathbf{B}\_{\text{rem}}.\tag{2}
$$

The general solution to Poisson's Equation (2) is determined in each annular region depending on their electromagnetic characteristic and can be expressed as

$$A\_z^{(v)}\left(r,\varphi\right) = \sum\_{k=1,3,5,\dots}^{\infty} A\_{z,k}^{(v)}\left(r\right) \sin\left(kp\varphi\_r\right),\tag{3a}$$

where:

$$
\hat{A}\_{z,k}^{(v)}(r) = \left(\mathbb{C}\_k^{(v)} r^{kp} + D\_k^{(v)} r^{-kp} + \hat{A}\_{part,z,k}^{(v)} r\right). \tag{3b}
$$

In the present investigation, where focus is given to characterize saturation in the RS region, PM magnetization pattern and topology does not need to be determined. These aspects are needed, however, to obtain the particular solution, *A*ˆ(*v*) *part*,*z*,*k*, in (3b). Further discussion on this particular solution from different magnetization profiles (e.g., parallel, radial, Halbach, etc.) can be found in [3,6,10,31,36]. Regions where no magnetic field sources are defined, namely air-gap and RS, the Poisson's equation (1) simplifies to the corresponding Laplace's equation. The solution is also obtained from Equations (3a) and (3b), where the particular solution is null, i.e., *A*ˆ(*v*) *part*,*z*,*k* = 0.

The magnetic vector potential in region *v*, *A*(*v*) *z* , is described through an infinite series of space harmonics *k*. The regions are defined according to Figure 1, where *v* = 1, 2 and 3 for PM, RS and air-gap, respectively. Rotor reference frame, *ϕr*, is defined from magnetization axis. Finally, the constants *C*(*v*) *k* and *D*(*v*) *k* are determined from the boundary conditions, which are defined by the continuity of the tangential field intensity, *Hϕ*, and the normal flux density, *Br* and expressed as

$$H\_{\varphi}^{(v+1)}\left(r\_{\upsilon},\varphi\right) - H\_{\varphi}^{(v)}\left(r\_{\upsilon},\varphi\right) = K\_z^{(v)}\left(\varphi\right),\tag{4}$$

$$B\_r^{(v+1)}\left(r\_{v\prime},\varphi\right) - B\_r^{(v)}\left(r\_{v\prime},\varphi\right) = 0.\tag{5}$$

The surface density current, *K*(*v*) *z* , is usually defined along the *z*-axis, being positive in the out-of-the-page direction. Equations (3a) and (3b)–(5) are the basis to the field solution procedure. The analysis is usually simplified, since space harmonics of different orders are orthogonal. Hence, the constants *C*(*v*) *k* and *D*(*v*) *k* are related only to the field source space harmonic component of the *k*-th order, being either from PM or current sources.

The present formulation is based on the reevaluation of such criteria. Reexamining the magnetic flux density tangential component boundary condition for the interfaces of the RS, Equation (4), with the proposed relative permeability in (1), that is

$$\frac{B\_{\varphi}^{(2)}\left(r\_{sl\prime}\,\,\,\rho\right)}{\mu\_{r}(\,\,\rho)\mu\_{0}} = \frac{B\_{\varphi}^{(3)}\left(r\_{sl\prime}\,\,\,\,\rho\right)}{\mu\_{0}},\tag{6}$$

$$\frac{B\_{\varphi}^{(1)}\left(r\_{m\prime}\,\,\,\rho\right) - B\_{\text{rem},\varphi}\left(\,\,\,\rho\right)}{\mu\_0} = \frac{B\_{\varphi}^{(2)}\left(r\_{m\prime}\,\,\,\,\rho\right)}{\mu\_r(\,\,\,\rho)\mu\_0}.\tag{7}$$

With the tangential component for magnetic flux density calculated from the magnetic vector potential Equations (3a) and (3b), Equation (6) can be rewritten as

$$\frac{1}{\mu\_r(\varphi)}\sum\_{k=1,3,5,\dots}^{\infty} \mathcal{B}\_{q,k}^{(2)}(r\_{sl})\sin\left(kp\varphi\_r\right) = \sum\_{k=1,3,5,\dots}^{\infty} \mathcal{B}\_{q,k}^{(3)}(r\_{sl})\sin\left(kp\varphi\_r\right),\tag{8}$$

where *B*ˆ(*v*) *ϕ*,*k* (*r*) = *kp C*(*v*) *k rkp*−1 − *D*(*v*) *k r*−*kp*−1 + *∂ A*ˆ(*v*) *part*,*z*,*kr ∂r* . The solution, with relative permeability defined with (1), includes multiple summation terms. For generalization purposes, these summation

indexes are defined with different variables, *k*, *m* and *n*. Thus, from (8), we have

$$\sum\_{k=1,3,5,\dots}^{\infty} \mathcal{B}\_{q,k}^{(2)}(r\_{sl}) \sin\left(kp \, q\_{l'}\right) = \sum\_{m=0,1,2,3,\dots}^{\infty} \mu\_{r,m} \cos\left(2mp \, q\_{l'}\right) \sum\_{n=1,3,5,\dots}^{\infty} \mathcal{B}\_{q,n}^{(3)}(r\_{sl}) \sin\left(np \, q\_{l'}\right) \,. \tag{9}$$

The calculation of the constants and proper interaction of the different space harmonic indexes are found from the definition of Fourier series. Therefore, if we establish

$$f\left(\boldsymbol{\varrho}\_{r}\right) = \sum\_{k=1,3,5,\dots}^{\infty} \mathcal{B}\_{\boldsymbol{\varrho},k}^{(2)} \sin\left(k \boldsymbol{\varrho} \boldsymbol{\varrho}\_{r}\right),\tag{10}$$

the constants in (10) are determined from

$$
\hat{B}\_{\varphi,k}^{(2)} = \frac{1}{\pi t} \int\_0^{2\pi} f\left(\varphi\_r\right) \sin\left(kp\varphi\_r\right) d\varphi. \tag{11}
$$

Rearranging the right-hand side of (11) with the aid of (9) yields

$$\hat{B}\_{q,\lambda}^{(2)}(r\_{sl}) = \sum\_{m=0,1,2,3,\dots}^{\infty} \sum\_{n=1,3,5,\dots}^{\infty} \frac{\mu\_{r,m} \mathcal{B}\_{q,n}^{(3)}(r\_{sl})}{\pi} \int\_{0}^{2\pi} \sin\left(np\,\rho\_{l}\right) \cos\left(2mp\,\rho\_{l}\right) \sin\left(kp\,\rho\_{l}\right) \,d\rho,\tag{12}$$

where the space harmonic index *k* can be thought as a reference index. The integration on the right-hand side of (12) will not be null for specific values of *k*, *m* and *n*.The same reasoning is used for (7), that is,

$$\begin{split} \mathcal{B}^{(2)}\_{\boldsymbol{\varrho},\boldsymbol{k}}(r\_{\boldsymbol{m}}) &= \sum\_{m=-0,1,2,3,\ldots}^{\infty} \sum\_{\substack{i=1,3,5,\ldots\\m\end{pmatrix}}^{\infty} \frac{\hat{\mu}\_{r,m} \mathcal{B}^{(1)}\_{\boldsymbol{\varrho},\boldsymbol{k}}(r\_{\boldsymbol{m}})}{\pi} \int\_{0}^{2\pi} \sin\left(ip\,\rho\_{\boldsymbol{\vargamma}}\right) \cos\left(2mp\,\rho\_{\boldsymbol{\vargamma}}\right) \sin\left(kp\,\rho\_{\boldsymbol{\vargamma}}\right) d\boldsymbol{\upgamma} \\ &- \sum\_{m=0,1,2,3,\ldots}^{\infty} \sum\_{\substack{i=1,3,5,\ldots\\m\end{pmatrix}}^{\infty} \frac{\hat{\mu}\_{r,m} \hat{\mu}\_{\boldsymbol{\upgamma},\boldsymbol{q},\boldsymbol{l}}}{\pi} \int\_{0}^{2\pi} \sin\left(kp\,\rho\_{\boldsymbol{\vargamma}}\right) \cos\left(2mp\,\rho\_{\boldsymbol{\vargamma}}\right) \sin\left(kp\,\rho\_{\boldsymbol{\vargamma}}\right) d\boldsymbol{\upgamma} .\end{split} \tag{13}$$

The definition of a relative permeability from Fourier series expansion results in an interaction of magnetic fields of different orders. The integrations in (12) and (13) will be different from zero only for certain values of *i*, *k*, *l*, *m* and *n*. The trigonometric transformations needed to evaluate integration terms in Equations (12) and (13) are well known and will not be reproduced here. Summarizing the results, Table 1 presents the relations between the indexes which results non-zero integration values in Equation (12). Though there are different indexes relations in Equation (13), the results can be easily extended by replacing *n* with *i* or *l*.

**Table 1.** Space harmonic index relations for non-zero integral solutions.

One first conclusion derived from Table 1 is that, assuming homogeneous and constant RS magnetic parameter, the solution is determined from interaction of magnetic fields with the same space harmonic order. In such case, for *m* = 0, the linear equation systems are obtained with *k* = *i* = *l* = *n*. In addition, it is important to note that the results from the index relations must be added up accordingly with the summation notations. Therefore, from the index relations in Table 1, the result for Equations (12) and (13) can be derived as

$$\mathcal{B}^{(2)}\_{\varphi,k}(r\_{sl}) = \sum\_{m=0,1,2,3,\dots}^{\infty} \frac{\hat{\mu}\_{r,m}}{2} \left[ \mathcal{B}^{(3)}\_{\varphi,k-2m}(r\_{sl}) + \mathcal{B}^{(3)}\_{\varphi,k+2m}(r\_{sl}) - \mathcal{B}^{(3)}\_{\varphi,2m-k}(r\_{sl}) \right],\tag{14}$$

$$\begin{split} \mathcal{B}\_{\boldsymbol{\varphi},k}^{(2)}(\boldsymbol{r}\_{m}) &= \sum\_{m=0,1,2,3,\dots}^{\infty} \frac{\hat{\mu}\_{\boldsymbol{r},m}}{2} \left[ \mathcal{B}\_{\boldsymbol{\varphi},k-2m}^{(1)}(\boldsymbol{r}\_{m}) + \mathcal{B}\_{\boldsymbol{\varphi},k+2m}^{(1)}(\boldsymbol{r}\_{m}) - \mathcal{B}\_{\boldsymbol{\varphi},2m-k}^{(1)}(\boldsymbol{r}\_{m}) \right] \\ &- \sum\_{m=0,1,2,3,\dots}^{\infty} \frac{\hat{\mu}\_{\boldsymbol{r},m}}{2} \left[ \mathcal{B}\_{\boldsymbol{\varphi}rm,\boldsymbol{\varphi},k-2m} + \mathcal{B}\_{\boldsymbol{\varphi}rm,\boldsymbol{\varphi},k+2m} - \mathcal{B}\_{\boldsymbol{\varphi}rm,\boldsymbol{\varphi},2m-k} \right]. \end{split} \tag{15}$$

From (14) and (15), it can be verified that the Fourier coefficients of the magnetic flux density in the non-homogeneous region, *v* = 2, are related to those of different space harmonic orders of the surrounding regions, i.e., *v* = 1 and 3. For example, *B*ˆ(2) *ϕ*,*k* is expressed in terms of *B*ˆ(3) *ϕ*,*k*−2*m*, *B*ˆ(3) *ϕ*,*k*+2*m* and *B*ˆ(3) *ϕ*,2*m*−*k*, which are the Fourier coefficients of space harmonic orders (*k* − 2*m*), (*k* + 2*m*) and (2*m* − *k*) for flux density in region *v* = 3 and where *m* is the space harmonic order of the relative permeability Fourier series. These relations will be further detailed in the application example (Section 4).

The mathematical formulation presented in this section is general and can be applied to several PM machine configurations. However, to help illustrate the proposed methodology and to show the importance of evaluating the RS saturation, the following sections will discuss its application to a HS machine prototype.

#### **3. Evaluated Prototype**

The machine used for verifying the proposed methodology is part of a micro-compressed air energy storage (micro-CAES) system [37]. It was built from a conventional turbocharger system adapted to operate only with air [37]. Due to its high speeds, the power density is considerably higher. It is suitable for small scale applications, such as operating with distributed generation to improve power supply reliability. The prototype is illustrated in Figure 2, where Figure 2a shows the dissembled system, and Figure 2b illustrates the mounted system.

**Figure 2.** Low cost high-speed micro-CAES system: (**a**) disassembled components; and (**b**) system in operation [37].

In such application, generator design is carried out based on some overall restrictions, such as reduced size and cost. A single high-speed shaft eliminates the need of gear-box, thus improving system reliability [2]. However, a high-speed drive system design is critical, with difficulties associated with generator topology, material selection, manufacturing procedure, bearings types and many others. All of these concerns for the evaluated machine are addressed by Maia et al. [37].

The constructed machine consists of a two-pole surface-mounted PM generator with parallel magnetization. In order to fulfill the banding requirement, a stainless steel AISI 310 sleeve was adopted for its additional properties. The high electrical conductivity is expected to help shielding the magnets and rotor iron from high-frequency magnetic fields [3], thus reducing rotor losses. Moreover, behaving like a damper winding, it could help improving transient stability. The stator, manufactured with electric steel with 0.35 mm lamination thickness, was designed with full-pitch concentrated winding, as pictured in Figure 3a. The assembled prototype is presented in Figure 3b and final machine characteristics are shown in Table 2.

**Figure 3.** PM generator: (**a**) schematic view; and (**b**) assembled prototype stator.

**Table 2.** Prototype PM machine characteristics.

Maia et al. [37,38] successfully evaluated the electromechanical design of a HS generator directly driven by a micro-CAES system. Additionally, experimental measurements carried out helped reevaluate the assumptions usually made to simplify the mathematical formulation. Using flux sensing coils, as indicated in Figure 3a, measurements pointed that the RS had higher magnetic permeability than expected from the material specification. The RS machining, needed to achieve proper mechanical tolerances, changed the relative permeability to approximately 300. With the adjusted values, numerical simulations using FEMM [39] were performed. The results, compared with measurements, are presented in Figure 4.

**Figure 4.** Measured and simulated results for induced voltage in exploring coils.

Discarding any saturation effects in the RS region, the higher relative permeability increases the pole-to-pole leakage flux, reducing significantly the no-load stator flux linkage and electromotive force [5]. Experimental no-load measurements did not match the results from FEM simulations which assumed linear RS relative permeability. Reevaluating the numerical simulation with nonlinear relative permeability for the RS region yields accurate results compared with measurements, as illustrated in Figure 5. The saturated RS provides a preferential path for the magnetic flux in the magnetization axis, diminishing the pole-to-pole leakage flux. The saturation curve for the stainless-steel RS was based on [40] and adjusted to provide the unsaturated relative permeability of approximately 300. The relative permeability profile in the RS for the PM generator at no-load condition is illustrated in Figure 6.

**Figure 5.** Measured and simulated results for no-load voltage.

**Figure 6.** RS space-varying relative permeability profile.

Another aspect observed from the manufactured system is that, though not measured, cogging torque was perceived by manually rotating the machine shaft. This is an interesting result because it also leads to the assumption that saturation was taking place in the RS. A linear magnetic parameter, with this type of PM excitation and slot geometry, would result in a constant torque, without any harmonic content.

The interaction between magnetic flux produced by the PM and the stator non-uniform structure governs the cogging torque phenomenon [29]. This relationship can be clearly visualized through the complex permeance modeling method [8,13,41,42]. From closed form solutions for calculating the cogging torque using conformal mapping [41], no cogging torque would be expected for the studied motor. The rotor was manufactured with a two-pole NdFeB ring shaped with parallel magnetization, which produces a sinusoidal air-gap magnetic flux density with no space harmonic content, and the slotting configuration results in a constant torque profile, with no oscillation. Using the revised relative permeability for the RS, simulations using FEMM [39] were performed to illustrate the effect of saturation in the cogging torque.

To allow the evaluation of the stator core saturation effect on the cogging torque, the simulation was made using different remanence values for the PM. Figure 7a,b illustrate the results for different remanence values, *B*ˆ*rem* = 1.2 T and *B*ˆ*rem* = 1.3 T, respectively. For each PM remanent magnetization value, simulation evaluated effects of using linear or nonlinear permeability of RS, stator and rotor core.

**Figure 7.** Cogging torque FEMM simulation with: (**a**) *B*ˆ*rem* = 1.2 T; and (**b**) *B*ˆ*rem* = 1.3 T.

#### **4. Application Example Results and Validation**

The validation was supported by problems simplification. Machine slotting effect was neglected and many magnetic parameters were approximated, as discussed in Section 2. This approach was chosen to help clarify some key aspects in the new formulation proposed. Furthermore, the relative permeability was defined only with a constant and one fundamental space-varying component, i.e., *μr*(*ϕ*) = *μ*ˆ*r*,0 + *μ*ˆ*r*,1 cos (2*ϕr*). This approximation is justified here to help evaluate the interaction between space harmonic of different orders, showing that even low order spatial harmonics in the RS can significantly affect the air-gap density flux distribution. In addition, results from the analytical method were validated from FEM results.

The geometrical parameters used in the mathematical formulation are based on the PM machine structure illustrated in Figure 1 with values given in Table 3. The rotor type allows one more simplification without loss of generality. The rotor iron core region was neglected and the rotor structure is composed of full cylindrical PM and the RS, as illustrated in Figure 8. This reduces one boundary condition and, as a consequence, one linear equation (i.e., *D*(1) *k* = 0, from (3b)). In addition, the machine PM configuration simplifies the Poisson's equation since the curl of the parallel magnetization vector is null. Therefore, the radial and tangential components of the magnetization

were defined only with a fundamental space harmonic, i.e., **B***rem* = *B*ˆ*rem* cos (*ϕr*) **a***r* − *B*ˆ*rem* sin (*ϕr*) **a***ϕ*, where **a***r* and **a***ϕ* are the unit vectors in the radial and tangential direction.

**Table 3.** Parameters for analytical model of slotless geometry.

**Figure 8.** A quarter cross section of the PM synchronous machine with a segmented RS region for FEMM simulation.

Combining Equations (3a) and (3b) with the boundary conditions (4) and (5) establishes a set of linear equations. However, the RS boundary conditions were reevaluated with the procedure developed in Section 2. First, for the studied machine with only fundamental magnetization vector, the index *l* in (13) is set as equal to one. Table 4 presents the space harmonic indexes *m* and *n*, as defined in (12), which satisfies the conditions in Table 1. In addition, the superscripts in Table 4 relate the conditions in Table 1 that are satisfied.

**Table 4.** Space harmonic index relation for an evaluated PM machine.

With the presented space harmonic indexes combinations, the system of linear equations can be defined. First, three boundary conditions used to derive the linear equation system remain unchanged. For the normal flux density continuity at *r* = *rm* and *r* = *rsl* and for the tangential field intensity at *r* = *rs*, the expressions can be derived for any space harmonic indexes as follows:

$$r\_s^{k-1} \mathbb{C}\_k^{(3)} - r\_s^{-k-1} D\_k^{(3)} = 0,\tag{16}$$

$$r\_{sl}^{k-1} \mathbf{C}\_k^{(3)} + r\_{sl}^{-k-1} D\_k^{(3)} - r\_{sl}^{k-1} \mathbf{C}\_k^{(2)} - r\_{sl}^{-k-1} D\_k^{(2)} = 0,\tag{17}$$

$$r\_m^{k-1} \mathbb{C}\_k^{(2)} + r\_m^{-k-1} D\_k^{(2)} - r\_m^{k-1} \mathbb{C}\_k^{(1)} = 0. \tag{18}$$

Furthermore, from the space harmonic relations defined in Table 4, the tangential field intensity continuity at *r* = *rm* and *r* = *rsl* can be expressed as

• For *k* = 1:

$$\begin{aligned} \left(\hat{\mu}\_{r,0} - \frac{1}{2}\hat{\mu}\_{r,1}\right)\mathcal{C}\_1^{(3)} - \left(\hat{\mu}\_{r,0} - \frac{1}{2}\hat{\mu}\_{r,1}\right)r\_{sl}^{-2}D\_1^{(3)} - \mathcal{C}\_1^{(2)} + r\_{sl}^{-2}D\_1^{(2)}\\ + \frac{3}{2}\hat{\mu}\_{r,1}r\_{sl}^2\mathcal{C}\_3^{(3)} - \frac{3}{2}\hat{\mu}\_{r,1}r\_{sl}^{-4}D\_3^{(3)}, \end{aligned} \tag{19}$$

$$\begin{split} \mathbf{C}\_{1}^{(2)} - r\_{m}^{-2} \mathbf{D}\_{1}^{(2)} - \left( \widehat{\mu}\_{r,0} - \frac{1}{2} \widehat{\mu}\_{r,1} \right) \mathbf{C}\_{1}^{(1)} - \frac{3}{2} \widehat{\mu}\_{r,1} r\_{m}^{2} \mathbf{C}\_{3}^{(1)} \\ = -\left( \widehat{\mu}\_{r,0} - \frac{1}{2} \widehat{\mu}\_{r,1} \right) \mathbf{B}\_{\text{rem}}. \end{split} \tag{20}$$

• For *k* = 3:

$$\begin{split} & \frac{1}{6} \hat{\mu}\_{r,1} \mathbb{C}\_{1}^{(3)} - \frac{1}{6} \hat{\mu}\_{r,1} r\_{sl}^{-2} D\_{1}^{(3)} + \hat{\mu}\_{r,0} r\_{sl}^{2} \mathbb{C}\_{3}^{(3)} - \hat{\mu}\_{r,0} r\_{sl}^{-4} D\_{3}^{(3)} \\ & - r\_{sl}^{2} \mathbb{C}\_{3}^{(2)} + r\_{sl}^{-4} D\_{3}^{(2)} + \frac{5}{6} \hat{\mu}\_{r,1} r\_{sl}^{4} \mathbb{C}\_{5}^{(3)} - \frac{5}{6} \hat{\mu}\_{r,1} r\_{sl}^{-6} D\_{5}^{(3)} = 0, \end{split} \tag{21}$$

$$\begin{split} -\frac{1}{6}\hat{\mu}\_{r,1}\mathbf{C}\_{1}^{(1)} + r\_{m}^{2}\mathbf{C}\_{3}^{(2)} - r\_{m}^{-4}D\_{3}^{(2)} - \hat{\mu}\_{r,0}r\_{m}^{2}\mathbf{C}\_{3}^{(1)} - \frac{5}{6}\hat{\mu}\_{r,1}r\_{m}^{4}\mathbf{C}\_{5}^{(1)} \\ = \frac{1}{6}\hat{\mu}\_{r,1}\hat{\mu}\_{\text{rem}}. \end{split} \tag{22}$$

• For *k* ≥ 5:

$$\begin{split} \frac{k-2}{2k} \hat{\mu}\_{r,1} r\_{sl}^{k-3} \mathbf{C}\_{k-2}^{(3)} - \frac{k-2}{2k} \hat{\mu}\_{r,1} r\_{sl}^{-k+1} \mathbf{D}\_{k-2}^{(3)} \\ + \hat{\mu}\_{r,0} r\_{sl}^{k-1} \mathbf{C}\_{k}^{(3)} - \hat{\mu}\_{r,0} r\_{sl}^{-k-1} \mathbf{D}\_{k}^{(3)} - r\_{sl}^{k-1} \mathbf{C}\_{k}^{(2)} + r\_{sl}^{-k-1} \mathbf{C}\_{k}^{(2)} \\ + \frac{k+2}{2k} \hat{\mu}\_{r,1} r\_{sl}^{k+1} \mathbf{C}\_{k+2}^{(3)} - \frac{k+2}{2k} \hat{\mu}\_{r,1} r\_{sl}^{k-3} \mathbf{D}\_{k+2}^{(3)} = 0, \end{split} \tag{23}$$

$$\begin{split} -\frac{k-2}{2k} \hat{\mu}\_{r,1} r\_{um}^{k-3} \mathbf{C}\_{k-2}^{(1)} + r\_{um}^{k-1} \mathbf{C}\_{k}^{(2)} - r\_{m}^{-k-1} D\_{k}^{(2)} - \hat{\mu}\_{r,0} r\_{m}^{k-1} \mathbf{C}\_{k}^{(1)} - \mathbf{c} \end{split} \tag{24}$$

$$\frac{k+2}{2k} \hat{\mu}\_{r,1} r\_{m}^{k+1} \mathbf{C}\_{k+2}^{(1)} = 0. \tag{25}$$

The proposed formulation requires different space harmonics orders to be solved integrally, which differs from the usual procedure. In addition, the analytical method was validated based on results from numerical simulations. The FEM analysis was set with the same premises used in the analytical model. For the space-varying relative permeability in the RS, the region was discretized. To obtain the needed accuracy, 80 regions were defined to characterize this space-varying parameter in the RS, as illustrated in Figure 8.

The simulations were carried out with several values for the RS relative permeability, considering both homogeneous and space-varying parameter. Table 5 summarizes the obtained result where the RS was modeled based on the evaluated PM machine. The results show the radial component of the magnetic flux density in the air-gap. Moreover, to illustrate the analytical model application with different permeability values, the model validation was done with similar parameters but with different magnitude orders, as illustrated in Figure 9 with the results presented in Tables 6 and 7.

**Table 5.** Simulation result for *B*ˆ(3) *r*,*k* (*T*): relative permeability based on measurements.

**Figure 9.** Different relative permeability profile for method validation.

**Table 6.** Simulation result for *B*ˆ(3) *r*,*k* (*T*): relative permeability based on measurements.

**Table 7.** Simulation result for *B*ˆ(3) *r*,*k* (*T*): relative permeability based on measurements.

The results show a very good agreement with those obtained from the FEM simulations. As expected from the derived equations, the results show great interaction among space harmonics of different orders. Tables 5–7 present the results with different harmonic spectrum considered for

the analytical method, ranging from *k* = 1 to *kmax*. In addition, the results are illustrated in Figure 10, which depicts the magnetic flux density waveforms from both methods in different regions (viz., the average radius of the air-gap, RS and PM). In these figures, results are for a maximum number of harmonics equal to *kmax* = 100.

**Figure 10.** Analytical and numarical simulation comparison at different regions: (**a**) *μr* = 16 + 15 cos(2*ϕr*); (**b**) *μr* = 151 + 150 cos(2*ϕr*) ;and (**c**) *μr* = 1501 + 1500 cos(2*ϕr*).

As can be appreciated from the results in Tables 5–7, the analyital method accuracy is dependent on the maximum order of space harmonics, i.e., *kmax*. For a higher relative permeability case presented in Table 6 and in Figure 10c, the analytical method converged (with an error smaller than 5%) only for space harmonic content evaluated up to the thirtieth order. In addition, the relative permeability evaluated with only one harmonic content leads to a air-gap magnetic flux density with significant spacial harmonic content.

Furthermore, if only the fundamental component would exist in a system, with *μr*,1 = 0 in the evaluated PM machine, no harmonic content would be observed using either techniques, FEM or analytical method. This result is expected since the PM topology with parallel magnetization would result in only fundamental magnetic fields. In addition, as expected from the no-load voltage measurements in the prototype in comparison with numerical results, higher relative permeability in the RS region increases the pole-to-pole leakage flux significantly. This behavior would diminish air-gap magnetic fields and stator linkage flux if no saturation effect took place in the RS. Therefore, the space-varying relative permeability helps to maintain the magnetization level.

#### **5. Conclusions**

In this work, a new analytical model is presented that incorporates non-homogeneous magnetic parameters. The novel procedure evaluates an annular region defined with a space-varying relative permeability, where iterative methods could be used to study saturation phenomenon. This new solution is needed because experimental analysis pointed out a higher order space harmonics spectrum for the calculated magnetic field. In the assessed PM machine, linear models provide only fundamental components of the magnetic fields. Such assumption would lead to erroneous conclusions.

The novel analytical method was validated with results obtained from finite-element analysis. Seeking a clearer and direct way to analyze the derived results, the PM machine model and parameters not related to the proposed solution were simplified. In this case, linear system equations were derived illustrating the proposed methodology. Comparison with numerical solutions for magnetic flux density at different PM machine regions showed very good agreement for both amplitude and phase.

The proposed technique is derived from well-established analytical solution methodology, where other researched improvements can be extended with it. For example, evaluation of slotting effect can be easily carried out with a subdomain model. This will help suppress known analytical method drawbacks, such as the need to define linear material properties to solve the problem. Thus, improved accuracy in magnetic field calculation, aligned with fast computation time, can extend analytical method application capabilities throughout PM machine and drive design procedures.

**Author Contributions:** Supervision, B.J.C.F.; Writing—original draft, G.A.M.; Writing—review & editing, T.A.C.M.

**Acknowledgments:** This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) - Finance Code 001, the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and the Universidade Federal de Minas Gerais (UFMG).

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).