Next Article in Journal
On the Mixed Dirichlet–Steklov-Type and Steklov-Type Biharmonic Problems in Weighted Spaces
Next Article in Special Issue
Study of a Hybrid Excitation Synchronous Machine: Modeling and Experimental Validation
Previous Article in Journal
An Elastodiffusive Orthotropic Euler–Bernoulli Beam Considering Diffusion Flux Relaxation
Previous Article in Special Issue
Magnetic Field Analytical Solution for Non-homogeneous Permeability in Retaining Sleeve of a High-Speed Permanent-Magnet Machine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exact Two-Dimensional Analytical Calculations for Magnetic Field, Electromagnetic Torque, UMF, Back-EMF, and Inductance of Outer Rotor Surface Inset Permanent Magnet Machines

Department of Electrical and Electronics Engineering, Shiraz University of Technology, Shiraz 13876-71557, Iran
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2019, 24(1), 24; https://0-doi-org.brum.beds.ac.uk/10.3390/mca24010024
Submission received: 7 January 2019 / Revised: 8 February 2019 / Accepted: 13 February 2019 / Published: 17 February 2019
(This article belongs to the Special Issue Mathematical Models for the Design of Electrical Machines)

Abstract

:
This paper presents a two-dimensional analytical model of outer rotor permanent magnet machines equipped with surface inset permanent magnets. To obtain the analytical model, the whole model is divided into the sub-domains, according to the magnetic properties and geometries. Maxwell equations in each sub-domain are expressed and analytically solved. By using the boundary/interface conditions between adjacent sub-regions, integral coefficients in the general solutions are obtained. At the end, the analytically calculated results of the air-gap magnetic flux density, electromagnetic torque, unbalanced magnetic force (UMF), back-electromotive force (EMF) and inductances are verified by comparing them with those obtained from finite element method (FEM). One of the merits of this method in comparison with the numerical model is the capability of rapid calculation with the highest precision, which made it suitable for optimization problems.

1. Introduction

The existence of different types of PM brushless machines (PMBLM) made them applicable for a wide range of applications. PMBLMs have superiorities in comparison with other rivals like induction machines or reluctance synchronous machines due to higher efficiency, high torque per volume, lower torque ripple, lower vibration, and lower acoustic noise.
PMBLM can be categorized in terms of various criteria such as the topology of PMs, the relative position of the rotor and stator, the slotted or slotless stator structure, etc.
Various PM topologies such as surface-mounted, surface-inset, and interior are used where each has its own advantages and disadvantages. Among these topologies, surface-inset can provide a compromise between the other two topologies.
Electric machines with single rotor and single stator can be either inner rotor or outer rotor. The outer rotor motors can develop more output torque than the inner ones for the same volume of the machine. Usually, inner rotor machines are used for applications, which need rapid acceleration and deceleration. Outer rotor machines usually are used for applications which need constant speed. Also, the mechanical robustness of the PMs in the outer rotor configuration is higher than the inner one.
In this paper because of the aforementioned advantages of the outer rotor machines and surface inset PMs, an exact two-dimensional electromagnetic model for this type of machines is extracted.
In the design procedure, the static model is normally considered. Numerous static models have been presented for electric machines, in which some of them are based on the analytical approaches [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52], and the others are based on numerical methods [53,54,55].
The presented analytical model for electric machines are based on permeance model [1] or magnetic equivalent circuit, also known as (a.k.a.) 0-D analytic model [2], or resolution of the Maxwell’s equations in 2-D plane (a.k.a. 2-D analytic model) [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48] or 3-D plane (a.k.a. 3-D analytic model) [49,50]. Also, other methods based on mapping techniques such as Schwarz-Christoffel have been used to extract the model of some machine with the analytical approach [51,52]. The most accurate presented models among the analytical methods are 3-D analytic and 2-D analytic. The 2-D analytic method can be used instead of 3-D when the model is symmetric and has no skewing. Two-dimensional analytical models are not only capable of considering a high number of harmonics which elevate the precision of the models, but also have less computational time in comparison with the numerical methods and made them appropriate for optimal design problems.
Two-dimensional analytical models are presented for the slotted [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,19,20,29,30,31,32,33,34,35,36,37,38,41,42,44,45] and slotless [21,22,23,24,25,26,27,28] machines, equipped with surface mounted [3,4,5,6,10,11,12,14,15,16,17,18,19,20,23,26,32,34,37,41,52] and surface inset [9,13,21,22,24,25,35,36] or spoke type magnets [7,8], in order to obtain the important quantities like magnetic flux density, electromagnetic torque, unbalanced magnetic force, back-EMF, and inductances. Also, the 2-D analytical model is used to calculate the eddy current effect [11,18,22,31,32,49,50] in electrical machines. Most of the abovementioned publications are focused on the inner rotor structure [5,6,7,8,9,10,11,12,13,14,15,16,17,18,22,23,24,25,26,27,31,32,33,34,35,36,37,38,42,44,45,46] and only a few of them present the 2-D model for outer rotor machines with surface mounted PMs [4,28,29,30,41]. Therefore, to the best knowledge of the authors, it is for the first time that 2-D analytical model of brushless PM machines with outer rotor and surface inset PM is presented using the subdomain technique.
Most of the developed 2-D or 3-D analytic model are assumed with the infinite permeability of the iron parts. New techniques to account for finite soft-magnetic permeability have been recently developed, i.e., the multi-layer model using the Cauchy’s product theorem is presented in [38], and the subdomain technique by applying the superposition principle in both directions is proposed in [39,40,41,42,43,44,45,46] which can be used to calculate the core losses and saturation phenomena.
An overview of the analytical models in the Maxwell‒Fourier method with a global or local saturation effect has been realized in [40]. According to [48], Dubas’ superposition technique [39,40] is very interesting since it enables the magnetic field calculation in the material of slotted geometries. This superposition technique has been implemented in radial-flux electrical machines with(out) PMs supplied by a direct or alternate current [44,45].
The presented technique in [39,40,41,42,43,44,45,46] is not only used to predict the magnetic field in all parts of the electrical machines, but also it is used to obtain a 2-D analytical model of the steady state heat transfer of the electrical machines by solving the heat equations [47].
The aim of this paper is to extract a 2-D analytical model of PM brushless outer rotor machines equipped with surface inset PMs. The model is used to analytically compute the electromagnetic torque, torque ripple, back-electromotive force (EMF), inductances and unbalanced magnetic forces (UMF).
Therefore, this paper is organized as follows. In Section 2 the procedure of extracting the 2-D model is explained. Section 3 is dedicated to the calculation of the electromagnetic quantities. In Section 4 the analytical results of the case study are presented and compared with those of the numerical method. In the final part this paper is concluded.

2. Extracting the Magnetic Model

2.1. Assumptions

Figure 1 shows the topology of an outer rotor surface inset brushless permanent magnet machine.
In order to make the problem solvable, below assumptions are made:
(a)
According to the geometry and the absence of skewing, the problem is solved in 2-D polar coordinates which means the end effect is neglected.
(b)
Magnetic vector potential has just axial component which is function of r and θ . Consequently, magnetic flux density has radial and tangential component; i.e., A = [ 0 , 0 , A z ] , B = [ B r , B θ , 0 ] .
(c)
All materials are isotropic.
(d)
Rotor and stator back iron have infinite permeability.
(e)
The edges of the slots and slot-openings have radial direction.
(f)
The eddy current effect is neglected.

2.2. Dividing Region into Sub-Regions

According to the shape and material characteristics, the whole domain is divided into a number of sub-domains. All the sub-domains are illustrated in Figure 1 and listed in Table 1 for a PMBLM with Q slots and p pole-pairs.
When the winding is single-layer or double-layer non-overlapping, as shown respectively in Figure 2a,b, each slot is a single subdomain; however, if the winding is two-layer overlapping, as shown in Figure 2c, each slot is divided into two subdomains.

2.3. Extracting the Magnetic Model

In this part, for each sub-domain a partial differential equation is extracted based on Maxwell’s equations.
Maxwell’s equations in quasi-static form are as follows:
B = 0
× H = J + J ^
where B is the magnetic flux density vector. Ampere’s law represents the relation between the magnetic field intensity vector ( H ), external current density vector ( J ), and current density vector in media ( J ^ ). In this investigation the current density vector in media is assumed to be negligible, i.e., J ^ = 0 .
The relation between the magnetic flux density vector and the magnetic field intensity vector in permanent magnet media with linear demagnetizing curve is as follows:
B = μ 0 μ r H + μ 0 M
where M is the magnetization vector.
Substituting (3) in (2) yields
× B = μ 0 μ r J + μ 0 × M
The magnetic flux density vector can be represented as the curl of the magnetic vector potential ( A ):
B = × A
Using (4) and (5) the following expression is obtained:
2 A = μ 0 μ r J μ 0 × M
For each sub-domain, Equation (6) results in Poisson equations for the magnet and slot regions, and Laplace equations for the air-gap and slot-opening sub-domains, as represented below.
2 A s l , j = μ 0 μ r J
2 A m , k = μ 0 × M
2 A χ = 0 ,   χ = { ( a ) , ( s o , j ) }
In 2-D polar coordinates, the magnetic vector potential and the current density vector have just a component along z , i.e., A = [ 0 , 0 , A z ( r , θ ) ] and J = [ 0 , 0 , J z ( θ , t ) ] . Also, the magnetic flux density vector and the magnetization vector have radial and tangential components as below:
B = [ B r ( r , θ ) , B θ ( r , θ ) , 0 ]
M = [ M r ( r , θ ) , M θ ( r , θ ) , 0 ]
Therefore, Equations (7)–(9) are rewritten as
1 r r ( r A z s l , j r ) + 1 r 2 2 A z s l , j θ 2 = μ 0 J z s l , j
1 r r ( r A z m , k r ) + 1 r 2 2 A z m , k θ 2 = μ 0 r ( M θ k M r k θ )
1 r r ( r A z χ r ) + 1 r 2 2 A z χ θ 2 = 0 ,   χ = { ( a ) , ( s o , j ) }

2.4. Boundary Conditions

The perpendicular magnetic flux density in two adjacent sub-domains must be equal as mathematically represented as follows:
n . ( B i B i + ) = 0
In this equation, B i is the magnetic flux density in sub-domain i , and B i + is the magnetic flux density in the sub-domain i + .
Also, if there is no current between the two adjacent sub-domains, the tangential components of the magnetic field intensity at the boundary of the two sub-domains are equal; this expression is shown mathematically by Relation (14).
n × ( H i H i + ) = 0
In this equation H i is the magnetic field intensity of the sub-domain i and H i + is the magnetic field intensity of sub-domain i + .
In both (13) and (14), n is the perpendicular unit vector to the interface between two adjacent sub-domains.
According to Figure 1, all boundary/interface conditions between sub-domains have been shown from (15) to (25) where α r , β , and δ are respectively the magnet arc per pole pitch ratio, the span angle of slot-openings, and the span angle of slots as shown in Figure 1.
Domain (i)Domain (i+)EquationBorder of the InterfaceLimit
MagnetRotor yoke H θ m , k ( r , θ ) = 0 r = R r | θ α k π / p | α r π / 2 p (15)
MagnetIron next to the PM pole H r m , k ( r , θ ) = 0 θ = α + k π / p ±   α r π / 2 p R r r R m (16)
Air-gapMagnet B r a ( r , θ ) = B r m , k ( r , θ ) r = R m | θ α k π / p | α r π / 2 p (17)
Air-gapMagnet H θ a ( r , θ ) = { k = 0 p 1 H θ m , k ( r , θ ) 0 r = R m { | θ α k π / p | α r π / 2 p elsewhere (18)
Slot-openingEdges of slot-opening H r s o , j ( r , θ ) = 0 θ = θ j ± β / 2 R s r R s o (19)
Air-gapSlot-opening B r a ( r , θ ) = B r s o , j ( r , θ ) r = R s | θ θ j | β 2 (20)
Air-gapSlot-opening H θ a ( r , θ ) = { j = 1 Q H θ s o , j ( r , θ ) 0 r = R s { | θ θ j | β 2 elsewhere (21)
SlotSlot-opening B r s l , j ( r , θ ) = B r s o , j ( r , θ ) r = R s o | θ θ j | β 2 (22)
SlotEdge of tooth
Slot-opening
Edge of tooth
H θ s l , j ( r , θ ) = { 0 H θ s o , j ( r , θ ) 0 r = R s o { θ j δ 2 θ < θ j β 2 θ j β 2 θ θ j + β 2 θ j + β 2 < θ θ j + δ 2 (23)
SlotEdges of the slot H r s l , j ( r , θ ) = 0 θ = θ j ± δ / 2 R s o r R s l (24)
SlotStator yoke H θ s l , j ( r , θ ) = 0 r = R s l | θ θ j | δ 2 (25)

2.5. Extracting the Fourier Series of the Armature Reaction

To consider the effect of the armature reaction, the current density of each slot should be represented as Fourier series. The current of each phase varies with time and can be represented as Equation (26).
I k ( t ) = v I v sin [ v ( p ω t γ k ) + θ v ] ,   k = 1 , 2 , , q
where q is the number of the phases, v shows the order of the harmonics, ω is the angular velocity of the rotor, p is the number of the pole-pairs, γ k = 2 π ( k 1 ) / q is the time offset of the phase k t h respect to the first phase. Also I v and θ v are the magnitude and phase offset of v t h harmonic, respectively.
It is obvious that the relation between the current density in each slot and phase current is dependent on the winding configuration. If the winding configuration is like Figure 2a,b, each slot is considered as one sub-region, like Figure 3a,b. But, if the configuration is similar to Figure 2c, each slot consists of two sub-regions (upper and lower sub-regions), as represented in Figure 3c.
The current density in each sub-region of a slot can be represented in Fourier series form as in Equation (27).
J z j ( θ , t ) = J 0 j ( t ) + v = 1 J v j ( t ) cos ( π v δ ( θ θ j + δ / 2 ) ) θ j δ / 2 θ θ j + δ / 2 .
where J 0 j ( t ) and J v j ( t ) are as follows:
J 0 j ( t ) = J j ( t ) + J r j ( t ) 2
J v j ( t ) = J j ( t ) J r j ( t ) π v / 2 sin ( π v / 2 )
In order to complete the Fourier series of each sub-region of a slot, it is necessary to obtain the current density of phases in each sub-region of a slot at a specific time by Equations (30) and (31):
J = I K f A s l o t
J = I K f A s l o t / 2
For instance, the current density in each sub-domain of a slot could be as Figure 4a,b.
If the winding configuration is as shown in Figure 2a,c, the figure of the current density in each sub-region of slot, at a time instant will be as Figure 4a, and if the winding configuration is as shown in Figure 2b, the figure of the current density in a time instant will be as Figure 4b. Also if the configuration of the winding is similar to Figure 2c, the current density in each sub-region of the slot will be as Figure 4a.

2.6. Extracting the Fourier Series of the Magnetization

In the 2-D polar coordinate system, the magnetization vector only has the radial and tangential components as Equation (32).
M = M r r + M θ θ
where r and θ are the radial and tangential unit vectors. M r and M θ are the components of magnetization vector which can be represented as Fourier series expansion of Equations (33) and (34).
M r ( θ r ) = n = 1 , 3 , 5 , M r n cos ( n p θ r )
M θ ( θ r ) = n = 1 , 3 , 5 , M θ n sin ( n p θ r )
where M r n and M θ n respectively are the radial and tangential coefficient of the Fourier series and will be determined according to the magnetization pattern (Table 2). In this paper, only the radial magnetization pattern has been used and represented in Figure 5.

2.7. Finding the General Solution

The overall format for the general solution in all sub-domains can be represented as Equation (35):
A ( r , θ ) = n = 1 ( A n r n + B n r n ) . ( C n cos n θ + D n sin n θ ) + ( A 0 lnr + B 0 ) ( C 0 θ + D 0 )
The general solution not only has the capability to satisfy the related PDE, but must satisfy the boundary conditions of the related sub-domain, especially Equations (16), (19), and (24). So the general solutions for sub-domains are as Equations (36)–(40).
The general solution of Poisson equation in slot sub-domain will be as follows:
A z s l , j ( r , θ ) = b 0 s l , j ln r + A p s l , j ( r , θ ) + v = 1 [ a v s l , j ( r R s o ) π v δ + b v s l , j ( R s l r ) π v δ ] × cos ( π v δ ( θ θ j + δ 2 ) )
The particular solution is as follows:
A p s l , j ( r , θ ) = μ 0 4 J 0 j r 2 + v = 1 μ 0 J v j r 2 ( π v δ ) 2 4 cos ( π v δ ( θ θ j + δ 2 ) )
Also the general solution for the slot-opening sub-domain is
A z s o , j ( r , θ ) = b 0 s o , j ln r + u = 1 [ a u s o , j ( r R s ) π u β + b u s o , j ( R s o r ) π u β ] × cos ( π u β ( θ θ j + β 2 ) )
The general solution for the air-gap sub-domain is
A z s o , j ( r , θ ) = b 0 s o , j ln r + u = 1 [ a u s o , j ( r R s ) π u β + b u s o , j ( R s o r ) π u β ] × cos ( π u β ( θ θ j + β 2 ) )
And finally, Poisson equation in the PMs sub-domain has the general solution as Equation (40).
A z m , k ( r , θ ) = b 0 m , k ln r + w = 1 [ a w m , k ( r R r ) w p α r + b w m , k ( R m r ) w p α r + k w k r ] × cos ( w p α r ( θ α k π p + α r π 2 p ) )
where
k w k = μ 0 χ w { w p α r M r w k M θ w k ( w p α r ) 2 1 w p α r M r w k M θ w k 2 ln r w p = α r
χ w = 1 ( 1 ) w 2
In order to simplify the general solution in PM sub-domain, boundary condition (15) has been implemented.
A z m , k ( r , θ ) = w = 1 { b w m , k [ ( R m R r ) w p α r ( r R r ) w p α r + ( R m r ) w p α r ] + R r ζ w 1 k ( r R r ) w p α r + k w k r } × cos ( w p α r ( θ α k π p + α r π 2 p ) )
where
ζ w 1 k = α r p w ( d k w k r d r | r = R r + μ 0 χ w M θ w k ) = μ 0 χ w { M r w k w p α r M θ w k ( w p α r ) 2 1 w p α r M r w k M θ w k 2 ( 1 + ln R r ) M θ w k w p = α r .
Also, by implementing boundary condition (25), the general solution in slot sub-domain will be simplified as in Equation (45).
A z s l , j ( r , θ ) = μ 0 4 J 0 s l , j ( 2 R s l 2 ln ( r ) r 2 ) + v = 1 V { b v s l , j [ ( r R s o ) π v δ + ( R s l r ) π v δ ( R s l R s o ) π v δ ] + μ 0 J v j ( π v δ ) 2 4 [ r 2 2 R s l π v δ ( R s l r ) π v δ ] } × cos ( π v δ ( θ θ j + δ 2 ) )

2.8. Obtaining Integral Coefficients

For implementing boundary condition (17), the correlation technique must be used [31].
By multiplying Equation (17) in 2 p α r π sin ( w p α r ( θ α k π p + α r π 2 p ) ) and integration over [ α + k π p π α r 2 p , α + k π p + π α r 2 p ] , Equation (46) will be obtained.
2 p α r π α + k π / p α r π / 2 p α + k π / p + α r π / 2 p B r a | r = R m sin ( w p α r ( θ α k π p + α r π 2 p ) ) d θ = 2 p α r π α + k π / p α r π / 2 p α + k π / p + α r π / 2 p B r m , k | r = R m sin ( w p α r ( θ α k π p + α r π 2 p ) ) d θ
From Equation (46), Equation (47) will be deduced.
w p α r [ 1 + ( R m R r ) 2 w p α r ] a w m , k n = 1 N n { [ a n a ( R s R m ) n + b n a ] σ s ( n , w , k ) [ c n a ( R s R m ) n + d n a ] σ c ( n , w , k ) } = R m w p α r [ ζ w 1 k ( R m R r ) w p α r + 1 + ζ w 2 k ]
where
ζ w 2 k = k w k | r = R m = μ 0 χ w { w p α r M r w k M θ w k ( w p α r ) 2 1 w p α r M r w k M θ w k 2 ln R m w p = α r
σ c ( n , w , k ) = 2 p α r π α + k π / p α r π / 2 p α + k π / p + α r π / 2 p cos ( n θ ) sin ( w p α r ( θ α k π p + α r π 2 p ) ) d θ
σ s ( n , w , k ) = 2 p α r π α + k π / p α r π / 2 p α + k π / p + α r π / 2 p sin ( n θ ) sin ( w p α r ( θ α k π p + α r π 2 p ) ) d θ
The solutions of the integrals have been given in the Appendix A.
By using correlation technique, Relation (18) must be multiplied by 1 π sin ( n θ ) and integration on the interval [ α π , α + π ] , Equation (51) will be obtained.
1 π α π α + π H θ a | r = R m sin ( n θ ) d θ = 1 π k = 0 2 p 1 α + k π / p α r π / 2 p α + k π / p + α r π / 2 p H θ m , k | r = R m sin ( n θ ) d θ
Equation (51) results in Equation (52):
k = 0 2 p 1 w = 1 W w p μ r α r [ ( R m R r ) 2 w p α r 1 ] ρ s ( n , w , k ) a w m , k + n [ c n a ( R s R m ) n d n a ]     = k = 0 2 p 1 w = 1 W w p R m μ r α r [ ζ w 1 k ( R m R r ) w p α r 1 + ζ w 3 k ] ρ s ( n , w , k )
where
ζ w 3 k = α r w p ( d k w k r d r | r = R m + μ 0 χ w M θ w k ) = μ 0 χ w { M r w k w p α r M θ w k ( w p α r ) 2 1 w p α r M r w k M θ w k 2 ( 1 + ln R m ) M θ w k w p = α r
ρ s ( n , w , k ) = 1 π α + k π / p α r π / 2 p α + k π / p + α r π / 2 p cos ( w p α r ( θ α k π p + α r π 2 p ) ) sin ( n θ )   d θ
The solution of the integral has been given in the Appendix A.
Again Equation (18) must be multiplied by 1 π cos ( n θ ) , then integration on the interval [ α π , α + π ] causes:
1 π α π α + π H θ a | r = R m cos ( n θ ) d θ = 1 π k = 0 2 p 1 α + k π / p α r π / 2 p α + k π / p + α r π / 2 p H θ m , k | r = R m cos ( n θ ) d θ
Equation (55) results in Equation (56).
k = 0 2 p 1 w = 1 W w p μ r α r [ ( R m R r ) 2 w p α r 1 ] ρ c ( n , w , k ) a w m , k + n [ a n a ( R s R m ) n b n a ]    = k = 0 2 p 1 w = 1 W w p R m μ r α r [ ζ w 1 k ( R m R r ) w p α r 1 + ζ w 3 k ] ρ c ( n , w , k )
where
ρ c ( n , w , k ) = 1 π α + k π / p α r π / 2 p α + k π / p + α r π / 2 p cos ( w p α r ( θ α k π p + α r π 2 p ) )   cos ( n θ )   d θ
The solution of the integral has been given in the Appendix A.
For implementing boundary condition (20), multiplying 2 β sin ( π u β ( θ θ j + β 2 ) ) to Equation (20) and integration over [ θ j β 2 , θ j + β 2 ] yields
2 β θ j β / 2 θ j + β / 2 B r s o , j ( r , θ ) | r = R s sin ( π u β ( θ θ j + β 2 ) ) d θ = 2 β θ j β / 2 θ j + β / 2 B r a ( r , θ ) | r = R s sin ( π u β ( θ θ j + β 2 ) ) d θ
Equation (59) is obtained via the simplification of Equation (58).
n = 1 N n { [ a n a + b n a ( R s R m ) n ] ε s ( n , u , j ) [ c n a + d n a ( R s R m ) n ] ε c ( n , u , j ) } + π u β ( R s o R s ) π u β a u s o , j + π u β b u s o , j = 0
where
ε s ( n , u , j ) = 2 β θ j β / 2 θ j + β / 2 sin ( n θ ) sin ( π u β ( θ θ j + β 2 ) ) d θ
ε c ( n , u , j ) = 2 β θ j β / 2 θ j + β / 2 cos ( n θ ) sin ( π u β ( θ θ j + β 2 ) ) d θ
The solution of the integral has been given in the Appendix A.
The correlation technique is used for boundary condition (21) and 1 π cos ( n θ ) is multiplied to Equation (21) and integration is taken over [ π , π ] .
1 π π π H θ a ( r , θ ) | r = R s cos ( n θ ) d θ = 1 π j = 1 Q θ j β / 2 θ j + β / 2 H θ s o , j ( r , θ ) | r = R s cos ( n θ ) d θ
Simplifying Equation (62) yields (63).
n [ a n a b n a ( R s R m ) n ] j = 1 Q u = 1 U π u β [ ( R s o R s ) π u β a u s o , j b u s o , j ] η c ( n , u , j ) j = 1 Q η c ( n , 0 , j ) b 0 s o , j = 0
where
η c ( n , u , j ) = 1 π θ j β / 2 θ j + β / 2 cos ( n θ ) cos ( π u β ( θ θ j + β 2 ) ) d θ
Also, by multiplying Equation (21) to 1 π sin ( n θ ) and integration over interval [ π , π ] the following expression is obtained:
1 π π π H θ a ( r , θ ) | r = R s sin ( n θ ) d θ = 1 π j = 1 Q θ j β / 2 θ j + β / 2 H θ s o , j ( r , θ ) | r = R s sin ( n θ ) d θ
Simplification of Equation (65) causes the formation of Equation (66).
n [ c n a ( R s R m ) n d n a ] j = 1 Q u = 1 U π u β [ ( R s o R s ) π u β a u s o , j b u s o , j ] η s ( n , u , j ) j = 1 Q η s ( n , 0 , j ) b 0 s o , j = 0
where
η s ( n , u , j ) = 1 π θ j β / 2 θ j + β / 2 sin ( n θ ) cos ( π u β ( θ θ j + β 2 ) ) d θ
The correlation technique is used and boundary condition (22) is multiplied by 2 β sin ( π u β ( θ θ j + β 2 ) ) . Integration over interval [ θ j β 2 , θ j + β 2 ] results in Equation (68).
2 β θ j β / 2 θ j + β / 2 B r s l , j ( r , θ ) | r = R s o sin ( π u β ( θ θ j + β 2 ) ) d θ = 2 β θ j β / 2 θ j + β / 2 B r s o , j ( r , θ ) | r = R s o sin ( π u β ( θ θ j + β 2 ) ) d θ
Simplifying (68) results in (69).
π u β [ a u s o , j + ( R s o R s ) π u β b u s o , j ] v = 1 V π v δ [ ( R s l R s o ) 2 π v δ + 1 ] γ s ( u , v ) b v s l , j = v = 1 V μ 0 J v j ( π v δ ) 2 4 [ π v δ R s l 2 2 R s o 2 ( R s l R s o ) π v δ ] γ s ( u , v )
where
γ s ( u , v ) = 2 β θ j β / 2 θ j + β / 2 sin ( π v δ ( θ θ j + δ 2 ) ) sin ( π u β ( θ θ j + β 2 ) ) d θ
Solution of the integral are given in the Appendix A.
Correlation technique is implemented for boundary condition (23) and it is multiplied by 2 δ cos ( π v δ ( θ θ j + δ 2 ) ) . Integration over interval [ θ j δ 2 , θ j + δ 2 ] yields Equation (71).
2 δ θ j δ / 2 θ j + δ / 2 H θ s l , j ( r , θ ) | r = R s o cos ( π v δ ( θ θ j + δ 2 ) ) d θ = 2 δ θ j β / 2 θ j + β / 2 H θ s o , j ( r , θ ) | r = R s o cos ( π v δ ( θ θ j + δ 2 ) ) d θ
Simplifying Equation (71) results in Equation (72).
u = 1 U π u β [ a u s o , j + ( R s o R s ) π u β b u s o , j ] γ c ( u , v ) γ c ( 0 , v ) b 0 s o , j + π v δ [ ( R s l R s o ) 2 π v δ 1 ] b v s l , j = 2 μ 0 J v j ( π v δ ) 2 4 [ R s l 2 R s o 2 ( R s l R s o ) π v δ ]
where
γ c ( u , v ) = 2 δ θ j β / 2 θ j + β / 2 cos ( π v δ ( θ θ j + δ 2 ) ) cos ( π u β ( θ θ j + β 2 ) ) d θ
Integration over interval [ θ j δ 2 , θ j + δ 2 ] on boundary condition (23) causes Equation (74).
θ j δ / 2 θ j + δ / 2 H θ s l , j ( r , θ ) | r = R s o d θ = θ j β / 2 θ j + β / 2 H θ s o , j ( r , θ ) | r = R s o d θ
From Equation (74), Equation (75) will be deduced.
b 0 s o , j = μ 0 J 0 j 2 ( R s o 2 R s l 2 ) δ β

2.9. Overlapping Winding

In overlapping winding, each slot is divided into two sub-domains as represented in Figure 3c. The upper and lower sub-domains respectively are indicated with s l b and s l t indices. Hence Equation (69) will be as follows:
π u β [ a u s o , j + ( R s o R s ) π u β b u s o , j ] v = 1 V π v δ [ ( R s l / R s o ) 2 π v δ + 1 ( R s l m / R s o ) π v δ ] γ s ( u , v ) b v s l b , j = 0
Also, Equations (72) and (75) will be modified as (77) and (78) respectively.
u = 1 U π u β [ a u s o , j + ( R s o R s ) π u β b u s o , j ] γ c ( u , v ) + π v δ [ ( R s l / R s o ) 2 π v δ 1 ( R s l m / R s o ) π v δ ] b v s l b , j γ c ( 0 , v ) b 0 s o , j = 0
b 0 s o , j = μ 0 2 [ J b 0 j ( R s l m 2 R s l 2 ) + J t 0 j ( R s o 2 R s l m 2 ) ] δ β
where R s l m = ( R s l 2 + R s o 2 ) / 2 is the radii of middle of the slot which divides it into two equal areas. J b 0 j and J t 0 j are current densities in the lower and upper sub-domains in a slot.

3. Quantities

3.1. Flux Density

The air-gap flux density vector is one of the most important quantities required for the calculation of other quantities. For obtaining the air-gap flux density, Equation (5) is expanded and Relations (78) and (79) in 2-D polar coordinates are deduced.
B r a ( r , θ ) = 1 r A z a θ = n = 1 N n { [ a n a R m ( r R m ) n 1 + b n a R s ( R s r ) n + 1 ] sin ( n θ ) [ c n a R m ( r R m ) n 1 + d n a R s ( R s r ) n + 1 ] cos ( n θ ) }
B θ a ( r , θ ) = A z a r = n = 1 N n { [ [ a n a R m ( r R m ) n 1 b n a R s ( R s r ) n + 1 ] ] cos ( n θ ) + [ c n a R m ( r R m ) n 1 d n a R s ( R s r ) n + 1 ] sin ( n θ ) }

3.2. Inductances

For calculation of the inductances, just the flux produced by armature is considered. Inductances between phase k and k are obtained as follows:
L k , k = j k   & j k λ j , j i j
where j = 1 , 2 , , Q and j = 1 , 2 , , Q are the indices of the coils. i j is the current of phase k and λ j , j is the flux linked by coil j , which is produced by coil j . If k = k , self-inductance of the phase is calculated.

3.3. Back-EMF

In order to calculate the no-load back-EMF of a phase, the permanent magnet flux linked by coils must be calculated by Equation (82).
φ = B d S
According to Faraday’s law, induced voltage in jth coil obtain by Equation (83).
E j = N t ω d φ j d α
where N t is the number of turns of the coil, ω and α are the angular velocity and rotor position respectively. The total back-EMF of a phase depends on coils connections.

3.4. Instantaneous Electromagnetic Torque

Instantaneous torque consists of cogging torque ( T c o g ), electromagnetic torque ( T e m ) and reluctance torque ( T r e l ).
T ( t ) = T c o g ( t ) + T e m ( t ) + T r e l ( t )
By using Maxwell stress tensor, the instantaneous electromagnetic torque can be obtained as follows:
T ( t ) =   1 μ 0 B r B θ d s
By expanding Equation (85), Relations (86) and (87) will be obtained.
T ( t ) = L s π π 1 μ 0 ( B r , P M a + B r , A R a ) ( B θ , P M a + B θ , A R a ) | r = R c R c 2 d θ
T ( t ) = L s R c 2 μ 0 π π ( B r , P M a B θ , P M a + B r , A R a B θ , P M a + B r , P M a B θ , A R a + B r , A R a B θ , A R a ) | r = R c d θ
where B r , P M a and B θ , P M a respectively are the radial and tangential flux density components in the air-gap, due to the PMs. Also B r , A R a and B θ , A R a are the radial and tangential magnetic flux density components due to the armature reaction in the air-gap. The parameter R c is the radius of the middle of inner air-gap.

3.5. Unbalanced Magnetic Force

The radial and tangential components of the local traction exerted on the rotor surface can be obtained by Maxwell stress tensor as follows:
f r = 1 2 μ 0 ( B r 2 B θ 2 )
f θ = 1 μ 0 B r B θ
By transforming these local tractions to the Cartesian plane, and summation of the same directions, Equations (90)–(93) will be obtained.
f x = f r cos θ f θ sin θ
f y = f r sin θ + f θ cos θ
F x ( t ) = L / 2 L / 2 π π f x r d θ d z = L π π f x r d θ
F y ( t ) = L / 2 L / 2 π π f y r d θ d z = L π π f y r d θ
Finally, the amplitude of the unbalanced magnetic force can be obtained by Equation (94).
F r = | F ( t ) | = F x 2 ( t ) + F y 2 ( t )

4. Results

In order to investigate the efficacy of the model, a case study with the parameters listed in Table 3 has been used.
Also, the winding configuration for this case study has been shown in Figure 5.

4.1. Flux Density

The radial and tangential magnetic flux density components due to the open-circuit and armature reaction are respectively depicted in Figure 6, Figure 7, Figure 8 and Figure 9. The numerical results obtained from FEM are shown respectively in each figure and confirm the accuracy of the proposed model.

4.2. Torque

Instantaneous electromagnetic torque, reluctance torque and cogging torque of the machine have been depicted in Figure 10, Figure 11 and Figure 12. Both analytic and numeric methods show good agreement which confirms the efficacy of the proposed model.

4.3. Back-EMF and Inductance

Results of the phase back-EMF and line back-EMF are shown in Figure 13. Again it is shown that both analytic and numeric results have good conformity.
Also, the self and mutual inductances have been depicted in Figure 14.

4.4. Unbalanced Magnetic Force (UMF)

Unbalanced magnetic forces due to the open-circuit, armature reaction, and both of them have been depicted in Figure 15, Figure 16 and Figure 17. As evident from these figures, unbalanced forces due to the armature reaction exert considerable forces compared to those of the open circuit.
Both analytic and numeric results are shown in Figure 15, Figure 16 and Figure 17 and confirm the correctness of the proposed model.

5. Conclusions

A 2-D analytical magnetic model is presented for brushless synchronous outer rotor machines with surface inset PMs. For this purpose, Maxwell’s equations in the form of the Laplace and Poisson equations are solved in predefined sub-domains of the 2-D polar coordinates. The general and particular solutions for each sub-region are presented so that they have the capability to satisfy the governing PDE and related boundary conditions. Finally, by imposing boundary conditions and solving simultaneous linear algebraic equations, all important quantities such as magnetic flux density, electromagnetic torque, UMF, back-EMF, and inductances are calculated and validated by those obtained by FEM. The results of the analytical model show the efficacy of the proposed approach.

Author Contributions

A.V. has done formal analysis, investigation, methodology, writing and editing the original draft. A.R. has done supervision, conceptualization, project administration and reviewing the original draft. H.M.-J. has done formal analysis, investigation, methodology, software, and validation. A.G. has done review and editing.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

To define the left and right side current density in each slot for the two layers non-overlapping concentrated winding topology, we have:
J r = i ( t ) K f A w / 2 C r
J l = i ( t ) K f A w / 2 C l
where J r = [ J r 1   J r j J r Q ] , J l = [ J l 1   J l j J l Q ] , i ( t ) = [ i a i b i c ] , K f is the filling factor and A w is the slot area. On the other hand, C r and C l are as follows:
C r ( i , j ) = { C ( i , j ) + | C ( i , j ) | 2 i f 1 < C ( i , j ) < 1 1 i f C ( i , j ) = 2 1 i f C ( i , j ) = 2
C l ( i , j ) = { C ( i , j ) | C ( i , j ) | 2 i f 1 < C ( i , j ) < 1 1 i f C ( i , j ) = 2 1 i f C ( i , j ) = 2
According to the winding topology, C is as follows:
C = [ 2 2 1 0 0 0 0 0 1 0 0 1 2 2 1 0 0 0 0 0 0 0 0 1 2 2 1 ] 3 × Q = 3 × 9
where C ( i , j ) = 2 or −2 means slot j accommodates two sides of two coils of phase i, which respectively carry positive or negative current. C ( i , j ) = 1 or −1 means slot j accommodates one side of one coil of phase i, which respectively carry positive or negative current. Also, 0 means there is no coil of the phase i th in the slot.
For α r n w p :
ρ s ( n , w , k ) = α r 2 π { cos ( w π + n π α r 2 p + n α + k n π p ) + cos ( n π α r 2 p n α k n π p ) α r n + w p     cos ( w π n π α r 2 p n α k n π p ) cos ( n π α r 2 p n α k n π p ) α r n w p }
ρ c ( n , w , k ) = α r 2 π { sin ( w π + n π α r 2 p + n α + k n π p ) + sin ( n π α r 2 p n α k n π p ) α r n + w p     sin ( w π n π α r 2 p n α k n π p ) sin ( n π α r 2 p n α k n π p ) α r n w p }
σ s ( n , w , k ) = p π   { sin ( w π + n π α r 2 p + n α + k n π p ) sin ( n π α r 2 p n α k n π p ) α r n + w p sin ( w π n π α r 2 p n α k n π p ) sin ( n π α r 2 p n α k n π p ) α r n w p }
σ c ( n , w , k ) = p π { cos ( w π + n π α r 2 p + n α + k n π p ) + cos ( n π α r 2 p n α k n π p ) α r n + w p + cos ( w π n π α r 2 p n α k n π p ) cos ( n π α r 2 p n α k n π p ) α r n w p }
For α r n = w p :
ρ s ( n , w , k ) = 1 4 n π [ cos ( 3 w π / 2 + n α + k n π p ) cos ( w π / 2 n α k n π p ) ] α r 2 p sin ( w π / 2 n α k n π p )
ρ c ( n , w , k ) = 1 4 n π [ sin ( 3 w π / 2 + n α + k n π p ) + sin ( w π / 2 n α k n π p ) ] + α r 2 p cos ( w π / 2 n α k n π p )
σ s ( n , w , k ) = 1 2 w π [ sin ( 3 w π / 2 + n α + k n π p ) + sin ( w π / 2 n α k n π p ) ] + cos ( w π / 2 n α k n π p )
σ c ( n , w , k ) = 1 2 w π [ cos ( 3 w π / 2 + n α + k n π p ) cos ( w π / 2 n α k n π p ) ] + sin ( w π / 2 n α k n π p )
For π u β n :
ε s ( n , u , j ) = 2 π u ( 1 ) u + 1 sin ( n ( θ j + β 2 ) ) + sin ( n ( θ j β 2 ) ) π 2 u 2 β 2 n 2
ε c ( n , u , j ) = 2 π u ( 1 ) u + 1 cos ( n ( θ j + β 2 ) ) + cos ( n ( θ j β 2 ) ) π 2 u 2 β 2 n 2
η s ( n , u , j ) = β 2 n π ( 1 ) u cos ( n ( θ j + β 2 ) ) cos ( n ( θ j β 2 ) ) π 2 u 2 β 2 n 2
η c ( n , u , j ) = β 2 n π ( 1 ) u + 1 sin ( n ( θ j + β 2 ) ) + sin ( n ( θ j β 2 ) ) π 2 u 2 β 2 n 2
For π u = β n :
ε s ( n , u , j ) = cos ( n ( θ j β 2 ) ) sin ( n ( θ j + 3 β 2 ) ) sin ( n ( θ j β 2 ) ) 2 n β
ε c ( n , u , j ) = sin ( n ( θ j β 2 ) ) cos ( n ( θ j + 3 β 2 ) ) cos ( n ( θ j β 2 ) ) 2 n β
η s ( n , u , j ) = sin ( n ( θ j β 2 ) ) 2 π / β cos ( n ( θ j + 3 β 2 ) ) cos ( n ( θ j β 2 ) ) 4 n π
η c ( n , u , j ) = cos ( n ( θ j β 2 ) ) 2 π / β   + sin ( n ( θ j + 3 β 2 ) ) sin ( n ( θ j β 2 ) ) 4 n π
For δ u β v :
γ s ( u , v ) = 2 δ 2 u π ( 1 ) u + 1 sin ( π v 2 δ ( δ + β ) ) + sin ( π v 2 δ ( δ β ) ) δ 2 u 2 β 2 v 2
γ c ( u , v ) = 2 β 2 v π ( 1 ) u + 1 sin ( π v 2 δ ( δ + β ) ) + sin ( π v 2 δ ( δ β ) ) δ 2 u 2 β 2 v 2
For δ u = β v :
γ s ( u , v ) = 2 π u cos ( π 2 ( u v ) ) sin ( π 2 ( 3 u + v ) ) sin ( π 2 ( u v ) ) 2 π u
γ c ( u , v ) = 2 π u cos ( π 2 ( u v ) ) + sin ( π 2 ( 3 u + v ) ) + sin ( π 2 ( u v ) ) 2 π v

References

  1. Zarko, D.; Ban, D.; Lipo, T.A. Analytical calculation of magnetic field distribution in the slotted air gap of a surface permanent-magnet motor using complex relative air-gap permeance. IEEE Trans. Magn. 2006, 42, 1828–1837. [Google Scholar] [CrossRef]
  2. Hur, J.; Yoon, S.; Hwang, D.; Hyun, D. Analysis of PMLSM using three dimensional equivalent magnetic circuit network method. IEEE Trans. Magn. 1997, 33, 4143–4145. [Google Scholar]
  3. Rahideh, A.; Vahaj, A.A.; Mardaneh, M.; Lubin, T. Two-Dimensional Analytical Investigation of the Parameters and the Effects of Magnetisation Patterns on the Performance of Coaxial Magnetic Gears. IET Electr. Syst. Transp. 2017, 7, 230–245. [Google Scholar] [CrossRef]
  4. Moayed-Jahromi, H.; Rahideh, A.; Mardaneh, M. 2-D Analytical Model for External Rotor Brushless PM Machines. IEEE Trans. Energy Convers. 2016, 31, 1100–1109. [Google Scholar] [CrossRef]
  5. Wu, L.J.; Zhu, Z.Q.; Staton, D.; Popescu, M.; Hawkins, D. An improved subdomain model for predicting magnetic field of surface-mounted permanent-magnet machines accounting for toothtips. IEEE Trans. Magn. 2011, 47, 1693–1704. [Google Scholar] [CrossRef]
  6. Rahideh, A.; Korakianitis, T. Subdomain analytical magnetic field prediction of slotted brushless machines with surface mounted magnets. Int. Rev. Electr. Eng. 2012, 7, 3891–3909. [Google Scholar]
  7. Pourahmadi-Nakhli, M.; Rahideh, A.; Mardaneh, M. Analytical 2-D model of slotted brushless machines with cubic spoke-type permanent magnets. IEEE Trans. Energy Convers. 2018, 33, 373–382. [Google Scholar] [CrossRef]
  8. Boughrara, K.; Ibtiouen, R.; Lubin, T. Analytical Prediction of Magnetic Field in Parallel Double Excitation and Spoke-Type Permanent-Magnet Machines Accounting for Tooth-Tips and Shape of Polar Pieces. IEEE Trans. Magn. 2012, 48, 2121–2137. [Google Scholar] [CrossRef]
  9. Rahideh, A.; Korakianitis, T. Analytical magnetic field calculation of slotted brushless PM machines with surface inset magnets. IEEE Trans. Magn. 2012, 48, 2633–2649. [Google Scholar] [CrossRef]
  10. Dubas, F.; Sari, A.; Kauffmann, J.M.; Espanet, C. Cogging torque evaluation through a magnetic field analytical computation in permanent magnet motors. In Proceedings of the 2009 International Conference on Electrical Machines and Systems, Tokyo, Japan, 15–18 November 2009; pp. 1–5. [Google Scholar]
  11. Dubas, F.; Espanet, C. Semi-analytical Solution of 2-D rotor eddy-current losses due to the slotting effect in SMPMM. In Proceedings of the 17th Conference on the Computation of Electromagnetic Fields COMPUMAG 2009, Florianopolis, Brasil, 22–26 November 2009; pp. 20–25. [Google Scholar]
  12. Lubin, T.; Mezani, S.; Rezzoug, A. 2-D exact analytical model for surface-mounted permanent magnet motors with semi-closed slots. IEEE Trans. Magn. 2011, 47, 479–492. [Google Scholar] [CrossRef]
  13. Lubin, T.; Mezani, S.; Rezzoug, A. Two-dimensional analytical calculation of magnetic field and electromagnetic torque for surface-inset permanent magnet motors. IEEE Trans. Magn. 2012, 48, 2080–2091. [Google Scholar] [CrossRef]
  14. Dubas, F.; Espanet, C. Analytical solution of the magnetic field in permanent-magnet motors taking into account slotting effect: No-load vector potential and flux density calculation. IEEE Trans. Magn. 2009, 45, 2097–2109. [Google Scholar] [CrossRef]
  15. Zhu, Z.Q.; Wu, L.J.; Xia, Z.P. An accurate subdomain model for magnetic field computation in slotted surface mounted permanent-magnet machines. IEEE Trans. Magn. 2010, 46, 1100–1115. [Google Scholar] [CrossRef]
  16. Boughrara, K.; Chikouche, B.L.; Ibtiouen, R.; Zarko, D.; Touhami, O. Analytical model of slotted air-gap surface mounted permanent-magnet synchronous motor with magnet bars magnetized in the shifting direction. IEEE Trans. Magn. 2009, 45, 747–758. [Google Scholar] [CrossRef]
  17. Bellara, A.; Amara, Y.; Barakat, G.; Dakyo, B. Two-dimensional exact analytical solution of armature reaction field in slotted surface mounted PM radial flux synchronous machines. IEEE Trans. Magn. 2009, 45, 4534–4538. [Google Scholar] [CrossRef]
  18. Dubas, F.; Espanet, C.; Miraoui, A. Field diffusion equation in high-speed surface mounted permanent magnet motors, parasitic eddy-current losses. In Proceedings of the 6th International Symposium on Advanced Electromechanical Motion Systems, Lausanne, Switzerland, 27–29 September 2005; pp. 1–6. [Google Scholar]
  19. Liu, Z.J.; Li, J.T. Analytical solution of air-gap field in permanent-magnet motors taking into account the effect of pole transition over slots. IEEE Trans. Magn. 2007, 43, 3872–3883. [Google Scholar] [CrossRef]
  20. Zhu, Z.Q.; Howe, D.; Chan, C.C. Improved analytical model for predicting the magnetic field distribution in brushless permanent-magnet machines. IEEE Trans. Magn. 2002, 38, 229–238. [Google Scholar] [CrossRef]
  21. Rahideh, A.; Korakianitis, T. Analytical magnetic field distribution of slotless brushless machines with inset permanent magnets. IEEE Trans. Magn. 2011, 47, 1763–1774. [Google Scholar] [CrossRef]
  22. Dubas, F.; Rahideh, A. 2-D analytical PM eddy-current loss calculations in slotless PMSM equipped with surface-inset magnets. IEEE Trans. Magn. 2014, 50, 6300320. [Google Scholar] [CrossRef]
  23. Rahideh, A.; Korakianitis, T. Analytical open-circuit magnetic field distribution of slotless brushless permanent magnet machines with rotor eccentricity. IEEE Trans. Magn. 2011, 47, 4791–4808. [Google Scholar] [CrossRef]
  24. Rahideh, A.; Mardaneh, M.; Korakianitis, T. Analytical 2-D calculations of torque, inductance, and back-EMF for brushless slotless machines with surface inset magnets. IEEE Trans. Magn. 2013, 49, 4873–4884. [Google Scholar] [CrossRef]
  25. Rahideh, A.; Korakianitis, T. Analytical, armature reaction field distribution of slotless brushless machines with inset permanent magnets. IEEE Trans. Magn. 2012, 48, 2178–2191. [Google Scholar] [CrossRef]
  26. Atallah, K.; Zhu, Z.Q.; Howe, D.; Birch, T.S. Armature reaction field and winding inductances of slotless permanent-magnet brushless machines. IEEE Trans. Magn. 1998, 34, 3737–3744. [Google Scholar] [CrossRef]
  27. Pfister, P.D.; Perriard, Y. Slotless permanent-magnet machines: General analytical magnetic field calculation. IEEE Trans. Magn. 2011, 47, 1739–1752. [Google Scholar] [CrossRef]
  28. Holm, S.R.; Polinder, H.; Ferreira, J.A. Analytical modeling of a permanent-magnet synchronous machine in a flywheel. IEEE Trans. Magn. 2007, 43, 1955–1967. [Google Scholar] [CrossRef]
  29. Liu, Z.J.; Li, J.T. Accurate prediction of magnetic field and magnetic forces in permanent magnet motors using an analytical solution. IEEE Trans. Energy Convers. 2008, 23, 717–726. [Google Scholar] [CrossRef]
  30. Liu, Z.J.; Li, J.T.; Jiang, Q. An improved analytical solution for predicting magnetic forces in permanent magnet motors. J. Appl. Phys. 2008, 103, 07F135. [Google Scholar] [CrossRef]
  31. Amara, Y.; Reghem, P.; Barakat, G. Analytical prediction of eddy-current loss in armature windings of permanent magnet brushless AC machines. IEEE Trans. Magn. 2010, 46, 3481–3484. [Google Scholar] [CrossRef]
  32. Pfister, P.; Yin, X.; Fang, Y. Slotted Permanent-Magnet Machines: General Analytical Model of Magnetic Fields, Torque, Eddy Currents, and Permanent-Magnet Power Losses Including the Diffusion Effect. IEEE Trans. Magn. 2016, 52, 1–13. [Google Scholar] [CrossRef]
  33. Wu, L.J.; Zhu, Z.Q.; Staton, D.; Popescu, M.; Hawkins, D. Subdomain model for predicting armature reaction field of surface-mounted permanent-magnet machines accounting for tooth-tips. IEEE Trans. Magn. 2011, 47, 812–822. [Google Scholar] [CrossRef]
  34. Gysen, B.L.J.; Meessen, K.J.; Paulides, J.J.H.; Lomonova, E.A. General formulation of the electromagnetic field distribution in machines and devices using fourier analysis. IEEE Trans. Magn. 2010, 46, 39–52. [Google Scholar] [CrossRef]
  35. Zhu, Z.Q.; Howe, D.; Xia, Z.P. Prediction of open-circuit airgap field distribution in brushless machines having an inset permanent magnet rotor topology. IEEE Trans. Magn. 1994, 30, 98–107. [Google Scholar] [CrossRef]
  36. Jian, L.; Chau, K.T.; Gong, Y.; Yu, C.; Li, W. Analytical calculation of magnetic field in surface-inset permanent-magnet motors. IEEE Trans. Magn. 2009, 45, 4688–4691. [Google Scholar] [CrossRef]
  37. Zhu, Z.Q.; Ishak, D.; Howe, D.; Chen, J. Unbalanced magnetic forces in permanent-magnet brushless machines with diametrically asymmetric phase windings. IEEE Trans. Ind. Appl. 2007, 43, 1544–1553. [Google Scholar] [CrossRef]
  38. Sprangers, R.L.J.; Paulides, J.J.H.; Gysen, B.L.J.; Lomonova, E.A. Magnetic Saturation in Semi-Analytical Harmonic Modeling for Electric Machine Analysis. IEEE Trans. Magn. 2016, 52, 1–10. [Google Scholar] [CrossRef]
  39. Dubas, F.; Boughrara, K. New scientific contribution on the 2-D subdomain technique in polar coordinates: Taking into account of iron parts. Math. Comput. Appl. 2017, 22, 42. [Google Scholar] [CrossRef]
  40. Dubas, F.; Boughrara, K. New scientific contribution on the 2-D subdomain technique in Cartesian coordinates: Taking into account of iron parts. Math. Comput. Appl. 2017, 22, 17. [Google Scholar] [CrossRef]
  41. Djelloul-Khedda, Z.; Boughrara, K.; Dubas, F.; Kechroud, A.; Tikellaline, A. Analytical Prediction of Iron-Core Losses in Flux-Modulated Permanent-Magnet Synchronous Machines. IEEE Trans. Magn. 2019, 55, 1–12. [Google Scholar] [CrossRef]
  42. Djelloul-Khedda, Z.; Boughrara, K.; Dubas, F.; Ibtiouen, R. Nonlinear Analytical Prediction of Magnetic Field and Electromagnetic Performances in Switched Reluctance Machines. IEEE Trans. Magn. 2017, 53, 1–11. [Google Scholar] [CrossRef]
  43. Djelloul-Khedda, Z.; Boughrara, K.; Dubas, F.; Kechroud, A.; Souleyman, B. Semi-Analytical Magnetic Field Predicting in Many Structures of Permanent-Magnet Synchronous Machines Considering the Iron Permeability. IEEE Trans. Magn. 2018, 54, 1–21. [Google Scholar] [CrossRef]
  44. Roubache, L.; Boughrara, K.; Dubas, F.; Ibtiouen, R. New Subdomain Technique for Electromagnetic Performances Calculation in Radial-Flux Electrical Machines Considering Finite Soft-Magnetic Material Permeability. IEEE Trans. Magn. 2018, 54, 1–15. [Google Scholar] [CrossRef]
  45. Ben Yahia, M.; Boughrara, K.; Dubas, F.; Roubache, L.; Ibtiouen, R. Two-Dimensional Exact Subdomain Technique of Switched Reluctance Machines with Sinusoidal Current Excitation. Math. Comput. Appl. 2018, 23, 59. [Google Scholar]
  46. Roubache, L.; Boughrara, K.; Dubas, F.; Ibtiouen, R. Elementary subdomain technique for magnetic field calculation in rotating electrical machines with local saturation effect. Int. J. Comput. Math. Electr. Electron. Eng. 2018. [Google Scholar] [CrossRef]
  47. Boughrara, K.; Dubas, F.; Ibtiouen, R. 2-D Exact Analytical Method for Steady-State Heat Transfer Prediction in Rotating Electrical Machines. IEEE Trans. Magn. 2018, 54, 1–19. [Google Scholar] [CrossRef]
  48. Hannon, B.; Sergeant, P.; Dupré, L. Two-dimensional Fourier-based modeling of electric machines. In Proceedings of the 2017 IEEE International Electric Machines and Drives Conference, Miami, FL, USA, 21–24 May 2017; pp. 1–8. [Google Scholar]
  49. Lubin, T.; Rezzoug, A. 3-D Analytical Model for Axial-Flux Eddy-Current Couplings and Brakes Under Steady-State Conditions. IEEE Trans. Magn. 2015, 51, 1–12. [Google Scholar] [CrossRef]
  50. Lubin, T.; Rezzoug, A. Improved 3-D Analytical Model for Axial-Flux Eddy-Current Couplings with Curvature Effects. IEEE Trans. Magn. 2017, 53, 1–9. [Google Scholar] [CrossRef]
  51. Zarko, D.; Ban, D.; Lipo, T.A. Analytical solution for cogging torque in surface permanent magnet motors using conformal mapping. IEEE Trans. Magn. 2008, 44, 52–65. [Google Scholar] [CrossRef]
  52. Boughrara, K.; Zarko, D.; Ibtiouen, R.; Touhami, O.; Rezzoug, A. Magnetic field analysis of inset and surface-mounted permanent-magnet synchronous motors using Schwarz–Christoffel transformation. IEEE Trans. Magn. 2009, 45, 3166–3178. [Google Scholar] [CrossRef]
  53. Clemens, M.; Lang, J.; Teleaga, D.; Wimmer, G. Transient 3D magnetic field simulations with combined space and time mesh adaptivity for lowest order Whitney finite element formulations. IET Sci. Meas. Technol. 2009, 3, 377–383. [Google Scholar] [CrossRef]
  54. Liang, Y.; Bian, X.; Yang, L.; Wu, L. Numerical calculation of circulating current losses in stator transposition bar of large hydro-generator. IET Sci. Meas. Technol. 2015, 9, 485–491. [Google Scholar] [CrossRef]
  55. Mohammed, O.A.; Liu, S.; Liu, Z. FE-based physical phase variable model of PM synchronous machines under stator winding short circuit faults. IET Sci. Meas. Technol. 2007, 1, 12–16. [Google Scholar] [CrossRef]
Figure 1. Outer rotor surface inset brushless permanent magnet machine illustration.
Figure 1. Outer rotor surface inset brushless permanent magnet machine illustration.
Mca 24 00024 g001
Figure 2. Winding topologies (a) single layer alternate teeth wound (b) double layer all teeth wound non-overlapping (c) double-layer overlapping.
Figure 2. Winding topologies (a) single layer alternate teeth wound (b) double layer all teeth wound non-overlapping (c) double-layer overlapping.
Mca 24 00024 g002
Figure 3. Sub-region division according to the winding configuration: (a) whole slot is considered as one region and belongs to one coil, (b) whole slot is considered as one region and left side and right side of the slot belongs to different coils, (c) whole slot divided into upper and lower sub-regions and each part belongs to one coil.
Figure 3. Sub-region division according to the winding configuration: (a) whole slot is considered as one region and belongs to one coil, (b) whole slot is considered as one region and left side and right side of the slot belongs to different coils, (c) whole slot divided into upper and lower sub-regions and each part belongs to one coil.
Mca 24 00024 g003
Figure 4. Value of the current density in each sub-region of slot, according to represented winding configurations.
Figure 4. Value of the current density in each sub-region of slot, according to represented winding configurations.
Mca 24 00024 g004
Figure 5. The machine topology and winding configuration.
Figure 5. The machine topology and winding configuration.
Mca 24 00024 g005
Figure 6. Radial magnetic flux density due to just PMs in the middle of the air-gap, when the rotor position is set to zero.
Figure 6. Radial magnetic flux density due to just PMs in the middle of the air-gap, when the rotor position is set to zero.
Mca 24 00024 g006
Figure 7. Tangential magnetic flux density due to just PMs in the middle of the air-gap, when the rotor position is set to zero.
Figure 7. Tangential magnetic flux density due to just PMs in the middle of the air-gap, when the rotor position is set to zero.
Mca 24 00024 g007
Figure 8. Radial magnetic flux density due to just armature winding in the middle of the air-gap, when the current density of phase A is zero, current density of phase B is 4.33 A/mm2 and phase C is 4.33 A/mm2.
Figure 8. Radial magnetic flux density due to just armature winding in the middle of the air-gap, when the current density of phase A is zero, current density of phase B is 4.33 A/mm2 and phase C is 4.33 A/mm2.
Mca 24 00024 g008
Figure 9. Tangential magnetic flux density due to just armature winding in the middle of the air-gap, when the current density of phase A is zero, current density of phase B is 4.33 A/mm2 and phase C is 4.33 A/mm2.
Figure 9. Tangential magnetic flux density due to just armature winding in the middle of the air-gap, when the current density of phase A is zero, current density of phase B is 4.33 A/mm2 and phase C is 4.33 A/mm2.
Mca 24 00024 g009
Figure 10. Electromagnetic torque vs. rotor position.
Figure 10. Electromagnetic torque vs. rotor position.
Mca 24 00024 g010
Figure 11. Reluctance torque vs. rotor position.
Figure 11. Reluctance torque vs. rotor position.
Mca 24 00024 g011
Figure 12. Cogging torque vs. rotor position.
Figure 12. Cogging torque vs. rotor position.
Mca 24 00024 g012
Figure 13. Back-electromotive force (EMF) of the first phase and first line.
Figure 13. Back-electromotive force (EMF) of the first phase and first line.
Mca 24 00024 g013
Figure 14. Self and mutual inductance of a phase.
Figure 14. Self and mutual inductance of a phase.
Mca 24 00024 g014
Figure 15. Unbalanced magnetic forces just due to the PMs.
Figure 15. Unbalanced magnetic forces just due to the PMs.
Mca 24 00024 g015
Figure 16. Unbalanced magnetic forces just due to the armature reaction.
Figure 16. Unbalanced magnetic forces just due to the armature reaction.
Mca 24 00024 g016
Figure 17. Total unbalanced magnetic forces exerted on the rotor.
Figure 17. Total unbalanced magnetic forces exerted on the rotor.
Mca 24 00024 g017
Table 1. The sub-domains and related symbols.
Table 1. The sub-domains and related symbols.
Sub-DomainSymbolNumber of Sub-Regions
Magnet m 1 , 2 , , 2 p
Air-gap a 1
Slot-openingso 1 , 2 , , Q
Slot s l 1 , 2 , , Q
Table 2. Radial magnetization pattern and its Fourier series components.
Table 2. Radial magnetization pattern and its Fourier series components.
Magnetization PatternIllustrationRadial Waveform ComponentTangential Waveform ComponentCoefficient of the Radial Component M r w k = ( 1 ) w 1 2 + k B r e m μ 0 × Coefficient of the Tangential Component M θ w k = ( 1 ) w 1 2 + k B r e m μ 0 ×
Radial Magnetization Mca 24 00024 i001 4 w π sin ( w π α p 2 α r ) 0
Table 3. Parameters used in the case study.
Table 3. Parameters used in the case study.
ParametersUnitSymbolValue
Number of phases q 3
Number of the pole-pair p 4
Number of slots Q 9
Outer radius of the slots(mm) R s l 18
Outer radius of the slot-opening(mm) R s o 29
Stator radius(mm) R s 31
Magnet radius(mm) R m 32
Radius of the rotor back iron(mm) R r 38
Axial length(m)Ls0.1
Span of the slot(rad) δ 0.6
Span of the slot-opening(rad) β 0.3
Pole arc to pole pith of the magnet α p 0.85
Ratio of the rotor back iron to the pole pitch α r 0.85
Remanence of magnet(T) B r e m 1
Relative permeability of the magnet μ r 1.05
Number of harmonics in each sub-domain N, U, V, W100

Share and Cite

MDPI and ACS Style

Vahaj, A.; Rahideh, A.; Moayed-Jahromi, H.; Ghaffari, A. Exact Two-Dimensional Analytical Calculations for Magnetic Field, Electromagnetic Torque, UMF, Back-EMF, and Inductance of Outer Rotor Surface Inset Permanent Magnet Machines. Math. Comput. Appl. 2019, 24, 24. https://0-doi-org.brum.beds.ac.uk/10.3390/mca24010024

AMA Style

Vahaj A, Rahideh A, Moayed-Jahromi H, Ghaffari A. Exact Two-Dimensional Analytical Calculations for Magnetic Field, Electromagnetic Torque, UMF, Back-EMF, and Inductance of Outer Rotor Surface Inset Permanent Magnet Machines. Mathematical and Computational Applications. 2019; 24(1):24. https://0-doi-org.brum.beds.ac.uk/10.3390/mca24010024

Chicago/Turabian Style

Vahaj, AmirAbbas, Akbar Rahideh, Hossein Moayed-Jahromi, and AliReza Ghaffari. 2019. "Exact Two-Dimensional Analytical Calculations for Magnetic Field, Electromagnetic Torque, UMF, Back-EMF, and Inductance of Outer Rotor Surface Inset Permanent Magnet Machines" Mathematical and Computational Applications 24, no. 1: 24. https://0-doi-org.brum.beds.ac.uk/10.3390/mca24010024

Article Metrics

Back to TopTop