Prediction and Validation of Load-Dependent Behavior of the Tibiofemoral and Patellofemoral Joints During Movement

The impact of PCA derived gait kinematic variations on estimated medial knee contact forces in a knee osteoarthritis population

Osteoarthritis (OA) is a prevalent musculoskeletal condition leading to functional limitations, especially among the elderly. Current treatments focus on pain relief and functional improvement, however there is a lack of approaches which slow disease progression. A promising approach focusses on reducing knee joint loading, as excessive loading contributes to knee OA progression. This study explores kinematic variations in the knee OA population, utilizing principal component analysis (PCA) to examine gait variations (primitives) in both healthy individuals and those with knee osteoarthritis (KOA) and their implications for knee joint loading. The KOA population exhibited 14 modes of variation representing 95% of the cumulative variance, compared to 20 in the healthy population, indicating lower variability with KOA. The relation between identified gait primitives and knee loading parameters, revealed complex relationships. Surprisingly, modes with the largest kinematic variations did not consistently correspond to the highest variations in knee loading parameters revealing degrees of freedom which may have a larger role in determining joint loading. Moreover, potential gait-retraining strategies for KOA, associating specific kinematic combinations with altered knee loading were identified. The results showed a good agreement with previously applied strategies. However, this study highlights the importance of analyzing whole-body kinematics for effective gait retraining, as opposed to focusing on one single joint variation. The study's insights contribute to understanding the intricate interplay between gait pattern variations and knee joint loading changes in healthy and KOA populations, offering practical applications for guiding interventions and estimating loading parameters.

Effect of uncertainties in musculoskeletal modeling inputs on sensitivity of knee joint finite element simulations

Musculoskeletal finite element modeling is used to estimate mechanical responses of knee joint tissues but involves uncertainties in muscle activations, marker locations, cartilage stiffness, maximum isometric forces, and gait parameter personalization. This study investigates how these uncertainties affect cartilage mechanical responses in knee joint finite element models during walking. We selected three subjects and constructed five musculoskeletal models for each, representing different variations of modeling assumptions, along with a reference model using conventional assumptions. We then ran finite element simulations of knee joints using both personalized gait inputs (motion and loading boundary conditions) and non-personalized gait inputs from literature. Our results demonstrated that varying modeling assumptions, such as optimization function for muscle activation patterns, knee marker position, knee cartilage stiffness, and maximum isometric force, produced highly subject-specific effects. Differences between the reference and altered models ranged from 3% to 30% in musculoskeletal modeling and from 1% to 61% in finite element modeling results. The largest effects occurred with non-personalized gait data, resulting in up to 6- and 2-fold changes in musculoskeletal and finite element modeling results, respectively. This study highlights the sensitivity of knee mechanics to different modeling assumptions and underscores the importance of applying personalized gait parameters for accurate finite element simulations.

Limb Underloading in Walking Transmits Less Dynamic Knee Joint Contact Forces after Anterior Cruciate Ligament Reconstruction

INTRODUCTION

Individuals with anterior cruciate ligament reconstruction (ACLR) often walk with a less dynamic vertical ground reaction force (vGRF), exemplified by a reduced first peak vGRF and elevated midstance vGRF compared with uninjured controls. However, the mechanism by which altered limb loading affects actual tibial plateau contact forces during walking remains unclear.

METHODS

Our purpose was to use musculoskeletal simulation to evaluate the effects of first peak vGRF biofeedback on bilateral tibiofemoral contact forces relevant to the development of post-traumatic osteoarthritis in 20 individuals with ACLR. We hypothesized that reduced first peak vGRF would produce less dynamic tibial plateau contact forces during walking in individuals with ACLR.

RESULTS

As the pivotal outcome from this study, and in support of our hypothesis, we found that less dynamic vGRF profiles in individuals with ACLR-observations that have associated in prior studies with more cartilage breakdown serum biomarkers and reduced proteoglycan density-are accompanied by less dynamic tibiofemoral joint contact forces during walking.

CONCLUSION

We conclude that more sustained limb-level loading, a phenotype that associates with worse knee joint health outcomes after ACLR and was prescribed herein using biofeedback, alters the loading profile and magnitude of force applied to tibiofemoral cartilage.

An Automated and Robust Tool for Musculoskeletal and Finite Element Modeling of the Knee Joint

OBJECTIVE

To develop and assess an automatic and robust knee musculoskeletal finite element (MSK-FE) modeling pipeline.

METHODS

Magnetic resonance images (MRIs) were used to train nnU-Net networks for auto-segmentation of knee bones (femur, tibia, patella, and fibula), cartilages (femur, tibia, and patella), menisci, and major knee ligaments. Two different MRI sequences were used to broaden applicability. Next, we created MSK-FE models of an unseen dataset using two MSK-FE modeling pipelines: template-based and auto-meshing. MSK models had personalized knee geometries with multi-degree-of-freedom elastic foundation contacts. FE models used fibril-reinforced poroviscoelastic swelling material models for cartilages and menisci.

RESULTS

Volumes of knee bones, cartilages, and menisci did not significantly differ (p>0.05) across MRI sequences. MSK models estimated secondary knee kinematics during passive knee flexion tests consistent with in vivo and simulation-based values from the literature. Between the template-based and auto-meshing FE models, estimated cartilage mechanics often differed significantly (p<0.05), though differences were <15% (considering peaks during walking), i.e., <1.5 MPa for maximum principal stress, <1 percentage point for collagen fibril strain, and <3 percentage points for maximum shear strain.

CONCLUSION

The template-based modeling provided a more rapid and robust tool than the auto-meshing approach, while the estimated knee biomechanics were comparable. Nonetheless, the auto-meshing approach might provide more accurate estimates in subjects with distinct knee irregularities, e.g., cartilage lesions.

SIGNIFICANCE

The MSK-FE modeling tool provides a rapid, easy-to-use, and robust approach for investigating task- and person-specific mechanical responses of the knee cartilage and menisci, holding significant promise, e.g., in personalized rehabilitation planning.

An open-source framework for the generation of OpenSim models with personalised knee joint geometries for the estimation of articular contact mechanics

Musculoskeletal modelling pipelines typically use generic models scaled to individual's anthropometry. The ability to represent variations in bone or joint geometry and alignment is highly limited. This may have a large effect, particularly when modelling contact between articular surfaces such as for the knee where articular contact mechanics are used to determine joint kinematics and the resulting cartilage contact pressures and locations. Here we describe a developed open-source framework for the personalisation of such models and compare dynamic simulation outputs. The framework involves three main steps: (1) positions personalised geometries from magnetic resonance imaging and replaces generic bone and contact geometries. (2) Repositions muscle and ligament attachments and via points and optimisation of wrapping surfaces to ensure physiological lengthening behaviour. Finally, (3) muscle and ligament properties are calibrated to ensure physiological behaviour. Following model creation, dynamic simulations from a single participant and gait trial were compared. Small changes in knee adduction/abduction and rotation angles were observed between models. Joint moment differences however were present in not only the knee but also hip and ankle joints. These differences resulted in changes in both the magnitude and location of knee joint contact pressure. The framework developed is automated and requires only minimal user interaction and is built using open-source software packages which can be freely downloaded and installed. The adoption of such personalised modelling approaches facilitates patient specific modelling and may provide more detailed information regarding disease progression, patient stratification and facilitate personalised rehabilitation and treatment planning.

On the consequences of intra-operative release versus over-tensioning of the posterior cruciate ligament in total knee arthroplasty

Intra-operative tensioning of the posterior cruciate ligament (PCL) in total knee arthroplasty (TKA) is commonly based on the surgeon's experience, resulting in a possibly loose or overly tight PCL. To date, the consequences of different PCL tensioning scenarios for the post-operative biomechanics of the knee remain unclear. Using a comprehensive musculoskeletal modelling approach that allows predictive joint kinematic and kinetic balance, we assessed variations in the movement and loading patterns of the knee as well as changes in ligament and muscle forces during walking in response to systematic variations in the PCL reference strain. The results indicate only small differences in the tibiofemoral and patellofemoral kinematics and kinetics for scenarios involving up to 10% release of the PCL (relative to the baseline reference scenario with 2% residual strain). These observations remain valid for simulations performed with high- as well as with low-conformity implant designs. However, over-tensioning of the ligament was found to considerably overload the tibiofemoral joint, including altered contact mechanics, and may therefore shorten the implant longevity. Finally, no meaningful impact of the PCL reference strain on the muscle force patterns was observed. This study therefore favours balancing the knee with a slightly loose rather than tense PCL, if appropriate intra-operative PCL tension cannot be objectively achieved.

Patient-specific gait pattern in individuals with patellofemoral instability reduces knee joint loads

Patellofemoral instability is influenced by morphological factors and associated with compensational alterations in gait pattern. Recent simulation studies investigated the impact of knee morphology on the stability and loading of the patellofemoral joint but neglected the patient-specific gait pattern. The aim of this study was to investigate the impact of patient-specific gait pattern on muscle forces and joint loading in individuals with patellofemoral instability. Musculoskeletal simulations with a model including a twelve degrees of freedom knee joint were performed based on three-dimensional motion capture data of 21 individuals with chronic patellofemoral instability and 17 healthy control participants. The patellofemoral instability group walked with a less flexed knee joint and reduced knee flexion and abduction moments compared to the control group, which required less quadriceps muscle forces. Lower quadriceps muscle forces resulted in a reduction of tibiofemoral and patellofemoral joint contact forces despite similar walking velocities between both groups. Furthermore, we observed decreased lateralizing patella forces in subjects with patella instability, which could potentially reduce the risk of patella dislocation. Our findings highlight the importance of accounting for the patient-specific gait pattern when analysing knee loads in individuals with patellofemoral instability.

Methodology to identify subject-specific dynamic laxity tests to stretch individual parts of knee ligaments

The mechanical properties of ligaments are important for multiple applications and are often estimated from laxity tests. However, the typical laxity tests are not optimized for this application and, a potential exists to develop better laxity tests in this respect. Therefore, the purpose of this study was to develop a methodology to identify optimal, dynamic laxity tests that isolate the stretch of the individual ligaments from each other. To this end, we applied an existing rigid body-based knee model and a dataset of ∼100.000 random samples of applied forces (0-150 N), moments (0-10 Nm) and knee flexion angles (0-90°) through Monte Carlo Simulations. For each modelled ligament bundle, we identified ten load cases; one producing the highest force and nine equally spaced between the maximal and zero force, where the maximal force in all other ligament bundles were minimized. We compared these novel laxity tests to standard internal/external and varus/valgus laxity tests using an isolation metric. We found that no laxity test could stretch the anterior part of the posterior cruciate and medial cruciate ligaments (PCL and MCL), whereas for all other ligaments, except the posterior PCL, the new laxity tests isolated the ligament stretch 28 % to 450 % better than standard tests. From our study, we conclude that it is possible to define better laxity tests than currently exist and these may be highly relevant for determination of mechanical properties of ligaments in vivo. Future studies should generalize our results and translate them to modern laxity measurements technologies.

Population-based in silico modeling of anatomical shape variation of the knee and its impact on joint loading in knee osteoarthritis

Anatomical knee joint features and osteoarthritis (OA) severity are associated, however confirming causals link to altered knee loading is challenging. This study leverages statistical shape models (SSM) to investigate the relationship between joint shape/alignment and knee loading during gait in knee OA (KOA) patients to understand their contribution to elevated medial knee loading in OA. Musculoskeletal (MSK) models were created for the mean as well as the first eight SSM principal modes of variation (-3,-2,-1, +1, +2, +3 standard deviations for each mode) and used as input to a MSK modeling framework. Using an identical KOA gait pattern (i.e., joint kinematics and ground reaction forces), we ran simulations for each MSK model and evaluated medial compartment loading magnitude and contact distribution at the instant of first and second peak of knee joint loading. An increase in external rotation, posterior tibia translation and a decrease in medial joint space and medial femoral condylar size predisposed the medial compartment knee joint to overloading during gait. This was coupled with an anterior and medial shift in contact location with increasing external rotated tibial position and increasing posterior tibial translation with respect to the femur. Next, results also highlighted a posterior shift of the medial compartment loading location with decreasing medial joint space. This study provides important population-based insights on how knee shape and alignment predispose individuals with KOA to elevated medial compartmental knee loading. This information can be crucial in assessing the risk for medial KOA development and progression.

Hindfoot kinematics and kinetics - A combined in vivo and in silico analysis approach

BACKGROUND

The complex anatomical structure of the foot-ankle imposes challenges to accurately quantify detailed hindfoot kinematics and estimate musculoskeletal loading parameters. Most systems used to capture or estimate dynamic joint function oversimplify the anatomical structure by reducing its complexity.

RESEARCH QUESTION

Can four dimensional computed tomography (4D CT) imaging in combination with an innovative foot manipulator capture in vivo hindfoot kinematics during a simulated stance phase of walking and can talocrural and subtalar articular joint mechanics be estimated based on a detailed in silico musculoskeletal foot-ankle model.

METHODS

A foot manipulator imposed plantar/dorsiflexion and inversion/eversion representing a healthy stance phase of gait in 12 healthy participants while simultaneously acquiring 4D CT images. Participant-specific 3D hindfoot rotations and translations were calculated based on bone-specific anatomical coordinate systems. Articular cartilage contact area and contact pressure of the talocrural and subtalar joints were estimated using an extended foot-ankle model updated with an elastic foundation contact model upon prescribing the participant-specific rotations measured in the 4D CT measurement.

RESULTS

Plantar/dorsiflexion predominantly occurred at the talocrural joint (RoM 15.9±3.9°), while inversion/eversion (RoM 5.9±3.9°) occurred mostly at the subtalar joint, with the contact area being larger at the subtalar than at the talocrural joint. Contact pressure was evenly distributed between the talocrural and subtalar joint at the beginning of the simulated stance phase but was then redistributed from the talocrural to the subtalar joint with increasing dorsiflexion.

SIGNIFICANCE

In a clinical case study, the healthy participants were compared with four patients after surgically treaded intra-articular calcaneal fracture. The proposed workflow was able to detect small but meaningful differences in hindfoot kinematics and kinetics, indicative of remaining hindfoot pathomechanics that may influence the onset and progression of degenerative joint diseases.

Knee kinematics are primarily determined by implant alignment but knee kinetics are mainly influenced by muscle coordination strategy

Implant malalignment has been reported to be a primary reason for revision total knee arthroplasty (TKA). In addition, altered muscle coordination patterns are commonly observed in TKA patients, which is thought to alter knee contact loads. A comprehensive understanding of the influence of surgical implantation and muscle recruitment strategies on joint contact mechanics is crucial to improve surgical techniques, increase implant longevity, and inform rehabilitation protocols. In this study, a detailed musculoskeletal model with a 12 degrees of freedom knee was developed to represent a TKA subject from the CAMS-Knee datasets. Using motion capture and ground reaction force data, a level walking cycle was simulated and the joint movement and loading patterns were estimated using a novel technique for concurrent optimization of muscle activations and joint kinematics. In addition, over 12'000 Monte Carlo simulations were performed to predict knee contact mechanics during walking, considering numerous combinations of implant alignment and muscle activation scenarios. Validation of our baseline simulation showed good agreement between the model kinematics and loading patterns against the in vivo data. Our analyses reveal a considerable impact of implant alignment on the joint kinematics, while variation in muscle activation strategies mainly affects knee contact loading. Moreover, our results indicate that high knee compressive forces do not necessarily originate from extreme kinematics and vice versa. This study provides an improved understanding of the complex inter-relationships between loading and movement patterns resulting from different surgical implantation and muscle coordination strategies and presents a validated framework towards population-based modelling in TKA.

Prediction of Model Generated Patellofemoral Joint Contact Forces Using Principal Component Prediction and Reconstruction

It is not currently possible to directly and noninvasively measure in vivo patellofemoral joint contact force during dynamic movement; therefore, indirect methods are required. Simple models may be inaccurate because patellofemoral contact forces vary for the same knee flexion angle, and the patellofemoral joint has substantial out-of-plane motion. More sophisticated models use 3-dimensional kinematics and kinetics coupled to a subject-specific anatomical model to predict contact forces; however, these models are time consuming and expensive. We applied a principal component analysis prediction and regression method to predict patellofemoral joint contact forces derived from a robust musculoskeletal model using exclusively optical motion capture kinematics (external approach), and with both patellofemoral and optical motion capture kinematics (internal approach). We tested this on a heterogeneous population of asymptomatic subjects (n = 8) during ground-level walking (n = 12). We developed equations that successfully capture subject-specific gait characteristics with the internal approach outperforming the external. These approaches were compared with a knee-flexion based model in literature (Brechter model). Both outperformed the Brechter model in interquartile range, limits of agreement, and the coefficient of determination. The equations generated by these approaches are less computationally demanding than a musculoskeletal model and may act as an effective tool in future rapid gait analysis and biofeedback applications.

Validation and evaluation of subject-specific finite element models of the pediatric knee

Finite element (FE) models have been widely used to investigate knee joint biomechanics. Most of these models have been developed to study adult knees, neglecting pediatric populations. In this study, an atlas-based approach was employed to develop subject-specific FE models of the knee for eight typically developing pediatric individuals. Initially, validation simulations were performed at four passive tibiofemoral joint (TFJ) flexion angles, and the resulting TFJ and patellofemoral joint (PFJ) kinematics were compared to corresponding patient-matched measurements derived from magnetic resonance imaging (MRI). A neuromusculoskeletal-(NMSK)-FE pipeline was then used to simulate knee biomechanics during stance phase of walking gait for each participant to evaluate model simulation of a common motor task. Validation simulations demonstrated minimal error and strong correlations between FE-predicted and MRI-measured TFJ and PFJ kinematics (ensemble average of root mean square errors < 5 mm for translations and < 4.1° for rotations). The FE-predicted kinematics were strongly correlated with published reports (ensemble average of Pearson's correlation coefficients (ρ) > 0.9 for translations and ρ > 0.8 for rotations), except for TFJ mediolateral translation and abduction/adduction rotation. For walking gait, NMSK-FE model-predicted knee kinematics, contact areas, and contact pressures were consistent with experimental reports from literature. The strong agreement between model predictions and experimental reports underscores the capability of sequentially linked NMSK-FE models to accurately predict pediatric knee kinematics, as well as complex contact pressure distributions across the TFJ articulations. These models hold promise as effective tools for parametric analyses, population-based clinical studies, and enhancing our understanding of various pediatric knee injury mechanisms. They also support intervention design and prediction of surgical outcomes in pediatric populations.

The History and Future of Neuromusculoskeletal Biomechanics

The Executive Council of the International Society of Biomechanics has initiated and overseen the commemorations of the Society's 50th Anniversary in 2023. This included multiple series of lectures at the ninth World Congress of Biomechanics in 2022 and XXIXth Congress of the International Society of Biomechanics in 2023, all linked to special issues of International Society of Biomechanics' affiliated journals. This special issue of the Journal of Applied Biomechanics is dedicated to the biomechanics of the neuromusculoskeletal system. The reader is encouraged to explore this special issue which comprises 6 papers exploring the current state-of the-art, and future directions and roles for neuromusculoskeletal biomechanics. This editorial presents a very brief history of the science of the neuromusculoskeletal system's 4 main components: the central nervous system, musculotendon units, the musculoskeletal system, and joints, and how they biomechanically integrate to enable an understanding of the generation and control of human movement. This also entails a quick exploration of contemporary neuromusculoskeletal biomechanics and its future with new fields of application.

Signatures of disease progression in knee osteoarthritis: insights from an integrated multi-scale modeling approach, a proof of concept

Introduction: Knee osteoarthritis (KOA) is characterized by articular cartilage degeneration. It has been widely accepted that the mechanical joint environment plays a significant role in the onset and progression of this disease. In silico models have been used to study the interplay between mechanical loading and cartilage degeneration, hereby relying mainly on two key mechanoregulatory factors indicative of collagen degradation and proteoglycans depletion. These factors are the strain in collagen fibril direction (SFD) and maximum shear strain (MSS) respectively. Methods: In this study, a multi-scale in silico modeling approach was used based on a synergy between musculoskeletal and finite element modeling to evaluate the SFD and MSS. These strains were evaluated during gait based on subject-specific gait analysis data collected at baseline (before a 2-year follow-up) for a healthy and progressive early-stage KOA subject with similar demographics. Results: The results show that both SFD and MSS factors allowed distinguishing between a healthy subject and a KOA subject, showing progression at 2 years follow-up, at the instance of peak contact force as well as during the stance phase of the gait cycle. At the peak of the stance phase, the SFD were found to be more elevated in the KOA patient with the median being 0.82% higher in the lateral and 0.4% higher in the medial compartment of the tibial cartilage compared to the healthy subject. Similarly, for the MSS, the median strains were found to be 3.6% higher in the lateral and 0.7% higher in the medial tibial compartment of the KOA patient compared to the healthy subject. Based on these intersubject SFD and MSS differences, we were additionally able to identify that the tibial compartment of the KOA subject at risk of progression. Conclusion/discussion: We confirmed the mechanoregulatory factors as potential biomarkers to discriminate patients at risk of disease progression. Future studies should evaluate the sensitivity of the mechanoregulatory factors calculated based on this multi-scale modeling workflow in larger patient and control cohorts.

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.

Peak Tibiofemoral Contact Forces Estimated Using IMU-Based Approaches Are Not Significantly Different from Motion Capture-Based Estimations in Patients with Knee Osteoarthritis

Altered tibiofemoral contact forces represent a risk factor for osteoarthritis onset and progression, making optimization of the knee force distribution a target of treatment strategies. Musculoskeletal model-based simulations are a state-of-the-art method to estimate joint contact forces, but they typically require laboratory-based input and skilled operators. To overcome these limitations, ambulatory methods, relying on inertial measurement units, have been proposed to estimated ground reaction forces and, consequently, knee contact forces out-of-the-lab. This study proposes the use of a full inertial-capture-based musculoskeletal modelling workflow with an underlying probabilistic principal component analysis model trained on 1787 gait cycles in patients with knee osteoarthritis. As validation, five patients with knee osteoarthritis were instrumented with 17 inertial measurement units and 76 opto-reflective markers. Participants performed multiple overground walking trials while motion and inertial capture methods were synchronously recorded. Moderate to strong correlations were found for the inertial capture-based knee contact forces compared to motion capture with root mean square error between 0.15 and 0.40 of body weight. The results show that our workflow can inform and potentially assist clinical practitioners to monitor knee joint loading in physical therapy sessions and eventually assess long-term therapeutic effects in a clinical context.

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.

Patella height influences patellofemoral contact and kinematics following cruciate-retaining total knee replacement

The role of patella height is discussed controversially in total knee arthroplasty (TKA). Therefore, this computational study aims to systematically analyze the biomechanical effect of different patella heights on patellofemoral (PF) forces and kinematics after cruciate-retaining (CR) TKA. We implemented a CR bicondylar TKA with a dome patellar button in a validated dynamic musculoskeletal multibody model of a male human knee joint. Retropatellar dynamics (contact force [N], shear force [N], patellar shift [mm], tilt [°], and rotation [°]) were evaluated during dual-limb squat motion (flexion from 0° to 90°) with simulated active muscle forces and the effects of different patella heights (Blackburne-Peel [BP] ratio of 0.39, 0.49, 0.65, 0.85, 1.01, and 1.1 were systematically examined). As active knee flexion increased, PF contact force also increased. Patella alta (BP = 1.1) resulted in higher PF contact forces compared to normal patella height (BP = 0.65) by up to 16%. Contrarily, patella baja was associated with decreased PF forces by 7%. Compared to patella baja (BP = 0.39), patella alta (BP = 1.1) considerably increased the contact force by up to 25%. Different patellar heights mainly affected PF shear forces during early knee flexion. Concerning PF kinematics, patella alta (BP = 1.1) yielded a greater lateral tilt of more than 4° and higher patellar rotation by up to 3° during deep knee flexion, compared to normal patella height (BP = 0.65). Our computational study indicates that patella alta is associated with the highest PF contact and shear force after the implantation of a CR bicondylar TKA. This should be considered in PF disorders following TKA.

Determining the relationship between tibiofemoral geometry and passive motion with partial least squares regression

Tibiofemoral geometry influences knee passive motion and understanding their relationship can provide insight into knee function and mechanisms of injury. However, the complexity of the geometric constraints has made characterizing the relationship challenging. The aim of this study was to determine the tibiofemoral bone geometries that explain the variation in passive motion using a partial least squares regression (PLSR) model. The PLSR model was developed for 29 healthy cadaver specimens (10 female, 19 male) with femur and tibia geometries retrieved from MRI images and six degree-of-freedom tibiofemoral kinematics determined during a flexion cycle with minimal medial pressure. The first 13 partial least squares (PLS) components explained 90% of the variation in the kinematics and accounted for 89% of the variation in geometry. The first three PLS components which shared geometric changes to particular surface congruencies of the tibial and femoral condyles explained the most amount of variation in the kinematics, primarily in anterior-posterior translation. Meanwhile, variations in femoral condyle width and the intercondylar space, tibia plateau size and conformity, and tibia eminences heights in PLS 2 and 4 explained the greatest amount of variation in internal-external rotation. PLS 4 exhibiting variation in overall size of the knee accounted for greatest amount of variation in geometry (50%) and had the greatest influence on the abduction-adduction motion and some on internal-external rotation but, overall, explained only a small proportion of the kinematics (10%). Elucidating the complex relationship between tibiofemoral bone geometry and passive kinematics may help personalize treatments for improved functional outcomes in patients.

Simulating Knee-Stress Distribution Using a Computed Tomography-Based Finite Element Model: A Case Study

This study aimed to evaluate the mechanism of progression involved in knee osteoarthritis (OA). We used the computed tomography-based finite element method (CT-FEM) of quantitative X-ray CT imaging to calculate and create a model of the load response phase, wherein the greatest burden is placed on the knee joint while walking. Weight gain was simulated by asking a male individual with a normal gait to carry sandbags on both shoulders. We developed a CT-FEM model that incorporated walking characteristics of individuals. Upon simulating changes owing to a weight gain of approximately 20%, the equivalent stress increased extensively in both medial and lower leg aspects of the femur and increased medio-posteriorly by approximately 230%. As the varus angle increased, stress on the surface of the femoral cartilage did not change significantly. However, the equivalent stress on the surface of the subchondral femur was distributed over a wider area, increasing by approximately 170% in the medio-posterior direction. The range of equivalent stress affecting the lower-leg end of the knee joint widened, and stress on the posterior medial side also increased significantly. It was reconfirmed that weight gain and varus enhancement increase knee-joint stress and cause the progression of OA.

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.

A Musculoskeletal Model for Estimating Hip Contact Pressure During Walking

Cartilage contact pressures are major factors in osteoarthritis etiology and are commonly estimated using finite element analysis (FEA). FEA models often include subject-specific joint geometry, but lack subject-specific joint kinematics and muscle forces. Musculoskeletal models use subject-specific kinematics and muscle forces but often lack methods for estimating cartilage contact pressures. Our objective was to adapt an elastic foundation (EF) contact model within OpenSim software to predict hip cartilage contact pressures and compare results to validated FEA models. EF and FEA models were built for five subjects. In the EF models, kinematics and muscle forces were applied and pressure was calculated as a function of cartilage overlap depth. Cartilage material properties were perturbed to find the best match to pressures from FEA. EF models with elastic modulus = 15 MPa and Poisson's ratio = 0.475 yielded results most comparable to FEA, with peak pressure differences of 4.34 ± 1.98 MPa (% difference = 39.96 ± 24.64) and contact area differences of 3.73 ± 2.92% (% difference = 13.4 ± 11.3). Peak pressure location matched between FEA and EF for 3 of 5 subjects, thus we do not recommend this model if the location of peak contact pressure is critically important to the research question. Contact area magnitudes and patterns matched reasonably between FEA and EF, suggesting that this model may be useful for questions related to those variables, especially if researchers desire inclusion of subject-specific geometry, kinematics, muscle forces, and dynamic motion in a computationally efficient framework.

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.

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

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

Effects of the Anterolateral Ligament and Anterior Cruciate Ligament on Knee Joint Mechanics: A Biomechanical Study Using Computational Modeling

Background:

Recent studies on lateral knee anatomy have reported the presence of a true ligament structure, the anterolateral ligament (ALL), in the anterolateral region of the knee joint. However, its biomechanical effects have not been fully elucidated.

Purpose:

To investigate, by using computer simulation, the association between the ALL and anterior cruciate ligament (ACL) under dynamic loading conditions.

Study Design:

Descriptive laboratory study; Level of evidence, 5.

Methods:

The authors combined medical imaging from 5 healthy participants with motion capture to create participant-specific knee models that simulated the entire 12 degrees of freedom of tibiofemoral (TF) and patellofemoral (PF) joint behaviors. These dynamic computational models were validated using electromyographic data, muscle activation data, and data from previous experimental studies. Forces exerted on the ALL with ACL deficiency and on the ACL with ALL deficiency, as well as TF and PF contact forces with deficiencies of the ACL, ALL, and the entire ligament structure, were evaluated under gait and squat loading. A single gait cycle and squat cycle were divided into 11 time points (periods 0.0-1.0). Simulated ligament forces and contact forces were compared using nonparametric repeated-measures Friedman tests.

Results:

Force exerted on the ALL significantly increased with ACL deficiency under both gait- and squat-loading conditions. With ACL deficiency, the mean force on the ALL increased by 129.7% under gait loading in the 0.4 period ( P < .05) and increased by 189% under high flexion during the entire cycle of squat loading ( P < .05). A similar trend of significantly increased force on the ACL was observed with ALL deficiency. Contact forces on the TF and PF joints with deficiencies of the ACL, ALL, and entire ligament structure showed a complicated pattern. However, contact force exerted on TF and PF joints with respect to deficiencies of ACL and ALL significantly increased under both gait- and squat-loading conditions.

Conclusion:

The results of this computer simulation study indicate that the ACL and the ALL of the lateral knee joint act as secondary stabilizers to each other under dynamic load conditions.

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.

Influence of Articular Geometry and Tibial Tubercle Location on Patellofemoral Kinematics and Contact Mechanics

Trochlear groove geometry and the location of the tibial tubercle, where the patellar tendon inserts, have both been associated with patellofemoral instability and can be modified surgically. Although their effects on patellofemoral biomechanics have been investigated individually, the interaction between the two is unclear. The authors' aim was to use statistical shape modeling and musculoskeletal simulation to examine the effect of patellofemoral geometry on the relationship between tibial tubercle location and patellofemoral function. A statistical shape model was used to generate new knee geometries with trochlear grooves ranging from shallow to deep. A Monte Carlo approach was used to create 750 knee models by randomly selecting a geometry and randomly translating the tibial tubercle medially/laterally and anteriorly. Each knee model was incorporated into a musculoskeletal model, and an overground walking trial was simulated. Knees with shallow trochlear geometry were more sensitive to tubercle medialization with greater changes in lateral patella position (-3.0 mm/cm medialization shallow vs -0.6 mm/cm deep) and cartilage contact pressure (-0.51 MPa/cm medialization shallow vs 0.04 MPa/cm deep). However, knees with deep trochlear geometry experienced greater increases in medial cartilage contact pressure with medialization. This modeling framework has the potential to aid in surgical decision making.

A Computationally Efficient Lower Limb Finite Element Musculoskeletal Framework Directly Driven Solely by Inertial Measurement Unit Sensors

Finite element musculoskeletal (FEMS) approaches using concurrent musculoskeletal and finite element models driven by motion data such as marker-based motion trajectory can provide insight into the interactions between the knee joint secondary kinematics, contact mechanics, and muscle forces in subject-specific biomechanical investigations. However, these data-driven FEMS systems have two major disadvantages that make them challenging to apply in clinical environments: they are computationally expensive and they require expensive and inconvenient equipment for data acquisition. In this study, we developed an FEMS model of the lower limb driven solely by inertial measurement unit sensors that includes the tissue geometries of the entire knee joint and combines muscle modeling and elastic foundation theory-based contact analysis of knee into a single framework. The model requires only the angular velocities and accelerations measured by the sensors as input, and the target outputs (knee contact mechanics, secondary kinematics, and muscle forces) are predicted from the convergence results of iterative calculations of muscle force optimization and knee contact mechanics. To evaluate its accuracy, the model was compared with in vivo experimental data during gait. The maximum contact pressure (12.6 MPa) in the rigid body contact analysis occurred on the medial side of the cartilage at the maximum loading response. The proposed computationally efficient framework drastically reduced the computational time (97.5% reduction) in comparison with the conventional deformable finite element analysis. The developed framework combines measurement convenience and computational efficiency and shows promise for clinical applications.

Magnetic Resonance Imaging-based biomechanical simulation of cartilage: A systematic review

MRI-based mathematical and computational modeling studies can contribute to a better understanding of the mechanisms governing cartilage's mechanical performance and cartilage disease. In addition, distinct modeling of cartilage is needed to optimize artificial cartilage production. These studies have opened up the prospect of further deepening our understanding of cartilage function. Furthermore, these studies reveal the initiation of an engineering-level approach to how cartilage disease affects material properties and cartilage function. Aimed at researchers in the field of MRI-based cartilage simulation, research articles pertinent to MRI-based cartilage modeling were identified, reviewed, and summarized systematically. Various MRI applications for cartilage modeling are highlighted, and the limitations of different constitutive models used are addressed. In addition, the clinical application of simulations and studied diseases are discussed. The paper's quality, based on the developed questionnaire, was assessed, and out of 79 reviewed papers, 34 papers were determined as high-quality. Due to the lack of the best constitutive models for various clinical conditions, researchers may consider the effect of constitutive material models on the cartilage disease simulation. In the future, research groups may incorporate various aspects of machine learning into constitutive models and MRI data extraction to further refine the study methodology. Moreover, researchers should strive for further reproducibility and rigorous model validation and verification, such as gait analysis.

Hamstrings Contraction Regulates the Magnitude and Timing of the Peak ACL Loading During the Drop Vertical Jump in Female Athletes

Background

Anterior cruciate ligament (ACL) injury reduction training has focused on lower body strengthening and landing stabilization. In vitro studies have shown that quadriceps forces increase ACL strain, and hamstring forces decrease ACL strain. However, the magnitude of the effect of the quadriceps and hamstrings forces on ACL loading and its timing during in vivo landings remains unclear.

Purpose

To investigate the effect and timing of knee muscle forces on ACL loading during landing.

Study Design

Descriptive laboratory study.

Methods

A total of 13 young female athletes performed drop vertical jump trials, and their movements were recorded with 3-dimensional motion capture. Lower limb joint motion and muscle forces were estimated with OpenSim and applied to a musculoskeletal finite element (FE) model to estimate ACL loading during landings. The FE simulations were performed with 5 different conditions that included/excluded kinematics, ground-reaction force (GRF), and muscle forces.

Results

Simulation of landing kinematics without GRF or muscle forces yielded an estimated median ACL strain and force of 5.1% and 282.6 N. Addition of GRF to kinematic simulations increased ACL strain and force to 6.8% and 418.4 N (P < .05). Addition of quadriceps force to kinematics + GRF simulations nonsignificantly increased ACL strain and force to 7.2% and 478.5 N. Addition of hamstrings force to kinematics + GRF simulations decreased ACL strain and force to 2.6% and 171.4 N (P < .001). Addition of all muscles to kinematics + GRF simulations decreased ACL strain and force to 3.3% and 195.1 N (P < .001). With hamstrings force, ACL loading decreased from initial contact (time of peak: 1-18 milliseconds) while ACL loading without hamstrings force peaked at 47 to 98 milliseconds after initial contact (P = .024-.001). The knee flexion angle increased from 20.9° to 73.1° within 100 milliseconds after initial contact.

Conclusion

Hamstrings activation had greater effect relative to GRF and quadriceps activation on ACL loading, which significantly decreased and regulated the magnitude and timing of ACL loading during in vivo landings.

Clinical Relevance

Clinical training should focus on strategies that influence increased hamstrings activation during landing to reduce ACL loads.

Knee extension moment arm variations relate to mechanical function in walking and running

The patellofemoral joint plays a crucial mechanical role during walking and running. It increases the knee extensor mechanism's moment arm and reduces the knee extension muscle forces required to generate the extension moment that supports body weight, prevents knee buckling and propels the centre of mass. However, the mechanical implications of moment arm variation caused by patellofemoral and tibiofemoral motion remain unclear. We used a data-driven musculoskeletal model with a 12-degree-of-freedom knee to simulate the knee extension moment arm during walking and running. Using a geometric method to calculate the moment arm, we found smaller moment arms during running than during walking in the swing phase. Overall, knee flexion causes differences between running and walking moment arms as increased flexion causes a posterior shift in the tibiofemoral rotation axis and patella articulation with the distal femur. Moment arms were also affected by knee motion direction and best predicted by separating by direction instead of across the entire gait cycle. Furthermore, we found high inter-subject variation in the moment arm that was largely explained by out-of-plane motion. Our results are consistent with the concept that shorter moment arms increase the effective mechanical advantage of the knee and may contribute to increased running velocity.

Patella Apex Influences Patellar Ligament Forces and Ratio

The relationship between three-dimensional shape and patellofemoral mechanics is complicated. The Wiberg patella classification is a method of distinguishing shape differences in the axial plane of the patella that can be used to connect shape differences to observed mechanics. This study uses the Wiberg patella classification to differentiate variance in a statistical shape model describing changes in patella morphology and height. We investigate how patella morphology influences force distribution within the patellofemoral joint. The Wiberg type I patella has a more symmetrical medial and lateral facet while the type III patella has a larger lateral facet compared to medial. The second principal component of the statistical shape model described shape variation that qualitatively resembled the different Wiberg patellas. We generated patellofemoral morphologies from the statistical shape model and integrated them into a musculoskeletal model with a twelve degrees-of-freedom knee. We simulated an overground walking trial with these morphologies and recorded patellofemoral mechanics and ligament forces. An increase in patellar ligament force corresponded with an increase in patella height. Wiberg type III patellas had a sharper patella apex which related to lower ratios of quadriceps tendon forces to patellar ligament forces. The change in pivot point of the patella affects the ratio of forces as well as the patellofemoral reaction force. This study provides a better understanding of how patella morphology affects fundamental patella mechanics which may help identify at-risk populations for pathology development.

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

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

A Methodology to Evaluate the Effects of Kinematic Measurement Uncertainties on Knee Ligament Properties Estimated From Laxity Measurements

Ligaments are important joint stabilizers but assessing their mechanical properties remain challenging. We developed a methodology to investigate the effects of kinematic measurement uncertainty during laxity tests on optimization-based estimation of ligament properties. We applied this methodology to a subject-specific knee model with known ligament properties as inputs and compared the estimated to the known knee ligament properties under the influence of noise. Four different sets of laxity tests were simulated with an increasing number of load cases, capturing anterior/posterior, varus/valgus, and internal/external rotation loads at 0 deg and 30 deg of knee flexion. 20 samples of uniform random noise ([-0.5,0.5] mm and degrees) were added to each set and fed into an optimization routine that subsequently estimated the ligament properties based on the noise targets. We found a large range of estimated ligament properties (stiffness ranges of 5.97 kN, 7.64 kN, 8.72 kN, and 3.86 kN; reference strain ranges of 3.11%, 2.53%, 1.88%, and 1.58% for anterior cruciate ligament (ACL), posterior cruciate ligament (PCL), medical collateral ligament (MCL), and lateral collateral ligament (LCL), respectively) for three sets of laxity tests, including up to 22 load cases. A set of laxity tests with 60 load cases kept the stiffness and reference strain ranges below 470 N per unit strain and 0.85%, respectively. These results illustrate that kinematic measurement noise have a large impact on estimated ligament properties and we recommend that future studies assess and report both the estimated ligament properties and the associated uncertainties due to kinematic measurement noise.

Moment arm of the knee-extensor mechanism measured in vivo across a range of daily activities

We measured the moment arm of the knee-extensor mechanism as ten healthy young individuals performed six functional activities: level walking, downhill walking, stair ascent, stair descent, open-chain (non-weight-bearing) knee flexion, and open-chain knee extension. The moment arm of the knee-extensor mechanism was described by the moment arm of the patellar-tendon force, which acts to rotate the tibia about the instantaneous axis of rotation (screw axis) of the knee. A mobile biplane X-ray imaging system enabled simultaneous measurements of the three-dimensional movements of the femur, tibia and patella during each activity, from which the position and orientation of the screw axis and the patellar-tendon moment arm (PTMA) were determined. Mean PTMA across all activities and all participants remained nearly constant (~46 mm) from 0° to 70° of knee flexion and decreased by no more than 20% at higher flexion angles. The magnitude of the PTMA varied more substantially across individuals than across activities, indicating that the moment arm is more heavily influenced by differences in knee-joint geometry than muscle loading. Hence, PTMA measurements obtained for a given activity performed by one individual may be used with good confidence to describe the PTMA for any other activity performed by the same individual. Caution is advised when using PTMA measurements obtained from one individual to describe the moment arm in another individual even once the data are normalized by knee bone size, as the PTMA varied by as much as 13% from the mean across individuals.

Development and evaluation of a new methodology for Soft Tissue Artifact compensation in the lower limb

Skin Marker (SM) based motion capture is the most widespread technique used for motion analysis. Yet, the accuracy is often hindered by Soft Tissue Artifact (STA). This is a major issue in clinical gait analysis where kinematic results are used for decision-making. It also has a considerable influence on the results of rigid body and Finite Element (FE) musculoskeletal models that rely on SM-based kinematics to estimate muscle, contact and ligament forces. Current techniques devised to compensate for STA, in particular multi-body optimization methods, often consider simplified joint models. Although joint personalization with anatomical constraints has improved kinematic estimation, these models yet don't represent a fully reliable solution to the STA problem, thus allowing us to envisage an alternative approach. In this perspective, we propose to develop a conceptual FE-based model of the lower limb for STA compensation and evaluate it for 66 healthy subjects under level walking motor task. Both hip and knee joint kinematics were analyzed, considering both rotational and translational joint motion. Results showed that STA caused underestimation of the hip joint kinematics (up to 2.2°) for all rotational DoF, and overestimation of knee joint kinematics (up to 12°) except in flexion/extension. Joint kinematics, in particular the knee joint, appeared to be sensitive to soft tissue stiffness parameters (rotational and translational mean difference up to 1.5° and 3.4 mm). Analysis of the results using alternative joint representations highlighted the versatility of the proposed modeling approach. This work paves the way for using personalized models to compensate for STA in healthy subjects and different activities.

Anterior Cruciate Ligament Graft Tunnel Placement and Graft Angle Are Primary Determinants of Internal Knee Mechanics After Reconstructive Surgery

BACKGROUND

Graft placement is a modifiable and often discussed surgical factor in anterior cruciate ligament (ACL) reconstruction (ACLR). However, the sensitivity of functional knee mechanics to variability in graft placement is not well understood.

PURPOSE

To (1) investigate the relationship of ACL graft tunnel location and graft angle with tibiofemoral kinematics in patients with ACLR, (2) compare experimentally measured relationships with those observed with a computational model to assess the predictive capabilities of the model, and (3) use the computational model to determine the effect of varying ACL graft tunnel placement on tibiofemoral joint mechanics during walking.

STUDY DESIGN

Controlled laboratory study.

METHODS

Eighteen participants who had undergone ACLR were tested. Bilateral ACL footprint location and graft angle were assessed using magnetic resonance imaging (MRI). Bilateral knee laxity was assessed at the completion of rehabilitation. Dynamic MRI was used to measure tibiofemoral kinematics and cartilage contact during active knee flexion-extension. Additionally, a total of 500 virtual ACLR models were created from a nominal computational knee model by varying ACL footprint locations, graft stiffness, and initial tension. Laxity tests, active knee extension, and walking were simulated with each virtual ACLR model. Linear regressions were performed between internal knee mechanics and ACL graft tunnel locations and angles for the patients with ACLR and the virtual ACLR models.

RESULTS

Static and dynamic MRI revealed that a more vertical graft in the sagittal plane was significantly related (P < .05) to a greater laxity compliance index (R² = 0.40) and greater anterior tibial translation and internal tibial rotation during active knee extension (R² = 0.22 and 0.23, respectively). Similarly, knee extension simulations with the virtual ACLR models revealed that a more vertical graft led to greater laxity compliance index, anterior translation, and internal rotation (R² = 0.56, 0.26, and 0.13). These effects extended to simulations of walking, with a more vertical ACL graft inducing greater anterior tibial translation, ACL loading, and posterior migration of contact on the tibial plateaus.

CONCLUSION

This study provides clinical evidence from patients who underwent ACLR and from complementary modeling that functional postoperative knee mechanics are sensitive to graft tunnel locations and graft angle. Of the factors studied, the sagittal angle of the ACL was particularly influential on knee mechanics.

CLINICAL RELEVANCE

Early-onset osteoarthritis from altered cartilage loading after ACLR is common. This study shows that postoperative cartilage loading is sensitive to graft angle. Therefore, variability in graft tunnel placement resulting in small deviations from the anatomic ACL angle might contribute to the elevated risk of osteoarthritis after ACLR.

12 Degrees of Freedom Muscle Force Driven Fibril-Reinforced Poroviscoelastic Finite Element Model of the Knee Joint

Accurate knowledge of the joint kinematics, kinetics, and soft tissue mechanical responses is essential in the evaluation of musculoskeletal (MS) disorders. Since in vivo measurement of these quantities requires invasive methods, musculoskeletal finite element (MSFE) models are widely used for simulations. There are, however, limitations in the current approaches. Sequentially linked MSFE models benefit from complex MS and FE models; however, MS model's outputs are independent of the FE model calculations. On the other hand, due to the computational burden, embedded (concurrent) MSFE models are limited to simple material models and cannot estimate detailed responses of the soft tissue. Thus, first we developed a MSFE model of the knee with a subject-specific MS model utilizing an embedded 12 degrees of freedom (DoFs) knee joint with elastic cartilages in which included both secondary kinematic and soft tissue deformations in the muscle force estimation (inverse dynamics). Then, a muscle-force-driven FE model with fibril-reinforced poroviscoelastic cartilages and fibril-reinforced poroelastic menisci was used in series to calculate detailed tissue mechanical responses (forward dynamics). Second, to demonstrate that our workflow improves the simulation results, outputs were compared to results from the same FE models which were driven by conventional MS models with a 1 DoF knee, with and without electromyography (EMG) assistance. The FE model driven by both the embedded and the EMG-assisted MS models estimated similar results and consistent with experiments from literature, compared to the results estimated by the FE model driven by the MS model with 1 DoF knee without EMG assistance.

ReadySim: A computational framework for building explicit finite element musculoskeletal simulations directly from motion laboratory data

Musculoskeletal modeling allows researchers insight into joint mechanics which might not otherwise be obtainable through in vivo or in vitro studies. Common musculoskeletal modeling techniques involve rigid body dynamics software which often employ simplified joint representations. These representations have proven useful but are limited in performing single-framework deformable analyzes in structures of interest. Musculoskeletal finite element (MSFE) analysis allows for representation of structures in sufficient detail to obtain accurate solutions of the internal stresses and strains including complex contact conditions and material representations. Studies which performed muscle force optimization directly in a finite element framework were often limited in complexity to minimize computational time. Recent advances in computational efficiency and control schemes for muscle force prediction have made these solutions more practical. Yet, the formulation of subject-specific simulations remains a challenging problem. The objectives of this work were to develop an open-source computational framework to build and run simulations which (a) scale the size of MSFE models and efficiently estimate (b) joint kinematics and (c) muscle forces from human motion data collected in a typical gait laboratory. A computational framework was built using MATLAB and Python to interface with model input and output files. The software uses laboratory marker data to scale model segment lengths and estimate joint kinematics. Concurrent muscle force and tissue strain estimations are performed based on the estimated kinematics and ground reaction forces. This software will improve the usability and consistency of single-framework MSFE simulations. Both software and template model are made freely available on SimTK.Novelty Statement Single framework musculoskeletal modeling directly in a finite element environment for muscle force estimation and tissue strain analysis. Open dissemination of unilateral musculoskeletal finite element model and software used in manuscript.

In vitro compressive properties of skeletal muscles and inverse finite element analysis: Comparison of human versus animals

Virtual finite element human body models have been widely used in biomedical engineering, traffic safety injury analysis, etc. Soft tissue modeling like skeletal muscle accounts for a large portion of a human body model establishment, and its modeling method is not enough explored. The present study aims to investigate the compressive properties of skeletal muscles due to different species, loading rates and fiber orientations, in order to obtain available parameters of specific material laws as references for building or improving the human body model concerning both modeling accuracy and computational cost. A series of compressive experiments of skeletal muscles were implemented for human gastrocnemius muscle, bovine and porcine hind leg muscle. To avoid long-time preservation effects, all experimental tests were carried out in 24 h after that the samples were harvested. Considering computational cost and generally used in the previous human body models, one-order hyperelastic Ogden model and three-term simplified viscoelastic quasi-linear viscoelastic (QLV) were selected for numerical analysis. Inverse finite element analysis was employed to obtain corresponding material parameters. With good fitting records, the simulation results presented available material parameters for human body model establishment, and also indicated significant differences of muscle compressive properties due to species, loading rates and fiber orientations. When considering one-order Ogden law, it is worthy of noting that the inversed material parameters of the porcine muscles are similar to those of the human gastrocnemius regardless of fiber orientations. In conclusion, the obtained material parameters in the present study can be references for global human body and body segment modeling.

Development and Evaluation of a Subject-Specific Lower Limb Model With an Eleven-Degrees-of-Freedom Natural Knee Model Using Magnetic Resonance and Biplanar X-Ray Imaging During a Quasi-Static Lunge

Musculoskeletal (MS) models can be used to study the muscle, ligament, and joint mechanics of natural knees. However, models that both capture subject-specific geometry and contain a detailed joint model do not currently exist. This study aims to first develop magnetic resonance image (MRI)-based subject-specific models with a detailed natural knee joint capable of simultaneously estimating in vivo ligament, muscle, tibiofemoral (TF), and patellofemoral (PF) joint contact forces and secondary joint kinematics. Then, to evaluate the models, the predicted secondary joint kinematics were compared to in vivo joint kinematics extracted from biplanar X-ray images (acquired using slot scanning technology) during a quasi-static lunge. To construct the models, bone, ligament, and cartilage structures were segmented from MRI scans of four subjects. The models were then used to simulate lunges based on motion capture and force place data. Accurate estimates of TF secondary joint kinematics and PF translations were found: translations were predicted with a mean difference (MD) and standard error (SE) of 2.13 ± 0.22 mm between all trials and measures, while rotations had a MD ± SE of 8.57 ± 0.63 deg. Ligament and contact forces were also reported. The presented modeling workflow and the resulting knee joint model have potential to aid in the understanding of subject-specific biomechanics and simulating the effects of surgical treatment and/or external devices on functional knee mechanics on an individual level.

Functional movement assessment by means of inertial sensor technology to discriminate between movement behaviour of healthy controls and persons with knee osteoarthritis

BACKGROUND

Apart from biomechanical alterations in movement patterns, it is known that movement limitations in persons with knee osteoarthritis (PwKOA) are related to an individual's perception and belief regarding pain and disability. To gain more insights into the functional movement behaviour of PwKOA in a clinical setting, inertial sensor technology can be applied. This study first aims to evaluate the ability of inertial sensors to discriminate between healthy controls (HC) and PwKOA. Secondly, this study aims to determine the relationship between movement behaviour, pain-related factors and disability scores.

METHODS

Twelve HC and 19 PwKOA were included. Five repetitions of six functional movement tasks (walking, forward lunge, sideward lunge, ascent and descent stairs, single leg squat and sit-to-stand) were simultaneously recorded by the inertial sensor system and a camera-based motion analysis system. Statistically significant differences in angular waveforms of the trunk, pelvis and lower limb joints between HC and PwKOA were determined using one-dimensional statistical parametric mapping (SPM1D). The Knee injury and Osteoarthritis Outcome Score and TAMPA scale for Kinesiophobia were used to evaluate the relationship between discriminating joint motion, pain-related factors and disability using spearman's correlation coefficients.

RESULTS

PwKOA had significantly less trunk rotation, internal pelvis rotation and knee flexion ROM during walking. Additionally, the reduced knee flexion (i.e. at the end of the stance phase and swing phase) was related to increased level of perceived pain. During the sideward lunge, PwKOA had significantly less knee flexion, ankle plantarflexion and hip abduction. This decreased hip abduction (i.e. during stance) was related to higher fear of movement. Finally, PwKOA had significantly less knee flexion during the forward lunge, single leg squat and during ascent and descent stairs. No significant correlations were observed with disability.

CONCLUSIONS

Inertial sensors were able to discriminate between movement characteristics of PwKOA and HC. Additionally, significant relationships were found between joint motion, perceived pain and fear of movement. Since inertial sensors can be used outside the laboratory setting, these results are promising as they indicate the ability to evaluate movement deviations. Further research is required to enable measurements of small movement deviations in clinically relevant tasks.

Discriminant validity of 3D joint kinematics and centre of mass displacement measured by inertial sensor technology during the unipodal stance task

Background

The unipodal stance task is a clinical task that quantifies postural stability and alignment of the lower limb joints, while weight bearing on one leg. As persons with knee osteoarthritis (KOA) have poor postural and knee joint stability, objective assessment of this task might be useful.

Objective

To investigate the discriminant validity of three-dimensional joint kinematics and centre of mass displacement (COM) between healthy controls and persons with knee KOA, during unipodal stance using inertial sensors. Additionally, the reliability, agreement and construct validity are assessed to determine the reproducibility and accuracy of the discriminating parameters.

Methods

Twenty healthy controls and 19 persons with unilateral severe KOA were included. Five repetitions of the unipodal stance task were simultaneously recorded by an inertial sensor system and a camera-based system (gold standard). Statistical significant differences in kinematic waveforms between healthy controls and persons with severe knee KOA were determined using one-dimensional statistical parametric mapping (SPM1D).

Results

Persons with severe knee KOA had more lateral trunk lean towards the contralateral leg, more hip flexion throughout the performance of the unipodal stance task, more pelvic obliquity and COM displacement towards the contralateral side. However, for the latter two parameters the minimum detectable change was greater than the difference between healthy controls and persons with severe knee KOA. The construct validity was good (coefficient of multiple correlation 0.75, 0.83 respectively) and the root mean squared error (RMSE) was low (RMSE <1.5°) for the discriminant parameters.

Conclusion

Inertial sensor based movement analysis can discriminate between healthy controls and persons with severe knee KOA for lateral trunk lean and hip flexion, but unfortunately not for the knee angles. Further research is required to improve the reproducibility and accuracy of the inertial sensor measurements before they can be used to assess differences in tasks with a small range of motion.

Patellofemoral Mechanics: a Review of Pathomechanics and Research Approaches

PURPOSE OF REVIEW

The patellofemoral joint is a complicated articulation of the patella and femur that is prone to pathologies. The purpose of this review is to report on the current methods of investigating patellofemoral mechanics, factors that affect joint function, and future directions in patellofemoral joint research with emerging technologies and techniques.

RECENT FINDINGS

While previous hypotheses have suggested that the patella is only a moment arm extender, recent literature has suggested that the patella influences the control of knee moments and forces acting on the tibia as well as contributes to various aspects of patellar function with minimal neural input. With advancements in simulating a six-degrees-of-freedom patellofemoral joint, we have gained a better understanding of patella motion and have shown that geometry and muscle activations impact patella mechanics. Research into influences on patella mechanics from other joints such as the hip and foot has become more prevalent. In this review, we report current in vivo, in vitro, and in silico approaches to studying the patellofemoral joint. Kinematic and anatomical factors that affect patellofemoral joint function such as patella alta and tilt or bone morphology and ligaments are discussed. Moving forward, we suggest that advanced in vivo dynamic imaging methods coupled to musculoskeletal simulation will provide further understanding of patellofemoral pathomechanics and allow engineers and clinicians to design interventions to mitigate or prevent patellofemoral pathologies.