Development of robust finite element models of porcine tibiofemoral joints loaded under varied flexion angles and tibial freedoms

Effect of meniscus modelling assumptions in a static tibiofemoral finite element model: importance of geometry over material

Finite element studies of the tibiofemoral joint have increased use in research, with attention often placed on the material models. Few studies assess the effect of meniscus modelling assumptions in image-based models on contact mechanics outcomes. This work aimed to assess the effect of modelling assumptions of the meniscus on knee contact mechanics and meniscus kinematics. A sensitivity analysis was performed using three specimen-specific tibiofemoral models and one generic knee model. The assumptions in representing the meniscus attachment on the tibia (shape of the roots and position of the attachment), the material properties of the meniscus, the shape of the meniscus and the alignment of the joint were evaluated, creating 40 model instances. The values of material parameters for the meniscus and the position of the root attachment had a small influence on the total contact area but not on the meniscus displacement or the force balance between condyles. Using 3D shapes to represent the roots instead of springs had a large influence in meniscus displacement but not in knee contact area. Changes in meniscus shape and in knee alignment had a significantly larger influence on all outcomes of interest, with differences two to six times larger than those due to material properties. The sensitivity study demonstrated the importance of meniscus shape and knee alignment on meniscus kinematics and knee contact mechanics, both being more important than the material properties or the position of the roots. It also showed that differences between knees were large, suggesting that clinical interpretations of modelling studies using single geometries should be avoided.

A Finite Element Model to Investigate the Stability of Osteochondral Grafts Within a Human Tibiofemoral Joint

Osteochondral grafting has demonstrated positive outcomes for treating articular cartilage defects by replacing the damaged region with a cylindrical graft consisting of bone with a layer of cartilage. However, factors that cause graft subsidence are not well understood. The aim of this study was to develop finite element (FE) models of osteochondral grafts within a tibiofemoral joint, suitable for an investigation of parameters affecting graft stability. Cadaveric femurs were used to experimentally calibrate the bone properties and graft-bone frictional forces for use in corresponding image-based FE models, generated from µCT scan data. Effects of cartilage defects and osteochondral graft repair were measured by examining contact pressure changes using further in vitro tests. Here, six defects were created in the femoral condyles, which were subsequently treated with osteochondral autografts or metal pins. Matching image-based FE models were created, and the contact patches were compared. The bone material properties and graft-bone frictional forces were successfully calibrated from the initial tests with good resulting levels of agreement (CCC = 0.87). The tibiofemoral joint experiment provided a range of cases that were accurately described in the resultant pressure maps and were well represented in the FE models. Cartilage defects and repair quality were experimentally measurable with good agreement in the FE model pressure maps. Model confidence was built through extensive validation and sensitivity testing. It was found that specimen-specific properties were required to accurately represent graft behaviour. The final models produced are suitable for a range of parametric testing to investigate immediate graft stability.

Computational modelling of articular joints with biphasic cartilage: recent advances, challenges and opportunities

Biphasic models have been widely used to simulate the time-dependent biomechanical response of soft tissues. Modelling techniques of joints with biphasic weight-bearing soft tissues have been markedly improved over the last decade, enhancing our understanding of the function, degenerative mechanism and outcomes of interventions of joints. This paper reviews the recent advances, challenges and opportunities in computational models of joints with biphasic weight-bearing soft tissues. The review begins with an introduction of the function and degeneration of joints from a biomechanical aspect. Different constitutive models of articular cartilage, in particular biphasic materials, are illustrated in the context of the study of contact mechanics in joints. Approaches, advances and major findings of biphasic models of the hip and knee are presented, followed by a discussion of the challenges awaiting to be addressed, including the convergence issue, high computational cost and inadequate validation. Finally, opportunities and clinical insights in the areas of subject-specific modeling and tissue engineering are provided and discussed.

Development of robust finite element models to investigate the stability of osteochondral grafts within porcine femoral condyles

Osteoarthritis (OA) is the most prevalent chronic rheumatic disease worldwide with knee OA having an estimated lifetime risk of approximately 14%. Autologous osteochondral grafting has demonstrated positive outcomes in some patients, however, understanding of the biomechanical function and how treatments can be optimised remains limited. Increased short-term stability of the grafts allows cartilage surfaces to remain congruent prior to graft integration. In this study methods for generating specimen specific finite element (FE) models of osteochondral grafts were developed, using parallel experimental data for calibration and validation. Experimental testing of the force required to displace osteochondral grafts by 2 mm was conducted on three porcine knees, each with four grafts. Specimen specific FE models of the hosts and grafts were created from registered μCT scans captured from each knee (pre- and post-test). Material properties were based on the μCT background with a conversion between μCT voxel brightness and Young's modulus. This conversion was based on the results of the separate testing of eight porcine condyles and optimization of specimen specific FE models. The comparison between the experimental and computational push-in forces gave a strong agreement with a concordance correlation coefficient (CCC) = 0.75, validating the modelling approach. The modelling process showed that homogenous material properties based on whole bone BV/TV calculations are insufficient for accurate modelling and that an intricate description of the density distribution is required. The robust methodology can provide a method of testing different treatment options and can be used to investigate graft stability in full tibiofemoral joints.