Compressive tibiofemoral force during crouch gait

Feasibility of predicting changes in gait biomechanics following muscle strength perturbations using optimal control in patients with transfemoral amputation

Bone-anchored limbs (BALs) are socket prosthesis alternatives, directly fixing to residual bone via osseointegrated implant. There is a need to quantify multi-level effects of rehabilitation for transfemoral BAL users (i.e. changes in joint loading and movement patterns). Our primary objective was determining feasibility of using optimal control to predict gait biomechanics compared to ground-truth experimental data from transfemoral BAL users. A secondary objective was examining biomechanical effects from estimated changes in hip abductor muscle strength. We developed and validated a workflow for predicting gait biomechanics in four transfemoral BAL users and investigated the biomechanical effects of altered hip abductor strengths.

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.

The effect of impaired unilateral ankle propulsion on contralateral knee joint loading

BACKGROUND

Impairments in unilateral ankle propulsion may result from restriction by an external device or pathology such as lower limb amputation. Models of gait suggest this reduction may lead to increased collisional force on the contralateral side, potentially increasing force through the knee and increasing the risk of knee pain or osteoarthritis.

RESEARCH QUESTION

How do restrictions in unilateral ankle propulsive force affect contralateral knee joint loading in otherwise healthy individuals?

METHODS

18 individuals without impairment walked on a treadmill at 1.5 m/s for two conditions: one free of restrictions, and one where a randomized limb's ankle propulsive force was restricted using an off-the-shelf ankle-foot orthosis (AFO). Ankle propulsive power, lower extremity joint work, and ground reaction force variables were calculated for the final 3 gait cycles of each condition. Tibiofemoral joint contact force (TJCF) for the limb contralateral to the AFO was calculated through a standard OpenSim workflow utilizing the gait2392 model. Intra-limb pair-wise comparisons were made between conditions.

RESULTS

Compared to walking unrestricted, the limb wearing the AFO demonstrated a significant reduction in peak ankle propulsive power and positive ankle work by approximately 50 % each (p<0.01). With ankle restriction, the ipsilateral knee significantly increased positive work (p<0.01). The overall propulsion produced by that limb did not change between conditions, demonstrated by a lack of change in anterior ground reaction force impulse (p=0.11). The knee of the limb contralateral to the AFO did not display differences in any TJCF variable between conditions (all p>0.07).

SIGNIFICANCE

These results suggest a unilateral deficit in ankle propulsion will not increase contralateral knee joint forces in individuals who are able to use other joints of the limb to compensate for the loss of ankle function. However, further research should investigate this relationship in those who display pathologies that may prevent more proximal compensations.

Predicting Knee Joint Contact Forces During Normal Walking Using Kinematic Inputs With a Long-Short Term Neural Network

Knee joint contact forces are commonly estimated via surrogate measures (i.e., external knee adduction moments or musculoskeletal modeling). Despite its capabilities, modeling is not optimal for clinicians or persons with limited experience. The purpose of this study was to design a novel prediction method for knee joint contact forces that is simplistic in terms of required inputs. This study included marker trajectories and instrumented knee forces during normal walking from the "Grand Challenge" (n = 6) and "CAMS" (n = 2) datasets. Inverse kinematics were used to derive stance phase hip (sagittal, frontal, transverse), knee (sagittal, frontal), ankle (sagittal), and trunk (frontal) kinematics. A long-short term memory network (LSTM) was created using matlab to predict medial and lateral knee force waveforms using combinations of the kinematics. The Grand Challenge and CAMS datasets trained and tested the network, respectively. Musculoskeletal modeling forces were derived using static optimization and joint reaction tools in OpenSim. Waveform accuracy was determined as the proportion of variance and root-mean-square error between network predictions and in vivo data. The LSTM network was highly accurate for medial forces (R2 = 0.77, RMSE = 0.27 BW) and required only frontal hip and knee and sagittal hip and ankle kinematics. Modeled medial force predictions were excellent (R2 = 0.77, RMSE = 0.33 BW). Lateral force predictions were poor for both methods (LSTM R2 = 0.18, RMSE = 0.08 BW; modeling R2 = 0.21, RMSE = 0.54 BW). The designed LSTM network outperformed most reports of musculoskeletal modeling, including those reached in this study, revealing knee joint forces can accurately be predicted by using only kinematic input variables.

Investigation of joint contact forces during walking in the subjects with toe in gait due to increasing in femoral head anteversion angle

Toe-in gait is a pathology in which the child walks and turns the foot inward instead of pointing straight ahead. The alignment of the lower limb structure changes in this disease, increasing the incidence of knee and hip osteoarthritis. This study aimed to determine the kinematic and joint loading in subjects walking with a toe-in gait pattern. This study selected two groups of subjects: normal subjects and those with toe-in gait due to an increased femoral head anteversion angle (each group consisted of 15 subjects). A Qualisys motion analysis system and a Kistler force plate were used to record the motions and forces applied to the leg while walking. OpenSim software (version 3.3) was used to analyze the range of motion, moments, muscle forces, and joint contact forces in both groups of subjects. The mean values of stride length for normal subjects (1.1 ± 0.141 m) and those with toe-in gait (0.94 ± 0.183 m) differed significantly. The mediolateral component of the ground reaction force decreased significantly in the toe-in gait group compared to normal subjects (p-value = 0.05). The peak force of most of the hip joint muscles increased significantly in those with toe-in gait compared to normal subjects (p-value < 0.05). The results of this study showed that those with toe-in gait, due to an increase in femoral head anteversion angle, only had a change in rotation of the pelvic and hip joint. There was no significant difference between walking speed and most ground reaction force components between normal subjects and those with toe-in gait. As the peaks of most of the hip joint muscles increased significantly in the toe-in gait group, this increased joint contact forces (especially the anteroposterior component of the hip joint and the mediolateral component of the knee joint), which may ultimately increase the incidence of hip and knee joint osteoarthritis.

Sensitivity of knee cartilage biomechanics in finite element analysis to selected Musculoskeletal models

Knee joint kinematics and kinetics analyzed by musculoskeletal (MS) modeling are often utilized in finite element (FE) models, estimating tissue-level mechanical responses. We compared knee cartilage stresses, strains, and centers of pressure of FE models driven by two widely used MS models, implemented in AnyBody and OpenSim. Minor discrepancies in the results were observed between the models. AnyBody-driven FE models showed slightly higher stresses in the medial tibial cartilage, while OpenSim-driven FE models estimated more anterior and lateral center of pressure. Recognizing these differences in the MS-FE models is important to ensure reliable analysis of cartilage mechanics and failure and simulation of rehabilitation.

Muscle synergy-informed neuromusculoskeletal modelling to estimate knee contact forces in children with cerebral palsy

Cerebral palsy (CP) includes a group of neurological conditions caused by damage to the developing brain, resulting in maladaptive alterations of muscle coordination and movement. Estimates of joint moments and contact forces during locomotion are important to establish the trajectory of disease progression and plan appropriate surgical interventions in children with CP. Joint moments and contact forces can be estimated using electromyogram (EMG)-informed neuromusculoskeletal models, but a reduced number of EMG sensors would facilitate translation of these computational methods to clinics. This study developed and evaluated a muscle synergy-informed neuromusculoskeletal modelling approach using EMG recordings from three to four muscles to estimate joint moments and knee contact forces of children with CP and typically developing (TD) children during walking. Using only three to four experimental EMG sensors attached to a single leg and leveraging an EMG database of walking data of TD children, the synergy-informed approach estimated total knee contact forces comparable to those estimated by EMG-assisted approaches that used 13 EMG sensors (children with CP, n = 3, R2 = 0.95 ± 0.01, RMSE = 0.40 ± 0.14 BW; TD controls, n = 3, R2 = 0.93 ± 0.07, RMSE = 0.19 ± 0.05 BW). The proposed synergy-informed neuromusculoskeletal modelling approach could enable rapid evaluation of joint biomechanics in children with unimpaired and impaired motor control within a clinical environment.

Hamstrings vibration reduces tibiofemoral compressive force following anterior cruciate ligament reconstruction

Individuals who have undergone anterior cruciate ligament reconstruction (ACLR) are at greater risk of developing knee osteoarthritis (OA). This elevated risk of knee OA is associated with high tibiofemoral (TF) compressive force, due to a combination of low knee flexion angles and increased co-contraction of the hamstrings and quadriceps during limb loading. Prolonged vibration of the hamstrings fatigues the intrafusal muscle fibers, which reduces autonomic reflexive excitation of the hamstrings and alleviates reciprocal inhibition to the quadriceps. The aim of this study was to examine the effect of prolonged hamstrings vibration on TF compressive force in individuals who have undergone ACL reconstruction. Fourteen participants with unilateral ACLR and 14 participants without knee injury performed a single-leg drop-land task before and after prolonged (20 min) vibration of the hamstrings. Peak TF compressive force, knee flexion angle, and hamstrings/quadriceps co-contraction were calculated during the deceleration phase of the drop-land task before and after vibration. The ACLR group experienced an 18% decrease in TF compressive force, a 32% increase in knee flexion angle, and a 38% decrease in hamstrings/quadriceps co-contraction after hamstrings vibration. There was no difference in any of the parameters in the noninjured group after vibration. These data suggest that acute prolonged hamstrings vibration has the potential to mitigate TF compressive force, which may protect the knee joint in the long term. Clinical significance: The results of this research are expected to lead to improved clinical care for ACLR patients because it holds promise for mitigating altered joint mechanics and perhaps slowing down the onset of posttraumatic knee osteoarthritis.

Chronic ankle instability affects the association between knee joint angle and loading: musculoskeletal simulation study

The purpose of this study was to calculate and compare (1) knee loads, (2) muscle-specific contributions to knee loads, and (3) effects of knee flexion angle on knee loads and muscle-specific load contributions during a forward jump-landing task in people with and without chronic ankle instability (CAI). Eight CAI patients and seven healthy controls performed a forward jump-landing task. We collected 3D kinematics, ground reaction force, and muscle activation and used musculoskeletal modeling. The results showed that only healthy controls exhibited an association between knee flexion angle and knee compressive impulse (r = 0.854, p = .014). The lack of association in CAI group may lead to knee instability and increase knee injury risk in people with CAI.

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.

Hip Contact Forces During Sprinting in Femoroacetabular Impingement Syndrome

PURPOSE

Sprinting often provokes hip pain in individuals with femoroacetabular impingement syndrome (FAIS). Asphericity of the femoral head-neck junction (cam morphology) characteristic of FAIS can increase the risk of anterior-superior acetabular cartilage damage. This study aimed to 1) compare hip contact forces (magnitude and direction) during sprinting between individuals with FAIS, asymptomatic cam morphology (CAM), and controls without cam morphology, and 2) identify the phases of sprinting with high levels of anteriorly directed hip contact forces.

METHODS

Forty-six recreationally active individuals with comparable levels of physical activity were divided into three groups (FAIS, 14; CAM, 15; control, 17) based on their history of hip/groin pain, results of clinical impingement tests, and presence of cam morphology (alpha angle >55°). Three-dimensional marker trajectories, ground reaction forces, and electromyograms from 12 lower-limb muscles were recorded during 10-m overground sprinting trials. A linearly scaled electromyogram-informed neuromusculoskeletal model was used to calculate hip contact force magnitude (resultant, anterior-posterior, inferior-superior, medio-lateral) and angle (sagittal and frontal planes). Between-group comparisons were made using two-sample t -tests via statistical parametric mapping ( P < 0.05).

RESULTS

No significant differences in magnitude or direction of hip contact forces were observed between FAIS and CAM or between FAIS and control groups during any phase of the sprint cycle. The highest anteriorly directed hip contact forces were observed during the initial swing phase of the sprint cycle.

CONCLUSIONS

Hip contact forces during sprinting do not differentiate recreationally active individuals with FAIS from asymptomatic individuals with and without cam morphology. Hip loading during early swing, where peak anterior loading occurs, may be a potential mechanism for cartilage damage during sprinting-related sports in individuals with FAIS and/or asymptomatic cam morphology.

Does Overhead Squat Performance Affect the Swing Kinematics and Lumbar Spine Loads during the Golf Downswing?

The performance of the overhead squat may affect the golf swing mechanics associated with golf-related low back pain. This study investigates the difference in lumbar kinematics and joint loads during the golf downswing between golfers with different overhead squat abilities. Based on the performance of the overhead squat test, 21 golfers aged 18 to 30 years were divided into the highest-scoring group (HS, N = 10, 1.61 ± 0.05 cm, and 68.06 ± 13.67 kg) and lowest-scoring group (LS, N = 11, 1.68 ± 0.10 cm, and 75.00 ± 14.37 kg). For data collection, a motion analysis system, two force plates, and TrackMan were used. OpenSim 4.3 software was used to simulate the joint loads for each lumbar joint. An independent t-test was used for statistical analysis. Compared to golfers demonstrating limitations in the overhead squat test, golfers with better performance in the overhead squat test demonstrated significantly greater angular extension displacement on the sagittal plane, smaller lumbar extension angular velocity, and smaller L4-S1 joint shear force. Consequently, the overhead squat test is a useful index to reflect lumbar kinematics and joint loading patterns during the downswing and provides a good training guide reference for reducing the risk of a golf-related lower back injury.

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.

Applying elastic resistance bands for gait training: A simulation-based study to determine how band configuration affects gait biomechanics and muscle activation

BACKGROUND

Wearable robotic exoskeletons and leg braces are desirable for gait rehabilitation because they can apply loads directly to an affected joint. Yet, they are not widely used in clinics because they are costly and complex to set up. Conversely, tethered devices, such as elastic resistance bands, are widely available in clinics, are low-cost, and are quick to set up. However, resistance bands will affect walking differently based on how they are configured to pull on the leg (e.g., pulling forward or backward).

RESEARCH QUESTION

How can a resistance band be configured to alter muscle activation and gait biomechanics based on the segment it is attached to and the angle with which it attaches?

METHODS

We used an open-source musculoskeletal modeling platform to emulate several configurations of an elastic band pulling on the ankle, calf, and thigh at various angles during non-pathological walking. We evaluated gait biomechanics and simulated muscle activation using computed muscle control (CMC) and identified a subset of four configurations with potential applications for gait training. Eight non-pathological participants then walked on a treadmill under these configurations to verify how these configurations altered muscle activation.

RESULTS

We found that muscle activity greatly varied based on the location where the elastic band is attached and the angle with which the elastic band pulls on the leg. Notably, specific angles can be used to pull on the legs to elicit an increase or decrease in muscle activation.

SIGNIFICANCE

This study provides insight into how tethered devices can be configured to provide assistance or resistance during gait training. This information can be applied when developing low-cost gait training solutions for addressing individuals' impairments.

Per-step and cumulative load at three common running injury locations: The effect of speed, surface gradient, and cadence

Understanding how loading and damage on common running injury locations changes across speeds, surface gradients, and step frequencies may inform training programs and help guide progression/rehabilitation after injuries. However, research investigating tissue loading and damage in running is limited and fragmented across different studies, thereby impairing comparison between conditions and injury locations. This study examined per-step peak load and impulse, cumulative impulse, and cumulative weighted impulse (hereafter referred to as cumulative damage) on three common injury locations (patellofemoral joint, tibia, and Achilles tendon) across different speeds, surface gradients, and cadences. We also explored how cumulative damage in the different tissues changed across conditions relative to each other. Nineteen runners ran at five speeds (2.78, 3.0, 3.33, 4.0, 5.0 m s-1 ), and four gradients (-6, -3, +3, +6°), and three cadences (preferred, ±10 steps min-1 ) each at one speed. Patellofemoral, tibial, and Achilles tendon loading and damage were estimated from kinematic and kinetic data and compared between conditions using a linear mixed model. Increases in running speed increased patellofemoral cumulative damage, with nonsignificant increases for the tibia and Achilles tendon. Increases in cadence reduced damage to all tissues. Uphill running increased tibial and Achilles tendon, but decreased patellofemoral damage, while downhill running showed the reverse pattern. Per-step and cumulative loading, and cumulative loading and cumulative damage indices diverged across conditions. Moreover, changes in running speed, surface gradient, and step frequency lead to disproportional changes in relative cumulative damage on different structures. Methodological and practical implications for researchers and practitioners are discussed.

Does patella lowering as part of multilevel surgery improve knee kinematics in children with cerebral palsy and crouch gait? A meta-analysis of comparative studies

Purpose

To evaluate differences in knee kinematic outcomes of patellar-lowering surgery, specifically patellar tendon advancement or patellar tendon shortening, compared with no-patellar-lowering surgery in multilevel surgery for children with cerebral palsy and crouch gait.

Methods

Four databases were searched to retrieve studies published from inception until 2023. Three reviewers independently screened for studies with observational or randomized control designs, comparing two groups of patients with cerebral palsy and crouch gait who underwent multilevel surgery (with patellar-lowering surgery versus no-patellar-lowering surgery), where various gait analysis outcomes were reported (CRD42023450692). The risk of bias was assessed with the Risk Of Bias In Non-randomised Studies - of Interventions (ROBINS-I) tool.

Results

Seven studies (249 patients and 368 limbs) met the eligibility criteria. Patients undergoing patellar-lowering surgery demonstrated statistically significant improvements in knee flexion at initial contact (mean difference = -6.39; 95% confidence interval = [-10.4, -2.75]; p = 0.0006; I2 = 84%), minimum knee flexion in stance (mean difference = -14.27; 95% confidence interval = [-18.31, -10.23]; p < 0.00001; I2 = 89%), and clinical knee flexion contracture (mean difference = -5.6; 95% confidence interval = [-9.59, -1.6]; p = 0.006; I2 = 95%), with a significant increase in anterior pelvic tilt (mean difference = 2.97; 95% confidence interval = [0.58, 5.36]; p = 0.01; I2 = 15%). However, improvements in gait deviation index and decrease in peak knee flexion in swing did not reach statistical significance. Subgroup analysis reduced heterogeneity and revealed (1) greater improvement using patellar tendon shortening versus patellar tendon advancement techniques; (2) lack of knee flexion contracture improvement in high-quality or longer-term studies; (3) longer-term improvement only in minimum knee flexion in stance, with a decrease in peak knee flexion in swing; and (4) an inability to assess the potential benefit of rectus femoris procedure and hamstring preservation.

Conclusions

Overall, the combination of patellar-lowering surgery with multilevel surgery demonstrated superior improvements in stance-phase knee kinematics compared with multilevel surgery alone, despite an increase in anterior pelvic tilt and a longer-term knee flexion reduction during the swing phase.

Level of evidence

Level III, Systematic review of level III studies.

Muscle contributions to reduced ankle joint contact force during drop vertical jumps in patients with chronic ankle instability

Chronic ankle instability is a condition linked to progressive early ankle joint degeneration. Patients with chronic ankle instability exhibit altered biomechanics during gait and jump landings and these alterations are believed to contribute to aberrant joint loading and subsequent joint degeneration. Musculoskeletal modeling has the capacity to estimate joint loads from individual muscle forces. However, the influence of chronic ankle instability on joint contact forces remains largely unknown. The objective of this study was to compare tri-axial (i.e., compressive, anterior-posterior, and medial-lateral) ankle joint contact forces between those with and without chronic ankle instability during the ground contact phase of a drop vertical jump. Fifteen individuals with and 15 individuals without chronic ankle instability completed drop vertical jump maneuvers in a research laboratory. We used those data to drive three-dimensional musculoskeletal simulations and estimate muscle forces and tri-axial joint contact force variables (i.e., peak and impulse). Compared to those without chronic ankle instability, the ankles of patients with chronic ankle instability underwent lower compressive ankle joint contact forces as well as lower anterior-posterior and medial-lateral shearing forces during the weight acceptance phase of landing (p <.05). These findings suggest that patients with chronic ankle instability exhibit lower ankle joint loading patterns than uninjured individuals during a drop vertical jump, which may be considered in rehabilitation to potentially reduce the risk of early onset of ankle joint degeneration.

Muscle activity during crouched walking

OBJECTIVES

Muscle activity during crouched walking has been previously studied in the context of the evolution of hominin bipedalism and human movement disorders. However, crouched walking could also be used in approach hunting where postural height (actual height of the body from the ground to the top of the head during locomotion) is the limiting factor. Here, we aim to analyze the relationship between relative postural height (%stature), kinematics, and muscle activity during crouched walking.

MATERIALS AND METHODS

Adult males (n = 19) walked with extended limbs and at three degrees of crouch while their 3D motion capture kinematics and lower limb muscle electromyography were recorded. We measured activation of tibialis anterior, soleus, gastrocnemius medialis, gastrocnemius lateralis, vastus lateralis, rectus femoris, biceps femoris, and gluteus maximus. We analyzed the effects of postural height on kinematics and muscle activation using linear mixed effects model.

RESULTS

Flexion angles, individual muscle activation (except for medial gastrocnemius), and total muscle activation were negatively related to relative postural height, that is, were greater at more crouched postures. Relative postural height had a stronger effect on the activation of the thigh and gluteal muscles compared to shank muscles.

DISCUSSION

General increase in lower limb muscle activation at lower postural heights suggests a negative relationship between relative postural height and fatigue, and may indicate a possible mechanism by which short stature could benefit the hunter in approach hunting. Greater activation of thigh and gluteal muscles relative to shank muscles may help to identify crouched walking in past human populations.

External weight mass and carrying position influence peak patellar tendon force and patellofemoral joint contact force independently during forward lunge

BACKGROUND

The forward lunge is a common exercise in the rehabilitation of patellar tendinopathy and patellofemoral pain syndrome. External weights are frequently used to increase the peak patellar tendon force and patellofemoral joint contact force during this exercise. The weight's position might influence this relationship. The objective of this study was to investigate the combined effect of an external weight's mass and carrying position on the peak patellar tendon force and patellofemoral joint contact force during a forward lunge.

METHODS

Ten healthy individuals performed forward lunges holding external weights between 0.1 and 0.3 times body mass either in one hand at the ipsilateral or contralateral side of the leading leg, or in two hands at the side or in front of the trunk. Three-dimensional kinematic data and ground reaction forces were collected and peak patellar tendon force and patellofemoral joint contact force were calculated using musculoskeletal modelling. Two-way repeated measures ANOVA's determined the main effects for the external weight's mass and position as well as their interaction effect.

FINDINGS

Increasing the mass of the external weights increased both the peak patellar tendon force and patellofemoral joint contact force linearly and at the same rate in all positions. Both peak forces were larger in the one-hand ipsilateral and two-hand side positions.

INTERPRETATION

An external weight's mass and position both influence the peak patellar tendon force and patellofemoral joint contact force during a forward lunge. The rate of increase in peak forces with increasing mass was similar for all weight-carrying positions.

Are Abnormal Muscle Biomechanics and Patient-reported Outcomes Associated in Patients With Hip Dysplasia?

BACKGROUND

Developmental dysplasia of the hip (DDH) is a major risk factor for the early development of hip osteoarthritis. Recent studies have demonstrated how DDH alters hip muscle moment arms and elevates muscle-induced biomechanical variables such as joint reaction forces and acetabular edge loads. Understanding the link between abnormal biomechanics and patient-reported outcome measures (PROMs) is important for evidence-based clinical interventions that improve patient symptoms and functional outcomes. To our knowledge, there are no reports of the relationships between muscle-induced biomechanics and PROMs.

QUESTIONS/PURPOSES

(1) Are there associations between PROMs and muscle-induced hip biomechanics during gait for patients with DDH and controls? (2) Are there associations among PROMs and separately among biomechanical variables?

METHODS

Participants in this prospective cross-sectional comparative study included 20 female patients with DDH who had no prior surgery or osteoarthritis and 15 female individuals with no evidence of hip pathology (controls) (age: median 23 years [range 16 to 39 years]; BMI: median 22 kg/m 2 [range 17 to 27 kg/m 2 ]). Muscle-induced biomechanical variables for this cohort were reported and had been calculated from patient-specific musculoskeletal models, motion data, and MRI. Biomechanical variables included joint reaction forces, acetabular edge loads, hip center lateralization, and gluteus medius muscle moment arm lengths. PROMs included the Hip Disability and Osteoarthritis Outcome Score (HOOS), the WOMAC, International Hip Outcome Tool-12, National Institutes of Health Patient-Reported Outcome Measure Information System (PROMIS) Pain Interference and Physical Function subscales, and University of California Los Angeles activity scale. Associations between PROMs and biomechanical variables were tested using Spearman rank-order correlations and corrected for multiple comparisons using the Benjamini-Yekutieli method. For this study, associations between variables were considered to exist when correlations were statistically significant (p < 0.05) and were either strong (ρ ≥ 0.60) or moderate (ρ = 0.40 to 0.59).

RESULTS

Acetabular edge load impulses (the cumulative acetabular edge load across the gait cycle), medially directed joint reaction forces, and hip center lateralization most commonly demonstrated moderate or strong associations with PROMs. The strongest associations were a negative correlation between acetabular edge load impulse on the superior acetabulum and the HOOS function in daily living subscale (ρ = -0.63; p = 0.001), followed by a negative correlation between hip center lateralization and the HOOS pain subscale (ρ = -0.6; p = 0.003), and a positive correlation between hip center lateralization and the PROMIS pain subscale (ρ = 0.62; p = 0.002). The University of California Los Angeles activity scale was the only PROM that did not demonstrate associations with any biomechanical variable. All PROMs, aside from the University of California Los Angeles activity scale, were associated with one another. Although most of the biomechanical variables were associated with one another, these relationships were not as consistent as those among PROMs.

CONCLUSION

The associations with PROMs detected in the current study suggest that muscle-induced biomechanics may have wide-reaching effects not only on loads within the hip, but also on patients' perceptions of their health and function. As the treatment of DDH evolves, patient-specific joint preservation strategies may benefit from targeting the underlying causes of biomechanical outcomes associated with PROMs.

LEVEL OF EVIDENCE

Level III, prognostic study.

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.

Influence of different calibration methods on surface electromyography-informed musculoskeletal models with few input signals

BACKGROUND

Although model personalization is critical when assessing individuals with morphological or neurological abnormalities, or even non-disabled subjects, its translation into routine clinical settings is hampered by the cumbersomeness of experimental data acquisition and lack of resources, which are linked to high costs and long processing pipelines. Quantifying the impact of neglecting subject-specific information in simulations that aim to estimate muscle forces with surface electromyography informed modeling approaches, can address their potential in relevant clinical questions. The present study investigates how different methods to fine-tune subject-specific neuromuscular parameters, reducing the number of electromyography input data, could affect the estimation of the unmeasured excitations and the musculotendon forces.

METHODS

Three-dimensional motion analysis was performed on 8 non-disabled adult subjects and 13 electromyographic signals captured. Four neuromusculoskeletal models were created for 8 participants: a reference model driven by a large set of sEMG signals; two models informed by four electromyographic signals but calibrated in different fashions; a model based on static optimization.

FINDINGS

The electromyography-informed models better predicted experimental excitations, including the unmeasured ones. The model based on static optimization obtained less reliable predictions of the experimental data. When comparing the different reduced models, no major differences were observed, suggesting that the less complex model may suffice for predicting muscle forces with a small set of input in clinical gait analysis tasks.

INTERPRETATION

Quantitative model performance evaluation in different conditions provides an objective indication of which method yields the most accurate prediction when a small set of electromyographic recordings is available.

Optimal degrees of freedom of the lower extremities for human walking and running

Determining the degrees of freedom (DOF) of the linked rigid-body model, representing a multi-body motion of the human lower extremity, is one of the most important procedures in locomotion analysis. However, a trade-off exists between the quality of data fitting and the generalizability of the model. This study aimed to determine the optimal DOF of the model for the lower extremities that balance the goodness-of-fit and generalizability of the model during walking and running using Akaike's information criterion (AIC). Empirically obtained kinematic data for the lower extremities during walking and running were fitted by models with 9, 18, or 22 DOF. The relative quality of these models was assessed using their bias-corrected AIC (cAIC) value. A significant simple main effect of the model was found on the cAIC value for both walking and running conditions. Pairwise comparisons revealed that the cAIC value of the 18-DOF model was significantly smaller than that of the 9-DOF (walking: p < 0.001, running: p = 0.010) and 22-DOF (walking: p < 0.001, running: p < 0.001) models. These findings suggest that the 18-DOF model is optimal for representing the lower extremities during walking and running, in terms of goodness-of-fit and generalizability.

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.

Comparison of biomechanical analysis results using different musculoskeletal models for children with cerebral palsy

Introduction: Musculoskeletal model-based simulations have gained popularity as a tool for analyzing human movement biomechanics. However, when examining the same gait, different models with varying anatomical data and assumptions may produce inconsistent biomechanical results. This inconsistency is particularly relevant for children with cerebral palsy, who often exhibit multiple pathological gait patterns that can impact model outputs. Methods: The aim of this study was to investigate the effect of selecting musculoskeletal models on the biomechanical analysis results in children with cerebral palsy. Gait data were collected from multiple participants at slow, medium, and fast velocities. Joint kinematics, joint dynamics, and muscle activation were calculated using six popular musculoskeletal models within a biomechanical simulation environment. Results: The degree of inconsistency, measured as the root-mean-square deviation, in kinematic and kinetic results produced by the different models ranged from 4% to 40% joint motion range and 0%-28% joint moment range, respectively. The correlation between the results of the different models (both kinematic and kinetic) was good (R> 0.85, P < 0.01), with a stronger correlation observed in the kinetic results. Four of the six models showed a positive correlation between the simulated muscle activation of rectus femoris and the surface EMG, while all models exhibited a positive correlation between the activation of medial gastrocnemius and the surface EMG (P < 0.01). Discussion: These results provide insights into the consistency of model results, factors influencing consistency, characteristics of each model's outputs, mechanisms underlying these characteristics, and feasible applications for each model. By elucidating the impact of model selection on biomechanical analysis outcomes, this study advances the field's understanding of musculoskeletal modeling and its implications for clinical gait analysis model decision-making in children with cerebral palsy.

Medial and lateral knee contact forces during walking, stair ascent and stair descent are more affected by contact locations than tibiofemoral alignment in knee osteoarthritis patients with varus malalignment

Introduction: Knee OA progression is related to medial knee contact forces, which can be altered by anatomical parameters of tibiofemoral alignment and contact point locations. There is limited and controversial literature on medial-lateral force distribution and the effect of anatomical parameters, especially in motor activities different from walking. We analyzed the effect of tibiofemoral alignment and contact point locations on knee contact forces, and the medial-lateral force distribution in knee OA subjects with varus malalignment during walking, stair ascending and stair descending. Methods: Fifty-one knee OA subjects with varus malalignment underwent weight-bearing radiographs and motion capture during walking, stair ascending and stair descending. We created a set of four musculoskeletal models per subject with increasing level of personalization, and calculated medial and lateral knee contact forces. To analyze the effect of the anatomical parameters, statistically-significant differences in knee contact forces among models were evaluated. Then, to analyze the force distribution, the medial-to-total contact force ratios were calculated from the fully-informed models. In addition, a multiple regression analysis was performed to evaluate correlations between forces and anatomical parameters. Results: The anatomical parameters significantly affected the knee contact forces. However, the contact points decreased medial forces and increased lateral forces and led to more marked variations compared to tibiofemoral alignment, which produced an opposite effect. The forces were less medially-distributed during stair negotiation, with medial-to-total ratios below 50% at force peaks. The anatomical parameters explained 30%-67% of the variability in the knee forces, where the medial contact points were the best predictors of medial contact forces. Discussion: Including personalized locations of contact points is crucial when analyzing knee contact forces in subjects with varus malalignment, and especially the medial contact points have a major effect on the forces rather than tibiofemoral alignment. Remarkably, the medial-lateral force distribution depends on the motor activity, where stair ascending and descending show increased lateral forces that lead to less medially-distributed loads compared to walking.

Increasing Step Frequency Reduces Patellofemoral Joint Stress and Patellar Tendon Force Impulse More at Low Running Speed

PURPOSE

Patellofemoral pain syndrome and patellar tendinopathy are important running-related overuse injuries. This study investigated the interaction of running speed and step frequency alterations on peak and cumulative patellofemoral joint stress (PFJS) and patellar tendon force (PTF) parameters.

METHODS

Twelve healthy individuals completed an incremental running speed protocol on a treadmill at habitual, increased and decreased step frequency. Peak PFJS and PTF, peak rate of PFJS and PTF development, and PFJS and PTF impulse per kilometer (km) were calculated using musculoskeletal modeling.

RESULTS

With increasing running speed, peak PFJS ( P < 0.001) and PTF ( P < 0.001) and peak rate of PFJS ( P < 0.001) and PTF ( P < 0.001) development increased, whereas PFJS ( P < 0.001) and PTF ( P < 0.001) impulse per km decreased. While increasing step frequency by 10%, the peak PFJS ( P < 0.001) and PTF ( P < 0.001) and the PFJS ( P < 0.001) and PTF ( P < 0.001) impulse per kilometer decreased. No significant effect of step frequency alteration was found for the peak rate of PFJS ( P = 0.008) and PTF ( P = 0.213) development. A significant interaction effect was found for PFJS ( P < 0.001) and PTF ( P < 0.001) impulse per km, suggesting that step frequency alteration was more effective at low running speed.

CONCLUSIONS

The effectiveness of step frequency alteration on PFJS and PTF impulse per km is dependent on the running speed. With regard to peak PFJS and PTF, step frequency alteration is equally effective at low and high running speeds. Step frequency alteration was not effective for peak rate of PFJS and PTF development. These findings can assist the optimization of patellofemoral joint and patellar tendon load management strategies.

Validation of a musculoskeletal model to investigate hip joint mechanics in response to dynamic multiplanar tasks

Existing hip-focused musculoskeletal (MSK) models are limited by the hip range of motion, hip musculature detail, or have only been qualitatively validated. The purposes of this study were to: i) modify the existing 2396Hip MSK model to simulate dynamic tasks with multiplanar hip joint motion; and ii) validate the modified MSK model quantitatively against experimental data. Experimental data was collected from five healthy adults (age = 25 [6] years, two females) during eight movement tasks. The motion and ground reaction force data were input into the MSK modeling software OpenSim to calculate muscle activations and hip contact forces (HCFs). The HCFs were compared to experimental HCFs previously measured in total hip arthroplasty (THA) patients using instrumented hip prostheses. A gait simulation was performed using data from one THA patient to directly assess the model's accuracy in estimating HCFs. The young adults' modeled and experimental muscle activations for seven muscles were compared using a cross-correlation function. The model only overestimated the peak resultant HCFs by 0.06-0.08 N/BW compared to the experimentally measured HCFs of the THA patient. The young adults' HCFs were over two standard deviations higher than previously measured in the THA patients, which is likely a result of different movement patterns. The correlation coefficients indicated strong correlations between experimental and modeled muscle activations in 50 of the 56 comparisons. The results of this study suggest the new MSK model is an appropriate method to quantify HCFs and muscle activations in response to dynamic, multiplanar tasks among young, healthy adults.

Do Interlimb Knee Joint Loading Asymmetries Persist throughout Stance during Uphill Walking Following Total Knee Arthroplasty?

The purpose of this study was to determine differences in total (TCF), medial compartment (MCF), and lateral compartment (LCF) tibiofemoral joint compressive forces and related muscle forces between replaced and non-replaced limbs during level and uphill walking at an incline of 10°. A musculoskeletal modeling and simulation approach using static optimization was used to determine the muscle forces and TCF, MCF, and LCF for 25 patients with primary TKA. A statistical parametric mapping repeated-measures ANOVA was conducted on knee compressive forces and muscle forces using statistical parametric mapping. Greater TCF, MCF, and LCF values were observed throughout the loading response, mid-stance, and late stance during uphill walking. During level walking, knee extensor muscle forces were greater throughout the first 50% of the stance during level walking, yet greater during uphill walking during the last 50% of the stance. Conversely, knee flexor muscle forces were greater through the loading response and push-off phases of the stance. No between-limb differences were observed for compressive or muscle forces, suggesting that uphill walking may promote a more balanced loading of replaced and non-replaced limbs. Additionally, patients with TKA appear to rely on the hamstrings muscle group during the late stance for knee joint control, thus supporting uphill walking as an effective exercise modality to improve posterior chain muscle strength.

Design and Implementation of a Portable Knee Actuator for the Improvement of Crouch Gait in Children with Cerebral Palsy

Common manifestation of spastic Cerebral Palsy (CP) are abnormal gait pathologies. These conditions require greater energy expenditure to successfully ambulate and are linked with significant deterioration in joint health and childhood musculoskeletal development. Crouch gait presents with knee hyperflexion throughout stance due to extensor muscle weakness and spasticity in flexor muscles stemming from neurological damage. The goal of this study was to develop a wearable cable-driven robotic system that applies controlled perturbation to the knee joint during overground walking in children with CP. Two children with spastic CP were recruited in this pilot study. They were tested in two conditions, i.e., applying knee resistance vs. knee assistance during overground walking. Kinematic and EMG data were recorded during overground walking. Data indicated that it was feasible to apply controlled knee perturbation torque during overground walking in children with crouch and preliminary results showed an improvement in crouch gait pattern in children with CP after one session of walking with the robotic system.Clinical Relevance- This study might have a potential clinical significance modifying neuromuscular control of CP patients with Crouch Gait.

The non-invasive evaluation technique of patellofemoral joint stress: a systematic literature review

Introduction: Patellofemoral joint stress (PFJS) is an important parameter for understanding the mechanism of patellofemoral joint pain, preventing patellofemoral joint injury, and evaluating the therapeutic efficacy of PFP rehabilitation programs. The purpose of this systematic review was to identify and categorize the non-invasive technique to evaluate the PFJS. Methods: Literature searches were conducted from January 2000 to October 2022 in electronic databases, namely, PubMed, Web of Science, and EBSCO (Medline, SPORTDiscus). This review includes studies that evaluated the patellofemoral joint reaction force (PJRF) or PFJS, with participants including both healthy individuals and those with patellofemoral joint pain, as well as cadavers with no organic changes. The study design includes cross-sectional studies, case-control studies, and randomized controlled trials. The JBI quality appraisal criteria tool was used to assess the risk of bias in the included studies. Results: In total, 5016 articles were identified in the database research and the citation network, and 69 studies were included in the review. Discussion: Researchers are still working to improve the accuracy of evaluation for PFJS by using a personalized model and optimizing quadriceps muscle strength calculations. In theory, the evaluation method of combining advanced computational and biplane fluoroscopy techniques has high accuracy in evaluating PFJS. The method should be further developed to establish the "gold standard" for PFJS evaluation. In practical applications, selecting appropriate methods and approaches based on theoretical considerations and ecological validity is essential.

Muscle force-length dynamics during walking over obstacles indicates delayed recovery and a shift towards more 'strut-like' function in birds with proprioceptive deficit

Recent studies of in vivo muscle function in guinea fowl revealed that distal leg muscles rapidly modulate force and work to stabilize running in uneven terrain. Previous studies focused on running only, and it remains unclear how muscular mechanisms for stability differ between walking and running. Here, we investigated in vivo function of the lateral gastrocnemius (LG) during walking over obstacles. We compared muscle function in birds with intact (iLG) versus self-reinnervated LG (rLG). Self-reinnervation results in proprioceptive feedback deficit due to loss of monosynaptic stretch reflex. We tested the hypothesis that proprioceptive deficit results in decreased modulation of EMG activity in response to obstacle contact, and a delayed obstacle recovery compared with that for iLG. We found that total myoelectric intensity (Etot) of iLG increased by 68% in obstacle strides (S 0) compared with level terrain, suggesting a substantial reflex-mediated response. In contrast, Etot of rLG increased by 31% in S 0 strides compared with level walking, but also increased by 43% in the first post-obstacle (S +1) stride. In iLG, muscle force and work differed significantly from level walking only in the S 0 stride, indicating a single-stride recovery. In rLG, force increased in S 0, S +1 and S +2 compared with level walking, indicating three-stride obstacle recovery. Interestingly, rLG showed little variation in work output and shortening velocity in obstacle terrain, indicating a shift towards near-isometric strut-like function. Reinnervated birds also adopted a more crouched posture across level and obstacle terrains compared with intact birds. These findings suggest gait-specific control mechanisms in walking and running.

Ten steps to becoming a musculoskeletal simulation expert: A half-century of progress and outlook for the future

Over the past half-century, musculoskeletal simulations have deepened our knowledge of human and animal movement. This article outlines ten steps to becoming a musculoskeletal simulation expert so you can contribute to the next half-century of technical innovation and scientific discovery. We advocate looking to the past, present, and future to harness the power of simulations that seek to understand and improve mobility. Instead of presenting a comprehensive literature review, we articulate a set of ideas intended to help researchers use simulations effectively and responsibly by understanding the work on which today's musculoskeletal simulations are built, following established modeling and simulation principles, and branching out in new directions.

Unilateral transfemoral osseointegrated prostheses improve joint loading during walking

People with unilateral transfemoral amputation using socket prostheses are at increased risk for developing osteoarthritis in both the residual hip and intact lower-limb joints. Osseointegrated prostheses are a surgical alternative to socket prostheses that directly attach to the residual femur via a bone-anchored implant, however their multi-joint loading effect is largely unknown. Our objective was to establish how osseointegrated prostheses influence joint loading during walking. Motion capture data (kinematics, ground reaction forces) were collected from 12 participants at baseline, with socket prostheses, and 12-months after prosthesis osseointegration during overground walking at self-selected speeds. Subject-specific musculoskeletal models were developed at each timepoint relative to osseointegration. Internal joint moments were calculated using inverse dynamics, muscle and joint reaction forces (JRFs) were estimated with static optimization. Changes in internal joint moments, JRFs, and joint loading-symmetry were compared using statistical parametric mapping (p≤ 0.05) before and after osseointegration. Amputated limb hip flexion moments and anterior JRFs decreased during terminal stance (p = 0.002, <0.001; respectively), while amputated limb hip abduction moments increased during mid-stance (p < 0.001), amputated hip rotation moment changed from internal to external throughout early stance (p < 0.001). Intact limb hip extension and knee flexion moments (p = 0.028, 0.032; respectively), superior and resultant knee JRFs (p = 0.046, 0.049; respectively) decreased during the loading response following prosthesis osseointegration. These results may indicate that the direct loading transmission of these novel prostheses create a more typical mechanical environment in bilateral joints, which is comparable with loading observed in able-bodied individuals and could decrease the risk of development or progression of osteoarthritis.

Monitoring Knee Contact Force with Force-Sensing Insoles

Numerous applications exist for monitoring knee contact force (KCF) throughout activities of daily living. However, the ability to estimate these forces is restricted to a laboratory setting. The purposes of this study are to develop KCF metric estimation models and explore the feasibility of monitoring KCF metrics via surrogate measures derived from force-sensing insole data. Nine healthy subjects (3F, age 27 ± 5 years, mass 74.8 ± 11.8 kg, height 1.7 ± 0.08 m) walked at multiple speeds (0.8-1.6 m/s) on an instrumented treadmill. Thirteen insole force features were calculated as potential predictors of peak KCF and KCF impulse per step, estimated with musculoskeletal modeling. The error was calculated with median symmetric accuracy. Pearson product-moment correlation coefficients defined the relationship between variables. Models develop per-limb demonstrated lower prediction error than those developed per-subject (KCF impulse: 2.2% vs 3.4%; peak KCF: 3.50% vs. 6.5%, respectively). Many insole features are moderately to strongly associated with peak KCF, but not KCF impulse across the group. We present methods to directly estimate and monitor changes in KCF using instrumented insoles. Our results carry promising implications for internal tissue loads monitoring outside of a laboratory with wearable sensors.

The influence of induced gait asymmetry on joint reaction forces

Chronic injury- or disease-induced joint impairments result in asymmetric gait deviations that may precipitate changes in joint loading associated with pain and osteoarthritis. Understanding the impact of gait deviations on joint reaction forces (JRFs) is challenging because of concurrent neurological and/or anatomical changes and because measuring JRFs requires medically invasive instrumented implants. Instead, we investigated the impact of joint motion limitations and induced asymmetry on JRFs by simulating data recorded as 8 unimpaired participants walked with bracing to unilaterally and bilaterally restrict ankle, knee, and simultaneous ankle + knee motion. Personalized models, calculated kinematics, and ground reaction forces (GRFs) were input into a computed muscle control tool to determine lower limb JRFs and simulated muscle activations guided by electromyography-driven timing constraints. Unilateral knee restriction increased GRF peak and loading rate ipsilaterally but peak values decreased contralaterally when compared to walking without joint restriction. GRF peak and loading rate increased with bilateral restriction compared to the contralateral limb of unilaterally restricted conditions. Despite changes in GRFs, JRFs were relatively unchanged due to reduced muscle forces during loading response. Thus, while joint restriction results in increased limb loading, reductions in muscle forces counteract changes in limb loading such that JRFs were relatively unchanged.

Can static optimization detect changes in peak medial knee contact forces induced by gait modifications?

Medial knee contact force (MCF) is related to the pathomechanics of medial knee osteoarthritis. However, MCF cannot be directly measured in the native knee, making it difficult for therapeutic gait modifications to target this metric. Static optimization, a musculoskeletal simulation technique, can estimate MCF, but there has been little work validating its ability to detect changes in MCF induced by gait modifications. In this study, we quantified the error in MCF estimates from static optimization compared to measurements from instrumented knee replacements during normal walking and seven different gait modifications. We then identified minimum magnitudes of simulated MCF changes for which static optimization correctly identified the direction of change (i.e., whether MCF increased or decreased) at least 70% of the time. A full-body musculoskeletal model with a multi-compartment knee and static optimization was used to estimate MCF. Simulations were evaluated using experimental data from three subjects with instrumented knee replacements who walked with various gait modifications for a total of 115 steps. Static optimization underpredicted the first peak (mean absolute error = 0.16 bodyweights) and overpredicted the second peak (mean absolute error = 0.31 bodyweights) of MCF. Average root mean square error in MCF over stance phase was 0.32 bodyweights. Static optimization detected the direction of change with at least 70% accuracy for early-stance reductions, late-stance reductions, and early-stance increases in peak MCF of at least 0.10 bodyweights. These results suggest that a static optimization approach accurately detects the direction of change in early-stance medial knee loading, potentially making it a valuable tool for evaluating the biomechanical efficacy of gait modifications for knee osteoarthritis.

Simulated Tibiofemoral Joint Reaction Forces for Three Previously Studied Gait Modifications in Healthy Controls

Gait modifications, such as lateral trunk lean (LTL), medial knee thrust (MKT), and toe-in gait (TIG), are frequently investigated interventions used to slow the progression of knee osteoarthritis. The Lerner knee model was developed to estimate the tibiofemoral joint reaction forces (JRF) in the medial and lateral compartments during gait. These models may be useful for estimating the effects on the JRF in the knee as a result of gait modifications. We hypothesized that all gait modifications would decrease the JRF compared to normal gait. Twenty healthy individuals volunteered for this study (26.7 ± 4.7 years, 1.75 ± 0.1 m, 73.4 ± 12.4 kg). Ten trials were collected for normal gait as well as for the three gait modifications: LTL, MKT, and TIG. The data were used to estimate the JRF in the first and second peaks for the medial and lateral compartments of the knee via opensim using the Lerner knee model. No significant difference from baseline was found for the first peak in the medial compartment. There was a decrease in JRF in the medial compartment during the loading phase of gait for TIG (6.6%) and LTL (4.9%) and an increasing JRF for MKT (2.6%). but none was statistically significant. A significant increase from baseline was found for TIG (5.8%) in the medial second peak. We found a large variation in individual responses to gait interventions, which may help explain the lack of statistically significant results. Possible factors influencing these wide ranges of responses to gait modifications include static alignment and the impacts of variation in muscle coordination strategies used, by participants, to implement gait modifications.

The impact of the parameters of the constitutive model on the distribution of strain in the femoral head

The rapid spread of the finite element method has caused that it has become, among other methods, the standard tool for pre-clinical estimates of bone properties. This paper presents an application of this method for the calculation and prediction of strain and stress fields in the femoral head. The aim of the work is to study the influence of the considered anisotropy and heterogeneity of the modeled bone on the mechanical fields during a typical gait cycle. Three material models were tested with different properties of porous bone carried out in literature: a homogeneous isotropic model, a heterogeneous isotropic model, and a heterogeneous anisotropic model. In three cases studied, the elastic properties of the bone were determined basing on the Zysset-Curnier approach. The tensor of elastic constants defining the local properties of porous bone is correlated with a local porosity and a second order fabric tensor describing the bone microstructure. In the calculations, a model of the femoral head generated from high-resolution tomographic scans was used. Experimental data were drawn from publicly available database "Osteoporotic Virtual Physiological Human Project." To realistically reflect the load on the femoral head, main muscles were considered, and their contraction forces were determined based on inverse kinematics. For this purpose, the results from OpenSim packet were used. The simulations demonstrated that differences between the results predicted by these material models are significant. Only the anisotropic model allowed for the plausible distribution of stresses along the main trabecular groups. The outcomes also showed that the precise evaluation of the mechanical fields is critical in the context of bone tissue remodeling under mechanical stimulations.

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.

Neuromusculoskeletal model calibration accounts for differences in electromechanical delay and maximum isometric muscle force

Electromechanical delay (EMD) and maximum isometric muscle force (F_o^M) are important parameters for joint contact force calculation with EMG-informed neuromusculoskeletal (NMS) models. These parameters can vary between tasks (EMD) and individuals (EMD and F_o^M), making it challenging to establish representative values. One promising approach is to personalise candidate parameters to the participant (e.g., F_o^M by regression equation) and then adjust all parameters within a calibration (i.e., numerical optimisation) to minimise error between corresponding pairs of experimental measures and model-predicted values. The purpose of this study was to determine whether calibration of an NMS model resulted in consistent joint contact forces, regardless of EMD value or personalisation of F_o^M. Hip, knee, and ankle contact forces were predicted for 28 participants using EMG-informed NMS models. Differences in joint contact forces with EMD were examined in six models, calibrated with EMD from 15 to 110 ms. Differences in joint contact forces with personalisation of F_o^M were examined in two models, both calibrated with the same initial EMD (50 ms), one with generic and one with personalised values for F_o^M. For all models, joint contact force peaks during the first and second halves of stance were extracted and compared using a repeated-measures analysis of variance. Calibrated models with EMD set between 35 and 70 ms produced similar magnitude and timing of peak joint contact forces. Compared with generic values, personalising and then calibrating F_o^M resulted in comparable peak contact forces at hip, but not knee or ankle, while also producing muscle-specific tensions similar to reported literature. Overall, EMD between 35 and 70 ms and personalised initial values of F_o^M before calibration are advised for EMG-informed NMS modelling.

Influence of muscle loading on early-stage bone fracture healing

Designing weight-bearing exercises for patients with lower-limb bone fractures is challenging and requires a systematic approach that accounts for patient-specific loading conditions. However, 'trial-and-error' approaches are commonplace in clinical settings due to the lack of a fundamental understanding of the effect of weight-bearing exercises on the bone healing process. Whilst computational modelling has the potential to assist clinicians in designing effective patient-specific weight-bearing exercises, current models do not explicitly account for the effects of muscle loading, which could play an important role in mediating the mechanical microenvironment of a fracture site. We combined a fracture healing model involving a tibial fracture stabilised with a locking compression plate (LCP) with a detailed musculoskeletal model of the lower limb to determine interfragmentary strains in the vicinity of the fracture site during both full weight-bearing (100% body weight) and partial weight-bearing (50% body weight) standing. We found that muscle loading significantly altered model predictions of interfragmentary strains. For a fractured bone with a standard LCP configuration (bone-plate distance = 2 mm, working length = 30 mm) subject to full weight-bearing, the predicted strains at the near and far cortices were 23% and 11% higher when muscle loading was included compared to the case when muscle loading was omitted. The knee and ankle muscles accounted for 38% of the contact force exerted at the knee joint during quiet standing and contributed significantly to the strains calculated at the fracture site. Thus, models of bone fracture healing ought to account explicitly for the effects of muscle loading. Furthermore, the study indicated that LCP configuration parameters play a crucial role in influencing the fracture site microenvironment. The results highlighted the dominance of working length over bone-plate distance in controlling the flexibility of fracture sites stabilised with LCP devices.

Tibiofemoral compressive force during downhill walking in patients with primary total knee arthroplasty: A statistical parametric mapping approach

BACKGROUND

Downhill walking is a necessary part of daily life and an effective activity in post-operative rehabilitation following total knee arthroplasty. The purpose of this study was to determine differences in the behavior of total, medial, and lateral tibiofemoral compressive forces as well as knee extensor and flexor muscle forces between different limbs of patients with total knee arthroplasty (replaced, non-replaced) during downhill and level walking.

METHODS

Musculoskeletal modeling and simulation were implemented to determine muscle forces and tibiofemoral compressive forces in 25 patients with total knee arthroplasty. A 2 × 2 [Limb (replaced, non-replaced) × Slope (0°, 10°)] Statistical parametric mapping repeated measures analysis of variance was conducted on selected variables.

FINDINGS

Statistical parametric mapping did not identify any between-limb differences for compressive or muscle forces. Differences in joint compressive and muscle forces persisted throughout different intervals of stance-phase between level and downhill walking. Knee extensor muscle forces were distinctly greater during level walking for nearly all of stance phase. Knee flexor muscle force was greater during downhill walking for >60% of stance. Statistical parametric mapping did identify regions of significance between level and downhill walking that coincided temporally (near loading response and push off) with peak joint moment and joint compressive forces traditionally reported using discrete variable analyses.

INTERPRETATION

Downhill walking may be a safe and useful rehabilitation tool for post-knee arthroplasty rehabilitation that will not disproportionally load either the replaced or the non-replaced joint and where the quadriceps muscles can be strengthened during a gait-specific task.

Elbow joint biomechanics during ADL focusing on total elbow arthroplasty - a scoping review

BACKGROUND

Overloading is hypothesized to be one of the failure mechanisms following total elbow arthroplasty (TEA). It is unclear whether the current post-operative loading instruction is compliant with reported failure mechanisms. Aim is therefore to evaluate the elbow joint load during activities of daily living (ADL) and compare these loads with reported failure limits from retrieval and finite element studies.

METHODS

A scoping review of studies until 23 November 2021 investigating elbow joint load during ADL were identified by searching PubMed/Medline and Web of Science. Studies were eligible when: (1) reporting on the elbow joint load in native elbows or elbows with an elbow arthroplasty in adults; (2) full-text article was available.

RESULTS

Twenty-eight studies with a total of 256 participants were included. Methodological quality was low in 3, moderate in 22 and high in 3 studies. Studies were categorized as 1) close to the body and 2) further away from the body. Tasks were then subdivided into: 1) cyclic flexion/extension, 2) push-up, 3) reaching, 4) self-care, 5) work. Mean flexion-extension joint load was 17 Nm, mean varus-valgus joint load 9 Nm, mean pronation-supination joint load 8 Nm and mean bone-on-bone contact force 337 N.

CONCLUSION

The results of our scoping review give a first overview of the current knowledge on elbow joint loads during ADL. Surprisingly, the current literature is not sufficient to formulate a postoperative instruction for elbow joint loading, which is compliant with failure limits of the prosthesis. In addition, our current instruction does not appear to be evidence-based. Our recommendations offer a starting point to assist clinicians in providing informed decisions about post-operative instructions for their patients.

Effects of Obesity on Medial Tibiofemoral Cartilage Mechanics in Females-An Exploration Using Musculoskeletal Simulation and Probabilistic Cartilage Failure Modelling

This study examined the effects of obesity on cartilage mechanics and longitudinal failure probability at the medial tibiofemoral compartment, using combined musculoskeletal simulation and probabilistic failure modelling approaches. The current investigation examined twenty obese females (BMI > 30.0 kg/m²) and 20 healthy weight (BMI < 25.0 kg/m²) females. Walking kinematics were obtained via an 8-camera optoelectric system, and a force plate was used to collect ground reaction forces. Musculoskeletal simulation and probabilistic failure modelling were utilized to explore medial tibiofemoral forces and cartilage probability. Comparisons between groups were undertaken using linear mixed-effects models. Net peak cartilage forces, stress and strain were significantly larger in the obese group (force = 2013.92 N, stress = 3.03 MPa & strain = 0.25), compared to health weight (force = 1493.21 N, stress 2.26 MPa & strain = 0.19). In addition, medial tibiofemoral cartilage failure probability was also significantly larger in the obese group (42.98%) compared to healthy weight (11.63%). The findings from the current investigation show that obesity has a profoundly negative influence on longitudinal medial knee cartilage health and strongly advocates for the implementation of effective weight management programs into long-term musculoskeletal management strategies.

The biomechanical fingerprint of hip and knee osteoarthritis patients during activities of daily living

BACKGROUND

Osteoarthritis is a highly prevalent disease affecting the hip and knee joint and is characterized by load-mediated pain and decreased quality of life. Dependent on involved joint, patients present antalgic movement compensations, aiming to decrease loading on the involved joint. However, the associated alterations in mechanical loading of the ipsi- and contra-lateral lower limb joints, are less documented. Here, we documented the biomechanical fingerprint of end-stage hip and knee osteoarthritis patients in terms of ipsilateral and contralateral hip and knee loading during walking and stair ambulation.

METHODS

Three-dimensional motion-analysis was performed in 20 hip, 18 knee osteoarthritis patients and 12 controls during level walking and stair ambulation. Joint contact forces were calculated using a standard musculoskeletal modelling workflow in Opensim. Involved and contralateral hip and knee joint loading was compared against healthy controls using independent t-tests (p < 0.05).

FINDINGS

Both hip and knee cohorts significantly decreased loading of the involved joint during gait and stair ambulation. Hip osteoarthritis patients presented no signs of ipsilateral knee nor contralateral leg overloading, during walking and stair ascending. However, knee osteoarthritis patients significantly increased loading at the ipsilateral hip, and contralateral hip and knee joints during stair ambulation compared to controls.

INTERPRETATION

The biomechanical fingerprint in knee and hip osteoarthritis patients confirmed antalgic movement strategies to unload the involved leg during gait. Only during stair ambulation in knee osteoarthritis patients, movement adaptations were confirmed that induced unbalanced intra- and inter-limb loading conditions, which are known risk factors for secondary osteoarthritis.

PATELLAR FRACTURE IN ANTERIOR CRUCIATE LIGAMENT RECONSTRUCTION: IN VITRO ANALYSIS

Objective

To determine, by biomechanical analysis, safe patellar cut limits in anterior cruciate ligament (ACL) reconstruction that minimize fracture risks.

Methods

From three-dimensional reconstruction, triangular cuts were made in the patella, with a depth of 6.5 mm and variable width and length (10 to 20 mm and 8 to 12 mm, respectively, both with an interval of 1 mm). The combinations of cuts constituted 55 models for tests, with five variations in width and 11 variations in length, tested with the finite element method (FEM).

Results

The mean of the localized principal maximum (traction force) values was 4.36 Pa (SD 0.87 ± 0.76) and the localized principal minimum (compression force) was -4.33 Pa (SD 1.05 ± 1.11). Comparing width and length to the tension force of the values of the main maximum, we found statistical significance from 11 mm for width and 13 mm for length.

Conclusion

In ACL reconstruction, the removal of the patellar bone fragment is safe for fragments smaller than 11 mm in width and 13 mm in length, which corresponds to 24% of the width and 28% of the length of the patella used. Level of Evidence II, Comparative Prospective Study.

Independently ambulatory children with spina bifida experience near-typical knee and ankle joint moments and forces during walking

BACKGROUND

Spina bifida, a neurological defect, can result in lower-limb muscle weakness. Altered ambulation and reduced musculoskeletal loading can yield decreased bone strength in individuals with spina bifida, yet individuals who remain ambulatory can exhibit normal bone outcomes.

RESEARCH QUESTION

During walking, how do lower-limb joint kinematics and moments and tibial forces in independently ambulatory children with spina bifida differ from those of children with typical development?

METHODS

We retrospectively analyzed data from 16 independently ambulatory children with spina bifida and 16 children with typical development and confirmed that tibial bone strength was similar between the two groups. Plantar flexor muscle strength was measured by manual muscle testing, and 14 of the children with spina bifida wore activity monitors for an average of 5 days. We estimated tibial forces at the knee and ankle using motion capture data and musculoskeletal simulations. We used Statistical Parametric Mapping t-tests to compare lower-limb joint kinematic and kinetic waveforms between the groups with spina bifida and typical development. Within the group with spina bifida, we examined relationships between plantar flexor muscle strength and peak tibial forces by calculating Spearman correlations.

RESULTS

Activity monitors from the children with spina bifida reported typical daily steps (9656 [SD 3095]). Despite slower walking speeds (p = 0.004) and altered lower-body kinematics (p < 0.001), children with spina bifida had knee and ankle joint moments and forces similar to those of children with typical development, with no detectable differences during stance. Plantar flexor muscle weakness was associated with increased compressive knee force (p = 0.002) and shear ankle force (p = 0.009).

SIGNIFICANCE

High-functioning, independently ambulatory children with spina bifida exhibited near-typical tibial bone strength and near-typical step counts and tibial load magnitudes. Our results suggest that the tibial forces in this group are of sufficient magnitudes to support the development of normal tibial bone strength.

An investigation into the hammer toe effects on the lower extremity mechanics and plantar fascia tension: A case for a vicious cycle and progressive damage

Hammer toes are one of the common deformities of the forefoot that can lead to compensatory changes during walking in individuals with this condition. Predicting the adverse effects of tissue damage on the performance of other limbs is very important in the prevention of progressive damage. Finite element (FE) and musculoskeletal modeling can be helpful by allowing such effects to be studied in a way where the internal stresses in the tissue could be investigated. Hence, this study aims to investigate the effects of the hammer toe deformity on the lower extremity, especially on the plantar fascia functions. To compare the joint reactions of the hammer toe foot (HTF) and healthy foot (HF), two musculoskeletal models (MSM) of the feet of a healthy individual and that of a participant with hammer toe foot were developed based on gait analysis. A previously validated 3D finite element model which was constructed using Magnetic Resonance Imaging (MRI) of the diabetic participant with the hammer toe deformity was processed at five different events during the stance phase of gait. It was found that the hammer toe deformity makes dorsiflexion of the toes and the windlass mechanism less effective during walking. Specifically, the FE analysis results showed that plantar fascia (PF) in HTF compared to HF played a less dominant role in load bearing with both medial and lateral parts of PF loaded. Also, the results indicated that the stored elastic energy in PF was less in HTF than the HF, which can indicate a higher metabolic cost during walking. Internal stress distribution shows that the majority of ground reaction forces are transmitted through the lateral metatarsals in hammer toe foot, and the probability of fifth metatarsal fracture and also progressive deformity was subsequently increased. The MSM results showed that the joint reaction forces and moments in the hammer toe foot have deviated from normal, where the metatarsophalangeal joint reactions in the hammer toe were less than the values in the healthy foot. This can indicate a vicious cycle of foot deformity, leading to changes in body weight force transmission line, and deviation of joint reactions and plantar fascia function from normal. These in turn lead to increased internal stress concentration, which in turn lead to further foot deformities. This vicious cycle cause progressive damage and can lead to an increase in the risk of ulceration in the diabetic foot.

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.