TIDUF67 Design guide

TIDUF67 April 2024 – December 2024

3.1.3.1.1 Design of ESMO for PMSM

The conventional PLL integrated into the SMO is shown in Figure 3-8.

TIDM-02018 Block Diagram of eSMO with PLL for a PMSM

Figure 3-8 Block Diagram of eSMO with PLL for a PMSM

The traditional reduced-order sliding mode observer is constructed, which mathematical model is shown in Equation 13 and the block diagram is shown in Figure 3-9.

Equation 13.

[\begin{matrix} {\dot{\hat{i}}}_{α} \\ {\dot{\hat{i}}}_{β} \end{matrix}] = \frac{1}{L_{d}} [\begin{matrix} {- R}_{s} & - {\hat{ω}}_{e} (L_{d} - L_{q}) \\ {\hat{ω}}_{e} (L_{d} - L_{q}) & {- R}_{s} \end{matrix}] [\begin{matrix} {\hat{i}}_{α} \\ {\hat{i}}_{β} \end{matrix}] + \frac{1}{L_{d}} [\begin{matrix} V_{α} - {\hat{e}}_{α} + z_{α} \\ V_{β} - {\hat{e}}_{β} + z_{β} \end{matrix}]

where z_α and z_β are sliding mode feedback components and are defined as:

Equation 14.

[\begin{matrix} z_{α} \\ z_{β} \end{matrix}] = [\begin{matrix} k_{α} sign ({\hat{i}}_{α} - i_{α}) \\ k_{β} sign ({\hat{i}}_{β} - i_{β}) \end{matrix}]

Where k_α and k_β are the constant sliding mode gain designed by Lyapunov stability analysis. If k_α and k_β are positive and significant enough to maintain the stable operation of the SMO, the k_α and k_β are usually large enough to hold Equation 15 and Equation 16.

Equation 15.

k_{α} > m a x (|e_{α}|)

Equation 16.

k_{β} > m a x (|e_{β}|)

TIDM-02018 Block Diagram of Traditional Sliding Mode Observer

Figure 3-9 Block Diagram of Traditional Sliding Mode Observer

The estimated value of EEMF in α-β axes (Equation 18, Equation 19) can be obtained by low-pass filter from the discontinuous switching signals z_α and z_β:

Equation 17.

[\begin{matrix} {\hat{e}}_{α} \\ {\hat{e}}_{β} \end{matrix}] = \frac{ω_{c}}{s + ω_{c}} [\begin{matrix} z_{α} \\ z_{β} \end{matrix}]

Equation 18.

ê_{α}

Equation 19.

ê_{β}

Where Equation 20 is the cutoff angular frequency of the LPF, which is usually selected according to the fundamental frequency of the stator current.

Equation 20.

ω_{c} = 2 {πf}_{c}

Therefore, the rotor position can be directly calculated from arc-tangent the back EMF, defined as follows

Equation 21.

{\hat{θ}}_{e} = - \tan^{- 1} (\frac{{\hat{e}}_{α}}{{\hat{e}}_{β}})

Low pass filter removes the high-frequency term of the sliding mode function, which leads to occur phase delay resulting. This can be compensated by the relationship between the cut-off frequency ω_c and back EMF frequency ω_e, which is defined as:

Equation 22.

∆ θ_{e} = - \tan^{- 1} (\frac{ω_{e}}{ω_{c}})

And then the estimated rotor position by using SMO method is:

Equation 23.

{\hat{θ}}_{e} = - \tan^{- 1} (\frac{{\hat{e}}_{α}}{{\hat{e}}_{β}}) + ∆ θ_{e}

In a digital control application, a time discrete equation of the SMO is needed. The Euler method is the appropriate way to transform to a time discrete observer. The time discrete system matrix of Equation 13 in α-β coordinates is given by Equation 24 as:

Equation 24.

[\begin{matrix} {\hat{\dot{i}}}_{α} (n + 1) \\ {\hat{\dot{i}}}_{β} (n + 1) \end{matrix}] = [\begin{matrix} F_{α} \\ F_{β} \end{matrix}] [\begin{matrix} {\hat{\dot{i}}}_{α} (n) \\ {\hat{\dot{i}}}_{β} (n) \end{matrix}] + [\begin{matrix} G_{α} \\ G_{β} \end{matrix}] [\begin{matrix} V_{α}^{*} (n) - {\hat{e}}_{α} (n) + z_{α} (n) \\ V_{β}^{*} (n) - {\hat{e}}_{β} (n) + z_{β} (n) \end{matrix}]

Where the matrix [F] and [G] are given by Equation 25 and Equation 26 as:

Equation 25.

[\begin{matrix} F_{α} \\ F_{β} \end{matrix}] = [\begin{matrix} e^{- \frac{R_{s}}{L_{d}}} \\ e^{- \frac{R_{s}}{L_{q}}} \end{matrix}]

Equation 26.

[\begin{matrix} G_{α} \\ G_{β} \end{matrix}] = \frac{1}{R_{s}} [\begin{matrix} {1 - e}^{- \frac{R_{s}}{L_{d}}} \\ 1 - e^{- \frac{R_{s}}{L_{q}}} \end{matrix}]

The time discrete form of Equation 17 is given by Equation 27 as:

Equation 27.

[\begin{matrix} {\hat{e}}_{α} (n + 1) \\ {\hat{e}}_{β} (n + 1) \end{matrix}] = [\begin{matrix} {\hat{e}}_{α} (n) \\ {\hat{e}}_{β} (n) \end{matrix}] + 2 {πf}_{c} [\begin{matrix} z_{α} (n) - {\hat{e}}_{α} (n) \\ z_{β} (n) - {\hat{e}}_{β} (n) \end{matrix}]