Sensitivity analysis of the estimated muscle forces during gait with respect to the musculoskeletal model parameters and dynamic simulation techniques

Predictive simulation of single-leg landing scenarios for ACL injury risk factors evaluation

The Anterior Cruciate Ligament (ACL) rupture is a very common knee injury during sport activities. Landing after jump is one of the most prominent human body movements that can lead to such an injury. The landing-related ACL injury risk factors have been in the spotlight of research interest. Over the years, researchers and clinicians acquire knowledge about human movement during daily-life activities by organizing complex in vivo studies that feature high complexity, costs and technical and most importantly physical challenges. In an attempt to overcome these limitations, this paper introduces a computational modeling and simulation pipeline that aims to predict and identify key parameters of interest that are related to ACL injury during single-leg landings. We examined the following conditions: a) landing height, b) hip internal and external rotation, c) lumbar forward and backward leaning, d) lumbar medial and lateral bending, e) muscle forces permutations and f) effort goal weight. Identified on related research studies, we evaluated the following risk factors: vertical Ground Reaction Force (vGRF), knee joint Anterior force (AF), Medial force (MF), Compressive force (CF), Abduction moment (AbdM), Internal rotation moment (IRM), quadricep and hamstring muscle forces and Quadriceps/Hamstrings force ratio (Q/H force ratio). Our study clearly demonstrated that ACL injury is a rather complicated mechanism with many associated risk factors which are evidently correlated. Nevertheless, the results were mostly in agreement with other research studies regarding the ACL risk factors. The presented pipeline showcased promising potential of predictive simulations to evaluate different aspects of complicated phenomena, such as the ACL injury.

A developed multibody knee model for unloading knee with cartilage penetration depth control

Unloader knee braces could relieve pain by decreasing the medial knee loading. Particularly for knee osteoarthritis (KOA) patients, this study investigates the relevance of the knee model after identifying the most influential parameter. Since KOA causes the cartilage in a joint to lose its elasticity and thickness, the lack of normal bone-to-bone separation can be painful. We believe that cartilage penetration depth control is an impactful strategy in this research. Moreover, the knee contact force in KOA is fewer than in healthy knees, confirming that the contact force control cannot be a straight factor. Therefore, a biomechanical human knee model is developed, and a generic procedure for dynamic analysis of contact problems in combination with the musculoskeletal model is introduced. The developed model includes the geometric expression of collision curves and an algorithm for determining collision points. This presentation addresses cartilage penetration depth and contact force calculation through nonlinear discontinuous contact law. In view of this, femur and tibia’s relative motion is analyzed through the combined collision reactions of cartilage and bone in the knee. In the simulation, maximum penetration depth in a healthy knee is reported to be 0.795 mm, while in a 75% KOA is 0.521 mm, including 0.5 mm cartilage-cartilage contact and 0.021 mm bone-bone contact. The top unloading 852 N is achieved, reducing penetration depth to 0.45 mm, avoiding bone-bone contact. This proposed procedure with low computation gives us a suitable analysis method for designing knee assistive devices.

Sensitivity analysis guided improvement of an electromyogram-driven lumped parameter musculoskeletal hand model

EMG-driven neuromusculoskeletal models have been used to study many impairments and hold great potential to facilitate human-machine interactions for rehabilitation. A challenge to successful clinical application is the need to optimize the model parameters to produce accurate kinematic predictions. In order to identify the key parameters, we used Monte-Carlo simulations to evaluate the sensitivities of wrist and metacarpophalangeal (MCP) flexion/extension prediction accuracies for an EMG-driven, lumped-parameter musculoskeletal model. Four muscles were modeled with 22 total optimizable parameters. Model predictions from EMG were compared with measured joint angles from 11 able-bodied subjects. While sensitivities varied by muscle, we determined muscle moment arms, maximum isometric force, and tendon slack length were highly influential, while passive stiffness and optimal fiber length were less influential. Removing the two least influential parameters from each muscle reduced the optimization search space from 22 to 14 parameters without significantly impacting prediction correlation (wrist: 0.90 ± 0.05 vs 0.90 ± 0.05, p = 0.96; MCP: 0.74 ± 0.20 vs 0.70 ± 0.23, p = 0.51) and normalized root mean square error (wrist: 0.18 ± 0.03 vs 0.19 ± 0.03, p = 0.16; MCP: 0.18 ± 0.06 vs 0.19 ± 0.06, p = 0.60). Additionally, we showed that wrist kinematic predictions were insensitive to parameters of the modeled MCP muscles. This allowed us to develop a novel optimization strategy that more reliably identified the optimal set of parameters for each subject (27.3 ± 19.5%) compared to the baseline optimization strategy (6.4 ± 8.1%; p = 0.004). This study demonstrated how sensitivity analyses can be used to guide model refinement and inform novel and improved optimization strategies, facilitating implementation of musculoskeletal models for clinical applications.

Uncertainty in Muscle–Tendon Parameters can Greatly Influence the Accuracy of Knee Contact Force Estimates of Musculoskeletal Models

Understanding the sources of error is critical before models of the musculoskeletal system can be usefully translated. Using in vivo measured tibiofemoral forces, the impact of uncertainty in muscle–tendon parameters on the accuracy of knee contact force estimates of a generic musculoskeletal model was investigated following a probabilistic approach. Population variability was introduced to the routine musculoskeletal modeling framework by perturbing input parameters of the lower limb muscles around their baseline values. Using ground reaction force and skin marker trajectory data collected from six subjects performing body-weight squat, the knee contact force was calculated for the perturbed models. The combined impact of input uncertainties resulted in a considerable variation in the knee contact force estimates (up to 2.1 BW change in the predicted force), especially at larger knee flexion angles, hence explaining up to 70% of the simulation error. Although individual muscle groups exhibited different contributions to the overall error, variation in the maximum isometric force and pathway of the muscles showed the highest impacts on the model outcomes. Importantly, this study highlights parameters that should be personalized in order to achieve the best possible predictions when using generic musculoskeletal models for activities involving deep knee flexion.

A review of musculoskeletal modelling of human locomotion

Locomotion through the environment is important because movement provides access to key resources, including food, shelter and mates. Central to many locomotion-focused questions is the need to understand internal forces, particularly muscle forces and joint reactions. Musculoskeletal modelling, which typically harnesses the power of inverse dynamics, unites experimental data that are collected on living subjects with virtual models of their morphology. The inputs required for producing good musculoskeletal models include body geometry, muscle parameters, motion variables and ground reaction forces. This methodological approach is critically informed by both biological anthropology, with its focus on variation in human form and function, and mechanical engineering, with a focus on the application of Newtonian mechanics to current problems. Here, we demonstrate the application of a musculoskeletal modelling approach to human walking using the data of a single male subject. Furthermore, we discuss the decisions required to build the model, including how to customize the musculoskeletal model, and suggest cautions that both biological anthropologists and engineers who are interested in this topic should consider.

Using Musculoskeletal Models to Estimate in vivo Total Knee Replacement Kinematics and Loads: Effect of Differences Between Models

Total knee replacement (TKR) is one of the most performed orthopedic surgeries to treat knee joint diseases in the elderly population. Although the survivorship of knee implants may extend beyond two decades, the poor outcome rate remains considerable. A recent computational approach used to better understand failure modes and improve TKR outcomes is based on the combination of musculoskeletal (MSK) and finite element models. This combined multiscale modeling approach is a promising strategy in the field of computational biomechanics; however, some critical aspects need to be investigated. In particular, the identification and quantification of the uncertainties related to the boundary conditions used as inputs to the finite element model due to a different definition of the MSK model are crucial. Therefore, the aim of this study is to investigate this problem, which is relevant for the model credibility assessment process. Three different generic MSK models available in the OpenSim platform were used to simulate gait, based on the experimental data from the fifth edition of the "Grand Challenge Competitions to Predict in vivo Knee Loads." The outputs of the MSK analyses were compared in terms of relative kinematics of the knee implant components and joint reaction (JR) forces and moments acting on the tibial insert. Additionally, the estimated knee JRs were compared with those measured by the instrumented knee implant so that the "global goodness of fit" was quantified for each model. Our results indicated that the different kinematic definitions of the knee joint and the muscle model implemented in the different MSK models influenced both the motion and the load history of the artificial joint. This study demonstrates the importance of examining the influence of the model assumptions on the output results and represents the first step for future studies that will investigate how the uncertainties in the MSK models propagate on disease-specific finite element model results.

An Anatomical Atlas-Based Scaling Study for Quantifying Muscle and Hip Joint Contact Forces in Above and Through-Knee Amputees Using Validated Musculoskeletal Modelling

OBJECTIVE

Customisation of musculoskeletal modelling using magnetic resonance imaging (MRI) significantly improves the model accuracy, but the process is time consuming and computationally intensive. This study hypothesizes that linear scaling to a lower limb amputee model with anthropometric similarity can accurately predict muscle and joint contact forces.

METHODS

An MRI-based anatomical atlas, comprising 18 trans-femoral and through-knee traumatic lower limb amputee models, is developed. Gait data, using a 10-camera motion capture system with two force plates, and surface electromyography (EMG) data were collected. Muscle and hip joint contact forces were quantified using musculoskeletal modelling. The predicted muscle activations from the subject-specific models were validated using EMG recordings. Anthropometry based multiple linear regression models, which minimize errors in force predictions, are presented.

RESULTS

All predictions showed excellent (error interval c = 0-0.15), very good (c = 0.15-0.30) or good (c = 0.30-0.45) similarity to the EMG data, demonstrating accurate computation of muscle activations. The primary predictors of discrepancies in force predictions were differences in pelvis width (p < 0.001), body mass index (BMI, p < 0.001) and stump length to pelvis width ratio (p < 0.001) between the respective individual and underlying dataset.

CONCLUSION

Linear scaling to a model with the most similar pelvis width, BMI and stump length to pelvis width ratio results in modelling outcomes with minimal errors.

SIGNIFICANCE

This study provides robust tools to perform accurate analyses of musculoskeletal mechanics for high-functioning lower limb military amputees, thus facilitating the further understanding and improvement of the amputee's function. The atlas is available in an open source repository.

A Conceptual Blueprint for Making Neuromusculoskeletal Models Clinically Useful

The ultimate goal of most neuromusculoskeletal modeling research is to improve the treatment of movement impairments. However, even though neuromusculoskeletal models have become more realistic anatomically, physiologically, and neurologically over the past 25 years, they have yet to make a positive impact on the design of clinical treatments for movement impairments. Such impairments are caused by common conditions such as stroke, osteoarthritis, Parkinson’s disease, spinal cord injury, cerebral palsy, limb amputation, and even cancer. The lack of clinical impact is somewhat surprising given that comparable computational technology has transformed the design of airplanes, automobiles, and other commercial products over the same time period. This paper provides the author’s personal perspective for how neuromusculoskeletal models can become clinically useful. First, the paper motivates the potential value of neuromusculoskeletal models for clinical treatment design. Next, it highlights five challenges to achieving clinical utility and provides suggestions for how to overcome them. After that, it describes clinical, technical, collaboration, and practical needs that must be addressed for neuromusculoskeletal models to fulfill their clinical potential, along with recommendations for meeting them. Finally, it discusses how more complex modeling and experimental methods could enhance neuromusculoskeletal model fidelity, personalization, and utilization. The author hopes that these ideas will provide a conceptual blueprint that will help the neuromusculoskeletal modeling research community work toward clinical utility.

Muscle Synergy Constraints Do Not Improve Estimates of Muscle Activity From Static Optimization During Gait for Unimpaired Children or Children With Cerebral Palsy

Neuromusculoskeletal simulation provides a promising platform to inform the design of assistive devices or inform rehabilitation. For these applications, a simulation must be able to accurately represent the person of interest, such as an individual with a neurologic injury. If a simulation fails to predict how an individual recruits and coordinates their muscles during movement, it will have limited utility for informing design or rehabilitation. While inverse dynamic simulations have previously been used to evaluate anticipated responses from interventions, like orthopedic surgery or orthoses, they frequently struggle to accurately estimate muscle activations, even for tasks like walking. The simulated muscle activity often fails to represent experimentally measured muscle activity from electromyographic (EMG) recordings. Research has theorized that the nervous system may simplify the range of possible activations used during dynamic tasks, by constraining activations to weighted groups of muscles, referred to as muscle synergies. Synergies are altered after neurological injury, such as stroke or cerebral palsy (CP), and may provide a method for improving subject-specific models of neuromuscular control. The aim of this study was to test whether constraining simulation to synergies could improve estimated muscle activations compared to EMG data. We evaluated modeled muscle activations during gait for six typically developing (TD) children and six children with CP. Muscle activations were estimated with: (1) static optimization (SO), minimizing muscle activations squared, and (2) synergy SO (SynSO), minimizing synergy activations squared using the weights identified from EMG data for two to five synergies. While SynSO caused changes in estimated activations compared to SO, the correlation to EMG data was not higher in SynSO than SO for either TD or CP groups. The correlations to EMG were higher in CP than TD for both SO (CP: 0.48, TD: 0.36) and SynSO (CP: 0.46, TD: 0.26 for five synergies). Constraining activations to SynSO caused the simulated muscle stress to increase compared to SO for all individuals, causing a 157% increase with two synergies. These results suggest that constraining simulated activations in inverse dynamic simulations to subject-specific synergies alone may not improve estimation of muscle activations during gait for generic musculoskeletal models.