Partitioning of knee joint internal forces in gait is dictated by the knee adduction angle and not by the knee adduction moment

Effects of overweight and obesity on lower limb walking characteristics from joint kinematics to muscle activations

BACKGROUND

Obesity is a crucial factor that increases the risk of initiating and advancing knee osteoarthritis. However, it remains unclear how obesity directly impacts the biomechanical experience of the lower limb joints, potentially triggering or exacerbating joint degeneration. This study investigated the interactive effects of BMI augmentation on lower limb kinematics, kinetics, and muscle activations during walking.

METHODOLOGY

A group of 60 participants underwent a three-dimensional gait analysis. These individuals were categorized into three groups based on their body mass index (BMI): those with a BMI below 25 were classified as having a healthy weight, those with a BMI between 25 and 30 were categorized as overweight, and those with a BMI exceeding 30 were considered obese. This study analyzed the gait of 60 participants categorized by BMI. During walking trials, they recorded ground reaction forces electromyography of leg muscles like the gastrocnemii, hamstrings, and quadriceps. Lower limb joint angles and net moments were also calculated. Statistical mapping identified variations in kinematic, kinetic, and muscle activation patterns across the stance phase between BMI groups.

RESULTS

The results displayed distinct biomechanical patterns in obese individuals. Notably, there was a significant increase in flexion observed in the hip and knee joints (P < 0.001) during the initial stance phase and an increase in hip and knee adduction angles and moments throughout the entire stance phase (P < 0.001). Additionally, muscle activations underwent significant changes (P < 0.01), with a positive correlation noted with the BMI factor. This correlation was most pronounced during the early stance phase for the quadriceps and hamstring muscles and the late stance phase for the gastrocnemius.

CONCLUSION

These findings represent a comprehensive picture that contributes to understanding how excess weight and obesity influence joint biomechanics, highlighting the associated risk of joint osteoarthritis.

The Estimation of Knee Medial Force with Substitution Parameters during Walking and Turning

PURPOSE

Knee adduction, flexion moment, and adduction angle are often used as surrogate parameters of knee medial force. To verify whether these parameters are suitable as surrogates under different walking states, we investigated the correlation between knee medial loading with the surrogates during walking and turning.

METHODS

Sixteen healthy subjects were recruited to complete straight walk (SW), step turn (ST), and crossover turn (CT). Knee joint moments were obtained using inverse dynamics, and knee medial force was computed using a previously validated musculoskeletal model, Freebody. Linear regression was used to predict the peak of knee medial force with the peaks of the surrogate parameters and walking speed.

RESULTS

There was no significant difference in walking speed among these three tasks. The peak knee adduction moment (pKAM) was a significant predictor of the peak knee medial force (pKMF) for SW, ST, and CT (p < 0.001), while the peak knee flexion moment (pKFM) was only a significant predictor of the pKMF for SW (p = 0.034). The statistical analysis showed that the pKMF increased, while the pKFM and the peak knee adduction angle (pKAA) decreased significantly during CT compared to those of SW and ST (p < 0.001). The correlation analysis indicated that the knee parameters during SW and ST were quite similar.

CONCLUSIONS

This study investigated the relationship between knee medial force and some surrogate parameters during walking and turning. KAM was still the best surrogate parameter for SW, ST, and CT. It is necessary to consider the type of movement when comparing the surrogate predictors of knee medial force, as the prediction equations differ significantly among movement types.

The effect of body weight on the knee joint biomechanics based on subject-specific finite element-musculoskeletal approach

Knee osteoarthritis (OA) and obesity are major public health concerns that are closely intertwined. This intimate relationship was documented by considering obesity as the most significant preventable risk factor associated with knee OA. To date, however, the effects of obesity on the knee joint's passive-active structure and cartilage loading have been inconclusive. Hence, this study investigates the intricate relationship between obesity and knee OA, centering on the biomechanical changes in knee joint active and passive reactions during the stance phase of gait. Using a subject-specific musculoskeletal and finite element approach, muscle forces, ligament stresses, and articular cartilage contact stresses were analyzed among 60 individuals with different body mass indices (BMI) classified under healthy weight, overweight, and obese categories. Our predicted results showed that obesity significantly influenced knee joint mechanical reaction, increasing muscle activations, ligament loading, and articular cartilage contact stresses, particularly during key instances of the gait cycle-first and second peak loading instances. The study underscores the critical role of excessive body weight in exacerbating knee joint stress distribution and cartilage damage. Hence, the insights gained provide a valuable biomechanical perspective on the interaction between body weight and knee joint health, offering a clinical utility in assessing the risks associated with obesity and knee OA.

Knee joint biomechanics and cartilage damage prediction during landing: A hybrid MD-FE-musculoskeletal modeling

Understanding the mechanics behind knee joint injuries and providing appropriate treatment is crucial for improving physical function, quality of life, and employability. In this study, we used a hybrid molecular dynamics-finite element-musculoskeletal model to determine the level of loads the knee can withstand when landing from different heights (20, 40, 60 cm), including the height at which cartilage damage occurs. The model was driven by kinematics-kinetics data of asymptomatic subjects at the peak loading instance of drop landing. Our analysis revealed that as landing height increased, the forces on the knee joint also increased, particularly in the vastus muscles and medial gastrocnemius. The patellar tendon experienced more stress than other ligaments, and the medial plateau supported most of the tibial cartilage contact forces and stresses. The load was mostly transmitted through cartilage-cartilage interaction and increased with landing height. The critical height of 126 cm, at which cartilage damage was initiated, was determined by extrapolating the collected data using an iterative approach. Damage initiation and propagation were mainly located in the superficial layers of the tibiofemoral and patellofemoral cartilage. Finally, this study provides valuable insights into the mechanisms of landing-associated cartilage damage and could help limit joint injuries and improve training programs.

High tibiofemoral contact and muscle forces during gait are associated with radiographic knee OA progression over 3 years

BACKGROUND

The objective of this study was to investigate differences in tibiofemoral joint contact forces between individuals with moderate medial OA who exhibit radiographic knee OA progression within 3 years versus those who do not, and to understand the relationship between model-predicted contact forces and net external moments for this population.

METHODS

27 individuals with moderate medial compartment knee OA underwent baseline instrumented gait analysis. OA progressors were defined as those who experienced at least a one grade increase in medial joint space narrowing at three years. An electromyography-driven musculoskeletal model was used to estimate muscle and tibiofemoral contact forces at baseline, which were compared between progressors and non-progressors using t-tests.

RESULTS

Seven individuals experienced radiographic OA progression by 3 years. Progressors walked with significantly higher peaks of medial and total tibiofemoral contact forces, and higher impulse of medial contact forces. Significant and high correlations were found between: first peaks of medial and total contact forces with first peak of the knee adduction moment (R² = 0.74; R² = 0.59); second peaks of medial and total knee contact forces with second peaks of knee adduction and flexion moments (R² = 0.71; R² = 0.68); medial knee contact force impulse with knee adduction moment impulse (R² = 0.76).

CONCLUSIONS

Higher tibiofemoral joint contact forces during walking were associated with three-year radiographic knee OA progression based on medial joint space narrowing. These results support the need for strategies that reduce compressive knee contact forces through the reduction of adduction and flexion moments during walking.

Using musculoskeletal modelling to estimate knee joint loading pre and post high tibial osteotomy

BACKGROUND

Both medial knee osteoarthritis and associated varus alignment have been proposed to alter knee joint loading and consequently overloading the medial compartment. Individuals with knee osteoarthritis and varus deformity are candidates for coronal plane corrective surgery, high tibial osteotomy. This study evaluated knee loading and contact location for a control group, a pre-surgery cohort and the same cohort 12 months post-surgery using a musculoskeletal modelling approach.

METHODS

Joint kinematics during gait were measured in 30 knee osteoarthritis patients, before and after high tibial osteotomy, and 28 healthy adults. Using a musculoskeletal model that incorporated patient-specific mechanical tibial femoral angle, the resulting muscle, ligament, and contact forces were calculated and the medial - lateral condyle load distribution was analysed.

FINDINGS

Surgery changed medial compartment contact force throughout stance relative to pre-surgery. This reduction in medial compartment contact force pre- vs post-HTO is observed despite a significant increase in post-surgery walking speed compared to pre-HTO, where increased speed is typically associated with increased joint loading.

INTERPRETATION

This study has estimated the effects of high tibial osteotomy on knee loading using a generic model that incorporates a detailed knee model to better understand tibiofemoral contact loading. The findings support the aim of surgery to unload the medial knee compartment and lateralise joint contact forces.

Relationships between tibial articular cartilage, in vivo external joint moments and static alignment in end-stage knee osteoarthritis: A micro-CT study

Biomechanical factors (e.g., joint loading) have a significant role in the progression of osteoarthritis (OA). However, some relationships between in vivo joint loading indices and tibial cartilage thickness are conflicting. This study investigated relationships between pre-operative in vivo external knee joint moments, joint alignment and regional tibial cartilage thickness using micro-CT in subjects with end-stage knee OA. Tibial plateaus from 25 patients that underwent knee replacement for OA were micro-CT scanned (17 µm/voxel). Prior to surgery, subjects underwent gait analysis to calculate external knee moments. The mechanical axis deviation (MAD) was obtained from pre-operative radiographs. Cartilage thickness (Cart.Th) was analyzed from micro-CT images, in anteromedial, anterolateral, posteromedial and posterolateral subregions of interest. Medial-to-lateral Cart.Th ratios were also explored. Relationships between Cart.Th and joint loading indices were examined using Pearson's correlations. Significant correlations were found between Cart.Th and joint loading indices, positive anteromedially with the first peak knee adduction moment (r = 0.55, p < 0.01) and external rotation moment (ERM; r = 0.52, p < 0.01), and negative with MAD (r = -0.76, p < 0.001). In the lateral regions, these correlations had opposite signs. The medial-to-lateral Cart.Th ratio correlated strongly with ERM (r = 0.63, p = 0.001) and MAD (r = -0.75, p < 0.001). Joint loading indices correlated with regional cartilage thickness values and their medial-to-lateral ratios in end-stage knee OA subjects, with higher regional loads corresponding to thinner cartilage. These relationships have the opposite sign compared to the subchondral bone microarchitecture found in our previous study on the same specimens, which may suggest a complementary bone-cartilage interplay in response to loading.

Towards Tailored Rehabilitation by Implementation of a Novel Musculoskeletal Finite Element Analysis Pipeline

Tissue-level mechanics (e.g., stress and strain) are important factors governing tissue remodeling and development of knee osteoarthritis (KOA), and hence, the success of physical rehabilitation. To date, no clinically feasible analysis toolbox has been introduced and used to inform clinical decision making with subject-specific in-depth joint mechanics of different activities. Herein, we utilized a rapid state-of-the-art electromyography-assisted musculoskeletal finite element analysis toolbox with fibril-reinforced poro(visco)elastic cartilages and menisci to investigate knee mechanics in different activities. Tissue mechanical responses, believed to govern collagen damage, cell death, and fixed charge density loss of proteoglycans, were characterized within 15 patients with KOA while various daily activities and rehabilitation exercises were performed. Results showed more inter-participant variation in joint mechanics during rehabilitation exercises compared to daily activities. Accordingly, the devised workflow may be used for designing subject-specific rehabilitation protocols. Further, results showed the potential to tailor rehabilitation exercises, or assess capacity for daily activity modifications, to optimally load knee tissue, especially when mechanically-induced cartilage degeneration and adaptation are of interest.

An EMG-assisted muscle-force driven finite element analysis pipeline to investigate joint- and tissue-level mechanical responses in functional activities: towards a rapid assessment toolbox

Joint tissue mechanics (e.g., stress and strain) are believed to have a major involvement in the onset and progression of musculoskeletal disorders, e.g., knee osteoarthritis (KOA). Accordingly, considerable efforts have been made to develop musculoskeletal finite element (MS-FE) models to estimate highly detailed tissue mechanics that predict cartilage degeneration. However, creating such models is time-consuming and requires advanced expertise. This limits these complex, yet promising MS-FE models to research applications with few participants and makes the models impractical for clinical assessments. Also, these previously developed MS-FE models have not been used to assess activities other than gait. This study introduces and verifies a semi-automated rapid state-of-the-art MS-FE modeling and simulation toolbox incorporating an electromyography- (EMG) assisted MS model and a muscle-force driven FE model of the knee with fibril-reinforced poro(visco)elastic cartilages and menisci. To showcase the usability of the pipeline, we estimated joint- and tissue-level knee mechanics in 15 KOA individuals performing different daily activities. The pipeline was verified by comparing the estimated muscle activations and joint mechanics to existing experimental data. To determine the importance of EMG-assisted MS approach, results were compared to those from the same FE models but driven by static-optimization-based MS models. The EMG-assisted MS-FE pipeline bore a closer resemblance to experiments compared to the static-optimization-based MS-FE pipeline. Importantly, the developed pipeline showed great potential as a rapid MS-FE analysis toolbox to investigate multiscale knee mechanics during different activities of individuals with KOA. The template FE model of the study is freely available here.

Biomechanical Characteristics of the Knee Joint during Gait in Obese versus Normal Subjects

Knee osteoarthritis (OA) is a growing source of pain and disability. Obesity is the most important avoidable risk factor underlying knee OA. The processes by which obesity impacts osteoarthritis are of tremendous interest to osteoarthritis researchers and physicians, where the joint mechanical load is one of the pathways generally thought to cause or intensify the disease process. In the current work, we developed a hybrid framework that simultaneously incorporates a detailed finite element model of the knee joint within a musculoskeletal model to compute lower extremity muscle forces and knee joint stresses in normal-weight (N) and obese (OB) subjects during the stance phase gait. This model accounts for the synergy between the active musculature and passive structures. In comparing OB subjects and normal ones, forces significantly increased in all muscle groups at most instances of stance. Mainly, much higher activation was computed with lateral hamstrings and medial gastrocnemius. Cartilage contact average pressure was mostly supported by the medial plateau and increased by 22%, with a larger portion of the load transmitted via menisci. This medial compartment experienced larger relative movement and cartilage stresses in the normal subjects and continued to do so with a higher level in the obese subjects. Finally, the developed bioengineering frame and the examined parameters during this investigation might be useful clinically in evaluating the initiation and propagation of knee OA.

Recent Advances in the Neural Control of Movements: Lessons for Functional Recovery

We review the current views on the control and coordination of movements following the traditions set by Nikolai Bernstein. In particular, we focus on the theory of neural control of effectors - from motor units to individual muscles, to joints, limbs, and to the whole body - with spatial referent coordinates organized into a hierarchy with multiple few-to-many mappings. Further, we discuss synergies ensuring stability of natural human movements within the uncontrolled manifold hypothesis. Synergies are organized within the neural control hierarchy based on the principle of motor abundance. Movement disorders are discussed as consequences of an inability to use the whole range of changes in referent coordinates (as in spasticity) and an inability to ensure controlled stability of salient variables as reflected in indices of multi-element synergies and their adjustments in preparation to actions (as in brain disorders, including Parkinson's disease, multiple-system atrophy, and stroke). At the end of the review, we discuss possible implications of this theoretical approach to peripheral disorders and their rehabilitations using, as an example, osteoarthritis. In particular, "joint stiffening" is viewed as a maladaptive strategy, which can compromise stability of salient variables during walking.

Tibial cartilage, subchondral bone plate and trabecular bone microarchitecture in varus- and valgus-osteoarthritis versus controls

This preliminary study quantified tibia cartilage thickness (Cart.Th), subchondral bone plate thickness (SBPl.Th) and subchondral trabecular bone (STB) microarchitecture in subjects with varus- or valgus- malaligned knees diagnosed with end-stage knee osteoarthritis (OA) and compared them to controls (non-OA). Tibial plateaus from 25 subjects with knee-OA (undergoing knee arthroplasty) and 15 cadavers (controls) were micro-CT scanned (17 µm/voxel). Joint alignment was classified radiographically for OA subjects (varus-aligned n = 18, valgus-aligned n = 7). Cart.Th, SBPl.Th, STB bone volume fraction (BV/TV) and their medial-to-lateral ratios were analyzed in anteromedial, anterolateral, posteromedial and posterolateral subregions. Varus-OA and valgus-OA were compared to controls. Compared to controls (1.19-1.54 mm), Cart.Th in varus-OA was significantly lower anteromedially (0.58 mm, -59%) and higher laterally (2.19-2.47 mm, +60-63%); in valgus-OA, Cart.Th was significantly higher posteromedially (1.86 mm, +56%). Control medial-to-lateral Cart.Th ratios were around unity (0.8-1.1), in varus-OA significantly below (0.2-0.6) and in valgus-OA slightly above (1.0-1.3) controls. SBPl.Th and BV/TV were significantly higher medially in varus-OA (0.58-0.72 mm and 37-44%, respectively) and laterally in valgus-OA (0.60-0.61 mm and 32-37%), compared to controls (0.26-0.47 mm and 18-37%). In varus-OA, the medial-to-lateral SBPl.Th and BV/TV ratios were above unity (1.4-2.4) and controls (0.8-2.1); in valgus-OA they were closer to unity (0.8-1.1) and below controls. Varus- and valgus-OA tibia differ significantly from controls in Cart.Th, SBPl.Th and STB microarchitecture depending on joint alignment, suggesting structural changes in OA may reflect differences in medial-to-lateral load distribution upon the tibial plateau. Here we identified an inverse relationship between cartilage thickness and underlying subchondral bone, suggesting a whole-joint response in OA to daily stimuli.

Effects of an Exercise Therapy Targeting Knee Kinetics on Pain, Function, and Gait Kinetics in Patients With Knee Osteoarthritis: A Randomized Clinical Trial

In this study, the effects of an exercise therapy comprising yoga exercises and medial-thrust gait (YogaMT) on lower-extremity kinetics, pain, and function in patients with medial knee osteoarthritis were investigated. Fifty-nine patients were randomly allocated to three treatment groups: (a) the YogaMT group practiced yoga exercises and medial thrust gait, (b) the knee-strengthening group performed quadriceps- and hamstring-strengthening exercises, and (c) the treadmill walking group practiced normal treadmill walking in 12 supervised sessions. The adduction and flexion moments of the hip, knee, and ankle; pain intensity; and 2-min walking test were assessed before and after treatment and at 1-month follow-up. The YogaMT group experienced a significant reduction in knee adduction moment. All groups showed significant improvement in pain and function. The YogaMT may reduce medial knee load in patients with knee osteoarthritis in the short term. A larger clinical trial is required to investigate the long-term outcomes of this intervention.

Effect of medial foot loading self-practice on lower limb kinematics in young individuals with asymptomatic varus knee alignment

BACKGROUND

Varus alignment of the knee is a risk factor for developing knee osteoarthritis. Recently, voluntary shifting the plantar pressure distribution medially (medial foot loading) during gait has been found to reduce knee adduction angle during stance, which may lower the joint load. However, it is not yet known whether such effect would persist after long-term self-practice. This study aimed to determine whether medial foot loading can be an effective self-care protocol for reducing the knee adduction angle.

METHODS

Eight subjects with asymptomatic varus knee alignment were trained on medial foot loading once in a laboratory, then carried out as self-practice for 8 weeks outside the laboratory. Spatiotemporal gait parameters and lower limb joint kinematics data were collected during natural walking prior to the training (baseline walking), during the practice session immediately after the initial training (trained walking), and during natural walking after the self-practice period (post-practice walking).

RESULTS

Participants walked significantly faster after the self-practice period with longer step length compared with the baseline. The knee adduction angle at initial contact, maximum angle during stance, and mean angle during a gait cycle were significantly decreased during both the trained and post-practice walking compared with baseline. The 8-week self-practice caused larger decrements in the three angles than the single training, but no significant differences were found between the two conditions.

CONCLUSIONS

Self-practice of medial foot loading walking could be an effective gait strategy to reduce the knee adduction angle. The effect could be sustained for individuals with asymptomatic varus knee alignment.

The effects of knee pain on knee contact force and external knee adduction moment in patients with knee osteoarthritis

Knee osteoarthritis (OA) is a major cause of knee pain, leading to physical dysfunction. External knee adduction moment (KAM), a surrogate measure of knee contact force (KCF) in the medial compartment, is related to knee pain, but the association between KCF and pain severity remains unclear. This study aimed to reveal the differences in KCF due to pain severity. Twenty-eight patients with knee OA were evaluated knee symptoms including pain severity via the Knee Society Score. Based on the median symptom score, 17 points in this study, subjects were classified as having Mild symptomatic OA (n = 15) and Severe symptomatic OA (n = 13). Subjects walked three times at a comfortable speed along a six-meter walkway, and we calculated KAM during the stance phase. KCF magnitude and distribution were also computed using the subject-specific musculoskeletal model, considering physical characteristics such as the femorotibial angle measured by X-ray. No differences in physical characteristics such as femorotibial angle and gait speed were found by symptom severity, whereas KAM and medial KCF at minimum and second peak in Severe symptomatic OA patients were significantly greater than those in Mild symptomatic OA. A significant medial shift of KCF in Severe symptomatic OA was also seen at first peak and minimum. Severe symptomatic OA had a greater medial KCF and medial shift of KCF. Detailed evaluations of KCF magnitude and distribution in addition to KAM would provide crucial information on knee contact force in relation to symptom severity.

Knee loading in OA subjects is correlated to flexion and adduction moments and to contact point locations

This study evaluated the association of contact point locations with the knee medial and lateral contact force (F_med, F_lat) alterations in OA and healthy subjects. A musculoskeletal model of the lower limb with subject-specific tibiofemoral contact point trajectories was used to estimate the F_med and F_lat in ten healthy and twelve OA subjects during treadmill gait. Regression analyses were performed to evaluate the correlation of the contact point locations, knee adduction moment (KAM), knee flexion moment (KFM), frontal plane alignment, and gait speed with the F_med and F_lat. Medial contact point locations in the medial–lateral direction showed a poor correlation with the F_med in OA (R² = 0.13, p = 0.01) and healthy (R² = 0.24, p = 0.001) subjects. Anterior–posterior location of the contact points also showed a poor correlation with the F_med of OA subjects (R² = 0.32, p < 0.001). Across all subjects, KAM and KFM remained the best predictors of the F_med and F_lat, respectively (R² between 0.62 and 0.69). Results suggest different mechanisms of contact force distribution in OA joints. The variations in the location of the contact points participate partially to explains the F_med variations in OA subjects together with the KFM and KAM.

EMG-Assisted Muscle Force Driven Finite Element Model of the Knee Joint with Fibril-Reinforced Poroelastic Cartilages and Menisci

Abnormal mechanical loading is essential in the onset and progression of knee osteoarthritis. Combined musculoskeletal (MS) and finite element (FE) modeling is a typical method to estimate load distribution and tissue responses in the knee joint. However, earlier combined models mostly utilize static-optimization based MS models and muscle force driven FE models typically use elastic materials for soft tissues or analyze specific time points of gait. Therefore, here we develop an electromyography-assisted muscle force driven FE model with fibril-reinforced poro(visco)elastic cartilages and menisci to analyze knee joint loading during the stance phase of gait. Moreover, since ligament pre-strains are one of the important uncertainties in joint modeling, we conducted a sensitivity analysis on the pre-strains of anterior and posterior cruciate ligaments (ACL and PCL) as well as medial and lateral collateral ligaments (MCL and LCL). The model produced kinematics and kinetics consistent with previous experimental data. Joint contact forces and contact areas were highly sensitive to ACL and PCL pre-strains, while those changed less cartilage stresses, fibril strains, and fluid pressures. The presented workflow could be used in a wide range of applications related to the aetiology of cartilage degeneration, optimization of rehabilitation exercises, and simulation of knee surgeries.

EMG-Assisted Muscle Force Driven Finite Element Model of the Knee Joint with Fibril-Reinforced Poroelastic Cartilages and Menisci

Abnormal mechanical loading is essential in the onset and progression of knee osteoarthritis. Combined musculoskeletal (MS) and finite element (FE) modeling is a typical method to estimate load distribution and tissue responses in the knee joint. However, earlier combined models mostly utilize static-optimization based MS models and muscle force driven FE models typically use elastic materials for soft tissues or analyze specific time points of gait. Therefore, here we develop an electromyography-assisted muscle force driven FE model with fibril-reinforced poro(visco)elastic cartilages and menisci to analyze knee joint loading during the stance phase of gait. Moreover, since ligament pre-strains are one of the important uncertainties in joint modeling, we conducted a sensitivity analysis on the pre-strains of anterior and posterior cruciate ligaments (ACL and PCL) as well as medial and lateral collateral ligaments (MCL and LCL). The model produced kinematics and kinetics consistent with previous experimental data. Joint contact forces and contact areas were highly sensitive to ACL and PCL pre-strains, while those changed less cartilage stresses, fibril strains, and fluid pressures. The presented workflow could be used in a wide range of applications related to the aetiology of cartilage degeneration, optimization of rehabilitation exercises, and simulation of knee surgeries.

The influence of knee joint geometry and alignment on the tibiofemoral load distribution: A computational study

BACKGROUND

Deviations in knee joint geometry and alignment were previously related to an increased risk for knee OA. These were hypothesized to influence the load distribution over the articular cartilage. Therefore, this study evaluated the effect of altered knee joint geometry and alignment in the coronal and transverse plane on the medial-lateral load distribution and ligament strain using a musculoskeletal modeling approach.

METHODS

Joint kinematics during gait were measured in 15 healthy adults. Using different musculoskeletal models with altered geometry of the tibia plateau or knee joint malalignment in the coronal and transverse plane, the resulting muscle, ligament and contact forces were calculated. Next, the distribution of the load over the medial and lateral condyle was analyzed and compared to the reference loading distribution, with neutral geometry and alignment, using repeated-measures ANOVA and individual t-tests, with a Bonferroni-corrected alpha level.

RESULTS

Coronal plane malalignment significantly affected the load distribution. Small changes in coronal tibial slope had less pronounced effects on the load distribution, but increased ligament strains. Transverse plane malalignment only minimally affected the load distribution.

CONCLUSION

Coronal plane knee malalignment affected knee loading, with increased varus alignment resulting in increased medial loading. This confirms a causal relation between coronal malalignment and increased medial compartment loading and suggests a potential role of aberrant coronal plane alignment on OA initiation. Altered coronal tibial slope induced increased ligament strains, potentially contributing to a cascade of knee laxity and subsequently more extreme knee malalignment.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Sensitivity of medial-lateral load sharing to changes in adduction moments or angles in an asymptomatic knee joint model during gait

BACKGROUND

Osteoarthritis (OA) of the knee joint is a common disease accompanied by pain and impaired mobility. Despite some recent concerns on the lack of correlation between the medial load and the knee adduction moment (KAM), KAM is routinely considered as a surrogate measure of medial load and hence a marker where its reduction is the main focus of preventive and treatment interventions.

RESEARCH QUESTION

Determine the relative sensitivity of the tibiofemoral medial-lateral contact load partitioning to changes in the knee adduction angle (KAA) versus KAM.

METHODS

Using a lower extremity hybrid musculoskeletal (MS) model driven by gait kinematics and kinetics, we compute here in asymptomatic subjects the sensitivity of the knee joint biomechanical response (muscle and ligament forces) in general and medial/lateral load partitioning in particular to the relative changes in the reported KAA versus changes in reported KAM (both by one standard deviation).

RESULTS

As KAA increased (at constant KAM), so did the passive moment resistance of the knee joint which as a result and at all stance periods substantially reduced forces in lateral hamstrings while increasing those in medial hamstrings. At 25% and 75% stance as two highly loaded periods of gait, the drop in KAA (from + SD to -SD while at constant KAM) drastically reduced the medial contact force by 44% and 30% and the medial over lateral contact load and area ratios by 92% and 79% as well as 64% and 51%, respectively. In contrast, the equivalent alterations in KAM (by ± SD at constant KAA) had lower and less consistent effects (<7%) showing much smaller sensitivity to changes in KAM alone. Ligament forces altered at various stance periods with inconsistent trends; peak values of 418 N in the anterior cruciate ligament (90% carried by the posterolateral bundle) and 1056 N in the patellar tendon were computed both at 25% stance and minimum KAA.

SIGNIFICANCE

These findings indicate a poor correlation between KAM and tibiofemoral load distribution suggesting instead that KAA and knee alignment should be in focus as the primary marker of knee joint load partitioning and associated prevention and treatment interventions.

Topographical Variation of Human Femoral Articular Cartilage Thickness, T1rho and T2 Relaxation Times Is Related to Local Loading during Walking

OBJECTIVE

Early detection of degenerative changes in the cartilage matrix composition is essential for evaluating early interventions that slow down osteoarthritis (OA) initiation. T1rho and T2 relaxation times were found to be effective for detecting early changes in proteoglycan and collagen content. To use these magnetic resonance imaging (MRI) methods, it is important to document the topographical variation in cartilage thickness, T1rho and T2 relaxation times in a healthy population. As OA is partially mechanically driven, the relation between these MRI-based parameters and localized mechanical loading during walking was investigated.

DESIGN

MR images were acquired in 14 healthy adults and cartilage thickness and T1rho and T2 relaxation times were determined. Experimental gait data was collected and processed using musculoskeletal modeling to identify weight-bearing zones and estimate the contact force impulse during gait. Variation of the cartilage properties (i.e., thickness, T1rho, and T2) over the femoral cartilage was analyzed and compared between the weight-bearing and non-weight-bearing zone of the medial and lateral condyle as well as the trochlea.

RESULTS

Medial condyle cartilage thickness was correlated to the contact force impulse ( r = 0.78). Lower T1rho, indicating increased proteoglycan content, was found in the medial weight-bearing zone. T2 was higher in all weight-bearing zones compared with the non-weight-bearing zones, indicating lower relative collagen content.

CONCLUSIONS

The current results suggest that medial condyle cartilage is adapted as a long-term protective response to localized loading during a frequently performed task and that the weight-bearing zone of the medial condyle has superior weight bearing capacities compared with the non-weight-bearing zones.

The effect of surface inclination and limb on knee loading measures in transtibial prosthesis users

Background

Osteoarthritis (OA) is a degenerative disease caused by the wearing of joint cartilage and bone. Literature has established that a prosthesis user’s intact limb is at greater risk of developing OA. This study analyzed the effect of commonly encountered surface inclinations on knee joint loading measures in able-bodied and transtibial prosthesis users.

Methods

12 transtibial prosthesis users and 12 able-bodied participants walked across level ground, up slope, down slope, and cross slope (further divided into top and bottom slope depending on the location of the limb being analyzed). First and second peak external knee adduction moment (KAM), external knee adduction moment rate, and external knee adduction moment impulse were extracted from the stance phase of gait. Mixed ANOVA statistics with Bonferonni post hoc analyses were performed.

Results

Significant limb differences were only found for KAM rate and first peak KAM. When compared to all other surfaces up slope had the significantly lowest KAM rate and was not significantly lower for all other tested variables. Down slope had significantly greater KAM rate than all surfaces except bottom slope. KAM second peak and KAM impulse analysis resulted in no significant differences.

Conclusions

Individuals at risk for developing, or currently dealing with, knee OA could avoid walking for extended periods on down slope. Walking up moderate slopes may be considered as a complementary activity to level walking for rehabilitation and delaying OA progression.

The lack of significant limb differences suggests that second peak KAM and KAM impulse may not be appropriate load-related indicators of OA initiation among prosthesis users without OA. KAM rate was the most sensitive joint loading variable and therefore should be investigated further as an appropriate variable for identifying OA risk in individuals with transtibial amputations.

The effect of articular geometry features identified using statistical shape modelling on knee biomechanics

Articular geometry in the knee varies widely among people which has implications for risk of injury and pathology. The goals of this work were to develop a framework to systematically vary geometry in a multibody knee model and to use this framework to investigate the effect of morphological features on dynamic knee kinematics and contact mechanics. A statistical shape model of the tibiofemoral and patellofemoral joints was created from magnetic resonance images of 14 asymptomatic knees. The shape model was then used to generate 37 unique multibody knee models based on -3 to +3 standard deviations of the scores for the first six principal components identified. Each multibody model was then incorporated into a lower extremity musculoskeletal model and the Concurrent Optimization of Muscle Activations and Kinematics (COMAK) routine was used to simulate knee mechanics for overground walking. Changes in articular geometry affected knee function, resulting in differences up to 17° in orientation, 8 mm in translation, 0.7 BW in contact force, and 2.0 MPa in mean cartilage contact pressure. Understanding the relationship between shape and function in a joint could provide insight into the mechanisms behind injury and pathology and the variability in response to treatment.

A lower extremity model for muscle-driven simulation of activity using explicit finite element modeling

A key strength of computational modeling is that it can provide estimates of muscle, ligament, and joint loads, stresses, and strains through non-invasive means. However, simulations that can predict the forces in the muscles during activity while maintaining sufficient complexity to realistically represent the muscles and joint structures can be computationally challenging. For this reason, the current state of the art is to apply separate rigid-body dynamic and finite-element (FE) analyses in series. However, the use of two or more disconnected models often fails to capture key interactions between the joint-level and whole-body scales. Single framework MSFE models have the potential to overcome the limitations associated with disconnected models in series. The objectives of the current study were to create a multi-scale FE model of the human lower extremity that combines optimization, dynamic muscle modeling, and structural FE analysis in a single framework and to apply this framework to evaluate the mechanics of healthy knee specimens during two activities. Two subject-specific FE models (Model 1, Model 2) of the lower extremity were developed in ABAQUS/Explicit including detailed representations of the muscles. Muscle forces, knee joint loading, and articular contact were calculated for two activities using an inverse dynamics approach and static optimization. Quadriceps muscle forces peaked at the onset of chair rise (2174 N, 1962 N) and in early stance phase (510 N, 525 N), while gait saw peak forces in the hamstrings (851 N, 868 N) in midstance. Joint forces were similar in magnitude to available telemetric patient data. This study demonstrates the feasibility of detailed quasi-static, muscle-driven simulations in an FE framework.

Subjects with medial and lateral tibiofemoral articular cartilage defects do not alter compartmental loading during walking

BACKGROUND

Healthy cartilage is essential for optimal joint function. Although, articular cartilage defects are highly prevalent in the active population and hamper joint function, the effect of articular cartilage defects on knee loading is not yet documented. Therefore, the present study compared knee contact forces and pressures between patients with tibiofemoral cartilage defects and healthy controls. Potentially this provides additional insights in movement adaptations and the role of altered loading in the progression from defect towards OA.

METHODS

Experimental gait data collected in 15 patients with isolated cartilage defects (8 medial involvement, 7 lateral-involvement) and 19 healthy asymptomatic controls was processed using a musculoskeletal model to calculate contact forces and pressures. Differences between two patient groups and controls were evaluated using Kruskal-Wallis tests and individually compared using Mann-Whitney-U tests (alpha <0.05).

FINDINGS

The patients with lateral involvement walked significantly slower compared to the healthy controls. No movement adaptations to decrease the loading on the injured condyle were observed. Additionally, the location of loading was not significantly affected.

INTERPRETATION

The current results suggest that isolated cartilage defects do not induce significant changes in the knee joint loading distribution. Consequently, the involved condyle will capture a physiological loading magnitude that should however be distributed over the cartilage surrounding the defect. This may cause local degenerative changes in the cartilage and in combination with inflammatory responses, might play a key role in the progression from articular cartilage defect to a more severe OA phenotype.

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

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

Frontal Plane Tibiofemoral Alignment is Strongly Related to Compartmental Knee Joint Contact Forces and Muscle Control Strategies During Stair Ascent

Static frontal plane tibiofemoral alignment is an important factor in dynamic knee alignment and knee adduction moments (KAMs). However, little is known about the relationship between alignment and compartment contact forces or muscle control strategies. The purpose of this study was to estimate medial (MCF) and lateral (LCF) compartment knee joint contact forces and muscle forces during stair ascent using a musculoskeletal model implementing subject-specific knee alignments. Kinematic and kinetic data from 20 healthy individuals with radiographically confirmed varus or valgus knee alignments were simulated using alignment specific models to predict MCFs and LCFs. Muscle forces were determined using static optimization. Independent samples t-tests compared contact and muscle forces between groups during weight acceptance and during pushoff. The varus group exhibited increased weight acceptance peak MCFs, while the valgus group exhibited increased pushoff peak LCFs. The varus group utilized increased vasti muscle forces during weight acceptance and adductor forces during pushoff. The valgus group utilized increased abductor forces during pushoff. The alignment-dependent contact forces provide evidence of the significance of frontal plane knee alignment in healthy individuals, which may be important in considering future knee joint health. The differing muscle control strategies between alignments detail-specific neuromuscular responses to control frontal plane knee loads.

Relationships between in vivo dynamic knee joint loading, static alignment and tibial subchondral bone microarchitecture in end-stage knee osteoarthritis

OBJECTIVE

To study, in end-stage knee osteoarthritis (OA) patients, relationships between indices of in vivo dynamic knee joint loads obtained pre-operatively using gait analysis, static knee alignment, and the subchondral trabecular bone (STB) microarchitecture of their excised tibial plateau quantified with 3D micro-CT.

DESIGN

Twenty-five knee OA patients scheduled for total knee arthroplasty underwent pre-operative gait analysis. Mechanical axis deviation (MAD) was determined radiographically. Following surgery, excised tibial plateaus were micro-CT-scanned and STB microarchitecture analysed in four subregions (anteromedial, posteromedial, anterolateral, posterolateral). Regional differences in STB microarchitecture and relationships between joint loading and microarchitecture were examined.

RESULTS

STB microarchitecture differed among subregions (P < 0.001), anteromedially exhibiting highest bone volume fraction (BV/TV) and lowest structure model index (SMI). Anteromedial BV/TV and SMI correlated strongest with the peak external rotation moment (ERM; r = -0.74, r = 0.67, P < 0.01), despite ERM being the lowest (by factor of 10) of the moments considered, with majority of ERM measures below accuracy thresholds; medial-to-lateral BV/TV ratios correlated with ERM, MAD, knee adduction moment (KAM) and internal rotation moment (|r|-range: 0.54-0.74). When controlling for walking speed, KAM and MAD, the ERM explained additional 11-30% of the variations in anteromedial BV/TV and medial-to-lateral BV/TV ratio (R² = 0.59, R² = 0.69, P < 0.01).

CONCLUSIONS

This preliminary study suggests significant associations between tibial plateau STB microarchitecture and knee joint loading indices in end-stage knee OA patients. Particularly, anteromedial BV/TV correlates strongest with ERM, whereas medial-to-lateral BV/TV ratio correlates strongest with indicators of medial-to-lateral joint loading (MAD, KAM) and rotational moments. However, associations with ERM should be interpreted with caution.

Knee Joint Loading in Healthy Adults During Functional Exercises: Implications for Rehabilitation Guidelines

Study Design Controlled laboratory study. Background The inclusion of specific exercises in rehabilitation after knee injury is currently expert based, as a thorough description of the knee contact forces during different exercises is lacking. Objective To quantify knee loading during frequently used activities such as squats, lunges, single-leg hops, walking stairs, standing up, and gait, and to grade knee joint loading during these activities. Methods Three-dimensional motion-analysis data of 15 healthy adults were acquired during 9 standardized activities used in rehabilitation. Experimental motion data were processed using musculoskeletal modeling to calculate contact and shear forces on the different knee compartments (tibiofemoral and patellofemoral). Using repeated-measures analyses of variance, contact and shear forces were compared between compartments and exercises, whereas muscle and average maximum femoral forces were compared only between exercises. Results With the exception of squats, all therapeutic exercises imposed higher forces to the tibiofemoral joint compared to gait. Likewise, patellofemoral forces were greater during all exercises when compared to gait. Greater compartmental contact forces were accompanied by greater compartmental shear forces. Furthermore, force distribution over the medial and lateral compartments varied between exercises. With increased knee flexion, more force was imposed on the posterior portion of the condyles. Conclusion These results suggest that with careful selection of exercises, forces on an injured zone of the joint can be reduced, as the force distribution differs strongly between exercises. Based on the results, a graded exercise program for progressive knee joint loading during rehabilitation can be conceptualized. J Orthop Sports Phys Ther 2018;48(3):162-173. Epub 6 Jan 2018. doi:10.2519/jospt.2018.7459.

Joint loading and proximal tibia subchondral trabecular bone microarchitecture differ with walking gait patterns in end-stage knee osteoarthritis

OBJECTIVES

To (1) stratify patient subgroups according to their distinct walking gait patterns in end-stage knee osteoarthritis (OA); (2) compare measures of joint loading and proximal tibia subchondral trabecular bone (STB) microarchitecture among these gait subgroups.

DESIGN

Twenty-five knee OA patients undergoing total knee arthroplasty (TKA) had pre-operative gait analysis. Following surgery, excised tibial plateaus were micro-CT-scanned and STB microarchitecture analysed in four tibial condylar regions of interest. Peak knee moments were input to k-means cluster analysis, to identify subgroups with homogeneous gait patterns. Joint loading and STB microarchitecture parameters were compared among gait subgroups (Kruskal-Wallis, Bonferroni-corrected Mann-Whitney U tests).

RESULTS

Three gait subgroups were revealed: biphasics (n = 7), flexors (n = 9), counter-rotators (n = 9). Peak knee adduction moment (KAM) and KAM impulse were significantly higher (P < 0.05) in biphasics than in flexors and counter-rotators (KAM = -0.65, -0.40 and -0.21 Nm/kg, respectively), suggesting a higher medial-to-lateral tibiofemoral load ratio in biphasics. Interestingly, STB medial-to-lateral bone volume fraction (BV/TV) ratio was also significantly higher (more than double) in biphasics and flexors than in counter-rotators (2.24, 2.00 and 1.00, respectively), whereas in biphasics it was only 10% higher than in flexors and not significantly so.

CONCLUSIONS

Within the confines of the limited sample size, data suggests that different mechanisms between the biphasic and flexor gait subroups may generate comparable loads upon the tibial plateau and corresponding bony responses, despite significantly lower KAM indices in flexors. Hence, in flexor gait OA patients, conservative treatments designed to reduce KAM, may not be appropriate. Understanding joint loading among walking gait patterns and relationships to bone microarchitecture may aid at identifying/improving management of persons at risk for developing knee OA.

Systematic mapping of the subchondral bone 3D microarchitecture in the human tibial plateau: Variations with joint alignment

Tibial subchondral bone plays an important role in knee osteoarthritis (OA). Microarchitectural characterization of subchondral bone plate (SBP), underlying subchondral trabecular bone (STB) and relationships between these compartments, however, is limited. The aim of this study was to characterize the spatial distribution of SBP thickness, SBP porosity and STB microarchitecture, and relationships among them, in OA tibiae of varying joint alignment. Twenty-five tibial plateaus from end-stage knee-OA patients, with varus (n = 17) or non-varus (n = 8) alignment were micro-CT scanned (17 μm/voxel). SBP and STB microarchitecture was quantified via a systematic mapping in 22 volumes of interest per knee (11 medial, 11 lateral). Significant within-condylar and between-condylar (medial vs. lateral) differences (p < 0.05) were found. In varus, STB bone volume fraction (BV/TV) was consistently high throughout the medial condyle, whereas in non-varus, medially, it was more heterogeneously distributed. Regions of high SBP thickness were co-located with regions of high STB BV/TV underneath. In varus, BV/TV was significantly higher medially than laterally, however, not so in non-varus. Moreover, region-specific significant associations between the SBP thickness and SBP porosity and the underlying STB microarchitecture were detected, which in general were not captured when considering the values averaged for each condyle. As subchondral bone changes reflect responses to local mechanical and biochemical factors within the joint, our results suggest that joint alignment influences both the medial-to-lateral and the within-condyle distribution of force across the tibia, generating corresponding local bony responses (adaptation) of both the subchondral bone plate and underlying subchondral trabecular bone microarchitecture. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:1927-1941, 2017.

Gait mechanics in those with/without medial compartment knee osteoarthritis 5 years after anterior cruciate ligament reconstruction

The objective of the study was to evaluate differences in gait mechanics 5 years after unilateral anterior cruciate ligament reconstruction surgery, for non-osteoarthritic (n = 24) versus osteoarthritic (n = 9) subjects. For the involved knee, the osteoarthritic group demonstrated significantly lower peak knee flexion angles (non-osteoarthritic = 24.3 ± 4.6°, osteoarthritic = 19.1 ± 2.9°, p = 0.01) and peak knee flexion moments (non-osteoarthritic = 5.3 ± 1.2% Body Weight × Height, osteoarthritic = 4.4 ± 1.2% Body Weight × Height, p = 0.05). Differences in peak knee adduction moment approached significance, with a higher magnitude for the osteoarthritic group (non-osteoarthritic = 2.4 ± 0.8% Body Weight × Height, osteoarthritic = 2.9 ± 0.5% Body Weight × Height, p = 0.09). Peak medial compartment joint load was evaluated using electromyography-informed neuromusculoskeletal modeling. Peak medial compartment joint load in the involved knee for the two groups was not different (non-osteoarthritic = 2.4 ± 0.4 Body Weight, osteoarthritic = 2.3 ± 0.6 Body Weight). The results suggest that subjects with dissimilar peak knee moments can have similar peak medial compartment joint load magnitudes. There was no evidence of inter-limb asymmetry for either group. Given the presence of inter-group differences (non-osteoarthritic vs. osteoarthritic) for the involved knee, but an absence of inter-limb asymmetry in either group, it may be necessary to evaluate how symmetry is achieved, over time, and to differentiate between good versus bad inter-limb symmetry, when evaluating knee gait parameters. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:625-633, 2017.

Musculoskeletal loading in the symptomatic and asymptomatic knees of middle-aged osteoarthritis patients

This study quantified the contributions by muscles, gravity, and inertia to the tibiofemoral compartment forces in the symptomatic (SYM) and asymptomatic (ASYM) limbs of varus mal-aligned medial knee osteoarthritis (OA) patients, and compared the results with healthy controls (CON). Muscle forces and tibiofemoral compartment loads were calculated using gait data from 39 OA patients and 15 controls aged 49 ± 7 years. Patients exhibited lower knee flexion angle, higher hip abduction, and knee adduction angles, lower internal knee flexion torque but higher external knee adduction moment. Muscle forces were highest in CON except hamstrings, which was highest in SYM. ASYM muscle forces were lowest for biceps femoris short head and gastrocnemius but otherwise intermediate between SYM and CON. In all subjects, vasti, hamstrings, gastrocnemius, soleus, gluteus medius, gluteus maximus, and gravity were the largest contributors to medial compartment force (MCF). Inertial contributions were negligible. Highest MCF was found in SYM throughout stance. Small increases in contributions from hamstrings, gluteus maximus, gastrocnemius, and gravity at the first peak; soleus and rectus femoris at the second peak; and soleus, gluteus maximus, gluteus medius, and gravity during mid-stance summed to produce significantly higher total MCF. Compared to CON, the ASYM limb exhibited similar peak MCF but higher mid-stance MCF. In patients, diminished non-knee-spanning muscle forces did not produce correspondingly diminished MCF contributions due to the influence of mal-alignment. Our findings emphasize consideration of muscle function, lower-limb alignment, and mid-stance loads in developing interventions for OA, and inclusion of the asymptomatic limb in clinical assessments. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:321-330, 2017.

Comparison of different material models of articular cartilage in 3D computational modeling of the knee: Data from the Osteoarthritis Initiative (OAI)

The intricate properties of articular cartilage and the complexity of the loading environment are some of the key challenges in developing models for biomechanical analysis of the knee joint. Fibril-reinforced poroelastic (FRPE) material models have been reported to accurately capture characteristic responses of cartilage during dynamic and static loadings. However, high computational and time costs associated with such advanced models limit applicability of FRPE models when multiple subjects need to be analyzed. If choosing simpler material models, it is important to show that they can still produce truthful predictions. Therefore, the aim of this study was to compare depth-dependent maximum principal stresses and strains within articular cartilage in the 3D knee joint between FRPE material models and simpler isotropic elastic (IE), isotropic poroelastic (IPE) and transversely isotropic poroelastic (TIPE) material models during simulated gait cycle. When cartilage-cartilage contact pressures were matched between the models (15% allowed difference), maximum principal stresses in the IE, IPE and TIPE models were substantially lower than those in the FRPE model (by more than 50%, TIPE model being closest to the FRPE model), and stresses occurred only in compression in the IE model. Additional simulations were performed to find material parameters for the TIPE model (due to its anisotropic nature) that would yield maximum principal stresses similar to the FRPE model. The modified homogeneous TIPE model was in a better agreement with the homogeneous FRPE model, and the average and maximum differences in maximum principal stresses throughout the depth of cartilage were 7% and 9%, respectively, in the lateral compartment and 9% and 11% in the medial compartment. This study revealed that it is possible to match simultaneously maximum principal stresses and strains of cartilage between non-fibril-reinforced and fibril-reinforced knee joint models during gait. Depending on the research question (such as analysis of fibril strain necessitates the use of fibril-reinforced material models) or clinical demand (fast simulations with simpler material models), the choice of the material model should be done carefully.

Wedged Insoles and Gait in Patients with Knee Osteoarthritis: A Biomechanical Review

The study of gait biomechanics in individuals with knee osteoarthritis has become widespread, especially in regards to the knee adduction moment-a variable commonly believed to be associated with knee osteoarthritis progression. Unfortunately, this variable is often studied clinically without considering how it is derived, or what it means in a mechanical context. The use of footwear for knee osteoarthritis management has received much attention as well. However, in many cases, footwear is studied without regard for the mechanical effects they actually induce on the patient. Therefore, this review aims to summarize the current state of knowledge in regards to knee osteoarthritis gait and footwear biomechanics, by taking a step back to review the foundations of these two research areas. First, an overview of the calculation of the knee adduction moment is provided, along with mechanical considerations. Then, this is used to discuss current evidence for wedged insoles and highlight knowledge gaps. The intent was to place this mechanical information in a clinically-oriented framework for approachability by scientists, engineers and clinicians alike. Based on this discussion, areas for future investigation are proposed.

Importance of Patella, Quadriceps Forces, and Depthwise Cartilage Structure on Knee Joint Motion and Cartilage Response During Gait

In finite-element (FE) models of the knee joint, patella is often omitted. We investigated the importance of patella and quadriceps forces on the knee joint motion by creating an FE model of the subject's knee. In addition, depthwise strains and stresses in patellar cartilage with different tissue properties were determined. An FE model was created from subject's magnetic resonance images. Knee rotations, moments, and translational forces during gait were recorded in a motion laboratory and used as an input for the model. Three material models were implemented into the patellar cartilage: (1) homogeneous model, (2) inhomogeneous (arcadelike fibrils), and (3) random fibrils at the superficial zone, mimicking early stages of osteoarthritis (OA). Implementation of patella and quadriceps forces into the model substantially reduced the internal–external femoral rotations (versus without patella). The simulated rotations in the model with the patella matched the measured rotations at its best. In the inhomogeneous model, maximum principal stresses increased substantially in the middle zone of the cartilage. The early OA model showed increased compressive strains in the superficial and middle zones of the cartilage and decreased stresses and fibril strains especially in the middle zone. The results suggest that patella and quadriceps forces should be included in moment- and force-driven FE knee joint models. The results indicate that the middle zone has a major role in resisting shear forces in the patellar cartilage. Also, early degenerative changes in the collagen network substantially affect the cartilage depthwise response in the patella during walking.

Osteoarthritis year in review 2015: mechanics

Motivated by the conceptual framework of multi-scale biomechanics, this narrative review highlights recent major advances with a focus on gait and joint kinematics, then tissue-level mechanics, cell mechanics and mechanotransduction, matrix mechanics, and finally the nanoscale mechanics of matrix macromolecules. A literature review was conducted from January 2014 to April 2015 using PubMed to identify major developments in mechanics related to osteoarthritis (OA). Studies of knee adduction, flexion, rotation, and contact mechanics have extended our understanding of medial compartment loading. In turn, advances in measurement methodologies have shown how injuries to both the meniscus and ligaments, together, can alter joint kinematics. At the tissue scale, novel findings have emerged regarding the mechanics of the meniscus as well as cartilage superficial zone. Moving to the cell level, poroelastic and poro-viscoelastic mechanisms underlying chondrocyte deformation have been reported, along with the response to osmotic stress. Further developments have emerged on the role of calcium signaling in chondrocyte mechanobiology, including exciting findings on the function of mechanically activated cation channels newly found to be expressed in chondrocytes. Finally, AFM-based nano-rheology systems have enabled studies of thin murine tissues and brush layers of matrix molecules over a wide range of loading rates including high rates corresponding to impact injury. With OA acknowledged to be a disease of the joint as an organ, understanding mechanical behavior at each length scale helps to elucidate the connections between cell biology, matrix biochemistry and tissue structure/function that may play a role in the pathomechanics of OA.

Alterations in knee contact forces and centers in stance phase of gait: A detailed lower extremity musculoskeletal model

Evaluation of contact forces-centers of the tibiofemoral joint in gait has crucial biomechanical and pathological consequences. It involves however difficulties and limitations in in vitro cadaver and in vivo imaging studies. The goal is to estimate total contact forces (CF) and location of contact centers (CC) on the medial and lateral plateaus using results computed by a validated finite element model simulating the stance phase of gait for normal as well as osteoarthritis, varus-valgus and posterior tibial slope altered subjects. Using foregoing contact results, six methods commonly used in the literature are also applied to estimate and compare locations of CC at 6 periods of stance phase (0%, 5%, 25%, 50%, 75% and 100%). TF joint contact forces are greater on the lateral plateau very early in stance and on the medial plateau thereafter during 25-100% stance periods. Large excursions in the location of CC (>17mm), especially on the medial plateau in the mediolateral direction, are computed. Various reported models estimate quite different CCs with much greater variations (~15mm) in the mediolateral direction on both plateaus. Compared to our accurately computed CCs taken as the gold standard, the centroid of contact area algorithm yielded least differences (except in the mediolateral direction on the medial plateau at ~5mm) whereas the contact point and weighted center of proximity algorithms resulted overall in greatest differences. Large movements in the location of CC should be considered when attempting to estimate TF compartmental contact forces in gait.

Calculation of external knee adduction moments: a comparison of an inverse dynamics approach and a simplified lever-arm approach

BACKGROUND

The external knee adduction moment (EKAM) is often studied in knee osteoarthritis research. This study compared EKAMs between two methods of calculation: a method that only requires ground reaction force and knee position data (i.e. lever-arm), and an inverse dynamics link-segment method.

METHODS

Sixteen participants walked while wearing a control shoe with and without a six millimeter lateral wedge insole. Peak EKAMs between the lever-arm and inverse dynamics methods were compared for the control condition, and the %change in moment induced by the lateral wedge was compared between methods.

RESULTS

When comparing EKAMs between methods, no correlation was found (r=0.24, p=0.36); peak EKAMs with the lever-arm method (26.0Nm) were significantly lower than EKAMs with the inverse dynamics method (40.2Nm, pb0.001); and Bland-Altman plots showed poor agreement between methods. When assessing the %change in moment with a lateral wedge, a moderate correlation was found (r=0.55, p=0.03) between methods; Bland-Altman plots showed moderate agreement between methods; and the lever-arm method (-6.4%) was not significantly different from the inverse dynamics method (-11.4%, p=0.09); however, the two methods produced opposite results 31% of the time.

CONCLUSION

The lever-arm method cannot estimate peak EKAMs, and can only approximate the %change in moment induced by a lateral wedge; however, the error rate was 31%. Therefore, the lever-arm method is not recommended for use in its current form.

CLINICAL RELEVANCE

This study may help guide the development of a fast and simple method for determining EKAMs for individuals with knee osteoarthritis.

A statistically-augmented computational platform for evaluating meniscal function

Meniscal implants have been developed in an attempt to provide pain relief and prevent pathological degeneration of articular cartilage. However, as yet there has been no systematic and comprehensive analysis of the effects of the meniscal design variables on meniscal function across a wide patient population, and there are no clear design criteria to ensure the functional performance of candidate meniscal implants. Our aim was to develop a statistically-augmented, experimentally-validated, computational platform to assess the effect of meniscal properties and patient variables on knee joint contact mechanics during the activity of walking. Our analysis used Finite Element Models (FEMs) that represented the geometry, kinematics as based on simulated gait and contact mechanics of three laboratory tested human cadaveric knees. The FEMs were subsequently programmed to represent prescribed meniscal variables (circumferential and radial/axial moduli-Ecm, Erm, stiffness of the meniscal attachments-Slpma, Slamp) and patient variables (varus/valgus alignment-VVA, and articular cartilage modulus-Ec). The contact mechanics data generated from the FEM runs were used as training data to a statistical interpolator which estimated joint contact data for untested configurations of input variables. Our data suggested that while Ecm and Erm of a meniscus are critical in determining knee joint mechanics in early and late stance (peak 1 and peak 3 of the gait cycle), for some knees that have greater laxity in the mid-stance phase of gait, the stiffness of the articular cartilage, Ec, can influence force distribution across the tibial plateau. We found that the medial meniscus plays a dominant load-carrying role in the early stance phase and less so in late stance, while the lateral meniscus distributes load throughout gait. Joint contact mechanics in the medial compartment are more sensitive to Ecm than those in the lateral compartment. Finally, throughout stance, varus-valgus alignment can overwhelm these relationships while the stiffness of meniscal attachments in the range studied have minimal effects on the knee joint mechanics. In summary, our statistically-augmented, computational platform allowed us to study how meniscal implant design variables (which can be controlled at the time of manufacture or implantation) interact with patient variables (which can be set in FEMs but cannot be controlled in patient studies) to affect joint contact mechanics during the activity of simulated walking.

Material models and properties in the finite element analysis of knee ligaments: a literature review

Knee ligaments are elastic bands of soft tissue with a complex microstructure and biomechanics, which are critical to determine the kinematics as well as the stress bearing behavior of the knee joint. Their correct implementation in terms of material models and properties is therefore necessary in the development of finite element models of the knee, which has been performed for decades for the investigation of both its basic biomechanics and the development of replacement implants and repair strategies for degenerative and traumatic pathologies. Indeed, a wide range of element types and material models has been used to represent knee ligaments, ranging from elastic unidimensional elements to complex hyperelastic three-dimensional structures with anatomically realistic shapes. This paper systematically reviews literature studies, which described finite element models of the knee, and summarizes the approaches, which have been used to model the ligaments highlighting their strengths and weaknesses.