Finite element simulation of adaptive bone remodelling

A stability criterion and a time stepping method

Timothy Harrigan, James J. Hamilton

Research output: Contribution to journalArticle

Abstract

Adaptive bone remodelling simulations which use finite element analysis can potentially aid in the design of orthopedic implants and can provide examples which test specific bone remodelling hypotheses in a quantitative manner. By concentrating on remodelling algorithms in which the geometry is fixed but the tissue stiffness changes based on strain energy density, we have predicted stability conditions for bone remodelling and we have tested the applicability of these conditions using numerical simulations. The stability requirements arrived at using a finite element formulation are similar to the requirements arrived at in an earlier analytical study. In order to test the stability conditions, we have developed an Euler backward time stepping technique which uses the derivation for stability. These simulations arrived at solutions which were impossible using Euler forward time stepping as applied in this study. Cases in which a simplified version of the derived Euler backward method are unstable or marginally stable have also been seen, but when the Euler backward method is applied using the full derived matrices, no instabilities are apparent. The results of the stability tests indicate that the converged density distributions in the examples studied are stable. Although a priori conditions which ensure stability are not found, a test for stability is provided, given an assumed density distribution.

Original languageEnglish (US)
Pages (from-to)837-854
Number of pages18
JournalInternational Journal for Numerical Methods in Engineering
Volume36
Issue number5
DOIs
StatePublished - Jan 1 1993
Externally publishedYes

Fingerprint

Bone Remodeling
Stability criteria
Time Stepping
Finite Element Simulation
Stability Criteria
Bone
Backward Euler Method
Stability Condition
Euler
Finite Element
Strain Energy Density
Remodeling
Implant
Requirements
Stiffness
Simulation
Unstable
Orthopedics
Strain energy
Numerical Simulation

ASJC Scopus subject areas

  • Numerical Analysis
  • Engineering(all)
  • Applied Mathematics

Cite this

Finite element simulation of adaptive bone remodelling : A stability criterion and a time stepping method. / Harrigan, Timothy; Hamilton, James J.

In: International Journal for Numerical Methods in Engineering, Vol. 36, No. 5, 01.01.1993, p. 837-854.

Research output: Contribution to journalArticle

@article{6cc162c638d94d5ab724d27d7f33dfed,
title = "Finite element simulation of adaptive bone remodelling: A stability criterion and a time stepping method",
abstract = "Adaptive bone remodelling simulations which use finite element analysis can potentially aid in the design of orthopedic implants and can provide examples which test specific bone remodelling hypotheses in a quantitative manner. By concentrating on remodelling algorithms in which the geometry is fixed but the tissue stiffness changes based on strain energy density, we have predicted stability conditions for bone remodelling and we have tested the applicability of these conditions using numerical simulations. The stability requirements arrived at using a finite element formulation are similar to the requirements arrived at in an earlier analytical study. In order to test the stability conditions, we have developed an Euler backward time stepping technique which uses the derivation for stability. These simulations arrived at solutions which were impossible using Euler forward time stepping as applied in this study. Cases in which a simplified version of the derived Euler backward method are unstable or marginally stable have also been seen, but when the Euler backward method is applied using the full derived matrices, no instabilities are apparent. The results of the stability tests indicate that the converged density distributions in the examples studied are stable. Although a priori conditions which ensure stability are not found, a test for stability is provided, given an assumed density distribution.",
author = "Timothy Harrigan and Hamilton, {James J.}",
year = "1993",
month = "1",
day = "1",
doi = "10.1002/nme.1620360508",
language = "English (US)",
volume = "36",
pages = "837--854",
journal = "International Journal for Numerical Methods in Engineering",
issn = "0029-5981",
publisher = "John Wiley and Sons Ltd",
number = "5",

}

TY - JOUR

T1 - Finite element simulation of adaptive bone remodelling

T2 - A stability criterion and a time stepping method

AU - Harrigan, Timothy

AU - Hamilton, James J.

PY - 1993/1/1

Y1 - 1993/1/1

N2 - Adaptive bone remodelling simulations which use finite element analysis can potentially aid in the design of orthopedic implants and can provide examples which test specific bone remodelling hypotheses in a quantitative manner. By concentrating on remodelling algorithms in which the geometry is fixed but the tissue stiffness changes based on strain energy density, we have predicted stability conditions for bone remodelling and we have tested the applicability of these conditions using numerical simulations. The stability requirements arrived at using a finite element formulation are similar to the requirements arrived at in an earlier analytical study. In order to test the stability conditions, we have developed an Euler backward time stepping technique which uses the derivation for stability. These simulations arrived at solutions which were impossible using Euler forward time stepping as applied in this study. Cases in which a simplified version of the derived Euler backward method are unstable or marginally stable have also been seen, but when the Euler backward method is applied using the full derived matrices, no instabilities are apparent. The results of the stability tests indicate that the converged density distributions in the examples studied are stable. Although a priori conditions which ensure stability are not found, a test for stability is provided, given an assumed density distribution.

AB - Adaptive bone remodelling simulations which use finite element analysis can potentially aid in the design of orthopedic implants and can provide examples which test specific bone remodelling hypotheses in a quantitative manner. By concentrating on remodelling algorithms in which the geometry is fixed but the tissue stiffness changes based on strain energy density, we have predicted stability conditions for bone remodelling and we have tested the applicability of these conditions using numerical simulations. The stability requirements arrived at using a finite element formulation are similar to the requirements arrived at in an earlier analytical study. In order to test the stability conditions, we have developed an Euler backward time stepping technique which uses the derivation for stability. These simulations arrived at solutions which were impossible using Euler forward time stepping as applied in this study. Cases in which a simplified version of the derived Euler backward method are unstable or marginally stable have also been seen, but when the Euler backward method is applied using the full derived matrices, no instabilities are apparent. The results of the stability tests indicate that the converged density distributions in the examples studied are stable. Although a priori conditions which ensure stability are not found, a test for stability is provided, given an assumed density distribution.

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

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

U2 - 10.1002/nme.1620360508

DO - 10.1002/nme.1620360508

M3 - Article

VL - 36

SP - 837

EP - 854

JO - International Journal for Numerical Methods in Engineering

JF - International Journal for Numerical Methods in Engineering

SN - 0029-5981

IS - 5

ER -