Knee medial and lateral contact forces in a musculoskeletal model with subject-specific contact point trajectories

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.

Can a Musculoskeletal Model Adapted to Knee Implant Geometry Improve Prediction of 3D Contact Forces and Moments?

Tibiofemoral contact loads are crucial parameters in the onset and progression of osteoarthrosis. While contact loads are frequently estimated from musculoskeletal models, their customization is often limited to scaling musculoskeletal geometry or adapting muscle lines. Moreover, studies have usually focused on superior-inferior contact force without investigating three-dimensional contact loads. Using experimental data from six patients with instrumented total knee arthroplasty (TKA), this study customized a lower limb musculoskeletal model to consider the positioning and the geometry of the implant at knee level. Static optimization was performed to estimate tibiofemoral contact forces and contact moments as well as musculotendinous forces. Predictions from both a generic and a customized model were compared to the instrumented implant measurements. Both models accurately predict superior-inferior (SI) force and abduction-adduction (AA) moment. Notably, the customization improves prediction of medial-lateral (ML) force and flexion-extension (FE) moments. However, there is subject-dependent variability in the prediction of anterior-posterior (AP) force. The customized models presented here predict loads on all joint axes and in most cases improve prediction. Unexpectedly, this improvement was more limited for patients with more rotated implants, suggesting a need for further model adaptations such as muscle wrapping or redefinition of hip and ankle joint centers and axes.

A subject-specific musculoskeletal model to predict the tibiofemoral contact forces during daily living activities

Accurate prediction of tibiofemoral contact force (TFCF) during daily living activities is significant for understanding the initiation, progression, and treatment of knee osteoarthritis (KOA). However, the diversity of target activities, prediction accuracy, and computational efficiency of the current musculoskeletal simulations need to be further improved. In this study, a subject-specific musculoskeletal model considered the tibiofemoral alignment, medial-lateral contact locations, secondary tibiofemoral and all patellofemoral motions, and knee ligaments was proposed to predict the TFCFs during the five daily activities (normal walking, sit-to-stand, stand-to-sit, stair ascent, and stair descent) in OpenSim software. The standing lower-limbs-full-length radiograph, local radiograph of knee joint, motion capture data, and force plate data of eighteen subjects were acquired as the input data of the musculoskeletal model. The results showed good agreements of TFCFs between the predictions based on our proposed musculoskeletal model and the in-vivo measurements based on instrumented knee implants during the five daily activities (RMSE: 0.16 ∼ 0.31 BW, R²: 0.88 ∼ 0.97, M: -0.11 ∼ -0.02, P: 0.03 ∼ 0.10, and C: 0.04 ∼ 0.14). Additionally, the order of the peak total and lateral TFCFs from low to high was normal walking, stair ascent and stand-to-sit, and stair descent and sit-to-stand (P < 0.05), and the peak medial TFCF was stand-to-sit, sit-to-stand, normal walking, stair ascent and stair descent (P < 0.05). The outcomes of this study are valuable for further understanding the knee biomechanics during daily living activities and providing theoretical guidance for the treatments of KOA.

An Efficient One-Step Moment Balancing Algorithm for Computing Medial and Lateral Knee Compartment Contact Forces

The knee adduction moment is associated with the progression of knee osteoarthritis (OA). The adduction moment reflects the net effect of muscles, passive tissues and bone-on-bone contact forces. Medial compartment OA is more common than lateral and therefore our ability to correctly partition bone-on-bones forces across the medial and lateral compartments is key to understanding mechanical factors associated with the onset and progression of knee OA. In this technical brief we present an efficient one-step moment balancing algorithm linking the sagittal and frontal planes in the determination of musculotendon forces. Novel to the one-step approach is the introduction of a penalty function limiting total compressive force from acting in the lateral compartment when the internal moment is net abduction (i.e., external knee adduction). Medial and lateral knee contact forces were computed using the one-step moment balancing algorithm for 10 subjects walking at a self-selected pace and compared to values determined using a well-established two-step frontal moment balancing approach. Overall, average peak differences in magnitude and timing were small and the ensemble-averaged contact force profiles were similar between methods. The only statistical difference was slightly larger (0.2 BWs) peak medial contact force for the one-step method during the first half of stance, however these differences are small relative to peak values and would not likely alter interpretation of the data. The 1-step moment balancing method is a more efficient methodology for computing medial and lateral knee contact forces that can be used in place of two-step frontal plane moment balancing.

Stresses and deformations of an osteosynthesis plate in a lateral tibia plateau fracture

This study has the aim to investigate the strain and stress in an anterolateral locking plate applied for the fixation of a lateral split fracture. To simulate a complex fracture situation, three segments were separated. With a FEM analysis, representative places for strain and stress measurement were determined. A locked osteosynthesis plate was instrumented with strain gauges and tested on a fractured and a non-fractured Saw Bone model. To simulate different loading situations, four different points of force application, from the center of the condyles to a 15 mm posterior position, were used with a medial-lateral load distribution of 60:40. The simulations as well as the biomechanical tests demonstrated that two deformations dominate the load on the plate: a bending into posterior direction and a bulging of the plate head. Shifting the point of application to the posterior direction resulted in increasing maximum stress, from 1.16 to 6.32 MPa (FEM analysis) and from 3.04 to 7.00 MPa (biomechanical study), respectively. Furthermore, the comparison of the non-fractured and fractured models showed an increase in maximum stress by the factor 2.06-2.2 (biomechanical investigation) and 1.5-3.3 (FEM analysis), respectively.

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.

Accuracy of the tibiofemoral contact forces estimated by a subject-specific musculoskeletal model with fluoroscopy-based contact point trajectories

Accurate estimation of the tibiofemoral contact forces relies on exact kinematics and joint geometry. Subject-specific kinematic constraints representing contact point trajectories derived from fluoroscopic measurements during lunge are introduced in a musculoskeletal model of the lower limb and compared to generic kinematic constraints. The medial, lateral, and total contact forces during gait and squat are validated using the data of four patients with an instrumented prosthesis. The accuracy of the estimated contact forces (both with subject-specific and generic kinematic constraints) remains close to the level reported in the literature. The mean root mean square errors range from 0.32 to 0.52 body weights for gait and from 0.27 to 0.72 body weights for squat. The impact of the subject-specific contact point trajectories is not found substantial or consistent between patients and tasks. Indeed, the kinematics of the total knee prostheses remains close to the kinematics of a hinge joint and the contact point locations remain generally centred at 20 mm from the tibia centreline (close to the constant value defined in the generic constraints). The contact point trajectories are also suspected to differ between tasks (lunge vs. gait and squat). While the contact point trajectories have been reported to be sensitive model parameters, no clear improvement of the contact force accuracy is demonstrated on patients with instrumented prosthesis. The introduction (as kinematic constraints) of fluoroscopy-based contact point trajectories may be considered in cases where these trajectories are significantly altered, as reported for osteoarthritis patients.

A method for quantitative evaluation of a valgus knee orthosis using biplane x-ray images

Knee orthoses are designed to reestablish the normal kinematics of the knee joint. However, the data on the effectiveness of them on modifying the internal joint kinematics are scarce. The aim of this study was to develop a method to allow accurate comparison of the knee contact kinematics in osteoarthritic (OA) subjects with and without wearing a valgus knee orthosis using imaging techniques. Biplane x-ray images of a subject (68 yrs., female, 1.70 m, 89 kg, left knee) was recorded during a weight-bearing squat at five positions. The same squat trial was repeated while wearing the orthosis. The 3D models of the knee were reconstructed from the biplane x-rays and the joint kinematics as well as the tibiofemoral contact point locations and bone-to-bone distance were compared at each posture. This could be seen as a proof of concept for the use of contact point locations as a parameter for evaluating the effectiveness of knee orthoses.Clinical Relevance- Joint kinematics derived from the skin markers suffer from low accuracy. The real impact of the knee orthoses on the skeleton takes vigorous techniques, which allows detecting the subtle kinematics changes directly at the joint level.

Reconstruction of the lower limb bones from digitised anatomical landmarks using statistical shape modelling

BACKGROUND

Bone shapes strongly influence force and moment predictions of kinematic and musculoskeletal models used in motion analysis. The precise determination of joint reference frames is essential for accurate predictions. Since clinical motion analysis typically does not include medical imaging, from which bone shapes may be obtained, scaling methods using reference subjects to create subject-specific bone geometries are widely used.

RESEARCH QUESTION

This study investigated if lower limb bone shape predictions from skin-based measurements, utilising an underlying statistical shape model (SSM) that corrects for soft tissue artefacts in digitisation, can be used to improve conventional linear scaling methods of bone geometries.

METHODS

SSMs created from 35 healthy adult femurs and tibiae/fibulae were used to reconstruct bone shapes by minimising the distance between anatomical landmarks on the models and those digitised in the motion laboratory or on medical images. Soft tissue artefacts were quantified from magnetic resonance images and then used to predict distances between landmarks digitised on the skin surface and bone. Reconstruction results were compared to linearly scaled models by measuring root mean squared distances to segmented surfaces, calculating differences of commonly used anatomical measures and the errors in the prediction of the hip joint centre.

RESULTS

SSM reconstructed surface predictions from varying landmark sets from skin and bone landmarks were more accurate compared to linear scaling methods (2.60-2.95 mm vs. 3.66-3.87 mm median error; p < 0.05). No significant differences were found between SSM reconstructions from bony landmarks and SSM reconstructions from digitised landmarks obtained in the motion lab and therefore reconstructions using skin landmarks are as accurate as reconstructions from landmarks obtained from medical images.

SIGNIFICANCE

These results indicate that SSM reconstructions can be used to increase the accuracy in obtaining bone shapes from surface digitised experimental data acquired in motion lab environments.