An evolutionary Wolff's law for trabecular architecture

A continuum model is proposed to describe the temporal evolution of both the density changes and the reorientation of the trabecular architecture given the applied stress state in the bone and certain material parameters of the bone. The data upon which the proposed model is to be based consist of experimentally determined remodeling rate coefficients and quantitative stereological and anisotropic elastic constant measurements of cancellous bone. The model shows that the system of differential equations governing the temporal changes in architecture is necessarily nonlinear. This nonlinearity is fundamental in that it stems from the fact that, during remodeling, the relationship between stress and strain is changing as the stress and strain variables themselves are changing. In order to preserve the remodeling property of the model, terms that are of the order strain times the changes in density and/or microstructural properties must be retained. If these terms were dropped, there would be no feedback mechanism for architectural adaptation and no adaptation of the trabecular architecture. There is, therefore, no linearized version of the model of the temporal evolution of trabecular architecture. An application of the model is illustrated by an example problem in which the temporal evolution of homogeneous trabecular architecture is predicted. A limitation of the proposed continuum model is the length scale below which it cannot be applied. The model cannot be applied in regions of cancellous bone where the trabecular bone architecture is relatively inhomogeneous or at a bone-implant interface.