Next Article in Journal
Hemodynamic Analysis on the Anastomosis Angle in Arteriovenous Graft Using Multiphase Blood Model
Previous Article in Journal
Blockchain-Based Information Management for Supply Chain Data-Platforms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Disturbance Observer-Based Chattering-Attenuated Terminal Sliding Mode Control for Nonlinear Systems Subject to Matched and Mismatched Disturbances

1
Graduate School of Mechanical Engineering, University of Ulsan, Ulsan 44610, Korea
2
Automatic Control Department, Ho Chi Minh City University of Technology and Education, Ho Chi Minh City 700000, Vietnam
3
School of Mechanical Engineering, University of Ulsan, Ulsan 44610, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(17), 8158; https://doi.org/10.3390/app11178158
Submission received: 7 August 2021 / Revised: 29 August 2021 / Accepted: 30 August 2021 / Published: 2 September 2021
(This article belongs to the Section Robotics and Automation)

Abstract

:
This work presents a new disturbance observer-based chattering-attenuated terminal sliding mode control for a class of nonlinear systems in the presence of both mismatched and matched disturbances. A nonlinear disturbance observer is typically employed to accurately estimate mismatched disturbances. In this study, a terminal sliding mode control was designed, based on the disturbance estimation results, to counter the effects of disturbances and ultimately stabilize the target system. The utilization of a chattering-attenuated full-order terminal sliding mode structure satisfactorily resolves both chattering and singularity problems in controller design. It was shown by theoretical analyses that both the disturbance estimation error and the system state converge to the equilibrium point in finite time. Two simulation studies, namely a numerical example and an application to an electro hydrostatic actuator system, were conducted to examine the characteristics and to verify the effectiveness of the proposed algorithm.

1. Introduction

Sliding mode control (SMC) is a popular control algorithm and possesses several desirable characteristics, such as design simplicity, ease of implementation, and robustness to external disturbances, model uncertainties, and parameter variations. Thus, SMC has received a significant amount of attention and is widely found in diverse practical applications, including robotics [1,2,3,4,5,6], electrical and mechanical systems [7,8,9,10,11,12], and aerospace engineering [13,14,15,16,17]. Nonetheless, the conventional SMC still exhibits several drawbacks, namely chattering, asymptotic convergence, and susceptibility to mismatched disturbances.
Chattering phenomena occur in the conventional SMC due to the use of a switching function, which leads to high oscillations in the control signal [18]. Significant efforts have been made in the literature to attenuate this problem, including the boundary layer [19], second and higher-order SMC [20,21,22], low-pass filtering [23,24], and the disturbance observer [25]. Another problem found in the conventional SMC is that system states can only asymptotically converge to the equilibrium point. By comparison, terminal sliding mode control (TSMC) can drive system states both to the sliding mode surface and to the equilibrium point in finite time [26]. Although the conventional TSMC inherits the robustness characteristic from its predecessor, it still suffers from the chattering phenomenon and is prone to the singularity problem as the control signal instantly jumps to infinity. A large number of studies have been devoted to solving this problem, e.g., [27,28,29]. Recently, a chattering-free full-order TSMC was developed in [30] that satisfactorily resolves both the chattering and singularity problems mentioned above.
A rather less common problem that SMC faces is mismatched disturbances, which act in a different channel than the control input. Some examples in which mismatched disturbances exist are magnetic levitation (MAGLEV) suspension systems, permanent magnet synchronous motor (PMSM), systems and flight control systems [31]. The various methods that contribute to tackling this problem include integral SMC [32,33,34] and adaptive SMC [35,36,37]. However, it is well known that the integral action may produce several undesirable effects, such as a large overshoot and a long settling time for the control system. Moreover, the nominal performance is normally sacrificed for robustness against disturbances in these methods. Disturbance observer-based SMC, in contrast, is a more promising approach to deal with mismatched disturbances [31,38,39,40,41,42]. In [38], a nonlinear disturbance observer is introduced to develop a new sliding mode surface in an attempt to attenuate mismatched disturbances. However, the mismatched disturbance considered in [38] must be constant rather than time varying, which is a highly restrictive constraint, because the disturbance can continuously change over time in practical applications; that is, unless the mismatched disturbance is time-invariant, the algorithm in [38] fails to stabilize the control system. In addition, the nonlinear disturbance observer utilized in [38] has been proven to be exponentially stable only in the case of slow time-varying mismatched disturbances [43]. A remedy for this constant restriction on the mismatched disturbance formulation is provided in [39], wherein the proposed extended disturbance observer can estimate nonlinear and time-varying mismatched disturbances. However, the problem of being Lyapunov stable without convergence remains, thus limiting the potential of the control algorithm. Another effort to stabilize nonlinear systems is found in [41], where mismatched disturbances are estimated and then suppressed by an extended-state observer-based sliding mode control. Both matched and mismatched disturbances are considered in the paper. Nonetheless, the estimation errors and the system state can only reach and vary in a small region around the equilibrium point, indicating that the observer and the closed-loop system are only Lyapunov stable without convergence. Consequently, the mismatched disturbance cannot be estimated precisely, and the system state cannot be stabilized in finite time.
Motivated by previous works, to accurately estimate and eventually counteract the effect of mismatched disturbances for a class of nonlinear systems, a new design of disturbance observer-based TSMC is introduced in this paper. The finite-time disturbance observer employed in this work originates from an arbitrary-order exact robust differentiator developed in [21]. It forms the basis for the development of this disturbance observer, which has been used by a large number of researchers in various applications. Some of the works are discussed as follows. In [44], the above-mentioned observer is implemented to tackle the tracking control problem of robotic manipulators. In [45], the influence of mismatched disturbances on a class of nonlinear systems is attenuated based on this observer. Simulation of a PMSM then proved the effectiveness of the proposed controller. The same disturbance observer is employed in [46] to address the problem of airgap control of a MAGLEV suspension vehicle. Another application of this disturbance observer can be found in [31], where the authors utilize it to control a DC-DC buck converter. Simulation and experiments on the converter have proven its disturbance rejection capability. These references convincingly demonstrated the effectiveness and the popularity of the disturbance observer considered in the current paper, for estimating and suppressing the influences of matched and mismatched disturbances in different highly nonlinear dynamic systems.
Major aspects of the proposed control algorithm are presented as follows. First, mismatched disturbances are accurately estimated by the nonlinear finite-time disturbance observer. Second, the estimation results are integrated into a new terminal sliding mode surface to enhance system robustness against external disturbances, uncertainties, and unmodeled dynamics. Third, the utilization of a chattering-attenuated full-order terminal sliding mode surface structure in controller design effectively suppresses chattering and singularity phenomena in the control signal. Fourth, Lyapunov stability analysis alongside the design step theoretically increases the reliability of the proposed control scheme. In addition, although the control algorithm is developed for a class of nonlinear second-order systems, the obtained results are extendable to the nth-order formulation. Finally, comparative simulation studies, including a numerical example and an application to an Electro Hydrostatic Actuator (EHA) system, demonstrate the effectiveness of the proposed algorithm.
This work can be perceived to be an extension and improvement of the references [30,41] for several reasons. First, in [30], rigorous research on chattering-free full-order terminal sliding mode control was presented, thus forming the backbone of the development of our work. However, because the control algorithm developed in [30] was only designated to tackle the problem of matched disturbances, it cannot cope with mismatched disturbances, which are more challenging. Second, in [41], an extended-state observer-based sliding mode control (ESOSMC) was introduced in an attempt to address the problem of mismatched disturbances. However, this approach nonetheless exhibited some drawbacks. Although it inherited the chattering attenuation characteristic from the work in [30], the algorithm proposed in [41] was only able to drive the disturbance estimation and system state to a small, bounded region around the equilibrium without actually converging to it. Moreover, the mismatched disturbance used in [41] gradually converged to a constant rather than being time varying throughout the system operation. Thus, the proposal was not practical and was unable to demonstrate the effectiveness of the controller. To address these problems, in this study, a new terminal sliding mode surface was developed. It is enhanced with disturbance rejection by a nonlinear disturbance observer. The proposed control algorithm is expected to suppress both time-varying matched and mismatched disturbances and provide a finite-time convergence characteristic. A brief comparison between the proposed controller and the works presented in [30] and [41] is shown in Table 1 below.
The remainder of the paper is organized as follows. Section 2 presents the problem formulation. Section 3 introduces the design procedure of the proposed algorithm, followed by stability analysis. Section 4 illustrates the effectiveness of the proposed approach via comparative simulation studies. Ultimately, the conclusion and future work are given in Section 5.

2. Problem Statement

Consider a class of second-order nonlinear systems that are susceptible to both matched and mismatched disturbances as follows [41]:
{ x ˙ 1 ( t ) = x 2 ( t ) + d 1 ( t ) x ˙ 2 ( t ) = f ( x , t ) + g ( x , t ) u ( t ) + d 2 ( t ) y ( t ) = x 1 ( t )
where x = [ x 1 , x 2 ] T is the system state and assumed to be available, f ( x , t ) and g ( x , t ) 0 are two known smooth nonlinear functions, u ( t ) is the control input, y ( t ) is the system output, and d i ( t ) , i = 1 , 2 are the time-varying mismatched and matched disturbances of the system.
Assumption 1.
Both mismatched and matched disturbances  d i ( t ) , i = 1 , 2 , and their derivatives, are bounded and satisfy:
| d i ( t ) | ϕ i , | d ˙ i ( t ) | ρ i , i = 1 , 2  
where  ϕ i > 0 and ρ i > 0 , ( i = 1 , 2 ) are positive constants. This assumption is reasonably realistic due to the physical limitation of machinery in practical applications.
This work aims to design a disturbance observer-based terminal sliding mode controller to drive the system state x 1 ( t ) to the origin in finite time in the presence of matched and mismatched disturbances. For this purpose, where the context is sufficiently explicit, the arguments of some functions might be omitted for brevity.
For the system (1), the chattering-free full-order TSMC developed in [30] is derived as below. First, a new terminal sliding mode surface is introduced with the form:
s = x ˙ 2 + c 2 | x 2 | α 2 sign ( x 2 ) + c 1 | x 1 | α 1 sign ( x 1 )
Then, a control law is designed accordingly as:
u = g 1 ( x , t ) ( u e q + u n )
where u e q denotes the equivalent control and u n denotes the switching term. The equivalent control u e q is designed as:
u e q = f ( x , t ) c 2 | x 2 | α 2 sign ( x 2 ) c 1 | x 1 | α 1 sign ( x 1 )
the switching term u n is defined and filtered by the following method:
u ˙ n + T u n = v
v = ( ρ 2 + k T + η ) sign ( s )
where η is a positive constant, T 0 , k T is chosen satisfying k T T ϕ 2 , c i and α i   ( i = 1 , 2 ) are positive constants. c i can be selected to ensure that the polynomial p 2 + c 2 p + c 1 is Hurwitz, and α i can be determined to satisfy [30,47]:
{ α 1 = α 2 2 α 2 α 2 = α , α ( 1 ε , 1 ) , ε ( 0 , 1 )
By substituting the system (1) and the control law (4) and (5) into the sliding surface (3), one can obtain:
s = x ˙ 2 + c 2 | x 2 | α 2 sign ( x 2 ) + c 1 | x 1 | α 1 sign ( x 1 ) = f ( x , t ) + g ( x , t ) u + d 2 ( t ) + c 2 | x 2 | α 2 sign ( x 2 ) + c 1 | x 1 | α 1 sign ( x 1 ) = d 2 ( t ) + u n
Consider a candidate Lyapunov function of the form: V ( s ) = s 2 / 2 . It was proven in [48] that the switching term satisfies | u n ( t ) | ϕ 2 . This statement is found to be crucial in proving the finite-time convergence of the sliding surface s . Thus, it is re-depicted here for clarity. Initially, u n ( t = 0 ) = 0 is taken, i.e., u n ( 0 ) [ ϕ 2 , ϕ 2 ] . Then, u n indicates that it cannot leave the interval [ ϕ 2 , ϕ 2 ] . Recall | d 2 | ϕ 2 , | d ˙ 2 | ρ 2 in (2), k T T ϕ 2 in the first line after (7), and denote β = ρ 2 + k T + η . Due to s = u n + d 2 in (9), (6) can be rewritten as:
u ˙ n + T u n = β sign ( s ) = β sign ( u n + d 2 )
If u n > ϕ 2 , then u n + T ϕ 2 > 0 , sign ( u n + d 2 ) = sign ( u n ) = 1 . It can be obtained from (10) that:
u ˙ n = T u n β < 0
Similarly, if u n < ϕ 2 , then sign ( u n + d 2 ) = sign ( u n ) = 1 . This leads to:
u ˙ n = T u n + β > 0
Thus, the inequality u n u ˙ n < 0 holds outside of [ ϕ 2 , ϕ 2 ] , and consequently u n cannot leave [ ϕ 2 , ϕ 2 ] , i.e., | u n | ϕ 2 and, accordingly, | d ˙ 2 ( t ) + T u n ( t ) | ρ 2 + T ϕ 2 are retained for any t 0 . Differentiating the Lyapunov function V ( s ) yields:
s s ˙ = s ( d ˙ 2 + u ˙ n ) = s ( d ˙ 2 + ( u ˙ n + T u n ) T u n ) = s ( d ˙ 2 + v T u n ) ( | d ˙ 2 T u n | ( ρ 2 + k T + η ) ) | s | = ( | d ˙ 2 T u n | ( ρ 2 + T ϕ 2 ) ( k T T ϕ 2 ) η ) | s | η | s |
The inequality (11) indicates that the system states will reach the sliding surface s = 0 in finite time. Once the desired sliding mode s = 0 is established, the dynamics of the system can be expressed as:
{ x ˙ 1 = x 2 + d 1 ( t ) x ˙ 2 = c 2 | x 2 | α 2 sign ( x 2 ) c 1 | x 1 | α 1 sign ( x 1 )
In the case in which there are no mismatched disturbances, (12) becomes:
{ x ˙ 1 = x 2 x ˙ 2 = c 2 | x 2 | α 2 sign ( x 2 ) c 1 | x 1 | α 1 sign ( x 1 )
From (13), if c i and α i , i = 1 , 2 are chosen as below (7) and as in (8), respectively, then the system states will converge to the equilibrium in finite time [47,49]. However, in the case in which mismatched disturbances exist, based on the equation set (12), one can obtain the below relationship:
x ¨ 1 + c 2 | x ˙ 1 d 1 | α 2 sign ( x ˙ 1 d 1 ) + c 1 | x 1 | α 1 sign ( x 1 ) = d ˙ 1
It is clear from the above equation that, despite the control effort, the mismatched disturbance still exists in the system dynamics. Thus, the resulting system state is strongly affected by the mismatched disturbance; that is, although the system state can be attracted to the sliding surface in finite time, it cannot converge to the equilibrium point due to the presence of the time-varying mismatched disturbance. As a result, the chattering-free TSMC developed in [30] can cope with matched disturbances, yet is unable to overcome the influence of mismatched disturbances.

3. Results

In this section, a new chattering-attenuated full-order terminal sliding mode control powered by a nonlinear disturbance observer is proposed for the system (1). First, the disturbance observer is employed to estimate the mismatched disturbance. Subsequently, a new terminal sliding surface and a control law are developed utilizing the estimation results to ensure that the system state converges to the origin in finite time.

3.1. Nonlinear Disturbance Observer

This section aims to employ a finite-time disturbance observer to estimate the time-varying mismatched disturbance existing in (1). Based on the research on globally convergent differentiators in [50], a nonlinear disturbance observer is designed as follows:
{ z ˙ 1 = v 1 + x 2 v 1 = λ 1 L 1 / 3 | z 1 x 1 | 2 / 3 sign ( z 1 x 1 ) μ 1 ( z 1 x 1 ) + z 2 z ˙ 2 = v 2 v 2 = λ 2 L 1 / 2 | z 2 v 1 | 1 / 2 sign ( z 2 v 1 ) μ 2 ( z 2 v 1 ) + z 3 z ˙ 3 = λ 3 L sign ( z 3 v 2 ) μ 3 ( z 3 v 2 )
where z 1 = x ^ 1 , z 2 = d ^ 1 , z 3 = d ˙ ^ 1 , and L > 0 , λ i , μ i > 0 ( i = 1 , 2 , 3 ) .
Assumption 2.
The time-varying mismatched disturbance d 1 ( t ) is assumed to be second-order differentiable and has a Lipschitz constant L .
Theorem 1.
Given that the disturbance observer is designed as in (15), if the design parameters are properly selected, the following results can be achieved in finite time:
z 1 = x 1 , z 2 = d 1 , z 3 = d ˙ 1
Proof of Theorem 1.
For the disturbance observer (15), estimation errors are defined as follows:
ε 1 = z 1 x 1 L , ε 2 = z 2 d 1 L , ε 3 = z 3 d ˙ 1 L
By taking the derivative of ε 1 , ε 2 and ε 3 , the observer error dynamics can be attained as:
{ ε ˙ 1 = λ 1 | ε 1 | 2 / 3 sign ( ε 1 ) μ 1 ε 1 + ε 2 ε ˙ 2 = λ 2 | ε 2 ε ˙ 1 | 1 / 2 sign ( ε 2 ε ˙ 1 ) μ 2 ( ε 2 ε ˙ 1 ) + ε 3 ε ˙ 3 = λ 3 sign ( ε 3 ε ˙ 2 ) μ 3 ( ε 3 ε ˙ 2 ) 1 L d ¨ 1
Then, according to Theorem 3.1, Lemma B.1, and their proofs in [51], the estimation errors (18) will converge to zero in finite time, indicating that the nonlinear disturbance observer can accurately estimate the system state x 1 ( t ) , the time-varying mismatched disturbance d 1 ( t ) , and its derivative d ˙ 1 ( t ) . This completes the proof of Theorem 1.  □
Remark 1.
Readers are strongly encouraged to inspect the works [21,50,51] and the references therein for better understanding of the rigorous stability analysis of the disturbance observer, and the appropriate selection of its design parameters. It is noted that the observer design parameters,  { λ i , μ i } ( i = 1 , 2 , 3 ) , are chosen recursively in such a manner that { λ i , μ i } provide for the convergence of the observer with the Lipschitz constant L indicated in Assumption 2. It is also mentioned by the authors that these parameters can be easily changed because the design is not highly sensitive to their values. Nonetheless, the tradeoff is as follows: generally, the greater the parameters, the faster the convergence of the disturbance observer, but the higher the sensitivity to input noises and the sampling step, and the larger the peaking phenomenon in the control signal, if the initial values of the system states and those of the observer states are different. Thus, in the controller design, particularly in the case of different initial values, a tradeoff should be made to ensure fast convergence yet an acceptable peaking phenomenon.

3.2. Disturbance Observer-Based Terminal Sliding Mode Control

In this section, based on the disturbance estimation, a new sliding surface and, accordingly, a new control law are proposed to drive the system state to the equilibrium point in finite time. This work can be considered as an extension to the control strategy developed in [30] and an improvement of the work [41].
A new disturbance observer (DOB)-based terminal sliding mode surface for the system (1) is proposed as:
σ = x ˙ 2 + c 2 | x 2 + z 2 | α 2 sign ( x 2 + z 2 ) + c 1 | x 1 | α 1 sign ( x 1 ) + z 3
where z 2 = d ^ 1 , z 3 = d ˙ ^ 1 are obtained from the disturbance observer, and c i and α i   ( i = 1 , 2 ) are positive constants. c i are chosen to guarantee that the polynomial p 2 + c 2 p + c 1 is Hurwitz, and α i are determined in the same manner as in (8).
Theorem 2.
Given the sliding surface (19), if the control law is developed as:
u = g 1 ( x , t ) ( u e q + u n )
with the equivalent control:
u e q = f ( x , t ) c 2 | x 2 + z 2 | α 2 sign ( x 2 + z 2 ) c 1 | x 1 | α 1 sign ( x 1 ) z 3
and the switching term filtered for chattering attenuation is:
u ˙ n + T u n = v
v = ( ρ 2 + k T + η ) sign ( σ )
where T , k T , η are defined as same as below (7), then the system (1) will reach the sliding mode surface σ = 0 in finite time, and the system state x 1 ( t ) will also converge to the origin in finite time.
Proof of Theorem 2.
By substituting the expression of x ˙ 2 ( t ) in (1) and the control law (20), (21) into the sliding surface (19), one can derive the following equality:
σ = x ˙ 2 + c 2 | x 2 + z 2 | α 2 sign ( x 2 + z 2 ) + c 1 | x 1 | α 1 sign ( x 1 ) + z 3 = [ f ( x , t ) + g ( x , t ) u + d 2 ( t ) ] + c 2 | x 2 + z 2 | α 2 sign ( x 2 + z 2 ) + c 1 | x 1 | α 1 sign ( x 1 ) + z 3 = d 2 ( t ) + u n
Consider a candidate Lyapunov function as V ( σ ) = σ 2 / 2 . Differentiating the Lyapunov function V ( σ ) and following the same procedure in (11) yields:
V ˙ ( σ ) = σ σ ˙ = σ [ d ˙ 2 ( t ) + u ˙ n ] η | σ |
which implies that the system will arrive at the sliding surface σ = 0 in finite time. Once the desired sliding mode occurs, combining system (1) and sliding surface (19) yields:
{ x ˙ 1 = x 2 + d 1 x ˙ 2 = c 2 | x 2 + z 2 | α 2 sign ( x 2 + z 2 ) c 1 | x 1 | α 1 sign ( x 1 ) z 3
From the equation set (26), the following relationship can be established:
x ¨ 1 = x ˙ 2 + d ˙ 1 = [ c 2 | x 2 + z 2 | α 2 sign ( x 2 + z 2 ) c 1 | x 1 | α 1 sign ( x 1 ) z 3 ] + d ˙ 1 = c 2 | ( x ˙ 1 d 1 ) + z 2 | α 2 sign ( ( x ˙ 1 d 1 ) + z 2 ) c 1 | x 1 | α 1 sign ( x 1 ) z 3 + d ˙ 1 = c 2 | x ˙ 1 + ( z 2 d 1 ) | α 2 sign ( x ˙ 1 + ( z 2 d 1 ) ) c 1 | x 1 | α 1 sign ( x 1 ) ( z 3 d ˙ 1 )
According to Theorem 1, z 2 = d 1 and z 3 = d ˙ 1 are achieved in finite time by the disturbance observer. Thus, after the convergence of the disturbance observer, the dynamics (27) become:
x ¨ 1 = c 2 | x ˙ 1 | α 2 sign ( x ˙ 1 ) c 1 | x 1 | α 1 sign ( x 1 )
Let w 1 = x 1 , w 2 = x ˙ 1 ; then, the dynamic Equation (28) can be transformed to:
{ w ˙ 1 = w 2 w ˙ 2 = u
where u = c 2 | w 2 | α 2 sign ( w 2 ) c 1 | w 1 | α 1 sign ( w 1 ) , and c i , α i   ( i = 1 , 2 ) are designed in the same manner as below equation (19) and equation (8). Then, according to Proposition 8.1 and its proof in [47], it can be concluded that the origin is a globally finite-time-stable equilibrium for the system (29), indicating that the state x 1 ( t ) of the system (1) will eventually converge to zero in finite time. This completes the proof for Theorem 2.  □
Remark 2.
As stated in [30], the control signal (22) acts in the same manner as a low-pass filter, wherein ω = T can be interpreted as the bandwidth of the filter. The non-smooth switching function u n is therefore smoothed to prevent the chattering phenomena from occurring, given that the design parameters have been properly chosen. Moreover, because it is a full-order terminal sliding mode surface, it allows the corresponding control law to be derived without differentiating the term c i | | α i sign ( ) , thus avoiding the singularity phenomena.
Remark 3.
It should be mentioned that the sliding surface (19) includes the “acceleration” signal x ˙ 2 ( t ) , which can be obtained by differentiating the “velocity” signal x 2 ( t ) . Nonetheless, this action is sensitive to noise and thus can degrade system performance. Note that the proposed control law only utilizes the sign of the sliding surface σ rather than its value. Based on this observation, the same idea suggested in [30] is re-depicted here to avoid the use of the “acceleration” signal. In this approach, sign ( σ ) can be alternatively determined as follows:
sign ( σ ) = sign ( h ( t ) h ( t τ ) )  
where h ( t ) is defined as:
h ( t ) = 0 t σ ( t ) d t = x 2 + 0 t ( c 2 | x 2 + z 2 | α 2 sign ( x 2 + z 2 ) + c 1 | x 1 | α 1 sign ( x 1 ) + z 3 ) d t
and τ is chosen to be the fundamental sampling time.
Alternatively, the nonlinear disturbance observer structure can be applied again to meet this need and is presented below:
{ z ˙ 4 = v 4 + x 3 v 4 = λ 4 L 2 1 / 3 | z 4 x 2 | 2 / 3 sign ( z 4 x 2 ) μ 4 ( z 4 x 2 ) + z 5 z ˙ 5 = v 5 v 5 = λ 5 L 2 1 / 2 | z 5 v 4 | 1 / 2 sign ( z 5 v 4 ) μ 5 ( z 5 v 4 ) + z 6 z ˙ 6 = λ 6 L 2 sign ( z 6 v 5 ) μ 6 ( z 6 v 5 )  
where x 3 = f ( x , t ) + g ( x , t ) u , z 4 = x ^ 2 , z 5 = d ^ 2 , z 6 = d ˙ ^ 2 . Because the disturbance observer is proven to be finite-time convergent, after a transient period, z ˙ 4 = x ^ ˙ 2 can be regarded as the “acceleration” signal. In addition, the matched disturbance d 2 ( t ) can be accurately estimated for further potential, such as reducing the switching gain in the control signal (23), thus decreasing the chance that chattering phenomena will occur. However, the use of another disturbance observer will increase the computational burden, and consequently, a longer response time is expected, which might represent a disadvantage in practice.
Remark 4.
Because the mismatched disturbance is absent from the control system, if the initial values of the observer states are selected as z 1 ( t 0 ) = x 1 ( t 0 ) , z 2 ( t 0 ) = z 3 ( t 0 ) = 0 , then it is derived from the observer formulation and observer error dynamics that ε 1 ( t ) = ε 2 ( t ) = ε 3 ( t ) = 0 and v 1 ( t ) = v 2 ( t ) = 0 during the absence of mismatched disturbances. In this case, the sliding surface (19) and the control law (20) reduce to those of the conventional TSMC. This indicates the nominal control performance preservation capability of the proposed controller.
Remark 5.
Because the formulation developed earlier in this paper is for second-order systems only, it is desirable to generalize it to an nth-order formulation for different applications. Therefore, the results obtained in this work are extended for a class of nonlinear nth-order systems of the form:
{ x ˙ i = x i + 1 + d i ( t ) , i = 1 , 2 , , n 1 x ˙ n = f ( x , t ) + g ( x , t ) u + d n ( t )  
To estimate the mismatched disturbances d i ( i = 1 , 2 , , n 1 ) , an nth-order disturbance observer is applied here as follows:
{ z ˙ 1 i = v 1 i + x i + 1 = φ 1 i ( z 1 i x i ) + z 2 i + x i + 1 z ˙ 2 i = v 2 i = φ 2 i ( z 2 i v 1 i ) + z 3 i z ˙ ( n + 1 ) i = φ ( n + 1 ) i ( z ( n + 1 ) i v n i )  
where φ j i are nonlinear functions of the form:
φ j i ( s ) = λ j i L i 1 / ( n j + 2 ) | s | ( n j + 1 ) / ( n j + 2 ) sign ( s ) μ j i ( s )  
and z 1 i = x ^ i , z 2 i = d ^ i , z 3 i = d ˙ ^ i ,…, z ( n + 1 ) i = d i ( n 1 ) ^ , L i > 0 , λ j i , μ j i > 0 , j = 1 , 2 , , n + 1 . Consequently, an nth-order disturbance observer-based TSMC is designed as:
σ = x ˙ n + k = 2 n ( c k | ψ k | α k sign ( ψ k ) ) + c 1 | x 1 | α 1 sign ( x 1 ) + χ z  
where ψ k = x k + z k 1 + z ( k 1 ) 2 + + z 2 ( k 1 ) , χ z = z 3 ( n 1 ) + z 4 ( n 2 ) + z ( n + 1 ) 1 . c k and α k are positive constants such that the polynomial p n + k = 2 n ( c k p k 1 ) + c 1 is Hurwitz and the following equation set is satisfied [30,47]:
{ α 1 = α , n = 1 α k 1 = α k α k + 1 2 α k + 1 α k , k = 2 , , n n 2  
where α n + 1 = 1 , α n = α , α ( 1 ε , 1 ) , ε ( 0 , 1 ) .

4. Simulations

In this section, two simulation studies are presented to verify the effectiveness of the control algorithm. Although a numerical example is provided to theoretically assess and discuss the characteristics and performance of the proposed controller, the task of controlling an electro hydrostatic actuator system partly reflects the practicality of the suggested control algorithm. These simulations were conducted in the MATLAB-Simulink environment, wherein the fixed-step ODE4 (Runge–Kutta) was chosen to be the solver. The sampling time was determined separately for each case study, and other settings remained as the defaults.

4.1. Numerical Example

To demonstrate the effectiveness of the proposed control algorithm, the succeeding nonlinear second-order system, which is subject to time-varying mismatched and matched disturbances, is considered:
{ x ˙ 1 = x 2 + d 1 ( t ) x ˙ 2 = 2 x 1 x 2 + e x 1 + u + d 2 ( t )
Here, the primary objective is to stabilize the above system, and the initial system state is selected away from the equilibrium as x = [ 1 , 1 ] T . Mismatched and matched disturbances are deliberately defined as d 1 ( t ) = 1 + 0.3 cos ( 2 t ) sin ( t ) and d 2 ( t ) = sin ( 2 t ) , respectively. The mismatched disturbance is selected such that as the running time increases and the system state x 1 ( t ) reaches its steady-state value. The mismatched disturbance is still time varying in a relatively large, bounded region rather than gradually converging to a constant as in the reference [41]. This is a decisive factor in demonstrating the superiority of the proposed control algorithm compared to strategies considered, which will be illustrated by later simulation results. Moreover, to illustrate the nominal control performance preservation characteristic of the proposed control algorithm, the mismatched disturbance is only active after t = 5 s .
For comparison purposes, the chattering-free full-order TSMC approach [30], denoted as TSMC, the extended-state observer-based chattering free SMC [41], denoted as ESOSMC, and the proposed control algorithm, denoted as DOBTSMC, are employed. It should be noted that both DOBTSMC and ESOSMC stem from TSMC, and that the proposed control algorithm, to a certain extent, can be interpreted as an extension and improvement of TSMC and ESOSMC in dealing with mismatched disturbances. Hence, the analogy is drawn here to demonstrate this argument.
Briefly, the ESOSMC proposed in [41] is constructed of a nonlinear extended-state observer (ESO) and an ESO-based SMC. The ESO is designed based on a hyperbolic tangent function of the form:
{ Ζ ˙ 1 = Ζ 2 a 1 ( Ζ 1 x 1 ) + x 2 Ζ ˙ 2 = a 2 tanh ( b ( Ζ 1 x 1 ) )
where Ζ 1 = x ^ 1 , Ζ 2 = d ^ 1 , tanh ( ) = ( e e ) / ( e + e ) , a i > 0 , i = 1 , 2 , b > 0 . The sliding mode surface integrated with the ESO is designed as:
s = x ˙ 2 + c 2 ( x 2 + Ζ 2 ) + c 1 x 1 + Ζ ˙ 2
where c i , i = 1 , 2 are chosen in the same manner as in previous sections.
The design control parameters are specifically selected to be the same for all control algorithms as c 1 = 5 , c 2 = 5 , T = 0.1 , ( ρ 2 + k T + η ) = 3 . In addition, for the TSMC and the DOBTSMC, following (8), α i ( i = 1 , 2 ) are designed as α 1 = 3 / 7 , α 2 = 3 / 5 . For the mismatched disturbance observers, ESO parameters are selected as a 1 = 20 , a 2 = 14 and b = 21 ; and DOB parameters are chosen as L = 3 , λ 1 = 6 , λ 2 = 11 , λ 3 = 6 , μ 1 = 18 , μ 2 = 6 , and μ 3 = 2 . The sampling time is set to be 0.1 ms .
The system state x 1 ( t ) delivered by all control algorithms is shown in Figure 1. For the first 5 s , all controllers are able to stabilize the target system. However, the convergence rate of ESOSMC is lower than that of the other two approaches. It takes ESOSMC about 3.4 s to reach the 0.01 mark (~1% of the initial value), whereas it takes TSMC and DOBTSMC only around 1.8 s to reach this mark. Moreover, because the trajectory of TSMC and DOBTSMC are identical in this period, it proves that the proposed control algorithm possesses the nominal performance preservation characteristic in the absence of mismatched disturbances. After t = 5 s , mismatched disturbance starts to affect the control system. It is noticeable that TSMC fails to stabilize the target system, which is already anticipated in the Problem Statement section. Consequently, TSMC signals are suppressed in the following figures for better visualization of other controllers. It can be seen in Figure 1 that, although ESOSMC can only drive the system state to a bounded region of about 0.08 around the equilibrium, DOBTSMC offers much higher accuracy, with a bound of [ 2.5 × 10 8 , 1.5 × 10 8 ] .
Figure 2 depicts the mismatched disturbance existing in the target system and the second state x 2 ( t ) response of ESOSMC and DOBTSMC. The figure illustrates that the mismatched disturbance is nonlinear and time varying in a bounded region. It can also be interpreted from the figure that, overall, the trajectory of the mismatched disturbance and that of the second state are symmetric about the horizontal axis 0. This phenomenon can be explained as follows. The system (38) formulation implies that, to stabilize the system state x 1 ( t ) to the equilibrium point, the following inequality must be fulfilled:
| x ˙ 1 | = | x 2 + d 1 | δ
where δ is a sufficiently small positive number and can represent the performance indicator for evaluation. Inequality (41) indicates that for the sake of stabilization, x 2 ( t ) must counteract the time-varying effect caused by the mismatched disturbance d 1 ( t ) . Consequently, x 2 ( t ) and d 1 ( t ) should be opposite. This also explains why the system state x 1 ( t ) converges to the equilibrium, but not x 2 ( t ) .
Figure 3a records the control signals generated by ESOSMC and the proposed algorithm. It clearly shows that both control signals are smooth without chattering and exhibit the same magnitude of peaking phenomenon at t = 5 s as the mismatched disturbance is suddenly initiated. The satisfactorily chattering-attenuated control signal results from an attenuation step depicted in Figure 3b, in which the chattering signal v is filtered to create the smooth switching control signal u n , as in (22). Therefore, it demonstrates the effectiveness of the proposed algorithm in attenuating chattering effects.
Figure 4 compares ESOSMC and the proposed control algorithm in mismatched disturbance estimation performance. Under the condition of time-varying mismatched disturbances, from t = 5 s onwards, the estimation error of ESOSMC varies inside a bounded region of 0.1 around the equilibrium. In comparison, DOBTSMC can deliver significantly higher precision, with an estimation error of approximately ± 5 × 10 7 . This subsequently contributes to the higher accuracy in the system state response discussed earlier and proves the effectiveness of the proposed controller in mismatched disturbance estimation.
Figure 5 illustrates the sliding surface trajectory of the proposed control algorithm, which oscillates inside the interval [ 0.005 ; 0.005 ] . Because the sliding surface σ , defined in (19), is a function of the second state x 2 , which varies with time to suppress the influence of the mismatched disturbance d 1 , it is reasonable that the sliding surface σ is also time varying. Nonetheless, the proposed control algorithm only utilizes the sign, but not the magnitude, of the sliding surface. Thus, the oscillating trajectory of the sliding surface is believed to have little or no effect on controller performance.
To illustrate the influence of different design parameter values on the performance of the disturbance observer, three simulations were conducted with different values of the Lipschitz constant, L = 1 , 3 , 5 . The simulations were carried out in a measurement noise environment with a mean of 0 and a variance of 5 × 10 3 , and the sampling time was set to t s = 0.01 s for practicality. Moreover, the initial values of the observer state and that of the system state were set to be different from each other, d ^ 1 ( 0 ) = z 2 ( 0 ) = 0.5 , x 1 ( 0 ) = 1 . The result obtained is shown in Figure 6, which strongly agrees with the statement noted in Remark 1. Specifically, it is evident from Figure 6 that ( L = 1 ) yields the longest convergence time, yet the smallest peaking phenomenon and is least susceptible to noise. In contrast, ( L = 5 ) takes the shortest amount of time for the estimation error to converge to zero but, at the same time, produces the largest peaking phenomenon and is most susceptible to noise. The most appropriate selection is ( L = 3 ) , which represents a tradeoff between convergence time and sensitivity to noise and different initial values. Thus, it is concluded here that the greater Lipschitz constant L leads to the faster convergence of the disturbance observer, yet more sensitivity to noise and larger peaking phenomenon in the transition period, especially in case where the initial values of the disturbance observer and the system are different from each other.

4.2. Application to Electro Hydrostatic Actuator System

As presented in this section, the proposed control algorithm developed previously was applied to stabilize an Electro Hydrostatic Actuator (EHA) system.
The EHA studied in this section is illustrated in Figure 7, and the major components are listed in Table 2. In the hydraulic circuit, cylinder (1) is driven by the valve system and the hydraulic pump (7), which is actuated by the AC motor (8). The flow discrepancy between the bore chamber and the rod chamber of the cylinder is automatically compensated by the check valves (4.1) and (4.2), with the support of the pilot-operated directional valve (5). In cooperation with two pilot-operated check valves (6), the hydraulic pump regulates the fluid flow inside the hydraulic circuit and consequently controls the movement of the cylinder.
Parameters used for mathematical modeling the EHA system are defined in Table 3 below.
Referring to [52], the dynamic model of the EHA system studied here consists of force dynamics and pressure dynamics, as follows. The force dynamics of the actuator is expressed as:
m x ¨ = P 1 A 1 P 2 A 2 + f F r + δ 2
where δ 2 represents a force disturbance that may include external forces and unmodeled force dynamics; f F r denotes a lumped friction force [53,54] and is approximated as a differential function of the form:
{ f F r = b 1 x ˙ b 2 tanh ( x ˙ ) tanh ( x ˙ ) = 1 e γ 2 x ˙ 1 + e γ 2 x ˙
Governed by the continuity law [55], the pressure dynamics of the actuator are:
{ V 1 t P ˙ 1 = β e ( η V D K d r u A 1 x ˙ C L i ( P 1 P 2 ) Q c v 1 ) + δ 31 V 2 t P ˙ 2 = β e ( η V D K d r u + A 2 x ˙ + C L i ( P 1 P 2 ) + Q c v 2 ) + δ 32
where Q c v 1 and Q c v 2 are, respectively, the flow in the check valves (4.1) and (4.2), and described as:
{ Q c v 1 = x ˙ ( A 1 A 2 ) sm ( x ˙ ) Q c v 2 = x ˙ ( A 1 A 2 ) sm ( x ˙ ) sm ( ) = ( 1 + e γ 3 ) 1
By defining the system state as [ x 1 , x 2 , x 3 ] T = [ x , x ˙ , ( P 1 A 1 P 2 A 2 ) / m ] T and combining equations from (42) to (45), the total dynamic system can be expressed in a canonical form as follows:
{ x ˙ 1 = x 2 x ˙ 2 = x 3 + d 1 x ˙ 3 = f + g u + d 2
where the dynamic functions are specifically defined as:
{ f = β e m [ A 1 V 1 t ( A 1 + ( A 1 A 2 ) s m ( x 2 ) ) A 2 V 2 t ( A 2 + ( A 1 A 2 ) s m ( x 2 ) ) ] x 2 g = β e η V D K d r m ( A 1 V 1 t + A 2 V 2 t ) d 1 = 1 m ( f F r + δ 2 ) d 2 = β e m C L i ( P 1 P 2 ) ( A 1 V 1 t + A 2 V 2 t ) + A 1 δ 31 m V 1 t A 2 δ 32 m V 2 t
It can be interpreted from the equation set (46) and (47) that d 1 , which consists of the lumped friction f F r and the force disturbance δ 2 , represents a mismatched disturbance, whereas d 2 , which includes internal leakages and pressure disturbances, is a matched disturbance. Thus, the EHA system considered here is subject to both mismatched and matched disturbances. It is assumed that these disturbances satisfy Assumptions 1 and 2.
To estimate the mismatched disturbance d 1 , one can apply (34) to yield a finite-time disturbance observer as follows:
{ z ˙ 12 = v 12 + x 3 v 12 = λ 12 L 2 1 / 4 | z 12 x 2 | 3 / 4 sign ( z 12 x 2 ) μ 12 ( z 12 x 2 ) + z 22 z ˙ 22 = v 22 v 22 = λ 22 L 2 1 / 3 | z 22 v 12 | 2 / 3 sign ( z 22 v 12 ) μ 22 ( z 22 v 12 ) + z 32 z ˙ 32 = v 32 v 32 = λ 32 L 2 1 / 2 | z 32 v 22 | 1 / 2 sign ( z 32 v 22 ) μ 32 ( z 32 v 22 ) + z 42 z ˙ 42 = λ 42 L 2 sign ( z 42 v 32 ) μ 42 ( z 42 v 32 )
where z 12 = x ^ 2 , z 22 = d ^ 1 , z 32 = d ˙ ^ 1 , z 42 = d ¨ ^ 1 , and L 2 > 0 , λ j 2 , μ j 2 > 0 ( j = 1 , , 4 ) . A chattering-attenuated full-order TSMC for the EHA system can then be designed based on (36) as:
σ = x ˙ 3 + c 3 | x 3 + z 22 | α 3 sign ( x 3 + z 22 ) + c 2 | x 2 | α 2 sign ( x 2 ) + c 1 | x 1 | α 1 sign ( x 1 ) + z 32
where c i ( i = 1 , 2 , 3 ) are positive constants satisfying the requirement that the polynomial p 3 + c 3 p 2 + c 2 p + c 1 is Hurwitz, and α i ( i = 1 , 2 , 3 ) are chosen as (37). Accordingly, the control law is developed with the form:
u = g 1 ( u e q + u n )
where the equivalent control u e q is designed as:
u e q = f c 3 | x 3 + z 22 | α 3 sign ( x 3 + z 22 ) c 2 | x 2 | α 2 sign ( x 2 ) c 1 | x 1 | α 1 sign ( x 1 ) z 32
and the switching term u n is defined as in (22) and (23).
The nominal values of the system parameters were determined as in Table 4 [56].
The cylinder stroke length is 0.3 m ; thus, the working range of the actuator is set to be x 1 [ 0.15 , 0.15 ] m . For safety purposes, the pressure threshold of the relief valve is set at 150 bar , meaning that P 1 , P 2 [ 0 , 150 × 10 5 ] Pa . To ensure the control task was challenging and to examine the effectiveness of control algorithms, both matched and mismatched disturbances were deliberately assumed to be time varying throughout the simulation period. Specifically, the force disturbance δ 2 and the pressure disturbances δ 31 , δ 32 were selected to be dependent on both time and system state as follows. δ 2 = 1000 ( 0.2 + 0.5 x 1 + x 2 + sin ( t ) ) , δ 31 = 1000 ( 0.1 + 3 x 2 + sin ( t ) ) , and δ 32 = 500 ( 0.2 + x 2 + 2 sin ( t ) ) . Moreover, for practicality, measurement noise and parameter uncertainty were added to the simulation. Specifically, the position and pressure measurements were subject to noise with a mean of 0 and a covariance of { 0.0005 m , 3500 Pa , 3500 Pa } [45]. Parameters used to design the controller have nominal values, whereas the parameters used to model the EHA system include variance (~ 5 % of nominal values), as indicated in Table 4. The sampling time was set to be 1 ms .
Again, the proposed control algorithm, DOBTSMC, was compared with TSMC and ESOSMC in this case study. The design control parameters of DOBTSMC were selected as c 1 = 315.6 , c 2 = 500 , c 3 = 270 ; α i ( i = 1 , 2 , 3 ) are designed as α 1 = 1 / 3 , α 2 = 1 / 5 , α 3 = 1 / 7 . The DOB parameters were chosen as L 2 = 500 , λ 12 = 30 , λ 22 = 65 , λ 32 = 7.5 , λ 42 = 5.5 , μ 12 = 6 , μ 22 = 13 , μ 32 = 1.5 , μ 42 = 1.1 .
Figure 8 presents the position response of the EHA system with the three considered control algorithms, namely, TSMC, ESOSMC, and DOBTSMC. Predictably, TSMC fails to overcome the effects of mismatched disturbances and thus cannot stabilize the target system. It is therefore excluded from the comparison. By comparison, both ESOSMC and DOBTSMC are able to force the system state to remain in the vicinity of the equilibrium. Specifically, the position state x 1 generated by ESOSMC varies in a small, bounded region of 3.5 mm around the zero axis. In contrast, the proposed controller, DOBTSMC, provides a much higher precision of about ± 0.3 mm , which is 10 times more accurate than that of ESOSMC.
The velocity trajectory of ESOSMC and the proposed controller are shown in Figure 9, which shows agreement with the position response in Figure 8. The precision of the DOBTSMC velocity signal, at ± 2.3 mm / s , is two-fold higher than that of the ESOSMC velocity signal, at ± 5 mm / s . Moreover, a larger peaking phenomenon in the transition period is seen for ESOSMC.
Figure 10 shows the trajectory of the mismatched disturbance and the third state of the control system. It is noticeable that the mismatched disturbance is time varying and that its trajectory and x 3 trajectories are symmetric. This phenomenon is consistent with the previous numerical study, where the system state x 3 is controlled to cancel out the influence of the mismatched disturbance d 1 . Thus, x 3 and d 1 are likely to be opposite in value.
The control signal generated by ESOSMC and the proposed control algorithm are recorded in Figure 11. Overall, both signals are smooth and show few chattering effects. The ESOSMC control signal oscillates inside ± 0.5 V , whereas the figure for DOBTSMC is slightly better, at ± 0.3 V . However, ESOSMC exhibits a larger peaking phenomenon than the proposed controller in the transition period.
Figure 12 illustrates the performance of the ESOSMC and that of the proposed controller in estimating mismatched disturbances. The estimation error of DOBTSMC, at about ± 0.6 m / s 2 , is approximately two-fold smaller than that of the ESOSMC, at about ± 1.2 m / s 2 . Both algorithms experience chattering phenomena in the trajectory. This can be explained as follows. Because both schemes utilize system states in observer design, measurement noise in system states has a direct impact on the performance of the observation process. Thus, chattering caused by noise appears and consequently degrades the estimation accuracy of both controllers.
Finally, working pressures of the proposed controller are shown in Figure 13, which clearly implies the nonlinearity of the EHA system. Specifically, the bore chamber area A 1 is larger than the rod chamber A 2 . Thus, to produce the same level of force, the bore chamber pressure P 1 must be accordingly smaller than the rod chamber pressure P 2 . This special condition makes the task of controlling a single rod cylinder much more challenging than the task of controlling a double rod cylinder, where the chamber areas are identical.

5. Conclusions

In this work, a new disturbance observer-based terminal sliding mode control was proposed to control a class of nonlinear systems under the influence of both matched and mismatched disturbances. Although matched disturbances were suppressed by the robust nature of the terminal sliding mode control, mismatched disturbances were compensated by the enhancement in disturbance observation. The employment of a nonlinear disturbance observer enabled the mismatched disturbances to be estimated, and then compensated, by the proposed control algorithm. A general guideline on the appropriate selection of parameters for this observer was provided, and the influence of parameter selection was later illustrated by simulations. Other aspects of the proposed controller, namely, chattering attenuation, sign function utilization, and nominal control performance preservation, were also discussed and demonstrated in this work. Although the proposed control scheme was first designed to stabilize a class of nonlinear second-order systems, it was also shown to be extendable to an nth-order formulation for wider applications. The formulation extension was subsequently used to control a third-order electro hydrostatic actuator system, which, along with a numerical example, revealed the characteristics and verified the effectiveness of the proposed algorithm. Comparative simulation results ultimately demonstrated the satisfactory performance of the proposed controller.
Throughout the study, system states were assumed to be available, and thus different sensors were required. However, there are numerous cases in practice in which only one sensor is implemented for several reasons, such as reducing the installation cost and space. Therefore, future research can be undertaken to design a disturbance-state observer to estimate both disturbances and unmeasured states in a fast-changing, noisy measurement environment.

Author Contributions

K.K.A. was the supervisor and administrator of the project and reviewed the manuscript. D.G.N. undertook the theoretical conceptualization, carried out simulation studies and examined the obtained results. D.T.T. checked the manuscript and provided valuable insights into the project. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Basic Science Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science and ICT, South Korea, under Grant NRF 2020R1A2B5B03001480.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
SMCSliding Mode Control
TSMCTerminal Sliding Mode Control
MAGLEVMagnetic Levitation
PMSMPermanent Magnet Synchronous Motor
EHAElectro Hydrostatic Actuator
DOBDisturbance Observer
ESOSMCExtended State Observer based Sliding Mode Control
DOBTSMCThe proposed controller

References

  1. Ginoya, D.; Shendge, P.; Phadke, S. Disturbance observer based sliding mode control of nonlinear mismatched uncertain systems. Commun. Nonlinear Sci. Numer. Simul. 2015, 26, 98–107. [Google Scholar] [CrossRef]
  2. Jing, C.; Xu, H.; Niu, X. Adaptive sliding mode disturbance rejection control with prescribed performance for robotic manipulators. ISA Trans. 2019, 91, 41–51. [Google Scholar] [CrossRef]
  3. Asl, R.M.; Hagh, Y.S.; Palm, R. Robust control by adaptive Non-singular Terminal Sliding Mode. Eng. Appl. Artif. Intell. 2017, 59, 205–217. [Google Scholar] [CrossRef]
  4. Ye, H.; Wang, S. Trajectory Tracking Control for Nonholonomic Wheeled Mobile Robots with External Disturbances and Parameter Uncertainties. Int. J. Control Autom. Syst. 2020, 18, 3015–3022. [Google Scholar] [CrossRef]
  5. Zhao, Y.; Zhang, Y.; Lee, J. Lyapunov and Sliding Mode Based Leader-follower Formation Control for Multiple Mobile Robots with an Augmented Distance-angle Strategy. Int. J. Control Autom. Syst. 2019, 17, 1314–1321. [Google Scholar] [CrossRef]
  6. Jin, Z.; Liang, Z.; Wang, X.; Zheng, M. Adaptive Backstepping Sliding Mode Control of Tractor-trailer System with Input Delay Based on RBF Neural Network. Int. J. Control Autom. Syst. 2020, 19, 76–87. [Google Scholar] [CrossRef]
  7. Napole, C.; Barambones, O.; Derbeli, M.; Calvo, I. Advanced Trajectory Control for Piezoelectric Actuators Based on Robust Control Combined with Artificial Neural Networks. Appl. Sci. 2021, 11, 7390. [Google Scholar] [CrossRef]
  8. Gao, Y.; Wu, Y.; Wang, X.; Chen, Q. Characteristic Model-based Adaptive Control with Genetic Algorithm Estimators for Four-PMSM Synchronization System. Int. J. Control Autom. Syst. 2020, 18, 1605–1616. [Google Scholar] [CrossRef]
  9. Wang, J.; Zhao, L.; Yu, L. Reduced-order Generalized Proportional Integral Observer Based Continuous Dynamic Sliding Mode Control for Magnetic Levitation System with Time-varying Disturbances. Int. J. Control Autom. Syst. 2020, 19, 439–448. [Google Scholar] [CrossRef]
  10. Din, S.U.; Khan, Q.; Rehman, F.U.; Akmeliawanti, R. A Comparative Experimental Study of Robust Sliding Mode Control Strategies for Underactuated Systems. IEEE Access 2017, 6, 1927–1939. [Google Scholar] [CrossRef]
  11. Derbeli, M.; Barambones, O.; Silaa, M.; Napole, C. Real-Time Implementation of a New MPPT Control Method for a DC-DC Boost Converter Used in a PEM Fuel Cell Power System. Actuators 2020, 9, 105. [Google Scholar] [CrossRef]
  12. Napole, C.; Derbeli, M.; Barambones, O. A global integral terminal sliding mode control based on a novel reaching law for a proton exchange membrane fuel cell system. Appl. Energy 2021, 301, 117473. [Google Scholar] [CrossRef]
  13. Chen, Z.; Tu, X.; Xing, L.; Fu, J.; Lozano, R.; Rogelio, L. A Special Kind of Sliding Mode Control for Nonlinear System with State Constraints. IEEE Access 2019, 7, 69998–70010. [Google Scholar] [CrossRef]
  14. Kamalesh, M.; Senthilnathan, N.; Bharatiraja, C. Design of a Novel Boomerang Trajectory for Sliding Mode Controller. Int. J. Control Autom. Syst. 2020, 18, 2917–2928. [Google Scholar] [CrossRef]
  15. Ye, D.; Zhang, H.; Tian, Y.; Zhao, Y.; Sun, Z. Fuzzy Sliding Mode Control of Nonparallel-ground-track Imaging Satellite with High Precision. Int. J. Control Autom. Syst. 2020, 18, 1617–1628. [Google Scholar] [CrossRef]
  16. Xu, S.S.-D.; Chen, C.-C.; Wu, Z.-L. Study of Nonsingular Fast Terminal Sliding-Mode Fault-Tolerant Control. IEEE Trans. Ind. Electron. 2015, 62, 3906–3913. [Google Scholar] [CrossRef]
  17. Feng, Y.; Zhou, M.; Zheng, X.; Han, F.; Yu, X. Full-order terminal sliding-mode control of MIMO systems with unmatched uncertainties. J. Frankl. Inst. 2018, 355, 653–674. [Google Scholar] [CrossRef]
  18. Levant, A. Chattering Analysis. IEEE Trans. Autom. Control 2010, 55, 1380–1389. [Google Scholar] [CrossRef]
  19. Kachroo, P.; Tomizuka, M. Chattering reduction and error convergence in the sliding-mode control of a class of nonlinear systems. IEEE Trans. Autom. Control 1996, 41, 1063–1068. [Google Scholar] [CrossRef] [Green Version]
  20. Bartolini, G.; Ferrara, A.; Usai, E. Chattering avoidance by second-order sliding mode control. IEEE Trans. Autom. Control 1998, 43, 241–246. [Google Scholar] [CrossRef]
  21. Levant, A. Higher-order sliding modes, differentiation and output-feedback control. Int. J. Control 2003, 76, 924–941. [Google Scholar] [CrossRef]
  22. Levant, A. Principles of 2-sliding mode design. Automatica 2007, 43, 576–586. [Google Scholar] [CrossRef] [Green Version]
  23. Xu, J.-X.; Pan, Y.-J.; Lee, T.-H. Sliding Mode Control with Closed-Loop Filtering Architecture for a Class of Nonlinear Systems. IEEE Trans. Circuits Syst. II Express Briefs 2004, 51, 168–173. [Google Scholar] [CrossRef]
  24. Tseng, M.-L.; Chen, M.-S. Chattering reduction of sliding mode control by low-pass filtering the control signal. Asian J. Control 2010, 12, 392–398. [Google Scholar] [CrossRef]
  25. Kawamura, A.; Itoh, H.; Sakamoto, K. Chattering reduction of disturbance observer based sliding mode control. IEEE Trans. Ind. Appl. 1994, 30, 456–461. [Google Scholar] [CrossRef]
  26. Mu, C.; He, H. Dynamic Behavior of Terminal Sliding Mode Control. IEEE Trans. Ind. Electron. 2017, 65, 3480–3490. [Google Scholar] [CrossRef]
  27. Corradini, M.L.; Cristofaro, A. Nonsingular terminal sliding-mode control of nonlinear planar systems with global fixed-time stability guarantees. Automatica 2018, 95, 561–565. [Google Scholar] [CrossRef]
  28. Rabiee, H.; Ataei, M.; Ekramian, M. Continuous nonsingular terminal sliding mode control based on adaptive sliding mode disturbance observer for uncertain nonlinear systems. Automatica 2019, 109, 108515. [Google Scholar] [CrossRef]
  29. Wang, Y.; Zhu, K.; Yan, F.; Chen, B. Adaptive super-twisting nonsingular fast terminal sliding mode control for cable-driven manipulators using time-delay estimation. Adv. Eng. Softw. 2018, 128, 113–124. [Google Scholar] [CrossRef]
  30. Feng, Y.; Han, F.; Yu, X. Chattering free full-order sliding-mode control. Automatica 2014, 50, 1310–1314. [Google Scholar] [CrossRef]
  31. Wang, J.; Li, S.; Yang, J.; Wu, B.; Li, Q. Finite-time disturbance observer based non-singular terminal sliding-mode control for pulse width modulation based DC–DC buck converters with mismatched load disturbances. IET Power Electron. 2016, 9, 1995–2002. [Google Scholar] [CrossRef]
  32. Utkin, V.; Shi, J. Integral sliding mode in systems operating under uncertainty conditions. In Proceedings of the 35th IEEE Conference on Decision and Control, Las Vegas, NV, USA, 13 December 2002; Volume 4, pp. 4591–4596. [Google Scholar] [CrossRef]
  33. Cao, W.-J.; Xu, J.-X. Nonlinear Integral-Type Sliding Surface for Both Matched and Unmatched Uncertain Systems. IEEE Trans. Autom. Control 2004, 49, 1355–1360. [Google Scholar] [CrossRef]
  34. Castanos, F.; Fridman, L. Analysis and Design of Integral Sliding Manifolds for Systems With Unmatched Perturbations. IEEE Trans. Autom. Control 2006, 51, 853–858. [Google Scholar] [CrossRef] [Green Version]
  35. Mondal, S.; Mahanta, C. Chattering free adaptive multivariable sliding mode controller for systems with matched and mismatched uncertainty. ISA Trans. 2013, 52, 335–341. [Google Scholar] [CrossRef]
  36. Kwan, C.-M. Sliding mode control of linear systems with mismatched uncertainties. Automatica 1995, 31, 303–307. [Google Scholar] [CrossRef]
  37. Wen, C.-C.; Cheng, C.-C. Design of sliding surface for mismatched uncertain systems to achieve asymptotical stability. J. Frankl. Inst. 2008, 345, 926–941. [Google Scholar] [CrossRef]
  38. Yang, J.; Li, S.; Yu, X. Sliding-Mode Control for Systems with Mismatched Uncertainties via a Disturbance Observer. IEEE Trans. Ind. Electron. 2012, 60, 160–169. [Google Scholar] [CrossRef]
  39. Ginoya, D.; Shendge, P.D.; Phadke, S. Sliding Mode Control for Mismatched Uncertain Systems Using an Extended Disturbance Observer. IEEE Trans. Ind. Electron. 2013, 61, 1983–1992. [Google Scholar] [CrossRef]
  40. Zhang, J.; Liu, X.; Xia, Y.; Zuo, Z.; Wang, Y. Disturbance Observer-Based Integral Sliding-Mode Control for Systems with Mismatched Disturbances. IEEE Trans. Ind. Electron. 2016, 63, 7040–7048. [Google Scholar] [CrossRef]
  41. Shi, S.-L.; Li, J.-X.; Fang, Y. Extended-State-Observer-Based Chattering Free Sliding Mode Control for Nonlinear Systems with Mismatched Disturbance. IEEE Access 2018, 6, 22952–22957. [Google Scholar] [CrossRef]
  42. Nguyen, N.P.; Oh, H.; Kim, Y.; Moon, J. Disturbance Observer-Based Continuous Finite-Time Sliding Mode Control against Matched and Mismatched Disturbances. Complexity 2020, 2020, 1–14. [Google Scholar] [CrossRef]
  43. Chen, W.-H. Nonlinear Disturbance Observer-Enhanced Dynamic Inversion Control of Missiles. J. Guid. Control Dyn. 2003, 26, 161–166. [Google Scholar] [CrossRef]
  44. Su, J.; Yang, J.; Li, S. Continuous finite-time anti-disturbance control for a class of uncertain nonlinear systems. Trans. Inst. Meas. Control 2013, 36, 300–311. [Google Scholar] [CrossRef]
  45. Yang, J.; Li, S.; Su, J.; Yu, X. Continuous nonsingular terminal sliding mode control for systems with mismatched disturbances. Automatica 2013, 49, 2287–2291. [Google Scholar] [CrossRef] [Green Version]
  46. Yang, J.; Su, J.; Li, S.; Yu, X. High-Order Mismatched Disturbance Compensation for Motion Control Systems Via a Continuous Dynamic Sliding-Mode Approach. IEEE Trans. Ind. Inform. 2013, 10, 604–614. [Google Scholar] [CrossRef] [Green Version]
  47. Bhat, S.P.; Bernstein, D.S. Geometric homogeneity with applications to finite-time stability. Math. Control Signals Syst. 2005, 17, 101–127. [Google Scholar] [CrossRef] [Green Version]
  48. Feng, Y.; Han, F.; Yu, X. Reply to “Comments on ‘Chattering free full-order sliding-mode control’ [Automatica 50 (2014) 1310–1314]”. Automatica 2016, 72, 255–256. [Google Scholar] [CrossRef]
  49. Bhat, S.; Bernstein, D. Finite-time stability of homogeneous systems. In Proceedings of the 1997 American Control Conference (Cat. No.97CH36041), Albuquerque, NM, USA, 6 June 1997. [Google Scholar] [CrossRef]
  50. Levant, A.; Livne, M. Globally convergent differentiators with variable gains. Int. J. Control 2018, 91, 1994–2008. [Google Scholar] [CrossRef]
  51. Levant, A. Homogeneity approach to high-order sliding mode design. Automatica 2005, 41, 823–830. [Google Scholar] [CrossRef]
  52. Ba, D.X.; Dinh, T.Q.; Bae, J.; Ahn, K.K. An Effective Disturbance-Observer-Based Nonlinear Controller for a Pump-Controlled Hydraulic System. IEEE/ASME Trans. Mechatron. 2019, 25, 32–43. [Google Scholar] [CrossRef]
  53. Yao, J.; Jiao, Z.; Ma, D.; Yan, L. High-Accuracy Tracking Control of Hydraulic Rotary Actuators with Modeling Uncertainties. IEEE/ASME Trans. Mechatron. 2013, 19, 633–641. [Google Scholar] [CrossRef]
  54. Won, D.; Kim, W.; Shin, D.; Chung, C.C. High-Gain Disturbance Observer-Based Backstepping Control With Output Tracking Error Constraint for Electro-Hydraulic Systems. IEEE Trans. Control Syst. Technol. 2014, 23, 787–795. [Google Scholar] [CrossRef]
  55. Merritt, H.E.; Pomper, V. Hydraulic control systems. J. Appl. Mech. 1968, 35, 200. [Google Scholar] [CrossRef]
  56. Ba, D.X.; Ahn, K.K.; Truong, D.Q.; Park, H.G. Integrated model-based backstepping control for an electro-hydraulic system. Int. J. Precis. Eng. Manuf. 2016, 17, 565–577. [Google Scholar] [CrossRef]
Figure 1. Response trajectory of the state x 1
Figure 1. Response trajectory of the state x 1
Applsci 11 08158 g001
Figure 2. Mismatched disturbance and second state response trajectories.
Figure 2. Mismatched disturbance and second state response trajectories.
Applsci 11 08158 g002
Figure 3. Total control signal (a) and switching control signal (b).
Figure 3. Total control signal (a) and switching control signal (b).
Applsci 11 08158 g003
Figure 4. Mismatched disturbance estimation error of ESOSMC and DOBTSMC.
Figure 4. Mismatched disturbance estimation error of ESOSMC and DOBTSMC.
Applsci 11 08158 g004
Figure 5. Sliding surface response of the proposed control algorithm.
Figure 5. Sliding surface response of the proposed control algorithm.
Applsci 11 08158 g005
Figure 6. Influence of design parameters on observer performance.
Figure 6. Influence of design parameters on observer performance.
Applsci 11 08158 g006
Figure 7. Schematic diagram of the studied EHA system.
Figure 7. Schematic diagram of the studied EHA system.
Applsci 11 08158 g007
Figure 8. Position response of the EHA system.
Figure 8. Position response of the EHA system.
Applsci 11 08158 g008
Figure 9. Velocity response of the EHA system.
Figure 9. Velocity response of the EHA system.
Applsci 11 08158 g009
Figure 10. Mismatched disturbance and acceleration signal of the EHA system.
Figure 10. Mismatched disturbance and acceleration signal of the EHA system.
Applsci 11 08158 g010
Figure 11. Control signal trajectory of ESOSMC and DOBTSMC.
Figure 11. Control signal trajectory of ESOSMC and DOBTSMC.
Applsci 11 08158 g011
Figure 12. Mismatched disturbance estimation error of ESOSMC and DOBTSMC.
Figure 12. Mismatched disturbance estimation error of ESOSMC and DOBTSMC.
Applsci 11 08158 g012
Figure 13. Working pressures of the proposed controller.
Figure 13. Working pressures of the proposed controller.
Applsci 11 08158 g013
Table 1. Comparison of the proposed controller and other works.
Table 1. Comparison of the proposed controller and other works.
Control
Algorithm
Matched
Disturbance
Compensation
Mismatched
Disturbance
Compensation
Finite-Time
Convergence
TSMC [30]YesNoYes
EOSMC [41]YesYesNo
The proposedYesYesYes
Table 2. Major components in the considered EHA system.
Table 2. Major components in the considered EHA system.
No.DeviceNo.Device
1Hydraulic cylinder5Directional valve
2Relief valves6Pilot check valve
3Reservoir7Hydraulic pump
4.1Check valve 4.18AC motor
4.2Check valve 4.29AC motor driver
Table 3. Major components in the considered EHA system.
Table 3. Major components in the considered EHA system.
SymbolQuantity and Unit
x Displacement of the cylinder. ( m )
P 1 , P 2 Pressures inside the cylinder chambers. ( Pa )
A 1 , A 2 Bore- and rod- side areas of the cylinder. ( m 2 )
m Relative mass of the motion system. ( kg )
V 10 , V 20 Initial volumes of chamber 1 and 2. ( m 3 )
V 1 t = V 10 + A 1 x Active volume of chamber 1. ( m 3 )
V 2 t = V 20 A 2 x Active volume of chamber 2. ( m 3 )
β e Effective bulk modulus of the fluid. ( Pa )
C L i Internal leakage coefficient. ( m 3 / ( sPa ) )
D Displacement of the pump. ( m 3 / rad )
K d r Control gain of the motor driver. ( rad / ( sV ) )
η V Volumetric efficiency of the pump.
u Driving voltage of the motor driver. ( V )
b 1 , b 2 , γ 2 , γ 3 Positive constants in the system model.
δ 2 Lumped disturbance in force dynamics. ( N )
δ 31 , δ 32 Lumped disturbances in pressure dynamics. ( N / s )
Table 4. Nominal values of the EHA system parameters.
Table 4. Nominal values of the EHA system parameters.
ParameterUnitNominal (+ Variance) Value
m kg 9.32 ( + 0.46 )
A 1 m 2 1.9635 × 10 3 ( + 0.0982 × 10 3 )
A 2 m 2 1.2566 × 10 3 ( + 0.0628 × 10 3 )
V 10 m 3 4.375 × 10 4 ( + 0.219 × 10 4 )
V 20 m 3 2.88 × 10 4 ( + 0.144 × 10 4 )
β e Pa 5.34 × 10 8
K d r rad / ( sV ) 10 π
D 1 η V D m 3 / rad 5.83 × 10 7 ( + 0.29 × 10 7 )
C L i m 3 / ( sPa ) 1.51 × 10 15
b 1 Ns / m 258
b 2 N 532
γ 2 10
γ 3 15 ( + 0.75 )
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nguyen, D.G.; Tran, D.T.; Ahn, K.K. Disturbance Observer-Based Chattering-Attenuated Terminal Sliding Mode Control for Nonlinear Systems Subject to Matched and Mismatched Disturbances. Appl. Sci. 2021, 11, 8158. https://doi.org/10.3390/app11178158

AMA Style

Nguyen DG, Tran DT, Ahn KK. Disturbance Observer-Based Chattering-Attenuated Terminal Sliding Mode Control for Nonlinear Systems Subject to Matched and Mismatched Disturbances. Applied Sciences. 2021; 11(17):8158. https://doi.org/10.3390/app11178158

Chicago/Turabian Style

Nguyen, Duc Giap, Duc Thien Tran, and Kyoung Kwan Ahn. 2021. "Disturbance Observer-Based Chattering-Attenuated Terminal Sliding Mode Control for Nonlinear Systems Subject to Matched and Mismatched Disturbances" Applied Sciences 11, no. 17: 8158. https://doi.org/10.3390/app11178158

APA Style

Nguyen, D. G., Tran, D. T., & Ahn, K. K. (2021). Disturbance Observer-Based Chattering-Attenuated Terminal Sliding Mode Control for Nonlinear Systems Subject to Matched and Mismatched Disturbances. Applied Sciences, 11(17), 8158. https://doi.org/10.3390/app11178158

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

Article Metrics

Back to TopTop