## Servicios Personalizados

## Articulo

## Indicadores

- Citado por SciELO

## Links relacionados

- Citado por Google
- Similares en SciELO
- Similares en Google

## Bookmark

## Latin American applied research

*versión impresa* ISSN 0327-0793

### Lat. Am. appl. res. v.36 n.2 Bahía Blanca abr./jun. 2006

**Teleoperation of mobile robots**

**E. Slawiñski, V. Mut and J. F. Postigo**

*Instituto de Automática (INAUT). Universidad Nacional de San Juan. Av. Libertador San Martín 1109 (oeste). J5400ARL. San Juan, Argentina. e-mail:{slawinski, vmut, jpostigo}@inaut.unsj.edu.ar*

*Abstract* — This paper proposes a stable control structure for bilateral teleoperation of mobile robots. The proposed control structure includes a time-delay compensation placed on both the local and remote sites of the teleoperation system. Teleoperation experiments through a simulated and real (using Internet) communication channel are presented to illustrate the performance and stability of the proposed control structure.

*Keywords* — Asymptotic Stability. Mobile Robots. Teleoperation. Time Varying Delay.

**I. INTRODUCTION**

Teleoperation systems have been developed to allow human operators to execute tasks in remote or hazardous environments, with a variety of applications ranging from space to underwater, nuclear plants, and so on (Sheridan, 1995). In general, the bilateral teleoperation systems of mobile robots are composed by a local site (where a human operator drives a hand-controller device); a remote site (where a mobile robot interacts with the physical world); and a communication channel that links both sites. The human operator generates velocity commands to the remote mobile robot, while the position of the mobile robot is back-fed to the human operator through the communication channel.

Perhaps, the most interesting case appears when there exits a distance between the local and remote sites of a teleoperation system. This generally introduces time varying delays adding distortion in the reference commands and feedback signals. The presence of time delay may induce instability or poor performance of a teleoperation system (Fiorini and Oboe, 1997; Richard, 2003). In general, in the design of teleoperation systems there is a trade-off between high transparency and sufficient stability margins (Lawrence, 1993). Main control strategies proposed for bilateral teleoperation of delayed systems are described in Anderson and Spong (1989), Niemeyer and Slotine (1991), Oboe and Fiorini (1998), Oboe (2003), Elhajj *et. al.* (2003) and Chopra and Spong (2003). In general, the proposed control structures keep the passivity or stability at the expenses of reducing the system transparency (Arcara and Melchiorri, 2002).

This paper proposes a stable control structure for bilateral teleoperation of mobile robots. The proposed control structure is based on combining the velocity command generated by the human operator in a delayed time instant, the received position information (which stimulates the operator) in such moment and the current position of the mobile robot to set the velocity reference of the mobile robot. The time proposed delay compensation is placed on the local and remote sites. Moreover, experiences of teleoperation of a mobile robot are shown to test the stability and performance of the designed teleoperation system.

The paper is organized as follows: Section II gives the notation used in this paper. In section III, some background material on delayed differential equations is introduced. Section IV presents the statement of the control problem. In Section V, a model of the human operator for motion control is presented. In Section VI, a stable control structure for bilateral teleoperation of mobile robots is proposed. In section VII, the stability and performance of the proposed control structure are analyzed through teleoperation experiments using a simulated and real (Internet) communication channel. Finally, the conclusions of this paper are given in Section VIII.

**II. NOTATION**

In this paper, the following notation is used: *h*(*t*) ∈ ℜ^{+} denotes the time delay. Here, *x*(*t*) ∈ ℜ^{n}, *x ^{T}* is the transpose of

*x*, |

*x*| is the Euclidean norm of

*x*,

*x*(for a given time instant

_{t}*t*) is the function defined by

*x*(θ) =

_{t}*x*(

*t*+ θ) for θ ∈ [-

*h*(

*t*), 0], for example:

*x*(0) =

_{t}*x*(

*t*),

*x*(-

_{t}*h*) =

*x*(

*t - h*); and the norm ||.|| is defined by ||

*x*|| = ·

_{t}*C*

_{n,h(t)}is the n-dimensional space of continuous functions on the interval [

*t*-

*h*(

*t*),

*t*] at any time

*t*, then the function

*x*∈

_{t}*C*

_{n,h(t)}. On the other hand, given a non-linear differentiable function (

*t*) =

*g*(

*x*(

*t*),

*x*(

*t - h*)), the incremental gain of

*g*is defined as |

*g*| = inf{γ:|

*g*(

*x*

_{1},

*x*

_{2}) -

*g*(

*y*

_{1},

*y*

_{2})| ≤ γ |[

*x*

_{1}

*x*

_{2}] - [

*y*

_{1}

*y*

_{2}]|} ∀

*x*

_{1},

*x*

_{1},

*y*

_{1},

*y*

_{2}∈ ℜ

^{n}.

**III. STABILITY OF DELAYED SYSTEMS**

The robot teleoperation systems are represented by delayed differential equations. In this section, we show standard definitions and facts in the theory of delayed functional differential equations (Krasovskii, 1963; Hale 1977; Kolmanovskii and Myshkis, 1999). In addition, we propose a stability condition for systems with time delay, which will be used in Section VI for the stability analysis of the proposed teleoperation system.

Let's consider the delayed functional differential equation given by,

(*t*) = *f*_{1}(*t*, *x _{t}*),

*f*

_{1}(

*t*, 0) = 0, ∀

*t*≥

*t*

_{0}(1)

where *x* ∈ ℜ^{n}, *x _{t}* ∈

*C*

_{n,h(t)},

*t*,

*t*

_{0}∈ ℜ

^{+}, and

*f*

_{1}:ℜ

^{+}×

*C*

_{n,h(t)}→ ℜ

^{n}. It is assumed that there exists a solution x(

*t*;

*t*

_{0}, ψ

_{0}) of (1) with initial data [

*t*

_{0}, ψ

_{0}], where ψ

_{0}=

*x*(

*t*

_{0}+ θ) for θ ∈ [-

*h*(

*t*

_{0}), 0] with ||ψ

_{0}|| <

*H*∈ ℜ

^{+}, which depends continuously on the initial data. From now on, we will denote the solution norm by |

*x*(

*t*;

*t*

_{0}, ψ

_{0})| = |

*x*(0)|.

_{t} **Definition 1.** The solution *x _{t}* = 0 of (1) is said to be asymptotically stable if,

a) For every ε > 0 and each *t*_{0} ≥ 0 there exists ρ = ρ(ε, *t*_{0}) such that ||ψ_{0}|| < ρ implies that |*x*(*t*; *t*_{0}, ψ_{0})| < ε for all *t*_{0} ≥ 0 .

b) For every *t*_{0} ≥ 0 there exists ε_{1} = ε_{1}(*t*_{0}) such that if ||ψ_{0}|| < ε_{1}, then |*x*(*t*; *t*_{0}, ψ_{0})| → 0 as *t* → ∞.

If ρ and ε_{1} are independent from the initial time *t*_{0}, then the zero-solution is uniformly asymptotically stable.

**Fact 1 **(Krasovskii, 1963). Let's suppose that the function *f*_{1}:ℜ^{+} × *C*_{n,h(t)} → ℜ^{n} maps bounded sets of *C*_{n,h(t)} in bounded sets of |*g*| , and that *u*(.), *v*(.) and *w*(.) are scalar, continuous, positive and non-decreasing functions. If there exists a continuous functional *V*:ℜ^{+} × *C*_{n,h(t)} → ℜ^{+}, and the following conditions hold:

*u*(|*x _{t}*(0)|) ≤

*V*(

*t*,

*x*) ≤

_{t}*v*(||

*x*||) (2)

_{t} (*t*, *x _{t}*) < -

*w*(|

*x*(0)|), (3)

_{t} where (*t*, *x _{t}*) is the time-derivative of

*V*(

*t*,

*x*) along the trajectories of (1); then the solution

_{t}*x*= 0 is uniformly asymptotically stable.

_{t}Now, let us consider a non-linear system with time varying delay described by,

(*t*) = *f*_{2}(*t*, *x*(*t*)) + *g*(*t*, *x*(*t*), *x*(*t - h*)), (4)

where 0 ≤ *h*(*t*) ≤ *h _{m}* and (

*t*) < τ < 1, with

*h*∈ ℜ

_{m}^{+},

*x*∈ ℜ

^{n},

*t*,

*t*

_{0}∈ ℜ

^{+},

*f*

_{2}:ℜ

^{+}× ℜ

^{n}→ ℜ

^{n}, and

*g*:ℜ

^{+}× ℜ

^{n}→ ℜ

^{n}. In addition, we assume that

*f*

_{2}(

*t*, 0) = 0 and

*g*(

*t*, 0, 0) = 0, for ∀

*t*≥

*t*

_{0}.

**Lemma 1.** If a system represented by = *f*_{2}(*x*) is exponentially stable, then there exist α, λ ∈ ℜ^{+} such that *x ^{T}f*

_{2}(

*x*) ≤ -λ

*x*, where |

^{T}x*x*| ≤ α|

*x*(0)|

*e*.

^{-λt} **Proof.** If the system = *f*_{2}(*x*) is exponentially stable, then it satisfies that |*x*| ≤ α|*x*(0)|*e ^{-λt}*, and therefore:

|*x*| ≤ α|*x*(0)|*e ^{-λt} I*

_{1}(5) .

where *I*_{1} ∈ ℜ^{n} with *I*_{1}(*i*) = 1 for 0 < *i* ≤ *n*. From (5), the evolution of verifies that,

|*x*| ≤ -λα|*x*(0)|*e ^{-λt} I*

_{1}(6) .

Using (5) and (6) on -λ*x ^{T}x*, the following can be expressed:

-λ*x ^{T}x* ≥ -λ(α|

*x*(0)|

*e*)

^{-λt}*I*

_{1}

^{T}I_{1}≥

*x*=

^{T}*x*

^{T}f_{2}(

*x*). (7)

Inequality (7) proves the proposed Lemma. Δ

**Theorem 1.** *Let us suppose that the subsystem = f_{2}(x) of the system *(4)

*is exponentially stable with rate λ, then the following condition ensures the asymptotic stability of the system*(4):

(8) |

*where λ, |g| ∈ ℜ ^{+} and (t) < τ < 1. The norm |g| is the incremental gain of the operator g*(

*.*)

*.*

**Proof.** A functional *V*:ℜ × *C*_{n,h(t)} → ℜ^{+} is proposed as follows,

(9) |

where the proposed functional incorporates information of the delayed dynamics (|*g*|) that will help to reach a stability condition that will directly depend on the time-derivative of the varying-time delay and the non-delayed dynamics of the delayed system.

From (9) and considering that the delay *h*(*t*) is bounded (*h*(*t*) ≤ *h _{m}*) and that

*x*= |

^{T}x*x*(0)|

_{t}^{2}≤ ||

*x*||

_{t}^{2}(by using norm properties), then the proposed functional

*V*(

*t*,

*x*) verifies condition (2) –given by

_{t}*Fact 1*-,

(10) |

The time-derivative of *V*(*t*, *x _{t}*) along the system trajectories (4) is,

(11) |

Now, the following inequalities are attained using norm properties,

(12) |

Putting (12) in (11), it yields,

(13) |

The third term of the right hand in (13) is negative definite because < τ < 1. By applying *Lemma 1* to (13)

and organizing terms, it yields,

(14) |

From (14), condition (3) –*Fact 1*- is satisfied if,

(15) |

Inequalities (10) and (14) verify the stability conditions given by *Fact 1* –inequalities (2) and (3)-. Then, the proposed *Theorem 1* is proven ensuring the asymptotic stability of the system (4). Δ

Figure 1 shows the effect of the maximum derivative of the time delay on the stability region -given by (15)- for three arbitrary values |*g*| = 0.5, |*g*| = 1 and |*g*| = 2.

Fig. 1. Stability region in function of τ.

The achieved stability condition is independent of the delay amplitude and it depends on three main factors: the exponential rate λ of the non-delayed system = *f*_{2}(*x*), the norm |*g*| of the delayed non-linear function *g*(*x*, *x*(*t*-*h*)) and the maximum time-derivative τ of the time delay. Moreover, the greater the temporal derivative of the time delay (τ), the stronger the stability of the non-delayed system (higher λ) to reach the stability of the system with time delay. In addition, if |*g*| → 0 then the proposed stability condition tends to the stability condition of a non-delayed system, this is: λ < 0.

**IV. STATEMENT OF THE CONTROL PROBLEM**

This section describes the analysed control problem on a bilateral teleoperation system of mobile robots. Figure 2 shows a general diagram of a teleoperation system.

Fig. 2. General block diagram of a teleoperation system of a mobile robot.

The human operator drives a mobile robot through a hand-controller generating velocity commands to send to the remote site, which will be executed by the mobile robot. The mobile robot and obstacles position is visually back-fed to the human operator. We suppose that the obstacles position generates a fictitious force, which depends on the distance between the mobile robot and the obstacle.

The main signals of the system are the position *x _{r}* and force

*f*on the remote site, the received position

_{r}*x*and force

_{l}*f*on the local site, the velocity command

_{l}*v*generated in the local site and the velocity reference

_{l}*v*applied to the mobile robot. On the other hand, the communication channel is represented by a time delay

_{r}*h*composed by a forward delay

*h*

_{2}(from the local site to the remote site) and a backward delay

*h*

_{1}(from the remote site to the local site), i.e.,

*h*(*t*) = *h*_{1}(*t*) + *h*_{2}(*t*). (16)

We will consider the mobile robot as a unicycle located at a non-zero distance from the objective frame <g>. In addition, attached to the robot there exists the frame <a>, as shown in Fig. 3.

Fig. 3. Position and orientation of a mobile robot.

We consider the vehicle position in Polar Coordinates, where the state variables are the polar coordinates ρ, α, θ measured between the frame <g> and the frame<a>. The kinematic equations can be written as,

(17) |

Where *v _{rv}*,

*v*are the linear and angular velocities of the mobile robot.

_{rω} The objective of the teleoperation system is that a human operator (placed on the local site) drives a mobile robot (placed on the remote site) to reach the frame <g> in spite of the time varying delay, this is, that the distance error (state) –in this case, without final orientation- x := [ρ, α] → 0 as *t* → ∞ starting from any non-zero distance from <g>.

**V. MODEL OF MOTION CONTROL OF THE HUMAN OPERATOR**

This section presents a model for the motion control of the human operator, which will be used later (Section VI) by the proposed delay compensation.

**A. Human operator's model for position control**

The kinematic model proposed for the position controller of the human operator, which generates velocity commands *v _{l}* = [

*v*,

_{lv}*v*], is the following (Slawiñski

_{lω}*et. al.*, 2005):

(18) |

where *k _{v}*,

*k*> 0. Introducing the human controller (18) into the kinematic equations of the mobile robot (17), we obtain the following closed loop equations for the state

_{ω}*x*:= [ρ, α]:

(19) |

**Lemma 2.** The non-delayed teleoperation system (given by (19)) of a mobile robot (17) driven by a human controller represented by (18) is exponentially stable with rate λ = min{*k _{v}*,

*k*}.

_{ω} **Proof.** The proposed Lyapunov candidate function is,

(20) |

The time derivative of *V*(ρ α) along the trajectories of the system (19) is,

= *-k _{v}*ρ

^{2}cos

^{2}α -

*k*

_{ω}α

^{2}. (21)

__Remark 1__: The time-derivative of the functional *V*(ρ, α) is negative definite (21), then the trivial solution is globally asymptotically stable.

__Remark 2__: From (19), the solution for α is . The initial condition has a range given by |α(0)| ≤ π.

The problem to establish an exponential response on ρ is resolved by steps:

a) If the initial condition is |α(0)| ≤ *d*π, where *d* is a positive arbitrary constant lower than 0.5-, then (21) can be expressed as,

< *-k _{v}*' ρ

^{2}-

*k*α

_{ω}^{2}, (22)

where *k _{v}*' =

*k*cos

_{v}^{2}(

*d*π) > 0. From (20) and (22) is simple to deduce that the non-delay system is exponentially locally stable with exponential rate

*k*= min{

_{m}*k*',

_{v}*k*}.

_{ω} b) If the initial condition is |α(0)| > *d*π, then there exists a finite time *T _{m}* defined by (from

*Remark 2*) which assures that α(

*T*) ≤

_{m}*d*π.

c) We propose that the response of ρ is bounded by,

. (23)

d) If *t* ≤ *T _{m}*, from (23) and considering that the system is globally asymptotically stable –

*Remark 1*-, it yields,

e) If *t* > *T _{m}*, then α(

*t*) <

*d*π -from (b)-. Then, from (a) the response of ρ is bounded by the exponential response given by (23).

__Remark 3__: From steps (a), (b), (c), (d), (e) and *Remark 2*, the equilibrium point *x* = [ρ α]^{T} = 0 is exponentially stable with rate λ = *k _{m}*. We choose

*d*near zero, such that

*k*' ≈

_{v}*k*> 0, then λ =

_{v}*k*= min{

_{m}*k*,

_{v}*k*}. Δ

_{ω} **B. Fictitious force**

The mechanical impedance regulation needs the feedback from the interaction force between the robot and its environment. The interaction forces imply physical contacts with the environment which, in the case of mobile robots, means a collision. To avoid obstacles, however, it's necessary to interact with the environment without causing any collision. In such case, the interaction force is represented by a fictitious force, which depends on the distance between the robot and the obstacle, as shown the Fig. 4.

Fig. 4. Impedance control with fictitious force.

The magnitude of the fictitious force *f* is computed as *f*(*t*) = *a - bd*(*t*), where *a*, *b* are positive constants such that *a - bd _{m}* = 0,

*d*is the robot-obstacle maximum distance, and

_{m}*d*(

*t*) is the robot-obstacle distance (0 ≤

*d*(

*t*) ≤

*d*), which is measured through ultrasonic or type-camera sensors. On the other hand, the angle of the fictitious force is β (see Fig. 4). The fictitious force on the remote site is

_{m}*f*:= [

_{r}*f*]

_{t}f_{n}^{T}(see Fig. 2).

**C. Reactive control mode with decision of the human operator**

The impedance model of the human operator is defined by *Z = Bs + K*, where *B*, *K* are positive constants; while the reference error is defined by = *Z*^{-1}*f _{t}*, where

*f*cos β is the component of

_{t}) f*f*on the robot motion direction. The reference error is transformed to a rotation angle ψ =

*D*(

*t*) applied on the position reference (Mut

*et. al.*, 2002); where

*f*sin β is the component of

_{n}= f*f*normal to the robot's motion direction and

*D*(

*t*) represents the human operator's decision. If the environment and the task are perfectly known, then the decision could be predicted. We assume that the task and environment are according to

*D = sign*(

*f*). When the fictitious force is zero, the reference error is zero too, and then the objective of the motion control is achieved.

_{n} **D. Reactive control mode of the human operator**

When the environment or the task are not known, then the fictitious force modifies the distance error ρ and the angular error α as: , where *f _{t}* is the component of

*f*on the robot motion direction,

*f*is the component of

_{n}*f*normal to the robot's motion direction, and the impedance model of the human operator is defined by

*K*

^{-1}=

*diag*[

*K*

_{ρ}

*K*

_{α}] ∈ ℜ

^{2×2}, where

*K*

_{ρ},

*K*

_{α}> 0 represent the human operator's elasticity in response to fictitious force generated by the distance robot-obstacle.

**E. Experimental validation of the human operator's model to drive the mobile robot**

Figure 5 shows the executed trajectories by the mobile robot driven by a human operator in two different experiences and also the trajectory using an automatic control, it is composed by both the position controller position described in sub-section A, and the impedance controller described in sub-section C. We conclude that the proposed model of the human operator is satisfactory. The impedance loop (the desired impedance is represented by a stable and proper strictly linear filter) only modifies the reference of the motion control. However, it won't be considered later to simplify the stability analysis of the bilateral teleoperation system (Section VI).

Fig. 5. Trajectories of the mobile robot using manual teleoperation and automatic control.

**VI. CONTROL STRUCTURE FOR BILATERAL TELEOPERATION OF MOBILE ROBOTS**

This section describes the proposed control structure applied to a bilateral teleoperation system of mobile robots.

The proposed delay compensation does not modify the feedback position from the remote site. In addition, the local site sends a signal *v _{l}*(

*t*) - Δ

*v*(

*t - h*

_{1}) to the remote site; this signal combines the velocity command generated by the human operator in a time instant and the received position information (which stimulates the operator) in such moment. In the remote site, the proposed delay compensation uses the current position of the mobile robot to modify the signal

*v*(

_{l}*t - h*

_{2}) - Δ

*v*(

*t*- (

*h*

_{1}+

*h*

_{2})) and to establish the velocity reference

*v*(

_{r}*t*). Figure 6 shows a block diagram of the delayed bilateral system introducing the proposed time-delay compensation.

Fig. 6. Block diagram of the teleoperation system with delay compensation.

The delay compensation is placed on both the local and remote sites and it is defined by an approximated model of the local site (Section V) as follows,

(24) |

where are the parameters of the delay compensation and the vector Δ*v* = [Δ*v _{v}*, Δ

*v*] is the output of the proposed delay compensation.

_{ω} Now, we analyze the system stability using the proposed delay compensation and also considering that the local site is represented by a time-invariant kinematic model. We computed the vector *v _{r}* = [

*v*,

_{rv}*v*] (Fig. 6), which is applied on the mobile robot (17) as follows,

_{rω}(25) |

From (17), the evolution of the state *x* := [ρ, α] of the delayed system is given by,

(26) |

We put (16), (18) and (24) in (25) to obtain a intermediate equation which is incorporated in (26) describing the delayed system as follows:

(27) |

where and .

If the delay compensation is an exact model of the local site, then and therefore |*g*| → 0 in (27). Then, from (19) and (27) the system will represent the non-delayed real system and from *Lemma 2*, the delayed system will be asymptotically stable.

On the other hand, if . To simplify this analysis, we suppose that , then from (27), the incremental gain is,

(28)

From *Theorem 1* (given by (8)), *Lemma 2* and (28), the stability condition is expressed as,

(29) |

The proposed control structure allows us to separate the delayed system into the non-delay real system (*f*_{2} in (27)) and a new delayed subsystem (delayed function *g* in (27)). In the general case (|*g*| ≠ 0), the proposed control strategy allows ensuring the system stability through a condition imposed on the non-delayed system (29) that depends on the maximum derivative τ of the delay and the gain |*g*| of the delayed non-linearity of the system.

**VII. EXPERIMENTAL RESULTS**

To illustrate the performance and stability of the proposed control structure for mobile robot teleoperation, experiments have been conducted on a Pioneer 2DX mobile robot through a simulated and real communication channel. The human operator receives visual feedback of position from a Logitech webcam placed on the remote site. The objective of control is to achieve the position reference avoiding a type cylinder or cube obstacle placed on the workspace of the remote site. It should be noticed that the impedance control loop is active when the mobile robot detects an obstacle at a distance less than 1.5[*m*] using ultrasonic sensors.

**A. Teleoperation with simulated delay**

The hand-controller used in this experiment is a Logitech Wingman joystick. The initial condition is ρ(θ) = 3.7[*m*], α(θ) = 0[*rad*] for θ ∈ [-*h*(*t*_{0}), 0]. The time delay is simulated by software. The used parameters for the delay compensation are: for the force compensation (sub-section C in section V) and for the position compensation (compensation of position and force in presence of time delay).

Figure 7 shows the executed trajectories by the Pioneer 2DX mobile robot for diverse delays. The advantage of using the delay compensation is clear.

Fig. 7. Trajectories of the mobile robot using the delay compensation for various delays.

Figure 8 shows the evolution of the state norm |x| = |[ρ, α]| of the delayed teleoperation system for various delays. On the other hand, the Fig. 9 shows the evolution of the linear velocity of the mobile robot. The maximum linear velocity varies between 0.4 and 0.5 for the diverse experiments.

Fig. 8. Evolution of the norm of the state *x*.

Fig. 9. Linear velocity of the mobile robot.

The response of the delayed teleoperation system using the delay compensation is similar to the manual teleoperation without time delay (reference response); therefore the performance of the system is good.

**B. Teleoperation through Internet between Brazil and Argentina**

Now, the performance of the proposed control structure for bilateral teleoperation of a mobile robot (Pioneer 2DX) driven by a human operator through Internet between San Juan (Argentina) and Vitoria (Espírito Santo, Brazil) is presented. The hand-controller used is

a commercial steering wheel with accelerator pedal. An obstacle type-cube is placed on the workspace of the remote site.

The used parameters for the delay compensation are: for the force compensation (sub-section D in section V) and for the position compensation (compensation of position and force in presence of time delay).

Figure 10 shows the evolution of the time delay *h* and (which is estimated using a differentiator filter). Figure 11 shows the trajectory of the mobile robot for this experiment. The human operator drives the mobile robot to reach the objective position avoiding the obstacle placed in the remote site.

Fig. 10. Time-varying delay for the experiment B using Internet to link the local and remote sites.

Fig. 11. Trajectory of the mobile robot teleoperated by a human operator through Internet.

On the other, Fig. 12 shows the temporal evolution of the distance error ρ, it tends to zero as *t* → ∞.

Fig. 12. Temporal evolution of the distance error ρ.

The response of the teleoperation system is satisfactory in spite of the time varying delay added by Internet.

**C. Stability**

Now, we analyse the stability of the teleoperation system. The maximum time-derivative of the time delay added by the simulated and real (Fig. 10) communication channel is approximately τ = 0.2. On the other hand, the exponential rate of the non-delayed system is λ = min{*k _{v}*,

*k*} (see

_{ω}*Lemma 2*). From (29), we can express the stability condition on as,

(30) |

From (30) and the value of *k*_{ω}, we can conclude that the model used by the delay compensation could have parametric errors to a percentage of 20%.

**VIII. CONCLUSIONS**

In this paper it has been proposed a stable control structure for bilateral teleoperation systems of mobile robots. The proposed strategy includes a delay compensation placed on the local and remote sites of the teleoperation system and it uses a model of the local site.

Several experiments have shown a stable response with good performance and transparency. In addition, a mobile robot was driven by a human operator with visual feedback through a simulated and real communication channel in a continuous way. From these results, we may conclude that the application of the proposed control structure on an industrial or commercial system is feasible.

The future work will be incorporating the dynamic model of the human operator on the proposed control structure. In addition, the parameters of the human operator will be identified to improve the performance of the teleoperation system.

**REFERENCES**

1. Anderson, R.J. and M. Spong. "Bilateral control of Teleoperators with time delay", *IEEE Transaction and Automatic Control*, **34, **494-501 (1989). [ Links ]

2. Arcara, P. and C. Melchiorri. "Control schemes for teloperation with time delay: A comparative study", *Robotics and Autonomous Systems*, **38**, 49-64 (2002). [ Links ]

3. Chopra N. and M. Spong. "Bilateral Teleoperation over the Internet: the Time Varying Delay Problem", *Proceedings of the American Control Conference*, Denver, Colorado, June 4-6, 155-160 (2003). [ Links ]

4. Elhajj, I., N. Xi, W. Fung and Y. Liu, "Supermedia-Enhanced Internet-Based Telerobotics", *Proceedings of the IEEE*, **91**, 396-421 (2003). [ Links ]

5. Fiorini, P. and R. Oboe, "Internet-Based Telerobotics: Problems and Approaches", *ICAR'97*, Monterey, CA, July 7-9, 765-770 (1997). [ Links ]

6. Hale, J.K., *Theory of Functional Differential Equations*, Springer Verlag, New York (1977). [ Links ]

7. Kolmanovskii, V.B. and A. Myshkis. *Introduction to the theory and applications of functional differential equations*. Dordrecht: Kluwer Academy (1999). [ Links ]

8. Krasovskii, N.N., *Stability of Motion*, Stanford University Press (1963). [ Links ]

9. Lawrence, A., "Stability and Transparency in Bilateral Teleoperation", *IEEE Transactions on Robotics and Automation*, **9**, 624-637 (1993). [ Links ]

10. Mut, V., J. Postigo, E. Slawiñski and B. Kuchen. "Bilateral teleoperation of mobile robots", *Robotica*, Cambrige University Press, **20**, 213-221 (2002). [ Links ]

11. Niemeyer, G. and J.J.E. Slotine. "Stable Adaptive Teleoperation", *IEEE Journal of Oceanic Engineering*, **16, **152-162 (1991). [ Links ]

12. Oboe, R. and P. Fiorini. "A Design and Control Environment for Internet-Based Telerobotics", *The International Journal of Robotics Research*, **17**, 433-449 (1998). [ Links ]

13. Oboe, R., "Force Reflecting Teleoperation Over the Internet: The JBIT Project", *Proceedings of the IEEE*, **91**, 449-462 (2003). [ Links ]

14. Richard, J.P., "Time-delay systems: an overview of some recent advances and open problems", *Automatica*, **39**, 1667-1694 (2003). [ Links ]

15. Sheridan, T.B., "Teleoperation, Telerobotics and Telepresence: A Progress report", *Control Eng. Practice*, **3**, 205-214 (1995). [ Links ]

16. Slawiñski E., V. Mut and J. Postigo, "Bilateral teleoperation of mobile robots with delay", *Proceeding of the IEEE International Conference on Mechatronics and Automation* (ICMA 2005), Niagara Falls, Canadá, July, 1672-1677 (2005). [ Links ]

**Received: September 21, 2005. Accepted for publication: March 3, 2006. Recommended by Guest Editors C. De Angelo, J. Figueroa, G. García and J. Solsona.**