The behavior of adaptive bone-remodeling simulation models

Abstract

The process of adaptive bone remodeling can be described mathematically and simulated in a computer model, integrated with the finite element method. In the model discussed here, cortical and trabecular bone are described as continuous materials with variable density. The remodeling rule applied to simulate the remodeling process in each element individually is, in fact, an objective function for an optimization process, relative to the external load. Its purpose is to obtain a constant, preset value for the strain energy per unit bone mass, by adapting the density. If an element in the structure cannot achieve that, it either turns to its maximal density (cortical bone) or resorbs completely. It is found that the solution obtained in generally a discontinuous patchwork. For a two-dimensional proximal femur model this patchwork shows a good resemblance with the density distribution of a real proximal femur. It is shown that the discontinuous end configuration is dictated by the nature of the differential equations describing the remodeling process. This process can be considered as a nonlinear dynamical system with many degrees of freedom, which behaves divergent relative to the objective, leading to many possible solutions. The precise solution is dependent on the parameters in the remodeling rule, the load and the initial conditions. The feedback mechanism in the process is self-enhancing, denser bone attracts more strain energy, whereby the bone becomes even more dense. It is suggested that this positive feedback of the attractor state (the strain energy field) creates order in the end configuration. In addition, the process ensures that the discontinuous end configuration is a structure with a relatively low mass, perhaps a minimal-mass structure, although this is no explicit objective in the optimization process. It is hypothesized that trabecular bone is a chaotically ordered structure which can be considered as a fractal with characteristics of optimal mechanical resistance and minimal mass, of which the actual morphology depends on the local (internal) loading characteristics, the sensor-cell density and the degree of mineralization.

Interaction of bone–dental implant with new ultra low modulus alloy using a numerical approach

“…Obtaining post-operative osseointegration is necessary in order to establish a solid and durable connection between the implant and the osseous structure. In agreement with the Wolff law, a process of osseous remodelling adapted to the stress level occurs after implantation [2][3][4]. This process is controlled by mechanical loads.…”

Untitled

“…How far back the phylogenetic tree should an ancestor be so that its descendent can be considered optimized with respect to it? (Computational analyses predicting optimization of, say, the femur suffer from an analogous difficulty because the assumed initial pattern largely determines the final outcome [60,61].) On the subject of optimization, it is perhaps salutary to relate the story of a former Senior Fellow of Trinity College Dublin, a one Rev.…”

Computational Modelling of Mechanical Behavior of Biological Tissues for Biomedical Applications

“…In this study, a cubic relationship between the Young's modulus and bone density was in-troduced for the bone-remodelling region according to the theory of Carter and Hayes [24]. Furthermore, based on the remodelling theory of Weinans et al [25], the density change in the bone tissue was in accordance with the strain-energy density. An ABAQUS user subroutine USDFLD was developed to implement the remodelling algorithm.…”