Numerical analysis of a strain-adaptive bone remodelling problem

In this paper we study from the numerical point of view a strain-adaptive bone remodelling that couples the displacements and the apparent density (the porosity) of the bone. The rate of this density at a particular location is described as an objective function, which depends on a particular stimulus at that location. The variational problem is written as a coupled system of a nonlinear variational equation for the displacement field and a nonlinear parabolic variational inequality for the apparent density. Then, fully discrete approximations are provided by using the finite element method to approximate the spatial variable and an Euler scheme to discretize the time derivatives. Error estimates are proved, from which, under adequate regularity conditions, the linear convergence of the algorithm is deduced. Finally, some numerical simulations involving one- and two-dimensional test examples are presented to demonstrate the accuracy of the approximation and the behaviour of the solution.

[1]  José R. Fernández,et al.  Analysis of a bone remodeling model , 2008 .

[2]  Isabel N. Figueiredo,et al.  Asymptotic Model of a Nonlinear Adaptive Elastic Rod , 2004 .

[3]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[4]  Isabel N. Figueiredo,et al.  Shape Analysis of an Adaptive Elastic Rod Model , 2005, SIAM J. Appl. Math..

[5]  L. Trabucho.  Non-linear bone remodelling: an existence and uniqueness result , 2000 .

[6]  Análisis de un método numérico con elementos finitos para problemas de contacto unilateral sin rozamiento en elasticidad. II parte: aproximación y resolución de los problemas discretos , 1986 .

[7]  Stephen C. Cowin,et al.  The Exact Stimulus for the Strain Adaptation of Bone Tissue is Unknown , 2006 .

[8]  H Weinans,et al.  A physiological approach to the simulation of bone remodeling as a self-organizational control process. , 1994, Journal of biomechanics.

[9]  H. Grootenboer,et al.  The behavior of adaptive bone-remodeling simulation models. , 1992, Journal of biomechanics.

[10]  Meir Shillor,et al.  Quasistatic evolution of damage in an elastic body: numerical analysis and computational experiments , 2007 .

[11]  S. Cowin,et al.  Bone remodeling I: theory of adaptive elasticity , 1976 .

[12]  Stephen C. Cowin,et al.  Bone remodeling II: small strain adaptive elasticity , 1976 .

[13]  W. Hayes,et al.  The compressive behavior of bone as a two-phase porous structure. , 1977, The Journal of bone and joint surgery. American volume.

[14]  H. M. Frost,et al.  Vital biomechanics: Proposed general concepts for skeletal adaptations to mechanical usage , 1988, Calcified Tissue International.

[15]  D. Carter,et al.  A unifying principle relating stress to trabecular bone morphology , 1986, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[16]  H Weinans,et al.  Adaptive bone remodeling and biomechanical design considerations for noncemented total hip arthroplasty. , 1989, Orthopedics.

[17]  W. Han,et al.  Quasistatic Contact Problems in Viscoelasticity and Viscoplasticity , 2002 .

[18]  H. Grootenboer,et al.  Adaptive bone-remodeling theory applied to prosthetic-design analysis. , 1987, Journal of biomechanics.

[19]  I. Figueiredo,et al.  Conical differentiability for bone remodeling contact rod models , 2005 .

[20]  J. Monnier,et al.  Existence and Uniqueness of a Solution to an Adaptive Elasticity Model , 1998 .

[21]  R. Glowinski,et al.  Numerical Methods for Nonlinear Variational Problems , 1985 .

[22]  M. Barboteu,et al.  A CLASS OF EVOLUTIONARY VARIATIONAL INEQUALITIES WITH APPLICATIONS IN VISCOELASTICITY , 2005 .

[23]  G S Beaupré,et al.  An approach for time‐dependent bone modeling and remodeling—theoretical development , 1990, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[24]  I. Figueiredo.  Approximation of bone remodeling models , 2005 .