An orthotropic continuum model with substructure evolution for describing bone remodeling: an interpretation of the primary mechanism behind Wolff’s law

We propose a variational approach that employs a generalized principle of virtual work to estimate both the mechanical response and the changes in living bone tissue during the remodeling process. This approach provides an explanation for the adaptive regulation of the bone substructure in the context of orthotropic material symmetry. We specifically focus upon the crucial gradual adjustment of bone tissue as a structural material that adapts its mechanical features, such as materials stiffnesses and microstructure, in response to the evolving loading conditions. We postulate that the evolution process relies on a feedback mechanism involving multiple stimulus signals. The mechanical and remodeling behavior of bone tissue is clearly a complex process that is difficult to describe within the framework of classical continuum theories. For this reason, a generalized continuum elastic theory is employed as a proper mathematical context for an adequate description of the examined phenomenon. To simplify the investigation, we considered a two-dimensional problem. Numerical simulations have been performed to illustrate bone evolution in a few significant cases: the bending of a rectangular cantilever plate and a three-point flexure test. The results are encouraging because they can replicate the optimization process observed in bone remodeling. The proposed model provides a likely distribution of stiffnesses and accurately represents the arrangement of trabeculae macroscopically described by the orthotropic symmetry directions, as supported by experimental evidence from the trajectorial theory.