1
|
Quexada D, Ramtani S, Trabelsi O, Marquez K, Marie-Christine, Linero Segrera DL, Duque-Daza C, Garzón Alvarado DA. A unified framework of cell population dynamics and mechanical stimulus using a discrete approach in bone remodelling. Comput Methods Biomech Biomed Engin 2023; 26:399-411. [PMID: 35587027 DOI: 10.1080/10255842.2022.2065201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
Abstract
Multiphysics models have become a key tool in understanding the way different phenomenon are related in bone remodeling and various approaches have been proposed, yet, to the best of the author's knowledge there is no model able to link a cell population model with a mechanical stimulus model using a discrete approach, which allows for an easy implementation. This article couples two classical models, the cell population model from Komarova and the Nackenhorst model in a 2D domain, where correlations between the mechanical loading and the cell population dynamics can be established, furthermore the effect of different paracrine and autocrine regulators is seen on the overall density of a portion of trabecular bone. A discretization is performed using frame 1D finite elements, representing the trabecular structure. The Nackenhorst model is implemented by using the finite element method to calculate the strain energy as the main mechanical stimulus that determines the bone mass density evolution in time. This density is normalized to be added to the bone mass percentage proposed by the Komarova model, where coupling terms have been added as well that guarantee a stable response. In the simulations, the equations were solved employing the finite element method with a user subroutine implemented in ABAQUS (2017) and by applying a direct formulation. The methodology presented can model the cell dynamics occurring in bone remodelling in accordance with the asynchronous nature of this process, yet allowing to differentiate zones with higher density, the main trabecular groups are obtained for the proximal femur. Finally, the model is tested in pathological cases, such as osteoporosis and osteopetrosis, yielding results similar to the pathology behavior. Furthermore, the discrete modelling technique is shown to be of use in this particular application.
Collapse
Affiliation(s)
- Diego Quexada
- Université de technologie de Compiègne, CNRS, Biomechanics and Bioengineering, Centre de Recherche Royallieu, CS 60 319 - 60 203 Compiègne Cedex, France.,Numerical Methods and Modeling Research Group (GNUM), Universidad Nacional de Colombia
| | - Salah Ramtani
- Université Sorbonne Paris Nord, Laboratoire CSPBAT, équipe LBPS, CNRS (UMR 7244), Institut Galilée, F93430, Villetaneuse, France
| | - Olfa Trabelsi
- Université de technologie de Compiègne, CNRS, Biomechanics and Bioengineering, Centre de Recherche Royallieu, CS 60 319 - 60 203 Compiègne Cedex, France
| | | | - Marie-Christine
- Université de technologie de Compiègne, CNRS, Biomechanics and Bioengineering, Centre de Recherche Royallieu, CS 60 319 - 60 203 Compiègne Cedex, France
| | | | - Carlos Duque-Daza
- Numerical Methods and Modeling Research Group (GNUM), Universidad Nacional de Colombia
| | | |
Collapse
|