Next Article in Journal
Fracture Behavior of Permeable Asphalt Mixtures with Steel Slag under Low Temperature Based on Acoustic Emission Technique
Previous Article in Journal
Dynamic Displacement Estimation for Long-Span Bridges Using Acceleration and Heuristically Enhanced Displacement Measurements of Real-Time Kinematic Global Navigation System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discrete Sliding Mode Control Design for Bilateral Teleoperation System via Adaptive Extended State Observer

1
Beijing Advanced Innovation Center for Biomedical Engineering, Beihang University, Beijing 100191, China
2
School of Biological Science and Medical Engineering, Beihang University, Beijing 100191, China
3
The Institute of Electrical Engineering, Yanshan University, Qinhuangdao 066004, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(18), 5091; https://doi.org/10.3390/s20185091
Submission received: 12 July 2020 / Revised: 23 August 2020 / Accepted: 2 September 2020 / Published: 7 September 2020
(This article belongs to the Section Sensors and Robotics)

Abstract

:
The goal of this paper is to improve the synchronization control performance of nonlinear teleoperation systems with system uncertainties in the presence of time delays. In view of the nonlinear discrete states of the teleoperation system in packet-switched communication networks, a new discrete sliding mode control (DSMC) strategy is performed via a new reaching law in task space. The new reaching law is designed to reduce the chattering and improve control performance. Moreover, an adaptive extended state observer (AESO) is used to estimate the total system disturbances. The additional gain of AESO is adjusted in time to decrease the estimation errors of both system states and disturbances automatically and improve the estimation performances of the AESO. Finally, the validity of the designed control strategy is demonstrated by both simulation and experiments. Furthermore, the experimental comparison results indicate that the improvement is achievable with the proposed AESO and DSMC.

1. Introduction

Nonlinear bilateral teleoperation systems can expand the range of perception, and enable humans to complete complex tasks in a remote operating environment. A representative nonlinear bilateral teleoperation system consists of the following five parts: human operator, master robot, the network communication channel, slave robot, and remote environment. In recent years, the potential applications of teleoperation systems are in the various fields, i.e., the remote handling of hazardous materials [1], underwater maintenance and repairing tasks [2], space exploration [3], telemedicine [4], and so on. In all of these applications, the tasks occur in long-distance and/or hazardous environments. However, due to the limitation of bandwidth for the communication, there will be inevitable time delays. As we all know, time delay is one of the factors influencing the stability of nonlinear bilateral teleoperation systems. In practice, because of the extremely complicated nonlinearity of the nonlinear bilateral teleoperation systems and the certain constraint conditions on their outputs or states, nonlinear bilateral teleoperation systems havedifficultly in performing ideal performances.
Nowadays, a series of effective control strategies have been performed to solve the stability problems of nonlinear bilateral teleoperation systems. In [5], a notion of wave variable had been developed to handle the time delay issue. In [6], the instability, caused by time delays, had been conquered with a new passivity and scattering theory. In [7], a robust adaptive control algorithm is developed to deal with the system’s uncertainties and to provide a smooth estimation of delayed reference signals. In [8], Baranitha tackled the time-varying delay problem for a single-master multi-slave teleoperation system by assuming an asymmetric and semi-Markovian jump protocol for communication of the slaves with the master. There were also the passivity-based assumptions employed to ensure the stability of the position and velocity of the teleoperation systems, such as [9, 10]. After that, Lu proposes a relative impedance-based force control method for a bimanual robot teleoperation system with varying time delays in [11]. Lu adopts a strategies that design a hybrid error item to minimize both the position and force tracking errors. Additionally, in order to reduce the tracking error and ensure the stability of the system, an asymmetrical wave variable compensation method is proposed [12], where the forward wave variable compensates the backward wave variable. However, the abovementioned research methods only work for joint space-based teleoperation systems. When the slave and the master are kinematically different; for instance, the slave robot is bigger/smaller than another one. Previous controllers in the joint space cannot bring about a satisfactory working performance. In order to make the robot end effector reach the desired position in the task space, Takegaki and Arimoto proposed a new position control method-based Jacobian task space transition matrix [13]. With the development of this study, many methodologies have been widely described in the literature to aim at task space-based teleoperation system control such as the new nonlinear saturated proportional derivative (SPD) control strategy based on approximate Jacobian matrix [14], a novel (nP+D)-like controller for task-space tracking performance [15], nonlinear adaptive saturation control strategy with limited amplitude [16], adaptive control schemes based on assumed parametric linearization of kinematics and dynamics [17], and so on. It is worth noting that there is no full discussion about the system uncertainties and external disturbances in the above literature, although the system stability and synchronization performances were achieved. Indeed, the above-mentioned consequences are performed under certain assumptions where the models for the communication network and the master/slave controllers have been supposed to be in the continuous-time domain. In reality, the communication media are composed of undependable communication networks that may have variable delays, lost or reordered packets. Therefore, there is an urgent need for nonlinear bilateral teleoperation in developing a discrete-time theory.
It is well known that sliding mode control (SMC) has effective control performance for both linear and nonlinear systems, and provides more noticeable robustness and simpler calculation than other robust control approaches [18]. However, the computation of the computer is based on discrete sample systems. If the continuous-time SMC algorithms are applied directly to discrete-time systems, it is of concern that there would be some indomitable problems such as chattering and discretization errors. Therefore, the design of discrete sliding mode control (DSMC) has attracted a lot of attention in recent years. There are also some contributions in the area of DSMC such as [19]. Better performance could be guaranteed by considering a sampling period in the design phase, even if the sampling period was quite long [20]. In [21], Ma developed a novel dead zone sliding mode reaching law with disturbance compensation for uncertain discrete-time systems. A new adaptive second order DSMC scheme is proposed for a class of uncertain nonlinear systems [22]; however, better robustness and trace performance come at the cost of slightly more complex control logic. In [23], an optimal sliding surface coupled with a delay predictor had been presented to construct a DSMC for overcoming system parametric uncertainty. Moreover, we concentrate particularly on the active disturbance rejection control (ADRC). In the ADRC strategy, the extended state observer (ESO) is adopted to obtain the real-time estimation of the lumped system uncertainties including both internal unmodeled dynamics and external disturbance in the system [24]. In [25], a third-order super-twisting extended state observer was proposed, which enhanced the estimation convergence and provided strong disturbance estimation against fast speed and load variation. Liu [26] proposed that an improved model predictive control (MPC) combined with extended state observe exhibits better control performance and faster dynamic response, where the ESO provides real-time disturbance compensation for the prediction control algorithm. Consequently, in this paper, we propose a reasonable method to effectively handle total system disturbances is by introducing ESOs for the uncertain teleoperation system.
Since the communication between the master and the slave is achieved through a packet-switched network, this paper focuses specifically on the synchronization control for discrete-time nonlinear bilateral teleoperation system in task space, which includes internal unmodeled dynamics, external disturbance and time-delays. Therefore, this paper is aimed at developing a new discrete sliding mode control algorithm to ensure the synchronization performance between the master and the slave via applying the adaptive extended state observer (AESO) to compensate total disturbances. The main accomplishments of the proposed strategy are summarized as follows: (i) with the aim of avoiding discretization after control design, a discrete sliding mode control algorithm is designed, in order to consider the discrete states caused by unreliable communication networks; (ii) a new reaching law of DSMC is developed to reduce the chattering while ensuring the tracking error quickly converges to zero domain; (iii) The parameter θ of AESO is designed so that the estimation errors quickly converge to smaller neighborhood and experimental comparisons demonstrate better. Finally, the proposed control method is simulated and tested by an experiment executed on a nonlinear bilateral teleoperation system composed of two Phantom Premium 1.5A robot manipulators. The test results reveal that the proposed control algorithm yields remarkable control performance.
The rest of this paper is organized as following. The concerned background is discussed in Section 2. In Section 3, we propose AESO-based DSMC and present the stability analyses of the nonlinear bilateral teleoperation system based on the designed controllers. Section 4 shows the simulation and experiment results. The concluding remarks are given in Section 5.

2. Problem Statement and Preliminaries

In this section, a n 1 -degree-of-freedom (DOF) master-slave nonlinear bilateral telerobotic system is considered as follows:
{ M q m ( q m ) q ¨ m + C q m ( q m , q ˙ m ) q ˙ m + g q m ( q m ) + f q m ( q ˙ m ) + B q m ( q m ) = τ m + J m T ( q m ) F h M q s ( q s ) q ¨ s + C q s ( q s , q ˙ s ) q ˙ s + g q s ( q s ) + f q s ( q ˙ s ) + B q s ( q s ) = τ s J s T ( q s ) F e
where m and s denote the master side and the slave side of the nonlinear bilateral teleoperation system, respectively, q i R n 1 with i = m , s is the joint position vector, q ˙ i R n 1 is the joint velocity vector, q ¨ i R n 1 is the joint accelerated velocity vector, M q i ( q i ) R n 1 × n 1 is the positive-definite inertia matrix, C q i ( q i , q ˙ i ) R n 1 × n 1 is the matrix of centripetal and Coriolis term, g q i ( q i ) R n 1 is the gravitational vector, f q i ( q ˙ i ) R n 1 is the viscous friction vector, B q i ( q i ) R n 1 denotes the unknown bounded external disturbance, τ i R n 1 is control torque and τ i = J i J ( q i ) u i , u i R n is applied input control vector, J i ( q i ) R n × n 1 is Jacobian matrix, and F h , F e R n denote master operator force and external environmental force, respectively.
We review the properties [27,28,29] for teleoperation system as follows:
Property 1.
The inertia matrix M ( q ) is a symmetric positive-definite matrix, and there are two positive constants m 1 and m 2 such that m 1 M ( q ) m 2 .
Property 2.
There exists a positive scalar b such that C ( q , x ) y b x y , with regard to all q , x , y R n .
Property 3.
If q ¨ and q ˙ are bounded, the time derivative of the term C ( q , q ˙ ) is also bounded.
Assumption 1.
In this article, the Jacobian matrix J i ( q i ) , i = m , s is supposed to be known and that the robot arms are working in a finite task space; in other words, the Jacobian matrix has full rank.
Let χ m , χ s R n represent the task coordinates of the task-space end effectors. The positional relationship between task space and joint space can be represented by the following relation
χ i = h i ( q i ) , χ ˙ i = J i ( q i ) q ˙ i
where i = m , s , h i ( ) : R n R n represents the mapping relationship from joint space to task space, and J i = h i ( q i ) / q i .
Then, q ˙ i , q ˙ i are expressed as follows:
q ˙ i = J i 1 ( q i ) χ ˙ i q i = J ˙ i 1 ( q i ) χ ˙ i + J i 1 ( q i ) χ ¨ i
After the model transformation, in the task space, the system (1) can be rewritten as:
M m χ ¨ m + C m χ ˙ m + g m + f m + B m = u m + F h M s χ ¨ s + C s χ ˙ s + g s + f s + B s = u s F e
where i = m , s ,
M i = ( J i T ( q i ) ) 1 M q i ( q i ) J i 1 ( q i ) C i = ( J i T ( q i ) ) 1 ( M q i ( q i ) J ˙ i 1 ( q i ) + C q i ( q i , q ˙ i ) J i 1 ( q i ) ) g i = ( J i T ( q i ) ) 1 g q i ( q i ) f i = ( J i T ( q i ) ) 1 f q i ( q i ) B i = ( J i T ( q i ) ) 1 B q i ( q i )
In most practical applications, the precise mode cannot be obtained directly, due to the noise, friction, viscous friction, uncertain disturbances and so on. Due to the existence of certain uncertainties in M i and C i , the system (4) is rewritten as follows:
M om χ ¨ m + C om χ ˙ m M om Θ m = u m M os χ ¨ s + C os χ ˙ s M os Θ s = u s
where M i = M o i + Δ M i , C i = C o i + Δ C i , Θ i = M o i 1 ( Δ M i χ ¨ i + Δ C i χ ˙ i + g i + f i + B i F i ) , i = m , s . M o i , C o i represent the nominal parts, while Δ M i and Δ C i represent the uncertainties. When i = m , F i = F h , otherwise, F i = F e . In this paper, Θ i denotes the lumped system uncertainty, and it is assumed to be bounded.
Introducing the state vector X i = [ X i 1 T X i 2 T ] T , i = m , s , let X m 1 = χ m , X m 2 = χ ˙ m , X s 1 = χ s , X s 2 = χ ˙ s , then the system (5) is transformed as follows:
X ˙ i 1 = X i 2 X ˙ i 2 = f i ( X i ) + Θ i + H i u i
Furthermore, the system (6) is rewritten as follows:
X ˙ i = A X i + B f i ( X i ) + B Θ i + B H i u i
where i = m , s , A = [ 0 I 0 0 ] , B = [ 0 I ] , H i = M o i 1 .
Notice that all the terms proposed in (6) and (7) can be easily calculated from (5). With the sampling time h , the discretization of the uncertain model equation is given as:
X i ( j h + h ) = A ¯ i X i ( j h ) + B ¯ i f i ( j h ) + B ¯ i Θ i ( j h ) + B ¯ i H i u i ( j h )
where
A ¯ i = exp ( [ A i ] h ) = [ I h I 0 I ]
B ¯ i = 0 h exp ( [ A i ] h ) d γ [ B i ] = [ h 2 2 ! I h I ]

3. Main Results

This section addresses a presentation of a proposed control scheme based on discrete-time nonlinear bilateral teleoperation system, where the new AESO is introduced to estimate and compensate the uncertainty. Meanwhile, the development of a discrete-time SMC algorithm is depicted, with a view to settle the synchronization problem of bilateral teleoperation system in task space.

3.1. Adaptive Extended State Observer

In this part, the adaptive extended state observer will be employed for the system (9). Firstly, for the convenience of further analysis, the following state measurement values are introduced:
Y i ( j h ) = X i 1 ( j h ) + n i j
where i = m , s X i 1 ( j h ) is the output to be controlled, and Y i ( j h ) R n is the one to be measured, which includes the measurement noise vector n i j R n . Then, the following assumptions for the observer are given as [30]:
Assumption 2.
{ n i j } 1 is a white random sequence and
E ( n i j n i j T ) R ¯ i
where R ¯ i is a known matrix.
Assumption 3.
E [ X i ( 0 ) X ^ i ( 0 ) Θ i ( 0 ) Θ ^ i ( 0 ) ] [ X i ( 0 ) X ^ i ( 0 ) Θ i ( 0 ) Θ ^ i ( 0 ) ] T P i 0
where [ X i ( 0 ) X ^ i ( 0 ) Θ i ( 0 ) Θ ^ i ( 0 ) ] is the estimation error of the AESO, and P i 0 is a known matrix.
Assumption 4.
E ( Θ i ( X i ( t ) , t ) Θ i ( X i ( j h ) , j h ) ) ( Θ i ( X i ( t ) , t ) Θ i ( X i ( j h ) , j h ) ) T Q ¯ i
where t [ j h , j h + h ) , j 0 , Θ i ( X i ( t ) , t ) Θ i ( X i ( j h ) , j h ) is the discretized error and Q ¯ i is a known diagonal matrix.
Therefore, the linear structure of AESO is designed as follows:
[ X ^ i ( j h + h ) Θ ^ i ( j h + h ) ] = A ˜ i [ X ^ i ( j h ) Θ ^ i ( j h ) ] + B ˜ i H i u i ( j h ) + B ˜ i f i ( j h ) + L i j ( Y i ( j h ) X ^ i 1 ( j h ) )
where A ˜ i = exp ( [ A i B i 0 0 ] h ) = [ I h I h 2 2 ! I 0 I h I 0 0 I ] , B ˜ i = 0 h exp ( [ A i B i 0 0 ] h ) d γ [ B i 0 ] = [ h 2 2 ! I h I 0 ] , and L i j is the gain of the discrete AESO, ensuring that X ^ i ( j h ) and Θ ^ i ( j h ) can be employed as the estimation of X i ( j h ) and Θ i ( j h ) , respectively. Then, the original values of AESO (13) are taken as:
X ^ i 1 ( 0 ) = Y i ( 0 ) , X ^ i k ( 0 ) = 0 , k 2 , Θ ^ i ( 0 ) = 0
Therefore, the gain L i j of the AESO is designed as follows:
L i j = A ˜ i ( I + θ i ) P i j C i ( C i T ( I + θ i ) P i j C i + R ¯ i ) 1
P i ( j + 1 ) = ( A ˜ i L i j C i T ) ( I + θ i ) P i j ( A ˜ i L i j C i T ) T + L i j R ¯ i L i j T + ( I + θ i 1 ) Q ¯ i
where
Q i = 3 n [ h 4 Q ¯ i 0 0 0 h 2 Q ¯ i 0 0 0 Q ¯ i ] 3 n × 3 n ,   C i = [ I 0 0 ] 3 ( n + 1 ) × n
And e i g ( θ i ) [ θ i 1 θ i 2 θ i ( 3 n ) ] T , θ i k ¯ = Q i ( k ¯ , k ¯ ) / P i 0 ( k ¯ , k ¯ ) , k ¯ = 1 , 2 , , 3 n
Define the estimation error ξ i ( j h ) as follows:
ξ i ( j h ) = [ X i ( j h ) X ^ i ( j h ) Θ i ( X i ( j h ) , j h ) Θ ^ i ( j h ) ]
then, we will provide the following Theorem 1 showing the property of the estimation error ξ i ( j h ) .
Theorem 1.
If Assumptions 2–4 hold, and there exist the positive real numbers a , c , p 1 , p 2 , ε 1 , ε 2 , θ ¯ , q and r as following:
A ˜ i = a θ i = θ ¯ C i = c Q i = q p 1 I P i j p 2 I R ¯ i = r ε 1 ξ i ( j h ) ε 2
so that the inequality is fulfilled:
ο = 1 p 1 ( a + ( 1 + θ ¯ ) a p 2 c 2 / r ) ( q + p 2 ) + q p 1 + ( 1 + θ ¯ ) 2 a 2 p 2 2 c 2 p 1 r < a ε 1 2 p 2
Then the estimation error { ξ i ( j h ) } j = 0 will be uniformly bounded in the mean square, if the initial estimation error ξ i ( 0 ) satisfies:
ξ i ( 0 ) ε 2
Before proof of Theorem 1, the following lemma is first discussed.
Lemma 1
([31]).In the case of Theorem 1, there is a real number 0 < α < 1 such that Π i j = P i j 1 meet the inequality:
( A ˜ i L i j C i T ) T Π i ( j + 1 ) ( A ˜ i L i j C i T ) ( 1 α ) Π i j
For j 0 with L i j given by (15).
Proof. 
From (16), we have:
P i ( j + 1 ) = ( A ˜ i L i j C i T ) ( I + θ i ) P i j ( A ˜ i L i j C i T ) T + L i j R ¯ i L i j T + ( I + θ i 1 ) Q ¯ i = ( A ˜ i L i j C i T ) [ ( I + θ i ) P i j + ( A ˜ i L i j C i T ) 1 ( L i j R ¯ i L i j T + ( I + θ i 1 ) Q ¯ i ) ( A ˜ i L i j C i T ) T ] ( A ˜ i L i j C i T ) T
From (19), and C i T ( I + θ i ) P i j C i > 0 , we have:
L i j ( 1 + θ ¯ ) a p 2 c 1 r
since, we get:
P i ( j + 1 ) ( A ˜ i L i j C i T ) [ ( I + θ i ) P i j + ( 1 + θ ¯ ) 2 a 2 p 2 2 c 2 / r + ( 1 + θ ¯ 1 ) q ( a + ( 1 + θ ¯ ) a p 2 c 2 / r ) 2 I ] ( A ˜ i L i j C i T ) T
Because of P i j p 1 I , invert both sides of this inequality (25), multiply both sides by ( A ˜ i L i j C i T ) T and ( A ˜ i L i j C i T ) , then we finally obtain:
( A ˜ i L i j C i T ) T Π i ( j + 1 ) ( A ˜ i L i j C i T ) [ 1 + θ ¯ + ( 1 + θ ¯ ) 2 a 2 p 2 2 c 2 / r + ( 1 + θ ¯ 1 ) q ( a + ( 1 + θ ¯ ) a p 2 c 2 / r ) 2 ] 1 Π i j
Therefore:
1 α = [ 1 + θ ¯ + ( 1 + θ ¯ ) 2 a 2 p 2 2 c 2 / r + ( 1 + θ ¯ 1 ) q ( a + ( 1 + θ ¯ ) a p 2 c 2 / r ) 2 ] 1
 □
Proof of Theorem 1.
From (7) and (13), we have:
ξ i ( j h + h ) = ( A ˜ i L i j C i T ) ξ i ( j h ) + L i j n i j + [ W i j G i j ]
where:
W i j = j h ( j + 1 ) h [ ( ( j + 1 ) h γ ) n 1 n 1 I ( ( j + 1 ) h γ ) 1 I I ] ( Θ i ( X i ( t ) , t ) Θ i ( X i ( j h ) , j h ) ) d γ
G i j = Θ i ( j h + h ) Θ i ( j h )
and the last term of Equation (28) satisfies:
E [ W i j G i j ] [ W i j G i j ] T Q i
Define Γ i j = E ( ξ i ( j h ) ξ i ( j h ) T ) , thus:
Γ i ( j + 1 ) = ( A ˜ i L i j C i T ) Γ i j ( A ˜ i L i j C i T ) T + E [ W i j G i j ] [ W i j G i j ] T + ( A ˜ i L i j C i T ) E ξ i ( j h ) ( L i j n i j + [ W i j G i j ] ) T + E ( L i j n i j + [ W i j G i j ] ) ξ i T ( j h ) ( A ˜ i L i j C i T ) T + L i j E n i j n i j T L i j T
With the above knowable, the measurement noise vector n i j and the estimation error ξ i ( j h ) are unrelated, so the following inequality is given as:
Γ i ( j + 1 ) ( A ˜ i L i j C i T ) Γ i j ( A ˜ i L i j C i T ) T + E [ W i j G i j ] [ W i j G i j ] T + ( A ˜ i L i j C i T ) E ( ξ i ( j h ) [ W i j G i j ] ) T + E ( [ W i j G i j ] ξ i T ( j h ) ) ( A ˜ i L i j C i T ) T + L i j R ¯ i L i j T
According to Assumption 4, the last third and second terms of (31) have the following upper bound:
( A ˜ i L i j C i T ) E ( ξ i ( j h ) [ W i j G i j ] ) T + E ( [ W i j G i j ] ξ i T ( j h ) ) ( A ˜ i L i j C i T ) T ( A ˜ i L i j C i T ) θ i Γ i j ( A ˜ i L i j C i T ) T + θ i 1 Q i
Note that the proper θ i can make the two sides of the inequality as close as possible. If and only if j = 1 , the equation of (32) can be achieved. In consequence, we get the following equation:
( A ˜ i L i 0 C i T ) θ i E ξ i ( 0 ) ξ i T ( 0 ) ( A ˜ i L i 0 C i T ) T = θ i 1 E [ W i 0 G i 0 ] [ W i 0 G i 0 ] T
where:
E ξ i ( 0 ) ξ i T ( 0 ) = P i 0 ,   E [ W i 0 G i 0 ] [ W i 0 G i 0 ] T = Q i
Because of A ˜ i L i 0 C i T I , the Equation (33) and Equation (34) indicate that:
θ i 2 = Q i P i 0 1
Since Q i and P i 0 are diagonal matrices, we set:
e i g ( θ i ) [ θ i 1 θ i 2 θ i ( 3 n ) ] T ,   θ i k ¯ = Q i ( k ¯ , k ¯ ) / P i 0 ( k ¯ , k ¯ ) ,   k ¯ = 1 , 2 , , 3 n
 □
Remark 1.
For the parameter θ i , if it is chosen as a constant, each item of gain L i j is iterated in the same way. Therefore, it is possible that the corresponding gain for the derivative of total disturbance estimation is large, which could easily lead to an overgrowth of the total disturbance estimation. Consequently, the gain of controller would be set sufficiently large to stabilize the closed-loop teleoperation system. Unfortunately, this will come with the actuator saturation problem. In order to solve this problem, the new parameter θ i in the form of a diagonal matrix, is proposed in this paper. In other words, parameters are adjusted, respectively, according to the variation of different system states in this paper. Therefore, each state could reach its ideal one as quickly as possible, then the stability will be guaranteed.
Next:
Γ i ( j + 1 ) ( A ˜ i L i j C i T ) ( I + θ i ) Γ i j ( A ˜ i L i j C i T ) T + L i j R ¯ i L i j T + ( I + θ i 1 ) Q i
Thus:
P i ( j + 1 ) = ( A ˜ i L i j C i T ) ( I + θ i ) Γ i j ( A ˜ i L i j C i T ) T + L i j R ¯ i L i j T + ( I + θ i 1 ) Q i
and P i 0 Γ i 0 , we get:
Γ i ( j + 1 ) P i ( j + 1 ) ,   j 0
Therefore, P i ( j + 1 ) is regarded as the upper bound of the covariance matrix of ξ i ( j h + h ) , and it is also recommended to be minimized by proposing L i j . Design
J i j = t r a c e ( P i ( j + 1 ) )
When the partial derivative of J i j with respect to L i j is zero, that is:
J i j L i j = 0
It is easy to know
L i j = A ˜ i ( I + θ i ) P i j C i ( C i T ( I + θ i ) P i j C i + R ¯ i ) 1
We design Lyapunov function as follows:
V i j = ξ i T ( j h ) Π i j ξ i ( j h )
From (28), we have:
V i ( j + 1 ) = [ ( A ˜ i L i j C i T ) ξ i ( j h ) + L i j n i j + [ W i j G i j ] ] T Π i ( j + 1 ) [ ( A ˜ i L i j C i T ) ξ i ( j h ) + L i j n i j + [ W i j G i j ] ]
and using Lemma 1 and Equation (43), the Equation (44) can be transformed as:
V i ( j + 1 ) ( 1 α ) V i j + [ W i j G i j ] T Π i ( j + 1 ) [ ( A ˜ i L i j C i T ) ξ i ( j h ) + L i j n i j ] + ξ i T ( j h ) ( A ˜ i L i j C i T ) T Π i ( j + 1 ) ( L i j n i j + [ W i j G i j ] ) + n i j T L i j T Π i ( j + 1 ) ( ( A ˜ i L i j C i T ) ξ i ( j h ) + [ W i j G i j ] ) + [ W i j G i j ] T Π i ( j + 1 ) [ W i j G i j ] + n i j T L i j T Π i ( j + 1 ) L i j n i j
Since ξ i ( j h ) is not related to n i j , we can get further:
V i ( j + 1 ) ( 1 α ) V i j + [ W i j G i j ] T Π i ( j + 1 ) ( A ˜ i L i j C i T ) ξ i ( j h ) + ξ i T ( j h ) ( A ˜ i L i j C i T ) T Π i ( j + 1 ) [ W i j G i j ] + [ W i j G i j ] T Π i ( j + 1 ) [ W i j G i j ] + n i j T L i j T Π i ( j + 1 ) L i j n i j
For the second and third terms of the (46), apply the (19) and we will have:
[ W i j G i j ] T Π i ( j + 1 ) ( A ˜ i L i j C i T ) ξ i ( j h ) + ξ i T ( j h ) ( A ˜ i L i j C i T ) T Π i ( j + 1 ) [ W i j G i j ] 1 p 1 ( a + ( 1 + θ ¯ ) a p 2 c 2 / r ) ( [ W i j G i j ] T ξ i ( j h ) + ξ i T ( j h ) [ W i j G i j ] ) 1 p 1 ( a + ( 1 + θ ¯ ) a p 2 c 2 / r ) ( [ W i j G i j ] T [ W i j G i j ] + ξ i T ( j h ) ξ i ( j h ) ) 1 p 1 ( a + ( 1 + θ ¯ ) a p 2 c 2 / r ) ( q + p 2 )
For the fourth term of the (46), apply the (19) and we will have:
[ W i j G i j ] T Π i ( j + 1 ) [ W i j G i j ] 1 p 1 [ W i j G i j ] T [ W i j G i j ] q p 1
and for the fifth term of the (46), apply the (19) and we will have
n i j T L i j T Π i ( j + 1 ) L i j n i j 1 p 1 n i j T L i j T L i j n i j
Since both sides of (49) are scalars, we can trace the right side of (49) and at the same time will not change its value, that is
n i j T L i j T Π i ( j + 1 ) L i j n i j 1 p 1 t r ( n i j T L i j T L i j n i j ) 1 p 1 t r ( L i j n i j n i j T L i j T ) ( 1 + θ ¯ ) 2 a 2 p 2 2 c p 1 r 2
Finally, synthesize the inequalities (46)–(50) above and we can get:
V i ( j + 1 ) V i j α V i j + ο
where ο = 1 p 1 ( a + ( 1 + θ ¯ ) a p 2 c 2 / r ) ( q + p 2 ) + q p 1 + ( 1 + θ ¯ ) 2 a 2 p 2 2 c p 1 r 2 . When ο < a ε 1 2 p 2 , using (19) and (43), we obtain:
V i ( j + 1 ) V i j 0
From discussion above, we may conclude that the appropriate positive parameters a , c , p 1 , p 2 , ε 1 , ε 2 , θ ¯ , q and r may be utilized to guarantee the stability of the error system (143), and the estimation error ξ i ( j h ) converges to a smaller zero domain.

3.2. Discrete Sliding Mode Surface

In the above subsection, the lumped system uncertainty of the bilateral teleoperation system is estimated by the new AESO (13). In this part, its estimation is employed as the compensation for the uncertainty, and DSMC is provided by the following main theorem. For the nonlinear bilateral teleoperation system, tracking trajectories of master and slave robots are X dm ( j h ) = X s ( j h T s ) and X ds ( j h ) = X m ( j h T m ) , respectively, where T i , i = m , s is the constant time delay. T m represents the time delay of signal transmission from the master robot to the slave robot and T s stands for the time delay of signal transmission from the slave robot to the master robot. So, the position synchronization errors between the master and the slave are defined as follows:
e m ( j h ) = X m ( j h ) X dm ( j h ) e s ( j h ) = X s ( j h ) X ds ( j h )
Next, a switching surface is constructed, and the dynamic equation of the ideal sliding mode is obtained by using the equivalent control method. Then, sufficient conditions for the asymptotic dynamic stability of the sliding mode are obtained by using Lyapunov stability theory and LMI (Linear Matrix Inequality) method.
Select the following switching function:
s i ( j h ) = Λ i e i ( j h )
with Λ i R n × 2 n , i = m , s , which is assumed that Λ i B ¯ i 0 . Then
s i ( j h + h ) = Λ i e i ( j h + h ) = Λ i A ¯ i X i ( j h ) + Λ i B ¯ i f i ( j h ) + Λ i B ¯ i Θ i ( j h ) + Λ i B ¯ i H i u i ( j h ) Λ i X d i ( j h + h )
According to the discrete sliding mode theory, the equation s i ( j h + h ) = s i ( j h ) = 0 is obtained when the system states reach the switching surface. Therefore, the ideal sliding mode equivalent control law can be obtained from formulas (8) and (55) as follows:
u i eq ( j h ) = ( Λ i B ¯ i H i ) 1 ( Λ i A ¯ i X i ( j h ) + Λ i B ¯ i f i ( j h ) + Λ i B ¯ i Θ i ( j h ) Λ i X d i ( j h + h ) )
Note 1.
The equivalent control law shown in Equation (57) is only used to analyze the stability of the sliding mode dynamic equation, and it is only used as a tool here. The actual sliding mode control law will be designed below.
By substituting the equivalent control law (56) into the system (8), the sliding mode dynamic equation of the uncertain discrete time-delay system under ideal conditions can be obtained as follows:
X i ( j h + h ) = ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) X i ( j h ) + ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) f i ( j h ) + ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) Θ i ( j h ) + B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i X d i ( j h + h )
The stability of the sliding mode dynamic equation (58) is further analyzed by using the LMI technique, and sufficient conditions are given to guarantee the asymptotic stability of the sliding mode dynamic equation.
Theorem 2.
For uncertain discrete time-delay system (8), switch function (54) is selected. If there exist positive definite matrix P i > 0 and scalar ϖ i 1 , ϖ i 2 , and the following LMI inequalities (58)-(60) are satisfied, then the sliding mode dynamic system (57) is asymptotically stable:
[ Ξ i 2 ϒ i T P i 8 ( ϒ i A ¯ i ) T P i 2 P i ϒ i P i 0 8 P i ϒ i A ¯ i 0 P i ] < 0
[ P P B ¯ i E E T B ¯ i T P ϖ i 1 I ] < 0
[ P i P i ϒ i B ¯ i E E T ( ϒ i B ¯ i ) T P i ϖ i 2 I ] < 0
where i = m , s , Ξ i = 8 A ¯ i T P i A ¯ i P i + 8 ϖ i 1 ( D f i T D f i + D Θ i T D Θ i ) + 8 ϖ i 2 ( D f i T D f i + D Θ i T D Θ i ) , ϒ i = B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i .
Assumption 5.
f i and Θ i represent the nonlinear part of the teleoperation system and the disturbance, so the following condition is satisfied:
[ f i ( j h ) Θ i ( j h ) ] = E F ( j h ) [ D f i D Θ i ]
where E , D f i , D Θ i and F ( j h ) are known constant matrices, and F ( j h ) satisfies F T ( j h ) F ( j h ) I .
Lemma 2
([32]). If E and D are the real matrixes with the appropriate dimensions, and F ( j h ) satisfies F T ( j h ) F ( j h ) I . Then, for any non-zero constant ϖ > 0 , the following inequality exists:
E F ( j h ) D + D T F T ( j h ) E T = ϖ 1 E T E + ϖ D T D
Lemma 3
([32]). If ν and υ are the real matrixes with the appropriate dimensions, for any non-zero matrix U > 0 , then
ν υ + ν T υ T = ν U ν T + υ T U 1 υ
Lemma 4
([33]).  (Schur’s theorem) For the following LMI inequality
[ Ω 11 Ω 12 Ω 12 Ω 22 ] < 0
where Ω 22 > 0 , Ω 11 + Ω 12 Ω 22 1 Ω 12 < 0 , so this is the same thing as Ω 11 = Ω 11 T , Ω 22 = Ω 22 T .
Proof. 
In this proof, to facilitate the proof process, the hypothesizes and lemmas above are given. □
Let us take the Lyapunov function as
V i ( j h ) = X i T ( j h ) P i X i ( j h )
Along the state trajectory of the system (57), it can be obtained that:
Δ V i ( j h ) = V i ( j h + h ) V i ( j h ) = Δ V i 1 ( j h ) + Δ V i 2 ( j h )
where
Δ V i 1 ( j h ) = X i T ( j h ) ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) T P i ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) X i ( j h ) + f i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) f i ( j h ) + Θ i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) Θ i ( j h ) + X d i T ( j h + h ) ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) T P i ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) X d i ( j h + h ) X i T ( j h ) P i X i ( j h )
Δ V i 2 ( j h ) = X i T ( j h ) ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) T P i ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) f i ( j h ) + f i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) X i ( j h ) + X i T ( j h ) ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) T P i ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) Θ i ( j h ) + Θ i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) X i ( j h ) + X i T ( j h ) ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) T P i ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) X d i ( j h ) + X d i T ( j h ) ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) T P i ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) X i ( j h ) + f i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) Θ i ( j h ) + Θ i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) f i ( j h ) + f i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) X d i ( j h ) + X d i T ( j h ) ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) T P i ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) f i ( j h ) + Θ i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) X d i ( j h ) + X d i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) Θ i ( j h )
For the first and second terms of the expression Δ V i 2 ( j h ) , it can be seen from the Lemma 3 that:
X i T ( j h ) ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) T P i ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) f i ( j h ) + f i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) X i ( j h ) X i T ( j h ) ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) T P ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) X i ( j h ) + f i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) f i ( j h )
Similarly, for other terms of the expression Δ V i 2 ( j h ) , we can get:
Δ V i 2 ( j h ) 3 X i T ( j h ) ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) T P i ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) X i ( j h ) + 3 f i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) f i ( j h ) + 3 Θ i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) Θ i ( j h ) + 3 X d i T ( j h ) ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) T P i ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) X d i ( j h )
Therefore,
Δ V i ( j h ) X i T ( j h ) P i X i ( j h ) + 4 X i T ( j h ) ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) T P i ( A ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) X i ( j h ) + 4 f i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) f i ( j h ) + 4 Θ i T ( j h ) ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) Θ i ( j h ) + 4 X d i T ( j h ) ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) T P i ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) X d i ( j h )
According to the Lemma 3, the following inequality is further obtained:
Δ V i ( j h ) X i T ( j h ) ( 8 A ¯ i T P i A ¯ i P i + ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) T P i ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i A ¯ i ) ) X i ( j h ) + f i T ( j h ) ( 8 B ¯ i T P i B ¯ i + ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) ) f i ( j h ) + Θ i T ( j h ) ( 8 B ¯ i T P i B ¯ i + ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) T P i ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i B ¯ i ) ) Θ i ( j h ) + 4 X d i T ( j h ) ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) T P i ( B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i ) X d i ( j h )
In addition, it can be obtained from Assumption 5 that
Δ V i ( j h ) X i T ( j h ) ( 8 A ¯ i T P i A ¯ i P i + ( ϒ i A ¯ i ) T P i ( ϒ i A ¯ i ) ) X i ( j h ) + X i T ( j h ) D f i T F T ( j h ) E T ( 8 B ¯ i T P i B ¯ i + ( ϒ i B ¯ i ) T P i ( ϒ i B ¯ i ) ) E F ( j h ) D f i X i ( j h ) + X i T ( j h ) D Θ i T F T ( j h ) E T ( 8 B ¯ i T P i B ¯ i + ( ϒ i B ¯ i ) T P i ( ϒ i B ¯ i ) ) E F ( j h ) D Θ i X i ( j h ) + 4 X d i T ( j h ) ( ϒ i ) T P i ( ϒ i ) X d i ( j h )
where i = m , s , ϒ i = B ¯ i H i ( Λ i B ¯ i H i ) 1 Λ i . So, here we have
Δ V i ( j h ) [ X i T ( j h ) X d i T ( j h ) ( ϒ i ) T ] T Σ i [ X i ( j h ) X d i ( j h ) ]
where:
Σ i = [ Ξ i 1 0 0 Ξ i 2 ]
Ξ i 1 = 8 A ¯ i T P i A ¯ i P i + ( ϒ i A ¯ i ) T P i ( ϒ i A ¯ i ) + D f i T F T ( j h ) E T ( 8 B ¯ i T P i B ¯ i + ( ϒ i B ¯ i ) T P i ( ϒ i B ¯ i ) ) E F ( j h ) D f i + D Θ i T F T ( j h ) E T ( 8 B ¯ i T P i B ¯ i + ( ϒ i B ¯ i ) T P i ( ϒ i B ¯ i ) ) E F ( j h ) D Θ i
Ξ i 2 = 4 ( ϒ i ) T P i ( ϒ i )
Therefore, when Σ i < 0 , Δ V i ( j h ) < 0 is true (If X i ( j h ) 0 ). According to the Lemma 4, the inequality Σ i < 0 can be equivalent to:
[ 8 A ¯ i T P i A ¯ i P i 2 ϒ i T P i 8 ( ϒ i A ¯ i ) T P i 8 D f i T F T ( j h ) E T B ¯ i T P i 2 ϒ i P i P i 0 0 8 P i ϒ i A ¯ i 0 P i 0 8 P i B ¯ i E F ( j h ) D f i 0 0 P i 0 8 P i B ¯ i E F ( j h ) D Θ i 0 0 8 P i ϒ i B ¯ i E F ( j h ) D f i 0 0 0 0 8 P i ϒ i B ¯ i E F ( j h ) D Θ i 0 0 0 8 D f i T F T ( j h ) E T B ¯ i T ϒ i T P i 0 8 D Θ i T F T ( j h ) E T B ¯ i T P i 0 8 D Θ i T F T ( j h ) E T B ¯ i T ϒ i T P 0 0 0 0 0 0 P i 0 0 0 P i 0 0 0 P i ] < 0
Furthermore, the inequality (73) is equivalent to
[ 8 A ¯ i T P i A ¯ i P i 2 ϒ i T P i 8 ( ϒ i A ¯ i ) T P i 0 2 ϒ i P i P i 0 0 8 P i ϒ i A ¯ i 0 P i 0 0 0 0 P i 0 0 0 0 8 P i ϒ i B ¯ i E F ( j h ) D f i 0 0 0 0 8 P i ϒ i B ¯ i E F ( j h ) D Θ i 0 0 0 8 D f i T F T ( j h ) E T B ¯ i T ϒ i T P i 0 0 0 8 D Θ i T F T ( j h ) E T B ¯ i T ϒ i T P 0 0 0 0 0 0 P i 0 0 0 P i 0 0 0 P i ] + W i 1 T [ F T ( j h ) 0 0 F T ( j h ) ] Y i 1 T + Y i 1 [ F ( j h ) 0 0 F ( j h ) ] W i 1 < 0
where
W i 1 T = [ 8 D f i T 0 0 0 0 0 0 0 8 D Θ i T 0 0 0 0 0 ] Y i 1 T = [ 0 0 0 E T B ¯ i T P i 0 0 0 0 0 0 0 E T B ¯ i T P i 0 0 ] T
According to the Lemma 2, we can get:
W i 1 T F T ( j h ) Y i 1 T + Y i 1 F ( j h ) W i 1 < ϖ i 1 1 Y i 1 T Y i 1 + ϖ i 1 W i 1 W i 1 T
Then, the following inequalities (76) and (77) can ensure that the (75) is true
[ Ξ i 3 2 ϒ i T P i 8 ( ϒ i A ¯ i ) T P i 2 ϒ i P i P i 0 8 P i ϒ i A ¯ i 0 P i 8 P i ϒ i B ¯ i E F ( j h ) D f i 0 0 0 8 P i ϒ i B ¯ i E F ( j h ) D Θ i 0 8 D f i T F T ( j h ) E T B ¯ i T ϒ i T P i 0 0 8 D Θ i T F T ( j h ) E T B ¯ i T ϒ i T P 0 0 P i 0 0 P i ] < 0
P i + ϖ i 1 1 P i B ¯ i E E T B ¯ i T P i < 0
where Ξ i 3 = 8 A ¯ i T P i A ¯ i P i + 8 ϖ i 1 ( D f i T D f i + D Θ i T D Θ i ) .
Similarly, the inequality (76) is equivalent to
W i 2 T [ F T ( j h ) 0 0 F T ( j h ) ] Y i 2 T + Y i 2 [ F ( j h ) 0 0 F ( j h ) ] W i 2 + [ Ξ i 3 2 ϒ i T P i 8 ( ϒ i A ¯ i ) T P i 2 ϒ i P i P i 0 8 P i ϒ i A ¯ i 0 P i 0 0 0 0 0 0 0 0 0 0 0 0 P i 0 0 P i ] < 0
where
W i 2 T = [ 8 D f i T 0 0 0 0 0 8 D Θ i T 0 0 0 ]
Y i 2 T = [ 0 0 0 E T B ¯ i T ϒ i T P i 0 0 0 0 0 E T B ¯ i T ϒ i T P i ] T
Then the following inequalities (79) and (80) can guarantee the existence of the (78):
[ Ξ i 2 ϒ i T P i 8 ( ϒ i A ¯ i ) T P i 2 ϒ i P i P i 0 8 P i ϒ i A ¯ i 0 P i ] < 0
P i + ϖ i 2 1 P i ϒ i B ¯ i E E T B ¯ i T ϒ i T P i < 0
where Ξ i = 8 A ¯ i T P i A ¯ i P i + 8 ϖ i 1 ( D f i T D f i + D Θ i T D Θ i ) + 8 ϖ i 2 ( D f i T D f i + D Θ i T D Θ i ) .
Therefore, according to the Lemma 4, if LMI inequalities (58)–(60) are true, then Σ i < 0 . Therefore, Δ V i ( j h ) < 0 . According to Lyapunov stability theory, the sliding mode dynamic system (57) is asymptotically stable. The proof process is over.

3.3. Discrete Sliding Mode Controller

This is an example of an equation: The control objective of this paper is that the system trajectory tracking error from any initial state can arrive at a switching surface s i ( j h ) = 0 and reach the origin along the sliding surface. So, define a switching belt encircling sliding surface as follows:
s i Δ ( j h ) = { s i ( j h ) | Δ i k s i k ( j h ) Δ i k , k = 1 , 2 , , n }
where i = m , s , k denotes the kth coordinate in task space, s i ( j h ) = [ s i 1 ( j h ) s i n ( j h ) ] T , Δ i = [ Δ i 1 Δ i n ] T . The 2 Δ i is the width of the switching belt. In order to reduce chattering and improve dynamic quality in the sliding stage, a new reaching law, which consists of the index term Φ i ( j h ) and the convergence parameter p i > 0 ( 1 p i h < 1 ) is designed by
s i ( j h + h ) = ( 1 p i h ) Φ i ( j h ) s i ( j h ) λ i Φ i ( j h ) tanh ( s i ( j h ) σ i )
where λ i , i = m , s , is a switching gain, σ i > 0 , and
Φ i ( j h ) = δ i + ( 1 δ i ) e φ i s i ( j h ) γ i { 0 < δ i < 1 ,   φ i > 0 γ i > 0 ,   a n d   γ i N
Remark 2.
In practical applications, the tanh function is advocated to replace the sign function in order to reduce chattering, and a mass of experimental data also verifies its effectiveness. However, no one has theoretically analyzed the feasibility for the tanh function. This paper will make an attempt in this area. Theoretical analysis will be presented in Theorem 3. According to expression (83), a conclusion can be drawn that if the system tracking error is far away from the discrete sliding surface, Φ i ( j h ) tends to δ i . Thus ( 1 p i h ) Φ i ( j h ) tends to ( 1 p i h ) δ i , which is less than ( 1 p i h ) , and λ i Φ i ( j h ) tends to λ i δ i , which is greater than λ i . On the contrary, when the system tracking error is close to the discrete sliding surface, Φ i ( j h ) tends to 1. Thus ( 1 p i h ) Φ i ( j h ) tends to ( 1 p i h ) , and λ i Φ i ( j h ) tends to λ i . Obviously, the index term Φ i ( j h ) is positive, therefore the stability of the nonlinear bilateral teleoperation system will be unaffected.
With consideration of the system (8) and the new reaching law (82), the proposed DSMC controller for the nonlinear bilateral teleoperation system can be given as
u i ( j h ) = ( Λ i B ¯ i H i ) 1 ( Λ i A ¯ i X i ( j h ) + Λ i B ¯ i f i ( j h ) Λ i X d i ( j h + h ) ( 1 p i h ) Φ i ( j h ) s i ( j h ) + λ i Φ i ( j h ) tanh ( s i ( j h ) σ i ) ) M o i Θ ^ i ( j h )
The control block diagram of nonlinear bilateral teleoperation system with controller (86) is displayed in Figure 1.
For the design of the discrete sliding mode controller, the following lemma and assumption will be utilized.
Lemma 5
([34]). For x R , x x tanh ( x δ ) 0.2785 δ , where δ > 0 R .
Assumption 6.
According to the Theory 1, the disturbanceestimation error Θ i ( j h ) Θ ^ i ( j h ) is bounded and converges to a very small range. Therefore, it is reasonable to assume that:
ς i k ( j h ) ς i k
where i = m , s , k = 1 , 2 , , n , ς i ( j h ) = Λ i B ¯ i ( Θ i ( j h ) Θ ^ i ( j h ) ) = [ ς i 1 ( j h ) ς i n ( j h ) ] T , ς i = [ ς i 1 ς i n ] T , and ς i k is the positive constant and denotes the upper bound of ς i k ( j h ) .
Theorem 3.
For the discrete nonlinear bilateral teleoperation system (89) based on the controller (84), under the Assumption 6, if the following condition is maintained:
ς i k λ i
(a). The system trajectory tracking error from any initial state will enter this switching belt s i Δ ( j h ) of DSMC defined by
s i Δ ( j h ) = { s i ( j h ) | s i k ( j h ) Δ i k = 0.2785 δ i λ i λ i ς i k Φ i ( j h ) , k = 1 , 2 , , n }
(b). Once the system trajectory tracking errors enter this switching belt s i Δ ( j h ) , they cannot escape from it.
In this paper, because the coupling relationship between the states of nonlinear bilateral teleoperation system is compensated as the uncertainty during the AESO design phase, with regard to k , the controllers are independent on each other and then the stability analysis can be demonstrated in the same way. Hereinafter, the stability analysis is discussed for only one.
Proof. 
In this proof, two cases will be considered, namely, the positive and negative values of s i k ( j h ) . The equivalent form of Theorem 3 is represented as follows:
0.2785 δ i λ i λ i ς i k Φ i ( j h ) < s i k ( j h ) < 0.2785 δ i λ i λ i ς i k Φ i ( j h )
where i = m , s , k = 1 , 2 , , n .
(a). Case 1: If s i k ( j h ) > 0 , due to 0 < Φ i ( j h ) < 1 , it can be obtained from (82) and (86) that:
s i k ( j h + h ) = ( 1 p i h ) Φ i ( j h ) s i k ( j h ) λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k ( j h ) s i k ( j h ) λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k
According to the Lemma 5, then
s i k ( j h + h ) s i k ( j h ) λ i Φ i ( j h ) ( 1 0.2785 δ i s i k ( j h ) ) + ς i k s i k ( j h )
Therefore, the sequence { s i k ( j h ) } is strictly monotonously decreasing when s i k ( j h ) > 0 . It is concluded that there must exist a positive integer j = j so that the following inequality (91) holds:
s i k ( j h + h ) = ( 1 p i h ) Φ i ( j h ) s i k ( j h ) λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k ( j h ) 0.2785 δ i λ i λ i ς i k Φ i ( j h )
when j j , the system trajectory tracking error enters the switching belt s i Δ ( j h ) of DSMC.
Case 2: If s i k ( j h ) < 0 , due to 0 < Φ i ( j h ) < 1 and the condition (86), it can be obtained that
s i k ( j h + h ) = ( 1 p i h ) Φ i ( j h ) s i k ( j h ) λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k ( j h ) s i k ( j h ) λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k s i k ( j h ) λ i Φ i ( j h ) ( 1 0.2785 δ i s i k ( j h ) ) + ς i k s i k ( j h )
Therefore, the sequence { s i k ( j h ) } is strictly monotonously increasing when s i k ( j h ) < 0 . It is concluded that there must exist a positive integer j = j so that the following inequality (93) holds:
s i k ( j h + h ) = ( 1 p i h ) Φ i ( j h ) s i k ( j h ) λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k ( j h ) > 0.2785 δ i λ i λ i ς i k Φ i ( j h )
when j j , the system trajectory tracking error enters the switching belt s i Δ ( j h ) of DSMC.
Therefore, in view of (91) and (93), it can be concluded that if s i k ( j h ) lies outside the switching belt s i Δ ( j h ) defined by (87), then the system trajectory tracking error from any initial state will enter into this switching belt of DSMC.
(b). Case 1: When the s i k ( j h ) enters the switching belt s i Δ ( j h ) , namely 0 < s i k ( j h ) < 0.2785 δ i λ i λ i ς i k Φ i ( j h ) , then:
s i k ( j h + h ) = ( 1 p i h ) Φ i ( j h ) s i k ( j h ) λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k ( j h ) λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k ( j h ) 0.2785 δ i λ i λ i ς i k Φ i ( j h )
Suppose s i k ( j h ) has escaped the switching belt s i Δ ( j h ) again, namely, s i k ( j h + h ) 0.2785 δ i λ i λ i ς i k Φ i ( j h ) s i k ( j h ) , then
s i k ( j h + h ) = ( 1 p i h ) Φ i ( j h ) s i k ( j h ) λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k ( j h ) < s i k ( j h ) λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k ( j h ) < 0.2785 δ i λ i λ i ς i k Φ i ( j h )
However, the result is inconsistent with the hypothesis s i k ( j h + h ) s i k ( j h ) , so that the hypothesis is not set up. Hence,
s i k ( j h + h ) < 0.2785 δ i λ i λ i ς i k Φ i ( j h )
Case 2: When the s i k ( j h ) enters the switching belt s i Δ ( j h ) defined by (87), namely, 0.2785 δ i λ i λ i ς i k Φ i ( j h ) < s i k ( j h ) < 0 , then
s i k ( j h + h ) = ( 1 p i h ) Φ i ( j h ) s i k ( j h ) λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k ( j h ) < λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k ( j h ) < 0.2785 δ i λ i λ i ς i k Φ i ( j h )
Suppose s i k ( j h ) has escaped the switching belt s i Δ ( j h ) again, namely, s i k ( j h + h ) 0.2785 δ i λ i λ i ς i k Φ i ( j h ) s i k ( j h ) , then
s i k ( j h + h ) = ( 1 p i h ) Φ i ( j h ) s i k ( j h ) λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k ( j h ) > s i k ( j h ) λ i Φ i ( j h ) tanh ( s i k ( j h ) σ i ) + ς i k ( j h ) > 0.2785 δ i λ i λ i ς i k Φ i ( j h )
However, the result is inconsistent with the hypothesis s i k ( j h + h ) s i k ( j h ) so that the hypothesis is not set up. Hence,
s i k ( j h + h ) 0.2785 δ i λ i λ i ς i k Φ i ( j h )
In short, s i k ( j h ) [ 0.2785 δ i λ i λ i ς i k Φ i ( j h ) , 0.2785 δ i λ i λ i ς i k Φ i ( j h ) ] , when s i k ( j h + h ) [ 0.2785 δ i λ i λ i ς i k Φ i ( j h ) , 0.2785 δ i λ i λ i ς i k Φ i ( j h ) ] . That is to say, once the system trajectory tracking error enters this switching belt s i Δ ( j h ) , they cannot escape from it. □

4. Simulation and Experiment

In this section, the simulated analyses and experimental results on the strength of the proposed control strategy are presented, with the purpose of verifying the effectiveness of the discrete SMC algorithm and the superiority of its control performance. Firstly, the simulations are implemented on two three-degree-of-freedom robot arms
{ M q m ( q m ) q ¨ m + C q m ( q m , q ˙ m ) q ˙ m + g q m ( q m ) + f q m ( q ˙ m ) + B q m ( q m ) = τ m + J m T ( q m ) F h M q s ( q s ) q ¨ s + C q s ( q s , q ˙ s ) q ˙ s + g q s ( q s ) + f q s ( q ˙ s ) + B q s ( q s ) = τ s J s T ( q s ) F e
The kinematic relationship between task space and joint space of the nonlinear bilateral teleoperation system is given by
χ i 1 = l i 1 cos ( q i 1 ) + l i 2 cos ( q i 1 + q i 2 ) + l i 3 cos ( q i 1 + q i 2 + q i 3 ) χ i 2 = l i 1 sin ( q i 1 ) + l i 2 sin ( q i 1 + q i 2 ) + l i 3 sin ( q i 1 + q i 2 + q i 3 ) χ i 3 = q i 1 + q i 2 + q i 3
where i = m , s , l i 1 , l i 2 , l i 3 represent the link lengths and q i 1 , q i 2 , q i 3 represent the joint angles of the robot arms.
For the simulation, the relevant parameter values are set as l i 1 = 2.05 , l i 2 = 2.05 , l i 3 = 0.5 , i = m , s . The human force F h is imposed on the master robot, which is depicted in Figure 2. While in the slave site, the force is kept to zero. The simulation results are utilized to prove the following: (1) When the master robot moves, does the slave robot follow the master? (2) When the human-input force vanishes, does the tracking error between the master and slave vanish as well? The aim of (1) and (2) is to present the stability of the nonlinear bilateral teleoperation system.
Now, we evaluate the validity of the proposed control strategy for the nonlinear bilateral teleoperation system (100) with existence of both time delays and flexible friction forces. First, time delays are adjusted to T m = T s = 300   ms . The sampling period is h = 0.001   s , and parameters of adaptive ESO are chosen as
e i g ( R ¯ m ) = e i g ( R ¯ s ) = [ 0.002 2 0.002 2 0.002 2 ] T
e i g ( Q ¯ m ) = e i g ( Q ¯ s ) = [ 0.1 2 0.1 2 0.1 2 ] T
e i g ( Q ¯ m ) = e i g ( Q ¯ s ) = [ 0.005 2 0.005 2 0.005 2 0.05 2 0.05 2 0.05 2 0.04 2 0.04 2 0.04 2 ] T
Parameters of discrete-time sliding mode surface are set as λ m = λ s = 0.15 , δ m = δ s = 0.25 , φ m = φ s = 20 , γ m = γ s = 10 , Λ m = Λ s = [ 0.3 0 0 0.05 0 0 0 0.2 0 0 0.02 0 0 0 0.2 0 0 0.02 ] , σ m = σ s = 1 , p m = p s = 350 .
Afterwards, the trajectory estimation errors of master and slave robots are exhibited in Figure 3 and Figure 4, respectively.
The disturbance estimation errors of the master and the slave are displayed in Figure 5 and Figure 6. As shown, estimation errors of both the motion trajectory and the disturbance converge to a very small range, and then approach rapidly to zero after F h = 0 . Therefore, accurate estimation of total system disturbances is provided by AESO. In order to show the superiority of the proposed controller, PD controller (Proportional derivative controller) is taken as the comparison term. Now, the motion trajectories of the master and the slave end effectors with proposed DSMC are shown in Figure 7. Obviously, the slave can accurately reproduce the trajectory of the master and the prominent synchronization performance is completed. On the contrary, in Figure 8, at the beginning of the movement, the chattering is quite serious, and the position tracking effect of the master robot and the slave robot is poor.
To further prove the superiority of the controller designed in this paper, the control algorithm is implemented on the teleoperation experiment platform built in the laboratory. The nonlinear bilateral teleoperation system model includes a couple of Phantom Premium 1.5A robot arms (SensAble Technologies, Inc.) to be performed, which is shown in Figure 9. In the experiment, the sampling period is h = 0.002   s , and the parameters of the discrete sliding mode control algorithm are selected as λ m = λ s = 0.4 , δ m = δ s = 0.25 , φ m = φ s = 20 , γ m = γ s = 10 , Λ m = Λ s = [ 0.3 0 0 0.08 0 0 0 0.3 0 0 0.05 0 0 0 0.3 0 0 0.05 ] , σ m = σ s = 1 , p m = p s = 350 .
Then the parameters of adaptive ESO are chosen as
e i g ( R ¯ m ) = e i g ( R ¯ s ) = [ 0.002 2 0.002 2 0.002 2 ] T
e i g ( Q ¯ m ) = e i g ( Q ¯ s ) = [ 0.1 2 0.1 2 0.1 2 ] T
e i g ( Q ¯ m ) = e i g ( Q ¯ s ) = [ 6 × 10 18 6 × 10 18 6 × 10 18 0.0245 2 0.0245 2 0.0245 2 12.247 2 12.247 2 12.247 2 ] T
Before the experiment, the information transmission delay between the master and the slave robots by Simulink module is set to T m = T s = 200   m s . With the aim of verifying the performance superiority of the AESO designed in this paper, the contrast experiment between the proposed AESO and another compared one with a constant parameter θ i is provided. The experiment results are displayed in Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16. As shown in Figure 10, Figure 11, Figure 12 and Figure 13, the motion trajectory estimation errors and the total disturbance estimation of the master and the slave based on this paper and the compared one is shown. It is noticed that the trajectory estimation errors based on this paper converge within 0.02 m, while others based on the compared one converge within 0.3 m. Comparing the total disturbance estimation, the AESO designed in this paper could effectively estimate the total disturbance; whereas, another one based on the compared one brings a hundredfold increase in total disturbance estimation and thus loses the estimation performance.
In addition, Figure 14 and Figure 15 display the control torques of the master and the slave robots, respectively. The motion trajectory tracking between the master and the slave is shown in Figure 16. We can see that the master robot stops moving after 50 s, and then the slave robot also promptly stops moving. Ultimately, the two robots stop in the same position.
Next, the communication time delay is increased to T m = T s = 500 ms. The trajectory estimations and the estimation errors for the master and the slave based on this paper are displayed in Figure 17, Figure 18, Figure 19 and Figure 20. The control torques of the master and the slave robots are shown in Figure 21 and Figure 22, respectively. Furthermore, Figure 23 presents the motion trajectory tracking between the master and the slave robots. From this experiment, although the slave robot responds slowly due to the increased time delay, it still follows the master motion. The experimental results and analyses further validate the superiority and validity of the proposed control strategy in this paper.

5. Conclusions

This paper discusses the synchronization control issue for asymmetric bilateral teleoperation systems with time delays. Due to the discrete states caused by unreliable communication networks, new task-space-based discrete sliding mode controllers are designed. The new reaching law of DSMC consisting of the index term and tanh function is developed to reduce the chattering significantly. In addition, the input controllers are formed by undelayed position signals, delayed reference signals, and AESO-based term to compensate the total disturbance associated with the master and slave robots model dynamic. For the sake of stabilizing AESO, a new parameter θ is introduced in the form of a matrix to adjust the corresponding gain for each state separately. It is shown that estimation errors of AESO and tracking errors of teleoperation system are bound in a certain range, which is not affected by time delays and flexible friction forces. Finally, the simulated and experimental results are provided to reveal the validity and superiority of the proposed strategy.

Author Contributions

Conceptualization, Y.Y. (Yongli Yan) and Y.Y. (Yana Yang); methodology, Y.Y. (Yongli Yan) and Y.Y. (Yana Yang); validation, L.D., F.L. and Y.Y. (Yana Yang); formal analysis, L.D. and F.L.; investigation, Y.Y. (Yongli Yan), L.D., F.L. and Y.Y. (Yana Yang); Writing—Original draft preparation, Y.Y. (Yongli Yan); Writing—Review and editing, Y.Y. (Yongli Yan) and Y.Y. (Yana Yang); supervision, L.D. and F.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China, grant number [Grant Nos. 2017YFC0110503, 2017YFC0110501] And The APC was funded by [Grant Nos. 2017YFC0110503, 2017YFC0110501].

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

{ } m , { } s Subscripts { m , s } denote the master and slave robots, respectively
{ } s Subscripts s denote the slave robots
q i Joint position of the master ( i = m ) /slave ( i = s ) robot in joint space
q ˙ i Joint velocity of the master ( i = m ) /slave ( i = s ) robot in joint space
q ¨ i Joint acceleration of the master ( i = m ) /slave ( i = s ) robot in joint space
χ i , χ ˙ i , χ ¨ i Position, velocity, acceleration of the master ( i = m ) /slave ( i = s ) task-space end effectors
M q i ( q i ) Symmetric positive definite inertia matrix of the master ( i = m ) /slave ( i = s ) robot in joint space
M i ( q i ) Symmetric positive definite inertia matrix of the master ( i = m ) /slave ( i = s ) robot in task space
C q i ( q i , q ˙ i ) q ˙ i Coriolis/centrifugal matrix of the master ( i = m ) /slave ( i = s ) robot in joint space
C i ( q i , q ˙ i ) q ˙ i Coriolis/centrifugal matrix of the master ( i = m ) /slave ( i = s ) robot in task space
g q i ( q i ) Gravitational torque of the master ( i = m ) /slave ( i = s ) robot in joint space
g i ( q i ) Gravitational torque of the master ( i = m ) /slave ( i = s ) robot in task space
f q i ( q ˙ i ) Viscous friction vector of the master ( i = m ) /slave ( i = s ) robot in joint space
f i ( q ˙ i ) Viscous friction vector of the master ( i = m ) /slave ( i = s ) robot in task space
B q i ( q i ) Disturbance vector of the master ( i = m ) /slave ( i = s ) robot in joint space
B i ( q i ) Disturbance vector of the master ( i = m ) /slave ( i = s ) robot in task space
F h , F e Forces exerted on the end-effectors of the master and slave robots by the human operator and environment in joint space, respectively
τ i ( q i ) Control torque of the master ( i = m ) /slave ( i = s ) robot in joint space
u i Input control torque of the master ( i = m ) /slave ( i = s ) robot in task space
J i ( q i ) Jacobian matrix of the master ( i = m ) /slave ( i = s ) robot
T i Time delays from the master robot to the slave robot and from the slave robot to the master robot in task space, respectively
X i Position and speed state vector set of the master ( i = m ) /slave ( i = s ) robot in task space
X d i Position tracking in task space
Θ i Total disturbance vector of the master ( i = m ) /slave ( i = s ) robot in task space
A State coefficient matrix in the system equation
B Unknown coefficient matrix in the system equation
{ ¯ } Discrete form of corresponding coefficient matrix
{ ˜ } The expanded form of the corresponding coefficient matrix
h Sampling time of the master ( i = m ) /slave ( i = s ) robot in task space
j Jth sampling moment of the master ( i = m ) /slave ( i = s ) robot in task space
R Real number vector set
R n N-dimensional real number vector set
R n × n N-row and n-column real matrix set
R n × 2 n N-row and 2n-column real matrix set
I Identity matrix with appropriate dimensions
{ } T Superscript T denotes transpose matrix
{ } 1 Superscript -1 denotes inverse matrix
{ ^ } Diacritical mark wedge denotes the estimation
L i j Observer gain at the jth sampling time of the master ( i = m ) /slave ( i = s ) robot in task space
α , β Lowercase Greek letters indicate positive definite constant values
e i Position synchronization error of the master ( i = m ) /slave ( i = s ) robot in task space
Λ i Sliding mode switching parameters of the master ( i = m ) /slave ( i = s ) robot in task space
k Kth joint in task space
s i Sliding mode surface of the master ( i = m ) /slave ( i = s ) robot in task space
s i Δ i Sliding surface switching band of the master ( i = m ) /slave ( i = s ) robot in task space
Δ i Bandwidth of sliding mode surface switching band of the master ( i = m ) /slave ( i = s ) robot in task space
λ i Switching gain of sliding surface of the master ( i = m ) /slave ( i = s ) robot in task space
Φ i Reaching law gain of sliding mode surface of the master ( i = m ) /slave ( i = s ) robot in task space
Euler norm

References

  1. Clement, G.; Vertut, J.; Fournier, R.; Espiau, B.; Andre, G. An overview of CAT control in nuclear services. IEEE Int. Conf. Robot. Autom. 1985, 2, 713–718. [Google Scholar]
  2. Yoerger, D.; Slotine, J.J. Supervisory control architecture for underwater teleoperation. ICRA 1987, 4, 2068–2073. [Google Scholar]
  3. Kobayashi, H.; Ohnishi, K. Realization of virtual force sensation through bilateral teleoperation. In Proceedings of the 9th IEEE International Workshop on Advanced Motion Control, Istanbul, Turkey, 27–29 March 2006. [Google Scholar]
  4. Diez, S.P.; Borghesan, G.; Joyeux, L.; Meuleman, C.; Deprest, J.; Stoyanov, D.; Ourselin, S.; Vercauteren, T.; Reynaerts, D.; Poorten, E.V. Evaluation of Haptic Feedback on Bimanually Teleoperated Laparoscopy for Endometriosis Surgery. IEEE. Trans. Biomed. Eng. 2019, 66, 1207–1221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Niemeyer, G.; Slotine, J.J.E. Stable adaptive teleoperation. Proc. Am. Control Conf. 2009, 16, 1186–1191. [Google Scholar]
  6. Anderson, R.J.; Spong, M.W. Bilateral control of teleoperators with time delay. IEEE Trans. Autom. Control 1989, 34, 494–501. [Google Scholar] [CrossRef]
  7. Kebria, P.M.; Khosravi, A.; Nahavandi, S.; Shi, P.; Alizadehsani, R. Robust Adaptive Control Scheme for Teleoperation Systems with Delay and Uncertainties. IEEE Trans. Cybern. 2019, 50, 3243–3253. [Google Scholar] [CrossRef]
  8. Baranitha, R.; Mohajerpoor, R.; Rakkiyappan, R. Bilateral Teleoperation of Single-Master Multislave Systems with Semi-Markovian Jump Stochastic Interval Time-Varying Delayed Communication Channels. IEEE Trans. Cybern. 2019, 1–11. [Google Scholar] [CrossRef]
  9. Hokayem, P.F.; Spong, M.W. Bilateral teleoperation: An historical survey. Automatica 2006, 42, 2035–2057. [Google Scholar] [CrossRef]
  10. Huang, J.; Shi, Y.; Wu, J. Transparent virtual coupler design for networked haptic systems with a mixed virtual wall. IEEE ASME Trans. Mechatron. 2012, 17, 480–487. [Google Scholar] [CrossRef]
  11. Lu, Z.; Huang, P.; Liu, Z. Relative Impedance-Based Internal Force Control for Bimanual Robot Teleoperation with Varying Time Delay. IEEE Trans. Ind. 2019, 67, 778–789. [Google Scholar] [CrossRef]
  12. Huang, P.; Dai, P.; Lu, Z.; Liu, Z. Asymmetric wave variable compensation method in dual-master-dual-slave multilateral teleoperation system. Mechatronics 2018, 49, 1–10. [Google Scholar] [CrossRef]
  13. Takegaki, M.; Arimoto, S. A new feedback method for dynamic control of manipulators. J. Dyn. Syst. Meas. Control 1981, 103, 119–125. [Google Scholar] [CrossRef]
  14. Cheah, C.C.; Hirano, M.; Kawamura, S.; Arimoto, S. Approximate jacobian control with task-space damping for robot manipulators. IEEE Trans. Autom. Control 2004, 49, 752–757. [Google Scholar] [CrossRef]
  15. Zakerimanesh, A.; Hashemzadeh, F.; Torabi, A.; Tavakoli, M. Controlled synchronization of nonlinear teleoperation in task-space with Timevarying Delays. Int. J. Control Autom. 2019, 17, 1875–1883. [Google Scholar] [CrossRef]
  16. Dixion, W.E. Adaptive regulation of amplitude limited robot manipulators with uncertain kinematics and dynamics. IEEE Trans. Autom. Control 2007, 52, 488–493. [Google Scholar] [CrossRef] [Green Version]
  17. Wang, H. Adaptive control of robot manipulators with uncertain kinematics and dynamics. IEEE Trans. Autom. Control 2017, 62, 948–954. [Google Scholar] [CrossRef] [Green Version]
  18. Cheng, Y.; Wen, G.; Du, H. Design of Robust Discretized Sliding Mode Controller: Analysis and Application to Buck Converters. IEEE Trans. Ind. 2020, 50, 3243–3253. [Google Scholar] [CrossRef]
  19. Kali, Y.; Saad, M.; Doval-Gandoy, J.; Rodas, J.; Benjelloun, K. Discrete Sliding Mode Control based on Exponential Reaching Law and Time Delay Estimation for an Asymmetrical Six-Phase Induction Machine Drive. IET Electr. Power Appl. 2019, 3, 1660–1671. [Google Scholar] [CrossRef]
  20. Golo, G.; Milosavljevi, C. Robust discrete-time chattering free sliding mode control. Syst. Control Lett. 2000, 41, 19–28. [Google Scholar] [CrossRef]
  21. Ma, H.; Li, Y. A Novel Dead Zone Reaching Law of Discrete-Time Sliding Mode Control with Disturbance Compensation. IEEE Trans. Ind. 2020, 67, 4815–4825. [Google Scholar] [CrossRef]
  22. Reza, A.M.; Mahdi, S.; Selina, P. Adaptive discrete second-order sliding mode control with application to nonlinear automotive systems. Dyn. Syst. Meas. Control 2018, 140, 121010. [Google Scholar]
  23. Khandekar, A.A.; Patre, B.M. Discrete sliding mode control for robust tracking of time-delay systems. Syst. Sci. Control Eng. 2014, 2, 457–464. [Google Scholar] [CrossRef]
  24. Zhuang, H.; Sun, Q.; Chen, Z.; Jiang, Y. Back-stepping sliding mode control for pressure regulation of oxygen mask based on an extended state observer. Automatica 2020, 119, 109106. [Google Scholar] [CrossRef]
  25. Zhang, T.; Xu, Z.; Li, J.; Zhang, H.; Gerada, C. A Third-Order Super-Twisting Extended State Observer for Dynamic Performance Enhancement of Sensorless IPMSM Drives. IEEE Trans. Ind. 2020, 67, 5948–5958. [Google Scholar] [CrossRef]
  26. Liu, Z.; Wang, Y.; Liu, S.; Li, Z.; Zhang, H.; Zhang, Z. An Approach to Suppress Low-frequency Oscillation by Combining Extended State Observer with Model Predictive Control of EMUs Rectifier. IEEE Trans. 2019, 34, 10282–10297. [Google Scholar] [CrossRef]
  27. Lee, D.; Spong, M.W. Passive bilateral teleoperation with constant time delay. IEEE Trans. Robot. 2006, 22, 269–281. [Google Scholar] [CrossRef]
  28. Yang, X.; Hua, C.; Yan, J. Synchronization analysis for nonlinear bilateral teleoperator with interval time-varying delay. Int. J. Robust Nonlinear Control 2014, 25, 2142–2161. [Google Scholar] [CrossRef]
  29. Yang, X.; Hua, C.; Yan, J. New stability criteria for networked teleoperation system. Inf. Sci. 2013, 233, 244–254. [Google Scholar] [CrossRef]
  30. Xue, W.; Bai, W.; Yang, S. ADRC with Adaptive Extended State Observer and Its Application to Air-fuel Ratio Control in Gasoline Engines. IEEE Trans. Ind. 2015, 62, 5847–5857. [Google Scholar] [CrossRef]
  31. Kim, K.; Jee, G.; Park, C.; Lee, J. The stability analysis of the adaptive fading extended Kalman filter using the innovation covariance. Int. J. Control Autom. 2009, 7, 49–56. [Google Scholar] [CrossRef]
  32. Petersen, I.R. A Stabilization Algorithm for a Class of Uncertain Linear Systems. Syst. Control Lett. 1987, 8, 351–357. [Google Scholar] [CrossRef]
  33. Reif, K.; Gunther, S.; Yaz, E. Stochastic stability of the discrete-time extended Kalman filter. IEEE Trans. Autom. Control 1999, 44, 714–728. [Google Scholar] [CrossRef]
  34. Chen, S.; Ho, D.W.; Li, L.; Liu, M. Fault-Tolerant Consensus of Multi-Agent System with Distributed Adaptive Protocol. IEEE Trans. Cybern. 2017, 45, 2142–2155. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Control block diagram of teleoperation system with adaptive extended state observer (AESO)-based discrete sliding mode control scheme.
Figure 1. Control block diagram of teleoperation system with adaptive extended state observer (AESO)-based discrete sliding mode control scheme.
Sensors 20 05091 g001
Figure 2. Human force of the master.
Figure 2. Human force of the master.
Sensors 20 05091 g002
Figure 3. Position estimation error of the master.
Figure 3. Position estimation error of the master.
Sensors 20 05091 g003
Figure 4. Position estimation error of the master.
Figure 4. Position estimation error of the master.
Sensors 20 05091 g004
Figure 5. Disturbance estimation errors of the master.
Figure 5. Disturbance estimation errors of the master.
Sensors 20 05091 g005
Figure 6. Disturbance estimation errors of the slave.
Figure 6. Disturbance estimation errors of the slave.
Sensors 20 05091 g006
Figure 7. Position tracking motion for the master and slave.
Figure 7. Position tracking motion for the master and slave.
Sensors 20 05091 g007
Figure 8. Position tracking motion for master and slave with PD controller.
Figure 8. Position tracking motion for master and slave with PD controller.
Sensors 20 05091 g008
Figure 9. The nonlinear bilateral teleoperation experiment platform.
Figure 9. The nonlinear bilateral teleoperation experiment platform.
Sensors 20 05091 g009
Figure 10. Master position estimation errors with T m = 200 ms.
Figure 10. Master position estimation errors with T m = 200 ms.
Sensors 20 05091 g010
Figure 11. Slave position estimation errors with T s = 200 ms.
Figure 11. Slave position estimation errors with T s = 200 ms.
Sensors 20 05091 g011
Figure 12. Master disturbance estimation with T m = 200 ms.
Figure 12. Master disturbance estimation with T m = 200 ms.
Sensors 20 05091 g012
Figure 13. Slave disturbance estimation with T s = 200 ms.
Figure 13. Slave disturbance estimation with T s = 200 ms.
Sensors 20 05091 g013
Figure 14. Control torque of the master with T m = 200   ms .
Figure 14. Control torque of the master with T m = 200   ms .
Sensors 20 05091 g014
Figure 15. Control torque of the slave with T s = 200   ms .
Figure 15. Control torque of the slave with T s = 200   ms .
Sensors 20 05091 g015
Figure 16. Position tracking motion of the master and the slave with T m = T s = 200   ms .
Figure 16. Position tracking motion of the master and the slave with T m = T s = 200   ms .
Sensors 20 05091 g016
Figure 17. Master position estimation with T m = 500   ms .
Figure 17. Master position estimation with T m = 500   ms .
Sensors 20 05091 g017
Figure 18. Slave position estimation with T s = 500   ms .
Figure 18. Slave position estimation with T s = 500   ms .
Sensors 20 05091 g018
Figure 19. Master position estimation errors with T m = 500   ms .
Figure 19. Master position estimation errors with T m = 500   ms .
Sensors 20 05091 g019
Figure 20. Slave position estimation errors with T s = 500   ms .
Figure 20. Slave position estimation errors with T s = 500   ms .
Sensors 20 05091 g020
Figure 21. Control torque of the master T m = 500   ms .
Figure 21. Control torque of the master T m = 500   ms .
Sensors 20 05091 g021
Figure 22. Control torque of the slave with T s = 500   ms .
Figure 22. Control torque of the slave with T s = 500   ms .
Sensors 20 05091 g022
Figure 23. Position tracking motion of the master and the slave with T m = T s = 500   ms .
Figure 23. Position tracking motion of the master and the slave with T m = T s = 500   ms .
Sensors 20 05091 g023

Share and Cite

MDPI and ACS Style

Yan, Y.; Ding, L.; Yang, Y.; Liu, F. Discrete Sliding Mode Control Design for Bilateral Teleoperation System via Adaptive Extended State Observer. Sensors 2020, 20, 5091. https://doi.org/10.3390/s20185091

AMA Style

Yan Y, Ding L, Yang Y, Liu F. Discrete Sliding Mode Control Design for Bilateral Teleoperation System via Adaptive Extended State Observer. Sensors. 2020; 20(18):5091. https://doi.org/10.3390/s20185091

Chicago/Turabian Style

Yan, Yongli, Li Ding, Yana Yang, and Fucai Liu. 2020. "Discrete Sliding Mode Control Design for Bilateral Teleoperation System via Adaptive Extended State Observer" Sensors 20, no. 18: 5091. https://doi.org/10.3390/s20185091

APA Style

Yan, Y., Ding, L., Yang, Y., & Liu, F. (2020). Discrete Sliding Mode Control Design for Bilateral Teleoperation System via Adaptive Extended State Observer. Sensors, 20(18), 5091. https://doi.org/10.3390/s20185091

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