Modélisation du transport couplé de fluide dans les milieux poreux : application aux os

Un modèle continu à deux échelles est développé afin de décrire le comportement chimio-electro-hydraulique du fluide interstitiel saturant un milieu poreux. Pour cela, le modèle couplé Moyne et Murad est modifié afin de tenir compte d'effets visqueux localisés à une échelle inférieure à la microéchelle. Le modèle à la microéchelle associe donc à une loi de Brinkman modifiée les équations de transport chimique et le bilan électrostatique. Après un changement de variables approprié, une procédure d'homogénéisation périodique est mise en \oe uvre pour dériver une description macroscopique. L'obtention d'une loi de Darcy couplée utilisable pour exprimer l'écoulement interstitiel dans des milieux à pores à géométrie cylindrique comme les os est finalement établie

Une approche multi-échelle pour la modélisation des fibrilles de l’os

Les propriétés mécaniques du tissu osseux sont dues à sa structure hiérarchique. À l’échelle de quelques micromètres, la structure que l’on peut observer est la fibrille de collagène minéralisée. Dans ce travail, nous présentons un modèle mécanique multi-échelle de la fibrille basé sur une technique d’homogénéisation variationnelle. Des simulations numériques montrent la dépendance des propriétés élastiques de la fibrille vis-à-vis de sa structure à l'échelle nanométrique

Influence of microcracks on interstitial fluid flows at the osteonal scale

Les microfissures sont connues comme un facteur important qui stimule le remodelage de l'os. Ce travail étudie les conséquences de l'apparition d'une microfissure sur le comportement hydraulique osteonal. Un modèle poroélastique de type Biot a été développé. Les résultats numériques obtenus montrent que les microfissures pourraient modifier significativement les champs de vitesses des écoulements à l'intérieur des ostéons avoisinants. Ceci permetrait d'expliquer le mécanisme de remodelage ciblé par lequel les microfissures sont naturellement détectées et réparées

Multiphysics modelling of transport phenomena within cortical bone treated as a biporous material

Interstitial fluid and ionic transport taking place in the fluid compartments of bone are thought to play a major role in bone mechanotransduction. In this study, we present a three-scale model of the multiphysical transport phenomena taking place within the vasculature porosity and the lacuno-canalicular network of cortical bone. These two fluid porosity levels exchange mass and ions through the permeable outer wall of the Haversian-Volkmann channels. Thus, coupled equations of electro-chemo-hydraulic transport are derived from the nanoscale of the canaliculi toward the cortical tissue, considering the intermediate scale of the intra-osteonal tissue. In particular, the Onsager reciprocity laws that govern the coupled transport are checked

De la détermination de la perméabilité lacuno-canaliculaire dans l'os

La perméabilité lacuno-canaliculaire régit l'écoulement interstitiel cortical, ce qui en fait un élément-clé quant à la compréhension du remodelage osseux. Dans cette étude est présentée une mise en perspective des écarts entre les prédictions théoriques (10e-20 m2) et les récentes mesures de ce paramètre (10e-24 m2) [1]. Pour illustrer notre propos, nous proposons une loi de Kozeny numérique donnant l'évolution de la perméabilité avec les paramètres de la microstructure osseuse. Cette loi pourra servir d'entrée dans un modèle de remodelage osseux. [1] Gardinier et al, 2010, Bone 46:107

A thermodynamics framework to describe bone remodeling: a 2D study

Introduction Bone is a living material which is continuously reorganized by bone cells in response to their mechanical and biochemical environment. This process, known as bone remodeling, is of major importance in everyday life, in case of fractures and to allow osseointegration phenomena around implants. Bone remodeling can be described as a stress- and chemistry-driven evolution of the mechanical properties of bone tissue. The interplay between mechanics and biochemistry, as well as the multiple scales involved in this process, represent serious challenges for the development of realistic bone adaptation models. This study describes a novel, thermodynamically sound model of bone remodeling which, while being mechanistic in nature, is able to account for the above issues.   Methods The theory of material remodeling [1] offers a firm basis to build a model of bone remodeling. Bone is described as an orthotropic elastic medium whose elastic properties evolve in time according to the prevailing mechanical and chemical stimuli. In this study, we focus on a special class of evolution, namely the stress-driven rotation of the elastic principal axes of the bony material. Our model is based on balance laws derived through a suitable statement of the virtual power principle, as well as the description of a constitutive theory. The latter is based on the definition of a strain-energy density depending only on the elastic strain, and on a formulation of the dissipation principle incorporating the dissipation due to remodeling. This modeling framework leads to a remodeling evolution law giving an explicit relationship coupling the dissipation related to the remodeling, the stress and strain tensors, the rotation of the material axes and its evolution. Remodeling equilibrium is achieved when material properties no longer evolve, corresponding to a stationary state of the rotation. It is worth noting that this model predicts the principal axes of the strain and stress tensors to be collinear at the remodeling equilibrium [2]. Thus, a physically sound condition for remodeling equilibrium [3] is recovered without any ad-hoc assumption.   Results and Discussion The model was studied in 2D where the rotation of the elastic principal axes is parameterized by an angle. In case of uniform stress/strain conditions, stable remodeling equilibrium states were found to correspond to strain energy minima (imposed displacements). Finite element simulations were also performed to study the evolution of the elastic principal axes resulting from different boundary conditions. The prediction of the rotation of the principal axes of the material in simple loading configurations is consistent with the superimposed boundary conditions confirming the ability of the model to simulate the material response to non-uniform stress configurations.   Conclusion A novel, thermodynamically sound model of bone remodeling was proposed. Preliminary results obtained in 2D are promising and show the potential of this approach. Model predictions for in vivo biomechanical loading configurations need to be further tested. Suitable experimental data will be identified. The particular case of tissue surrounding an implant will be also studied. Our model can also integrate the mechanobiological phenomena regulating bone remodeling. However, this would require a reliable description of the biochemical stimuli of bone remodeling. This matter is out of the scope of this paper and will be addressed in future works.   Acknowledgement This project has received funding from the European Research Council under the European Union's Horizon 2020 research and innovation program (grant agreement No 682001, project ERC Consolidator Grant 2015 BoneImplant).   References [1] DiCarlo A, Naili S, and Quiligotti S (2006) C. R. Mecanique, 334:651-661 [2] Sansalone V, Naili S, and Di Carlo A (2011) Comput Methods Biomech Biomed Engin, 14(s1):203-204 [3] Cowin S (1986) J Biomech Eng, 108:83-8

Influence des flux de calcium sur la contrainte de cisaillement agissant sur les ostéocytes dans l'os cortical

Nous avons modélisé l'écoulement du fluide encerclant les cellules osseuses mécano-sensibles (ostéocytes). Le but de cette étude est d'améliorer nos modèles précédents en incluant des flux de calcium apparaissant lors de la dissolution ou de la précipitation de la matrice osseuse. Même si ces flux ne semblent pas altérer de manière significative la vitesse du fluide, ils peuvent changer le cisaillement ressenti par les ostéocytes. Par conséquent, nous avons examiné la façon dont de tels échanges chimiques affectent l'écoulement du fluide et donc la mécano-sensibilité des ostéocytes