Speed_Sensorless_Control_of_PMSM_Using_An_Improved
Speed_Sensorless_Control_of_PMSM_Using_An_Improved
net/publication/242571707
CITATIONS READS
21 523
2 authors, including:
Mongkol Konghirun
King Mongkut's University of Technology Thonburi
65 PUBLICATIONS 636 CITATIONS
SEE PROFILE
All content following this page was uploaded by Mongkol Konghirun on 02 December 2016.
z α = eα + ∆α
zβ = eβ + ∆β (4)
Where ∆α , ∆β are the errors determined by the
distortions of both slow and fast components of the
discontinuous filter inputs.
For high performance applications, the bang-bang
control is replaced by saturation function. Moreover, Fig.2: Conventional sliding mode observer
back emf observer is used to design for better filter-
ing and estimating the rotor speed simultaneously as
shown in Fig.3. This observer is formed as follows [7]:
∗
êα = ω̂e êβ − k2 (êα + zα )
∗
êβ = ω̂e êα − k2 (êβ + zβ ) (5)
∗
ω̂e = (êβ − zβ )êα − (êα − zα )êβ
Where k2 is positive gain. Its proper value is de-
rived by using Lyapunov stability.
Fig.3: Sliding mode observer with saturation func-
The proposed sliding mode observer can be de-
tion
picted in Fig.4. It only consists of the three parts,
i.e., current observer, continuous control function or where a is a positive constant that can be adjusted
sigmoid function and rotor position and speed calcu- the slope of the sigmoid function. The i¯α = îα − iα
lator, which explains as follows: and i¯β = îβ − iβ are the error between estimated
Current observer model: phase current and actual one in the α− and β−axis,
∗ respectively.
R 1 1
îα = − îα + uα − k1 H(îα − iα ) The sliding mode surface is defined as
L L L
T
∗
R 1 1 sn = [sα sβ ] (8)
îβ = − îβ + uβ − k1 H(îβ − iβ ) (6)
L L L where sα = i¯α and sβ = i¯β . When the estimation
where superscript “ ∧ ” represents the estimated error trajectories reach the sliding surface, i.e.,sn = 0.
quantities, k1 is constant current observer gain, and It is obvious that the observed currents will converge
Speed Sensorless Control of PMSM Using An Improved Sliding Mode Observer With Sigmoid Function 53
4. EXPERIMENTAL RESULTS
The hardware configuration of a DSP-based exper-
imental system is shown in Fig.5. The process of the
system is a fixed-point TMS320F2812 DSP with a
clock frequency of 150 MHz. The PMSM is supplied
by a three phase voltage source PWM inverter with
a switching frequency of 20 kHz. A Quadratic En-
coder Pulse (QEP) sensor is used to detect the rotor
position and merely compare with the estimated ro-
Fig.4: Improved sliding mode observer tor position. Two motor phase currents are sensed,
rescaled, and converted to digital values by on-chip
to the actual ones, i.e., îα = iα and îβ = iβ . It is ADC with 12-bit resolution. The PWM gate firing
important that on the sliding surface, the observer signals from the desired phase voltage commands are
will not be affected by any system parameters or any generated by means of the space vector modulation
disturbances. technique. The sampling period of the control system
∗
The sliding mode exists when sn sn < 0, that is is set at 50 µs.
sn → 0 as t → ∞ . Let’s design the Lyapunov func- Fig.6 shows the speed response of PI controller sys-
tion to find such condition of sliding mode existence. tem when using the improved sliding mode observer
at no load condition for a step change of speed from
1 T 1 0.4 PU (360 rpm) to 0.7 PU (630 rpm) at time of 0.18
V = sn sn = (s2α + s2β ) (9)
2 2 s. It is seen that the motor speed can be successfully
The error equation is obtained by subtracting (6) converged to the speed command.
from (1) as Next, the steady-state estimated rotor position re-
∗
sponses and actual rotor position (from QEP) at a
∗ ∗ ∗ R 1 1
s¯α = i¯α = iˆα − iα = − i¯α + eα − k1 H(i¯α ) controlled speed of 0.5 PU (450 rpm) for three differ-
L L L ent structures of sliding mode observers are compared
∗ ∗ ∗ ∗ R 1 1 and shown in Fig.7 through 9. In conventional one
s¯β = i¯β = iˆβ − iβ = − i¯β + eβ − k1 H(i¯β ) (10) and one with saturation function, the rotor position
L L L
∗
errors are significant in the same magnitude, as seen
∗
The sliding mode exists V = sTn sn < 0 [3], i.e., in Figs 7(b) and 8(b). This is due to the effect of de-
lay time of low pass filter. Thus, the estimated rotor
∗ R ³ ¯ 2 ¯ 2´ 1 position should be further compensated by adding an
sTn sn = − iα + iβ + (eα i¯α − k1 i¯α H(i¯α ))
L L offset according to operating speeds. When using the
improved sliding mode observer, the slow component
1 could be extracted directly from the sigmoid function
+ (eβ i¯β − k1 i¯β H(i¯β )) < 0 (11)
L without low pass filter. Thus, it could represent the
As a result, back emf. Then the resulting rotor position can be
k1 ≥ max (|eα |, |eβ |) (12) computed without angle compensation due to delay
time of filter. The rotor position error is greatly re-
Once the sliding mode occurs when selecting a duced as shown in Fig 9(b).
large enough k1 , then the sliding surface becomes as Fig.10 shows the steady-state slow components in
h∗ ∗i the stationary reference frame (α-axis) obtained by
T
sα sβ = [sα sβ ] = [0 0] (13) these sliding mode observers at the same controlled
speed. In conventional one and one with saturation
and (13) can be rewritten as
function, the slow components are obtained from the
low pass filters. It is obvious that the magnitude
(k1 H (i¯α ))eq = eα
of slow component is significantly reduced when the
(k1 H (i¯β ))eq = eβ (14) continuous control (i.e., saturation and sigmoid func-
The chattering problem could be rejected by the tions) is being used rather than discontinuous con-
sigmoid function which is the continuous control. trol (i.e. bang-bang control). This implies that the
Back emf from equation (14) can be used to calcu- high oscillation, causing chattering problem, on the
late the rotor speed and position directly as follows. observed back emf is reduced.
µ ¶
−1 eα
θ̂e = − tan (15)
eβ
dθ̂e
ω̂e = (16)
dt
54 ECTI TRANSACTIONS ON ELECTRICAL ENG., ELECTRONICS, AND COMMUNICATIONS VOL.5, NO.1 February 2007