Towards Human-like Walking with Biomechanical and Neuromuscular Control Features: Personalized Attachment Point Optimization Method of Cable-Driven Exoskeleton

The cable-driven exoskeleton can avoid joint misalignment, and is substantial alterations in the pattern of muscle synergy coordination, which arouse more attention in recent years to facilitate exercise for older adults and improve their overall quality of life. This study leverages principles from neuroscience and biomechanical analysis to select attachment points for cable-driven soft exoskeletons. By extracting key features of human movement, the objective is to develop a subject-specific design methodology that provides precise and personalized support in the attachment points optimization of cable-driven exoskeleton to achieve natural gait, energy efficiency, and muscle coordination controllable in the domain of human mobility and rehabilitation. To achieve this, the study first analyzes human walking experimental data and extracts biomechanical features. These features are then used to generate trajectories, allowing better natural movement under complete cable-driven exoskeleton control. Next, a genetic algorithm-based method is employed to minimize energy consumption and optimize the attachment points of the cable-driven system. This process identifies connections that are better suited for the human model, leading to improved efficiency and natural movement. By comparing the calculated elderly human model driven by exoskeleton with experimental subject in terms of joint angles, joint torques and muscle forces, the human model can successfully replicate subject movement and the cable output forces can mimic human muscle coordination. The optimized cable attachment points facilitate more natural and efficient collaboration between humans and the exoskeleton, making significant contributions to the field of assisting the elderly in rehabilitation.

1 Introduction

Due to advances in medical technology, improvements in living conditions, increased health awareness, and declining birth rates, the global population is experiencing rapid aging (Organization et al., 2020). As individuals age, they often encounter a gradual decrease in muscle mass and the quantity of muscle fibers, a condition referred to as muscle atrophy. Concurrently, the incidence of stroke, particularly among the elderly, rises significantly, with rates as high as 65% (Lo Coco et al., 2016). This muscle atrophy is typically attributed to various factors, including alterations in the nervous system, hormonal fluctuations, and reduced physical activity levels (Ziaaldini et al., 2017). These factors collectively contribute to the weakening of lower limb strength in older adults, thereby elevating the risk of falls, disability, and mortality (Thompson et al., 2021).

Existing research has demonstrated the high effectiveness of regular exercise in preventing and treating age-related skeletal muscle loss symptoms (Frizziero et al., 2016; Yoo et al., 2018; Shur et al., 2021). Furthermore, physical exercise shows promise in enhancing functional motor recovery by promoting neuroplasticity at various cellular and molecular levels (Pin-Barre and Laurin, 2015; Liang et al., 2021). To facilitate exercise for older adults, and improve their overall quality of life, robot technologies have found applications in the field of elderly care (Alnajjar et al., 2019), including the use of active exoskeleton robots (Li et al., 2021; Xiong et al., 2022) and assistive robots (Pino et al., 2015). To achieve the goal of “rehabilitation training in daily life” and meet the requirements for joint driving capability and movement precision, many internationally developed wearable rehabilitation exoskeletons rely on rigid components such as motors and gear reducers (Carnevale et al., 2018; Nakashima et al., 2020). However, these traditional rigid exoskeleton designs tend to be heavy, increasing the mass and inertia of the user's distal body parts and consequently elevating metabolic costs during movement. Notably, one significant drawback of these rigid exoskeleton designs is the issue of joint misalignment (Yeem et al., 2018). Additionally, wearing such exoskeletons often results in substantial alterations in the average pattern of muscle synergy coordination among participants (Li et al., 2019).

Inspired by human muscles and tendons, flexible drives systems such as pneumatics (Prasad et al., 2022), hydraulics (Glowinski et al., 2020), and cables (Prasad et al., 2022) offer distinct advantages in the design of exoskeletons, including high power density ratio, wide contraction range, rapid response, flexible structure, and safety. Among these options, lighter and simpler cable drives exhibit particular advantages in exoskeletons, as they efficiently replicate the biomimetic properties of human tendon-joints (Xiloyannis et al., 2021). An illustrative example is the cable-driven exoskeleton garment, known as an Exosuit, which employs electric motors to drive hip and ankle joint movements through Bowden cables, enabling flexible actuation for assisting in horizontal walking (De Rossi et al., 2022). Another innovative development is the HitExosuit, which features twisted cable actuators and artificial muscles based on steel cables, supplying essential force to each hip joint during stair climbing, providing assistance torque to the knee joints, which are most crucial during this activity (Zhao et al., 2019). In the context of stair ascent and descent, a soft exoskeleton system with an additional knee joint effectively reduces muscle activity, thus lightening the user's load (Lee et al., 2020). A wearable soft robotic exosuit incorporating a portable cable actuation system is introduced, applying hip abduction torque to potentially mitigate knee osteoarthritis progression by reducing external knee adduction moment (Yang et al., 2022). A portable waist-loaded soft exosuit, concentrated at the waist, effectively aids hip flexion during running, demonstrated through motion flexibility experiments (Chen et al., 2022). In Biao et al. (2023), a novel flexible Bowden cable-driven exoskeleton is introduced, with each joint controlled by a pair of antagonistic muscles. This design utilizes direct lower limb motion feedback to achieve a natural gait during walking.

However, biomechanical investigations have highlighted the importance of biarticular elements in enhancing energy efficiency during gravity-intensive daily activities, signifying the crucial role of biarticular muscles in the future of exoskeleton-driven systems (Sharbafi et al., 2018). An exemplary flexible wearable device in this context is the exoskeleton robot Myosuit, which employs biarticular and coordinated control methods, making it suitable for gravity-intensive daily activities like sitting transfers (Schmidt et al., 2017; Ganesan and Gupta, 2021). BATEX, a new exosuit, is inspired by human musculoskeletal systems and neural control, to enhance metabolic efficiency, adaptability, and overall effectiveness (Firouzi et al., 2021). Additionally, exosuits utilize flexible drive units with Bowden cables to connect the ankle joint and thigh, delivering active assistance during ankle plantar flexion, consequently reducing metabolic consumption and muscle output (He et al., 2023). These exosuits contribute significantly to assisting the elderly in various tasks across different scenarios. Nevertheless, it is worth noting that the attachment points or routing positions of most exoskeleton suits are typically determined based on experiential and intuitive knowledge, with only a limited number of exosuits being optimized for factors such as stiffness (Wu et al., 2019), cable tension (Grosu et al., 2018), joint torque (Prasad et al., 2023), or metabolic factors (Bardi et al., 2022). In Prasad et al. (2022) a promising design option is proposed for a lower limb rehabilitation exoskeleton. Its identification is based on tracking performance, model requirements, and forces exerted on the limb by its components.

On the other hand, the biomechanical characteristics of the human body encompass not only the interaction of biarticular muscles but also the coordination between monoarticular and biarticular muscles (Junius et al., 2017). Moreover, aging constitutes a systemic reorganization of the neuro-musculo-skeletal system, extending beyond the mere degeneration of muscles (Sleimen-Malkoun et al., 2014). As individuals age, changes in age-related neural circuits can lead to alterations in muscle coordination, prompting compensatory changes to adapt to new conditions (Guo et al., 2022). Fortunately, the aging brain retains the capacity to reorganize activation patterns between neural circuits, accommodating anatomical and physiological changes (Vernooij et al., 2016). As motor tasks become more complex, this expanded activation is typically more pronounced, resulting in the emergence of new muscle synergy patterns, which are closely associated with changes in white matter and neural plasticity within the brain (Singh et al., 2018). Therefore, rehabilitation strategies founded on muscle synergistic interactions hold great promise in inducing genuine recovery and enhancing neuro-muscular coordination (Hong et al., 2021).

In light of their ability to mimic human muscle coordination and facilitate its achievement, exoskeletons hold immense significance in assisting older adults in their rehabilitation journey. Additionally, the selection of attachment points for soft exoskeletons significantly impacts their comfort, energy efficiency, and wearability (Asbeck et al., 2014; Tonazzini et al., 2018). The activation patterns of lower limb muscles exhibit variations across different walking speeds (Arnold et al., 2013). The cable arrangement and attachment points should also be finely tuned to meet the demands of various Activities of Daily Living (ADL) (Joshi et al., 2022). Exosuit-assisted strategies that rely on individual measurements of muscle dynamics in specific tasks, coupled with personalized approaches for achieving natural and complete gait restoration, along with targeted feedback, enhance exosuit performance. Furthermore, these strategies promote the establishment of new neural control centers and stimulate the activation of the entire lower limb muscle group, thereby providing valuable feedback to the cerebral cortex (Biao et al., 2023).

To address the aforementioned requirements effectively, before the implementation of cable-driven exoskeleton system, three key problems need to be tackled:

① The method to extract the biomechanical and neuro-control features from human motion experimental data, providing a more natural and effective motion reference for the cable-driven exoskeleton.

② The cable-driven exoskeleton modeling method considering the dynamics of the elderly human model and cable-driven system as well.

③ The optimization method of cable attachment points to make the cable-driven exoskeleton realize human-like walking motion with biomechanical features (e.g., ZMP, COM, joint angles and torques trajectories) and neuromuscular control features (e.g., muscle output forces and muscle synergy patterns).

With these challenges in mind, this article proposes a subject-specific method for determining the attachment points of a cable-driven exoskeleton. This approach allows the generated forces of the cables to adapt to the biological synergy during motion, thereby facilitating neuroplasticity in patients with neurological impairments and facilitating the restoration of normal walking ability. Additionally, this method holds promise for the development of humanoid musculoskeletal assistive robots, as older adults tend to accept humanoid robots more readily (Choudhury et al., 2018; Feingold Polak et al., 2018). Musculoskeletal robots can assess the muscle strength of older adults during movement and obtain precise data throughout the rehabilitation process, enabling adjustments to rehabilitation plans based on individual circumstances and needs (Farhat et al., 2022).

2 Materials and methods

With the continuous advancement of science and technology, exoskeleton robotic systems have made significant breakthroughs in improving the quality of life and providing assistance. This research aims to further drive the design of subject-specific cable-driven exoskeleton before its mechanical implementation to provide more human-like nature motion to support the fields of elderly human mobility and rehabilitation. The flowchart of the proposed approach is shown in Figure 1.

Figure 1

Figure 1. Subject-specific method for attachment points optimization of cable-driven exoskeleton.

Firstly, subject walking experiments are performed, and the experimental data is analyzed to extract biomechanical features. These features are applied to generate trajectories for the elderly human model, aiming to achieve similar movements under complete cable-driven exoskeleton control. Additionally, the subject's neural-control features such as muscle forces and muscle activation order are calculated using electromyography (EMG) and dynamic data.

Next, based on neuroscience, an analysis of the roles and contributions of lower limb muscles during movement is conducted. Using elderly subject data, a human model is established, and the cable-driven system is determined. Separately, the lower limb biomechanics model for the elderly human model, the dynamics model for the cable-driven system, and the ground reaction force model are constructed to derive the dynamic model of cable-driven exoskeleton. The calculation of cable forces are calculated based on the specific application scenarios in which the human motion are totally activated by the cable-driven system.

Finally, the genetic algorithm (GA) is employed to optimize the attachment points of the cable-driven system, seeking attachment points that are well-suited for the human model, in order to make the system generate human-like walking motion with biomechanical features and neural-control features as well. In the following sections, the detailed materials and methods used in our research are delved into.

2.1 Modeling

Modeling forms the foundation of cable-driven exoskeleton optimization. In the following section, we will provide an in-depth examination of the modeling techniques and methods employed in our research.

2.1.1 Biomechanical and neuromuscular control analysis of human lower limbs' motion

The lower limbs serve as the foundational support structure of the human body, maintaining an upright posture and facilitating movement. These lower limbs consist of three major joints: the hip joint, knee joint, and ankle joint. However, designing a lower limb exoskeleton that strictly adheres to the 14 degrees of freedom found in the human lower limb system is not only challenging but also demands extensive computational resources. According to Park's (2008) research, movements of the human upper limbs, including the arms, have minimal impact on sagittal dynamics. As a result, the human leg is simplified to incorporate 6 degrees of freedom in the sagittal plane, while the upper limbs, trunk, and head are collectively simplified as a single entity.

Human motion involves a complex interplay of more than 40 muscles that are relevant to various movements. However, attempting to incorporate all of these muscles into an exoskeleton system is not a practical endeavor. As previously mentioned, human joint movements are achieved through the complementary functions of monoarticular and biarticular muscles. Monoarticular muscles primarily generate unidirectional forces in the form of tension and are capable of performing pulling actions. To enable bidirectional motion, two monoarticular muscles are typically arranged as flexors and extensors. On the other hand, biarticular muscles, which span two joints, play a crucial role in human walking. They serve several vital functions, including the reduction of limb inertia during fast movements. They also contribute to gait stability, extension of the working range of motion, and optimization of the working point of monoarticular muscles. Additionally, biarticular muscles assist in maintaining limb balance, contribute to the control of leg swings, and have a positive impact on the development of running exercises. The consideration of neuroscience in this context adds to our understanding of the intricate neuromuscular control involved in human motion.

The results of neuroscience research indicate that in the sagittal plane, the knee and hip joints are driven by seven major muscles, which include four monoarticular muscles: iliacus (ILPSO) for hip flexion, gluteus maximus (GMAX) for hip extension, biceps femoris (BFS) for knee flexion, and vastus (VAS) for knee extension. There are also three biarticular muscles involved: rectus femoris (RF) for hip extension and knee flexion, hamstring (HAMS) for hip extension and knee flexion, and gastrocnemius (GAS) for knee flexion and ankle extension, as documented in the research by Hwang et al. (2018). Furthermore, two monoarticular muscles that drive the ankle joint have been selected for inclusion: the soleus (SOL) for ankle plantar flexion and the tibialis anterior (TA) for ankle dorsiflexion. These muscles play critical roles in controlling foot movements and maintaining balance during human locomotion, as outlined in the study by Di Giulio et al. (2009). These nine major muscles as the mimicry targets for the cable-driven exoskeleton are selected to drive the human model.

2.1.2 Lower limb modeling of elderly: an example based on average data

As a result of the aging process, the bodies of elderly individuals undergo regressive changes, leading to noticeable alterations in their body size, particularly a decrease in height compared to their younger years. These aging-related changes also affect various other body dimensions. Consequently, the height dimensions of elderly individuals cannot be simply equated to those of the local average young population. For this study, we have utilized an example of average dimensions of male elderly individuals, which have been sourced from the research conducted in Hu et al. (2007). These dimensions include an average height of 1,657 mm, a weight of 67 kg, a thigh length of 448 mm, a calf length of 401 mm, an ankle height of 66 mm, a foot length of 241 mm, and a foot width of. These measurements provide a more accurate representation of the physical characteristics of elderly individuals. To further enhance the accuracy of our analysis, the estimation of the center of mass and moments of inertia for each body segment is based on Cheng et al. (2000). For a comprehensive breakdown of these parameters, please refer to Table 1.

Table 1

Table 1. Parameters of human body model for elderly people.

Once the initial human body model has been established, the next step involves performing a kinematic analysis of the human body. In this analysis, the human model is regarded as a free-system configuration, and a set of generalized coordinates is defined as shown in Equation (1):

$\begin{array}{ll}\text{q}={\left[x,y,q_1,q_2,q_3,q_4,q_5,q_6,q_7\right]}^T& (1)\end{array}$

In this set of coordinates, x,y,q₁~q₇ represent generalized coordinates as shown in Figure 2A. Additionally, there are parameters denoted as l₁, l₂, l₃, l₄, l₅ that correspond to specific lengths in the human model: the length of the trunk, the length of the thigh, the length of the shank, the height from the foot to the ankle joint, the length of the heel, and the length of the toe. These parameters are illustrated in Figure 2A. The DH parameters are further detailed in Table 2. With these DH parameters and the known joint angles, it becomes feasible to calculate the kinematics of the human model, including the position, velocity, and acceleration of each component of the human body. Subsequently, the dynamic model will be analyzed in further detail, building upon the obtained kinematic information.

Figure 2

Figure 2. Human model and simplified cable-driven exoskeleton model. (A) Coordination of human model. (B) Simplified model of cable-driven exoskeleton.

Table 2

Table 2. D-H parameters in left leg.

2.1.3 Cable-driven system modeling

According to the previous analysis, it is understood that nine major muscles of human leg play major role in the motion of sagittal plane. Therefore, in this chapter, nine cables have been set up, each designed to mimic the functions of these muscles, in order to drive the elderly human model to move like a healthy individual. The cable distribution is depicted in Figure 2B. The tension of the i_th cable is denoted as f_ci, which is a scalar value. We set the attachment point closer to the trunk as the starting point, and the other as the ending point. The starting and ending points of the i_th cable on the human model are set as r_si and r_ei. The forces applied at r_si and r_ei are represented as f_si, f_ei. The cables are configured as a massless rigid cable model with infinite stiffness. The properties to be defined for this model include a set of attachment locations and the minimum and maximum force bounds for the cables (Eden et al., 2018).

Consistent with the prior analysis, the resolution of cable tensions is profoundly affected by the selection of attachment points. Prior to identifying the optimal attachment points, we initially established the initial positions of cable attachment points based on the attachment points of lower limb muscles using a linear mapping. The initial attachment points are determined by linearly mapping human muscle attachment points to the attachment points in the human model. Muscle origin and insertion data are collected from eight males with an average height of 175 (±8.36) cm and an average age of 38.75 (±10.61) years (Shan, 2003). The distance between the hip joint and the knee joint averages 40.82 (±2.35) cm, while the distance between the knee joint and the ankle joint averages 38.00 (±2.37) cm.

Due to the fact that cables can only provide tension, the minimum force is set to 0. The maximum force values are determined by the driver. In this paper, we plan to use MYO-muscle to mimic the functions of muscles through the cables (Marques et al., 2013; Eden et al., 2018). Other structures have relatively small masses compared to the human model, and since this paper primarily focuses on cable forces and attachment point optimization, the modeling of other mechanical structures in the exoskeleton is intentionally disregarded. To calculate the cable forces in the exoskeleton, the cable dynamic model will be discussed in detail in the next section.

2.1.4 Dynamic model

The cable-driven exoskeleton including human model is a system with multiple degrees of freedom, and it is crucial to establish its dynamic model and conduct relevant analyses. The Lagrange equation is employed to establish the system's dynamics equations:

$\begin{array}{ll}{F}_i=\frac{d}{dt}\frac{\partial L}{\dot{\partial q_i}}-\frac{\partial L}{\partial q_i}(i=1,2,3...n)& (2)\end{array}$

where q_i represents the generalized coordinates of joint i shown in Figure 2A, $\dot{q_i}$ is the corresponding velocity, F_i is the nonconservative force acting on the joint i of human model, and n is the number of links of human model. The Lagrangian function L is defined as the difference between the kinetic energy K and the potential energy V of the cable-driven exoskeleton including the human model: L = K−V.

The center of mass of the i_th link of human model in the generalized coordinate system is represented as $p_i={\left[x_i(q),y_i(q)\right]}^T$, which are obtained through direct kinematics. The velocity and angular velocity of the i_th link in the generalized coordinate system are represented as: $\dot{p_i}=J_{v_i}\dot{q}$ and ${\omega }_i=J_{{\omega }_i}\dot{q}$, J_vi and J_ωi are, respectively, the velocity Jacobian matrix and angular velocity Jacobian matrix. So the kinetic enery and potential energy can be calculated as $K=\sum _{i=1}^n(\frac{1}{2}{\dot{p_i}}^T{m}_i\dot{p_i}+\frac{1}{2}{{\omega }_i}^T{R}_i{I}_i{R}_i^T{\omega }_i)=\frac{1}{2}\dot{q^T}{\sum }_{i=1}^n(m_i{J}_{v_i}^T{J}_{v_i}+J_{{\omega }_i}^T{R}_i{I}_i{R}_i^T{J}_{{\omega }_i})\dot{q}$, the potential energy $V=\sum _{i=1}^n{m}_{i}g{y}_i$. I_i is moments of inertia and R_i is the attitude transformation matrix of the link coordinate system relative to the world coordinate system. Clearly, K is a function of q and $\dot{q}$, while V is a function of q. Then the Equation (2) can be deduced as

$\begin{array}{ll}{F}_i=\frac{d}{dt}\frac{\partial K}{\dot{\partial q_i}}-\frac{\partial K}{\partial q_i}+\frac{\partial V}{\partial q_i}(i=1,2,3...n)& (3)\end{array}$

The first term in the Equation (3) can be expanded as

$\begin{array}{ll}\frac{d}{dt}\frac{\partial K}{\dot{\partial q_i}}=\frac{\partial }{\partial \dot{q_i}}(\frac{\partial K}{\partial \dot{q_i}})\ddot{q_i}+\frac{\partial }{\partial \dot{q_i}}(\frac{\partial K}{\partial q_i})\dot{q_i}& (4)\end{array}$

We set $D=\sum _{i=1}^n(m_i{J}_{v_i}^T{J}_{v_i}+J_{w_i}{R}_i{I}_i{R}_i^T{J}_{w_i})$, the Equation (4) can be further written as shown in Equation (5)

$\begin{array}{ll}\frac{d}{dt}\frac{\partial K}{\stackrel{.}{\partial q_i}}={{\displaystyle \sum }}_{j=1}^n{D}_{ij}{\ddot{q}}_i+{{\displaystyle \sum }}_{j,k=1}^n\frac{\partial D_{kj}}{\partial q_i}{\dot{q}}_j{\dot{q}}_k& (5)\end{array}$

The second term in the Equation (3) can be expanded as shown in Equation (6)

$\begin{array}{ll}\frac{\partial K}{\partial q_i}=\frac{1}{2}{{\displaystyle \sum }}_{j,k=1}^n\frac{\partial D_{kj}}{\partial q_i}{\dot{q}}_j{\dot{q}}_k& (6)\end{array}$

Set ${\Gamma }_{ijk}=\frac{1}{2}(\frac{\partial D_{ij}}{\partial k}+\frac{\partial D_{ik}}{\partial j}+\frac{\partial D_{kj}}{\partial i})$

The analytical expression of the dynamic equation is shown in Equation (7)

$\begin{array}{ll}D(q)\ddot{q}+C(q,\dot{q})\dot{q}+G(q)=F& (7)\end{array}$

D represents the mass-inertia matrix, while $C_{i,j}=\sum _{k=1}^n{\Gamma }_{ijk}\dot{q_k}$ and $G=\frac{\partial V}{\partial q}$ represent the centrifugal-coriolis and gravitational terms, respectively. The nonconservative force F can be expressed as shown in Equation (8):

$\begin{array}{ll}F=\tau +J_{grf}^T·f_{grf}+J_{\text{c}}^T·f_c& (8)\end{array}$

τ represents the human joint drive torques. In our application case of the exoskeleton, we suppose that elder patient's lower limb movement is entirely driven by cable-driven exoskeleton, the patient joint drive torques τ can be set to zero. f_grf represents the ground reaction force, J_grf is the Jacobian matrix that relates the ground reaction force to the joints. f_c denotes the cable force vector, which should be positive. Additionally, J_c is the Jacobian matrix that establishes the relationship between the cable force vector and the joints. f_grf, J_grf, f_c, and J_c will be elaborated in detail in the following sections.

2.1.4.1 Dynamic model of cable-driven system

To calculate the Jacobian matrix, J_c, the following steps are taken. As shown in Figure 2B, r_si represents the attachment starting point of the i_th cable, and f_si represents the tension of the i_th cable acting on the attachment starting point. Similarly, r_ei represents the attachment ending point of the i_th cable, and f_ei represents the tension of the i_th cable acting on the attachment ending point. The unit vectors of the i_th cable's attachment points are defined by r_si and r_ei. Additionally, the cable tension is represented as f_ci, which is a scalar. The forces acting on r_si and r_ei can be expressed as $f_{si}=f_{ci}·\widehat{I_i}$ and $f_{ei}=-f_{ci}·\widehat{I_i}$, respectively, where $\widehat{I_i}$ is a unit direction vector as shown in Equation (9).

$\begin{array}{ll}{\widehat{\text{I}}}_i=\frac{r_{\text{ei}}-r_{si}}{\sqrt{{(r_{ei}-r_{si})}^T(r_{ei}-r_{si})}}& (9)\end{array}$

Component force F_ij of f_si and f_ei acting on j_th joint are expressed as shown in Equation (10):

$\begin{array}{ll}{F}_{ij}={(\frac{\partial r_{si}}{\partial q_j})}^T{f}_{si}+{(\frac{\partial r_{ei}}{\partial q_j})}^T{f}_{ei}={\left(\frac{\partial r_{si}}{\partial q_j}-\frac{\partial r_{ei}}{\partial q_j}\right)}^T·{\widehat{\text{I}}}_i·f_{ci}& (10)\end{array}$

So the Jacobin matrix ${J_c}^T$ is expressed as shown in Equation (11)

$\begin{array}{ll}{J}_c^T=\left[\begin{array}{ccc}{\left(\frac{\partial r_{e1}-\partial r_{s1}}{\partial q_1}\right)}^T·{\widehat{\text{I}}}_1& \cdots & {\left(\frac{\partial r_{ei}-\partial r_{si}}{\partial q_1}\right)}^T·{\widehat{\text{I}}}_i\\ ⋮\\ {\left(\frac{\partial r_{e1}-\partial r_{s1}}{\partial q_j}\right)}^T·{\widehat{\text{I}}}_1& \cdots & {\left(\frac{\partial r_{ei}-\partial r_{si}}{\partial q_j}\right)}^T·{\widehat{\text{I}}}_i\end{array}\right]& (11)\end{array}$

2.1.4.2 Estimation of ground reaction force

The ground reaction force on the standing leg during the single supporting phase is calculated based on both inertial force and gravity as shown in Equation (12),

$\begin{array}{ll}\{\begin{array}{ll}{F}_{grf}\hfill & ={\displaystyle \sum (m_i{\ddot{p}}_i-m_{i}g)}\hfill \\ M_{grf}\hfill & ={\displaystyle \sum (\frac{d}{dt}(}I\omega )-(m_i{\ddot{p}}_i-m_{i}g)\times r_i)\hfill \end{array}& (12)\end{array}$

In the double-supporting phase, the ground reaction forces are distributed between both legs. The distribution of this load is approximated using a weighted linear relationship between the two contact points on the foot and the position of the Zero Moment Point (ZMP) (Xiang et al., 2009). As shown in Figure 3, F_fore−grf, M_fore−grf and F_rear−grf, M_rear−grf represent the ground reaction forces on the forefoot and rearfoot, respectively. d₁ and d₂ represent the vector from the forefoot and rearfoot contact points to the ZMP, respectively. The trajectory of ZMP is determined based on human motion, which will be elaborated on in the following section.

$\begin{array}{l}{F}_{fore-grf}=\frac{\left|\right|d_2\left|\right|}{\left|\right|d_1\left|\right|+\left|\right|d_2\left|\right|}{F}_{grf}\\ M_{fore-grf}=\frac{\left|\right|d_2\left|\right|}{\left|\right|d_1\left|\right|+\left|\right|d_2\left|\right|}{M}_{grf}+d_1\times F_{fore-grf}\\ F_{rear-grf}=\frac{\left|\right|d_1\left|\right|}{\left|\right|d_1\left|\right|+\left|\right|d_2\left|\right|}{F}_{grf}\\ M_{rear-grf}=\frac{\left|\right|d_1\left|\right|}{\left|\right|d_1\left|\right|+\left|\right|d_2\left|\right|}{M}_{grf}+d_2\times F_{rear-grf}& (13)\end{array}$

Figure 3

Figure 3. Partition of ground reaction forces.

The interaction force between the human model and the ground can be obtained according to the force and moment on each leg of the human model in Equation (13) (Xiang et al., 2009);

During the single supporting phase, the contact point between the foot and the ground is modeled by Φ(q) as shown in Equation (14).

$\begin{array}{ll}{J}_{grf}=\frac{\partial \Phi (q)}{\partial q}& (14)\end{array}$

During the double supporting phase, the contact points between the foot and the ground are Φ_fore−foot(q) and Φ_rear−foot(q).

2.1.5 Calculation of cable forces

The resulting dynamic model becomes a redundant dynamical system, leading to multiple solutions during the calculation of cable force. Consequently, the inverse dynamic problem is typically formulated as an optimization problem. Taking into account the constraints of positive cable force, the optimal set of cable forces can be determined by solving Equation (15).

$\begin{array}{l}{f}_c^{*}=\mathrm{argmin}Q(f_c)\\ \text{ s.t. }D(q)\ddot{q}+C(q,\stackrel{.}{q})\stackrel{.}{q}+G(q)-J_{grf}^T{f}_{grf}=J_{\text{c}}^T·f_c\\ f_{\text{min}}\le f_c\le f_{\text{max}}& (15)\end{array}$

The minimum and maximum bounds on the cable force are represented by the vectors f_min and f_max, respectively. A typical objective function $Q(f_c)=\sum _m^i{f}_i^2$ is selected to achieve a desired goal (Lau et al., 2016). The optimization is solved in interior point method. In next section, the human walking study is conducted in order to complete model analysis.

2.2 Human walking experiment and biomechanics features extraction

The daily activities require a wide range of lower limb movements, including walking, climbing stairs, bending, standing and sitting, walking, and changing clothes. Among these activities, balance-recovery stepping is a particularly common movement, and is considered a typical daily activity (Liu and Lockhart, 2009). Therefore, this chapter primarily focuses on generating motion data for the human model, specifically pertaining to the balance-recovery stepping motion observed in healthy individuals. Generating biomimetic gait is a complex proposition, involving not only kinematics but also dynamics and biomechanics. Many researchers have established a strong correlation between specific characteristics of human gait and aspects of walking stability and energy efficiency (Chen et al., 2017). Consequently, experiments on human walking was been designed to explore these characteristics through motion analysis. In addition, the muscle forces of the subject was also be analyzed to serve as an evaluation criterion for the performance of the human model driven by the exoskeleton.

2.2.1 Human walking experiment

The experimental environment is depicted in Figure 4A. It includes three force plates (40 × 60 cm, AMTI, USA) embedded in the floor to record ground reaction forces. Over 30 reflective skin markers (9-mm diameter) were strategically placed at joint locations and other significant anatomical landmarks. The selection of marker locations aligned with those in the Gait2354 Oensim model and includes markers at the acromion process of the scapula, anterior superior iliac spine, front thigh, lateral thigh, medial knee (medial femoral condyle), lateral knee (lateral femoral condyle), front shank, lateral shank, medial malleolus, lateral malleolus, heel (calcaneus), head of the fifth metatarsal, head of the first metatarsal, and tip of the first toe. These markers were placed on both the right and left sides of the body, as well as on the top of the head, top of the sacrum, and sternum shown in Figure 4B.

Figure 4

Figure 4. Experimental setup. (A) Motion capture system. (B) The position of the reflective markers.

To collect kinematic data, a twelve-camera motion capture system (VICON T40S, Oxford, UK) was employed to capture data at a 100 Hz rate, while force-plate data was collected at 1,000 Hz. The analysis of the acquired experimental data was conducted using Opensim. In addition to the motion capture experiment, surface electromyography (EMG) data were also collected. As described in the previous section, six muscle signals were chosen for EMG recording, specifically the medial femoral muscle, rectus femoris, biceps femoris long head, gastrocnemius Lateralis, soleus muscle, and tibialis anterior muscle. The placement of electrodes follows the guidelines established in Stegeman and Hermens (2007). An eight-channel EMG acquisition box (Biovison, 8 Chanel, Germany), was utilized to capture EMG signals, with two of its channels synchronized with the VICON system at a sampling rate of 1,000 Hz.

To enhance the resemblance of the exoskeleton-driven human's motion to that of healthy individual, an experiment was conducted. The experiment involves a healthy adult male subject, with height of 1.78 m and weight of 56 kg, and without lower limb injury or illness. The subject was instructed to assume a double-stance position, with one foot placed in front of the other, each foot in contact with a force plate. The subject's arms hanged naturally at their sides. Upon receiving an “all set” signal, the participant initiated the gait and proceeded to walk straight ahead to the end of the walkway. The speed, height, and step length were determined under the most comfortable conditions.

The subject was fully informed about the experiment, and the experiment was being conducted following the review by the Institutional Ethics Review Board at Northwestern Polytechnical University (No. 202302009).

2.2.2 Motion biomechanics analysis

After obtaining experimental data, the next step is to analyze the motion data to extract biomechanical features. The subject's joint angles are obtained through kinematic simulation using motion capture data from markers and plantar force data. Inverse dynamics simulation is then performed to calculate joint torques, and based on the plantar force data, the ZMP trajectory can be derived.

Figure 5A illustrates some key frames during the motion, such as toe off, maximum altitude of heel, maximum altitude of toe, and heel on. At time t1, the toe-off event occurs, and the motion before t1 is in the double support phase (DSP). At time t2, the heel makes contact with the ground, and this phase is also considered as the double support phase (DSP). At the conclusion of the single support phase (SSP), the supporting foot pivots around its toe, while the swinging foot rotates around its heel upon contact with the ground. It's important to note that the heel of the swinging foot touches the ground before the heel of the supporting foot lifts off. The phase between t1 and t2 is the single support phase (SSP).

Figure 5

Figure 5. Motion analysis for the walking experiment. (A) Some key frames during the motion. (B) Snapshot of the motion. (C) ZMP trajectory of the three experiments.

The snapshot of the motion is presented in Figure 5B, with green indicating the stance leg (left leg in this case) and black representing the swing leg (right leg in this case). There are four trajectories shown in the figure, with the ankle joint of the stance leg not displayed as its range of motion is limited. Evidently during walking, the hip joint undergoes a subtle, periodic up-and-down motion in the vertical direction. This periodic motion in the upper body leverages gravity to provide forward momentum for the body while reducing excessive energy consumption caused by knee joint bending. This results in an overall decrease in energy expenditure during walking (Zhu et al., 2016).

The ZMP trajectory is illustrated in Figure 5C, a vital indicator of human balance. During the initiation of forward walking, the ZMP initially moves backward, and the center of gravity shifts rearward to store energy. Subsequently, the torso swings forward, and the center of gravity shifts forward, initiating forward movement. During the single support phase, the Zero Moment Point (ZMP) safely progresses along the direction of the foot. In the double support phase, the ZMP moves in a straight line from the rear foot to the front foot. In the process of generating humanoid gait, it's essential to consider key features of human motion, such as the heel strike and toe-off of the feet, the movement of the center of gravity, and the trajectory of the ZMP.

2.2.3 Human-like trajectory generation based on biomechanics features

This section is dedicated to the generation of motion trajectory for the human model driven by expskeleton. The process of generating humanoid trajectories is based on the methodology outlined in Boutin et al. (2011).

Initially, the joint angles obtained from OpenSim are directly mapped to the human model. However, this process can generate incorrect positions due to kinematic mismatches between the experimental subject and the elderly human model. Consequently, it becomes necessary to adjust the incorrect foot positions. The principles guiding this adjustment process are as follows: ensure that the human model's feet do not intersect with the ground and that the model does not slide during motion.

After the adjustment, key phases in the sequence are extracted, such as the moments when the heel strikes the ground and when the toe leaves the floor as shown in Figure 5A. The generated foot trajectories closely resemble those of the human body, ensuring that the angles of toe-off and heel strike are similar to those of the subject. The trajectories of the feet at these moments are linearly fitted to maintain the continuity of position, velocity, and acceleration. While these generated trajectories maintain a similar gait frequency to that of subject, the walking step length is adjusted based on the size of the human model.

To ensure model balance, it is crucial to control the ZMP trajectory. As previously analyzed, the ZMP remains inside the supporting leg and advances in the direction of the foot during the single stance phase. In the double stance phase, the ZMP moves in a straight line from the rearfoot to the forefoot. Based on this analysis and the previously generated feet trajectories, the ZMP trajectory of human model is planned. Although it is challenging to derive joint angles directly from ZMP trajectories, it is possible to calculate joint angles from the Center of Mass (CoM) using inverse kinematics simulation. Human motion analysis indicates that the CoM moves up and down during movement, and this information can be used in human-like motion generation for exoskeleton. An inverted pendulum model can be employed to find a CoM trajectory that ensures the ZMP follows the reference trajectory (Fayong et al., 2014).

After obtaining these values, optimization is performed using a damped least squares method to adjust the joint angles while taking into account the constraints related to the Center of Mass (CoM) and foot motion. This optimization process allows us to derive the desired trajectory for the human model to be driven by the cable-driven exoskeleton.

2.2.4 Neuromuscular control analysis

In the realm of cable-driven exoskeleton robot system, a key focus lies in understanding and validating the neuromuscular control of human muscle forces during different activities. The human lower limb's muscular and skeletal system exhibits a high degree of redundancy, meaning that various muscle combinations and joint movement patterns can accomplish similar tasks. In such scenario, the use of Neuro-Cybernetics control methods can more effectively simulate and comprehend how the human biological nervous system coordinates muscles to perform motor tasks (Wong, 2023). To validate the biomechanics of these systems and understand muscle activation during human motion, the CEINMS toolbox is employed, as detailed in Pizzolato et al. (2015).

To achieve this, electromyogram (EMG) signals collected from the subject play a pivotal role. These signals are processed meticulously, following a structured approach. Initially, a zero-lag fourth-order recursive Butterworth high-pass filter is applied, efficiently eliminating noise from the raw EMG data. Subsequent steps involve full-wave rectification and low-pass filtering using a Butterworth filter with 6 Hz cutoff frequency, adhering to established techniques (Lloyd and Besier, 2003). These processed EMG signals are then mapped onto Musculotendon Units (MTUs), effectively associating specific muscle groups with their corresponding activation levels. This mapping is a vital bridge connecting EMG data to the intricacies of muscle dynamics. Further insights into the biomechanics are gleaned by calculating joint angles, Musculotendon Unit (MTU) lengths, and moment arms. Inverse kinematics aids in determining joint angles, while inverse dynamics calculations provide the necessary MTU parameters. OpenSim's muscle analysis tools prove instrumental in these computations. Crucially, the human body exhibits variations in muscle parameters across individuals, necessitating a meticulous calibration process. This step aims to optimize parameters within the Neuromusculoskeletal (NMS) model. Employing a calibration loop, the process minimizes discrepancies between estimated joint moments, derived from the NMS model, and experimental joint moments, directly measured from the subject. To initiate this calibration, initial parameters for the MTUs are drawn from a standardized model, such as the Gait2354 model within OpenSim. With successfully calibrated parameters, the estimation of Musculotendon Unit (MTU) forces becomes possible. This estimation considers calibrated MTUs, muscle excitations (activation levels), MTU kinematics (lengths and moment arms), and external joint moments. The culmination of these intricate steps provides valuable insights into muscle activation and forces during various human activities.

These insights, founded on rigorous calibration and meticulous analysis of EMG signals, serve as pivotal neuromuscular control validation data for cable-driven exoskeleton robot systems, enriching our understanding of human locomotion and facilitating the development of advanced robotic systems.

2.3 Optimization method of cable attachment points of exoskeleton

In the previous section, the cable attachment points of the cable-driven system were directly mapped from the human muscle attachment points data and need to be adjusted to better suit our requirements. To optimize these attachment points for improved cable performance that closely simulates the muscles of a healthy individual, a classic intelligent optimization algorithm, the genetic algorithm, was employed to determine the most suitable cable attachment points under conditions resembling human motion.

Within the genetic algorithm, the design variables pertain to the coordinates of the cable attachment points. As the muscle origin and insertion points vary among individuals, the cable attachment points will also vary within a certain interval (Shan, 2003). The upper and lower limits of these design variables are determined by appropriately expanding this interval. To enhance performance, it is crucial to define an appropriate evaluation function. For this purpose, we chose energy efficiency, which is a critical aspect of human dynamic walking and a significant indicator for exoskeleton systems (Roberts et al., 2017).

In this study, to calculate the energy consumption of Cable-driven exoskeleton System, we assume that cable tension is achieved through the use of MYO-muscles. The MYO-muscle is driven by a Maxon brushless motor with a power rating of 100 watts and is equipped with a Hall sensor. However, in the real system implementation, if different type of cable-driven artifical muscle is selected, the corresponding parameters in Equation (16) should be adjusted accordingly. The total energy consumption (E) is calculated as follows (Zhu et al., 2016):

$\begin{array}{l}E={\displaystyle \int P_{\text{motor }}}dt={\displaystyle \int u^T}idt\cong {\displaystyle \int {\displaystyle \sum _{i\in \{L,R\}}(C_m{\tau }_i^T{\dot{q}}_i+C_t{\tau }_i^T{\tau }_i+}}\\ C_v{\dot{q}}_i^T{\dot{q}}_i)C_m=\frac{(2RF+BK)}{{(\gamma K)}^2},C_t=\frac{R}{{(n\gamma K)}^2},C_v=\frac{n^2(FR+BK)F}{K^2}\end{array}$

R, B, K, F, γ, and n represent the armature resistance, back EMF constant, torque constant, viscous friction coefficient, gear efficiency, and gear ratio of the motor of the joint, respectively. These values can be obtained from the Maxon official website.

Taking the energy consumption as the optimization target, the Genetic Algorithm (GA) is selected to optimize the cable attachment points of exoskeleton system. Finally, the parameters for the Genetic Algorithm are configured as follows: Population size: 1,000, Number of terminating evolutionary generations: 100, Selection Probability: 0.7, Probability of crossover: 0.8, Probability of mutation: 0.08. These parameters define how the genetic algorithm operates during the optimization process to find the most suitable cable attachment points for the cable-driven exoskeleton system generating human-like motion.

3 Results and analysis

The exoskeleton dynamic simulation process including the elderly human lower limb modeling (Section 2.1.2), cable driven system modeling (Sections 2.1.3 Section 2.1.4), cable force calculations (Section 2.1.5), human-like walking trajectories generation (Section 2.2.3), and cable attachment points optimization (Section 2.3), were performed using MATLAB. The genetic algorithm utilized parallel computing to accelerate the calculations and took ~42,101 s to complete. The calculation of the subject's muscle forces was conducted with the CEINMS toolbox and took ~10 h.

In Figure 6, the red lines represent the subject's COM trajectories and joint trajectories in the walking experiment, while the black lines depict the motion of the elderly human model driven by the cable exoskeleton, closely resembling human-like movements. From the dynamic simulation results, it's evident that the generated COM and joint trajectories of human model closely matches that of the experimental subject. The motion trends of the human model's torso closely follows the subject's trajectory but with smoother pattern. Notably, the entire trajectory of the human model's left leg aligns more closely with the subject's trajectory compared to the right leg.

Figure 6

Figure 6. COM and joint trajectories of elderly human model driven by cable exoskeleton compared with experimental subject's data.

As shown in Figure 7, the red lines represent the subject's joint torques calculated by Opensim, while the black lines represent the joint angular torques of the human model driven by the exoskeleton during the same motion. All torques have been normalized based on body weight. The figure illustrates that the trends in joint torques for the human model and the subject are generally similar. Specifically, the joint torque trends at the right knee, left hip, left knee, and left ankle are very closely aligned. Prior to t < 0.75, the joint torques of the right leg are notably higher than those of the subject, whereas the joint torques of the left leg are lower. Subsequently, after t > 0.75s, the joint torque trends for both the left and right legs of the human model closely resemble the subject's joint torque changes.

Figure 7

Figure 7. Joint torque of elderly human model driven by cable exoskeleton compared with experimental subject.

The average fitness of the genetic algorithm is presented in Figure 8A. It is noticeable that the average fitness converges at around the 60th generation. Figure 8B showcases the initial cable attachment point distribution, while Figure 8C presents the attachment point distribution after optimization. The blue outlines in Figure 8B represent various parts of the human model, including the trunk, thighs, lower legs, and feet, while the red outlines indicate that cable 1 to cable 9 are configured according to Figure 2B. The content in Figures 8B, C is identical, except for the specific cable attachment points' positions. In these figures, the thick black line corresponds to the left leg of elderly human model, the thick green line represents the right leg, the thin red lines symbolize the cables attached to the left leg, and the thin blue lines denote the cables attached to the right leg. It is evident that a significant difference exists between the initial attachment point distribution and the distribution after optimization. After optimization, the cable attachment points of the hip and knee joint extend to some degree in the direction along the lower limb, with minimal changes in the radial direction of the lower limb. Notably, the attachment points of the cable corresponding to the calf muscle after optimization are positioned closer to the ankle joint, whereas the attachment point of the cable corresponding to TA muscle moves farther away from the ankle joint.

Figure 8

Figure 8. Cable attachment points optimization results. (A) Average fitness of the genetic algorithm. (B) Initial attachment points of cable-driven system. (C) Optimized attachment points of cable-driven system.

In Figure 9A, the cable forces imitating specific muscle of right leg in various conditions are depicted. In this figure, the black line represents cable forces under the initial cable configuration, the blue line represents cable forces after optimization, and the red line represents the subject's muscle forces calculated using the CEIMS toolbox based on experimental data. All forces have been normalized based on body weight. It is apparent that the overall trends in forces of the cable corresponding to specific muscle before and after optimization closely resemble the subject's muscle force trends. Specifically, the cable force trends after optimization, particularly for the VAS muscle, align more closely with the subject's muscle activation level and timing compared to the initial cable force trends.

Figure 9

Figure 9. Corresponding cable forces before and after optimization, compared with experimental subject's muscle forces. (A) Corresponding cable forces of right leg compared with subject's muscle forces. (B) Corresponding cable forces of the left leg compared with subject's muscle forces.

At the hip joint after optimization, the forces of cable simulating single-joint muscles ILPSO and GMAX exhibit slightly larger amplitude changes compared to the cable forces under the initial configuration. Conversely, the forces of cable corresponding to RF muscle displays smaller amplitude changes after optimization compared to the initial configuration. There is not much change in the amplitude of cable force corresponding to the HAMS muscle before and after optimization. After optimization, the amplitude change of cable force imitating the VAS muscle exatation is smaller than that under the original configuration, while the force of cable simulating the BFS muscle activation exhibits greater amplitude of change. The forces of cable corresponding to GAS and SOL muscles after optimization closely resemble the initial muscle forces in the first half, with significantly smaller amplitude of change afterward compared to the forces under the original configuration. The force of cable imitating TA muscle after optimization is greater than the force under the original configuration.

The cable forces of the elderly human model's left leg is presented in Figure 9B. The overall similarity in the trend of cable force changes before and after optimization in the left leg with respect to the subject is not as pronounced as in the right leg. Specifically, the trends in forces of cables corresponding to ILPSO, RF, HAMS, GAS, TA, and SOL muscles closely resemble the subject's muscle activation level and timing. However, for the GMAX, VAS, and BFS muscles, there is a noticeable phase difference in the cable force trends compared to the subject.

In terms of force changes, except for the cable corresponding to HAMS, BFS, and SOL muscles, the amplitudes of cable forces before and after optimization are smaller than those of the subject's muscle forces. Especially for ILPSO, RF, HAMS, GAS and SOL muscles, the trends and magnitudes of cable force changes after optimization are relatively close to muscles activation level and timing.

4 Discussion

While achieving a state in older adults that is similar to that of healthy adults individuals may not be realistic, some biomechanical data from individuals still holds certain reference value. In the design of human model, its height based on average elderly human data is shorter compared to the subject's height. Consequently, the forward speed, displacement, and the amplitude of the center of gravity's vertical oscillation are relatively smaller. In the planned motion, the left leg serves as the stance leg, while the right leg transitions from the stance leg to the swing leg with a larger range of motion. Therefore, the movement trends of the human model's left leg is closer to that of the subject. The joint angles of both ankles are directly mapped and remain unadjusted, hence the ankle joint angle trajectories are more akin to those of the subject.

In the Double Support Phase (DSP) of one walking gait, the Ground Reaction Forces (GRF) acts on both legs. In this paper, GRF is estimated based on the Zero Moment Point (ZMP) of subject, which is not entirely consistent with the human model, leading to not totally motion consistencies between the human model and the subject. However, when the human model's right leg leaves the ground, transitioning from DSP to Single Support Phase (SSP), the left leg remains in a stance position while the right leg swings. During this phase, GRF is entirely applied to the left leg. Consequently, the motion similarity between the human model and the subject is significantly higher.

Because the left leg remains the stance leg throughout and bears a greater load, the overall muscle force on the left leg tends to be larger. The calculation process for the subject's muscle forces in the right leg involves EMG signals from the muscles of the right leg, whereas the subject's muscle force in the left leg is entirely determined through biomechanical calculations. It is evident that the trends in cable force changes for the right leg aligns more closely with the human experiment data. Therefore, in subsequent work, increasing the number of surface electromyography (EMG) sensors can be considered to obtain more accurate muscle force data and better optimized cable forces.

When comparing the human model to the subject, it is observed that the human model has a slightly lower height but greater mass. When both entities move at similar speeds, the overall force experienced by the human model is slightly higher compared to that experienced by the subject. Notably, this difference is more pronounced in left leg. Furthermore, OpenSim modeling includes a rich set of muscles at the hip and knee joints, resulting in muscle force trends that do not entirely match. In contrast, at the ankle joint, where there are fewer muscles involved, the results tend to be more similar. It's worth noting that due to the optimization objective of minimizing energy consumption, cables corresponding to muscles like RF and HAMS, which cross both the hip and knee joints and have a large range of motion and force generation, experience significant reductions in forces after optimization. This also makes the output forces of the actuators be easily realized.

In future work, additional optimization objectives such as joint stiffness can be introduced into the optimization process for cable attachment points. Furthermore, more data from healthy elderly individuals will be considered to generate human-like motion trajectories. An analysis of the characteristics of muscle force calculations during human motion will be conducted, with the aim of enabling single-step calculations to provide more biomechanical support for the cable forces. Additionally, there are also plans to construct a lower limb exoskeleton prototype.

Despite some limitations in the method, the final results still indicate that the subject-specific method for determining the cable attachment points of a exoskeleton contributes to the generation of human-like motion with biomechanical and neuromuscular-control features.

5 Conclusion

This study aimed to investigate a personalized attachment point optimization method for cable-driven exoskeleton. The results indicate that, the elderly human model driven by the cable system exhibits cable forces similar to the muscle activation of healthy individual during similar movements. This similarity can facilitate muscle coordination and, ultimately, neurological recovery. Based on the analysis and results from the previous chapters, the following conclusions can be drawn.

The personalized attachment point optimization method successfully emulated the motion trajectories of experimental subjects, resulting in the exoskeleton device closely aligning with natural human coordination during motion. Estimation of joint torques in the exoskeleton also benefited from the application of personalized attachment points. Torque trends throughout motion closely matched those of experimental subject, particularly during transitions between double-leg support and single-leg support phases. Performance of the exoskeleton in terms of cable output forces and coordination improved through the optimization of personalized cable attachment points. Trends in the forces of cables corresponding to the key muscle groups closely aligned with those of experimental subject, especially in the optimized design where cable force changes better adhered to human motion neuromuscular-control features, such as muscle output forces and muscle activation timing sequence.

While this method still faces certain limitations, such as computational complexity and simplification of biomechanical models, these findings underscore the importance of personalized methods in the design of cable-driven exoskeletons. This approach holds promising prospects for applications in assisting older adults rehabilitation, providing valuable guidance for the cable-driven exoskeleton prototype realization.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by Institutional Ethics Review Board, Northwestern Polytechnical University.

Author contributions

YC: Conceptualization, Data curation, Formal analysis, Methodology, Software, Visualization, Writing—original draft. WY: Conceptualization, Funding acquisition, Methodology, Project administration, Supervision, Writing—review & editing. AB: Conceptualization, Supervision, Writing—review & editing. DL: Data curation, Formal analysis, Investigation, Software, Writing—original draft. SK: Formal Analysis, Writing—original draft. RW: Conceptualization, Funding acquisition, Supervision, Writing—review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was funded by High-end Foreign Experts Introduction Plan (Grant No. G2023183021L) and National Natural Science Foundation of China (Grant No. 61603302).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2024.1327397/full#supplementary-material

References

Alnajjar, F., Khalid, S., Vogan, A. A., Shimoda, S., Nouchi, R., and Kawashima, R. (2019). Emerging cognitive intervention technologies to meet the needs of an aging population: a systematic review. Front. Aging Neurosci. 11:291. doi: 10.3389/fnagi.2019.00291

Crossref Full Text | Google Scholar

Arnold, E. M., Hamner, S. R., Seth, A., Millard, M., and Delp, S. L. (2013). How muscle fiber lengths and velocities affect muscle force generation as humans walk and run at different speeds. J. Exp. Biol. 216, 2150–2160. doi: 10.1242/jeb.075697

PubMed Abstract | Crossref Full Text | Google Scholar

Asbeck, A. T., De Rossi, S. M., Galiana, I., Ding, Y., and Walsh, C. J. (2014). Stronger, smarter, softer: next-generation wearable robots. IEEE Roboti. Automat. Mag. 21, 22–33. doi: 10.1109/MRA.2014.2360283

Crossref Full Text | Google Scholar

Bardi, E., Ambrosini, E., Pirelli, A., Pedrocchi, A., Braghin, F., Covarrubias, M., et al. (2022). “Upper limb exosuit cable routing optimization,” in 2022 International Conference on Rehabilitation Robotics (ICORR) (Rotterdam: IEEE), 1–6.

PubMed Abstract | Google Scholar

Biao, L., Youwei, L., Xiaoming, X., Haoyi, W., and Longhan, X. (2023). Design and control of a flexible exoskeleton to generate a natural full gait for lower-limb rehabilitation. J. Mech. Robot. 15, 011005. doi: 10.1115/1.4054248

Crossref Full Text | Google Scholar

Boutin, L., Eon, A., Zeghloul, S., and Lacouture, P. (2011). From human motion capture to humanoid locomotion imitation application to the robots hrp-2 and hoap-3. Robotica 29, 325–334. doi: 10.1017/S0263574710000172

Crossref Full Text | Google Scholar

Carnevale, L., Calabrò, R. S., Celesti, A., Leo, A., Fazio, M., Bramanti, P., et al. (2018). Toward improving robotic-assisted gait training: can big data analysis help us? IEEE Int. Things J. 6, 1419–1426. doi: 10.1109/JIOT.2018.2855937

Crossref Full Text | Google Scholar

Chen, L., Chen, C., Ye, X., Wang, Z., Liu, Y., Cao, W., et al. (2022). A portable waist-loaded soft exosuit for hip flexion assistance with running. Micromachines 13, 157. doi: 10.3390/mi13020157

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, X., Yu, Z., Zhang, W., Zheng, Y., Huang, Q., and Ming, A. (2017). Bioinspired control of walking with toe-off, heel-strike, and disturbance rejection for a biped robot. IEEE Transact. Ind. Elect. 64, 7962–7971. doi: 10.1109/TIE.2017.2698361

Crossref Full Text | Google Scholar

Cheng, C.-K., Chen, H.-H., Chen, C.-S., Lee, C.-L., and Chen, C.-Y. (2000). Segment inertial properties of chinese adults determined from magnetic resonance imaging. Clin. Biomech. 15, 559–566. doi: 10.1016/S0268-0033(00)00016-4

PubMed Abstract | Crossref Full Text | Google Scholar

Choudhury, A., Li, H., Greene, C., and Perumalla, S. (2018). Humanoid robot-application and influence. arXiv [preprint]. doi: 10.26502/acbr.50170059

Crossref Full Text | Google Scholar

De Rossi, S. M. M., O'donnell, K. E., Bae, J., Asbeck, A. T., Holt, K. G., and Walsh, C. J. (2022). Assistive Flexible Suits, Flexible Suit sysTems, and Methods for Making and Control Thereof to Assist Human Mobility. US Patent 11,324,655. Washington, DC: U. S. Patent and Trademark Office.

Google Scholar

Di Giulio, I., Maganaris, C. N., Baltzopoulos, V., and Loram, I. D. (2009). The proprioceptive and agonist roles of gastrocnemius, soleus and tibialis anterior muscles in maintaining human upright posture. J. Physiol. 587, 2399–2416. doi: 10.1113/jphysiol.2009.168690

PubMed Abstract | Crossref Full Text | Google Scholar

Eden, J., Song, C., Tan, Y., Oetomo, D., and Lau, D. (2018). “Caspr-ros: a generalised cable robot software in ros for hardware,” in Cable-Driven Parallel Robots: Proceedings of the Third International Conference on Cable-Driven Parallel Robots (Cham: Springer), 50–61.

Google Scholar

Farhat, N., Zamora, P., Reichert, D., Mata, V., Page, A., and Valera, A. (2022). A computationally efficient musculoskeletal model of the lower limb for the control of rehabilitation robots: assumptions and validation. Appl. Sci. 12, 2654. doi: 10.3390/app12052654

Crossref Full Text | Google Scholar

Fayong, G., Jianghai, Z., Tao, M., Minzhou, L., and Xiaobo, S. (2014). “A modified gait generator for humanoid robots based on height compensation of center of mass,” in 2014 IEEE International Conference on Robotics and Biomimetics (ROBIO 2014) (Bali: IEEE), 1278–1283.

Google Scholar

Feingold Polak, R., Elishay, A., Shachar, Y., Stein, M., Edan, Y., and Levy Tzedek, S. (2018). “Differences between young and old users when interacting with a humanoid robot: a qualitative usability study,” in Companion of the 2018 ACM/IEEE International Conference on Human-Robot Interaction (Chicago, IL), 107–108.

Google Scholar

Firouzi, V., Davoodi, A., Bahrami, F., and Sharbafi, M. A. (2021). From a biological template model to gait assistance with an exosuit. Bioinspir. Biomimet. 16, 066024. doi: 10.1088/1748-3190/ac2e0d

PubMed Abstract | Crossref Full Text | Google Scholar

Frizziero, A., Salamanna, F., Della Bella, E., Vittadini, F., Gasparre, G., Nicoli Aldini, N., et al. (2016). The role of detraining in tendon mechanobiology. Front. Aging Neurosci. 8, 43. doi: 10.3389/fnagi.2016.00043

PubMed Abstract | Crossref Full Text | Google Scholar

Ganesan, K., and Gupta, A. (2021). “Human-centric optimal design of biomimetic exosuit for loaded walking: a simulation study,” in Mechanism and Machine Science: Select Proceedings of Asian MMS 2018 (Singapore: Springer), 497–510.

Google Scholar

Glowinski, S., Krzyzynski, T., Bryndal, A., and Maciejewski, I. (2020). A kinematic model of a humanoid lower limb exoskeleton with hydraulic actuators. Sensors 20, 6116. doi: 10.3390/s20216116

PubMed Abstract | Crossref Full Text | Google Scholar

Grosu, S., De Rijcke, L., Grosu, V., Geeroms, J., Vanderboght, B., Lefeber, D., et al. (2018). Driving robotic exoskeletons using cable-based transmissions: a qualitative analysis and overview. Appl. Mech. Rev. 70:60801. doi: 10.1115/1.4042399

Crossref Full Text | Google Scholar

Guo, X., He, B., Lau, K. Y., Chan, P. P., Liu, R., Xie, J. J., et al. (2022). Age-related modifications of muscle synergies and their temporal activations for overground walking. IEEE Transact. Neural Syst. Rehabil. Eng. 30, 2700–2709. doi: 10.1109/TNSRE.2022.3206887

PubMed Abstract | Crossref Full Text | Google Scholar

He, L., Xu, C., and Guan, X. (2023). Design methodology and experimental study of a lower extremity soft exosuit. Electronics 12, 2502. doi: 10.3390/electronics12112502

Crossref Full Text | Google Scholar

Hong, Y. N. G., Ballekere, A. N., Fregly, B. J., and Roh, J. (2021). Are muscle synergies useful for stroke rehabilitation? Curr. Opin. Biomed. Eng. 19:100315. doi: 10.1016/j.cobme.2021.100315

Crossref Full Text | Google Scholar

Hu, H., Li, Z., Yan, J., Wang, X., Xiao, H., Duan, J., et al. (2007). Anthropometric measurement of the chinese elderly living in the beijing area. Int. J. Ind. Ergon. 37, 303–311. doi: 10.1016/j.ergon.2006.11.006

Crossref Full Text | Google Scholar

Hwang, B., Oh, B.-M., and Jeon, D. (2018). An optimal method of training the specific lower limb muscle group using an exoskeletal robot. IEEE Transact. Neural Syst. Rehabil. Eng. 26, 830–838. doi: 10.1109/TNSRE.2018.2816961

PubMed Abstract | Crossref Full Text | Google Scholar

Joshi, S., Beck, I., Seth, A., and Della Santina, C. (2022). “Minimalistic soft exosuit for assisting the shoulder via biomechanics-aware optimization,” in 2022 IEEE-RAS 21st International Conference on Humanoid Robots (Humanoids) (Ginowan: IEEE), 667–673.

Google Scholar

Junius, K., Moltedo, M., Cherelle, P., Rodriguez-Guerrero, C., Vanderborght, B., and Lefeber, D. (2017). Biarticular elements as a contributor to energy efficiency: biomechanical review and application in bio-inspired robotics. Bioinspir. Biomimet. 12, 061001. doi: 10.1088/1748-3190/aa806e

PubMed Abstract | Crossref Full Text | Google Scholar

Lau, D., Eden, J., Tan, Y., and Oetomo, D. (2016). “Caspr: a comprehensive cable-robot analysis and simulation platform for the research of cable-driven parallel robots,” in 2016 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS) (Daejeon: IEEE), 3004–3011.

Google Scholar

Lee, H. D., Park, H., Seongho, B., and Kang, T. H. (2020). Development of a soft exosuit system for walking assistance during stair ascent and descent. Int. J. Control Automat. Syst. 18, 2678–2686. doi: 10.1007/s12555-019-0584-5

Crossref Full Text | Google Scholar

Li, Y., Fan, T., Qi, Q., Wang, J., Qiu, H., Zhang, L., et al. (2021). Efficacy of a novel exoskeletal robot for locomotor rehabilitation in stroke patients: a multi-center, non-inferiority, randomized controlled trial. Front. Aging Neurosci. 13:706569. doi: 10.3389/fnagi.2021.706569

PubMed Abstract | Crossref Full Text | Google Scholar

Li, Z., Liu, H., Yin, Z., and Chen, K. (2019). Muscle synergy alteration of human during walking with lower limb exoskeleton. Front. Neurosci. 12:1050. doi: 10.3389/fnins.2018.01050

PubMed Abstract | Crossref Full Text | Google Scholar

Liang, J., Wang, H., Zeng, Y., Qu, Y., Liu, Q., Zhao, F., et al. (2021). Physical exercise promotes brain remodeling by regulating epigenetics, neuroplasticity and neurotrophins. Rev. Neurosci. 32, 615–629. doi: 10.1515/revneuro-2020-0099

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, J., and Lockhart, T. E. (2009). Age-related joint moment characteristics during normal gait and successful reactive-recovery from unexpected slip perturbations. Gait Post. 30, 276–281. doi: 10.1016/j.gaitpost.2009.04.005

PubMed Abstract | Crossref Full Text | Google Scholar

Lloyd, D. G., and Besier, T. F. (2003). An emg-driven musculoskeletal model to estimate muscle forces and knee joint moments in vivo. J. Biomech. 36, 765–776. doi: 10.1016/S0021-9290(03)00010-1

Crossref Full Text | Google Scholar

Lo Coco, D., Lopez, G., and Corrao, S. (2016). Cognitive impairment and stroke in elderly patients. Vasc. Health Risk Manag. 105–116. doi: 10.2147/VHRM.S75306

PubMed Abstract | Crossref Full Text | Google Scholar

Marques, H. G., Christophe, M., Lenz, A., Dalamagkidis, K., Culha, U., Siee, M., et al. (2013). “Myorobotics: a modular toolkit for legged locomotion research using musculoskeletal designs,” in 6th International Symposium on Adaptive Motion of Animals and Machines, Darmstadt, Germany, Vol. 11 (Darmstadt).

Google Scholar

Nakashima, I., Imoto, D., Hirano, S., Mukaino, M., Imaida, M., Saitoh, E., et al. (2020). “Development of an abnormal gait analysis system in gait exercise assist robot “welwalk”? for hemiplegic stroke patients,” in 2020 8th IEEE RAS/EMBS International Conference for Biomedical Robotics and Biomechatronics (BioRob) (New York, NY: IEEE), 1030–1035.

Google Scholar

Park, J. (2008). Synthesis of natural arm swing motion in human bipedal walking. J. Biomech. 41, 1417–1426. doi: 10.1016/j.jbiomech.2008.02.031

PubMed Abstract | Crossref Full Text | Google Scholar

Pin-Barre, C., and Laurin, J. (2015). Physical exercise as a diagnostic, rehabilitation, and preventive tool: influence on neuroplasticity and motor recovery after stroke. Neural Plast. 2015:608581. doi: 10.1155/2015/608581

PubMed Abstract | Crossref Full Text | Google Scholar

Pino, M., Boulay, M., Jouen, F., and Rigaud, A.-S. (2015). “Are we ready for robots that care for us??? attitudes and opinions of older adults toward socially assistive robots. Front. Aging Neurosci. 7:141. doi: 10.3389/fnagi.2015.00141

PubMed Abstract | Crossref Full Text | Google Scholar

Pizzolato, C., Lloyd, D. G., Sartori, M., Ceseracciu, E., Besier, T. F., Fregly, B. J., et al. (2015). Ceinms: a toolbox to investigate the influence of different neural control solutions on the prediction of muscle excitation and joint moments during dynamic motor tasks. J. Biomech. 48, 3929–3936. doi: 10.1016/j.jbiomech.2015.09.021

PubMed Abstract | Crossref Full Text | Google Scholar

Prasad, R., El-Rich, M., Awad, M. I., Agrawal, S. K., and Khalaf, K. (2023). Bi-planar trajectory tracking with a novel 3dof cable driven lower limb rehabilitation exoskeleton (c-lrex). Sensors 23, 1677. doi: 10.3390/s23031677

PubMed Abstract | Crossref Full Text | Google Scholar

Prasad, R., El-Rich, M., Awad, M. I., Hussain, I., Jelinek, H., Huzaifa, U., et al. (2022). A framework for determining the performance and requirements of cable-driven mobile lower limb rehabilitation exoskeletons. Front. Bioeng. Biotechnol. 10:920462. doi: 10.3389/fbioe.2022.920462

PubMed Abstract | Crossref Full Text | Google Scholar

Roberts, D., Quacinella, J., and Kim, J. H. (2017). Energy expenditure of a biped walking robot: instantaneous and degree-of-freedom-based instrumentation with human gait implications. Robotica 35, 1054–1071. doi: 10.1017/S0263574715000983

Crossref Full Text | Google Scholar

Schmidt, K., Duarte, J. E., Grimmer, M., Sancho-Puchades, A., Wei, H., Easthope, C. S., et al. (2017). The myosuit: bi-articular anti-gravity exosuit that reduces hip extensor activity in sitting transfers. Front. Neurorobot. 11:57. doi: 10.3389/fnbot.2017.00057

PubMed Abstract | Crossref Full Text | Google Scholar

Shan, D.-M. (2003). Study on the muscular function model and its application of human lower extremity. Sports Sci. Res. 24, 16–20.

Google Scholar

Sharbafi, M. A., Barazesh, H., Iranikhah, M., and Seyfarth, A. (2018). Leg force control through biarticular muscles for human walking assistance. Front. Neurorobot. 12:39. doi: 10.3389/fnbot.2018.00039

PubMed Abstract | Crossref Full Text | Google Scholar

Shur, N., Creedon, L., Skirrow, S., Atherton, P., MacDonald, I., Lund, J., et al. (2021). Age-related changes in muscle architecture and metabolism in humans: the likely contribution of physical inactivity to age-related functional decline. Ageing Res. Rev. 68:101344. doi: 10.1016/j.arr.2021.101344

PubMed Abstract | Crossref Full Text | Google Scholar

Singh, R. E., Iqbal, K., White, G., and Hutchinson, T. E. (2018). A systematic review on muscle synergies: from building blocks of motor behavior to a neurorehabilitation tool. Appl. Bionics Biomech. 2018:3615368. doi: 10.1155/2018/3615368

PubMed Abstract | Crossref Full Text | Google Scholar

Sleimen-Malkoun, R., Temprado, J.-J., and Hong, S. L. (2014). Aging induced loss of complexity and dedifferentiation: consequences for coordination dynamics within and between brain, muscular and behavioral levels. Front. Aging Neurosci. 6:140. doi: 10.3389/fnagi.2014.00140

PubMed Abstract | Crossref Full Text | Google Scholar

Stegeman, D., and Hermens, H. (2007). Standards for surface electromyography: the european project surface emg for non-invasive assessment of muscles (seniam). Enschede 10, 8–12.

Google Scholar

Thompson, M. Q., Yu, S., Tucker, G. R., Adams, R. J., Cesari, M., Theou, O., et al. (2021). Frailty and sarcopenia in combination are more predictive of mortality than either condition alone. Maturitas 144, 102–107. doi: 10.1016/j.maturitas.2020.11.009

PubMed Abstract | Crossref Full Text | Google Scholar

Tonazzini, A., Shintake, J., Rognon, C., Ramachandran, V., Mintchev, S., and Floreano, D. (2018). “Variable stiffness strip with strain sensing for wearable robotics,” in 2018 IEEE International Conference on Soft Robotics (RoboSoft) (Livorno: IEEE), 485–490.

Google Scholar

Vernooij, C. A., Rao, G., Berton, E., Retornaz, F., and Temprado, J.-J. (2016). The effect of aging on muscular dynamics underlying movement patterns changes. Front. Aging Neurosci. 8:309. doi: 10.3389/fnagi.2016.00309

PubMed Abstract | Crossref Full Text | Google Scholar

Wong, K. K. (2023). Cybernetical Intelligence: Engineering Cybernetics with Machine Intelligence. Hoboken, NJ: John Wiley & Sons.

Google Scholar

World Health Organization (2020). Decade of Healthy Ageing: Baseline Report. Geneva.

Google Scholar

Wu, Q., Chen, B., and Wu, H. (2019). Neural-network-enhanced torque estimation control of a soft wearable exoskeleton for elbow assistance. Mechatronics 63:102279. doi: 10.1016/j.mechatronics.2019.102279

Crossref Full Text | Google Scholar

Xiang, Y., Arora, J. S., Rahmatalla, S., and Abdel-Malek, K. (2009). Optimization-based dynamic human walking prediction: one step formulation. Int. J. Numer. Methods Eng. 79, 667–695. doi: 10.1002/nme.2575

Crossref Full Text | Google Scholar

Xiloyannis, M., Alicea, R., Georgarakis, A.-M., Haufe, F. L., Wolf, P., Masia, L., et al. (2021). Soft robotic suits: state of the art, core technologies, and open challenges. IEEE Transact. Robot. 38, 1343–1362. doi: 10.1109/TRO.2021.3084466

Crossref Full Text | Google Scholar

Xiong, F., Liao, X., Xiao, J., Bai, X., Huang, J., Zhang, B., et al. (2022). Emerging limb rehabilitation therapy after post-stroke motor recovery. Front. Aging Neurosci. 14:863379. doi: 10.3389/fnagi.2022.863379

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, H. D., Cooper, M., Eckert-Erdheim, A., Orzel, D., and Walsh, C. J. (2022). A soft exosuit assisting hip abduction for knee adduction moment reduction during walking. IEEE Robot. Automat. Lett. 7, 7439–7446. doi: 10.1109/LRA.2022.3182106

Crossref Full Text | Google Scholar

Yeem, S., Heo, J., Kim, H., and Kwon, Y. (2018). Technical analysis of exoskeleton robot. World J. Eng. Technol. 7, 68–79. doi: 10.4236/wjet.2019.71004

Crossref Full Text | Google Scholar

Yoo, S.-Z., No, M.-H., Heo, J.-W., Park, D.-H., Kang, J.-H., Kim, S. H., et al. (2018). Role of exercise in age-related sarcopenia. J. Exerc. Rehabil. 14, 551. doi: 10.12965/jer.1836268.134

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, S., Yang, Y., Gao, Y., Zhang, Z., Zheng, T., and Zhu, Y. (2019). “Development of a soft knee exosuit with twisted string actuators for stair climbing assistance,” in 2019 IEEE international conference on robotics and biomimetics (ROBIO) (Dali: IEEE), 2541–2546.

Google Scholar

Zhu, H., Luo, M., Mei, T., Zhao, J., Li, T., and Guo, F. (2016). Energy-efficient bio-inspired gait planning and control for biped robot based on human locomotion analysis. J. Bionic Eng. 13, 271–282. doi: 10.1016/S1672-6529(16)60300-1

Crossref Full Text | Google Scholar

Ziaaldini, M. M., Marzetti, E., Picca, A., and Murlasits, Z. (2017). Biochemical pathways of sarcopenia and their modulation by physical exercise: a narrative review. Front. Med. 4:167. doi: 10.3389/fmed.2017.00167

Crossref Full Text | Google Scholar