Tapping Into Skeletal Muscle Biomechanics for Design and Control of Lower Limb Exoskeletons: A Narrative Review.

Lower limb exoskeletons and exosuits ("exos") are traditionally designed with a strong focus on mechatronics and actuation, whereas the "human side" is often disregarded or minimally modeled. Muscle biomechanics principles and skeletal muscle response to robot-delivered loads should be incorporated in design/control of exos. In this narrative review, we summarize the advances in literature with respect to the fusion of muscle biomechanics and lower limb exoskeletons. We report methods to measure muscle biomechanics directly and indirectly and summarize the studies that have incorporated muscle measures for improved design and control of intuitive lower limb exos. Finally, we delve into articles that have studied how the human-exo interaction influences muscle biomechanics during locomotion. To support neurorehabilitation and facilitate everyday use of wearable assistive technologies, we believe that future studies should investigate and predict how exoskeleton assistance strategies would structurally remodel skeletal muscle over time. Real-time mapping of the neuromechanical origin and generation of muscle force resulting in joint torques should be combined with musculoskeletal models to address time-varying parameters such as adaptation to exos and fatigue. Development of smarter predictive controllers that steer rather than assist biological components could result in a synchronized human-machine system that optimizes the biological and electromechanical performance of the combined system.

A review of user needs to drive the development of lower limb prostheses.

BACKGROUND: The development of bionic legs has seen substantial improvements in the past years but people with lower-limb amputation still suffer from impairments in mobility (e.g., altered balance and gait control) due to significant limitations of the contemporary prostheses. Approaching the problem from a human-centered perspective by focusing on user-specific needs can allow identifying critical improvements that can increase the quality of life. While there are several reviews of user needs regarding upper limb prostheses, a comprehensive summary of such needs for those affected by lower limb loss does not exist. METHODS: We have conducted a systematic review of the literature to extract important needs of the users of lower-limb prostheses. The review included 56 articles in which a need (desire, wish) was reported explicitly by the recruited people with lower limb amputation (N = 8149). RESULTS: An exhaustive list of user needs was collected and subdivided into functional, psychological, cognitive, ergonomics, and other domain. Where appropriate, we have also briefly discussed the developments in prosthetic devices that are related to or could have an impact on those needs. In summary, the users would like to lead an independent life and reintegrate into society by coming back to work and participating in social and leisure activities. Efficient, versatile, and stable gait, but also support to other activities (e.g., sit to stand), contribute to safety and confidence, while appearance and comfort are important for the body image. However, the relation between specific needs, objective measures of performance, and overall satisfaction and quality of life is still an open question. CONCLUSIONS: Identifying user needs is a critical step for the development of new generation lower limb prostheses that aim to improve the quality of life of their users. However, this is not a simple task, as the needs interact with each other and depend on multiple factors (e.g., mobility level, age, gender), while evolving in time with the use of the device. Hence, novel assessment methods are required that can evaluate the impact of the system from a holistic perspective, capturing objective outcomes but also overall user experience and satisfaction in the relevant environment (daily life).

The Sensor-Based Biomechanical Risk Assessment at the Base of the Need for Revising of Standards for Human Ergonomics.

Due to the epochal changes introduced by "Industry 4.0", it is getting harder to apply the varying approaches for biomechanical risk assessment of manual handling tasks used to prevent work-related musculoskeletal disorders (WMDs) considered within the International Standards for ergonomics. In fact, the innovative human-robot collaboration (HRC) systems are widening the number of work motor tasks that cannot be assessed. On the other hand, new sensor-based tools for biomechanical risk assessment could be used for both quantitative "direct instrumental evaluations" and "rating of standard methods", allowing certain improvements over traditional methods. In this light, this Letter aims at detecting the need for revising the standards for human ergonomics and biomechanical risk assessment by analyzing the WMDs prevalence and incidence; additionally, the strengths and weaknesses of traditional methods listed within the International Standards for manual handling activities and the next challenges needed for their revision are considered. As a representative example, the discussion is referred to the lifting of heavy loads where the revision should include the use of sensor-based tools for biomechanical risk assessment during lifting performed with the use of exoskeletons, by more than one person (team lifting) and when the traditional methods cannot be applied. The wearability of sensing and feedback sensors in addition to human augmentation technologies allows for increasing workers' awareness about possible risks and enhance the effectiveness and safety during the execution of in many manual handling activities.

Voluntary control of wearable robotic exoskeletons by patients with paresis via neuromechanical modeling.

BACKGROUND: Research efforts in neurorehabilitation technologies have been directed towards creating robotic exoskeletons to restore motor function in impaired individuals. However, despite advances in mechatronics and bioelectrical signal processing, current robotic exoskeletons have had only modest clinical impact. A major limitation is the inability to enable exoskeleton voluntary control in neurologically impaired individuals. This hinders the possibility of optimally inducing the activity-driven neuroplastic changes that are required for recovery. METHODS: We have developed a patient-specific computational model of the human musculoskeletal system controlled via neural surrogates, i.e., electromyography-derived neural activations to muscles. The electromyography-driven musculoskeletal model was synthesized into a human-machine interface (HMI) that enabled poststroke and incomplete spinal cord injury patients to voluntarily control multiple joints in a multifunctional robotic exoskeleton in real time. RESULTS: We demonstrated patients' control accuracy across a wide range of lower-extremity motor tasks. Remarkably, an increased level of exoskeleton assistance always resulted in a reduction in both amplitude and variability in muscle activations as well as in the mechanical moments required to perform a motor task. Since small discrepancies in onset time between human limb movement and that of the parallel exoskeleton would potentially increase human neuromuscular effort, these results demonstrate that the developed HMI precisely synchronizes the device actuation with residual voluntary muscle contraction capacity in neurologically impaired patients. CONCLUSIONS: Continuous voluntary control of robotic exoskeletons (i.e. event-free and task-independent) has never been demonstrated before in populations with paretic and spastic-like muscle activity, such as those investigated in this study. Our proposed methodology may open new avenues for harnessing residual neuromuscular function in neurologically impaired individuals via symbiotic wearable robots.

Modeling and simulating the neuromuscular mechanisms regulating ankle and knee joint stiffness during human locomotion.

This work presents an electrophysiologically and dynamically consistent musculoskeletal model to predict stiffness in the human ankle and knee joints as derived from the joints constituent biological tissues (i.e., the spanning musculotendon units). The modeling method we propose uses electromyography (EMG) recordings from 13 muscle groups to drive forward dynamic simulations of the human leg in five healthy subjects during overground walking and running. The EMG-driven musculoskeletal model estimates musculotendon and resulting joint stiffness that is consistent with experimental EMG data as well as with the experimental joint moments. This provides a framework that allows for the first time observing 1) the elastic interplay between the knee and ankle joints, 2) the individual muscle contribution to joint stiffness, and 3) the underlying co-contraction strategies. It provides a theoretical description of how stiffness modulates as a function of muscle activation, fiber contraction, and interacting tendon dynamics. Furthermore, it describes how this differs from currently available stiffness definitions, including quasi-stiffness and short-range stiffness. This work offers a theoretical and computational basis for describing and investigating the neuromuscular mechanisms underlying human locomotion.

Electromyography-driven musculoskeletal models with time-varying fatigue dynamics improve lumbosacral joint moments during lifting.

Muscle fatigue is prevalent across different aspects of daily life. Tracking muscle fatigue is useful to understand muscle overuse and possible risk of injury leading to musculoskeletal disorders. Current fatigue models are not suitable for real-world settings as they are either validated using simulations or non-functional tasks. Moreover, models that capture the changes to muscle activity due to fatigue either assume a linear relationship between muscle activity and muscle force or utilize a simple muscle model. Personalised electromygraphy (EMG)-driven musculoskeletal models (pEMS) offer person-specific approaches to model muscle and joint kinetics during a wide repertoire of daily life tasks. These models utilize EMG, thus capturing central fatigue-dependent changes in multi-muscle bio-electrical activity. However, the peripheral muscle force decay is missing in these models. Thus, we studied the influence of fatigue on a large scale pEMS of the trunk. Eleven healthy participants performed functional asymmetric lifting task. Average peak body-weight normalized lumbosacral moments (BW-LM) were estimated to be 2.55 ± 0.26 Nm/kg by reference inverse dynamics. After complete exhaustion of the lower back, the pEMS overestimated the peak BW-LM by 0.64 ± 0.37 Nm/kg. Then, we developed a time-varying muscle force decay model resulting in a time-varying pEMS (t-pEMS). This reduced the difference between BW-LM estimated by the t-pEMS and reference to 0.49 ± 0.14 Nm/kg. We also showed that five fatiguing contractions are sufficient to calibrate the t-pEMS. Thus, this study presents a person and muscle specific model to track fatigue during functional tasks.

The Simultaneous Model-Based Estimation of Joint, Muscle, and Tendon Stiffness is Highly Sensitive to the Tendon Force-Strain Relationship.

OBJECTIVE: Accurate estimation of stiffness across anatomical levels (i.e., joint, muscle, and tendon) in vivo has long been a challenge in biomechanics. Recent advances in electromyography (EMG)-driven musculoskeletal modeling have allowed the non-invasive estimation of stiffness during dynamic joint rotations. Nevertheless, validation has been limited to the joint level due to a lack of simultaneous in vivo experimental measurements of muscle and tendon stiffness. METHODS: With a focus on the triceps surae, we employed a novel perturbation-based experimental technique informed by dynamometry and ultrasonography to derive reference stiffness at the joint, muscle, and tendon levels simultaneously. Here, we propose a new EMG-driven model-based approach that does not require external joint perturbation, nor ultrasonography, to estimate multi-level stiffness. We present a novel set of closed-form equations that enables the person-specific tuning of musculoskeletal parameters dictating biological stiffness, including passive force-length relationships in modeled muscles and tendons. RESULTS: Calibrated EMG-driven musculoskeletal models estimated the reference data with average normalized root-mean-square error ≈ 20%. Moreover, only when calibrated tendons were approximately four times more compliant than typically modeled, our approach could estimate multi-level reference stiffness. CONCLUSION: EMG-driven musculoskeletal models can be calibrated on a larger set of reference data to provide more realistic values for the biomechanical variables across multiple anatomical levels. Moreover, the tendon models that are typically used in musculoskeletal modeling are too stiff. SIGNIFICANCE: Calibrated musculoskeletal models informed by experimental measurements give access to an augmented range of biomechanical variables that might not be easily measured with sensors alone.

Benchmarking commercially available soft and rigid passive back exoskeletons for an industrial workplace.

Low-back pain is a common occupational hazard for industrial workers. Several studies show the advantages of using rigid and soft back-support passive exoskeletons and exosuits (exos) to reduce the low-back loading and risk of injury. However, benefits of using these exos have been shown to be task-specific. Therefore, in this study, we developed a benchmarking approach to assess exos for an industrial workplace at Hankamp Gears B.V. We assessed two rigid (Laevo Flex, Paexo back) and two soft (Auxivo Liftsuit 1.0, and Darwing Hakobelude) exos for tasks resembling the workplace. We measured the assistive moment provided by each exo and their respective influence on muscle activity as well as the user's perception of comfort and exertion. Ten participants performed four lifting tasks (Static hold, Asymmetric, Squat, and Stoop), while their electromyography and subjective measures were collected. The two rigid exos provided the largest assistance during the Dynamic tasks. Reductions in erector spinae activity were seen to be task-specific, with larger reductions for the two rigid exos. Overall, Laevo Flex offered a good balance between assistive moments, reductions in muscle activity, as well as user comfort and reductions in perceived exertion. Thus, we recommend benchmarking exos for intended use in the industrial workplace. This will hopefully result in a better adoption of the back-support exoskeletons in the workplace and help reduce low-back pain.

A Wearable Real-time Kinematic and Kinetic Measurement Sensor Setup for Human Locomotion.

Current laboratory-based setups (optical marker cameras + force plates) for human motion measurement require participants to stay in a constrained capture region which forbids rich movement types. This study established a fully wearable system, based on commercially available sensors (inertial measurement units + pressure insoles) that can measure both kinematic and kinetic motion data simultaneously and support wireless frame-by-frame streaming. In addition, its capability and accuracy were tested against a conventional laboratory-based setup. An experiment was conducted, with 9 participants wearing the wearable measurement system and performing 13 daily motion activities, from slow walking to fast running, together with vertical jump, squat, lunge and single-leg landing, inside the capture space of the laboratory-based motion capture system. The recorded sensor data were post-processed to obtain joint angles, ground reaction forces (GRFs), and joint torques (via multi-body inverse dynamics). Compared to the laboratory-based system, the established wearable measurement system can measure accurate information of all lower limb joint angles (Pearson's r = 0.929), vertical GRFs (Pearson's r = 0.954), and ankle joint torques (Pearson's r = 0.917). Center of pressure (CoP) in the anterior-posterior direction and knee joint torques were fairly matched (Pearson's r = 0.683 and 0.612, respectively). Calculated hip joint torques and measured medial-lateral CoP did not match with the laboratory-based system (Pearson's r = 0.21 and 0.47, respectively). Furthermore, both raw and processed datasets are openly accessible (https://doi.org/10.5281/zenodo.6457662). Documentation, data processing codes, and guidelines to establish the real-time wearable kinetic measurement system are also shared (https://github.com/HuaweiWang/WearableMeasurementSystem).

Person-specific biophysical modelling of alpha-motoneuron pools driven by in vivo decoded neural synaptic input.

Interfacing with alpha-motoneurons (MNs) is key to understand and control motor impairment and neurorehabilitation technologies. Depending on the neurophysiological condition of each individual, MN pools exhibit distinct neuro-anatomical properties and firing behaviors. Hence, the ability to assess subject-specific characteristics of MN pools is essential for unravelling the neural mechanisms and adaptations underlying motor control, both in healthy and impaired individuals. However, measuring in vivo the properties of complete human MN pools remains an open challenge. Therefore, this work proposes a novel approach based on decoding neural discharges from human MNs in vivo for driving the metaheuristic optimization of biophysically realistic MN models. First, we show that this framework provides subject-specific estimates of MN pool properties from the tibialis anterior muscle on five healthy individuals. Second, we propose a methodology to create complete pools of in silico MNs for each subject. Lastly, we show that neural-data driven complete in silico MN pools reproduce in vivo MN firing characteristics and muscle activation profiles during force-tracking tasks involving isometric ankle dorsi-flexion, at different levels of amplitude. This approach can open new avenues for understanding human neuro-mechanics and, particularly, MN pool dynamics, in a person-specific way. Thereby enabling the development of personalized neurorehabilitation and motor restoring technologies.

Predictive Control of Peak Achilles Tendon Force in a Simulated System of the Human Ankle Joint with a Parallel Artificial Actuator During Hopping.

Latest advances in wearable exoskeletons for the human lower extremity predominantly focus on minimising metabolic cost of walking. However, there currently is no robotic exoskeleton that gains control on the mechanics of biological tissues such as biological muscles or series-elastic tendons. Achieving robotic control of biological tissue mechanics would enable prevention of musculoskeletal injuries or the personalization of rehabilitation treatments following injury with levels of precisions not attained before. In this paper, we introduce a new framework that uses nonlinear model predictive control (NMPC) for the closed-loop control of peak tendon force in a simulated system of the human ankle joint with parallel exoskeletal actuation. We propose a computationally efficient NMPC's inner model consisting of explicit, closed-form equations of muscle-tendon dynamics along with those of the ankle joint with parallel actuation. The proposed formulation is tested and verified on movement data collected during dynamic ankle dorsiflexion/plantarflexion rotations executed on a dynamometer as well as during walking and running on a treadmill. The framework designed using the NMPC controller showed a promising performance in keeping the Achilles tendon force under a predefined threshold. Results indicated that our proposed model was generalizable to different muscles and gaits and suitable for real-time applications due to its low computational time.

Characterization of Motor Unit Firing and Twitch Properties for Decoding Musculoskeletal Force in the Human Ankle Joint In Vivo.

Understanding how motor units (MUs) contribute to skeletal mechanical force is crucial for unraveling the underlying mechanism of human movement. Alterations in MU firing, contractile and force-generating properties emerge in response to physical training, aging or injury. However, how changes in MU firing and twitch properties dictate skeletal muscle force generation in healthy and impaired individuals remains an open question. In this work, we present a MU-specific approach to identify firing and twitch properties of MU samples and employ them to decode musculoskeletal function in vivo. First, MU firing events were decomposed offline from high-density electromyography (HD-EMG) of six lower leg muscles involved in ankle plantar-dorsi flexion. We characterized their twitch responses based on the statistical distributions of their firing properties and employed them to compute MU-specific activation dynamics. Subsequently, we decoded ankle joint moments by linking our framework to a subject-specific musculoskeletal model. We validated our approach at different ankle positions and levels of activation and compared it with traditional EMG-driven models. Our proposed MU-specific formulation achieves higher generalization across conditions than the EMG-driven models, with significantly lower coefficients of variation in torque predictions. Furthermore, our approach shows distinct neural strategies across a large repertoire of contractile conditions in different muscles. Our proposed approach may open new avenues for characterizing the relationship between MU firing and twitch properties and their influence on force capacity. This can facilitate the development of targeted rehabilitation strategies tailored to individuals with specific neuromuscular conditions.

A Neuromechanical Model-Based Strategy to Estimate the Operator's Payload in Industrial Lifting Tasks.

One of the main technological barriers hindering the development of active industrial exoskeleton is today represented by the lack of suitable payload estimation algorithms characterized by high accuracy and low calibration time. The knowledge of the payload enables exoskeletons to dynamically provide the required assistance to the user. This work proposes a payload estimation methodology based on personalized Electromyography-driven musculoskeletal models (pEMS) combined with a payload estimation method we called "delta torque" that allows the decoupling of payload dynamical properties from human dynamical properties. The contribution of this work lies in the conceptualization of such methodology and its validation considering human operators during industrial lifting tasks. With respect to existing solutions often based on machine learning, our methodology requires smaller training datasets and can better generalize across different payloads and tasks. The proposed payload estimation methodology has been validated on lifting tasks with 0kg, 5kg, 10kg and 15kg, resulting in an average MAE of about 1.4 Kg. Even if 5kg and 10Kg lifting tasks were out of the training set, the MAE related to these tasks are 1.6 kg and 1.1 kg, respectively, demonstrating the generalizing property of the proposed methodology. To the best of the authors' knowledge, this is the first time that an EMG-driven model-based approach is proposed for human payload estimation.

Automated spatial localization of ankle muscle sites and model-based estimation of joint torque post-stroke via a wearable sensorised leg garment.

Assessing a patient's musculoskeletal function during over-ground walking is a primary objective in post-stroke rehabilitation, due to the importance of walking recovery for everyday life. However, the quantitative assessment of musculoskeletal function currently requires lab-constrained equipment, and labor-intensive analyses, which hampers assessment in standard clinical settings. The development of fully wearable systems for the online estimation of muscle-tendon forces and resulting joint torque would aid clinical assessment of motor recovery, it would enhance the detection of neuro-muscular anomalies and it would consequently enable highly personalized treatments. Here, we present a wearable technology that combines (1) a soft garment for the human leg sensorized with 64 flexible and dry electromyography (EMG) electrodes, (2) a generalized and automated algorithm for the localization of leg muscle sites, and (3) an EMG-driven musculoskeletal modeling framework for the estimation of ankle dorsi-plantar flexion torques. Our results showed that the automated clustering algorithm could detect muscle locations in both healthy and post-stroke individuals. The estimated muscle-specific EMG envelopes could be used to drive forward person-specific musculoskeletal models and estimate resulting joint torques accurately across all healthy and post-stroke individuals and across different walking speeds (R2  > 0.82 and RMSD  < 0.16). The technology we proposed opens new avenues for automated muscle localization and quantitative musculoskeletal function assessment during gait in both healthy and neurologically impaired individuals.

Real-time myoelectric control of wrist/hand motion in Duchenne muscular dystrophy: A case study.

Introduction: Duchenne muscular dystrophy (DMD) is a genetic disorder that induces progressive muscular degeneration. Currently, the increase in DMD individuals' life expectancy is not being matched by an increase in quality of life. The functioning of the hand and wrist is central for performing daily activities and for providing a higher degree of independence. Active exoskeletons can assist this functioning but require the accurate decoding of the users' motor intention. These methods have, however, never been systematically analyzed in the context of DMD. Methods: This case study evaluated direct control (DC) and pattern recognition (PR), combined with an admittance model. This enabled customization of myoelectric controllers to one DMD individual and to a control population of ten healthy participants during a target-reaching task in 1- and 2- degrees of freedom (DOF). We quantified real-time myocontrol performance using target reaching times and compared the differences between the healthy individuals and the DMD individual. Results and Discussion: Our findings suggest that despite the muscle tissue degeneration, the myocontrol performance of the DMD individual was comparable to that of the healthy individuals in both DOFs and with both control approaches. It was also evident that PR control performed better for the 2-DOF tasks for both DMD and healthy participants, while DC performed better for the 1-DOF tasks. The insights gained from this study can lead to further developments for the intuitive multi-DOF myoelectric control of active hand exoskeletons for individuals with DMD.

Synergy-Driven Musculoskeletal Modeling to Estimate Muscle Excitations and Joint Moments at Different Walking Speeds in Individuals with Transtibial Amputation.

The main requirement for an amputee is to regain the function of the lost limb. In order to fully benefit from powered prosthetic legs, the user must rely on the dynamic control of the device. Progress in high-level control for powered prosthetic legs is currently challenged by the inability of current control schemes to generalize to large repertoires of movements as well as adapting to external mechanical demands. This ultimately leads the user to adopt compensatory movements, lack of comfort, higher energy requirements during walking and standing. This study uses a feedforward model of muscle activation and force generation that applies mathematical formulations of muscle synergies to generate synthetic activation profiles underlying walking across different speeds. Estimated activation profiles are used to drive forward subject-specific numerical models of the lower extremity musculoskeletal system. The model was validated on one individual with uni-lateral transtibial amputation and its predictions were compared to experimental torques from inverse dynamic calculations. Results showed that a generic muscle synergy driven personalized musculoskeletal model can fit the ankle torques of the intact limb of a person with transtibial amputation (RMSD = 0.1329±0.02). The estimated moments might be suitable as the control signal to drive powered prostheses to ultimately improve physical interaction between the user and a powered prostheses during dynamic motor tasks.

Trunk extensor muscle endurance and its relationship to action potential conduction velocity and spectral parameters estimated using high-density electromyography.

Trunk extensor muscle fatigue typically manifests as a decline in spectral content of surface electromyography. However, previous research on the relationship of this decline with trunk extensor muscle endurance have shown inconsistent results. The decline of spectral content mainly reflects the decrease in average motor unit action potential conduction velocity (CV). We evaluated whether the rate of change in CV, as well as two approaches employing the change in spectral content, are related to trunk extensor muscle endurance. Fourteen healthy male participants without a low-back pain history performed a non-strictly controlled static forward trunk bending trial until exhaustion while standing. For 13 participants, physiologically plausible CV estimates were obtained from high-density surface electromyography bilaterally from T6 to L5. Laterally between L1 and L2, the linear rate of CV change was strongly correlated to endurance time (R2 = 0.79), whereas analyses involving the linear rate of change in spectral measures showed a lower (R2 = 0.38) or no correlation. For medial electrode locations, estimating CV and its relationship with endurance time was less successful, while the linear rate of change in spectral measures correlated moderately to endurance time (R2 = 0.44; R2 = 0.56). This study provides guidance on monitoring trunk extensor muscle fatigue development using electromyography.

Does a Soft Actuated Back Exosuit Influence Multimodal Physiological Measurements and User Perception During an Industry Inspired Task?

Back support soft exosuits are promising solutions to reduce risk of musculoskeletal injuries at workplaces resulting from physically demanding and repetitive lifting tasks. Design of novel active exosuits address the impact on the muscle activity and metabolic costs but do not consider other critical aspects such as comfort and user perception during the intended tasks. Thus, in this study, we describe a novel soft active exosuit in line with its impact on physiological and subjective measures during lifting. We tested four healthy participants who performed repetitive lifting tasks with and without this exosuit. The exosuit provided assistance proportional to the lumbar flexion angle measured using an inertial measurement unit. We measured the participant's multimodal physiological measures including surface electromyography, metabolic cost, heart rate, and skin temperature. We also measured subjective scores on user exertion, task load, and device acceptability. All participants perceived a reduction in task load when using the exosuit. Three participants showed reduction of muscle activity for the erector spinae muscles. The metabolic costs and heart rate reserve reduced for two participants, with similar trends for skin temperature. For future development of workplace exosuits, we recommend incorporating assessments of both physiological and subjective measures, considering the user-dependent response to the exosuit.

Automated estimation of ankle muscle EMG envelopes and resulting plantar-dorsi flexion torque from 64 garment-embedded electrodes uniformly distributed around the human leg.

The design of personalized movement training and rehabilitation pipelines relies on the ability of assessing the activation of individual muscles concurrently with the resulting joint torques exerted during functional movements. Despite advances in motion capturing, force sensing and bio-electrical recording technologies, the estimation of muscle activation and resulting force still relies on lengthy experimental and computational procedures that are not clinically viable. This work proposes a wearable technology for the rapid, yet quantitative, assessment of musculoskeletal function. It comprises of (1) a soft leg garment sensorized with 64 uniformly distributed electromyography (EMG) electrodes, (2) an algorithm that automatically groups electrodes into seven muscle-specific clusters, and (3) a EMG-driven musculoskeletal model that estimates the resulting force and torque produced about the ankle joint sagittal plane. Our results show the ability of the proposed technology to automatically select a sub-set of muscle-specific electrodes that enabled accurate estimation of muscle excitations and resulting joint torques across a large range of biomechanically diverse movements, underlying different excitation patterns, in a group of eight healthy individuals. This may substantially decrease time needed for localization of muscle sites and electrode placement procedures, thereby facilitating applicability of EMG-driven modelling pipelines in standard clinical protocols.

Electromyography-driven model-based estimation of ankle torque and stiffness during dynamic joint rotations in perturbed and unperturbed conditions.

The simultaneous modulation of joint torque and stiffness enables humans to perform large repertoires of movements, while versatilely adapting to external mechanical demands. Multi-muscle force control is key for joint torque and stiffness modulation. However, the inability to directly measure muscle force in the intact moving human prevents understanding how muscle force causally links to joint torque and stiffness. Joint stiffness is predominantly estimated via joint perturbation-based experiments in combination with system identification techniques. However, these techniques provide joint-level stiffness estimations with no causal link to the underlying muscle forces. Moreover, the need for joint perturbations limits the generalizability and applicability to study natural movements. Here, we present an electromyography (EMG)-driven musculoskeletal modeling framework that can be calibrated to match reference joint torque and stiffness profiles simultaneously via a multi-term objective function. EMG-driven models calibrated on <2 s of reference torque and stiffness data could blindly estimate reference profiles across 100 s of data not used for calibration. Model calibrations using an objective function comprising torque and stiffness terms always provided less feasible solutions than an objective function comprising solely a torque term, thereby reducing the space of feasible muscle-tendon parameters. Results also showed the proposed framework's ability to estimate joint stiffness in unperturbed conditions, while capturing differences against stiffness profiles derived during perturbed conditions. The proposed framework may provide new ways for studying causal relationships between muscle force and joint torque and stiffness during movements in interaction with the environment, with broad implications across biomechanics, rehabilitation and robotics.