Next Article in Journal
Visualization Study on Thermo-Hydrodynamic Behaviors of a Flat Two-Phase Thermosyphon
Next Article in Special Issue
On Field Weakening Performance of a Brushless Direct Current Motor with Higher Winding Inductance: Why Does Design Matter?
Previous Article in Journal
A Hybrid Framework for Short Term Multi-Step Wind Speed Forecasting Based on Variational Model Decomposition and Convolutional Neural Network
Previous Article in Special Issue
Quantitative Comparisons of Six-Phase Outer-Rotor Permanent-Magnet Brushless Machines for Electric Vehicles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Node Mapping Criterion for Highly Saturated Interior PMSMs Using Magnetic Reluctance Network

by
Damian Caballero
1,2,*,
Borja Prieto
1,2,
Gurutz Artetxe
1,2,
Ibon Elosegui
1,2 and
Miguel Martinez-Iturralde
1,2
1
Ceit, Manuel Lardizabal 15, 20018 Donostia/San Sebastian, Spain
2
Universidad de Navarra, Tecnun, Manuel Lardizabal 13, 20018 Donostia/San Sebastian, Spain
*
Author to whom correspondence should be addressed.
Energies 2018, 11(9), 2294; https://doi.org/10.3390/en11092294
Submission received: 19 July 2018 / Revised: 24 August 2018 / Accepted: 27 August 2018 / Published: 31 August 2018
(This article belongs to the Special Issue Permanent Magnet Synchronous Machines)

Abstract

:
Interior Permanent Magnet Synchronous Machine (IPMSM) are high torque density machines that usually work under heavy load conditions, becoming magnetically saturated. To obtain properly their performance, this paper presents a node mapping criterion that ensure accurate results when calculating the performance of a highly saturated IPMSM via a novel magnetic reluctance network approach. For this purpose, a Magnetic Circuit Model (MCM) with variable discretization levels for the different geometrical domains is developed. The proposed MCM caters to V-shaped IPMSMs with variable magnet depth and angle between magnets. Its structure allows static and dynamic time stepping simulations to be performed by taking into account complex phenomena such as magnetic saturation, cross-coupling saturation effect and stator slotting effect. The results of the proposed model are compared to those obtained by Finite Element Method (FEM) for a number of IPMSMs obtaining excellent results. Finally, its accuracy is validated comparing the calculated performance with experimental results on a real prototype.

1. Introduction

The demand for Permanent Magnet Synchronous Machines (PMSM) is rapidly increasing in high-performance applications, such as electric vehicles [1,2,3,4], due to their high power density. In particular, the configuration of buried magnets inside the rotor is becoming very popular, because of the additional torque made available to saliency, their wide constant power speed region, and their high demagnetization withstand capability, among others [5,6].
The design process of a PMSM frequently involves the aid of software based on Finite Elements Methods (FEM) [7,8,9,10,11,12,13]. Although accuracy of the results is very high, it requires a high computational cost, together with an elevated amount of time to define the problem. This makes FEM more suitable for validation purposes rather than for preliminary machine design by iterative process. Consequently, the use of analytical design tools that rely on magnetic circuits instead of FEM is becoming widespread.
Different authors have presented simple magnetic circuits for different PMSM topologies [14,15,16]. The results are acceptable as a first estimation of the machine performance, but they are not comprehensive enough if a transient analysis or a deeper study is required. For this purpose, complex magnetic circuit models based on magnetic reluctance network, known as Magnetic Circuit Model (MCM), have been proposed [17]. The methodology is very similar to simple magnetic equivalent circuits; the main difference lies in the larger number of elements that the machine’s geometry is discretized.
In the literature, different MCMs have been proposed [18,19,20,21,22,23]. However, a clear meshing criterion that guarantees accurate results regardless of the size or geometry of the machine, i.e., a node mapping criterion for different machine regions, is not provided. In general, the reluctance element’s distribution is set according to the main flux paths [20,24]. However, in machines that work at heavy load conditions, the flux paths are in most cases unpredictable, especially in IPMSMs, where the complex geometry of the rotor makes difficult their analytical modelling [8,25,26]. In addition, the fact that the magnetic flux can only pass through an element in a unique defined direction makes it advisable to establish a general node mapping criteria for IPMSMs.
In this paper, a general node mapping criterion for IPMSMs of any geometry and size is presented. To that end, MCM with variable discretization levels for the different geometric domains is developed. The proposed MCM models V-shaped IPMSMs, with variable magnet depth and angle between magnets. A suitable MCM structure composed of generic cells, named nodal elements, as modelling element for any part of the machine is proposed. The nodal elements contain geometric and electromagnetic information of the modelled physical domain. The model contemplates rotor motion, allowing dynamic analysis, i.e., time stepping simulations. The MCM takes into account the effects of stator slotting, the airgap magnetomotive force (MMF) waveform due to armature current, and cross-magnetizing saturation effects due to an improved magnetic flux path definition. This, in turn, leads to more accurate results than other simpler magnetic circuit models shown in the bibliography.
This paper is organized as follows. First, the MCM developed for IPMSMs is described. Next, the magnetic phase flux linkage at heavy load conditions for different V-shaped IPMSMs is calculated. Then, an analysis of the accuracy of the results from using different node mappings in the MCM is conducted, and the results are compared to those obtained by the FEM software. This analysis provides a node mapping criterion to ensure sufficient accuracy when calculating the performance of an IPMSM with a magnetic reluctance network. Finally, once the node mapping criterion is described, the presented methodology is validated by comparing the calculated performance to tests on a real traction IPMSM working under heavy load conditions.

2. Definition of the Magnetic Circuit Model

In this section, the MCM developed for IPMSM is described. The MCM allows the number of nodes the machine is discretized into to be selected. The proposed MCM is divided into three regions: stator, rotor and airgap. A generic nodal element, shown in Figure 1, is proposed.
The node is set in the center of the nodal element. Each nodal element has its own coordinates (i,j), so it may be stored in a nodal elements matrix. Appropriate row and column are assigned to each nodal element according to its location on the associated discretized geometrical domain being modelled. Information regarding the associated geometry (identifier, row, column, region, position, etc.) is stored as well. Each nodal element has four sub-elements that capture both the radial and orthoradial magnetic fluxes and, therefore, possible cross-coupling effects [27,28]. The sub-element geometrical domain can be rectangular, trapezoidal or even triangular. Sub-element information that is held is fundamental for the solving process: orientation length L, transversal section A i , permeance, permeability μ , MMF source, etc. Magnetic permeance, P i , o r i e n t a t i o n , is employed for convenience in the solving process; it is the inverse of magnetic reluctance, i , o r i e n t a t i o n , which is calculated by Equation (1) [29].
i , o r i e n t a t i o n = 0 L d x μ ( x ) · A i , o r i e n t a t i o n ( x )
To reduce the computational effort, if periodicities exist within the machine, just part of it is modelled. The number of periodicities is obtained according the number of pole pairs and slots (Equation (2)).
N s i m = G C D ( Q s , p )
where Q s is the number of stator slots, p the number of pole pairs and GCD stands for Greatest Common Divisor.
The geometric positioning of the nodal elements that model different regions of the machine, is established according to the existing main magnetic flux paths. In addition, it facilitates the connection between nodal elements, and therefore, the different regions of the machine can be accurately meshed.

2.1. Stator

The stator yoke, slots, slot-openings, teeth, and tooth-tip are modelled separately. The MCM parameters that define the stator node mapping are presented in Table 1. As an illustrative example, the values given in the table result in the node mapping shown in Figure 2.
Presented parameters in Table 1 define the number of node columns and rows in stator regions, as can be observed in Figure 2. For open slot machines, N c o l , t 0 is null. The armature current MMF sources are located at the north and south sub-elements of each element belonging to the teeth and are calculated by Equation (3).
M M F w = Z s l o t , p h , s N w , p a r a l , s · N r o w , t , s · 2 · [ W s ] · [ I U V W ]
where [ W s ] is the winding sequence matrix that relates the teeth wound by each coil to the corresponding phase winding [30]. [ I U V W ] is the phase current vector. Z s l o t , p h , s is the number of turns per slot and layer. Finally, N w , p a r a l , s is the number of parallel connected winding groups.

2.2. Rotor

The parameters that define the reluctance network for IPMSM rotors are presented in Table 2. Due to the existing symmetry, only the mapping for half a pole needs to be defined. To ease the comprehension, the values given in the table result in the node mapping shown in Figure 3.
In Figure 4, the whole pole node mapping scheme for V-shaped IPMSMs is shown.
Regarding the rotor nodal mapping, some nuances must be taken into consideration:
  • N r o w , P M also define the number of node rows in the non-magnetic material and in the q-axis magnetic bridge.
  • N c o l , P M :
    -
    It defines rotor yoke’s node mapping: over and below PM side. These node distributions are “triangular”, with N c o l , P M columns and as we move away from d-axis, the number of rows decreases from N c o l , P M to one.
    -
    It defines with other parameters as N c o l , b r i d g e and N c o l , b r i d g e q and one necessary extra node (for node connection purposes), the number of columns in the N r o w , b r i d g e upper rotor nodal rows, belonging to an arc whose thickness is the bridge height.
It is remarkable that the same nodal mapping defined by Table 2 may be established for extreme cases of IPMSM. As an example, V-shaped in Figure 5 and embedded in Figure 6 are shown.
The magnets, the rotor yoke, the magnetic bridge and the rotor slot non-magnetic material, which is responsible for preventing a large magnetic flux leakage, are modelled by a selectable number of elements. The MMF contribution is computed by Equation (4) and assigned to the north and south sub-elements corresponding to permanent magnets.
M M F P M = h P M · H c
h P M is the magnet height, and H c the coercive field strength.

2.3. Airgap

This region is the most significant part of the MCM and the overall precision of the model depends on its modelling. Moreover, it is the link between the stator and rotor models, thus a correct modelling of the airgap is necessary. Therefore, the airgap model is built by taking into account the relative position between the stator and the rotor. Since machine rotation is considered, airgap node mapping is varied for each time step. An airgap element is placed at the same angular position of each stator and rotor nodal element in contact with the airgap, as shown in Figure 7. It is important to note that more finely discretized stator and/or rotor models entail a more detailed airgap region node mapping.
The number of airgap row elements is determined by the parameter N r o w , a g a p . In Figure 7, the proposed airgap node mapping is visually displayed for a N r o w , a g a p equal to 2.
To link the airgap with stator and rotor entities, auxiliary nodes are placed in the airgap region boundaries, as shown in Figure 7 with small yellow circles. They only have north and south sub-elements, and an infinite permeance. Therefore, it is a mathematical element that is introduced in the circuit matrix system. Its necessity is due to the inter-entities borders, where the airgap sub-elements are connected to the nearest auxiliary node (Figure 8). If they are connected simultaneously to various airgap nodes, the existence of these auxiliary nodes guarantees that each nodal element belonging to any machine entity is connected to other nodes by no more than four sub-elements. Thus, it is possible to generate the flux and permeanace matrix, according to the defined nodal element structure (Figure 1).

2.4. Solving Process

Once all the nodes are linked and the branches defined, the MCM is completely set. Owing to the non-linear magnetic behavior of the stator and rotor materials, whose B-H curve are imported from materials database and used as a lookup table, the MCM needs to be solved iteratively. Firstly, permeance and magnetic source flux matrices, P and ϕ , are calculated for each sub-element. The matrix circuit is solved in terms of Kirchoff’s Voltage Law (KVL), and the scalar magnetic potential at each node is obtained by:
V m = ϕ · P 1
The magnetic flux that crosses each corresponding branch from node i to node j is calculated by:
ϕ m i , j = V m i V m j + M M F i , j · P i , j
where MMF i , j is the addition of different MMF sources at each branch. A weighted average of the updated and previous iteration magnetic permeability is used in the following iteration step by the corresponding permeance sub-element, rewriting data in P and recalculating Equations (5) and (6) [20].
The iteration process is finished when the following convergence criterion is met:
μ e r r o r , k μ e r r o r , k 1 ϵ
where μ e r r o r , k is the committed mean error in sub-element permeability at iteration k.

2.5. Data Post-Processing

A key parameter for determining machine behavior is the phase magnetic flux linkage, Ψ p h [31]. The magnetic flux linked by each phase winding can be computed by Equation (8).
Ψ p h = N s i m · Z s l o t , p h , s N w , p a r a l · Φ t , s · W s
where N s i m is the number of symmetries which the problem has been divided due to the existing geometry periodicity, Z s l o t , p h , s is the number of turns per slot and layer, and N w , p a r a l is the number of pole groups in parallel. The matrix [ Φ t , s ] is the magnetic flux that crosses each pair of slot-teeth. This is easily obtained once the MCM is solved.
Depending on whether the MCM is solved for a load operating point or at no-load condition, the phase magnetic flux linkage is denoted as Ψ l o a d , p h and Ψ P M , p h , respectively. Using Park’s transformation matrix, it is possible to work at d-q rotational reference system, whose main advantage is the fact that the different variables are time invariant. The electromagnetic torque is computed by Equation (9).
T e m = 3 2 · p Ψ d · I q Ψ q · I d
The phase back EMF (Equation (10)) and the phase voltage (Equation (11)) are obtained by the time derivatives of the no-load and load fluxes:
E p h ( t ) = d Ψ P M , p h ( t ) d t
U p h ( t ) = I p h ( t ) · R p h + d Ψ l o a d , p h ( t ) d t
where I p h is the phase current and R p h is the stator phase winding resistance.

3. Results Using MCM with Different Node Mapping

The described MCM is implemented in MATLAB. To obtain a MCM node mapping criterion that guarantee an acceptable balance between results accuracy and computational costs, an analysis of the influence of the discretization level obtained by the model was conducted.
Furthermore, to assess the validity of the proposed MCM, three different size IPMSMs, whose characteristics are presented in Table 3, were evaluated.
The analysed machines have very different geometries to validate the proposed MCM and establish a node mapping criterion, regardless of the size of the machine.
The load conditions for each machine are presented in Table 4.
As Table 5 and Figure 9 illustrate, the three machines operate under heavy load conditions for the considered operating points. In the case of Motor C (Figure 9c), although the stator is not highly saturated, it can be observed that the armature current MMF is comparable to permanent magnet MMF, thereby establishing the magnetic axis almost in the q-axis. Besides, the machine comprises wide slots that behaves as barriers to the magnetic flux and, therefore, increasing the magnetic saturation at the rotor bridges. Altogether, this makes not only the magnetic flux paths more unpredictable but also the cross-coupling effect appreciable.
To establish a general node mapping criterion, the presented motor geometries have been evaluated with the aid of the developed MCM and the discretization level presented in Table 6. Each column defines a whole machine node mapping and is defined so that a specific machine’s modelled part is studied, i.e., in the case of Stator X, the node mapping parameters that define the stator are modified while the rest of parameters are kept at their lowest possible value. The parametric analysis carried out allows to separately analyse the influence of each model. The Xs in Table 6 refer to values that have been evaluated in each analysis, ranging from 1 to 6. For the parameter controlling the rotor discretization (Rotor X), the Xs values are evaluated in threes due to the complexity of the IPMSM rotor geometry and the difficulties to predict the magnetic flux paths at heavy load conditions.
As explained in Section 2.5, the key parameter to determine the machine performance is the magnetic flux linkage, Ψ l o a d , p h . Thus, it is carried out a further analysis on this parameter obtaining.
In Table 7, Ψ l o a d , p h results obtained by FEM are presented.
In Table 8, the relative errors between the results of the MCM and FEM are presented for the first and third time harmonics of the Ψ l o a d , p h waveforms. Each row corresponds to results obtained when a MCM node mapping defined in Table 6 is employed with a determined X value. For example, Stator 1 means one value for all the stator parameters, and the other parameters are set to the default value.
In Figure 10, the error tendency for the first harmonic is graphically plotted.
The discrepancies between the results obtained by the MCM and the FEM are fairly low, especially taking into account the deep saturation of the iron parts for the considered operating points.
To show graphically the great accuracy obtained by the developed MCM, some important magnetic calculated characteristics are presented and compared with obtained by FEM simulation.
In Figure 11, the Ψ l o a d , p h corresponding to Motor A is presented. Its time distribution at an electric cycle for a given stator phase is shown.
In Figure 12, the spatial distribution of airgap magnetic flux density along a pole for Motor A at no load operation point is shown. It reflects the accuracy obtained in terms of taking into account the slotting effect, noticing the well calculated magnetic flux density ripple, and also, in terms of the well calculated beginning and ending of the mentioned waveform, which means that the bridge’s model works correctly.
The MCM node mapping used is the one termed “ S t a t o r 6 ”, due to the excellent results obtained.
Employing MCM magnetic information, the main electromagnetic characteristics at load conditions have been calculated as explained in Section 2 for the three machines. These are presented in Table 9.

4. Results Discussion and Node Mapping Criterion

When a machine is modelled by a magnetic reluctance network, it is usually not clear which node mapping configuration offers the optimal balance between accuracy and simulation time. With the aim of establishing a node mapping criterion, the results collected in Table 8 and plotted in Figure 10 are examined in greater detail in the following paragraphs:
Stator (Figure 10 a)
As can be observed, the higher is the number of stator nodal elements, the better is the accuracy of the results. The main reason is that, when the number of stator elements increases, so does the number of airgap elements.
Rotor (Figure 10 b)
Although the rotor has a considerable number of nodal elements, the stator is poorly meshed. This means that the linking between airgap and stator is not good enough, especially in terms of taking the slotting effect into account.
For Motor A and Motor B, the number of airgap nodes does not increase as much when N c o l , P M is increased, unlike the case where the stator is more finely discretized. Nevertheless, for Motor C, the number of airgap nodes is significantly increased, because the number of slots per pole is low, and therefore the effect of stator discretization in the airgap meshing is lower.
Bridge (Figure 10 c)
Because of the high magnetic saturation in the bridge, the number of nodes does not affect the results.
Airgap (Figure 10 d)
The airgap meshing is controlled by the number of existing nodal elements at the stator–airgap border, and the rotor–airgap border. A different number of N r o w , a g a p , does not involve a thorough whole machine physical domain modelling. Therefore, the obtained results are invariant.
Additionally, the computational time was measured for the different mapping configurations. In Figure 13, the normalized simulation time is displayed for Motor B using different node densities.
For the analysis called BridgeX, the defined MCM has a slightly higher computational cost due to the increasing number of deeply saturated permeance sub-elements. In addition, it should be noted that the high computational cost for the RotorX analysis is due to the fact that the number of rotor elements increases three times faster than parameter X.
The reflected data in Table 10 show the required computation time using FEM and MCM corresponding to Stator1 node mapping. The time data relate to one time step solving process.
It should be highlighted that the required computation time for solving MCM shown in Table 10 is up to seven times less than the required time for solving FEM. This together with the fact that, in contrast to FEM, MCM is instantaneously generated, makes our proposed MCM a suitable and fast designing tool.
Based on the accuracy of the results and the required solving time, it can be stated that a good node mapping choice is in the range of configurations S t a t o r 3 to S t a t o r 6 .
In Table 11, the airgap node mapping densities are presented for all three machine configurations under study. The information corresponds to a single airgap row. Various airgap node row details are given: number of nodes per pole pair, average distance between two consecutive nodes, and number of nodes per slot pitch.
Next, an analysis to obtain a general node mapping criterion was carried out.
  • The number of nodes per pole pair is increased because of the higher value of X. Additionally, the difference between different geometries is due to the fact that the number of slots per pole are different, and the stator mapping is modelled for each slot-tooth pair. In this regard, setting a number of nodes per pole is not sufficient to establish a general node mapping criterion.
  • The average distance between consecutive nodes depends on the number of airgap nodes and especially on the airgap diameter, i.e., motor size. As can be seen in Motor C, despite having very low values, it does not guarantee excellent results.
  • The number of nodes per slot pitch for Motor A and Motor B with at least a value of 5 is sufficient for obtaining acceptable results. Nevertheless, when analysing rotor mappings, although this ratio is reached, the stator magnetic flux paths are not very well defined. In the case of Motor C, due to the lower number of nodes per pole pair, a value of at least 10 (Stator3) is needed. This value is also reached at rotor analysis, but, as explained previously, the stator node mapping is not very well defined.
The analysis presented here leads to the next node mapping criterion for modelling IPMSM by MCM.
  • The key point is that the number of nodes per slot in the airgap must be at least seven.
    As stated in Table 8 and the airgap region nodal information in Table 11, specifically those displayed in the third column (defined with the heading Nodes per slot and pitch), the best results are obtained for MCM node mappings with a ratio value of at least seven nodes per slot in the airgap region.
  • There must be at least three equidistant nodes in the stator slot and the tooth span.
    To consider the slotting effect and hence obtain accurate results, according to Table 8, the nodal configuration named Stator 3, that guaranties having at least three equidistant nodes in the stator slot and the tooth span, is enough.
  • At least three rows of nodal elements show be employed at both rotor yokes.
    In consonance with Figure 10, to increase the stator nodal density (Figure 10a), it is fundamental to obtain accurate results, even with the minimum number of rotor nodal rows ( N c o l , P M ), which in this case has been established as three.
By applying these criteria, an adequate connection between the stator and rotor models can be generated, leading to accurate results.

5. Experimental Validation

With the aim of validating the proposed methodology, a real traction motor with embedded magnets was tested. Its main characteristics are shown in Table 12.
The prototype machine has been previously simulated via FEM at rated load. In Figure 14, it can be appreciated that the motor is magnetically saturated, reaching more than 2 T at some points.
After applying the node mapping criteria, the main characteristics were obtained and compared with FEM simulations, as shown in Table 13. The operation point is set by defining the same current values for both numerical models.
The prototype performance was obtained by standard IEC testing in the test bench shown in Figure 15.
In Table 14, the experimental results at rated load conditions are presented.To allow for a fair comparison, the MCM was evaluated at an operating point corresponding to the rated power.
Finally, to visually compare the carried out measurements, in Figure 16, the obtained Back electromotive force waveform at tests and MCM is shown.
As it can be observed in Table 13, and especially in Table 14 and Figure 16, both experimental and MCM results match largely, showing that the presented approach is suitable for designing IPMSM and predicting its performance even under heavy saturation.

6. Conclusions

Nowadays, the use of PMSMs is exponentially increasing, with a specific interest in its use at very demanding conditions. Consequently, taking into account the different electromagnetic phenomena that take place inside the machine is crucial to predict in a very precise manner machine performance.
To this aim, in this paper, a general node mapping criterion for modelling highly saturated IPMSM using a magnetic reluctance network has been proposed. A MCM model, based on a magnetic reluctance network, is developed to model V-shaped interior mounted magnet rotors. The proposed MCM model allows selecting the discretization levels for the different machine parts. Furthermore, it not only allows simulating rotor motion but also considers the magnetic cross-coupling effect, the slotting effect and the iron magnetic saturation.
To validate the MCM model, Several V-shaped machines of different sizes and geometries were used together with FEM simulations. The results from these validations were remarkably accurate and efficient, requiring less time to complete the process. Finally, aiming to validate the MCM model with the suggested node mapping criterion, a real IPMSM prototype was tested. The comparison of the obtained results shows a great correspondence, proving the validity of the proposed method to determine highly saturated Interior PMSMs performance.
Overall, the present study expands the field of magnetic circuit models for highly saturated machines. Specifically, it provides the means to evolve within this area towards using the node mapping criterion to apply it to other highly demanded PMSM rotor solutions such as the Spoke type or the multilayer IPMSM. In addition, it would be interesting to analyse the way to optimize the implemented software code, as well as to study the use of other numerical analysis methods for solving nonlinear systems since it could lead to the development of an even faster machine designing tool.

Author Contributions

The presented work was carried out through the cooperation of all authors. D.C., B.P. and G.A. conducted the research and wrote the paper. I.E. and M.M.I. edited the manuscript and supervised the study.

Acknowledgments

This work was financially supported by the Basque Country Government Economic Development and Infraestructure Department, by means of grants program “ELKARTEK” (project KK20170095 Electromagnetic, mechanical and thermal study of light motors-MOTLIG).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
IPMSMInterior mounted Permanent Magnet Synchronous Machine
MCMMagnetic Circuit Model
FEMFinite Elements Methods
EMFElectromotive Force

References

  1. Burress, T.; Campbell, S. Benchmarking EV and HEV power electronics and electric machines. In Proceedings of the 2013 Transportation Electrification Conference and Expo (ITEC), Detroit, MI, USA, 16–19 June 2013. [Google Scholar]
  2. Sun, L.; Cheng, M.; Jia, H. Analysis of a Novel Magnetic-Geared Dual-Rotor Motor With Complementary Structure. IEEE Trans. Ind. Electron. 2015, 62, 6737–6747. [Google Scholar] [CrossRef]
  3. Zheng, P.; Wang, W.; Wang, M.; Liu, Y.; Fu, Z. Investigation of the Magnetic Circuit and Performance of Less-Rare-Earth Interior Permanent-Magnet Synchronous Machines Used for Electric Vehicles. Energies 2017, 10, 2173. [Google Scholar] [CrossRef]
  4. Gu, W.; Zhu, X.; Quan, L.; Du, Y. Design and optimization of permanent magnet brushless machines for electric vehicle applications. Energies 2015, 8, 13996–14008. [Google Scholar] [CrossRef]
  5. Wu, W.; Zhu, X.; Quan, L.; Du, Y.; Xiang, Z.; Zhu, X. Design and Analysis of a Hybrid Permanent Magnet Assisted Synchronous Reluctance Motor Considering Magnetic Saliency and PM Usage. IEEE Trans. Appl. Supercond. 2018, 28, 1–6. [Google Scholar] [CrossRef]
  6. Yue, L.; Yulong, P.; Yanjun, Y.; Yanwen, S.; Feng, C. Increasing the saliency ratio of fractional slot concentrated winding interior permanent magnet synchronous motors. Electr. Power Appl. IET 2015, 9, 439–448. [Google Scholar] [CrossRef]
  7. Bai, J.; Zheng, P.; Tong, C.; Song, Z.; Zhao, Q. Characteristic Analysis and Verification of the Magnetic-Field- Modulated Brushless Double-Rotor Machine. IEEE Trans. Ind. Electron. 2015, 62, 4023–4033. [Google Scholar] [CrossRef]
  8. Sizov, G.; Ionel, D.; Demerdash, N. Modeling and Parametric Design of Permanent-Magnet AC Machines Using Computationally Efficient Finite-Element Analysis. IEEE Trans. Ind. Electron. 2012, 59, 2403–2413. [Google Scholar] [CrossRef]
  9. Parasiliti, F.; Villani, M.; Lucidi, S.; Rinaldi, F. Finite-Element-Based Multiobjective Design Optimization Procedure of Interior Permanent Magnet Synchronous Motors for Wide Constant-Power Region Operation. IEEE Trans. Ind. Electron. 2012, 59, 2503–2514. [Google Scholar] [CrossRef]
  10. Ruuskanen, V.; Nerg, J.; Pyrhonen, J.; Ruotsalainen, S.; Kennel, R. Drive Cycle Analysis of a Permanent-Magnet Traction Motor Based on Magnetostatic Finite-Element Analysis. Veh. Technol. IEEE Trans. 2015, 64, 1249–1254. [Google Scholar] [CrossRef]
  11. Zheng, P.; Zhao, J.; Liu, R.; Tong, C.; Wu, Q. Magnetic characteristics investigation of an axial-axial flux compound-structure PMSM used for HEVs. Magn. IEEE Trans. 2010, 46, 2191–2194. [Google Scholar] [CrossRef]
  12. Chong, L.; Rahman, M. Saliency ratio derivation and optimisation for an interior permanent magnet machine with concentrated windings using finite-element analysis. IET Electr. Power Appl. 2010, 4, 249–258. [Google Scholar] [CrossRef]
  13. Cavagnino, A.; Bramerdorfer, G.; Tapia, J.A. Optimization of Electric Machine Designs—Part I. IEEE Trans. Ind. Electron. 2017, 64, 9716–9720. [Google Scholar] [CrossRef]
  14. Hwang, C.C.; Cho, Y. Effects of leakage flux on magnetic fields of interior permanent magnet synchronous motors. IEEE Trans. Magn. 2001, 37, 3021–3024. [Google Scholar] [CrossRef]
  15. Zhu, L.; Jiang, S.; Zhu, Z.; Chan, C. Analytical Modeling of Open-Circuit Air-Gap Field Distributions in Multisegment and Multilayer Interior Permanent-Magnet Machines. IEEE Trans. Magn. 2009, 45, 3121–3130. [Google Scholar] [CrossRef] [Green Version]
  16. Chen, Q.; Liu, G.; Zhao, W.; Shao, M. Nonlinear adaptive lumped parameter magnetic circuit analysis for spoke-type fault-tolerant permanent-magnet motors. Magn. IEEE Trans. 2013, 49, 5150–5157. [Google Scholar] [CrossRef]
  17. Perho, J. Reluctance Network for Analysing Induction Machines; Helsinki University of Technology: Otaniemi, Finland, 2002. [Google Scholar]
  18. Vincent, R.; Emmanuel, V.; Lauric, G.; Laurent, G. Optimal sizing of an electrical machine using a magnetic circuit model: application to a hybrid electrical vehicle. IET Electr. Syst. Transp. 2015, 6, 27–33. [Google Scholar] [CrossRef]
  19. Raminosoa, T.; Rasoanarivo, I.; Meibody-Tabar, F.; Sargos, F.M. Time-Stepping Simulation of Synchronous Reluctance Motors Using a Nonlinear Reluctance Network Method. IEEE Trans. Magn. 2008, 44, 4618–4625. [Google Scholar] [CrossRef]
  20. Tangudu, J.K.; Jahns, T.M.; El-Refaie, A.; Zhu, Z. Lumped parameter magnetic circuit model for fractional-slot concentrated-winding interior permanent magnet machines. In Proceedings of the Energy Conversion Congress and Exposition, San Jose, CA, USA, 20–24 September 2009; pp. 2423–2430. [Google Scholar]
  21. Farooq, J.; Srairi, S.; Djerdir, A.; Miraoui, A. Use of permeance network method in the demagnetization phenomenon modeling in a permanent magnet motor. Magn. IEEE Trans. 2006, 42, 1295–1298. [Google Scholar] [CrossRef]
  22. Kuttler, S.; Benkara, K.; Friedrich, G.; Vangraefschepe, F.; Abdelli, A. Analytical model taking into account the cross saturation for the optimal sizing of IPMSM. In Proceedings of the 2012 XXth International Conference on Electrical Machines (ICEM), Marseille, France, 2–5 September 2012; pp. 2779–2785. [Google Scholar]
  23. Aden, A.; Amara, Y.; Barakat, G.; Hlioui, S.; De La Barriere, O.; Gabsi, M. Modeling of a radial flux PM rotating machine using a new hybrid analytical model. In Proceedings of the 2014 International Conference on Electrical Sciences and Technologies in Maghreb (CISTEM), Tunis, Tunisia, 3–6 November 2014; pp. 1–5. [Google Scholar]
  24. Amrhein, M.; Krein, P. Induction Machine Modeling Approach Based on 3-D Magnetic Equivalent Circuit Framework. IEEE Trans. Energy Convers. 2010, 25, 339–347. [Google Scholar] [CrossRef]
  25. Rasmussen, C.; Ritchie, E. A magnetic equivalent circuit approach for predicting PM motor performance. In Proceedings of the Conference Record of the 1997 IEEE Industry Applications Conference Thirty-Second IAS Annual Meeting, New Orleans, LA, USA, 5–9 October 1997; pp. 10–17. [Google Scholar]
  26. Lovelace, E.; Jahns, T.; Lang, J.H. A saturating lumped-parameter model for an interior PM synchronous machine. IEEE Trans. Ind. Appl. 2002, 38, 645–650. [Google Scholar] [CrossRef] [Green Version]
  27. Yamazaki, K.; Kumagai, M. Torque Analysis of Interior Permanent-Magnet Synchronous Motors by Considering Cross-Magnetization: Variation in Torque Components With Permanent-Magnet Configurations. IEEE Trans. Ind. Electron. 2014, 61, 3192–3201. [Google Scholar] [CrossRef]
  28. Lee, S.; Jeong, Y.S.; Kim, Y.J.; Jung, S.Y. Novel Analysis and Design Methodology of Interior Permanent-Magnet Synchronous Motor Using Newly Adopted Synthetic Flux Linkage. IEEE Trans. Ind. Electron. 2011, 58, 3806–3814. [Google Scholar]
  29. Ostovic, V. Dynamics of Saturated Electric Machines; Springer: New York, NY, USA, 1989. [Google Scholar]
  30. Han, S.H.; Jahns, T.; Soong, W. A Magnetic Circuit Model for an IPM Synchronous Machine Incorporating Moving Airgap and Cross-Coupled Saturation Effects. In Proceedings of the Electric Machines & Drives 2007 International Conference, Antalya, Turkey, 3–5 May 2007; pp. 21–26. [Google Scholar]
  31. Ong, C.M. Dynamic Simulation of Electric Machinery: Using MATLAB/SIMULINK; Upper Saddle River: Bergen, NJ, USA, 1998. [Google Scholar]
Figure 1. Nodal element structure.
Figure 1. Nodal element structure.
Energies 11 02294 g001
Figure 2. Stator nodal model.
Figure 2. Stator nodal model.
Energies 11 02294 g002
Figure 3. Half pole rotor nodal model for V-shaped PMSMs.
Figure 3. Half pole rotor nodal model for V-shaped PMSMs.
Energies 11 02294 g003
Figure 4. Whole pole rotor nodal model for V-shaped PMSMs.
Figure 4. Whole pole rotor nodal model for V-shaped PMSMs.
Energies 11 02294 g004
Figure 5. Rotor nodal model for V-shaped PMSMs.
Figure 5. Rotor nodal model for V-shaped PMSMs.
Energies 11 02294 g005
Figure 6. Rotor nodal model for embedded PMSMs.
Figure 6. Rotor nodal model for embedded PMSMs.
Energies 11 02294 g006
Figure 7. Airgap nodal model.
Figure 7. Airgap nodal model.
Energies 11 02294 g007
Figure 8. Airgap multinode linking.
Figure 8. Airgap multinode linking.
Energies 11 02294 g008
Figure 9. Magnetic flux density distribution for the three motors.
Figure 9. Magnetic flux density distribution for the three motors.
Energies 11 02294 g009
Figure 10. MCM node mapping Ψ l o a d , p h relative error.
Figure 10. MCM node mapping Ψ l o a d , p h relative error.
Energies 11 02294 g010
Figure 11. Phase magnetic flux linkage (Motor A).
Figure 11. Phase magnetic flux linkage (Motor A).
Energies 11 02294 g011
Figure 12. Airgap magnetic flux density at no load condition (Motor A).
Figure 12. Airgap magnetic flux density at no load condition (Motor A).
Energies 11 02294 g012
Figure 13. Normalized time simulation vs MCM node map.
Figure 13. Normalized time simulation vs MCM node map.
Energies 11 02294 g013
Figure 14. Prototype’s Magnetic flux density distribution.
Figure 14. Prototype’s Magnetic flux density distribution.
Energies 11 02294 g014
Figure 15. IPMSM Prototype at test bench.
Figure 15. IPMSM Prototype at test bench.
Energies 11 02294 g015
Figure 16. Protype’s Back EMF at 50% rated speed.
Figure 16. Protype’s Back EMF at 50% rated speed.
Energies 11 02294 g016
Table 1. Stator reluctance network parameters.
Table 1. Stator reluctance network parameters.
DefinitionParameterValue (See Figure 2)
Yoke rows N r o w , y , s 1
Tooth rows N r o w , t , s 4
Tooth-tip rows N r o w , s 0 1
Tooth columns N c o l , t , s 3
Tooth-tip columns (One side) N c o l , t 0 2
Slot-opening columns N c o l , s 0 2
Table 2. IPMSM rotor reluctance network parameters.
Table 2. IPMSM rotor reluctance network parameters.
DefinitionParameterValue (See Figure 3)
PM rows N r o w , P M 1
Bridge rows N r o w , b r i d g e 1
Bridge columns N c o l , b r i d g e 1
Bridge inter pole columns N c o l , b r i d g e , q 1
PM columns N c o l , P M 3
Table 3. V-shaped IPMSM characteristics.
Table 3. V-shaped IPMSM characteristics.
ParameterMotor AMotor BMotor C
Number of poles468
Number of slots485424
Airgap length [mm]221
Stator outer diameter [mm]820370150
Rotor outer diameter [mm]570276108
Stack length [mm]820404300
V-shaped PM depth [mm]1002410
V-shaped PM height [mm]953.5
Stator yoke height [mm]642810
Stator slot width [mm]156.56
Magnet remanence [T]1.3041.3041.13
Magnet relative permeability1.061.061.04
Table 4. Load operation point.
Table 4. Load operation point.
ParameterMotor AMotor BMotor C
Output power [kW]14718050
Speed [r/min]70012502500
Phase current RMS [A]160400130
q-component current [A]148392125
d-component current [A]−60−79−34
Table 5. FEM measured B m a x .
Table 5. FEM measured B m a x .
RegionMotor AMotor BMotor C
Stator Tooth [T]1.82.01.5
Stator Core [T]1.91.91.6
Rotor Bridge [T]2.32.52.6
Rotor Core [T]1.61.51.4
Table 6. MCM node mapping employed in the analysis.
Table 6. MCM node mapping employed in the analysis.
ParameterStatorXAirgapXBridgeXRotorX
N r o w , s 0 X111
N c o l , s 0 X111
N r o w , t s X111
N c o l , t s X111
N r o w , y s X111
N r o w , b r i d g e 11X1
N r o w , P M 2222
N c o l , P M 333 3 · X
N c o l , b r i d g e 11X1
N c o l , b r i d g e q 11X1
N r o w , a g a p 2X22
Table 7. Ψ l o a d , p h FEM results.
Table 7. Ψ l o a d , p h FEM results.
1st Harmonic3rd Harmonic
Motor A2988 mWb250 mWb
Motor B702 mWb29.4 mWb
Motor C246 mWb14 mWb
Table 8. MCM node mapping Ψ l o a d , p h relative error values.
Table 8. MCM node mapping Ψ l o a d , p h relative error values.
1st Harmonic [%]3rd Harmonic [%]
MotorABCABC
Stator15.39.321.79.112.215.1
Stator22.53.013.57.712.911.8
Stator31.00.84.22.615.37.4
Stator41.10.43.20.813.78.4
Stator50.50.51.02.115.26.3
Stator60.60.40.40.413.87.2
Rotor35.49.321.89.112.615.0
Rotor65.27.018.87.225.514.2
Rotor95.38.614.51.315.89.5
Rotor125.29.012.310.022.98.1
Rotor154.98.811.012.327.212.6
Rotor184.69.410.113.928.16.1
Bridge15.49.321.89.212.615.0
Bridge24.86.222.516.713.515.0
Bridge34.68.023.012.37.215.1
Airgap15.49.421.89.112.014.7
Airgap25.39.321.89.212.215.0
Airgap35.39.321.79.312.215.0
Table 9. Load performance.
Table 9. Load performance.
ValueRelative Error [%]
MotorABCABC
Torque [N· m]200513701730.2−0.5−4.8
U p h [V]1st harm434.0285.8285.2−0.73.66.3
3rd harm103.138.942.8−7.15.0−2.7
Table 10. Computational time MCM vs FEM.
Table 10. Computational time MCM vs FEM.
Number of ElementsComputation Time [s]
MCMFEMMCMFEM
Motor A424531541.39.1
Motor A364211491.23.6
Motor A208250671.45.1
Table 11. Airgap region nodal information.
Table 11. Airgap region nodal information.
Node MapNodes per Pole PairAverage Nodes Distance [mm]Nodes per Slot Pitch
ABCABCABC
Stator172603612.44.82.433.36
Stator212096487.53.01.855.38
Stator3168132605.32.21.477.310
Stator4216168724.21.71.299.312
Stator5264204843.41.41.01111.314
Stator6312240962.91.20.91313.316
Rotor372603612.44.82.433.36
Rotor684724810.64.01.83.548
Rotor99684609.33.51.444.610
Rotor1210896728.33.01.24.55.312
Rotor15120108847.52.71.05614
Rotor18132120966.82.40.95.56.616
Table 12. Tested IPMSM characteristics.
Table 12. Tested IPMSM characteristics.
ParameterValue
Output power [kW]225
Voltage [V]318
Stator outer diameter [mm]420
Rotor outer diameter [mm]317
Stack length [mm]380
Magnet qualityN40UH
Table 13. Prototype’s performance calculated via proposed MCM and FEM.
Table 13. Prototype’s performance calculated via proposed MCM and FEM.
ParameterMCMFEMRelative Error [%]
Back EMF [V]3293154.4
Electromagnetic Power [kW]2302242.7
Phase current [A]4344340.0
Power factor0.990.990.9
Line Voltage [V]3163053.3
D-axis load Magnetic flux linkage [mWb]3163005.4
Q-axis load Magnetic flux linkage [mWb]2712730.8
Stator Tooth load Magnetic flux density [T]2.22.16.3
Stator Core load Magnetic flux density [T]1.41.42.9
Rotor Bridge load Magnetic flux density [T]2.52.810.7
Rotor Core load Magnetic flux density [T]2.02.13.4
Table 14. Prototype measured performance at test bench.
Table 14. Prototype measured performance at test bench.
ParameterMCMTestRelative Error [%]
Back EMF [V]3293028.9
Mechanical Power [kW]2252230.9
Phase current [A]4344252.1
Power factor0.990.990.8
Line Voltage [V]3183180.0

Share and Cite

MDPI and ACS Style

Caballero, D.; Prieto, B.; Artetxe, G.; Elosegui, I.; Martinez-Iturralde, M. Node Mapping Criterion for Highly Saturated Interior PMSMs Using Magnetic Reluctance Network. Energies 2018, 11, 2294. https://doi.org/10.3390/en11092294

AMA Style

Caballero D, Prieto B, Artetxe G, Elosegui I, Martinez-Iturralde M. Node Mapping Criterion for Highly Saturated Interior PMSMs Using Magnetic Reluctance Network. Energies. 2018; 11(9):2294. https://doi.org/10.3390/en11092294

Chicago/Turabian Style

Caballero, Damian, Borja Prieto, Gurutz Artetxe, Ibon Elosegui, and Miguel Martinez-Iturralde. 2018. "Node Mapping Criterion for Highly Saturated Interior PMSMs Using Magnetic Reluctance Network" Energies 11, no. 9: 2294. https://doi.org/10.3390/en11092294

APA Style

Caballero, D., Prieto, B., Artetxe, G., Elosegui, I., & Martinez-Iturralde, M. (2018). Node Mapping Criterion for Highly Saturated Interior PMSMs Using Magnetic Reluctance Network. Energies, 11(9), 2294. https://doi.org/10.3390/en11092294

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop