Research on the Axial Force of Conical-Rotor Permanent Magnet Synchronous Motors with Turbines

: The general method to suppress the axial force of the permanent magnet synchronous motor (PMSM) direct-drive turbine is to increase the number of balance devices, such as balance disks and special bearings, to counteract its inﬂuence, but this also leads to complex system structure and higher mechanical losses. Aiming to solve the above issue, this paper presents a novel PMSM structure with a conical-rotor (CR). Due to its adaptive equilibrium of axial force and simple structure of rotor with turbine, the CR-PMSM can help improve the system efﬁciency. Both surface-type and interior-type motors are analyzed, and the axial magnetic force of CR-PMSM is studied in detail. The 3-D ﬁnite-element method (FEM) is used to model and simulate the machine, and the magnetic-ﬁeld distribution, axial magnetic force and driving performance are obtained. Also, the control rule of d- axis current is analyzed to achieve the adaptive equilibrium of axial force. A 2.0 kW, 6000 r/min prototype motor is fabricated and tested to validate the theory.


Introduction
Thanks to their high energy transfer efficiency and high torque density, PMSM direct-drive high-speed turbines (such as turbo-compressors and turbo-expanders, etc.) have been widely accepted. However, to ensure the safe and reliable operation of the system, the axial force of the turbine under the impact of high-speed working-fluid flow (which could be as large as 145 N for a 1.5 kW, 6000 r/min small turbo-compressor, and 3100 N for a 350 kW, 8000 r/min large turbo-expander) should be balanced, which is usually accomplished by using mechanical devices such as balance disks, mechanical linkages, sliders and special bearings, which significantly increase the system complexity and bring more losses. Moreover, the axial force changes aperiodically due to the influence of temperature, pressure and flow of the working fluid, so the impeller blades will collide with the volute and experience vibration and friction, causing abnormal wear and tear of the bearings, balance plates and springs. As a result, air gap eccentricity faults often occur, which cause the distortion of the air gap magnetic flux, and then affect the spectral content of both the stator currents and the zero-sequence voltage component [1,2]. Especially when under load, the vibration caused by the rotor eccentricity can seriously affect the stable operation of the motor [3], so it is necessary to diagnose the air gap eccentricity, resulting in additional work efforts. Meanwhile, for designers, not only should the axial force be calculated accurately, but also the balanced device should be designed and selected strictly, which significantly increases the manufacturing cost and design difficulty. Practically, this traditional method cannot obtain a satisfactory balancing result. Similar issues can be found in different types of flow turbines such as turbomachinery, turbopumps, turbochargers and turbogenerators, etc. [4][5][6][7]. The common problem is the inability of balancing the device to quickly and adaptively suppress the axial forces under variable conditions, which results in shorter operating cycles.
Thus, this paper presents an axial force suppression method for motor-turbine systems by using a CR-PMSM instead of increasing the number of balance devices. The rotor of the CR-PSMSM can be moved outside in an axial direction which decreases the air-gap flux density and increases the effective space between the stator and the rotor, so the axial magnetic force of the motor can be adjusted by changing the current or the the air-gap length to balance the axial force of the turbine.
As a novel PMSM with special structure, the CR-PMSM has received the attention of researchers in recent years. For example, in [8] a 3-D finite element model was used to simulate the magnetic-field distribution of the air gap, flux weakening, and driving performance in CR-PMSM according to rotor axial displacement, instead of changing the current angle. The decrease of flux linkage can be achieved by moving the conical-rotor in the axial direction, however, the axial magnetic force was not analyzed. Reference [9] introduced a PMSM consisting of two conical air gap bearingless PM half-motors. The two half-motors are mounted on a single shaft, and each conical stator has two separate three-phase distributed windings with different magnetic pole pair counts. A 1 kW, 18,000 r/min prototype has been designed and built. However, the effect of q-axis current on axial force was not analyzed.
It must be also noted that the conical structure brings some additional serious challenges to the designer because the CR dimensions change axially. Research indicates that the value of the axial force depends on the PM air gap flux density and the d-axis current, besides the geometry. Equivalent magnetic circuit models under open circuit condition (OCC) are often used to investigate electric machines [10][11][12][13][14].
As there is almost no precedent for the use of conical-rotor machines in turbine systems, the study in this paper may provide a theoretical guide for the application of CR-PMSM in turbomachinery. In this paper, the performance of the CR-PMSM is studied in detail. First of all, the basic features of the motor, including surface-type and interior-type, are analyzed. In order to determine the specific method of the axial force suppression, an axial magnetic force changing ratio parameter is defined by theoretical analysis. Then, 3-D FEM is adopted to model the motor and analyze the electromagnetic field, axial magnetic force and driving performance with the axial displacement of the rotor. Finally, a 2.0 kW, 6000 r/min prototype CR-PMSM with a cone angle of 6.0 • is designed and built, and tests of the motor-turbine system are carried out to validate the theory.

Theoretical Analysis of CR-PMSM
The stator winding structure of the CR-PMSM is almost the same as that of the conventional PMSM. However, the two ends of the CR have different diameters, and both the yoke height and tooth width along the axial direction change linearly with the rotor diameter [15].
A schematic drawing of CR-PMSM with the cone angle α is shown in Figure 1, in which L s and R in_av are defined as the length of iron core and the average radius of stator, respectively. The CR is supported by cylindrical roller bearings and can be easily moved in the axial direction. The air gap is minimal (shown by dashed lines in Figure 1 when the conical rotor is in alignment with the stator (∆z = 0). When the axial displacement of a CR is produced (shown by solid lines in Figure 1), the stator and the rotor are misplaced (∆z > 0), and the effective air-gap length (∆δ) can be calculated by: ∆δ(α, ∆z) = δ + ∆z sin α (1) where, δ is the actual air-gap length, α is the cone angle of the rotor.

Theoretical Calculation for Magnetic Circuit
In order to obtain simple design principles, the following assumptions are proposed: Equation (1) the relative permeability of the iron cores is infinite and Equation (2) the leakage flux and the cogging effect are ignored. The air gap magnetic field is synthesized by PM magnetic field and armature reaction magnetic field if the saturation is not considered. According to the double reaction theory of synchronous motor, armature reaction can be equivalent to the superposition of the effect of the d-axis and the q-axis component [16].
According to the magnetic circuit method, the flux density Bδ under OCC can be expressed as: where, Bδ0, B1d; B1q, θ and Z are the no-load air-gap flux density, the d-axis armature reaction flux density, the q-axis armature reaction flux density, the circumferential angle and the axial position of the rotor, respectively. At a fixed axial position, the specific expressions of the flux densities are as follows: d dp 1d 0 1q q dp 1q ( ) sin( ) + + sin where, μ0 is the permeability of vacuum (μ0 = 4π × 10 −7 H/m); Br, hm, bδ and bm are the PM remanence, the PM thickness (length of the magnetization direction), the PM width and the width of air-gap per pole, respectively; Kdp, Ns, Id, Iq, F1d, F1q and p are the winding coefficient, turn numbers per phase, d-axis current, q-axis current, d-axis magnetic potential, q-axis magnetic potential and pole pairs, respectively. When the motor carries load, B1d is of the same phase as Bδ0, enhancing or weakening the magnetic field. While, B1q is in the orthogonal phase of Bδ0, enhancing the magnetic field and causing waveform change, and no-load of CR-PMSM, Equation (3) could be approximately simplified to be: where, , Fm is the magneto-motive force (MMF) in a closed magnetic circuit, Rm and Rδ are the reluctance in the PM and air gap, respectively.

Theoretical Calculation for Magnetic Circuit
In order to obtain simple design principles, the following assumptions are proposed: Equation (1) the relative permeability of the iron cores is infinite and Equation (2) the leakage flux and the cogging effect are ignored. The air gap magnetic field is synthesized by PM magnetic field and armature reaction magnetic field if the saturation is not considered. According to the double reaction theory of synchronous motor, armature reaction can be equivalent to the superposition of the effect of the d-axis and the q-axis component [16].
According to the magnetic circuit method, the flux density B δ under OCC can be expressed as: where, B δ0 , B 1d ; B 1q , θ and Z are the no-load air-gap flux density, the d-axis armature reaction flux density, the q-axis armature reaction flux density, the circumferential angle and the axial position of the rotor, respectively. At a fixed axial position, the specific expressions of the flux densities are as follows: where, µ 0 is the permeability of vacuum (µ 0 = 4π × 10 −7 H/m); B r , h m , b δ and b m are the PM remanence, the PM thickness (length of the magnetization direction), the PM width and the width of air-gap per pole, respectively; K dp , N s , I d , I q , F 1d , F 1q and p are the winding coefficient, turn numbers per phase, d-axis current, q-axis current, d-axis magnetic potential, q-axis magnetic potential and pole pairs, respectively. When the motor carries load, B 1d is of the same phase as B δ0 , enhancing or weakening the magnetic field. While, B 1q is in the orthogonal phase of B δ0 , enhancing the magnetic field and causing waveform change, and no-load of CR-PMSM, Equation (3) could be approximately simplified to be: where, Φ δ0 = F m R m +R δ , F m is the magneto-motive force (MMF) in a closed magnetic circuit, R m and R δ are the reluctance in the PM and air gap, respectively.
Among various types of the rotor topologies, surface-type and interior-type CR-PMSM, as shown in Figure 2, are studied and compared in this paper. Among various types of the rotor topologies, surface-type and interior-type CR-PMSM, as shown in Figure 2, are studied and compared in this paper. For the surface-type CR-PMSM as shown in Figure 2a, the width of PM is equal to the width of air-gap per pole (bδ/bm = 1), and the flux passes through the PM twice in a closed magnetic circuit. For the interior-type CR-PMSM as shown in Figure 2b, the width of PM can be selected without regard of the air-gap width, and the flux passes through the PM once. If the width of the PM is bigger than that of air gap (bδ/bm < 1), the same rotor displacement will have a greater influence on the air-gap flux density, as shown in Figure 3. In order to keep the stator slot width bslot and height hslot dimension the same in the axial direction (z axis), and a constant the radius of stator Rout. So, through the magnetic circuit method, the maximum and minimum radius of the stator Rin_m can be calculated by: where, Rin_min take the value of "−" and Rin_max take the value of "+". Then, the stator slot width bslot, tooth width btooth_m and yoke height hyoke_m can be expressed as: For the surface-type CR-PMSM as shown in Figure 2a, the width of PM is equal to the width of air-gap per pole (b δ /b m = 1), and the flux passes through the PM twice in a closed magnetic circuit. For the interior-type CR-PMSM as shown in Figure 2b, the width of PM can be selected without regard of the air-gap width, and the flux passes through the PM once. If the width of the PM is bigger than that of air gap (b δ /b m < 1), the same rotor displacement will have a greater influence on the air-gap flux density, as shown in Figure 3. Among various types of the rotor topologies, surface-type and interior-type CR-PMSM, as shown in Figure 2, are studied and compared in this paper. For the surface-type CR-PMSM as shown in Figure 2a, the width of PM is equal to the width of air-gap per pole (bδ/bm = 1), and the flux passes through the PM twice in a closed magnetic circuit. For the interior-type CR-PMSM as shown in Figure 2b, the width of PM can be selected without regard of the air-gap width, and the flux passes through the PM once. If the width of the PM is bigger than that of air gap (bδ/bm < 1), the same rotor displacement will have a greater influence on the air-gap flux density, as shown in Figure 3. In order to keep the stator slot width bslot and height hslot dimension the same in the axial direction (z axis), and a constant the radius of stator Rout. So, through the magnetic circuit method, the maximum and minimum radius of the stator Rin_m can be calculated by: where, Rin_min take the value of "−" and Rin_max take the value of "+". Then, the stator slot width bslot, tooth width btooth_m and yoke height hyoke_m can be expressed as: In order to keep the stator slot width b slot and height h slot dimension the same in the axial direction (z axis), and a constant the radius of stator R out . So, through the magnetic circuit method, the maximum and minimum radius of the stator R in_m can be calculated by: where, R in_min take the value of "−" and R in_max take the value of "+". Then, the stator slot width b slot , tooth width b tooth_m and yoke height h yoke_m can be expressed as: where, k s and Q s are the ratio of tooth width to pitch and the stator slot number, respectively. So, the flux density of stator tooth B tooth_m and yoke B yoke_m can be expressed as: As shown in Equation (7), the flux density of the stator tooth B tooth_m and yoke B yoke_m varies with the axial position and the cone angle of the rotor. Inevitably, the torque density of CR-PMSM will be reduced slightly indicated by the parameter (L s tan α). The flux density of yoke B yoke_m can be adjusted by the outer radius R out independently. The relationship between the flux density of tooth B tooth_m and the length-to-diameter is: In order to maintain a rational magnetic loading (such as B δ = 0.7 T, B tooth_max = 1.54 T, and k s = 0.4), the ratio of length to diameter should meet the condition that L s tanα/R in_av < 0.15.

Modeling of Axial Magnetic Force Produced for CR-PMSM
In general, the flux lines in the air gap length δ of a low saturated PMSM are nearly perpendicular to the surface of the stator [17]. Unlike general cylindrical motors, the flux density B δ in the conical air gap can be decomposed into a radial and an axial component. If no rotor eccentricity occurs, the sum of the radial forces acting between the stator and the rotor is zero. Only the axial force is acting on the stator and the rotor, respectively [18]. The axial magnetic force can be calculated by the virtual displacement method [19] or the Maxwell stress tensor method [20]. The latter is employed in this paper for investigating the magnetic force, and the force components per air gap area element can be determined by: where f r , f t and f z are respectively the radial, tangential and axial component of the force, of which f r and f z compose to the normal component; B δn and B δt are the normal and tangential component of the air gap flux density. As can be seen from Figure 4, the unit area dS on the rotor surface can be modeled and analyzed by:  In a low saturated PM machine, the flux lines in the air gap are nearly perpendicular to the outer surface of the rotor, which implies that Bδt << Bδn, so the force densities could be further simplified to: Substituting Equations (2), (3) and (10) into (11), yields: Terms with '·' on head in Equation (12) mean that they are vector quantities. Since the q-axis flux is in theory perpendicular to the d-axis flux, Equation (12) could be approximately simplified to be According to Equations (12) and (13), once the CR-PMSM rotates, an axial magnetic force will be produced, and it will change with the stator current.
When the CR-PMSM under no-load, Equation (13) could be approximately simplified to be: In conclusion, the axial magnetic force, no-load back EMF and electromagnetic torque are mainly relevant to the air gap magnetic field, air gap length, outer diameter of the CR, core length and cone angle, etc. [21] From Equations (2) and (13), once the geometry of the CR-PMSM is fixed, the axial force Fz can be magnetically enhanced or weakened by changing the amplitudes of the dq-axis currents [22,23]. Thus, the analysis of Bδ will be necessary for the calculation of Fz.
In order to determine the axial magnetic force performance of CR-PMSM, when the CR moved outside in axial direction Δz and enhanced (or weakened) by changing the amplitudes of the d-axis current, a parameter of axial magnetic force ratio FΔz should be expressed as: In a low saturated PM machine, the flux lines in the air gap are nearly perpendicular to the outer surface of the rotor, which implies that B δt << B δn , so the force densities could be further simplified to: Substituting Equations (2), (3) and (10) into (11), yields: Terms with '·' on head in Equation (12) mean that they are vector quantities. Since the q-axis flux is in theory perpendicular to the d-axis flux, Equation (12) could be approximately simplified to be According to Equations (12) and (13), once the CR-PMSM rotates, an axial magnetic force will be produced, and it will change with the stator current.
When the CR-PMSM under no-load, Equation (13) could be approximately simplified to be: In conclusion, the axial magnetic force, no-load back EMF and electromagnetic torque are mainly relevant to the air gap magnetic field, air gap length, outer diameter of the CR, core length and cone angle, etc. [21] From Equations (2) and (13), once the geometry of the CR-PMSM is fixed, the axial force F z can be magnetically enhanced or weakened by changing the amplitudes of the dq-axis currents [22,23]. Thus, the analysis of B δ will be necessary for the calculation of F z .
In order to determine the axial magnetic force performance of CR-PMSM, when the CR moved outside in axial direction ∆z and enhanced (or weakened) by changing the amplitudes of the d-axis current, a parameter of axial magnetic force ratio F ∆z should be expressed as: where, F ∆z = 0 and F ∆z_max are the no-load axial magnetic force at the rated speed when ∆z = 0 and ∆z = max, F id = 0 and F id_max are the load axial magnetic force at the rated speed when i d = 0 and i d = max.
As can be seen from Equation (15), when the cone angle of CR-PMSM is set, the axial magnetic force is only related to the axial displacement ∆z of the CR and the i d current of the stator, so Equation (15) could be approximately simplified to be: The axial magnetic force changing ratio F ∆z,id can be decomposed to the air gap flux density changing ratio C B and the effective air gap space decrease ratio C S : Equation (17) shows the relationship between the axial magnetic force changing ratio F ∆z,id and the CR axial displacement ∆z under the i d current of the stator. Meanwhile, 6 • is an appropriate cone angle, which not only keeps F ∆z,id an approximately linear relation with wide range of ∆z, but also results in the balance point of good f axial magnetic force changing ratio and not quite a large displacement.

3-D FE Analysis of CR-PMSM
The axial magnetic force can be obtained by 3-D FEM when the armature winding is with different d, q-axis currents [18]. The main geometry and electrical parameters of the prototype machine are provided in Table 1. In order to save calculation time, according to the structural symmetry of the motor, a quarter of 3-D model is established by manual modeling using ANSYS software (ANSYS 16.0, Ansys Inc., Canonsburg, PA, USA), as shown in Figure 5. For distinction, the end surface with the minimum rotor diameter of the machine is defined as side A, and the other end surface of the machine (with the maximum rotor diameter) is defined as side B.

Magnetic Field Distribution
By using 3-D FE simulation, the magnetic field of the motor with cone angle from 4°-8° is analyzed, as shown in Figure 6a-d.

Magnetic Field Distribution
By using 3-D FE simulation, the magnetic field of the motor with cone angle from 4 • -8 • is analyzed, as shown in Figure 6a-d.

Magnetic Field Distribution
By using 3-D FE simulation, the magnetic field of the motor with cone angle from 4°-8° is analyzed, as shown in Figure 6a-d.  It can be seen from Figure 6 that the presence of the cone angle will make the end CR magnetically dense. When the cone angle is 8 • , the degree of stator saturation is obvious, so it is necessary to properly design the cone angle to prevent the CR-PMSM from operating in an over-saturated state. When the cone angle is fixed, the longer the axial length of the core, the larger the diameter of the stator and rotor changes, and the more prone to saturation. Therefore, the core should not be designed too long. Figure 7 shows the 3-D magnetic-field distribution of the air-gap magnetic density B δ , which is decomposed into the radial component B δr and the axial component B δz . As the cone angle increases, the radial magnetic density amplitude decreases slightly, and the axial magnetic density rises slightly. This tendency can be clearly seen when the cone angle is 6 • and the cone angle is 8 • .
Energies 2017, 10, x FOR PEER REVIEW 9 of 20 It can be seen from Figure 6 that the presence of the cone angle will make the end CR magnetically dense. When the cone angle is 8°, the degree of stator saturation is obvious, so it is necessary to properly design the cone angle to prevent the CR-PMSM from operating in an over-saturated state. When the cone angle is fixed, the longer the axial length of the core, the larger the diameter of the stator and rotor changes, and the more prone to saturation. Therefore, the core should not be designed too long. Figure 7 shows the 3-D magnetic-field distribution of the air-gap magnetic density Bδ, which is decomposed into the radial component Bδr and the axial component Bδz. As the cone angle increases, the radial magnetic density amplitude decreases slightly, and the axial magnetic density rises slightly. This tendency can be clearly seen when the cone angle is 6° and the cone angle is 8°. The flux density distribution with different axial displacement of CR is shown in Figure 8. At the larger end of the stator bore, the magnetic density on the top of the stator tooth appears to be partially saturated, and the magnetic density on the top of the stator tooth is obviously reduced at the smaller end of the stator core.
The flux density distribution with different axial displacement of CR is shown in Figure 8. At the larger end of the stator bore, the magnetic density on the top of the stator tooth appears to be partially saturated, and the magnetic density on the top of the stator tooth is obviously reduced at the smaller end of the stator core. Energies 2017, 10, x FOR PEER REVIEW 10 of 20        The amplitudes of air-gap flux density B δ along the axial direction are shown in Figure 10. It can be seen that the amplitude of B δ will be decreased accordingly as the CR axial displacement, but the change gradient is smaller and smaller. Although the amplitude of the axial component B δz is small. Thus, the calculation of B δ will be necessary for the calculation of the magnetic axial force. The amplitudes of air-gap flux density Bδ along the axial direction are shown in Figure 10. It can be seen that the amplitude of Bδ will be decreased accordingly as the CR axial displacement, but the change gradient is smaller and smaller. Although the amplitude of the axial component Bδz is small. Thus, the calculation of Bδ will be necessary for the calculation of the magnetic axial force.

Analysis of the Axial Magnetic Force
It can be seen from Figure 11 that when α is small, FZ increases slowly with a small axial displacement, and when α is larger, FZ decreases linearly with the axial displacement. As can be seen from Figure 12, the axial magnetic force has only a quite slight increase with the increase of q-axis current, which confirms the previous analysis and conclusion in Equations (12) and (13). In addition, the magnetic axial force will increase with id.

Analysis of the Axial Magnetic Force
It can be seen from Figure 11 that when α is small, F Z increases slowly with a small axial displacement, and when α is larger, F Z decreases linearly with the axial displacement. The amplitudes of air-gap flux density Bδ along the axial direction are shown in Figure 10. It can be seen that the amplitude of Bδ will be decreased accordingly as the CR axial displacement, but the change gradient is smaller and smaller. Although the amplitude of the axial component Bδz is small. Thus, the calculation of Bδ will be necessary for the calculation of the magnetic axial force.

Analysis of the Axial Magnetic Force
It can be seen from Figure 11 that when α is small, FZ increases slowly with a small axial displacement, and when α is larger, FZ decreases linearly with the axial displacement. As can be seen from Figure 12, the axial magnetic force has only a quite slight increase with the increase of q-axis current, which confirms the previous analysis and conclusion in Equations (12) and (13). In addition, the magnetic axial force will increase with id.  As can be seen from Figure 12, the axial magnetic force has only a quite slight increase with the increase of q-axis current, which confirms the previous analysis and conclusion in Equations (12) and (13). In addition, the magnetic axial force will increase with i d . The amplitudes of air-gap flux density Bδ along the axial direction are shown in Figure 10. It can be seen that the amplitude of Bδ will be decreased accordingly as the CR axial displacement, but the change gradient is smaller and smaller. Although the amplitude of the axial component Bδz is small. Thus, the calculation of Bδ will be necessary for the calculation of the magnetic axial force.

Analysis of the Axial Magnetic Force
It can be seen from Figure 11 that when α is small, FZ increases slowly with a small axial displacement, and when α is larger, FZ decreases linearly with the axial displacement. As can be seen from Figure 12, the axial magnetic force has only a quite slight increase with the increase of q-axis current, which confirms the previous analysis and conclusion in Equations (12) and (13). In addition, the magnetic axial force will increase with id.  Based on the simulation model, the operating performance curves of the CR-PMSM can be obtained, as shown in Figure 13. The x-coordinate is the ratio of output power (P out ) to rated power (P N ), and the y-coordinate is normalized, including the power factor (cosϕ), efficiency (η), the ratio of current (I) and the ratio of torque (T e ). Based on the simulation model, the operating performance curves of the CR-PMSM can be obtained, as shown in Figure 13. The x-coordinate is the ratio of output power (Pout) to rated power (PN), and the y-coordinate is normalized, including the power factor (cosφ), efficiency (η), the ratio of current (I) and the ratio of torque (Te). More simulation results are presented in Section 4, comparing with the experimental results, to verify the design rationality.

Experimental Validation and Analysis
The technique flowchart for design, analysis and validation of the proposed CR-PMSM is described in Figure 14, which contains four steps. Steps 1 to 3 mainly focus on the design and simulation of the motor, which all belong to theoretical research and have been discussed in detail in the sections above. In this section, the operation of the prototype experiment and the related analysis results in detail will be introduced.

Test Rig of Prototype Machine
In order to evaluate the performance of developed CR-PMSM, experimental system is designed and built. As shown in Figure 15, the parameters that can be adjusted on experimental platform are the flow and the pressure of compressed air. The inlet of the turbocompressor is connected to the air source, and the flow rate of normal air into the turbine-compressor can be controlled by adjusting the valve. The inlet pipeline is equipped with a thermal resistance and mass flowmeter to measure the temperature and flow parameters of the air at the inlet. While, the outlet pipeline is equipped with a thermal resistance and a high-precision pressure gauge to measure the temperature and pressure of the compressed air. The compressed air is injected into the pressure tank, and under rated condition the output flow rate is 250 kg/h and the pressure is 0.15 MPa. The turbocompressor is directly connected to the shaft of the conical motor, driven by a special frequency converter, and the other end of the motor shaft is equipped with a rotor axial displacement device. More simulation results are presented in Section 4, comparing with the experimental results, to verify the design rationality.

Experimental Validation and Analysis
The technique flowchart for design, analysis and validation of the proposed CR-PMSM is described in Figure 14, which contains four steps. Steps 1 to 3 mainly focus on the design and simulation of the motor, which all belong to theoretical research and have been discussed in detail in the sections above. In this section, the operation of the prototype experiment and the related analysis results in detail will be introduced.

Test Rig of Prototype Machine
In order to evaluate the performance of developed CR-PMSM, experimental system is designed and built. As shown in Figure 15, the parameters that can be adjusted on experimental platform are the flow and the pressure of compressed air. The inlet of the turbocompressor is connected to the air source, and the flow rate of normal air into the turbine-compressor can be controlled by adjusting the valve. The inlet pipeline is equipped with a thermal resistance and mass flowmeter to measure the temperature and flow parameters of the air at the inlet. While, the outlet pipeline is equipped with a thermal resistance and a high-precision pressure gauge to measure the temperature and pressure of the compressed air. The compressed air is injected into the pressure tank, and under rated condition the output flow rate is 250 kg/h and the pressure is 0.15 MPa. The turbocompressor is directly connected to the shaft of the conical motor, driven by a special frequency converter, and the other end of the motor shaft is equipped with a rotor axial displacement device.     A 2.0 kW, 6000 r/min prototype machine with a cone angle of 6 • was fabricated. The whole prototype test rig, which is composed of a CR-PMSM prototype, an oscilloscope, a turbo-compressor, a volute and a controller is shown in Figure 16. Firstly, the rotor of the CR-PMSM will rotate, and due to the connection of the shaft and the turbine, the turbo-compressor starts to rotate, converting the electricity power into mechanical energy. According to the flowchart shown in Figure 14, when the motor-turbine system runs at a steady speed, the valve size is adjusted to change the load, which will cause the changes of the axial force of the turbine; at this time, the axial position of the rotor and the stator current should be controlled to change the axial electromagnetic force of the CR-PMSM, thereby suppressing the axial force of the turbine as the pressure changes.
Energies 2017, 10, x FOR PEER REVIEW 14 of 20 A 2.0 kW, 6000 r/min prototype machine with a cone angle of 6° was fabricated. The whole prototype test rig, which is composed of a CR-PMSM prototype, an oscilloscope, a turbo-compressor, a volute and a controller is shown in Figure 16. Firstly, the rotor of the CR-PMSM will rotate, and due to the connection of the shaft and the turbine, the turbo-compressor starts to rotate, converting the electricity power into mechanical energy. According to the flowchart shown in Figure 14, when the motor-turbine system runs at a steady speed, the valve size is adjusted to change the load, which will cause the changes of the axial force of the turbine; at this time, the axial position of the rotor and the stator current should be controlled to change the axial electromagnetic force of the CR-PMSM, thereby suppressing the axial force of the turbine as the pressure changes. It should be noted that the thimble of the axial displacement device is closely attached to the shaft end of the motor. Restricted by the experimental conditions, the axial length of the thimble can only be adjusted by manual manipulation of the rotary handle, to drive the motor to move axially and obtain the required axial position of the rotor. Meanwhile, the maximum axial displacement is 4.5 mm, which is measured by the displacement sensor.

Analysis of Results
The no-load back EMF is measured and depicted in Figure 17a, of which the RMS value is about 238.2 V. Compared with the 3-D FEM result, it has a quite slight difference (3%). As shown in Figure  17b, with the increase of CR axial displacement, the no-load back EMF of the motor basically decreases linearly.  It should be noted that the thimble of the axial displacement device is closely attached to the shaft end of the motor. Restricted by the experimental conditions, the axial length of the thimble can only be adjusted by manual manipulation of the rotary handle, to drive the motor to move axially and obtain the required axial position of the rotor. Meanwhile, the maximum axial displacement is 4.5 mm, which is measured by the displacement sensor.

Analysis of Results
The no-load back EMF is measured and depicted in Figure 17a, of which the RMS value is about 238.2 V. Compared with the 3-D FEM result, it has a quite slight difference (3%). As shown in Figure 17b, with the increase of CR axial displacement, the no-load back EMF of the motor basically decreases linearly.
Energies 2017, 10, x FOR PEER REVIEW 14 of 20 A 2.0 kW, 6000 r/min prototype machine with a cone angle of 6° was fabricated. The whole prototype test rig, which is composed of a CR-PMSM prototype, an oscilloscope, a turbo-compressor, a volute and a controller is shown in Figure 16. Firstly, the rotor of the CR-PMSM will rotate, and due to the connection of the shaft and the turbine, the turbo-compressor starts to rotate, converting the electricity power into mechanical energy. According to the flowchart shown in Figure 14, when the motor-turbine system runs at a steady speed, the valve size is adjusted to change the load, which will cause the changes of the axial force of the turbine; at this time, the axial position of the rotor and the stator current should be controlled to change the axial electromagnetic force of the CR-PMSM, thereby suppressing the axial force of the turbine as the pressure changes. It should be noted that the thimble of the axial displacement device is closely attached to the shaft end of the motor. Restricted by the experimental conditions, the axial length of the thimble can only be adjusted by manual manipulation of the rotary handle, to drive the motor to move axially and obtain the required axial position of the rotor. Meanwhile, the maximum axial displacement is 4.5 mm, which is measured by the displacement sensor.

Analysis of Results
The no-load back EMF is measured and depicted in Figure 17a, of which the RMS value is about 238.2 V. Compared with the 3-D FEM result, it has a quite slight difference (3%). As shown in Figure  17b, with the increase of CR axial displacement, the no-load back EMF of the motor basically decreases linearly.  The back EMF by theoretical analysis with 3-D FEM is also given in Figure 16. Because of the material parameter differences between the actual prototype and the simulation, and the influence of the manufacturing process, the calculated back EMF is larger, especially at the initial stage of the CR axial displacement. The match among the three methods is nevertheless quite good. Note that the decreasing slopes of the back EMF for all three methods are the same, this proving the exactness and validity of the analytic calculation.
During the test process, the outlet pressure should be set to adapt to 50% (half load), 100% (full load) and 150% (over load) of rated load, respectively, so the output power curve should be set and adjusted according to the 3-D FE simulation curve shown in Figure 18 to carry out the corresponding load test verification, but of course, the actual tested output power curve deviates slightly from the simulation curve, because the response to load changes is somewhat delayed, and the curve shows transient fluctuations due to disturbances caused by the load change. The back EMF by theoretical analysis with 3-D FEM is also given in Figure 16. Because of the material parameter differences between the actual prototype and the simulation, and the influence of the manufacturing process, the calculated back EMF is larger, especially at the initial stage of the CR axial displacement. The match among the three methods is nevertheless quite good. Note that the decreasing slopes of the back EMF for all three methods are the same, this proving the exactness and validity of the analytic calculation.
During the test process, the outlet pressure should be set to adapt to 50% (half load), 100% (full load) and 150% (over load) of rated load, respectively, so the output power curve should be set and adjusted according to the 3-D FE simulation curve shown in Figure 18 to carry out the corresponding load test verification, but of course, the actual tested output power curve deviates slightly from the simulation curve, because the response to load changes is somewhat delayed, and the curve shows transient fluctuations due to disturbances caused by the load change. According to Figure 12, the d-axis current has a greater influence on the axial magnetic force, while the q-axis current mainly affects the output power. In order to verify the influence of current on axial force, adjust the input current of the motor while keeping the rotor axial displacement unchanged (Δz = 0), and the d, q-axis currents and the axial magnetic force can be measured. On this basis, the current curves and axial force curve can be obtained, as shown in Figures 19 and 20, which are basically consistent with the simulation results.  According to Figure 12, the d-axis current has a greater influence on the axial magnetic force, while the q-axis current mainly affects the output power. In order to verify the influence of current on axial force, adjust the input current of the motor while keeping the rotor axial displacement unchanged (∆z = 0), and the d, q-axis currents and the axial magnetic force can be measured. On this basis, the current curves and axial force curve can be obtained, as shown in Figures 19 and 20, which are basically consistent with the simulation results. The back EMF by theoretical analysis with 3-D FEM is also given in Figure 16. Because of the material parameter differences between the actual prototype and the simulation, and the influence of the manufacturing process, the calculated back EMF is larger, especially at the initial stage of the CR axial displacement. The match among the three methods is nevertheless quite good. Note that the decreasing slopes of the back EMF for all three methods are the same, this proving the exactness and validity of the analytic calculation.
During the test process, the outlet pressure should be set to adapt to 50% (half load), 100% (full load) and 150% (over load) of rated load, respectively, so the output power curve should be set and adjusted according to the 3-D FE simulation curve shown in Figure 18 to carry out the corresponding load test verification, but of course, the actual tested output power curve deviates slightly from the simulation curve, because the response to load changes is somewhat delayed, and the curve shows transient fluctuations due to disturbances caused by the load change. According to Figure 12, the d-axis current has a greater influence on the axial magnetic force, while the q-axis current mainly affects the output power. In order to verify the influence of current on axial force, adjust the input current of the motor while keeping the rotor axial displacement unchanged (Δz = 0), and the d, q-axis currents and the axial magnetic force can be measured. On this basis, the current curves and axial force curve can be obtained, as shown in Figures 19 and 20, which are basically consistent with the simulation results.  Tests of the axial magnetic force are also carried out when both the rotor axial displacement and the input current are changed. The test results are organized into a 3-D curved surface, as shown in Figure 21b, which is similar to the simulated results shown in Figure 21a. It can be seen that when the d-axis current is enhancing the flux (id > 0), the greater the axial magnetic force, when the straight axis current is demagnetization (id < 0), the reverse is the case. The overall change trend of the torque-angular (β) characteristics of the prototype is shown in Figure 22. It can be seen that the torque decreases monotonically with the increase of the axial position of the rotor under the condition of a certain rotor angle. Compared with the 3-D FE method result, it has a quite slight difference (2%), which confirms the accuracy of developed 3-D model.  Tests of the axial magnetic force are also carried out when both the rotor axial displacement and the input current are changed. The test results are organized into a 3-D curved surface, as shown in Figure 21b, which is similar to the simulated results shown in Figure 21a. It can be seen that when the d-axis current is enhancing the flux (i d > 0), the greater the axial magnetic force, when the straight axis current is demagnetization (i d < 0), the reverse is the case. Tests of the axial magnetic force are also carried out when both the rotor axial displacement and the input current are changed. The test results are organized into a 3-D curved surface, as shown in Figure 21b, which is similar to the simulated results shown in Figure 21a. It can be seen that when the d-axis current is enhancing the flux (id > 0), the greater the axial magnetic force, when the straight axis current is demagnetization (id < 0), the reverse is the case. The overall change trend of the torque-angular (β) characteristics of the prototype is shown in Figure 22. It can be seen that the torque decreases monotonically with the increase of the axial position of the rotor under the condition of a certain rotor angle. Compared with the 3-D FE method result, it has a quite slight difference (2%), which confirms the accuracy of developed 3-D model.  The overall change trend of the torque-angular (β) characteristics of the prototype is shown in Figure 22. It can be seen that the torque decreases monotonically with the increase of the axial position of the rotor under the condition of a certain rotor angle. Compared with the 3-D FE method result, it has a quite slight difference (2%), which confirms the accuracy of developed 3-D model. Tests of the axial magnetic force are also carried out when both the rotor axial displacement and the input current are changed. The test results are organized into a 3-D curved surface, as shown in Figure 21b, which is similar to the simulated results shown in Figure 21a. It can be seen that when the d-axis current is enhancing the flux (id > 0), the greater the axial magnetic force, when the straight axis current is demagnetization (id < 0), the reverse is the case. The overall change trend of the torque-angular (β) characteristics of the prototype is shown in Figure 22. It can be seen that the torque decreases monotonically with the increase of the axial position of the rotor under the condition of a certain rotor angle. Compared with the 3-D FE method result, it has a quite slight difference (2%), which confirms the accuracy of developed 3-D model.  It should be noted that since the rotor axial displacement is manually adjusted, the continuity of the test data is not good, thus the graphics shown in Figures 21b and 22b are not smooth. The axial force of the rotor shaft with a turbine (axial force resultant) is tested, and the result curve is shown in Figure 23, where it can be seen that the axial force of the turbine is well counteracted by the axial magnetic force generated by the CR-PMSM, and the average value of total axial force is strictly constrained around zero, but unlike the ideal curve obtained by simulation, which is always zero, there are small disturbances in the test curve, which occur at those load change moments. It should be noted that since the rotor axial displacement is manually adjusted, the continuity of the test data is not good, thus the graphics shown in Figures 21b and 22b are not smooth. The axial force of the rotor shaft with a turbine (axial force resultant) is tested, and the result curve is shown in Figure 23, where it can be seen that the axial force of the turbine is well counteracted by the axial magnetic force generated by the CR-PMSM, and the average value of total axial force is strictly constrained around zero, but unlike the ideal curve obtained by simulation, which is always zero, there are small disturbances in the test curve, which occur at those load change moments. The efficiency of the CR-PMSM direct-drive turbine system is calculated by using the measured data including the input power, outlet pressure and flow rate of the turbine. Based on the turbine load characteristic, the efficiency curve varying with different load conditions can be obtained, as shown in Figure 24. A conventional PMSM direct-drive turbine system with the same power is also tested for comparison. It can be seen that the efficiency of the CR-PMSM direct-drive turbine system is about 1.5 percentage points higher at different load conditions. The main comparison parameters of the two motors are listed in Table 2. The motor constant (CA) of the CR-PMSM is basically close to but slightly lower than that of the conventional PMSM. Meanwhile, because the inner diameter of each stamping sheet of the stator core is different and the outer diameter of each stamping sheet of the rotor core is also different, the stamping dies and the laminating dies are more complicated than that of conventional machines. So, the manufacturing cost of the CR-PMSM is higher than that of conventional PMSM.
For this prototype, the manufacturing cost of the main components is listed also in Table 2, which indicates that manufacturing cost of this prototype CR-PMSM is about 16.7% higher than that of a conventional PMSM with the same power level and almost the same volume. Of course, if The efficiency of the CR-PMSM direct-drive turbine system is calculated by using the measured data including the input power, outlet pressure and flow rate of the turbine. Based on the turbine load characteristic, the efficiency curve varying with different load conditions can be obtained, as shown in Figure 24. A conventional PMSM direct-drive turbine system with the same power is also tested for comparison. It can be seen that the efficiency of the CR-PMSM direct-drive turbine system is about 1.5 percentage points higher at different load conditions. It should be noted that since the rotor axial displacement is manually adjusted, the continuity of the test data is not good, thus the graphics shown in Figures 21b and 22b are not smooth. The axial force of the rotor shaft with a turbine (axial force resultant) is tested, and the result curve is shown in Figure 23, where it can be seen that the axial force of the turbine is well counteracted by the axial magnetic force generated by the CR-PMSM, and the average value of total axial force is strictly constrained around zero, but unlike the ideal curve obtained by simulation, which is always zero, there are small disturbances in the test curve, which occur at those load change moments. The efficiency of the CR-PMSM direct-drive turbine system is calculated by using the measured data including the input power, outlet pressure and flow rate of the turbine. Based on the turbine load characteristic, the efficiency curve varying with different load conditions can be obtained, as shown in Figure 24. A conventional PMSM direct-drive turbine system with the same power is also tested for comparison. It can be seen that the efficiency of the CR-PMSM direct-drive turbine system is about 1.5 percentage points higher at different load conditions. The main comparison parameters of the two motors are listed in Table 2. The motor constant (CA) of the CR-PMSM is basically close to but slightly lower than that of the conventional PMSM. Meanwhile, because the inner diameter of each stamping sheet of the stator core is different and the outer diameter of each stamping sheet of the rotor core is also different, the stamping dies and the laminating dies are more complicated than that of conventional machines. So, the manufacturing cost of the CR-PMSM is higher than that of conventional PMSM.
For this prototype, the manufacturing cost of the main components is listed also in Table 2, which indicates that manufacturing cost of this prototype CR-PMSM is about 16.7% higher than that of a conventional PMSM with the same power level and almost the same volume. Of course, if The main comparison parameters of the two motors are listed in Table 2. The motor constant (C A ) of the CR-PMSM is basically close to but slightly lower than that of the conventional PMSM. Meanwhile, because the inner diameter of each stamping sheet of the stator core is different and the outer diameter of each stamping sheet of the rotor core is also different, the stamping dies and the laminating dies are more complicated than that of conventional machines. So, the manufacturing cost of the CR-PMSM is higher than that of conventional PMSM.
For this prototype, the manufacturing cost of the main components is listed also in Table 2, which indicates that manufacturing cost of this prototype CR-PMSM is about 16.7% higher than that of a conventional PMSM with the same power level and almost the same volume. Of course, if large-scale production is carried out, the cost of each motor will be reduced. How much reduction is related to the scale of production and the power of the motor, cannot be generalized.  Figure 25 shows the measured line voltage (U a ) and current (I a and I b ) at half load, full load and over load conditions, and the waveform presents a quite good sinusoidal. large-scale production is carried out, the cost of each motor will be reduced. How much reduction is related to the scale of production and the power of the motor, cannot be generalized.  Figure 25 shows the measured line voltage (Ua) and current (Ia and Ib) at half load, full load and over load conditions, and the waveform presents a quite good sinusoidal. The output performances at different load conditions is listed in Table 3, which shows there is a quite small difference (≤5%) between the simulated and measured results. Based on the above analysis, the axial magnetic force only has a greater increase with the increase of d-axis current, which confirms the previous analysis and conclusions in Equations (12) and (13). As shown in Figure 17 and Table 3, the air gap flux density changing ratio CB and the effective air-gap space decrease ratio CS is listed in Table 4. The output performances at different load conditions is listed in Table 3, which shows there is a quite small difference (≤5%) between the simulated and measured results. Based on the above analysis, the axial magnetic force only has a greater increase with the increase of d-axis current, which confirms the previous analysis and conclusions in Equations (12) and (13).
As shown in Figure 17 and Table 3, the air gap flux density changing ratio C B and the effective air-gap space decrease ratio C S is listed in Table 4.

Conclusions
In this paper, the CR-PMSM is proposed to realize the axial force suppression of the motor-turbine system, without balance devices. The magnetic field characteristics, no-load back EMF and axial magnetic force of the motor were analyzed. The change regulations were obtained, which indicates that the air gap magnetic field weakens with the increase of the CR axial displacement and the d-axis current demagnetization current has a greater influence on the axial magnetic force. The parameter of the axial magnetic force changing ratio F ∆z,id was defined to calculate the axial magnetic force performance. The results of theoretical calculation and experiment research show that the maximum axial displacement of the CR ∆z_ max and d-axis current are the key parameters which influence the flux weakening ratio F ∆z,id . The above conclusions can be applied to both surface-type and interior-type PMSMs. In the motor design process, an interior-type rotor and a small length to diameter ratio are appropriate properties for CR-PMSM. However, there are still some additional research work to be done and perfected in the future. It is quite challenging to derive the expression of axial force with respect to d-axis current and cone angle, as well as to analyze the efficiency in a wide speed range, which will be the priority of future work.