Design of Lumped Disturbance Observer in Current Loop of IPMSM Based on Recursive Integral Sliding Mode Surface

Teng, Yongxiao; Gao, Qiang; Chen, Xuehan; Xu, Dianguo

doi:10.3390/en17040836

Open AccessArticle

Design of Lumped Disturbance Observer in Current Loop of IPMSM Based on Recursive Integral Sliding Mode Surface

by

Yongxiao Teng

,

Qiang Gao

^*,

Xuehan Chen

and

Dianguo Xu

School of Electrical Engineering and Automation, Harbin Institute of Technology, Harbin 150001, China

^*

Author to whom correspondence should be addressed.

Energies 2024, 17(4), 836; https://doi.org/10.3390/en17040836

Submission received: 10 January 2024 / Revised: 4 February 2024 / Accepted: 8 February 2024 / Published: 9 February 2024

(This article belongs to the Special Issue Advanced Modeling and Optimization of Electrical Drives Technology)

Download

Browse Figures

Versions Notes

Abstract

:

To overcome the problem of current control effect being reduced by unideal factors in a motor control system, such as motor parameter variation, inverter dead time, nonlinearity of the system, etc., a sliding mode disturbance observer for an interior permanent magnet synchronous motor is proposed in this paper. The model of an interior permanent magnet synchronous motor with unideal factors is designed, and the unideal factors are unified into lumped disturbances of motor stator voltage. Then, the observer for lumped disturbance is designed. A recursive integral sliding surface is used to replace the terminal sliding surface to avoid the noise sensitivity and singularity problem of the traditional terminal sliding mode observer. The observer can estimate the lumped disturbance of the current loop without relying on the accurate system model in finite time. Moreover, the structure of the current loop does not need to be adjusted while using the observer to observe and compensate for disturbances. Experiments are carried out to verify the effectiveness of the proposed observer.

Keywords:

disturbance estimate; sliding mode disturbance observer; recursive integral sliding surface; interior permanent magnet synchronous motor

1. Introduction

The interior permanent magnet synchronous motor (IPMSM) has been widely utilized due to its simple structure, high efficiency, high power density, fast response and other excellent features [1,2]. Among the IPMSM control methods, the field-oriented control (FOC) is most widely applied due to the features of low torque ripple, fast response and high control precision. At the same time, the double-loop structure of FOC can control the torque and magnetic flux independently [3,4,5]. The inner current loop determines the dynamic and steady-state performance of stator currents [6,7]. However, IPMSM is a nonlinear and strong coupling system and contains unavoidable and unmeasured disturbances such as possible parameter variations in operating [8,9,10]. Moreover, the inverter includes nonlinearity such as dead time, voltage drop and other factors [11,12]. These unideal factors, which are usually ignored in motor control system modeling, will affect the current control accuracy and cause controller parameter mismatch. It is difficult to meet the high performance requirements of the IPMSM control system using traditional control methods.

With the development of power electronics and digital signal processing, many methods to improve current control accuracy have been proposed. In general, these methods can be divided into two groups. One is to improve the current control performance by introducing new control laws. Study [13] presents a method based on an adaptive linear neuron to suppress the current harmonics by a self-tuning voltage feedforward, which does not need to extract the current harmonics. Study [14] proposes a model predictive current control (MPCC) method with adaptive-adjusting timescales based on the motor operation conditions, which can improve the system dynamic response and the current prediction accuracy. Study [15] presents a fixed switching period sliding mode control (SMC) method which combines the advantages of FOC and direct torque control (DTC). This method can allow the control law to be transferred to the inverter directly by a novel decoupling algorithm to find a set of decoupled sliding surfaces. Study [16] analyzes the relationship between the low-pass filter in the feedback loop and the controller parameters, and proposes a modified model compensation active disturbance rejection control (ADRC) method, which can make the current response closer to that of the ideal closed-loop control system. Study [17] proposes a complex vector decoupling method considering the time delay caused by sampling, calculation and PWM modulation, which quantifies and compensates the phase delay to improve the dynamic performance.

Although the above methods improved the current control performance by introducing new control laws, the parameters’ perturbation and other unideal factors are ignored, and the system disturbances have not been eliminated yet.

The other methods obtained system information through sensors or observers and improve the controller [18,19,20,21,22,23,24,25,26]. Among them, disturbance observer (DO) gains the favor of many scholars and becomes a research hotspot because of its simple structure, lack of dependence on system model accuracy and lack of an additional hardware cost. At the same time, compensation based on the disturbance observer does not change the original structure of the controller. In [21], a multiple-models adaptive disturbance observer is proposed to solve the steady-state current error by estimating the lumped disturbances, which can improve the robustness to parameter perturbations and model mismatch. Study [22] proposes an adaptive observer based on the model reference adaptive system (MRAS), which improves the control performance, anti-interference ability and robustness of the control system. Study [23] presents a novel super-twisting algorithm associated with extended state observer (ESO) for FOC, which can estimate system disturbances, eliminate chattering, improve the system anti-interference ability, and obtain a fast dynamic response. Study [24] introduces an adaptive internal model observer (AIMO) to observe system disturbances quickly without static difference, which avoid disturbances caused by parameter variations. Study [25] proposes a sliding-mode observer based on adaptive super-twisting algorithm (ASTA-SMO) to improve the performance of deadbeat predictive current control (DPCC), which can predict future values of stator current and track system disturbance caused by parameter mismatch. Study [26] presents an adaptive stator current disturbance observer considering the discrete characteristics of the current prediction error, which structures a parallel strategy, observes the discrete disturbance related to each voltage vector separately, and then eliminates the prediction error.

The above disturbance observation methods can effectively estimate and compensate for system disturbances. Among them, the SMO has simple structure and is easy to implement. Other methods have disadvantages such as difficulty with observer parameter adjustment, a poor observation effect on nonlinear disturbances, etc. However, due to the discontinuous switching characteristics, the ordinary sliding mode method has a chattering problem [27]. The terminal sliding mode method introduces nonlinear terms into the sliding mode surface, which can make the observation error converge to zero in finite time, and its dynamic performance is better than the ordinary sliding mode method with state asymptotic convergence under linear sliding mode surface conditions [28]. And the terminal sliding mode method does not have switching items, which can effectively eliminate chattering. However, the introduction of nonlinear terms results in the existence of negative exponential term in the sliding mode control function which will bring about singularity problems [29].

In this research, a terminal sliding mode disturbance observer (SMDO) based on a recursive integral sliding mode surface for lumped disturbances of IPMSM is proposed. The observer can realize the lumped estimation of the electromagnetic parameter variation and external disturbances of the control system without relying on the accurate system model in finite time. At the same time, the noise sensitivity caused by the forward Euler derivative and the singularity problem of the traditional terminal sliding mode observer are avoided. In this paper, a novel IPMSM model is presented, and all the system disturbances as the lumped stator voltage variation are considered in Section 2. The recursive integral sliding mode disturbance observer is analyzed and designed in Section 3. The experimental results are given to verify the feasibility of the proposed observer in Section 4.

2. Modeling of IPMSM with Parameter Variation and Unknown Disturbance

According to study [2], the state equation of the IPMSM dynamic model in a synchronously rotating reference frame can be denoted as (1):

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

(1)

where u_d, u_q, i_d and i_q are the stator voltage and current components in d- and q-axis, respectively; R_s is the nominal value of stator resistance; L_d and L_q are the nominal values of the inductance in the d- and q-axis, respectively; and ψ_f is the nominal value of permanent magnet flux linkage.

The system state variables are defined as x = [i_d i_q]^T, the input variables are defined as u = [u_d u_q]^T, then (1) can be expressed as (2):

\dot{x} = A x + B u + C

(2)

where

A = [\begin{matrix} - \frac{R_{s}}{L_{d}} & \frac{ω_{e} L_{q}}{L_{d}} \\ - \frac{ω_{e} L_{d}}{L_{q}} & - \frac{R_{s}}{L_{q}} \end{matrix}], B = [\begin{matrix} \frac{1}{L_{d}} & 0 \\ 0 & \frac{1}{L_{q}} \end{matrix}], C = [\begin{matrix} 0 \\ - \frac{ω_{e} ψ_{f}}{L_{q}} \end{matrix}] .

The motor parameters including stator resistance, stator inductance, and permanent magnet flux linkage vary with the varying operating condition, then the controller performance can be affected by the varying motor parameters. Meanwhile, the system nonlinearity, dead time and other factors of the inverter can also affect the motor voltage and current. After considering the above, the system state equation can be expressed as [2]:

\{\begin{matrix} \frac{d i_{d}}{d t} = - \frac{R_{s} + Δ R_{s}}{L_{d} + Δ L_{d}} i_{d} + \frac{L_{q} + Δ L_{q}}{L_{d} + Δ L_{d}} ω_{e} i_{q} + \frac{1}{L_{d} + Δ L_{d}} (u_{d} + Δ u_{d}) \\ \frac{d i_{q}}{d t} = - \frac{R_{s} + Δ R_{s}}{L_{q} + Δ L_{q}} i_{q} - \frac{L_{d} + Δ L_{d}}{L_{q} + Δ L_{q}} ω_{e} i_{d} + \frac{1}{L_{q} + Δ L_{q}} (u_{q} + Δ u_{q}) - \frac{ψ_{f} + Δ ψ_{f}}{L_{q} + Δ L_{q}} ω_{e} \end{matrix}

(3)

where ΔR, ΔL_d, ΔL_q, Δψ_f are variations of R_s, L_d, L_q, and ψ_f, respectively; Δu_d and Δu_q are the voltage variations caused by the other unideal factors. After considering the above issues, (2) can be obtained:

\dot{x} = (A + Δ A) x + (B + Δ B) (u + Δ u) + (C + Δ C)

(4)

From (4), the variation in current response will induce additional errors, and cause parameter mismatch of the current controller, resulting in a decrease in the current control effect.

Considering the effects above, all the disturbances can be regarded as the lumped disturbance f = [f_d f_q]^T in the d- and q-axis stator voltage.

f = B^{- 1} (Δ A x + B Δ u + u Δ B + Δ B Δ u + Δ C)

(5)

Then (4) can be expressed as:

\dot{x} = A x + B (u + f) + C

(6)

3. Design and Analysis of Recursive Integral Sliding Mode Disturbance Observer

Disturbance observer is an effective method to solve the current controller mismatch. After the lumped disturbance f is estimated by the disturbance observer, disturbance compensation controller can eliminate the disturbance influence efficaciously. The block diagram of the IPMSM control scheme with disturbance compensation is shown in Figure 1, where ω_e* is the reference value of the rotating speed speed, x* is the reference value of current,

\hat{f}

is the estimated value of the lumped disturbances f, and f_c is the disturbance compensation value obtained by

\hat{f}

.

The controlled system with the disturbance observation and compensation can be equivalent to Figure 2, then the system state equation can be shown as (7):

\dot{x} = A x + B (u + f - f_{c}) + C

(7)

According to (7), when f − f_c = 0, the effects of system disturbances can be compensated and the actual system state can be equivalent to the ideal state without disturbances. Then, the current loop parameters can be designed according to (2) under ideal conditions.

3.1. Design of Recursive Integral SMDO

The SMDO is constructed based on (6), and the observation equation is as follows:

\dot{\hat{x}} = A \hat{x} + B u + C + u_{smo}

(8)

where

\dot{\hat{x}}

= [

{\hat{i}}_{d}

{\hat{i}}_{q}

]^T,

{\hat{i}}_{d}

and

{\hat{i}}_{q}

are the observed values of i_d and i_q, respectively; u_smo = [u_smod u_smoq]^T, u_smod and u_smoq are the d- and q-axis control functions, respectively.

The current observation error e = [e_id e_iq]^T is defined as (9):

e = x - \hat{x}

(9)

By substituting (6) and (8) into (9), the state equation of the current observation error is derived as:

\dot{e} = A e + B f - u_{smo}

(10)

To ensure the observation error can converge to 0 in finite time, and the observer control functions u_smod and u_smoq are non-singular, an integral fast terminal sliding mode surface is selected as [27]:

s_{i} = [\begin{matrix} s_{id} \\ s_{iq} \end{matrix}] = [\begin{matrix} e_{id} + \int_{0}^{t} [α_{id} e_{id} + β_{id} {|e_{id}|}^{γ_{id}} sign (e_{id})] d τ \\ e_{iq} + \int_{0}^{t} [α_{iq} e_{iq} + β_{iq} {|e_{iq}|}^{γ_{iq}} sign (e_{iq})] d τ \end{matrix}]

(11)

where α_id, α_iq, β_id, β_iq, γ_id, and γ_iq are sliding mode observer coefficients; in general, α_id > 0, α_iq > 0, β_id > 0, β_iq > 0, 0 < γ_id < 1, and 0 < γ_iq < 1.

When the control functions u_smo satisfy the SMDO reachability condition and the observer is in sliding mode, the first-order derivative of the sliding mode surface is 0:

{\dot{s}}_{i} = [\begin{matrix} {\dot{s}}_{id} \\ {\dot{s}}_{iq} \end{matrix}] = [\begin{matrix} {\dot{e}}_{id} + α_{id} e_{id} + β_{id} {|e_{id}|}^{γ_{id}} sign (e_{id}) \\ {\dot{e}}_{iq} + α_{iq} e_{iq} + β_{iq} {|e_{iq}|}^{γ_{iq}} sign (e_{iq}) \end{matrix}] = [\begin{matrix} 0 \\ 0 \end{matrix}]

(12)

The equivalent observer control function u_eq can be obtained according to (10) and (12):

u_{eq} = A e + [\begin{matrix} α_{id} e_{id} + β_{id} {|e_{id}|}^{γ_{id}} sign (e_{id}) \\ α_{iq} e_{iq} + β_{iq} {|e_{iq}|}^{γ_{iq}} sign (e_{iq}) \end{matrix}]

(13)

Since there are nonlinear and unmeasurable factors in the lumped disturbances f, and the observer cannot obtain the estimated values

\hat{f}

directly by only using the equivalent control functions u_eq, it is necessary to design the disturbance tracking functions u_f = [u_fd u_fq]^T to estimate the lumped disturbances. According to (6) and (10), (14) can be obtained as follows:

\dot{\hat{x}} = A \hat{x} + B u + C + B u_{f} + u_{eq}

(14)

The state equation of error e is rewritten as:

\dot{e} = A e + B (f - u_{f}) - u_{eq}

(15)

On the basis of an integral fast terminal sliding mode surface s_i, the recursive sliding mode surface s_f is adopted to realize the estimation of lumped disturbances f, and the SMDO based on the recursive terminal sliding mode surface is obtained. The s_f is as (16):

s_{f} = [\begin{matrix} s_{fd} \\ s_{fq} \end{matrix}] = [\begin{matrix} {\dot{s}}_{id} + α_{fd} s_{id} + β_{fd} {|s_{id}|}^{γ_{fd}} sign (s_{id}) \\ {\dot{s}}_{iq} + α_{fq} s_{iq} + β_{fq} {|s_{iq}|}^{γ_{fq}} sign (s_{iq}) \end{matrix}]

(16)

where α_fd, α_fq, β_fd, β_fq, γ_fd, and γ_fq are sliding mode observer coefficients of s_f; in general, α_fd > 0, α_fq > 0, β_fd > 0, β_fq > 0, 0 < γ_fd < 1, and 0 < γ_fq < 1.

The first-order derivative of s_f is given in (17):

{\dot{s}}_{f} = [\begin{matrix} {\dot{s}}_{fd} \\ {\dot{s}}_{fq} \end{matrix}] = [\begin{matrix} {\ddot{s}}_{id} + [α_{fd} s_{id} + β_{fd} {|s_{id}|}^{γ_{fd}} sign (s_{id})]^{'} \\ {\ddot{s}}_{iq} + [α_{fq} s_{iq} + β_{fq} {|s_{iq}|}^{γ_{fq}} sign (s_{iq})]^{'} \end{matrix}]

(17)

In order to solve the disturbance tracking function u_f in SMDO, (15) is substituted into (11) to obtain the first and second order derivative of s_i:

{\dot{s}}_{i} = B (f - u_{f})

(18)

{\ddot{s}}_{i} = B (\dot{f} - {\dot{u}}_{f})

(19)

The disturbance tracking function u_f is selected as (20):

u_{f} = [\begin{matrix} L_{d} [α_{fd} s_{id} + β_{fd} {|s_{id}|}^{γ_{fd}} sign (s_{id})] \\ L_{q} [α_{fq} s_{iq} + β_{fq} {|s_{iq}|}^{γ_{fq}} sign (s_{iq})] \end{matrix}]

(20)

According to Equations (13), (14) and (20), the observer control function u_smo can be obtained as (21):

\begin{matrix} u_{smo} = u_{eq} + B u_{f} \\ = A e + [\begin{matrix} α_{id} e_{id} + β_{id} {|e_{id}|}^{γ_{id}} sign (e_{id}) \\ α_{iq} e_{iq} + β_{iq} {|e_{iq}|}^{γ_{iq}} sign (e_{iq}) \end{matrix}] + [\begin{matrix} α_{fd} s_{id} + β_{fd} {|s_{id}|}^{γ_{fd}} sign (s_{id}) \\ α_{fq} s_{iq} + β_{fq} {|s_{iq}|}^{γ_{fq}} sign (s_{iq}) \end{matrix}] \end{matrix}

(21)

When the observer is in the balance state e = 0, the equivalent observer control function u_eq in (21) is equal to zero and the disturbance tracking function u_f in (21) is equal to the lumped disturbance f, so u_f is the observed value

\hat{f}

of the lumped disturbance.

The block diagrams of the disturbance observer are shown in Figure 3.

While using the traditional terminal sliding mode method to observe the system states, the system state equation is usually in the form of an integral chain [28]. At this time, the system state equation as follows:

\{\begin{matrix} \dot{e} = e_{f} \\ {\dot{e}}_{f} = \dot{f} - u_{smoT} \end{matrix}

(22)

where e_f = f −

\hat{f}

= [e_fd e_fq]^T, and u_smoT is the sliding mode control function and is also the observed value of

\hat{f}

.

The sliding mode surface is selected as (23) according to study [29], and the recursive sliding mode surfaces is selected as (24):

s_{1} = [\begin{matrix} s_{1 d} \\ s_{1 q} \end{matrix}] = [\begin{matrix} {\dot{e}}_{id} + α_{1 d} e_{id} + β_{1 d} {e_{id}}^{γ_{1 d}} \\ {\dot{e}}_{iq} + α_{1 q} e_{iq} + β_{1 q} {e_{iq}}^{γ_{iq}} \end{matrix}]

(23)

s_{2} = [\begin{matrix} s_{fd} \\ s_{fq} \end{matrix}] = [\begin{matrix} {\dot{s}}_{1 d} + α_{2 d} s_{1 d} + β_{2 d} {s_{1 d}}^{γ_{2 d}} \\ {\dot{s}}_{1 q} + α_{2 q} s_{1 q} + β_{2 q} {s_{1 q}}^{γ_{2 q}} \end{matrix}]

(24)

where α_ij > 0, β_ij > 0, and 0 < γ_ij < 1 (i =1, 2; j = d, q) are sliding mode coefficients.

According to (22)–(24), the control function u_smoT is selected as:

u_{smoT} = [\begin{matrix} α_{1 d} e_{fd} + β_{1 d} γ_{1 d} {e_{id}}^{γ_{1 d} - 1} e_{fd} + α_{2 d} s_{1 d} + β_{2 d} {s_{1 d}}^{γ_{2 d}} \\ α_{1 q} e_{fq} + β_{1 q} γ_{1 q} {e_{iq}}^{γ_{1 q} - 1} e_{fq} + α_{2 q} s_{1 q} + β_{2 q} {s_{1 q}}^{γ_{2 q}} \end{matrix}]

(25)

According to (25), the control function u_smoT contains the current observation error e and the disturbance observation error e_f. e can be obtained directly based on the actual current of the system and observed current of the observer, but e_f needs to be obtained by derivation of e. In common, e_f is calculated as

e_{f} (k) = \dot{e} (k) \approx [e (k) - e (k - 1)] / T_{c}

, which is sensitive to noise. T_c is the control period of the current loop. At the same time, since γ_1d − 1 < 0 and γ_1q − 1 < 0, the observer has a singular problem, while e = 0 and e_f ≠ 0.

Compared with traditional terminal SMDO, the proposed recursive integral SMDO considers lumped disturbance as a motor stator voltage variation. The state equation of observer is the first-order differential equation shown in (10). The observer control function u_smo as shown in (21) only requires the current observation error e and does not contain the derivative s_i or e. And the observer control function u_smo does not contain the negative exponential term, there is no singular problem when the observer is in sliding mode. Therefore, the proposed recursive integral SMDO can avoid the noise sensitivity and singularity problem while tracking the system lumped disturbance in finite time.

The SMDO only demands the input and output information of a system, and the lumped disturbance can be eliminated by disturbance observation and compensation.

3.2. Stability Analysis

The Lyapunov function V_i is defined as (26), and its first-order differential is given as (27):

V_{i} = \frac{1}{2} {\dot{s}}_{i}^{T} {\dot{s}}_{i}

(26)

\begin{array}{l} {\dot{V}}_{i} = {\dot{s}}_{i}^{T} {\ddot{s}}_{i} & = \frac{1}{L_{d}} {\dot{s}}_{id} ({\dot{f}}_{d} - α_{fd} {\dot{s}}_{id} - β_{fd} γ_{fd} {|s_{id}|}^{γ_{fd} - 1} {\dot{s}}_{id}) + \frac{1}{L_{q}} {\dot{s}}_{iq} ({\dot{f}}_{q} - α_{fq} {\dot{s}}_{iq} - β_{fq} γ_{fq} {|s_{iq}|}^{γ_{fq} - 1} {\dot{s}}_{iq}) \\ \leq \frac{1}{L_{d}} |{\dot{s}}_{id}| (|{\dot{f}}_{d}| - α_{fd} |{\dot{s}}_{id}| - β_{fd} γ_{fd} {|s_{id}|}^{γ_{fd} - 1} |{\dot{s}}_{id}|) + \frac{1}{L_{d}} |{\dot{s}}_{id}| (|{\dot{f}}_{d}| - α_{fd} |{\dot{s}}_{id}| - β_{fd} γ_{fd} {|s_{id}|}^{γ_{fd} - 1} |{\dot{s}}_{id}|) \end{array}

(27)

When the control period of the current loop is short enough,

{\dot{f}}_{d}

and

{\dot{f}}_{q}

can be considered as zero.

{\dot{V}}_{i}

< 0 is drawn while the observer parameter ranges are α_fd > 0, α_fq > 0, β_fd > 0, β_fq > 0, 1 > γ_fd > 0, and 1 > γ_fq > 0. It can be derived that

{\dot{s}}_{i}

is asymptotically stable at itsorigin, which implies that the current observation error e can converge to 0 in finite time with the control function u_smo.

Define the disturbance observation error e_f = f − u_f, and define the Lyapunov function V_f as follows:

V_{f} = \frac{1}{2} e_{f}^{T} e_{f}

(28)

According to (18) and (19), the first-order derivative of V_f is given as:

{\dot{V}}_{f} = e_{f}^{T} {\dot{e}}_{f} = {(B^{- 1} {\dot{s}}_{i})}^{T} (B^{- 1} {\ddot{s}}_{i})

(29)

It can be derived that e_f is also asymptotically stable at its origin. Thus, the current observation results

\dot{\hat{x}}

of SMDO can converge to the actual currents, and the disturbance tracking functions u_f converge to the lumped disturbance f in finite time.

The trajectory of the observation error e is denoted as (30):

\{\begin{matrix} {\dot{e}}_{id} + α_{id} e_{id} + β_{id} {|e_{id}|}^{γ_{id}} sign (e_{id}) = 0 \\ {\dot{e}}_{iq} + α_{iq} e_{iq} + β_{iq} {|e_{iq}|}^{γ_{iq}} sign (e_{iq}) = 0 \end{matrix}

(30)

When the current observation error is close to the balance state e = 0, the convergence speed is mainly determined by (31):

\{\begin{matrix} {\dot{e}}_{id} = - β_{id} {|e_{id}|}^{γ_{id}} sign (e_{id}) \\ {\dot{e}}_{iq} = - β_{iq} {|e_{iq}|}^{γ_{iq}} sign (e_{iq}) \end{matrix}

(31)

When the current observation error is far from its origin, the convergence speed is mainly determined by (32):

\{\begin{matrix} {\dot{e}}_{id} = - α_{id} e_{id} \\ {\dot{e}}_{iq} = - α_{iq} e_{iq} \end{matrix}

(32)

By solving (30), the time to arrive at e = 0 can be derived as:

\{\begin{matrix} t_{id} = \frac{1}{α_{id} (1 - γ_{id})} \ln \frac{α_{id} {|e_{id} (0)|}^{1 - γ_{id}} + β_{id}}{β_{id}} \\ t_{iq} = \frac{1}{α_{iq} (1 - γ_{iq})} \ln \frac{α_{iq} {|e_{iq} (0)|}^{1 - γ_{iq}} + β_{iq}}{β_{iq}} \end{matrix}

(33)

4. Experimental Results

An experimental setup was built to validate the proposed algorithm, as shown in Figure 4. An IPMSM, an asynchronous induction motor (ACIM) and a torque sensor are connected via diaphragm couplings, where the IPMSM is used to test the algorithm proposed in this paper, and the ACIM is used to provide the load torque. The parameters of IPMSM are listed in Table 1.

The control algorithm of IPMSM is implemented through DSP TI TMS320F28377D. The CPU1 of DSP communicates with the PC while running the basic motor control algorithm. The CPU2 is responsible for realizing the disturbance observer. The interfaces of multiple peripherals including resolver-to-digital converter (R/D Converter), analog–digital converter (ADC), and digital–analog converter (DAC), the external I/O interface, and the hardware protection logic are implemented in FPGA Xilinx XC6SLX16. Hall sensors are used to obtain voltage and current information of IPMSM. A rotating transformer is used to measure the speed and position of IPMSM. The pulse width modulation (PWM) switching frequency is 20 kHz, and the DC bus voltage is 200 V. An oscilloscope (DLM5054), a recorder (DL850) and a power analyzer (WT5000) manufactured by Yokogawa were used to record the system data during the experiments. The experimental setup is shown in Figure 5.

Experiments are carried out to validate the proposed recursive integral SMDO. During the experiments, the same parameters are applied to the speed loop and current loop. After IPMSM reaches 1400 rpm, the load torque increases from 0 N·m to 20 N·m, then decreases to 0 N·m.

Figure 6 shows the experimental results of currents without disturbance compensation. In Figure 6a,b, the top waveforms in each figure are the currents i_dn and i_qn, calculated by DSP based on the current motor currents i_d and i_q and the voltage given values u_d and u_q of the current loop through the state Equation (1) under ideal conditions. The middle waveforms in each figure are the actual i_d and i_q obtained by park transformation of the measured phase currents. The waveforms at the bottom of each figure are the current observation values

{\hat{i}}_{d}

and

{\hat{i}}_{q}

of the observer. i_d varies between −1.4 A and −5.7 A, and i_q varies between 34.2 A and 38.2 A when speed is 1400 rpm and torque is 20 N·m.

Due to the disturbances shown in (3), the actual currents i_d and i_q are not equal to the theoretical currents i_dn and i_qn calculated by DSP. The difference between the actual currents and the theoretical currents error_id= i_dn − i_d and error_iq = i_qn − i_q are shown in Figure 7. The absolute values of error_id and error_iq increase with speed and load increasing. The error_id varies between −3.8 A and −1.8 A, and the error_iq varies between −0.2 A and −0.9 A when speed is 1400 rpm and torque is 20 N·m.

Figure 8 shows the difference between the observed currents and the actual currents, that is, the current observation error e in (9). Although the fluctuation range of observation error increases with the increase in speed and load, it still remains within a small range. Both e_id and e_iq fluctuate between −0.2 and 0.2.

Figure 9 shows the derivative of the current observation error calculated by DSP based on the formula

\dot{e} (k) \approx [e (k) - e (k - 1)] / T_{c}

. The derivation based on the backward Euler method will amplify the quantization errors and noise in the error e, and also bring about truncation error. These factors will affect the accuracy of disturbance observation while using the traditional terminal SMDO.

The observed lumped disturbance

\hat{f}

is shown in Figure 10. The absolute values of observed disturbance values also increase with speed and load increasing.

{\hat{f}}_{d}

varies between 8.1 V and 13.4 V, and

{\hat{f}}_{q}

varies between 2.2 V and 5.5 V when speed is 1400 rpm and torque is 20 N·m.

Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15 show the experimental results when the observed lumped disturbance

\hat{f}

is utilized in the system.

Figure 11 shows the currents with the disturbance compensation. i_d varies between −1.8 A and −5.2 A, and i_q varies between −34.6 A and 38.1 A under 1400 rpm and 20 N·m. Compared with the current waveforms in Figure 6, the fluctuations in the motor currents i_d and i_q become smaller.

error_id and error_iq after adding the disturbance compensation are shown in Figure 12, which indicates that the absolute values of steady-state differences no longer increase with the increase in speed and load, but remain close to zero. The error_id varies between −0.9 A and 0.8 A, and the error_iq varies between −0.3 A and −0.3 A under 1400 rpm and 20 N·m.

Figure 13 shows the current observation error e. e_id and e_iq still fluctuate between −0.2 and 0.2. The introduction of the disturbance compensation does not affect the accuracy of current observation.

The observed lumped disturbance value

\hat{f}

with the disturbance compensation is shown in Figure 14. At this time,

\hat{f}

is no longer the observed values of f, but f − f_c in (7), that is, the remaining disturbances after compensation. The steady-state of

\hat{f}

also become close to zero.

{\hat{f}}_{d}

varies between −2.2 V and 2.1 V, and

{\hat{f}}_{d}

varies between −1.2 V and 1.2 V when speed is 1400 rpm and torque is 20 N·m.

Figure 15 shows the compensation value f_c of the lumped disturbance. The absolute value of f_c also increases with the increase in speed and load torque. f_cd varies between 9.8 V and 12.0 V, and f_cd varies between 3.1 V and 4.5 V when speed is 1400 rpm and torque is 20 N·m. Compared with Figure 10, f_c is close to the observed value

\hat{f}

when the system is not compensated.

Table 2 shows the comparison of the peak-to-peak values of system states when the lumped disturbances are compensated and uncompensated under 1400 rpm and 20 N·m where the current fluctuation and control error are at the maximum. It can be seen from Table 2 that the fluctuation in the system state becomes smaller after compensation, which means that the current loop control precision is higher and the motor output torque is more stable.

5. Conclusions and Future Studies

This article proposed a sliding mode disturbance observer for IPMSM. An IPMSM model which unifies the unideal factors in the current loop into the lumped disturbance of the motor stator voltage is considered. Furthermore, a terminal SMDO based on the recursive integral sliding mode surface is used to observe the lumped disturbances of the system. The integral sliding mode surface avoids the noise sensitivity and singularity problems of the traditional terminal sliding mode surface, and the introduction of the recursive sliding mode surface enables the estimation of lumped disturbances. The stability analysis shows that the recursive integral SMDO proposed in this article can realize the estimated values of currents and lumped disturbances to converge to the real values in finite time. Compensation based on observation values of lumped disturbance can improve current control accuracy and eliminate controller current loop parameter mismatch without changing the structure of the current loop. The experimental results are presented to validate the proposed observer. With the introduction of disturbance observation and compensation, the accuracy current control is improved, the current control errors are maintained near 0, and the fluctuations in error are reduced by about 15%. At the same time, the fluctuation in current in the d-axis is reduced by 20.9%, and the fluctuation in current in the q-axis is reduced by 12.5%. The proposed terminal sliding mode SMDO only considers the voltage equation of IPMSM to observe lumped disturbance of the current loop, and considers the motor speed to be a constant rather than a variable. Future work will combine the IPMSM torque equation to estimate the disturbance of speed loop disturbance while observing the lumped disturbance of the current loop, and improve the performance of the IPMSM control system by compensating for disturbances in the current loop and speed loop.

Author Contributions

Methodology, Y.T.; Software, X.C.; Resources, Q.G.; Supervision, D.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

Wu, G.; Huang, S.; Wu, Q.; Rong, F.; Zhang, C.; Liao, W. Robust Predictive Torque Control of N*3-Phase PMSM for High-Power Traction Application. IEEE Trans. Power Electron. 2020, 35, 10799–10809. [Google Scholar] [CrossRef]
Yang, J.; Chen, W.-H.; Li, S.; Guo, L.; Yan, Y. Disturbance/Uncertainty Estimation and Attenuation Techniques in PMSM Drives—A Survey. IEEE Trans. Ind. Electron. 2017, 64, 3273–3285. [Google Scholar] [CrossRef]
Rauth, S.S.; Samanta, B. Comparative Analysis of IM/BLDC/PMSM Drives for Electric Vehicle Traction Applications Using ANN-Based FOC. In Proceedings of the 2020 IEEE 17th India Council International Conference (INDICON), New Delhi, India, 10–13 December 2020; pp. 1–8. [Google Scholar]
Gupta, S.; George, S.; Awate, V. PI Controller Design and Application for SVPWM Switching Technique Based FOC of PMSM. In Proceedings of the 2023 Second International Conference on Trends in Electrical, Electronics, and Computer Engineering (TEECCON), Bangalore, India, 26–27 August 2023; pp. 172–177. [Google Scholar]
Zhang, Y.; Shen, W.; Yang, H. An Improved Deadbeat Predictive Current Control of PMSM Drives Based on the Ultra-Local Model. Chin. J. Electr. Eng. 2023, 9, 27–37. [Google Scholar] [CrossRef]
Xu, C.; Han, Z.; Lu, S. Deadbeat Predictive Current Control for Permanent Magnet Synchronous Machines with Closed-Form Error Compensation. IEEE Trans. Power Electron. 2020, 35, 5018–5030. [Google Scholar] [CrossRef]
Yang, R.; Sun, T.; Feng, W.; He, S.; Zhu, S.; Chen, X. Accurate Online MTPA Control of IPMSM Considering Derivative Terms. Chin. J. Electr. Eng. 2021, 7, 100–110. [Google Scholar] [CrossRef]
Agoro, S.; Husain, I. High-Fidelity Nonlinear Modeling of an Asymmetrical Dual Three-Phase PMSM. In Proceedings of the 2023 IEEE International Electric Machines & Drives Conference (IEMDC), San Francisco, CA, USA, 15–18 May 2023; pp. 1–7. [Google Scholar]
Li, Z.; Li, H. MTPA Control of PMSM System Considering Saturation and Cross-Coupling. In Proceedings of the 2012 15th International Conference on Electrical Machines and Systems (ICEMS), Sapporo, Japan, 21–24 October 2012; pp. 1–5. [Google Scholar]
Xia, J.; Li, Z.; Yu, D.; Guo, Y.; Zhang, X. Robust Speed and Current Control with Parametric Adaptation for Surface-Mounted PMSM Considering System Perturbations. IEEE J. Emerg. Sel. Top. Power Electron. 2021, 9, 2807–2817. [Google Scholar] [CrossRef]
Zhou, G.; Wang, B.; Xu, K.; Liu, Q. An Inverter Nonlinearity Compensation Method for PMSM Drives Based on Dead-Time Elimination. In Proceedings of the 2023 26th International Conference on Electrical Machines and Systems (ICEMS), Zhuhai, China, 5–8 November 2023; pp. 662–667. [Google Scholar]
Zhu, S.; Huang, W.; Yan, Y. A Simple Inverter Nonlinearity Compensation Method Using On-Line Voltage Error Observer. In Proceedings of the 2019 22nd International Conference on Electrical Machines and Systems (ICEMS), Harbin, China, 11 August 2019; pp. 1–6. [Google Scholar]
Wang, L.; Zhu, Z.Q.; Bin, H.; Gong, L.M. Current Harmonics Suppression Strategy for PMSM with Nonsinusoidal Back-EMF Based on Adaptive Linear Neuron Method. IEEE Trans. Ind. Electron. 2020, 67, 9164–9173. [Google Scholar] [CrossRef]
Niu, F.; Chen, X.; Huang, S.; Huang, X.; Wu, L.; Li, K.; Fang, Y. Model Predictive Current Control with Adaptive-Adjusting Timescales for PMSMs. CES Trans. Electr. Mach. Syst. 2021, 5, 108–117. [Google Scholar] [CrossRef]
Repecho, V.; Biel, D.; Arias, A. Fixed Switching Period Discrete-Time Sliding Mode Current Control of a PMSM. IEEE Trans. Ind. Electron. 2018, 65, 2039–2048. [Google Scholar] [CrossRef]
Gao, Y.; Zhao, H.; Ma, K.; Huo, X. Control Design for Current Loop of PMSM Using A Modified Model-Compensation ADRC Controller. In Proceedings of the 2020 39th Chinese Control Conference (CCC), Shenyang, China, 27–29 July 2020; pp. 2545–2550. [Google Scholar]
Ding, Y.; Sun, L. Complex Vector Decoupling and Delay Compensation for PMSM at Low Carrier Ratio. In Proceedings of the 2023 IEEE 6th International Electrical and Energy Conference (CIEEC), Hefei, China, 12–14 May 2023; pp. 1025–1031. [Google Scholar]
Ma, X.; Bi, C. A Technology for Online Parameter Identification of Permanent Magnet Synchronous Motor. CES Trans. Electr. Mach. Syst. 2020, 4, 237–242. [Google Scholar] [CrossRef]
Liang, G.; Huang, S.; Li, M.; Wu, X. A high-order fast terminal sliding-mode disturbance observer based on mechanical parameter identification for PMSM. Trans. China Electrotech. Soc. 2020, 35, 395–403. [Google Scholar] [CrossRef]
Liu, Y.; Wang, X.; Zhou, K. Current deviation decoupling control with a sliding mode observer for permanent magnet synchronous motor. Trans. China Electrotech. Soc. 2020, 35, 1642–1652. [Google Scholar] [CrossRef]
Bai, C.; Yin, Z.; Zhang, Y.; Liu, J. Multiple-Models Adaptive Disturbance Observer-Based Predictive Control for Linear Permanent-Magnet Synchronous Motor Vector Drive. IEEE Trans. Power Electron. 2022, 37, 9596–9611. [Google Scholar] [CrossRef]
Li, Z.; Zhang, Z.; Feng, S.; Wang, J.; Guo, X.; Sun, H. Design of Model-Free Speed Regulation System for Permanent Magnet Synchronous Linear Motor Based on Adaptive Observer. IEEE Access 2022, 10, 68545–68556. [Google Scholar] [CrossRef]
Le, T.-L.; Hsieh, M.-F.; Nguyen, M.-T. Robust Speed Control with Disturbance Rejection for Surface-Mounted Permanent Magnet Synchronous Motor. In Proceedings of the 2022 6th International Conference on Green Technology and Sustainable Development (GTSD), Hefei, China, 12–14 July 2022; pp. 459–463. [Google Scholar]
Chen, M.; Wang, F.; He, L.; Ke, D.; Zuo, K.; Rodriguez, J. Predictive Current Control of Permanent Magnet Synchronous Motor Based on An Adaptive Internal Model Observer. In Proceedings of the 2020 IEEE 9th International Power Electronics and Motion Control Conference (IPEMC2020-ECCE Asia), Nanjing, China, 29 November–2 December 2020; pp. 2567–2571. [Google Scholar]
Zhao, F.; Lian, G.; Ban, F.; Hou, Z. Improved Deadbeat Predictive Current Control Combined Adaptive Super-Twisting Algorithm Based Sliding-Mode Observer for Permanent-Magnet Synchronous Motors Drives. In Proceedings of the 2020 IEEE 5th Information Technology and Mechatronics Engineering Conference (ITOEC), Chongqing, China, 12–14 June 2020; pp. 1749–1754. [Google Scholar]
Wang, F.; Zuo, K.; Lin, G.; He, L.; Rodriguez, J.; Garcia, C. Adaptive Stator Current Disturbance Observer Based on the Predictive Current Control for PMSM. In Proceedings of the IECON 2020 The 46th Annual Conference of the IEEE Industrial Electronics Society, Singapore, 18–21 October 2020; pp. 4612–4618. [Google Scholar]
Jie, Y.; Qing-Lin, W.; Yuan, L. Sliding Mode Variable Structure Control Theory: A Survey. In Proceedings of the 31st Chinese Control Conference, Hefei, China, 25–27 July 2012; pp. 3197–3202. [Google Scholar]
Du, X.; Fang, X.; Liu, F. Continuous Full-Order Nonsingular Terminal Sliding Mode Control for Systems with Matched and Mismatched Disturbances. IEEE Access 2019, 7, 130970–130976. [Google Scholar] [CrossRef]
Zhao, Y. A Novel Piece-Wise Non-Singular Fast Terminal Sliding Mode Control for Second-Order Nonlinear System. In Proceedings of the 2022 World Automation Congress (WAC), San Antonio, TX, USA, 11–15 October 2022; pp. 18–22. [Google Scholar]

Figure 1. Control scheme of IPMSM with disturbance compensation.

Figure 2. Equivalent model of controlled system.

Figure 3. Block diagrams of SMDO (a) traditional terminal SMDO (b) recursive integral SMDO.

Figure 4. Block diagram of test platform.

Figure 5. Experimental setup.

Figure 6. Waveforms of currents without disturbance compensation. (a) d-axis currents (b) q-axis currents.

Figure 7. error_id and error_iq without disturbance compensation.

Figure 8. Current observation error without disturbance compensation.

Figure 9. Derivative of current observation error.

Figure 10. Observation values of lumped disturbance.

Figure 11. Waveforms of currents with disturbance compensation. (a) d-axis currents (b) q-axis currents.

Figure 12. error_id and error_iq with disturbance compensation.

Figure 13. Current observation error with disturbance compensation.

Figure 14. Observation values of remaining lumped disturbances.

Figure 15. Disturbance compensation values.

Table 1. Parameters of IPMSM.

Descriptions	Parameters	Values
Stator resistance	R_s	12.4 mΩ
d-axis inductance	L_d	190 μH
q-axis inductance	L_q	400 μH
Flux	ψ_f	0.0712 Wb
Pole pairs	p	6
Inertia	J	0.09615 kg·m²

Table 2. Comparison of system states.

System States	Peak-to-Peak Value without Compensation	Peak-to-Peak Value with Compensation	Reduced Proportion
i_d	4.3 A	3.4 A	20.9%
i_q	4.0 A	3.5 A	12.5%
error_id	2.0 A	1.7 A	15.0%
error_iq	0.7 A	0.6 A	14.2%
${\hat{f}}_{d}$	5.3 V	4.3 V	18.9%
${\hat{f}}_{q}$	3.3 V	2.4 V	27.3%

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

Teng, Y.; Gao, Q.; Chen, X.; Xu, D. Design of Lumped Disturbance Observer in Current Loop of IPMSM Based on Recursive Integral Sliding Mode Surface. Energies 2024, 17, 836. https://doi.org/10.3390/en17040836

AMA Style

Teng Y, Gao Q, Chen X, Xu D. Design of Lumped Disturbance Observer in Current Loop of IPMSM Based on Recursive Integral Sliding Mode Surface. Energies. 2024; 17(4):836. https://doi.org/10.3390/en17040836

Chicago/Turabian Style

Teng, Yongxiao, Qiang Gao, Xuehan Chen, and Dianguo Xu. 2024. "Design of Lumped Disturbance Observer in Current Loop of IPMSM Based on Recursive Integral Sliding Mode Surface" Energies 17, no. 4: 836. https://doi.org/10.3390/en17040836

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

Article Menu

Design of Lumped Disturbance Observer in Current Loop of IPMSM Based on Recursive Integral Sliding Mode Surface

Abstract

1. Introduction

2. Modeling of IPMSM with Parameter Variation and Unknown Disturbance

3. Design and Analysis of Recursive Integral Sliding Mode Disturbance Observer

3.1. Design of Recursive Integral SMDO

3.2. Stability Analysis

4. Experimental Results

5. Conclusions and Future Studies

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