Adaptive bone remodeling incorporating simultaneous density and anisotropy considerations

doi:10.1016/S0021-9290(96)00189-3

Journal of Biomechanics

Volume 30, Issue 6, June 1997, Pages 603-613

https://doi.org/10.1016/S0021-9290(96)00189-3 Get rights and content

Abstract

Over 100 years ago, Wolff hypothesized that cancellous bone altered both its apparent density and trabecular orientation in response to mechanical loads. A mathematical counterpart of this principle is derived by adding a remodeling rule for the rate-of-change of the full anisotropic stiffness tensor (all 21 independent terms) to the density rate-of-change rule adapted from an existing isotropic theory. As a result, anisotropy and density patterns develop such that the local stiffness tensor is optimal for the given series of applied loadings. The method does not rely on additional morphological measures of trabecular orientation. Furthermore, assumptions of material symmetry are not required, and any observed regions of orthotropy, transverse isotropy, or isotropy are a result entirely of the functional adaptation of the bone and not the consequence of a modeling assumption. This approach has been implemented with the finite element method and applied to a two-dimensional model of the proximal femur with encouraging results.

References (43)

D.R. Carter
Mechanical loading history and skeletal biology
J. Biomechanics
(1987)
D.R. Carter et al.
Trabecular bone density and loading history: regulation of connective tissue biology by mechanical energy
J. Biomechanics
(1987)
D.R. Carter et al.
Relationships between loading history and femoral cancellous bone architecture
J. Biomechanics
(1989)
S.A. Goldstein
The mechanical properties of trabecular bone: dependence on anatomic location and function
J. Biomechanics
(1987)
T.P. Harrigan et al.
An analytical and numerical study of the stability of bone remodelling theories: dependence on microstructural stimulus
J. Biomechanics
(1992)
T.P. Harrigan et al.
Limitations of the continuum assumption in cancellous bone
J. Biomechanics
(1988)
R. Huiskes et al.
Adaptive bone-remodeling theory applied to prosthetic-design analysis
J. Biomechanics
(1987)
M.E. Levenston et al.
The role of loading memory in bone adaptation simulations
Bone
(1994)
H. Roesler
The history of some fundamental concepts in bone biomechanics
J. Biomechanics
(1987)
N. Sasaki et al.
Orientation of bone mineral and its role in the anisotropic mechanical properties of bone—transverse anisotropy
J. Biomechanics
(1989)

J.C. Simo et al.

Strain and stress based continuum damage models—I. Formulation

Int. J. Solids Struct

(1987)

C.H. Turner

Homeostatic control of bone structure: an application of feedback theory

Bone

(1991)

C.H. Turner

Functional determinants of bone structure: beyond Wolff's law of bone transformation

Bone

(1992)

C.H. Turner et al.

Prediction of proximal femoral architecture using the isotropic strain criterion

J. Biomechanics

(1993)

H. Weinans et al.

The behavior of adaptive bone-remodeling simulation models

J. Biomechanics

(1992)

G.E. Bacon et al.

A neutron diffraction study of the bones of the foot

J. Anat.

(1984)

G.S. Beaupré et al.

An approach for time-dependent bone modeling and remodeling-theoretical development

J. Orthop. Res.

(1990)

G.S. Beaupré et al.

An approach for time-dependent bone modeling and remodeling-application: a preliminary remodeling simulation

J. Orthop. Res.

(1990)

J.M. Bourgery

Traité complet de l'anatomie de l'homme. I.

Osteologic

(1832)

P.G. Ciarlet

S.C. Cowin et al.

Bone remodeling I: a theory of adaptive elasticity

J. Elasticity

(1976)

Cited by (185)

Interstitial fluid transport in cortical bone porosities: Effects of blood pressure and mass exchange using porous media theory
2024, Mechanics of Materials
This research focuses on the importance of understanding bone adaptation for gaining insights into diseases such as osteoporosis. The study investigates the role of fluid flow in cortical bone adaptation, with a specific focus on the network of osteocyte canaliculi and the chemical reactions involved in bone resorption. Prior studies have implemented methods such as poroelasticity theory and dual porosity, but these models did not account for mass exchange between the solid and fluid phases within the bone or the pressure exerted by arteries. The movement of interstitial fluid between the vascular porosity and lacuno-canalicular porosity in cortical bone, as well as the velocity of the fluid, are affected by both the pressure exerted by arteries and the exchange of mass between the solid and fluid components of the bone. To incorporate these effects in the model, the authors propose a dual porosity model based on the theory of porous media. The model considers the transfer of mass between vascular and lacunar-canalicular porosities, as well as between solid and fluid phases, while also accounting for blood pressure variations. To verify the proposed model, a numerical example is solved and compared with the results of previous studies. The model enables the measurement and analysis of the interstitial fluid velocity and the pressure in the respective porosities. The results reveal that accounting for blood pressure pulses and mass exchange between solid and fluid phases can alter the velocity and pressure fields of interstitial fluid within cortical bone. In addition, this study investigates transient examples and various scenarios of the bone remodeling. Overall, this research provides a more accurate understanding of the bone remodeling and the effects of blood pressure variation and mass exchange between solid and fluid phases of cortical bone.
Approximated decompositions for computational continuum mechanics
2023, Journal of Computational Physics
Solid and fluid mechanics problems often involve the computation of eigenvalues and eigenvectors. One has short algebraic expressions for the Euler equations and the perfect gas equation of state. However, algebraic expressions for solid mechanics models can become substantial. Computing the eigenvalues and eigenvectors algebraically for different models can become cumbersome. Also, larger formulae may result in increased overall computational time. Internal energy plays a significant role in modeling the medium, and thermodynamic consistency is indispensable. The models should correctly capture heat transfer, mechanical response, and fluid flow, e.g., in mechanocaloric and bone remodeling (BR) problems. In mechanocaloric problems, temperature variation occurs because of load and unload stresses. In BR problems, several mathematical models were developed for modeling its behavior. Flexibility for testing different models is an essential requirement. Hence, our goal is to propose and test approximated eigensystem and singular value (SVD) decompositions for continuum mechanics problems. We tested the resolution and accuracy of our propositions in benchmarking CFD problems. For the solid mechanics, we solved a few classical problems, verified the accuracy, and solved mechanocaloric tests and bar under tension. Although more expensive, the approximated eigensystem decomposition is more flexible. It is high-order and high-resolution, can capture heat transfer, fluid flow, discontinuities, and mechanical response, and can model real-world applications. The SVD is a viable possibility and can be used to adjust the computational time of severe situations, but it needs to be improved.
Simple anisotropic model of Bone Adaptation - SAMBA
2022, Journal of the Mechanical Behavior of Biomedical Materials
Bone presents the ability to adapt itself to the evolving mechanical environment. A simple anisotropic model for bone adaptation allowing reproducing the evolution of the elastic properties and the reorientation of the anisotropy frame is proposed is this work. The elastic properties are related to the value of the bone apparent density. The evolution law of the density is described via two functions reflecting the activities of the osteoclast and osteoblast cells. The anisotropy of the elastic properties of the bone is assumed evolving continuously between those of trabecular and compact tissues. The existence of a target material frame is assumed to describe its reorientation. The rate of rotation or spin of the material frame is supposed to be proportional to the target orientation angle and to the activity of the osteoclasts and osteoblasts. The mechanical stimulus governing the evolution of the apparent density is defined as the ratio between the current strain energy density and its critical value corresponding to the initiation of bone damage. The simulations showed that this simple model can reproduce some essential phenomena observed during bone adaptation process.
Digital twins for understanding the mechanical adaptation of bone in disease and postsurgery
2022, Digital Human Modeling and Medicine: The Digital Twin
Bone is a dynamic tissue, constantly undergoing remodeling in response to mechanical and chemical signals. The concerted action of different bone cells and the role of mechanobiological pathways in bone remodeling has been studied extensively in the past few decades, but many aspects of this intricate physiological process still remain unknown. In silico modeling provides a powerful tool, complementing in vivo and in vitro experiments to gain new insights into mechano-adaptive remodeling of bone tissue in healthy, diseased and post-operative states. Therefore, this chapter will review the application of digital twins for studying mechanically-induced bone remodeling at the cellular and tissue scales.
Techniques for mitigating the checkerboard formation: application in bone remodeling simulations
2022, Medical Engineering and Physics
Citation Excerpt :
The phenomenological model proposed by Jacobs et al. [4] describes the behavior of bone tissue. We developed a 2D geometric model of the human femur (Fig. 1) similarly to presented in Jacobs et al. [4,36]. Mesh has 11,038 constant strain triangle (CST) elements (5,692 nodes).
The checkerboard phenomenon is a numerical problem that can occur in bone remodeling simulations. This phenomenon degrades the quality of the results and emphasizes a strong mesh dependency. The node-based approach is a well-used method for solving this problem. However, this method does not provide a stable solution, and long-time simulations do not reach steady-state. This paper aims to discuss techniques for mitigating the checkerboard formation and obtaining bone characterization. For this, we apply the simple nodal average technique on stress, density, or both simultaneously. Then, we compare the element-based and node-based approaches with and without smoothing. Also, we introduced a quantification parameter called standard density variation. This parameter measures the intensity of the checkerboard formation and analyzes the evolution of the bone tissue behavior. The application of the smoothing techniques mitigates the checkerboard formation. Also, the application allows obtaining a stable behavior for long-time simulations. The standard density variation allowed visualizing this stable behavior. In particular, the stress field smoothing technique provides the best results. Besides the stabilization and checkerboard control, the technique allows obtaining a physiological density distribution.
Bone remodelling in implanted proximal femur using topology optimization and parameterized cellular model
2022, Journal of the Mechanical Behavior of Biomedical Materials
The clinical relevance of bone remodelling predictions calls for accurate finite element (FE) modelling of implant-bone structure and musculoskeletal loading conditions. However, simplifications in muscle loading, material properties, has often been used in FE simulations. Bone adaptation induces changes in bone apparent density and its microstructure. Multiscale simulations, involving optimization methods and biomimetic microstructural models, have proven to be promising for predicting changes in bone morphology. The objective of the study is to develop a novel computational framework to predict bone remodelling around an uncemented femoral implant, using multiscale topology optimization and a parameterized cellular model. The efficacy of the scheme was evaluated by comparing the remodelling predictions with those of isotropic strain energy density (SED) and orthotropy based formulations. The characteristic functional groups and low-density regions of Ward's triangle, predicted by the optimization scheme, were comparable to micro-CT images of the proximal femur. Although the optimization scheme predicted well comparable material distribution in the 2D femur models, the obscured material orientations in some planes of the 3D model indicate the need for a more robust modelling of the boundary conditions. Regression analysis revealed a higher correlation (0.6472) between the topology optimization and SED models than the orthotropic predictions (0.4219). Despite higher bone apposition of 10–20% around the distal tip of the implant, the bone density distributions were well comparable to clinical observations towards the proximal femur. The proposed computational scheme appears to be a viable method for including bone anisotropy in the remodelling formulation.

View all citing articles on Scopus

View full text

ESB research award 1996Adaptive bone remodeling incorporating simultaneous density and anisotropy considerations

Abstract

J. Biomechanics

J. Biomechanics

J. Biomechanics

J. Biomechanics

J. Biomechanics

J. Biomechanics

J. Biomechanics

Bone

J. Biomechanics

J. Biomechanics

Int. J. Solids Struct

Bone

Bone

J. Biomechanics

J. Biomechanics

A neutron diffraction study of the bones of the foot

J. Anat.

An approach for time-dependent bone modeling and remodeling-theoretical development

J. Orthop. Res.

An approach for time-dependent bone modeling and remodeling-application: a preliminary remodeling simulation

J. Orthop. Res.

Traité complet de l'anatomie de l'homme. I.

Osteologic

Bone remodeling I: a theory of adaptive elasticity

J. Elasticity

ESB research award 1996
Adaptive bone remodeling incorporating simultaneous density and anisotropy considerations