Compliant-Control-Based Assisted Walking with Mobile Manipulator

In this paper, a new approach involving the use of a mobile manipulator to assist humans with mobility impairments to walk is proposed. First, in order to achieve flexible interaction between humans and mobile manipulators, we propose a variable admittance controller that can adaptively regulate the virtual mass and damping parameters based on the interaction forces and the human motion intention predicted using the fuzzy theory. Moreover, a feedforward velocity compensator based on a designed state observer is proposed to decrease the inertia resistance of the manipulator, effectively enhancing the compliance of the human–robot interaction. Then, the configuration of the mobile manipulator is optimized based on a null-space approach by considering the singularity, force capacity, and deformation induced by gravity. Finally, the proposed assisted walking approach for the mobile manipulator is implemented using the human–robot interaction controller and the null-space controller. The validity of the proposed controllers and the feasibility of assisted human walking are verified by conducting a set of tests involving different human volunteers.


Introduction
In recent years, the elderly population has globally increased rapidly, and there is a great need for assisted walking robots to help the elderly or disabled people walk.Among the multitude of devices and robots, especially walking assistants, wheelchairs have long served as a good option for this group of people.However, wheelchairs can lead to muscle weakness and joint stiffness, which can hamper the recovery of such people's athletic abilities.Therefore, many scholars have developed and studied various types of devices for providing assistance to the elderly [1][2][3][4][5][6], such as powered lower-limb orthoses, wearable exoskeletons, intelligent wheelchairs, walkers, and assistive robots.Robot-assisted exercise training can provide sustainable intensive rehabilitation treatment to the elderly and help them recover their motion functions [7].Traditional walking aids include the U-Step walking stabilizer, whose velocity can be adjusted by cranking or height changes.However, the use of such assistants increases the risk of falling [8].Smart walkers can provide superior walking assistance, but several of these devices are extremely complex to use [4].
Mobile manipulators have been widely used in walking assistance robots, which generally consist of a mobile base and a manipulator.Herein, the main challenge is realizing cooperative motion control between the mobile manipulator and the human.Mobile manipulators have been used to assist humans with opening doors and carrying heavy objects [9][10][11].To obtain the desired human-robot commands, several controllers have been proposed, including compliant control, human motion intention recognition, and redundant solutions.
In the task of assisted walking, the interaction force from the human hand should be decreased by using compliant control.This allows humans to interact effortlessly with the mobile manipulator.
In general, compliant control can be divided into active and passive compliance [12].Active compliance can be further subdivided into impedance control, admittance control, and hybrid force/position control [13].Passive compliance has a critical disadvantage in that it cannot be used to vary the compliance effect depending on the external environment, and it is effective only in a few simple applications [14].
Compared with passive compliance, with active compliance, it is easier to realize and modify the compliance effect of the robot.In manipulators, admittance controllers are widely used to achieve compliance control [15].Thus far, active compliance has invariably been implemented by regulating the virtual damping, virtual stiffness [16], virtual mass, or simultaneously changing the virtual damping parameter and virtual mass.For example, Jiang et al. [17] proposed a variable admittance algorithm based on Fuzzy Sarsa (λ)-learning to regulate the virtual damping parameters.Kuang et al. [18] proposed a variable admittance controller based on fuzzy reference learning control by changing the virtual damping parameters in the admittance model to ensure that the manipulator moves according to the surgeon's motion intention, thus exhibiting a certain level of compliance.However, these methods can provide good results only on specific tasks, and they have relatively strong limitations.In addition, radial basis function (RBF) neural networks have been widely used to adjust the damping parameters of admittance controllers [19].Sharkawy et al. [20] proposed a method to modify only the virtual mass parameter through online training of a multilayer feedforward neural network.Thereafter, they proposed another method to simultaneously regulate the virtual mass and virtual damping by using a Jordan recurrent neural network [21].
In the task of assisted walking, it is also necessary to recognize the human motion intention to realize forward, backward, and steering operations, ensuring that mobile manipulators can effectively track human motion and provide assistance.
To identify human motion intention, video-based image recognition methods such as the OpenPose algorithm [22,23] are widely used to identify human hands and limbs.However, during the walking process, due to the camera view, it is usually difficult to obtain all the image information of the limbs.
In addition, sensor-based human behavior recognition methods [24] have emerged recently, and they can be subdivided into methods based on wearable and non-wearable sensors.The wearable sensors used in such methods include electroencephalography (EEG) collectors, electrooculography (EOG) collectors [25,26] and electromyography (EMG) collectors [27].However, it is common for signals to be easily interfered with, difficult to measure, and the wearing process to be complex, which brings significant inconvenience to daily care.
Therefore, many scholars have studied intention-recognition methods based on nonwearable sensors.In [28], the use of high-order Markov chains to gauge human motion intentions was proposed.Tomislav et al. [29] proposed a method to predict human behavioral intentions by using a hidden Markov model.Palm et al. [30] proposed a collisionavoidance method based on fuzzy control rules and verified its effectiveness by conducting a simulation and an experiment.However, the identification of human motion intention during walking is fraught with several novel challenges, like how to accurately identify the intention of human movement without wearing any sensors in the contact interaction.In such cases, information pertaining to interaction forces and velocities effectively serves as the source data.Lecours et al. [31] proposed a variable admittance controller that can infer human intentions by using velocity and acceleration data.
In the task of assisted walking, the mobile manipulator should also be able to move reasonably.Generally, a mobile manipulator is a system that possesses redundant degrees of freedom.For solving such redundancy, usually the extended Jacobi method [32] and the gradient projection method [33] are employed.Cetin et al. [34] proposed an extended Jacobi matrix formulation for a fixed laser or optic camera mounted on the first joint to accurately track the end-effector.Boundec et al. [35] employed the gradient projection method to realize obstacle avoidance.Under this method, they developed three objective functions to achieve obstacle avoidance, help the robot avoid exceeding its physical limits, and prevent the robot from colliding with its own base.
Therefore, we propose an approach by employing a mobile manipulator to assist human walking in this study, which involves having an experienced nurse help a patient walk.The mobile manipulator used in this work consists of a four-wheeled mobile robot and a six-DOF (degree of freedom) manipulator.Assisted walking based on mobile manipulators can provide sustainable, intensive rehabilitation for patients with mobility impairments.The main challenges of assisted walking depend on the flexibility of human-robot interaction, the configuration optimization of the mobile manipulator, and the smoothness and robustness of the human-robot collaborative controller.
To implement compliant-control-based walking assistance by using a mobile manipulator, we propose a human motion intention recognition method and a variable admittance controller based on the fuzzy theory, which can accurately predict human motion intention.Simultaneously, we eliminate the inertial resistance of the manipulator by means of velocity feedforward compensation provided using a developed state observer, which can effectively enhance the compliance of the interaction between the human and the mobile manipulator.Furthermore, to solve the redundant degrees of freedom of the mobile manipulator, a null-space controller based on the gradient projection method is proposed.This controller can identify the optimal manipulator configuration.To guarantee that the position and posture of the end-effector remain unchanged during assisted walking, a position and posture controller for the end-effector is proposed based on the proportional control algorithm.The main contributions of this paper are as follows: 1.
A variable admittance controller is proposed to simultaneously adjust virtual mass and virtual damping based on the fuzzy theory, and a state-observer-based velocity compensator is also proposed.

2.
A null-space controller is presented to optimize the mobile manipulator's configuration so that it can avoid singularities, enhance force capacity in the vertical direction, and optimize gravity-induced deformation.3.
By using the aforementioned controllers, the mobile manipulator is able to successfully provide walking assistance to humans.
The remainder of this paper is structured as follows: In Section 2, the kinematics model of the mobile manipulator is established.In Section 3, a human-robot-compliant interaction method is proposed based on the fuzzy theory.In Section 4, we propose a method to optimize the configuration of the mobile manipulator based on a null-space controller.Section 5 describes a series of experiments that verify the effectiveness and correctness of the methods proposed herein.Our concluding remarks are presented in Section 6.

Kinematics Model of Mobile Manipulator
In this paper, a mobile manipulator consisting of a four-wheeled differential mobile robot and a 6-DOF manipulator is researched, as illustrated in Figure 1.
To coordinate the kinematic relationship between the mobile platform and the manipulator joints, we first develop the integrated kinematics of the mobile manipulator.The world coordinate system, mobile platform coordinate system, manipulator base coordinate system, and manipulator end-effector coordinate system are represented by Σ w , Σ b , Σ 0 , and Σ ee , as depicted in Figure 1.To coordinate the kinematic relationship between the mobile platform and the manipulator joints, we first develop the integrated kinematics of the mobile manipulator.The world coordinate system, mobile platform coordinate system, manipulator base coordinate system, and manipulator end-effector coordinate system are represented by Σw, Σb, Σ0, and Σee, as depicted in Figure 1.
Herein, the Denavit-Hartenberg (D-H) method is used for kinematic modeling of the mobile manipulator.qb = [xb yb φ] T represents the position and orientation of the mobile base coordinate system in the world coordinate system, where xb, yb, and φ are respectively the x-coordinate, y-coordinate, and rotation angle of the mobile platform in the world coordinate system.qm = [q1 q2 q3 q4 q5 q6] T is the joint angle vector of the manipulator, and qj is the angle of the j-th (j = 1, 2, 3, 4, 5, 6) joint.Therefore, the transformation matrix from Σb to Σw can be obtained as follows: ( ) where P(qb) and ub are the kinematic constraint and velocity vector of the mobile base, respectively.vb and ωb are the longitudinal and heading angle speeds of the mobile base, respectively.
In addition, the transformation matrix from the end-effector to the manipulator base is defined as 0 6 T , and the transformation matrix from the manipulator base to the mobile platform is defined as 0 b T , which are both the constant matrixes.Then, the transformation matrix from the end-effector to the world coordinate system can be expressed as Herein, the Denavit-Hartenberg (D-H) method is used for kinematic modeling of the mobile manipulator.q b = [x b y b φ] T represents the position and orientation of the mobile base coordinate system in the world coordinate system, where x b , y b , and φ are respectively the x-coordinate, y-coordinate, and rotation angle of the mobile platform in the world coordinate system.q m = [q 1 q 2 q 3 q 4 q 5 q 6 ] T is the joint angle vector of the manipulator, and q j is the angle of the j-th (j = 1, 2, 3, 4, 5, 6) joint.Therefore, the transformation matrix from Σ b to Σ w can be obtained as follows: (1) where P(q b ) and u b are the kinematic constraint and velocity vector of the mobile base, respectively.v b and ω b are the longitudinal and heading angle speeds of the mobile base, respectively.
In addition, the transformation matrix from the end-effector to the manipulator base is defined as 0 6 T, and the transformation matrix from the manipulator base to the mobile platform is defined as b 0 T, which are both the constant matrixes.Then, the transformation matrix from the end-effector to the world coordinate system can be expressed as w e T = w b T b 0 T 0 6 T. By using the matrix w e T, the Jacobian matrix of the mobile manipulator can be constructed, including the differential kinematic model of the mobile platform.Herein, we assume that the DOF of the manipulator in the task space is r, the DOF in the joint space is m, the DOF of the mobile robot is n b , and its degree of mobility under the nonholonomic constraints is b.The total DOF of the mobile manipulator is n = n b + m.Then, the kinematic model of the mobile manipulator can be expressed as follows: where r = m = 6, J ∈ R r×n is the Jacobian matrix of the mobile manipulator without considering the nonholonomic constraints; J b ∈ R r×n b is the Jacobian matrix of the mobile platform, and it is related to the matrix in (1); J m (q m ) ∈ R r×m is the Jacobian matrix of the manipulator; .
x ∈ R r is the velocity vector of the end-effector; . q b ∈ R n b is the velocity vector of the mobile platform; u m ∈ R m is the velocity vector of the manipulator joints; and r×(b+m) is the Jacobian matrix of the mobile manipulator considering the nonholonomic constraints.
By using the mobile manipulator model (3), we can reasonably allocate u b and u m .

Human-Robot Compliant Interaction Based on Fuzzy Control
To realize compliant interaction between the manipulator and the human user in the process of providing walking assistance, the admittance control method is employed in this study.A virtual "Inertia-Damper-Spring" model is established between the manipulator and the human hand.The target velocity, or angular velocity, of the end-effector is calculated based on the interaction force or torque exerted on the manipulator by the human hand.This method mainly realizes the compliant interaction between the mobile manipulator and the human hand in the horizontal direction and provides adequate support to the human in the vertical direction.
Herein, we focus on the performance of force interaction between the human and the manipulator while ignoring the interaction torque and spring effect.The admittance model used in this study can be expressed as follows: where M v ∈ R 3×3 is a virtual mass matrix; M x , M y , and M z are virtual masses in the x, y, and z-directions, respectively; C v ∈ R 3×3 is a virtual damping matrix; C vx , C vy , and C vz denote virtual damping in the x, y, and z-directions, respectively; and F h ∈ R 3×1 is the human-robot interaction force.For the single-DOF admittance model, the transfer function between the end-effector velocity and the interaction force can be expressed as follows: where β is the ratio of virtual mass to virtual damping.C v and β affect the steady-state response and poles of the interaction system, thereby affecting the dynamic response characteristics of the system.If β is small, the response speed of the system will be faster, but the system will be prone to overshoot and vibration.By contrast, if β is large, the system will be more stable, but its response time will be longer.
In walking-assistance technology, when the mobile manipulator moves rapidly with the human, the interaction system needs a small damping coefficient and offers superior compliance.When the human user decelerates or stops moving, the interactive system needs a large damping coefficient to offer superior stability and accuracy.However, it is difficult to implement such a regulation by using the traditional admittance control scheme with fixed parameters, and for this reason, a variable-structure admittance controller is proposed herein.
In [17,36], the variable structure admittance controller only adjusts C v while keeping M v unchanged, resulting in the change of β.However, M v also affects the magnitude of F h in the process of acceleration and deceleration.To further improve the performance of the interaction system, we propose that M v should be varied with C v .Therefore, according to [31], we define the virtual mass as a piecewise function of virtual damping, as expressed in the following equation: A schematic diagram of piecewise function ( 6) is shown in Figure 2.
where C vi is the virtual damping coefficient in the i direction, and i represents the x, y, and z-directions, which is a scalar.β is the ratio of steady-state virtual mass to virtual damping, and β > 0; ε is the adjustment factor of steady-state virtual mass to the virtual damping ratio, and it satisfies 0 < ε < 1; µ is the smoothness factor, and it satisfies µ > 0; i represents the x, y, and z-directions.C 0 represents the virtual damping coefficient in steady state, and C 0 = 50 Ns/m in this paper.
[31], we define the virtual mass as a piecewise function of virtual damping, as expressed in the following equation: A schematic diagram of piecewise function ( 6) is shown in Figure 2.
where   is the virtual damping coefficient in the i direction, and i represents the x, y, and z-directions, which is a scalar.β is the ratio of steady-state virtual mass to virtual damping, and β > 0; ε is the adjustment factor of steady-state virtual mass to the virtual damping ratio, and it satisfies 0 < ε < 1; μ is the smoothness factor, and it satisfies μ > 0; i represents the x, y, and z-directions.C0 represents the virtual damping coefficient in steady state, and C0 = 50 Ns/m in this paper.To realize the adaptive variable structure admittance control, we adjust the virtual damping according to the end-effector velocity and the interaction force between the manipulator and the human.In view of the robustness of fuzzy control, we build a fuzzy function to describe the relationship between   ,  ℎ , and Fhi. .
where ( ) is a fuzzy function of  ℎ and Fhi, which will be covered in the next section.Cmax and Cmin are, respectively, the upper and lower limits of the virtual damping adjustment, which corresponds to the upper and lower limits of fuzzy function values; Fhi is the interaction force; and  ℎ is the end-effector velocity.In this paper, the minimum and maximum values of Cvi are set to be Cmin = 30 Ns/m and Cmax = 70 Ns/m.The above equation can be used to effectively regulate the relationship between the virtual mass and virtual damping to achieve superior human-robot interaction performance in the case of assisted walking.Then, the desired end-effector velocity vd is determined using the modified admittance controller.By (3), we can obtain the desired ub and um, and the movement of the mobile manipulator according to the human user's motion intention during the process of walking assistance is completed simultaneously.

Prediction of Human Motion Intention Using Fuzzy Theory
In this section, the fuzzy theory is employed as it can well describe the complex relationship between the required damping coefficients and human-robot interaction motion and forces and has good robustness for the uncertainty of human-robot interaction.Based To realize the adaptive variable structure admittance control, we adjust the virtual damping according to the end-effector velocity and the interaction force between the manipulator and the human.In view of the robustness of fuzzy control, we build a fuzzy function to describe the relationship between C vi , v hi , and F hi .
where f (v hi , F hi ) is a fuzzy function of v hi and F hi , which will be covered in the next section.C max and C min are, respectively, the upper and lower limits of the virtual damping adjustment, which corresponds to the upper and lower limits of fuzzy function values; F hi is the interaction force; and v hi is the end-effector velocity.In this paper, the minimum and maximum values of C vi are set to be C min = 30 Ns/m and C max = 70 Ns/m.The above equation can be used to effectively regulate the relationship between the virtual mass and virtual damping to achieve superior human-robot interaction performance in the case of assisted walking.Then, the desired end-effector velocity v d is determined using the modified admittance controller.By (3), we can obtain the desired u b and u m , and the movement of the mobile manipulator according to the human user's motion intention during the process of walking assistance is completed simultaneously.

Prediction of Human Motion Intention Using Fuzzy Theory
In this section, the fuzzy theory is employed as it can well describe the complex relationship between the required damping coefficients and human-robot interaction motion and forces and has good robustness for the uncertainty of human-robot interaction.Based on a series of pretests, the interaction force and velocity of the end-effector can be used to represent the human motion intention, and the suitable damping is further analyzed to improve human comfort.Further, the fuzzy theory is employed to establish the relationship between virtual damping and human-robot interaction.
(a) When the interaction force F h or the end-effector velocity v h increases gradually from 0, it can be judged that the human user intends to start moving, and accordingly, the damping coefficient should be decreased gradually.(b) When the directions of F h and v h are the same and their magnitudes are large, it can be judged that the human user intends to accelerate significantly, and accordingly, a small damping coefficient is required.(c) When the directions of F h and v h are not the same but their magnitudes are large, it can be judged that the human user intends to decelerate significantly, and accordingly, a large damping coefficient is required.
(d) When F h and v h are extremely small, it can be judged that the human user intends to maintain the current position or move with greater precision, and accordingly, a large damping coefficient is required.
Then, the virtual damping C v in (7) can be regulated based on the prediction results.Thus, the inputs to the proposed fuzzy system are F h and v h , and its output is the virtual damping C v .Assuming that the maximum and minimum values of F h in the three directions are identical as F max and F min , respectively, and the maximum and minimum values of v h in the three directions are defined to be equal as v max and v min , respectively, the domains of F h and v h are F h :[F min F max ] and v h :[v min v max ], respectively.Considering the directionality of force and velocity, F h can be divided into seven fuzzy subsets {NB (Negative Big), NM (Negative Medium), NS (Negative Small), ZO (Zero), PS (Positive Small), PM (Positive Medium), and PB (Positive Big)}, and the end-effector velocity v h can be divided into seven fuzzy subsets {NB (Negative Big), NM (Negative Medium), NS (Negative Small), ZO (Zero), PS (Positive Small), PM (Positive Middle), PB (Positive Big)}.Similarly, the virtual damping C v can be divided into five fuzzy subsets: {VS (Very Small), S (Small), M (Medium), B (Big), VB (Very Big)}.
The triangular membership function is employed to fuzzify F h , v h , and C v .The parameters in the membership function are obtained experimentally, and the distances between the peaks of F h and v h are, respectively, 3 and 0.04.The designed membership functions of F h , v h , and C v are depicted in Figures 3 and 4. A fuzzy rule table is then developed, as shown in Table 1.
damping coefficient should be decreased gradually.(b) When the directions of Fh and vh are the same and their magnitudes are large, it can be judged that the human user intends to accelerate significantly, and accordingly, a small damping coefficient is required.(c) When the directions of Fh and vh are not the same but their magnitudes are large, it can be judged that the human user intends to decelerate significantly, and accordingly, a large damping coefficient is required.(d) When Fh and vh are extremely small, it can be judged that the human user intends to maintain the current position or move with greater precision, and accordingly, a large damping coefficient is required.
Then, the virtual damping Cv in ( 7) can be regulated based on the prediction results.Thus, the inputs to the proposed fuzzy system are Fh and vh, and its output is the virtual damping Cv.Assuming that the maximum and minimum values of Fh in the three directions are identical as Fmax and Fmin, respectively, and the maximum and minimum values of vh in the three directions are defined to be equal as vmax and vmin, respectively, the domains of Fh and vh are Fh:[Fmin Fmax] and vh:[vmin vmax], respectively.Considering the directionality of force and velocity, Fh can be divided into seven fuzzy subsets {NB (Negative Big), NM (Negative Medium), NS (Negative Small), ZO (Zero), PS (Positive Small), PM (Positive Medium), and PB (Positive Big)}, and the end-effector velocity vh can be divided into seven fuzzy subsets {NB (Negative Big), NM (Negative Medium), NS (Negative Small), ZO (Zero), PS (Positive Small), PM (Positive Middle), PB (Positive Big)}.Similarly, the virtual damping Cv can be divided into five fuzzy subsets: {VS (Very Small), S (Small), M (Medium), B (Big), VB (Very Big)}.
The triangular membership function is employed to fuzzify Fh, vh, and Cv.The parameters in the membership function are obtained experimentally, and the distances between the peaks of Fh and vh are, respectively, 3 and 0.04.The designed membership functions of Fh, vh, and Cv are depicted in Figures 3 and 4. A fuzzy rule table is then developed, as shown in Table 1.The virtual damping is obtained through defuzzification by using the discrete centroid method: The virtual damping is obtained through defuzzification by using the discrete centroid method: where c pk denotes the horizontal coordinates of the peak of the membership function with virtual damping, k = 1, 2, 3, 4, 5.By using (8), the human motion intention can be predicted in terms of the interaction force and end-effector velocity, and then the effective virtual damping can be obtained using (7).Therefore, variable admittance control is realized by adjusting C v based on the fuzzy theory and M v by using ( 6) to obtain the designed end-effector velocity with good compliance performance.

Velocity Feedforward Compensation Based on the State Observer
In the human-robot interaction process, the manipulator reacts according to the interaction force, and the robot's inertia inevitably offers resistance to the human hand, resulting in delayed responses.Therefore, a state observer is designed to predict the joint torque due to the interaction force by considering the current motion state, and a velocity feedforward controller is proposed to append a compensatory velocity to the end-effector and improve system compliance.
During assisted human walking with a mobile manipulator, the desired velocity induced by the interaction force from the human hand is allocated to the manipulator and the mobile platform, respectively.Compliance is mainly reflected in the interaction between the manipulator and the human user.Thus, the effect of the mobile platform on the compliance of the interaction system can be neglected.
Based on the above assumptions, a state observer is designed for the manipulator to predict the joint torque induced by the interaction force.The dynamic model of the manipulator can be written as follows: where M(q m ) ∈ R m×m is the inertia term of the manipulator, C q m , .q m ∈ R m×m is the centrifugal and Coriolis term, G(q m ) ∈ R m×1 is the gravity term, F f q m , .q m ∈ R m×1 is the friction term, τ m ∈ R m×1 is the joint driving torque of the manipulator, and τ ext ∈ R m×1 is the joint torque generated due to the external force acting on the end-effector.
According to (9), when the manipulator is subjected to the interaction force exerted by the human user, the joint torque caused by the external force can be expressed as follows: Based on (10), a nonlinear disturbance state observer can be constructed to estimate the joint torque generated by the external force.The estimated values of M(q m ) and N q m , .q m are M(q m ) and N q m , .q m , respectively, and their uncertainty factors are M, N, and M(q m ) is a symmetric positive-definite matrix as well.We can obtain A disturbance vector τ d ∈ R m×1 is designed as follows: If there are adequate and reasonable estimates of M(q m ) and N q m , .q m , the joint torque τ ext of the manipulator due to the external force is almost equal to the disturbance vector τ d .By combining ( 10)-( 12), we can obtain: Thus, the following state observer is proposed to observe the joint torque τ ext . .
where L is the observer gain matrix.
In practice, to obtain an accurate .. q m , the first-order difference of .q m and the Kalman filter are employed here.
Then, to verify the stability of the proposed state observer, its error vector e dis ∈ R m×1 is defined as follows: The change rate of this error is as follows: . e dis = .
During assisted walking, the mobile manipulator moves at a low velocity, and the disturbance vector varies slowly.In [37], the designed observer in ( 14) can be regarded as a non-dynamic nonlinear disturbance observer (NDOB).Therefore, it is reasonably assumed that: .
By ensuring that the proposed observer is asymptotically stable over a broad range, we define L as follows: where ν is a positive constant.Therefore, the proposed state observer (14), constrained by (20), is stable and effective.Then, we can obtain a reasonable τ ext by employing the observer (14).
Based on ( 14), the following feedforward velocity is designed to eliminate the robot's inertia. .
τd is the feedforward gain matrix, and T is the control period.
From (21), we can obtain .q f with the help of τd , .
τd , which can be used to improve the compliance of interactive systems.
However, upon the direct use of ( 21), M and N induced notable manipulator jitter in the experiments.Because .. q m is obtained after Kalman filtering, M has little effect on this phenomenon, and the N is the main driver of this phenomenon, possibly because of the friction torque and joint vibration.In addition, the position and velocity of the end-effector are mainly determined by the first three joints of the manipulator, and the latter three joints mainly affect the posture of the end-effector.Additionally, these noises have stronger effects on the latter three joints of the manipulator because of their small inertia.Therefore, to avoid unreasonable movements of the latter three joints and prevent manipulator jitter, the feedforward velocity compensation is applied only to the first three joints of the manipulator in this study.Herein, we can define K f τd , .τd in (21) as where τdl is the joint torque threshold to activate the process of velocity feedforward; . τdl , .τdu are the thresholds of the change rate for the joint torques to initiate the process of velocity feedforward; s a and s b are smoothing factors; and j denotes the j-th joint of the manipulator.
According to (22), the process of velocity feedforward is activated only when τd is relatively large or changes rapidly.Thus, the adverse effects of uncertain factors on manipulator motion can be avoided.
According to (3) and ( 21), . q f is mapped to ∆v f in the Cartesian space as follows: To guarantee that the mobile manipulator can better assist human walking, the endeffector of the manipulator should track the human user's motion, such that the interactive system simultaneously has good compliance and motion-tracking performance.Moreover, it should also provide a certain vertical force to support the human hand while walking.Therefore, the magnitudes of v d in the x and y-directions are computed using the proposed variable admittance controller (4), while the velocity in the z-direction is computed using the following constant admittance controller: Finally, the proposed variable admittance controller can be expressed as follows: The desired end-effector velocity required for the mobile manipulator is then modified as Then, the desired v r can be expressed as follows: Therefore, M x and M y and C x and C y in (25) are revised based on ( 6)- (8).Then, v r is projected onto u b and u m by (3).In this manner, compliant interaction between the manipulator and the human user is realized, and the mobile manipulator can move according to the human motion intention.

Configuration Optimization of Mobile Manipulator
Given that the total DOF of the mobile manipulator when considering the nonholonomic constraints is b + r = 8, its Jacobian matrix J u (6 × 8) is non-square.Because of the existing redundant DOFs of the mobile manipulator, there may exist infinite sets of joint velocities for a desired end-effector velocity in Cartesian space.
To guarantee that the mobile manipulator has the optimal configuration when tracking the desired end-effector commands (v r ), it is necessary to implement de-redundancy to realize motion control of the mobile manipulator.
In this study, a null-space controller based on the projected gradient method is proposed.

Null-Space Controller
Considering v r as the desired command for the mobile manipulator, the secondary task is to implement an optimal configuration that can be defined on the basis of the null space of the mobile manipulator.The redundant solution of the mobile manipulator in (3) can be expressed as follows: where J † u is the pseudo-inverse of the Jacobian matrix of the mobile manipulator, I − J † u J u is the orthographic projection operator representing the null space N (J u ) of the mobile manipulator, and u n is the joint velocity vector required to perform secondary tasks.
The goal of the null-space control herein is to ensure that the mobile manipulator maintains its optimal configuration.Therefore, ( 27) is directly used to optimize the configuration of the mobile manipulator under null-space control.
(1) Before the manipulator reaches the optimal configuration, the desired velocity u rn is converted into the commands u b and u m for the mobile robot and the manipulator through J u (q) −1 .The mobile robot and the manipulator move, respectively, with the decomposed commands decided by ( 3) and ( 27). ( 2) After reaching the optimal configuration, u m along the horizontal direction is mostly distributed to the moving platform; the velocity in the vertical direction should be close to zero to ensure the provision of vertical support forces.

Null-Space Tasks
In this section, the projected gradient method is employed to solve the null-space tasks in (27) by designing an optimization function F(u).Then, with u n = ∇ u F(u), we can implement the required secondary tasks.
As a secondary task, u n is constantly projected onto the null space N (J u ) of J u , and it does not affect the designed task v r in the task space.Here, F(u) should be a differentiable function.
To design a feasible F(u), we mainly consider the following three aspects: (1) The singular configuration cannot appear when the mobile manipulator flexibly interacts with the human user, meaning that the manipulator should be as far away from the singularity as possible.
(2) The manipulator can adequately withstand external forces from the human user's hand.
(3) The effect of the manipulator's gravity on itself should be as small as possible to ensure that the manipulator can adequately support the human hand in the z-direction.
For the first target, the joint velocity of the manipulator can be considered as the following velocity ellipsoid [38]: The velocity ellipsoid can be the operational evaluation index of the manipulator.The larger the value of E 1 , the wider the non-singular configuration range of the manipulator.
For convenient calculation, the square of the velocity ellipsoid is taken as the first objective function: For the second target, the external force can be considered from the manipulator force ellipsoid [38].However, the force ellipsoid and the velocity ellipsoid are reciprocals of each other [39].If the capacity of the end-effector to withstand external forces is increased, the corresponding improvement in the force ellipsoid will cause the manipulator to move closer to the singularity configuration.Therefore, it would be impossible to improve the capacity of the end-effector to withstand external forces as the manipulator moves farther away from the singular configuration.Thus, the force ellipsoid and the velocity ellipsoid should be balanced.
To ensure that the mobile manipulator can assist human walking, the interaction force between the end-effector and the human hand should be as small as possible in the horizontal direction, and it should perform a supporting role in the vertical direction.Therefore, we only need to increase the force capacity in the vertical direction.
The normalized vector of the manipulator joint torque is defined as follows: where τ norm is the normalized joint torque vector, W lim is the proportional matrix considering the joint torque limits, and τ mlimj is the limit torque of joint j.
With the definition of force ellipsoid in [38], the constraint for the manipulator joint torques is defined as follows: When the end-effector of the manipulator is subjected to an external force, F ext , the joint torque is mainly induced by F ext and the manipulator's own gravity, as follows: Then, we can project the gravity component G(q m ) in (32) onto its Cartesian space.
By substituting (33) in (31), we have According to (34), if τ T norm τ norm is constant, meaning that the manipulator joint torque is constant, a reduction in J m W lim W lim J T m can lead to an increase in the external force that the manipulator can bear.However, directly considering J m W lim W lim J T m reduction as the optimization goal will increase the force capacity of the manipulator in all directions [40], which may weaken the force application ability in the vertical direction under the constraints of joint maximum torques.Moreover, there are different requirements for force capacity in the horizontal and vertical directions during assisted walking.
Therefore, it is assumed that the end-effector has the equivalent stiffness K opt ∈ R 3×3 , and the stiffness K opt of the end-effector is obtained through the pretests to guarantee it can supply enough vertical force under a permissible vertical deformation.The virtual deformation of the end-effector due to F ext and G m are ∆x ext = F ext /K opt and ∆x g = G m /K opt , respectively.Then, by substituting these expressions into (34), we obtain According to (35), by reducing K opt J m W lim W lim J T m K opt , we can increase the maximum deformation of the end-effector by considering the external forces along different directions, thereby increasing the capacity of the end-effector to withstand external forces.Then, the optimization function of the second objective along the vertical direction that is needed to increase the capacity of the end-effector to withstand external forces can be expressed as follows: where δ ∈ R 3×1 is the specified optimized direction vector for the enhanced force capacity, and δ = [0, 0, 1] T .According to (35), the self-gravity of the end-effector affects its capacity to withstand external forces in the vertical direction.Therefore, to optimize the gravity-induced deformation, it is reasonable to set the third objective function as follows: Thus, we can decrease the effect of G(q m ) in ( 32) by reducing the objective function F 3 in (37), meaning that we can reduce the joint torques induced by G(q m ) by reducing ∆x g .
During optimization, the three functions in ( 29), (36), and (37) should have the same order of magnitude; then, the overall optimization function can be expressed as follows: where w 1 , w 2 , and w 3 are the weights of the three optimization functions; α 1 , α 2 , and α 3 are the normalized factors; τ jmin , τ jmax are the minimum and maximum torque of the j-th joint; and q jmin and q jmax are the minimum and maximum angle of the j-th joint.
The velocity tasks defined by the null space can be expressed as follows: where λ is the control gain and ∇ q m F is the gradient of F(q m ) in (38).By substituting (39) in (27), we obtain the following null-space controller:

Position and Posture Controller for the End-Effector
During assisted walking, the position and posture of the end-effector should be maintained at their optimal configuration.Therefore, we apply a proportional controller to the end-effector.
The quaternion in ideal position and posture of the end-effector Q de can be obtained from w e R in w e T. The desired end-effector position in Σ b is set as x e .Assuming that the current end-effector position is x c , the error between x e and x c is x err = x e − x c .Therefore, the desired linear velocity v de of the end-effector can be expressed as follows: where K v is the proportional coefficient.
The current posture of the end-effector can be expressed as Q c = w c + u c , and the error between Q de and Q c can be written as: where Then, the angle errors can be obtained from Q err , as follows: Therefore, the desired angular velocity of the end-effector can be obtained.
where K w is the proportional coefficient.
The controller for the end-effector can be summarized as Then, v e is projected onto the mobile platform and the manipulator, separately, by using J u (q) in (3): where u be and u me are the velocity of the mobile platform and manipulator.

Integrated Controller for Assisted Walking
From ( 40) and (46), the velocity relationship between the mobile platform and the manipulator can be inferred clearly.
(1) In the optimization process, there is no interaction between the human user and the mobile manipulator during optimization.Therefore, the velocities of the mobile platform and manipulator can be computed as follows: (2) During assisted walking, the velocities of the mobile platform and manipulator can be obtained using the total velocity controller expressed as (47): The assisted walking system comprising the mobile manipulator is developed, as depicted in Figure 5.In the following section of null-space optimization, the desired velocity of the end-effector can be obtained by (47), and the mobile platform and the manipulator move together.After null-space optimization, the end-effector maintains its posture by the proposed controller in (45), and therefore, the desired velocity of the end-effector in the horizontal direction is allocated to the mobile platform while that in the vertical direction is set to be 0 to provide adequate assisted walking support to the user.The mobile manipulator walking-assistance technology based on the proposed velocity feedforward compensator in (21) and variable admittance controller in (25) can be implemented.It is worth noting that the configuration optimization is carried out first by the objective function (38) and the null-space controller (40), and the human-robot interaction begins after the optimization is completed.Then, the optimized configuration can be maintained by (40) in the process of assisted walking.

Experimental Verification
In order to measure the compliance, stability, and comfort of the human-robot interaction, we will collect the velocity, the interaction force, the trajectory of the end-effector, and the joint angle of the manipulator in the interaction process shown in Video S1.Due to the limitations of the laboratory conditions, we selected three volunteers to test the different requirements of walking assistance.

Experimental Setup
The experimental platform consists of a UR3e manipulator (Universal Robots Inc., Aarhus, Denmark) and a four-wheeled mobile robot named SCOUT (AgileX Inc., Dongguan, China), as shown in Figure 6, and the parameters are summarized in Table 2. Herein, L, W, and H are the length, width, and height of the SCOUT robot; l0, w0, and h0 are the installation positions of the manipulator; and l1 and l2 are the wheelbase and wheeltrack.The proposed assisted walking system is implemented in ROS (Robot Operation System) Kinetic, while the mobile robot owns an odometer, the manipulator owns a sixaxis force sensor at the end-effector, and the motion for its joints can be obtained in real time.The global controller for the mobile robot and the manipulator are integrated on an upper computer based on Figure 5, which communicates with the lower control units for the mobile robot and the manipulator through CAN (Controller Area Network) and LAN (Local Area Network), respectively.The calculation frequency of the whole system is 500 Hz.The six-axis force sensor shown in Figure 6 is employed to obtain the interacting forces for the assisted walking system shown in Figure 5.The experiments can be seen in Video S1.

Experimental Verification
In order to measure the compliance, stability, and comfort of the human-robot interaction, we will collect the velocity, the interaction force, the trajectory of the end-effector, and the joint angle of the manipulator in the interaction process shown in Video S1.Due to the limitations of the laboratory conditions, we selected three volunteers to test the different requirements of walking assistance.

Experimental Setup
The experimental platform consists of a UR3e manipulator (Universal Robots Inc., Aarhus, Denmark) and a four-wheeled mobile robot named SCOUT (AgileX Inc., Dongguan, China), as shown in Figure 6, and the parameters are summarized in Table 2. Herein, L, W, and H are the length, width, and height of the SCOUT robot; l 0 , w 0 , and h 0 are the installation positions of the manipulator; and l 1 and l 2 are the wheelbase and wheel-track.The proposed assisted walking system is implemented in ROS (Robot Operation System) Kinetic, while the mobile robot owns an odometer, the manipulator owns a six-axis force sensor at the end-effector, and the motion for its joints can be obtained in real time.The global controller for the mobile robot and the manipulator are integrated on an upper computer based on Figure 5, which communicates with the lower control units for the mobile robot and the manipulator through CAN (Controller Area Network) and LAN (Local Area Network), respectively.The calculation frequency of the whole system is 500 Hz.The six-axis force sensor shown in Figure 6 is employed to obtain the interacting forces for the assisted walking system shown in Figure 5.The experiments can be seen in Video S1.First, a human-machine flexible interaction experiment and the following null-space control experiment are conducted.The effectiveness of the proposed controllers in (38) and ( 40) is reflected by the optimization functions F1, F2, F3, and F and the joint torques of the manipulator before and after optimization.
Second, a compliant interactive experiment is conducted to verify the validity of the proposed state observer (14), the velocity feedforward compensator (21), and the variable admittance controller (25) based on the fuzzy theory.The effectiveness of the above control algorithm is reflected by the real and predictive values of joint torques, the interaction forces, and the velocities of the end-effector and mobile platform in the process of compliant interaction.
Finally, different volunteers conducted the experiments to verify the stability and robustness of the assisted walking technology based on the compliant interaction approach proposed herein.The effectiveness of the approach is reflected in the interaction forces reported by the various volunteers and the end-effector velocities recorded during assisted walking.

Human-Machine Flexible Interaction Experiment
In this paper, the proposed human-robot flexible interaction system is verified by experiments.The specific test is to let the operator hold the end-effector, and then the operator pulls the manipulator to carry out a series of complex actions.The traditional admittance control (Method I [41]), admittance control (Method II [42]) with fuzzy adjustment of virtual damping, admittance control (Method III [31]) with fuzzy adjustment of virtual damping, and adjustment of virtual mass are tested together with the flexible  First, a human-machine flexible interaction experiment and the following null-space control experiment are conducted.The effectiveness of the proposed controllers in (38) and ( 40) is reflected by the optimization functions F 1 , F 2 , F 3 , and F and the joint torques of the manipulator before and after optimization.
Second, a compliant interactive experiment is conducted to verify the validity of the proposed state observer (14), the velocity feedforward compensator (21), and the variable admittance controller (25) based on the fuzzy theory.The effectiveness of the above control algorithm is reflected by the real and predictive values of joint torques, the interaction forces, and the velocities of the end-effector and mobile platform in the process of compliant interaction.
Finally, different volunteers conducted the experiments to verify the stability and robustness of the assisted walking technology based on the compliant interaction approach proposed herein.The effectiveness of the approach is reflected in the interaction forces reported by the various volunteers and the end-effector velocities recorded during assisted walking.

Human-Machine Flexible Interaction Experiment
In this paper, the proposed human-robot flexible interaction system is verified by experiments.The specific test is to let the operator hold the end-effector, and then the operator pulls the manipulator to carry out a series of complex actions.The traditional admittance control (Method I [41]), admittance control (Method II [42]) with fuzzy adjustment of virtual damping, admittance control (Method III [31]) with fuzzy adjustment of virtual damping, and adjustment of virtual mass are tested together with the flexible interactive control method (Method IV) proposed in this paper, and among them, Methods I, II, and III do not have feedforward velocity control.The parameters of the fuzzy inference system and virtual mass control law used in Methods II and III are consistent with those of Method IV.The experimental results are shown in Figures 7-9   As shown in Table 3, the interaction forces in Method IV are much smaller, both at maximum and average, than in the previous three methods.Although the average value of the interaction force in Method II is slightly higher than that in Method I, the velocity of the end-effector is also higher than that in Method I, and it is more obvious.In addition, interactive control method (Method IV) proposed in this paper, and among them, Methods I, II, and III do not have feedforward velocity control.The parameters of the fuzzy inference system and virtual mass control law used in Methods II and III are consistent with those of Method IV.The experimental results are shown in Figures 7-9 and Table 1.As shown in Table 3, the interaction forces in Method IV are much smaller, both at maximum and average, than in the previous three methods.Although the average value of the interaction force in Method II is slightly higher than that in Method I, the velocity of the end-effector is also higher than that in Method I, and it is more obvious.In addition, interactive control method (Method IV) proposed in this paper, and among them, Methods I, II, and III do not have feedforward velocity control.The parameters of the fuzzy inference system and virtual mass control law used in Methods II and III are consistent with those of Method IV.The experimental results are shown in Figures 7-9 and Table 1.As shown in Table 3, the interaction forces in Method IV are much smaller, both at maximum and average, than in the previous three methods.Although the average value of the interaction force in Method II is slightly higher than that in Method I, the velocity of the end-effector is also higher than that in Method I, and it is more obvious.In addition, As shown in Table 3, the interaction forces in Method IV are much smaller, both at maximum and average, than in the previous three methods.Although the average value of the interaction force in Method II is slightly higher than that in Method I, the velocity of the end-effector is also higher than that in Method I, and it is more obvious.In addition, the maximum interaction force in Method II is smaller than that in Method I, while the maximum speed of the end-effector is larger.It can be seen that the human motion intention prediction method proposed in this paper can improve the compliance of the human-robot interaction system.The maximum value and average value of interaction force in Method III are smaller than those in Method II, and the end-effector movement speed is faster, which proves that maintaining the dynamic response characteristics of the system can further improve the compliance of the interactive system.In addition, although the maximum and average velocity of the end-effector in Method IV are 4.62% and 1.75% lower than that in Method III, the maximum and average interaction force are 5.83% and 11.34% lower than that in Method III, respectively.Method IV can make the manipulator end effector produce a speed similar to Method III with less force due to the use of the feedforward velocity control.In summary, the Method IV proposed in this paper has good compliance and comfort among the four methods.

Null-Space Control Experiment
To verify whether the proposed objective function (38) can find an optimized configuration for the mobile manipulator while the end-effector is maintained at a given point.
Figure 9 shows the configurations of the mobile manipulator before and after optimization.According to Figure 10, in the initial posture, the three objective functions are close.As the optimization proceeds, F 2 in (36) and F 1 in (29) exhibit two opposite trends, indicating that enhancement of the force capacity of the end-effector and maintenance of distance from the singular configuration cannot be achieved simultaneously.However, the F 3 in (37) decreases rapidly and then increases marginally to reach a small steady-state value, indicating that the gravity effect of the manipulator is effectively suppressed, and the overall objective function F (38) gradually converges to a stable value.The minimum value of F1 is greater than 0, indicating it can avoid singularity.The minimum value of F3 is smaller than its stable value, indicating that if the force capacity of the end-effector in the vertical direction is improved, the influence of gravity can be decreased, but it cannot be eliminated.The minimum value of F 1 is greater than 0, indicating it can avoid singularity.The minimum value of F 3 is smaller than its stable value, indicating that if the force capacity of the end-effector in the vertical direction is improved, the influence of gravity can be decreased, but it cannot be eliminated.
As shown in Figure 11, the pose of the end-effector may deviate marginally from the desired pose.According to Figure 11a, the end-effector has different degrees of deviation in three directions, of which the position error in the x-direction is the greatest, and its maximum value is 1.492 cm.The obvious deviation in the x-direction may have been caused by the mobile platform, which cannot implement the desired small velocity commands.The end posture errors shown in Figure 11b are within an extremely small range (approximately 2 × 10 −4 rad).Therefore, it can be concluded that the mobile manipulator can well maintain the desired pose of the end-effector during the optimization, indicating that the null-space approach to finding an optimized configuration is effective.The minimum value of F1 is greater than 0, indicating it can avoid singularity.The minimum value of F3 is smaller than its stable value, indicating that if the force capacity of the end-effector in the vertical direction is improved, the influence of gravity can be decreased, but it cannot be eliminated.
As shown in Figure 11, the pose of the end-effector may deviate marginally from the desired pose.According to Figure 11a, the end-effector has different degrees of deviation in three directions, of which the position error in the x-direction is the greatest, and its maximum value is 1.492 cm.The obvious deviation in the x-direction may have been caused by the mobile platform, which cannot implement the desired small velocity commands.The end posture errors shown in Figure 11b are within an extremely small range (approximately 2 × 10 −4 rad).Therefore, it can be concluded that the mobile manipulator can well maintain the desired pose of the end-effector during the optimization, indicating that the null-space approach to finding an optimized configuration is effective.To verify whether the proposed null-space controller (40) can improve the force capacity of the end-effector in the vertical direction while maintaining the end-effector pose, two comparison tests are conducted by applying the same load to the end-effector before and after optimization.An end gripper is used to grab a bottle of water (approximately 550 mL), meaning that a vertical load of approximately 6 N is applied to the end-effector, as shown in Figure 12.To verify whether the proposed null-space controller (40) can improve the force capacity of the end-effector in the vertical direction while maintaining the end-effector pose, two comparison tests are conducted by applying the same load to the end-effector before and after optimization.An end gripper is used to grab a bottle of water (approximately 550 mL), meaning that a vertical load of approximately 6 N is applied to the end-effector, as shown in Figure 12.As shown in Figure 13, before the optimization, the heavy load is grabbed by the manipulator gripper at approximately 35 s, and after the optimization, it is grabbed at approximately 27 s.Regardless of whether the end-effector is subjected to an external force, after optimization, the torques of joints 2 and 3 are significantly smaller than those before optimization, meaning that the proposed objective functions F2 (36) and F3 (37) are effective.Moreover, the torques of the other joints are always small because these joints are not used to balance the vertical load against gravity.The torque details of the joints are summarized in Table 4.As shown in Figure 13, before the optimization, the heavy load is grabbed by the manipulator gripper at approximately 35 s, and after the optimization, it is grabbed at approximately 27 s.Regardless of whether the end-effector is subjected to an external force, after optimization, the torques of joints 2 and 3 are significantly smaller than those before optimization, meaning that the proposed objective functions F 2 (36) and F 3 (37) are effective.Moreover, the torques of the other joints are always small because these joints are not used to balance the vertical load against gravity.The torque details of the joints are summarized in Table 4.
As shown in Figure 13, before the optimization, the heavy load is grabbed by the manipulator gripper at approximately 35 s, and after the optimization, it is grabbed at approximately 27 s.Regardless of whether the end-effector is subjected to an external force, after optimization, the torques of joints 2 and 3 are significantly smaller than those before optimization, meaning that the proposed objective functions F2 (36) and F3 (37) are effective.Moreover, the torques of the other joints are always small because these joints are not used to balance the vertical load against gravity.The torque details of the joints are summarized in Table 4.According to Table 4, when the manipulator grasps the load and remains static, the torques of joints 2, 3, and 4 vary substantially, and the 2-norm of the torque vector increases accordingly.The torques of joints 2, 3, and 4 are mainly used to balance the vertical load against gravity.In the absence of external load, the torques of joints 2, 3, and 4 decrease by 50.78%, 28.11%, and 55.02%, respectively, after optimization, meaning that the gravity effect of the manipulator itself is effectively weakened.Under a load, the torques of joints 2 and 3 are 63.07%and 14.18% smaller, respectively, than those before optimization, and the torque of joint 4 is close to that before optimization.In general, the 2-norm of the torque vector after optimization is 44.34% and 41.77% smaller than that before optimization, respectively.This indicates that the null-space controller (40) can significantly improve the force capacity of the end-effector in the vertical direction.According to Table 4, when the manipulator grasps the load and remains static, the torques of joints 2, 3, and 4 vary substantially, and the 2-norm of the torque vector increases accordingly.The torques of joints 2, 3, and 4 are mainly used to balance the vertical load against gravity.In the absence of external load, the torques of joints 2, 3, and 4 decrease by 50.78%, 28.11%, and 55.02%, respectively, after optimization, meaning that the gravity effect of the manipulator itself is effectively weakened.Under a load, the torques of joints 2 and 3 are 63.07%and 14.18% smaller, respectively, than those before optimization, and the torque of joint 4 is close to that before optimization.In general, the 2-norm of the torque vector after optimization is 44.34% and 41.77% smaller than that before optimization, respectively.This indicates that the null-space controller (40) can significantly improve the force capacity of the end-effector in the vertical direction.

Assisted Walking Test
In the assisted walking test, the manipulator position is set to q1 = [0, −2π/3, −π/2, −π/3, −π/2, 0] T , and the end-effector position is x1 = [0.713,−0.131, 0.348] T .First, the nullspace controller is used to ensure that the mobile manipulator reaches the optimal configuration.Then, the manipulator position is q2 = [0, −89π/200, −387π/500, −281π/1000, −π/2, 0] T , and the end-effector position is x2 = [0.698,−0.133, 0.348] T .Next, the experimenter's arm is clamped using the gripper, and compliant force interaction is performed along the x-, y-, and z-directions of the mobile manipulator.The experimental process is shown in Figure 14. Figure 15 shows the trajectory of the end-effector in the world coordinate system during assisted walking.The trajectory is smooth and consistent with the human motion intention.A complete circle is drawn in the horizontal plane, which verifies that the compliant interaction controller proposed in this paper can not only perceive and execute the Figure 15 shows the trajectory of the end-effector in the world coordinate system during assisted walking.The trajectory is smooth and consistent with the human motion intention.A complete circle is drawn in the horizontal plane, which verifies that the compliant interaction controller proposed in this paper can not only perceive and execute the human user's forward or backward motion intention but also comply with the human user's steering intention.Figure 15 shows the trajectory of the end-effector in the world coordinate system during assisted walking.The trajectory is smooth and consistent with the human motion intention.A complete circle is drawn in the horizontal plane, which verifies that the compliant interaction controller proposed in this paper can not only perceive and execute the human user's forward or backward motion intention but also comply with the human user's steering intention.Figure 16 presents a comparison between the real (values observed using sensors) and predicted (values obtained using the state observer ( 14)) values of the joint torques induced by the human arm on the end-effector during assisted walking.Herein, the proposed velocity feedforward compensator ( 21) is only applied to the first three joints, and the root mean square errors of joints 1, 2, and 3 are 0.244, 0.171, and 0.154, respectively, and the mean square errors of the latter three joints are 0.236, 0.267, and 0.178, respectively.In conclusion, the errors between the predicted and real values are small for the first three joints, and none of the predicted values are greater than the corresponding real values, thereby guaranteeing that the manipulator moves within the allowable and safe range.Figure 16 presents a comparison between the real (values observed using sensors) and predicted (values obtained using the state observer ( 14)) values of the joint torques induced by the human arm on the end-effector during assisted walking.Herein, the proposed velocity feedforward compensator ( 21) is only applied to the first three joints, and the root mean square errors of joints 1, 2, and 3 are 0.244, 0.171, and 0.154, respectively, and the mean square errors of the latter three joints are 0.236, 0.267, and 0.178, respectively.In conclusion, the errors between the predicted and real values are small for the first three joints, and none of the predicted values are greater than the corresponding real values, thereby guaranteeing that the manipulator moves within the allowable and safe range.
These findings validate the effectiveness of the proposed velocity feedforward controller (21) and the proposed state observer (14).Moreover, the predictions for joint torque are accurate and obtained in real time, and they can be used to further improve the performance of motion-tracking and the compliance between the human and the robot.
Figures 17 and 18 present the human-robot interaction forces and end-effector velocity during assisted walking.The interaction forces in the x-and y-directions are less than 15 N, and the end-effector velocity can even exceed 0.3 m/s, meaning that the compliant interaction method proposed herein offers good compliance in the horizontal direction.At 65 s, a vertical load (of approximately 30 N) is applied to the end-effector, resulting in a small vertical velocity (about 0.02 m/s), which validates that the mobile manipulator equipped with the designed controller can provide reliable support to the human user's hand.
As shown in Figure 19, when the end-effector is subjected to an interaction force in the x-direction (from 5 s to 22 s), the mobile platform tracks the desired velocity u b in (48), which is forward or backward.The mobile platform then rotates to track the desired lateral motion or turning motion when the end-effector is subjected to an interaction force in the y-direction (from 22 s to 65 s).This is because the mobile platform used in this study is constrained by non-holonomic kinematics.
As shown in Figure 20, the joint angles of the manipulator vary in a small range because the proposed approach can maintain the optimal configuration by the controllers (40) and (45).This guarantees that the robot always has good force capacity and stability.
Figure 21 shows the position errors of the end-effector along the x, y, and z axes.The maximum position error along the x-direction (approximately 1.5 cm) is greater than that along the y-direction (approximately 0.2 cm) and z-direction (approximately 0.7 cm) from 5 s to 22 s.This is because the maximum interaction force along the x-direction (approximately 12 N) is greater than that along the y-direction (approximately 5 N), while a stiffness coefficient is applied along the z-direction in (25).From 22 s to 65 s, the proposed approach performs similarly in the y-direction.At 65 s, the end-effector is subjected to a large vertical load (approximately 30 N), and the displacement along the z-direction becomes approximately 1.7 cm.The displacement in the z-direction can generate an elastic force that can support the load in the vertical direction.However, the deviation varies within a small range (0-2 cm), which verifies the effectiveness of the proposed position and posture controller for the end-effector (45) in terms of controlling its position.These findings validate the effectiveness of the proposed velocity feedforward controller (21) and the proposed state observer (14).Moreover, the predictions for joint torque are accurate and obtained in real time, and they can be used to further improve the performance of motion-tracking and the compliance between the human and the robot.
Figures 17 and 18 present the human-robot interaction forces and end-effector velocity during assisted walking.The interaction forces in the x-and y-directions are less than 15 N, and the end-effector velocity can even exceed 0.3 m/s, meaning that the compliant interaction method proposed herein offers good compliance in the horizontal direction.At 65 s, a vertical load (of approximately 30 N) is applied to the end-effector, resulting in a small vertical velocity (about 0.02 m/s), which validates that the mobile manipulator equipped with the designed controller can provide reliable support to the human user's hand.These findings validate the effectiveness of the proposed velocity feedforward controller ( 21) and the proposed state observer (14).Moreover, the predictions for joint torque are accurate and obtained in real time, and they can be used to further improve the performance of motion-tracking and the compliance between the human and the robot.
Figures 17 and 18 present the human-robot interaction forces and end-effector velocity during assisted walking.The interaction forces in the x-and y-directions are less than 15 N, and the end-effector velocity can even exceed 0.3 m/s, meaning that the compliant interaction method proposed herein offers good compliance in the horizontal direction.At 65 s, a vertical load (of approximately 30 N) is applied to the end-effector, resulting in a small vertical velocity (about 0.02 m/s), which validates that the mobile manipulator equipped with the designed controller can provide reliable support to the human user's hand.As shown in Figure 19, when the end-effector is subjected to an interaction force in the x-direction (from 5 s to 22 s), the mobile platform tracks the desired velocity ub in (48), which is forward or backward.The mobile platform then rotates to track the desired lateral motion or turning motion when the end-effector is subjected to an interaction force in the y-direction (from 22 s to 65 s).This is because the mobile platform used in this study is constrained by non-holonomic kinematics.As shown in Figure 20, the joint angles of the manipulator vary in a small range because the proposed approach can maintain the optimal configuration by the controllers (40) and (45).This guarantees that the robot always has good force capacity and stability.As shown in Figure 19, when the end-effector is subjected to an interaction force in the x-direction (from 5 s to 22 s), the mobile platform tracks the desired velocity ub in (48), which is forward or backward.The mobile platform then rotates to track the desired lateral motion or turning motion when the end-effector is subjected to an interaction force in the y-direction (from 22 s to 65 s).This is because the mobile platform used in this study is constrained by non-holonomic kinematics.As shown in Figure 20, the joint angles of the manipulator vary in a small range because the proposed approach can maintain the optimal configuration by the controllers (40) and (45).This guarantees that the robot always has good force capacity and stability.As shown in Figure 19, when the end-effector is subjected to an interaction force in the x-direction (from 5 s to 22 s), the mobile platform tracks the desired velocity ub in (48), which is forward or backward.The mobile platform then rotates to track the desired lateral motion or turning motion when the end-effector is subjected to an interaction force in the y-direction (from 22 s to 65 s).This is because the mobile platform used in this study is constrained by non-holonomic kinematics.As shown in Figure 20, the joint angles of the manipulator vary in a small range because the proposed approach can maintain the optimal configuration by the controllers ( 40) and ( 45).This guarantees that the robot always has good force capacity and stability.Figure 21 shows the position errors of the end-effector along the x, y, and z axes.The maximum position error along the x-direction (approximately 1.5 cm) is greater than that along the y-direction (approximately 0.2 cm) and z-direction (approximately 0.7 cm) from 5 s to 22 s.This is because the maximum interaction force along the x-direction (approximately 12 N) is greater than that along the y-direction (approximately 5 N), while a stiffness coefficient is applied along the z-direction in (25).From 22 s to 65 s, the proposed approach performs similarly in the y-direction.At 65 s, the end-effector is subjected to a large vertical load (approximately 30 N), and the displacement along the z-direction becomes approximately 1.7 cm.The displacement in the z-direction can generate an elastic force that can support the load in the vertical direction.However, the deviation varies within a small range (0-2 cm), which verifies the effectiveness of the proposed position and posture controller for the end-effector (45) in terms of controlling its position.Figure 22 presents the posture errors of the end-effector relative to the mobile platform.The angle error around the z axis is approximately 0.04 rad, while others are about 0. This is because when the mobile manipulator is subjected to the interaction force in the y-direction, the mobile platform rotates around the z axis.Due to this force, an obvious  The angle error around the z axis is approximately 0.04 rad, while others are about 0. This is because when the mobile manipulator is subjected to the interaction force in the y-direction, the mobile platform rotates around the z axis.Due to this force, an obvious error appears around the z-axis.However, the posture errors vary within a small range (0-0.04 rad), which verifies the validity of the proposed position and posture controller for the end-effector (45) in terms of controlling its posture.Figure 22 presents the posture errors of the end-effector relative to the mobile platform.The angle error around the z axis is approximately 0.04 rad, while others are about 0. This is because when the mobile manipulator is subjected to the interaction force in the y-direction, the mobile platform rotates around the z axis.Due to this force, an obvious error appears around the z-axis.However, the posture errors vary within a small range (0-0.04 rad), which verifies the validity of the proposed position and posture controller for the end-effector (45) in terms of controlling its posture.

Assisted Walking Tests with Three Volunteers
To further verify the adaptability of the assisted walking approach proposed herein, comparison tests are conducted.Herein, we select three volunteers who have different

Assisted Walking Tests with Three Volunteers
To further verify the adaptability of the assisted walking approach proposed herein, comparison tests are conducted.Herein, we select three volunteers who have different physical characteristics and can represent different groups of people.These volunteers are A (male), B (female), and C (male).Their heights and weights are 185 cm, 165 cm, and 175 cm, and 75 kg, 45 kg, and 60 kg, respectively.In the following experiments, the reference path is similar to a sine curve, as shown in Figure 23.This path is only a guideline for the volunteers, and it is not strictly tracked by humans.The experimental results are presented in Figures 24-26, and the performances of the three volunteers are comparatively summarized in Table 5.Moreover, we interviewed three volunteers to evaluate the comfort, stability, and compliance of the assisted walking process, and the score range was 0-10 points, as shown in Table 6.  5.Moreover, we interviewed three volunteers to evaluate the comfort, stability, and compliance of the assisted walking process, and the score range was 0-10 points, as shown in Table 6.Through analyzing the experimental data of three different volunteers in Figures 24-26, the velocity of the end-effector all presents the same change trend with the interaction force.In addition, by analyzing the data in Table 5, it can be seen that the absolute maximum force between three different volunteers and the mobile manipulator in the x-direction is 16.82 N, 15.48 N, and 13.51 N, respectively.The maximum absolute values of the endeffector velocity in the x-direction are 0.45 m/s, 0.37 m/s, and 0.32 m/s, respectively.It can be concluded that the bigger the interaction force is in a certain direction of the end-effector, the bigger the velocity of the end-effector is in that direction.It can be seen that the proposed controller (48) can produce a corresponding speed according to the interaction forces exerted by different users and obediently track the movement of the users.It indicates that the proposed assisted walking controller (48) is robust for users with different gaits.
In the vertical direction, the maximum absolute values of the interaction forces are 24.41N, 27.32 N, and 28.69 N, respectively, and the corresponding maximum absolute values of the end-effector velocities in the z-direction are 0.0051 m/s, 0.0093 m/s, and 0.0096 m/s, respectively.Moreover, the velocities vary within a small range, as shown in Figures 24b, 25b and 26b.Therefore, during assisted walking, even if the end-effector bears a large vertical load, its deformation is small, meaning that the robot can provide adequate support force to the human user in the vertical direction.
In Table 6, compared with the other two volunteers, volunteer C had the best experience during the assisted walking process and gave a high evaluation.Volunteer B, the only female, only had a high evaluation on stability.Volunteer A, who is the highest among the three volunteers, also had a good evaluation of the compliance, stability, and compliance of the assisted walking system.Based on their feedback, we can know that the assisted walking system has good compliance, stability, and robustness for different humans.However, they all mentioned that the comfort and compliance of the walking assistance system when turning is not as good as walking in a straight line, and we will further improve the comfort and compliance of assisted walking when turning in the future.
The results of the three abovementioned experiments involving different volunteers indicate that the proposed assisted walking approach with the mobile manipulator can accurately identify human motion intentions and implement compliant tracking of human motion across individuals.In conclusion, the proposed approach has good stability, high accuracy, and a high level of robustness.Although the amplitudes of the interaction forces shown in  are small enough, this may not be suitable for a small number of people, and we will make further research on the amplitude control of the interactive forces.

Conclusions
In this paper, a novel approach to human-robot-compliant interaction for mobile manipulators is proposed.This approach can be used to assist human walking.To realize compliant interaction between a human user and the mobile manipulator, a method for predicting the human user's motion intention based on the fuzzy theory is proposed, and it is used to regulate the mass and damping coefficients.Moreover, a velocity feedforward compensator with a state observer is used to reduce the resistance caused by the inertia of the mobile manipulator.Furthermore, a null-space controller based on the projection gradient method is proposed to maintain the optimal manipulator configuration during assisted walking.The experimental results obtained through experiments involving different volunteers indicated that the proposed approach has good stability, high accuracy, and a high level of robustness.
At present, the aging society requires more assisted walking equipment for rehabilitation and daily nursing.As a multi-functional robot, the mobile manipulator can not only meet other daily services but also assist humans in walking, according to the proposed approach in this paper.In the future, we will seek cooperation from different kinds of people to extend the application of the proposed assisted walking system in different environments, including elderly people and hospital patients, and further research on the anti-fall safety strategies and improve its performance, which can make it safer and more intelligent and can be widely applied for elderly or disabled persons.

Figure 1 .
Figure 1.Schematic diagram of mobile manipulator system.

Figure 1 .
Figure 1.Schematic diagram of mobile manipulator system.

Figure 3 .
Figure 3. Membership function of interaction force and end-effector velocity.Figure 3. Membership function of interaction force and end-effector velocity.

Figure 3 .Figure 4 .
Figure 3. Membership function of interaction force and end-effector velocity.Figure 3. Membership function of interaction force and end-effector velocity.Biomimetics 2024, 9, x FOR PEER REVIEW 8 of 30

Figure 4 .
Figure 4. Membership function of virtual damping.

Figure 5 .
Figure 5. Assisted walking system with mobile manipulator.

Figure 5 .
Figure 5. Assisted walking system with mobile manipulator.

Figure 9 .
Figure 9.Comparison of robot configuration before and after optimization.(a) Before optimization; (b) After optimization.

Figure 11 .
Figure 11.Pose errors of the end-effector during optimization.(a) Position errors; (b) Posture errors.

Figure 11 .
Figure 11.Pose errors of the end-effector during optimization.(a) Position errors; (b) Posture errors.

Figure 12 .
Figure 12.Comparison of force capacity before and after optimization.(a) Before optimization; (b) After optimization.

Figure 12 .
Figure 12.Comparison of force capacity before and after optimization.(a) Before optimization; (b) After optimization.

Figure 16 .
Figure 16.Real and predicted values of joint torque.

Figure 16 .Figure 16 .
Figure 16.Real and predicted values of joint torque.

Figure 17 .
Figure 17.Interaction forces of end-effector (from 65 s, a vertical load (of approximately 30 N) is applied to the end-effector, resulting in a small vertical velocity (about 0.02 m/s in Figure18)).

Figure 17 .
Figure17.Interaction forces of end-effector (from 65 s, a vertical load (of approximately 30 N) is applied to the end-effector, resulting in a small vertical velocity (about 0.02 m/s in Figure18)).

Figure 20 .
Figure 20.Angle of each joint.Figure 20.Angle of each joint.

Figure 20 .
Figure 20.Angle of each joint.Figure 20.Angle of each joint.

Figure 21 .
Figure 21.Position errors of the end-effector relative to the mobile platform.

Figure 21 .
Figure 21.Position errors of the end-effector relative to the mobile platform.

Figure 22
Figure22presents the posture errors of the end-effector relative to the mobile platform.The angle error around the z axis is approximately 0.04 rad, while others are about 0. This is because when the mobile manipulator is subjected to the interaction force in the y-direction, the mobile platform rotates around the z axis.Due to this force, an obvious error appears around the z-axis.However, the posture errors vary within a small range (0-0.04 rad), which verifies the validity of the proposed position and posture controller for the end-effector (45) in terms of controlling its posture.

Figure 21 .
Figure 21.Position errors of the end-effector relative to the mobile platform.

Figure 22 .
Figure 22.Posture errors of the end-effector relative to the mobile platform.

Figure 22 .
Figure 22.Posture errors of the end-effector relative to the mobile platform.

Biomimetics 2024, 9 ,
x FOR PEER REVIEW 26 of 30 physical characteristics and can represent different groups of people.These volunteers are A (male), B (female), and C (male).Their heights and weights are 185 cm, 165 cm, and 175 cm, and 75 kg, 45 kg, and 60 kg, respectively.In the following experiments, the reference path is similar to a sine curve, as shown in Figure23.This path is only a guideline for the volunteers, and it is not strictly tracked by humans.The experimental results are presented in Figures24-26, and the performances of the three volunteers are comparatively summarized in Table

Figure 24 .
Figure 24.Test results of volunteer A. (a) Interactive forces; (b) End effector velocity.

Figure 24 .
Figure 24.Test results of volunteer A. (a) Interactive forces; (b) End effector velocity.

Figure 24 .
Figure 24.Test results of volunteer A. (a) Interactive forces; (b) End effector velocity.

Figure 25 .
Figure 25.Test results of volunteer B. (a) Interactive forces; (b) End effector velocity.Figure 25.Test results of volunteer B. (a) Interactive forces; (b) End effector velocity.

Figure 26 .
Figure 26.Test results of volunteer C. (a) Interactive forces; (b) End effector velocity.

Table 2 .
Control algorithm parameters and geometric parameters of moving manipulator.

Table 2 .
Control algorithm parameters and geometric parameters of moving manipulator.

Table 1 .
and Biomimetics 2024, 9, x FOR PEER REVIEW 18 of 30 interactive control method (Method IV) proposed in this paper, and among them, Methods I, II, and III do not have feedforward velocity control.The parameters of the fuzzy inference system and virtual mass control law used in Methods II and III are consistent with those of Method IV.The experimental results are shown in Figures 7-9 and Table 1.

Table 3 .
Specific statistics of the interaction force and end effector velocity.

Table 4 .
Specific joint torque data.

Table 4 .
Specific joint torque data.

Table 5 .
Performance comparison of volunteers A, B, and C.

Table 6 .
Evaluation of volunteers A, B, and C.

Table 5 .
Performance comparison of volunteers A, B, and C.

Table 6 .
Evaluation of volunteers A, B, and C.