Influence of posterior cruciate ligament tension on tibiofemoral and patellofemoral joint contact mechanics in cruciate-retaining total knee replacement: a combined musculoskeletal multibody and finite-element simulation

The influence of posterior cruciate ligament (PCL) tension on the clinical outcome of cruciate-retaining total knee replacement (CR-TKR) remains controversial. Various numerical approaches have been used to study this influence systematically, but the models used are limited by certain assumptions and simplifications. Therefore, the objective of this computational study was to develop a combined musculoskeletal multibody and finite-element simulation during a squat motion to 90° knee flexion with a CR-TKR design to overcome previous limitations regarding model inputs. In addition, different PCL tensions (tight, lax, resected) were modeled and the influence on tibiofemoral and resurfaced patellofemoral joint dynamics and contact stresses was evaluated. The effect of the PCL on knee joint dynamics and contact stresses was more pronounced at higher flexion angles. Tibiofemoral joint dynamics were influenced and a tight PCL induced increased posterior femoral translation during flexion. The maximum contact stress in the tibial insert increased from 20.6 MPa to 22.5 MPa for the resected and tightest PCL at 90° knee flexion. Patellofemoral joint dynamics were only slightly affected by PCL tension. However, the maximum contact stress in the patellar component decreased from 58.0 MPa to 53.7 MPa for the resected and tightest PCL at 90° knee flexion. The combination of musculoskeletal multibody and finite-element simulation is a sufficient method to comprehensively investigate knee joint dynamics and contact stresses in CR-TKR. The PCL tension after CR-TKR affects joint dynamics and contact stresses at the articulating implant surfaces.

A systematic comparison between FEBio and PolyFEM for biomechanical systems

BACKGROUND AND OBJECTIVES

Finite element simulations are widely employed as a non-invasive and cost-effective approach for predicting outcomes in biomechanical simulations. However, traditional finite element software, primarily designed for engineering materials, often encountered limitations in contact detection and enforcement, leading to simulation failure when dealing with complex biomechanical configurations. Currently, a lot of model tuning is required to get physically accurate finite element simulations without failures. This adds significant human interaction to each iteration of a biomechanical model. This study addressed these issues by introducing PolyFEM, a novel finite element solver that guarantees inversion- and intersection-free solutions with completely automatic collision detection. The objective of this research is to validate PolyFEM's capabilities by comparing its results with those obtained from a well-established finite element solver, FEBio.

METHODS

To achieve this goal, five comparison scenarios were formulated to assess and validate PolyFEM's performance. The simulations were reproduced using both PolyFEM and FEBio, and the final results were compared. The five comparison scenarios included: (1) reproducing simulations from the FEBio test suite, consisting of static, dynamic, and contact-driven simulations; (2) replicating simulations from the verification paper published alongside the original release of FEBio; (3) a biomechanically based contact problem; (4) creating a custom simulation involving high-energy collisions between soft materials to highlight the difference in collision methods between the two solvers; and (5) performing biomechanical simulations of biting and quasi-stance.

RESULTS

We found that PolyFEM was capable of replicating all simulations previously conducted in FEBio. Particularly noteworthy is PolyFEM's superiority in high-energy contact simulations, where FEBio fell short, unable to complete over half of the simulations in Scenario 4. Although some of the simulations required significantly more simulation time in PolyFEM compared to FEBio, it is important to highlight that PolyFEM achieved these results without the need for any additional model tuning or contact declaration.

DISCUSSION

Despite being in the early stages of development, PolyFEM currently provides verified solutions for hyperelastic materials that are consistent with FEBio, both in previously published workflows and novel finite element scenarios. PolyFEM exhibited the ability to tackle challenging biomechanical problems where other solvers fell short, thus offering the potential to enhance the accuracy and realism of future finite element analyses.

Endoscopic versus laparoscopic bariatric procedures: A computational biomechanical study through a patient-specific approach

BACKGROUND AND OBJECTIVES

Within the framework of computational biomechanics, finite element models of the gastric district could be seen as a potential clinical tool not only to study the effects apported by bariatric surgery, but also to compare different surgical techniques such as the new emerging Endoscopic Sleeve Gastroplasty (ESG) with respect to well-established ones (such as the Laparoscopic Sleeve Gastrectomy, LSG).

METHODS

This work realized a fully computational comparison between the outcomes obtained from 10 patient-specific stomach models, which were used to simulate ESG, and the complementary results obtained from models representing the post-LSG of the same subjects. Specifically, once the ESG was simulated, a mechanical stimulus was applied by increasing an intragastric pressure up to a maximum of 5 kPa, in order to replicate the process of food intake, as well as for post-LSG models.

RESULTS

Results revealed non negligible differences between the techniques also within the same subject. In particular, not only LSG could lead to a greater reduction in the stomach volume (about 77 % at baseline, which is strictly linked to weight loss), but also influence the gastric distension (12 % less than pre-operative models). On the contrary, if ESG would be performed, a more similar pre-operative mechanical stimulation of the gastric walls may be seen (difference of about 1 %), thus preserving the mechanosensation, but the detriment of the volume reduction (about 56 % at baseline, and even decreases with increasing pressure). Moreover, since results suggested ESG may be more influenced by the pre-operative gastric cavity than LSG, a predictive model was proposed to support the surgical planning and the estimation of the volume reduction after ESG.

CONCLUSIONS

ESG and LSG have substantial differences in their protocols and post-surgical effects. This work pointed out that variations between the two procedures may be observed also from a computational point of view, especially when including patient-specific geometries. These insights support gastric modelling as a valuable tool to evaluate, design and critically compare emerging bariatric surgical procedures, not only from empirical aspects and clinical outcomes, but also from a mechanical point of view.

Mechanical behavior of infrapatellar fat pad of patients affected by osteoarthritis

The infrapatellar fat pad (IFP) is an adipose tissue present in the knee that lies between the patella, femur, meniscus and tibia, filling the space between these structures. IFP facilitates the distribution of the synovial fluid and may act to absorb impulsive actions generated through the joint. IFP in osteoarthritis (OA) pathology undergoes structural changes characterized by inflammation, hypertrophy and fibrosis. The aim of the present study is to analyze the mechanical behavior of the IFP in patients affected by end-stage OA. A specific test fixture was designed and indentation tests were performed on IFP specimens harvested from OA patients who underwent total knee arthroplasty. Experiments allowed to assess the typical features of mechanical response, such as non-linear stress-strain behavior and time-dependent effects. Results from mechanical experimentations were implemented within the framework of a visco-hyperelastic constitutive theory, with the aim to provide data for computational modelling of OA IFP role in knee mechanics. Initial and final indentation stiffness were calculated for all subjects and statistical results reveled that OA IFP mechanics was not significantly influenced by gender, BMI and sample preparation. OA IFP mechanical behavior was also compared to that of other adipose tissues. OA IFP appeared to be a stiffer adipose tissue compared to subcutaneous, visceral adipose tissues and heel fat pads. It is reasonable that fibrosis induces a modification of the tissue destabilizing the normal distribution of forces in the joint during movement, causing a worsening of the disease.

Multiscale finite element musculoskeletal model for intact knee dynamics

BACKGROUND AND OBJECTIVE

The dynamic characteristics of the intact knee joint are valuable for treating knee osteoarthritis and designing knee prostheses. However, it remains a challenge to elucidate the detailed dynamics of the knee due to its complexity of anatomical structure and complex interaction with body dynamics.

METHODS

In this study, a unique subject-specific musculoskeletal model with a concurrent high-accuracy intact finite element knee model was created and used to simultaneously evaluate the kinematics and mechanics of an intact knee joint during the gait cycle.

RESULTS

A medial pivot motion with external rotation, and a large parallel anterior translation were observed in the stance and swing phases, respectively, which is consistent with the in vivo fluoroscopy measurements. The maximum axial contact force on the knee joint, observed at 45% of the gait cycle, is approximately 2.89 times the body weight. The medial cartilage bears 65.7% of the total axial contact force. The results demonstrate that the cartilage-cartilage contact bears most of the joint load (62.5%) compared to the cartilage-meniscus-cartilage contact (37.5%). Regarding contact mechanics, the maximum contact pressure on both sides of the tibial cartilage (8.2 MPa) is almost similar to the first axial loading peak (14%) of the gait cycle. Additionally, the maximum contact pressure (6.01 MPa) was observed during the stance phase of the gait cycle on the patellofemoral joint.

CONCLUSIONS

The predicted results on the tibiofemoral and patellofemoral joints provide a theoretical basis for the treatment of knee joint diseases and knee prosthesis design. Moreover, this approach presents a comprehensive tool to evaluate the mechanics at both the body and tissue levels. Therefore, it has a high potential for application in human biomechanics.

Using Musculoskeletal Models to Estimate in vivo Total Knee Replacement Kinematics and Loads: Effect of Differences Between Models

Total knee replacement (TKR) is one of the most performed orthopedic surgeries to treat knee joint diseases in the elderly population. Although the survivorship of knee implants may extend beyond two decades, the poor outcome rate remains considerable. A recent computational approach used to better understand failure modes and improve TKR outcomes is based on the combination of musculoskeletal (MSK) and finite element models. This combined multiscale modeling approach is a promising strategy in the field of computational biomechanics; however, some critical aspects need to be investigated. In particular, the identification and quantification of the uncertainties related to the boundary conditions used as inputs to the finite element model due to a different definition of the MSK model are crucial. Therefore, the aim of this study is to investigate this problem, which is relevant for the model credibility assessment process. Three different generic MSK models available in the OpenSim platform were used to simulate gait, based on the experimental data from the fifth edition of the "Grand Challenge Competitions to Predict in vivo Knee Loads." The outputs of the MSK analyses were compared in terms of relative kinematics of the knee implant components and joint reaction (JR) forces and moments acting on the tibial insert. Additionally, the estimated knee JRs were compared with those measured by the instrumented knee implant so that the "global goodness of fit" was quantified for each model. Our results indicated that the different kinematic definitions of the knee joint and the muscle model implemented in the different MSK models influenced both the motion and the load history of the artificial joint. This study demonstrates the importance of examining the influence of the model assumptions on the output results and represents the first step for future studies that will investigate how the uncertainties in the MSK models propagate on disease-specific finite element model results.

Comparison of Kinematics and Contact Mechanics in Normal Knee and Total Knee Replacements: A Computational Investigation

An objective of total knee replacement (TKR) is to restore the mechanical function of a normal knee. Joint kinematics and contact mechanics performance are two of the primary indices that indicate the success of TKR devices. The aim of this study was to compare the kinematics and contact mechanics of TKR and normal knee joints. An experimentally evaluated finite-element (FE) knee model was developed and used to investigate the performance of four TKR designs (fixed cruciate-retaining (CR), mobile CR, posterior-stabilized (PS), medial pivot design (MP)) and the normal knee joint during a gait cycle. The predicted kinematic results showed that the MP design presented similar kinematics to those of the normal knee joint and did not demonstrate paradoxical motion of the femur. A considerably larger contact area and lower contact pressure were found on the normal knee joint (1315 mm², and 14.8 MPa, respectively) than on the TKRs, which was consistent with the previous in-vivo fluoroscopic investigation. The mobile CR and PS designs exhibited the smallest and greatest contact pressures of the four TKR designs, respectively. The results of the present study help to understand the kinematics and contact mechanics in the TKR during the gait cycle, and provide comprehensive information about the performance of the normal knee joint.

In vivo kinematical validated knee model for preclinical testing of total knee replacement

BACKGROUND AND OBJECTIVE

A computational knee model facilitates efficient component design evaluations and preclinical testing under various dynamic loadings. However, the development of a highly mimicked dynamic whole knee model with specified ligament constraints that provides high predictive accuracy with in-vivo experiments remains a challenge.

METHODS

In the present study, a musculoskeletal integrated force-driven explicit finite-element knee model with tibiofemoral and patellofemoral joints constrained with detailed soft tissue was developed. A proportional-integral-derivative controller was concurrently added to the knee model to track the boundary conditions. The actuations of the quadriceps and hamstrings were predicted via a subject-specific musculoskeletal model and matched with electromyography results.

RESULTS

Compared to in-vivo fluoroscopic results in a gait cycle, the predicted results of the kinematics of the tibiofemoral joint exhibited an agreement in terms of tendency and magnitude (anterior-posterior translation: RMSE = 1.1 mm, r² = 0.87; inferior-superior translation: RMSE = 0.83 mm, r² = 0.84; medial-lateral translation: RMSE = 0.82 mm, r² = 0.05; flexion-extension rotation: RMSE = 0.23°, r² = 1; internal-external rotation: RMSE = 1.85°, r² = 0.65; varus-valgus rotation: RMSE = 1.39°, r² = 0.08). Contact mechanics, including the contact area, pressure, and stress, were synchronously simulated on the tibiofemoral and patellofemoral joints.

CONCLUSIONS

The study provides a calibrated knee model and a kinematical validation approach that can be widely used in preclinical testing and knee prosthesis design.

Enhanced In-Silico Polyethylene Wear Simulation of Total Knee Replacements During Daily Activities

A computational wear simulator is an efficient tool for evaluating the wear of artificial knee joints. The classical Archard's wear law-based simulator has questionable accuracy and is focused on walking. In this study, an in silico polyethylene wear simulation of total knee replacements was developed considering the various highly demanding daily activities. A good predicted accuracy (error = 8.1%) was found through comparison of the experimental results. A relatively larger averaged wear loss was found under the loading condition (1.53 mg/mc) of daily activities compared with the walking condition (1.32 mg/mc). The squatting movement (2.57 mg/mc) produces the highest overall wear rate. In addition, a relatively larger amount of wear was found on the medial side knee prosthesis than that on the lateral side. The enhanced in silico polyethylene wear simulator provides an accurate and comprehensive tool for wear prediction in preclinical wear testing.