An Analysis and Comparison of Convergence and Uniqueness of Time-Independent Bone Adaptation Models

Several stimuli are proposed in the bone remodeling theory. It is not clear, if a unique solution exists and if the result is convergent using a certain stimulus. In this study, the strain stimulus, strain energy stimulus and the von Mises stress stimulus for bone remodeling are compared and applied to a square plate model using the finite element method. In the plane stress state, the remodeling equilibrium equations are transformed into functions of only the principal strains and the graphs of these functions are drawn in a diagram using the principal strains as the variables of two coordinate axes. The equation of the sum of principal strain squared equal to a constant is a circle in the diagram. The remodeling equilibrium equation of the strain stimulus is a quadrangle fitting into the circle, the remodeling equilibrium equation of the strain energy stimulus is an ellipse and the remodeling equilibrium equation of the von Mises stress stimulus is also an ellipse close to the principal strains circle when we take the same constants in the above equations. Using the finite element method, two models are performed with the uniform initial elastic properties and with the semi-random initial distribution of the elastic properties. The principal strains as the final finite element results converge within 2% of the objective constant for all the different stimuli. The obtained Young's moduli of two models as the adaptation object are different but in equilibrium, i.e. the equilibrium solution of adaptation model is not unique. The principal strains can not be used to examine the uniqueness of solution, since two different solutions can have the same results of principal strains. Using a certain stimulus, certain initial properties and a certain iterative equation, the solution is unique in equilibrium. The results using the model in this study show also that the same results can be obtained using any of the three stimuli when a proper constant in each remodeling equilibrium equation is chosen.