Peak stress levels predicted in finite element analysis (FEA) usually depend on mesh density, due to singular points in the model. In an earlier study, an FEA algorithm was developed to simulate the damage accumulation process in the cement mantle around total hip replacement (THR) implants. It allows cement crack formation to be predicted, as a function of the local cement stress levels. As the simulation is driven by mesh-dependent peak stresses, predicted crack formation rates are also likely to be mesh dependent. The aim of this study was to evaluate the mesh dependence of the predicted crack formation process, and to present a method to reduce the mesh dependence. Crack-propagation experiments were simulated. Experimental specimens, re...