Next Article in Journal
On a Certain Functional Equation and Its Application to the Schwarz Problem
Previous Article in Journal
An Ensemble Deep Learning Model for Provincial Load Forecasting Based on Reduced Dimensional Clustering and Decomposition Strategies
Previous Article in Special Issue
LCANet: A Lightweight Context-Aware Network for Bladder Tumor Segmentation in MRI Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Finite-Time Sliding Mode Control Approach for Constrained Euler–Lagrange System

School of Astronautics, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(12), 2788; https://doi.org/10.3390/math11122788
Submission received: 12 May 2023 / Revised: 15 June 2023 / Accepted: 19 June 2023 / Published: 20 June 2023
(This article belongs to the Special Issue Advanced Control Theory with Applications)

Abstract

:
This paper investigates a general control strategy to track the reference trajectory for the constrained Euler–Lagrange system with model uncertainties and unknown external disturbances. Unlike the disturbances assumed to be upper-bounded by a constant in other papers, we consider the disturbances to be bounded by a function of the system states, which are more realistic. First, the nominal controller was designed based on the nonsingular fast terminal sliding mode control, and global fast finite-time convergence to the sliding surface was guaranteed. As the system is state-constrained in this paper, we introduce the control barrier function approach to formulate the constraints and ensure the system does not break the restrictions. The proposed control strategy was numerically assessed on a two-link robot manipulator system, and the simulation results illustrate the effectiveness of the proposed control strategy.

1. Introduction

The Euler–Lagrange system, as a mechanical system that is widely used in robotics, such as unmanned ground vehicles [1,2], unmanned air vehicles [3,4], robot manipulators [5,6], satellites [7], and so on, has been of broad interest in recent years. The system model is derived from the kinetic equations and the Euler–Lagrange equation. Consequently, it is worthy of further research, especially in practical applications. Moreover, as a typical motion control system, controller design based on the Euler–Lagrange system model is more versatile and dumbs down the designing process. Nonetheless, the sophistication of the dynamics makes it challenging to identify the system parameters for a specific robot. Given the vital importance of the linear in parameters (LIP) property [8] of the Euler–Lagrange system, estimating the uncertainties is achievable via an appropriate observer [1,7]. In addition, the Euler–Lagrange system is highly nonlinear and constantly affected by various external unknown disturbances. Hence, for the trajectory tracking control of the Euler–Lagrange system, researchers have introduced various advanced control strategies to solve the above issues in the actual control scenario, such as model predictive control [9,10] and sliding mode control [11].
As more attention has been paid to swarm behavior in insects and birds, approaches involving individual behavior have certain limitations, especially in complex task scenarios. Through the cooperation between multiple agents, work efficiency and reliability can be improved. Even if a single or part of the agents fails, the system can remain functional through appropriate reactions. Most agents can be modeled as a Euler–Lagrange system, even for a heterogeneous multiagent system. Recently, researchers have increasingly focused on multiagent cooperative control based on the Euler–Lagrange system [12,13].
Unlike classical control theory, which is drawn from engineering practice and focuses on the system’s dynamic performance, the modern control approach is based on stability analysis theory, i.e., designing the controller to guarantee the system converges to the equilibrium point. In practice, the rapid convergence of the system has a very significant impact on the control performance, especially in environments that require high dynamics.
Sliding mode control (SMC) is one of the most effective control approaches as a result of its insensitivity to certain external disturbances and model uncertainties. The basic idea of SMC is to steer the trajectory of the system to the sliding surface and converge to the equilibrium point. However, as a linear sliding surface is adopted, it can only guarantee the asymptotic stability of the system. Compared with the conventional sliding mode control, terminal sliding mode control (TSMC) can drive the system to the equilibrium point in finite time via a nonlinear function in the sliding surface. In [14], a continuous TSMC was tested for mechanical servo systems, which achieved robustness and finite-time convergence. Although TSMC could significantly improve the dynamic performance of the system, the singularity problem remained. A new terminal sliding surface [15] was proposed to eliminate the singularity problem. For further improvement on the convergence rate of the system, considering the different behavior of linear and nonlinear functions in different intervals, researchers proposed many improvement methods for the sliding mode surface. A robust ATNTSM control approach was introduced in [16] for robotic manipulators with finite-time convergence guaranteed based on the arctan function. It was shown that the convergence rate is apparently faster than the traditional TSMC and NTSMC. An adaptive nonsingular fast TSMC was proposed in [17] to control both the position and attitude of a quadrotor subject to model uncertainties and external disturbances.
In addition, there always exist various constraints for the agent to achieve complex tasks working in a restricted and practical environment, such as adaptive cruise control [18] and collision avoidance [19]. Model predictive control [20] is widely used for systems with constraints. However, model predictive control is a model-based optimal control approach, which means the model accuracy will affect the control performance directly. Due to the high calculation burden, more computing resources are required to solve the nonlinear optimization problem, especially in high dynamics environments. Another control approach is based on the barrier Lyapunov function. Unlike quadratic functions, the barrier will force the system to back away as it approaches the state’s boundary. In [21], an adaptive fuzzy controller was proposed based on a tan-type barrier Lyapunov function to achieve trajectory tracking with output constraints. Similarly, a tan-type barrier Lyapunov-based optimized backstepping controller was proposed in [22]. While the desired signal does not break the constraints, it is hard to say whether the state constraints can be satisfied in practice. A mixed H / H fault-detection filtering was proposed in [23] to improve the safety and reliability of industrial processes. An observer-based fault estimation and fault-tolerant control of the linear discrete time-varying stochastic system was presented in [24].
In recent years, a safety-critical control strategy was proposed in [25] based on the set invariance theory to drive the system states into a predefined safe set. This work unifies the control Lyapunov function (CLF) with the control barrier function (CBF) through an online quadratic programming (QP) process for adaptive cruise controller design. It guarantees that the car can always go forward with the desired velocity unless the front vehicle is too close to maintain a safe distance. It has been shown that the online QP can be performed with sample rates of 200 Hz to 1 kHz in real time, which means the CLF-CBF-QP framework is an effective way to stabilize the system with additional constraints, such as state and output constraints. Similar frameworks have also been extended to various attractive applications in robotics, such as collision-free behaviors in multirobot systems [26,27], dynamic locomotion over rough terrain for quadruped robots [28], and bipedal robotic walking [29].
In this paper, we tackle the tracking problem of the constrained Euler–Lagrange system with model uncertainties and unknown external disturbances. The main contributions of this paper are summed up as follows: (i) To improve the control performance of the conventional NTSMC, we propose a novel nonsingular fast terminal sliding mode controller (NFTSMC). Compared with the algorithm in [16], the NFTSMC shows a faster convergence rate and better robustness; (ii) the state constraints are formulated as a function-defined safe set. The tracking problem is solved using a CBF-QP, and the simulation results illustrate the effectiveness of tracking a reference signal that sometimes breaks the limits on states.
The rest of this paper is organized as follows. The Euler–Lagrange system model and some preliminaries are given in Section 2. The NFTSMC and CBF-QP control strategy is presented in Section 3, and simulation results illustrate the effectiveness of the proposed control strategy in Section 4. The paper is concluded in Section 5.

2. Dynamic Model and Preliminaries

2.1. Dynamic Model

Consider a typical Euler–Lagrange system
M ( q ¨ ) + C ( q , q ˙ ) q ˙ + G ( q ) = τ + τ d ,
where q , q ˙ , q ¨ R n denote the generalized position, velocity, and acceleration of the system, respectively. M ( q ) R n × n , which is symmetric positive definite denotes the inertia matrix and C ( q , q ˙ ) R n × n , is the centrifugal and Coriolis matrix. G ( q ) R n is a generalized force caused by the potential energy. τ , τ d R n are generalized control input and external disturbances applied to the system, respectively. For the convenience, the matrices are abbreviated as M , C , G . In fact, we cannot obtain the system parameters precisely, i.e., there exists uncertain terms δ M , δ C , and δ G , which can be expressed as
M = M 0 + δ M , C = C 0 + δ C , G = G 0 + δ G ,
where M 0 , C 0 , and G 0 are estimated system parameters. Then, system (1) can be reformed as
M 0 q ¨ + C 0 q ˙ + G 0 = τ + f ,
where f = τ d δ M q ¨ δ C q ˙ δ G is the total disturbances including model uncertainties and external disturbances.
Let x 1 = q , x 2 = q ˙ , then system (2) can be rewritten as
x ˙ 1 = x 2 , x ˙ 2 = F ( x 1 , x 2 ) + M 0 1 ( x ) τ + d ( t , x ) ,
where
F ( x 1 , x 2 ) = M 0 1 ( C 0 q ˙ + G 0 ) , d ( t , x ) = M 0 1 f ,
The external disturbances such as friction and wind drag are bounded [15] and can be expressed as a function of the states. Then, we have the following assumption:
Assumption 1.
The disturbance d ( t , x ) is bounded by
d ( t , x ) M 0 1 b 0 + b 1 q + b 2 q ˙
where b 0 , b 1 , and b 2 are positive constants. Let ρ = M 0 1 b 0 + b 1 q + b 2 q ˙ . It is trivial that ρ > 0 .
The following lemma is given to analyze the stability of the system:
Lemma 1.
The states of the system will converge to the origin within finite time if there exists a continuous positive function V ( x ) : R n R for any given initial condition x ( 0 ) = x 0 , such that [30]
V ˙ ( x ) α 1 V ( x ) α 2 V α ( x ) ,
where α 1 , α 2 > 0 , 0 < α < 1 .

2.2. Conventional Nonsingular TSMC

Consider the following second-order nonlinear system with uncertainties:
x ˙ 1 = x 2 , x ˙ 2 = f ( x ) + g ( x ) u + d ( x , t ) ,
where d ( x , t ) denotes the uncertainties and external disturbances, which is bounded by d ( x , t ) D , g ( x ) 0 , and D is a positive constant. Let x = [ x 1 , x 2 ] T be the state vector with x 1 R and x 2 R . The nonsingular terminal sliding surface is chosen as
s = x 1 + 1 β x 2 p q ,
where β > 0 and p, q are positive odd integers satisfying p > q . For system (6), the nonsingular terminal sliding mode controller is designed as
u = g 1 ( x ) f ( x ) + β q p x 2 2 p q + ( D + η ) sign ( s ) ,
where 1 < p / q < 2 , η > 0 . For convenience, let λ = 1 / β , γ = q / p . Then, when the sliding mode motion is achieved, the states will converge to the equilibrium point within
T n = 0 | x 1 ( 0 ) | λ 1 / γ x 1 1 / γ d x 1 = λ 1 / γ · | x 1 ( 0 ) | 1 1 / γ 1 1 / γ .

2.3. Control Barrier Function

The basic idea of the control barrier function method is to design a barrier function that satisfies some specific conditions and keeps the safe set forward invariant.
Definition 1.
The set C is the aforementioned forward invariant for system (10) if, from any initial state x ( t 0 ) C , the solutions x ( t ) remain in that set for all t t 0 .
Consider a general control affine system,
x ˙ = f ( x ) + g ( x ) u
with f and g being locally Lipschitz, and x X R n being the state and u U R m the control input. A safe set C is defined as
C = { x R n : h ( x ) 0 } ,
where h ( x ) : R n R is a continuously differentiable function and will be defined later. The safe set C is forward invariant if and only if there exists a control input u, such that
h ˙ ( x , u ) α ( h ( x ) ) ,
where α ( · ) is an extended class K function that is strictly increasing and α ( 0 ) = 0 . There are many other equivalent conditions; for more details see [31] and the references therein. Then, we can have the following definition of control barrier function:
Definition 2.
The function h ( x ) is said to be a zero control barrier function if there exists an extended class K function α such that [32]
sup u R m ( L f h ( x ) + L g h ( x ) u + α ( h ( x ) ) ) 0 ,
where L f h ( x ) and L g h ( x ) denote the Lie derivatives of h ( x ) as
L f h ( x ) = h ( x ) x f ( x ) , L g h ( x ) = h ( x ) x g ( x ) .
One of the benefits of CBF is that we can design the controller separately. Given the nominal input u norm , which is intended to stabilize or tracks a reference signal for system (10), a CBF-QP problem can be expressed as
u * ( x ) = arg min u R m 1 2 u u norm 2 s . t . L f h ( x ) + L g h ( x ) u + α ( h ( x ) ) 0 .
The applied control input u is equivalent to the nominal one if the safe condition is satisfied. Once the system become potentially unsafe, the CBF-QP guarantees that the system will not break the safe set and that the applied control input is as close as possible to the nominal control.

3. Controller Design

In this section, an NFTSMC controller design method is given such that the Euler–Lagrange system can track a given reference trajectory within finite time with a strong antidisturbance capacity. To overcome the singularity problem in conventional TSMC design and further improve the convergence rate, inspired by [33], for the cases of relative degree higher than 1, an exponential CBF (ECBF) is used to ensure the safe set C forward invariant.

3.1. NFTSM Controller Design

To achieve fast convergence and disturbance rejection, the nonsingular terminal sliding surface was designed as
s = x 1 + Λ 1 sign Γ 1 ( x 1 ) + Λ 2 sign Γ 2 ( x 2 ) ,
where
Γ 1 = diag ( γ 11 , γ 12 , , γ 1 n ) , Γ 2 = diag ( γ 21 , γ 22 , , γ 2 n ) Λ 1 = diag ( λ 11 , λ 12 , , λ 1 n ) , Λ 2 = diag ( λ 21 , λ 22 , , λ 2 n )
with γ 1 i > γ 2 i , 1 < γ 2 i < 2 and λ 1 i , λ 2 i > 0 for i = 1 , 2 , , n . Note that sign Γ ( x ) denotes
sign Γ ( x ) = diag ( sign ( x ) ) · | x | Γ = [ | x 1 | γ 1 sign ( x 1 ) | x n | γ n sign ( x n ) ] T ,
where sign ( x ) = [ sign ( x 1 ) sign ( x n ) ] T , | x | Γ = [ | x 1 | γ 1 | x n | γ n ] T and it is trivial that
d d t sign Γ ( x ) = Γ diag ( | x | Γ I n ) · x ˙ ,
where I n is an n-dimensional identity matrix. Then, we have the time derivative of (15) as
s ˙ = x 2 + Λ 1 Γ 1 diag ( | x 1 | Γ 1 I n ) x 2 + Λ 2 Γ 2 diag ( | x 2 | Γ 2 I n ) x ˙ 2 .
Theorem 1.
Consider the rewritten Euler–Lagrange system (3) with uncertainties and disturbances satisfying Assumption 1, the system state will converge to the origin in finite time via the following control law:
τ = M 0 ( M 2 s + ( ρ + M 1 ) s s M 0 1 ( C 0 x 2 + G 0 ) + Λ 2 1 Γ 2 1 ( I n + Λ 1 Γ 1 diag ( | x 1 | Γ 1 I n ) ) sign 2 I n Γ 2 ( x 2 ) ) ,
where M 1 and M 2 are positive constants.
Proof. 
Consider the following Lyapunov function
V = 1 2 s T s .
Then, the derivative of V becomes
V ˙ = s T s = s T ( x 2 + Λ 1 Γ 1 diag ( | x 1 | Γ 1 I n ) x 2 + Λ 2 Γ 2 diag ( | x 2 | Γ 2 I n ) ( F ( x 1 , x 2 ) + g ( x ) τ + d ( t , x ) ) ) .
Substituting (17) into (19), it follows that
V ˙ = s T ( x ˙ 1 + Λ 1 Γ 1 diag ( | x 1 | Γ 1 I n ) x ˙ 1 + Λ 2 Γ 2 diag ( | x 2 | Γ 2 I n ) ( d M 2 s ( ρ + M 1 ) s s Λ 2 1 Γ 2 1 I n + Λ 1 Γ 1 diag ( | x 1 | Γ 1 I n ) sign 2 I n Γ 2 ( x 2 ) ) ) = s T ( x ˙ 1 + Λ 1 Γ 1 diag ( | x 1 | Γ 1 I n ) x ˙ 1 + Λ 2 Γ 2 diag ( | x 2 | Γ 2 I n ) ( d M 2 s ( ρ + M 1 ) s s ) I n + Λ 1 Γ 1 diag ( | x 1 | Γ 1 I n ) x ˙ 1 ) = s T Λ 2 Γ 2 diag ( | x 2 | Γ 2 I n ) ( ρ + M 1 ) s s + M 2 s d .
Notice that
( s T ρ s s s T d ) s ρ + s · d s ρ + s · ρ = 0 .
Therefore, it is obviously that
V ˙ s T Λ 2 Γ 2 diag ( | x 2 | Γ 2 I n ) M 1 s + M 2 s min λ 2 i γ 2 i | x 2 i | γ 2 i 1 M 1 s + M 2 s T s .
Let
ρ 1 ( x 2 ) = min λ 2 i γ 2 i | x 2 i | γ 2 i 1 · 2 M 2 , ρ 2 ( x 2 ) = min λ 2 i γ 2 i | x 2 i | γ 2 i 1 · 2 M 1 .
Then, it follows that
V ˙ ρ 1 V ρ 2 V . ,
When x 2 i 0 , it is trivial that ρ 1 > 0 and ρ 2 > 0 . By Lemma 1, it can be confirmed that the states of the system will converge to the sliding surface s = 0 within finite time. Then, we consider the case that x 2 i = 0 . Substituting (17) into (3), we have
x ˙ 2 = ( C 0 q ˙ + G 0 ) M 0 ( M 2 s + ( ρ + M 1 ) s s M 0 1 ( C 0 x 2 + G 0 ) + Λ 2 1 Γ 2 1 ( I n + Λ 1 Γ 1 diag ( | x 1 | Γ 1 I n ) ) sign 2 I n Γ 2 ( x 2 ) ) + d ( t , x ) = M 2 s M 1 s s + ρ s s + d ( t , x ) M 2 s M 1 s s
which suggests that x ˙ 2 i < M 1 and x ˙ 2 i > M 1 for the cases s i > 0 and s i < 0 , respectively. This indicates that the states of the system will not remain at the points x 2 i = 0 . It is shown in Figure 1 that the controller will drive the states to the sliding surface s = 0 when x 2 = 0 . □
Remark 1.
Both NFTSM and conventional TSM share a nonlinear term of tracking error to achieve fast convergence when the states are far from the origin. However, the convergence rate changes slowly when the system approaches the sliding surface. Unlike TSM, a linear term of x 1 is introduced into the sliding surface to ensure a fast rate while approaching the nearby origin. See [15,16] for more details on how the states converge to zero when the sliding surface is reached.
Proposition 1.
If the following conditions are satisfied:
λ i = λ 2 i λ 1 i γ i γ 2 i , γ i = γ 2 i γ 1 i ,
then the states converge to the equilibrium point controlled using the NFTSMC faster than with the conventional nonsingular TSMC for any nonzero initial state x 1 ( 0 ) . λ i and γ i are parameters of the conventional nonsingular TSMC given by (6).
Proof. 
It is trivial that
( x 1 i + λ 1 i x 1 i γ 1 i ) 1 γ 2 i > λ 1 i 1 γ 2 i x 1 i γ 1 i γ 2 i ,
which follows that
T n f = 0 | x 1 i ( 0 ) | λ 2 i 1 γ 2 i ( x 1 i + λ 1 i x 1 i γ 1 i ) 1 γ 2 i d x 1 i < 0 | x 1 i ( 0 ) | ( λ 2 i / λ 1 i ) 1 / γ 2 i x 1 i γ 1 i / γ 2 i d x 1 i = T n ,
where T n f is the convergence time of the proposed controller. Therefore, the NFTSMC converges to the equilibrium faster than the conventional nonsingular TSMC. This completes the proof. □
Remark 2.
To reduce the chattering caused by switching, we use a continuous function in control law (17) to replace s / s ,
Θ = s s + δ ,
where δ is a positive small constant. As δ is selected to be small enough, Θ s / s .

3.2. High-Order CBF Design

Without loss of generality, consider the control affine system (10). Assume that the barrier function h ( x ) is a continuously differentiable function and satisfies
h ( r ) + k r 1 h ( r 1 ) + + k 1 h ˙ + k 0 h 0 ,
where k 0 , , k r 1 are real constants such that the roots of the following polynomial:
p ( λ ) = λ r + k r 1 λ r 1 + + k 1 λ + k 0
are all negative real, i.e., λ 1 , , λ r ( λ i > 0 , 1 i r ). We define
s 0 ( t , x ) = h ( x ) , s k = ( d d t + λ k ) s k 1 , 1 k r .
It is obvious from (26) that (24) is equivalent to s r ( t , x ) 0 . Denote s k ( 0 , x ( 0 ) ) as s k ( 0 ) , then we have the following lemma:
Lemma 2.
Given a continuously differentiable function h ( x ) and a set of real constant k 0 , , k r 1 such that (24) is satisfied and the roots of polynomial (25) are all negative real [34]. If s i ( 0 ) 0 , i = 0 , 1 , , r 1 , then h ( x ) 0 for any t 0 .
Proof. 
As s r ( t , x ) 0 and (26), we have
d d t e λ r t s r 1 ( t , x ( t ) ) 0 .
Integrate both sides on [ 0 , t ] and we have
s r 1 ( t , x ( t ) ) s r 1 ( 0 ) e λ r t .
Similarly, multiply both sides by e λ r 1 t and from (26), we have
d d t e λ r 1 t s r 2 ( t , x ( t ) ) s r 1 ( 0 ) e ( λ r 1 λ r ) t .
After integrating, we have
s r 2 ( t , x ( t ) ) s r 1 ( 0 ) e λ r 1 t 0 t e ( λ r 1 λ r ) τ 1 d τ 1 + s r 2 ( 0 ) e λ r 1 t .
Repeat the above iterative process and it becomes
s 0 ( t , x ( t ) ) s 0 ( 0 ) e λ 1 t + k = 1 r 1 s k ( 0 ) e λ 1 t 0 t e λ 1 λ 2 τ k 0 τ k e λ 2 λ 3 τ k 1 0 τ 2 e λ k λ k + 1 τ 1 d τ 1 d τ k 1 d τ k .
It is easy to check that the integral item of the above inequality is positive and, as t , it converges to zero. Since s i ( 0 ) 0 , it is trivial that s 0 ( t , x ( t ) ) 0 . Therefore, h ( x ) 0 for any t 0 . □
Then, we can have the following definition:
Definition 3.
Given control affine system (10), the smooth function h ( x ) : X R with relative degree r is defined as an exponential control barrier function if there exists K R r , such that x C
sup u U L f r h ( x ) + L g L f r 1 h ( x ) u + K T H 0 ,
where K = [ k 0 k 1 k r 1 ] T can be decided by the pole placement approach and
H = h ( x ) L f h ( x ) L f r 1 h ( x ) .
For all x X , define the set
K ecbf ( x ) = { u U | L f r h ( x ) + L g L f r 1 h ( x ) u + K T H 0 } .
According to Lemma 2, we can have the following proposition:
Proposition 2.
Given control affine system (10), h ( x ) is an ECBF with relative degree r. Then, any controller u K ecbf will render the safe set C defined by the (11) forward invariant.
Proposition 2 can be easily proved by Lemma 2. As defined above, the set K ecbf provides a family of applicable control inputs that render the system inside the safe set. Therefore, we only need to find a proper control input inside K ecbf to guarantee that the tracking performance is as accurate as possible.
We consider the state constraints as safety conditions to ensure that the Euler–Lagrange system can track the reference trajectory within state limits. Inspired by the work in [27], we use a super-ellipsoid safety region to ensure the safety of the system. We only consider the position constraints applied to the systems, and the safe set is defined as
C = { x R 2 n : h ( x ) 0 } ,
h ( x ) = 1 x 11 c 1 p 1 4 x 1 n c n p n 4 ,
where x i j is jth component of x i , c i is the center of the ellipse, and p i is the bound on the state. By taking the derivatives of h ( x ) along x , we have
h ( x ) x = 4 · x 11 c 1 p 1 3 x 1 n c n p n 3 0 1 × n , L f h ( x ) = h ( x ) x x 2 T F ( x 1 , x 2 ) T T = 4 x 11 c 1 p 1 3 x 21 4 x 12 c 1 p 1 3 x 22 4 x 1 n c 1 p 1 3 x 2 n ,
L g h ( x ) = h ( x ) x 0 n × n ( M 0 1 ) T T = 0 1 × n .
Apparently, function h ( x ) has relative degree two, which means the control input cannot be expressed explicitly, and the CBF-QP (14) of relative degree one is no longer applicable. Continue to take the derivative of h ( x ) and we have
L f h ( x ) = L f h ( x ) x , L f 2 h ( x ) = L f h ( x ) · x 2 T F ( x 1 , x 2 ) T T ,
L g L f h ( x ) = L f h ( x ) · 0 n × n ( M 0 1 ) T T .

3.3. ECBF-Based QP Design

We can formulate the ECBF condition into a quadratic program to refine the nominal control law (17) in real time. To minimize the control deviation from the nominal one, we take the control input error as
e u 2 = ( u τ ) T ( u τ ) = u T u 2 τ T u + τ T τ .
Then, we can reformulate the CBF-QP problem (14) as
u * ( x ) = arg min u R m u T u 2 τ T u s . t . L f 2 h ( x ) + L g L f h ( x ) u + K T H 0 ,
where K = [ k 0 k 1 ] T , H = [ h ( x ) L f h ( x ) ] T and L f h ( x ) , L f 2 h ( x ) , L g L f h ( x ) are defined by (30)–(33), respectively.
Remark 3.
Safety is always required to be satisfied even though the states of the system break the limit. When the constraints hold, the tracking error cannot converge to the origin, which means the system is unstable from this perspective. Therefore, this is a compromise between safety and performance.

4. Numerical Simulation

In this section, numerical simulations exemplify the validity and performance of the proposed control strategy. Consider a two-link robotic manipulator system [15], as shown in Figure 2, whose dynamics equation is described as an Euler–Lagrange system:
M ( q ¨ ) + C ( q , q ˙ ) q ˙ + G ( q ) = τ + τ d .
In this case, q i denotes the angular position of the ith arm, r i the length of the arm, m i the mass, and J i the moment of inertial. Denote q = [ q 1 q 2 ] T as the state of the system and the relative matrices are defined as
M ( q ¨ ) = M 11 M 12 M 21 M 22 ,
C ( q , q ˙ ) = β 12 q 2 q ˙ 1 2 2 β 12 q 2 q ˙ 1 q ˙ 2 β 12 q 2 q ˙ 2 2 ,
G ( q ) = m 1 + m 2 g r 1 cos q 2 + m 2 g r 2 cos q 1 + q 2 m 2 g r 2 cos q 1 + q 2 ,
where
M 11 = m 1 + m 2 r 1 2 + m 2 r 2 2 + 2 m 2 r 1 r 2 cos q 2 + J 1 , M 12 = m 2 r 2 2 + m 2 r 1 r 2 cos q 2 , M 21 = m 2 r 2 2 + m 2 r 1 r 2 cos q 2 , M 22 = m 2 r 2 2 + J 2 , β 12 ( q 2 ) = m 2 r 1 r 2 sin q 2
with parameters given as r 1 = 1   m , r 2 = 0.8   m , J 1 = 5   k g   m 2 , J 2 = 5   k g   m 2 , m 1 =   0.5   k g , m 2 = 1.5   k g , and g = 9.81   m / s 2 . As mentioned above, the system parameters sometimes cannot be acquired accurately and we assume that the mass is uncertain. Instead, we assume that the nominal values of m 1 and m 2 are m ^ 1 = 0.4   k g and m ^ 2 = 1.2   k g , respectively. The external disturbances are given by τ d = [ 2 sin ( t ) + 0.5 sin ( 200 π t ) , cos ( 2 t ) + 0.5 sin ( 200 π t ) ] T . Select the following initial conditions: q 10 = 1 , q 20 = 0 , and the desired trajectory is given by q d = [ 1.5 sin ( 0.8 t ) , 2.5 cos ( t ) ] T . The uncertainty term is assumed to satisfy Assumption 1 where the parameters are given by b 0 = 12 , b 1 = 2.2 , b 2 = 2.8 .
Define the tracking error ε 1 = q q d and ε 2 = q ˙ q ˙ d . Therefore, the tracking error model of system (3) can be expressed as
ε ˙ 1 = ε 2 , ε ˙ 2 = F e ( ε ) + M 0 1 ( q ) τ + d ( t , x ) ,
where F e ( ε ) = F ( x 1 , x 2 ) q ¨ d .

4.1. Controllers for Comparison

To demonstrate the effectiveness and robustness of the proposed NFTSM controller, we used the ATNTSM control scheme proposed in [16] with the same system parameters for comparison. In addition, the ATNTSM surface was chosen as
s = σ + C a ε 2 α / β ,
and the control input τ can be expressed as
τ = τ eq + τ sw ,
where
τ e q = M 0 ( q ) ( Π + F ( ε ) ) , τ s w = M 0 ( q ) a 0 + a 1 q + a 2 q ˙ 2 sign ( Λ s ) k M 0 ( q ) s ,
and k > 0 is an appropriately chosen constant. The remaining parameters of the controllers are given in Table 1.
The tracking performance of the given two-link robotic manipulator system is illustrated in Figure 3 and Figure 4. As shown in Figure 3, both controllers were able to track the reference trajectory, and the NFTSM controller performed better when the system was disturbed with a relatively large amplitude signal. To further illustrate the angular position tracking performance of the controllers, the position tracking error is shown in Figure 5. The NFTSM controller ensured that the position tracking error converged to zero, while the error caused by the ANTSM controller fluctuated in a small range. As can be seen in Figure 4, the ANTSM controller also fluctuated to a certain extent during speed tracking. Although the saturation function was utilized instead of the sign function, the ANTSM control input still changed relatively sharply. Figure 6 illustrates the control signal of each joint, and it shows that no singularity indeed occurred. Figure 7 shows that the sliding surfaces converged to zero within 0.5   s , much faster than the compared ANTSMC algorithm, and finite-time convergence was achieved. A comparison of the two controllers revealed that the NFTSM controller guarantees faster convergence and has stronger robustness.

4.2. States Constrained

Considering the states were constrained, we formulated the constraints on the angular position of each joint as safe constraints, we defined the safe set by (28), and h ( x ) was defined as
h ( x ) = 1 q 1 p 1 4 q 2 p 2 4 .
To illustrate the effectiveness of our proposed approach, we assumed the constraints on the angular position of each joint to be 1 q 1 1 and 2 q 2 2 , i.e., p 1 = 1 , p 2 = 2 . We still took q d = [ 1.5 sin ( 0.8 t ) , 2.5 cos ( t ) ] T as the reference trajectory, which sometimes broke the state constraints. As shown in Figure 8, when q i was about to break the limits, the controller responded quickly through the CBF-QP process to guarantee the states remained in the predefined safe set. When the system was operating safely, the zero tracking error with high robustness and finite-time convergence was guaranteed. The control inputs are shown in Figure 9.

4.3. Safety Behavior

We considered physical limitations on the robot manipulator system for more practical scenarios. We assumed that there was a wall above, and the joints should not hit the wall. The control purpose was to enforce the joints to the desired place without collisions. The initial position was q 0 = [ 0 , 0 ] T and desired position was q d = [ π , π ] T . The wall was about 1.5   m high, i.e.,
r 1 sin q 1 + r 2 sin q 2 1.5 .
Thus, we defined the control barrier function as
h ( x ) = 1.5 r 1 sin q 1 r 2 sin q 2 .
We applied the CBF-QP controller to this two-link robot manipulator with the aforementioned safe constraints. Figure 10 and Figure 11 show the trajectory of the angular position and velocity of each joint. It is shown that the robot reached the desired position in about 3   s . For a more precise illustration, the time evolution of the robot manipulator is provided in Figure 12. We can see that the joints were controlled correctly to avoid collision, while the joints controlled by the NFTSM-only controller hit the wall.

5. Conclusions

In this paper, a novel terminal sliding surface is introduced to improve the convergence rate of the conventional nonsingular terminal sliding mode control. We illustrate the effect of model uncertainties and external disturbances of the Euler–Lagrange system, and the high control performance and robustness of the proposed controller are guaranteed. Furthermore, the CBF-QP control strategy ensures that the system remains safe with the formulated state constraints. The simulation results illustrate the control performance and effectiveness of the proposed control strategy. This control strategy provides a general and practical approach for any mechanical system formulated as a Euler–Lagrange system.

Author Contributions

Methodology, G.S.; writing—original draft, G.S.; supervision, Q.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data is contained within the article. The data presented in this study are available in Section 4.

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
ATNTSMThree letter acronym
NFTSMCNonsingular fast terminal sliding mode control
QPQuadratic programming
CBFControl barrier function

References

  1. Cai, X.; de Queiroz, M. Adaptive rigidity-based formation control for multirobotic vehicles with dynamics. IEEE Trans. Control Syst. Technol. 2015, 23, 389–396. [Google Scholar] [CrossRef]
  2. Roy, S.; Roy, S.B.; Kar, I.N. Adaptive-Robust Control of Euler-Lagrange Systems with Linearly Parametrizable Uncertainty Bound. IEEE Trans. Control Syst. Technol. 2018, 26, 1842–1850. [Google Scholar] [CrossRef] [Green Version]
  3. Shao, X.; Sun, G.; Yao, W.; Liu, J.; Wu, L. Adaptive Sliding Mode Control for Quadrotor UAVs with Input Saturation. IEEE/ASME Trans. Mechatron. 2022, 27, 1498–1509. [Google Scholar]
  4. De Marinis, A.; Iavernaro, F.; Mazzia, F. A minimum-time obstacle-avoidance path planning algorithm for unmanned aerial vehicles. Numer. Algorithms 2022, 89, 1639–1661. [Google Scholar]
  5. Shakouri, A.; Shakouri, A. Prescribed-Time Control for Perturbed Euler-Lagrange Systems with Obstacle Avoidance. IEEE Trans. Autom. Control 2022, 67, 3754–3761. [Google Scholar] [CrossRef]
  6. Mobayen, S.; Alattas, K.A.; Assawinchaichote, W. Adaptive continuous barrier function terminal sliding mode control technique for disturbed robotic manipulator. IEEE Trans. Circuits Syst. I Regul. Pap. 2021, 68, 4403–4412. [Google Scholar]
  7. Chen, L.; Guo, Y.; Li, C.; Huang, J. Satellite Formationcontainment flying control with collision avoidance. AIAA J. Aerosp. Inf. Syst. 2018, 15, 253–270. [Google Scholar]
  8. Spong, M.W.; Hutchinson, S.; Vidyasagar, M. Robot Dynamics and Control; John Wiley & Sons: Hoboken, NJ, USA, 2008. [Google Scholar]
  9. Andrade, R.; Raffo, G.V.; Normey-Rico, J.E. Model predictive control of a tilt-rotor UAV for load transportation. In Proceedings of the 2016 European Control Conference (ECC), Aalborg, Denmark, 19 June–1 July 2016; pp. 2165–2170. [Google Scholar]
  10. Scaglioni, B.; Previtera, L.; Martin, J.; Norton, J.; Obstein, K.L.; Valdastri, P. Explicit Model Predictive Control of a Magnetic Flexible Endoscope. IEEE Robot. Autom. Lett. 2019, 4, 716–723. [Google Scholar] [PubMed] [Green Version]
  11. Song, Z.; Sun, K. Adaptive Sliding Mode Tracking Control for Uncertain Euler-Lagrange System. IEEE Access 2019, 7, 56817–56825. [Google Scholar]
  12. Li, H.; Li, X. Distributed Model Predictive Consensus of Heterogeneous Time-Varying Multi-Agent Systems: With and Without Self-Triggered Mechanism. IEEE Trans. Circuits Syst. I Regul. Pap. 2020, 67, 5358–5368. [Google Scholar]
  13. Mehrabian, A.; Khorasani, K. Distributed Formation Recovery Control of Heterogeneous Multiagent Euler–Lagrange Systems Subject to Network Switching and Diagnostic Imperfections. IEEE Trans. Control Syst. Technol. 2016, 24, 2158–2166. [Google Scholar] [CrossRef]
  14. Hou, H.; Yu, X.; Xu, L.; Rsetam, K.; Cao, Z. Finite-Time Continuous Terminal Sliding Mode Control of Servo Motor Systems. IEEE Trans. Ind. Electron. 2020, 67, 5647–5656. [Google Scholar]
  15. Feng, Y.; Yu, X.; Man, Z. Non-singular terminal sliding mode control of rigid manipulators. Automatica 2002, 38, 2159–2167. [Google Scholar]
  16. Zhai, J.; Xu, G. A novel non-singular terminal sliding mode trajectory tracking control for robotic manipulators. IEEE Trans. Circuits Syst. II Express Briefs 2021, 68, 391–395. [Google Scholar] [CrossRef]
  17. Labbadi, M.; Cherkaoui, M. Robust adaptive nonsingular fast terminal sliding-mode tracking control for an uncertain quadrotor UAV subjected to disturbances. ISA Trans. 2020, 99, 290–304. [Google Scholar] [PubMed]
  18. Moser, D.; Schmied, R.; Waschl, H.; del Re, L. Flexible Spacing Adaptive Cruise Control Using Stochastic Model Predictive Control. IEEE Trans. Control Syst. Technol. 2018, 26, 114–127. [Google Scholar]
  19. Hu, J.; Zhang, H.; Liu, L.; Zhu, X.; Zhao, C.; Pan, Q. Convergent Multiagent Formation Control With Collision Avoidance. IEEE Trans. Robot. 2020, 36, 1805–1818. [Google Scholar]
  20. Lee, J.H. Model predictive control: Review of the three decades of development. Int. J. Control Autom. Syst. 2011, 9, 415–424. [Google Scholar]
  21. Zhang, Z.; Wu, Y. Adaptive Fuzzy Tracking Control of Autonomous Underwater Vehicles with Output Constraints. IEEE Trans. Fuzzy Syst. 2021, 29, 1311–1319. [Google Scholar] [CrossRef]
  22. Zhang, J.; Li, K.; Li, Y. Output-Feedback Based Simplified Optimized Backstepping Control for Strict-Feedback Systems with Input and State Constraints. IEEE/CAA J. Autom. Sin. 2021, 8, 1119–1132. [Google Scholar]
  23. Jiang, X.; Zhao, D. Event-triggered fault detection for nonlinear discrete-time switched stochastic systems: A convex function method. Sci. China Inf. Sci. 2021, 64, 200204. [Google Scholar] [CrossRef]
  24. Zhang, T.; Deng, F.; Sun, Y.; Shi, P. Fault estimation and fault-tolerant control for linear discrete time-varying stochastic systems. Sci. China Inf. Sci. 2021, 60, 200201. [Google Scholar]
  25. Ames, A.D.; Grizzle, J.W.; Tabuada, P. Control Barrier Function based Quadratic Programs with Application to Adaptive Cruise Control. In Proceedings of the 53rd IEEE Conference on Decision and Control, Los Angeles, CA, USA, 15–17 December 2014; pp. 6271–6278. [Google Scholar]
  26. Wang, L.; Ames, A.D.; Egerstedt, M. Safety Barrier Certificates for Collisions-Free Multirobot Systems. IEEE Trans. Robot. 2017, 33, 661–674. [Google Scholar]
  27. Wang, L.; Ames, A.D.; Egerstedt, M. Safe certificate-based maneuvers for teams of quadrotors using differential flatness. In Proceedings of the 2017 IEEE International Conference on Robotics and Automation (ICRA), Singapore, 29 May–3 June 2017; pp. 3293–3298. [Google Scholar]
  28. Grandia, R.; Taylor, A.J.; Ames, A.D.; Hutter, M. Multi-layered safety for legged robots via control barrier functions and model predictive control. In Proceedings of the 2021 IEEE International Conference on Robotics and Automation (ICRA), Xi’an, China, 30 May–5 June 2021; pp. 8352–8358. [Google Scholar]
  29. Hsu, S.C.; Xu, X.; Ames, A.D. Control barrier function based quadratic programs with application to bipedal robotic walking. In Proceedings of the 2015 American Control Conference (ACC), Chicago, IL, USA, 1–3 July 2015; pp. 4542–4548. [Google Scholar]
  30. Yu, S.; Yu, X.; Shirinzadeh, B.; Man, Z. Continuous finite-time control for robotic manipulators with terminal sliding mode. Automatica 2005, 41, 1957–1964. [Google Scholar] [CrossRef]
  31. Ames, A.D.; Coogan, S.; Egerstedt, M.; Notomista, G.; Sreenath, K.; Tabuada, P. Control barrier functions: Theory and applications. In Proceedings of the 2019 18th European Control Conference (ECC), Napoli, Italy, 25–28 June 2019; pp. 3420–3431. [Google Scholar]
  32. Ames, A.D.; Xu, X.; Grizzle, J.W.; Tabuada, P. Control Barrier Function Based Quadratic Programs for Safety Critical Systems. IEEE Trans. Autom. Control 2017, 62, 3861–3876. [Google Scholar]
  33. Nguyen, Q.; Sreenath, K. Exponential control barrier functions for enforcing high relative-degree safety-critical constraints. In Proceedings of the 2016 American Control Conference (ACC), Boston, MA, USA, 6–8 July 2016; pp. 322–328. [Google Scholar]
  34. Xu, X. Constrained control of input-output linearizable systems using control sharing barrier functions. Automatica 2018, 87, 195–201. [Google Scholar]
Figure 1. The phase trajectory of the system.
Figure 1. The phase trajectory of the system.
Mathematics 11 02788 g001
Figure 2. 2-DOF robot manipulator.
Figure 2. 2-DOF robot manipulator.
Mathematics 11 02788 g002
Figure 3. Angular position of joint–1 (upper) and joint–2 (lower) tracking responses to a sin-type desired trajectory with unknown external disturbances.
Figure 3. Angular position of joint–1 (upper) and joint–2 (lower) tracking responses to a sin-type desired trajectory with unknown external disturbances.
Mathematics 11 02788 g003
Figure 4. Angular velocity of joint–1 (upper) and joint–2 (lower) tracking responses.
Figure 4. Angular velocity of joint–1 (upper) and joint–2 (lower) tracking responses.
Mathematics 11 02788 g004
Figure 5. Angular position tracking error of joint–1 (upper) and joint–2 (lower).
Figure 5. Angular position tracking error of joint–1 (upper) and joint–2 (lower).
Mathematics 11 02788 g005
Figure 6. Control responses to the desired trajectory of joint–1 (upper) and joint–2 (lower).
Figure 6. Control responses to the desired trajectory of joint–1 (upper) and joint–2 (lower).
Mathematics 11 02788 g006
Figure 7. The time evolution of sliding surface of each joint.
Figure 7. The time evolution of sliding surface of each joint.
Mathematics 11 02788 g007
Figure 8. Trajectory tracking of the angular position with safe constraints.
Figure 8. Trajectory tracking of the angular position with safe constraints.
Mathematics 11 02788 g008
Figure 9. Trajectory tracking of the angular position with safe constraints.
Figure 9. Trajectory tracking of the angular position with safe constraints.
Mathematics 11 02788 g009
Figure 10. Angular position trajectory of the joints.
Figure 10. Angular position trajectory of the joints.
Mathematics 11 02788 g010
Figure 11. Angular velocity trajectory of the joints.
Figure 11. Angular velocity trajectory of the joints.
Mathematics 11 02788 g011
Figure 12. Two-link robot manipulator with physical constraints. (a) NFTSMC only; (b) the CBF-QP was utilized based on the nominal NFTSM controller. The dashed line is state constraint. The blue and orange line are joint–1 and joint–2, respectively.
Figure 12. Two-link robot manipulator with physical constraints. (a) NFTSMC only; (b) the CBF-QP was utilized based on the nominal NFTSM controller. The dashed line is state constraint. The blue and orange line are joint–1 and joint–2, respectively.
Mathematics 11 02788 g012
Table 1. Controller parameters.
Table 1. Controller parameters.
VariablesValuesVariablesValues
M 1 1 k 0 600
M 2 2 k 1 50
Γ 1 diag ( [ 2 2 ] ) p 1 1
Γ 2 diag ( [ 5 / 3 5 / 3 ] ) p 2 2
Λ 1 diag ( [ 2 2 ] ) α 5
Λ 2 diag ( [ 1 1 ] ) β 3
C a diag ( [ 0.3 0.4 ] ) k2.5
δ 1 , δ 2 3.5 × 10 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.

Share and Cite

MDPI and ACS Style

Sun, G.; Zeng, Q. A Finite-Time Sliding Mode Control Approach for Constrained Euler–Lagrange System. Mathematics 2023, 11, 2788. https://doi.org/10.3390/math11122788

AMA Style

Sun G, Zeng Q. A Finite-Time Sliding Mode Control Approach for Constrained Euler–Lagrange System. Mathematics. 2023; 11(12):2788. https://doi.org/10.3390/math11122788

Chicago/Turabian Style

Sun, Guhao, and Qingshuang Zeng. 2023. "A Finite-Time Sliding Mode Control Approach for Constrained Euler–Lagrange System" Mathematics 11, no. 12: 2788. https://doi.org/10.3390/math11122788

APA Style

Sun, G., & Zeng, Q. (2023). A Finite-Time Sliding Mode Control Approach for Constrained Euler–Lagrange System. Mathematics, 11(12), 2788. https://doi.org/10.3390/math11122788

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