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.
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
, where
p is the number of pole pairs. The relative permeability distribution is defined in terms of a Fourier series as
where
is the space-varying relative permeability.
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
The general solution to Poisson’s Equation (
2) is determined in each annular region depending on their electromagnetic characteristic and can be expressed as
where:
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,
, 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.,
.
The magnetic vector potential in region
v,
, is described through an infinite series of space harmonics
k. The regions are defined according to
Figure 1, where
and 3 for PM, RS and air-gap, respectively. Rotor reference frame,
, is defined from magnetization axis. Finally, the constants
and
are determined from the boundary conditions, which are defined by the continuity of the tangential field intensity,
, and the normal flux density,
and expressed as
The surface density current,
, 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
and
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
With the tangential component for magnetic flux density calculated from the magnetic vector potential Equations (
3a) and (
3b), Equation (
6) can be rewritten as
where
. 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
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
the constants in (
10) are determined from
Rearranging the right-hand side of (
11) with the aid of (
9) yields
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,
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.
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
, the linear equation systems are obtained with
. 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
From (
14) and (
15), it can be verified that the Fourier coefficients of the magnetic flux density in the non-homogeneous region,
, are related to those of different space harmonic orders of the surrounding regions, i.e.,
and 3. For example,
is expressed in terms of
,
and
, which are the Fourier coefficients of space harmonic orders
,
and
for flux density in region
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.