I Introduction
Generally speaking, exoskeletons are humanrobot interfaces (HRIs) that ease interaction and communication between human operators and robots. In the context of teleoperation systems, a human operator needs to have an adequate sense of the remote site and control of the slave robot to accomplish a successful teleoperation task [1]. As a reasonable solution, an exoskeleton transfers human body movements to the slave robot and provides a sense of the remote site through haptic technology, allowing the operator to execute dexterous tasks.
The aim of our project is to study two upperlimb 7DoF exoskeletons as master robots and two heavyduty hydraulic mobile manipulators as slave robots within the concept of multimaster multislave (MMMS) teleoperation systems [2]. Like in our earlier work [3], we considered realworld dissimilar master and slave robots, where the range of motion and force amplification was 1:4 for motion and 1:800 for force amplification. This required us to consider the problem of motion/force scaling which resulted complicated stability analysis tackled in [3]. However, in [3], a commercial 6DoF desktop haptic manipulator was utilized for controlling a stationary heavyduty hydraulic slave robot, whereas our current target is to utilize wearable 7DoF exoskeletons to control multiple mobile slave robots in humanlike object manipulation scenarios. Therefore, it is essential to design a control strategy with the capability of handling these highly complex MMMS systems for the required high performance control objectives, such as high transparency as in [3]. For this ambitious objective, in this paper, decentralized nonlinear control problem of the 7DoF righthand upperlimb exoskeleton is examined. In the future, this work will be expanded to the control of MMMS teleoperation systems.
The control problem of onearmed upperlimb exoskeletons has been widely examined. In [4][5], adaptive integral terminal sliding mode control (SMC) and integral secondorder terminal SMC have been employed to control upperlimb exoskeletons, respectively. In [6], an integral secondorder terminal SMC with quasitime delay estimation (QTDE) was employed for force and position control of an upperlimb exoskeleton. Adaptive impedance control with a nonlinear disturbance observer with estimation of humandesired intended motion has been considered in [7] and [8]. Other control schemes, such as backstepping [9], backsteppingSMC [10], and intelligent approaches like fuzzy SMC [11] have also been applied to the exoskeleton control problem. Although the abovementioned control schemes have shown good performance, they are not perfect choices for MMMS systems, especially when MMMS teleoperation systems work to manipulate a common object. In such a scenario, the levels of nonlinearity and complexity of the system are intensely high, and the aforementioned approaches can face difficulties in handling them. For example, backstepping can end up with an explosion of complexity, even for a simple system [12]. Moreover, most of the mentioned works are about the applications of exoskeletons in rehabilitation with one arm, while our goal is to use exoskeletons in MMMS teleoperation systems. However, MMMS systems encompass rehabilitation and telerehabilitation.
The virtual decomposition control (VDC) approach is a decentralized control scheme that breaks the entire complex system into multiple subsystems that are easy to deal with and design control action [13]. Virtual cutting points (VCPs) are points at which a system is broken down into subsystems. With the introduction of virtual power flows (VPFs) at VCPs as stability connectors, the stability of each subsystem, which is examined distinctly at the subsystem level, is expanded to the whole system and culminates in asymptotic convergence. Such an approach of modulation makes VDC a much more stable control algorithm in a way that altering the control action of one subsystem does not affect the others. Therefore, VDC is able to handle extremely high DoF dissimilar teleoperation systems and perfect choice for our future aims.
So far, the VDC method has been utilized in fields such as hydraulic manipulators [14][16], mobile manipulators [17], and underwater applications [18]. Nonetheless, for the exoskeleton application, only [19]
has employed an adaptive VDC scheme to control a 7DoF exoskeleton with estimation of the joint’s moment of inertia and friction. In contrast to
[19], which examined rehabilitation without motion/force scaling, our aim is to use the exoskeleton in dissimilar MMMS teleoperation systems with high motion/force scaling. However, because adaptive VDC utilizes a projection function [13] to estimate unknown parameters of the rigid body, it has two major drawbacks: ) it requires lower and upper bounds of each parameter, and ) it designates an adaptation gain to each parameter. On the other hand, the inertial parameter vector of the VDCbased dynamic equation encompasses 13 parameters, which in general, the unique inertial parameter vector must contain 10 parameters [20]. This means that 13 adaptation gains and 26 lower and upper bounds are required for adaptive VDC, where is the number of rigid bodies. Therefore, as the DoF of a system increases, like in [19] and MMMS teleoperation systems, preparing and tuning all these parameters can make VDC timeconsuming, despite its benefit of being decentralized.The aim of this paper is to eliminate the abovementioned problems by proposing a novel adaptive VDC approach that utilizes NAL [21] for the estimation of unknown rigid body parameters. The NAL function is defined based on the Riemannian metric and is coordinate invariant with respect to coordinate frames or physical unit selection. This adaptation law only requires one tuning gain for the entire chain of rigid bodies with no need for upper and lower bounds. The contributions of this paper are as follows:

Redesigning dynamic equations of VDC into the 10form inertial parameter vectors.

Proposing an asymptoticstability guaranteed, novel VDCNAL approach. By doing so, all the 26 lower and upper bounds are eliminated and 13 adaptation gains are reduced to 1, while the physical consistency of the estimated parameters is ensured.

Applying the presented approach to control a 7DoF exoskeleton. The commercial exoskeleton ABLE by Haption is used in the experiment to examine the performance of the modified controller.
Ii Proposed VDC Structure
In this section, the prevailing form of VDC is explained and essential theorems and definitions are provided to prove stability of the system in the sense of virtual stability. Then, the dynamic equations of the system are reformulated to suit the presented adaptation law.
Iia Virtual Decomposition Control Approach
Consider {A} as a frame that is attached to a rigid body. Then, the 6D linear/angular velocity vector and force/moment vector can be expressed as follows [13]:
(1) 
where and are the linear and angular velocities of frame {A}, and and are the force and moment expressed in frame {A}, respectively. The transformation matrix that transforms force/moment vectors and velocity vectors between frame {A} and {B}, where {B} is another frame attached to the rigid body, is [13]
(2) 
where is a rotation matrix between frame {A} and {B}, and denotes a vector from the origin of frame {A} to the origin of frame {B}, expressed in {A}. Based on the , the force/moment and velocity vectors can be transformed between frames as [13]
(3) 
Then, the dynamic equation of the free rigid body expressed in frame {A} can be derived as follows [13]:
(4) 
where is the mass matrix, is the centrifugal and coriolis matrix, is the gravity vector, and is the net force/moment vector applied to the rigid body. These matrices are [13]
(5) 
where denotes the mass of the rigid body, is a 3
where denotes the moment of inertia matrix around the center of mass, and is the gravity vector.The dynamic equation (4) can be rewritten as a linearinparameter expression:
(6) 
where is the regressor matrix generated from , is the required velocity vector, and is the rigid body inertial parameter vector. Furthermore, the required net force/moment vector can be designed as
(7) 
where is the estimation of the unknown inertial parameter vector, and is the symmetric positivedefinite matrix. To estimate the inertial parameter vector of each rigid body, the projection function is utilized as follows [13]:
(8) 
with
(9) 
and
(10) 
where is the adaptation gain, and and are the lower and upper limits of the parameter, respectively.
Each adjacent subsystem interacts with each other at VCPs with governing dynamics of (4) and a control law of (7). Such interactions are designated as VPFs. Considering the fixed frame {A}, VPF can be defined as follows [13]:
(11) 
where is the required force. VPFs act as stability connectors in the concept of virtual stability and extend the stability of each subsystem to the stability of entire system. The following lemma expresses the virtual stability concept of VDC:
Definition 1.[13] A single subsystem that is virtually decomposed from complex system with dynamics (4), control law (7) with known , and affiliated function can be said to be virtually stable if and only if there exists a nonnegative accompanying function
(12) 
in a way that
(13) 
where and are two blockdiagonal positivedefinite matrices, and and are two adjacent neighbors of A.
Theorem 1.[13] Consider a complex system that is virtually decomposed into subsystems. If all the decomposed subsystems are virtually stable in the sense of definition 1, then the entire system is stable.
As can be seen from (8) and (9), 13 adaptation gains must be tuned and 13 lowerbound and 13 upperbound parameters must be designated for a single subsystem, which for a system with rigid bodies, they become 13 and 26. The aim of this paper is to reduce the amount of required parameters by replacing the projection function with the NAL adaptation function. This replacement eliminates all the aforementioned burdens by reducing 13 gains to 1 and eliminating 26 inertial parameter bounds. However, as can be seen from (6), contains 13 parameters, whereas NAL requires to be . Consequently, the next step is to reformulate the prevailing equations.
IiB Redesigned Dynamic Equations
In this part, a new structure for the VDC scheme is proposed. This structure equips VDC with 10 unique representations of the rigid body in space with respect to frame {A} and puts it in a natural form, making it compatible with NAL. By defining , we can rewrite (5) as
(14) 
Details are presented in Appendix A. It must be mentioned that gravity vector does not need to be reformulated. Consequently, the linearinparameter form of new matrices can be written as
(15) 
where is the new regressor matrix and is the 10form inertial parameter vector. Details can be found in Appendix B. It is worth mentioning that the reqressor matrix represented in (B.2) has a more compact and simple form than the one in [13]. Now, by using (15), it is possible to merge VDC and NAL.
Iii Natural Adaptation Law
In the previous section, we derived a 10form inertial parameter vector for the VDC context. This vector, (detailed in Appendix B), is a unique representation of a rigid body in space that contains mass and inertia represented in the frame {A}. However, to have such a unique representation just a subset of can be assigned to . Such a subset must satisfy physical consistency conditions to be a unique representation of the rigid body.
Proposition 1.[22] Consider a nonnegative mass density function , such that
(16) 
where is a symmetric matrix, is a second moment matrix, and is the space of realsymmetric matrices. The rigid body inertial parameter vector is physically consistent if there exists a onetoone linear map ,
(17) 
(18) 
such that is positivedefinite (i.e., ). is the trace operator and .
Based on proposition 1, the set of physically consistent inertial parameter vectors for a rigid body can be defined on manifold as
(19) 
where is the space of all realsymmetric, positivedefinite matrices.
Iiia Riemannian Geometry of and Bregman Divergence
is the tangent space for a given . Therefore, Riemannian metric invariant under group action , where is any nonsingular matrix, for given and can be defined as
(20) 
Then, the geodesic distance between two arbitrary points can be written in the sense of (20)
(21) 
where
are the eigenvalues of
. Thereafter, a distance metric on can be defined as(22) 
One other possible way to define distance metric is Bregman divergence. Let be a continuousdifferentiable, strictly convex function defined on a closed convex set . The Bregman distance is defined as
(23) 
By considering and assigning logdet function to (i.e., for ), the resulting Bregman divergence is,
(24) 
where are the eigenvalues of . Now, by utilizing a onetoone mapping from to , a distance metric on can be expressed as
(25) 
IiiB Extracted NAL
In this part, the NAL function is embedded into VDC and proves that asymptotic stability is ensured. Consider a nonnegative accompanying function
(26) 
where is in the sense of (12), and is the accompanying function of the adaptive term. By taking the derivative of (26) and utilizing (15) and the control law in the sense of (7), the following is achieved:
(27) 
where and is defined in (10). Most importantly, must possess the feature in the form of
(28) 
The common quadratic function, which results in an adaptation function like the projection function, embodies such a feature. Another possible candidate for can be geodesic distance in the form of (21). However, due to the nonlinearity of (21), its time derivative can not be written in the form of (28). Therefore, stability of the overall system cannot be guaranteed. Another choice is to consider as the Bregman divergence distance metric[21]:
(29) 
where . From (24) and according to the fact that , it can be concluded that (29) is a valid accompanying function. Taking the derivative of (29) leads to
(30) 
where .
Proposition 2. Suppose that the adaptive VDC control law in the sense of (7) results in a time derivative of the accompanying function in the form of (27). Considering as defined in (29) and its time derivative in (30), the NAL is defined as[21]
(31) 
and results in asymptotic stability of the system. is a unique symmetric matrix that satisfies .
Proof. Substituting (30) in (27) and using (31) leads to
which is in the sense of definition 1. Therefore, by using theorem 1, asymptotic stability of the whole system is established.
As demonstrated in (31), there is only one tuning gain for the estimation of rigid bodies in the absence of lower and upper limits for the inertial parameters vector. In addition, the NAL function ensures that all the estimated parameters will remain physically consistent by ascertaining that will remain positivedefinite[21].
Iv 7 DOF Exoskeleton Control
In this section, the modified VDC is employed to control the 7DoF commercial ABLE exoskeleton. Assigned VDC frames are demonstrated in Fig. 1. Using these frames, the entire robot configuration is decomposed into subsystems of joints and rigid bodies. It is assumed that joints have negligible mass.
Iva Rigid Body Control Equations
The velocity vector can be computed as
(32) 
(33) 
where , is the velocity vector of the ground, and is the velocity vector of the endeffector, and are the joint angles and velocities, respectively, for , for , and for , which , , and . Moreover, can be computed with (2) by replacing and with and , respectively. The required velocities also can be expressed as
(34) 
(35) 
where is the required velocity of the ground, and is the required velocity vector of the endeffector, and can be defined as follows:
(36) 
where is the Jacobian matrix, and are the desired velocity and pose of the endeffector in the Cartesian space, respectively, is the diagonal, positivedefinite matrix, and is the computed pose of the endeffector. Here, the defined endeffector’s pose error is
(37) 
where is the common position error with , and is the quaternionbased orientation error, defined as
(38) 
where and are unit quaternions computed from the rotation matrix, and are desired unit quaternions computed from the desired Euler angles of the endeffector with angles related to convention, and is a matrix operator. Now, the net force/moment vectors can be computed for each rigid body by replacing with in the (4) for . Then, the force/moment vectors can be expressed as
(39) 
(40) 
with required terms as
(41) 
(42) 
where is the interaction force between the endeffector and environment, and is the applied force from joint 1 to the ground. is the required net force/moment vector that is defined in (7).
IvB Virtual stability
In this section, the virtual stability of the rigid body dynamics (15) under the control action of (43) is established.
Lemma 1. If the accompanying function for links with the control of (43) and the adaptation law of (31) is chosen as
(44) 
it results in
(45) 
which shows the virtual stability of the link’s dynamics in the sense of definition 1.
Proof. Proof can be established by means of proposition 2 and [13].
IvC Joint Control Equations
The dynamic equations of joints in the absence of friction can be expressed as
(46) 
for , are the moment of inertia of the joint, denotes the net torque applied to the joint, is the control torque of the joint, and . The required net torque for each joint can be defined as
(47) 
where , , are estimations of inertial parameter vectors for each motor with the update function of
(48) 
where is corresponding pseudo inertia matrix. Therefore, the control torque can be designed as
(49) 
where .
IvD Virtual stability
Lemma 2. If the accompanying function for the joint with control action of (49) and adaptation law of (48) is chosen as
(50) 
it results in
(51) 
which is virtually stable in the sense of definition 1.
Proof. Proof is the same as lemma 1.
Theorem 2. If the accompanying function of the entire system is considered as
(52) 
it follows,
(53) 
This shows that the entire system is asymptotically stable in the sense of theorem 1.
Proof. The procedure is the same as Appendix C of [23].
V Experimental Results
In this section, the experimental results are presented to exhibit the performance of the new adaptive VDC method. The commercial upperlimb exoskeleton ABLE is utilized for the experiments. The haptic technology that this robot contains makes it a perfect choice for bilateral teleoperation, which is the goal of our future work. All highlevel processes are accomplished in a host computer, and computed torque signals are transmitted to the robot using SIMULINK and the Virtuose interface, provided by Haption, with a sample time of 1 ms. Experiments are performed through the following scenario, as shown in Fig. 2. First, the desired Cartesian space targets are given to the order trajectory generator [24], which generates smooth acceleration trajectories in Cartesian space. Then, the required joint velocities are derived by employing the ClosedLoop Inverse Kinematic (CLICK) algorithm [25], which is an online inverse kinematic solver. The computed required joint velocities are then fed to the proposed adaptive VDC method to compute the required net torques to execute Cartesianspace tasks.
To show that the new adaptive VDC has the same performance as the original one, simulation results related to the joint control of a 3DoF robot utilizing both the original and proposed adaptive VDC are demonstrated in Fig. 3, where the solid lines are the tracking error of each joint with the projection function and the dashed lines correspond to the NAL function. The control objective is to track the desire trajectory for each joint, starting from the zero initial condition. All the adaptation and control gains are considered the same. As can be seen, the performance of the proposed adaptive VDC is as good as the original one, and the tracking error is almost zero. The difference is that the NAL function requires only one gain for tuning, but the original adaptive VDC requires 39 adaptation gains and 78 upper and lower limits. This shows that applying the proposed adaptive VDC to highDoF systems can be highly timeefficient.
In the following part, the experimental results of the pose control of the 7DoF redundant exoskeleton are presented. It must be mentioned that to activate joint motors, the human operator must grasp the handle of the exoskeleton, which, as can be seen in Fig. 1, is the endeffector. Therefore, the intention force of the operator, which can be considered an unknown timevarying disturbance, affects the performance of the controller in reaching the desired pose. Moreover, the utilized exoskeleton is a redundant robot that can have different inverse kinematic solutions each time the robot performs a task, especially for different initial values of robot configuration. To have reliable results for the performance of the proposed controller, the experiment is repeated 6 times, and the means and standard deviations for position and orientation errors are computed. The utilized control and adaptation gains are as follows:
. The initial condition for the adaptation function of (31) is chosen in a way that satisfies the physical consistency condition and . DH parameters and the Jacobian matrix are derived based on the convention explained in [26].Fig. 4(a) demonstrates the mean and standard deviation of the position errors of the endeffector in each direction, and Fig. 4(b) shows a sample of desired path tracking. As can be seen in Fig. 4, the proposed adaptive VDC has excellent performance in tracking the desired path in Cartesian space. The results demonstrate that the adaptive VDC is robust against unknown timevarying disturbances and stably maintains the system throughout the working time. The average Root Mean Square Error (RMSE) of position errors for all experiments is . Moreover, Fig. 5 shows the mean and standard deviation of the quaternion orientation errors of the endeffector computed with (38). It can be seen that in the presence of unknown external timevarying disturbances, which are directly applied to the endeffector, adaptive VDC exhibits perfect tracking of the desired orientation. The average RMSE of orientation errors for all experiments is .
Vi Conclusions
In this paper, we proposed a new adaptive VDC approach that uses NAL to estimate unknown inertial parameters. The original adaptive VDC requires tuning of 13 adaptation gains and assigning 26 upper and lower parameter limits, which makes the process excessively tedious. However, the proposed adaptive VDC removes all the mentioned burdens and requires only one adaptation gain for all the rigid bodies. First, we reformulated the dynamic equations of the VDC context to put them into 10 unique inertial parameter vectors, making them compatible with the NAL function. Then, it was proved that employing the NAL function as the adaptation function for adaptive VDC culminates in asymptotic stability. Finally, the experimental results of applying the new adaptive VDC to the redundant upperlimb exoskeleton demonstrated the excellent performance of the VDC approach.
Our future work target is the control problem of two master exoskeleton robots that command heavyduty hydraulic mobile manipulators as slaves for handling common objects in humanlike manipulation scenarios with high scaling. As we showed in our earlier work [3], the high motion/force scaling must be tackled in the controller design and stability analysis, which makes analysis much more complicated. Therefore, a novel adaptive VDC scheme developed in this paper will greatly lessen the implementation burden of this highly complex MMMS system.
References
 [1] Sheridan, Thomas B. Telerobotics, automation, and human supervisory control. MIT press, 1992.
 [2] Azimifar, Farhad, et al. ”Transparency performance improvement for multimaster multislave teleoperation systems with external force estimation.” Transactions of the Institute of Measurement and Control 40.13 (2018): 38513859.
 [3] Lampinen, Santeri, et al. ”ForceSensorLess Bilateral Teleoperation Control of Dissimilar MasterSlave System With Arbitrary Scaling.” IEEE Transactions on Control Systems Technology (2021).
 [4] Riani, A., et al. ”Adaptive integral terminal sliding mode control for upperlimb rehabilitation exoskeleton.” Control Engineering Practice 75 (2018): 108117.
 [5] Brahmi, Brahim, et al. ”A new integral secondorder terminal sliding mode control with time delay estimation for an exoskeleton robot with dynamics uncertainties.” Proceedings of the 5th International Conference of Control, Dynamic Systems, and Robotics (CDSR’18). Vol. 113. 2018.
 [6] Brahmi, Brahim, et al. ”Adaptive force and position control based on quasitime delay estimation of exoskeleton robot for rehabilitation.” IEEE Transactions on Control Systems Technology 28.6 (2019): 21522163.
 [7] Brahmi, Brahim, et al. ”Novel adaptive impedance control for exoskeleton robot for rehabilitation using a nonlinear timedelay disturbance observer.” ISA transactions 108 (2021): 381392.
 [8] Brahmi, Brahim, et al. ”Impedance learning control for physical humanrobot cooperative interaction.” Mathematics and Computers in Simulation 190 (2021): 12241242.
 [9] Zhang, GaoWei, et al. ”Integrated observerbased fixedtime control with backstepping method for exoskeleton robot.” International Journal of Automation and Computing 17.1 (2020): 7182.
 [10] Wu, Qingcong, Bai Chen, and Hongtao Wu. ”Rbfnbased adaptive backstepping sliding mode control of an upperlimb exoskeleton with dynamic uncertainties.” IEEE Access 7 (2019): 134635134646.
 [11] Rahmani, Mehran, and Mohammad Habibur Rahman. ”An upperlimb exoskeleton robot control using a novel fast fuzzy sliding mode control.” Journal of Intelligent & Fuzzy Systems 36.3 (2019): 25812592.
 [12] Yip, P. Patrick, and J. Karl Hedrick. ”Adaptive dynamic surface control: a simplified algorithm for adaptive backstepping control of nonlinear systems.” International Journal of Control 71.5 (1998): 959979.
 [13] Zhu, WenHong. Virtual decomposition control: toward hyper degrees of freedom robots. Vol. 60. Springer Science & Business Media, 2010.
 [14] Shen, Youhao, et al. ”An Adaptive Control Method for Electrohydrostatic Actuator Based on Virtual Decomposition Control.” 2020 AsiaPacific International Symposium on Advanced Reliability and Maintenance Modeling (APARM). IEEE, 2020.
 [15] Koivumäki, Janne, and Jouni Mattila. ”Stabilityguaranteed forcesensorless contact force/motion control of heavyduty hydraulic manipulators.” IEEE Transactions on Robotics 31.4 (2015): 918935.
 [16] Humaloja, JukkaPekka, et al. ”Decentralized Observer Design for Virtual Decomposition Control.” IEEE Transactions on Automatic Control 67.5 (2021): 25292536.
 [17] Brahmi, Abdelkrim, et al. ”Adaptive Control of Mobile Manipulator Robot based on Virtual Decomposition Approach.” ICINCO (2). 2016.
 [18] Antonelli, Gianluca, Fabrizio Caccavale, and Stefano Chiaverini. ”Adaptive tracking control of underwater vehiclemanipulator systems based on the virtual decomposition approach.” IEEE Transactions on Robotics and Automation 20.3 (2004): 594602.
 [19] Luna, Cristobal Ochoa, et al. ”Virtual decomposition control of an exoskeleton robot arm.” Robotica 34.7 (2016): 15871609.
 [20] Wensing, Patrick M., Günter Niemeyer, and JeanJacques E. Slotine. ”Observability in inertial parameter identification.” arXiv preprint arXiv:1711.03896 (2017).
 [21] Lee, Taeyoon, Jaewoon Kwon, and Frank C. Park. ”A natural adaptive control law for robot manipulators.” 2018 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS). IEEE, 2018.
 [22] Wensing, Patrick M., Sangbae Kim, and JeanJacques E. Slotine. ”Linear matrix inequalities for physically consistent inertial parameter identification: A statistical perspective on the mass distribution.” IEEE Robotics and Automation Letters 3.1 (2017): 6067.
 [23] Petrović, Goran R., and Jouni Mattila. ”Mathematical modelling and virtual decomposition control of heavyduty parallel–serial hydraulic manipulators.” Mechanism and Machine Theory 170 (2022): 104680.
 [24] Jazar, Reza N. Theory of applied robotics. Springer Science+ Business Media, LLC, 2010.
 [25] Chiaverini, Stefano, and Bruno Siciliano. ”The unit quaternion: A useful tool for inverse kinematics of robot manipulators.” Systems Analysis Modelling Simulation 35.1 (1999): 4560.
 [26] Sciavicco, Lorenzo, and Bruno Siciliano. Modelling and control of robot manipulators. Springer Science & Business Media, 2001.
Appendix
Via Restructured Dynamic Equations
To rewrite the matrix, it is enough to rewrite . By expanding , we get to