1
|
Giorgio I, dell'Isola F, Andreaus U, Misra A. An orthotropic continuum model with substructure evolution for describing bone remodeling: an interpretation of the primary mechanism behind Wolff's law. Biomech Model Mechanobiol 2023; 22:2135-2152. [PMID: 37542620 PMCID: PMC10613191 DOI: 10.1007/s10237-023-01755-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2023] [Accepted: 07/16/2023] [Indexed: 08/07/2023]
Abstract
We propose a variational approach that employs a generalized principle of virtual work to estimate both the mechanical response and the changes in living bone tissue during the remodeling process. This approach provides an explanation for the adaptive regulation of the bone substructure in the context of orthotropic material symmetry. We specifically focus upon the crucial gradual adjustment of bone tissue as a structural material that adapts its mechanical features, such as materials stiffnesses and microstructure, in response to the evolving loading conditions. We postulate that the evolution process relies on a feedback mechanism involving multiple stimulus signals. The mechanical and remodeling behavior of bone tissue is clearly a complex process that is difficult to describe within the framework of classical continuum theories. For this reason, a generalized continuum elastic theory is employed as a proper mathematical context for an adequate description of the examined phenomenon. To simplify the investigation, we considered a two-dimensional problem. Numerical simulations have been performed to illustrate bone evolution in a few significant cases: the bending of a rectangular cantilever plate and a three-point flexure test. The results are encouraging because they can replicate the optimization process observed in bone remodeling. The proposed model provides a likely distribution of stiffnesses and accurately represents the arrangement of trabeculae macroscopically described by the orthotropic symmetry directions, as supported by experimental evidence from the trajectorial theory.
Collapse
Affiliation(s)
- Ivan Giorgio
- Department of Civil, Construction-Architectural and Environmental Engineering (DICEAA), University of L'Aquila, 1, P.zza Ernesto Pontieri, Monteluco di Roio, L'Aquila, 67100, Italy.
- International Research Center for the Mathematics and Mechanics of Complex Systems (M &MoCS), University of L'Aquila, 1, P.zza Ernesto Pontieri, Monteluco di Roio, L'Aquila, 67100, Italy.
| | - Francesco dell'Isola
- Department of Civil, Construction-Architectural and Environmental Engineering (DICEAA), University of L'Aquila, 1, P.zza Ernesto Pontieri, Monteluco di Roio, L'Aquila, 67100, Italy
- International Research Center for the Mathematics and Mechanics of Complex Systems (M &MoCS), University of L'Aquila, 1, P.zza Ernesto Pontieri, Monteluco di Roio, L'Aquila, 67100, Italy
- Faculty of Mechanical and Industrial Engineering, Warsaw University of Technology, ul. Narbutta 85, Warsaw, 02-524, Poland
- CNRS Fellow, ENS Paris-Saclay, 4, avenue des Sciences, Gif-sur-Yvette, 91190, France
| | - Ugo Andreaus
- Department of Structural and Geotechnical Engineering (DISG), Università di Roma La Sapienza, 18, Via Eudossiana, Rome, 00184, Italy
| | - Anil Misra
- International Research Center for the Mathematics and Mechanics of Complex Systems (M &MoCS), University of L'Aquila, 1, P.zza Ernesto Pontieri, Monteluco di Roio, L'Aquila, 67100, Italy
- Civil, Environmental and Architectural Engineering Department (CEAE), The University of Kansas, 1530 W. 15th Street, Learned Hall, Lawrence, 66045-7609, Kansas, USA
| |
Collapse
|
2
|
Hayashi K, Yanagisawa T, Kishida R, Tsuchiya A, Ishikawa K. Gear-shaped carbonate apatite granules with a hexagonal macropore for rapid bone regeneration. Comput Struct Biotechnol J 2023; 21:2514-2523. [PMID: 37077175 PMCID: PMC10106487 DOI: 10.1016/j.csbj.2023.03.053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2022] [Revised: 03/31/2023] [Accepted: 03/31/2023] [Indexed: 04/08/2023] Open
Abstract
Synthetic bone grafts are in high demand owing to increased age-related bone disorders in the global aging population. Here, we report fabrication of gear-shaped granules (G-GRNs) for rapid bone healing. G-GRNs possessed six protrusions and a hexagonal macropore in the granular center. These were composed of carbonate apatite, i.e., bone mineral, microspheres with ∼1-μm micropores in the spaces between the microspheres. G-GRNs formed new bone and blood vessels (both on the granular surface and within the macropores) 4 weeks after implantation in the rabbit femur defects. The formed bone structure was similar to that of cancellous bone. The bone percentage in the defect recovered to that in a normal rabbit femur at week-4 post-implantation, and the bone percentage remained constant for the following 8 weeks. Throughout the entire period, the bone percentage in the G-GRN-implanted group was ∼10% higher than that of the group implanted with conventional carbonate apatite granules. Furthermore, a portion of the G-GRNs resorbed at week-4, and resorption continued for the following 8 weeks. Thus, G-GRNs are involved in bone remodeling and are gradually replaced with new bone while maintaining a suitable bone level. These findings provide a basis for the design and fabrication of synthetic bone grafts for achieving rapid bone regeneration.
Collapse
|
3
|
Sen A, Follet H, Sornay-Rendu E, Rémond Y, George D. Prediction of osteoporotic degradation of tibia human bone at trabecular scale. J Mech Behav Biomed Mater 2023; 139:105650. [PMID: 36657191 DOI: 10.1016/j.jmbbm.2023.105650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2022] [Revised: 11/18/2022] [Accepted: 01/01/2023] [Indexed: 01/11/2023]
Abstract
A theoretical numerical model is proposed to predict patient dependent osteoporotic bone degradation. The model parameters are identified through a particle swarm optimization algorithm and based on individual patient high resolution peripherical quantitative computer tomography (HRpQCT) scan data. The degradation model is based on cellular activity initiated by the elastic strain energy developed in the bone microstructure through patient's body weight. The macro (organ scale) and meso (trabecular scale) scale analyses are carried out and predicted bone volume fraction and microstructure evolution are compared with in-vivo experimental bone degradation for four elderly women over a period of 10 years. A significant correlation (r > 0.9) is observed between the model predictions and in-vivo experiments in all cases with an average deviation error of 1.46%. The model can easily be extended to other patients and provide good predictions for different population categories such as ethnicity, gender, age, etc.
Collapse
Affiliation(s)
- Ahmet Sen
- University of Strasbourg, CNRS, ICUBE Laboratory, Strasbourg, France
| | - Hélène Follet
- University Claude Bernard Lyon 1, INSERM, LYOS UMR 1033, 69008, Lyon, France.
| | - Elisabeth Sornay-Rendu
- University Claude Bernard Lyon 1, INSERM, LYOS UMR 1033, 69008, Lyon, France; Edouard Heriot Hospital, Hospices Civils of Lyon, 69437, Lyon, France
| | - Yves Rémond
- University of Strasbourg, CNRS, ICUBE Laboratory, Strasbourg, France
| | - Daniel George
- University of Strasbourg, CNRS, ICUBE Laboratory, Strasbourg, France.
| |
Collapse
|
4
|
A Proposal for a Novel Formulation Based on the Hyperbolic Cattaneo’s Equation to Describe the Mechano-Transduction Process Occurring in Bone Remodeling. Symmetry (Basel) 2022. [DOI: 10.3390/sym14112436] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
In this paper, we propose a model for the mechanical stimulus involved in the process of bone remodeling together with its evolution over time. Accumulated evidence suggests that bone remodeling could be interpreted as a feedback control process in which the mechanical state of the bone tissue is monitored, then appropriate signals are derived from the daily mechanical usage of the bone, these signals are transmitted into the surrounding region, and then they are detected by other agents whose purpose is to adapt the bone mass to the mechanical requirements of the environment. Therefore, we employ the diffusion equation for mass transport which is improved with Cattaneo’s correction to model the stimulus. This last improvement considers the effects of relaxation and non-locality, which we believe play essential roles in signaling messengers transport phenomena and are essential to match the evidence that suggests time-dependent excitations provide a more significant response at specific frequencies. To illustrate this particular behavior, numerical simulations have been performed in a 2D framework. The results fit the central aspect addressed, related to the dependency of the time of the adaptive process of bone, suggesting that our model is promising and deserves further investigation, both theoretical and experimental.
Collapse
|
5
|
Perier-Metz C, Corté L, Allena R, Checa S. A 3D in Silico Multi-Tissue Evolution Model Highlights the Relevance of Local Strain Accumulation in Bone Fracture Remodeling. Front Bioeng Biotechnol 2022; 10:835094. [PMID: 35433640 PMCID: PMC9008279 DOI: 10.3389/fbioe.2022.835094] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Accepted: 03/08/2022] [Indexed: 12/03/2022] Open
Abstract
Since 5–10% of all bone fractures result in non-healing situations, a thorough understanding of the various bone fracture healing phases is necessary to propose adequate therapeutic strategies. In silico models have greatly contributed to the understanding of the influence of mechanics on tissue formation and resorption during the soft and hard callus phases. However, the late-stage remodeling phase has not been investigated from a mechanobiological viewpoint so far. Here, we propose an in silico multi-tissue evolution model based on mechanical strain accumulation to investigate the mechanobiological regulation of bone remodeling during the late phase of healing. Computer model predictions are compared to histological data of two different pre-clinical studies of bone healing. The model predicted the bone marrow cavity re-opening and the resorption of the external callus. Our results suggest that the local strain accumulation can explain the fracture remodeling process and that this mechanobiological response is conserved among different mammal species. Our study paves the way for further understanding of non-healing situations that could help adapting therapeutic strategies to foster bone healing.
Collapse
Affiliation(s)
- Camille Perier-Metz
- Julius Wolff Institute, Berlin Institute of Health at Charité—Universitätsmedizin Berlin, Berlin, Germany
- Centre des Matériaux, MINES Paris–PSL, Paris, France
| | - Laurent Corté
- Centre des Matériaux, MINES Paris–PSL, Paris, France
- Chimie Moléculaire, Macromoléculaire et Matériaux, ESPCI Paris–PSL, Paris, France
| | - Rachele Allena
- Laboratoire Mathématiques and Interactions J. A. Dieudonné, UMR 7351 CNRS, Université Côte d’Azur, Nice, France
| | - Sara Checa
- Julius Wolff Institute, Berlin Institute of Health at Charité—Universitätsmedizin Berlin, Berlin, Germany
- *Correspondence: Sara Checa,
| |
Collapse
|
6
|
Mathai B, Dhara S, Gupta S. Orthotropic bone remodelling around uncemented femoral implant: a comparison with isotropic formulation. Biomech Model Mechanobiol 2021; 20:1115-1134. [PMID: 33768358 DOI: 10.1007/s10237-021-01436-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2020] [Accepted: 02/11/2021] [Indexed: 11/25/2022]
Abstract
Peri-prosthetic bone adaptation has usually been predicted using subject-specific finite element analysis in combination with remodelling algorithms and assuming isotropic bone material property. The objective of the study is to develop an orthotropic bone remodelling algorithm for evaluation of peri-prosthetic bone adaptation in the uncemented implanted femur. The simulations considered loading conditions from a variety of daily activities. The orthotropic algorithm was tested on 2D and 3D models of the intact femur for verification of predicted results. The predicted orthotropic directionality, based on principal stress directions, was in agreement with the trabecular orientation in a micro-CT data of proximal femur. The validity of the proposed strain-based algorithm was assessed by comparing the predicted results of the orthotropic model with those of the strain-energy-density-based isotropic formulation. Despite agreement in cortical densities [Formula: see text], the isotropic remodelling algorithm tends to predict relatively higher values around the distal tip of the implant as compared to the orthotropic model. Both formulations predicted 4-8% bone resorption in the proximal femur. A linear regression analysis revealed a significant correlation [Formula: see text] between the stresses and strains on the cortex of the proximal femur, predicted by the isotropic and orthotropic formulations. Despite reasonable agreement in peri-prosthetic bone density distributions, the quantitative differences with isotropic model predictions highlight the combined influences of bone orthotropy and mechanical stimulus in the adaptation process.
Collapse
Affiliation(s)
- Basil Mathai
- Department of Mechanical Engineering, Indian Institute of Technology Kharagpur, Kharagpur, West Bengal, 721 302, India
| | - Santanu Dhara
- School of Medical Science and Technology, Indian Institute of Technology Kharagpur, Kharagpur, West Bengal, 721 302, India
| | - Sanjay Gupta
- Department of Mechanical Engineering, Indian Institute of Technology Kharagpur, Kharagpur, West Bengal, 721 302, India.
| |
Collapse
|
7
|
Scerrato D, Bersani AM, Giorgio I. Bio-Inspired Design of a Porous Resorbable Scaffold for Bone Reconstruction: A Preliminary Study. Biomimetics (Basel) 2021; 6:18. [PMID: 33802227 PMCID: PMC8006156 DOI: 10.3390/biomimetics6010018] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2021] [Revised: 03/06/2021] [Accepted: 03/08/2021] [Indexed: 02/07/2023] Open
Abstract
The study and imitation of the biological and mechanical systems present in nature and living beings always have been sources of inspiration for improving existent technologies and establishing new ones. Pursuing this line of thought, we consider an artificial graft typical in the bone reconstruction surgery with the same microstructure of the bone living tissue and examine the interaction between these two phases, namely bone and the graft material. Specifically, a visco-poroelastic second gradient model is adopted for the bone-graft composite system to describe it at a macroscopic level of observation. The second gradient formulation is employed to consider possibly size effects and as a macroscopic source of interstitial fluid flow, which is usually regarded as a key factor in bone remodeling. With the help of the proposed formulation and via a simple example, we show that the model can be used as a graft design tool. As a matter of fact, an optimization of the characteristics of the implant can be carried out by numerical investigations. In this paper, we observe that the size of the graft considerably influences the interaction between bone tissue and artificial bio-resorbable material and the possibility that the bone tissue might substitute more or less partially the foreign graft for better bone healing.
Collapse
Affiliation(s)
- Daria Scerrato
- Dipartimento di Scienze di Base ed Applicate per l’Ingegneria (SBAI), University of Rome La Sapienza, 00161 Roma, Italy;
- International Research Center for the Mathematics and Mechanics of Complex Systems (M&MoCS), University of L’Aquila, 67100 L’Aquila, Italy;
| | - Alberto Maria Bersani
- International Research Center for the Mathematics and Mechanics of Complex Systems (M&MoCS), University of L’Aquila, 67100 L’Aquila, Italy;
- Dipartimento di Ingegneria Meccanica e Aerospaziale (DIMA), University of Rome La Sapienza, 00184 Roma, Italy
| | - Ivan Giorgio
- International Research Center for the Mathematics and Mechanics of Complex Systems (M&MoCS), University of L’Aquila, 67100 L’Aquila, Italy;
- Dipartimento di Ingegneria Civile, Edile-Architettura e Ambientale (DICEAA), University of L’Aquila, 67100 L’Aquila, Italy
| |
Collapse
|
8
|
Reisinger AG, Frank M, Thurner PJ, Pahr DH. A two-layer elasto-visco-plastic rheological model for the material parameter identification of bone tissue. Biomech Model Mechanobiol 2020; 19:2149-2162. [PMID: 32377934 PMCID: PMC7603462 DOI: 10.1007/s10237-020-01329-0] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2019] [Accepted: 04/13/2020] [Indexed: 11/29/2022]
Abstract
The ability to measure bone tissue material properties plays a major role in diagnosis of diseases and material modeling. Bone's response to loading is complex and shows a viscous contribution to stiffness, yield and failure. It is also ductile and damaging and exhibits plastic hardening until failure. When performing mechanical tests on bone tissue, these constitutive effects are difficult to quantify, as only their combination is visible in resulting stress-strain data. In this study, a methodology for the identification of stiffness, damping, yield stress and hardening coefficients of bone from a single cyclic tensile test is proposed. The method is based on a two-layer elasto-visco-plastic rheological model that is capable of reproducing the specimens' pre- and postyield response. The model's structure enables for capturing the viscously induced increase in stiffness, yield, and ultimate stress and for a direct computation of the loss tangent. Material parameters are obtained in an inverse approach by optimizing the model response to fit the experimental data. The proposed approach is demonstrated by identifying material properties of individual bone trabeculae that were tested under wet conditions. The mechanical tests were conducted according to an already published methodology for tensile experiments on single trabeculae. As a result, long-term and instantaneous Young's moduli were obtained, which were on average 3.64 GPa and 5.61 GPa, respectively. The found yield stress of 16.89 MPa was lower than previous studies suggest, while the loss tangent of 0.04 is in good agreement. In general, the two-layer model was able to reproduce the cyclic mechanical test data of single trabeculae with an root-mean-square error of 2.91 ± 1.77 MPa. The results show that inverse rheological modeling can be of great advantage when multiple constitutive contributions shall be quantified based on a single mechanical measurement.
Collapse
Affiliation(s)
- Andreas G Reisinger
- Division Biomechanics, Department of Anatomy and Biomechanics, Karl Landsteiner University of Health Sciences, Krems an der Donau, Austria.
- Institute of Lightweight Design and Structural Biomechanics, Vienna University of Technology, Vienna, Austria.
| | - Martin Frank
- Institute of Lightweight Design and Structural Biomechanics, Vienna University of Technology, Vienna, Austria
| | - Philipp J Thurner
- Institute of Lightweight Design and Structural Biomechanics, Vienna University of Technology, Vienna, Austria
| | - Dieter H Pahr
- Division Biomechanics, Department of Anatomy and Biomechanics, Karl Landsteiner University of Health Sciences, Krems an der Donau, Austria
- Institute of Lightweight Design and Structural Biomechanics, Vienna University of Technology, Vienna, Austria
| |
Collapse
|
9
|
FERNÁNDEZ MPEÑA, WITTE F, TOZZI G. Applications of X‐ray computed tomography for the evaluation of biomaterial‐mediated bone regeneration in critical‐sized defects. J Microsc 2020; 277:179-196. [DOI: 10.1111/jmi.12844] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2019] [Revised: 10/06/2019] [Accepted: 11/04/2019] [Indexed: 12/16/2022]
Affiliation(s)
- M. PEÑA FERNÁNDEZ
- Zeiss Global Centre, School of Mechanical and Design EngineeringUniversity of Portsmouth Portsmouth UK
| | - F. WITTE
- Biotrics Bioimplants GmbH Berlin Germany
| | - G. TOZZI
- Zeiss Global Centre, School of Mechanical and Design EngineeringUniversity of Portsmouth Portsmouth UK
| |
Collapse
|
10
|
Delpuech B, Nicolle S, Confavreux CB, Bouazza L, Clezardin P, Mitton D, Follet H. Failure Prediction of Tumoral Bone with Osteolytic Lesion in Mice. ADVANCED STRUCTURED MATERIALS 2020. [DOI: 10.1007/978-3-030-50464-9_2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
|
11
|
Della Corte A, Giorgio I, Scerrato D. A review of recent developments in mathematical modeling of bone remodeling. Proc Inst Mech Eng H 2019; 234:273-281. [DOI: 10.1177/0954411919857599] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
In this article, we summarize the developments in the mathematical modeling of the mechanics of bone and related biological phenomena. We will devote special attention to the results of the last 10–15 years, although we will cover some relevant classical work to better frame the more recent researches. We will propose a division of the literature based on the main aim of the model (mechanical/biomathematical) and the type of biological phenomena considered (stimulus, growth, cell population dynamics). Finally, we will suggest some possible directions for future investigations.
Collapse
Affiliation(s)
- Alessandro Della Corte
- International Research Center for the Mathematics and Mechanics of Complex Systems, University of L’Aquila, L’Aquila, Italy
| | - Ivan Giorgio
- International Research Center for the Mathematics and Mechanics of Complex Systems, University of L’Aquila, L’Aquila, Italy
- Department of Structural and Geotechnical Engineering, SAPIENZA Università di Roma, Rome, Italy
| | - Daria Scerrato
- International Research Center for the Mathematics and Mechanics of Complex Systems, University of L’Aquila, L’Aquila, Italy
| |
Collapse
|
12
|
Giorgio I, dell’Isola F, Andreaus U, Alzahrani F, Hayat T, Lekszycki T. On mechanically driven biological stimulus for bone remodeling as a diffusive phenomenon. Biomech Model Mechanobiol 2019; 18:1639-1663. [DOI: 10.1007/s10237-019-01166-w] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2018] [Accepted: 05/08/2019] [Indexed: 10/26/2022]
|
13
|
Othman MIA, Abd-Elaziz EM. Effect of initial stress and Hall current on a magneto-thermoelastic porous medium with microtemperatures. INDIAN JOURNAL OF PHYSICS 2019; 93:475-485. [DOI: 10.1007/s12648-018-1313-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/11/2017] [Accepted: 07/12/2018] [Indexed: 09/02/2023]
|
14
|
George D, Allena R, Rémond Y. Cell nutriments and motility for mechanobiological bone remodeling in the context of orthodontic periodontal ligament deformation. ACTA ACUST UNITED AC 2018. [DOI: 10.1016/j.jocit.2018.09.006] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
15
|
Kim JJ, Nam J, Jang IG. Computational study of estimating 3D trabecular bone microstructure for the volume of interest from CT scan data. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2018; 34:e2950. [PMID: 29218827 DOI: 10.1002/cnm.2950] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/16/2017] [Revised: 11/15/2017] [Accepted: 11/25/2017] [Indexed: 06/07/2023]
Abstract
Inspired by the self-optimizing capabilities of bone, a new concept of bone microstructure reconstruction has been recently introduced by using 2D synthetic skeletal images. As a preliminary clinical study, this paper proposes a topology optimization-based method that can estimate 3D trabecular bone microstructure for the volume of interest (VOI) from 3D computed tomography (CT) scan data with enhanced computational efficiency and phenomenological accuracy. For this purpose, a localized finite element (FE) model is constructed by segmenting a target bone from CT scan data and determining the physiological local loads for the VOI. Then, topology optimization is conducted with multiresolution bone mineral density (BMD) deviation constraints to preserve the patient-specific spatial bone distribution obtained from the CT scan data. For the first time, to our knowledge, this study has demonstrated that 60-μm resolution trabecular bone images can be reconstructed from 600-μm resolution CT scan data (a 62-year-old woman with no metabolic bone disorder) for the 4 VOIs in the proximal femur. The reconstructed trabecular bone includes the characteristic trabecular patterns and has morphometric indices that are in good agreement with the anatomical data in the literature. As for computational efficiency, the localization for the VOI reduces the number of FEs by 99%, compared with that of the full FE model. Compared with the previous single-resolution BMD deviation constraint, the proposed multiresolution BMD deviation constraints enable at least 65% and 47% reductions in the number of iterations and computing time, respectively. These results demonstrate the clinical feasibility and potential of the proposed method.
Collapse
Affiliation(s)
- Jung Jin Kim
- The Cho Chun Shik Graduate School of Green Transportation, Korea Advanced Institute of Science and Technology, 291, Daehak-ro, Yuseong-gu, Daejeon, 34141, South Korea
| | - Jimin Nam
- The Cho Chun Shik Graduate School of Green Transportation, Korea Advanced Institute of Science and Technology, 291, Daehak-ro, Yuseong-gu, Daejeon, 34141, South Korea
| | - In Gwun Jang
- The Cho Chun Shik Graduate School of Green Transportation, Korea Advanced Institute of Science and Technology, 291, Daehak-ro, Yuseong-gu, Daejeon, 34141, South Korea
| |
Collapse
|
16
|
Computational Investigation on the Biomechanical Responses of the Osteocytes to the Compressive Stimulus: A Poroelastic Model. BIOMED RESEARCH INTERNATIONAL 2018; 2018:4071356. [PMID: 29581973 PMCID: PMC5822791 DOI: 10.1155/2018/4071356] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/29/2017] [Revised: 12/02/2017] [Accepted: 12/19/2017] [Indexed: 11/17/2022]
Abstract
Osteocytes, the major type of bone cells embedded in the bone matrix and surrounded by the lacunar and canalicular system, can serve as biomechanosensors and biomechanotranducers of the bone. Theoretical analytical methods have been employed to investigate the biomechanical responses of osteocytes in vivo; the poroelastic properties have not been taken into consideration in the three-dimensional (3D) finite element model. In this study, a 3D poroelastic idealized finite element model was developed and was used to predict biomechanical behaviours (maximal principal strain, pore pressure, and fluid velocity) of the osteocyte-lacunar-canalicular system under 150-, 1000-, 3000-, and 5000-microstrain compressive loads, respectively, representing disuse, physiological, overuse, and pathological overload loading stimuli. The highest local strain, pore pressure, and fluid velocity were found to be highest at the proximal region of cell processes. These data suggest that the strain, pore pressure, and fluid velocity of the osteocyte-lacunar-canalicular system increase with the global loading and that the poroelastic material property affects the biomechanical responses to the compressive stimulus. This new model can be used to predict the mechanobiological behaviours of osteocytes under the four different compressive loadings and may provide an insight into the mechanisms of mechanosensation and mechanotransduction of the bone.
Collapse
|
17
|
George D, Allena R, Rémond Y. Mechanobiological stimuli for bone remodeling: mechanical energy, cell nutriments and mobility. Comput Methods Biomech Biomed Engin 2017; 20:91-92. [DOI: 10.1080/10255842.2017.1382876] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- D. George
- ICUBE, CNRS, University of Strasbourg, 67000 Strasbourg, France
| | - R. Allena
- Arts et Métiers ParisTech, LBM/Institut de Biomécanique Humaine Georges Charpak, 151 bd de l’ Hôpital 75013 Paris, France
| | - Y. Rémond
- ICUBE, CNRS, University of Strasbourg, 67000 Strasbourg, France
| |
Collapse
|
18
|
A Review of Mixture Theory for Deformable Porous Media and Applications. APPLIED SCIENCES-BASEL 2017. [DOI: 10.3390/app7090917] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
19
|
Spingarn C, Wagner D, Rémond Y, George D. Multiphysics of bone remodeling: A 2D mesoscale activation simulation. Biomed Mater Eng 2017; 28:S153-S158. [PMID: 28372290 DOI: 10.3233/bme-171636] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
In this work, we present an evolutive trabecular model for bone remodeling based on a boundary detection algorithm accounting for both biology and applied mechanical forces, known to be an important factor in bone evolution. A finite element (FE) numerical model using the Abaqus/Standard® software was used with a UMAT subroutine to solve the governing coupled mechanical-biological non-linear differential equations of the bone evolution model. The simulations present cell activation on a simplified trabeculae configuration organization with trabecular thickness of 200µm. For this activation process, the results confirm that the trabeculae are mainly oriented in the active direction of the principal mechanical stresses and according to the principal applied mechanical load directions. The trabeculae surface activation is clearly identified and can provide understanding of the different bone cell activations in more complex geometries and load conditions.
Collapse
Affiliation(s)
- C Spingarn
- Laboratoire des Sciences de l'Ingenieur, de l'Informatique et de l'Imagerie (Icube), Université de Strasbourg, CNRS, 67000 Strasbourg, France
| | - D Wagner
- Laboratoire des Sciences de l'Ingenieur, de l'Informatique et de l'Imagerie (Icube), Université de Strasbourg, CNRS, 67000 Strasbourg, France.,Faculté de Chirurgie Dentaire, Hôpitaux Universitaires de Strasbourg, 67000 Strasbourg, France
| | - Y Rémond
- Laboratoire des Sciences de l'Ingenieur, de l'Informatique et de l'Imagerie (Icube), Université de Strasbourg, CNRS, 67000 Strasbourg, France
| | - D George
- Laboratoire des Sciences de l'Ingenieur, de l'Informatique et de l'Imagerie (Icube), Université de Strasbourg, CNRS, 67000 Strasbourg, France
| |
Collapse
|
20
|
George D, Spingarn C, Dissaux C, Nierenberger M, Rahman RA, Rémond Y. Examples of multiscale and multiphysics numerical modeling of biological tissues. Biomed Mater Eng 2017; 28:S15-S27. [DOI: 10.3233/bme-171621] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Affiliation(s)
- Daniel George
- ICube Laboratory, University of Strasbourg, CNRS, 2 rue Boussingault, 67000 Strasbourg, France
- International Center MEMOCS, University del Aquila, Palazzo Caetani, Cisterna di Latina, Italy
| | - Camille Spingarn
- ICube Laboratory, University of Strasbourg, CNRS, 2 rue Boussingault, 67000 Strasbourg, France
| | - Caroline Dissaux
- Service of Maxillo-Facial Surgery and Stomatology, Strasbourg University Hospital, 1 place de l’Hôpital, 67000 Strasbourg, France
| | - Mathieu Nierenberger
- ICube Laboratory, University of Strasbourg, CNRS, 2 rue Boussingault, 67000 Strasbourg, France
| | | | - Yves Rémond
- ICube Laboratory, University of Strasbourg, CNRS, 2 rue Boussingault, 67000 Strasbourg, France
- International Center MEMOCS, University del Aquila, Palazzo Caetani, Cisterna di Latina, Italy
| |
Collapse
|