Numerical analysis of a diffusive strain-adaptive bone remodelling theory

详细信息	查看全文 | 推荐本文 |

• 作者：J.R. Ferná ; ndez^a ; ^{jose.fernandez@uvigo.es} ; J.M. Garcí ; a-Aznar^b ; R. Martí ; nez^c
• 关键词：Bone remodelling ; Diffusive strain-adaptivity ; Error estimates ; Finite elements ; Numerical simulations
• 刊名：International Journal of Solids and Structures
• 出版年：2012
• 期刊代码：93_00207683
• 类别：et
• 出版时间：August, 2012
• 卷：49
• 期：15-16
• 页码：2085-2093
• 文件大小：1065 K

摘要

In this paper we revisit a strain-adaptive bone remodelling problem, assuming that the rate of the apparent density at a particular location is described as a local objective function and depending on a particular stimulus at that location. Normally, continuum mathematical descriptions of adaptive bone remodelling can lead to discontinuous solutions on the global apparent density distribution. To improve this numerical solution, in this work, as the main novelty, we include the diffusion of the bone remodelling into the model. The variational problem is written as a coupled system of a nonlinear variational equation for the displacement field and a parabolic elliptic variational inequality for the apparent density. An existence and uniqueness result is stated. Then, a fully discrete problem is introduced by using the finite element method to approximate the spatial variable and an Euler scheme to discretize the time derivatives. A priori error estimates are proved from which, under adequate additional regularity conditions, the linear convergence of the algorithm is deduced. Finally, some numerical simulations are presented to demonstrate the accuracy of the approximation and the behaviour of the solution.