Active Loading Control Design for a Wearable Exoskeleton with a Bowden Cable for Transmission

: Exoskeletons with a Bowden cable for power transmission have the advantages of a concentrated mass and ﬂexible movement. However, their integrated motor is disturbed by the Bowden cable’s friction, which limits the performance of the force loading response. In this paper, we solve this problem by designing an outer-loop feedforward-feedback proportion-differentiation controller based on an inner loop disturbance observer. Firstly, the inner loop’s dynamic performance is equivalent to the designed nominal model using the proposed disturbance observer, which effectively compensates for the parameter perturbation and friction disturbance. Secondly, based on an analysis of the stability of the inner loop controller, we obtain the stability condition and discuss the inﬂuence of modeling errors on the inner loop’s dynamic performance. Thirdly, to avoid excessive noise from the force sensors being introduced into the designed disturbance observer, we propose the feedforward-feedback proportion-differentiation controller based on the nominal model and pole conﬁguration, which improves the outer loop’s force loading performance. Experiments are conducted, which verify the effectiveness of the proposed methods. output shaft and the knee joint sheave. Two lightweight carbon ﬁberboards (LCFB) form a rotating joint, and the hinge mechanism forms a rotating joint on the coronal plane of the knee joint. The bafﬂes are connected to the thigh and calve with a few ﬂexible straps. Since the assistance needs to be performed based on the human body’s active movement, an encoder at the knee joint is installed. Moreover, a force loading control framework based on a speed loop is proposed. loop nominal experimental results show phase


Introduction
Lower-limb exoskeletons can be divided into three types, according to their different auxiliary tasks. The first one enhances human strength and endurance [1,2]. The second focuses on assisting patients or elderly people with a partial decline in their motion performance [3,4]. The third is an exoskeleton [5,6] that is suitable for the rehabilitation training of patients with a near-loss walking ability. To reduce the impact of heavy loads on humans, weight-bearing is usually transferred to the ground [2]. The body balance of patients also needs to be supported by the rehabilitation of their exoskeleton [5]. The structures of these exoskeletons are usually rigid and run parallel to the limbs. Traditional rigid exoskeletons facilitate sensor integration and installation, and a dynamic model of these exoskeletons [7,8] has been established to realize zero-moment. However, for people with a normal motion capacity, the flexibility and comfort of rigid exoskeletons are still limited.
A wearable soft exosuit has been proposed [9], in which a Bowden cable is used for transmission. A Bowden cable can transmit displacement through the relative movement of an internal cable relative to an external sheath, which has the advantages of a flexible layout and long-distance transmission. Therefore, a wearable soft exosuit significantly enhances movement flexibility. Some studies have also applied a Bowden cable to the upper limb [10,11] and hands [12,13]. However, friction is produced during transmission loop. However, to overcome the nonlinear and unknown disturbance, the sliding-mode control will produce high-frequency chattering. When the inner loop adopts the DOB, designing the controller based on the force error signal can introduce noise from the force sensor into the inner loop. Therefore, a feedforward-feedback proportion-differentiation controller (FFPD) based on a nominal model is proposed in this paper.
The main contribution of this paper is the design of an active loading control strategy, with a Bowden cable for transmission. Firstly, the control framework of the outer force loop based on the inner velocity loop is established. To improve the anti-interference ability of the inner loop, this paper proposes a DOB. The inner loop is equivalent to the designed nominal model, which effectively compensates for the parameter perturbation of the inner loop and the friction disturbance caused by the Bowden cable. The influence of modeling errors, with an increasing frequency, on the inner loop's dynamic performance is further discussed. Moreover, the stability condition of the inner loop is obtained. On this basis, the FFPD is proposed, which can avoid excessive force sensor noise being introduced into the inner loop and improve the force loading performance of the outer loop through pole configuration.
This paper is organized as follows. The mathematical modeling of the system is established in Section 2. The inner loop and outer-loop controllers are elaborated in Section 3. The effectiveness of the method proposed in this paper is verified in Section 4. The conclusion is given in Section 5.

Mathematical Modeling of the System
A Bowden cable can transmit displacement through the relative movement of an internal cable relative to an external sheath, which has the advantages of a flexible layout and long-distance transmission. Therefore, a wearable exoskeleton with a Bowden cable for transmission is introduced in this paper. A power supply, motors, controllers, and other components are integrated near the body's center of gravity. The power is transmitted to the lower limb fixation mechanism of the knee joint through the Bowden cable.
The fixation mechanism of the knee joint is shown in Figure 1. The cable sheath is connected to the motor base and the support module. Moreover, the inner cable is connected to the motor output shaft and the knee joint sheave. Two lightweight carbon fiberboards (LCFB) form a rotating joint, and the hinge mechanism forms a rotating joint on the coronal plane of the knee joint. The baffles are connected to the thigh and calve with a few flexible straps. Since the assistance needs to be performed based on the human body's active movement, an encoder at the knee joint is installed. Moreover, a force loading control framework based on a speed loop is proposed.
inner loop adopts the DOB, designing the controller based on the force error signal c introduce noise from the force sensor into the inner loop. Therefore, a feedforwa feedback proportion-differentiation controller (FFPD) based on a nominal model proposed in this paper.
The main contribution of this paper is the design of an active loading control strate with a Bowden cable for transmission. Firstly, the control framework of the outer fo loop based on the inner velocity loop is established. To improve the anti-interferen ability of the inner loop, this paper proposes a DOB. The inner loop is equivalent to designed nominal model, which effectively compensates for the parameter perturbati of the inner loop and the friction disturbance caused by the Bowden cable. The influen of modeling errors, with an increasing frequency, on the inner loop's dynam performance is further discussed. Moreover, the stability condition of the inner loop obtained. On this basis, the FFPD is proposed, which can avoid excessive force sen noise being introduced into the inner loop and improve the force loading performance the outer loop through pole configuration.
This paper is organized as follows. The mathematical modeling of the system established in Section 2. The inner loop and outer-loop controllers are elaborated Section 3. The effectiveness of the method proposed in this paper is verified in Section The conclusion is given in Section 5.

Mathematical Modeling of the System
A Bowden cable can transmit displacement through the relative movement of internal cable relative to an external sheath, which has the advantages of a flexible layo and long-distance transmission. Therefore, a wearable exoskeleton with a Bowden ca for transmission is introduced in this paper. A power supply, motors, controllers, a other components are integrated near the body's center of gravity. The power transmitted to the lower limb fixation mechanism of the knee joint through the Bowd cable.
The fixation mechanism of the knee joint is shown in Figure 1. The cable sheath connected to the motor base and the support module. Moreover, the inner cable connected to the motor output shaft and the knee joint sheave. Two lightweight carb fiberboards (LCFB) form a rotating joint, and the hinge mechanism forms a rotating jo on the coronal plane of the knee joint. The baffles are connected to the thigh and ca with a few flexible straps. Since the assistance needs to be performed based on the hum body's active movement, an encoder at the knee joint is installed. Moreover, a fo loading control framework based on a speed loop is proposed. However, friction will be produced during transmission through the Bowden cable. More than half of the motor torque is used to overcome the Bowden cable's friction [27]. In the process of human movement, the Bowden cable's angle changes with the swing of the thigh. The friction between the cable sheath and the inner cable will also change constantly. These factors all increase the difficulty of friction modeling. On the other hand, due to the small power of the integrated motor, the anti-interference ability is insufficient. Under the interference of the Bowden cable's friction, the parameter perturbation of the motor will be increased, which will further deteriorate its dynamic performance. Therefore, due to the limited number of sensors and the actuator's poor anti-interference ability, the controller design should be used to improve the dynamic performance of force loading.
This article pays more attention to solving the following problems: compensating for the Bowden cable's friction and improving the dynamic performance of force loading. Therefore, the exoskeleton and human movement are not involved in this paper. The end of the inner cable is connected to the fixed base. The entire system can be regarded as active force loading, with a Bowden cable for transmission. The system's mathematical model only considers the driver, integrated servo motor, Bowden cable, spring, and force sensor, as shown in Figure 2. F d is the force command. F is the force response.  However, friction will be produced during transmission through the Bowden More than half of the motor torque is used to overcome the Bowden cable's friction In the process of human movement, the Bowden cable's angle changes with the sw the thigh. The friction between the cable sheath and the inner cable will also c constantly. These factors all increase the difficulty of friction modeling. On the other due to the small power of the integrated motor, the anti-interference ability is insuff Under the interference of the Bowden cable's friction, the parameter perturbation motor will be increased, which will further deteriorate its dynamic perform Therefore, due to the limited number of sensors and the actuator's poor anti-interfe ability, the controller design should be used to improve the dynamic performance of loading.
This article pays more attention to solving the following problems: compensati the Bowden cable's friction and improving the dynamic performance of force loa Therefore, the exoskeleton and human movement are not involved in this paper. Th of the inner cable is connected to the fixed base. The entire system can be regard active force loading, with a Bowden cable for transmission. The system's mathem model only considers the driver, integrated servo motor, Bowden cable, spring, and sensor, as shown in Figure 2. d F is the force command. F is the force response.
where e J is the moment of inertia equivalent to the motor shaft, L is the induc coefficient, i K is the current feedback coefficient, R is the resistance coefficient, the amplification factor of the velocity loop, u k is the amplification factor, e C motor back EMF coefficient, e b is the equivalent viscous damping, m k is the t coefficient, and v k is the feedback coefficient of the velocity loop.
The moment balance equation is: where J e is the moment of inertia equivalent to the motor shaft, L is the inductance coefficient, K i is the current feedback coefficient, R is the resistance coefficient, k d is the amplification factor of the velocity loop, k u is the amplification factor, C e is the motor back EMF coefficient, b e is the equivalent viscous damping, k m is the torque coefficient, and k v is the feedback coefficient of the velocity loop. The moment balance equation is: where τ m is the output torque of the current loop, τ sm is the equivalent torque from the interaction force to the motor shaft, τ f is the equivalent torque from the Bowden cable friction to the motor shaft, s is a Laplacian operator, .
θ is the angular velocity of the motor shaft, n is the reduction ratio, F is the interaction force, and r m is the radius of the motor sheave.
Combined with Equations (1)-(4), the block diagram of the inner loop can be obtained, as shown in Figure 3. k is the spring stiffness, and f ( . θ) is a function of the Bowden cable's friction. It can be seen, from Figure 3, that the block diagram of the inner loop is composed of a servo motor current loop, speed loop, and mechanism expression. When the motor is set as the current loop, friction torque disturbance is not included in the current loop. Therefore, it can only be compensated for by the Bowden cable friction model. When the motor is set as the velocity loop, the disturbance caused by friction torque is included. On the one hand, the velocity loop can decrease the disturbance. On the other hand, the disturbance can be compensated for by the DOB. Moreover, the number of sensors can be reduced.
interaction force to the motor shaft,  f is the equivalent torque from the Bowden ca friction to the motor shaft, s is a Laplacian operator,  is the angular velocity of t motor shaft, n is the reduction ratio, F is the interaction force, and m r is the radius the motor sheave.
Combined with Equations (1)-(4), the block diagram of the inner loop can obtained, as shown in Figure 3. k is the spring stiffness, and   f  is a function of t Bowden cable's friction. It can be seen, from Figure 3, that the block diagram of the inn loop is composed of a servo motor current loop, speed loop, and mechanism expressio When the motor is set as the current loop, friction torque disturbance is not included the current loop. Therefore, it can only be compensated for by the Bowden cable fricti model. When the motor is set as the velocity loop, the disturbance caused by fricti torque is included. On the one hand, the velocity loop can decrease the disturbance. O the other hand, the disturbance can be compensated for by the DOB. Moreover, t number of sensors can be reduced. Figure 3. Block diagram of the inner loop.

Design of the Inner-Loop Controller
In a gait cycle, there is only one peak in the unidirectional assist curves of the hip a ankle joints. The force command can be set to 0 during the rest of a gait cycle. There a three peaks in the flexion direction of the knee joint. Moreover, the knee's maximu flexion angular velocity is twice that of the hip joint and ankle joint. Since the inner lo needs to overcome the interference of the friction and interaction force, the output pow is enormous. The internal parameters of the motor have a certain nonlinearity and tim variance. In this case, it is not ideal to enhance the dynamic response using t feedforward model. Generally, although the inner loop's nominal model is not the sam as the actual model, the transfer function is of the same type. Therefore, the inner loo nominal model can be expressed as follows: Gs as bs c .
The physical meaning of each nominal parameter in Equation (5) is the same as th of each nominal parameter in Equation (1). The DOB can be designed as follows:

Design of the Inner-Loop Controller
In a gait cycle, there is only one peak in the unidirectional assist curves of the hip and ankle joints. The force command can be set to 0 during the rest of a gait cycle. There are three peaks in the flexion direction of the knee joint. Moreover, the knee's maximum flexion angular velocity is twice that of the hip joint and ankle joint. Since the inner loop needs to overcome the interference of the friction and interaction force, the output power is enormous. The internal parameters of the motor have a certain nonlinearity and time variance. In this case, it is not ideal to enhance the dynamic response using the feedforward model. Generally, although the inner loop's nominal model is not the same as the actual model, the transfer function is of the same type. Therefore, the inner loop's nominal model can be expressed as follows:Ĝ 2 (s) =ĝ as 2 +bs +ĉ .
The physical meaning of each nominal parameter in Equation (5) is the same as that of each nominal parameter in Equation (1). The DOB can be designed as follows: whered is the equivalent disturbance expressed by the DOB, e is the velocity error, and d is the equivalent disturbance caused by τ sm and τ f , the formula of which is: Combining Equations (5)-(7), the control block diagram of the inner loop based on the DOB can be simplified, as shown in Figure 4. and d is the equivalent disturbance caused by  sm and  f , the formula of which is: Combining Equations (5)- (7), the control block diagram of the inner loop based on the DOB can be simplified, as shown in Figure 4. Gs , which will affect the stability of the system. Therefore, it is necessary to design a low-pass filter

 
Qs. According to Figure 4, the transfer functions are as follows: The error between the nominal and the actual models is small in the low-frequency range, as shown in Equation (11) In the assistance process, the inner loop needs to overcome the friction disturbance torque and interactive torque. Therefore, the actual model of the inner loop has parameter perturbation. Combining Figure 4 and Equation (8), Equation (12) can be obtained as follows: In Figure 4, ξ is the high-frequency noise of the speed signal. Q(s) is a low-pass filter. The requirement of keeping Q(s)Ĝ −1 2 (s) regularization should be satisfied. As shown in Figure 4, ξ can be amplified byĜ −1 2 (s), which will affect the stability of the system. Therefore, it is necessary to design a low-pass filter Q(s). According to Figure 4, the transfer functions are as follows: .
. θ ξ The error between the nominal and the actual models is small in the low-frequency range, as shown in Equation (11) G 2 (s) ≈Ĝ 2 (s). (11) In the assistance process, the inner loop needs to overcome the friction disturbance torque and interactive torque. Therefore, the actual model of the inner loop has parameter perturbation. Combining Figure 4 and Equation (8), Equation (12) can be obtained as follows: where ∆G 2 (s) is the perturbation of the actual model relative to the nominal model in the presence of disturbance. Combining Equations (8), (11) and (12), the sensitivity function can be obtained as follows: Then, the complementary sensitivity function can be written as: According to the robust stability theorem, the stability condition of the inner loop can be obtained as follows: The sinusoidal sweep signal, with an amplitude of 45 r/min and frequency of 1-7 Hz, is adopted as the speed command of the inner loop. The tracking curve in the 3-7 Hz frequency band of the sweep signal is intercepted to observe the tracking effect more clearly. The dynamic response curve of the inner loop without a load is shown in Figure 5. It can be seen, from Figure 5, that the amplitude of the inner loop has almost no attenuation. Only the phase lags by 15 ms.
According to the robust stability theorem, the stability condition of the inner loop can be obtained as follows: The sinusoidal sweep signal, with an amplitude of 45 r/min and frequency of 1-7 Hz, is adopted as the speed command of the inner loop. The tracking curve in the 3-7 Hz frequency band of the sweep signal is intercepted to observe the tracking effect more clearly. The dynamic response curve of the inner loop without a load is shown in Figure  5. It can be seen, from Figure 5, that the amplitude of the inner loop has almost no attenuation. Only the phase lags by 15 ms.
According to Equation (16), the design of the filter should include the ability to suppress disturbance and noise sensitivity. In the low-frequency band, the motor's response is consistent with the nominal model established in this paper. The DOB can effectively compensate for the friction disturbance torque and interaction torque. Moreover, the highfrequency noise has little effect on the output of the motor. In the high-frequency band, Q(s) = 0. Substituting it into Equations (8)-(10), the results are as follows: In the high-frequency band, the motor performance shows its actual dynamic response. The influence of friction torque and interaction torque cannot be suppressed. High-frequency noise has little effect on the output speed. Since human motion is in the low-frequency band, the dynamic response of the velocity loop in the high-frequency range can be ignored. Therefore, within the bandwidth of Q(s), the torque disturbance, the modeling perturbance between the actual system and the nominal model will be compensated for. The dynamic performance of the inner loop can be regarded asĜ 2 (s). Combined with Equations (8), (16) and (17), Equation (18) can be rewritten as: According to Equations (15) and (18), the inner loop's anti-interference ability is enhanced in the low-frequency band. The conditions for improving the stability of the inner loop are obtained considering the modeling error. The accuracy of the nominal model should be improved as much as possible, on the one hand, to ensure the inner loop's stability, and on the other, to reduce the operating range of the inner loop in order to minimize the perturbation. Under the premise of the inner loop stability, the broader the bandwidth of the filter is, the stronger the anti-interference ability in the low-frequency band is. Consequently, the potential force loading performance can be further enhanced.

Design of the Oute-Loop Controller
The block diagram of the outer loop is shown in Figure 6. F d is the force command, and F is the force response.
According to Equations (15) and (18), the inner loop's anti-interference ability is enhanced in the low-frequency band. The conditions for improving the stability of the inner loop are obtained considering the modeling error. The accuracy of the nominal model should be improved as much as possible, on the one hand, to ensure the inner loop's stability, and on the other, to reduce the operating range of the inner loop in order to minimize the perturbation. Under the premise of the inner loop stability, the broader the bandwidth of the filter is, the stronger the anti-interference ability in the lowfrequency band is. Consequently, the potential force loading performance can be further enhanced.

Design of the Oute-Loop Controller
The block diagram of the outer loop is shown in Figure 6. d F is the force command, and F is the force response.   According to Equation (18), the dynamic performance of the inner loop in the lowfrequency band is consistent with the nominal model. In other words, the DOB control method proposed in the inner loop can effectively compensate for not only the interference of friction with the inner loop, but also the nonlinearity and modeling errors caused by the parameter perturbation. Therefore, the outer-loop controller can also be designed based on G 2 (s), which can effectively avoid the model parameter perturbation and nonlinearity of the inner loop.
In the low-frequency band, combined with Equations (5), (18) and Figure 6, the outer loop model is shown in the following: where b =ˆbâ , c =ĉâ , g =ĝ r m an .
The force error of the outer loop can be expressed as: By substituting Equation (21) into Equation (19) and combining it with the PD control, the control rate can be obtained as: The equation of the force error can be obtained by substituting Equation (22) into Equation (22): Using the Laplace transform, Equation (23) can be rewritten as: Equation (24) is a characteristic equation concerning force error. The solution of Equation (24) can represent the rapidity and stability of the outer loop. The dynamic performance in the outer loop can be changed by adjusting k p and k d . To ensure that the eigenvalues are located on the left side of the virtual axis, the desired characteristic equation can be written as: The coefficients of Equations (24) and (25) are made to be equal, that is: where x 1 , x 2 , and x 3 are the expected poles. Hence, k p and k d can be obtained according to the expected x 1 , x 2 , and x 3 .

Test Setup
As a preliminary detection, the corresponding force loading control strategy, with a Bowden cable for transmission, is proposed. The scheme process is shown in Figure 7, including the control system, mechanical structure, and sensors.

Test Setup
As a preliminary detection, the corresponding force loading control strategy, with a Bowden cable for transmission, is proposed. The scheme process is shown in Figure 7, including the control system, mechanical structure, and sensors.
The control system is a single task mode, which includes a host computer and a target computer. The PC can be used as a host computer, and the LabVIEW software can be run on it. The control program can be built and compiled in the LabVIEW software. Then, the compiled program can be deployed to the target computer through a shared local area network (LAN). An industrial computer can be used as the target computer. The target computer has a high-speed real-time computing environment, which can collect data and send control signals through the program deployed therein. To improve the integration, the Sbrio-9636 is selected as the target computer. The vendor of the servo motor is the Hai Tai Electromechanical Equipment in YiWu, China. The model used in this paper is the HT-02. The servo motor integrates a 6:1 reduction ratio, and the total weight is 480 g. The rated torque and rated speed after the reducer are 7.6 Nm and 250 RPM, respectively. The mechanical structure is composed of the servo motor, Bowden cable, spring, and three Lshaped supports. The motor and three L-shaped supports are fixed on a steel plate. The cable sheath is connected to a 1# L-shaped support and 2# L-shaped support. The inner cable is connected to the motor output shaft and the spring. The inner cable and the cable sheath move relative to each other and are driven by the servo motor.  The torque curve in the knee joint's flexion direction has three peaks, and its frequency gradually increases. The curve selected in this paper can basically cover the force loading frequency range of human motion. Therefore, the average force curve of the knee joint flexion direction is obtained by the Vicon system when the human body is The control system is a single task mode, which includes a host computer and a target computer. The PC can be used as a host computer, and the LabVIEW software can be run on it. The control program can be built and compiled in the LabVIEW software. Then, the compiled program can be deployed to the target computer through a shared local area network (LAN). An industrial computer can be used as the target computer. The target computer has a high-speed real-time computing environment, which can collect data and send control signals through the program deployed therein. To improve the integration, the Sbrio-9636 is selected as the target computer. The vendor of the servo motor is the Hai Tai Electromechanical Equipment in YiWu, China. The model used in this paper is the HT-02. The servo motor integrates a 6:1 reduction ratio, and the total weight is 480 g. The rated torque and rated speed after the reducer are 7.6 Nm and 250 RPM, respectively. The mechanical structure is composed of the servo motor, Bowden cable, spring, and three L-shaped supports. The motor and three L-shaped supports are fixed on a steel plate. The cable sheath is connected to a 1# L-shaped support and 2# L-shaped support. The inner cable is connected to the motor output shaft and the spring. The inner cable and the cable sheath move relative to each other and are driven by the servo motor.
The torque curve in the knee joint's flexion direction has three peaks, and its frequency gradually increases. The curve selected in this paper can basically cover the force loading frequency range of human motion. Therefore, the average force curve of the knee joint flexion direction is obtained by the Vicon system when the human body is walking at 4.5 km/h. To verify the effectiveness of the method proposed in this paper, the dynamic performance of the system using the different control methods is comprehensively compared through experiments. Firstly, both the inner loop and outer loop adopt PID control to verify the necessity of the control method proposed in this article. Secondly, the experimental results in Sections 4.2 and 4.3 were compared to further verify the influence of the inner loop performance on the outer loop. Finally, springs with different stiffnesses are used to verify the influence of model errors on the dynamic performance of the system.
Because of the high flexibility, the inner cable may not be tensioned before the experiment. Therefore, at the beginning of the experiment, the motor rotates at the speed of 0.2 r/min until the force sensor value is greater than 5 N.

Friction Verification of the Bowden Cable
In this section, the tension relationship of the inner cable at the 1# and 3# L-shaped supports is used to verify the hysteresis phenomenon of the Bowden cable's friction. The proportional control is adopted, and the force response curve is shown in Figure 8. F i and F are the force of the inner cable at the 1# and 3# L-shaped supports, respectively. It can be seen, from Figure 8, that 43% of F i is used to overcome the Bowden cable's friction. It also can be seen, from Figure 9, that the friction of the Bowden cable has a strong hysteresis. Therefore, the Bowden cable's friction cannot be ignored or simply suppressed by closed-loop control during the force loading. Because of the high flexibility, the inner cable may not be tensioned be experiment. Therefore, at the beginning of the experiment, the motor rotates at t of 0.2 r/min until the force sensor value is greater than 5 N.

Friction Verification of the Bowden Cable
In this section, the tension relationship of the inner cable at the 1# and 3# L supports is used to verify the hysteresis phenomenon of the Bowden cable's fric proportional control is adopted, and the force response curve is shown in Fig   and F are the force of the inner cable at the 1# and 3# L-shaped supports, respec can be seen, from Figure 8, that 43% of i F is used to overcome the Bowden cable' It also can be seen, from Figure 9, that the friction of the Bowden cable has hysteresis. Therefore, the Bowden cable's friction cannot be ignored or simply su by closed-loop control during the force loading.

PID Control in Both the Inner and Outer Loop
Firstly, PID control is adopted in both the inner loop and outer loop, the responses of which are shown in Figures 10-12. As shown in Figure 10, the phase lag from 0 s to 0.4 s is about 0.1 s, and there is almost no amplitude attenuation. However, within 0.8-1.1 s, the frequency of the force curve increases. The force response is almost constant during the rapid change in the force command, indicating that the dynamic performance of the outer loop is poor. Figure 12 shows the integral curve of the velocity response. Moreover, the force response of Figure 10 is scaled to the same order of magnitude as the integral curve of the velocity response. As shown in Figure 12, F m is the force response curve after scaling. It can be seen that the displacement response of the Bowden cable close to the motor is significantly earlier than the force response, which shows that the friction hinders the Bowden cable's transmission. It can also be seen, from Figure 12, that when the loading force is small, the flexible connection between the Bowden cable and the spring has a more significant impact on the force loading. Therefore, it is difficult for PID control to meet the force loading performance. ators 2021, 10, x FOR PEER REVIEW Firstly, PID control is adopted in both the inner loop and ou which are shown in Figures 10-12. As shown in Figure 10, the p is about 0.1 s, and there is almost no amplitude attenuation. H the frequency of the force curve increases. The force response i the rapid change in the force command, indicating that the dyn outer loop is poor. Figure 12 shows the integral curve of the velo the force response of Figure 10 is scaled to the same order of m curve of the velocity response. As shown in Figure 12, m F is the f scaling. It can be seen that the displacement response of the B motor is significantly earlier than the force response, which show the Bowden cable's transmission. It can also be seen, from Figure force is small, the flexible connection between the Bowden cable significant impact on the force loading. Therefore, it is difficult fo force loading performance.

The FFPD in the Outer Loop
The sinusoidal sweep signal, with an amplitude of 45 r/min and frequency of 1-7 Hz, is adopted as the speed command in the inner loop. The tracking curve in the 3-7 Hz range is intercepted to observe the tracking effect more clearly. The identified model without a load is as follows: The dynamic response of the nominal model (NM) and actual model without a load (WAM) is shown in Figure 13. The nominal model shown in Equation (27) reduces the time delay by 5 ms, compared to that of the actual model. The dynamic performance of the nominal model is basically consistent with that of the actual model. Therefore, the feedforward compensation and the DOB can be designed based on Equation (27).

The FFPD in the Outer Loop
The sinusoidal sweep signal, with an amplitude of 45 r/min and frequency of 1-7 Hz, is adopted as the speed command in the inner loop. The tracking curve in the 3-7 Hz range is intercepted to observe the tracking effect more clearly. The identified model without a load is as follows:Ĝ 2 (s) = 3915 s 2 + 198.1s + 23, 520 .
The dynamic response of the nominal model (NM) and actual model without a load (WAM) is shown in Figure 13. The nominal model shown in Equation (27) reduces the time delay by 5 ms, compared to that of the actual model. The dynamic performance of the nominal model is basically consistent with that of the actual model. Therefore, the feedforward compensation and the DOB can be designed based on Equation (27).
The FFPD is adopted to further verify the influence of the inner loop performance on the outer loop. According to Equation (27), the dynamic responses under the FFPD are obtained, as shown in Figures 14-16. In Figure 14, the phase lag from 0-0.4 s is about 0.06 s, which reduces the time delay by 40%, compared with that of PD control. Therefore, it can be proved that the phase lag can be compensated for by the feedforward control. From 0.8-1.1 s, the tracking performance of the first peak is greatly improved, compared with PD control. Moreover, the phase of the first peak is advanced by 25%. However, the tracking performance of the second peak is poor. The dynamic response of the inner loop becomes worse in Figure 15, which further limits the improvement of the dynamic performance in the outer loop. Therefore, it is necessary to compensate for the Bowden cable's friction and improve the inner loop's anti-interference ability. As shown in both Figures 12 and 16, the interaction force does not change proportionally with the displacement change. This also proves the phenomenon of friction hysteresis.
is intercepted to observe the tracking effect more clearly. The identified model witho load is as follows: The dynamic response of the nominal model (NM) and actual model without a (WAM) is shown in Figure 13. The nominal model shown in Equation (27) reduce time delay by 5 ms, compared to that of the actual model. The dynamic performan the nominal model is basically consistent with that of the actual model. Therefore feedforward compensation and the DOB can be designed based on Equation (27). The FFPD is adopted to further verify the influence of the inner loop performanc the outer loop. According to Equation (27), the dynamic responses under the FFPD obtained, as shown in Figures 14-16. In Figure 14, the phase lag from 0-0.4 s is about s, which reduces the time delay by 40%, compared with that of PD control. Therefo can be proved that the phase lag can be compensated for by the feedforward control. F 0.8-1.1 s, the tracking performance of the first peak is greatly improved, compared PD control. Moreover, the phase of the first peak is advanced by 25%. However tracking performance of the second peak is poor. The dynamic response of the inner becomes worse in Figure 15, which further limits the improvement of the dynamic performance in the outer loop. Therefore, it is necessary to compensate for the Bowden cable's friction and improve the inner loop's anti-interference ability. As shown in both Figure 12 and Figure 16, the interaction force does not change proportionally with the displacement change. This also proves the phenomenon of friction hysteresis.   becomes worse in Figure 15, which further limits the improvement of the dynamic performance in the outer loop. Therefore, it is necessary to compensate for the Bowden cable's friction and improve the inner loop's anti-interference ability. As shown in both Figure 12 and Figure 16, the interaction force does not change proportionally with the displacement change. This also proves the phenomenon of friction hysteresis.     Feedforward control can improve the dynamic performance of force loading to a certain extent. However, due to the limitation of the inner loop, the force loading performance still cannot be further improved. Feedforward control can improve the dynamic performance of force loading to a certain extent. However, due to the limitation of the inner loop, the force loading performance still cannot be further improved.

The FFPD in the Outer Loop and the DOB in the Inner Loop
The following conclusions are obtained through theoretical analysis. The enhancement of the dynamic performance in the inner loop can further improve the dynamic response in the outer loop. Moreover, the stability of the inner loop depends on the perturbation of the actual model, according to Equation (15). With the increase of the frequency, the inner loop's perturbation is stronger, and the anti-interference ability becomes worse. To verify the correctness of the above analysis, the spring stiffness k 2 is increased and compared with the dynamic performance of k 1 . The values of the spring stiffnesses k 1 and k 2 are 10, 000 N/m and 30, 000 N/m, respectively. The dynamic response is shown in Figures 17-19, when the outer-loop FFPD controller, based on an inner-loop DOB, is adopted. ators 2021, 10, x FOR PEER REVIEW

The FFPD in the Outer Loop and the DOB in the Inner Loop
The following conclusions are obtained through th enhancement of the dynamic performance in the inner loop dynamic response in the outer loop. Moreover, the stability of th the perturbation of the actual model, according to Equation (15) frequency, the inner loop's perturbation is stronger, and the becomes worse. To verify the correctness of the above analysis, increased and compared with the dynamic performance of 1 k . stiffnesses 1 k and 2 k are 10,000 N m and 30,000 N m , res response is shown in Figures 17-19, when the outer-loop FFPD inner-loop DOB, is adopted.     According to Figure 17, the phase lag of the force curve is about 0.07 s in 0-0.4 s, which reduces the time delay by 30%, compared with that of PD control. The response amplitude of the second peak increases from 7 N to 80 N and 66 N, respectively. Compared with the feedforward controller, the resulting amplitude is much improved. Comparing Figure 15 with Figure 18, the command amplitude and frequency of the inner loop have both increased. This shows that the DOB controller can not only effectively compensate for the Bowden cable's friction, but also effectively enhance the antiinterference ability in the inner loop. Moreover, the dynamic performance is approximately the nominal model. Designing the feedforward control in the outer loop based on the nominal model is also more conducive to improving the force loading performance.
The force response curve under 1 k has an overshoot of 8 N at 1.2 s. Correspondingly, the force response curve under 2 k has no overshoot. Besides, the force response curve with 2 k is ahead of that with 1 k in phase. It can also be verified, from Figure 18 and Figure 19, that when the stiffness is 1 k , the bandwidth required by the inner loop increases. Therefore, it is verified that the larger the model error, the worse the system stability.
The experimental results show that the higher the spring stiffness, the better the force loading response. Comparing Figure 14 and Figure 17, it can be concluded that only by improving the inner loop's dynamic response can the outer loop's dynamic performance be further improved. While there is still a phase lag of 0.06 s, the number of sensors required is reduced. Based on the above analysis, the inner loop adopts the DOB, and the outer loop adopts the FFPD, which can meet the force loading performance of the knee joint. According to Figure 17, the phase lag of the force curve is about 0.07 s in 0-0.4 s, which reduces the time delay by 30%, compared with that of PD control. The response amplitude of the second peak increases from 7 N to 80 N and 66 N, respectively. Compared with the feedforward controller, the resulting amplitude is much improved. Comparing Figure 15 with Figure 18, the command amplitude and frequency of the inner loop have both increased. This shows that the DOB controller can not only effectively compensate for the Bowden cable's friction, but also effectively enhance the anti-interference ability in the inner loop. Moreover, the dynamic performance is approximately the nominal model. Designing the feedforward control in the outer loop based on the nominal model is also more conducive to improving the force loading performance.

Conclusions
The force response curve under k 1 has an overshoot of 8 N at 1.2 s. Correspondingly, the force response curve under k 2 has no overshoot. Besides, the force response curve with k 2 is ahead of that with k 1 in phase. It can also be verified, from Figures 18 and 19, that when the stiffness is k 1 , the bandwidth required by the inner loop increases. Therefore, it is verified that the larger the model error, the worse the system stability.
The experimental results show that the higher the spring stiffness, the better the force loading response. Comparing Figures 14 and 17, it can be concluded that only by improving the inner loop's dynamic response can the outer loop's dynamic performance be further improved. While there is still a phase lag of 0.06 s, the number of sensors required is reduced. Based on the above analysis, the inner loop adopts the DOB, and the outer loop adopts the FFPD, which can meet the force loading performance of the knee joint.

Conclusions
In this paper, an outer-loop FFPD controller with a DOB in the inner loop is proposed. Besides, the force loading performance with the Bowden cable for transmission is effectively