Non-invasive regional parameter identification of degenerated human meniscus

Abstract

Accurate identification of local changes in the biomechanical properties of the normal and degenerative meniscus is critical to better understand knee joint osteoarthritis onset and progression. Ex-vivo material characterization is typically performed on specimens obtained from different locations, compromising the tissue's structural integrity and thus altering its mechanical behavior. Therefore, the aim of this in-silico study was to establish a non-invasive method to determine the region-specific material properties of the degenerated human meniscus. In a previous experimental magnetic resonance imaging (MRI) study, the spatial displacement of the meniscus and its root attachments in mildly degenerated (n = 12) and severely degenerated (n = 12) cadaveric knee joints was determined under controlled subject-specific axial joint loading. To simulate the experimental response of the lateral and medial menisci, individual finite element models were created utilizing a transverse isotropic hyper-poroelastic constitutive material formulation. The superficial displacements were applied to the individual models to calculate the femoral reaction force in an inverse finite element analysis. During particle swarm optimization, the four most sensitive material parameters were varied to minimize the error between the femoral reaction force and the force applied in the MRI loading experiment. Individual global and regional parameter sets were identified. In addition to in-depth model verification, prediction errors were determined to quantify the reliability of the identified parameter sets. Both compressibility of the solid meniscus matrix (+141 %, p ≤ 0.04) and hydraulic permeability (+53 %, p ≤ 0.04) were significantly increased in the menisci of severely degenerated knees compared to mildly degenerated knees, irrespective of the meniscus region. By contrast, tensile and shear properties were unaffected by progressive knee joint degeneration. Overall, the optimization procedure resulted in reliable and robust parameter sets, as evidenced by mean prediction errors of <1 %. In conclusion, the proposed approach demonstrated high potential for application in clinical practice, where it might provide a non-invasive diagnostic tool for the early detection of osteoarthritic changes within the knee joint.