Abstract

A procedure for the guaranteed cost control design of delayed bilateral teleoperation systems with nonlinear external forces is proposed. The assumption that the external forces are nonlinear functions of velocities and/or positions of local devices and that one part of these forces satisfies sector condition has been made. A stability criterion is formulated firstly and then the optimal guaranteed cost controller is obtained by solving the equivalent convex optimization problem in the form of linear matrix inequalities (LMIs). The controller preserves closed-loop stability regardless of the delay length. The behavior of the resulting teleoperation system is illustrated in simulations.

1. Introduction

Teleoperation systems enable humans to extend their capacity to manipulate remote interfaces with better safety, at less cost, and even with better accuracy. Their rich applications vary from nuclear operations and space and underwater exploration to medical surgery; see [1, 2] and references therein. In such an application, a human operator uses an interface referred to as a “master device” to conduct tasks by remotely controlling a “slave device” at the distant or hazardous environment. The general bilateral system setup is depicted in Figure 1. The main objectives of the control design for bilateral teleoperation systems are position coordination between the master and the slave and accurate haptic display of environment forces to the human operator. However, due to the very nature of teleoperation, time delays associated with communication between the local and the remote sites are inevitable. Handling the delays is especially problematic in the case of bilateral teleoperation, where the measurements are communicated in both directions to allow haptic feedback. The commonly proposed approaches to deal with bilateral teleoperation with time delays in the communication channel are mainly based on scattering theory formalism [3] or wave variable concept [4], which are called passivity-based control approaches. These controllers render the communication link passive and thus guarantee stable bilateral teleoperation. More recent research focused on improving tracking performance for constant and time varying delays and increasing the transparency of the teleoperator system and also the extensions to task space and network system control. For example, constant delays in the communication channel were considered in [2, 5], while [68] dealed with the communication channel with time varying time delays. Data loss and packet-switched network were considered in [9], where a passivity-based framework was proposed, and extension to the case where communication channel is a packet-switched network was also provided. Reference [10] used a modified Smith predictor, a kalman filter, and an energy regulator to control a teleoperator over Internet-based communication. Both quantization and time delay issues for teleoperation systems over networks were studied in [11]. Control design for bilateral teleoperation systems in task space was considered in [12, 13], where the task space position tracking was guaranteed. Another research related to task space control design of teleoperation systems can be found in [14, 15]. In addition to the above mentioned results, passivity-based controllers relying on damping injection [5, 16], passive output synchronization [16, 17], and adaptive control [18] were also developed, recently. Several other approaches were proposed to deal with delays and are reviewed in [19].

Passivity-based controllers have the advantage of being able to guarantee stability independent of the user and environment dynamics; however, when controlling a real teleoperator, it is desirable to design a controller which not only guarantees the stability of the closed-loop system but also guarantees adequate level of performance. In this work a guaranteed cost control design is proposed for teleoperation systems with delays to guarantee a specific level of performance. We suppose that the external forces are nonlinear, while one part of them satisfies sector condition. Note that only one part of these forces satisfies a sector condition with respect to the output, which means that other parts may not be passive. Thus passivity is actually not required in our case. The position/velocity tracking performance is optimized by using the information of sector bounds. Constant and symmetric communication time delays are assumed to simplify the analysis; extensions for time varying delays are similar and easy to be realized.

The rest of the paper is organized as follows. In Section 2, the teleoperation model is generalized and then the problem formulation is given in Section 3. The controller design and its main results are given in Section 4, following which the simulation result is given to show the effectiveness of the proposed method in Section 5. Finally, the conclusion is included in Section 6.

The notation used throughout the paper is standard. denotes the set of -dimensional real vectors. For a real symmetric matrix , the notation of is used to denote positive (negative) definiteness. is used to denote the identity matrix of appropriate dimensions, is the zero matrix of appropriate dimension, and is used to represent a block matrix which is readily referred to by symmetry. For a matrix , denotes its transpose. For any function , the -norm is defined as and the -norm is defined as . and spaces are defined as the sets and , respectively.

2. Preliminaries: Bilateral System Modeling

Consider a simple teleoperation system as follows:where , , are the joint positions, velocities, and accelerations of the master and slave devices with or representing master or slave robot manipulators, respectively. Similarly, , , are the effective mass, damping, and spring coefficients of the master and slave devices. External forces applied to the devices by the human operator and the environment are represented by , , respectively, while and are the Jacobian matrices of the master and the slave, respectively. , stand for the control signals.

Defining , for each , and , a minimal state-space realization of system (1) is reformulated aswhere

3. Problem Formulation

The problem formulation based on the following assumptions is considered:(A1)The forward/backward time delays in the communication channels are the same and are denoted by .(A2)The positions and velocities are available for measurement.

Remark 1. Assumptions A1 and A2 are very general in the control of teleoperation systems. Assumption A1 is for the simplification of the analysis, while the results in this paper can be easily extended to the case with asymmetric time varying delays.

The states and the input of the teleoperation system are defined as and the output is defined asand the evolution of the states is governed bywith , where Note that in dynamics (6) is included here for the convenience of stability analysis, and the details will be found in the following sections.

The following assumption is furthermore given:(A3) can be modeled aswhere is continuous uniformly bounded function satisfying and is a nonlinear function of and it is piecewise continuous in , globally Lipschitz in with , while satisfying the following sector condition (a) or (b) for . Consider(a)(b)where are constant real matrices of appropriate dimensions.

Remark 2. Equations (11) and (12) are special cases of the prevailing passivity-type assumption. The sector bounds and are introduced to provide bounding information of the nonlinearities in the external forces, which also make the following optimization steps possible. Note that external forces do not satisfy passivity assumption due to the time varying . However, when , Assumption A3 reduces to the prevailing passivity-type assumption.

The state feedback controllers for the teleoperation system (1) or (6)–(8) are proposed as follows:where are matrices to be designed. Hence the control law can be rewritten aswhere .

Remark 3. Note that pretended time delays are intentionally introduced in system output (5) and controllers (13) and (14); by doing this, it avoids the case that the resulting setup follows into a group of controller design problems with intrinsically quadratically invariant constraints on the controller parameters; see [20] for more discussions. The basic idea is to design the guaranteed cost controller for system (6)–(8) under the worst case; that is, there are delays both in the local and in the remote variables; then better performance can be achieved when local variables actually have no time delays. For a practical teleoperation system of which local state variables do not have time delays, the stability and performance will be guaranteed as well. In the simulation, the time delays at the local sides were removed to show the stability.

To introduce the main objective of this paper, a cost function associated with system (6)–(8) is introduced:where are given positive-definite symmetric matrices.

Hence, the main objective of this paper is to design a guaranteed cost control law (14) for the teleoperation system (1) with the human force and environmental force satisfying Assumption A3 such that the closed-loop system is stable with an associated cost less than a guaranteed cost . The problem setup is shown in Figure 2.

4. Controller Design

The case when the nonlinear function satisfies (11) is firstly considered and the following result is obtained. After this, the result will be extended to the general case in which (12) is satisfied.

Theorem 4. For the teleoperation system (6)–(8) with external forces satisfying Assumption A3 with the simple sector condition (11), the guaranteed cost controller (14) exists if there exist symmetric positive-definite matrices and scalars such that

Proof. Choose a Lyapunov-Krasovskii functional candidate aswhere . Taking derivative of with respect to along the trajectory of system (6)–(8) yieldsDenote , and one hasFrom (11), one hasEmploy the inequality obtained by completion of squares:where , and then from (19), (20), and (21) one obtains that hence (11) implies thatBy integrating both sides of inequality (23) from to , one hasNow, in view of (16) and (24), one can conclude that is negative if , and thus this implies that , , and hence , . Furthermore, one could easily get that This completes the proof.

Based on the stability criterion above, the original problem of designing a guaranteed cost controller (14) can be transformed to a convex optimization problem, which is presented below.

Theorem 5. For the teleoperation system (6)–(8) with external forces satisfying Assumption A3 with the simple sector condition (11), the control law (14) is an optimal memoryless guaranteed cost control law which ensures the minimization of the guaranteed cost (15) if the following convex optimization problem,subject to (1)(2)the matrices which are symmetric, and(3) (4)has a solution . Here, denotes the trace of the matrix and . In addition, the controller gains in (14) are given as .

Proof. It follows from the Schur complement that (16) is equivalent toPre- and postmultiplying both sides of (31) by yield LMI (27), where . By applying Theorem 3 in [21], convex optimization problem (26) was obtained. It completes the proof.

For the nonlinearity satisfying the more general sector condition (12), by applying the idea known as loop transformation [22], one can conclude that the absolute stability of system (6)–(8) with the nonlinearity within the sector is equivalent to that of the following system:where satisfies the sector condition for . Hence by Theorem 4 the following result is obtained.

Proposition 6. For the teleoperation system (6)–(8) with external forces satisfying Assumption A3 with the sector condition (12), the guaranteed cost controller (14) exists if there exist symmetric such thatwhere

Now, following Theorem 5, the guaranteed cost controller (14) for system (6)–(8) with nonlinearities (12) can be obtained by solving the optimization problem (26) with constraints (28), (29), (30), and (35).

Remark 7. The main results in this paper provide a delay-independent stability criterion for the nonlinear system (6)–(8). If the Lyapunov-Krasovskii functional is chosen as , where and is defined in (17), a delay-dependent stability criterion for the system can be obtained.

5. Simulations and Results

Consider a master-slave teleoperation system (1) with the effective mass Kg, the damping  Ns/m, and the spring coefficients  N/m.

Hence it can be described by (6)–(8) with and the initial states are chosen as with the time delay  s. Here the master and the slave with different initial positions and velocities are chosen to see the tracking performance.

To access the teleoperation system behavior, the human force and the environment force are firstly chosen to be linear; that is, with as the external input which is chosen to be the rectangle signal depicted in Figure 3, and , . The coefficients are chosen as , , , and . Hence, the nonlinearity satisfies (11) with The cost function (15) is chosen with

By Theorem 5 and solving the corresponding optimization problem by using YALMIP toolbox [23], the optimal guaranteed cost control law with the solution is obtained.

Noting the structure of the control gains in (42), one can rewrite (14) with parameter (42) asIt can be seen that the optimal controllers are Proportional-Derivative (PD) controllers. Note that the time delays in the local variables are omitted since in the simulations there is no time delay in the local variables. Applying controller (43) to the described system, the simulation results shown in Figures 4(a) and 4(b) are obtained. It can be seen that the master and the slave response stably. The positions and the velocities of the slave follow the master’s motion. Specifically, the master and the slave move forward and backward with the increase and decrease of the human force, respectively. Let , the value of the cost can also be tracked in Figure 4(c), and the best cost is .

Now nonlinearities into the external forces are introduced; that is,where is the inverse function of and , .

By Theorem 5 the optimal guaranteed cost control law with the solutionis obtained and thus the equivalent controllers in the traditional form are It can be seen that the optimal controllers are Proportional (P) controllers. Simulation results for the system with nonlinear external forces (44) are presented in Figures 5(a) and 5(b). There is still no time delay existing in the local variables, which is more realistic. It can be seen that the simulation results for the case when the external forces are nonlinear are similar to the results for the case when the external results are linear. This implies that the controller is valid for the system with nonlinear external forces. The value of the cost is also tracked in Figure 5(c), and the best cost is .

6. Conclusion

In this paper, a guaranteed cost control design procedure for delayed bilateral teleoperation systems is proposed. The assumption that there exist time delays in the local devices has been made to obtain the desired controllers under the worst case, and thus the quadratically invariant problem setup is avoided. Furthermore, only one part of the external forces is assumed to satisfy some section conditions, and thus the guaranteed cost controller is introduced. To obtain the guaranteed cost control gains, we firstly formulate a stability criterion and then convert it to a convex optimization problem which can be solved by linear matrix inequality (LMI) technique. Simulations for a simple example show that the stability is guaranteed with very good performance when the designed controller is applied to the systems which do not actually have delays in the local devices. It is also found that the actual controllers obtained from our method are PD controllers or P controllers, but there is no need to pick any PD or P parameters in our method while guaranteed level of performance is satisfied. It is very beneficial for the control of teleoperation systems with high degrees of freedom in practical applications. Further study will focus on the bilateral teleoperation systems with complicated, nonlinear manipulator models, varying time delays, and parameter uncertainties.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work was supported by Key Program of National Natural Science Foundation of China (Grant no. 61333002).