Numerical modeling of long bone adaptation due to mechanical loading: Correlation with experiments

Natarajan Chennimalai Kumar, Jonathan A. Dantzig, Iwona M. Jasiuk, Alexander Robling, Charles H. Turner

Research output: Contribution to journalArticle

25 Citations (Scopus)

Abstract

The process of external bone adaptation in cortical bone is modeled mathematically using finite element (FE) stress analysis coupled with an evolution model, in which adaptation response is triggered by mechanical stimulus represented by strain energy density. The model is applied to experiments in which a rat ulna is subjected to cyclic loading, and the results demonstrate the ability of the model to predict the bone adaptation response. The FE mesh is generated from micro-computed tomography (μCT) images of the rat ulna, and the stress analysis is carried out using boundary and loading conditions on the rat ulna obtained from the experiments [Robling, A. G., F. M. Hinant, D. B. Burr, and C. H. Turner. J. Bone Miner. Res. 17:1545-1554, 2002]. The external adaptation process is implemented in the model by moving the surface nodes of the FE mesh based on an evolution law characterized by two parameters: one that captures the rate of the adaptation process (referred to as gain); and the other characterizing the threshold value of the mechanical stimulus required for adaptation (referred to as threshold-sensitivity). A parametric study is carried out to evaluate the effect of these two parameters on the adaptation response. We show, following comparison of results from the simulations to the experimental observations of Robling et al. (J. Bone Miner. Res. 17:1545-1554, 2002), that splitting the loading cycles into different number of bouts affects the threshold-sensitivity but not the rate of adaptation. We also show that the threshold-sensitivity parameter can quantify the mechanosensitivity of the osteocytes.

Original languageEnglish (US)
Pages (from-to)594-604
Number of pages11
JournalAnnals of Biomedical Engineering
Volume38
Issue number3
DOIs
StatePublished - Mar 2010

Fingerprint

Bone
Rats
Miners
Stress analysis
Experiments
Strain energy
Tomography

Keywords

  • Cortical bone adaptation
  • Evolution law
  • Finite element modeling
  • Loading bouts
  • Rat ulna

ASJC Scopus subject areas

  • Biomedical Engineering

Cite this

Numerical modeling of long bone adaptation due to mechanical loading : Correlation with experiments. / Chennimalai Kumar, Natarajan; Dantzig, Jonathan A.; Jasiuk, Iwona M.; Robling, Alexander; Turner, Charles H.

In: Annals of Biomedical Engineering, Vol. 38, No. 3, 03.2010, p. 594-604.

Research output: Contribution to journalArticle

Chennimalai Kumar, Natarajan ; Dantzig, Jonathan A. ; Jasiuk, Iwona M. ; Robling, Alexander ; Turner, Charles H. / Numerical modeling of long bone adaptation due to mechanical loading : Correlation with experiments. In: Annals of Biomedical Engineering. 2010 ; Vol. 38, No. 3. pp. 594-604.
@article{a4f7219cc4684486b9e9c7d2e6ecf766,
title = "Numerical modeling of long bone adaptation due to mechanical loading: Correlation with experiments",
abstract = "The process of external bone adaptation in cortical bone is modeled mathematically using finite element (FE) stress analysis coupled with an evolution model, in which adaptation response is triggered by mechanical stimulus represented by strain energy density. The model is applied to experiments in which a rat ulna is subjected to cyclic loading, and the results demonstrate the ability of the model to predict the bone adaptation response. The FE mesh is generated from micro-computed tomography (μCT) images of the rat ulna, and the stress analysis is carried out using boundary and loading conditions on the rat ulna obtained from the experiments [Robling, A. G., F. M. Hinant, D. B. Burr, and C. H. Turner. J. Bone Miner. Res. 17:1545-1554, 2002]. The external adaptation process is implemented in the model by moving the surface nodes of the FE mesh based on an evolution law characterized by two parameters: one that captures the rate of the adaptation process (referred to as gain); and the other characterizing the threshold value of the mechanical stimulus required for adaptation (referred to as threshold-sensitivity). A parametric study is carried out to evaluate the effect of these two parameters on the adaptation response. We show, following comparison of results from the simulations to the experimental observations of Robling et al. (J. Bone Miner. Res. 17:1545-1554, 2002), that splitting the loading cycles into different number of bouts affects the threshold-sensitivity but not the rate of adaptation. We also show that the threshold-sensitivity parameter can quantify the mechanosensitivity of the osteocytes.",
keywords = "Cortical bone adaptation, Evolution law, Finite element modeling, Loading bouts, Rat ulna",
author = "{Chennimalai Kumar}, Natarajan and Dantzig, {Jonathan A.} and Jasiuk, {Iwona M.} and Alexander Robling and Turner, {Charles H.}",
year = "2010",
month = "3",
doi = "10.1007/s10439-009-9861-4",
language = "English (US)",
volume = "38",
pages = "594--604",
journal = "Annals of Biomedical Engineering",
issn = "0090-6964",
publisher = "Springer Netherlands",
number = "3",

}

TY - JOUR

T1 - Numerical modeling of long bone adaptation due to mechanical loading

T2 - Correlation with experiments

AU - Chennimalai Kumar, Natarajan

AU - Dantzig, Jonathan A.

AU - Jasiuk, Iwona M.

AU - Robling, Alexander

AU - Turner, Charles H.

PY - 2010/3

Y1 - 2010/3

N2 - The process of external bone adaptation in cortical bone is modeled mathematically using finite element (FE) stress analysis coupled with an evolution model, in which adaptation response is triggered by mechanical stimulus represented by strain energy density. The model is applied to experiments in which a rat ulna is subjected to cyclic loading, and the results demonstrate the ability of the model to predict the bone adaptation response. The FE mesh is generated from micro-computed tomography (μCT) images of the rat ulna, and the stress analysis is carried out using boundary and loading conditions on the rat ulna obtained from the experiments [Robling, A. G., F. M. Hinant, D. B. Burr, and C. H. Turner. J. Bone Miner. Res. 17:1545-1554, 2002]. The external adaptation process is implemented in the model by moving the surface nodes of the FE mesh based on an evolution law characterized by two parameters: one that captures the rate of the adaptation process (referred to as gain); and the other characterizing the threshold value of the mechanical stimulus required for adaptation (referred to as threshold-sensitivity). A parametric study is carried out to evaluate the effect of these two parameters on the adaptation response. We show, following comparison of results from the simulations to the experimental observations of Robling et al. (J. Bone Miner. Res. 17:1545-1554, 2002), that splitting the loading cycles into different number of bouts affects the threshold-sensitivity but not the rate of adaptation. We also show that the threshold-sensitivity parameter can quantify the mechanosensitivity of the osteocytes.

AB - The process of external bone adaptation in cortical bone is modeled mathematically using finite element (FE) stress analysis coupled with an evolution model, in which adaptation response is triggered by mechanical stimulus represented by strain energy density. The model is applied to experiments in which a rat ulna is subjected to cyclic loading, and the results demonstrate the ability of the model to predict the bone adaptation response. The FE mesh is generated from micro-computed tomography (μCT) images of the rat ulna, and the stress analysis is carried out using boundary and loading conditions on the rat ulna obtained from the experiments [Robling, A. G., F. M. Hinant, D. B. Burr, and C. H. Turner. J. Bone Miner. Res. 17:1545-1554, 2002]. The external adaptation process is implemented in the model by moving the surface nodes of the FE mesh based on an evolution law characterized by two parameters: one that captures the rate of the adaptation process (referred to as gain); and the other characterizing the threshold value of the mechanical stimulus required for adaptation (referred to as threshold-sensitivity). A parametric study is carried out to evaluate the effect of these two parameters on the adaptation response. We show, following comparison of results from the simulations to the experimental observations of Robling et al. (J. Bone Miner. Res. 17:1545-1554, 2002), that splitting the loading cycles into different number of bouts affects the threshold-sensitivity but not the rate of adaptation. We also show that the threshold-sensitivity parameter can quantify the mechanosensitivity of the osteocytes.

KW - Cortical bone adaptation

KW - Evolution law

KW - Finite element modeling

KW - Loading bouts

KW - Rat ulna

UR - http://www.scopus.com/inward/record.url?scp=77952010862&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=77952010862&partnerID=8YFLogxK

U2 - 10.1007/s10439-009-9861-4

DO - 10.1007/s10439-009-9861-4

M3 - Article

C2 - 20013156

AN - SCOPUS:77952010862

VL - 38

SP - 594

EP - 604

JO - Annals of Biomedical Engineering

JF - Annals of Biomedical Engineering

SN - 0090-6964

IS - 3

ER -