S Majid Nazemi1, Morteza Amini2, Saija A Kontulainen3, Jaques S Milner4, David W Holdsworth4, Bassam A Masri5, David R Wilson5, James D Johnston6. 1. Department of Mechanical Engineering, University of Saskatchewan, Saskatoon, Canada. Electronic address: majid.nazemi@gmail.com. 2. Department of Mechanical Engineering, University of Saskatchewan, Saskatoon, Canada. 3. College of Kinesiology, University of Saskatchewan, Saskatoon, Canada. 4. Robarts Research Institute, Western University, London, Canada. 5. Department of Orthopaedics and Centre for Hip Health and Mobility, University of British Columbia, Vancouver, BC, Canada. 6. Department of Mechanical Engineering, University of Saskatchewan, Saskatoon, Canada. Electronic address: jd.johnston@usask.ca.
Abstract
BACKGROUND: Quantitative computed tomography based subject-specific finite element modeling has potential to clarify the role of subchondral bone alterations in knee osteoarthritis initiation, progression, and pain initiation. Calculation of bone elastic moduli from image data is a basic step when constructing finite element models. However, different relationships between elastic moduli and imaged density (known as density-modulus relationships) have been reported in the literature. The objective of this study was to apply seven different trabecular-specific and two cortical-specific density-modulus relationships from the literature to finite element models of proximal tibia subchondral bone, and identify the relationship(s) that best predicted experimentally measured local subchondral structural stiffness with highest explained variance and least error. METHODS: Thirteen proximal tibial compartments were imaged via quantitative computed tomography. Imaged bone mineral density was converted to elastic moduli using published density-modulus relationships and mapped to corresponding finite element models. Proximal tibial structural stiffness values were compared to experimentally measured stiffness values from in-situ macro-indentation testing directly on the subchondral bone surface (47 indentation points). FINDINGS: Regression lines between experimentally measured and finite element calculated stiffness had R(2) values ranging from 0.56 to 0.77. Normalized root mean squared error varied from 16.6% to 337.6%. INTERPRETATION: Of the 21 evaluated density-modulus relationships in this study, Goulet combined with Snyder and Schneider or Rho appeared most appropriate for finite element modeling of local subchondral bone structural stiffness. Though, further studies are needed to optimize density-modulus relationships and improve finite element estimates of local subchondral bone structural stiffness.
BACKGROUND: Quantitative computed tomography based subject-specific finite element modeling has potential to clarify the role of subchondral bone alterations in knee osteoarthritis initiation, progression, and pain initiation. Calculation of bone elastic moduli from image data is a basic step when constructing finite element models. However, different relationships between elastic moduli and imaged density (known as density-modulus relationships) have been reported in the literature. The objective of this study was to apply seven different trabecular-specific and two cortical-specific density-modulus relationships from the literature to finite element models of proximal tibia subchondral bone, and identify the relationship(s) that best predicted experimentally measured local subchondral structural stiffness with highest explained variance and least error. METHODS: Thirteen proximal tibial compartments were imaged via quantitative computed tomography. Imaged bone mineral density was converted to elastic moduli using published density-modulus relationships and mapped to corresponding finite element models. Proximal tibial structural stiffness values were compared to experimentally measured stiffness values from in-situ macro-indentation testing directly on the subchondral bone surface (47 indentation points). FINDINGS: Regression lines between experimentally measured and finite element calculated stiffness had R(2) values ranging from 0.56 to 0.77. Normalized root mean squared error varied from 16.6% to 337.6%. INTERPRETATION: Of the 21 evaluated density-modulus relationships in this study, Goulet combined with Snyder and Schneider or Rho appeared most appropriate for finite element modeling of local subchondral bone structural stiffness. Though, further studies are needed to optimize density-modulus relationships and improve finite element estimates of local subchondral bone structural stiffness.
Authors: Amadeus C S de Alcântara; Israel Assis; Daniel Prada; Konrad Mehle; Stefan Schwan; Lucia Costa-Paiva; Munir S Skaf; Luiz C Wrobel; Paulo Sollero Journal: Materials (Basel) Date: 2019-12-24 Impact factor: 3.623
Authors: Alisdair R MacLeod; Nicholas Peckham; Gil Serrancolí; Ines Rombach; Patrick Hourigan; Vipul I Mandalia; Andrew D Toms; Benjamin J Fregly; Harinderjit S Gill Journal: Commun Med (Lond) Date: 2021-06-30
Authors: Hanieh Arjmand; Majid Nazemi; Saija A Kontulainen; Christine E McLennan; David J Hunter; David R Wilson; James D Johnston Journal: Sci Rep Date: 2018-07-31 Impact factor: 4.379