Next Article in Journal
The Influence of Speed Ratio on the Nonlinear Dynamics of a Magnetic Suspended Dual-Rotor System with a Fixed-Point Rubbing
Next Article in Special Issue
Extended State Observer-Based Sliding Mode Control Design of Two-DOF Lower Limb Exoskeleton
Previous Article in Journal
Research on Speed Control Methods and Energy-Saving for High-Voltage Transmission Line Inspection Robots along Cable Downhill
Previous Article in Special Issue
Research on Machine Vision-Based Control System for Cold Storage Warehouse Robots
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling and RBF Control of Low-Limb Swinging Dynamics of a Human–Exoskeleton System

School of Aeronautics and Astronautics, University of Electronic Science and Technology of China, Chengdu 611731, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Actuators 2023, 12(9), 353; https://doi.org/10.3390/act12090353
Submission received: 10 August 2023 / Revised: 29 August 2023 / Accepted: 3 September 2023 / Published: 6 September 2023

Abstract

:
With the increase in the elderly population in China and the growing number of individuals who are unable to walk normally, research on lower limb exoskeletons is becoming increasingly important. This study proposes a complete dynamic model parameter identification scheme for the human–machine coupling model of lower limb exoskeletons. Firstly, based on the coupling model, the excitation trajectory is optimized, data collection experiments are conducted, and the dynamic parameter vector of the system is identified using the least squares method. Secondly, this lays the foundation for designing adaptive control based on RBF neural network approximation. Thirdly, the Lyapunov function is used to prove that the RBF neural network adaptive controller can achieve stable tracking of the lower limb exoskeleton. Finally, simulation analysis reveals that increasing the gains of the RBF controllers effectively reduces tracking errors. Furthermore, the tracking errors and control torques show that adaptive control based on the RBF neural network approximation works well.

1. Introduction

According to some references, by 2030, 18.2% of Chinese people will be over 65 years old [1]. Additionally, the number of people with physical disabilities will reach 24.12 million, accounting for 29.07% of the total number of disabled individuals, among which there are approximately 1.58 million people with lower limb paralysis [2]. Currently, China is becoming an aging society, and with the growth of the elderly population and the occurrence of various accidents, the number of people who have difficulty in walking is increasing year by year. However, due to various reasons, only a small portion of the population can receive timely and effective rehabilitation treatment [3]. This can cause significant potential harm to the lives of the patients and their families [4]. At this point, the emergence of lower limb exoskeletons brings significant help to patients and their families. Lower limb exoskeletons can assist patients in effective rehabilitation training, improve their lower limb motor function, and alleviate the burden on the patients’ caregivers. The emergence of lower limb exoskeletons brings new hope to patients and improves their quality of life [5]. Generally speaking, lower limb exoskeletons have two different objectives: assisting patients in rehabilitation training, and assisting in human working activities [6]. Robotic technologies that assist physical health and provide support for the elderly are rapidly developing. In the past three decades, lower limb exoskeleton robotics and assistive technologies have been a focus of attention in the industry [7]. As a typical application of human–machine interaction devices, the perception and control of lower limb exoskeleton robots will significantly affect the actual wearing effects of users in lower limb assistance or augmentation [8]. Therefore, in order to achieve the swinging control of the lower limb exoskeleton, it is essential to design an appropriate swing controller during the development stage [9].
Currently, an increasing number of researchers aim to enhance the development and improvement of lower limb exoskeletons through control, thereby improving their performance [10]. With the rise of model-based control techniques, knowledge of the precise dynamic mathematical model of lower limb exoskeletons is required. However, in practical systems, various uncertainties exist that have a significant impact on the model-based controllers, thereby severely affecting their performance. Therefore, designing a reasonable control algorithm to address the uncertainties in lower limb exoskeleton systems has become a research hotspot [11]. Currently, optimization control strategies for lower limb exoskeletons need to meet the criteria of safety, stability, effective rehabilitation, assistance, and efficient control simultaneously [12]. Furthermore, for lower limb exoskeletons, controllers need to demonstrate good tracking performance and stability while minimizing tracking errors [13]. Therefore, appropriate control algorithms are required to ensure the stability and robustness of lower limb exoskeletons during motion [14]. Currently, the main control strategies for lower limb exoskeletons include position tracking control, force impedance control, biological signal control, and more. Among them, position tracking control serves as the foundation for other control methods [15]. Additionally, the control modes of exoskeletons primarily involve active control and passive control. In active control, the exoskeleton provides necessary assistance and closely follows human motion. In passive control, the human body acts as the load for the exoskeleton, and the exoskeleton is entirely driven by external forces. The key aspect of passive control lies in studying the dynamic model of the exoskeleton [16]. From this perspective, the development of control systems is one of the most critical parts of exoskeleton systems. RBF neural networks are chosen due to their excellent generalization ability, simple structure, and the ability to approximate any nonlinear function with arbitrary precision. They can be used to estimate the unknown parts of the hyper-local model [17]. Reference [18] considers the exoskeleton as a nonlinear motion system and applies RBF neural networks for real-time identification of the system. The current hip and ankle joint angles of the exoskeleton, along with the previous knee joint angle, are defined as the network inputs at the current time. A general error function is defined, which includes tracking errors of the exoskeleton and interaction torques between the pilot and the exoskeleton to ensure convergence. Subsequently, a mathematical model of the human–machine system is established. It can be seen that the RBF neural network controller performs well in controlling the lower limb exoskeleton robot through simulation experiments.
Furthermore, through extensive literature review, it is evident that determining the inertial parameters of the robot is necessary to develop advanced control algorithms. The design of nonlinear robot controllers typically relies on the robot model, and their performance directly depends on accurate inertial parameters. Due to the complexity of certain robot structures and the nonlinear nature of loads, determining dynamic parameters can be challenging, and parameter identification is the only effective method to obtain precise inertial parameters [19]. As model-based control is crucial, the identification of inertial parameters has gained extensive attention from researchers [20]. Identifying the inertial parameters first requires establishing a mathematical model and linearizing it to obtain a coefficient matrix for parameter identification [21]. In the offline dynamic identification stage the coefficient matrix is utilized, and the parameters are calculated using the least squares method [22].
It can be seen that the control of lower limb exoskeletons is crucial in their research. Model-based control methods have been used to achieve a desirable control performance. This study focuses on lower limb exoskeleton robots. Considering the individual differences among people, this work firstly determines the inertial parameters of the human–machine coupling model through experiments to establish a foundation for control simulations, rather than using average anatomical values as the model’s inertial parameters. In the control simulation stage, assuming the exoskeleton operates in passive control mode, an adaptive controller based on RBF neural network approximation is designed to track the desired trajectory. Simulations are conducted using MATLAB to indicate the control effects of the RBF neural network controller.

2. Materials and Methods

2.1. Coupled Dynamics Model of 2-DOF Lower Extremity Exoskeleton and Human Body

In this paper, a Lagrange modelling approach is utilized to research the dynamics of a lower limb exoskeleton robot and the human body. The Lagrange method firstly finds the total kinetic energy E k and the potential energy E p of the model, and then puts E k and E p into the Lagrange function to get L a = E k E p . Then, the Lagrange function equation can be derived from its partial derivation to find out the magnitude of the actuating torque needed to rotate the corresponding joints in the model.
Figure 1 shows the schematic diagram of a 2-DOF lower extremity exoskeleton coupled with the human body. The corresponding symbols for the physical parameters in the figure are shown in Table 1. We treat the human body and lower limb exoskeleton as a whole in this section.
According to the kinetic energy formula, the translational kinetic energy and the rotational kinetic energy of the center of mass, respectively, make up the kinetic energy of the thigh E k 1 and the shank E k 2 :
E k 1 = 1 2 m 1 c 1 2 θ ˙ 1 2 + 1 2 I 1 θ ˙ 1 2
E k 2 = 1 2 m 2 ( l 1 2 θ ˙ 1 2 + c 2 2 θ ˙ 1 2 + c 2 2 θ ˙ 2 2 + 2 L 1 c 2 θ ˙ 1 2 cos θ 2 + 2 L 1 c 2 θ ˙ 1 θ ˙ 2 cos θ 2 + 2 c 2 2 θ ˙ 1 θ ˙ 2 ) + 1 2 I 2 ( θ ˙ 1 + θ ˙ 2 ) 2
From the potential energy equation, the potential energy of the thigh E p 1 and the shank E p 2 can be obtained:
E p 1 = c 1 m 1 g cos θ 1
E p 2 = m 2 g [ L 1 cos θ 1 + c 2 cos ( θ 1 + θ 2 ) ]
where gravitational acceleration g = 9.8   m / s 2 .
Then the Lagrange function L a can be written as:
L a ( θ 1 , θ 2 , θ ˙ 1 , θ ˙ 2 ) = E k 1 + E k 2 E p 1 E p 2
Lagrange dynamic method is used to solve the system’s equations, and the Lagrange equation is:
T = t L a θ ˙ L a θ
where θ = [ θ 1 , θ 2 ] T 2 represents the dual angle vector representing the coupled model, T represents the joint generalized vector.
According to the above derivation, the dynamic model of the coupled human lower extremity exoskeleton can be represented as:
M ( θ ) θ ¨ + C ( θ , θ ˙ ) θ ˙ + G ( θ ) + τ f i c = τ
In Equation (7), M ( θ ) , C ( θ , θ ˙ ) and G ( θ ) represent the Inertia matrix, Coriolis matrix and Gravity matrix, τ f i c 2 represent the dual joint friction torque, τ 2 represent the dual joint control torque.
Convert the dynamic model into a linear form:
M ( θ ) θ ¨ + C ( θ , θ ˙ ) θ ˙ + G ( θ ) + τ f i c = A ( θ , θ ˙ , θ ¨ ) X = τ
In Equation (8), X 8 represents the parameter vector to be identified, A ( θ , θ ˙ , θ ¨ ) 2 × 8 represents the regression matrix composed of the angle, angular velocity, and angular acceleration of the hip and knee joints. The specific forms of the parameter vector to be identified and the regression matrix are provided in Appendix A.
Assuming that there are N sets of sampling points in the experiment, the sampling regression matrix A ¯ 2 N × 8 and the sampling torque vector τ ¯ 2 N × 1 can be represented as:
A ¯ = [ A ( 1 ) A ( 2 ) A ( N ) ] , τ ¯ = [ τ ( 1 ) τ ( 2 ) τ ( N ) ]
The parameter vector to be identified can be obtained using the least square method:
X ^ = ( A ¯ T A ¯ ) 1 A ¯ T τ ¯
Due to the presence of signification noise in the sampled data, designing a suitable excitation trajectory can help suppress the noise and thereby improve the accuracy of the parameter vector to be identified.
The excitation trajectory can be represented using a second-order Fourier series:
θ d = θ 0 + k = 0 n ¯ ( a i , k sin ( k ω t ) + b i , k cos ( k ω t ) ) θ ˙ d = k = 0 n ¯ ( a i , k k ω cos ( k ω t ) b i , k k ω sin ( k ω t ) ) θ ¨ d = k = 0 n ¯ ( a i , k ( k ω ) 2 sin ( k ω t ) b i , k ( k ω ) 2 cos ( k ω t ) )
In Equation (11), i = 1 , 2 , k = 1 , , n ¯ , t [ 0 , T ] and T represents the duration of the data collection experiment for design, ω represents the set base frequency, k represents the set sampling parameter, θ i , 0 represents the initial bias angle of the exoskeleton hip and knee joints, and a i , k and b i , k represent the parameters to be optimized. Considering the joint angle limitation range during human movement [23] as well as our hardware device restriction, the constraint intervals for the excitation trajectory are shown in Table 2.
In MATLAB, the Equation (11) is sampled to obtain a sampled regression matrix. By using the particle swarm algorithm to optimize the obtained matrix, the parameters a i , k and b i , k to be optimized can be obtained. The condition number of the sampling regression matrix is selected as the fitness function ( C o n d ( A ¯ ) ) for the optimization process. Meanwhile, the set excitation trajectory should follow the constraint conditions shown in Table 2. Based on the final set, the excitation trajectory is as follows:
F i t ( A ¯ ) = { C o n d ( A ¯ ) if { 0 . 0872 < θ 1 < 1 . 2217 1 . 6232 < θ 2 < 0 . 2269 1 . 6449 < θ ˙ , θ ¨ < 1.6449 I n f Otherwise

2.2. Adaptive Control of Exoskeleton Based on RBF Neural Network Approximation

The most classic control algorithms for exoskeletons include PID control, neural network control, fuzzy control, iterative learning control, and robot inversion control [24]. The PD controller is the most widely used control algorithm, especially in a nonlinear control system. The PD controller’s leading characteristics can be utilized to enhance the dynamic performance and robustness of the control system [25]. RBF neural networks, due to their good performance and simple structure, can avoid unnecessary and complex calculations that are commonly used in a nonlinear system [26]. Compared to the PD controller, the RBF neural network can eliminate the error caused by disturbance and accelerate convergence speed [27].
In practical application engineering, there may be unknown external disturbance acting on the torque τ d applied to the exoskeleton [28]. When an unknown external disturbance is added to system, the mathematical expression for the model of the limb exoskeleton is:
M θ ¨ + C θ ˙ + G = τ τ d τ f r i c
The M, C, G matrices represent the inertia matrix, Coriolis acceleration, and gravity matric of the identified exoskeleton and human body system.
The input of the system θ d ( t ) is the desired angle vector for the exoskeleton’s hip and knee joints, and the output θ ( t ) is the actual tracked angle vector for the exoskeleton’s hip and knee joints. Therefore, the error of the system is given by:
e ( t ) = θ d ( t ) θ ( t )
In order to achieve a better control performance, it is necessary to compensate for the unknown disturbances. Because the RBF neural network has a good error compensate capability, it is used to compensate for the unknown external disturbances in the exoskeleton. The error function is defined as:
q = e ˙ + Λ e
In Equation (15) Λ = Λ T > 0 , which represents the magnification. By combining Equations (14) and (15), we obtain:
θ ˙ = q + θ ˙ d + Λ e
Taking the derivative of Equation (16) with regards to time and left-multiplying by the matrix M , we obtain:
M q ˙ = M ( θ ¨ d θ ¨ + Λ e ˙ ) = M ( θ ¨ d + Λ e ˙ ) M θ ¨
Substituting Equations (13) to (17) into the previous equation, we obtain:
M q ˙ = C q τ + f + τ d
The expression for f in Equation (18) is:
f = M ( θ ¨ d + Λ e ) + C ( θ d + Λ e ) + G + τ f r i c
The unknown f represents the model uncertainty in practical engineering applications, and it needs to be approximated. Equation (19) indicates that this uncertainty is a nonlinear function, which can be approximated using an RBF neural network. The chosen RBF neural network is a three-layer feedforward network with a single hidden layer, which is more straightforward and computationally straightforward than higher level neural networks, while still satisfying our model’s control requirements. The transformation from the input layer to the hidden layer is nonlinear. The input signals of the network are denoted as x = [ e T e ˙ T q d T q ˙ d T q ¨ d T ] . The transformation function in the hidden layer is the radial basis function, represented as hj, j = 1, 2, …, m. The radial basis function is used in the Gaussian kernel function, and its expression is:
h j ( x ) = exp [ x c j 2 b j 2 ]
In Equation (20), c j represents the center of the j-th radial basis function, and b j represents the width of the j-th radial basis function.
The transformation from the hidden layer to the output layer in the RBF neural network is a linear transformation. The output expression is as follows:
f ^ = W ^ T H ( x )
Let
W ˜ = W W ^
In Equation (22), W F W m a x and W ^ T represent the weight matrix of the RBF neural network, while H ( x ) represents the output matrix of the radial basis functions.
Since the role of the neural network is to compensate for the model uncertainty, the control law of the designed RBF neural network controller is as follows:
τ = W ^ T H ( x ) = K v q v
In Equation (23), K v represents the amplification coefficient of the error, while v represents the robust term used to overcome approximation errors in the neural network.

2.3. Adaptive Analysis of Neural Network Stablility and Converagence

By substituting the control law (23) into Equation (18) and simplifying, we obtain:
M q ˙ = ( K v + C ) q + ξ 1
where ξ 1 = W ˜ T H ( x ) + ( ε + τ d ) + v .
Designing the robust term v as follows:
v = ( ε N + τ d ) sgn ( q )
In Equation (25), ε N and b d represent the upper bounds of ε and τ d .
The expression for the Lyapunov function is defined as:
L = 1 2 q T M q + 1 2 t r ( W ˜ T F 1 W ˜ )
In Equation (26), F is a positive definite matrix, and taking the derivative of Equation (26), we obtain:
L ˙ = q T M q ˙ + 1 2 q T M ˙ q + t r ( W ˜ T F 1 W ˜ ˙ )
By substituting Equation (22) into the previous equation, we obtain:
L ˙ = q T K v q + 1 2 q T ( M ˙ 2 C ) q + t r W ˜ T ( F 1 W ˜ ˙ + H q T ) + q T ( ε + τ d + v )
Taking into account the characteristics of the exoskeleton, we choose the adaptive law of the neural network as:
W ˜ ˙ = F H q T
Therefore:
L ˙ = q T K v q + q T ( ε + τ d + v )
Due to:
q T ( ε + τ d + v ) = q T ( ε + τ d ) + q T v = q T ( ε + τ d ) q ( ε N + τ d ) 0
Based on the above analysis, we can conclude that: L ˙ 0 . When L ˙ 0 , q 0 , according to LaSalle’s invariance principle, the closed-loop system is asymptotically stable, when t , q 0 and the tracking error e corresponds to zero.

3. Results

3.1. Identification of Human–Machine Coupling Parameters

The Particle Swarm Optimization (PSO) algorithm is used to optimize the excitation trajectory, with the chosen parameters of ω = 0.2 π , n ¯ = 2 and θ 0 = [ 0.611 , 0.925 ] T . To sample Equation (11) in MATLAB, the sampling regression matrix A ¯ is obtained. The condition number of the sampling regression matrix is calculated. By utilizing the PSO algorithm, the parameters a i , k and b i , k can be optimized to obtain the desired results. The optimization results based on the PSO are as follows: a 1 , 1 = 0.0229 , b 1 , 1 = 0.0006 , a 1 , 2 = 24 . 0422 , b 1 , 2 = 17 . 9290 , a 2 , 1 = 13 . 0688 , b 2 , 1 = 22.6565 , a 2 , 2 = 3 . 0457 , b 2 , 2 = 15 . 1847 , and the fitness at this moment is 11.6060.
In the data sampling experiment, a model-free PD controller and the UEXO-I Lower Limb Exoskeleton Prototype is used. The optimized excitation trajectory is input into the computer to enable the exoskeleton to swing according to the optimized trajectory. In the first step of the experiment, the exoskeleton swings alone under the control of the PD controller, without any participation of experimental personnel. This step aims to separate the friction parameters. As seen in Figure 2, two servo motor actuators (GDM1-100N2/120N2) and two motor drivers (Elmo-G-SOLHOR15/100EE) drove the exoskeleton’s thigh and shank. Two absolute encoders (INC-4-150 and INC-3-125) were used to measure the rotational angles of the robot’s hip and knee, and four 3D force sensors (JNSH-2-10kg-BSQ-12) were used to detect the coupling forces.
In the second step, the experimental subject, an 88 kg/1.73 m healthy adult, is driven by the exoskeleton without actively exerting force. The subject was standing during the trial, enabling the entire leg to swing and the knee to be in a natural posture. This step is used to determine the inertia parameters of the human body and exoskeleton. Ethical statements: The subject gave informed consent for inclusion before their participation in the study. The study was conducted in accordance with the Declaration of Helsinki, and protocol was approved by the Ethics Committee of the University of Electronic Science and Technology of China (106142023021725090). The parameters obtained from the human–machine identification experiment are shown in Table 3.
The root mean square error between the actual sampled torque of the hip joint and the estimated torque obtained by substituting the identified parameters into Equation (2) is 6.9395. Similarly, the root mean square for the knee joint is 10.2776.

3.2. Simulation Experiment of RBF Control

The simulation of RBF control was conducted using Simulink. The required parameters for the neural network were obtained from reference [28]. The desired trajectory settings were obtained from reference [29]. Furthermore, considering that commercially available motors generally provide a maximum torque of exceeding 200 (N·m), the maximum torque output of the controller was set to 200 (N·m).
In the control of the RBF neural network, the parameter Kv amplifies the error function, and the selection of an appropriate Kv directly affects the control effect. The influence of Kv on the maximum absolute error of the system after stabilization is shown in Figure 3.
In Figure 3, it can be observed that increasing Kv initially leads to a rapid decrease in the maximum absolute error of the hip and knee joints. It can be observed that the error decreases rapidly when Kv is between 100 and 200. However, once Kv increases to 300, the rate of error reduction becomes very slow. Increasing Kv further from 300 to 500 does not result in an error reduction of less than 0.01 (rad). However, once Kv exceeds 300, the error reduction becomes very slow. Additionally, in the simulation experiment it was discovered that when Kv is too large the system’s noise is amplified, resulting in poor system performance. When selecting Kv, it is important to ensure that the system’s error does not decrease sharply with increasing Kv and that the system is not affected by noise. Therefore, Kv = 300 was chosen.
Suitable gains were selected for RBF control and simulation verification was concluded in MATLAB. The comparison between the desired trajectory and the actual tracking trajectory of the system under the RBF neural network control is shown in Figure 4.
In Figure 4, it can be observed that the RBF neural network control can track the desired trajectory of the system. From the graph, it is evident that the RBF control can quickly track the desired trajectory in less than 1 s at the initial stage. The errors between the actual trajectory and the desired trajectory under the RBF neural network control are shown in Figure 5.
In Figure 5, it can be observed that the control performance of the hip joint is superior to that of the knee joint. In the control of the hip joint, the RBF neural network control shows almost no error after reaching a steady state, with the error approaching zero. This is because in the RBF neural network control, the network adaptive law is updated at each sampling instance, enabling dynamic adjustment of the system and effectively compensating for uncertainties in the exoskeleton model. The unknown uncertainties—including noise—that can be approximated by the RBF network come from the actual environment [28]. Therefore, the RBF neural network control demonstrates superior performance in tracking the desired trajectory.
The torque required for tracking the desired trajectory in the RBF neural network control is shown in Figure 6. At the beginning stage, in order to track the desired trajectory quickly, the controller requires a large amount of torque. From the graph, it can be observed that the maximum torque in the RBF control reaches the set upper limit of 200 (N·m). If no limiting control is applied at this point, the excessively large pulses may pose certain risks to the human body and the exoskeleton.

4. Discussion and Conclusions

This study focuses on the exoskeleton human–machine system. To enhance the control effectiveness of the model-based control, an identification experiment is used to obtain the inertial parameters of the human–machine system. Based on this, a simulation is made on the control effects of the RBF adaptive control. In the identification stage, firstly, the exoskeleton human–machine mathematical model is linearized and the excitation trajectory is optimized to use the particle swarm optimization algorithm. Secondly, experimental data is collected to obtained the identification parameters required for the mathematical model. Then, the inertial parameters of the human–machine model are obtained through the least square method, which helps determine the parameters for subsequent simulation experiments. Subsequently, MATLAB is used to simulate and compare the exoskeleton human–machine model. The relationship between the maximum error after stabilization and the control parameters are analyzed to select to the control gain parameters for RBF neural network control, optimizing the control for the control method. The simulation is conducted using MATLAB’s S-function, and the comparative graphs of the desired and actual trajectories of the exoskeleton’s hip and knee joints, error comparison, and torque comparison are obtained. Based on the dynamic compensation ability of RBF in handing system uncertainties, where control is based on the system error and its derivative, comprehensive analysis of the simulation results indicates that RBF provides desirable trajectory tracking control for the exoskeleton when an appropriate parameter is selected. The error decreases when increasing the RBF parameter Kv within a specified range. However, this reducing impact becomes smaller as the parameter values rise. From the results, we can judge that RBF achieve controlling consequences with fast convergence as well as small errors. In addition, the control performance of the hip joint is superior to that of the knee joint.
Substantial progress has been made in the research of lower limb exoskeleton robot control over the past few decades. However, whether the control strategies of lower limb exoskeletons can more effectively stimulate the functional recovery of patients remains an unresolved question. Future research should focus on structured and standardized studies aimed at identifying the relationship between control strategies and a set of core clinical outcome measures, taking into account the impact of participants’ initial impairment level and training intensity [30]. In this study, we assumed that the patient does not exert any force and is completely driven by the exoskeleton, without considering the patient’s intentions. In future research, we will take the patient’s intentions into account and conduct in-depth investigations into the influencing factors. We plan to carry out further research on active control in the future.

Author Contributions

Conceptualization, X.P. and Y.Y.; methodology, X.P.; software, S.Z.; formal analysis, M.C.; writing—original draft preparation, X.P. and S.Z.; writing—review and editing, Y.Y.; visualization, Y.Y.; supervision, Y.Y.; project administration, Y.Y.; funding acquisition, Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grants No. 12072068, 11872147, 11932015, and 52175046), and the Sichuan Science and Technology Program (Grant No. 2023YFG0050).

Data Availability Statement

The numerical and experimental data sets generated and analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The expression for M ( θ ) , C ( θ , θ ˙ ) , G ( θ ) and τ f i c in the dynamic model for human–machine coupling are as follows:
M = [ M 11 M 12 M 21 M 22 ] , C = [ C 11 C 12 C 21 C 22 ] , G = [ G 1 G 2 ] , M = [ X 1 + 2 X 3 L 1 cos ( θ 2 ) X 2 + X 3 L 1 cos ( θ 2 ) X 2 + X 3 L 1 cos ( θ 2 ) X 2 ] C = [ 2 X 3 L 1 θ ˙ 2 sin ( θ 2 ) X 3 L 1 θ ˙ 2 sin ( θ 2 ) X 3 L 1 θ ˙ 1 sin ( θ 2 ) 0 ] G = [ X 3 g sin ( θ 1 + θ 2 ) + X 4 g sin ( θ 2 ) X 3 g sin ( θ 1 + θ 2 ) ] τ f i c = [ K S 1 sgn ( θ ˙ 1 ) + K m 1 θ ˙ 1 K S 2 sgn ( θ ˙ 2 ) + K m 2 θ ˙ 2 ] X 1 = m 1 c 1 2 + I 1 + m 2 L 1 2 + m 1 c 2 2 + I 2 X 2 = m 2 c 2 2 + I 2 X 3 = m 2 c 2 X 4 = m 1 c 1 + m 2 L 1
In the above equation, subscript 1 represents the parameters of the thigh and subscript 2 represents the parameters of the shank. By measuring the length of the exoskeleton’s thigh, L1 = 0.38952. In this context, g represents the gravitational constant.
For the linear form of the dynamic model, the specific expression of the regression matrix A can be represented as:
A = [ A 11 A 12 A 13 A 14 A 15 A 16 A 17 A 18 A 21 A 21 A 22 A 23 A 24 A 25 A 26 A 27 ] A 11 = θ ¨ 1 A 12 = θ ¨ 2 A 13 = L 1 ( 2 θ ¨ 1 cos ( θ 2 ) + θ ¨ 2 cos ( θ 2 ) 2 θ ˙ 1 θ ˙ 2 sin ( θ 2 ) θ ˙ 2 2 sin ( θ 2 ) ) + g sin ( θ 1 + θ 2 ) A 14 = g sin θ 1 A 15 = sgn ( θ ˙ 1 ) A 16 = θ ˙ 1 A 17 = A 18 = A 21 = A 24 = A 25 = A 26 = 0 A 22 = θ ¨ 1 + θ ¨ 2 A 23 = L 1 ( θ ¨ 1 cos ( θ 2 ) + θ ˙ 1 2 sin ( θ 2 ) ) + g sin ( θ 1 + θ 2 ) A 27 = sgn ( θ ˙ 2 ) A 28 = θ ˙ 2
Similarly, the specific form of the parameter vector to be identified can be expressed as:
X = [ X 1 X 2 X 3 X 4 K S 1 K m 1 K S 2 K m 2 ] T

References

  1. Shi, D.; Zhang, W.; Zhang, W.; Ding, X. A Review on Lower Limb Rehabilitation Exoskeleton Robots. J. Chin. Mech. Eng. 2019, 32, 74. [Google Scholar] [CrossRef]
  2. Bao, R.R.; Feng, J.; Lu, Y.N.; Shen, Q.J.; Deng, D.Z. Research progress of lower limb dynamic exoskeleton robot. Sci. Technol. Innov. Her. 2019, 16, 76–77. [Google Scholar]
  3. Lv, Y.; Fang, H.B.; Xu, J.; Ma, J.; Wang, Q.; Zhang, X. Dynamic modeling and analysis of the lower limb prosthesis four-bar linkage prosthetic knee. Chin. J. Theor. Appl. Mech. 2020, 52, 1157–1173. [Google Scholar]
  4. Kalita, B.; Narayan, J.; Dwivedy, S.K. Development of active lower limb robotic-based orthosis and exoskeleton devices: A systematic review. Int. J. Soc. Robot. 2021, 13, 775–793. [Google Scholar] [CrossRef]
  5. Zhang, X.; Zan, Y.; Jing, W. Robotics in Lower-Limb Rehabilitation after Stroke. Behav. Neurol. 2017, 2017, 3731802. [Google Scholar] [CrossRef]
  6. Pamungkas, D.S.; Caesarendra, W.; Soebakti, H.; Analia, R.; Susanto, S. Overview: Types of Lower Limb Exoskeletons. Electronics 2019, 8, 1283. [Google Scholar] [CrossRef]
  7. Hussain, F.; Goecke, R.; Mohammadian, M. Exoskeleton robots for lower limb assistance: A review of materials, actuation, and manufacturing methods. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2021, 235, 1375–1385. [Google Scholar] [CrossRef]
  8. Sun, Y.; Tang, Y.; Zheng, J.; Dong, D.; Chen, X.; Bai, L. From sensing to control of lower limb exoskeleton: A systematic review. Annu. Rev. Control 2022, 53, 83–96. [Google Scholar] [CrossRef]
  9. Choi, B.; Seo, C.; Lee, S.; Kim, B.; Kim, D. Swing Control of a Lower Extremity Exoskeleton Using Echo State Networks. Ifac Pap. 2017, 50, 1328–1333. [Google Scholar] [CrossRef]
  10. Baud, R.; Manzoori, A.R.; Ijspeert, A.; Bouri, M. Review of control strategies for lower-limb exoskeletons to assist gait. J. Neuro Eng. Rehabil. 2021, 18, 119. [Google Scholar] [CrossRef]
  11. Chen, J.; Fan, Y.; Sheng, M.; Zhu, M. Optimized Control for Exoskeleton for Lower Limb Rehabilitation with Uncertainty. In Proceedings of the 2019 31ST Chinese Control and Decision Conference (CCDC2019), Nanchang, China, 3–5 June 2019. [Google Scholar]
  12. Tijjani, I.; Kumar, S.; Boukheddimi, M. A Survey on Design and Control of Lower Extremity Exoskeletons for Bipedal Walking. Appl. Sci. 2022, 12, 2395. [Google Scholar] [CrossRef]
  13. Huang, P.; Yuan, W.; Li, Q.; Feng, Y. Neural network-based optimal control of a lower-limb exoskeleton robot. In Proceedings of the 2021 6th IEEE International Conference on Advanced Robotics and Mechatronics (ICARM), Chongqing, China, 31 December 2021. [Google Scholar]
  14. Zhang, J.Y.; Cao, D.G.; Cao, J.X.; Zhang, T.C.; Wang, J.S.; Cheng, X. Neural network sliding mode control for trajectory tracking of lower limb exoskeleton robot. Process Autom. Instrum. 2021, 42, 47–52. [Google Scholar]
  15. Luo, D.J.; Gao, X.S.; Li, J.; Zhang, P.F. Trajectory tracking control of lower limb exoskeleton rehabilitation robot. Ordnance Ind. Autom. 2020, 39, 87–91. [Google Scholar]
  16. Gilbert, M.; Zhang, X.; Yin, G. Modeling and design on control system of lower limb rehabilitation exoskeleton robot. In Proceedings of the International Conference on Ubiquitous Robots and Ambient Intelligence, Xian, China, 19–22 August 2016. [Google Scholar]
  17. Chen, Y.; Liu, J.; Wang, H.; Pan, Z.; Han, S. Model-free based adaptive RBF neural network control for a rehabilitation exoskeleton. In Proceedings of the 2019 Chinese Control and Decision Conference (CCDC), Nanchang, China, 3–5 June 2019. [Google Scholar]
  18. Liu, D.; Tang, Z.; Pei, Z. The motion control of lower extremity exoskeleton based on RBF neural network identification. In Proceedings of the IEEE International Conference on Information and Automation 2015, Lijiang, China, 8–10 August 2015. [Google Scholar]
  19. Bingül, Z.; Karahan, O. Dynamic identification of Staubli RX-60 robot using PSO and LS methods. Expert Syst. Appl. 2011, 38, 4136–4149. [Google Scholar] [CrossRef]
  20. Wu, J.; Wang, J.; You, Z. An overview of dynamic parameter identification of robots. Robot. Comput. Integr. Manuf. 2010, 26, 414–419. [Google Scholar] [CrossRef]
  21. Gao, G.; Sun, G.; Na, J.; Guo, Y.; Wu, X. Structural parameter identification for 6 DOF industrial robots. Mech. Syst. Signal Process. 2017, 113, 145–155. [Google Scholar] [CrossRef]
  22. Maxime, G.; Briot, S. Dynamic Parameter Identification of a 6 DOF Industrial Robot using Power Model. In Proceedings of the IEEE International Conference on Robotics and Automation ICRA, Karlsruhe, Germany, 6–10 May 2013. [Google Scholar]
  23. Kapandji, A.I. The Physiology of the Joints; Churchill Livingstone: Edinburgh, UK, 2011. [Google Scholar]
  24. Lei, L.; Li, J.; Wu, Q. Gait Tracking of Lower Limb Exoskeleton Based on RBF Neural Network Adaptive Control. J. Guangxi Univ. Sci. Technol. 2021, 32, 7. [Google Scholar]
  25. Hu, H.; Hu, L.; Liu, Y.; Cao, W.; Chen, C. Research on a Flexible Lower Limb Exoskeleton Control Strategy. J. Instrum. Meas. 2020, 41, 184–191. [Google Scholar]
  26. Song, S.; Zhang, X.; Tan, Z. RBF Neural Network Based Sliding Mode Control of a Lower Limb Exoskeleton Suit. J. Mech. Eng. 2014, 60, 437–446. [Google Scholar] [CrossRef]
  27. Sun, Z.; Li, F.; Wang, G.; Liu, Y.; Lian, Y.; Liu, K. A novel RBF neural network-based iterative learning control for lower limb rehabilitation robot with strong robustness. In Proceedings of the Chinese Control Conference, Guangzhou, China, 27–30 July 2019. [Google Scholar]
  28. Liu, J. Robot Control System Design and MATLAB Simulation; Tsinghua University Press: Beijing, China, 2008; pp. 27–31. [Google Scholar]
  29. Chen, Z.; Guo, Q.; Xiong, H.; Jiang, D.; Yan, Y. Control and Implementation of 2-DOF Lower Limb Exoskeleton Experiment Platform. Chin. J. Mech. Eng. Engl. Ed. 2021, 34, 17. [Google Scholar] [CrossRef]
  30. de Miguel-Fernández, J.; Lobo-Prat, J.; Prinsen, E.; Font-Llagunes, J.M.; Marchal-Crespo, L. Control strategies used in lower limb exoskeletons for gait rehabilitation after brain injury: A systematic review and analysis of clinical effectivenes. J. NeuroEngineering Rehabil. 2023, 20, 23. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Model of Lower Extremity Exoskeleton and Human Body.
Figure 1. Model of Lower Extremity Exoskeleton and Human Body.
Actuators 12 00353 g001
Figure 2. The exoskeleton used for experiment.
Figure 2. The exoskeleton used for experiment.
Actuators 12 00353 g002
Figure 3. Relationship between Kv and the maximum absolute value.
Figure 3. Relationship between Kv and the maximum absolute value.
Actuators 12 00353 g003
Figure 4. Comparison of trajectories between the desired trajectory and RBF neural network.
Figure 4. Comparison of trajectories between the desired trajectory and RBF neural network.
Actuators 12 00353 g004
Figure 5. Error in RBF neural network control.
Figure 5. Error in RBF neural network control.
Actuators 12 00353 g005
Figure 6. Control torque on RBF neural network control.
Figure 6. Control torque on RBF neural network control.
Actuators 12 00353 g006
Table 1. Symbolic representation of physical parameters for the coupled model.
Table 1. Symbolic representation of physical parameters for the coupled model.
Meaning of ParametersNotation
Hip joint angle/Knee joint angle θ 1 / θ 2
Thigh length/Shank lengthL1/L2
Thigh mass/Shank massm1/m2
Thigh center of mass/Shank center of massc1/c2
Thigh moment of inertia/Shank moment of inertiaI1/I2
Table 2. The constraint intervals for the excitation trajectory.
Table 2. The constraint intervals for the excitation trajectory.
Constraint TermsConstraint Intervals
Hip joint angle/ Knee joint angle (rad)[0.0872, 1.2217]/[−1.6232, −0.2269]
Hip/Knee angular velocities (rad/s)[−1.6449, 1.6449]
Hip/Knee angular accelerations (rad/s2)[−5.1677, 5.1677]
Table 3. Identification parameters of human body and exoskeleton system.
Table 3. Identification parameters of human body and exoskeleton system.
ParameterValueParameterValue
X110.3016Hip joint’s coefficient of static friction (Ks1)2.6596
X23.0157Hip joint’s coefficient of viscous friction (Km1)4.14895
X32.2664Knee joint’s coefficient of static friction (Ks2)4.7197
X45.5354Knee joint’s coefficient of viscous friction (Km2)7.3602
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Peng, X.; Zhang, S.; Cai, M.; Yan, Y. Modelling and RBF Control of Low-Limb Swinging Dynamics of a Human–Exoskeleton System. Actuators 2023, 12, 353. https://doi.org/10.3390/act12090353

AMA Style

Peng X, Zhang S, Cai M, Yan Y. Modelling and RBF Control of Low-Limb Swinging Dynamics of a Human–Exoskeleton System. Actuators. 2023; 12(9):353. https://doi.org/10.3390/act12090353

Chicago/Turabian Style

Peng, Xinyu, Shujun Zhang, Mengling Cai, and Yao Yan. 2023. "Modelling and RBF Control of Low-Limb Swinging Dynamics of a Human–Exoskeleton System" Actuators 12, no. 9: 353. https://doi.org/10.3390/act12090353

APA Style

Peng, X., Zhang, S., Cai, M., & Yan, Y. (2023). Modelling and RBF Control of Low-Limb Swinging Dynamics of a Human–Exoskeleton System. Actuators, 12(9), 353. https://doi.org/10.3390/act12090353

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop