Adaptive control of time delay teleoperation system with uncertain dynamics

A bilateral adaptive control method based on PEB control structure is designed for a class of time-delay force feedback teleoperation system without external interference and internal friction to study the uncertainty of dynamic parameters and time delay. The stability and tracking performances of the closed-loop constant time delay teleoperation system are analyzed by Lyapunov stability theory. Finally, the controller designed in this paper is successfully applied to the teleoperation system composed of a two-degree of freedom rotating manipulator as the master robot and the slave robot. The simulation is carried out in no operator and environment force or with operator and environment force. The adaptive bilateral control method's control performance is compared with that of the traditional time-delay teleoperation system. Finally, it is verified that the method has good control performance.


Introduction
After decades of research and rapid development of teleoperation robot systems, they have been applied in many fields, such as unmanned submersibles (Sayers and Paul, 1994), space robots (Bejczy, 1994), remote surgery robots (Wright et al., 2006;Su et al., 2020c;Tang et al., 2020), and teleoperation mobile robots (DiMaio et al., 2011;Su et al., 2020a,b,d). In a word, through the research of teleoperation robot systems, human intelligence and robotics can be combined to improve efficiency and reduce cost when performing tasks extensively. Therefore, it has a broad application prospect and rich practical significance.
Usually, in a teleoperation robot system, the operator controls the equipment at the main end to carry out a particular action and then transmits the equipment's action signal at the main end to the controller of the slave end through the communication channel. It then predicts and post-processes the device's control signal at the slave end by calculation (Zhang et al., 2022;Zheng et al., 2022), commands the equipment at the slave end, realizes the human's work skills, and completes the corresponding work tasks. The available remote operation robot system includes a master module, operator module, and the Master Controller, communication channel, Slave Controller, Slave, and .
/fnbot. . environment. Cho and Park (2005) designed an impedance controller for teleoperation systems with communication delay for varying damping. Do and Namerikawa (2009) proposed an impedance control method for the time-delay force feedback teleoperation system based on IOS small gain theory. However, the transparency of the teleoperation system with impedance control structure is not ideal. Position error-based structure (PEB), is a control structure that uses the position signals of the master robot and the slave robot as its own desired signals, and then designs a controller for the generated error signals so as to achieve the goal of system stability and improve tracking performance. When the slave robot contacts the working environment in the teleoperation system, the weighted position error can be fed back to the operator so that the position error presents the interaction force, and then the control torque can be adjusted to reduce the tracking error between the master robot and the slave robot. Nuño et al. (2008) designed a globally stable PD bilateral controller based on the bilateral teleoperation system position error structure. Nuño et al. (2008) proposed an adaptive bilateral control method based on the non-linear teleoperation system's position error structure. Joinié-Maurin et al. (2011) proposed a force feedback teleoperation system based on position error structure. They can compensate for external interference of the system and solve the influence of interference on system stability. Franken et al. (2012) proposed a stable position-based teleoperation system by introducing damping compensation. Kim and Ahn (2013) used a composite adaptive controller to design a teleoperation system based on position differences. Lawrence (1994), Yokokohji and Yoshikawa (1994), and Hastrudi-Zaad and Salcudean (1999) proposed a remote operation system based on four communication channels (4-Channel, 4-Ch). The teleoperation system has been proved to have good force tracking and transparency. Liu and Tavakoli (2011) proposed an adaptive inverse dynamics control method for a non-linear uncertain teleoperation system based on a 4-Ch structure. Dehghan et al. (2016) proposed a non-linear bilateral control method based on a 4-Ch teleoperation system using an adaptive force estimator.
Most of the teleoperation systems studied are based on the position-position type and 4-Ch control structure because these two structures have better characteristics.
The control signals between the master robot and the slave robot communicate through the forward communication channel and the reverse communication channel. Usually, the distance between the system's communication network channels is relatively long, which will cause delays. As should go without saying, the existence of a time delay will lead to 'poor system performance. It will bring the problem of destroying the stability of the teleoperation system. The masterslave robot's dynamic mathematical model is established in the teleoperation system by analyzing its motion characteristics. The structure and mechanical parameters of the master and slave robots are involved in the mathematical model Zhang et al., 2022). In general, most of the teleoperation system's control methods are realized under the condition of the master-slave robot's precise structure and mechanical parameters. However, in the actual teleoperation mechanical system, it is difficult to obtain the robot's precise mechanical parameters, such as mass, length, the center of mass, and moment of inertia. As a result, the system dynamic parameters (inertia vector matrix, centrifugal force matrix, and gravity term matrix) are not accurate, which is common in robot workspace control. In this case, it is impossible to establish an accurate mathematical model of the teleoperation system, which will reduce the system's transparency and performance and even cause the whole system instability. After the above analysis, considering the influence of dynamic parameter uncertainty and time delay, we mainly solve the two problems of time delay and dynamic parameter uncertainty in the teleoperation system.
For the above two problems to be solved in this paper, from the mentioned teleoperation system's research status, we can conclude that there are mainly adaptive control methods to solve communication channel delay and system dynamic parameter uncertainty (Polushin et al., 2010;Haddadi et al., 2015;Wang et al., 2021).
Moreover, most of the existing bilateral adaptive control methods for time-delay teleoperation systems do not consider operator and environmental dynamic parameter uncertainty. They usually choose the adaptive law proposed by Slotine and Li (1991) to estimate the system dynamic parameter vector. However, Slotine and Li (1991) also pointed out that if the reference signal does not meet the persistent excitation condition, the estimated value's accuracy cannot be guaranteed.
Based on the above discussion, in this paper, based on the position error control structure, a new adaptive bilateral teleoperation control method is designed. This new method solves the problems of constant time delay and uncertain dynamic parameters in the teleoperation system. It will ensure the closed-loop constant time-delay system's stability and improve the system's tracking control (Xu et al., 2020;Zhang et al., 2021) performance.

Method and experiments
Mathematical model of teleoperation system The spatial dynamic model of master-slave robot joint Without considering the joint friction and external interference, the general dynamic equations of the master robot and the slave robot of the teleoperation system can be expressed by the following Euler Lagrange formula (Nuño et al., 2008): M q s q s q s + C q s q s ,q s q s + G q s (q s ) = τ s − J T s q s F e (2) Frontiers in Neurorobotics frontiersin.org . /fnbot. .
In order to simplify the description, the subscripts i (i = m, s) are defined to represent the master and slave robots. Then q i ∈ R n×1 ,q i ∈ R n×1 , andq i ∈ R n×1 are divided into the joint angular position, angular velocity, and angular acceleration of the master robot and the slave robot. M q i q i ∈ R n×n is the inertia matrix of the master robot and the slave robot.
C q i q i ,q i ∈ R n×n is the Coriolis force and centripetal force matrix of the master robot and the slave robot. τ i ∈ R n×1 is the joint control torque of the master robot and the slave robot.
G q i q i ∈ R n×1 is the gravity matrix of the slave robot. J i q i ∈ R n×n represents the Jacobian matrix of the master robot and the slave robot. J T i q i is the transposition of the Jacobian matrix. F h ∈ R n×1 represents the force exerted by the operator on the master robot and F e ∈ R n×1 is the interaction force between the slave robot and the environment module.
The dynamic equations of the master robot and slave robot in the teleoperation system, namely formula (1) and formula (2), have the following properties: Property (1): the inertia matrix, M q i q i , is symmetric and positive definite, with maximum and minimum values.
Property (2): the matrix of Coriolis and centrifugal force C q i q i ,q i and satisfies: Equivalently, we can get ζ T 1 2Ṁ q i q i − C q i (q i ,q i ) ζ = 0, ∀ζ ∈ R n× 1 Property (3): We linearly transform the terms on the left of the dynamic formulas (1) and (2) of the master robot and the slave robot of the teleoperation system, and define the unknown constant parameter vector of the robot as θ d = θ d1 , · · ·, θ dr T , then θ d is linear, and: Among them, Y d q i ,q i ,q i ∈ R n×r is called the dynamic regression matrix, which is the known function matrix about the robot's joint variables.

Kinematics model of master-slave robot in workspace
Through geometric analysis of the structure of the robot, we can establish that there is a conversion relationship between the joint position of the robot and the end position of the robot actuator x m , x s ∈ R n×1 can be expressed as Among them, h i q i ∈ R n → R n is a non-linear transformation used to describe the relationship between the master-slave robot actuator's end position and the joint space angle position. Thus, the end velocity of the actuator in the workspace of the master robot and the slave robot is: x m = J m q m q m ,ẋ s = J s q s q s (7) Among them, J m q m , J s q s are the Jacobian matrix of the master robot and the slave robot, respectively. For the sake of simplicity, it is abbreviated as J m , J s . Derivation of formula (7) to time, the acceleration of the executive end of the master-slave robot in the workspace is obtained as follows:ẍ m =J mqm + J mqm (8) x s =J sqs + J sqs (9) Therefore, the above formula's joint space models (1) and (2) can be transformed into the workspace's motion model. It can better describe the robot's contact and the environment and more intuitively relate the robot's end-effector's velocity and acceleration with force acting on the environment's end.

The dynamic model of operator and environment workspace
The operator and environment dynamic models in formulas (1) and (2) of the teleoperation system are usually described by the workspace robot's end actuator position. The expressions of the force exerted by the operator of the master robot's actuator and the interaction between the slave robot and the environment are as follows: F e = f * e + e (x s ,ẋ s ,ẍ s ) Among them, f * h represents the external force exerted by the operator on the main robot, meeting the requirements represents the external force applied by the working environment module to the slave robot.
h (x m ,ẋ m ,ẍ m ) and e (x s ,ẋ s ,ẍ s ) represent the operator's inertial and viscous characteristics and the working environment model, respectively, and can be linear or non-linear passive dynamic functions.
The inertia and viscosity functions in the dynamic model of the operator and environment workspace are represented as h/e (x i ,ẋ i ,ẍ i ) , i = m, s. According to the literature (Malysz and Sirouspour, 2009;Haddadi et al., 2015), h/e (x i ,ẋ i ,ẍ i ) can include the following three situations: Non-linear model:

The formula of non-linear viscous function combined with linear damping:
is a positive definite constant diagonal matrix corresponding to the mass, damping, and elastic coefficient matrices of the operator and the environment.
For different application scenarios, the dynamic models of the operator and environment workspace are different. Because the inertia and viscosity characteristics of the workspace dynamic model of the operator and the environment are very complex, this paper selects the second-order decoupling lineartime-invariant (LTI) model in (3) to approximately establish the dynamic model of the operator and environment workspace in the teleoperation system.
To sum up, we can observe the dynamic model of the operator and environment workspace in the teleoperation system. In practical application, the teleoperation system's dynamic parameters' uncertainty includes master and slave robots' mechanical parameters and the operator and environment models' mechanical parameters, such as mass coefficient, damping coefficient, and elastic coefficient.

Joint space dynamic model of combined teleoperation system
By establishing the master robot module's mathematical models, the slave robot module, the operator module, and the environment module in the teleoperation system, we can observe the master and slave robots' dynamic models in the joint space. The dynamic model of the operator and the environment is in the working space. Thus, the teleoperation system's dynamic model cannot be unified, and the dynamic model is increased. It is challenging to design a bilateral controller. Therefore, it is necessary to use the master robot and slave-robot workspace kinematics model to transform the operator module's dynamic model and environment module workspace into joint space. Finally, the joint space dynamic model of the master robot and slave robot is sorted out. Then the simplified joint space model of the teleoperation system is obtained. By substituting formula (6) to formula (9) into formula (10) and formula (11), the joint space dynamic model of operator and environment is obtained as follows: The two sides of formula (12) and formula (13) are multiplied by J T m and substituting into formula (1) and formula (2), respectively, the simplified joint space model of the teleoperation system is obtained: Among which: After unifying the dynamics of each module in the teleoperation system into the joint space, according to property (1) to property (3), we can deduce the mathematical models of the combined teleoperation system described in formula (14) and formula (15) have the following new properties. for all i = m, s they, respectively represent the master and the slave. Property (4): the inertial matrix M i q i is symmetric and positive definite, with maximum and minimum values.
Property (5): For ∀ξ ∈ R n×1 , the Coriolis matrixṀ i q i and the centrifugal force matrix C i q i ,q i satisfy: Property (6): unify each module's mathematical models in the teleoperation system into the joint space and substitute them into the joint space mathematical models of the master robot and the slave robot. After sorting out the teleoperation system's dynamic models, the items on the left of formula (14) and formula (15) are obtained. The unknown constant parameter vector of the robot can also be defined as θ z = [θ z1 , · · ·, θ zr ] T . After linear transformation, it is obtained that the parameter vector θ z of the robot is linear.
Where i = m, s, Y z q i ,q i ,q i ∈ R n×r is called the dynamic regression matrix, which is the known function matrix about the robot's joint variables.

FIGURE
Control structure diagram of teleoperation robot system based on the position error.

Control structure based on the position error
The control structure of the teleoperation robot system based on position error (Bessa et al., 2010;Polushin et al., 2010) is shown in Figure 1. The control structure connects the master-slave robot's joint position and speed signals through communication channels to form a network. In this network, the master-slave robot's desired position and speed signals are the joint position and speed signals of the slave robot through the reverse channel. The expected position and speed signals of the slave robot are the joint position and speed signals of the master robot through the forward channel, meeting the following requirements: Where q d m (t),q d s (t),q d s (t),q d m (t) , respectively represent the expected position and speed signals of the master-slave robot. q m (t − T m ), q s (t − T s ),q m (t − T m ), andq s (t − T s ) , respectively represent the position and speed signals of the master-slave robot joints passing through the communication channel, T m , T s , respectively represent the delay of the forward communication channel and the delay of the reverse communication channel.
As shown in Figure 1, this structure uses the position signals of the master robot and the slave robot as their expected signals. It then designs the controller for the generated error signals to achieve system stability and improve tracking performance. In the teleoperation system, when the slave robot and the working environment contact each other, the weighted position error can be fed back to the operator so that the interaction force can be presented through the position error. The control torque can then be adjusted to reduce the tracking error between the master robot and the slave robot. Therefore, the teleoperation system based on position error control structure is simple.

The proposed teleoperation system with adaptive control
This paper is devoted to solving constant communication delay and dynamic parameter uncertainty in the teleoperation system. For a teleoperation system without considering the master-slave robot mechanism's internal friction and external interference, the joint space dynamic model can be expressed by formulas (14) and (15). An adaptive control method is designed. We apply the designed bilateral controller to the teleoperation system of the master robot and the slave robot. They are composed of two degrees of freedom and two link manipulators.
In the combined teleoperation system introduced in Section Joint space dynamic model of combined teleoperation system, there is a constant communication delay between the forward channel and the reverse channel. Therefore, the combined teleoperation system's adaptive control structure based on the position error structure is shown in Figure 2. The control structure connects the position and velocity signals of the master-slave robot's end arm through forward and reverses channels to form a network. In this network, the desired position and speed signals of the master robot are the end position and speed signals of the slave robot actuator through the reverse channel. The desired position and speed signals of the slave robot are the end position, and speed signals of the main robot actuator passing through the forward channel, that is: According to the self-adaptive control structure block diagram of the teleoperation system based on the position error structure, the robot's end position's tracking error at the master end and the slave end can be obtained. Thus, the end position's tracking error is obtained through the relationship between the end position of the master robot and the slave robot and the joint angle. That is, the position signal is converted into the joint signal, and the joint . /fnbot. .

FIGURE
Adaptive control block diagram of teleoperation system based on position error structure.
position tracking error of the master robot and the slave robot is obtained as follows: The mathematical models of each module in the teleoperation system are unified into the joint space. They have substituted into the joint space mathematical models of the master robot and the slave robot.
M s q s q s +Ĉ s q s ,q s q s +Ĝ s (q s ) =τ s = Y zs q s ,q s ,q s θ zs Y z q i ,q i ,q i are the dynamic regression function matrix of the known variables of the master and slave robot joints.
It can be seen from formula (28) and formula (29) that the dynamic uncertainty parameter vector θ zm obtained by a linear transformation of the master robot dynamic combination model is composed of the dynamic uncertain parameters of the master robot and the operator. Similarly, the dynamic uncertainty parameter vector obtained from the linear transformation of the robot dynamics combination model consists of the uncertain parameters of the slave robot and the environment dynamics form. In this way, we include the external force f * h of the operator and the external force f * e of the environment into the dynamic uncertainty vector θ zm , θ zs of the system.

Controller design and stability analysis of teleoperation system
The control objective of the teleoperation system based on the position error control structure in this section is to design an adaptive controller at the master end and an adaptive controller at the slave end for the master robot and the slave robot, respectively, to solve the problem that the dynamic parameters of the master robot and the slave robot are uncertain. The communication channel has a constant time delay, ensuring the system's stability and making the system have stable performance. It has good instantaneous characteristics and transparency, thus improving the tracking performance of the system .
Therefore, based on the position error control structure, two-sided adaptive controllers are designed for the master robot and the slave robot, respectively. The stability and position tracking of the teleoperation system's performance is analyzed using high calculus Lyapunov Krasovski functional.

Controller design
Now we start to design non-linear bilateral adaptive controllers. Firstly, to stabilize the motion of the master robot and the slave robot in the teleoperation system, t → ∞, wheṅ q i (t) → 0. The following auxiliary variables need to be defined: Where λ i is the diagonal matrix of positive definite constant and represents the weight of tracking error. Since the acceleration signal cannot be used in a teleoperation system, a regression matrixȲ i q i ,q i , e i ,ė i about tracking error is defined as follows: (33) Formula (33) shows that when e i → 0, there werē Thus, the adaptive controller of the master robot and the adaptive controller of the slave robot is designed as: τ s = −Ȳ s q s ,q s , e s ,ė s θ zs − k s s s − β sės Where, for convenience, here i = m, s, k i is a positive definite constant diagonal matrix, which represents the proportional gain. The termȲ i q i ,q i , e i ,ė i θ i is a model-based term, that is used to compensate for the error.θ i is the estimated value of θ i . The proportion term for auxiliary variables s i is the additional term of error compensation. The other term β iėi is used to eliminate the additional energy dissipation caused by a time delay, where β i is the diagonal matrix of positive definite constant, representing the dissipation coefficient. By substituting formula (34) and formula (35) into formula (14) and formula (15), the following results are obtained: Among them, θ zm , θ zs , respectively represents the estimated error of θ zm , θ zs , which satisfies θ zi = θ zi −θ zi . Usually, the adaptive law proposed by Slotine and Li (1991) Among them, i = m, s, Ŵ i is the positive definite diagonal constant matrix gain. The above estimation method has been widely used in adaptive control applications. However, Slotine and Li (1991) also pointed out that if the reference signal does not meet the persistent excitation conditions, the above estimation accuracy cannot be guaranteed. Therefore, to improve the adaptive estimation value's accuracy and convergence rate, the above formula's estimation method is improved using the input torque estimation error. The adaptive estimation law of uncertain parameters θ zi is defined as follows: Among them, for the sake of convenience i = m, s, ε i = Y zi θ zi − Y ziθzi = Y zi θ zi denotes the estimation error of calculating input torque used for correcting the estimated value. Here θ i represents the estimation error. Ŵ i (t) = diag γ 1 (t) , · · ·, γ p (t) is the adaptive gain matrix, satisfying where γ n (t) = a n e − t 0 f n (υ)dυ + b n . a n . b n are the constants, and f n (υ) ≥ 0, n = 1, · · ·, p . p i (t) is the low-pass filtered signal of Y T zi Y zi which satisfies the formula where µ and η are constant and ξ 1 and ξ 2 are constant too. According to the robust control technique in Nuño et al. (2010), the coefficients ξ 1 can be defined as: where α, δ is the constant. λ Y T zi Y zi min is the minimum eigenvalue of Y T zi Y zi .

Analysis of system stability and tracking performance
In this paper, the stability of closed-loop systems with time delay is described in formulas (14) and (15), adaptive bilateral control laws in formula (36) and formula (37), and adaptive laws shown in formula (39) are analyzed. The joint position tracking errorė m ,ė s and the estimation error of the adaptive law θ zm , θ zs of the teleoperation system are bounded. The velocityq m ,q s and tracking errorė m ,ė s converge to zero.: For a teleoperation robot system, a high calculus Lyapunov Krasovski candidate function V (t) is defined: From property (4), we know that M m q m , M s q s are symmetric positive definite matrices, Ŵ i (t), k i , β i are positive definite matrices. Therefore, V (t) is also positive definite matrices.
By deriving the two sides of the candidate Lyapunov function formula (41), the following results are obtained: According to property (5): Because θ zi = θ zi −θ zi , so there are: By substituting formulas (43) and (44) into formula (42), the results are as follows: (45) From the previous definition, Ŵ i (t) > 0 andŴ i (t) ≤ 0 are negative semidefinite, theṅ Where, B i , k i , β i is the positive definite diagonal constant, Therefore, V (t) ≥ 0 whileV (t) ≤ 0. The auxiliary variables s i of the system, the tracking error e i , and the estimation error θ zi of the adaptive law are bounded. Then we use the theorem of Barbalat to know that theV (t) asymptotic approach to 0, then at that time t → ∞, there are s i → 0,ė i → 0.
Thus, lim t→∞ (s m + s s ) = 0, define r (t) = x m +x s . We can get the following conclusions: It can be inferred that lim t→∞ s 2 r (s) = 0. that isẋ m +ẋ s → 0 when t → ∞.
Similarly, define r (t) = x m − x s , getẋ m −ẋ s → 0 as t → ∞. Therefore, it is proved that the velocity of the master robot and the slave robot approaches zero asymptotically, that isẋ m ,ẋ s → 0, as t → ∞. Considering the auxiliary variables s i =ẋ i + λ i e i → 0, the position tracking error converges to 0.

Simulation verification
The master robot and the slave robot are simulated, respectively, in free motion and in contact with the operator and the environment to verify the above control algorithm's effectiveness. In these two cases, the position tracking and force tracking of the traditional adaptive control algorithm in Nuño et al. (2010) and the designed adaptive control method are compared.
In this paper's simulation experiment, the teleoperation system structure is selected, as shown in Figure 3. The master and slave robots of the teleoperation system with time-delay force feedback are selected as the two degrees of freedom two-link rotary manipulator robots. Its dynamic mathematical model is formula (48) ∼ formula (52). As shown in Figure 3, the kinematic Jacobian matrices of the master and slave robots of the closed-loop delay teleoperation system are: Frontiers in Neurorobotics frontiersin.org . /fnbot. .

FIGURE
Schematic diagram of teleoperation system selected by simulation.

Simulation of free motion
Firstly, in the teleoperation system based on the position error structure shown in Figure 2, the effectiveness of the adaptive control method proposed in this chapter is verified by simulation in the case of the free motion of the master robot and the slave robot. There is no contact force in the process of the

Name Symbol Value Unit
Mass of rod 1 m 1 1 kg Length of rod 1 l 1 1 m Distance from the centroid of rod 1 to joint 0 l c1 1/2 m Moment of inertia of rod 1 I 1 1/12 kg × m 2 Mass of rod 2 m 2 3 kg Length of rod 2 l 2 2 m Distance from the centroid of rod 2 to joint 1 l c2 1 m Moment of inertia of rod 2 I 2 2/5 kg × m 2 master robot and the robot, so the expression formula is F h = 0 and F e = 0. The simulation experiment is carried out by using the MATLAB Simulink tool. The experiment proves that the bilateral teleoperation system with the dynamic model of formula (14) and formula (15) is stable under the control of adaptive controller type formula (34) and formula (35) and adaptive law formula (39). The position tracking error is bounded and can converge to near 0.
In the experiment, it is assumed that the dynamic parameters of the master robot and the slave robot are unknown. Let the dynamic parameter vectors of the master robot and slave robot be: Here = m 1 l 2 c1 + m 2 l 2 1 + I 1 + m 2 l 2 c2 + I 2 , ν = m 2 l 2 c2 + I 2 , and κ = m c l 1 l c2 . The initial parameters of the master robot and the slave robot of the teleoperation system are set as initial position parameters x m (0) = [0.184, 2.6389] and x s (0) = [3, 0]. The initial dynamic uncertainty parameter vector is θ m (0) = θ s (0) = [4.1, 1.9, 1.7] T . The gain parameters of the master controller and slave controller are set k m = diag (30, 30) and k s = diag (30, 30), dissipation coefficient β m = β s = diag (0.8, 0.6), and λ m = λ s = diag (5, 5). The adaptive gain matrix of master adaptive law and slave adaptive law is as follows: And the constant parameters of the adaptive law and the slave adaptive law are set to α m = α s = 1, δ m = δ s = 1, ξ 2m = ξ 2s = 0.8, µ m = µ s = 1, and η m = η s = 1. Acceleration of gravity g = 9.81. In the teleoperation system, the forward channel delay and the reverse channel delay are T m = T s = 0.5 s. The simulation results are shown in Figures 4, 5.
The traditional adaptive two-sided control method of the time-delay teleoperation system in Nuño et al. (2010) .
/fnbot. .  to illustrate the advantages of the two-sided controller are compared and analyzed. The control law is as follows: The adaptive law is as follows: Where r i =ė i + α i e i , α i , K i , L id are constant diagonal matrix.
In the simulation experiment, after repeated debugging, we can analyze and set the parameters of bilateral controller and adaptive law as α m = α s = 0.25, K m = K s = diag (100, 100), and L md = L sd = diag (40, 40).

Simulation with contact
In the case of contact motion, the operator contacts the master robot. From the contact between the robot and the environment, the effectiveness of the adaptive bilateral controller scheme designed in this chapter for the time-delay force feedback teleoperation system is verified.
This section also uses MATLAB SIMULINK for simulation verification. In the simulation teleoperation system and the environment space model is set. The environmental dynamic parameters are set as m e = 0.1 kg, where the mass coefficient In the experiment, it is considered that the dynamic parameters of the master robot and the slave robot are unknown, and the bilateral adaptive controller can solve this problem. According to properties (6), the dynamic models of the master robot and the slave robot are transformed linearly. The dynamic uncertain parameter vectors are obtained as follows: θ s = l 2 c2 m 2 + l 2 2 M e , l 1 l c2 m 2 , l 1 l 2 m e , l 1 l 2 b e , l 2 2 b e , l 1 l 2 k e , l 2 1 m 2 + I 1 + I 2 + l 2 c1 m 1 + l 2 1 m e , l 2 1 b e T By introducing the parameters in Table 1 into formula (59) and formula (60), the actual dynamic parameter vector can be calculated as follows: 3.8, 3, 0.4, 100, 200, 200, 3.933, 50, 25 (1 − cos (πt)) , 50 (1 − cos (πt))] T θ s = [3.4, 3, 0.2, 100, 200, 200, 3.833, 50] Table 1. Acceleration of gravity g = 9.81. Set the initial dynamic parameter vector as follows: 1.9, 1.7, 0.5, 90, 190, 180, 2.0, 40, 20, 45] T θ s (0) = [1. 9, 1.7, 0.1, 90, 190, 190, 1.9, 40]  After repeated debugging, the design parameters of the master controller and slave controller formula (33) and formula (34) are set as follows k m = diag (50, 50) and k s = diag (50, 50). gain parameter is β m = β s = diag (0.8, 0.6) and λ m = λ s = diag (5, 5). Dissipation coefficient is β m = β s = diag (0.8, 0.6). The adaptive law of the master robot and the slave robot is formula (3-14). After repeated debugging, the adaptive gain matrix of the master robot and the slave robot is: Ŵ m (t) = e −t + 2 I 8×8 and Ŵ s (t) = e −t + 2 I 6×6 . The constant parameters of the adaptive law and the slave adaptive law are set to α m = α s = 1, δ m = δ s = 1, ξ 2m = ξ 2s = 0.8, µ m = µ s = 1, and η m = η s = 1. In the teleoperation system, the forward channel delay and the reverse channel delay are T m = T s = 0.5s. The experimental results are shown in Figure 7. The adaptive two-sided controller described in formula (57) and adaptive law described in formula (58) in Nuño et al. (2010) are compared under the same environmental conditions to illustrate the advantages of the two-sided controller. The parameters of the controller in Nuño et al. (2010) are set to α m = α s = 5, K m = K s = 200I, L md = L sd = 40I.

Results and findings Simulation results of free motion
This section verifies the proposed adaptive control method's effectiveness while the master robot and slave robot are in free motion. The proposed method is compared with the method mentioned in Nuño et al. (2010), which could not be fully discussed and explained in this paper. In the simulation experiment, when the teleoperation system moves freely, the master robot and the slave robot are not affected by external forces. Their position tracking trajectory is shown in Figure 4. We can tell from Figure 4 that the slave robot of the teleoperation system can track the position of the master robot within 1.8 s and ensure the stability of the closed-loop system with time-delay force feedback.
The simulation results of the comparative experiment are shown in Figures 5, 6. In the figures, ACIT (adaptive bilateral control with improved tracking performance) represents the bilateral adaptive control method designed. TAC (traditional adaptive bilateral controller) represents the traditional bilateral adaptive control method in Nuño et al. (2010). Figure 5 shows the end position tracking error curves of the master robot and the slave robot of the teleoperation system with time-delay force feedback under the two control methods. From this graph, we can observe that the time of position tracking error convergence to zero in the two-sided adaptive control method is about 1.8 s. In comparison, the traditional adaptive control method is about 2.8 s. Figures 6A,B show the estimated values of the dynamic uncertain parameter vector of the slave robot of the teleoperation system under the traditional bilateral adaptive control method and the bilateral adaptive control method designed in this chapter, respectively. For the time-delay teleoperation system with uncertain dynamic parameters, the adaptive bilateral control method designed in this chapter can quickly converge to the real value in 0.5 s compared with the traditional two-sided adaptive control method in Nuño et al. (2010). However, in the traditional adaptive control method, it takes 1.8 s to converge to the region near the real value, which is unstable and fluctuates greatly. Its accuracy is relatively low.

Simulation results with contact
In the case of contact motion, the operator contacts the master robot. From the contact between the robot and the environment, the effectiveness of the adaptive bilateral controller scheme designed in this chapter for the time-delay force feedback teleoperation system is verified. The experimental results are shown in Figure 7. Figure 7A shows the position tracking the trajectory of the mechanical end of the master robot and the slave robot of the teleoperation system. We can observe that the position curves of the master robot and the slave robot almost coincide in about 3.5 s. In Figure 7B, The input torques of joint 1 and joint 2 of the master robot and slave robot are displayed. The slave robot of the teleoperation system can track the   position of the upper master robot under the contact condition, and the time-delay force feedback teleoperation system can maintain stability.
The comparison of experimental results shows that the estimation error of uncertain vector of dynamic parameters of the teleoperation system with time-delay force feedback under the two control methods and the tracking error curve of the master robot and the slave robot end position is shown in Figures 8, 9. In the figure, TAC is used to represent the method in Nuño et al. (2010). ACIT represents the method designed in this paper. Figure 8A shows the estimation error of the uncertainty vector θ s of the robot dynamic parameters for the teleoperation system under the traditional adaptive control method in Nuño et al. (2010). We can observe that the position tracking convergence speed of the teleoperation system under the traditional adaptive control method is slow, and the tracking error accuracy is not high. At the same time, its adaptive estimation converges slowly and fluctuates greatly, which will affect the control performance of the system in practical application. Figure 8B shows the estimated error of the teleoperation system from the uncertain vector θ s of robot dynamic parameters under the control method designed in this paper. The estimated vector value from the uncertainty of robot dynamic parameters can quickly converge to the true value in 2 s. Figure 9 compares the traditional adaptive control method and the control method designed in this chapter. We can observe that the adaptive control method in this paper has better tracking performance.

Conclusions
Aiming at the uncertainty of dynamic parameters and time delay in the system, a bilateral adaptive control method based on PEB control structure is designed for a class of time-delay force feedback teleoperation systems without external interference and internal friction. The stability and tracking performance of the closed-loop constant time delay teleoperation system is analyzed by Lyapunov stability theory. Finally, the controller designed in this paper is successfully applied to the teleoperation system composed of a two-degree freedom rotating manipulator as a master robot and a slave robot. The advantages of this method are verified by simulation under the conditions of free motion and environmental forces and with contact and environmental forces, respectively. The control performance of this method is compared with that of the traditional adaptive bilateral control method for the timedelay teleoperation system, and the effectiveness of this method is analyzed by experiments. Finally, the method is proved to have good control performance.
The bilateral controller designed in this paper has the following advantages: (1) Based on the position error structure, the operator and environment's dynamic parameters can be included in the unknown vector of the whole system for estimation to avoid the influence of uncertain parameters in the operator and environment model on the system. (2) The estimation error is used to compensate for the adaptive estimation law to improve the system's tracking performance to improve the adaptive estimation accuracy and convergence speed.
(3) By compensating for the dissipation caused by a time delay, the controller's performance is improved, thus ensuring the system's stability and improving the system's tracking performance.
There are still many possible improvements and experiments that could be performed based on this paper. For example, for more complex cases, such as non-linear models and timevarying cases, the control of teleoperation systems needs to be further studied. Besides, This paper only analyzes and verifies the control method of the time-delay force feedback teleoperation system in this article from two aspects of theory and simulation. There may be some problems that have not been considered in the actual system. Therefore, in future research work, the above method should be applied to the real teleoperation system for verification and improvement, so as to make it more practical.