Is my model good enough? Best practices for verification and validation of musculoskeletal models and simulations of movement

Lower Extremity Joint Contributions to Trunk Control During Walking in Persons with Transtibial Amputation

Controlled trunk motion is crucial for balance and stability during walking. Persons with lower extremity amputation often exhibit abnormal trunk motion, yet underlying mechanisms are not well understood nor have optimal clinical interventions been established. The aim of this work was to characterize associations between altered lower extremity joint moments and altered trunk dynamics in persons with unilateral, transtibial amputation (TTA). Full-body gait data were collected from 10 persons with TTA and 10 uninjured persons walking overground (~1.4 m/s). Experimentally-measured trunk angular accelerations were decomposed into constituent accelerations caused by net joint moments throughout the body using an induced acceleration analysis. Results showed persons with TTA had similar ankle moment magnitude relative to uninjured persons (P > 0.05), but greater trunk angular acceleration induced by the prosthetic ankle which acted to lean the trunk ipsilaterally (P = 0.003). Additionally, persons with TTA had a reduced knee extensor moment relative to uninjured persons (P < 0.001), resulting in lesser sagittal and frontal induced trunk angular accelerations (P < 0.001). These data indicate kinetic compensations at joints other than the lumbar and hip contribute to altered trunk dynamics in persons with a unilateral TTA. Findings may inform development of new clinical strategies to modify problematic trunk motion.

Musculoskeletal models of a human and bonobo finger: parameter identification and comparison to in vitro experiments

Introduction

Knowledge of internal finger loading during human and non-human primate activities such as tool use or knuckle-walking has become increasingly important to reconstruct the behaviour of fossil hominins based on bone morphology. Musculoskeletal models have proven useful for predicting these internal loads during human activities, but load predictions for non-human primate activities are missing due to a lack of suitable finger models. The main goal of this study was to implement both a human and a representative non-human primate finger model to facilitate comparative studies on metacarpal bone loading. To ensure that the model predictions are sufficiently accurate, the specific goals were: (1) to identify species-specific model parameters based on in vitro measured fingertip forces resulting from single tendon loading and (2) to evaluate the model accuracy of predicted fingertip forces and net metacarpal bone loading in a different loading scenario.

Materials & Methods

Three human and one bonobo (Pan paniscus) fingers were tested in vitro using a previously developed experimental setup. The cadaveric fingers were positioned in four static postures and load was applied by attaching weights to the tendons of the finger muscles. For parameter identification, fingertip forces were measured by loading each tendon individually in each posture. For the evaluation of model accuracy, the extrinsic flexor muscles were loaded simultaneously and both the fingertip force and net metacarpal bone force were measured. The finger models were implemented using custom Python scripts. Initial parameters were taken from literature for the human model and own dissection data for the bonobo model. Optimized model parameters were identified by minimizing the error between predicted and experimentally measured fingertip forces. Fingertip forces and net metacarpal bone loading in the combined loading scenario were predicted using the optimized models and the remaining error with respect to the experimental data was evaluated.

Results

The parameter identification procedure led to minor model adjustments but considerably reduced the error in the predicted fingertip forces (root mean square error reduced from 0.53/0.69 N to 0.11/0.20 N for the human/bonobo model). Both models remained physiologically plausible after the parameter identification. In the combined loading scenario, fingertip and net metacarpal forces were predicted with average directional errors below 6° and magnitude errors below 12%.

Conclusions

This study presents the first attempt to implement both a human and non-human primate finger model for comparative palaeoanthropological studies. The good agreement between predicted and experimental forces involving the action of extrinsic flexors—which are most relevant for forceful grasping—shows that the models are likely sufficiently accurate for comparisons of internal loads occurring during human and non-human primate manual activities.

Episiotomy: the biomechanical impact of multiple small incisions during a normal vaginal delivery

Childbirth-related injuries are one of the main causes of pelvic floor dysfunction. To attempt to avoid serious tears during delivery, an episiotomy can be performed. In this study, we intended to investigate the biomechanical performance of the pelvic floor muscles after performing different episiotomies using a physics-based computational model which includes the pelvic floor muscles and the fetus. Previous biomechanical studies have analysed the mechanical effects of single incisions of different lengths; in this study, we intend to analyse the implications of multiple small incisions, evaluating the reaction forces, the stress on the muscles and the loss of tissue integrity sustained by the pelvic floor. The obtained results predict that an episiotomy delivery reduces the likelihood of macroscopic levator trauma by decreasing the stress on the region of insertion of the rectal area of the levator ani in the symphysis pubis. From the mechanical point of view, multiple incisions do not bring benefits compared to larger incisions. However, nothing can be ascertained about the clinical benefit of such an approach.

Musculoskeletal full-body models including a detailed thoracolumbar spine for children and adolescents aged 6-18 years

Currently available musculoskeletal inverse-dynamics thoracolumbar spine models are entirely based on data from adults and might therefore not be applicable for simulations in children and adolescents. In addition, these models lack lower extremities, which are required for comprehensive evaluations of functional activities or therapeutic exercises. We therefore created OpenSim-based musculoskeletal full-body models including a detailed thoracolumbar spine for children and adolescents aged 6-18 years and validated by comparing model predictions to in vivo data. After combining our recently developed adult thoracolumbar spine model with a lower extremity model, children and adolescent models were created for each year of age by adjusting segmental length and mass distribution, center of mass positions and moments of inertia of the major body segments as well as sagittal pelvis and spine alignment based on literature data. Similarly, muscle strength properties were adjusted based on CT-derived cross-sectional area measurements. Simulations were conducted from in vivo studies reported in the literature involving children and adolescents evaluating maximum trunk muscle strength (MTMS), lumbar disc compressibility (LDC), intradiscal pressure (IDP) and trunk muscle activity (MA). Model predictions correlated highly with in vivo data (MTMS: r ≥ 0.82, p ≤ 0.03; LDC: r = 0.77, p < 0.001; IDP: r ≥ 0.78, p < 0.001; MA: r ≥ 0.90, p < 0.001), indicating suitability for the reasonably accurate prediction of maximal trunk muscle strength, segmental loading and trunk muscle activity in children and adolescents. When aiming at investigating children or adolescents with pathologies such as idiopathic scoliosis, our models can serve as a basis for the creation of deformed spine models and for comparative purposes.

SimCP: A Simulation Platform to Predict Gait Performance Following Orthopedic Intervention in Children With Cerebral Palsy

Gait deficits in cerebral palsy (CP) are often treated with a single-event multi-level surgery (SEMLS). Selecting the treatment options (combination of bony and soft tissue corrections) for a specific patient is a complex endeavor and very often treatment outcome is not satisfying. A deterioration in 22.8% of the parameters describing gait performance has been reported and there is need for additional surgery in 11% of the patients. Computational simulations based on musculoskeletal models that allow clinicians to test the effects of different treatment options before surgery have the potential to drastically improve treatment outcome. However, to date, no such simulation and modeling method is available. Two important challenges are the development of methods to include patient-specific neuromechanical impairments into the models and to simulate the effect of different surgical procedures on post-operative gait performance. Therefore, we developed the SimCP framework that allows the evaluation of the effect of different simulated surgeries on gait performance of a specific patient and includes a graphical user interface (GUI) that enables performing virtual surgery on the models. We demonstrated the potential of our framework for two case studies. Models reflecting the patient-specific musculoskeletal geometry and muscle properties are generated based solely on data collected before the treatment. The patient's motor control is described based on muscle synergies derived from pre-operative EMG. The GUI is then used to modify the musculoskeletal properties according to the surgical plan. Since SEMLS does not affect motor control, the same motor control model is used to define gait performance pre- and post-operative. We use the capability gap (CG), i.e., the difference between the joint moments needed to perform healthy walking and the joint moments the personalized model can generate, to quantify gait performance. In both cases, the CG was smaller post- then pre-operative and this was in accordance with the measured change in gait kinematics after treatment.

An Optimization Method Tracking EMG, Ground Reactions Forces, and Marker Trajectories for Musculo-Tendon Forces Estimation in Equinus Gait

In the context of neuro-orthopedic pathologies affecting walking and thus patients' quality of life, understanding the mechanisms of gait deviations and identifying the causal motor impairments is of primary importance. Beside other approaches, neuromusculoskeletal simulations may be used to provide insight into this matter. To the best of our knowledge, no computational framework exists in the literature that allows for predictive simulations featuring muscle co-contractions, and the introduction of various types of perturbations during both healthy and pathological gait types. The aim of this preliminary study was to adapt a recently proposed EMG-marker tracking optimization process to a lower limb musculoskeletal model during equinus gait, a multiphase problem with contact forces. The resulting optimization method tracking EMG, ground reactions forces, and marker trajectories allowed an accurate reproduction of joint kinematics (average error of 5.4 ± 3.3 mm for pelvis translations, and 1.9 ± 1.3° for pelvis rotation and joint angles) and ensured good temporal agreement in muscle activity (the concordance between estimated and measured excitations was 76.8 ± 5.3 %) in a relatively fast process (3.88 ± 1.04 h). We have also highlighted that the tracking of ground reaction forces was possible and accurate (average error of 17.3 ± 5.5 N), even without the use of a complex foot-ground contact model.

Multibody Kinematics Optimization for the Estimation of Upper and Lower Limb Human Joint Kinematics: A Systematized Methodological Review

Multibody kinematics optimization (MKO) aims to reduce soft tissue artefact (STA) and is a key step in musculoskeletal modeling. The objective of this review was to identify the numerical methods, their validation and performance for the estimation of the human joint kinematics using MKO. Seventy-four papers were extracted from a systematized search in five databases and cross-referencing. Model-derived kinematics were obtained using either constrained optimization or Kalman filtering to minimize the difference between measured (i.e., by skin markers, electromagnetic or inertial sensors) and model-derived positions and/or orientations. While hinge, universal, and spherical joints prevail, advanced models (e.g., parallel and four-bar mechanisms, elastic joint) have been introduced, mainly for the knee and shoulder joints. Models and methods were evaluated using: (i) simulated data based, however, on oversimplified STA and joint models; (ii) reconstruction residual errors, ranging from 4 mm to 40 mm; (iii) sensitivity analyses which highlighted the effect (up to 36 deg and 12 mm) of model geometrical parameters, joint models, and computational methods; (iv) comparison with other approaches (i.e., single body kinematics optimization and nonoptimized kinematics); (v) repeatability studies that showed low intra- and inter-observer variability; and (vi) validation against ground-truth bone kinematics (with errors between 1 deg and 22 deg for tibiofemoral rotations and between 3 deg and 10 deg for glenohumeral rotations). Moreover, MKO was applied to various movements (e.g., walking, running, arm elevation). Additional validations, especially for the upper limb, should be undertaken and we recommend a more systematic approach for the evaluation of MKO. In addition, further model development, scaling, and personalization methods are required to better estimate the secondary degrees-of-freedom (DoF).

Comprehensive Uncertainty Quantification and Sensitivity Analysis for Cardiac Action Potential Models

Recent efforts to ensure the reliability of computational model-based predictions in healthcare, such as the ASME V&V40 Standard, emphasize the importance of uncertainty quantification (UQ) and sensitivity analysis (SA) when evaluating computational models. UQ involves empirically determining the uncertainty in model inputs-typically resulting from natural variability or measurement error-and then calculating the resultant uncertainty in model outputs. SA involves calculating how uncertainty in model outputs can be apportioned to input uncertainty. Rigorous comprehensive UQ/SA provides confidence that model-based decisions are robust to underlying uncertainties. However, comprehensive UQ/SA is not currently feasible for whole heart models, due to numerous factors including model complexity and difficulty in measuring variability in the many parameters. Here, we present a significant step to developing a framework to overcome these limitations. We: (i) developed a novel action potential (AP) model of moderate complexity (six currents, seven variables, 36 parameters); (ii) prescribed input variability for all parameters (not empirically derived); (iii) used a single "hyper-parameter" to study increasing levels of parameter uncertainty; (iv) performed UQ and SA for a range of model-derived quantities with physiological relevance; and (v) present quantitative and qualitative ways to analyze different behaviors that occur under parameter uncertainty, including "model failure". This is the first time uncertainty in every parameter (including conductances, steady-state parameters, and time constant parameters) of every ionic current in a cardiac model has been studied. This approach allowed us to demonstrate that, for this model, the simulated AP is fully robust to low levels of parameter uncertainty - to our knowledge the first time this has been shown of any cardiac model. A range of dynamics was observed at larger parameter uncertainty (e.g., oscillatory dynamics); analysis revealed that five parameters were highly influential in these dynamics. Overall, we demonstrate feasibility of performing comprehensive UQ/SA for cardiac cell models and demonstrate how to assess robustness and overcome model failure when performing cardiac UQ analyses. The approach presented here represents an important and significant step toward the development of model-based clinical tools which are demonstrably robust to all underlying uncertainties and therefore more reliable in safety-critical decision-making.

An actuated dissipative spring-mass walking model: Predicting human-like ground reaction forces and the effects of model parameters

Simple models are widely used to understand the mechanics of human walking. The optimization-based minimal biped model and spring-loaded-inverted-pendulum (SLIP) model are two popular models that can achieve human-like walking patterns. However, ground reaction forces (GRF) from these two models still deviate from experimental data. In this paper, we proposed an actuated dissipative spring-mass model by integrating these two models to realize more human-like GRF patterns. We first explored the function of stiffness, damping, and weights of both energy cost and force cost in the objective function and found that these parameters have distinctly different influences on the optimized gait and GRF profiles. The stiffness and objective weight affect the number and size of peaks in the vertical GRF and stance time. The damping changes the relative size of the peaks but has little influence on stance time. Based on these observations, these parameters were manually tuned at three different speeds to approach experimentally measured vertical GRF and the highest correlation coefficient can reach 0.983. These results indicate that the stiffness, damping, and proper objective functions are all important factors in achieving human-like motion for this simple walking model. These findings can facilitate the understanding of human walking dynamics and may be applied in future biped models.

Assessment of the energy-related cost function over a range of walking speeds

Cost funtions are needed for calculation of muscle forces in musculoskeletal models. The behavior of the energy-related cost function, proposed by Praagman et al. (J Biomech 39(4):758-765, 2006. https://doi.org/10.1016/j.jbiomech.2004.11.034 ) (CFP), can be used as an optimization criteria in musculoskeletal models for studying gait. In particular, in this work, its performance is compared against two empirical phenomenological models at different walking speed conditions. Also, the sensitivity of the CFP function to model parameters, such as muscle mass, maximal isometric muscle force, optimal muscle fiber length and maximum muscle velocity of the contractile element, was analyzed. The obtained results showed that CFP presents different behavior (in terms of the normalized root-mean-squared deviation (NRMSD) and the coefficient of multiple correlation (CMC)) for different muscles. Also, it provided estimates with median of NRMSD between 0.176 and 0.299 and median of CMC between 0.703 and 0.865 both metrics for slow, free and fast walking speed, which could be considered as acceptable results. Furthermore, the results indicated that CFP is insensitive to changes in muscle mass and relatively sensitive to maximal isometric muscle force. However, CFP presented a noisy behavior on estimations of muscle energy rate for some muscle as compared to phenomenological models. Finally, estimations by CFP during gait are within the values obtained by the empirical phenomenological models.

Contributions of muscles and external forces to medial knee load reduction due to osteoarthritis braces

BACKGROUND

Braces for medial knee osteoarthritis can reduce medial joint loads through a combination of three mechanisms: application of an external brace abduction moment, alteration of gait dynamics, and reduced activation of antagonistic muscles. Although the effect of knee bracing has been reported independently for each of these parameters, no previous study has quantified their relative contributions to reducing medial knee loads.

METHODS

In this study, we used a detailed musculoskeletal model to investigate immediate changes in medial and lateral loads caused by two different knee braces: OA Assist and OA Adjuster 3 (DJO Global). Seventeen osteoarthritis subjects and eighteen healthy controls performed overground gait trials in unbraced and braced conditions.

RESULTS

Across all subjects, bracing reduced medial loads by 0.1 to 0.3 times bodyweight (BW), or roughly 10%, and increased lateral loads by 0.03 to 0.2 BW. Changes in gait kinematics due to bracing were subtle, and had little effect on medial and lateral joint loads. The knee adduction moment was unaltered unless the brace moment was included in its computation. Only one muscle, biceps femoris, showed a significant change in EMG with bracing, but this did not contribute to altered peak medial contact loads.

CONCLUSIONS

Knee braces reduced medial tibiofemoral loads primarily by applying a direct, and substantial, abduction moment to each subject's knee. To further enhance brace effectiveness, future brace designs should seek to enhance the magnitude of this unloader moment, and possibly exploit additional kinematic or neuromuscular gait modifications.

Mass distribution of wheelchair athletes assessed using DXA scans and biomechanical simulations

The anthropometries of elite wheelchair racing athletes differ to the generic, able-bodied anthropometries commonly used in computational biomechanical simulations. The impact of using able-bodied parameters on the accuracy of simulations involving wheelchair racing is currently unknown. In this study, athlete-specific mass segment inertial parameters of five elite wheelchair athletes were calculated using dual-energy X-ray absorptiometry scans. These were compared against commonly used anthropometrics parameters of data presented in the literature. A computational biomechanical simulation of wheelchair propulsion assessed the sensitivity of athlete-specific mass parameters using Kruskal-Wallis analysis, Mann-Whitney U analysis and Spearman correlations. Substantial between-athlete body mass distribution variances (thigh mass < 14.6% total body mass), and between-limb asymmetries (<62.4%; 3.1 kg) were observed. Compared to non-athletic able-bodied anthropometric data, wheelchair racing athletes demonstrated greater mass in the upper extremities (up to 3.8% total body mass), and less in the lower extremities (up to 9.8% total body mass). Computational simulations were sensitive to individual body mass distribution, with simulation outputs increasing by up to 12.5% when measured segment masses were 14.3% greater than the generic counterpart. These data suggest non-athletic, able-bodied mass segment inertial parameters are inappropriate for analysing elite wheelchair racing motion.

Musculoskeletal modelling of the human cervical spine for the investigation of injury mechanisms during axial impacts

Head collisions in sport can result in catastrophic injuries to the cervical spine. Musculoskeletal modelling can help analyse the relationship between motion, external forces and internal loads that lead to injury. However, impact specific musculoskeletal models are lacking as current viscoelastic values used to describe cervical spine joint dynamics have been obtained from unrepresentative quasi-static or static experiments. The aim of this study was to develop and validate a cervical spine musculoskeletal model for use in axial impacts. Cervical spine specimens (C2-C6) were tested under measured sub-catastrophic loads and the resulting 3D motion of the vertebrae was measured. Specimen specific musculoskeletal models were then created and used to estimate the axial and shear viscoelastic (stiffness and damping) properties of the joints through an optimisation algorithm that minimised tracking errors between measured and simulated kinematics. A five-fold cross validation and a Monte Carlo sensitivity analysis were conducted to assess the performance of the newly estimated parameters. The impact-specific parameters were integrated in a population specific musculoskeletal model and used to assess cervical spine loads measured from Rugby union impacts compared to available models. Results of the optimisation showed a larger increase of axial joint stiffness compared to axial damping and shear viscoelastic parameters for all models. The sensitivity analysis revealed that lower values of axial stiffness and shear damping reduced the models performance considerably compared to other degrees of freedom. The impact-specific parameters integrated in the population specific model estimated more appropriate joint displacements for axial head impacts compared to available models and are therefore more suited for injury mechanism analysis.

Neuro-musculoskeletal flexible multibody simulation yields a framework for efficient bone failure risk assessment

Fragility fractures are a major socioeconomic problem. A non-invasive, computationally-efficient method for the identification of fracture risk scenarios under the representation of neuro-musculoskeletal dynamics does not exist. We introduce a computational workflow that integrates modally-reduced, quantitative CT-based finite-element models into neuro-musculoskeletal flexible multibody simulation (NfMBS) for early bone fracture risk assessment. Our workflow quantifies the bone strength via the osteogenic stresses and strains that arise due to the physiological-like loading of the bone under the representation of patient-specific neuro-musculoskeletal dynamics. This allows for non-invasive, computationally-efficient dynamic analysis over the enormous parameter space of fracture risk scenarios, while requiring only sparse clinical data. Experimental validation on a fresh human femur specimen together with femur strength computations that were consistent with literature findings provide confidence in the workflow: The simulation of an entire squat took only 38 s CPU-time. Owing to the loss (16% cortical, 33% trabecular) of bone mineral density (BMD), the strain measure that is associated with bone fracture increased by 31.4%; and yielded an elevated risk of a femoral hip fracture. Our novel workflow could offer clinicians with decision-making guidance by enabling the first combined in-silico analysis tool using NfMBS and BMD measurements for optimized bone fracture risk assessment.

Gait Analysis of Patients Subjected to the Atrophic Mandible Augmentation with Iliac Bone Graft

In this study, we aimed to quantitatively monitor and describe the gait functions of patients, who underwent iliac crest bone grafting in atrophic jaw augmentation operation, by taking into account the alterations of gait parameters and muscle forces in the early recovery course. To do so, temporospatial and kinematic gait parameters of ten patients during pre- and postoperative periods were recorded, and forces of the gluteus medius, gluteus maximus, and iliacus muscles were calculated. Three postoperative periods were specified as one week (post-op1), two weeks (post-op2), and three weeks (post-op3) after the surgery. Restoring process of the gait patterns was comparatively evaluated by analyzing the gait parameters and muscle forces for pre- and postoperative periods. Temporospatial and kinematic parameters of post-op3 were closer to those obtained in pre-op than those in post-op1 and post-op2 (p < 0.05). Muscle forces calculated in post-op3 showed the best agreement with those in pre-op among the postoperative periods in terms of both magnitude and correlation (p < 0.05). In conclusion, the patients began to regain their preoperative gait characteristics from the second week after surgery, but complete recovery in gait was observed three weeks after the surgery.

Neuromusculoskeletal Simulation Reveals Abnormal Rectus Femoris-Gluteus Medius Coupling in Post-stroke Gait

Post-stroke gait is often accompanied by muscle impairments that result in adaptations such as hip circumduction to compensate for lack of knee flexion. Our previous work robotically enhanced knee flexion in individuals post-stroke with Stiff-Knee Gait (SKG), however, this resulted in greater circumduction, suggesting the existence of abnormal coordination in SKG. The purpose of this work is to investigate two possible mechanisms of the abnormal coordination: (1) a reflex coupling between stretched quadriceps and abductors, and (2) a coupling between volitionally activated knee flexors and abductors. We used previously collected kinematic, kinetic and EMG measures from nine participants with chronic stroke and five healthy controls during walking with and without the applied knee flexion torque perturbations in the pre-swing phase of gait in the neuromusculoskeletal simulation. The measured muscle activity was supplemented by simulated muscle activations to estimate the muscle states of the quadriceps, hamstrings and hip abductors. We used linear mixed models to investigate two hypotheses: (H1) association between quadriceps and abductor activation during an involuntary period (reflex latency) following the perturbation and (H2) association between hamstrings and abductor activation after the perturbation was removed. We observed significantly higher rectus femoris (RF) activation in stroke participants compared to healthy controls within the involuntary response period following the perturbation based on both measured (H1, p < 0.001) and simulated (H1, p = 0.022) activity. Simulated RF and gluteus medius (GMed) activations were correlated only in those with SKG, which was significantly higher compared to healthy controls (H1, p = 0.030). There was no evidence of synergistic coupling between any combination of hamstrings and hip abductors (H2, p > 0.05) when the perturbation was removed. The RF-GMed coupling suggests an underlying abnormal coordination pattern in post-stroke SKG, likely reflexive in origin. These results challenge earlier assumptions that hip circumduction in stroke is simply a kinematic adaptation due to reduced toe clearance. Instead, abnormal coordination may underlie circumduction, illustrating the deleterious role of abnormal coordination in post-stroke gait.

Best Practices to Maximize the Use and Reuse of Quantitative and Systems Pharmacology Models: Recommendations From the United Kingdom Quantitative and Systems Pharmacology Network

The lack of standardization in the way that quantitative and systems pharmacology (QSP) models are developed, tested, and documented hinders their reproducibility, reusability, and expansion or reduction to alternative contexts. This in turn undermines the potential impact of QSP in academic, industrial, and regulatory frameworks. This article presents a minimum set of recommendations from the UK Quantitative and Systems Pharmacology Network (UK QSP Network) to guide QSP practitioners seeking to maximize their impact, and stakeholders considering the use of QSP models in their environment.

Spatial Dependency of Glenohumeral Joint Stability during Dynamic Unimanual and Bimanual Pushing and Pulling

Degenerative wear to the glenoid from repetitive loading can reduce effective concavity depth and lead to future instability. Workspace design should consider glenohumeral stability to prevent initial wear. While stability has been previously explored for activities of daily living including push-pull tasks, whether stability is spatially dependent is unexplored. We simulated bimanual and unimanual push-pull tasks to 4 horizontal targets (planes of elevation: 0º, 45º, 90º, and 135º) at 90º thoracohumeral elevation and 3 elevation targets (thoracohumeral elevations: 20º, 90º, 170º) at 90º plane of elevation. The 45º horizontal target was most stable regardless of exertion type and would be the ideal target placement when considering stability. This target is likely more stable because the applied load acts perpendicular to the glenoid, limiting shear force production. The 135º horizontal target was particularly unstable for unimanual pushing (143% less stable than the 45º target), and the applied force acts parallel to the glenoid, likely creating shear forces. Pushing was less stable than pulling (all targets except sagittal 170º for both task types and horizontal 45º for bimanual) (p<0.01), which is consistent with prior reports. For example, unimanual pushing at the 90º horizontal target was 197% less stable than unimanual pulling. There were limited stability benefits to task placement for pushing, and larger stability benefits may be seen from converting pushing to pulling rather than optimizing task layout. There was no difference in stability between bimanual and unimanual tasks, suggesting no stability benefit to bimanual operation.

In vivo prediction of temporomandibular joint disc thickness and position changes for different jaw positions

Temporomandibular joint disorders (TMD) are common dysfunctions of the masticatory region and are often linked to dislocation or changes of the temporomandibular joint (TMJ) disc. Magnetic resonance imaging (MRI) is the gold standard for TMJ imaging but standard clinical sequences do not deliver a sufficient resolution and contrast for the creation of detailed meshes of the TMJ disc. Additionally, bony structures cannot be captured appropriately using standard MRI sequences due to their low signal intensity. The objective of this study was to enable researchers to create high resolution representations of all structures of the TMJ and consequently investigate morphological as well as positional changes of the masticatory system. To create meshes of the bony structures, a single computed tomography (CT) scan was acquired. In addition, a high-resolution MRI sequence was produced, which is used to collect the thickness and position change of the disc for various static postures using bite blocks. Changes in thickness of the TMJ disc as well as disc translation were measured. The newly developed workflow successfully allows researchers to create high resolution models of all structures of the TMJ for various static positions, enabling the investigation of TMJ disc translation and deformation. Discs were thinnest in the lateral part and moved mainly anteriorly and slightly medially. The procedure offers the most comprehensive picture of disc positioning and thickness changes reported to date. The presented data can be used for the development of a biomechanical computer model of TMJ anatomy and to investigate dynamic and static loads on the components of the system, which could be useful for the prediction of TMD onset.

Musculoskeletal models with generic and subject-specific geometry estimate different joint biomechanics in dysplastic hips

Optimizing the geometric complexity of musculoskeletal models is important for reliable yet feasible estimation of joint biomechanics. This study investigated the effects of subject-specific model geometry on hip joint reaction forces (JRFs) and muscle forces in patients with developmental dysplasia of the hip (DDH) and healthy controls. For nine DDH and nine control subjects, three models were created with increasingly subject-specific pelvis geometry, hip joint center locations and muscle attachments. Hip JRFs and muscle forces during a gait cycle were compared among the models. For DDH subjects, resultant JRFs from highly specific models including subject-specific pelvis geometry, joint locations and muscle attachments were not significantly different compared to models using generic geometry in early stance, but were significantly higher in late stance (p = 0.03). Estimates from moderately specific models using CT-informed scaling of generic pelvis geometry were not significantly different from low specificity models using generic geometry scaled with skin markers. For controls, resultant JRFs in early stance from highly specific models were significantly lower than moderate and low specificity models (p ≤ 0.02) with no significant differences in late stance. Inter-model JRF differences were larger for DDH subjects than controls. Inter-model differences for JRF components and muscle forces were similar to resultant JRFs. Incorporating subject-specific pelvis geometry significantly affects JRF and muscle force estimates in both DDH and control groups, which may be especially important for reliable estimation of pathomechanics in dysplastic hips.

Hip joint muscle forces during gait in patients with femoroacetabular impingement syndrome are associated with patient reported outcomes and cartilage composition

Femoroacetabular impingement syndrome (FAIS) consists of abnormal hip joint morphology and pain during activities of daily living. Abnormal gait mechanics and potentially abnormal muscle forces within FAI patients leads to articular cartilage damage. Therefore, there is a necessity to understand the effects of FAI on hip joint muscle forces during gait and the link between muscle forces, patient reported outcomes (PRO) and articular cartilage health. The purposes of this study were to assess: (1) hip muscle forces between FAI patients and healthy controls and (2) the associations between hip muscle forces with PRO and cartilage composition (T_1ρ/T₂ mapping) within FAI patients. Musculoskeletal simulations were used to estimate peak muscle forces during the stance phase of gait in 24 FAI patients and 24 healthy controls. Compared to controls, FAI patients ambulated with lower vasti (30% body-weight, p = 0.01) and higher sartorius (4.0% body-weight, p < 0.01) forces. Within FAI patients, lower peak gluteus medius, gluteus minimus, sartorius and iliopsoas forces were associated with worse hip joint pain and function (R = 0.43-0.70, p = 0-0.04), while lower muscle forces were associated with increased T_1ρ and T₂ values (i.e. altered cartilage composition) within the hip joint cartilage (R = -0.44 to -0.58, p = 0.006-0.05). Although FAI patients demonstrate abnormal muscle forces, it is unknown whether or not these altered muscle force patterns are associated with pain avoidance or weak musculature. Further investigation is required in order to better understand the effects of FAI on hip joint muscle forces and the associations with hip joint cartilage degeneration.

OpenSim Versus Human Body Model: A Comparison Study for the Lower Limbs During Gait

Musculoskeletal modeling and simulations have become popular tools for analyzing human movements. However, end users are often not aware of underlying modeling and computational assumptions. This study investigates how these assumptions affect biomechanical gait analysis outcomes performed with Human Body Model and the OpenSim gait2392 model. The authors compared joint kinematics, kinetics, and muscle forces resulting from processing data from 7 healthy adults with both models. Although outcome variables had similar patterns, there were statistically significant differences in joint kinematics (maximal difference: 9.8° [1.5°] in sagittal plane hip rotation), kinetics (maximal difference: 0.36 [0.10] N·m/kg in sagittal plane hip moment), and muscle forces (maximal difference: 8.51 [1.80] N/kg for psoas). These differences might be explained by differences in hip and knee joint center locations up to 2.4 (0.5) and 1.9 (0.2) cm in the posteroanterior and inferosuperior directions, respectively, and by the offset in pelvic reference frames of about 10° around the mediolateral axis. The choice of model may not influence the conclusions in clinical settings, where the focus is on interpreting deviations from the reference data, but it will affect the conclusions of mechanical analyses in which the goal is to obtain accurate estimates of kinematics and loading.

Comparison between kinetic and kinetic-kinematic driven knee joint finite element models

Use of knee joint finite element models for diagnostic purposes is challenging due to their complexity. Therefore, simpler models are needed for studies where a high number of patients need to be analyzed, without compromising the results of the model. In this study, more complex, kinetic (forces and moments) and simpler, kinetic-kinematic (forces and angles) driven finite element models were compared during the stance phase of gait. Patella and tendons were included in the most complex model, while they were absent in the simplest model. The greatest difference between the most complex and simplest models was observed in the internal-external rotation and axial joint reaction force, while all other rotations, translations and joint reaction forces were similar to one another. In terms of cartilage stresses and strains, the simpler models behaved similarly with the more complex models in the lateral joint compartment, while minor differences were observed in the medial compartment at the beginning of the stance phase. We suggest that it is feasible to use kinetic-kinematic driven knee joint models with a simpler geometry in studies with a large cohort size, particularly when analyzing cartilage responses and failures related to potential overloads.

Limb Kinematics, Kinetics and Muscle Dynamics During the Sit-to-Stand Transition in Greyhounds

Standing up from a prone position is a critical daily activity for animals: failing to do so effectively may cause an injurious fall or increase predation susceptibility. This sit-to-stand behaviour (StS) is biomechanically interesting because it necessitates transitioning through near-maximal joint motion ranges from a crouched (i.e., poor mechanical advantage) to a more upright posture. Such large joint excursions should require large length changes of muscle-tendon units. Here we integrate experimental and musculoskeletal simulation methods to quantify the joint motions, limb forces, and muscle fibre forces, activations and length changes during StS in an extreme athlete-the greyhound-which has large hindlimb muscles bearing short-fibred distal muscles and long tendons. Study results indicate that hindlimb anti-gravity muscle fibres operate near their ~50% limits of length change during StS; mostly by starting at highly lengthened positions. StS also requires high muscle activations (>50%), in part due to non-sagittal motions. Finally, StS movements require passive non-muscular support in the distal hindlimb where short-fibred muscles are incapable of sustaining StS themselves. Non-locomotor behaviours like StS likely impose important trade-offs between muscle fibre force capacity and length changes, as well as active and passive mechanisms of support, that have been neglected in locomotor biomechanics studies.

Impact of antagonistic muscle co-contraction on in vivo knee contact forces

Background

The onset and progression of osteoarthritis, but also the wear and loosening of the components of an artificial joint, are commonly associated with mechanical overloading of the structures. Knowledge of the mechanical forces acting at the joints, together with an understanding of the key factors that can alter them, are critical to develop effective treatments for restoring joint function. While static anatomy is usually the clinical focus, less is known about the impact of dynamic factors, such as individual muscle recruitment, on joint contact forces.

Methods

In this study, instrumented knee implants provided accurate in vivo tibio-femoral contact forces in a unique cohort of 9 patients, which were used as input for subject specific musculoskeletal models, to quantify the individual muscle forces during walking and stair negotiation.

Results

Even between patients with a very similar self-selected gait speed, the total tibio-femoral peak forces varied 1.7-fold, but had only weak correlation with static alignment (varus/valgus). In some patients, muscle co-contraction of quadriceps and gastrocnemii during walking added up to 1 bodyweight (~ 50%) to the peak tibio-femoral contact force during late stance. The greatest impact of co-contraction was observed in the late stance phase of stair ascent, with an increase of the peak tibio-femoral contact force by up to 1.7 bodyweight (66%).

Conclusions

Treatment of diseased and failed joints should therefore not only be restricted to anatomical reconstruction of static limb axes alignment. The dynamic activation of muscles, as a key modifier of lower limb biomechanics, should also be taken into account and thus also represents a promising target for restoring function, patient mobility, and preventing future joint failure.

Trial registration

German Clinical Trials Register: ID: DRKS00000606, date: 05.11.2010.

Electronic supplementary material

The online version of this article (10.1186/s12984-018-0434-3) contains supplementary material, which is available to authorized users.

Building a Bird: Musculoskeletal Modeling and Simulation of Wing-Assisted Incline Running During Avian Ontogeny

Flapping flight is the most power-demanding mode of locomotion, associated with a suite of anatomical specializations in extant adult birds. In contrast, many developing birds use their forelimbs to negotiate environments long before acquiring "flight adaptations," recruiting their developing wings to continuously enhance leg performance and, in some cases, fly. How does anatomical development influence these locomotor behaviors? Isolating morphological contributions to wing performance is extremely challenging using purely empirical approaches. However, musculoskeletal modeling and simulation techniques can incorporate empirical data to explicitly examine the functional consequences of changing morphology by manipulating anatomical parameters individually and estimating their effects on locomotion. To assess how ontogenetic changes in anatomy affect locomotor capacity, we combined existing empirical data on muscle morphology, skeletal kinematics, and aerodynamic force production with advanced biomechanical modeling and simulation techniques to analyze the ontogeny of pectoral limb function in a precocial ground bird (Alectoris chukar). Simulations of wing-assisted incline running (WAIR) using these newly developed musculoskeletal models collectively suggest that immature birds have excess muscle capacity and are limited more by feather morphology, possibly because feathers grow more quickly and have a different style of growth than bones and muscles. These results provide critical information about the ontogeny and evolution of avian locomotion by (i) establishing how muscular and aerodynamic forces interface with the skeletal system to generate movement in morphing juvenile birds, and (ii) providing a benchmark to inform biomechanical modeling and simulation of other locomotor behaviors, both across extant species and among extinct theropod dinosaurs.

A subject-specific integrative biomechanical framework of the pelvis for gait analysis

Analysis of the human locomotor system using rigid-body musculoskeletal models has increased in the biomechanical community with the objective of studying muscle activations of different movements. Simultaneously, the finite element method has emerged as a complementary approach for analyzing the mechanical behavior of tissues. This study presents an integrative biomechanical framework for gait analysis by linking a musculoskeletal model and a subject-specific finite element model of the pelvis. To investigate its performance, a convergence study was performed and its sensitivity to the use of non-subject-specific material properties was studied. The total hip joint force estimated by the rigid musculoskeletal model and by the finite element model showed good agreement, suggesting that the integrative approach estimates adequately (in shape and magnitude) the hip total contact force. Previous studies found movements of up to 1.4 mm in the anterior-posterior direction, for single leg stance. These results are comparable with the displacement values found in this study: 0-0.5 mm in the sagittal axis. Maximum von Mises stress values of approximately 17 MPa were found in the pelvic bone. Comparing this results with a previous study of our group, the new findings show that the introduction of muscular boundary conditions and the flexion-extension movement of the hip reduce the regions of high stress and distributes more uniformly the stress across the pelvic bone. Thus, it is thought that muscle force has a relevant impact in reducing stresses in pelvic bone during walking of the finite element model proposed in this study. Future work will focus on including other deformable structures, such as the femur and the tibia, and subject-specific material properties.

Exploring Validation and Verification: How they Different and What They Mean to Healthcare Simulation

STATEMENT

The healthcare simulation (HCS) community recognizes the importance of quality management because many novel simulation devices and techniques include some sort of description of how they tested and assured their simulation's quality. Verification and validation play a key role in quality management; however, literature published on HCS has many different interpretations of what these terms mean and how to accomplish them. The varied use of these terms leads to varied interpretations of how verification process is different from validation process. We set out to explore the concepts of verification and validation in this article by reviewing current psychometric science description of the concepts and exploring how other communities relevant to HCS, such as medical device manufacturing, aviation simulation, and the fields of software and engineering, which are building blocks of technology-enhanced HCS, use the terms, with the focus of trying to clarify the process of verification. We also review current literature available on verification, as compared with validation in HCS and, finally, offer a working definition and concept for each of these terms with hopes to facilitate improved communication within, and with colleagues outside, the HCS community.

Construction and evaluation of a model for wheelchair propulsion in an individual with tetraplegia

Upper limb overuse injuries are common in manual wheelchair users with spinal cord injury. Patient-specific in silico models enhance experimental biomechanical analyses by estimating in vivo shoulder muscle and joint contact forces. Current models exclude deep shoulder muscles that have important roles in wheelchair propulsion. Freely accessible patient-specific models have not been generated for persons with tetraplegia, who have a greater risk for shoulder pain and injury. The objectives of this work were to (i) construct a freely accessible, in silico, musculoskeletal model capable of generating patient-specific dynamic simulations of wheelchair propulsion and (ii) establish proof-of-concept with data obtained from an individual with tetraplegia. Constructed with OpenSim, the model features muscles excluded in existing models. Shoulder muscle forces and activations were estimated via inverse dynamics. Mean absolute error of estimated muscle activations and fine-wire electromyography (EMG) recordings was computed. Mean muscle activation for five consecutive stroke cycles demonstrated good correlation (0.15-0.17) with fine-wire EMG. These findings, comparable to other studies, suggest that the model is capable of estimating shoulder muscle forces during wheelchair propulsion. The additional muscles may provide a greater understanding of shoulder muscle contribution to wheelchair propulsion. The model may ultimately serve as a powerful clinical tool. Graphical abstract ᅟ.

A Monte Carlo analysis of muscle force estimation sensitivity to muscle-tendon properties using a Hill-based muscle model

Surface electromyography driven models are desirable for estimating subject-specific muscle forces. However, these models include parameters that come from an array of sources, thus creating uncertainty in the model-estimated force. In this study, we used Monte-Carlo simulations to evaluate the sensitivity of Hill-based model muscle forces to changes in 11 parameters in the muscle-tendon unit morphological properties and in the model force-length and force-velocity relationships. We decomposed the force variability and ranked the sensitivity of the model to the underlying parameters using the Variogram Analysis of Response Surfaces. For the analyzed running experiments and the adopted Hill model structure, our results show that the parameters are separable into four groups, where the parameters in each group have a synergistic contribution to the model global sensitivity. The first group consists of the maximum isometric force and the pennation angle. The second group contains the optimal fiber length, the tendon slack length, the tendon reference strain and the tendon shape factor. The third group contains the width and shape at the extremities of the active contractile element, along with the maximum contraction velocity and the curvature constant in the force-velocity curve. The fourth group consisted only of the force enhancement during eccentric contraction. The first two groups revealed the largest influence on the output force sensitivity. As many input parameters are difficult to measure and impact estimated forces, we propose that model estimates be presented with confidence intervals as well as inter-parameter relationships, to encourage users to explicitly consider the model uncertainty.

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

The purpose of the current study was to investigate the robustness of dynamic simulation results in the presence of uncertainties resulting from application of a scaled-generic musculoskeletal model instead of a subject-specific model as well as the effect of the choice of simulation method on the obtained muscle forces. The performed sensitivity analysis consisted of the following multibody parameter modifications: maximum isometric muscle forces, number of muscles, the hip joint centre location, segment masses as well as different dynamic simulation methods, namely static optimization with three different criteria and a computed muscle control algorithm (hybrid approach combining forward and inverse dynamics). Twenty-four different models and fifty-five resultant dynamic simulation data sets were analysed. The effects of model perturbation on the magnitude and profile of muscle forces were compared. It has been shown that estimated muscle forces are very sensitive to model parameters. The greatest impact was observed in the case of the force magnitude of the muscles generating high forces during gait (regardless of the modification introduced). However, the force profiles of those muscles were preserved. Relatively large differences in muscle forces were observed for different simulation techniques, which included both magnitude and profile of muscle forces. Personalization of model parameters would affect the resultant muscle forces and seems to be necessary to improve general accuracy of the estimated parameters. However, personalization alone will not ensure high accuracy due to the still unresolved muscle force sharing problem.

Estimation of muscle activation during different walking speeds with two mathematical approaches compared to surface EMG

BACKGROUND

Muscle force estimation could improve clinical gait analysis by enhancing insight into causes of impairments and informing targeted treatments. However, it is not currently standard practice to use muscle force models to augment clinical gait analysis, partly, because robust validations of estimated muscle activations, underpinning force modelling processes, against recorded electromyography (EMG) are lacking.

RESEARCH QUESTION

Therefore, in order to facilitate future clinical use, this study sought to validate estimated lower limb muscle activation using two mathematical models (static optimisation SO, computed muscle control CMC) against recorded muscle activations of ten healthy participants.

METHODS

Participants walked at five speeds. Visual agreement in activation onset and offset as well as linear correlation (r) and mean absolute error (MAE) between models and EMG were evaluated.

RESULTS

MAE between measured and recorded activations were variable across speeds (SO vs EMG 15-68%, CMC vs EMG 13-69%). Slower speeds resulted in smaller deviations (mean MAE < 30%) than faster speeds. Correlation was high (r > 0.5) for only 11/40 (CMC) and 6/40 (SO) conditions (muscles X speeds) compared to EMG.

SIGNIFICANCE

Modelling approaches do not yet show sufficient consistency of agreement between estimated and recorded muscle activation to support recommending immediate clinical adoption of muscle force modelling. This may be because assumptions underlying muscle activation estimations (e.g. muscles' anatomy and maximum voluntary contraction) are not yet sufficiently individualizable. Future research needs to find timely and cost efficient ways to scale musculoskeletal models for better individualisation to facilitate future clinical implementation.

The inclusion of hyoid muscles improve moment generating capacity and dynamic simulations in musculoskeletal models of the head and neck

OpenSim musculoskeletal models of the head and neck can provide information about muscle activity and the response of the head and neck to a variety of situations. Previous models report weak flexion strength, which is partially due to lacking moment generating capacity in the upper cervical spine. Previous models have also lacked realistic hyoid muscles, which have the capability to improve flexion strength and control in the upper cervical spine. Suprahyoid and infrahyoid muscles were incorporated in an OpenSim musculoskeletal model of the head and neck. This model was based on previous OpenSim models, and now includes hyoid muscles and passive elements. The moment generating capacity of the model was tested by simulating physical experiments in the OpenSim environment. The flexor and extensor muscle strengths were scaled to match static experimental results. Models with and without hyoid muscles were used to simulate experimentally captured motions, and the need for reserve actuators was evaluated. The addition of hyoid muscles greatly increased flexion strength, and the model is the first of its kind to have realistic strength values in all directions. Less reserve actuator moment was required to simulate real motions with the addition of hyoid muscles. Several additional ways of improving flexion strength were investigated. Hyoid muscles add control and strength to OpenSim musculoskeletal models of the head and neck and improve simulations of head and neck movements.

A Dynamic Simulation of Musculoskeletal Function in the Mouse Hindlimb During Trotting Locomotion

Mice are often used as animal models of various human neuromuscular diseases, and analysis of these models often requires detailed gait analysis. However, little is known of the dynamics of the mouse musculoskeletal system during locomotion. In this study, we used computer optimization procedures to create a simulation of trotting in a mouse, using a previously developed mouse hindlimb musculoskeletal model in conjunction with new experimental data, allowing muscle forces, activation patterns, and levels of mechanical work to be estimated. Analyzing musculotendon unit (MTU) mechanical work throughout the stride allowed a deeper understanding of their respective functions, with the rectus femoris MTU dominating the generation of positive and negative mechanical work during the swing and stance phases. This analysis also tested previous functional inferences of the mouse hindlimb made from anatomical data alone, such as the existence of a proximo-distal gradient of muscle function, thought to reflect adaptations for energy-efficient locomotion. The results do not strongly support the presence of this gradient within the mouse musculoskeletal system, particularly given relatively high negative net work output from the ankle plantarflexor MTUs, although more detailed simulations could test this further. This modeling analysis lays a foundation for future studies of the control of vertebrate movement through the development of neuromechanical simulations.

Patient-specific musculoskeletal modeling of the hip joint for preoperative planning of total hip arthroplasty: A validation study based on in vivo measurements

Validation of musculoskeletal models for application in preoperative planning is still a challenging task. Ideally, the simulation results of a patient-specific musculoskeletal model are compared to corresponding in vivo measurements. Currently, the only possibility to measure in vivo joint forces is to implant an instrumented prosthesis in patients undergoing a total joint replacement. In this study, a musculoskeletal model of the AnyBody Modeling System was adapted patient-specifically and validated against the in vivo hip joint force measurements of ten subjects performing one-leg stance and level walking. The impact of four model parameters was evaluated; hip joint width, muscle strength, muscle recruitment, and type of muscle model. The smallest difference between simulated and in vivo hip joint force was achieved by using the hip joint width measured in computed tomography images, a muscle strength of 90 N/cm², a third order polynomial muscle recruitment, and a simple muscle model. This parameter combination reached mean deviations between simulation and in vivo measurement during the peak force phase of 12% ± 14% in magnitude and 11° ± 5° in orientation for one-leg stance and 8% ± 6% in magnitude and 10° ± 5° in orientation for level walking.

Assessing rugby place kick performance from initial ball flight kinematics: development, validation and application of a new measure

The appropriate determination of performance outcome is critical when appraising a performer's technique. Previous studies of rugby place kicking technique have typically assessed performance based on ball velocity, but this is not the sole requirement. Therefore, a mathematical model of rugby place kick ball flight was developed to yield a single measure more representative of true performance. The model, which requires only initial ball flight kinematics, was calibrated and validated using empirical place kick data, and found to predict ball position with a mean error of 4.0% after 22 m of ball flight. The model was then applied to the performances of 33 place kickers. The predicted maximum distance, a single performance measure which accounted for initial ball velocity magnitude and direction, and spin, was determined using the model and was compared against ball velocity magnitude. A moderate association in the rank-order of the kicks between these two measures (ρ = 0.52) revealed that the relative success of the kicks would be assessed differently with each measure. The developed model provides a representative measure of place kick performance that is understandable for coaches, and can be used to predict changes in performance outcome under different ball launch or environmental conditions.