Computational modelling of bone cement polymerization: Temperature and residual stresses

The two major concerns associated with the use of bone cement are the generation of residual stresses and possible thermal necrosis of surrounding bone. An accurate modelling of these two factors could be a helpful tool to improve cemented hip designs. Therefore, a computational methodology based on previous published works is presented in this paper combining a kinetic and an energy balance equation. New assumptions are that both the elasticity modulus and the thermal expansion coefficient depend on the bone cement polymerization fraction. This model allows to estimate the thermal distribution in the cement which is later used to predict the stress-locking effect, and to also estimate the cement residual stresses. In order to validate the model, computational results are compared with experiments performed on an idealized cemented femoral implant. It will be shown that the use of the standard finite element approach cannot predict the exact temporal evolution of the temperature nor the residual stresses, underestimating and overestimating their value, respectively. However, this standard approach can estimate the peak and long-term values of temperature and residual stresses within acceptable limits of measured values. Therefore, this approach is adequate to evaluate residual stresses for the mechanical design of cemented implants. In conclusion, new numerical techniques should be proposed in order to achieve accurate simulations of the problem involved in cemented hip replacements.