Voltage Sources in 2 D Fourier-Based Analytical Models of Electric Machines

The importance of extensive optimizations during the design of electric machines entails a need for fast and accurate simulation tools. For that reason, Fourier-based analytical models have gained a lot of popularity. The problem, however, is that these models typically require a current density as input.This is in contrast with the fact that the greatmajority ofmodern drive trains are powered with the help of a pulse-width modulated voltage-source inverter. To overcome that mismatch, this paper presents a coupling of classical Fourier-based models with the equation for the terminal voltage of an electric machine, a technique that is well known in finite-element modeling but has not yet been translated to Fourier-based analytical models. Both a very general discussion of the technique and a specific example are discussed. The presented work is validated with the help of a finite-element model. A very good accuracy is obtained.


Introduction
With evermore strict demands on the performance of electric machines, the importance of optimizations during early design stages is growing.Typically, such optimizations have a very large design space.In order to limit the associated computational burden, fast machine models are required.Therefore, often very simple models are used [1,2].However, typically, these models require a lot of simplifications or experimental parameters.One class of models that combines a high level of accuracy and a low computational time is the class of Fourier-based analytical models.Moreover, these Fourier-based analytical models do not require any experimental parameters.It is therefore no surprise that the interest in such models is very high [3][4][5][6][7][8][9][10][11].
The Fourier-based analytical models that are presented in literature require a current density as input.The problem is that nowadays most electric drives are powered with the help of a pulse-width modulated voltage-source inverter.To overcome that mismatch, this work extends the magnetic calculations of Fourier-based models with the equation for the terminal voltage of an electric machine.By doing so, it is possible to directly account for a voltage source instead of the classical approach of imposing a current density.The technique of combining magnetic calculations with the equation for the terminal voltage has already proven its worth in finite-element models [12].However, to date, it has not yet been translated to Fourier-based analytical models.Moreover, despite its apparent simplicity, a general formulation of the technique is not obvious.For those reasons, this work presents both a general discussion and a specific example on how to couple classical Fourier-based models with the equation for the terminal voltage.
The presented work consists of three major parts.Firstly, a very general discussion on the technique of coupling magnetic calculations with the equation for the terminal voltage of electric machines is presented.This discussion is spread over Sections 2-5.Section 2 introduces the magnetic calculations of classic Fourier-based models.In Sections 3 and 4, the equation for the terminal voltage is rewritten so that it can be used in Fourier-based analytical models.These two aspects are combined in Section 5, resulting in a model that directly accounts for the terminal voltage of electric machines.
In the second part of this work, the theoretical discussion of part one is concretized with the help of an example.
The machine that is considered is a rotating, radial-flux permanent-magnet synchronous machine (PMSM) for highspeed applications.This machine is discussed in Section 6.
The same machine topology is used in the third part to validate the presented work with the help of a finite-element model.This is done in Section 7. Section 8 concludes the work.

Fourier-Based Analytical Models
By means of introduction, this section briefly discusses Fourier-based analytical models.The goal is to present a set of equations that defines the magnetic field as a function of the geometry, the remanent magnetic flux density, and the current density.
Typically, Fourier-based models use a potential formulation, that is, the magnetic vector potential (MVP) or the Magnetic Scalar Potential (MSP), to rewrite Maxwell's equations in the form of a second-order partial-differential equation.Because of the fact that it is the most generally valid formulation, the focus in this work is on the MVP.
The magnetic vector potential, which is denoted by A, is defined through its curl: where B is the magnetic flux density.
The magnetic field and therefore the MVP depend on time and space.Fourier-based models rely on periodicities in these time and spatial dependencies to formulate a solution for A. While time is simply denoted by , a coordinate system has to be chosen to describe the spatial dependency.This choice usually depends on the shape of the studied problem.For the sake of generality, an arbitrary coordinate system (, , ) is defined so that the spatial periodicity is along the -axis.The -direction is referred to as the parallel direction.The -direction is referred to as the normal direction and the -direction is the longitudinal direction.Note that, in practice, most authors use a Cartesian coordinate system ( = ,  = ,  = ), a cylindrical coordinate system ( = ,  = ,  = ), or an axisymmetric coordinate system ( = ,  = ,  = ).More information on the use of these coordinate systems can be found in [6].
For simplicity, in this work, it is assumed that the studied problem is invariant along the longitudinal axis; that is, the problem can be regarded in two dimensions.This approximation is very generally used and implies that the MVP's only nonzero component is the one along the -direction [3][4][5][6][7][8][9].Evidently, this nonzero component is independent from the invariant direction.Mathematically, this assumption implies the following: In rotational machines with a radial flux, for example, usually a cylindrical system with the -direction along the machine's axis is used.The problem is then often assumed to be invariant along the machine's axis, that is, along the -direction.

Governing Equation.
As already mentioned, the MVP is used to rewrite Maxwell's equations in the form of a differential equation.Assuming quasistatic conditions, this governing equation is written as [3,4] where the time-derivative term accounts for eddy-currents and the two terms in the right-hand side of (3) account for the sources, that is, externally imposed current densities and remanent magnetic flux densities. with where  and  are the time-and spatial-harmonic orders,  (]) 0 is the starting angle of subdomain ],  (]) , and  (]) , are integration constants,  (]) , () and  (]) , () are -dependent functions that are defined by the governing equation, and   and   are the time and spatial periods.ℎ (]) , () is the particular solution of the governing equation, which is either zero or dependent on the source terms in subdomain ] (J (]) ext and B (]) rem ).Note that in practice the infinite summations have to be truncated.The highest time-and spatial-harmonic orders that are considered are ℎ  and ℎ  , respectively.

Boundary Conditions
. The next step is to link the solutions in the different subdomains back together.This is done by imposing physical boundary conditions, that is, Ampère's law and Gauss law for magnetism.On the boundary between subdomains ] and ] + 1, these conditions are written as follows: where H (]) is the magnetic field strength in subdomain ] and K (]) is the current density on the boundary between subdomains ] and ] + 1.Using the definition of the MVP and the magnetic constitutive relation, these boundary conditions can be written in terms of the MVP [3-9].

System of Equations.
In a final step, a system of equations is constructed by combining the equation for the MVP (4) with the boundary conditions (6a) and (6b).As shown in [3], the system can be solved separately for every time-harmonic order ().In matrix form, the result is then written as follows: where [  ] is a 2(2ℎ  + 1) ] × 1 matrix, every row of which refers to an integration constant.Analogously, [  ] is a (2ℎ  + 1) ] × 1 matrix that contains all of the spatial-harmonic coefficients of the remanent magnetic induction in each of the subdomains.[  ] contains the current densities in each of the subdomains.Although some authors account for the spatial dependency of the current density in the subdomains [11,13], usually the current density is assumed to be constant in a subdomain.For simplicity reasons, this assumption is also adopted in this work.[  ] is then a  ] × 1 matrix, with each row referring to the current density in a subdomain. [ and 2(2ℎ  + 1) ] × (2ℎ  + 1) ] .These matrices depend on the machine's geometry and are therefore not regarded here.
The above implies that the system contains 2(2ℎ  + 1) ] unknown integration constants and an equal amount of boundary conditions.
By solving (7) for every time-harmonic order, the MVP is uniquely defined in every subdomain.This implies that the magnetic field in the studied machine is known.However, if the current densities, that is, [  ], are unknown, the system is underdetermined and extra  ] equations are required.In the following sections, these equations will be derived from the equation for the terminal voltage of an electric machine.
2.6.Assumptions.The analytical approach described in the above requires some basic assumptions.Primarily, the machine is assumed to operate in steady state.This assumption is necessary in order to impose a time periodicity, as explained in Section 2. A second assumption is that the problem can be regarded in 2 dimensions.Although this assumption is not strictly necessary, it greatly simplifies the calculus.Thirdly, the situation is assumed to be quasistatic.Again, this approximation results in a reduced computational complexity.For that same reason, the externally imposed current density is assumed to be spatially constant in every current-carrying region, for example, a slot.Finally, the machine's soft-magnetic material is assumed to be infinitely permeable.This assumption is mandatory to analytically solve the governing equation.
These five assumptions are listed as follows: (i) Steady-state operation.
(ii) 2D approximation of the problem.
(iv) Uniform current density in every subdomain.
(v) Infinite permeability of the iron.

Terminal Voltage and Current Density
The previous section introduced the calculation of the magnetic field of an electric machine with the help of the Fourierbased modeling technique.As can be seen from ( 3), one of the inputs of such a model is the current density.In order to input the terminal voltages instead, this section links the current densities of every subdomain to the applied terminal voltages.The discussion starts from the classical equation for the terminal voltage of an arbitrary coil  in an electric machine: The flux coupled with a coil  ( () tot ()) can be divided in a component related to the active part of the coil ( () ()) and a component related to the end-windings ( ()  ew ()).The above then results in the following: Note that it is assumed that every coil has the same ohmic resistance  and the same end-windings inductance  ew .These values can be obtained with classical formulas, such as the ones found in [14].The functions in (9) can be written in terms of their Fourier series: where  is the machine's mechanical pulsation, which equals 2/  .The above can be rewritten for every time-harmonic order  separately: Imply that the current's th harmonic order can be calculated as follows: The amplitude of the th harmonic order of the current density in a subdomain ], for example, a slot, can now be calculated as follows: In (13),  (,]) is the amount of conductors that coil  has in subdomain ] and   is the amount of coils in the machine.
(]) is the surface of subdomain ].It can easily be reasoned that substitution of (13) in the governing equation will allow for directly accounting for the terminal voltage as a source.However, Ψ ()  has to be calculated from the magnetic field, that is, from the solution of the governing equation.This prevents a direct coupling between the calculation of the magnetic field and the equation for the terminal voltage.To overcome that problem, Ψ ()   is rewritten in terms of the magnetic vector potential in the following section.

Flux Linkage
The goal of this section is to express the flux related to the active part of an arbitrary coil  as a function of the magnetic vector potential.In a first step, the flux coupled with a single turn of the coil is derived; in a second step, the flux coupled with the entire coil is regarded.

Flux Coupled with a Single
Turn.The physical flux through a single turn  of coil  is calculated as the integration of the flux density over a surface spanned by that turn: where the definition of the MVP and Stokes theorem were used.  is the surface of the turn and   is the boundary of that surface.Since the magnetic vector potential is assumed to only have an -component, the integration of A along the turn's contour will only be nonzero along the -direction.Noting that the MVP is independent of , this implies that the integration along the contour of the turn can be rewritten as follows: where   is the longitudinal length of the studied machine.] + is the subdomain in which the direction of the integration is along the positive -axis; A and s then have the same direction and sense.] − represents the subdomain in which the coil returns; the integration direction is opposed to the -axis and thus to A.
+ and   − are the normal positions of turn  and   + and   − are the tangential positions of turn .
The direction of the integration is chosen so that it corresponds to the reference direction of the current.This implies that the going conductor of turn  is located in subdomain ] + and the returning conductor is located in subdomain ] − .Equation (15) can now be rewritten as follows: where  (])  is 1 in the subdomain that contains the going conductor of , −1 in the subdomain that contains the returning conductor, and 0 in the other subdomains.
Usually, the exact position of the turn cannot be determined; that is, ( (])  ,  (])  ) is unknown.For that reason, most authors either choose an arbitrary position or use the average MVP in the considered subdomain.Doing so results in one MVP value for every subdomain; this value will be referred to as  (]) () in the following.The flux coupled with a single turn  of coil  can then be written as follows: (]) () . (17)

Flux
Coupled with a Coil.The flux coupled with coil  is calculated by summing the fluxes coupled with each of its turns.From (17), it can be written that where which implies that  (,]) equals  (,]) if subdomain ] contains going conductors of ,  (,]) equals − (,]) if subdomain ] contains returning conductors of , and  (,]) equals 0 if subdomain ] does not contain any conductors of coil .
The th time-harmonic coefficient of () () can now be written as follows: where  (])  is the th time-harmonic coefficient of  (]) ().Since  (]) () is a direct function of the MVP in subdomain ], which in turn is determined by the integration constants and the source terms, (20) can be written in matrix form as follows: where [ () ] is a 1 ×  ] matrix describing the winding configuration of coil ; that is, [ () ] 1,] =  (,]) .[  ] greatly depends on the studied geometry.For that reason, it will not be regarded here.However, [ 4   ], [ 5   ], and [ 6   ] can directly be derived from the MVP.It can easily be seen that the above effectively expresses the flux linkage of coil  in terms of the machine's geometry, the integration constant, and the classical source terms (J ext and B rem ).

A New System of Equations
In (7), the system of a traditional Fourier-based analytical model is introduced.However, as mentioned in Section 2, that system is underdetermined if the current densities are not known.In Section 3, an equation for the current density was proposed (13); combining this equation with the equation for the flux coupled with a coil (20) gives the following: Considering ( 21), this can be written in matrix notation as follows: where [  ] is   × 1 matrix containing the terminal voltages of every coil.
[] is the   ×  ] winding matrix of the machine.This means that [] ,] =  (,]) .[ (]) ] in turn is a 1×  matrix, the th element of which is calculated as follows: The above equation for the current density is valid in each of the  ] subdomains.A matrix notation for the resulting set of equations can be found: where Rearranging gives the following: where [ where   ] is the identity matrix of size  ] .
The above implies that the combination of ( 7) and (26) is a system of equations that uniquely defines both the integration constants and the current densities in each of the subdomains.This system is written in matrix form as follows: Figure 1: Geometry and subdomains of the studied machine.
Solving this system for every time-harmonic order will uniquely define the MVP and the current density in each of the subdomains.It can readily be seen that ( 28) offers to directly impose a voltage signal to the coils instead of the classical approach of imposing a current density to the subdomains.

Example
By means of example, this section discusses the determination of the coefficient matrices for a permanent-magnet synchronous machine.The goal is to clarify the above theory.
The same PMSM topology will be used in Section 7 to validate the presented work.
6.1.Geometry.The machine that is studied in this section is a high-speed, radial-flux PMSM with a Shielding Cylinder (SC).The latter is a conductive sleeve that is wrapped around the magnets.Shielding Cylinders are often used in high-speed PMSMs to reduce the rotor eddy-current losses, thereby reducing the risk of permanent demagnetization of the magnets.The studied machine is shown in Figure 1.Its geometrical parameters include outer radius of the rotor yoke ( 1 ), outer radius of the magnets ( 2 ), outer radius of the SC ( 3 ), outer radius of the air gap ( 4 ), outer radius of the slots ( 5 ), and outer radius of the stator yoke ( 6 ).Each magnet spans an angle of   radians.The slots, which are indicated with a slot number , have a starting angle   and an opening angle .The machine has a three-phase, single-layer, and distributed winding with one slot per pole and per phase.Every phase consists of 2 parallel coils, which in turn consist of   turns.
As already mentioned, the studied geometry is divided into a number of subdomains.Those subdomains are shown in Figure 1.Note that every slot is a separate subdomain, indicated with an index 4 where  is the slot number.Evidently, the machine's 6 coils are distributed solely over the slots.The externally imposed current densities in subdomains 1, 2, and 3 are thus a priori known; they are zero.This means that the dimension of [  ] is reduced to the amount of current-carrying subdomains, that is, the 12 slots in this example.Moreover, the dimension of the electric problem ( 26) is reduced to the 12 slots as well.

Magnetic Calculations.
A Fourier-based analytical model for the machine depicted in Figure 1 was introduced by the authors in [3].Using a cylindrical coordinate system fixed to the rotor (, , ), these authors presented solutions for the MVP and a set of boundary conditions.This implies that the results of [3] allow for easily determining [ 1   ], [ 2  ], and [ 3  ].For that reason, the determination of these matrices will not be discussed in detail here.
Similarly, [] can be obtained from Figure 1 and (24).Its first row, corresponding to the rightmost slot in Figure 1, is as follows: where  slot is the surface of a slot and every column corresponds to a coil, the first one being the first coil of phase .Note that (29) and (30) regard all 12 slots and all 6 coils.Because of the machine's periodicity, it is possible to only consider 6 slots and 3 coils.This of course reduces the amount of unknowns.However, for the sake of generality, all 12 slots and all 6 coils are regarded here.From (20) and ( 21), it can be seen that matrices [ 4  ], [ 5  ], and [ 6   ] are determined by  (]) (), that is, by the expression for the MVP in each of the subdomains.In [3], an equation for the MVP in slot  was calculated as follows: where As mentioned in Section 4, there are several ways to compute  (]) (), that is, to approximate  (]) ( (])  ,  (])  , ).In this work, an average value is computed by integrating the MVP over the slot and dividing the result by the surface of the slot: where which implies that the th time-harmonic order of  (4) () can be calculated as follows: Combining the above with (20) and ( 21), it can be concluded that every row  of [ 4  ] contains only one nonzero element, that is, the element corresponding to  (4) ,0 .This element simply equals 1. [ 5   ] turns out to be a scalar matrix of size 12, whose scalar is .[ 6   ] is a zero matrix.The latter could indeed be expected as there is no remanent magnetic flux density in the slots.
As  ] can easily be calculated from (27a), (27b), (27c), and (27d).This implies that the complete electromagnetical system (28) is now defined in terms of the machine's geometrical parameters, its winding distribution, and its mechanical pulsation.

Validation
The goal of this section is to validate the work that was presented in Sections 2-6.To do so, a voltage-powered PMSM, as the one presented in Figure 1, is studied with the Fourier-based analytical model.The results are compared to results obtained from a finite-element model.
The parameters of the studied machine are shown in Table 1.
A sinusoidal voltage with an amplitude of 150 V and a frequency of 1000 Hz is applied to the machine.Because of this very high frequency, Litz wire is used for the windings.This implies that the ohmic resistance of the coils can be computed with the help of Pouillet's law.The inductance of the endwindings was calculated using the formulas presented in the book of Pyrhönen et al. [14].The highest time-and spatialharmonic orders are 50; that is, ℎ  = ℎ  = 50.
Note that the interest here is to study the accuracy of the armature-reaction field.Indeed, as saturation is neglected, applying a voltage instead of a current has an effect on the armature-reaction situation, not on the no-load situation.For that reason, the remanent magnetic flux of the magnets (B rem ) is set to 0 T.
The calculated radial and tangential components of B in the middle of the air gap are shown in Figures 2(a  respectively.The agreement between the results obtained from the analytical model and the results obtained from the finite-element model is very good.This proves that the presented theory is indeed effective.

Conclusion
This work presents a technique to directly impose the terminal voltage in Fourier-based analytical models for electric machines.The idea is to combine the magnetic calculations of classical Fourier-based analytical models with the equation for the terminal voltage of an electric machine, a technique that is well known in finite-element modeling but was not yet translated to Fourier-based analytical models.Firstly, a very general discussion on the technique and its implementation is presented.This discussion is then illustrated by means of an example.Finally, the work is validated with the help of a finite-element model.The accuracy was found to be very good.
It can be concluded that this work provides an extension to the existing Fourier-based analytical models for electric machines.It presents both a very general discussion and an example on how to directly account for a voltage source instead of having to impose current sources.As the great majority of modern drive trains are powered with a voltage source, the work offers more realistic analytical modeling of electric machines.This conclusion is of great significance for machine designers who require fast and accurate modeling tools to cope with large design spaces.
Fourier-based models are a great tool because of their accurate and fast calculations.However, they are relatively complex when compared to more traditional analytical models.Although the work presented in this paper allows for even better optimization procedures, it adds to the models complexity.For that reason, machine designers should make a well-considered choice between a more realistic model with voltage sources or a simpler model with current sources.It would therefore be interesting to compare both approaches in future work.

Figure 2 :
Figure 2: Analytical and finite-element results for the magnetic flux density in the middle of the air gap at armature-reaction conditions.

Table 1 :
Parameters of the studied machine.