Design_of_robust_adaptive_fuzzy_control_for_uncert
Design_of_robust_adaptive_fuzzy_control_for_uncert
DOI: 10.1049/cth2.12407
ORIGINAL RESEARCH
1
Department of Computer Engineering, College of Abstract
Computer Engineering and Sciences, Prince Sattam
In this study, a novel method based on a robust adaptive fuzzy control approach is devel-
bin Abdulaziz University, Al-Kharj, Saudi Arabia
2
oped for nonlinear teleoperation systems. Its main objectives are to ensure system stability
Control and Energy Management Laboratory
(CEM Lab), Ecole Nationale d Ingenieurs de Sfax
and properly mitigating parametric uncertainties stemming from external disturbances and
(ENIS), Institut Superieur de Biotechnologie de Sfax un-modelled dynamics. For the communication channel, instead of the direct transmission
(ISBS), University of Sfax, Sfax, Tunisia of environmental torque signals, the approximated environmental parameters by the fuzzy
3
Department of Computer Science and Artificial system are transmitted to the master side for the prediction of environmental torque, thus
Intelligence, College of Computer Science and successfully avoiding the transmission of the power signals in the delayed communication
Engineering, University of Jeddah, Jeddah, Saudi
Arabia channel and solving the passivity problem in the teleoperation system. Besides, a trajec-
4
tory generator is employed in the master side, whereas a trajectory smoothing is provided
College of Science, Shanghai Institute of
Technology, Shanghai, China in the slave side. Theoretically, it was proven that both position tracking and force feed-
5
Department of Electrical and Computer
back problems are attained. Using Lyapunov stability analysis, this work illustrates that the
Engineering, University of Louisiana at Lafayette, robust adaptive fuzzy controller based on the backstepping approach guarantees the sys-
Lafayette, Louisiana, USA tem’s asymptotic stability. Simulation results confirm the efficiency of the suggested control
6
Graduate School of Intelligent Data Science, technique in achieving the stability and tracking objectives of the uncertain nonlinear
National Yunlin University of Science and teleoperation system.
Technology, Douliou, Yunlin, Taiwan
Correspondence
Saleh Mobayen, Graduate School of Intelligent Data
Science, National Yunlin University of Science and
Technology, 123 University Road, Section 3,
Douliou, Yunlin 640301, Taiwan.
Email: [email protected]
This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is
properly cited.
© 2023 The Authors. IET Control Theory & Applications published by John Wiley & Sons Ltd on behalf of The Institution of Engineering and Technology.
difficult, deteriorate the performance of the feedback-based in the presence of environmental forces a challenging prob-
control function and can even lead to instability [23, 24]. Due lem. In [45], a variable structure control approach based on
to their ability to stabilize structured and unstructured uncertain the neural network and optimized fractional-order selection
systems, fuzzy control and neural network approaches are con- law has been suggested for the sensor-less tele-robotic system
sidered as powerful instruments for the control of the robotic in the presence of uncertain time-delays, model uncertainties,
systems [25–29]. In [11], an adaptive controller is introduced fractional calculus numerical approximation bias, and external
which is capable of synchronization in the presence of dynamic disturbances. According to the integer-order calculus based on
uncertainties without the appearance of any delay in infor- the definition of Caputo and comparison principle, the Rie-
mation from the communication channel. The fuzzy control mann Liouville uniformly-ultimately-bounded (UUB) stability
approach in [30] has been proposed to ensure the sustainabil- is presented in order to synthesize the SMC technique with
ity of non-linear systems; however, these designs have mostly varying order. Combined with the mixed-type error, the pro-
focused on mobile robots and rarely addressed delayed non- posed hybrid integer and fractional-order analysis techniques
linear teleoperation systems. The teleoperation systems with are able to guarantee the Riemann–Liouville-based UUB sta-
constant known delays have been studied in [31, 32]. The bility of tracking errors in the presence of time delays. Based
time delays of these systems are practically variable. Actually, on the Euler difference, the recommended synthesis technique
the major object of teleoperation systems is to achieve robust suppresses the combined boundedness error of approximating
stability on the work environment [33, 34]. Because of some double Euler difference errors.
inherent fundamental problems of these systems, achieving Many of the previous approaches applied for master-slave
this aim seems to be complex. One of the reasons mentioned robot control are mainly based on robust adaptive fuzzy control
for such problems is the delay in the master-slave communi- of teleoperation systems [50]. Due to the simplicity of teleopera-
cation channel. This factor even results in instability of the tion systems’ dynamic equations, controlling master-slave robot
control system [35–37]. Another method for counteracting the control will be complex accordingly [51]. Besides, the use of pas-
inherent delay in such systems is the use of Smith’s estimator sivity conditions [52] in proving controllers’ stability will result
configuration [38]. In [39], the system has become stable by in some limitations.
removing delays from equation components by means of the The most important challenge in the bilateral teleoperation
Smith estimator. Anyway, there is no exact model for such non- systems is to achieve the operator’s true sense of the envi-
linear systems; however, this kind of modelling is not much ronment and guarantee the stability of the system. Until now,
reliable because of system non-linearities. The kind of system various approaches have been developed to design such systems
examined in [40] is non-linear with variable delays and sys- in which the passivity condition for environmental forces and
tem stability was examined through passivity condition and operators have been used to prove the stability of the Lyapunov
Lyapunov–Krasovskii function. Though some attempts have function. These conditions can impose serious restrictions,
been made to find a bound for the time delay, one of the weak- especially in the field of medical applications on teleoperation
nesses of this method, however, lies in the use of passivity systems. The existence of time-varying delay in the master-slave
condition which has been considered as a system limitation in communication channel and also uncertainties (external dis-
[41]. Accordingly, none of these resources have investigated the turbances and un-modelled dynamics), which sometimes may
issue of model uncertainty. destabilize the control system, are important factors and under-
Because of lengthy information communication, the inher- lie in attracting researchers’ attention. Considering the above
ent time delay within the communication channel is negligible challenges and the previous literature, a robust backstepping
[42, 43]. Master-slave synchronization [44] and stability analy- fuzzy adaptive controller for the control of uncertain non-linear
sis of the teleoperation systems under various time delays such bilateral teleoperation system in the presence of time delay and
as constant delay [11], time-variable delays [45, 46] and random parametric uncertainty. The main contributions of this paper are
delays [47] are considered as interesting topics in the field of threefold:
teleoperation systems. In [48], a delay compensator has been
proposed for teleoperation systems in the presence of feed- ∙ Introduction of a robustifying term to overcome external
back forces for control of the workspace framework. Since the disturbances and fuzzy approximation errors.
exact measurement of both kinematic and dynamic parameters ∙ A design that successfully avoids the transmission of
is considered to be a difficult task, multiple control algorithms the power signals in the delayed communication chan-
have been proposed to counteract these two difficulties. By nel and solves the passivity problem in the teleoperation
combining the adaptive and robust control techniques, the sys- system.
tem will have the advantages of both types of controllers. The ∙ A design that successfully alleviates the limitations stemming
robust adaptive control technique is useful for tracking control from the passivity problem.
of robots with unstructured dynamics. However, the control
approach is accompanied by a high computational load [49]. The remainder of the paper is structured as follows. Section 2
Another non-linear control method is the backstepping con- elaborates upon the research hypotheses and lemmas. Section 3
trol approach. Contrary to sliding mode control (SMC), it is details the proposed control approach and provides the stability
not required for the equations to be in the portable form in analysis. The simulation results are given in Section 4. Finally,
the backstepping control [1]; thus, making the robot control some conclusions are provided in Section 5.
802 BOUTERAA ET AL.
2 PROBLEM DESCRIPTION AND environmental forces. In this method, the bound for all
NECESSARY CONCEPTS uncertainties will be approximated using an adaptive rule.
The teleoperation system’s model will be as follows:
In order to design the proposed control approach, the following ( ) ( ) ( )
properties and lemmas are considered: Mm qm q̈ m + Cm qm , q̇ m q̇ m + gm qm = 𝜏m + Dm + 𝜏h
( ) ( ) ( )
Property 1. [53]: The mass inertia matrices Mm and Ms are symmetric Ms qs q̈ s + Cs qs , q̇ s q̇ s + gs qs = 𝜏s + Ds − 𝜏e (2)
positive-definite and the upper and lower bounds are defined as 0 < 𝜇m I ≤
Mm ≤ 𝜇M I and 0 < 𝜇s I ≤ Ms ≤ 𝜇S I , where 𝜇S and 𝜇m are constant such that m and s indices are master and slave robots, corre-
scalar. spondingly, qm and qs are joints’ location vectors, q̇ m and q̇ s are
joints; speed vectors, Mm (qm ) and Ms (qs ) are definite positive
Property 2. [54]: The Coriolis matrices have the following features: inertia matrices, Cm (qm , q̇ m ) and Cs (qs , q̇ s ) are Coriolis torques,
(1) The dynamic equation of the Coriolis inertia matrices gm (qm ) and gs (qs ) are gravity torques, 𝜏h and 𝜏e are human and
is considered as Ṁ i (q, q)
̇ = Ci (q, q)
̇ + CiT (q, q),̇ ∀q, q̇ ∈ Rn×1 environment torques, and finally 𝜏m and 𝜏s are control signals of
(i = m, s), where Ci is Coriolis inertia matrix. torque material. Besides, Dm and Ds are the total time-varying
(2) The conditions q̇ mT (Ṁ m − 2Cm )q̇ m = 0 and values related to un-modelled dynamics and external distur-
T ̇
q̇ s (Ms − 2Cs )q̇ s = 0 are satisfied, where the asymmetric bance. Considering the order of the teleoperation system model
property of the master-slave system (Ṁ i − 2Ci = 0) is obtained. discussed above, the following two systems can be considered
for error dynamics:
Lemma 1. [55]: For any continuous function f(x), there is a fuzzy
system as 𝜃T 𝜁(x) which holds sup | f (x) − 𝜃T 𝜁(x)| ≤ w, ∀w > Zm1 = qm − qS (t − Ts (t ))
x∈Ω
0, where 𝜃T 𝜁(x) will be defined in the following. Zs1 = qs − qm (t − Tm (t )) (3)
ẋ 1 = x2 1 T 1
V1 = Zm1 Zm1 + Zs1T Zs1 (5)
2 2
ẋ 2 = Mi−1 (−Ci x2 − gi + 𝜏i + Di + 𝜏ext ) (1)
where the time-derivative of V1 can be computed as:
where x1 = qi and x2 = q̇ i (i = m, s). Considering the above V̇ 1 = Zm1
T ̇
Zm1 + Zs1T Ż s1
state space, backstepping control is planned in two stages which ( ( ) )
require the recommendation of primary Lyapunov function = Zm1
T
q̇ m − 1 − Ṫ s (t ) q̇ s (t − Ts (t ))
(first stage) and final Lyapunov function (second stage). ( ( ) )
Note 2: To approximate the indefinite terms, fuzzy approx- + Z T q̇ s − 1 − Ṫ m (t ) q̇ m (t − Tm (t ))
s1 (6)
imator is used which is based on the single point fuzzification,
product inference machine, and centric mean defuzzification. It is noteworthy that in the communication channel, the delay
Note 3: In this study, the robustifying term has been used to term and its derivative are bounded; that is, Tm (t ) ≤ d1 , Ts (t ) ≤
overcome external disturbances and fuzzy approximation error, d2 , Ṫ m (t ) ≤ 𝜇1 ≤ 1 and Ṫ s (t ) ≤ 𝜇2 ≤ 1; where d1 and d2 are pos-
which is considered as one of the strengths of work in the field itive scalar values. Now, by substituting the time derivative of
of teleoperation systems. Equation (4) into (6), the following equation can be obtained:
( ( ) )
V̇ 1 = Zm1T
Zm2 + 𝛼m1 − 1 − Ṫ s (t ) q̇ s (t − Ts (t ))
3 RECOMMENDED CONTROL ( ( ) )
METHOD + Zs1T Zs2 + 𝛼s1 − 1 − Ṫ m (t ) q̇ m (t − Tm (t )) (7)
where Km1 and Ks1 are positive constant values. Then, by sub- equation can be re-written as:
stituting Equation (8) into (7), the following equation is found:
V̇ 2 = −Zm1
T
Km1 Zm1 − Zs1T Ks1 Zs1 + Zm1
T
Zm2 + Zs1T Zs2
( ( ) ( )
V̇ 1 = −Zm1
T
Km1 Zm1 − Zs1T Ks1 Zs1 + Zm1
T
Zm2 + Zs1T Zs2 (9) + Zm2
T
− Cm qm , q̇ m q̇ m − gm qm + 𝜏m + Dm
) ( ( ) ( )
On the other hand, based on Equation (5), we have: + 𝜏h − Mm 𝛼̇ m1 + Zs2T − Cs qs , q̇ s q̇ s − gs qs
( )( ( ) ( ) ) )
q̈ m = Mm−1 qm −Cm qm , q̇ m q̇ m − gm qm + 𝜏m + Dm + 𝜏h + 𝜏s + Ds + 𝜏e − Ms 𝛼̇ s1 + Zm2 T
Cm Zm2 + Zs2T Cs Zs2
( )( ( ) ( ) ) + 𝜃̃mT Γ−1 ̃̇ ̃ T −1 ̃̇ ̃ T −1 ̃̇
m 𝜃 m + 𝜃s Γ s 𝜃 s + 𝜃e Γ e 𝜃 e
q̈ s = Ms−1 qs −Cs qs , q̇ s q̇ s − gs qs + 𝜏s + Ds + 𝜏e (10)
+ 𝛿̃mT 𝜇m−1 𝛿̃̇ m + 𝛿̃sT 𝜇s−1 𝛿̃̇ s (16)
The error dynamics derivatives in the second design step will
be obtained as follows:
On the other hand, Equation (16) can be expressed as:
( )( ( ) ( )
Ż m2 = q̈ m − 𝛼̇ m1 = Mm−1 qm −Cm qm , q̇ m q̇ m gm qm
V̇ 2 = −Zm1
T
Km1 Zm1 − Zs1T Ks1 Zs1 + Zm1
T
Zm2 + Zs1T Zs2
+ 𝜏m + Dm + 𝜏h ) − 𝛼̇ m1 (11) ( ( )
+ Zm2
T
− Cm q̇ mCm Zm2 − gm qm + 𝜏m + Dm + 𝜏h
( )( ( ) ( ) ) ( ( )
Ż s2 = q̈ s − 𝛼̇ s1 = Ms−1 qs −Cs qs , q̇ s q̇ s − gs qs − Mm 𝛼̇ m1 + Zs2T − Cs q̇ s + Cs q̇ s + Cs Zs2 − gs qs + 𝜏s
+ 𝜏s + Ds + 𝜏e ) − 𝛼̇ s1 (12) )
+ Ds + 𝜏e − Ms 𝛼̇ s1 + 𝜃̃mT Γ−1 ̃̇ ̃ T −1 ̃̇ ̃ T −1 ̃̇
m 𝜃 m + 𝜃s Γ s 𝜃 s + 𝜃e Γ e 𝜃 e
The system’s final Lyapunov function will be recommended
as + 𝛿̃mT 𝜇m−1 𝜃̃̇ m + 𝛿̃sT 𝜇s−1 𝜃̃̇ s (17)
1 T 1 1
V2 = V1 + Zm2 Mm Zm2 + Zs2T Ms Zs2 + 𝜃̃mT Γ−1 ̃
m 𝜃m
Considering Equation (4), the −Ci q̇ i + Ci Zi2 = −Ci 𝛼i1
2 2 2 holds. On the other hand, the following expressions will be
1 1 ̃ T −1 ̃ 1 ̃ T −1 ̃ 1 ̃ T −1 ̃ considered as an indefinite part to be approximated by the fuzzy.
+ 𝜃̃sT Γ−1 ̃
s 𝜃s + 𝜃e Γe 𝜃e + 𝛿m 𝜇m 𝛿m + 𝛿s 𝜇s 𝛿s
2 2 2 2
(13) −Cm 𝛼m1 − gm − Mm 𝛼̇ m1 = 𝜙m = 𝜃mT 𝜁m − Cs 𝛼s1 − gs − Ms 𝛼̇ s1
= 𝜙s = 𝜃sT 𝜁s (18)
where 𝜇m , Γe , Γs , Γm and 𝜇s are adjustment parameters. Besides,
we have:
Generally, the if-so fuzzy rules describing output behaviour
𝜃̃i = 𝜃i∗ − 𝜃̂i , (i = m, s ) can be written as:
𝜃̃e = 𝜃e∗ − 𝜃̂e (14) where Ad1 is the fuzzy set, ȳ d is a constant value, and d =
1, … , M is the number of fuzzy rules, x(1, 2, … , n) and y are
where 𝜃̃i , 𝜃̃e and 𝛿̃i are master, slave, environment, and indef- inputs and outputs of the fuzzy system, respectively. Therefore,
inite bound parameters’ approximation errors. Finally, the the output of the fuzzy system can be defined as (20):
derivative of the system’s final Lyapunov function are obtained
as: ∑M [∏ ]
n
d =1 𝜃d w=1 𝜇Adw
(xw )
V̇ 2 = −Zm1
T
Km1 Zm1 − Zs1T Ks1 Zs1 + Zm1
T
Zm2 + Zs1T Zs2 y (xw ) = ∑ [∏ ] = 𝜁 (xw ) 𝜃T (20)
M n
w=1 𝜇Adw
d =1
(xw )
1 T ̇ 1
+ Zm2
T
Mm Ż m2 + Zs2T Ms Ż s2 + Zm2 Mm Zm2 + Zs2T Ṁ s Zs2
2 2 T
where ` = [̄y1 , … , ȳd ] is the vector for output fuzzy member-
+ 𝜃̃mT Γ−1 ̃̇ ̃ T −1 ̃̇ ̃ T −1 ̃̇ ̃ T −1 ̃̇ ship function centers with input membership function of 𝜇Adw .
m 𝜃 m + 𝜃s Γ s 𝜃 s + 𝜃e Γ e 𝜃 e + 𝛿 m 𝜇m 𝛿 m
Besides, vector 𝜁(xw ) can be defined as Equation (21).
+ 𝛿̃sT 𝜇s−1 𝛿̃̇ s (15) ∏n
w=1 𝜇Adw
(xw )
Considering the asymmetric property (Ṁ − 2C = 0) of the 𝜁d (xw ) = ∑ [∏ ] (21)
M n
robot and substituting Equations (11) and (12), the above d =1 𝜇 d (
w=1 Aw w x )
804 BOUTERAA ET AL.
Considering Note 2, the corresponding fuzzy system for The following equation can be obtained:
master and slave robot can be written as:
∑M [∏ ] 𝜙 j − 𝜙̂ j = 𝜙̂ j (X |𝜃∗j ) − 𝜙̂ j + 𝜔 j = 𝜃̃ Tj 𝜁 j + 𝜔 j , ( j = m, s, e)
n
d =1 𝜃 ikd 𝜇 d
w=1 Aiw
( x iw )
𝜙ik (xiw ) = ∑ [∏ ] = 𝜁kT (xiw ) 𝜃ik (22) (32)
M n
w=1 𝜇Ad xiw
d =1
( )
iw Therefore, Equation (30) is rewritten in the following
where = 1, 2, … u, u is the degree of freedom for master and form:
slave robot (i = m, s). then,
V̇ 2 = −Zm1
T
Km1 Zm1 − Zs1T Ks1 Zs1 − Zm2
T
Km2 Zm2 − Zs2T Ks2 Zs2
[ T
𝜙i = 𝜙i1 𝜙i2 ⋯ 𝜙iu ] = 𝜁 T (xiw ) 𝜃i (23) T ̃T
+ Zm2 𝜃m 𝜁m + Zm2
T
(𝜔m + Dm − urm ) + Zs2T 𝜃̃sT 𝜁s
then,
+ Zs2T (𝜔s + 𝜔e + Ds − urm ) + Zs2T 𝜃̃eT 𝜁e + 𝜃̃mT Γ−1 ̃̇
m 𝜃m
[ ]
𝜁T (xiw ) = dig 𝜁1T (xiw )𝜁2T (xiw ) ⋯ 𝜁uT (xiw ) ,
[ + 𝜃̃sT Γ−1 ̃̇ ̃ T −1 ̃̇ ̃ T −1 ̃̇ ̃ T −1 ̃̇
s 𝜃 s + 𝜃e Γe 𝜃 e + 𝛿m 𝜇m 𝛿 m + 𝛿s 𝜇s 𝛿 s (33)
T
𝜃i = 𝜃i1 𝜃i2 … 𝜃iu ] (24)
It is noteworthy that uncertainties bound will be considered
Optimized parameters can be defined as (l = m, s, e): as |𝜔m + Dm |max = 𝛿m∗ and |𝜔s + 𝜔e + Ds |max = 𝛿s∗ , where the
adaptive rules are recommended by this way:
𝜃l∗ = arg min𝜃l 𝜖Ω [supx(lw) 𝜖Rn |𝜁 T (xiw )𝜃l∗ − 𝜁 T (xiw )𝜃l |] (25)
Ω = {𝜃 ∶ m ≤ 𝜃 ≤ m}
−
̄ (26) 𝜃̂̇ s = Γs Zs2T 𝜁s
𝜃̂̇ e = Γe Zs2T 𝜁e
where m and m are design parameters. If the control signals will
(34)
be recommended as follows:
where the environment’s fuzzy regressors are considered as
𝜏m = −Km2 Zm2 − Zm1 − 𝜏h − 𝜙̂ m − urm
𝜁e ∝ (Zs1 , Zs2 , Zm1 , Zm2 ). By employing the adaptive laws (34)
𝜏s = −Ks2 Zs2 − Zs1 − 𝜏̂ e − 𝜙̂ s − urs (27) and control torques (27), all signals of the master side and the
slave of the bilateral teleoperation system are guaranteed to be
such that uri is the robustifying term of the master-slave control bounded and the tracking errors Zm and Zs converge to zero
signals and 𝜙̂ i is an approximation of 𝜙i . The control signal of as t → ∞, which can be small by reducing the values of the
the modeled environment using fuzzy logic can be obtained as fuzzy approximation error and the upper bound of exterior
follows [1–5]: disturbances.
Thus, by substituting Equation (34) into (33) and by
𝜏e = 𝜁 T (xew ) 𝜃e (28) considering the uncertainties bound, we have:
where 𝜃e is the unknown environmental parameter and 𝜏e is the V̇ 2 ≤ −Zm1
T
Km1 Zm1 − Zs1T Ks1 Zs1 − Zm2
T
Km2 Zm2 − Zs2T Ks2 Zs2
output of the fuzzy system. Besides, xew (qs , q̇ s , q̈ s ) is considered
T ̃ T ̂
as the input for the fuzzy system. Because of the inherent delay + Zm2 𝛿m + Zm2 𝛿m − Zm2
T
urm + Zs2T 𝛿̃s + Zs2T 𝛿̂s
within the communication channel, the environment estimator
in the master side can be defined as [1–5]: − Zs2T urs + 𝛿̃mT 𝜇m−1 𝛿̃̇ m + 𝛿̃sT 𝜇s−1 𝛿̃̇ s (35)
V̇ 2 = −Zm1
T
Km1 Zm1 − Zs1T Ks1 Zs1 − Zm2
T
Km2 Zm2 − Zs2T Ks2 Zs2
𝛿̃̇ s = 𝜇s Zs2 → urs = 𝛿̂s Zs2 (36)
+ T
Zm2 (𝜙m − 𝜃̂m ) + Zm2
T
(Dm − urm ) + Zs2T (𝜙s − 𝜃̂s )
then, we will have:
.
+ Zs2T (Ds − urs ) + Zs2T (𝜏e − 𝜏̂ e ) + 𝜃̃mT Γ−1 ̃
m 𝜃m V̇ 2 ≤ −Zm1
T
Km1 Zm1 − Zs1T Ks1 Zs1 − Zm2
T
Km2 Zm2 − Zs2T Ks2 Zs2
. . . .
+ 𝜃̃sT Γ−1 ̃ ̃ T −1 ̃ ̃ T −1 ̃ ̃ T −1 ̃ ≤0 (37)
s 𝜃 s + 𝜃e Γe 𝜃 e + 𝛿m 𝜇m 𝛿 m + 𝛿s 𝜇s 𝛿 s (30)
If we consider fuzzy approximation error as this: Thus, the master-slave tracking error will asymptotically con-
verge to zero and the stability of the control system will be
𝜔i = 𝜙i − 𝜙̂ i (X |𝜃i∗ ) (31) ensured. It is noteworthy that the recommended method does
BOUTERAA ET AL. 805
I = m, s:
[ ]
Mi11(qi ) Mi12(qi )
Mi(qi ) =
Mi21(qi ) Mi22(qi )
[ ]
Ci11(qi ,q̇ i ) Ci12(qi ,q̇ i )
Ci(qi ,q̇ i ) =
Ci21(qi ,q̇ i ) Ci22(qi ,q̇ i )
[ ]
gi1
gi = (38)
gi2
With Mi11(qi ) = li22 mi12 + li12 Mi + 2li1 li2 mi2 ci2 , Mi22(qi ) = li22 mi2 ,
Mi12(qi ) = Mi21(qi ) = li22 mi2 + li1 li2 mi2 ci2 , Ci11(qi ,q̇ i ) = −li1 li2 mi2
FIGURE 2 Two-link robot manipulator
si2 q̇ i2 , Ci12(qi ,q̇ i ) = −li1 li2 mi2 si2 (q̇ i1 + q̇ i2 ), Ci21(qi ,q̇ i ) =
li1 li2 mi2 si2 q̇ i1 , Ci22(qi ,q̇ i ) = 0, gi1 (qi ) = gli2 mi2 ci12 + li1 Mi ci1 ,
gi2 (qi ) = gli2 mi2 ci12 , where mik and lik are the mass and length
not apply passivity condition to prove its stability and is consid- of kth link of the robot respectively (k ∈ {1, 2}). The Mi =
ered as a powerful method in non-linear teleoperation systems. mi1 + mi2 expression is the sum of all robot links, ci12 ′ ci2 ′ ci1 and
The block diagram of the recommended control method is si2 are the abbreviations for cos(qi1 + qi2 ), cos(qi2 ), cos(qi1 ) and
illustrated in Figure 1. sin(qi2 ), respectively. The initially given gravitational accelera-
tion to the controllers has been considered as equal to gi = 9.81
for both master and slave sides. The robot, controller and time
4 SIMULATION RESULTS delay parameters are illustrated in Table 1. The recommended
human-based controller in Figure 3 is applied to the master
In order to evaluate the performance of the proposed approach, robot along y vertical axis. However, because of using torque
two examples of two and three link robot manipulator are equations in dynamical model of the system, the human force
considered. should be transformed into torque forces: 𝜏h = JmT × [0 1] × Fh .
Case 1: Two-link robot manipulator It is noteworthy that any of the fuzzy inputs used in this sim-
For simulation, two robots with two links with round joints ulation is based on the membership functions presented by
(Figure 2) in master and slave sides which have been modelled in Mamdani fuzzy system of Figure 4. The following three mem-
Equation (2) are used [56]. Ci, Mi and gi matrices are as follows bership functions have been considered for the fuzzy system
806 BOUTERAA ET AL.
methods, with the recommended method having a lower track- FIGURE 9 Human-environment torques
ing error rate. Figure 9 explains that the environmental torque
has been applied on the slave robot just in the opposite direction
of human torque; however, tracking (as it has been discussed
earlier) has been done well. These results illustrate the proper
performance of the proposed control approach.
Furthermore, generally in comparison with the reference
method [1], the three following points are drawn:
(a) In the recommended plan, proof of stability can be done
simultaneously for the teleoperation system; while, in reference
[1], the stability analysis is done on a separate basis for both
master and slave robots.
(b) In the recommended approach, the fuzzy approximation
error is considered and the robustifying term is in charge of its
elimination; however, in the reference given, the approximation
error is neglected.
(c) The time delay in this plan is time-variable; however, the
delay has been considered as constant in the reference approach.
Case 2: Three-link robot manipulator
A three-link planar robot manipulator is illustrated in
Figure 10, with link lengths denoted by L1 , L2 and L3 , respec-
tively. The parameters of the planar articulated three-link robot
manipulator considered for simulation are shown in Table 2. FIGURE 10 Three-link robot manipulator with passive second joint [51]
808 BOUTERAA ET AL.
Length of links l1i = 0.3 m, l2i = 0.3 m, l3i = 0.3 m Robustifying parameters 𝜇s = 0.2, 𝜇m = 0.2
Mass of links m1m = 1 kg, m2m = 1 kg, m3m = 1 kg Time-varying values Dm = Ds = [0.1 0.1]
Initial conditions 0 Delay on master robot dm (t ) = 0.2 × sin(t )
Control gains Km1 = 0.02, Km2 = 0.2, Ks1 = 50, Ks2 = 80 Environment parameters Γe = 0.1, Γs = 0.5, Γm = 0.1
Delay on slave robot ds (t ) = 0.2sin(t )
The desired configuration of the system should become: (Ms (qs ) + ΔMs (qs ))̈qs + (Cs (qs , q̇ s ) + ΔCs (qs , q̇ s ))̈qs
q1 = 3 rad, q1 = 1.16 rad, q1 = 6 rad. It can be observed from
Figure 11 that the active joint 1 is capable of reaching its desired
positions in 15 ms. Moreover, the passive joint 2 reaches its + gs (qs ) = 𝜏s + Ds − 𝜏e (40)
desired position in 17 ms, before getting locked as demonstrated
in Figure 12. The active joint 3 reaches the desired position in In order to evaluate the performance of the controller in the
20 ms, as displayed in Figure 13. presence of external disturbances and parametric uncertainties
As can be seen from Figures 13 through 15, by applying vari- in the system, these perturbations have been considered in this
able inputs to the mechanical arm system of the two-axis robot section. The amplitude of the error increases to a maximum of
modelled in this research, the outputs of the system converge to 0.1% following the occurrence of uncertainty in the system but
their permanent value after a number of oscillations. Figure 14 quickly settles to zero due to the ability of the control system
shows the control efforts applied by the controller to stabilize to mitigate system uncertainties and disturbances. The position
the output variables of the mechanical arm system of the robot. error in tracking the reference path in the presence of uncer-
Figure 15 shows the dynamics of the tracking position error in tainties is shown in Figure 15. Note from the dynamics of the
the presence of uncertainties. position error that the proposed controller is quite successful at
The dynamic equations of the robot in the presence of eliminating the external disturbances and parameter uncertain-
parametric uncertainties are written as follows: ties. In order to further validate the performance of the control
approach, we compared the results obtained in this study with
those obtained when implementing the approach proposed in
(Mm (qm ) + ΔMm (qm ))̈qm + (Cm (qm , q̇ m ) + ΔCm (qm , q̇ m ))q̇ m [51]. The position control of the mechanical arm of the three-
link robot using the proposed controller and the controller in
+ gm (qm ) = 𝜏m + Dm + 𝜏h [51] is highlighted in Figure 16.
BOUTERAA ET AL. 809
23. Lu, K., Liu, Z., Wang, Y., Chen, C.P.: Fixed-time adaptive fuzzy control for 40. Hua, C.-C., Liu, X.P.: Delay-dependent stability criteria of teleoperation
uncertain nonlinear systems. IEEE Trans. Fuzzy Syst. 29(12), 3769-3781 systems with asymmetric time-varying delays. IEEE Trans. Rob., 26, (5),
(2020) 925–932 (2010)
24. Capriglione, D., Cerro, G., Ferrigno, L., Milano, F., Moschitta, A.: A multi- 41. Hua, C.-C., Liu, X.P.: A new coordinated slave torque feedback control
frequency approach to mitigate the performance degradation of a magnetic algorithm for network-based teleoperation systems. IEEE/ASME Trans.
positioning system under Cw disturbance conditions. Measurement 161, Mechatron., 18, (2), 764–774 (2012)
107842, (2020) 42. Liu, K., Seuret, A., Xia, Y.: Stability analysis of systems with time-varying
25. Yang, C., Ge, S.S., Xiang, C., Chai, T., Lee, T.H.: Output feedback NN delays via the second-order Bessel–Legendre inequality. Automatica 76,
control for two classes of discrete-time systems with unknown control 138–142 (2017)
directions in a unified approach. IEEE Trans. Neural Networks 19(11), 43. Nemati, A., Peimani, M., Mobayen, S., Sayyedfattahi, S.: Adaptive
1873–1886 (2008) non-singular finite time control of nonlinear disturbed cyber-physical sys-
26. Liu, Y.-J., Wang, W., Tong, S.-C., Liu, Y.-S.: Robust adaptive tracking control tems with actuator cyber-attacks and time-varying delays. Inf. Sci. 612,
for nonlinear systems based on bounds of fuzzy approximation parame- 1111–1126 (2022)
ters. IEEE Trans. Syst. Man Cybern. Part A Syst. Humans 40(1), 170–184 44. Wu, Z.-G., Xu, Y., Lu, R., Wu, Y., Huang, T.: Event-triggered control for
(2009) consensus of multiagent systems with fixed/switching topologies. IEEE
27. Lima, B.N., Balducci, P., Passos, R.P., Novelli, C., Fileni, C.H.P., Vieira, Trans. Syst. Man Cybern. Syst. 48(10), 1736–1746 (2017)
F., de Camargo, L.B., Junior, G.B.V.: Artificial intelligence based on fuzzy 45. Li, Y., Johansson, R., Liu, K., Yin, Y.: Guaranteed cost control design
logic for the analysis of human movement in healthy people: A systematic for delayed teleoperation systems. J. Franklin Inst. 352, (11), 5085–5105
review. Artif. Intell. Rev. 54, 1507–1523 (2021) (2015),
28. Aly, A.A., Vu, M., El-Sousy, F.F., Alotaibi, A., Mousa, G., Le, D.-N., 46. Yang, H., Liu, L., Wang, Y.: Observer-based sliding mode control for bilat-
Mobayen, S.: Fuzzy-based fixed-time nonsingular tracker of exoskeleton eral teleoperation with time-varying delays. Control Eng. Pract. 91, 104097,
robots for disabilities using sliding mode state observer. Mathematics (2019)
10(17), 3147, (2022) 47. Kang, Y., Li, Z., Cao, X., Zhai, D.: Robust control of motion/force for
29. Vu, M.T., Alattas, K.A., Bouteraa, Y., Rahmani, R., Fekih, A., Mobayen, S., robotic manipulators with random time delays. IEEE Trans. Control Syst.
Assawinchaichote, W.: Optimized fuzzy enhanced robust control design Technol. 21, (5), 1708–1718 (2012)
for a Stewart parallel robot. Mathematics 10(11), 1917, (2022) 48. Abidi, K., Yildiz, Y., Korpe, B.: Explicit time-delay compensation in tele-
30. Hwang, C.-L., Chang, N.-W.: Fuzzy decentralized sliding-mode control of operation: An adaptive control approach. Int. J. Robust Nonlinear Control
a car-like mobile robot in distributed sensor-network spaces. IEEE Trans. 26(15), 3388–3403 (2016)
Fuzzy Syst. 16(1), 97–109 (2008) 49. Jing, B., Na, J., Gao, G., Yang, C.: Robust adaptive control for bilat-
31. Santacruz-Reyes, H., Garcia-Valdovinos, L.G., Jiménez-Hernández, H., eral teleoperation systems with guaranteed parameter estimation. In: 2016
Salgado-Jiménez, T., Garcia-Zarco, L.A.: Higher order sliding mode based International Conference on Advanced Robotics and Mechatronics (ICARM). (2016)
impedance control for dual-user bilateral teleoperation under unknown 50. Lu, Z., Huang, P., Liu, Z.: Predictive approach for sensorless bimanual tele-
constant time delay. 2015 IEEE/RSJ International Conference on Intelligent operation under random time delays with adaptive fuzzy control. IEEE
Robots and Systems (IROS), Hamburg, Germany, 5209–5215, (2015) https:// Trans. Ind. Electron. 65, (3), 2439–2448 (2017)
doi.org/10.1109/IROS.2015.7354111 51. Bhave, M., Janardhanan, S., Dewan, L.: Configuration control of pla-
32. Shahdi, A., Sirouspour, S.: Adaptive/robust control for time-delay nar underactuated robotic manipulator using terminal sliding mode.
teleoperation. IEEE Trans. Rob. 25(1), 196–205 (2009) IFAC-PapersOnLine (2016), 49, (1), 148–153
33. Mohammadi, L., Alfi, A., Xu, B.: Robust bilateral control for state con- 52. Li, H., Kawashima, K.: Bilateral teleoperation with delayed force feed-
vergence in uncertain teleoperation systems with time-varying delay: A back using time domain passivity controller. Rob. Comput. Integr. Manuf.
guaranteed cost control design. Nonlinear Dyn. 88(2), 1413–1426 (2017) (2016), 37, 188–196
34. Andaluz, V.H., Quevedo, W.X., Chicaiza, F.A., Varela, J., Gallardo, C., 53. Spong, M.W., Hutchinson, S., Vidyasagar, M.: Robot Modeling and
Sánchez, J.S., Arteaga, O.: Transparency of a bilateral tele-operation scheme Control. (Vol. 3, pp. 75–118). New York: Wiley (2006)
of a mobile manipulator robot. In: International Conference on Augmented 54. Kelly, R., Davila, V.S., Perez, J.A.L.: Control of Robot Manipulators in Joint
Reality, Virtual Reality and Computer Graphics. (2016) Space. Springer Science & Business Media, Berlin/Heidelberg, Germany
35. Tian, D., Yashiro, D., Ohnishi, K.: Haptic transmission by weighting con- (2006)
trol under time-varying communication delay. IET Control Theory Appl. 55. Wang, L.X., Ying, H.: Adaptive fuzzy systems and control: Design and
6(3), 420–429 (2012) stability analysis. J. Intell. Fuzzy Syst. (1995), 3, (2), 187
36. Ghaffari, V., Mobayen, S., Din, S., Rojsiraphisal, T., Vu, M.T.: Robust 56. Slotine, J.-J.E., Li, W.: Applied Nonlinear Control. Prentice Hall, Engle-
tracking composite nonlinear feedback controller design for time-delay wood Cliffs, NJ (1991)
uncertain systems in the presence of input saturation. ISA Trans. 129(B),
88–99 (2022)
37. Alattas, K.A., Mohammadzadeh, A., Mobayen, S., Abo-Dief, H.M.,
Alanazi, A.K., Vu, M.T., Chang, A.: Automatic control for time delay
How to cite this article: Bouteraa, Y., Alattas, K.A.,
markov jump systems under polytopic uncertainties. Mathematics 10(2),
187, (2022) Peng, T., Fekih, A., Rahmani, R., Mobayen, S.: Design of
38. De Rossi, G., Muradore, R.: A bilateral teleoperation architecture using robust adaptive fuzzy control for uncertain bilateral
smith predictor and adaptive network buffering. IFAC-PapersOnLine teleoperation systems based on backstepping approach.
50(1), 11421–11426 (2017) IET Control Theory Appl. 17, 800–813 (2023).
39. Hou, W., Fu, M., Zhang, H.: Distributed consensus of third-order multi-
https://doi.org/10.1049/cth2.12407
agent systems with communication delay. Asian J. Control 20(2), 956–961
(2018)