I. INTRODUCTION
Interior permanent magnet motors (IPMSMs) have been widely used in a variety of industrial applications due to their high power density, high efficiency, high reliability and wide constant power operating range [1]-[3]. With the development of modern power electronics conversion technologies and computer control theory, multiphase motor driving systems have experienced rapid development. When compared with three-phase interior permanent magnet synchronous motors, multiphase motors have higher power density, reliability and fault-tolerant performance [4]. In addition, they are a better fit for aerospace, ship electric propulsion, electric vehicles, etc.
To achieve better performance for motor driving, maximum torque per ampere (MTPA) control strategies are often used. The MTPA control strategy was initially proposed in field-oriented-control (FOC) with the goal of delivering the same electromagnetic torque with the minimum current amplitude. The existing MTPA control strategies can be divided into three main types, mathematical model based methods, lookup-table (LUT) based methods and extremum seeking methods. In [5]-[7], the MTPA points were calculated directly with constant parameters. However, these parameters may change due to the influence of cross-saturation, magnetic saturation and temperature in practical application [8], [9]. To solve this problem, researchers adopted online parameter estimation algorithms in [10], [11] to improve the accuracy of MTPA tracking. However, estimation accuracy affect the MTPA tracking performances. In addition, the cross-saturation effect is not considered. LUT-based methods are employed in [12]-[14] for MTPA tracking. In this method, LUT data can be acquired by simulations or numerical analysis. However, its accuracy is affected by manufacturing error, variations of material characteristics and temperature in real application. To improve the accuracy of the algorithm, LUT data can be obtained by a series of experiments. However, this method is quite time consuming and resource consuming. In addition, it is impractical to test these data in every different motors. Therefore, the application scope of LUT-based methods are greatly limited.
Extremum seeking control strategies, which can also be called online-search-based MTPA control strategies [15], including the signal injection based MTPA control schemes [16]-[18], reach MTPA point by adjusting current angle through disturbance. These methods are parameter independent. However, the dynamic performance of search algorithm is relatively poor because the convergence rate is slow. Additionally, the signal-injection-based MTPA control strategies result in additional copper losses and torque ripple due to the real signal injected into the motor. The recently proposed virtual-signal-injection-based MTPA control strategy [19] can realize the online search of the MTPA point by virtually injecting a high frequency into the torque equation mathematically. Therefore, it does not cause additional copper losses or torque ripple. However, this method is only realized under FOC.
When compared with FOC, direct torque control (DTC) has better torque dynamic performance. Therefore, it is therefore widely applied in industry servo applications. DTC was first proposed in three-phase induction motor drives [20]. However, it has since been extended to permanent magnet synchronous motors (PMSMs) and multiphase motors. However, research on using the MTPA control strategy for DTC drives is not as common as that for FOC drives. The main reasons are that DTC drives focus more on fast dynamic performance which fits industry servo applications better. In addition, the current is hardly controlled in DTC drives. However, many emerging applications (like mining and steel industries) require more and more energy source. The MTPA control strategy, which can improve system efficiency, is of great demand.
Conventional DTC drives usually set the reference flux linkage to a fixed value, such as the permanent magnetic flux. To maintain a constant value, the stator current often contains more reactive power components, which increases the copper loss and reduces the motor efficiency. In [21], the bestefficiency stator flux linkage reference is found offline, which can give the optimal flux linkage via a different torque and speed. However, this LUT-based method requires a massive amount of offline work and considerable hardware resources. An MTPA control strategy based on signal injection was proposed in [16]. This method can realize the online optimization of the reference flux linkage and motor parameters are not required. However, it causes extra copper losses and torque ripple due to the injected signal.
In this paper, a new virtual signal injected MTPA control strategy is proposed for five-phase IPMSM drives. The high frequency signal is injected into the flux linkage amplitude mathematically and the resonant current is calculated and analyzed. This method does not cause additional copper losses or the resonant problems associated with real signal injection since it does not injects any real signals into the motor.
The organization of this paper is as follows. The relationship between the stator current amplitude and the stator flux linkage amplitude is analyzed in Section II. The proposed virtual signal injected MTPA algorithm is introduced in Section III. Simulation verifications are presented in Section IV. Experiments are implemented to verify the validity of the proposed VSI-MTPA algorithm in Section V. Finally, some conclusions are presented in Section VI.
II. RELATIONSHIP BETWEEN THE STATOR CURRENT AMPLITUDE AND THE STATOR FLUX LINKAGE AMPLITUDE
Fig. 1 shows a five-phase IPMSM prototype with 40 slots and 8 poles. Unequal pole arcs and shifting pole-pairs are applied to lower the torque ripple [22].
Fig. 1. Prototype of a five-phase IPMSM.
For a five-phase IPMSM, torque expression without considering magnetic saturation is as follow:
\(T_{e}=\frac{5 p}{2}\left[\psi_{m} i_{q}+\left(L_{d}-L_{q}\right) i_{d} i_{q}\right]\) (1)
\(\left\{\begin{array}{l} i_{d}=\left|i_{s}\right| \cos (\beta) \\ i_{q}=\left|i_{s}\right| \sin (\beta) \end{array}\right.\) (2)
Where id and iq are the d-axis and q-axis currents, Ld and Lq are the d-axis and q-axis stator inductances, ψm is the flux linkage of the permanent magnet, p is the number of pole pairs, |is| is the amplitude of the stator current, and β is the current vector angle.
From expression (1), the electromagnetic torque includes two parts. The former that is proportional to the q-axis current is the permanent magnet torque. The latter that is proportional to the product of the d-q axis currents is the reluctance torque. For the FOC, the MTPA point corresponds to the current angle β or d-axis current that can achieve the maximum torque for a given magnitude of current vector.
In the most common online-search-based MTPA algorithms (including signal-injection-based methods), the corresponding relation between the stator current magnitude and the current vector angle shown in Fig. 2 is used to judge the working condition of motors. As shown in Fig. 2, for current curves under a constant torque, the value of d|i|/dβ is zero at the lowest point of the curve, which corresponds to the MTPA point. When the current angle β is less than βMTPA, the value of d|i|/dβ is negative, or it is positive. This feature can be used to estimate the working condition of the motor and to search for the MTPA point online.
Fig. 2. Current amplitude locus with current angle under constant torque.
However, in DTC drive systems, the current angles are not directly regulated. Instead, the flux linkage amplitudes are used to control the whole system. Therefore, the feature of the current β cannot be directly applied in DTC drive systems
For the DTC, the torque expression (1) can be expressed as (3) in terms of the flux linkage with respect to the d-q reference frame.
\(T_{e}=\frac{5 p}{2}\left(\psi_{d} i_{q}-\psi_{q} i_{d}\right)\) (3)
\(\left\{\begin{array}{l} \psi_{d}=\left|\psi_{s}\right| \cos (\delta)=L_{d} i_{d}+\psi_{m} \\ \psi_{q}=\left|\psi_{s}\right| \sin (\delta)=L_{q} i_{q} \end{array}\right.\) (4)
Where ψd and ψq are the d-axis and q-axis flux linkage amplitudes in the d-q reference frame, respectively, |ψs| represents the amplitude of the stator flux linkage, and δ represents the angle of the flux linkage.
The stator current amplitude |is| can be expressed by the flux linkage amplitude |ψs| and the angle δ based on equation (4) as follows.
\(\left|i_{s}\right|=\sqrt{\left(\cos ^{2} \delta / L_{d}^{2}+\sin ^{2} \delta / L_{q}^{2}\right) \psi_{s}-2 \psi_{m} \cos \delta \psi_{s} / L_{d}^{2}+\psi_{m}^{2} / L_{d}^{2}}\) (5)
The relationship between the stator current and the stator flux linkage under the constant torque constrain can be obtained based on equations (1), (2) and (4). The result has been shown in Fig. 3. It can be seen that there exists a minimum in the objective function |is| for a constant torque which corresponds to the MTPA point. The flux linkage amplitude opposite the MTPA point is the MTPA flux linkage amplitude ψs_MTPA. The value of d|is|/d|ψs| is zero in the MTPA point. Therefore, the MTPA flux linkage amplitude ψs_MTPA can be calculated as follow by setting the value of d|is|/d|ψs| to zero. The result is as follow:
\(\psi_{s_{-} M T P A}=\psi_{m} \cos \delta /\left(\cos ^{2} \delta+\sin ^{2} \delta\left(L_{d} / L_{q}\right)^{2}\right)\) (6)
Fig. 3. Current amplitude locus with flux linkage amplitude under constant torque.
Equation (6) shows the optimal flux linkage amplitude which can achieve the minimum current amplitude for a constant torque, namely the MTPA point. However, the equation for ψs_MTPA contains three motor parameters. These parameters are the permanent magnet flux linkage ψm, and the inductances of d-axis and q-axis Ld and Lq. In real industrial applications, the motor parameters can change due to the change of motor working conditions. Therefore, the flux linkage amplitude provided by equation (6) may not be in accordance with the actual MTPA flux linkage amplitude.
As can be seen in Fig. 3, the corresponding relation between |i| and |ψs| is clear. The stator current amplitude |i| decreases until it reaches the bottom. Then it increases with an increasing of the flux linkage amplitude |ψs|. That is to say, the value of d|is|/d|ψs| is negative if |ψs| is less than ψs_MTPA. Meanwhile, it is positive if |ψs| is greater than ψs_MTPA. The slope, which is represented by d|is|/d|ψs|, approaches zero around the MTPA point.
These features can be used for MTPA tracking. The specific process of MTPA implementation will be introduced in the next section.
III. PROPOSED MTPA PROCEDURE
A. Method to Obtain Current Variation Information
For a DTC drive system, the electromagnetic torque can be expressed as follows by substituting (4) into (3):
\(T_{e}=\frac{5 p}{2}\left[\frac{1}{2}\left(\frac{1}{L_{q}}-\frac{1}{L_{d}}\right)\left|\psi_{s}\right| \sin (2 \delta)+\frac{\psi_{m}}{L_{d}} \sin (\delta)\right]\left|\psi_{s}\right|\) (7)
Equation (7) represents the base equation of a DTC algorithm for a five-phase IPMSM. As can be seen from the equation, the torque is mainly decided by two variables, the flux linkage amplitude |ψs| and the flux linkage angle δ in the d-q reference frame. In addition, the change of torque can be obtained by changing either of the two variables. A fast torque change can be obtained by varying the flux linkage angle δ, while variation of the flux linkage amplitude |ψs| is usually used to change the operating region of the motor or to reach particular working conditions [16]. That is to say, the torque is mainly decided by the flux linkage angle δ, and the flux linkage amplitude |ψs| has little effect. If a small high-frequency sinusoidal signal is added to the stator flux linkage amplitude, the torque fluctuates. However, the average torque is maintained as a stable value. Therefore, the torque can be assumed to be constant for a small high-frequency sinusoidal flux linkage amplitude variation.
The expression of the d-q axis current amplitude can be expressed by flux linkage amplitude in the d-q reference frame based on equation (8).
\(\left\{\begin{array}{l} i_{d}=\left(\psi_{d}-\psi_{m}\right) / L_{d} \\ i_{q}=\psi_{q} / L_{q} \end{array}\right.\) (8)
The square of the stator current amplitude can be calculated with the d-q axis current amplitude as follow:
\(\left|i_{s}\right|^{2}=i_{d}^{2}+i_{q}^{2}\) (9)
By injecting a small signal into the stator flux linkage amplitude, the results of d-axis and q-axis flux linkage amplitude containing perturbation are as follow:
\(\left\{\begin{array}{l} \psi_{d}^{h}=\left(\psi_{s}+\Delta \psi_{s}\right) \cos \delta \\ \psi_{q}^{h}=\left(\psi_{s}+\Delta \psi_{s}\right) \sin \delta \end{array}\right.\) (10)
Submitting (10) into (8), the corresponding d-axis and q-axis current after signal injection are as follows:
\(\left\{\begin{array}{l} i_{d}^{h}=\left(\psi_{d}^{h}-\psi_{d}\right) / L_{d}+i_{d} \\ i_{q}^{h}=\left(\psi_{d}^{h} i_{q}\right) / \psi_{q} \end{array}\right.\) (11)
Combining (9) and (11), the resultant current amplitude variation can be expressed as follows:
\(\left|i_{s}^{h}\right|=\sqrt{\left[i_{d}+\left(\psi_{d}^{h}-\psi_{d}\right) / L_{d}\right]^{2}+\left[\left(\psi_{q}^{h} i_{q}\right) / \psi_{q}\right]^{2}}\) (12)
From equation (12), the sum current magnitude variation in (5) due to the injected signal can be acquired from the measured d-q axis current, d-q axis flux linkage amplitude and the d-q axis flux linkage amplitude containing perturbation given by (10). The d-axis inductance is also necessary. When compared to the three-parameter dependent equation (6), there is only one parameter Ld that is needed. The values of Ld can be acquired from corresponding LUTs or nominal values since the value of the d-axis inductance remains relatively steady during operation [14]. It can be found that the current amplitude variation can be extracted by mathematical signal injection instead of real signal injection. Therefore, the proposed method can be classified as a virtual-signal-injection based method.
Fig. 4 shows the specific flow of the current amplitude calculation processing based on virtual signal injection. The final current amplitude variation can be calculated based on ψαβ, iαβ and θr, which represent the α-β-axis flux linkage amplitude, α-β-axis current amplitude in stationary coordinates, and rotor position, respectively. The current and flux linkage of d-axis and q-axis components can be obtained with a rotation transformation. Both the current and flux linkage are low-pass filtered to cut down undesirable noise. The magnitude and angle of the flux linkage can be calculated with ψdq, and the high-frequency signal is then directly injected into the flux linkage amplitude. The flux linkage of d-axis and q-axis containing a high-frequency component can be calculated through equation (10) The final current amplitude variation can be obtained based on equation (12). It should be noticed that the rotor positions are used only in the virtual signal injection block and are not involved in the DTC loop. Therefore, the virtual signal injection block does not influence the DTC system and does not affect the inherent superiority of the DTC. For example, it does not affect the fast response of the torque.
Fig. 4. Schematic of current amplitude variation calculation.
From Section II, it is found that the MTPA point can be obtained by comparing the variations between |ψs| and |is|. First, a sinusoidal signal is injected into the flux linkage amplitude |ψs| shown in Fig. 4. The reference flux linkage amplitude is gradually increased from values less than the optimal flux linkage amplitude to values greater than the optimal flux linkage amplitude. Meanwhile, the current magnitude variation and the perturbation signal, Asin(ωh)t, are required. In addition, they are filtered by the band-pass filter to extract the given frequency component. The results are shown in Fig. 5. At first, when the reference flux linkage amplitude is less than the optimal flux linkage amplitude, the resonant current is out of phase with the perturbation flux linkage. The resonant current becomes less relevant to the perturbation of the flux linkage when the reference flux linkage approaches the optimal one. Meanwhile, the phase relationship gets to be in phase when the reference flux linkage amplitude becomes greater than the optimal flux linkage amplitude. The current amplitude variation around the MPTA point is less than that of the other two cases. The injected signal frequency is 300 Hz. All of the above correspond to what was discussed in Section II, which proves the rationality of the proposed calculating method of current amplitude.
Fig. 5. Phase relationship between the |is| and |ψs|.
B. Signal Processing for the Extraction of the MTPA Criterion
When the perturbation is added into the flux linkage amplitude |ψs|, current amplitude with perturbation is as follow:
\(\left|i_{s}^{h}\right|=\left|i_{s}\right|+A \frac{d\left|i_{s}\right|}{d\left|\psi_{s}\right|} \sin \left(\omega_{h} t\right)+A^{2} \frac{d}{d\left|\psi_{s}\right|}\left(\frac{d\left|i_{s}\right|}{d\left|\psi_{s}\right|}\right) \sin ^{2}\left(\omega_{h} t\right)+\cdots\) (13)
d|is|/d|ψs| can be extracted from current amplitude containing a high-frequency |is h | with the signal processing method shown in Fig. 6. A band-pass filter with a center frequency equal to the high-frequency injected is utilized to extract a specific frequency ωh signal in |is h |. Then the signal output of the band-pass filter is multiplied with sin(ωht).
\(\begin{aligned} K \frac{d\left|i_{s}\right|}{d\left|\psi_{s}\right|} A \sin ^{2}\left(\omega_{h} t\right) &=K \frac{d\left|i_{s}\right|}{d\left|\psi_{s}\right|} A\left\{\frac{1}{2}\left[\cos (0)-\cos \left(2 \omega_{h} t\right)\right]\right\} \\ &=\frac{1}{2} K A \frac{d\left|i_{s}\right|}{d\left|\psi_{s}\right|}-\frac{1}{2} K A \frac{d\left|i_{s}\right|}{d\left|\psi_{s}\right|} \cos \left(2 \omega_{h} t\right) \end{aligned}\) (14)
d|is|/d|ψs| can be extracted by a first-order low-pass filter whose cutoff frequency is below the frequency ωh of the virtually injected signal. The output ε of the low-pass filter, which is proportional to d|is|/d|ψs|, is input to an integrator to obtain the compensation flux linkage amplitude of the reference value of the flux linkage amplitude.
Fig. 6. Schematic of the signal processing block.
Simulations on the change of the MTPA criterion value have also been implemented. The obtained results are shown in Fig. 7. At the same time, the current and flux amplitude waveforms are given. In Fig. 7, the dotted line represents the minimum current moment. The MTPA criterion value corresponding to the dotted line is about zero at the same time. It can also be seen from Fig. 7 that when the reference flux linkage amplitude is less than the MTPA flux linkage amplitude, the criterion value is negative. Meanwhile, it is positive when the reference flux linkage amplitude is greater than the MTPA flux linkage amplitude. This is consistent with the conclusions shown in Section II and the results shown in Fig. 5.
Fig. 7. Current magnitude, MTPA criterion value and reference flux linkage amplitude at a 6 Nm load torque.
C. Control System Configuration
When compared to three-phase IPMSMs, there exists high order harmonics in five-phase IPMSMs for it owns three dimensional subspaces [23]. To restrain high order harmonics, switch tables containing both α-β and α3-β3 reference frames have been used to improve the performance of IPMSMs [24].
Fig. 8 shows five-phase IPMSM space voltage vectors in both the α-β and α3-β3 frames. It can be found from Fig. 8 that the large voltage vector in the α-β frame becomes a small vector in the α3-β3 frame, the medium vector is still a medium vector in the α3-β3 frame, and the small vector becomes a large vector in the α3-β3 frame. In addition, the medium vector in the α-β frame is opposite that of the large vector and the small vector in the α3-β3 frame. This means that the effect on the flux linkage of the medium vector is opposite that of the large vector and the small vector in the α3-β3 frame. Therefore, appropriate selection of the medium voltage vector can guarantee both the controlled state of the torque and the flux linkage, and effectively suppress the magnitude of the flux linkage in the α3-β3 frame, which reduces the high order harmonics. According to the above analysis, new switch tables that consider both the α-β and α3-β3 frames have been proposed. These switch tables are shown in Table I and Table II.
Fig. 8. Voltage vector distributions. (a) Voltage vectors in the α-β frame. (b) Voltage vectors in the α3-β3 frame.
TABLE I IMPROVED SWITCH TABLE IN THE Α-Β FRAME
TABLE II SWITCH TABLE IN THE Α3-Β3 FRAME
Fig. 9 shows the control block diagram of a five-phase IPMSM drive containing the proposed VSI-MTPA algorithm. The fundamental torque control and flux linkage amplitude control loop are just the same as the conventional methods. The reference torque can be regulated by speed based on a PI controller just like the traditional one. The torque and the flux linkage estimator are shown in (15) and (16), respectively.
\(T_{e}=\frac{5}{2} p\left(\psi_{a} i_{\beta}-\psi_{\beta} i_{\alpha}\right)\) (15)
\(\left\{\begin{array}{l} \psi_{\alpha}=\int\left(u_{\alpha}-R i_{\alpha}\right) \\ \psi_{\beta}=\int\left(u_{\beta}-R i_{\beta}\right) \end{array}\right.\) (16)
Fig. 9. Control block diagram of the proposed virtual signal injection.
Instead of giving a designated value, the reference flux linkage amplitude is compensated by a VSI block. The details of the VSI block are shown in Fig.4 and Fig.6, and discussed above.
It should be noticed that the convergence rate of the VSI block does no influence the fast response of DTC drives. In addition, the virtually injected signal does not affect the torque or current because it does not inject any real signals. Therefore, the frequency of the injected signal is easy to choose. In addition, the integral parameter, shown in Fig. 6, that decides the convergence rate of the flux linkage amplitude can also be set to a higher value to improve the dynamic performance.
IV. SIMULATION VERIFICATION
Simulations have been implemented and compared with a prototype five-phase IPMSM drive system to demonstrate the effectiveness of the proposed VSI-MTPA algorithm. The motor parameters are shown in Table III.
TABLE III PARAMETERS OF THE IPMSM
In the simulation, the sampling time of the DTC drive system is 1μs, the load torque is 4 Nm, and the speed is set to 300r/min. First, the MTPA flux linkage amplitude is tested and the final value is 0.132Wb. As shown in Fig. 10(a), the referenced flux linkage amplitude is set to 0.132Wb. Then, it steps to 0.111Wb at 0.4s. After 0.6s, the proposed VSI-MTPA control strategy is adopted. After the proposed VSI-MTPA control strategy is used, the stator flux linkage increases gradually and finally approaches 0.132Wb, which is the MTPA flux linkage amplitude under the given load torque. As shown in Fig. 10(b), the current amplitude steps from 3.6A to 4.5A when the flux linkage amplitude decreases from 0.132Wb to 0.111Wb. After the proposed MTPA strategy is adopted, the current amplitude gradually decreases and finally reaches its initial value. The torque shown in Fig. 10(c) fluctuates slightly at 0.4s due to current changes. The final torque becomes as stable as that before 0.4s. It can also be observed that the torque fluctuation is less than 10%. Moreover, there is no additional torque ripple after 0.6s since the proposed method is based on virtual signal injection. Variations of the MTPA criterion value shown in Fig. 10(d), correspond to what was discussed in Section III. The magnitude of Δψs, shown in Fig. 10(e), corresponds to the compensation required for the stator flux linkage magnitude in Fig. 10(a). The speed performance is shown in Fig. 10(f). As shown in Fig. 10(f), the speed fluctuation is less than 1%. The speed fluctuates slightly at 0.4s due to the torque wave. Then it quickly recovers to the normal condition. Therefore, the algorithm can be shown to have no effect on the rotational speed performance. The simulation result in Fig. 10 shows that the proposed MTPA tracking method can quickly optimize the flux linkage amplitude to the MTPA value and that it does not cause extra torque ripple. The proposed method is shown to have excellent steady state performance.
Fig. 10. Simulated results. (a) Stator flux linkage amplitude and reference flux linkage amplitude tracking performance. (b) Variation of the stator current amplitude. (c) Variation of the electromagnetic torque. (d) Variation of the MTPA criterion value. (e) Variation of Δψs. (f) Locus of the speed.
As mentioned before, the proposed method is virtualsignal-injected based. In order to show this performance, FFT analyses of the flux linkage in the α-β coordinates for the constant flux linkage amplitude algorithm and the VSI-MTPA algorithm have been provided in Fig. 11. The fundamental value is equal to one. The fundamental frequency is 20Hz, and the injected signal frequency is 150Hz. The THD is 1.23% for the constant flux linkage amplitude algorithm and 0.99% for the VSI-MTPA algorithm. From Fig. 11, it can be seen that the harmonic proportion of the VSI-MTPA algorithm is basically the same as that of the constant flux amplitude algorithm, which demonstrates that the proposed VSI-MTPA algorithm does not cause extra harmonics.
Fig. 11. Comparison of the α-β flux linkage harmonic proportion by the constant flux amplitude algorithm and the VSI MTPA algorithm.
Fig. 12 shows simulation results of the torque and flux linkage amplitude for load torque step. The ideal MTPA flux linkage amplitude value at 6 Nm is shown with dotted line. From Fig. 12, the proposed VSI-MTPA algorithm show good tracking performance for load step. It can be concluded that the proposed method has excellent dynamic performance. Detailed flux linkages in the d-q coordinate are shown in Fig. 13. It can be seen that the d-q flux linkages have a fast response to load steps. Fig. 14 shows the trajectories of currents in the d-q coordinates. The dotted line represents the MTPA d-axis current at 6Nm. It can be seen that the currents in the d-q coordinates can response to the load step and show good accuracy in tracking the theoretical MTPA d-axis current.
Fig. 12. Stator flux linkage and electromagnetic torque tracking response to the load step.
Fig. 13. Flux linkages in the d-q coordinate tracking response to the load step.
Fig. 14. Currents in the d-q coordinate tracking response to the load step.
V. EXPERIMENT VERIFICATION
To verify the feasibility of the proposed VSI-MTPA algorithm for DTC drive systems, corresponding test platform has been built. The platform is shown in Fig. 15. The experimental platform includes a five-phase IPMSM, a five-phase inverter, and a magnetic powder brake. A dSPACE1005 controller is used to implement the whole algorithm. The torque is measured by a high precision torque sensor (YH502/20Nm). The sample rate of the sampling circuit is 10 kHz, and the given DC voltage is 80 V. In this experiment, the frequency of the injected signal in the VSI block is 150 Hz.
Fig. 15. Experimental platform.
Fig. 16 shows torque and current waveforms from the constant flux linkage amplitude algorithm to the VSI based MTPA algorithm. The dotted line represents the algorithm switching point. The current waveform is a bit distorted due to high order harmonics. Although a harmonic suppression method is employed, there are still some other high order harmonics, which deteriorate the current performance. However, the MTPA method itself effectively reduces the current amplitude after the algorithm switching. Meanwhile, it also keeps the average value of the torque constant at 4 Nm before and after switching. The value of the current amplitude remains at 4.6A with the constant flux linkage amplitude algorithm and decreases to 3.7A after the algorithm is switched to the proposed VSI-MTPA method. The current amplitude declined by 0.9A, while the torque remained constant, which demonstrates the function of the proposed method.
Fig. 16. Torque and currents waveforms from the constant flux linkage amplitude algorithm to the VSI-MTPA operation.
Fig. 17 shows waveform of the speed and the d-q axis currents for algorithm switching. From Fig. 17, the speed fluctuation is 2.7%. The speed of the motor is kept stable during the algorithm switch period, and the d-q axis currents change quickly to track to the MTPA trajectory.
Fig. 17. Trajectories of the speed and d-q axis currents for algorithm switching.
Fig. 18 shows the d-q axis flux linkage amplitudes trajectories for algorithm switching. The d-q axis flux linkages change quickly after the algorithm is switched from the constant flux linkage magnitude to the VSI-MTPA algorithm. Fig. 17 and Fig. 18 prove that the proposed method is effective for MTPA tracking and that it does not effect the speed performance.
Fig. 18. Trajectories of the d-q axis flux linkage amplitudes for algorithm switching.
Experiment for the comparisons of the α-β flux linkage harmonic proportion by the constant flux amplitude algorithm and the VSI-MTPA algorithm have been implemented. The results are shown in Fig. 19. The fundamental value is equal to one. The fundamental frequency is 20Hz, and the frequency of the injected signal is 150Hz. The THD is 2.76% for the constant flux linkage amplitude algorithm, and 2.05% for the VSI-MTPA algorithm. From Fig. 19, the harmonic proportion of the VSI-MTPA algorithm is basically the same as that of constant flux amplitude algorithm. Therefore, it can be demonstrated that no actual injection signal occurs in the proposed method.
Fig. 19. Comparison of the α-β flux linkage harmonic proportion by the constant flux amplitude algorithm and the VSI-MTPA algorithm.
Experiments on the corresponding relations for the current magnitude, the stator flux linkage amplitude and the MTPA criterion value have been implemented. The results are shown in Fig. 20. When the reference flux linkage amplitude varies from 0.111Wb to 0.17Wb, the minimum values of the current amplitude correspond to the zero point of ε, which is in accordance with the content discussed in Section III. The feasibility of the algorithm is verified in principle by the experiment. After 30s, the proposed algorithm is used. The values of |i|, ε and |ψs| step to 5.2A, 0 and 0.145Wb, respectively. All three of these parameters step quickly because the compensation of the flux linkage amplitude reached the MTPA compensation value before 30s. Then the values of the three parameters step according to the MTPA value when the algorithm is used. From the current amplitude |i| and the ε trajectory, it is clear that the final current amplitude adjusted by the proposed VSI-MTPA is equal to the minimum values from before 30s. In addition, ε fluctuates around zero which means that it is a MTPA point. Therefore, the proposed method can be shown to accurately track the MTPA flux linkage amplitude. The experimental results are basically the same as the simulated ones. There are some small errors in the numerical values due to the simulated model errors, and manufacturing and measuring errors of the motor. However, both the simulations and experimental results demonstrate the accuracy of the proposed algorithm.
Fig. 20. Corresponding relations for the current amplitude, MTPA criterion value and flux linkage amplitude at a 6 Nm load torque.
Fig. 21 shows waveforms of the torque and currents corresponding to a load step for the VSI-MTPA algorithm. It can be seen that both the electromagnetic torque and current of phase-A change quickly to a step change in the torque, which shows the good dynamic performance of the proposed algorithm. In addition, both the torque and current can maintain relative stability during the load step, which demonstrates the stability of the transient state of the VSI-MTPA algorithm.
Fig. 21. Torque and currents waveforms during a load step for the VSI based MTPA algorithm.
Fig. 22 shows the trajectories of electromagnetic torque and speed for a step change in torque. In Fig. 22, the electromagnetic torque changes quickly with the load step. The speed only fluctuates when the load steps, and it is quickly restored to its initial stable state within 2s. Fig. 23 shows the flux linkage tracking response to a step change in torque. The dotted line represents the theoretic MTPA flux linkage amplitude. It can be seen that the stator flux linkage amplitude changes quickly with the load step and that it accurately tracks the theoretical MTPA value. Fig. 24 and Fig. 25 show the trajectories of the d-q axis currents and the d-q axis flux linkage amplitudes to a step change in torque, where both the d-q axis currents and the d-q axis flux linkage amplitudes respond well to the load step. The obtained experimental results are also consistent with the simulated results. Although there are some small errors in the numerical values due to manufacturing and measuring errors, both the simulations and experimental results demonstrate the dynamic performance of the proposed algorithm. Fig. 26 shows the stator flux linkage amplitude tracking response to a step reduction in torque. It is obvious that the stator flux amplitude can still accurately track the theoretical MTPA value for load reduction. From Fig. 21 to Fig. 26, the proposed algorithm has been shown to possess excellent dynamic performance and good tracking accuracy.
Fig. 22. Trajectories of electromagnetic torque and speed for step change in torque.
Fig. 23. Stator flux linkage amplitude tracking response to step change in torque.
Fig. 24. Trajectories of d-q axis currents to step change in torque.
Fig. 25. Trajectories of the d-q axis flux linkages to a step change in torque.
Fig. 26. Stator flux linkage amplitude tracking response to a step reduction in torque.
Fig. 27 shows the trajectories of the speed and the d-q axis currents to a step change in speed. It is clear that the speed have a good performance for speed step. The d-q axis currents fluctuate when the speed steps and is quickly restored to its original value. Fig. 28 shows the trajectories of the speed and the stator flux linkage amplitude corresponding to a step change in speed. The stator flux linkage amplitude fluctuates when the speed steps and is restored to its original value from before the speed step. Both Fig. 27 and Fig. 28 show that the algorithm has a good speed dynamic performance and that the tracking of the currents or stator flux linkage amplitudes are independent of speed.
Fig. 27. Trajectories of the speed and d-q axis currents to a step change in speed.
Fig. 28. Trajectories of speed and stator flux linkage to a step change in speed.
Fig. 29 shows the starting performance of the proposed VSI-MTPA algorithm. The referenced speed is set as 300r/min, and the load torque is set as 2Nm. It can be seen that the proposed MTPA algorithm can quickly and accurately realize the load starting and track the MTPA point. As can be seen in both the torque and speed step experiments, all of the variables stay relative stable during the transient state. It can be concluded that the proposed method has a good stability for the transient state.
Fig. 29. Starting performance of the proposed MTPA algorithm at a 2Nm load torque.
From both the simulations and experimental results, the proposed method is shown to be effective for MTPA tracking and has a good accuracy for the seeking of the MTPA stator flux linkage amplitude. It should be noted that the VSI-MTPA algorithm has universal applicability since the formulas used to derive the algorithm can be used on all motors. Therefore, the algorithm can also be applied to other motors such as three-phase motors.
Although the experiment results have demonstrated the validity of the proposed method, it has some limitations. The value of the d-axis inductance is needed for the MTPA calculations. In addition, there are some fluctuations in the speed, flux and currents. However, these fluctuations can be reduced by new modulation methods, e.g. SVPWM [25].
VI. CONCLUSION
In this paper, an online MTPA tracking algorithm for DTC-based five-phase IPMSM drives has been proposed. The proposed method is VSI based and it mathematically injects a high frequency signal into the flux linkage amplitude of a DTC drive. MTPA points are obtained by using the relationship between the current magnitude and the stator flux linkage amplitude. Since this method does not inject any real signals, it does not result in additional copper loss or torque ripple when compared to real-signal-injection based methods. Simulations and experiments have been implemented to verify the proposed algorithm. These results demonstrate that the proposed method have a good performance for MTPA tracking.
참고문헌
- G. Wang, L. Yang, B. Yuan, B. Wang, G. Zhang, and D. Xu, “Pseudo-random high-frequency square-wave voltage injection based sensorless control of IPMSM drives for audible noise reduction,” IEEE Trans. Ind. Electron., Vol. 63, No. 12, pp. 7423-7433, Dec. 2016. https://doi.org/10.1109/TIE.2016.2594171
- Y. Kano, “Torque ripple reduction of saliency-based sensorless drive concentrated-winding IPMSM using novel flux barrier,” IEEE Trans. Ind. Appl., Vol. 51, No. 4, pp. 2905-2916, Jul./Aug. 2015. https://doi.org/10.1109/TIA.2015.2391431
- M. M. I. Chy and M. N. Uddin, “Development and implementation of a new adaptive intelligent speed controller for IPMSM drive,” IEEE Trans. Ind. Appl., Vol. 45, No. 3, pp. 1106-1115, May/Jun. 2009. https://doi.org/10.1109/TIA.2009.2018918
- Q. Chen, W. Zhao G. Liu, and Z. Lin, "Extension of virtual-signal-injection-based MTPA control for five-phase IPMSM into fault-tolerant operation" IEEE Trans. Ind. Electron., Vol: 66, No. 2, pp. 944-955, Feb. 2019. https://doi.org/10.1109/tie.2018.2826473
- T. M. Jahns, G. B. Kliman, and T. W. Neumann, “Interior permanent magnet synchronous motors for adjustablespeed drives,” IEEE Trans. Ind. Appl., Vol. IA-22, No. 4, pp. 738-747, Jul. 1986. https://doi.org/10.1109/TIA.1986.4504786
- S. Morimoto, M. Sanada, and Y. Takeda, “Wide-speed operation of interior permanent magnet synchronous motors with high-performance current regulator,” IEEE Trans. Ind. Appl., Vol. 30, No. 4, pp. 920-926, Jul./Aug. 1994. https://doi.org/10.1109/28.297908
- M. N. Uddin, T. S. Radwan, and M. A. Rahman, “Performance of interior permanent magnet motor drive over wide speed range,” IEEE Trans. Energy Convers., Vol. 17, No. 1, pp. 79-84, Mar. 2002. https://doi.org/10.1109/60.986441
- H. W. Lee, K.-D. Lee, W.-H. Kim, I.-S. Jang, M.-J. Kim, J.-J. Lee, and J. Lee, “Parameter design of IPMSM with concentrated winding considering partial magnetic saturation,” IEEE Trans. Magn., Vol. 47, No. 10, pp. 3653-3656, Oct. 2011. https://doi.org/10.1109/TMAG.2011.2157957
- X. Liu, H. Chen, J. Zhao, and A. Belahcen, “Research on the performances and parameters of interior PMSM used for electric vehicles,” IEEE Trans. Ind. Electron., Vol. 63, No. 6, pp. 3533-3545, Jun. 2016. https://doi.org/10.1109/TIE.2016.2524415
- R. Ni, D. Xu, G. Wang, L. Ding, G. Zhang, and L. Qu, “Maximum efficiency per ampere control of permanentmagnet synchronous machines,” IEEE Trans. Ind. Electron., Vol. 62, No. 4, pp. 2135-2143, Apr. 2015. https://doi.org/10.1109/TIE.2014.2354238
- P. Niazi, H. A. Toliyat, and A. Goodarzi, “Robust maximum torque per ampere (MTPA) control of PM-assisted SynRM for traction applications,” IEEE Trans. Veh. Technol., Vol. 56, No. 4, pp. 1538-1545, Jul. 2007. https://doi.org/10.1109/TVT.2007.896974
- T. M. Jahns, G. B. Kliman, and T. W. Neumann, “Interior permanent-magnet synchronous motors for adjustablespeed drives,” IEEE Trans. Ind. Appl., Vol. IA-22, No. 4, pp. 738-747, Jul. 1986. https://doi.org/10.1109/TIA.1986.4504786
- A. Consoli, G. Scarcella, G. Scelba, and A. Testa, “Steadystate and transient operation of IPMSMs under maximumtorque-per-ampere control,” IEEE Trans. Ind. Appl., Vol. 46, No. 1, pp. 121-129, Jan./Feb. 2010. https://doi.org/10.1109/TIA.2009.2036665
- S. Y. Jung, J. Hong, and K. Nam, “Current minimizing torque control of the IPMSM using Ferrari’s method,” IEEE Trans. Power Electron., Vol. 28, No. 12, pp. 5603-5617, Dec. 2013. https://doi.org/10.1109/TPEL.2013.2245920
- R. S. Colby and D. W. Novotny, “An efficiency-optimizing permanent magnet synchronous motor drive,” IEEE Trans. Ind. Appl., Vol. 24, No. 3, pp. 462-469, May/Jun. 1988. https://doi.org/10.1109/28.2897
- S. Bolognani, L. Peretti, and M. Zigliotto, “Online MTPA control strategy for DTC synchronous-reluctance-motor drives,” IEEE Trans. Power Electron., Vol. 26, No. 1, pp. 20-28, Jan. 2011. https://doi.org/10.1109/TPEL.2010.2050493
- G. Liu, J. Wang, W. Zhao, and Q. Chen, “A novel MTPA control strategy for IPMSM drives by space vector signal injection,” IEEE Trans. Ind. Electron., Vol. 64, No. 12, pp. 9243-9252, Dec. 2017. https://doi.org/10.1109/TIE.2017.2711507
- R. Antonello, M. Carraro, and M. Zigliotto, “Maximumtorque-per-ampere operation of anisotropic synchronous permanent-magnet motors based on extremum seeking control,” IEEE Trans. Ind. Electron., Vol. 61, No. 9, pp. 5086-5093, Sep. 2014. https://doi.org/10.1109/TIE.2013.2278518
- T. Sun, J. Wang, and X. Chen, “Maximum torque per ampere (MTPA) control for interior permanent magnet synchronous machine drives based on virtual signal injection,” IEEE Trans. Power Electron., Vol. 30, No. 9, pp. 5036-5045, Sep. 2015. https://doi.org/10.1109/TPEL.2014.2365814
- I. Takahashi and T. Noguchi, “A new quick-response and high-efficiency control strategy of an induction motor,” IEEE Trans. Ind. Appl., Vol. IA-22, No. 5, pp. 820-827, Sep. 1986. https://doi.org/10.1109/TIA.1986.4504799
- J. Habibi and S. Vaez.-Zadeh, "Efficiency-optimizing direct torque control of permanent magnet synchronous machines," in Proc. IEEE Power Electron. Spec. Conf. (PESC), pp. 759-764, 2005.
- Q. Chen, G. Xu, W. Zhao, L. Liu, and Z. Lin, "Torque ripple reduction in five-phase IPM motors by lowering interactional MMF" IEEE Trans. Ind. Electron., Vol. 65, No. 11, pp. 8520-8531, Nov. 2018 https://doi.org/10.1109/TIE.2018.2807392
- L. Parsa and H. A. Toliyat, “Sensorless Direct torque control of five-phase interior permanent-magnet motor drives,” IEEE Trans. Ind. Appl., Vol. 43, No. 4, pp. 952-959, Jul. 2007. https://doi.org/10.1109/TIA.2007.900444
- X. Huang, G. Liu, H. Zhou, and J. Ji, "Direct thrust control for five-phase tubular linear PM motor based on thirdharmonic current suppression" in Proc. Int. Conf. Electr. Mach. Syst. (ICEMS), pp. 1-4, 2017.
- F. Yu, X. Zhang, M. Qiao, and C. Du, "The direct torque control of multiphase permanent magnet synchronous motor based on low harmonic space vector PWM," in Proc. IEEE Int. Conf. Ind. Technol., pp. 1-5, 2008.