A Finite-Set Integral Sliding Modes Predictive Control for a Permanent Magnet Synchronous Motor Drive System

Hidalgo, Hector; Orosco, Rodolfo; Huerta, Hector; Vazquez, Nimrod; Estrada, Leonel; Pinto, Sergio; de Castro, Angel

doi:10.3390/wevj15070277

Open AccessArticle

A Finite-Set Integral Sliding Modes Predictive Control for a Permanent Magnet Synchronous Motor Drive System

by

Hector Hidalgo

¹

,

Rodolfo Orosco

^2,*,

Hector Huerta

³,

Nimrod Vazquez

^2,*

,

Leonel Estrada

⁴

,

Sergio Pinto

^5,6

and

Angel de Castro

⁷

¹

Mechatronics Department, TecNM/Instituto Tecnológico Superior de Villa La Venta, Huimanguillo 86410, Mexico

²

Electronics Department, TecNM/Instituto Tecnológico de Celaya, Celaya 38010, Mexico

³

Department of Computational Sciences and Engineering, Universidad de Guadalajara/Centro Universitario de los Valles, Ameca 46600, Mexico

⁴

Electronics Department, TecNM/Instituto Tecnológico Superior del Sur de Guanajuato, Benito Juárez, Guanajuato 38980, Mexico

⁵

School of Digital Innovation, Instituto Técnico Superior Especializado (ITSE) de Panamá, Tocumen, Avenida Domingo Díaz, Panama City 07202, Panama

⁶

Faculty of Informatics, Electronics, and Communications, Central Campus, Universidad de Panama, Panama City 3366, Panama

⁷

Electronics and Communications Technology Department, Universidad Autonoma de Madrid, 28049 Madrid, Spain

^*

Authors to whom correspondence should be addressed.

World Electr. Veh. J. 2024, 15(7), 277; https://doi.org/10.3390/wevj15070277

Submission received: 22 May 2024 / Revised: 17 June 2024 / Accepted: 19 June 2024 / Published: 21 June 2024

(This article belongs to the Special Issue Electric Vehicle Technology Development, Energy and Environmental Implications, and Decarbonization)

Download

Browse Figures

Versions Notes

Abstract

:

Finite-set model predictive control (FS-MPC) is an easy and intuitive control technique. However, parametric uncertainties reduce the accuracy of the prediction. Classical MPC requires many calculations; therefore, the calculation time generates a considerable time delay in the actuation. This delay deteriorates the performance of the system and generates a significant current ripple. This paper proposes a finite-set integral sliding modes predictive control (FS-ISMPC) for a permanent magnet synchronous motor (PMSM). The conventional decision function is replaced by an integral sliding cost function, which has several advantages, such as robustness to parameter uncertainties, and convergence in finite time. The proposed decision function does not require the inductance and resistance parameters of the motor. In addition, the proposal includes compensation for the calculation delay of the control vector. The proposed control strategy was compared with traditional predictive control with delay compensation using a real-time hardware-in-the-loop (HIL) simulation. The results obtained from the comparison indicated that the proposed controller has a lower THD and computational burden.

Keywords:

computational burden; model predictive control; real-time simulation; HIL testing; permanent magnet synchronous motor

1. Introduction

The permanent magnet synchronous motor (PMSM) offers high efficiency, power density, and the ability to maintain constant power over a wide range of speeds. PMSM is widely used in industrial applications and electric vehicle traction systems. PMSM can be controlled using both trapezoidal and sinusoidal PWM techniques for speed and position. In [1], a genetic algorithm-based controller for speed and position is presented. The study compares the performance of trapezoidal and sinusoidal PWM commutation techniques, providing insights for optimizing motor control systems. The classic control methods used in PMSM drives are current vector control [2] and direct torque control [3]. However, traditional controllers have poor dynamic performance. Some of the problems include overshoot and torque fluctuations at low speeds. Finite-set model predictive control (FS-MPC) has gained high relevance in the last decade due to its implementation in digital systems. FSMPC does not require modulation, and nonlinear constraints can be easily incorporated [4].

Controllers based on FS-MPC for PMSM suffer from several drawbacks. The online control vector calculation time affects the performance of the controller [5]. Computation time includes filtering, sampling, and other factors, which generate a delay in the control signal. The state of the switches changes once the calculation is completed. Consequently, the measured currents present a ripple of significant amplitude [6,7]. Therefore, it is necessary to add compensation for the delay and improve the performance of the control system. The most widely used method for delay compensation is to increase the prediction and control horizon [8,9,10]. To increase the performance of the compensation, the control vector of the next sampling is calculated. This process is repeated due to the receding horizon principle. However, this compensation method depends on the accuracy of the model. On the other hand, parameter uncertainty is a challenge for FS-MPC, since parameters change due to PMSM conditions. The saturation effect of the magnetic field changes the inductance of the motor under load conditions. Temperature changes due to mechanical stress decrease the value of the constant flux linkage. In addition, the rising temperature increases the resistance of the winding [11]. Several solutions are presented in the literature to improve the robustness of the prediction. In practical applications, there is not only parametric uncertainty but also disturbances and nonlinearities. Therefore, disturbance estimation and compensation techniques are required [12,13,14]. In [15], a predictive current control based on a composite observer is proposed. The estimation method eliminates the influence of motor parameter variation and permanent magnet demagnetization on the current vector. The proposed observer is designed using sliding modes (SM) and Luenberger approaches. In [16], a moving horizon estimator is utilized; the effect of constant flux linkage is eliminated, and voltage errors caused by resistance and inductance are compensated.

The recent trend to solve the parameter mismatch has led to the study of model-free control (MFC) methods. In [17], a nonparametric predictive current control (NPCC) has been proposed that can predict future current behavior using the measured data instead of motor parameters. However, the performance of NPCC is deficient if two consecutive switching states are identical, causing a stagnation problem [18,19]. MFC can be performed using the ultra-local model (ULM) approach. The ULM technique can be implemented by utilizing the input and output measurement data of the system without prior knowledge of the system parameters. The mathematical model is replaced by a differential equation of order selected by the practitioner [20]. ULM has been successfully applied to MPCs for PMSM drives [21,22]. In [23], a model of extended affine ULM is built containing a two-order term based on affine arithmetic. The accuracy of the model is improved using extended affine ULM, even in the presence of nonlinear terms. ULM provides a detailed representation of the system behavior in a small region of operation, allowing for more accurate predictions and control strategies tailored to that specific region [24].

Several robust control techniques exist in the field of AC drive systems. SM control is a special nonlinear control strategy, which considers robustness as a part of the design process. This control strategy can also be adapted to power electronics due to the switching nature of the control law. A stability analysis, using Lyapunov theory, is the main advantage [25]. However, sliding chattering is related to the robustness of the overall control system. In [26], an adaptive integral sliding mode predictive control is proposed to solve the chattering problem and improve the dynamics of the system. Recent research has presented a study of MPC using the SM cost function [27]. In [28], a different cost function using SM theory is proposed for a current controller of a grid-connected three-phase bidirectional power inverter. The sliding cost function presents a fast dynamic response and reduces computational burden. In [29], a finite-set sliding modes predictive current control (FS-SMPC) is applied to a PMSM.

Apart from the control technique, its testing and validation have critical importance. HIL is a real-time simulation technique. HIL verifies control algorithms and plant models, saving cost and time with high fidelity. PMSM and power converters have been successfully simulated in real time. Moreover, parameters can be changed during the simulation [30]. In [31], HIL implementation using digital signal processors (DSPs) and field programmable gate arrays (FPGA) for electric machine drives has been investigated. In addition, the drawbacks of both methods, such as complexity in development and verification, are presented. In [32], a high-fidelity real-time simulation using HIL is presented; the simulation is implemented on an FPGA platform; and HIL testing is validated against finite element analysis and experimental results.

In this paper, a robust model predictive control based on sliding mode theory is proposed for the current control of the PMSM drive system, instead of the traditional FS-MPC, which considers an error square cost function. The chattering effect is reduced by an integral sliding surface and an extended control set. To eliminate the delay due to calculation time, a one-step-ahead control horizon is implemented. The contributions of the paper are: (1) A robust current controller for PMSM is designed using an integral sliding cost function. (2) The computational burden of the proposed controller is low compared to traditional FS-MPC. (3) The designed controller is simple, and real-time simulation demonstrates the parameter robustness. This paper is organized as follows: In Section 2, the mathematical model of PMSM and traditional FS-MPC is introduced. In Section 3, the proposed FS-ISMPC strategy is described. Real-time simulation is demonstrated to verify the effectiveness of the proposed method in Section 4, which is followed by a discussion in Section 5. Finally, the paper is concluded in Section 6.

2. Mathematical Model of PMSM and FS-Strategies

2.1. The Mathematical Model of PMSM

In this study, the mathematical model of a PMSM in a dq synchronous reference frame is described as follows [33]:

\{\begin{matrix} \frac{d i_{d}}{d t} = - \frac{R}{L_{d}} i_{d} + \frac{L_{q}}{L_{d}} i_{q} ω_{e} + \frac{1}{L_{d}} u_{d} + f_{d} \\ \frac{d i_{q}}{d t} = - \frac{R}{L_{q}} i_{q} - \frac{L_{d}}{L_{q}} i_{d} ω_{e} - \frac{ω_{e} ψ_{f}}{L_{q}} + \frac{1}{L_{q}} u_{q} + f_{q} . \end{matrix}

(1)

In this model,

u_{d}

and

u_{q}

, represent the d and q axis voltages, respectively;

i_{d}

and

i_{q}

, represent the d and q axis currents, respectively;

L_{d}

and

L_{q}

are the stator inductance;

R

is the winding resistance;

ψ_{f}

is the flux linkage of permanent magnets;

ω_{e}

is the electrical angular velocity;

f_{d}

and

f_{q}

are defined to represent the disturbance caused by the parameter variation. The parameters of the motor are mainly influenced by the temperature and magnetic saturation. Therefore, the deviation of the internal parameters is ranged in a finite scale. The mechanical motion equation and the electromagnetic torque of PMSM are written as follows:

\{\begin{matrix} \frac{d ω_{m}}{d t} = \frac{T_{e}}{J} - \frac{T_{L}}{J} - \frac{B}{J} ω_{m} \\ T_{e} = \frac{3}{2} n_{p} [ψ_{f} + (L_{d} - L_{q}) i_{d}] i_{q}, \end{matrix}

(2)

where

T_{e}

and

T_{L}

, represent the electromagnetic torque and load torque, respectively;

B

represents the viscous friction;

ω_{m}

is the mechanical angular velocity;

J

is the rotor inertia;

n_{p}

is the number of pole pairs. Moreover, for surface-mounted PMSM inductance

L_{d}

is equal to

L_{q}

, therefore electromagnetic torque can be simplified as follows:

T_{e} = K_{t} i_{q},

(3)

where

K_{t}

is the torque constant, assumed to be equal to

(3 / 2) n_{p} ψ_{f}

.

2.2. Traditional FS-MPC Strategy

In power electronics, MPC can be implemented considering a finite number of control actions or switching states. FS-MPC is intuitive, the first step is obtaining the dynamic model of the system. According to the Euler equation and neglecting the perturbation terms, (1) can be written in discrete form as follows:

i_{d q}^{p} (k + 1) = A (k) i_{d q} (k) + B u_{d q} (k) + C (k),

(4)

where

k

and

k + 1

represent the value of the current moment and the next moment, respectively,

T_{s}

is the sampling period, and written as follows:

\begin{matrix} i_{d q}^{p} (k + 1) = [\begin{matrix} i_{d}^{p} (k + 1) \\ i_{q}^{p} (k + 1) \end{matrix}]; i_{d q} (k) = [\begin{matrix} i_{d} (k) \\ i_{q} (k) \end{matrix}]; \\ u_{d q} (k) = [\begin{matrix} u_{d} (k) \\ u_{q} (k) \end{matrix}]; A (k) = [\begin{matrix} 1 - \frac{T_{s} R}{L_{d}} & ω_{e} (k) T_{s} \\ - ω_{e} (k) T_{s} & 1 - \frac{T_{s} R}{L_{q}} \end{matrix}]; \\ B = [\begin{matrix} \frac{T_{s}}{L_{d}} & 0 \\ 0 & \frac{T_{s}}{L_{q}} \end{matrix}]; C (k) = [\begin{matrix} 0 \\ - \frac{ω_{e} (k) ψ_{f} T_{s}}{L_{q}} \end{matrix}] . \end{matrix}

The PMSM is driven by a two-level three-phase voltage source inverter. The switching states can be summarized in a single complex switching function given by:

S = \frac{2}{3} (S_{a} + S_{b} a + S_{c} a^{2}),

(5)

where

S_{a}

,

S_{b}

, and

S_{c}

are the switching functions of the inverter,

V_{d c}

is the DC source voltage and

a

is a unitary vector equal to

e^{j 2 π / 3}

. Considering all possible combinations, there are eight different voltage vectors, as depicted in Figure 1. To ensure that the output current follows a specified trajectory, a cost function is defined as follows [3]:

J_{k} = {[i_{d}^{r e f} - i_{d}^{p} (k + 1)]}^{2} + {[i_{q}^{r e f} - i_{q}^{p} (k + 1)]}^{2},

(6)

where

i_{d}^{r e f}

is the d axis reference current;

i_{q}^{r e f}

is the q axis reference current. When field-oriented control is adopted

i_{d}^{r e f}

= 0, and

i_{q}^{r e f}

is obtained from the speed controller.

Considering the digital delay in the control action, a one-step-ahead control horizon is used. The cost function (6) evaluates the error two sampling periods ahead. The traditional cost function with delay compensation is defined as follows [34]:

J_{k + 1} = {[i_{d}^{r e f} - i_{d}^{p} (k + 2)]}^{2} + {[i_{q}^{r e f} - i_{q}^{p} (k + 2)]}^{2} .

(7)

2.3. FS-SMPC Strategy

MPC is applied in a way that an optimal control sequence is given every sampling period. Moreover, the prediction is performed in an open loop. Therefore, stability is not guaranteed [35]. In sliding mode control theory, the stability of the closed-loop system is proved using a Lyapunov function [25].

The space vector representation of (1) in the

α β

stationary reference frame is given by:

\frac{d i_{s}}{d t} = - \frac{R}{L_{s}} i_{s} - \frac{1}{L_{s}} ξ + \frac{1}{L_{s}} u_{s} + f_{s},

(8)

where

u_{s}

is the stator voltage,

f_{s}

is the disturbance,

i_{s}

is the stator current,

θ_{e}

is the electrical rotor position,

L_{s} = L_{d} = L_{q}

, and

\frac{d i_{s}}{d t} = [\begin{matrix} \frac{d i_{\propto}}{d t} \\ \frac{d i_{β}}{d t} \end{matrix}]; i_{s} = [\begin{matrix} i_{\propto} \\ i_{β} \end{matrix}]; u_{s} = [\begin{matrix} u_{\propto} \\ u_{β} \end{matrix}]; ξ = j ω_{e} ψ_{f} e^{j θ_{e}} .

A sliding surface is designed as the difference between the measured current and the current reference. The sliding surface can be written as follows:

σ = i_{s} - i_{s}^{r e f} .

(9)

The high-frequency switching ripple is directly related to the stator current vector

i_{s}

and the reference current

i_{r e f}

is a pure sinusoidal waveform. Taking the first-time derivative of the sliding surface, and considering (8), this yields the following:

\dot{σ} \approx \frac{d i_{s}}{d t} .

(10)

This leads to the following cost function and optimization problem:

\begin{matrix} J_{o p t} = \min \{σ^{T} \dot{σ}\} < 0 \\ s . t . u_{s} \in V_{d c} S . \end{matrix}

(11)

The above expression can be simplified since

L_{s}

is always positive and

R i_{s} ≪ u_{s}

. The cost function (11) can be written in discrete form as follows [26,27]:

J_{k}^{s m} = {(i_{s} (k) - i_{s}^{r e f} (k))}^{T} (- ξ (k) + u_{s} (k)) .

(12)

3. Proposed FS-ISMPC

3.1. Integral Sliding Mode Predictive Control

The integral sliding mode control algorithm is an attractive alternative. The order of the motion equation in integral sliding mode is equal to the order of the original system [26]. Therefore, the system robustness can be ensured by an overall response of the system starting from the initial time instance. Let

e_{s} = i_{s} - i_{s}^{r e f}

be the tracking error vector, and

η

be the control gain. Then, the integral sliding mode surface (ISMS) is adopted as follows [26]:

σ = e_{s} + η \int_{0}^{t} e_{s} d τ .

(13)

Replacing (8) into the first-time derivative of (13), and neglecting the perturbation term, yields the following:

\dot{σ} = - \frac{R}{L_{s}} i_{s} - \frac{1}{L_{s}} ξ + \frac{1}{L_{s}} u_{s} + η e_{s} .

(14)

ISMS is similar to a proportional-integral controller. In this work, the discrete integral term

e_{i}

is expressed as [36]:

e_{i} (k) = e_{s} (k) + e_{i} (k - 1) .

(15)

The discrete-time version of (13) is given by the following equation:

σ (k) = e_{s} (k) + η e_{i} (k) .

(16)

The cost function can be acquired by applying (11), and the following discrete expression is obtained using the following equation:

J_{k}^{i s m} = σ^{T} (k) (- \frac{R}{L_{s}} i_{s} (k) - \frac{1}{L_{s}} ξ (k) + \frac{1}{L_{s}} u_{s} (k) + η e_{s} (k)) .

(17)

From (17), it can be seen that the control input is optimally calculated based on the SM existence condition. Moreover, the same online optimization mechanism of typical FS-MPC is employed. The more negative the value of (17) is, the faster it reaches the sliding surface. This aspect is used in the cost function definition of the proposed predictive controller.

3.2. Stability Analysis

The stability is compared by introducing a continuous approximation of the control input. The cost function (17) selects the amplitude and the sign of the optimal vector. In SM, the control input has two parts, an equivalent control input part, and a discontinuous control input part [37]. According to [25], the optimal voltage vector of FS-MPC can be expressed by its continuous set optimal voltage vector and the quantization error. In order to explain the robustness of the proposed controller, the required optimal voltage must be as follows:

u_{o p t} = u_{e q} - L_{s} ε s i g n (σ),

(18)

where:

u_{o p t}

is the optimal voltage vector of FS-ISMPC,

u_{e q}

is the equivalent control input, and

ε

is the quantization error. In the case of FS-MPC, the continuous-set voltage vector is obtained based on the deadbeat solution of (7). However, the equivalent control is added to cancel known terms. The control input in (17) presents the same property of cancelation. From (14), the equivalent control can be obtained by applying the invariance condition:

u_{e q} = R i_{s} + ξ - {L_{s} k}_{i} e_{s} .

(19)

Consider the following Lyapunov function candidate:

V = \frac{1}{2} σ^{T} σ .

(20)

The system is stable if the variation of

V

is decreasing:

\dot{V} = σ^{T} \dot{σ} .

(21)

Then, using (8) and (21), the following expression is obtained:

\dot{V} = σ^{T} (- ε s i g n (σ) + f_{s}) < 0 .

(22)

According to the above analysis, the sufficient condition for (21) to decrease is

‖ε‖ > 0

. This is met due to the negative value of the cost function, which is mentioned in (17). Under the condition

‖ε‖ > ‖f_{s}‖

, the Equation (22) is negative and the surface (13) converges to zero in a finite time.

3.3. Delay Compensation

The discrete current predictive model of (8) can be written as follows:

i_{s}^{p} (k + 1) = (1 - \frac{T_{s} R}{L_{s}}) i_{s} (k) + \frac{1}{L_{s}} [- ξ (k) + u_{s} (k)] .

(23)

Considering the digital delay in the control input calculation, one step ahead should be shifted. Assuming that

ξ (k + 1) \approx ξ (k)

, (17) is modified as follows:

J_{k + 1}^{i s m} = σ^{T} (k + 1) (- \frac{R}{L_{s}} i_{s} (k + 1) - \frac{1}{L_{s}} ξ (k) + \frac{1}{L_{s}} u_{s} (k + 1) + η e_{s} (k + 1)) .

(24)

Similarly, (16) is modified as the following:

σ (k + 1) = e_{s} (k + 1) + η [e_{s} (k + 1) + e_{i} (k)] .

(25)

The above expression requires of one-step ahead error vector. The reference current is pure sinusoidal. Assuming that

i_{s}^{r e f} (k + 1) \approx i_{s}^{r e f} (k)

, one-step-ahead error can be calculated as follows:

e_{s} (k + 1) = i_{s}^{p} (k + 1) - i_{s}^{r e f} (k) .

(26)

According to [26], (24) can be simplified since the product

R i_{s}

is negligible. Moreover,

L_{s}

is always positive. The cost function (24) is modified as follows:

J_{k + 1}^{i s m} = σ^{T} (k + 1) [- ξ (k) + u_{s} (k + 1) + η e_{s} (k + 1)] .

(27)

3.4. Extended Control Set

FS-MPC utilizes two null vectors and six active vectors to control the system. The output of the controller is a switching pattern. Therefore, a modulation technique is not required. The transient performance of the controller is faster than a classical controller. However, the tracking current presents a large ripple. To improve the steady-state performance of the two-level power inverter, an additional six intermediate vectors are added to the existing set of voltage vectors, as shown in Figure 2.

The cost function (27) indicates the direction of the current change. Therefore, the decision function selects a voltage vector with magnitude

2 V_{d c} / 3

. This can be solved by adding a soft constraint on the rate of voltage change. This type of functionality cannot be included in a traditional SM theory. An additional term can be included in the cost function. The resulting cost function is expressed as follows:

J_{k + 1}^{T} = J_{k + 1}^{i s m} + λ ‖u_{s} (k + 1) - u_{s} (k)‖ .

(28)

where

J_{k + 1}^{T}

is the cost function with constraints, and

λ

is the weighting factor. The voltage vectors in the αβ frame and the corresponding duty cycles are shown in Table 1. These voltage vectors are stored in a look-up table avoiding the necessity for an online modulation algorithm. The duty cycle is compared with a triangular carrier signal and the intersections define the switching instants. The combination of switching states is applied for a corresponding period. Figure 3 shows the switching pattern of each phase over a sampling period. The carrier signal limits the maximum switching frequency and alternates the zero vectors included in

{\vec{V}}_{0}

. However, the variable switching frequency is the main disadvantage of this control method. The overall flowchart of the proposed FS-ISMPC is depicted in Figure 4. The control diagram is shown in Figure 5.

4. Performance and Evaluation

4.1. Simulation Test

The proposed FS-ISMPC was compared to FS-PMC and FS-SMPC via simulations using MATLAB/Simulink models. The simulation parameters of the PMSM are listed in Table 2. The sampling period

T_{s}

was set as 20 μs, the coefficient

η

= 0.12, and

λ

= 0.5. The PI speed controller parameters were set to

K_{p}

= 2 and

K_{i}

= 0.015.

Figure 6 shows the tracking current comparisons to evaluate the integral term in the cost function. It can be observed that FS-ISMPC presents the smallest steady-state tracking error. In this simulation, the parameter values of the PMSM were not deviated. The rotor speed was 1000 rpm and the load torque was 16 N∙m. In terms of the current ripple in the αβ frame, the extended control set was reflected in the current quality. Therefore, the proposed method has an apparent improvement in steady-state performance.

To verify the robustness of the proposed current controller, a surface plot comparison between FS-MPC, FS-SMPC, and FS-ISMPC is presented in Figure 7. Two different cases are included in this plot: the static current tracking root mean square error (RMSE) by changing

L_{s}

and

ψ_{f}

in the PMSM model. The simulated results show that the sensitivity to parameter variations is similar for both FS-MPC and FS-SMPC, whereas the RMSE of the proposed method is not significant. Compared with the above parameters, the resistance value presents a negligible impact on the performance of all three methods. On the other hand, it is well known that variations in the system parameter values can affect the performance of classical FS-MPC in terms of RMSE.

Table 3 compares and summarizes the proposed and existing strategies based on SM. It can be seen that the similarity between all three methods is that no modulator is required. Meanwhile, the variable switching frequency is the main disadvantage of the three strategies. However, the actuation constraint of the proposed strategy forces the switching frequency to remain in the vicinity of its maximum frequency. The proposed FS-ISMPC can be considered as an optimized discrete SMC control strategy with a multi-objective optimization and constraint-handling capability. Moreover, the proposed strategy requires weighting factors (WF) and tunable coefficients, which can be tuned by using heuristic methods. A large value of

η

in the predictive surface generates a larger transient time in the current response, as depicted in Figure 8. A large value of

λ

in the cost function decreases the robustness.

4.2. HIL Simulation

Discrete-time equations of PMSM were required for the real-time simulation. The difference equations were programmed using fixed-point in LabVIEW predefined DSP blocks. The PMSM model and proposed FS-ISMPC were implemented on an NI cRIO-9067 using a 40 MHz clock. In the case of the proposed FS-ISMPC, six PWM gate drive signals were generated by a digital I/O module NI 9401. An analog input module NI 9291 was used for stator current acquisition. In addition, an isolated frequency input module NI 9326 was used for encoder signal acquisition. On the other hand, the PMSM model required a DAC module NI 9262 to generate the stator current. The encoder signal was generated by a digital I/O module NI 9401. The schematic diagram of overall system interconnections is shown in Figure 9.

Table 4 compares the FPGA resource usage for the control strategies. The most demanding strategy in terms of resources is the classical FS-MPC due to the use of the mathematical model of PMSM in the optimization stage. The proposal presents a significant usage of available DSP48 blocks in comparison with FS-SMPC. The integral term requires arithmetic calculations based on adders/multipliers and accumulators. The calculation time of each strategy is listed in Table 5. The computational burden of the proposed FS-ISMPC is slightly lower than the classical FS-MPC.

4.3. Real-Time Simulation Results

To validate the effectiveness of the proposed strategy, the classical FS-MPC, FS-SMPC, and the proposed FS-ISMPC were compared. The PMSM model implemented in HIL was (8), whereas the controllers equations implemented were (7) for FS-MPC, (12) for FS-SMPC, and (28) for proposed FS-ISMPC. A speed reference step change was applied, and the rotor speed was changed from 500 to 1000 rpm, as can be observed in Figure 10, the dynamic response is very similar for all three strategies. The settling time is about 30 ms, which is a fast response for PMSM. The performance of all three control strategies is shown in Figure 11, with a step change in load torque

T_{L}

from 8 Nm to 16 Nm, while maintaining a constant reference speed of 1000 rpm. In this case, the response of the three strategies was similar, and the setting time was 50 ms. Additionally, according to the load step test, the rotor speed undershoot was about 4 rpm. The torque ripple was compared using the RMS value. According to the torque fluctuation analysis, the value of torque fluctuation of the FSMPC was 0.380 Nm, the value of FS-SMPC was 0.618 Nm, and the value of torque fluctuation of the proposed FS-ISMPC was 0.342 Nm.

Figure 12 illustrates the total harmonic distortion (THD) of

i_{α}

under different load torque values. The THD was tested using nominal parameter values at 1000 rpm. The FS-MPC and the proposed FS-ISMPC present similar values when the load torque is greater than 50%. The THD of the three strategies increased when the load torque was low. The full-load THD of FS-MPC was 2.8%. In the case of FS-SMPC, it was 3.48%. Finally, the THD of the proposed FS-ISMPC was 2.1%. Figure 13 shows the current tracking RMSE to evaluate the FS-ISMPC robustness under parameter mismatches using HIL simulation. The RMSE is tested at 1000 rpm and 8 N.m. As shown, the proposed FS-ISMPC exhibits less sensitivity to parameter variations compared to the FS-SMPC. Figure 14 shows the current harmonic spectrum using the classical FS-MPC, FS-SMPC, and the proposed FS-ISMPC. The experimental data are acquired by the oscilloscope and are processed in the FFT module of MATLAB. The THD was calculated in MATLAB following the methodology proposed in [38]. In the case of the proposal, a visible harmonic at 50 kHz exists due to the use of carrier signal and the extended control set.

5. Discussion

The HIL simulation results indicate that the amplitude of the torque ripple depends on the sampling time. Notably, these findings corroborate with previous research, as evidenced by agreement with [10]. In the case of the proposed controller, a larger sampling time can decrease the effects of calculation delay; however, it also increases the amplitude of torque ripple. Increasing the sampling time could decrease the computational burden and potentially eliminate the necessity for compensation Equations (23) and (25).

The effects of parameter variation are negligible due to the uncertainty considered in this work. The robustness of the controller depends on the sliding surface, which aligns with [27,28]. However, in practice, there may be conditions that deviate the parameters outside the bounds, decreasing prediction accuracy. On the other hand, the tracking current error is diminished due to integral action on the sliding surface. The RMSE results indicate that the proposed controller exhibits improved performance similar to that reported by [26].

A system stability test similar to [27] was not directly performed. However, the analogy proposed in (18) suggests that the robustness of the controller depends on the quantization error; the larger the error, the greater the robustness.

6. Conclusions and Future Works

In this paper, an FS-ISMPC is presented for controlling a PMSM. The cost function is replaced by an SM existence condition. The stability and robustness of the proposed FSISMPC are theoretically proven, but also via MATLAB and real-time HIL simulations. The inductance and resistance values are not required in the cost function since the integral term reduces the RMSE. The proposed controller offers remarkable accuracy, robustness, low computational burden, and easy tuning for implementation. The THD of the proposal is lower than classical FS-MPC, even at low load torque. According to the torque fluctuation analysis, the value of the torque fluctuation of the FS-MPC is 0.380 Nm, while the value of torque fluctuation of the FS-ISMPC is 0.342 Nm; consequently, the torque fluctuation is reduced by 11%.

In future works, it might be possible to study alternative methods to enhance rejection for more complex loads, such as time-varying torque. As an alternative, a fixed-switching method can be included to increase global efficiency. Additionally, another type of nonlinear sliding surface should be explored as a case study. Also, the implementation of the online system identification approach may enhance the robustness of the proposed current controller.

Author Contributions

Conceptualization, H.H. (Hector Hidalgo), L.E. and N.V.; investigation, H.H. (Hector Hidalgo) and L.E.; methodology, N.V. and H.H. (Hector Hidalgo); supervision, N.V., R.O., H.H. (Hector Huerta), A.d.C. and S.P.; writing—original draft, H.H. (Hector Hidalgo); writing—review and editing, N.V., S.P. and A.d.C.; sponsorship, L.E., N.V. and A.d.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been partially funded by the Instituto de Innovación, Ciencia y Emprendimiento para la Competitividad para el Estado de Guanajuato “IDEA GTO”, under grant IDEAGTO/CONV/146/2023; by TecNM, under grant 19549.24-P; and the project PID2022-137593OB-I00, financed by the MCIN/AEI/10.13039/501100011033/FEDER, UE, and FSE+.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

Bekiroglu, E.; Dalkin, A. Comparison of Trapezoidal and Sinusoidal PWM Techniques for Speed and Position Control of PMSM. Adv. Electr. Electron. Eng. 2020, 18, 207–216. [Google Scholar] [CrossRef]
Hu, H.; Wang, H.; Liu, K.; Wei, J.; Shen, X. A Simplified Space Vector Pulse Width Modulation Algorithm of a High-Speed Permanent Magnet Synchronous Machine Drive for a Flywheel Energy Storage System. Energies 2022, 15, 4065. [Google Scholar] [CrossRef]
Bao, G.; Qi, W.; He, T. Direct Torque Control of PMSM with Modified Finite Set Model Predictive Control. Energies 2020, 13, 234. [Google Scholar] [CrossRef]
Peng, J.; Yao, M. Overview of Predictive Control Technology for Permanent Magnet Synchronous Motor Systems. Appl. Sci. 2023, 13, 6255. [Google Scholar] [CrossRef]
Gao, J.; Gong, C.; Li, W.; Liu, J. Novel Compensation Strategy for Calculation Delay of Finite Control Set Model Predictive Current Control in PMSM. IEEE Trans. Ind. Electron. 2020, 67, 5816–5819. [Google Scholar] [CrossRef]
Fang, S.; Meng, J.; Wang, W.; Meng, Y.; Wang, Y.; Huang, D. Compensation Strategy of PMSM Predictive Control with Reduced Parameter Disturbance. Sustainability 2022, 14, 9868. [Google Scholar] [CrossRef]
Yang, Y.; Wen, H.; Li, D. A Fast and Fixed Switching Frequency Model Predictive Control With Delay Compensation for Three-Phase Inverters. IEEE Access 2017, 5, 17904–17913. [Google Scholar] [CrossRef]
Luo, Y.; Liu, C. Model Predictive Control for a Six-Phase PMSM Motor With a Reduced-Dimension Cost Function. IEEE Trans. Ind. Electron. 2020, 67, 969–979. [Google Scholar] [CrossRef]
Zhang, Y.; Xu, D.; Liu, J.; Gao, S.; Xu, W. Performance Improvement of Model-Predictive Current Control of Permanent Magnet Synchronous Motor Drives. IEEE Trans. Ind. Appl. 2017, 53, 3683–3695. [Google Scholar] [CrossRef]
Zhang, X.; Zhang, L.; Zhang, Y. Model Predictive Current Control for PMSM Drives With Parameter Robustness Improvement. IEEE Trans. Power Electron. 2018, 34, 1645–1657. [Google Scholar] [CrossRef]
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]
Sun, X.; Zhang, Y.; Cai, Y.; Tian, X. Compensated Deadbeat Predictive Current Control Considering Disturbance and VSI Nonlinearity for In-Wheel PMSMs. IEEE/ASME Trans. Mechatron. 2022, 27, 3536–3547. [Google Scholar] [CrossRef]
Sun, X.; Zhang, Y.; Lei, G.; Guo, Y.; Zhu, J. An Improved Deadbeat Predictive Stator Flux Control with Reduced-Order Disturbance Observer for In-Wheel PMSMs. IEEE/ASME Trans. Mechatron. 2021, 27, 690–700. [Google Scholar] [CrossRef]
Niu, S.; Luo, Y.; Fu, W.; Zhang, X. Robust Model Predictive Control for a Three-Phase PMSM Motor With Improved Control Precision. IEEE Trans. Ind. Electron. 2020, 68, 838–849. [Google Scholar] [CrossRef]
Zhang, C.; Wu, G.; Rong, F.; Feng, J.; Jia, L.; He, J.; Huang, S. Robust Fault-Tolerant Predictive Current Control for Permanent Magnet Synchronous Motors Considering Demagnetization Fault. IEEE Trans. Ind. Electron. 2017, 65, 5324–5334. [Google Scholar] [CrossRef]
Pei, G.; Li, L.; Gao, X.; Liu, J.; Kennel, R. Predictive Current Trajectory Control for PMSM at Voltage Limit. IEEE Access 2020, 8, 1670–1679. [Google Scholar] [CrossRef]
Yuan, X.; Zhang, S.; Zhang, C. Nonparametric Predictive Current Control for PMSM. IEEE Trans. Power Electron. 2020, 35, 9332–9341. [Google Scholar] [CrossRef]
Agustin, C.A.; Yu, J.-T.; Cheng, Y.-S.; Lin, C.-K.; Yi, Y.-W. A Synchronized Current Difference Updating Technique for Model-Free Predictive Current Control of PMSM Drives. IEEE Access 2021, 9, 63306–63318. [Google Scholar] [CrossRef]
Lin, C.; Yu, J.; Lai, Y.; Yu, H.; Lin, Y.; Chen, F. Simplified model-free predictive current control for interior permanent magnet synchronous motors. Electron. Lett. 2016, 52, 49–50. [Google Scholar] [CrossRef]
Fliess, M.; Join, C. Model-free control. Int. J. Control. 2013, 86, 2228–2252. [Google Scholar] [CrossRef]
Zhang, Y.; Jin, J.; Huang, L. Model-Free Predictive Current Control of PMSM Drives Based on Extended State Observer Using Ultralocal Model. IEEE Trans. Ind. Electron. 2021, 68, 993–1003. [Google Scholar] [CrossRef]
Li, X.; Wang, Y.; Guo, X.; Cui, X.; Zhang, S.; Li, Y. An Improved Model-Free Current Predictive Control Method for SPMSM Drives. IEEE Access 2021, 9, 134672–134681. [Google Scholar] [CrossRef]
Wei, Y.; Young, H.; Ke, D.; Wang, F.; Rodríguez, J. Model-Free Predictive Current Control Using Extended Affine Ultralocal for PMSM Drives. IEEE Trans. Ind. Electron. 2023, 71, 6719–6729. [Google Scholar] [CrossRef]
Mousavi, M.S.; Davari, S.A.; Nekoukar, V.; Garcia, C.; Rodriguez, J. Integral Sliding Mode Observer-Based Ultralocal Model for Finite-Set Model Predictive Current Control of Induction Motor. IEEE J. Emerg. Sel. Top. Power Electron. 2022, 10, 2912–2922. [Google Scholar] [CrossRef]
Qu, L.; Qiao, W.; Qu, L. Active-Disturbance-Rejection-Based Sliding-Mode Current Control for Permanent-Magnet Synchronous Motors. IEEE Trans. Power Electron. 2020, 36, 751–760. [Google Scholar] [CrossRef]
Li, Z.; Wang, F.; Ke, D.; Li, J.; Zhang, W. Robust Continuous Model Predictive Speed and Current Control for PMSM With Adaptive Integral Sliding-Mode Approach. IEEE Trans. Power Electron. 2021, 36, 14398–14408. [Google Scholar] [CrossRef]
Zheng, C.; Gong, Z.; Wu, X.; Dragicevic, T.; Rodriguez, J.; Blaabjerg, F. Finite-Set Quasi-Sliding Mode Predictive Control of $ LC$-Filtered Voltage Source Inverters. IEEE Trans. Ind. Electron. 2021, 69, 11968–11978. [Google Scholar] [CrossRef]
Estrada, L.; Vazquez, N.; Vaquero, J.; Hernandez, C.; Arau, J.; Huerta, H. Finite Control Set—Model Predictive Control Based On Sliding Mode For Bidirectional Power Inverter. IEEE Trans. Energy Convers. 2021, 36, 2814–2824. [Google Scholar] [CrossRef]
Hidalgo, H.; Vázquez, N.; Orosco, R.; Hernandez, C.; López, H.; Huerta, H.; Pinto, S. Robust Modulated Model Predictive Control for a PMSM Using Sliding Mode Cost Function. In Proceedings of the 2023 IECON—49th Annual Conference of the IEEE Industrial Electronics Society, Singapore, 16–19 October 2023; pp. 1–6. [Google Scholar]
Abdelrahman, A.S.; Algarny, K.S.; Youssef, M.Z. A Novel Platform for Powertrain Modeling of Electric Cars With Experimental Validation Using Real-Time Hardware in the Loop (HIL): A Case Study of GM Second Generation Chevrolet Volt. IEEE Trans. Power Electron. 2018, 33, 9762–9771. [Google Scholar] [CrossRef]
Lee, J.S.; Choi, G. Modeling and hardware-in-the-loop system realization of electric machine drives—A review. CES Trans. Electr. Mach. Syst. 2021, 5, 194–201. [Google Scholar] [CrossRef]
Alvarez-Gonzalez, F.; Griffo, A.; Sen, B.; Wang, J. Real-Time Hardware-in-the-Loop Simulation of Permanent-Magnet Synchronous Motor Drives Under Stator Faults. IEEE Trans. Ind. Electron. 2017, 64, 6960–6969. [Google Scholar] [CrossRef]
Liu, X.; Yu, H.; Yu, J.; Zhao, L. Combined Speed and Current Terminal Sliding Mode Control With Nonlinear Disturbance Observer for PMSM Drive. IEEE Access 2018, 6, 29594–29601. [Google Scholar] [CrossRef]
Cortes, P.; Rodriguez, J.; Silva, C.; Flores, A. Delay Compensation in Model Predictive Current Control of a Three-Phase Inverter. IEEE Trans. Ind. Electron. 2011, 59, 1323–1325. [Google Scholar] [CrossRef]
Li, T.; Sun, X.; Lei, G.; Guo, Y.; Yang, Z.; Zhu, J. Finite-Control-Set Model Predictive Control of Permanent Magnet Synchronous Motor Drive Systems-An Overview. IEEE/CAA J. Autom. Sin. 2022, 9, 2087–2105. [Google Scholar] [CrossRef]
Lynch, K.M.; Marchuk, N.; Elwin, M.L. PID Feedback Control. In Embedded Computing and Mechatronics with the PIC32; Elsevier: Amsterdam, The Netherlands, 2016; pp. 375–385. [Google Scholar] [CrossRef]
Huerta, H.; Loukianov, A.G.; Canedo, J.M. Passivity Sliding Mode Control of Large-Scale Power Systems. IEEE Trans. Control. Syst. Technol. 2018, 27, 1219–1227. [Google Scholar] [CrossRef]
Panchbhai, A.; Parmar, S.; Vaghasiya, C. THD comparison for power quality (MATLAB simulation & hardware). In Proceedings of the 2016 IEEE 1st International Conference on Power Electronics, Intelligent Control and Energy Systems (ICPEICES), Delhi, India, 4–6 July 2016; pp. 1–4. [Google Scholar]

Figure 1. Voltage vectors of the inverter.

Figure 2. Extended control set.

Figure 3. FS-ISMPC switching pattern for

{\vec{V}}_{7}

and

{\vec{V}}_{9}

.

Figure 3. FS-ISMPC switching pattern for

{\vec{V}}_{7}

and

{\vec{V}}_{9}

.

Figure 4. Flowchart of proposed FS-ISMPC.

Figure 5. Block diagram of proposed FS-ISMPC.

Figure 6. Simulated steady state performance evaluation in the αβ frame: (a) classical FS-MPC; (b) FS-SMPC; and (c) proposed FS-ISMPC.

Figure 7. Simulated performance with mismatches

{\tilde{L}}_{s}

and

{\tilde{ψ}}_{f}

: (a) classical FS-MPC; (b) FS-SMPC; and (c) proposed FS-ISMPC.

Figure 7. Simulated performance with mismatches

{\tilde{L}}_{s}

and

{\tilde{ψ}}_{f}

: (a) classical FS-MPC; (b) FS-SMPC; and (c) proposed FS-ISMPC.

Figure 8. State trajectory and sliding surface.

Figure 9. HIL simulation implementation block diagram.

Figure 10. Response of the PMSM under a speed step change: (a) classical FS-MPC; (b) FS-SMPC; and (c) proposed FS-ISMPC.

Figure 11. Response of the PMSM under a load step change: (a) classical FS-MPC; (b) FS-SMPC; and (c) proposed FS-ISMPC.

Figure 12. THD of stator current under different load conditions.

Figure 13. Control strategy performance: (a) inductance mismatch; and (b) permanent magnet flux mismatch.

Figure 14. HIL simulation current harmonic spectra: (a) classical FS-MPC; (b) FS-SMPC; and (c) proposed FS-ISMPC.

Table 1. Voltage vector extended control set.

Vector	Rectangular Form	Duty Cycle
Vector	Rectangular Form	D₀	D₁	D₂
${\vec{V}}_{0}$	$0$	0.5	0.5	0.5
${\vec{V}}_{1}$	$2 V_{d c} / 3$	1	0	0
${\vec{V}}_{2}$	$V_{d c} / 3 + j \sqrt{3} V_{d c} / 3$	1	1	0
${\vec{V}}_{3}$	$- V_{d c} / 3 + j \sqrt{3} V_{d c} / 3$	0	1	0
${\vec{V}}_{4}$	$- 2 V_{d c} / 3$	0	1	1
${\vec{V}}_{5}$	$- V_{d c} / 3 - j \sqrt{3} V_{d c} / 3$	0	0	1
${\vec{V}}_{6}$	$V_{d c} / 3 - j \sqrt{3} V_{d c} / 3$	1	0	1
${\vec{V}}_{7}$	$V_{d c} / 2 + j \sqrt{3} V_{d c} / 6$	1	0.5	0
${\vec{V}}_{8}$	$j \sqrt{3} V_{d c} / 3$	0.5	1	0
${\vec{V}}_{9}$	$- V_{d c} / 2 + j \sqrt{3} V_{d c} / 6$	0	1	0.5
${\vec{V}}_{10}$	$- V_{d c} / 2 - j \sqrt{3} V_{d c} / 6$	0	0.5	1
${\vec{V}}_{11}$	$- j \sqrt{3} V_{d c} / 3$	0.5	0	1
${\vec{V}}_{12}$	$V_{d c} / 2 - j \sqrt{3} V_{d c} / 6$	1	0	0.5

Table 2. Parameters of the PMSM.

Parameter	Value
Rated torque	20 Nm
Rated current	24.5 A
Pole pair ( $n_{p}$ )	4
Stator Resistance ( $R$ )	0.203 Ω
Inductance ( $L_{s}$ )	2.1 mH
Inertia ( $J$ )	0.0048 kg∙m²
Permanent magnet flux ( $ψ_{f}$ )	0.123 Wb
Viscous friction coefficient (B)	0.001 Nm∙s/rad

Table 3. Comparison between existing strategies.

Item	[10]	[28]	Proposed FS-ISMPC
Modulator	No	No	No
Cost function	${‖e‖}^{2}$	$σ^{T} \dot{σ}$ $σ = e$	$σ^{T} \dot{σ}$ $σ = e + η \int e d t$
Multi-objective optimization	Yes	Yes	Yes
Coefficient/WF	No	No	$λ, η$
Robustness	Low	Moderate	High
Constraint handling	Yes	Yes	Yes

Table 4. Hardware resources employed in control strategies.

Resources	Total	FS-MPC [10]	FS-SMPC [28]	Proposed
Total slides	13,300	4945 37.2%	3618 27.2%	4913 36.9%
Slide registers	106,400	11,225 10.5%	9024 8.5%	11,181 10.5%
Slide LUTs	53,200	12,824 24.1%	9568 18%	12,644 23.8%
Block RAMs	140	4 2.9%	4 2.9%	4 2.9%
DSP48s	220	108 49.1%	8 3.6%	90 40.9%

Table 5. Computational burden comparison.

Strategy	Clock Cycles	Execution Time
FS-MPC [10]	212	5.30 μs
FS-SMPC [28]	137	3.43 μs
Proposed	194	4.85 μs

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Hidalgo, H.; Orosco, R.; Huerta, H.; Vazquez, N.; Estrada, L.; Pinto, S.; de Castro, A. A Finite-Set Integral Sliding Modes Predictive Control for a Permanent Magnet Synchronous Motor Drive System. World Electr. Veh. J. 2024, 15, 277. https://doi.org/10.3390/wevj15070277

AMA Style

Hidalgo H, Orosco R, Huerta H, Vazquez N, Estrada L, Pinto S, de Castro A. A Finite-Set Integral Sliding Modes Predictive Control for a Permanent Magnet Synchronous Motor Drive System. World Electric Vehicle Journal. 2024; 15(7):277. https://doi.org/10.3390/wevj15070277

Chicago/Turabian Style

Hidalgo, Hector, Rodolfo Orosco, Hector Huerta, Nimrod Vazquez, Leonel Estrada, Sergio Pinto, and Angel de Castro. 2024. "A Finite-Set Integral Sliding Modes Predictive Control for a Permanent Magnet Synchronous Motor Drive System" World Electric Vehicle Journal 15, no. 7: 277. https://doi.org/10.3390/wevj15070277

Article Menu

A Finite-Set Integral Sliding Modes Predictive Control for a Permanent Magnet Synchronous Motor Drive System

Abstract

1. Introduction

2. Mathematical Model of PMSM and FS-Strategies

2.1. The Mathematical Model of PMSM

2.2. Traditional FS-MPC Strategy

2.3. FS-SMPC Strategy

3. Proposed FS-ISMPC

3.1. Integral Sliding Mode Predictive Control

3.2. Stability Analysis

3.3. Delay Compensation

3.4. Extended Control Set

4. Performance and Evaluation

4.1. Simulation Test

4.2. HIL Simulation

4.3. Real-Time Simulation Results

5. Discussion

6. Conclusions and Future Works

Author Contributions

Funding

Data Availability Statement

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI