Continuous Finite-Time TSM Control For Electronic
Continuous Finite-Time TSM Control For Electronic
E-mail: [email protected]
Abstract: This article examines the position tracking difficulties of the electronic throttle (ET) system and presents a continuous
finite-time terminal sliding mode (TSM) control method. The development of this control method involves two procedures: (i)
designing a global finite-time observer to estimate the derivative and higher order derivatives of the ET system output, and the
total disturbances; and (ii) construction of a continuous finite-time TSM controller based on observer estimations, to ensure the
realisation of fast finite-time convergence of system output with a comparatively smooth control action. Comprehensive stability
proof, simulation study, and experimental implementation are provided to affirm the applicability of the proposed method.
J. Eng. 1
This is an open access article published by the IET under the Creative Commons Attribution-NonCommercial-NoDerivs License
(http://creativecommons.org/licenses/by-nc-nd/3.0/)
where Ra is motor resistance, La is motor inductance, ia is armature
current, ke is back-emf coefficient, θm is motor angular
displacement, and θ̇m is motor speed. A low-power DC motor with
small inductance La is often used to drive the throttle plate, and the
current response time constant is much less than the sampling time,
so the inductance dynamics can be ignored [21]. By setting La = 0,
motor armature current can be approximated, resulting into a
simplified form of (1) given as:
Fig. 1 Schematic structure of the electronic throttle system
u − keθ̇m
ia = (2)
Ra Table 1 ET system parameter values
Parameters Definition Values Units
The motor torque equation is expressed as: ks0 spring coefficient 0.1 N m rad−1
N gear transmission ratio 40 —
Jmθ̈m + Bmθ̇m + T m = ktia (3)
Ra0 armature resistance 2.01 Ω
where Bm is motor damping coefficient, Jm is motor inertia, T m is km0 back emf coefficient 0.0217 N m A−1
gear input torque, and kt is torque constant. By combining (2) and Jm0 motor inertia 3 × 10 −6
kg m2
(3), the standard DC motor dynamic equation can be obtained as: Jt 0 throttle torque inertia 2 × 106 kg m2
u − keθ̇m
Jmθ̈m + Bmθ̇m + T m = kt (4)
Ra
3 Control design
The throttle dynamic equation is expressed as follows: A global finite-time observer is implemented to estimate the
derivatives of the ET system output position and the total
Jtθ̈t + Btθ̇t + T sp + T f + T l = T out (5) disturbances, which are then used in the design of the CFTSMC
method. Let θtr be the desired output, then the tracking error is
where T sp, T f , and T out are the non-linear spring torque, coulomb defined as:
friction torque, and gearbox output torque, respectively. T l is load
x1 = θt − θtr
torque, Jt is throttle torque inertia, Bt is throttle damping (10)
coefficient, θt is throttle angular displacement, and θ̇t is throttle x2 = θ̇t − θ̇tr
speed. The non-linear spring torque T sp is expressed as:
Then, system (8) can be rewritten as:
T sp = ks(θt − θt0) (6)
ẋ1 = x2
(11)
and gearbox output torque ẋ2 = Gu − ψ θ(x1 + θtr) − ψ ω(x2 + θ̇tr) + D(t) − θ̈tr
T out = NT m + P(T m) (7) Equation (11) can be described by a general class of second-order
system given below:
where P(T m) is a bounded non-linear function of T m, and θt0 is the
default opening angle of the throttle plate called the limp-home ẋ1 = x2
(12)
position. By combining (4), (5), and (7), and taking into ẋ2 = Gu + ψ(x1, x2) + D(t)
consideration nominal parameter variations, the ET system-state
model can finally be obtained as: where ψ(x1, x2) = − ψ θ(x1 + θtr) − ψ ω(x2 + θ̇tr) − θ̈tr.
2 J. Eng.
This is an open access article published by the IET under the Creative Commons Attribution-NonCommercial-NoDerivs License
(http://creativecommons.org/licenses/by-nc-nd/3.0/)
ż1 = z2 + E1⌈x1 − z1⌋ 1
γ σ = − un (19)
γ2
ż2 = Gu + ψ(x1, z2) + z3 + E2⌈x1 − z1⌋ Taking the derivative of the sliding surface
ż3 = z4 + E3⌈x1 − z1⌋γ3
(13) λ3sgn(σ)
… σ̇ = − (20)
λ1 + αλ2 σ α − 1
γ
żn − 1 = zn + En − 1⌈x1 − z1⌋ n − 1
γ Considering a bounded candidate Lyapunov function V 0(σ, x1, x2)
żn = En⌈x1 − z1⌋ n defined as:
γ
where ⌈x1 − z1⌋γi denotes ⌈x1 − z1⌋γi = x1 − z1 i sign (x1 − z1), 1 1 1
V 0 = σ 2 + x12 + x22 (21)
γi = 1 + iτ, τ ∈ ( − (1/n), 0), (i = 1, 2, …, n) are the observer 2 2 2
coefficients; z1, z2, z3, z4,......, zn − 1 and zn are estimation of x1, x2,
Taking the time derivative of V 0, one can obtain:
D(t), D(i)(t),…, D(n − 3)(t) and D(n − 2)(t); and E1, E2, E3, …, En − 1, En
are the observer gains. Let the estimation error be defined as V̇ 0 = σσ̇ + x1 ẋ1 + x2 ẋ2 (22)
e1 = x1 − z1, e2 = x2 − z2, and e3 = D(t) − z3. Then, the error
dynamics can be written as: ẋ2 can be obtained by combining (15) and (12) as follows:
γ1
ė1 = e2 − E1⌈e1⌋ ẋ2 = Gu + ψ(x1, x^ 2) − [ψ(x1, x^ 2) − ψ(x1, x2)]
γ2
ė2 = e3 − ψ ωe2 − E2⌈e1⌋ ^ ^
+D − [D − D(t)]
(23)
ė3 = e4 − E3⌈e1⌋γ3 = σ + e3 − ψ ωe2 − k2 x^ 2
β2
sgn(x^ 2) − k1 x1
β1
sgn(x1)
(14)
… β2 β1
≤ σ + e3 + ψ ω e2 + k2 x^ 2 + k1 x1
γn − 1
ėn − 1 = en − En − 1⌈e1⌋
γ Since 0 < β < 1, then
ėn = D(n − 2) − En⌈e1⌋ n
β2
According to [23], there exists a finite time T e, such that for x^ 2 < 1 + β2 x^ 2 < 1 + x^ 2
(24)
∀te > T e, ei (i = 1, 2, …, n) will converge to zero in finite time. β1
x1 < 1 + β1 x1 < 1 + x1
Thus, z2 = x2 and z3 = D(t) will be realised.
Based on the definition of the observer estimation error
3.2 Controller design
x^ 2 = x2 + e2 ≤ x2 + e2 (25)
3.2.1 Theorem 3.1: From system (12), based on the observer
estimations, the sliding surface is designed as: α
Also, note that σ < 1 + σ . Substituting (23), (24) and (25) into
^ β2 β1 (22), one obtains
σ = Gu + ψ(x1, x2) + D + k2⌈x2⌋ + k1⌈x1⌋
^ ^
(15)
β
V̇ 0 ≤ x1 x2 + x2σ + x2e3 + ψ ω x2e2
where ⌈x⌋ β = x sgn (x), ki is selected such that the polynomial
+k2(1 + e2 ) x2 + k1(1 + x1 ) x2
n x12 + x22 x22 + e32 x22 + σ 2 x2 + e22
≤ + + + ψω 2
P(λ) = λ + n
∑ kiλ i−1
(16) 2 2 2 2 (26)
i=1 2 2 2 2 2 2
1+e x +e
2 1+x
2 x +x
2 2 1 2
+k2 + k2 + k1 + k1
is Hurwitz. The value of βi is derived based on the following 2 2 2 2
conditions: ≤ Lv1V 0 + Lv2
J. Eng. 3
This is an open access article published by the IET under the Creative Commons Attribution-NonCommercial-NoDerivs License
(http://creativecommons.org/licenses/by-nc-nd/3.0/)
where Q = λ3 /ℓ, ℓ = (λ1 σ 1 − α + αλ2). For σ ≠ 0, α ∈ (0, 1), and
λ1, λ2, λ3 > 0, then ℓ ≠ 0. According to the following inequality
V̇ 1 ≤ − qV 1r (28)
ẋ1 = x2 − e2
β β
ẋ2 = σ − k2⌈x2 − e2⌋ 2 − k1⌈x1⌋ 1 + e3 − ψ ωe2
Γ: σ = − un (29)
λ3sgn(σ)
u̇n =
λ1 + αλ2 σ α − 1
It has been discussed that at time t > T e, the observer errors will
converge to zero in finite time. Also, it has been proven in (26)
based on the proposition in [25], that σ, x1, and x2 will not escape to Fig. 2 Block diagram of the CFTSMC closed-loop control system
infinity in finite time. After t > T e, the closed-loop system
becomes:
ẋ1 = x2
β β1
ẋ2 = σ − k2⌈x2⌋ 2 − k1⌈x1⌋
Γ: σ = − un (30)
λ3sgn(σ)
u̇n =
λ1 + αλ2 σ α − 1
4 J. Eng.
This is an open access article published by the IET under the Creative Commons Attribution-NonCommercial-NoDerivs License
(http://creativecommons.org/licenses/by-nc-nd/3.0/)
Fig. 4 Time history of control input (simulation study)
(a) SMC control input, (b) FOTSMC control input, (c) Proposed control input
Fig. 5 Transient performance comparisons under constant and sinusoidal input voltage disturbance (simulation study)
(a) Constant disturbance, (b) 4 Hz frequency sinusoidal disturbance, (c) 8 Hz frequency sinusoidal disturbance
Fig. 7 Transient performance comparisons for step reference command (experimental case 1)
(a) Transient response, (b) Tracking error, (c) Control input
control. It is as well obvious in Fig. 4 that chattering, which exist in saturated within the input limit of the DC motor. The sampling
the SMC controller, is reduced to a large extent in both the period is selected as 1 ms.
FOTSMC controller and the proposed controller. The degree of
reduction is higher in the proposed controller than in the FOTSMC 4.2.2 Experimental results: Experimental case 1: Nominal
controller as seen from the zoom-in plot. This is due to the effect of performance evaluation
the designed continuous control law and low control gain. Fig. 5 The outcome of the first experimental case is plotted in Figs. 7
shows the performance of the three controllers, subjected to and 8. The proposed controller and FOTSMC controller have faster
constant and sinusoidal input voltage disturbances. The SMC is transient response. For the 4 Hz sinusoidal input command, the
invariant to matched uncertainties, provided the control gain is steady-state error of each controller is obtained as SMC = 1.57 deg,
higher than the upper bound of the disturbance. However, the high FOTSMC = 1.03 deg, and proposed = 0.34 deg. This result reveals
control gain needed to remove the offset increases chattering. The the superiority of the proposed control method with respect to
FOTSMC control law also requires high control gain to remove the tracking accuracy. It is also obvious from the plots that control
offset caused by the disturbances, which could result into steady- chattering is highly reduced in the FOTSMC and proposed control,
state fluctuation. On the contrary, the proposed CFTSMC control with the proposed control having the least chattering effect (Fig. 9).
law has the ability to remove the offset caused by the disturbances,
with small control gain. Experimental case 2: Disturbance rejection ability
In the second experimental test case, the disturbance rejection
4.2 Experimental verification ability of each controllers is verified in the presence of constant
input voltage and sinusoidal input voltage disturbances:
4.2.1 Experimental setup: The experimental substantiation of the
proposed control method is demonstrated in the setup configuration 12V(Nominal), for t ∈ [0, 2]s,
shown in Fig. 6. The set-up consists of a DS1104 dSPACE data
acquisition system and digital signal processor, power supply 7V(Decrease), for t ∈ [2, 8]s,
E= (37)
system, a H-bridge PWM driver, a throttle valve, and a set of high 12V(Increase), for t ∈ [8, 15]s,
processor speed computer on which the dSPACE ControlDesk 4sin(4t), for t ∈ [15, 20]s,
software and Matlab softwares are installed. The input range of the
throttle valve DC motor is ( ± 12V), so the control action is
J. Eng. 5
This is an open access article published by the IET under the Creative Commons Attribution-NonCommercial-NoDerivs License
(http://creativecommons.org/licenses/by-nc-nd/3.0/)
Fig. 8 Transient performance comparisons for 4 Hz frequency sinusoidal reference command (experimental case 1)
(a) Transient response, (b) Tracking error, (c) Control input
Fig. 9 Transient performance comparisons under constant and sinusoidal input voltage disturbance (experimental case 2)
(a) Transient response, (b) Tracking error, (c) Control input
Table 2 Transient performance comparisons under constant and sinusoidal input voltage disturbances (experimental case 2)
Disturbances Algorithms Max(e), deg Recovery time, s RMS(e), deg
constant SMC 2.93 — 1.33
FOTSMC 1.64 >3.0 0.59
proposed 1.59 0.18 0.16
sinusoidal SMC 2.30 — 0.87
FOTSMC 1.43 — 0.79
proposed 0.49 0.18 0.11
6 J. Eng.
This is an open access article published by the IET under the Creative Commons Attribution-NonCommercial-NoDerivs License
(http://creativecommons.org/licenses/by-nc-nd/3.0/)
[16] Wang, H., Shi, L., Man, Z., et al.: ‘Continuous fast nonsingular terminal [22] Du, H., Qian, C., Yang, S., et al.: ‘Recursive design of finite-time convergent
sliding mode control of automotive electronic throttle systems using finite- observers for a class of time-varying nonlinear systems’, Autom., 2013, 49,
time exact observer’, IEEE Trans. Ind. Electron., 2018, 65, (9), pp. 7160– (2), pp. 601–609
7172 [23] Wang, Z., Li, S., Yang, J., et al.: ‘Current sensorless finite-time control for
[17] Fridman, L., Moreno, J.A., Bandyopadhyay, B., et al.: ‘Continuous nested buck converters with time-varying disturbances’, Control Eng. Pract., 2018,
algorithms: the fifth generation of sliding mode controllers’ (Springer 77, pp. 127–137
International Publishing, Berlin Heidelberg, Germany, 2015) [24] Bhat, S.P., Bernstein, D.S.: ‘Geometric homogeneity with applications to
[18] Emel'Yanov, S.V., Korovin, S.K., Levant, A.: ‘High-order sliding modes in finite-time stability’, Math. Control Signals Syst., 2005, 17, (2), pp. 101–127
control systems’, Comput. Math. Model., 1996, 7, (3), pp. 294–318 [25] Li, S., Tian, Y.P.: ‘Finite-time stability of cascaded time-varying systems’, Int.
[19] Levant, A.: ‘Sliding order and sliding accuracy in sliding mode control’, Int. J. Control, 2007, 80, (4), pp. 646–657
J. Control, 1993, 58, (6), pp. 1247–1263 [26] Bhat, S.P., Bernstein, D.S.: ‘Finite-time stability of continuous autonomous
[20] Feng, Y., Han, F., Yu, X.: ‘Chattering free full-order sliding-mode control’, systems’, SIAM J. Control Optim., 2000, 38, (3), pp. 751–766
Autom., 2014, 50, (4), pp. 1310–1314
[21] Jiao, X., Zhang, J., Shen, T.: ‘An adaptive servo control strategy for
automotive electronic throttle and experimental validation’, IEEE Trans. Ind.
Electron., 2014, 61, (11), pp. 6275–6284
J. Eng. 7
This is an open access article published by the IET under the Creative Commons Attribution-NonCommercial-NoDerivs License
(http://creativecommons.org/licenses/by-nc-nd/3.0/)