Mathematical development of a novel discrete hip deformation algorithm for the in silico elasto-hydrodynamic lubrication modelling of total hip replacements

Abstract

In this paper, the procedure to achieve an accurate deformation model of a total hip replacement (THR) was proposed with the aim to obtain a numerical tool to be simply merged into THR elasto-hydrodynamic computational synovial lubrication algorithms. The approach was based on the Finite Element Method (FEM) and was developed in a Matlab code, allowing the definition of the influence matrix and of a boundary conditions vector. It works with linear tetrahedra and performs the displacement calculation for both the acetabular cup and the femoral head, taking into account the anatomical hip relative motion, by coupling them with a cubic interpolation matrix. Two simulations were conducted in order to validate the algorithm and the results were compared with the ones obtained by the commercial software Ansys. The comparison provides a satisfactory agreement in terms of surface deformation, Von Mises stress and strain energy, proving the reliability of the model and the possibility to use the model in the in silico prostheses tribological simulations, avoiding the complexity and the high computational resource requirement coming from the coupling between complex lubrication algorithms and FEM commercial software, and with the possibility to directly act on many key parameter characteristics of the investigated problem.