Micromechanical modeling of the contact stiffness of an osseointegrated bone–implant interface

AbstractBackgroundThe surgical success of cementless implants is determined by the evolution of the biomechanical properties of the bone–implant interface (BII). One difficulty to model the biomechanical behavior of the BII comes from the implant surface roughness and from the partial contact between bone tissue and the implant. The determination of the constitutive law of the BII would be of interest in the context of implant finite element (FE) modeling to take into account the imperfect characteristics of the BII. The aim of the present study is to determine an effective contact stiffness$$\left( {K_{c}^{\text{FEM}} } \right)$$KcFEMof an osseointegrated BII accounting for its micromechanical features such as surface roughness, bone–implant contact ratio (BIC) and periprosthetic bone properties. To do so, a 2D FE model of the BII under normal contact conditions was developed and was used to determine the behavior of$$K_{c}^{\text{FEM}}$$KcFEM.ResultsThe model is validated by comparison with three analytical schemes based on micromechanical homogenization including two Lekesiz’s models (considering interacting and non-interacting micro-cracks) and a Kachanov’s model.$$K_{c}^{\text{FEM}}$$KcFEMis found to be comprised between 1013and 1015 N/m3according to the properties of the BII.$$K_{c}^{\text{FEM}}$$KcFEMis shown to increase nonlinearly as a function of the BIC and to decrease as a function of the roughness amplitude for high BIC values (above around 20%). Moreover,$$K_{c}^{\text{FEM}}$$KcFEMdecreases as a function of the roughness wavelength and increases linearly as a function of the Young’s modulus of periprosthetic bone tissue.ConclusionsThese results open new paths in implant biomechanical modeling since this model may be used in future macroscopic finite element models modeling the bone–implant system to replace perfectly rigid BII conditions.