- Research
- Open Access
On a new law of bone remodeling based on damage elasticity: a thermodynamic approach
- Ahmed Idhammad^{1}Email author and
- Abdelmounaïm Abdali^{1}
https://doi.org/10.1186/1742-4682-9-51
© Idhammad and Abdali; licensee BioMed Central Ltd. 2012
- Received: 26 September 2012
- Accepted: 18 November 2012
- Published: 29 November 2012
Abstract
Background
Bone tissue is the main element of the human skeleton and is a dynamic tissue that is continuously renewed by bone-resorbing osteoclasts and bone-forming osteoblasts.
The bone is also capable of repairing itself and adapting its structure to changes in its load environment through the process of bone remodeling.
Therefore, this phenomenon has been gaining increasing interest in the last years and many laws have been developed in order to simulate this process.
Results
In this paper, we develop a new law of bone remodeling in the context of damaged elastic by applying the thermodynamic approach in the case of small perturbations.
The model is solved numerically by a finite difference method in the one-dimensional bone structure of a n-unit elements model.
Conclusion
In addition, several numerical simulations are presented that confirm the accuracy and effectiveness of the model.
Keywords
- Numerical simulation
- Thermodynamic approach
- Small perturbations hypothesis
- Bone remodeling
- Bone density
- Damage
- Fatigue
- Osteocyte
- Elasticity
- N-unit elements
Introduction
Bone is a living material that constantly replaces old tissue with new in a process called remodeling. It is also able to respond adaptively to its environment [1, 2].
This process is assumed to repair the microdamage and maintain bone quality; and also occurs continuously with each cycle lasting 4 to 7 months [6].
Over the past, the progress made in understanding bone remodeling, through two groups: phenomenological approach and thermodynamic approach, has been truly spectacular.
The thermodynamic approach was initiated first by chemists and was applied to continuum mechanics by Eckart and Biot around 1950. Furthermore, this approach was adopted by introducing state variables [7, 8] and thermodynamic potential which allows to define associated variables chosen for the study of the phenomenon [9].
Many theoretical and experimental works have been published emphasizing the importance of applying this thermodynamical approach to a bone structure [10–13].
In this study, we adopt thermodynamic approach of irreversible processes to get a new law of bone remodeling taking into account the bone density described by the law of Mullender et al. [14] and the damage evolution in the bone proposed by Martin [15].
The found equation is solved by the finite difference method (FDM) in the one-dimensional bone structure of a n-unit elements model.
Finally, we present some examples of numerical simulation results.
Fundamental assumptions
- 1)
The bone is considered as a linear-elastic, isotropic and inhomogeneous material.
- 2)
The external mechanical loading acts as a stimulus for bone remodeling.
- 3)
The n-unit elements model is applied to the damaged-bone structure in the one-dimensional case.
- 4)
The small perturbations hypothesis (displacements and their time and spatial variations are small).
- 5)
The state coupling of damage with elastic strain.
- 6)
The remodeling processes can be considered isothermal, adiabatic and without internal generation of heat.
- 7)The decoupling of the thermodynamic potential is assumed, such that:$\psi =\psi \left({\epsilon}^{e},D,\varphi \right)={\psi}_{e}\left({\epsilon}^{e},D\right)+{\psi}_{r}\left(\varphi \right)$
With:
ψ_{ e }(ε^{ e }, D) is the thermodynamic potential depending on the elastic strain tensor ε^{ e } and the damage variable D.
ψ_{ r }(ϕ) is the thermodynamic potential depending on the bone density ϕ.
Thermodynamic approach
The thermodynamics of irreversible processes allows the modeling of different materials behavior. This is accomplished by defining the state variables and the state potential and also the dissipation potential [16].
The general theory of adaptive damaged-elastic materials and general framework of continuum thermodynamics is considered to find a new law of bone remodeling [10, 17].
Thermodynamic variables
State variables | Associated variables | |
---|---|---|
Observable | Internal | |
T | s | |
ε | σ | |
V _{ k } | A_{ k } |
Within the hypothesis of small strains and small displacements, the state variables, observable and internal, are chosen in accordance with the physical mechanisms of deformation and degradation of the bone as follows [18].
Observable variables:
– ε is the total strain tensor associated with the stress tensor σ.
– T is the temperature associated with the specific entropy s.
Internal variables:
– ε^{e} is the elastic strain tensor associated with the stress tensor σ.
– D is the damage associated with a variable $\overline{\text{Y}}$.
– ϕ is the bone density associated with the bone remodeling variable R.
Chart of thermodynamic variables
State variables | Associated variables | |
---|---|---|
Observable | Internal | |
T | s | |
ε | σ | |
ε ^{ e } | σ | |
D | $\overline{\text{Y}}$ | |
ϕ | - R |
We assume a bone remodeling variable, which is characterized by:
R > 0 in the case of the formation phase.
R = 0 in the case of the equilibrium phase.
R < 0 in the case of the resorption phase.
We postulate the existence of a thermodynamic potential from which the state laws can be derived [16, 18, 19].
The state potential: ψ = ψ(ε^{ e }, D, ϕ)
With:
ψ_{ e }(ε^{ e }, D) is the thermodynamic potential depending on the elastic strain tensor and the damage variable.
ψ_{ r }(ϕ) is the thermodynamic potential depending on the bone density.
The development of the Clausius-Duhem inequality
We note that:
$\rho \xb7\stackrel{\xb7}{\mathit{\psi}}=\rho \frac{\partial {\psi}_{e}}{\partial {\epsilon}^{e}}:{\stackrel{\xb7}{\mathit{\epsilon}}}^{e}+\rho \frac{\partial {\psi}_{e}}{\partial D}.\stackrel{\xb7}{D}+\rho \frac{\partial {\psi}_{r}}{\partial \varphi}.\stackrel{\xb7}{\mathit{\varphi}}$ with ε=ε^{ e }
(the thermo-elasticity law) [17]
Discussion
In the case of constant damage: $\stackrel{\xb7}{D}=0$
The inequality (5) gives $+\text{R}.\stackrel{\xb7}{\varphi}\ge 0$
In the resorption area:
In the formation area:
Equilibrium area (dead zone):
The conservation of energy equation
With:
– e is the specific internal energy.
– r is the internal heat source.
– q is the heat flux.
We introduce the specific heat (capacity) defined by:
$\text{C}=\text{T}\frac{\partial \text{s}}{\partial \text{T}}$[16, 17] and taking into account Fourier's law $\overrightarrow{q}=-k.g\overrightarrow{ra}dT$[17]
Assumptions
The classical heat equation corresponds to a process: [17]
Without internal generation of heat created by the external sources: r=0.
With adiabatic evolution: k. ΔT = 0.
With isothermal transformation: $\stackrel{\xb7}{T}=0$ Therefore $s=\frac{C}{T}.\stackrel{\xb7}{T}=0$ Then s=0.
We have: $\overline{\text{Y}}=-Y$
(the strain energy release rate) [17]
We have also: $\rho .{\psi}_{e}=\frac{1}{2}.\left(1-D\right).E.{\stackrel{\xb7}{\mathit{\epsilon}}}^{e}:{\stackrel{\xb7}{\mathit{\epsilon}}}^{e}$ (the strain energy) [17] and the equivalent constraint σ_{ eq } is written by ${\sigma}_{\mathit{eq}}=\rho \frac{\partial {\psi}_{e}}{\partial {\epsilon}^{e}}=\left(1-D\right).E.{\stackrel{\xb7}{\mathit{\epsilon}}}^{e}$ (The thermo-elasticity law)
With: σ_{ eq } = σ in the one-dimensional case.
The Young's modulus of the bone which is an isotropic material and inhomogeneous is expressed as:
E = (1 − D).E_{0} with E_{0} = c.ϕ^{ α }[22, 23]
c=100 and α =3 are two constants characteristic of the bone
This equation represents the new law of bone remodeling developed by applying the thermodynamic approach in the context of damaged elastic.
In this study, we introduce the law of damage developed by Martin [15] which shows that the damage in human cortical bone can grow exponentially until the fatigue failure [15, 24, 25].
The evolution law for the damage is expressed as:
$D={D}_{0}.{e}^{{f}_{d}t}$ Then, $\frac{\partial D}{\partial t}=\stackrel{\xb7}{D}={f}_{d}D$
With:
f_{ d } : the fatigue life of the bone devoid of the remodeling [26]
D_{0} : the initial damage
t : the time
Numerical resolution
- 1.
a constant $\frac{{\sigma}^{2}{f}_{d}}{2c}$
- 2.
a damage function $\frac{D}{{\left(1-D\right)}^{3}}$
- 3.
a function of bone density $\frac{1}{{\varphi}^{\alpha}\stackrel{\xb7}{\varphi}}$
With :
– ϕ_{ min } ≤ ϕ ≤ ϕ_{ max }
– ϕ_{ min } is the density of completely resorbed bone
– ϕ_{ max } is the maximum density defined for a compact bone
– τ is a positive constant related to the reaction time of bone tissue (constant of bone remodeling)
– 1≤ i ≤n
– ϕ_{ i } density of bone tissue of element i
– m (m ≤ n) is the total number of osteocytes in the solid
– I_{ k } (1 ≤ k ≤m) corresponds to the series of numbers of the elements containing an osteocyte
– S_{ k } represents the density of deformation energy in I_{ k }
– S_{ ref } reference stimulus value
– β is a parameter reflecting the intensity of the stimulus cell
– d is the normalization factor limiting the area of influence of osteocyte
– d(i,I_{ k }) is the distance between the centers of geometric element i and the element I_{ k }
To solve the previous equation (16), we use the method of finite difference with an implicit scheme and the fixed point method [23, 27, 28].
Results and discussion
Values of the parameters used during the numerical simulations
Data | Symbol | Values | Units |
---|---|---|---|
Maximum density | ϕ_{max} | 1.75 | g/cm^{3} |
Minimal density | ϕ_{min} | 0.01 | g/cm^{3} |
Initial density | ϕ_{0} | 0.6 | g/cm^{3} |
The step of time | Δt | 5.10^{-3} | UT |
The total force | F | 10 | N |
The distance between 2 centers | d | 25 | mm |
Reference stimulus value | S_{ref} | 0.04 | MPA |
The fatigue life of the bone | f_{d} | 3 | years |
n-unit elements of the bone fragment | 50 |
Figure 4 shows the temporal evolution of the variable of bone remodeling in the case of a uniform distribution of the osteocyte cells (n=m=50), and of another heterogeneous case (n≠m with m=30 in the central package of the osteocytes).
The curves consists of three key periods. The first period of the curves corresponds to the resorption phase, where the variable of the bone remodeling was negative. The second period exhibits the formation phase, where the variable of the bone remodeling was positive. The third period defined as the interval between the resorption phase and the formation phase, which the variable of the bone remodeling reached zero.
The resorption phase takes approximately 18 days, which is then followed by an equilibrium phase that can last for up to 10 days and finally by the formation phase from 17 to 35 days. This is in agreement with results from the literature [6, 29, 30].
By comparing the curves in Figure 4 to the graph proposed by Terrier et al. [20], we see a good agreement. Furthermore, the curves that are found have a nonlinear shape.
Conclusion
In this paper, we proposed a thermodynamic approach in small perturbations for bone remodeling process.
The adopted model takes into consideration both the bone density and the damage and gives a new law of bone remodeling. Then, the governing equation of the process was solved by the finite difference method in the one-dimensional bone structure with n-unit elements model.
The numerical results obtained are in accordance with the experimental results found in the literature.
Authors’ information
Ahmed Idhammad - He received his engineering degree status at the National School of Mineral Industry in Rabat, Morocco. He is curzently a Doctoral student at the Faculty of Sciences and Technics in Marrakech, Morocco. His research interests are numerical simulation, biomechanics, bone remodeling, thermodynamic, fatigue and damage.
Abdelmounaïm Abdali - PhD in Solid Mechanics and Structures in University of Amiens in 1996, France. He is a Professor in computer science at the University Cadi Ayyad, Faculty of Sciences and Technics, Marrakech, Morocco. Member at the Laboratory of Applied Mathematics and Computer Science (LAMAI) Marrakech, Morocco. His research interests are numerical simulation, biomechanics, bone remodeling and damage, computer science, DTN Network.
Declarations
Acknowledgements
We are particularly grateful to the editor in chief Dr. Paul S. Agutter for his thoughtful suggestions for improving the language and style of the manuscript.
We equally thankful to all members of our laboratory team for stimulating discussions.
Authors’ Affiliations
References
- Rüberg T: Computer Simulation of Adaptive Bone Remodeling. 2003, Spain: Master’s thesis, Technische Universität Braunschweig, Centro Politécnico Superior, Universidad de ZaragozaGoogle Scholar
- García-Aznar JM, Rueberg T, Doblaré M: A bone remodeling model coupling microdamage growth and repair by 3D BMU-activity. Biomech Model Mechanobiol. 2005, 4: 147-167. 10.1007/s10237-005-0067-x.View ArticlePubMedGoogle Scholar
- Comfort P, Abrahamson E: Sports Rehabilitation and Injury Prevention. 2010, UK: Wiley-Blackwell Publishers, John Wiley & Sons Ltd, 114-115. 1stView ArticleGoogle Scholar
- Mohri T, Hanada K, Ozawa H: Coupling of resorption and formation on bone remodeling sequence in orthodontic tooth movement: A histochemical study. J Bone Miner Metab. 1991, 9: 57-69. 10.1007/BF02377987.View ArticleGoogle Scholar
- Raggatt LJ, Partridge NC: Cellular and Molecular Mechanisms of Bone Remodeling. J Biol Chem. 2010, 285: 25103-25108. 10.1074/jbc.R109.041087.PubMed CentralView ArticlePubMedGoogle Scholar
- Riggs BL, Parfitt AM: Drugs Used to Treat Osteoporosis: The Critical Need for a Uniform Nomenclature Based on Their Action on Bone Remodeling. J Bone Miner Res. 2005, 20: 177-184.View ArticlePubMedGoogle Scholar
- Coleman BD, Gurtin ME: Thermodynamics with internal state variables. J Chem Phys. 1967, 47: 597-613. 10.1063/1.1711937.View ArticleGoogle Scholar
- Coleman BD, Noll W: The thermodynamic of elastic materials with heat conduction and viscosity. Arch Ration Mech Anal. 1963, 13: 167-178. 10.1007/BF01262690.View ArticleGoogle Scholar
- Lemaitre J, Chaboche JL: Mechanics of solid materials. 1990, UK: Cambridge University PressView ArticleGoogle Scholar
- Ramtani S, Zidi M: A theoretical model of the effect of continuum damage on a bone adaptation model. J Biomech. 2001, 34: 471-479. 10.1016/S0021-9290(00)00215-3.View ArticlePubMedGoogle Scholar
- Doblaré M, Garcıa JM, Gomez MJ: Modelling bone tissue fracture and healing: a review. Eng Fract Mech. 2004, 71: 1809-1840. 10.1016/j.engfracmech.2003.08.003.View ArticleGoogle Scholar
- Kuhl E, Steinmann P: Theory and numerics of geometrically non-linear open system mechanics. Int J Numer Meth Eng. 2003, 58: 1593-1615. 10.1002/nme.827.View ArticleGoogle Scholar
- Hoger A, Lubarda VA: On the mechanics of solids with a growing mass. Int J Solids Struct. 2002, 39: 4627-4664. 10.1016/S0020-7683(02)00352-9.View ArticleGoogle Scholar
- Mullender M, Huiskes R, Weinans H: A physiological approach to the simulation of bone remodeling as self organizational control process. J Biomech. 1994, 27: 1389-1394. 10.1016/0021-9290(94)90049-3.View ArticlePubMedGoogle Scholar
- Martin RB: Fatigue damage, remodeling, and the minimization of skeletal weight. J Theor Biol. 2003, 220: 271-276. 10.1006/jtbi.2003.3148.View ArticlePubMedGoogle Scholar
- Lemaitre J, Desmorat R: Engineering Damage Mechanics. 2005, Berlin: SpringerGoogle Scholar
- Lemaitre J, Chaboche JL, Benallal A, Desmorat R: Mécanique des matériaux solides. 2009, Paris: DunodGoogle Scholar
- Lemaitre J: A Course on Damage Mechanics. 1992, NY: Springer-VerlagView ArticleGoogle Scholar
- Barbero EJ, Greco F, Lonetti P: Continuum Damage-healing Mechanics with Application to Self-healing Composites. Int J Damage Mech. 2005, 14: 51-81. 10.1177/1056789505045928.View ArticleGoogle Scholar
- Terrier A, Rakotomanana RL, Ramaniraka AN, Leyvraz PF: Adaptation models of anisotropic bone. Comput Meth Biomech Biomed Eng. 1997, 1: 47-59. 10.1080/01495739708936694.View ArticleGoogle Scholar
- Wriggers P: Computational Contact Mechanics. 2006, NY: Springer-Verlag, 2ndView ArticleGoogle Scholar
- Currey JD: The effect of porosity and mineral content on theYoung’s modulus elasticity of compact bone. J Biomech. 1988, 21: 131-139. 10.1016/0021-9290(88)90006-1.View ArticlePubMedGoogle Scholar
- Idhammad A, Abdali A, Bussy P: Numerical simulation of the process of bone remodeling in the context of damaged elastic. Int J Adv Sci Tech. 2011, 37: 87-98.Google Scholar
- Carter DR, Hayes WC, Schurman DJ: Fatigue life of compact bone-II. Effects of microstructure and density. J Biomech. 1976, 9: 211-214. 10.1016/0021-9290(76)90006-3.View ArticlePubMedGoogle Scholar
- Pattin CA, Caler WE, Carter DR: Cyclic mechanical property degradation during fatigue loading on cortical bone. J Biomech. 1996, 29: 69-79. 10.1016/0021-9290(94)00156-1.View ArticlePubMedGoogle Scholar
- Taylor D, Kuiper JH: The prediction of stress fractures using a stressed volume concept. J Orthop Res. 2001, 19: 919-926. 10.1016/S0736-0266(01)00009-2.View ArticlePubMedGoogle Scholar
- Abdali A, Almoatassime H, Errafay A: Extension of the law of Mullender within the viscoelastic framework for the digital simulation of the process of bone remodeling. Int J Math Stat. 2010, 6: 1-10. 10.3844/jmssp.2010.1.3.View ArticleGoogle Scholar
- Idhammad A, Abdali A: Numerical simulation by finite difference of the problem of bone remodeling: case of elasticity with damage. Proceedings of the International Congress on Numerical Analysis and Scientific Computing with Applications in Sciences and Engineering: 19–20 April 2011. 2011, Morocco: Settat, 183-186.Google Scholar
- Lee TC, Staines A, Taylor D: Bone adaptation to load: microdamage as a stimulus for bone remodelling. J Anat. 2002, 201: 437-446. 10.1046/j.1469-7580.2002.00123.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Komarova SV, Smith RJ, Dixon SJ, Sims SM, Wahl LM: Mathematical model predicts a critical role for osteoclast autocrine regulation in the control of bone remodeling. Bone. 2003, 33: 206-215. 10.1016/S8756-3282(03)00157-1.View ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.