Mathematical Problems in Engineering

Volume 2015, Article ID 195410, 8 pages

http://dx.doi.org/10.1155/2015/195410

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

^{1}Department of IT&C, Electrical Energy Research Group, Ghent University, V. Vaerwyckweg 1, 9000 Gent, Belgium^{2}Department of EESA, Electrical Energy Laboratory, Ghent University, Technologiepark-Zwijnaarde 913, 9052 Gent, Belgium

Received 29 June 2015; Accepted 19 September 2015

Academic Editor: Jamal Berakdar

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

#### Abstract

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 great majority of modern 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.

#### 1. 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–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 high-speed 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.

#### 2. 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 , is defined through its curl:where 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 . 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–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.

##### 2.1. 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.

##### 2.2. Subdomains

In order to simplify the governing equation, the problem is now divided in regions with constant electromagnetic properties, that is, constants and . Such regions are called subdomains; they are denoted with index . Usually, the governing equation will reduce to a Laplace, a Poisson, or a Helmholtz equation in each of the subdomains.

##### 2.3. Separation of Variables

Using separation of variables to solve (3) in subdomain , can be written as follows:withwhere and are the time- and spatial-harmonic orders, 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 .

Note that in practice the infinite summations have to be truncated. The highest time- and spatial-harmonic orders that are considered are and , respectively.

##### 2.4. 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 , these conditions are written as follows:where is the magnetic field strength in subdomain and is the current density on the boundary between subdomains and . Using the definition of the MVP and the magnetic constitutive relation, these boundary conditions can be written in terms of the MVP [3–9].

##### 2.5. 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 matrix, every row of which refers to an integration constant. Analogously, is a 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 matrix, with each row referring to the current density in a subdomain. , , and are coefficient matrices with respective size of , , and . These matrices depend on the machine’s geometry and are therefore not regarded here.

The above implies that the system contains 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.(iii)Quasistatic situation.(iv)Uniform current density in every subdomain.(v)Infinite permeability of the iron.

#### 3. Terminal Voltage and Current Density

The previous section introduced the calculation of the magnetic field of an electric machine with the help of the Fourier-based 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 can be divided in a component related to the active part of the coil and a component related to the end-windings . 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 . 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 .

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.

#### 4. 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.

##### 4.1. 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 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; and 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 .

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 in the subdomain that contains the going conductor of , in the subdomain that contains the returning conductor, and 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:

##### 4.2. 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 thatwherewhich implies that equals if subdomain contains going conductors of , equals if subdomain contains returning conductors of , and equals 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 matrix describing the winding configuration of coil ; that is, . , , and are coefficient matrices with respective sizes of , , and . The content of , , and greatly depends on the studied geometry. For that reason, it will not be regarded here. However, , , and 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 .

#### 5. 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 matrix containing the terminal voltages of every coil. is the winding matrix of the machine. This means that . in turn is a 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 is matrix whose th row equals . Rearranging gives the following:where , , , and are matrices with respective sizes of , , , and . They are calculated as follows: 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: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.

#### 6. 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 , outer radius of the magnets , outer radius of the SC , outer radius of the air gap , outer radius of the slots , and outer radius of the stator yoke . 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.