Computational simulation of the multiphasic degeneration of the bone-cartilage unit during osteoarthritis via indentation and unconfined compression tests.

It has been experimentally proposed that the discrete regions of articular cartilage, along with different subchondral bone tissues, known as the bone-cartilage unit, are biomechanically altered during osteoarthritis degeneration. However, a computational framework capturing all of the dominant changes in the multiphasic parameters has not yet been developed. This article proposes a new finite element model of the bone-cartilage unit by combining several validated, nonlinear, depth-dependent, fibril-reinforced, and swelling models, which can computationally simulate the variations in the dominant parameters during osteoarthritis degeneration by indentation and unconfined compression tests. The mentioned dominant parameters include the proteoglycan depletion, collagen fibrillar softening, permeability, and fluid fraction increase for approximately non-advanced osteoarthritis. The results depict the importance of subchondral bone tissues in fluid distribution within the bone-cartilage units by decreasing the fluid permeation and pressure (up to a maximum of 100 kPa) during osteoarthritis, supporting the notion that subchondral bones might play a role in the pathogenesis of osteoarthritis. Furthermore, the osteoarthritis composition-based studies shed light on the significant biomechanical role of the calcified cartilage, which experienced a maximum change of 70 kPa in stress, together with relative load contributions of articular cartilage constituents during osteoarthritis, in which the osmotic pressure bore around 70% of the loads after degeneration. To conclude, the new insights provided by the results reveal the significance of the multiphasic osteoarthritis simulation and demonstrate the functionality of the proposed bone-cartilage unit model.