A multi-scale modelling framework combining musculoskeletal rigid-body simulations with adaptive finite element analyses, to evaluate the impact of femoral geometry on hip joint contact forces and femoral bone growth

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.

Joint contact forces during semi-recumbent seated cycling

Semi-recumbent cycling performed from a wheelchair is a popular rehabilitation exercise following spinal cord injury (SCI) and is often paired with functional electrical stimulation. However, biomechanical assessment of this cycling modality is lacking, even in unimpaired populations, hindering the development of personalised and safe rehabilitation programs for those with SCI. This study developed a computational pipeline to determine lower limb kinematics, kinetics, and joint contact forces (JCF) in 11 unimpaired participants during voluntary semi-recumbent cycling using a rehabilitation ergometer. Two cadences (40 and 60 revolutions per minute) and three crank powers (15 W, 30 W, and 45 W) were assessed. A rigid body model of a rehabilitation ergometer was combined with a calibrated electromyogram-informed neuromusculoskeletal model to determine JCF at the hip, knee, and ankle. Joint excursions remained consistent across all cadence and powers, but joint moments and JCF differed between 40 and 60 revolutions per minute, with peak JCF force significantly greater at 40 compared to 60 revolutions per minute for all crank powers. Poor correlations were found between mean crank power and peak JCF across all joints. This study provides foundation data and computational methods to enable further evaluation and optimisation of semi-recumbent cycling for application in rehabilitation after SCI and other neurological disorders.

Understanding the role of head size and neck length in micromotion generation at the taper junction in total hip arthroplasty

Modular hip implants allow intra-operative adjustments for patient-specific customization and targeted replacement of damaged elements without full implant extraction. However, challenges arise from relative micromotions between components, potentially leading to implant failure due to cytotoxic metal debris. In this study magnitude and directions of micromotions at the taper junction were estimated, aiming to understand the effect of variations in head size and neck length. Starting from a reference configuration adhering to the 12/14 taper standard, six additional implant configurations were generated by varying the head size and/or neck length. A musculoskeletal multibody model of a prothesized lower limb was developed to estimate hip contact force and location during a normal walking task. Following the implant assembly, the multibody-derived loads were imposed as boundary conditions in a finite element analysis to compute the taper junction micromotions as the relative slip between the contacting surfaces. Results highlighted the L-size head as the most critical configuration, indicating a 2.81 μm relative slip at the mid-stance phase. The proposed approach enables the investigation of geometric variations in implants under accurate load conditions, providing valuable insights for designing less risky prostheses and informing clinical decision-making processes.

Combined musculoskeletal finite element modeling of femur stress during reactive balance training

The purpose of this study was to determine the material stresses experienced in the femoral neck during the stepping phase of recovery from a forward loss of balance achieved both using release from a static forward lean and rapid treadmill accelerations in 8 older adults. A scalable musculoskeletal model with 23 degrees of freedom and 92 force actuators was used to calculate joint reaction forces. A finite element model of the femur used joint reaction forces calculated by the musculoskeletal model to calculate the material stresses during stepping. Balance recovery from a static forward lean angle had a greater joint contact force and greater maximum tensile stress than a recovery from treadmill induced perturbations both before and after a training session. Hip joint contact loads were found to be large in magnitude, however, all stresses experienced by the bone are less than critical yield stresses for trabecular bone. We suggest that stepping balance recovery is safe for older adults with no obvious loss of bone density or strength and that analyses such as finite element analysis are necessary to understand stresses in the material at the joint level.

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.

Comparison of Locked Plates and Blade Plates for Varus Osteotomy of the Proximal Femur by the Finite Element Method

Objective:  The present study compared the difference in load and pressure distribution behavior of the blade plate and locked plate for varus osteotomy of the proximal femur per the finite element method. Methods:  Modeling was performed by scanning a medium-sized left femur with medial valgus deformity made of polyurethane. Results:  The stiffness of the locked plate is higher compared with that of the blade plate. However, this difference was not significant. In addition, the locked plate has proximal locking screws to ensure that the bending moments on the screws are smaller during loading. Conclusion:  In summary, both plates are well-established and effective. However, the study using the finite element method plays a fundamental role in understanding the load and pressure distribution of the implant. Moreover, it opens up new possibilities for further studies, including surgical proposals and customized implant materials.

Derotational femoral osteotomy locations and their influence on joint reaction forces in dysplastic hips

Femoral version (FV) deformities are common in patients with developmental dysplasia of the hip (DDH) and may contribute to cartilage damage due to abnormal joint loading. Derotational femoral osteotomy (DFO) surgery corrects FV deformities. However there is little consensus about the femoral transection location for DFO, and its influence on joint loads is unknown. The purpose of this study was to compare the effects of two common DFO locations on muscle forces and hip joint reaction forces (JRFs) in patients with DDH. DFO was simulated in nine patients with DDH and abnormal FV using patient-specific musculoskeletal models. Femoral transection for DFO was separately simulated proximal and distal to the lesser trochanter and FV values were corrected to an idealized 15°. JRFs during early and late stance of gait were compared between the two simulated transection locations. Most changes to JRFs were similar between proximal and distal DFO, however, statistically significant differences were found for the medial JRF component during late stance among patients with femoral anteversion (p = 0.01). Force changes from five hip muscles were significantly different between DFO locations, however, changes were minimal. Most changes after DFO in patients with femoral retroversion were opposite of those with femoral anteversion, with anteroposterior and superior JRFs increasing after retroversion correction. After DFO correction, superior and medial JRFs in DDH patients remained elevated compared to controls. Understanding the influence of DFO location on muscle-generated hip forces can help surgeons justify decisions and potentially standardize surgical correction of FV deformities in patients with DDH.

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.

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.

Can pelvic tilt cause cam morphology? A computational model of proximal femur development mechanobiology

Cam deformity of the proximal femur is a risk factor for early osteoarthritis. While cam morphology is related to mechanical force at a formative time in skeletal growth, the specific problematic forces contributing to the development of cam morphology remain unknown. Individuals with femoroacetabular impingement syndrome exhibit an increased anterior pelvic tilt during walking, which alters their hip joint forces. This study aims to investigate the influence of altered joint force caused by anterior pelvic tilt on proximal femur epiphyseal growth and the potential association between increased anterior pelvic tilt and the development of cam morphology. A computational model is utilized to simulate the endochondral ossification in the proximal femur and predict cam formation. Cartilage growth and ossification patterns for a gait cycle with and without anterior pelvic tilt were modeled. The simulated growth results indicated an increased alpha angle (53° for typically developing to 68° for anterior pelvic tilt) and aspherical femoral head in the model with anterior pelvic tilt. We conclude that anterior pelvic tilt may be sufficient to cause the formation of the cam morphology. Identifying the critical mechanical conditions that increase the risk of cam deformity could help prevent this condition by adjusting the physical activities before skeletal maturity.

Differences in hip musculoskeletal loads between limbs during daily activities in patients with 3D-printed hemipelvic reconstructions following tumor surgery

BACKGROUND

Anatomical custom-made prostheses, thanks to computer-aided design and 3D-printing technology, help improve osseointegration and reduce mechanical complications in bone reconstructions following bone tumors. A recent quantitative analysis of long-term recovery in patients with 3D-printed reconstructions following pelvic tumor surgery showed asymmetries in ground reaction forces between limbs during different motor activities, while standing very good motor performance and quality of life.

RESEARCH QUESTION

We analyzed hip contact forces and muscle forces in that cohort of six patients with an innovative custom-made reconstruction of the hemipelvis, and we tested the hypothesis that asymmetries in ground reaction forces would result in more marked differences in musculoskeletal forces.

METHODS

State-of-the-art musculoskeletal modeling in an optimization-based inverse-dynamics workflow was used to calculate hip contact forces and muscle forces during five motor activities, and the differences between limbs were statistically evaluated across the motor activity cycles and on the force peaks.

RESULTS

The musculoskeletal loads were found to be not symmetric, as hip loads were generally higher in the contralateral limb. We found significant differences in considerable portions of the motor activities cycles except squat, load symmetry indices indicating a load increase (median up to 25%) on the contralateral limb, especially during stair descent and chair rise/sit, and significantly higher values in the contralateral limb at force peaks.

SIGNIFICANCE

We confirmed the hypothesis that residual asymmetries found in ground reaction forces were amplified when hip musculoskeletal loads were investigated, reflecting a shift of the loads toward the intact limb. Despite the general trend of higher loads found in the contralateral hip, this cannot be considered a risk of overloading, as both hips supported loads in a physiological range or lower, indicating a likely optimal recovery.

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 functional role of hip muscles during gait in patients with increased femoral anteversion

BACKGROUND

Femoral anteversion affects the lever arm and moment-generating capacity of the hip abductors, while an increased hip internal rotation during walking was proposed to be a compensatory mechanism to restore the abductive lever arm. Children with isolated increased femoral anteversion, however, do not always present a deficit in the net hip abduction moment during gait, suggesting that a more comprehensive understanding of the effect of morphology and motion on muscle forces and moments is needed to aid clinical decision making.

RESEARCH QUESTION

Are muscle contributions to hip joint moments and muscle forces altered in patients with increased femoral anteversion and internally rotated gait pattern compared to a control group of typically developing children? And how would the functional role of the muscle be altered if the patients walked straight?

METHODS

This follow-up study compared patients with increased femoral anteversion (n = 42, 12.8 ± 1.9 years, femoral anteversion: 39.6 ± 6.9°) to controls (n = 9, 12.0 ± 3.0 years, femoral anteversion: 18.7 ± 4.1°). Muscle forces and moment contributions were calculated using personalized musculoskeletal models. Additionally, a hypothetical scenario, in which the gait of the controls was modelled with an anteverted femoral morphology, was used to understand what would happen if the patients walked straight.

RESULTS

Gluteus medius abductive contribution was lower in patients compared to controls, despite a comparable net abduction moment around the hip. Patients presented lower muscle forces. However, if modelled to walk straight, they would require higher forces as well as a larger co-contraction of both hip internal and external rotators in the transversal plane.

SIGNIFICANCE

This study suggests that patients with increased femoral anteversion walking with an internally rotated gait pattern present lower muscle forces, but when modelled to walk straight muscle forces increase. The current results provide important information to better understand this condition and improve treatment recommendations in these patients.

Using a finite element model of the thumb to study Trapeziometacarpal joint contact during lateral pinch

BACKGROUND

Finite element (FE) analysis is widely used in different fields of orthopaedic surgery, however, its application to the trapeziometacarpal joint has been limited due to the small size, complex biconcave-convex joint geometry, and complex musculature. The goal of this study was to improve upon existing models by creating a muscle-driven FE thumb model and use the model to simulate the biomechanical effect of hand therapy exercises and ligament reconstructive surgeries.

METHODS

Bone and cartilage geometry were based on a CT dataset of a subject performing a static lateral pinch task. A previously validated musculoskeletal model was utilized to extract electromyography (EMG)-driven muscle forces. Five ligaments with biomechanical significance were modeled as springs using literature values and attached according to their anatomical landmarks.

FINDINGS

The biomechanical consequence of various interventions was proxied as a change in the maximum cartilage stress. The result shows tightening the dorsal ligament complex (dorsal radial ligament, dorsal central ligament, posterior oblique ligament) is the most effective, achieving a stress reduction of 4.8%. Five exercises used in hand therapies were modeled, among which thenar eminence strengthening showed the most prominent stress reduction of 4.0%. Four ligament reconstructive surgeries were modeled, with Eaton-Littler reconstruction showed the most significant stress reduction of 25.0%.

INTERPRETATION

Among the routinely utilized treatment options for early thumb osteoarthritis, we found that three methods: dorsal ligament imbrication, thenar eminence exercise, and the Eaton-Littler method may confer biomechanical advantages cartilage loading. These advantages align with the clinically observed favorable outcomes.

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.

The role of computational models in mechanobiology of growing bone

Endochondral ossification, the process by which long bones grow in length, is regulated by mechanical forces. Computational models, specifically finite element models, have been used for decades to understand the role of mechanical loading on endochondral ossification. This perspective outlines the stages of model development in which models are used to: 1) explore phenomena, 2) explain pathologies, 3) predict clinical outcomes, and 4) design therapies. As the models progress through the stages, they increase in specificity and biofidelity. We give specific examples of models of endochondral ossification and expect models of other mechanobiological systems to follow similar development stages.

Increased Femoral Anteversion Does Not Lead to Increased Joint Forces During Gait in a Cohort of Adolescent Patients

Orthopedic complications were previously reported for patients with increased femoral anteversion. A more comprehensive analysis of the influence of increased femoral anteversion on joint loading in these patients is required to better understand the pathology and its clinical management. Therefore, the aim was to investigate lower-limb kinematics, joint moments and forces during gait in adolescent patients with increased, isolated femoral anteversion compared to typically developing controls. Secondly, relationships between the joint loads experienced by the patients and different morphological and kinematic features were investigated. Patients with increased femoral anteversion (n = 42, 12.8 ± 1.9 years, femoral anteversion: 39.6 ± 6.9°) were compared to typically developing controls (n = 9, 12.0 ± 3.0 years, femoral anteversion: 18.7 ± 4.1°). Hip and knee joint kinematics and kinetics were calculated using subject-specific musculoskeletal models. Differences between patients and controls in the investigated outcome variables (joint kinematics, moments, and forces) were evaluated through statistical parametric mapping with Hotelling T2 and t-tests (α = 0.05). Canonical correlation analyses (CCAs) and regression analyses were used to evaluate within the patients’ cohort the effect of different morphological and kinematic predictors on the outcome variables. Predicted compressive proximo-distal loads in both hip and knee joints were significantly reduced in patients compared to controls. A gait pattern characterized by increased knee flexion during terminal stance (KneeFlex_tSt) was significantly correlated with hip and knee forces, as well as with the resultant force exerted by the quadriceps on the patella. On the other hand, hip internal rotation and in-toeing, did not affect the loads in the joints. Based on the finding of the CCAs and linear regression analyses, patients were further divided into two subgroups based KneeFlex_tSt. Patients with excessive KneeFlex_tSt presented a significantly higher femoral anteversion than those with normal KneeFlex_tSt. Patients with excessive KneeFlex_tSt presented significantly larger quadriceps forces on the patella and a larger posteriorly-oriented shear force at the knee, compared to patients with normal KneeFlex_tSt, but both patients’ subgroups presented only limited differences in terms of joint loading compared to controls. This study showed that an altered femoral morphology does not necessarily lead to an increased risk of joint overloading, but instead patient-specific kinematics should be considered.

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

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

Altered Muscle Contributions are Required to Support the Stance Limb During Voluntary Toe-Walking

Toe-walking characterizes several neuromuscular conditions and is associated with a reduction in gait stability and efficiency, as well as in life quality. The optimal choice of treatment depends on a correct understanding of the underlying pathology and on the individual biomechanics of walking. The objective of this study was to describe gait deviations occurring in a cohort of healthy adult subjects when mimicking a unilateral toe-walking pattern compared to their normal heel-to-toe gait pattern. The focus was to characterize the functional adaptations of the major lower-limb muscles which are required in order to toe walk. Musculoskeletal modeling was used to estimate the required muscle contributions to the joint sagittal moments. The support moment, defined as the sum of the sagittal extensive moments at the ankle, knee, and hip joints, was used to evaluate the overall muscular effort necessary to maintain stance limb stability and prevent the collapse of the knee. Compared to a normal heel-to-toe gait pattern, toe-walking was characterized by significantly different lower-limb kinematics and kinetics. The altered kinetic demands at each joint translated into different necessary moment contributions from most muscles. In particular, an earlier and prolonged ankle plantarflexion contribution was required from the soleus and gastrocnemius during most of the stance phase. The hip extensors had to provide a higher extensive moment during loading response, while a significantly higher knee extension contribution from the vasti was necessary during mid-stance. Compensatory muscular activations are therefore functionally required at every joint level in order to toe walk. A higher support moment during toe-walking indicates an overall higher muscular effort necessary to maintain stance limb stability and prevent the collapse of the knee. Higher muscular demands during gait may lead to fatigue, pain, and reduced quality of life. Toe-walking is indeed associated with significantly larger muscle forces exerted by the quadriceps to the patella and prolonged force transmission through the Achilles tendon during stance phase. Optimal treatment options should therefore account for muscular demands and potential overloads associated with specific compensatory mechanisms.

Impaired bone quality in the superolateral femoral neck occurs independent of hip geometry and bone mineral density

Skeletal adaptation is substantially influenced by mechanical loads. Osteocytes and their lacuno-canalicular network have been identified as a key player in load sensation and bone quality regulation. In the femoral neck, one of the most common fracture sites, a complex loading pattern with lower habitual loading in the superolateral neck and higher compressive stresses in the inferomedial neck is present. Variations in the femoral neck-shaft angle (NSA), i.e., coxa vara or coxa valga, provide the opportunity to examine the influence of loading patterns on bone quality. We obtained femoral neck specimens of 28 osteoarthritic human subjects with coxa vara, coxa norma and coxa valga during total hip arthroplasty. Bone mineral density (BMD) was assessed preoperatively by dual energy X-ray absorptiometry (DXA). Cortical and trabecular microstructure and three-dimensional osteocyte lacunar characteristics were assessed in the superolateral and inferomedial neck using ex vivo high resolution micro-computed tomography. Additionally, BMD distribution and osteocyte lacunar characteristics were analyzed by quantitative backscattered electron imaging (qBEI). All groups presented thicker inferomedial than superolateral cortices. Furthermore, the superolateral site exhibited a lower osteocyte lacunar density along with lower lacunar sphericity than the inferomedial site, independent of NSA. Importantly, BMD and corresponding T-scores correlated with microstructural parameters at the inferomedial but not superolateral neck. In conclusion, we provide micromorphological evidence for fracture vulnerability of the superolateral neck, which is independent of NSA and BMD. The presented bone qualitative data provide an explanation why DXA may be insufficient to predict a substantial proportion of femoral neck fractures. STATEMENT OF SIGNIFICANCE: The femoral neck, one of the most common fracture sites, is subject to a complex loading pattern. Site-specific differences (i.e., superolateral vs. inferomedial) in bone quality influence fracture risk, but it is unclear how this relates to hip geometry and bone mineral density (BMD) measurements in vivo. Here, we examine femoral neck specimens using a variety of high-resolution imaging techniques and demonstrate impaired bone quality in the superolateral compared to the inferomedial neck. Specifically, we found impaired cortical and trabecular microarchitecture, mineralization, and osteocyte properties, regardless of neck-shaft angle. Since BMD correlated with bone quality of the inferomedial but not the superolateral neck, our results illustrate why bone densitometry may not predict a substantial proportion of femoral neck fractures.

Femoral version deformities alter joint reaction forces in dysplastic hips during gait

Developmental dysplasia of the hip (DDH) causes hip instability and early-onset osteoarthritis. The focus on pathomechanics in DDH has centered on the shallow acetabulum, however there is growing awareness of the role of femoral deformities in joint damage. The objective of this study was to determine the influence of femoral version (FV) on the muscle and joint reaction forces (JRFs) of dysplastic hips during gait. Magnetic resonance images, in-vivo gait data, and musculoskeletal models were used to calculate JRFs and simulate changes due to varying FV deformities. Rotation about the long axis of the femur was added in the musculoskeletal models to simulate FV values from -5° (relative retroversion) to + 35° (increased anteversion). In our simulations, FV deformities caused the largest changes to the anteroposterior and resultant JRFs. From a normal FV of 15°, a 15° increase in femoral anteversion caused JRFs to be less posterior in early stance (Δ = 0.43 ± 0.22 xbodyweight) and more anterior in late stance (Δ = 0.60 ± 14 xbodyweight). Relative retroversion caused anteroposterior changes that were similar to anteversion in early stance but opposite in late stance. Resultant JRFs experienced the largest changes during late stance where anteversion raised the peak by 0.48 ± 0.15 xbodyweight and relative retroversion lowered the peak by 0.32 ± 0.30 xbodyweight. Increasing anteversion increased hip flexor and abductor muscle forces, which caused the changes in JRFs. Identifying how FV deformities influence hip joint loading can elucidate their role in the mechanisms of hip degeneration in patients with DDH.

A Direct Comparison of Node and Element-Based Finite Element Modeling Approaches to Study Tissue Growth

Finite element analysis is a useful tool to model growth of biological tissues and predict how growth can be impacted by stimuli. Previous work has simulated growth using node-based or element-based approaches, and this implementation choice may influence predicted growth, irrespective of the applied growth model. This study directly compared node-based and element-based approaches to understand the isolated impact of implementation method on growth predictions by simulating growth of a bone rudiment geometry, and determined what conditions produce similar results between the approaches. We used a previously reported node-based approach implemented via thermal expansion and an element-based approach implemented via osmotic swelling, and we derived a mathematical relationship to relate the growth resulting from these approaches. We found that material properties (modulus) affected growth in the element-based approach, with growth completely restricted for high modulus values relative to the growth stimulus, and no restriction for low modulus values. The node-based approach was unaffected by modulus. Node- and element-based approaches matched marginally better when the conversion coefficient to relate the approaches was optimized based on the results of initial simulations, rather than using the theoretically predicted conversion coefficient (median difference in node position 0.042 cm versus 0.052 cm, respectively). In summary, we illustrate here the importance of the choice of implementation approach for modeling growth, provide a framework for converting models between implementation approaches, and highlight important considerations for comparing results in prior work and developing new models of tissue growth.

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.

Subject-Specific Modeling of Femoral Torsion Influences the Prediction of Hip Loading During Gait in Asymptomatic Adults

Hip osteoarthritis may be caused by increased or abnormal intra-articular forces, which are known to be related to structural articular cartilage damage. Femoral torsional deformities have previously been correlated with hip pain and labral damage, and they may contribute to the onset of hip osteoarthritis by exacerbating the effects of existing pathoanatomies, such as cam and pincer morphologies. A comprehensive understanding of the influence of femoral morphotypes on hip joint loading requires subject-specific morphometric and biomechanical data on the movement characteristics of individuals exhibiting varying degrees of femoral torsion. The aim of this study was to evaluate hip kinematics and kinetics as well as muscle and joint loads during gait in a group of adult subjects presenting a heterogeneous range of femoral torsion by means of personalized musculoskeletal models. Thirty-seven healthy volunteers underwent a 3D gait analysis at a self-selected walking speed. Femoral torsion was evaluated with low-dosage biplanar radiography. The collected motion capture data were used as input for an inverse dynamics analysis. Personalized musculoskeletal models were created by including femoral geometries that matched each subject’s radiographically measured femoral torsion. Correlations between femoral torsion and hip kinematics and kinetics, hip contact forces (HCFs), and muscle forces were analyzed. Within the investigated cohort, higher femoral antetorsion led to significantly higher anteromedial HCFs during gait (medial during loaded stance phase and anterior during swing phase). Most of the loads during gait are transmitted through the anterior/superolateral quadrant of the acetabulum. Correlations with hip kinematics and muscle forces were also observed. Femoral antetorsion, through altered kinematic strategies and different muscle activations and forces, may therefore lead to altered joint mechanics and pose a risk for articular damage. The method proposed in this study, which accounts for both morphological and kinematic characteristics, might help in identifying in a clinical setting patients who, as a consequence of altered femoral torsional alignment, present more severe functional impairments and altered joint mechanics and are therefore at a higher risk for cartilage damage and early onset of hip osteoarthritis.

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.

ESB Clinical Biomechanics Award 2020: Pelvis and hip movement strategies discriminate typical and pathological femoral growth - Insights gained from a multi-scale mechanobiological modelling framework

BACKGROUND

Many children with cerebral palsy (CP) develop skeletal deformities during childhood. So far, it is unknown why some children with CP develop bony deformities whereas others do not. The aims of this study were to (i) investigate what loading characteristics lead to typical and pathological femoral growth, and (ii) evaluate why some children with CP develop femoral deformities whereas other do not.

METHODS

A multi-scale mechanobiological modelling workflow was used to simulate femoral growth based on three-dimensional motion capture data of six typically developing children and 16 children with CP. Based on the growth results, the participants with CP were divided into two groups: typical growth group and pathological growth group. Gait kinematics and femoral loading were compared between simulations resulting in typical growth and those resulting in pathologic growth.

FINDINGS

Hip joint contact forces were less posteriorly-oriented in the pathological growth simulations compared to the typical ones. Compared to the typically developing participants, the CP group with pathological femoral growth presented increased knee flexion and no hip extension. The CP group with simulated typical growth presented similar sagittal plane joint kinematics but differed in the frontal plane pelvic and hip movement strategy, which normalized the hip joint contact force and therefore contributed to typical femoral growth trends.

INTERPRETATION

Our simulation results identified specific gait features, which may contribute to pathological femoral growth. Furthermore, the hip joint contact force orientation in the sagittal plane seems to be the dominant factor for determining femoral growth simulations.

Dependency of lower limb joint reaction forces on femoral version

BACKGROUND

Musculoskeletal (MSK) models based on literature data are meant to represent a generic anatomy and are a popular tool employed by biomechanists to estimate the internal loads occurring in the lower limb joints, such as joint reaction forces (JRFs). However, since these models are normally just linearly scaled to an individual's anthropometry, it is unclear how their estimations would be affected by the personalization of key features of the MSK anatomy, one of which is the femoral version angle.

RESEARCH QUESTION

How are the lower limb JRF magnitudes computed through a generic MSK model affected by changes in the femoral version?

METHODS

We developed a bone-deformation tool in MATLAB (shared at https://simtk.org/projects/bone_deformity) and used it to create a set of seven OpenSim models spanning from 2˚ femoral retroversion to 40˚ anteversion. We used these models to simulate the gait of an elderly individual with an instrumented prosthesis implanted at their knee joint (5th Grand Challenge dataset) and quantified both the changes in JRFs magnitude due to varying the skeletal anatomy and their accuracy against the correspondent in vivo measurements at the knee joint.

RESULTS

Hip and knee JRF magnitudes were affected by the femoral version with variations from the unmodified generic model up to 17.9 ± 4.5% at the hip and 43.4 ± 27.1% at the knee joint. The ankle joint was unaffected by the femoral geometry. The MSK models providing the most accurate knee JRFs (root mean squared error: 0.370 ± 0.068 body weight, coefficient of determination: 0.757 ± 0.104, peak error range: 0.09-0.42 body weight) were those with femoral anteversion angle closer to that measured on the segmented bone of the individual.

SIGNIFICANCE

Femoral version substantially affects hip and knee JRFs estimated with generic MSK models, suggesting that personalizing key MSK anatomical features might be necessary for accurate estimation of JRFs with these models.

Torsion Tool: An automated tool for personalising femoral and tibial geometries in OpenSim musculoskeletal models

Common practice in musculoskeletal modelling is to use scaled musculoskeletal models based on a healthy adult, but this does not consider subject-specific geometry, such as tibial torsion and femoral neck-shaft and anteversion angles (NSA and AVA). The aims of this study were to (1) develop an automated tool for creating OpenSim models with subject-specific tibial torsion and femoral NSA and AVA, (2) evaluate the femoral component, and (3) release the tool open-source. The Torsion Tool (https://simtk.org/projects/torsiontool) is a MATLAB-based tool that requires an individual's tibial torsion, NSA and AVA estimates as input and rotates corresponding bones and associated muscle points of a generic musculoskeletal model. Performance of the Torsion Tool was evaluated comparing femur bones as personalised with the Torsion Tool and scaled generic femurs with manually segmented bones as golden standard for six typically developing children and thirteen children with cerebral palsy. The tool generated femur geometries closer to the segmentations, with lower maximum (-19%) and root mean square (-18%) errors and higher Jaccard indices (+9%) compared to generic femurs. Furthermore, the tool resulted in larger improvements for participants with higher NSA and AVA deviations. The Torsion Tool allows an automatic, fast, and user-friendly way of personalising femoral and tibial geometry in an OpenSim musculoskeletal model. Personalisation is expected to be particularly relevant in pathological populations, as will be further investigated by evaluating the effects on simulation outcomes.

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.

Influence of simulated hip muscle weakness on hip joint forces during deep squatting

This study aimed to determine the effects of simulated hip muscle weakness on changes in hip joint forces during deep squat motion. Ten healthy individuals performed squat motion at three different positions (0° foot angle [N-squat], 10° toe-in [IN-squat], and 30° toe-out [OUT-squat]). A scaled musculoskeletal model for each participant was used to calculate the muscle and hip joint forces. For each hip muscle, models of full strength, mild muscle weakness (15% decrease), and severe muscle weakness (30% decrease) were created. The muscles affecting the hip joint forces were identified, and the rate of change in the joint forces was compared among the three squat conditions. The anterior hip joint force was increased in the muscle weakness models of the inferior gluteus maximus (iGlutMax) and iGlutMax+deep external rotator (ExtRot) muscles. With 30% muscle weakness of these muscles, statistically significant differences in the rate of increase in the anterior joint force were observed in the following order: IN-squat (iGlutMax, 29.5%; iGlutMax+ExtRot, 41.4%), N-squat (iGlutMax, 18.3%; iGlutMax+ExtRot, 27.8%), and OUT-squat (iGlutMax, 5.6%; iGlutMax+ExtRot, 9.3%). OUT-squat may be recommended to minimize the increase in hip joint forces if accompanied by hip muscle weakness.

Paggiosi M, McCloskey EV, Viceconti M, Mazzà C, Li X. Femoral neck strain prediction during level walking using a combined musculoskeletal and finite element model approach

Recently, coupled musculoskeletal-finite element modelling approaches have emerged as a way to investigate femoral neck loading during various daily activities. Combining personalised gait data with finite element models will not only allow us to study changes in motion/movement, but also their effects on critical internal structures, such as the femur. However, previous studies have been hampered by the small sample size and the lack of fully personalised data in order to construct the coupled model. Therefore, the aim of this study was to build a pipeline for a fully personalised multiscale (body-organ level) model to investigate the strain levels at the femoral neck during a normal gait cycle. Five postmenopausal women were included in this study. The CT and MRI scans of the lower limb, and gait data were collected for all participants. Muscle forces derived from the body level musculoskeletal models were used as boundary constraints on the finite element femur models. Principal strains were estimated at the femoral neck region during a full gait cycle. Considerable variation was found in the predicted peak strain among individuals with mean peak first principal strain of 0.24% ± 0.11% and mean third principal strain of -0.29% ± 0.24%. For four individuals, two overall peaks of the maximum strains were found to occur when both feet were in contact with the floor, while one individual had one peak at the toe-off phase. Both the joint contact forces and the muscular forces were found to substantially influence the loading at the femoral neck. A higher correlation was found between the predicted peak strains and the gluteus medius (R² ranged between 0.95 and 0.99) than the hip joint contact forces (R² ranged between 0.63 and 0.96). Therefore, the current findings suggest that personal variations are substantial, and hence it is important to consider multiple subjects before deriving general conclusions for a target population.

The effect of hip muscle weakness and femoral bony deformities on gait performance

BACKGROUND

Children with cerebral palsy (CP) present with a pathological gait pattern due to musculoskeletal impairments, such as muscle weakness and altered bony geometry. However, the effect of these impairments on gait performance is still unknown. Research aim:This study aimed to explore the effect of hip muscle weakness and femoral deformities on the gait performance of CP and typical developing (TD) subjects.

METHODS

6400 musculoskeletal models were created by weakening the hip extensors, abductors, adductors and flexors from 0% to 75 % and increasing the femoral anteversion angle (FAA) and neck shaft angle (NSA) from 20° to 60° and 120° to 160°, respectively. One TD and five CP gait patterns were imposed to each model and muscle forces were calculated. The effect of weakness and bony deformities on the capability gap (CG) at the hip, i.e. the lack in hip moment generating capacity to perform the gait pattern, was investigated using regression analysis.

RESULTS

The CG of apparent equinus, stiff knee gait, TD gait, jump gait and true equinus increased with 0.080, 0.038, 0.015, 0.023 and 0.005 Nm/kg per 10 percent hip abductor weakness increase, with 0.211, 0.130, 0.056, 0.045 and 0.011 Nm/kg per 10 degrees FAA increase and with 0.163, 0.080, 0.036, 0.043 and 0.011 Nm/kg per 10 degrees NSA increase, respectively. Combined weakness and bony deformities explained 96 %, 85 %, 82 %, 65 %, 40 % and 35 % of the variance in the CG of apparent equinus, TD gait, stiff knee gait, jump gait, true equinus and crouch gait, respectively.

SIGNIFICANCE

The results suggest that surgical correction of femoral deformities is more likely to be effective than strength training of hip muscles in enhancing CP gait performance. Jump gait, true equinus and especially crouch were more robust, while apparent equinus and stiff knee gait were limited by hip weakness and femoral deformities.