Numerical analysis of a poroelastic cartilage model: Investigating the influence of changing material properties in osteoarthritis

Abstract

Several changes occur in both the cartilage's material properties and anatomical structure as osteoarthritis progresses. Unlike most numerical studies that solely consider individual changes, our study aimed to understand the impact on cartilage mechanics by considering the combined effect of material properties and cartilage thickness varied with osteoarthritis progression. In total, 3 three-dimensional finite element models, representing the intact, early, and late osteoarthritis conditions, were developed to simulate a load-bearing area in the knee. The articular cartilage was modelled as fluid-saturated linear biphasic poroelastic to incorporate solid-fluid interaction. All models underwent prolonged creep (50 N) and relaxation (0.3 mm) analyses for 600 s. In the early stage of osteoarthritis, the tibial cartilage demonstrated an overall stiffer behaviour attributed to cartilage swelling despite decreased stiffness at the material level. On the other hand, in the late stage of osteoarthritis, the decrease in cartilage thickness led to increased knee deformation. Additionally, increased permeability resulted in accelerated fluid exudation across all osteoarthritis models, and the elevation in void ratio further intensified fluid pressure within the cartilage to a higher magnitude. Furthermore, these changes collectively influenced both the magnitude and distribution of the outcomes. A holistic understanding of the material properties altered in osteoarthritis may contribute to a better understanding of the mechanical performance of cartilage during disease progression.