Next Article in Journal
Effect of the Dendrimer Generation Used in the Synthesis of Pt-Ru Nanoparticles Supported on Carbon Nanofibers on the Catalytic Activity towards Methanol Oxidation
Previous Article in Journal
The Impact of Shale Gas on the Cost and Feasibility of Meeting Climate Targets—A Global Energy System Model Analysis and an Exploration of Uncertainties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Torque Estimation for an IPMSM with Cross-Coupling and Parameter Variations

1
Center, Dongwoon Anatech Co. Ltd., Seoul 05029, Korea
2
Department of Electrical Engineering, Konkuk University, Seoul 05029, Korea
3
Department of Energy Mechanical Engineering, Gyeonggi College of Science and Technology, Gyeonggi 15073, Korea
4
Department of Creative IT Engineering, Phohang University of Science and Technology, Gyeongbuk 37673, Korea
*
Author to whom correspondence should be addressed.
Energies 2017, 10(2), 167; https://doi.org/10.3390/en10020167
Submission received: 30 September 2016 / Accepted: 23 January 2017 / Published: 27 January 2017

Abstract

:
This paper presents a new adaptive torque estimation algorithm for an interior permanent magnet synchronous motor (IPMSM) with parameter variations and cross-coupling between d- and q-axis dynamics. All cross-coupled, time-varying, or uncertain terms that are not part of the nominal flux equations are included in two equivalent mutual inductances, which are described using the equivalent d- and q-axis back electromotive forces (EMFs). The proposed algorithm estimates the equivalent d- and q-axis back EMFs in a recursive and stability-guaranteed manner, in order to compute the equivalent mutual inductances between the d- and q-axes. Then, it provides a more accurate and adaptive torque equation by adding the correction terms obtained from the computed equivalent mutual inductances. Simulations and experiments demonstrate that torque estimation errors are remarkably reduced by capturing and compensating for the inherent cross-coupling effects and parameter variations adaptively, using the proposed algorithm.

1. Introduction

In recent years, interior permanent-magnet synchronous motors (IPMSMs) have emerged as superior solutions for use in high-performance drive applications, owing to their advantageous features such as high efficiency, high power density, robust and wide-speed operation, and low noise [1,2,3]. The benefits of using the IPMSM drives have been established well in [4,5].
In order to make good use of the aforementioned advantageous features of an IPMSM, and to achieve fast and accurate responses, it is very important to obtain the exact d- and q-axis reactance parameters of the IPMSM and to use them to estimate the applied torque with high accuracy [6]. Such accurate torque estimation is essential to produce the desired torque. In the electric-vehicle industry, it is believed that the torque tracking errors should be less than 5% in order to provide torque with tolerable precision with respect to the accelerator pedal position [7]. Widely used ideal models based on the orthogonal or decoupled d- and q-axes may be different from an actual system with cross-coupling effects and various parameter variations, as the operation of the IPMSM is strongly affected by the rotor magnetic saliency [8], core saturation [9], temperature [10], and armature reaction effects [11]. For this reason, the torque equation derived from an ideal model does not hold when applied in real cases [12,13,14]. For accurate torque estimation, a more sophisticated model is required, which should account for the cross-coupling and parameter variations.
Two representative approaches have been used to overcome the aforementioned inaccurate torque estimation issues arising from the inherent physical limitations: a lookup table method and a direct torque control (DTC) scheme. Lookup tables, as widely used methods for estimating the torque of an IPMSM, have been employed to tabulate the optimal current commands computed a priori with respect to various operation conditions such as the DC link voltage, rotor speed, magnet temperature, etc. [15,16]. The current commands corresponding to the desired torque references are drawn from these tables based on the operation conditions. Constructing such lookup tables requires time and effort, as the same experiments should be repeated under different conditions. Furthermore, minor parameter variations from aging effects may not be reflected immediately on an already-determined lookup table. In order to reduce the size of the tables, analytical approaches have been trying to solve the torque equations with some voltage constraints, at the expense of real-time computation [1,17]. However, since these analytical approaches are strongly dependent on mathematical models, they tend to be less robust and more sensitive to uncertainties in situations of unexpected phenomena, which are not reflected in the models. The DTC scheme for IPMSM drives, which is another approach widely used to overcome the inaccurate torque estimation, has been explored since the 1990s [18,19]. DTC does not require an accurate model or parameters, except for the armature resistances. However, DTC exhibits problems such as unfixed switching frequency and large torque and stator flux ripples due to less dependence on system parameters. These drawbacks may result in reduced system efficiency and a short life span [20]. In addition to the lookup table method and the DTC approach, the estimation of some parameters of IPMSMs have been attempted by using flux observers [21,22], by setting other parameters to their nominal [23,24] or measured [25,26] values, or by injecting perturbation signals [27,28,29] and DC offset voltages [30]. However, these existing estimation schemes often demonstrate mismatches between the nominal parameter values and the actual ones, as nominal models are used for estimation or unmodelled uncertainties are not considered. Especially, machine parameters can vary with the injected signals; this is difficult to describe in a mathematical fashion, and, as a result, often remains unmodelled. To the best of the authors’ knowledge, there has been no report of a real-time adaptive torque estimation algorithm that deals with a variety of uncertainties including cross-coupling effects and parameter variations, and that provides a more accurate torque equation. It would be very useful to develop such an adaptive torque estimation algorithm for IPMSMs.
In this paper, we present a new adaptive torque estimation algorithm for an IPMSM with cross-coupling and time-varying d- and q-axis reactance parameters. It is observed through experiments that the inductances are influenced by the magnitudes of both the corresponding current and the orthogonal current, and have their own nonlinearities and cross-coupling effects between the d- and q-axis components. Such nonlinear interactive inductances motivated us to introduce equivalent mutual inductances. All cross-coupling, time-varying, or uncertain terms that are not present in the nominal flux equations are placed into two equivalent mutual inductances described using the so-called equivalent d- and q-axis back electromotive forces (EMFs). The proposed algorithm estimates the equivalent d- and q-axis back EMFs in a recursive and stability-guaranteed way, in order to compute the equivalent mutual inductances between the d- and q-axes, and provides a more accurate and adaptive torque equation by adding the correction terms obtained from the computed equivalent mutual inductances. It is shown through simulations and experiments that the proposed torque estimation algorithm greatly reduces the torque estimation errors by capturing and compensating for the cross-coupling effects and parameter variations due to the changes in operating conditions.
The rest of this paper is organized as follows: Section 2 describes a torque equation derived from a more elaborate model of an IPMSM. In Section 3, the torque estimation algorithm is proposed to obtain the equivalent mutual inductance and finally produce an accurate torque equation. The performance of the proposed torque estimation algorithm is illustrated through simulations and experiments in Section 4 and Section 5, respectively. Finally, the proposed work is summarized and the conclusions are presented in Section 6.

2. Real IPMSM Model and Its Corresponding Torque Equation

A precise model of the IPMSM is very important for high-performance vector-controlled servo drive systems. However, for simplicity, the IPMSM is often represented in terms of an ideal model based on the following assumptions:
  • The three-phase stator windings are symmetrical.
  • The slot effect is neglected and the back EMF is exactly sinusoidal.
  • Saturation, eddy currents, and hysteresis losses are neglected.
  • No damping exists in the rotor.
With these assumptions, the flux equations of the IPMSM in the synchronous reference frame are given by
λ d = L d i d + λ m ,
λ q = L q i q ,
where i d and i q are the d- and q-axis stator currents, L d and L q are the d- and q-axis stator inductances, λ d and λ q are the d- and q-axis stator magnetic fluxes, and λ m is the rotor flux. The electrical torque T e based on the ideal model (1) and (2) is expressed as
T e = P 2 3 2 ( λ d i q λ q i d ) = P 2 3 2 ( λ m i q + ( L d L q ) i d i q ) ,
where P is the number of poles.
As seen in Equations (1) and (2), the ideal model assumes that the d- and q-axis stator inductances L d and L q are constant, and hence, independent of any currents, which is not the case in reality. Considering the core saturation effect and the spatial unbalance of fluxes, the d- and q-axis stator inductances are the functions of the d- and q-axis currents, and, furthermore, the cross-coupling phenomena cannot be neglected. Therefore, such current-dependent inductances and the mutual interactions between the d- and q-axes should be reflected for accurate torque estimation. In other words, it is necessary to modify the IPMSM model in Equations (1) and (2) and the corresponding torque Equation (3) into more realistic versions.
The coupling effect is illustrated through an experiment, which is shown in Figure 1. The dq axis output voltages were obtained by changing the dq axis currents, while keeping the motor speed constant at the rated speed. The value of the flux for the d–q current was calculated from the motor voltages and currents in the steady state. It is observed that the d-axis flux is affected by the q-axis current and the q-axis flux by the d-axis current. The ideal model (1) and (2) does not reflect such coupling. By using inductances varied by the d- and q-axis currents, and adding supplementary terms for reflecting the coupling effects, the ideal IPMSM model (1) and (2) can be rewritten as
λ d = L d ( i d , i q ) i d + λ m + L d q i q ,
λ q = L q ( i d , i q ) i q + L q d i d ,
where L d q and L q d are the mutual inductances of the d- and q-axis windings with respect to the q- and d-axis currents, respectively. Furthermore, considering the uncertainties of the inductances and the rotor flux yields:
λ d = ( L d + L d 0 ) i d + ( λ m + λ m 0 ) + ( L d q + L d q 0 ) i q ,
λ q = ( L q + L q 0 ) i q + ( L q d + L q d 0 ) i d ,
where L d 0 , L q 0 , λ m 0 , L d q 0 , and L q d 0 are constant nominal values and their variations are denoted by L d , L q , λ m , L d q , and L q d . Rearranging the terms in Equations (6) and (7), we obtain
λ d = L d 0 i d + λ m 0 + L ε d i q ,
λ q = L q 0 i q + L ε q i d ,
where L ε d and L ε q are given by
L ε d = L d q + L d q 0 + L d i d + λ m i q ,
L ε q = L q d + L q d 0 + L q 0 i q i d ,
and are called equivalent mutual inductances, in this paper. It is noted that the two mutual inductances— L ε d and L ε q —are the functions of the d- and q-axis currents and inductance uncertainties. For simple notations, L ε d and L ε q will be used without arguments throughout this paper. Considering the mutual inductances and uncertainties, or substituting the fluxes (8) and (9) into the torque Equation (3), we have
T e = P 2 3 2 ( λ m 0 i q + ( L d 0 L q 0 ) i d i q L ε q i d 2 + L ε d i q 2 ) .
It is observed from the comparison with Equation (3) that the torque components arising from the coupling effect and inductance variations are added in the third and fourth terms of Equation (12). Therefore, for more accurate torque estimation, the third and fourth terms in Equation (12) should be considered. In the following section, the estimation of L ε d and L ε q will be discussed. Figure 2 shows the efficiency map of the IPMSM at various operating points. It is calculated based on the above assumptions in which the core loss is not taken into account.

3. Algorithm for Estimating Mutual Inductances

In order to estimate L ε d and L ε q in Equations (10) and (11), we shall exploit the following voltage equations in the rotating dq frame:
V d = R i d + d λ d d t ω r λ q ,
V q = R i q + d λ q d t + ω r λ d ,
where V d and V q are the stator voltages, R is the stator resistance, λ d and λ q are the d- and q-axis stator magnetic fluxes, and ω r is the electrical rotor speed. Substituting Equations (8) and (9) into Equations (13) and (14), and using the matrix-vector form, we can write Equations (13) and (14) in a compact manner as follows:
[ V d V q ]   = [ R + L d 0 ρ ω r L ε q L ε d ρ ω r L q 0 L ε q ρ + ω r L d 0 R + L q 0 ρ + ω r L ε d ] [ i d i q ] + [ 0 ω r λ m 0 ] ,
where ρ denotes a differential operator. Representing the voltage Equations (15) in a simple form without the dq axis coupling effects, we have
[ V d V q ]   = [ R + L d 0 ρ ω r L q 0 ω r L d 0 R + L q 0 ρ ] [ i d i q ] + [ E x d E x q ] ,
where E x d and E x q are given by
E x d = ω r L ε q i d + L ε d ρ i q ,
E x q = ω r λ m 0 + ω r L ε d i q + L ε q ρ i d ,
and are called the equivalent d- and q-axis back EMFs.
It is noted that all nonlinear, time-varying, or uncertain terms not present in the nominal stator voltage equations are included into the equivalent back EMFs E x d and E x q in Equations (17)–(18). The voltage equations in (16) have the same form as the nominal stator voltage Equations (13) and (14) only if the equivalent back EMFs are replaced with [ 0 ω r λ m 0 ] T , and the mutual inductances L ε d and L ε q are set to zero.
From Equations (17) and (18), the mutual inductances can be obtained as
L ε d = E x d ρ i d + ω r i d E x q ω r 2 i d λ m 0 ω r 2 i d i q + ρ i d ρ i q ,
L ε q = ω r i q E x d E x q ρ i q + ω r λ m 0 ρ i q ω r 2 i d i q + ρ i d ρ i q ,
which means that L ε d and L ε q are available only if E x d and E x q are known.
By setting the estimate of E x d as E ^ x d , and combining Equations (8), (13), and (17), we can design an estimator as follows:
i d = ( V d + ω r L q 0 i q E x d ) 1 L d 0 s + R ,
i ^ d = ( V d + ω r L q 0 i q E ^ x d ) 1 L d 0 s + R ,
E ^ x d = ( i d i ^ d ) ( K P d + K I d s ) ,
where K P d and K I d are constant gains and i ^ d can be considered as an estimate of the real current i d . The estimators (21)–(23) have the physical meaning that E ^ x d can be obtained by reducing the difference between a real current i d and its estimate i ^ d . The mathematical expression (23) is similar to that of a proportional-integral (PI) control that makes i ^ d follow i d well. The estimators (21)–(23) are depicted in Figure 2, together with a typical PI plus forward current control. By substituting Equations (8) and (9) into Equations (13) and (14), we can observe that the forward current controls V d _ f f * and V q _ f f * in Figure 2 can be represented as
V d _ f f * = ω r λ q = ω r L q i q
and
V q _ f f * = ω r λ d = ω r ( L d i d + λ m ) ,
respectively. V d _ f b * and V q _ f b * denote the voltages commanded through feedback current controls.
In order to guarantee the convergence of E x d , we obtain the transfer function from E x d to E ^ x d .
Combining Equations (21)–(23) yields
E ^ x d E x d = K P d s + K I d L s 2 + ( R + K P d ) s + K I d .
In order to obtain a stable estimator, or to move all poles to the left half-plane, K P d > R   and K I d > 0   should be met. K P d can be chosen to be positive in order to avoid the non-minimum phase. It is noted that a unit DC gain is guaranteed in Equation (26) since the following relation holds:
l i m s 0 s E ^ x d E x d 1 s = 1 .
In a similar way, to estimate E x d , E x q can be estimated from i q and i ^ q . The bandwidth of the current controller and back EMF estimator in Figure 3 is selected as ω c c = 3600   [ r a d / s ] . By the pole–zero cancellation method of the open-loop transfer function, the gains of the d-axis current controller are determined as K P d s = L q 0 ω c c ,   K I d s = R ω c c . The gains of the q-axis current controller are also determined as K P q s = L d 0 ω c c ,   K I q s = R ω c c . The gains of the back EMF estimators are selected as K P d = L d 0 / R ,   K I d s = ω c c , K P q = L q 0 / R ,   K I q s = ω c c . The overall view of a control system including the estimators E ^ x d and E ^ x q is depicted in Figure 4. The equivalent mutual inductances in Equations (19) and (20) are obtained from the estimated equivalent back EMFs E ^ x d and E ^ x q . In Figure 4, the gains of the flux weakening controller, K P f w and K I f w are chosen as small as possible by the trial-and-error method because it is difficult to obtain an accurate model.
The proposed torque estimation algorithm uses the back EMF and the current value estimated using the PI state filter as shown in Figure 3. If there are motor resistance error ( R ) and voltage error ( Δ V q , Δ V d ) due to dead time effects between the voltage reference and the inverter output voltage, the estimated torque error can be expressed by Equation (28):
Δ T = - Δ R ( i d 2 ω r + i q 2 ω r ) ( i d 2 Δ V q ω r i q + i q 2 Δ V d ω r i d ) ( i d 2 Δ V q ω r i q + i q 2 Δ V d ω r i d ) .
In the high speed operation, the dead-time effect is not so serious because the relative voltage error becomes small. The estimation error of the torque due to the voltage error is insignificant. In the low speed region, the voltage error due to dead time effect is relatively large and the speed is minimal. Therefore, the torque estimation error has a larger value than its value in the high speed region. In this paper, the dead time is set to 3 μs. To reduce the voltage error, the error caused by the inverter is identified and is used in the voltage compensator, where the voltage error is cancelled [31].

4. Simulation

In order to illustrate the performance of the proposed torque estimation algorithm, simulations are carried out by perturbing inductances and fluxes. The accuracy of the torque estimates computed from the proposed torque equation is evaluated, and is compared with that from the existing equation.
The parameter and values for the IPMSM model used in the simulations are shown in Table 1; they are obtained from the real IPMSM employed in the experiments in this paper. The sampling frequency for control action is set at 10 kHz. A q-axis current of 130 A is applied and the corresponding d-axis current is determined in order to achieve maximum torque per ampere (MTPA) [2]:
i q = λ m 0 2 ( L q 0 L d 0 ) λ m 0 2 4 ( L q 0 L d 0 ) 2 + i q 2 .
It is noted that L d 0 and L q 0 are the nominal inductances in Equations (6) and (7). For implementation, care should be taken for the zero current in the steady state. Since the mutual inductances in Equations (19) and (20) reduce to
L ε d = E x q ω r λ m 0 ω r i q
and
L ε q = E x q ω r i d
in the steady state, the zero current can make the mutual inductances diverge. The proposed torque estimation scheme works well only if nonzero currents are applied.
For performance comparison, the percentage torque estimation error is computed as follows:
T m T e T m × 100 ,
where T m is the measured torque and T e is the torque estimated using Equations (3) or (12). The percentage torque estimation error is obtained by sweeping the inductances and the fluxes from 55% to 145% of the nominal values.
A more accurate and complex IPMSM model is employed for the simulation, which takes into consideration the magnetic saturation and cross-coupling effects [32]. Specially, the fluxes are computed from:
λ d ( i d , i q ) = K L d ( i d + I 0 ) 1 + K S d | i d + I 0 | + K S d q | i q | + λ 0 ,
λ q ( i d , i q ) = K L q i q 1 + K S q d | i d + I 0 | + K S q | i q | ,
where I 0 and λ 0 are the initial current and flux, respectively, and K L d , K L q , K S d , K S q , K S d q , K S q d , I 0 , and λ 0 are determined to be 0.000385987, 0.0003585, 0.00208, 0.00154, 0.005, 0.001298, 40, and 0.03363, respectively. All parameter values in Equations (33) and (34) are chosen to minimize the difference between the closed-form Expressions (33) and (34) and the experimental data obtained from Figure 1. To this end, a nonlinear least-square method is employed in this paper.
Figure 5 shows the flux computed from Equations (33) and (34), together with the one estimated from Equations (8) and (9). The current is applied from 0.15 s. As mentioned earlier, the mutual inductances are not well-defined, or diverging, for zero currents. Therefore, the dq fluxes are chosen to be λ m and 0 at the initial time when no current is applied. These choices make sense since the d and q fluxes corresponding to zero current are computed to be λ m and 0, from an ideal flux model (1) and (2). It is observed that significant flux estimation errors are present from 0 s to 0. 15 s. After 0.15 s, the estimated flux from the proposed algorithm follows the real one from a more realistic model (33) and (34) well. It can be seen that the estimated flux from an ideal model has more biased and larger estimation errors than that from the proposed algorithm. In particular, more significant performance improvements are obtained in estimating the d-axis flux.
Figure 6 shows the percentage torque estimation errors due to the inductance and flux variations. L d 0 , L q 0 , and λ m 0 are swept from 55% to 145% of the nominal values in steps of 15%. The conventional torque equation provides estimation errors between −18% and 3% due to L d 0 and L q 0 variations. The percentage torque estimation error due to the λ m 0 variation is more significant; it varies from −49% to 40%. The proposed torque equation provides a small estimation error between −2% and 0.9% due to   L d 0 , L q 0 , and λ m 0 variations. Even though the parameter values are quite different from the nominal ones, the estimation errors obtained from the proposed torque equation are still around 0%. The torque estimation scheme based on the proposed torque equation can be said to be adaptive and robust to various model uncertainties. Compared with the existing torque Equation (3), the proposed one (12) demonstrates superior performance while estimating the torque.

5. Experiment

Figure 7 depicts the experimental set up. All parameters of the IPMSM used for the experiment are the same as the ones used for simulation. A torque sensor is employed to measure the real torque and compare it with the estimated one computed from the torque equation. As in the simulation, we sweep   L d 0 , L q 0 , and λ m 0 from 55% to 145% of the nominal values in steps of 15%. For each inductance and flux, the experiment is carried out as shown in Figure 8. ①, ②, and ③ correspond to the torque values based on the conventional torque equation, the proposed one, and the torque sensor. The experiment is repeated for all parameter variations.
Figure 9 shows the percentage torque estimation errors due to the inductance and flux variations. The existing torque equation provides estimation errors ranging from −4% to 21%, −7% to 25.9%, and −28% to 46% due to the L d 0 , L q 0 , and λ m 0 variations, respectively. As in the case of the simulation, the torque estimation error due to the λ m 0 variation is the most sensitive. The proposed torque equation provides a small estimation error of −0.3% to 0.7% due to the L d 0 , L q 0 , and λ m 0 variations, which is consistent with the simulation results. Overall, the experiment illustrates the good performance of the proposed torque estimation scheme better than the simulation.
In order to verify the proposed torque equation, the experiment is carried out for various speeds and applied currents. The speed takes values of 500 rpm, 1000 rpm, and 1500 rpm. The current is chosen to be 10 A, 30 A, 50 A, 70 A, 100 A, 120 A, and 150 A. Figure 10 shows the torques and their percentage estimation errors against the q-axis currents for the three different speeds. The largest torque estimation error occurs for a current of 30 A at a speed of 500 rpm. As a whole, the estimation error comes within 5%.
Figure 11 shows torque estimation results when MTPA operation is performed while increasing the q-axis current from 0 A to 130 A in ramp form at a constant speed of 1500 rpm. It is possible to confirm that the proposed torque estimation method is more accurate than the conventional method in the steady state as well as in the transient state by following the actual torque better than the conventional torque calculation result.

6. Conclusions

This paper presented an adaptive torque estimation algorithm for IPMSMs with parameter variations including coupling effects and time-varying inductances. The key observations and conclusions are summarized below:
  • The d and q coupling effects were shown through experiments.
  • Such coupling effects and inductance uncertainties involved in the operation were represented using two equivalent mutual inductances and were estimated online.
  • A more elaborate torque equation was derived from the estimated inductances mentioned above.
For achieving more efficient operation over a wide range of speeds, future work is directed towards field weakening operations and MTPA algorithms, with consideration of the mutual inductances.

Author Contributions

Dooyoung Yang, Hyungsoo Mok, Juseok Lee and Soohee. Han conceived and designed the experiments; Dooyoung Yang performed the experiments; Dooyoung Yang, Hyungsoo Mok and Juseok Lee analyzed the data; and Dooyoung Yang, Soohee Han and Hyungsoo Mok wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jung, S.-Y.; Hong, J.; Nam, K. Current minimizing torque control of the IPMSM using Ferrari’s method. IEEE Trans. Power Electron. 2013, 28, 5603–5617. [Google Scholar] [CrossRef]
  2. Butt, C.B.; Hoque, M.A.; Rahman, M.A. Simplified fuzzy-logic-based MTPA speed control of IPMSM drive. IEEE Trans. Ind. Appl. 2004, 40, 1529–1535. [Google Scholar] [CrossRef]
  3. Do, T.D.; Kwak, S.; Choi, H.H.; Jung, J.-W. Suboptimal control scheme design for interior permanent-magnet synchronous motors: An SDRE-based approach. IEEE Trans. Power Electron. 2014, 29, 3020–3031. [Google Scholar] [CrossRef]
  4. Krishnan, R. Selection criteria for servo motor drives. IEEE Trans. Ind. Appl. 1987, 23, 270–275. [Google Scholar] [CrossRef]
  5. Rahman, M.A.; Zhou, P. Analysis of brushless permanent magnet synchronous motors. IEEE Trans. Ind. Electron. 1996, 43, 256–267. [Google Scholar] [CrossRef]
  6. Jung, S.-Y.; Mi, C.C.; Nam, K. Torque control of IPMSM in the field-weakening region with improved DC-link voltage utilization. IEEE Trans. Ind. Electron. 2015, 62, 3380–3387. [Google Scholar] [CrossRef]
  7. Kim, Y.-S.; Sul, S.-K. Torque control strategy of an IPMSM considering the flux variation of the permanent magnet. In Proceedings of the 42nd IAS Annual Meeting, Conference Record of the 2007 IEEE Industry Applications Conference, New Orleans, LA, USA, 23–27 September 2007; pp. 1301–1307. [CrossRef]
  8. Corley, M.J.; Lorenz, R.D. Rotor position and velocity estimation for a permanent magnet synchronous machine at standstill and high speeds. In Proceedings of the 31th IAS Annual Meeting (IAS ‘96), Conference Record of the 1996 IEEE Industry Applications Conference, San Diego, CA, USA, 6–10 October 1996; pp. 36–41. [CrossRef]
  9. Chalmers, B.J. Influence of saturation in brushless permanent-magnet motor drives. IEE Proc. B—Electr. Power Appl. 1992, 139, 51–52. [Google Scholar] [CrossRef]
  10. Choi, C.; Lee, W.; Kwon, S.-O.; Hong, J.-P. Experimental estimation of inductance for interior permanent magnet synchronous machine considering temperature distribution. IEEE Trans. Magn. 2013, 49, 2990–2996. [Google Scholar] [CrossRef]
  11. Morimoto, S.; Sanada, M.; Takeda, Y. Effects and compensation of magnetic saturation in flux-weakening controlled permanent magnet synchronous motor drives. IEEE Trans. Ind. Appl. 1994, 30, 1632–1637. [Google Scholar] [CrossRef]
  12. Nalepa, R.; Orlowska-Kowalska, T. Optimum trajectory control of the current vector of a nonsalient-pole PMSM in the field-weakening region. IEEE Trans. Ind. Electron. 2012, 59, 2867–2876. [Google Scholar] [CrossRef]
  13. Rang, G.; Lim, J.; Nam, K.; Ihm, H.-B.; Kim, H.-G. A MTPA control scheme for an IPM synchronous motor considering magnet flux variation caused by temperature. In Proceedings of the 19th Annual IEEE Applied Power Electronics Conference and Exposition, 2004. (APEC ’04), Anaheim, CA, USA, 22–26 February 2004; pp. 1617–1621. [CrossRef]
  14. Monajemy, R.; Krishnan, R. Implementation strategies for concurrent flux weakening and torque control of the PM synchronous motor. In Proceedings of the 30th IAS Annual Meeting, (IAS ’95), Conference Record of the 1995 IEEE Industry Applications Conference, Orlando, FL, USA, 8–12 October 1995; pp. 238–245. [CrossRef]
  15. Cheng, B.; Tesch, T.R. Torque feedforward control technique for permanent-magnet synchronous motors. IEEE Trans. Ind. Electron. 2010, 57, 969–974. [Google Scholar] [CrossRef]
  16. Rahman, M.F.; Zhong, L.; Lim, K.W. A direct torque-controlled interior permanent magnet synchronous motor drive incorporating field weakening. IEEE Trans. Ind. Appl. 1998, 34, 1246–1253. [Google Scholar] [CrossRef]
  17. Jeong, I.; Nam, K. Analytic expressions of torque and inductances via polynomial approximations of flux linkages. IEEE Trans. Magn. 2015, 51. [Google Scholar] [CrossRef]
  18. Zhong, L.; Rahman, M.F.; Hu, W.Y.; Lim, K.W. Analysis of direct torque control in permanent magnet synchronous motor drives. IEEE Trans. Power Electron. 1997, 12, 528–536. [Google Scholar] [CrossRef]
  19. Do, T.D.; Choi, H.H.; Jung, J.-W. Nonlinear optimal DTC design and stability analysis for interior permanent magnet synchronous motor drives. IEEE/ASME Trans. Mechatron. 2015, 20, 2716–2725. [Google Scholar] [CrossRef]
  20. Choi, Y.-S.; Choi, H.H.; Jung, J.-W. Feedback linearization direct torque control with reduced torque and flux ripples for IPMSM drives. IEEE Trans. Power Electron. 2016, 31, 3728–3737. [Google Scholar] [CrossRef]
  21. Yu, C.-Y.; Fukushige, T.; Limsuwan, N.; Kato, T.; Reigosa, D.D.; Lorenz, R.D. Variable-flux machine torque estimation and pulsating torque mitigation during magnetization state manipulation. IEEE Trans. Ind. Appl. 2014, 50, 3414–3422. [Google Scholar] [CrossRef]
  22. Xu, W.; Lorenz, R.D. Reduced parameter sensitivity stator flux linkage observer in deadbeat-direct torque and flux control for IPMSMs. IEEE Trans. Ind. Appl. 2014, 50, 2626–2636. [Google Scholar] [CrossRef]
  23. Ichikawa, S.; Tomita, M.; Doki, S.; Okuma, S. Sensorless control of permanent-magnet synchronous motors using online parameter identification based on system identification theory. IEEE Trans. Ind. Electron. 2006, 53, 363–372. [Google Scholar] [CrossRef]
  24. Liu, K.; Zhu, Z.Q. Mechanical parameter estimation of permanent-magnet synchronous machines with aiding from estimation of rotor PM flux linkage. IEEE Trans. Ind. Appl. 2015, 51, 3115–3125. [Google Scholar] [CrossRef]
  25. Yang, S.-M.; Lin, J.-D. Observer-based automatic control loop tuning for servo motor drives. In Proceedings of the 2013 IEEE 10th International Conference on Power Electronics and Drive Systems (PEDS), Kitakyushu, Japan, 22–25 April 2013; pp. 302–305. [CrossRef]
  26. Liu, K.; Zhu, Z.Q. Online estimation of the rotor flux linkage and voltage-source inverter nonlinearity in permanent magnet synchronous machine drives. IEEE Trans. Power Electron. 2014, 29, 418–427. [Google Scholar] [CrossRef]
  27. Underwood, S.J.; Husain, I. Online parameter estimation and adaptive control of permanent-magnet synchronous machines. IEEE Trans. Ind. Electron. 2010, 57, 2435–2443. [Google Scholar] [CrossRef]
  28. Wang, X.; Xie, W.; Dajaku, G.; Kennel, R.M.; Gerling, D.; Lorenz, R.D. Position self-sensing evaluation of novel CW-IPMSMs with an HF injection method. IEEE Trans. Ind. Appl. 2014, 50, 3325–3334. [Google Scholar] [CrossRef]
  29. Xu, W.; Lorenz, R.D. High-frequency injection-based stator flux linkage and torque estimation for DB-DTFC implementation on IPMSMs considering cross-saturation effects. IEEE Trans. Ind. Appl. 2014, 50, 3805–3815. [Google Scholar] [CrossRef]
  30. Lee, S.-B.; Habetler, T.G. A remote and sensorless thermal protection scheme for small line-connected ac machines. IEEE Trans. Ind. Appl. 2003, 39, 1323–1332. [Google Scholar] [CrossRef]
  31. Yukimori, I.; Koji, Y.; Shigeo, M.; Masayuki, S. Improvement of position estimation accuracy in IPMSM sensorless drives with on-line parameter identification. IEEJ Trans. Ind. Appl. 2010, 130, 867–873. [Google Scholar] [CrossRef]
  32. Nakatsugawa, J.; Iwasaki, N.; Nagura, H.; Iwaji, Y. Proposal of mathematical models taking into consideration magnetic saturation and cross-coupling effects in permanent magnet synchronous motors. IEEJ Trans. Ind. Appl. 2008, 128, 1212–1220. [Google Scholar] [CrossRef]
Figure 1. d- and q-axis flux variations due to the d- and q-axis currents: (a) d-axis flux; (b) q-axis flux.
Figure 1. d- and q-axis flux variations due to the d- and q-axis currents: (a) d-axis flux; (b) q-axis flux.
Energies 10 00167 g001
Figure 2. Efficiency map of an Interior Type Permanent Magnet Synchronous Motor (IPMSM) at various operating points.
Figure 2. Efficiency map of an Interior Type Permanent Magnet Synchronous Motor (IPMSM) at various operating points.
Energies 10 00167 g002
Figure 3. Back electromotive force (EMF) estimation block diagram with feedforward and feedback current controls ( K * denote control gains).
Figure 3. Back electromotive force (EMF) estimation block diagram with feedforward and feedback current controls ( K * denote control gains).
Energies 10 00167 g003
Figure 4. Overall control block diagram.
Figure 4. Overall control block diagram.
Energies 10 00167 g004
Figure 5. d- and q-axis fluxes (simulation): (a) d-axis flux; (b) q-axis flux.
Figure 5. d- and q-axis fluxes (simulation): (a) d-axis flux; (b) q-axis flux.
Energies 10 00167 g005
Figure 6. Percentage torque estimation errors due to the parameter variations: (a) existing torque Equation (3); (b) proposed torque Equation (12).
Figure 6. Percentage torque estimation errors due to the parameter variations: (a) existing torque Equation (3); (b) proposed torque Equation (12).
Energies 10 00167 g006
Figure 7. Experimental set up.
Figure 7. Experimental set up.
Energies 10 00167 g007
Figure 8. Comparison of real and estimated torques. (①: conventional torque equation, ②: proposed torque equation, ③: real torque by torque sensor).
Figure 8. Comparison of real and estimated torques. (①: conventional torque equation, ②: proposed torque equation, ③: real torque by torque sensor).
Energies 10 00167 g008
Figure 9. Percentage torque estimation errors due to the parameter variations (experiment): (a) conventional torque Equation (3); (b) proposed torque Equation (12).
Figure 9. Percentage torque estimation errors due to the parameter variations (experiment): (a) conventional torque Equation (3); (b) proposed torque Equation (12).
Energies 10 00167 g009
Figure 10. Percentage torque error rates for 500, 1000, 1500 rpm.
Figure 10. Percentage torque error rates for 500, 1000, 1500 rpm.
Energies 10 00167 g010aEnergies 10 00167 g010b
Figure 11. Dynamic performance of the proposed torque estimation algorithm.
Figure 11. Dynamic performance of the proposed torque estimation algorithm.
Energies 10 00167 g011
Table 1. Parameter values.
Table 1. Parameter values.
Rated Power15 [kW]
V d c 135 [V]
I max 250 [A]
Max speed4500 [rpm]
Rated speed1500 [rpm]
Rated torque70 [Nm]
The number of poles16
R 0.0128 [Ω]
L d 0 0.22 [mH]
L q 0 0.28 [mH]
λ m 0 0.0442 [V/rad/s]

Share and Cite

MDPI and ACS Style

Yang, D.; Mok, H.; Lee, J.; Han, S. Adaptive Torque Estimation for an IPMSM with Cross-Coupling and Parameter Variations. Energies 2017, 10, 167. https://doi.org/10.3390/en10020167

AMA Style

Yang D, Mok H, Lee J, Han S. Adaptive Torque Estimation for an IPMSM with Cross-Coupling and Parameter Variations. Energies. 2017; 10(2):167. https://doi.org/10.3390/en10020167

Chicago/Turabian Style

Yang, Dooyoung, Hyungsoo Mok, Jusuk Lee, and Soohee Han. 2017. "Adaptive Torque Estimation for an IPMSM with Cross-Coupling and Parameter Variations" Energies 10, no. 2: 167. https://doi.org/10.3390/en10020167

APA Style

Yang, D., Mok, H., Lee, J., & Han, S. (2017). Adaptive Torque Estimation for an IPMSM with Cross-Coupling and Parameter Variations. Energies, 10(2), 167. https://doi.org/10.3390/en10020167

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