# Distributed Synchronization Control to Trajectory Tracking of Multiple Robot Manipulators

Authors: Yassine Bouteraa; Jawhar Ghommam; Gérard Poisson; Nabil Derbel
Journal: Journal of Robotics
 (2011)

---

## Abstract

This paper investigates the issue of designing decentralized control
laws to cooperatively command a team of general fully actuated
manipulators. The purpose is to synchronize their movements
while tracking a common desired trajectory. Based on
the well-known consensus algorithm, the control strategy consists
in synchronizing the joint position and the velocity of each
robot in the network with respect to neighboring robots' joints
and velocities. Modeled by an undirected graph, the cooperative
robot network requires just local neighbor-to-neighbor information
exchange between manipulators. So, it does not assume
the existence of an explicit leader in the team. Based above all
on combination of Lyapunov direct method and cross-coupling
strategy, the proposed decentralized control law is extended to an
adaptive synchronization control taking into account parameter
uncertainties. To address the time delay problems in the network
communication channels, the suggested synchronization control
law robustly synchronizes robots to track a given trajectory. To
this end, Krasovskii functional method has been used to deal
with the delay-dependent stability problem. A real-time software
simulator is developed to visualize the robot manipulators
coordination.


---

## Body

## 1. Introduction

There has been a great research effort for synchronization problems, in distributed cooperative control strategies where robot control laws are coupled and each robot control is updated using local rule based on its own sensors and the states of its neighbors [1, 2]. Consensus algorithms, distributed coordination, and passivity-based output synchronization for networked Euler-Lagrange systems have been studied in [3–5], respectively. In this context, one recent representative work [6] shows that we can synchronize the multicomposed system when only position measurements are available. Synchronization framework that can be directly applied to cooperative control of multiagent systems in robotic manipulation and teleoperation has been presented in [1]. Specifically, the design based on the graph theory and the Laplacian matrix produces interesting results [7–9]. Adaptive control is an effective strategy used to address the synchronization problem [10, 11].In the literature, most of earlier works on multiagent coordination and consensus [2, 8, 9, 12] mainly deal with very simple dynamic models such as linear systems and focuses on algorithm taking the form of first-order dynamics or second-order dynamics without nonlinear inertia matrices [13–15]. Most previous works on consensus and coordination of multiagent systems using the graph theory and the Laplacian [2, 8, 9, 16, 17] have presented a synchronization to the weighted average of initial conditions but they do not consider multiagent systems where there is a desired path to follow.The objective of this paper is to develop a synchronized trajectory-tracking control of multiple robot manipulators. The proposed controller relies principally on a consensus algorithm for systems modeled by nonlinear second-order dynamics and applies the algorithm for the synchronization control problem by choosing appropriately information states on which consensus is reached. The concept key of the new synchronizing controller is the introduction of a state vector that quantifies the coordination degree between robot manipulator positions and different positions of its neighbors, using cross-coupling technique. Robot manipulators are widely used in production processes. In tasks that can not be fulfilled by a single robot, either because of the complexity of the task or the spatial and temporal limits of the robot, the use of cooperative robots proved to be a good determination [18]. The tasks being executed by each robot may terminate in different time due to nonhomogenity of robots or even electronics that drives motors of robots’ joints. Despite the common trajectory robots have to track, it is not almost sure that the final running time of each robot will definitely be the same. So, in order to ensure this same final running time, robots have to have image of the current state of the neighboring robots as such to synchronize their movement together.This paper presents a distributed control strategy. Through local interactions, the proposed approach can achieve more efficient performances, in particular in the presence of external disturbances. The proposed method based on a combination of Lyapunov direct method and concepts of cross-coupling, through a convenient mathematical manipulation, has converted the motion control problem of multiple systems into a stabilization problem for one single system. Compared to existing approaches, main contributions of this paper can be stated as follows.(1)
The present work deals with highly nonlinear systems. While [13–15], to name a few, deal with simple dynamic models.(2)
In contrast to [2, 8, 9, 16, 17] the proposed approach achieves not only global asymptotical synchronization of the configuration variables, but also global asymptotical convergence to the desired trajectory.(3)
The theory is extended to adaptive control and time delay control.(4)
In [10, 11, 19], the position synchronization error of each robot is defined as the differential position error between this robot and its two adjacent robots. However, the proposed design allows interconnections between all robots, such that all robots have direct influence in the combined dynamics, and consequently the synchronous behavior is the result of interactions between all robots. Moreover, it is straightforward to show that the particular choice of the tracking error surface in the proposed approach provides an exponential convergence of qi to qd if we have a constant tracking error surface.(5)
Using the Virtual Reality Modeling Language (VRML) a virtual world was developed to simulate the robot synchronization application in 3D scenes.


## 2. Preliminaries

The dynamic equation of a general rigid link manipulator havingn degrees of freedom in a free space can be written as (1)Mi(qi)q̈i+Ci(qi,q̇i)q̇i+gi(qi)=τi,
where i(1≤i≤p) denotes the ith robot index in the network and p is the total number of the individual elements. In addition, qi∈ℝn denotes the vector of generalized displacements of the ith robot coordinates, and τi∈ℝn denotes the vector of generalized control input torques in robot coordinates; Mi(qi)∈ℝn×n inertia matrix which is symmetric uniformly bounded and positive definite, Ci(qi,q̇i)q̇i∈ℝn is a vector function containing Coriolis and centrifugal forces, and gi(qi)∈ℝn is a vector function consisting of gravitational forces. According to [20, 21], we have some fundamental properties of motion equations.(i) The inertia matrixMi(qi) is symmetric, positive definite, and uniformly bounded: (2)Mi(qi)=MiT(qi)>0.(ii) Using a proper definition,Ṁi(qi)-2Ci(qi,q̇i) is a skew symmetric matrix, satisfying (3)XT(Ṁi(qi)-2Ci(q,q))X=0,
where XT is the transpose of a vector X∈ℝn.(iii) The Euler-Lagrange equation (1) is linear with respect to the structural parameter θ, hence, (4)Mi(qi)q̈i+Ci(qi,q̇i)q̇i+gi(qi)=Y(qi,q̇i,q̈i)θi,
where Y∈ℝn×a is the regressor matrix composed of known functions of q, q̇ and q̈, θi∈ℝa is the vector of structural parameters of the manipulator, and a is the number of unknown parameters. In the present topology, the edge represents bidirectional communication links. This consists of a group of p manipulators interchanging information that can be viewed as an undirected graph (Figure 1).Figure 1
Multirobot system under mutual synchronization scheme.


## 3. Controller Design

We design decentralized control laws forp robot manipulators such that all joint positions mutually synchronize and track a common desired trajectory. The control objective of the proposed synchronization controller scheme is to synchronize the ith-joint position and velocity qi, q̇i to the state of any manipulator qj, q̇j. Besides the controller is required to regulate the joint position qi to track a desired trajectory qd. Specifically, the control torque for the ith-robot is to control the tracking error to converge to zero and at the same time to synchronize its motion with respect to motions of the p-1 robots in the network, so that the synchronization error converges to zero. To this end, we define the tracking error surface of the ith manipulator as (5)ε1i(t)=qi(t)-qd(t)+∫t0tΛi[qi(λ)-qd(λ)]dλ,
where Λi is a diagonal positive-definite matrix. Information on the vector ε1i will give insight on the convergence of the joint positions to the desired trajectory. However, it is straightforward to show that this particular choice of this tracking error surface provides an exponential convergence of qi to qd if we have a constant tracking error surface. But the defined sliding error only guarantees the trajectory tracking. Nevertheless, it is required to know the performance of the controller, that is, to know how the trajectory of each robot manipulator converges with respect to each other. There are various ways to choose the synchronization error. For example in [6], authors include the error information of all systems involved in the synchronization. In [10], authors use the cross-coupling concept to solve the synchronization problem. Our approach will make use of the cross-coupling to propose a feasible and efficient synchronization error, which consists on a measure of the synchronization for robot manipulators as defined as follows: (6)ε2i(t)=∑j≠ipKij[(qi-qd)-(qj-qd)]=∑j≠ipKij(qi-qj),
where Kij is a symmetric positive-definite matrix that gives insight on the communication quality between the ith and jth robot manipulators. We note that Kij≠0 if and only if there is information exchange between robot manipulators i and j. Consequently, each robot is not necessarily aware of all other robots. Therefore, we define the synchronizing tracking error surface which encompasses both synchronization error and trajectory tracking error for manipulator i as: (7)ei=ε1i+∫t0tε2i(λ)dλ.
The objective is to design a control law such that coupling errors, that is, position errors, velocity errors, and synchronization errors, all converge to zero. For each manipulator, the control law τi is defined as follows: (8)τi=Ci(qi,q̇i)q̇i+gi(qi)+Mi(qi)[q̈d-Kpiei-Kdiėi-Λi(q̇i-q̇d)]+Mi(qi)[∑j≠iKij(ε2i-ε2j)+Λiε2i],
where qd is a common trajectory reference to be tracked, which is a smooth time-varying trajectory and for which the first and the second derivative exist for all t≥0. Kdi and Kpi are symmetric positive-definite matrices.Theorem 1.
IfKdi>2∑j≠iKij, the proposed synchronization tracking controller (8) guarantees asymptotic convergence to zero of position errors, velocity errors, and synchronization errors, that is,  ei→0 as  t→∞.Remark 1.
Note that the synchronization controller only requires local information of the joint positionqi, the joint velocity and the information about the desired trajectory.   Also in order to ensure synchronization with other robots it needs information on the joint position from neighbors. It has to be noticed that the joint position of the neighboring robots are lumped into the error variable ε2j.Proof.
Notice that the controller law (8) contains PD controller terms. However, these terms are premultiplied by the inertia matrix Mi(qi). Therefore it is clear that this is not a linear controller as the PD control law, since the position and velocity gains are not constant but they depend explicitly on the position error.
Substituting (8) in (1) yields:
(9)Mi(qi)q̈i=Mi(qi)[q̈d-Kpiei-Kdiėi-Λi(q̇i-q̇d)]+Mi(qi)[∑j≠iKij(ε2i-ε2j)+Λiε2i].
Multiplying by Mi-1 in both sides yields
(10)q̈i=q̈d-Kpiei-Kdiėi-Λi(q̇i-q̇d)+∑j≠iKij(ε2i-ε2j)+Λiε2i.
Adding ε̇2i in both sides yields
(11)q̈i-q̈d+Λi(q̇i-q̇d)+ε̇2i=-Kpiei-Kdiėi+∑j≠iKij(ε2i-ε2j)+Λiε2i+ε̇2i.
This results in
(12)ëi=-Kpiei-Kdiėi+∑j≠iKij(ε2i-ε2j)+Λiε2i+ε̇2i.
Using the expression of the synchronization error ε2i and its first derivative gives
(13)ëi=-Kpiei-Kdiėi+∑j≠iKij(ε2i-ε2j)+∑j≠i[Kij(q̇i-q̇d)-(q̇j-q̇d)]+∑j≠iΛiKij[(qi-qd)-(qj-qd)].
Further calculation will result in
(14)ëi=-Kpiei-Kdiėi+∑j≠iKij(ėi-ėj).
Equation (14) represents the closed-loop synchronized system for the ith manipulator. In the sequel, we proceed to analyze the stability properties of the proposed synchronized control scheme and ultimately show that the control goals: the position error, velocity error, and synchronization error, all converge to zero. To prove stability of the overall synchronized system, we define eT=[e1T,e2T,…,enT]T.
Using (14) we obtain the synchronized error dynamics:
(15)ë=-Kpe-(Kd-Kc)ė,
where Kp=diag(Kpi),Kd=diag(Kdi), and Kc is given by
(16)(∑j≠1K1j⋯-K1j⋯-K1n⋯⋯⋯⋯⋯-Ki1⋯∑j≠iKij⋯-Kin⋯⋯⋯⋯⋯-Kn1⋯-Knj⋯∑j≠nKnj).
Note that Kc is symmetric and positive semidefinite matrix.
The synchronized error dynamics (15) is a linear time invariant system described by a second-order linear differential equation. A sufficient condition for the error dynamics to be exponentially stable is that the matrices Kp and Kd-Kc are positive definite. In particular, matrices Kdi can be diagonal satisfying Kdi>2∑i≠jKij.
To analyze the stability properties of the closed-loop error dynamics (15), we take the following definite and radially unbounded Lyapunov function candidate:
(17)V=ėTė+eTKpe.
Its derivative with respect to time can be expressed as
(18)V̇=-2ėT(Kd-Kc)ė≤0.
It follows by direct application of LaSalle's invariance principle that the origin (e,ė)=(0,0) is globally asymptotically stable and limė→0 for t→∞.
Referring to the expression of the global error (7)
(19)ei=qi-qd+∫t0tΛ[qi(λ)-qd(λ)]dλ+∫t0t∑j≠iKij(qi-qj)
as ė=0 we have
(20)q̇i-q̇d=-Λi(qi-qd)-∑j≠iKij(qi-qj).
We set εi=qi-qd. Then (20) can be written as
(21)ε̇i=-Λiεi-∑j≠iKij(εi-εj).
Our objective is to show that limεi converges to zero as t↦∞. To this end, we define ε=[ε1T⋯εiT⋯εnT]T and Λ=diag(Λ1⋯Λi⋯Λn). Equation (21) can be written as
(22)ε̇=A⋅ε,
where matrix A is given by
(23)(-Λ1-∑j≠1K1j⋯K1j⋯K1n……………Ki1⋯-Λi-∑j≠iKij⋯Kin⋯⋯⋯⋯⋯Kn1⋯Knj⋯-Λn-∑j≠nKnj).
Consider the following Lyapunov function:
(24)v(t)=εTε.
Differentiating v(t) with respect to time yields
(25)v̇=2∑i=1nεiTε̇i=2∑i=1nεiT(-Λiεi-∑j≠iKij(εi-εj))=-2∑i=1nεiTΛiεi-2∑i=1n∑j≠iεiTKij(εi-εj)=-2∑i=1nεiTΛiεi-2∑i=1n∑j≠iεiTKijεi+2∑i=1n∑j≠iεiTKijεj.
Knowing that
(26)∑i=1n∑j≠iεiTKijεi=∑j=1n∑i≠jεjTKjiεj,
consequently,
(27)v̇=-2∑i=1nεiTΛiεi-∑i=1n∑j≠i(εi-εj)TKij(εi-εj)≤0.
It follows by direct application of LaSalle's invariance that the origin is globally asymptotically stable. Consequently we obtain limεi(t)→0 for t→∞. Then qi→qd and q̇i→q̇d for t→∞.
Referring to (19), we show that qi→qj for t→∞.


## 4. Adaptive Synchronization

In this section, we consider an uncertainty in the model parameters. We propose to extend the previous decentralized control law to an adaptive version. For that, we consider the following adaptive control law, which has the similar local coupling structure as the proposed control law in (8): (28)τi=Ĉi(qi,q̇i)+ĝi(qi)+M̂i(qi)×[q̈d-Kpiei-Kdiėi-Λi(q̇i-q̇d)]+M̂i(qi)[∑j≠iKij(ε2i-ε2j)+Λiε2i].
We recall that a similar relation to (1) holds when the estimate parameter vector θ̂ is used to replace the exact parameter vector θ: (29)M̂i(qi)q̈i+Ĉi(qi,q̇i)q̇i+ĝi(qi)=Yi(qi,q̇i,q̈i)θ̂i,
where θi∈ℝa is the vector of structural parameters of the manipulator and a is the number of unknown parameters. The estimated parameter θ̂i is subject to the adaptation law: (30)θ̂̇i=-Γi-1(M̂i-1Yi)Tėi,
where Γi is a diagonal positive-definite matrix. Since the value of the dynamic parameter θi is hard to be known exactly in practice, one defines θ̂i(t) as the estimate of θi. M̂i, Ĉi, and ĝi are estimates of Mi, Ci, and gi, respectively. Yi(qi,q̇i,q̈i) denotes a regression matrix. Define θ̃i=θi-θ̂i as a vector containing model estimation errors. Then the adaptation law (30) can be written as (31)θ̃̇i=Γi-1(M̂i-1Yi)Tėi.Theorem 2.
IfKdi>2∑j≠iKij, the proposed adaptive coupling controllers (28) and (31) guarantee the asymptotic convergence to zero of the joint position, velocity, and synchronization errors.Proof.
Substituting (28) into the dynamic model (1) leads to the following closed-loop dynamics:
(32)Mi(qi)q̈i+Ci(qi,q̇i)q̇i+gi(qi)=M̂i(qi)[∑j≠iKijq̈d-Kpiei-Kdiėi-Λi(q̇i-q̇d)+∑j≠iKij(ε2i-ε2j)+Λiε2i]+Ĉi(qi,q̇i)q̇i+ĝi(qi).
We define
(33)M̃i(qi)=Mi(qi)-M̂i(qi),C̃i(qi,q̇i)=Ci(qi,q̇i)-Ĉi(qi,q̇i),g̃i(qi)=gi(qi)-ĝi(qi).
Taking into account (33), (32) gives:
(34)M̃i(qi)q̈i+C̃i(qi,q̇i)q̇i+g̃i(qi)+M̂i(qi)(q̈i-q̈d)=M̂i(qi)[∑j≠iKij-Kpiei-Kdiėi-Λi(q̇i-q̇d)+∑j≠iKij(ε2i-ε2j)+Λiε2i].
Using (29), the dynamic (34) is written as
(35)Yi(qi,q̇i,q̈i)θ̃i+M̂i(qi)(q̈i-q̈d)=M̂i(qi)[∑j≠iKij-Kpiei-Kdiėi-Λi(q̇i-q̇d)+∑j≠iKij(ε2i-ε2j)+Λiε2i].
Consequently,
(36)M̂i(qi)ëi+Yi(qi,q̇i,q̈i)θ̃i=M̂i(qi)∑j≠iKij[-Kpiei-Kdiėi+∑j≠iKij(ėi-ėj)].
To prove the stability of the overall synchronized system, we define YT=[Y1T,…,YiT,…,YnT]T, θT=[θ1T,…,θiT,…,θnT]T, Γ=diag(Γ1T,…,ΓiT,…,ΓnT)Kd=diag(Kdi), Kp=diag(Kpi), M̂=diag(M̂i).
Writing (36) in a compact form gives the following:
(37)M̂[ë+Kpe+Kdė-Kcė]=-Y(qi,q̇i,q̈i)θ̃.
Multiplying by M̂-1 both sides yields
(38)ë=-Kpe-Kdė+Kcė-M̂-1Y(qi,q̇i,q̈i)θ̃.
Define the Lyapunov function candidate as
(39)v=12(ėTė+eTKpe+θ̃TΓθ̃).
Differentiating v(t) with respect to time yields
(40)v̇=ëTė+eTKpė+θ̃TΓθ̃̇.
Substituting (38) into (40)
(41)v̇=(-Kpe-Kdė+Kcė-M̂-1Yθ̃)Tė+eTkpė+θ̃TΓθ̃̇,
consequently
(42)v̇=-ėT(Kd-Kc)ė+θ̃T[Γθ̃̇-(M̂-1Y)Tė].
Using the adaptation law (31) gives
(43)v̇=-ėT(Kd-Kc)ė.
Since Kdi>∑j≠iKij, we have v̇(t)<0, and this yields that v(t)<v(0), which gives that e,ė, and M̂-1Yθ̃T are bounded.
Differentiatingv̇(t) with respect to time yields
(44)v̈=2(Kp+(Kd-Kc)ė+M̂-1Yθ̃)T(Kd-Kc)ė.
Using Barbalat's lemma, v̈ is bounded because e, ė, and M̂-1Yθ̃ are bounded. This implies v̇(t)→0 for t→∞ and hence ė→0 for t→∞. The proof pursued the same line reasoning as the proof of Section 3. Consequently, we show that position errors, velocity errors and synchronization errors, converge asymptotically to zero.


## 5. Synchronization with Time Delays

In this section, we study the coordination control problem taking into account time delays of communication channels. As a first assumption, we suppose that these delays can be justified by the fact that data information sent by the neighboring vehiclesj≠i reaches vehicle i after a certain time delay due to the short-range communication channels. To take into account time delays produced during the communication among robots, we introduce, in the coordination error expression, a term τ which represents the same time delay due to the short-range communication channels. Therefore, a coordination error, in the time delay context, will be presented as the well-known classical time delay model of multiagent network [9]: (45)ε2i(t)=∑j≠iKij[qi(t-τ)-qj(t-τ)].
It will be shown that the behavior of the coordinated system under the effect of time delay changes significantly. Consequently, the controller implanted in each lagrangian system among the network takes the following expression.Theorem 3.
Implemented controller taking the following expressions(i=1,…,p):
(46)τi=Ci(qi,q̇i)q̇i+gi(qi)+Mi(qi)×[q̈d-Kpiei(t)-Kdiėi(t)-Λi(q̇i-q̇d)]+Mi(qi)[∑j≠iKij[ε2i(t-τ)-ε2j(t-τ)]+Λiε2i(t-τ)∑j≠iKij[ε2i(t-τ)-ε2j(t-τ)]]
stabilize the behavior of the robot network.Proof.
Substituting (46) into (1) yields
(47)Mi(qi)q̈i=Mi(qi)[q̈d-Kpiei-Kdiėi-Λi(q̇i-q̇d)]+Mi(qi)[∑j≠iKij[ε2i(t-τ)-ε2j(t-τ)]+Λiε2i(t-τ)∑j≠iKij[ε2i(t-τ)-ε2j(t-τ)]].
By further calculation, we obtain the synchronized error dynamics:
(48)ë=-Kpe-Kdė+Kcė(t-τ),
where Kp, Kd, and Kc are the same matrices already defined (see Section 3). By the Leibnitz formula, we have
(49)ė-ė(t-τ)=∫t-τtë(λ)dλ.
Substituting (49) into (48) leads to
(50)ë=-Kpe-Kdė+Kc(ė-∫t-τtë(λ)dλ).
Setting ẽ=[eT,ėT]T. Therefore (50) can be written as
(51)ẽ̇=(0I-Kp-Kd+Kc)ẽ-(000Kc)∫t-τtẽ̇(λ)dλ.
This yields the following form:
(52)ẽ̇=β0ẽ-β1∫t-τtẽ̇(λ)dλ
with β0=(0I-Kp-Kd+Kc) and β1=(000Kc).
To analyze the stability of the global system, we consider the following Lyapunov-Krasovskii functional (LKF):(53)v(t)=v1(t)+v2(t)+v3(t),v1(t)=ẽT(t)Pẽ(t),v2(t)=∫t-τtẽT(λ)Rẽ(λ)dλ,v3(t)=∫-τ0(∫t+stẽ̇T(α)Zẽ̇(α)dα)dλ,
where P=PT>0, R=RT>0, Z=ZT>0 are weighting matrices of appropriate dimensions. A straightforward computation gives the time derivative of v(t) along the solution of (53) as
(54)v̇(t)=δT[2NTPM+NTRN-QTRQ+τMTZM]δ-∫t-τtẽ̇T(λ)Zẽ̇(λ)dλ,
where δ=[ẽT(t),ẽT(t-τ)]T, N=[β0,β1], M=[I,0], Q=[0,I]. The Jensen's inequality gives a suitable bound for the last term of (54):
(55)-∫t-τtẽ̇T(λ)Zẽ̇(λ)dλ≤∫t-τtẽ̇T(λ)dλ(Zτ)∫t-τtẽ̇(λ)dλ≤-ẽTTT(Zτ)Tẽ
with T=[I,-I]. Thus, the time derivative of the LKF (53) can be bounded by v̇(t)≤δTξδ, where
(56)ξ=2NTPM+NTRN-QTRQ+τMTZM-1τTTZT.
Then if the LMI ξ<0 is satisfied, the derivative of the Lyapunov-Krasovskii functional is negative definite. To ensure matrix ξ is negative definite, we select appropriate control gains Kp>Kp* and Kd-Kc>K* throught processing MATLAB's LMI solver such that
(57)2NTPM+NTRN-QTRQ+τMTZM-1τTTZT<0.
Then, if the LMI ξ<0 is satisfied, the derivative of the Lyapunov-Krasovskii functional is therefore negative definite. In consequence, the origin ẽ=0 is asymptotically stable.
This results ine→0 for t→∞ and ė→0 for t→∞. The proof for the asymptotic convergence of the coordinated tracking error ẽ is not sufficient to prove the convergence to zero of both errors e1 and e2. Now, our concern is is to show that coordination is successfully realized for a specific time delay τc.
The proof pursued the same line reasoning as the proof of Section3. Consequently, we obtain the following equation derived from the global error expression:
(58)ε̇i=-Λiεi-∑j≠iKij[εi(t-τ)-εj(t-τ)].
Rewriting all states of (58) into a compact representation and applying the Laplace transform lead to
(59)sε(s)-ε(0)=-Λε(s)-e-τsKcε(s).
This can be written as
(60)ε(s)=(sI+Λ+e-τsKc)-1ε(0).
If the characteristic equation P(s,τ)= det|sI+Λ+e-τsKc|=0 has all its zeros in the left half complex plan, then the system is stable and one can easily conclude about the convergence of qi to qd. Since the ordinal system, free from time delay (i.e, τ=0), is a continuous function of τ, then using the D-decomposition, the minimal positive solution to the following equation:
(61)τi=1(Kci2-Λi2)arccos(-ΛiKci)
which would make all the zeros of the characteristic equation in the left half complex plane. Therefore if we select τ∈[0,τc], where τc=sup{τi,∀1≤i≤p}, therefore solutions of (60) converge to zero and consequently qi→qd, q̇i→q̇d, and qi→qj for t→∞.


## 6. Simulation Results

To show the effectiveness of the proposed synchronizing controller, we run two types of simulation.Simulation 1
Simulations are performed on MATLAB/Simulink. These simulations were performed on identical three manipulators with two degrees of freedom. Robots do not have the same starting positions, and they start their motion att0=0. First, we assume that no parameter uncertainty is present and that the transmission delay can be neglected. Figure 2 illustrates the robots synchronization tracking a common trajectory. This proves that tracking and synchronization objectives are attained by the proposed controller which explains how robots, while tracking the desired trajectory, synchronize their position with others. Next, we assume that the length and the mass of manipulators are unknown. Consequently, mass, Coriolis, and gravity matrices are ill known. Figure 3 plots an adaptive position tracking while in Figures 4 and 5 tracking errors and synchronizing errors are shown, respectively. Finally, we consider the time delay in communication. Figures 6 and 7 illustrate that the behavior of the coordinated system changes significantly, under the effect of time delays in communication channels. The speed for achieving an agreement depends essentially on the time delay of communication channels.Figure 2
Synchronization of robots without uncertain parameters.Figure 3
Synchronization of robots in the presence of uncertain parameters.Figure 4
Position errors: adaptive synchronization.Figure 5
Synchronization errors: adaptive synchronization.Figure 6
Robots synchronization without time delay.Figure 7
Robots synchronization in presence of time delay.Simulation 2
One effective solution is to develop a platform to simulate the complex system design before hardware implementation. To design this virtual environment, we apply several tools, such as the modeling package Cinema 4D, MATLAB VRML (Virtual Reality Modeling Language) toolbox, and MATLAB/Simulink. Robot CAD (Concept Architecture Design) models created using 4D cinema are extracted in VRML97 (.wrl) format. The VRML file is edited by V-Realm Builder 2.0 before being imported into MATLAB-Simulink by the VR (Virtual Reality) Toolbox. A VRML file uses a standard text format which can be read with any text editor. For MATLAB, the Virtual Reality Toolbox software includes the Ligos V-Realm Builder application as a native 3D editor. But the 3D editing tools remain the most efficient. In fact, the CAD packages offer the power and versatility to create many types of practical models and techniques. In the developed real-time software simulator we used the 4D Cinema as CAD package where the robots are designed. Next, the 3D objects are imported. Simulink model generates signal data which is used to control and animate the virtual world. Reality Toolbox block is added to the Simulink model in order to ensure communication with the virtual world. The 3D virtual reality environment is presented in Figure8.The 3D virtual environment.
(a)(b)(c)(d)(e)


## 7. Conclusions

This paper has considered the synchronization problem in distributed Multirobot systems under cooperative schemes. The aim of this work is to find out a decentralized controller, applied to each manipulator; the synchronization is therefore met. It has been shown that the proposed strategy can coordinate manipulator articulations to track a given time-varying trajectory. As an extension, we moved to the step which entails the fact that mass, Coriolis, and gravity matrices are badly known. In this case, the proposed controller has encompassed an adaptive version. The obtained simulation results from Multirobot motion control system demonstrate the effectiveness of the synchronization approach. To deal with time delays problem in communication between robots, the proposed decentralized control guarantees that information variables of each robot reach agreement even in the presence of communication channels delays. A new package for simulation of coordinated robot manipulators is developed. Simulations which have been applied to an illustrative example have shown the effectiveness of the described strategy.



---

*Source: 652785-2011-08-23.xml*