# An Analytical Solution to the Electromagnetic Performance of an Ultra-High-Speed PM Motor

Authors: Wenjie Cheng; Zhikai Deng; Ling Xiao; Bin Zhong; Yanhua Sun
Journal: Journal of Electrical and Computer Engineering
 (2019)

---

## Abstract

Bipolar rotors with cylindrical permanent magnets (PMs), which are magnetized in parallel, are generally used in ultra-high-speed permanent magnet synchronous motors (PMSMs). Aiming for this type of PMSM, two methods, namely, the Poynting theorem and the Maxwell stress tensor, were used to deduce the analytic expression of the electromagnetic driving torque. A locked-rotor simulation by the finite element (FE) method verifies the analytical solution. In addition, the cogging torque, back-electromotive force (emf), armature reaction inductance, etc., were calculated, compared, and discussed. The no-load back-emf was checked by a 60,000 rpm free coast experiment. The proposed analytical solution is very convenient to guide the design of the ultra-high-speed PMSM and check its performance.


---

## Body

## 1. Introduction

With advancements in new materials, electrical and electronics, control technology, etc., the speed limits of rotating machines are being pushed higher, and ultra-high-speed motors are increasingly used in industry applications [1]. Particularly in recent times, the high-speed PMSM has attracted more attention due to its high efficiency, simple structure, high power density, and small size [2, 3]; for example, it can be used as the driving motor for the fuel cell air blower [4, 5], as an air centrifugal compressor [6], or as the generator of a micro-gas turbine [7].It is very significant to calculate the electromagnetic torque, because it can be used to check the motor power, estimate noise and vibration, and provide further guidance for the motor control. It is well known that the electromagnetic torque is generated from the interaction of the PM field and the armature reaction field; therefore, how to determine the above two magnetic fields becomes key. PMSMs often adopt a surface-mounted PM rotor, and this kind of rotor, as seen in most relevant literatures [8, 9], has a common characteristic, which is that the magnetic pole arc is less than a pole pitch (the pole embrace is less than 1), which means the air space between the magnets is occupied by unmagnetized materials, such as air or aluminum alloy. Rotors of the ultra-high-speed PMSM are always subject to extreme conditions, such as high speed and high temperature. To improve the bending stiffness and strength of the rotor, ultra-high-speed PMSMs generally use a special rotor construction, in which a cylindrical PM magnetized in parallel is interference-fitted into a metallic sleeve; therefore, the rotor has single pole pairs, and each half of the magnet has a 180° electric pole arc (the pole embrace is equal to 1), as shown in Figure 1. As mentioned in the literature [10], existing analytical methods are successful only in certain circumstances and are likely to fail in others, for example, when the magnets have a 180° electric pole arc. In the literature [11] study, an open-circuit magnetostatic field solution to a 2-pole rotor with a diametrically magnetized sintered NdFeB magnet was given for a 20,000 rpm, 3-phase brushless PM dc motor. Recently, the literature [12] gave analytical solutions to the PM field of bipolar rotors with ring or cylindrical PMs that are magnetized in parallel. Nevertheless, most relevant literature studies on analytical models for predicting the armature reaction field distribution are successful for the motor, owing to the high permeability shaft [13]; however, few studies involve the hollow stator bore as shown in Figure 2(b).Figure 1
Cross section of the ultra-high-speed PMSM.Figure 2
Models of the electromagnetic field of the PMSMs. (a) Open-circuit PM field. (b) Armature reaction field.
(a)
(b)For ultra-high-speed PMSMs including a stator with overlapping windings and 2-pole parallel magnetized cylindrical PMs, analytical solutions to the magnetic field and corresponding electromagnetic torque are rarely reported at present, which has resulted in a prolonged motor design cycle to evaluating the motor performance. Due to the need of cooling and reduction of pulse torque, the ultra-high-speed PMSM generally adopts a large air gap; hereupon, the slot effect is very weak. Therefore, the assumption of “slotless stator + equivalent current sheet” can be introduced, and some convenient formulas of estimating motor performance are expected to be obtained.Hence, this paper will first deduce the magnetic field, then obtain the electromagnetic driving torque by using the Poynting theorem and Maxwell stress tensor, and then verify the correctness by using the FE method. In addition, the cogging torque, the back-electromotive force (emf), and armature reaction inductance are calculated, compared, and discussed. A 10 kW, 120,000 rpm prototype is designed and manufactured, and the no-load back-emf is checked by a 60,000 rpm free coast experiment.


## 2. Analytical Model

The cross section of the ultra-high-speed PMSMs is shown in Figure1, in which slotted stator and distribution winding are used, and the rotor is composed by installing a cylindrical parallel polarized PM into the metallic sleeve.For the sintered NdFeB magnet, the demagnetization curve is approximately a straight line, and the relative recoil permeabilityμr is approximate to that of air. μr was chosen from 1.029–1.100 in most literatures. To simplify the field calculation, the following assumptions are made:(1)
The PM has a linear second-quadrant demagnetization curve, and the permeability of the stator iron is infinite(2)
The winding end effect and rotor eddy current are neglected(3)
The PM field is calculated with a 2D solution by assuming a smooth stator surface; that is, the slotting is neglected, as shown in Figure2(a)(4)
Conductors in the stator slot can be equivalent to the current sheet distributing evenly along the slot opening, as shown in Figure2(b)In the calculation of the PM field, the slot-opening effect is neglected; therefore, the calculation model of the PM field is shown in Figure2(a), where R1 and R2 are the outer radius and the inner radius of the PM, and I and II represent the PM and the air region, respectively.The sleeve applies to a low permeability alloy, which can be physically equivalent to the air; therefore, in the calculation model of the armature reaction field, the stator bore seems to be full of air, as shown in Figure2(b), where b0 is the width of the slot opening.The loaded magnetic field is the superposition of the open-circuit PM and the armature reaction field, and the loaded air-gap flux density is defined as follows:(1)BIIrload=BIIrW+BIIrPM,BIIθload=BIIθW+BIIθPM.In equation (1), load, PM, and W represent the loaded, PM, and winding magnetic fields, respectively. r and θ represent the radial and circumferential directions, respectively.


### 2.1. Open-Circuit Magnetic Field

The literature [12] has deduced the solution to the PM field shown in Figure 2(a), which is given as follows:(2)BIIrPM=−μ0MrR121+R22r−2R12μr−1−R22μr+1cosθ,BIIθPM=μ0MrR12−1+R22r−2R12μr−1−R22μr+1sinθ,where Mr is the amplitude of the residual magnetization of the PM, μr is the relative recoil permeability of the PM, μ0 is the vacuum permeability, and r and θ represent the radial and circumferential positions, respectively.


### 2.2. Armature Reaction Field

The governing equation of the air-gap region is as follows:(3)∂ϕ2∂2r+∂ϕr∂r+∂ϕ2r2∂θ2=0,where ϕ is a scalar magnetic potential, and the relation of ϕ and H (magnetic field intensity) is expressed as follows:(4)Hθ=−∂ϕr∂θ,Hr=−∂ϕ∂r.In the inner radius of the stator,(5)Hθr=R2=−J,where J is the Fourier series expansion of the current sheet. At the center of the stator, the H should meet the following:(6)Hθr⟶0≠∞,Hrr⟶0≠∞.Taking into account the symmetry of the magnetic field, the single conductor field can be solved according to equations (3)–(6). The method to solve the armature reaction field is as follows: the single conductor field is solved first, and then the single-turn coil field can be obtained by adding two single conductor fields at a polar distance apart from each other. Whereupon, the single-phase winding field is superpositioned onto fields that are excited by the multiple single-turn coils, and finally, the armature reaction field can be calculated by superimposing the fields of three-phase winding.For distributed winding, multiple single-turn coil fields that are a pitch apart from each other can be equivalent to a concentrated winding by taking into account the distribution coefficient and the pitch-shortening coefficient. The above-mentioned treatment will lightly change higher harmonic components of the armature reaction field but vastly simplify the solution of the single-phase winding field with good accuracy. In addition, it is assumed that instantaneous three-phase currentsia,ib,ic only involve fundamental components, and they can be written as follows:(7)ia=ImcosPωt+β,ib=ImcosPωt+β−2π3,ic=ImcosPωt+β−4π3,where Im is the current amplitude of a single conductor, P is the number of pole pairs, ω is the synchronous angular velocity, and β is the current phase.The armature reaction field in Figure2(b) is expressed as follows:(8)BIIrW=−2qNsμ0π∑v=1∞Kdpv2v−1Ksov2v−1⋅Fo2v−1,r⋅−1v⋅iacos2v−1α−γ+ibcos2v−1α−γ−2π3+iccos2v−1α−γ−4π3,BIIαW=−2qNsμ0π∑v=1∞Kdpv2v−1Ksov2v−1⋅Fo2v−1,r⋅−1v−1⋅iasin2v−1α−γ+ibsin2v−1α−γ−2π3+icsin2v−1α−γ−4π3,where Ksovv=sinvb0/2R2/vb0/2R2 is the slot-opening coefficient, which is a function of the spatial harmonic order v.For,v=rv−1/R2v reflects the influence of the air-gap thickness on the radial/circumferential flux density of the v order harmonic wave. α is the angle relative to the horizontal direction, as shown in Figure 2(b), and is positive in the counter clockwise direction. γ is an angle between the axis of winding A and the d-axis. q is the number of slots per phase per pole, and Ns is the number of conductors per slot.


### 2.3. Electromagnetic Torque and Cogging Torque

In this part, two methods, namely, the Poynting theorem and the Maxwell stress tensor, were used to deduce the analytic expression of the electromagnetic torque. For the 2D field, the electric field strength vectorE in the air-gap region only includes the axial component Ez; therefore, the Poynting vector is given by the following:(9)∏r=−EzHθ,∏θ=EzHr,∏z=0,where Ez=μ0∫∂Hθ/∂tdr. The electromagnetic power transmitting from the stator to the rotor is mainly dependent on the radial component of the Poynting vector Πr which can be calculated from equation (10) as follows:(10)PM=Temω=lef∫02π∏rrdθ,where lef is the length of the stator iron. Taking the horizontal axis in Figure 2(b) as the reference position, the PM field is changed as follows:(11)BIIrPM=−μ0MrR121+R22r−2R12μr−1−R22μr+1⋅cosφ,BIIθPM=μ0MrR12−1+R22r−2R12μr−1−R22μr+1⋅sinφ,where φ=θ−π/2=α. It is assumed that the electromagnetic power is mainly transmitted by the fundamental magnetic field of armature reaction field. When the number of pole pairs equals 1, by substituting equation (7) into equation (8), the armature reaction field will become as follows, using the properties of trigonometric functions:(12)BIIrW=μ0G1cosωt−φ+δ,BIIθW=μ0G1sinωt−φ+δ,where δ=β+γ represents the angle between the axis of armature reaction field and the axis of PM field and G1=3qNsKdp1Kso1Im/πR2. Equation (12) shows that the armature reaction field is a travelling wave, but the PM field is a standing wave shown by equation (11). To make the PM field rotate counterclockwise synchronously relative to the armature reaction field, equation (11) can written as the following:(13)BIIrPM=μ0−G2r−2−G3⋅cosφ−ωt,BIIθPM=μ0G2r−2−G3⋅sinφ−ωt,where G2=MrR12R22/R12μr−1−R22μr+1 and G3=MrR12/R12μr−1−R22μr+1. Therefore, the loaded gap magnetic intensity can be obtained as follows:(14)HIIr=G1cosωt−φ+δ+−G2r−2−G3⋅cosφ−ωt,HIIθ=G1sinωt−φ+δ+G2r−2−G3⋅sinφ−ωt.Furthermore,Ez can be obtained by HIIθ and is expressed as the following:(15)Ez=μ0ωrG1cosωt−φ+δ+G2r−2+G3cosφ−ωt.By substituting equation (15) into equation (9), Πr will be calculated, and according to equation (10), the electromagnetic torque can be obtained as follows:(16)Tem=−2πμ0lefG1G2sinδ,δ∈0,π.In the 2D electromagnetic field,fθ=BrBθ/μ0, and fθ is the circumferential electromagnetic force density acting on the stator or rotor. Therefore, the electromagnetic torque of the rotor, which is shown in Figure 1, can be written as the following:(17)Tem=2μ0lef∫0πHIIrHIIθr2dθ.By substituting equation (14) into equation (17), the expression of the electromagnetic torque Tem is as follows:(18)Tem=−2πμ0lefG1G3r2sinδ,δ∈0,π,≈6πμ0qKdp1Kso1NsIm2πR2R1R22Mrlefr2sinδ.In equation (18), when r = R2, equation (18) is equal to equation (16), which indicates that the torque acting on the rotor surface generated by the Maxwell stress equals the torque acting on the current sheet generated by the ampere force. The results show that the electromagnetic torque is proportional to the length of stator iron lef, the ampere turns per unit length of the stator inner surface NsIm/2πR2, and the magnetization Mr, while it is proportional to the square of the ratio R1/R22 (PM radius to stator inner radius).The cogging torque arises from the interaction of the magnetomotive force (mmf) harmonics and the airgap permeance harmonics; hence, it exists in almost all types of motors in which the airgap permeance is not constant. The cogging torque will cause undesirable speed pulsations, vibrations, possible resonance, and acoustic noise [10]. For the motor shown in Figure (1), because the radial and tangential components of the magnetic flux density vector, as presented in equation (11), are pure sinusoidal, the mmf harmonics only contain the fundamental component, and the cogging torque is caused by stator slots. In fact, the cogging torque can also be calculated by equation (17); however, HIIr and HIIθ should be replaced by HIIrPM and HIIθPM, respectively. In the paper, HIIrPM and HIIθPM are obtained by the hypothesis of a smooth stator surface; hereupon, they cannot reveal the slot effects. However, slot effects can be considered analytically by the subdomain technique. Works in the literature [10] have proposed a simple method to predict the cogging torque based on the analytical calculation of the air-gap field distribution and the net lateral force acting on the teeth.


### 2.4. Back-Electromotive Force

The flux through a phase winding excited by PM and by winding, respectively, can be written as follows:(19)Φpm=2BIIrPMmaxπτplef,ΦW=2BIIrWmaxπτplef,where τp is the pole pitch and lef is the length of stator iron. The back-emf amplitude in each phase winding, assuming a 120-degree coil span, is then obtained as follows:(20)E=4.44fkdpvNΦload,where f is the angular frequency of the rotor, kdpv is the winding coefficient, N is the number of turns per phase, and Φload=Φpm+ΦW. When ΦW = 0, it corresponds to an open-circuit situation, and the back-emf becomes the no-load back-emf E0.


### 2.5. Armature Reaction Inductance

The constant-power transformation from the three-phase coordinates to the d-q coordinate is given by the following:(21)idiqi0=23cosθcosθ−2π3cosθ+2π3−sinθ−sinθ−2π3−sinθ+2π3121212iAiBiC,where θ=∫ωdt represents the rotor position, which is an angle between the d-axis and winding axis of phase A. When the above two axes coincide, θ=0. Whereupon, according to equation (21), there are(22a)idiq=CiAiBiC,(22b)iAiBiC=CTidiq,where C=2/31−1/21/203/2−3/2 and CT and C−1 are transposed matrix and inverse matrix of C, respectively. Notably, there is CT=C−1. According to equation (22a), we have(23)id=23iA−12iB−12iC,iq=2332iB−32iC.θ
=
0, which means Pωt=0 in equation (7). Substituting (7) into equation (23) yields(24)id=3IRmscosβ,iq=3IRmssinβ,where IRms is the RMS of phase current, IRms=Im/2, and β denotes the electrical angle of the phase current vector beyond the d-axis. Let β = 0; by substituting equation (24) into equation (22b) to get the three-phase AC iA,iB,iC, when the three-phase AC flows into the stator windings, it will generate the d-axis armature reaction field. Obviously, the stator iron saturation is most likely to occur when the direction of the armature reaction field is in line with that of the PM field; accordingly, we only consider this situation.Denoting the inductance matrix in theabc coordinate system by(25)LABC=LAALABLACLBALBBLBCLCALCBLCC.The flux linkage in theabc coordinate system can be represented as(26)λAλBλC=LABCiAiBiC.The flux linkage in thedq coordinate system is also applied to equation (22a), that is,(27)λdλq=CλAλBλC.Substituting equation (26) and equation (22b) into equation (27) yields the following:(28)λdλq=CLABCCTidiq=Ldqidiq,where Ldq=LdLmdLmqLq=CLABCCT; it represents the inductance matrix in the dq coordinate system. If the end leakage is not considered and other leakage factors such as higher harmonics and slot are also ignored, then Ld and Lq will represent the d-axis and q-axis armature reaction inductances, namely, Lad and Laq, respectively. From the FEM, LABC will be obtained, and then it will turn into Ldq by using the above transformation.L
ad and Laq can also be calculated by analytical method, and the process is as follows.The load flux linkage of phase A is calculated as(29)λA=NR2lef∫−αy/2αy/2BIIrloadR2,θdθ,where N is the number of turns per coil of one phase winding and αy is the coil pitch angle. The load flux linkage of phase B and C is expressed as(30)λB=λAcos2π3,λC=λAcos4π3.According to the above-mentionedλA,λB,λC, from equation (27), we will get λd,λq. And then, there is(31)Ld=λdid.Similarly, neglecting some leakage factors,Ld calculated by equation (31), corresponds to Lad.


## 2.1. Open-Circuit Magnetic Field

The literature [12] has deduced the solution to the PM field shown in Figure 2(a), which is given as follows:(2)BIIrPM=−μ0MrR121+R22r−2R12μr−1−R22μr+1cosθ,BIIθPM=μ0MrR12−1+R22r−2R12μr−1−R22μr+1sinθ,where Mr is the amplitude of the residual magnetization of the PM, μr is the relative recoil permeability of the PM, μ0 is the vacuum permeability, and r and θ represent the radial and circumferential positions, respectively.


## 2.2. Armature Reaction Field

The governing equation of the air-gap region is as follows:(3)∂ϕ2∂2r+∂ϕr∂r+∂ϕ2r2∂θ2=0,where ϕ is a scalar magnetic potential, and the relation of ϕ and H (magnetic field intensity) is expressed as follows:(4)Hθ=−∂ϕr∂θ,Hr=−∂ϕ∂r.In the inner radius of the stator,(5)Hθr=R2=−J,where J is the Fourier series expansion of the current sheet. At the center of the stator, the H should meet the following:(6)Hθr⟶0≠∞,Hrr⟶0≠∞.Taking into account the symmetry of the magnetic field, the single conductor field can be solved according to equations (3)–(6). The method to solve the armature reaction field is as follows: the single conductor field is solved first, and then the single-turn coil field can be obtained by adding two single conductor fields at a polar distance apart from each other. Whereupon, the single-phase winding field is superpositioned onto fields that are excited by the multiple single-turn coils, and finally, the armature reaction field can be calculated by superimposing the fields of three-phase winding.For distributed winding, multiple single-turn coil fields that are a pitch apart from each other can be equivalent to a concentrated winding by taking into account the distribution coefficient and the pitch-shortening coefficient. The above-mentioned treatment will lightly change higher harmonic components of the armature reaction field but vastly simplify the solution of the single-phase winding field with good accuracy. In addition, it is assumed that instantaneous three-phase currentsia,ib,ic only involve fundamental components, and they can be written as follows:(7)ia=ImcosPωt+β,ib=ImcosPωt+β−2π3,ic=ImcosPωt+β−4π3,where Im is the current amplitude of a single conductor, P is the number of pole pairs, ω is the synchronous angular velocity, and β is the current phase.The armature reaction field in Figure2(b) is expressed as follows:(8)BIIrW=−2qNsμ0π∑v=1∞Kdpv2v−1Ksov2v−1⋅Fo2v−1,r⋅−1v⋅iacos2v−1α−γ+ibcos2v−1α−γ−2π3+iccos2v−1α−γ−4π3,BIIαW=−2qNsμ0π∑v=1∞Kdpv2v−1Ksov2v−1⋅Fo2v−1,r⋅−1v−1⋅iasin2v−1α−γ+ibsin2v−1α−γ−2π3+icsin2v−1α−γ−4π3,where Ksovv=sinvb0/2R2/vb0/2R2 is the slot-opening coefficient, which is a function of the spatial harmonic order v.For,v=rv−1/R2v reflects the influence of the air-gap thickness on the radial/circumferential flux density of the v order harmonic wave. α is the angle relative to the horizontal direction, as shown in Figure 2(b), and is positive in the counter clockwise direction. γ is an angle between the axis of winding A and the d-axis. q is the number of slots per phase per pole, and Ns is the number of conductors per slot.


## 2.3. Electromagnetic Torque and Cogging Torque

In this part, two methods, namely, the Poynting theorem and the Maxwell stress tensor, were used to deduce the analytic expression of the electromagnetic torque. For the 2D field, the electric field strength vectorE in the air-gap region only includes the axial component Ez; therefore, the Poynting vector is given by the following:(9)∏r=−EzHθ,∏θ=EzHr,∏z=0,where Ez=μ0∫∂Hθ/∂tdr. The electromagnetic power transmitting from the stator to the rotor is mainly dependent on the radial component of the Poynting vector Πr which can be calculated from equation (10) as follows:(10)PM=Temω=lef∫02π∏rrdθ,where lef is the length of the stator iron. Taking the horizontal axis in Figure 2(b) as the reference position, the PM field is changed as follows:(11)BIIrPM=−μ0MrR121+R22r−2R12μr−1−R22μr+1⋅cosφ,BIIθPM=μ0MrR12−1+R22r−2R12μr−1−R22μr+1⋅sinφ,where φ=θ−π/2=α. It is assumed that the electromagnetic power is mainly transmitted by the fundamental magnetic field of armature reaction field. When the number of pole pairs equals 1, by substituting equation (7) into equation (8), the armature reaction field will become as follows, using the properties of trigonometric functions:(12)BIIrW=μ0G1cosωt−φ+δ,BIIθW=μ0G1sinωt−φ+δ,where δ=β+γ represents the angle between the axis of armature reaction field and the axis of PM field and G1=3qNsKdp1Kso1Im/πR2. Equation (12) shows that the armature reaction field is a travelling wave, but the PM field is a standing wave shown by equation (11). To make the PM field rotate counterclockwise synchronously relative to the armature reaction field, equation (11) can written as the following:(13)BIIrPM=μ0−G2r−2−G3⋅cosφ−ωt,BIIθPM=μ0G2r−2−G3⋅sinφ−ωt,where G2=MrR12R22/R12μr−1−R22μr+1 and G3=MrR12/R12μr−1−R22μr+1. Therefore, the loaded gap magnetic intensity can be obtained as follows:(14)HIIr=G1cosωt−φ+δ+−G2r−2−G3⋅cosφ−ωt,HIIθ=G1sinωt−φ+δ+G2r−2−G3⋅sinφ−ωt.Furthermore,Ez can be obtained by HIIθ and is expressed as the following:(15)Ez=μ0ωrG1cosωt−φ+δ+G2r−2+G3cosφ−ωt.By substituting equation (15) into equation (9), Πr will be calculated, and according to equation (10), the electromagnetic torque can be obtained as follows:(16)Tem=−2πμ0lefG1G2sinδ,δ∈0,π.In the 2D electromagnetic field,fθ=BrBθ/μ0, and fθ is the circumferential electromagnetic force density acting on the stator or rotor. Therefore, the electromagnetic torque of the rotor, which is shown in Figure 1, can be written as the following:(17)Tem=2μ0lef∫0πHIIrHIIθr2dθ.By substituting equation (14) into equation (17), the expression of the electromagnetic torque Tem is as follows:(18)Tem=−2πμ0lefG1G3r2sinδ,δ∈0,π,≈6πμ0qKdp1Kso1NsIm2πR2R1R22Mrlefr2sinδ.In equation (18), when r = R2, equation (18) is equal to equation (16), which indicates that the torque acting on the rotor surface generated by the Maxwell stress equals the torque acting on the current sheet generated by the ampere force. The results show that the electromagnetic torque is proportional to the length of stator iron lef, the ampere turns per unit length of the stator inner surface NsIm/2πR2, and the magnetization Mr, while it is proportional to the square of the ratio R1/R22 (PM radius to stator inner radius).The cogging torque arises from the interaction of the magnetomotive force (mmf) harmonics and the airgap permeance harmonics; hence, it exists in almost all types of motors in which the airgap permeance is not constant. The cogging torque will cause undesirable speed pulsations, vibrations, possible resonance, and acoustic noise [10]. For the motor shown in Figure (1), because the radial and tangential components of the magnetic flux density vector, as presented in equation (11), are pure sinusoidal, the mmf harmonics only contain the fundamental component, and the cogging torque is caused by stator slots. In fact, the cogging torque can also be calculated by equation (17); however, HIIr and HIIθ should be replaced by HIIrPM and HIIθPM, respectively. In the paper, HIIrPM and HIIθPM are obtained by the hypothesis of a smooth stator surface; hereupon, they cannot reveal the slot effects. However, slot effects can be considered analytically by the subdomain technique. Works in the literature [10] have proposed a simple method to predict the cogging torque based on the analytical calculation of the air-gap field distribution and the net lateral force acting on the teeth.


## 2.4. Back-Electromotive Force

The flux through a phase winding excited by PM and by winding, respectively, can be written as follows:(19)Φpm=2BIIrPMmaxπτplef,ΦW=2BIIrWmaxπτplef,where τp is the pole pitch and lef is the length of stator iron. The back-emf amplitude in each phase winding, assuming a 120-degree coil span, is then obtained as follows:(20)E=4.44fkdpvNΦload,where f is the angular frequency of the rotor, kdpv is the winding coefficient, N is the number of turns per phase, and Φload=Φpm+ΦW. When ΦW = 0, it corresponds to an open-circuit situation, and the back-emf becomes the no-load back-emf E0.


## 2.5. Armature Reaction Inductance

The constant-power transformation from the three-phase coordinates to the d-q coordinate is given by the following:(21)idiqi0=23cosθcosθ−2π3cosθ+2π3−sinθ−sinθ−2π3−sinθ+2π3121212iAiBiC,where θ=∫ωdt represents the rotor position, which is an angle between the d-axis and winding axis of phase A. When the above two axes coincide, θ=0. Whereupon, according to equation (21), there are(22a)idiq=CiAiBiC,(22b)iAiBiC=CTidiq,where C=2/31−1/21/203/2−3/2 and CT and C−1 are transposed matrix and inverse matrix of C, respectively. Notably, there is CT=C−1. According to equation (22a), we have(23)id=23iA−12iB−12iC,iq=2332iB−32iC.θ
=
0, which means Pωt=0 in equation (7). Substituting (7) into equation (23) yields(24)id=3IRmscosβ,iq=3IRmssinβ,where IRms is the RMS of phase current, IRms=Im/2, and β denotes the electrical angle of the phase current vector beyond the d-axis. Let β = 0; by substituting equation (24) into equation (22b) to get the three-phase AC iA,iB,iC, when the three-phase AC flows into the stator windings, it will generate the d-axis armature reaction field. Obviously, the stator iron saturation is most likely to occur when the direction of the armature reaction field is in line with that of the PM field; accordingly, we only consider this situation.Denoting the inductance matrix in theabc coordinate system by(25)LABC=LAALABLACLBALBBLBCLCALCBLCC.The flux linkage in theabc coordinate system can be represented as(26)λAλBλC=LABCiAiBiC.The flux linkage in thedq coordinate system is also applied to equation (22a), that is,(27)λdλq=CλAλBλC.Substituting equation (26) and equation (22b) into equation (27) yields the following:(28)λdλq=CLABCCTidiq=Ldqidiq,where Ldq=LdLmdLmqLq=CLABCCT; it represents the inductance matrix in the dq coordinate system. If the end leakage is not considered and other leakage factors such as higher harmonics and slot are also ignored, then Ld and Lq will represent the d-axis and q-axis armature reaction inductances, namely, Lad and Laq, respectively. From the FEM, LABC will be obtained, and then it will turn into Ldq by using the above transformation.L
ad and Laq can also be calculated by analytical method, and the process is as follows.The load flux linkage of phase A is calculated as(29)λA=NR2lef∫−αy/2αy/2BIIrloadR2,θdθ,where N is the number of turns per coil of one phase winding and αy is the coil pitch angle. The load flux linkage of phase B and C is expressed as(30)λB=λAcos2π3,λC=λAcos4π3.According to the above-mentionedλA,λB,λC, from equation (27), we will get λd,λq. And then, there is(31)Ld=λdid.Similarly, neglecting some leakage factors,Ld calculated by equation (31), corresponds to Lad.


## 3. Example Validation

To verify the analytic solution to the electromagnetic torque, a locked-rotor simulation is carried out by the FEM method. The motor parameters are listed in Table1.Table 1
Parameters of the high-speed PMSM.
Symbol
Quantity
Values
R
1
PM radius
12 × 10−3 m
R
2
Inner radius of stator
17.4 × 10−3 m
Q
Slot numbers
12
l
ef
Core length
55 × 10−3 m
M
Magnetization
820000 A/m
ω
Angular velocity
12560 rad/s
N
s
Conductors per slot
20
μ
r
Relative recoil permeability
1.038
P
Number of pole pairs
1In the simulation, the windings are fed by a symmetrical three-phase alternating current to produce a rotating armature reaction field, but the rotor speed is set to zero; that is, the rotor is locked. When the armature reaction field rotates relative to the rotor with the angular velocityω, the locked-rotor torque will generate and vary periodically with δ (an angle between two axes of the armature reaction field and PM field, and it changes periodically from 0 to 360 degrees). The locked-rotor phenomenon can be applied to verify the electromagnetic torque of equation (16) and equation (18).When the motor is fed with currents as shown in equation (7), the locked-rotor torques are as shown in Figure 3. Three kinds of input currents of 2000 Hz with different amplitudes, such as 10 A, 30 A, and 50 A are considered. Figure 3(a) shows the three-phase currents with a 10 A peak value in FEM. The locked-rotor torque curves obtained by FEM are shown in Figure 3(b), which shows that the locked-rotor torque is also sinusoidal with time, and its peak value is increased with the increasing amplitude of the phase current. A comparison of the results of FEM and the analytical model is illustrated in Figure 3(c), in which it can be found that the analytical solution is in good agreement with the FEM solution.Figure 3
The results of the locked-rotor. (a) The three-phase alternating currents. (b) The torque vs time. (c) Result comparison of FEM and the analytical model.
(a)
(b)
(c)Peak torques obtained by the FE model are compared with those from the derived analytical equations. The results, which are listed in Table2, show that the errors are less than 5%, which proves the correctness of equation (16) and equation (18).Table 2
Peak torque comparison between FE and the theory.
I
m (A)
FEM (N·m)
Analytical (N·m)
Error (%)
10
0.266
0.271
1.87
30
0.802
0.815
1.62
50
1.333
1.357
1.76The magnetic field results that were calculated by FEM are shown in Figure4, where red and blue conductors indicate the current flows out of and into the cross section, respectively. In addition, the PM has a parallel magnetization (the direction of the magnetization is vertically upward).Figure 4
Flux density vector field and current vector field of the FE model. (a) Fields at 0.0002 s. (b) Fields at 0.0004 s. (c) Fields at 0.0005 s.
(a)
(b)
(c)Because the PM field is much larger than the armature reaction field, the magnetic fields of the motor in Figure4 at different times are approximately the same. At 0.2 ms, the locked-rotor torque in Figure 3(b) is negative, and the corresponding motor magnetic field is shown in Figure 4(a), in which the armature reaction field lags behind the PM field by approximately 135 degrees, which just causes a negative torque. Moreover, the armature reaction field weakens the PM field, the maximum flux density occurs in the stator tooth, and its value is 1.38 T. At 0.4 ms, the locked-rotor torque equals zero, as shown in Figure 3(b), and the corresponding motor magnetic field is shown in Figure 4(b), where the armature reaction field coincides with the PM field, which verifies the null torque. At 0.4 ms, the maximum flux density of the stator tooth is 1.70 T. At 0.5 ms, the locked-rotor torque approaches the maximum, as shown in Figure 3(b), and the corresponding motor magnetic field is shown in Figure 4(c). It can be seen that the armature reaction field direction is left (90 degrees ahead of the PM field), which verifies the maximum torque. At 0.5 ms, the maximum flux density of the stator tooth is 1.60 T.Three methods, namely, the FM method, the air-gap permeance method [10], and the subdomain method, were used to estimate the cogging torque, and the results are shown in Figure 5. The air-gap permeance method is based on the assumption that the flux crosses the magnet and air gap in a straight line wherever a magnet faces a tooth. However, the above assumption will no longer be suitable for the motor mentioned in the paper, because the flux will bend in the air gap, as shown in Figure 5(a); the cogging torque calculated by the air-gap permeance method is small relative to that from FEM and the subdomain method, as shown in Figure 5(b). The cogging torque peak is approximately 370 μNm and is 0.05% of the rated torque (0.8 Nm), which is too small to be ignored. In fact, some degrees of harmonic components will exist in the polarized PM; thus, the cogging torque peak will increase.Figure 5
PM field and cogging torque. (a) PM field. (b) Comparison of the cogging torque of FEM and the analytical model.
(a)
(b)The stator tooth tips in relation to the width of the slot openings also influence the back-emf waveform. When the height of the tooth tips is small, magnetic saturation exists and the emf waveform tends to be a trapezoid. However, when the tooth tip height is too large, the emf waveform will deteriorate [11]. The motor consists of a stator having 12 teeth, which carry overlapping winding. The no-load back-emf waveform obtained from the FEM is displayed in Figure 6(a), and the emf waveform is almost purely sinusoidal, which is ideal for ultra-high-speed PMSMs. As shown in Figure 5(a), a stator tooth has the maximum flux density, 1.542 T, which does not reach the saturation region of the material. The comparison of the no-load back-emf calculated by equation (20) and that of FEM is shown in Figure 6(b), which manifests a good agreement of both methods. Figure 6(c) shows the prototype of the PMSM, in which the rotor is supported by gas foil bearings, and the experimental results are shown in Figure 6(d). A coast-down test under 60,000 r/min is carried out; that is, the frequency converter will shut down when the rotor is stable at 60,000 r/min, and then the rotor speed will decrease to zero gradually by overcoming air friction. In this process, the stator has a back-emf, but no current; thereupon, the rotor is not subjected to the electromagnetic force. Hereupon, the curve of a no-load back-emf vs speed can be achieved. The experimental and design values of the no-load back-emf are 57.0 V and 59.6 V at 746 Hz, respectively, and the difference between them is within 4.6%.Figure 6
Results of the no-load back-emf. (a) No-load back-emf of FEM at 120,000 rpm. (b) Comparison of no-load back-emf at 120,000 rpm. (c) Prototype of the ultra-high-speed PMSM. (d) Experimental curves of the phase voltage and frequency.
(a)
(b)
(c)
(d)The results ofLad and Laq are shown in Figure 7. It can be found that Lad and Laq are same when current is small, and both of them decrease with the increase of the current. Because the d-axial magnetic circuit is more saturated, Lad descends faster. The magnetic saturation is neglected in the analytical method, whereas Lad not varies with the current.Figure 7
Curves ofLad and Laq vs the current.WhenIRms = 10 A and 100 A, the errors of the two methods are 9% and 2%, respectively. It shows that the results of the analytical method and FEM are very close. In the analytical method, only the fundamental harmonic is used to calculate the flux linkage. However, in the FEM, except for the fundamental harmonic, high-order harmonics are also taken into account in the flux linkage. Therefore, the result of the analytical method is smaller than that of FEM.It is worth noting that the stator tooth has the most serious saturation as shown in Figure8(a), and the flux densities reach 2.086 T and 1.880 T, when IRms = 100 A (2.8IN) and 50 A (1.4IN), respectively. IN denotes the rated current of the motor. The working points of stator tooth corresponding to the both working condition are illustrated in Figure 8(b). It shows that the working point is located at the end of the linear segment of the B-H curve when IRms = 1.4IN. It indicates that the magnetic saturation of the stator iron is not prominent when the prototype works under rated working conditions, and it applies to the proposed analytical method in this paper.Figure 8
Magnetic saturation condition of the stator iron. (a) Distribution of flux density forIRms = 100 A. (b) Working points of the stator tooth.
(a)
(b)In addition, in order to enhance the cooling of the rotor, the air gap should be chosen as large as possible, so that magnetic saturation in the tooth tips rarely occurs. Furthermore, the teeth in a 12-slot stator are diametrically symmetrical, the air-gap field distribution is symmetrical, and no unbalanced magnetic pull will exist if the rotor and stator axes are coincident. Although the rotor rotates eccentrically with respect to the stator, in fact, the eccentricity is too tiny due to the constraint of bearings to cause a notable unbalanced magnetic pull.


## 4. Conclusion

For the ultra-high-speed PMSM, including a stator with overlapping windings and the 2-pole parallel magnetized cylindrical PM, electromagnetic torques are rarely reported at present, which leads to a prolonged motor design cycle to evaluating motor performance. Therefore, the work of this paper is as follows:(1)
This paper deduces the analytical solution to the electromagnetic driving torque of ultra-high-speed PMSMs, which is verified by the FE method according to locked-rotor simulation. The results show that the electromagnetic torque is proportional to the length of stator iron, the ampere turns per unit length of the stator inner surface and the magnetization, while it is proportional to the square of the ratio, i.e., PM radius to stator inner radius.(2)
The cogging torque, back-emf, armature reaction inductance, and so on are calculated, compared, and discussed. For the PMSM studied in the paper, the air-gap permeance method was shown not to be applicable to the estimation of the cogging torque due to the fact that the flux bends in the air gap, and thus, a high-accuracy method (e.g., the subdomain method) should be adopted.(3)
A prototype is designed and manufactured, and the no-load back-emf is checked by a 60,000 rpm free coast experiment. The experimental and design values of the no-load back-emf are 57.0 V and 59.6 V at 746 Hz, respectively, and the difference between them is within 4.6%.



---

*Source: 1469637-2019-09-12.xml*