Accuracy and Reliability of Marker-Based Approaches to Scale the Pelvis, Thigh, and Shank Segments in Musculoskeletal Models.

Gait analysis together with musculoskeletal modeling is widely used for research. In the absence of medical images, surface marker locations are used to scale a generic model to the individual's anthropometry. Studies evaluating the accuracy and reliability of different scaling approaches in a pediatric and/or clinical population have not yet been conducted and, therefore, formed the aim of this study. Magnetic resonance images (MRI) and motion capture data were collected from 12 participants with cerebral palsy and 6 typically developed participants. Accuracy was assessed by comparing the scaled model's segment measures to the corresponding MRI measures, whereas reliability was assessed by comparing the model's segments scaled with the experimental marker locations from the first and second motion capture session. The inclusion of joint centers into the scaling process significantly increased the accuracy of thigh and shank segment length estimates compared to scaling with markers alone. Pelvis scaling approaches which included the pelvis depth measure led to the highest errors compared to the MRI measures. Reliability was similar between scaling approaches with mean ICC of 0.97. The pelvis should be scaled using pelvic width and height and the thigh and shank segment should be scaled using the proximal and distal joint centers.

Calcium phosphate cement augmentation after volar locking plating of distal radius fracture significantly increases stability.

PURPOSE: Distal radius fractures represent the most common fractures in adults. Volar locking plating to correct unstable fractures has become increasingly popular. Although reasonable primary reduction is possible in most cases, maintenance of reduction until the fracture is healed is often problematic in osteoporotic bone. To our knowledge, no biomechanical studies have compared the effect of enhancement with biomaterial on two different volar fixed-angle plates. METHODS: Human fresh-frozen cadaver pairs of radii were used to simulate an AO/OTA 23-A3 fracture. In a total of four groups (n = 7 for each group), two volar fixed-angle plates (Aptus 2.5 mm locking fracture plate, Medartis, Switzerland and VA-LCP two-column distal radius plate 2.4, volar, Synthes, Switzerland) with or without an additional injection of a biomaterial (Hydroset Injectable HA Bone Substitute, Stryker, Switzerland) into the dorsal comminution zone were used to fix the distal metaphyseal fragment. Each specimen was tested load-controlled under cyclic loading with a servo-hydraulic material testing machine. Displacement, stiffness, dissipated work and failure mode were recorded. RESULTS: Improved mechanical properties (decreased displacement, increased stiffness, decreased dissipated work) were found in both plates if the biomaterial was additionally injected. Improvement of mechanical parameters after biomaterial injection was more evident in the Synthes plate compared to the Aptus plate. Pushing out of the screws was noticed as a failure mode only in samples lacking supplementary biomaterial. CONCLUSIONS: Injection of a biomaterial into the dorsal comminution zone increases stability after volar locking plating of distal radius fractures in vitro.

Repeatability and minimal detectable change including clothing effects for smartphone-based 3D markerless motion capture.

OpenCap, a smartphone- and web-based markerless system, has shown acceptable accuracy compared to marker-based systems, but lacks information on repeatability. This study fills this gap by evaluating the intersession repeatability of OpenCap and investigating the effects of clothing on gait kinematics. Twenty healthy volunteers participated in a test-retest study, performing walking and sit-to-stand tasks with minimal clothing and regular street wear. Segment lengths and lower-limb kinematics were compared between both sessions and for both clothing conditions using the root-mean-square-deviation (RMSD) for entire waveforms and the standard error of measurement (SEM) and minimal detectable change (MDC) for discrete kinematic parameters. In general, the RMSD test-retest values were 2.8 degrees (SD: 1.0) for walking and 3.3 degrees (SD: 1.2) for sit-to-stand. The highest intersession variability was observed in the trunk, pelvis, and hip kinematics of the sagittal plane. SEM and MDC values were on average 2.2 and 6.0 degrees, respectively, for walking, and 2.4 and 6.5 degrees for sit-to-stand. Clothing had minimal effects on kinematics by adding on average less than one degree to the RMSD values for most variables. The segment lengths showed moderate to excellent agreement between both sessions and poor to moderate agreement between clothing conditions. The study highlights the reliability of OpenCap for markerless motion capture, emphasizing its potential for large-scale field studies. However, some variables showed high MDC values above 5 degrees and thus warrant further enhancement of the technology. Although clothing had minimal effects, it is still recommended to maintain consistent clothing to minimize overall variability.

Muscle synergies are shared across fundamental subtasks in complex movements of skateboarding.

A common theory of motor control posits that movement is controlled by muscle synergies. However, the behavior of these synergies during highly complex movements remains largely unexplored. Skateboarding is a hardly researched sport that requires rapid motor control to perform tricks. The objectives of this study were to investigate three key areas: (i) whether motor complexity differs between skateboard tricks, (ii) the inter-participant variability in synergies, and (iii) whether synergies are shared between different tricks. Electromyography data from eight muscles per leg were collected from seven experienced skateboarders performing three different tricks (Ollie, Kickflip, 360°-flip). Synergies were extracted using non-negative matrix factorization. The number of synergies (NoS) was determined using two criteria based on the total variance accounted for (tVAF > 90% and adding an additional synergy does not increase tVAF > 1%). In summary: (i) NoS and tVAF did not significantly differ between tricks, indicating similar motor complexity. (ii) High inter-participant variability exists across participants, potentially caused by the low number of constraints given to perform the tricks. (iii) Shared synergies were observed in every comparison of two tricks. Furthermore, each participant exhibited at least one synergy vector, which corresponds to the fundamental 'jumping' task, that was shared through all three tricks.

Increased trial-to-trial similarity and reduced temporal overlap of muscle synergy activation coefficients manifest during learning and with increasing movement proficiency.

Muscle synergy analyses are used to enhance our understanding of motor control. Spatially fixed synergy weights coordinate multiple co-active muscles through activation commands, known as activation coefficients. To gain a more comprehensive understanding of motor learning, it is essential to understand how activation coefficients vary during a learning task and at different levels of movement proficiency. Participants walked on a line, a beam, and learned to walk on a tightrope-tasks that represent different levels of proficiency. Muscle synergies were extracted from electromyography signals across all conditions and the number of synergies was determined by the knee-point of the total variance accounted for (tVAF) curve. The results indicated that the tVAF of one synergy decreased with task proficiency, with the tightrope task resulting in the highest tVAF compared to the line and beam tasks. Furthermore, with increasing proficiency and after a learning process, trial-to-trial similarity increased and temporal overlap of synergy activation coefficients decreased. Consequently, we propose that precise adjustment and refinement of synergy activation coefficients play a pivotal role in motor learning.

Demonstrating the utility of Instrumented Gait Analysis in the treatment of children with cerebral palsy.

BACKGROUND: Instrumented gait analysis (IGA) has been around for a long time but has never been shown to be useful for improving patient outcomes. In this study we demonstrate the potential utility of IGA by showing that machine learning models are better able to estimate treatment outcomes when they include both IGA and clinical (CLI) features compared to when they include CLI features alone. DESIGN: We carried out a retrospective analysis of data from ambulatory children diagnosed with cerebral palsy who were seen at least twice at our gait analysis center. Individuals underwent a variety of treatments (including no treatment) between sequential gait analyses. We fit Bayesian Additive Regression Tree (BART) models that estimated outcomes for mean stance foot progression to demonstrate the approach. We built two models: one using CLI features only, and one using CLI and IGA features. We then compared the models' performance in detail. We performed similar, but less detailed, analyses for a number of other outcomes. All results were based on independent test data from a 70%/30% training/testing split. RESULTS: The IGA model was more accurate than the CLI model for mean stance-phase foot progression outcomes (RMSEIGA = 11∘, RMSECLI = 13∘) and explained more than 1.5 × as much of the variance (R2IGA = .45, R2CLI = .28). The IGA model outperformed the CLI model for every level of treatment complexity, as measured by number of simultaneous surgeries. The IGA model also exhibited superior performance for estimating outcomes of mean stance-phase knee flexion, mean stance-phase ankle dorsiflexion, maximum swing-phase knee flexion, gait deviation index (GDI), and dimensionless speed. INTERPRETATION: The results show that IGA has the potential to be useful in the treatment planning process for ambulatory children diagnosed with cerebral palsy. We propose that the results of machine learning outcome estimators-including estimates of uncertainty-become the primary IGA tool utilized in the clinical process, complementing the standard medical practice of conducting a through patient history and physical exam, eliciting patient goals, reviewing relevant imaging data, and so on.

ESMAC Best Paper Award 2023: Increased knee flexion in participants with cerebral palsy results in altered stresses at the distal femoral growth plate compared to a typically developing cohort.

INTRODUCTION: Femoral deformities are highly prevalent in children with cerebral palsy (CP) and can have a severe impact on patients' gait abilities. While the mechanical stress regime within the distal femoral growth plate remains underexplored, understanding it is crucial given bone's adaptive response to mechanical stimuli. We quantified stresses at the distal femoral growth plate to deepen our understanding of the relationship between healthy and pathological gait patterns, internal loading, and femoral growth patterns. METHODS: This study included three-dimensional motion capture data and magnetic resonance images of 13 typically developing children and twelve participants with cerebral palsy. Employing a multi-scale mechanobiological approach, integrating musculoskeletal simulations and subject-specific finite element analysis, we investigated the orientation of the distal femoral growth plate and the stresses within it. Limbs of participants with CP were grouped depending on their knee flexion kinematics during stance phase as this potentially changes the stresses induced by knee and patellofemoral joint contact forces. RESULTS: Despite similar growth plate orientation across groups, significant differences were observed in the shape and distribution of growth values. Higher growth rates were noted in the anterior compartment in CP limbs with high knee flexion while CP limbs with normal knee flexion showed high similarity to the group of healthy participants. DISCUSSION: Results indicate that the knee flexion angle during the stance phase is of high relevance for typical bone growth at the distal femur. The evaluated growth rates reveal plausible results, as long-term promoted growth in the anterior compartment leads to anterior bending of the femur which was confirmed for the group with high knee flexion through analyses of the femoral geometry. The framework for these multi-scale simulations has been made accessible on GitHub, empowering peers to conduct similar mechanobiological studies. Advancing our understanding of femoral bone development could ultimately support clinical decision-making.

Efficacy of Real-Time Feedback Exercise Therapy in Patients Following Total Hip Arthroplasty: Protocol for a Pilot Cluster-Randomized Controlled Trial.

BACKGROUND: Osteoarthritis of the hip joint is an increasing functional and health-related problem. The most common surgical treatment is hip replacement to reduce pain and improve function. Rehabilitation after total hip arthroplasty (THA) is not regulated in Austria and mostly depends on the patient's own initiative and possibilities. Functional deficits, such as valgus thrust of the leg, functional Trendelenburg gait, or Duchenne limp, are characteristic symptoms before and, due to the performance learning effect prior to surgery, also after the operation. Addressing these deficits is possible through neuromuscular-focused exercise therapy. The efficacy of such therapy relies significantly on the quality of performance, the frequency of exercise, and the duration of engagement. Enhancing sustainability is achievable through increased motivation and real-time feedback (RTF) on exercise execution facilitated by digital feedback systems. OBJECTIVE: This study will be performed to quantify the medium-term effectiveness of digital home exercise feedback systems on functional performance following THA. METHODS: A clinical trial with a cluster-randomized, 2-arm, parallel-group design with an 8-week intervention phase and subsequent follow-ups at 3 and 6 months postsurgery will be conducted. Feedback during exercising will be provided through a blended-care program, combining a supervised group exercise program with a self-developed digital feedback system for home exercise. In total, 70 patients will be recruited for baseline. The primary outcome parameters will be the frontal knee range of motion, pelvic obliquity, and lateral trunk lean. Secondary outcomes will be the sum scores of patient-reported outcomes and relevant kinematic, kinetic, and spatiotemporal parameters. RESULTS: The trial started in January 2024, and the first results are anticipated to be published by June 2025. RTF-supported home exercise is expected to improve exercise execution quality and therapeutic adherence compared to using paper instructions for excise guidance. CONCLUSIONS: The anticipated findings of this study aim to offer new insights into the effect of a blended-care program incorporating digital RTF on exercise therapy after unilateral THA, in addition to knowledge on the functional status 3 and 6 months postsurgery, for further improvement in the development of rehabilitation guidelines following THA. TRIAL REGISTRATION: ClinicalTrials.gov: NCT06161194; https://clinicaltrials.gov/study/NCT06161194. INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID): PRR1-10.2196/59755.

A framework based on subject-specific musculoskeletal models and Monte Carlo simulations to personalize muscle coordination retraining.

Excessive loads at lower limb joints can lead to pain and degenerative diseases. Altering joint loads with muscle coordination retraining might help to treat or prevent clinical symptoms in a non-invasive way. Knowing how much muscle coordination retraining can reduce joint loads and which muscles have the biggest impact on joint loads is crucial for personalized gait retraining. We introduced a simulation framework to quantify the potential of muscle coordination retraining to reduce joint loads for an individuum. Furthermore, the proposed framework enables to pinpoint muscles, which alterations have the highest likelihood to reduce joint loads. Simulations were performed based on three-dimensional motion capture data of five healthy adolescents (femoral torsion 10°-29°, tibial torsion 19°-38°) and five patients with idiopathic torsional deformities at the femur and/or tibia (femoral torsion 18°-52°, tibial torsion 3°-50°). For each participant, a musculoskeletal model was modified to match the femoral and tibial geometry obtained from magnetic resonance images. Each participant's model and the corresponding motion capture data were used as input for a Monte Carlo analysis to investigate how different muscle coordination strategies influence joint loads. OpenSim was used to run 10,000 simulations for each participant. Root-mean-square of muscle forces and peak joint contact forces were compared between simulations. Depending on the participant, altering muscle coordination led to a maximum reduction in hip, knee, patellofemoral and ankle joint loads between 5 and 18%, 4% and 45%, 16% and 36%, and 2% and 6%, respectively. In some but not all participants reducing joint loads at one joint increased joint loads at other joints. The required alteration in muscle forces to achieve a reduction in joint loads showed a large variability between participants. The potential of muscle coordination retraining to reduce joint loads depends on the person's musculoskeletal geometry and gait pattern and therefore showed a large variability between participants, which highlights the usefulness and importance of the proposed framework to personalize gait retraining.

The influence of patellar resurfacing on patellar kinetics and retropatellar contact characteristics.

BACKGROUND: Femoropatellar complications are one of the most common problems after total knee arthroplasty (TKA). However, the question of whether to resurface the patella remains controversial. Therefore, we evaluated the kinetics and the retropatellar contact characteristics of patella resurfacing with fixed and gliding surfaces. METHODS: Eight Thiel-embalmed cadaver knees were tested--first intact, then after TKA without patellar resurfacing, and finally with additional patellar resurfacing--while flexing the knee from 0° to 100°. We tested a fixed as well as a gliding patella surface. During the examination, quadriceps and hamstring forces were applied. The retropatellar pressure was determined with a special patella sensor, and the patellar kinetics were measured using an optical three-dimensional motion analysis system. RESULTS: Resurfacing the patella caused a significant increase in retropatellar pressure and a significant decrease in retropatellar contact area. Using a fixed patella, the retropatellar pressure nearly quadrupled in higher flexion compared to the native patella. Furthermore, the lateral movement of the patella increased after TKA, especially after additional patellar resurfacing. CONCLUSIONS: Resurfacing the patella routinely is not advised. When osteoarthritis of the patella is found, the gliding patella should be preferred.

Imaging-based musculoskeletal models alter muscle and joint contact forces but do not improve the agreement with experimentally measured electromyography signals in children with cerebral palsy.

BACKGROUND: Musculoskeletal simulations are used to estimate muscle-tendon and joint contact forces (JCF). Personalizing the model's femoral geometry has been shown to improve the accuracy of JCF calculations. It is, however, unknown if the personalized geometry improves the agreement between estimated muscle activations and experimentally measured electromyography (EMG) signals. RESEARCH QUESTION: Does personalizing the musculoskeletal geometry improve the agreement between estimated muscle activations and EMG signals in terms of timing? METHODS: We retrospectively analysed data from Bosmans et al. [5], which included three-dimensional motion capture data, EMG signals of eight lower limb muscles on each leg, and magnetic resonance imaging (MRI) data from seven children with cerebral palsy. For each patient we created a generic-scaled model and MRI-based model, which accounted for the subject-specific musculoskeletal geometry. We calculated muscle activations, muscle-tendon forces and JCF. Muscle activations were compared to the EMG signals using coefficient of determination and cosines similarity. RESULTS: MRI-based models altered the magnitude of muscle activations and had a large impact on JCF but did not change the muscle activations profiles and therefore did not improve the agreement with EMG signals. SIGNIFICANCE: MRI-based models do not alter the shape of muscle activations. Hence, if detailed muscle activations are a desired output of the simulations, EMG-informed modeling approaches should be used for musculoskeletal simulation in children with cerebral palsy. Furthermore, our study highlighted that altered JCF does not necessarily mean accurate muscle activations. To improve patient-specific simulations, future work should focus on developing methods to estimate cost functions representative for the neural control of children with cerebral palsy.

The gait pattern and not the femoral morphology is the main contributor to asymmetric hip joint loading.

Gait asymmetry and skeletal deformities are common in many children with cerebral palsy (CP). Changes of the hip joint loading, i.e. hip joint contact force (HJCF), can lead to pathological femoral growth. A child's gait pattern and femoral morphology affect HJCFs. The twofold aim of this study was to (1) evaluate if the asymmetry in HJCFs is higher in children with CP compared to typically developing (TD) children and (2) identify if the bony morphology or the subject-specific gait pattern is the main contributor to asymmetric HJCFs. Magnetic resonance images (MRI) and three-dimensional gait analysis data of twelve children with CP and fifteen TD children were used to create subject-specific musculoskeletal models and calculate HJCF using OpenSim. Root-mean-square-differences between left and right HJCF magnitude and orientation were computed and compared between participant groups (CP versus TD). Additionally, the influence on HJCF asymmetries solely due to the femoral morphology and solely due to the gait pattern was quantified. Our findings demonstrate that the gait pattern is the main contributor to asymmetric HJCFs in CP and TD children. Children with CP have higher HJCF asymmetries which is probably the result of larger asymmetries in their gait pattern compared to TD children. The gained insights from our study highlight that clinical interventions should focus on normalizing the gait pattern and therefore the hip joint loading to avoid the development of femoral deformities.

Influence of femoral anteversion angle and neck-shaft angle on muscle forces and joint loading during walking.

Femoral deformities, e.g. increased or decreased femoral anteversion (AVA) and neck-shaft angle (NSA), can lead to pathological gait patterns, altered joint loads, and degenerative joint diseases. The mechanism how femoral geometry influences muscle forces and joint load during walking is still not fully understood. The objective of our study was to investigate the influence of femoral AVA and NSA on muscle forces and joint loads during walking. We conducted a comprehensive musculoskeletal modelling study based on three-dimensional motion capture data of a healthy person with a typical gait pattern. We created 25 musculoskeletal models with a variety of NSA (93°-153°) and AVA (-12°-48°). For each model we calculated moment arms, muscle forces, muscle moments, co-contraction indices and joint loads using OpenSim. Multiple regression analyses were used to predict muscle activations, muscle moments, co-contraction indices, and joint contact forces based on the femoral geometry. We found a significant increase in co-contraction of hip and knee joint spanning muscles in models with increasing AVA and NSA, which led to a substantial increase in hip and knee joint contact forces. Decreased AVA and NSA had a minor impact on muscle and joint contact forces. Large AVA lead to increases in both knee and hip contact forces. Large NSA (153°) combined with large AVA (48°) led to increases in hip joint contact forces by five times body weight. Low NSA (108° and 93°) combined with large AVA (48°) led to two-fold increases in the second peak of the knee contact forces. Increased joint contact forces in models with increased AVA and NSA were linked to changes in hip muscle moment arms and compensatory increases in hip and knee muscle forces. Knowing the influence of femoral geometry on muscle forces and joint loads can help clinicians to improve treatment strategies in patients with femoral deformities.

Intra- and inter-subject variability of femoral growth plate stresses in typically developing children and children with cerebral palsy.

Little is known about the influence of mechanical loading on growth plate stresses and femoral growth. A multi-scale workflow based on musculoskeletal simulations and mechanobiological finite element (FE) analysis can be used to estimate growth plate loading and femoral growth trends. Personalizing the model in this workflow is time-consuming and therefore previous studies included small sample sizes (N < 4) or generic finite element models. The aim of this study was to develop a semi-automated toolbox to perform this workflow and to quantify intra-subject variability in growth plate stresses in 13 typically developing (TD) children and 12 children with cerebral palsy (CP). Additionally, we investigated the influence of the musculoskeletal model and the chosen material properties on the simulation results. Intra-subject variability in growth plate stresses was higher in cerebral palsy than in typically developing children. The highest osteogenic index (OI) was observed in the posterior region in 62% of the TD femurs while in children with CP the lateral region was the most common (50%). A representative reference osteogenic index distribution heatmap generated from data of 26 TD children's femurs showed a ring shape with low values in the center region and high values at the border of the growth plate. Our simulation results can be used as reference values for further investigations. Furthermore, the code of the developed GP-Tool ("Growth Prediction-Tool") is freely available on GitHub (https://github.com/WilliKoller/GP-Tool) to enable peers to conduct mechanobiological growth studies with larger sample sizes to improve our understanding of femoral growth and to support clinical decision making in the near future.

Neuromechanics of finger hangs with arm lock-offs: analyzing joint moments and muscle activations to improve practice guidelines for climbing.

Introduction: Climbing imposes substantial demands on the upper limbs and understanding the mechanical loads experienced by the joints during climbing movements is crucial for injury prevention and optimizing training protocols. This study aimed to quantify and compare upper limb joint loads and muscle activations during isometric finger hanging exercises with different arm lock-off positions. Methods: Seventeen recreational climbers performed six finger dead hangs with arm lock-offs at 90° and 135° of elbow flexion, as well as arms fully extended. Upper limb joint moments were calculated using personalized models in OpenSim, based on three-dimensional motion capture data and forces measured on an instrumented hang board. Muscle activations of upper limb muscles were recorded with surface electromyography electrodes. Results: Results revealed that the shoulder exhibited higher flexion moments during arm lock-offs at 90° compared to full extension (p = 0.006). The adduction moment was higher at 135° and 90° compared to full extension (p < 0.001), as well as the rotation moments (p < 0.001). The elbows exhibited increasing flexion moments with the increase in the arm lock-off angle (p < 0.001). Muscle activations varied across conditions for biceps brachii (p < 0.001), trapezius (p < 0.001), and latissimus dorsi, except for the finger flexors (p = 0.15). Discussion: Our findings indicate that isometric finger dead hangs with arms fully extended are effective for training forearm force capacities while minimizing stress on the elbow and shoulder joints. These findings have important implications for injury prevention and optimizing training strategies in climbing.

Influence of total knee arthroplasty on patellar kinematics and contact characteristics.

PURPOSE: Femoro-patellar complications are one of the most common problems after total knee arthroplasty (TKA). TKA components that reduce patellar loads and preserve physiological patellar kinematics should reduce these problems. Therefore, we evaluated the patellar kinematics and the retro-patellar contact characteristics in both the intact knee and in the TKA-knee. METHODS: Eight Thiel-embalmed cadaver knees were tested first intact and then after TKA using rotating as well as gliding inlay and with additional patellar resurfacing while flexing the knee from 0° to 100°. During the examination quadriceps and hamstring forces were applied. RESULTS: TKA with additional patellar resurfacing led to an increased retro-patellar pressure, a decreased contact area and an increased lateral movement. Although patellar kinematics could not be changed by using a gliding inlay compared to a rotating inlay, the gliding inlay improved retro-patellar contact characteristics by reducing the pressure and increasing the contact area, especially in higher flexion. CONCLUSIONS: The increased retro-patellar pressure together with the increased lateral movement of the patella after TKA may be one important cause for anterior knee pain appearing after TKA. In view of the improved retro-patellar contact characteristics using a gliding inlay this inlay should be preferred, providing that the posterior cruciate ligament is intact.

The importance of a consistent workflow to estimate muscle-tendon lengths based on joint angles from the conventional gait model.

BACKGROUND: Musculoskeletal models enable us to estimate muscle-tendon length, which has been shown to improve clinical decision-making and outcomes in children with cerebral palsy. Most clinical gait analysis services, however, do not include muscle-tendon length estimation in their clinical routine. This is due, in part, to a lack of knowledge and trust in the musculoskeletal models, and to the complexity involved in the workflow to obtain the muscle-tendon length. RESEARCH QUESTION: Can the joint angles obtained with the conventional gait model (CGM) be used to generate accurate muscle-tendon length estimates? METHODS: Three-dimensional motion capture data of 15 children with cerebral palsy and 15 typically developing children were retrospectively analyzed and used to estimate muscle-tendon length with the following four modelling frameworks: (1) 2392-OSM-IK-angles: standard OpenSim workflow including scaling, inverse kinematics and muscle analysis; (2) 2392-OSM-CGM-angle: generic 2392-OpenSim model driven with joint angles from the CGM; (3) modif-OSM-IK-angles: standard OpenSim workflow including inverse kinematics and a modified model with segment coordinate systems and joint degrees-of-freedom similar to the CGM; (4) modif-OSM-CGM-angles: modified model driven with joint angles from the CGM. Joint kinematics and muscle-tendon length were compared between the different modelling frameworks. RESULTS: Large differences in hip joint kinematics were observed between the CGM and the 2392-OpenSim model. The modif-OSM showed similar kinematics as the CGM. Muscle-tendon length obtained with modif-OSM-IK-angles and modif-OSM-CGM-angles were similar, whereas large differences in some muscle-tendon length were observed between 2392-OSM-IK-angles and 2392-OSM-CGM-angles. SIGNIFICANCE: The modif-OSM-CGM-angles framework enabled us to estimate muscle-tendon lengths without the need for scaling a musculoskeletal model and running inverse kinematics. Hence, muscle-tendon length estimates can be obtained simply, without the need for the complexity, knowledge and time required for musculoskeletal modeling and associated software. An instruction showing how the framework can be used in a clinical setting is provided on https://github.com/HansUniVie/MuscleLength.

Impact of scaling errors of the thigh and shank segments on musculoskeletal simulation results.

BACKGROUND: Musculoskeletal simulations are widely used in the research community. The locations of surface markers are mostly used to scale a generic model to the participant's anthropometry. Marker-based scaling approaches include errors due to inaccuracies in marker placements. RESEARCH QUESTION: How do scaling errors of the thigh and shank segments influence simulation results? METHODS: Motion capture data and magnetic resonance images from a child with cerebral palsy and a typically developing child were used to create a subject-specific reference model for each child. These reference models were modified to mimic scaling errors due to inaccurately placed lateral epicondyle markers, which are frequently used to scale the thigh and shank segments. The thigh length was altered in 1 % steps from the original length and the shank length was accordingly adjusted to keep the total leg length constant. Thirty additional models were created, which included models with an altered thigh length of ±15 %. Subsequently, musculoskeletal simulations with OpenSim were performed with all models. Joint kinematics, joint kinetics, muscle forces and joint contact forces (JCF) were compared between the reference and altered models. RESULTS: The investigated scaling error influenced joint kinematics and joint kinetics by up to 9.4° (hip flexion angle) and 0.15 Nm/kg (knee flexion moment), respectively. Maximum muscle and JCF differences of 46 % (medial gastrocnemius) and 72 % (hip JCF) bodyweight, respectively, were observed between the reference and altered models. Scaling errors mainly changed the magnitude but not the shape of most analyzed parameters. The influence of scaling errors on simulation results were similar in both participants. SIGNIFICANCE: Scaling errors of the thigh segment influence simulation results at all joints due to the global optimization approach used in musculoskeletal simulations. Our findings can be used to estimate potential errors due to marker-based scaling approaches in previous and future studies.

Generic scaled versus subject-specific models for the calculation of musculoskeletal loading in cerebral palsy gait: Effect of personalized musculoskeletal geometry outweighs the effect of personalized neural control.

BACKGROUND: Musculoskeletal modelling is used to assess musculoskeletal loading during gait. Linear scaling methods are used to personalize generic models to each participant's anthropometry. This approach introduces simplifications, especially when used in paediatric and/or pathological populations. This study aimed to compare results from musculoskeletal simulations using various models ranging from linear scaled to highly subject-specific models, i.e., including the participant's musculoskeletal geometry and electromyography data. METHODS: Magnetic resonance images (MRI) and gait data of one typically developing child and three children with cerebral palsy were analysed. Musculoskeletal simulations were performed to calculate joint kinematics, joint kinetics, muscle forces and joint contact forces using four modelling frameworks: 1) Generic-scaled model with static optimization, 2) Generic-scaled model with an electromyography-informed approach, 3) MRI-based model with static optimization, and 4) MRI-based model with an electromyography-informed approach. FINDINGS: Root-mean-square-differences in joint kinematics and kinetics between generic-scaled and MRI-based models were below 5° and 0.15 Nm/kg, respectively. Root-mean-square-differences over all muscles was below 0.2 body weight for every participant. Root-mean-square-differences in joint contact forces between the different modelling frameworks were up to 2.2 body weight. Comparing the simulation results from the typically developing child with the results from the children with cerebral palsy showed similar root-mean-square-differences for all modelling frameworks. INTERPRETATION: In our participants, the impact of MRI-based models on joint contact forces was higher than the impact of including electromyography. Clinical reasoning based on overall root-mean-square-differences in musculoskeletal simulation results between healthy and pathological participants are unlikely to be affected by the modelling choice.

Quantifying Muscle Forces and Joint Loading During Hip Exercises Performed With and Without an Elastic Resistance Band.

An increase in hip joint contact forces (HJCFs) is one of the main contributing mechanical causes of hip joint pathologies, such as hip osteoarthritis, and its progression. The strengthening of the surrounding muscles of the joint is a way to increase joint stability, which results in the reduction of HJCF. Most of the exercise recommendations are based on expert opinions instead of evidence-based facts. This study aimed to quantify muscle forces and joint loading during rehabilitative exercises using an elastic resistance band (ERB). Hip exercise movements of 16 healthy volunteers were recorded using a three-dimensional motion capture system and two force plates. All exercises were performed without and with an ERB and two execution velocities. Hip joint kinematics, kinetics, muscle forces, and HJCF were calculated based on the musculoskeletal simulations in OpenSim. Time-normalized waveforms of the different exercise modalities were compared with each other and with reference values found during walking. The results showed that training with an ERB increases both target muscle forces and HJCF. Furthermore, the ERB reduced the hip joint range of motion during the exercises. The type of ERB used (soft vs. stiff ERB) and the execution velocity of the exercise had a minor impact on the peak muscle forces and HJCF. The velocity of exercise execution, however, had an influence on the total required muscle force. Performing hip exercises without an ERB resulted in similar or lower peak HJCF and lower muscle forces than those found during walking. Adding an ERB during hip exercises increased the peak muscle and HJCF but the values remained below those found during walking. Our workflow and findings can be used in conjunction with future studies to support exercise design.