1
|
Stolle J, Harper CM, Voegele KK, Najafi AR, Taheri M, MacBain J, Siegler S. A statistical analysis of human talar shape and bone density distribution. J Orthop Res 2024; 42:1780-1790. [PMID: 38483072 DOI: 10.1002/jor.25835] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Revised: 02/16/2024] [Accepted: 02/26/2024] [Indexed: 07/04/2024]
Abstract
The shape of the talus, its internal structure, and its mechanical properties are important in determining talar behavior during loading, which may be significant for the design of surgical tools and implants. Although recent studies using statistical shape modeling have described quantitative talar external shape variation, no similar quantitative study exists to describe the density distribution of internal talar structure. The goal of this study is to quantify statistical variation in talar shape and density to benefit the design of talar implants. To this end, weight-bearing computed tomography (CT) scans of the ankle were collected in neutral, bilateral standing posture, and three-dimensional models were generated for each talus. Local density derived from the Hounsfield unit of each CT voxel was extracted. A weighted spherical harmonic analysis was performed to quantify the talar external shape. One hundred and seventy-nine volumes of interest were placed in the same relative position within each talus to quantify the talar density. Additionally, a finite element analysis (FEA) was conducted on a talus with both heterogeneous and homogeneous material properties to observe the effect of these properties on the stress and strain response. Significant differences were found in the talar density in sex and age, as well as in the stress and strain response between homogeneous and heterogeneous FEA. These differences show the importance of considering heterogeneity when examining the load response of tarsal bones.
Collapse
Affiliation(s)
- Jordan Stolle
- Department of Mechanical Engineering and Mechanics, Drexel University, Philadelphia, Pennsylvania, USA
| | - Christine M Harper
- Department of Biomedical Sciences, Cooper Medical School of Rowan University, Camden, New Jersey, USA
| | - Kristyn K Voegele
- Department of Geology, Rowan University School of Earth and Environment, Glassboro, New Jersey, USA
| | - Ahmad R Najafi
- Department of Mechanical Engineering and Mechanics, Drexel University, Philadelphia, Pennsylvania, USA
| | - Mehrangiz Taheri
- Department of Mechanical Engineering and Mechanics, Drexel University, Philadelphia, Pennsylvania, USA
| | - Joshua MacBain
- Department of Mechanical Engineering and Mechanics, Drexel University, Philadelphia, Pennsylvania, USA
| | - Sorin Siegler
- Department of Mechanical Engineering and Mechanics, Drexel University, Philadelphia, Pennsylvania, USA
| |
Collapse
|
2
|
Bregoli C, Biffi CA, Tuissi A, Buccino F. Effect of trabecular architectures on the mechanical response in osteoporotic and healthy human bone. Med Biol Eng Comput 2024:10.1007/s11517-024-03134-8. [PMID: 38822996 DOI: 10.1007/s11517-024-03134-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2023] [Accepted: 05/18/2024] [Indexed: 06/03/2024]
Abstract
Research at the mesoscale bone trabeculae arrangement yields intriguing results that, due to their clinical resolution, can be applied in clinical field, contributing significantly to the diagnosis of bone-related diseases. While the literature offers quantitative morphometric parameters for a thorough characterization of the mesoscale bone network, there is a gap in understanding relationships among them, particularly in the context of various bone pathologies. This research aims to bridge these gaps by offering a quantitative evaluation of the interplay among morphometric parameters and mechanical response at mesoscale in osteoporotic and non-osteoporotic bones. Bone mechanical response, dependent on trabecular arrangement, is defined by apparent stiffness, computationally calculated using the Gibson-Ashby model. Key findings indicate that: (i) in addition to bone density, measured using X-ray absorptiometry, trabecular connectivity density, trabecular spacing and degree of anisotropy are crucial parameters for characterize osteoporosis state; (ii) apparent stiffness values exhibit strong correlations with bone density and connectivity density; (iii) connectivity density and degree of anisotropy result the best predictors of mechanical response. Despite the inherent heterogeneity in bone structure, suggesting the potential benefit of a larger sample size in the future, this approach presents a valuable method to enhance discrimination between osteoporotic and non-osteoporotic samples.
Collapse
Affiliation(s)
- Chiara Bregoli
- National Research Council, CNR-ICMTE, Lecco, Italy.
- Mechanical Engineering Department, Politecnico Di Milano, Milano, Italy.
| | | | | | - Federica Buccino
- Mechanical Engineering Department, Politecnico Di Milano, Milano, Italy
- IRCCS Istituto Ortopedico Galeazzi, Milan, Italy
| |
Collapse
|
3
|
Yang J, Pei Q, Wu X, Dai X, Li X, Pan J, Wang B. Stress reduction through cortical bone thickening improves bone mechanical behavior in adult female Beclin-1 +/- mice. Front Bioeng Biotechnol 2024; 12:1357686. [PMID: 38600946 PMCID: PMC11004267 DOI: 10.3389/fbioe.2024.1357686] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2023] [Accepted: 03/13/2024] [Indexed: 04/12/2024] Open
Abstract
Fragility fractures, which are more prevalent in women, may be significantly influenced by autophagy due to altered bone turnover. As an essential mediator of autophagy, Beclin-1 modulates bone homeostasis by regulating osteoclast and chondrocyte differentiation, however, the alteration in the local bone mechanical environment in female Beclin-1+/- mice remains unclear. In this study, our aim is to investigate the biomechanical behavior of femurs from seven-month-old female wild-type (WT) and Beclin-1+/- mice under peak physiological load, using finite element analysis on micro-CT images. Micro-CT imaging analyses revealed femoral cortical thickening in Beclin-1+/- female mice compared to WT. Three-point bending test demonstrated a 63.94% increase in whole-bone strength and a 61.18% increase in stiffness for female Beclin-1+/- murine femurs, indicating improved biomechanical integrity. After conducting finite element analysis, Beclin-1+/- mice exhibited a 26.99% reduction in von Mises stress and a 31.62% reduction in maximum principal strain in the femoral midshaft, as well as a 36.64% decrease of von Mises stress in the distal femurs, compared to WT mice. Subsequently, the strength-safety factor was determined using an empirical formula, revealing that Beclin-1+/- mice exhibited significantly higher minimum safety factors in both the midshaft and distal regions compared to WT mice. In summary, considering the increased response of bone adaptation to mechanical loading in female Beclin-1+/- mice, our findings indicate that increasing cortical bone thickness significantly improves bone biomechanical behavior by effectively reducing stress and strain within the femoral shaft.
Collapse
Affiliation(s)
- Jiaojiao Yang
- Key Laboratory for Biorheological Science and Technology of Ministry of Education, College of Bioengineering, Chongqing University, Chongqing, China
- Institute of Life Sciences, College of Basic Medicine, Chongqing Medical University, Chongqing, China
| | - Qilin Pei
- Institute of Life Sciences, College of Basic Medicine, Chongqing Medical University, Chongqing, China
- Department of Biomedical Engineering, Fourth Military Medical University, Xi’an, China
| | - Xingfan Wu
- Institute of Life Sciences, College of Basic Medicine, Chongqing Medical University, Chongqing, China
| | - Xin Dai
- Institute of Life Sciences, College of Basic Medicine, Chongqing Medical University, Chongqing, China
| | - Xi Li
- Institute of Life Sciences, College of Basic Medicine, Chongqing Medical University, Chongqing, China
| | - Jun Pan
- Key Laboratory for Biorheological Science and Technology of Ministry of Education, College of Bioengineering, Chongqing University, Chongqing, China
| | - Bin Wang
- Institute of Life Sciences, College of Basic Medicine, Chongqing Medical University, Chongqing, China
| |
Collapse
|
4
|
Verbruggen ASK, McNamara LM. Mechanoregulation may drive osteolysis during bone metastasis: A finite element analysis of the mechanical environment within bone tissue during bone metastasis and osteolytic resorption. J Mech Behav Biomed Mater 2023; 138:105662. [PMID: 36630755 DOI: 10.1016/j.jmbbm.2023.105662] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Revised: 12/22/2022] [Accepted: 01/04/2023] [Indexed: 01/09/2023]
Abstract
Metastatic bone disease occurs in 70-80% of advanced breast cancer patients and bone tissue is accepted to have attractive physical properties that facilitate cancer cell attraction, adhesion, and invasion. Bone cells also facilitate tumour invasion by biochemical signalling and through resorption of the bone matrix (osteolysis), which releases factors that further stimulate tumour cell activity. The evolving mechanical environment during tumour invasion might play an important role in these processes, as the activity of both bone and cancer cells is regulated by mechanical cues. In particular bone loss and altered mineralisation have been reported, yet how these alter the mechanical environment local to bone and tumour cells is unknown. The objective of this study is to quantify changes in the mechanical environment within bone tissue, during bone metastasis and osteolytic resorption, using finite element analysis (FEA) models reconstructed from high-resolution μCT images of metastatic mouse bone. In particular, we quantify time-dependent changes in mechanical stimuli, local to and distant from an invading tumour mass, to investigate putative mechanobiological cues for osteolysis during bone metastasis. We report here that in early metastasis (3 weeks after tumour inoculation), there was a decrease in strain distribution within the proximal femur trabecular and distal cortical bone tissue. These changes in the mechanical environment preceded extensive osteolytic destruction, but coincided with the onset of early osteolysis, cortical thickening and mineralisation of proximal and distal femur bone. We propose that early changes in the mechanical environment within bone tissue may activate resorption by osteoclast cells and thereby contribute to the extensive osteolytic bone loss at later stage (6 weeks) bone metastasis.
Collapse
Affiliation(s)
- Anneke S K Verbruggen
- Mechanobiology and Medical Device Research Group (MMDRG), Biomedical Engineering, College of Science and Engineering, University of Galway, Ireland
| | - Laoise M McNamara
- Mechanobiology and Medical Device Research Group (MMDRG), Biomedical Engineering, College of Science and Engineering, University of Galway, Ireland.
| |
Collapse
|
5
|
Computational modeling for osteogenic potential assessment of physical exercises based on loading-induced mechanobiological environments in cortical bone remodeling. Biomech Model Mechanobiol 2023; 22:281-295. [PMID: 36305993 DOI: 10.1007/s10237-022-01647-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2022] [Accepted: 10/08/2022] [Indexed: 11/02/2022]
Abstract
Osteoporosis and disuse can cause bone loss which reduces the weight-bearing strength of long bones. Physical exercise or mechanical loading prevents bone loss as it promotes bone modeling through osteogenesis, i.e., new bone formation. Several studies have observed distinct bone remodeling responses to physical exercises; nevertheless, the underlying mechanism behind such responses is not well established. Loading-induced pore-pressure and fluid motion act as mechanobiological stimuli to bone cells namely osteocytes which further initiate osteoactivities. The shape of loading waveforms also affects the poromechanical environment of bone. Accordingly, the present study hypothesizes that loading waveforms associated with physiological exercises may expose the bone to different mechanobiological stimuli resulting in distinct bone remodeling. A poromechanical finite element model is developed to compute pore-pressure and interstitial fluid velocity in femoral cortical bone tissue (healthy and osteoporotic) subjected to loading waveforms of three physiological exercises namely walking, running, and jumping. The model also computes the mechanobiological stimulus as a function of fluid velocity. The outcomes indicate that pore-pressure and fluid velocity decrease significantly in osteoporotic bone tissue in comparison with healthy tissue. Jumping and running both improve pore-pressure and fluid velocity in healthy and osteoporotic tissues, whereas running significantly enhances mechanobiological stimulus in both the tissues which indicates a possible explanation for distinct bone remodeling to different physical exercises. The present work also suggests that running may be recommended as a potential biomechanical therapeutic to prevent bone loss. Overall, the present work contributes to the area of orthopedic research to develop effective designs of prophylactic exercises to improve bone health.
Collapse
|
6
|
Ruf P, Orassi V, Fischer H, Steffen C, Duda GN, Heiland M, Kreutzer K, Checa S, Rendenbach C. Towards mechanobiologically optimized mandible reconstruction: CAD/CAM miniplates vs. reconstruction plates for fibula free flap fixation: A finite element study. Front Bioeng Biotechnol 2022; 10:1005022. [PMID: 36466355 PMCID: PMC9712730 DOI: 10.3389/fbioe.2022.1005022] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2022] [Accepted: 10/13/2022] [Indexed: 07/28/2023] Open
Abstract
Due to their advantages in applicability, patient-specific (CAD/CAM) reconstruction plates are increasingly used in fibula free flap mandible reconstruction. In addition, recently, CAD/CAM miniplates, with further advantages in postoperative management, have been introduced. However, biomechanical conditions induced by CAD/CAM systems remain partially unknown. This study aimed to evaluate the primary fixation stability of CAD/CAM fixators. For a patient-specific scenario, the biomechanical conditions induced in a one segmental fibula free flap stabilized using either a CAD/CAM reconstruction plate or CAD/CAM miniplates were determined using finite element analysis. The main output parameters were the strains between intersegmental bone surfaces and stresses in the fixation systems due to different biting scenarios. CAD/CAM miniplates resulted in higher mechanical strains in the mesial interosseous gap, whereas CAD/CAM reconstruction plate fixation resulted in higher strains in the distal interosseous gap. For all investigated fixation systems, stresses in the fixation systems were below the material yield stress and thus material failure would not be expected. While the use of CAD/CAM miniplates resulted in strain values considered adequate to promote bone healing in the mesial interosseous gap, in the distal interosseous gap CAD/CAM reconstruction plate fixation might result in more beneficial tissue straining. A mechanical failure of the fixation systems would not be expected.
Collapse
Affiliation(s)
- Philipp Ruf
- Department of Oral and Maxillofacial Surgery, Charité—Universitätsmedizin, Corporate member of Freie Universität Berlin, Humboldt-Universität zu Berlin and Berlin Institute of Health, Berlin, Germany
- Julius Wolff Institute, Berlin Institute of Health at Charité—Universitätsmedizin Berlin, Berlin, Germany
| | - Vincenzo Orassi
- Julius Wolff Institute, Berlin Institute of Health at Charité—Universitätsmedizin Berlin, Berlin, Germany
| | - Heilwig Fischer
- Department of Oral and Maxillofacial Surgery, Charité—Universitätsmedizin, Corporate member of Freie Universität Berlin, Humboldt-Universität zu Berlin and Berlin Institute of Health, Berlin, Germany
- Julius Wolff Institute, Berlin Institute of Health at Charité—Universitätsmedizin Berlin, Berlin, Germany
| | - Claudius Steffen
- Department of Oral and Maxillofacial Surgery, Charité—Universitätsmedizin, Corporate member of Freie Universität Berlin, Humboldt-Universität zu Berlin and Berlin Institute of Health, Berlin, Germany
| | - Georg N. Duda
- Julius Wolff Institute, Berlin Institute of Health at Charité—Universitätsmedizin Berlin, Berlin, Germany
| | - Max Heiland
- Department of Oral and Maxillofacial Surgery, Charité—Universitätsmedizin, Corporate member of Freie Universität Berlin, Humboldt-Universität zu Berlin and Berlin Institute of Health, Berlin, Germany
| | - Kilian Kreutzer
- Department of Oral and Maxillofacial Surgery, Charité—Universitätsmedizin, Corporate member of Freie Universität Berlin, Humboldt-Universität zu Berlin and Berlin Institute of Health, Berlin, Germany
| | - Sara Checa
- Julius Wolff Institute, Berlin Institute of Health at Charité—Universitätsmedizin Berlin, Berlin, Germany
| | - Carsten Rendenbach
- Department of Oral and Maxillofacial Surgery, Charité—Universitätsmedizin, Corporate member of Freie Universität Berlin, Humboldt-Universität zu Berlin and Berlin Institute of Health, Berlin, Germany
| |
Collapse
|
7
|
Henyš P, Kuchař M, Hájek P, Hammer N. Mechanical metric for skeletal biomechanics derived from spectral analysis of stiffness matrix. Sci Rep 2021; 11:15690. [PMID: 34344907 PMCID: PMC8333423 DOI: 10.1038/s41598-021-94998-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2021] [Accepted: 07/19/2021] [Indexed: 02/07/2023] Open
Abstract
A new metric for the quantitative and qualitative evaluation of bone stiffness is introduced. It is based on the spectral decomposition of stiffness matrix computed with finite element method. The here proposed metric is defined as an amplitude rescaled eigenvalues of stiffness matrix. The metric contains unique information on the principal stiffness of bone and reflects both bone shape and material properties. The metric was compared with anthropometrical measures and was tested for sex sensitivity on pelvis bone. Further, the smallest stiffness of pelvis was computed under a certain loading condition and analyzed with respect to sex and direction. The metric complements anthropometrical measures and provides a unique information about the smallest bone stiffness independent from the loading configuration and can be easily computed by state-of-the-art subject specified finite element algorithms.
Collapse
Affiliation(s)
- Petr Henyš
- grid.6912.c0000000110151740Institute of New Technologies and Applied Informatics, Faculty of Mechatronics, Informatics and Interdisciplinary Studies, Technical University of Liberec, Studentská 1402/2, 461 17 Liberec, Czech Republic
| | - Michal Kuchař
- grid.4491.80000 0004 1937 116XDepartment of Anatomy, Faculty of Medicine in Hradec Králové, Charles University, Šimkova 870, 500 03 Hradec Králové, Czech Republic
| | - Petr Hájek
- grid.4491.80000 0004 1937 116XDepartment of Anatomy, Faculty of Medicine in Hradec Králové, Charles University, Šimkova 870, 500 03 Hradec Králové, Czech Republic
| | - Niels Hammer
- grid.11598.340000 0000 8988 2476Department of Macroscopic and Clinical Anatomy, Medical University of Graz, Auenbruggerpl. 2, 8036 Graz, Austria ,grid.9647.c0000 0004 7669 9786Department of Orthopedic and Trauma Surgery, University of Leipzig, Leipzig, Germany ,grid.461651.10000 0004 0574 2038Fraunhofer Institute for Machine Tools and Forming Technology IWU, Nöthnitzer Straße 44, 01187 Dresden, Germany
| |
Collapse
|
8
|
Idkaidek A, Schwarcz H, Jasiuk I. Modeling of bending and torsional stiffnesses of bone at sub-microscale: Effect of curved mineral lamellae. J Biomech 2021; 123:110531. [PMID: 34051614 DOI: 10.1016/j.jbiomech.2021.110531] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2020] [Revised: 05/09/2021] [Accepted: 05/10/2021] [Indexed: 10/21/2022]
Abstract
Recent transmission electron microscopy images of transverse sections of human cortical bone showed that mineral lamellae (polycrystalline sheets of apatite crystals) form arcuate multi-radius patterns around collagen fibrils. The 3-6 nm thick mineral lamellae are arranged in stacks of 3-20 layers and curve around individual fibrils, few fibrils, and higher numbers of collagen fibrils. We evaluate the effect of these stacked mineral lamellae with various radius of curvature patterns on the elastic bending and torsional responses of bone at the sub-microscale using a finite element method. We find that the curved multi-radius stack patterns increased the bending and torsional stiffnesses by 7% and 23%, respectively, compared to when the stacks of mineral lamellae only encircle individual fibrils for the idealized geometric models considered. This study provides new insights into the structure-property relations for the bone ultrastructure.
Collapse
Affiliation(s)
- Ashraf Idkaidek
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, 1206 West Green Street, Urbana, IL 61801, USA
| | - Henry Schwarcz
- School of Earth, Environment & Society, McMaster University, Hamilton, ON, Canada
| | - Iwona Jasiuk
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, 1206 West Green Street, Urbana, IL 61801, USA.
| |
Collapse
|
9
|
Deng C, Gillette JC, Derrick TR. Finite element analysis of femoral neck strains during stair ascent and descent. Sci Rep 2021; 11:9183. [PMID: 33911099 PMCID: PMC8080611 DOI: 10.1038/s41598-021-87936-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2020] [Accepted: 03/23/2021] [Indexed: 11/30/2022] Open
Abstract
For older population, a better understanding of the hip joint loading environment is needed for the prevention of hip pain, and the reduction of the stress fractures and fall risks. Using the motion analysis and inverse dynamics methods, combined with musculoskeletal modelling, static optimization, and finite element (FE) femur model, the difference of femoral neck strains between stair ascent vs. descent, young vs. older populations was compared. A two-way repeated-measures MANOVA was applied to test the effect of age and stair direction on the femoral neck strains. The strains at the femoral neck cross-section were greater for stair descent than ascent for both age groups (mostly P = 0.001 to 0006) but there was no difference between age groups. In this study, femoral neck strains represented bone loading environment in more direct ways than joint reaction forces/moments or joint contact forces, the risk of hip pain, falls and stress fractures might be greater during stair descent than ascent. Possible preventative methods to reduce these risks should be developed in the future studies.
Collapse
Affiliation(s)
- Chen Deng
- Division of Sport Biomechanics, School of Sport Science, Beijing Sport University, Beijing, 100084, People's Republic of China.
| | - Jason C Gillette
- Department of Kinesiology, Iowa State University, Ames, IA, 50010, USA
| | - Timothy R Derrick
- Department of Kinesiology, Iowa State University, Ames, IA, 50010, USA
| |
Collapse
|
10
|
Brémaud L, Cai X, Brenner R, Grimal Q. Maximum effect of the heterogeneity of tissue mineralization on the effective cortical bone elastic properties. Biomech Model Mechanobiol 2021; 20:1509-1518. [PMID: 33884512 DOI: 10.1007/s10237-021-01459-z] [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: 12/07/2020] [Accepted: 04/07/2021] [Indexed: 11/29/2022]
Abstract
The mineralization level is heterogeneous in cortical bone extracellular matrix as a consequence of remodeling. Models of the effective elastic properties at the millimeter scale have been developed based on idealizations of the vascular pore network and matrix properties. Some popular models do not take into account the heterogeneity of the matrix. However, the errors on the predicted elasticity when the difference in elastic properties between osteonal and interstitial tissues is not modeled have not been quantified. This work provides an estimation of the maximum error. We compare the effective elasticity of a representative volume element (RVE) assuming (1) different elastic properties in osteonal and interstitial tissues vs. (2) average matrix properties. In order to account for the variability of bone microstructure, we use a collection of high resolution images of the pore network to build RVEs. In each RVE we assumed a constant osteonal wall thickness and we artificially varied this thickness between 35 and 140 [Formula: see text]m to create RVEs with different amounts of osteonal tissue. The homogenization problem was solved with a fast Fourier transform (FFT)-based numerical scheme. We found that the error depends on pore volume fraction and varies on average from 1 to [Formula: see text] depending on the assumed diameter of the osteons. The results suggest that matrix heterogeneity may be disregarded in cortical bone models in most practical cases.
Collapse
Affiliation(s)
- Luc Brémaud
- Sorbonne Université, INSERM, CNRS, Laboratoire d'Imagerie Biomédicale, 75006, Paris, France.,Sorbonne Université, CNRS, Institut Jean Le Rond d'Alembert, 75005, Paris, France
| | - Xiran Cai
- Sorbonne Université, INSERM, CNRS, Laboratoire d'Imagerie Biomédicale, 75006, Paris, France.,School of Information Science and Technology, ShanghaiTech University, Pudong District, 201210, Shanghai, China
| | - Renald Brenner
- Sorbonne Université, CNRS, Institut Jean Le Rond d'Alembert, 75005, Paris, France
| | - Quentin Grimal
- Sorbonne Université, INSERM, CNRS, Laboratoire d'Imagerie Biomédicale, 75006, Paris, France.
| |
Collapse
|
11
|
Biomechanical design prognosis of two extramedullary fixation devices for subtrochanteric femur fracture: a finite element study. Med Biol Eng Comput 2021; 59:271-285. [PMID: 33417126 DOI: 10.1007/s11517-020-02306-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: 09/15/2020] [Accepted: 12/26/2020] [Indexed: 10/22/2022]
Abstract
The design rationale of extramedullary fixation for femur fracture has remained a matter of debate in the orthopaedic community. The present work provides a comparative preclinical assessment between two standard fracture fixation techniques: dynamic hip screw (DHS) and proximal femoral locking plate (PFLP), by employing finite element (FE)-based in silico models. The study attempts to evaluate and compare the two implants on following biomechanical behaviours: (1) stress variation on the femur and implant, (2) axial displacement of the fixated femur constructs, (3) postoperative stress shielding and longer term external remodelling of the host bone. We hypothesised that, of the two implants, PFLP has better biomechanical characteristics when used for subtrochanteric femoral fracture (SFF) fixation considering long-term adaptation. A comminuted fracture, simulated as two-part fracture gap of 20 mm, was created in the subtrochanteric region of a femur CAD model. Non-uniform physiological load cases were considered. External bone adaptation was modelled mathematically using stress analysis coupled with a growth model, in which strain energy density (SED) acted as feedback control variable. The computational results predicted lower stress shielding (by ~ 6%) and relatively less cortical thinning beneath the plate for PFLP as compared to DHS. DHS-fixated femur, on the other hand, predicted superior postoperative rigidity. Graphical Abstract FE-based comparative assessment between two extramedullary femur fixation devices-dynamic hip screw (DHS) and proximal femoral locking plate (PFLP).
Collapse
|
12
|
Peyroteo MMA, Belinha J, Natal Jorge RM. A mathematical biomechanical model for bone remodeling integrated with a radial point interpolating meshless method. Comput Biol Med 2020; 129:104170. [PMID: 33352308 DOI: 10.1016/j.compbiomed.2020.104170] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2020] [Revised: 12/03/2020] [Accepted: 12/04/2020] [Indexed: 11/18/2022]
Abstract
Bone remodeling is a highly complex process, in which bone cells interact and regulate bone's apparent density as a response to several external and internal stimuli. In this work, this process is numerically described using a novel 2D biomechanical model. Some of the new features in this model are (i) the mathematical parameters used to determine bone's apparent density and cellular density; (ii) an automatic boundary recognition step to spatially control bone remodeling and (iii) an approach to mimic the mechanical transduction to osteoclasts and osteoblasts. Moreover, this model is combined with a meshless approach - the Radial Point Interpolation Method (RPIM). The use of RPIM is an asset for this application, especially in the construction of the boundary maps. This work studies bone's adaptation to a certain loading regime through bone resorption. The signaling pathways of bone cells are dependent on the level of strain energy density (SED) in bone. So, when SED changes, bone cells' functioning is affected, causing also changes on bone's apparent density. With this model, bone is able to achieve an equilibrium state, optimizing its structure to withstand the applied loads. Results suggest that this model has the potential to provide high quality solutions while being a simpler alternative to more complex bone remodeling models in the literature.
Collapse
Affiliation(s)
- M M A Peyroteo
- INEGI, Institute of Science and Innovation in Mechanical and Industrial Engineering, Rua Dr. Roberto Frias, 400, 4200-465, Porto, Portugal.
| | - J Belinha
- School of Engineering, Polytechnic of Porto (ISEP), Mechanical Engineering Department, Rua Dr. António Bernardino de Almeida, 431, 4200-072, Porto, Portugal.
| | - R M Natal Jorge
- Faculty of Engineering of the University of Porto, Mechanical Engineering Department, FEUP, Rua Dr. Roberto Frias, S/N, 4200-465, Porto, Portugal.
| |
Collapse
|
13
|
Ganesh T, Laughrey LE, Niroobakhsh M, Lara-Castillo N. Multiscale finite element modeling of mechanical strains and fluid flow in osteocyte lacunocanalicular system. Bone 2020; 137:115328. [PMID: 32201360 PMCID: PMC7354216 DOI: 10.1016/j.bone.2020.115328] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/11/2020] [Accepted: 03/17/2020] [Indexed: 12/20/2022]
Abstract
Osteocytes form over 90% of the bone cells and are postulated to be mechanosensors responsible for regulating the function of osteoclasts and osteoblasts in bone modeling and remodeling. Physical activity results in mechanical loading on the bones. Osteocytes are thought to be the main mechanosensory cells in bone. Upon load osteocytes secrete key factors initiating downstream signaling pathways that regulate skeletal metabolism including the Wnt/β-catenin signaling pathway. Osteocytes have dendritic structures and are housed in the lacunae and canaliculi within the bone matrix. Mechanical loading is known to have two primary effects, namely a mechanical strain (membrane disruption by stretching) on the lacunae/cells, and fluid flow, in the form of fluid flow shear stress (FFSS), in the space between the cell membranes and the lacuna-canalicular walls. In response, osteocytes get activated via a process called mechanotransduction in which mechanical signals are transduced to biological responses. The study of mechanotransduction is a complex subject involving principles of engineering mechanics as well as biological signaling pathway studies. Several length scales are involved as the mechanical loading on macro sized bones are converted to strain and FFSS responses at the micro-cellular level. Experimental measurements of strain and FFSS at the cellular level are very difficult and correlating them to specific biological activity makes this a very challenging task. One of the methods commonly adopted is a multi-scale approach that combines biological and mechanical experimentation with in silico numerical modeling of the engineering aspects of the problem. Finite element analysis along with fluid-structure interaction methodologies are used to compute the mechanical strain and FFSS. These types of analyses often involve a multi-length scale approach where models of both the macro bone structure and micro structure at the cellular length scale are used. Imaging modalities play a crucial role in the development of the models and present their own challenges. This paper reviews the efforts of various research groups in addressing this problem and presents the work in our research group. A clear understanding of how mechanical stimuli affect the lacunae and perilacunar tissue strains and shear stresses on the cellular membranes may ultimately lead to a better understanding of the process of osteocyte activation.
Collapse
Affiliation(s)
- Thiagarajan Ganesh
- Department of Civil and Mechanical Engineering, University of Missouri-Kansas City, 350L Flarsheim Hall, 5100 Rockhill Road, Kansas City, MO 64110, United States of America.
| | - Loretta E Laughrey
- Department of Civil and Mechanical Engineering, University of Missouri-Kansas City, 350L Flarsheim Hall, 5100 Rockhill Road, Kansas City, MO 64110, United States of America
| | - Mohammadmehdi Niroobakhsh
- Department of Civil and Mechanical Engineering, University of Missouri-Kansas City, 350L Flarsheim Hall, 5100 Rockhill Road, Kansas City, MO 64110, United States of America
| | - Nuria Lara-Castillo
- Department of Oral and Craniofacial Sciences, School of Dentistry, University of Missouri-Kansas City, 650 E 25th Street, Kansas City, MO 64108, United States of America
| |
Collapse
|
14
|
Wang L, Shi Q, Cai Y, Chen Q, Guo X, Li Z. Mechanical–chemical coupled modeling of bone regeneration within a biodegradable polymer scaffold loaded with VEGF. Biomech Model Mechanobiol 2020; 19:2285-2306. [DOI: 10.1007/s10237-020-01339-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2019] [Accepted: 05/12/2020] [Indexed: 10/24/2022]
|
15
|
Li Z, Betts D, Kuhn G, Schirmer M, Müller R, Ruffoni D. Mechanical regulation of bone formation and resorption around implants in a mouse model of osteopenic bone. J R Soc Interface 2020; 16:20180667. [PMID: 30890053 DOI: 10.1098/rsif.2018.0667] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Although mechanical stimulation is considered a promising approach to accelerate implant integration, our understanding of load-driven bone formation and resorption around implants is still limited. This lack of knowledge may delay the development of effective loading protocols to prevent implant loosening, especially in osteoporosis. In healthy bone, formation and resorption are mechanoregulated processes. In the intricate context of peri-implant bone regeneration, it is not clear whether bone (re)modelling can still be load-driven. Here, we investigated the mechanical control of peri-implant bone (re)modelling with a well-controlled mechanobiological experiment. We applied cyclic mechanical loading after implant insertion in tail vertebrae of oestrogen depleted mice and we monitored peri-implant bone response by in vivo micro-CT. Experimental data were combined with micro-finite element simulations to estimate local tissue strains in (re)modelling locations. We demonstrated that a substantial increase in bone mass around the implant could be obtained by loading the entire bone. This augmentation could be attributed to a large reduction in bone resorption rather than to an increase in bone formation. We also showed that following implantation, mechanical regulation of bone (re)modelling was transiently lost. Our findings should help to clarify the role of mechanical stimulation on the maintenance of peri-implant bone mass.
Collapse
Affiliation(s)
- Zihui Li
- 1 Institute for Biomechanics, ETH Zurich , Zurich , Switzerland
| | - Duncan Betts
- 1 Institute for Biomechanics, ETH Zurich , Zurich , Switzerland
| | - Gisela Kuhn
- 1 Institute for Biomechanics, ETH Zurich , Zurich , Switzerland
| | | | - Ralph Müller
- 1 Institute for Biomechanics, ETH Zurich , Zurich , Switzerland
| | - Davide Ruffoni
- 1 Institute for Biomechanics, ETH Zurich , Zurich , Switzerland.,3 Mechanics of Biological and Bioinspired Materials Laboratory, Department of Aerospace and Mechanical Engineering, University of Liège , Liège , Belgium
| |
Collapse
|
16
|
Zhang B, Guo L, Chen H, Ventikos Y, Narayan RJ, Huang J. Finite element evaluations of the mechanical properties of polycaprolactone/hydroxyapatite scaffolds by direct ink writing: Effects of pore geometry. J Mech Behav Biomed Mater 2020; 104:103665. [DOI: 10.1016/j.jmbbm.2020.103665] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2019] [Revised: 01/20/2020] [Accepted: 01/27/2020] [Indexed: 12/13/2022]
|
17
|
Kalajahi SMH, Nazemi SM, Johnston JD. An exclusion approach for addressing partial volume artifacts with quantititive computed tomography-based finite element modeling of the proximal tibia. Med Eng Phys 2019; 76:95-100. [PMID: 31870545 DOI: 10.1016/j.medengphy.2019.10.013] [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: 10/25/2018] [Revised: 10/14/2019] [Accepted: 10/20/2019] [Indexed: 10/25/2022]
Abstract
INTRODUCTION Quantitative computed tomography based finite element modeling (QCT-FE) has potential to clarify the role of subchondral bone stiffness in osteoarthritis. The limited spatial resolution of clinical QCT systems, however, results in partial volume (PV) artifacts and low contrast between cortical and trabecular bone, which adversely affects the accuracy of QCT-FE models. The objective of this research was to evaluate the agreement between stiffness predictions offered by QCT-FE models of proximal tibial subchondral bone (constructed with and without a new voxel-exclusion algorithm) with experimentally-derived local subchondral bone structural stiffness. METHODS Thirteen proximal tibial compartments were obtained and imaged using QCT. Two types of QCT-FE models were developed: (1) standard model, which employed the standard procedure for QCT-FE modeling; and (2) "voxel exclusion (VE)" model, which addressed PV artifacts by excluding low density voxels during the material mapping stage of construction. We assessed agreement between QCT-FE stiffness estimates (using standard and VE approaches) with experimental stiffness by reporting predicted variance from linear regression and mean bias with 95% Limits of Agreement (LOA). RESULTS The standard and VE models explained 81% and 84% of the variance in experimentally measured stiffness, respectively. The standard model showed a mean bias of -268 N/mm (LOA -1210 to 679 N/mm); the VE model showed a mean bias of +59 N/mm (LOA -762 to 910 N/mm). INTERPRETATION The VE model explained more variance in subchondral bone stiffness with less bias. Our findings indicate that the VE method has potential to improve QCT-FE models of bone affected by PV artifacts.
Collapse
Affiliation(s)
| | - S Majid Nazemi
- Biomechanics Research Unit, GIGA In Silico Medicine, University of Liège, Liège, Belgium
| | - James D Johnston
- Department of Mechanical Engineering, University of Saskatchewan, Saskatoon, Canada.
| |
Collapse
|
18
|
Alcântara ACS, Assis I, Prada D, Mehle K, Schwan S, Costa-Paiva L, Skaf MS, Wrobel LC, Sollero P. Patient-Specific Bone Multiscale Modelling, Fracture Simulation and Risk Analysis-A Survey. MATERIALS (BASEL, SWITZERLAND) 2019; 13:E106. [PMID: 31878356 PMCID: PMC6981613 DOI: 10.3390/ma13010106] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/13/2019] [Revised: 12/16/2019] [Accepted: 12/17/2019] [Indexed: 12/26/2022]
Abstract
This paper provides a starting point for researchers and practitioners from biology, medicine, physics and engineering who can benefit from an up-to-date literature survey on patient-specific bone fracture modelling, simulation and risk analysis. This survey hints at a framework for devising realistic patient-specific bone fracture simulations. This paper has 18 sections: Section 1 presents the main interested parties; Section 2 explains the organzation of the text; Section 3 motivates further work on patient-specific bone fracture simulation; Section 4 motivates this survey; Section 5 concerns the collection of bibliographical references; Section 6 motivates the physico-mathematical approach to bone fracture; Section 7 presents the modelling of bone as a continuum; Section 8 categorizes the surveyed literature into a continuum mechanics framework; Section 9 concerns the computational modelling of bone geometry; Section 10 concerns the estimation of bone mechanical properties; Section 11 concerns the selection of boundary conditions representative of bone trauma; Section 12 concerns bone fracture simulation; Section 13 presents the multiscale structure of bone; Section 14 concerns the multiscale mathematical modelling of bone; Section 15 concerns the experimental validation of bone fracture simulations; Section 16 concerns bone fracture risk assessment. Lastly, glossaries for symbols, acronyms, and physico-mathematical terms are provided.
Collapse
Affiliation(s)
- Amadeus C. S. Alcântara
- Department of Computational Mechanics, School of Mechanical Engineering, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil; (A.C.S.A.); (D.P.)
| | - Israel Assis
- Department of Integrated Systems, School of Mechanical Engineering, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil;
| | - Daniel Prada
- Department of Computational Mechanics, School of Mechanical Engineering, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil; (A.C.S.A.); (D.P.)
| | - Konrad Mehle
- Department of Engineering and Natural Sciences, University of Applied Sciences Merseburg, 06217 Merseburg, Germany;
| | - Stefan Schwan
- Fraunhofer Institute for Microstructure of Materials and Systems IMWS, 06120 Halle/Saale, Germany;
| | - Lúcia Costa-Paiva
- Department of Obstetrics and Gynecology, School of Medical Sciences, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-887, Brazil;
| | - Munir S. Skaf
- Institute of Chemistry and Center for Computing in Engineering and Sciences, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil;
| | - Luiz C. Wrobel
- Institute of Materials and Manufacturing, Brunel University London, Uxbridge UB8 3PH, UK;
- Department of Civil and Environmental Engineering, Pontifical Catholic University of Rio de Janeiro, Rio de Janeiro 22451-900, Brazil
| | - Paulo Sollero
- Department of Computational Mechanics, School of Mechanical Engineering, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil; (A.C.S.A.); (D.P.)
| |
Collapse
|
19
|
BARKAOUI ABDELWAHED, KAHLA RABEBBEN, MERZOUKI TAREK, HAMBLI RIDHA. NUMERICAL SIMULATION OF APPARENT DENSITY EVOLUTION OF TRABECULAR BONE UNDER FATIGUE LOADING: EFFECT OF BONE INITIAL PROPERTIES. J MECH MED BIOL 2019. [DOI: 10.1142/s0219519419500416] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Bone remodeling is a physiological phenomenon coupling resorption and formation processes that are mainly mediated by osteoclasts and osteoblasts, in response to mechanical stimuli transduced by osteocytes to biochemical signals activating the bone multicellular unit. Under normal loading conditions, bone resorption and formation are balanced by a homeostasis process. When bone is subjected to overstress, microdamaging occurs, which induces a modification of the structural integrity and microarchitecture. This has drawn significant attention to the mechanical properties of bone. In this context, the current study has been carried out with the aim of numerically investigating the impact of the mechanical properties on the remodeling process of the trabecular bone under cyclic loading, highlighting the effects of different values of the mineral density and the Young’s modulus. This was performed using a mechanobiological model, coupling mechanical and biological approaches, allowing to numerically simulate the effect of the selected parameters for a 20-year-period of cyclic loading for 2D and 3D models of a human femur head. The current work is an explorative numerical study, and the obtained results revealed the changes in the overall stiffness of the bone according to the mechanical properties.
Collapse
Affiliation(s)
- ABDELWAHED BARKAOUI
- Laboratoire des Energies Renouvelables et, Matériaux Avancés (LERMA), Ecole Supérieure de, l’Ingénierie de l’Energie, Université Internationale, de Rabat, Rocade Rabat-Salé, 11100, Rabat-Sala El Jadida, Maroc
| | - RABEB BEN KAHLA
- Laboratoire de Systèmes et de Mécanique, Appliquée (Lasmap-EPT), Ecole Polytechnique de Tunis, Université de Carthage, 2078 La Marsa, Tunisie
| | - TAREK MERZOUKI
- Laboratoire Ingénierie des Systèmes de Versailles, Université de Versailles St Quentin en Yvelines, 10 avenue de l’Europe, 78140 Velizy, France
| | - RIDHA HAMBLI
- Laboratoire Prisme, EA4229, Ecole d’Ingénieurs, Polytech Orléans, Université d’Orléans, 8 Rue Léonard de Vinci, 45072 Orléans, France
| |
Collapse
|
20
|
Cai X, Follet H, Peralta L, Gardegaront M, Farlay D, Gauthier R, Yu B, Gineyts E, Olivier C, Langer M, Gourrier A, Mitton D, Peyrin F, Grimal Q, Laugier P. Anisotropic elastic properties of human femoral cortical bone and relationships with composition and microstructure in elderly. Acta Biomater 2019; 90:254-266. [PMID: 30922952 DOI: 10.1016/j.actbio.2019.03.043] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2018] [Revised: 03/19/2019] [Accepted: 03/20/2019] [Indexed: 02/08/2023]
Abstract
The strong dependence between cortical bone elasticity at the millimetre-scale (mesoscale) and cortical porosity has been evidenced by previous studies. However, bone is an anisotropic composite material made by mineral, proteins and water assembled in a hierarchical structure. Whether the variations of structural and compositional properties of bone affect the different elastic coefficients at the mesoscale is not clear. Aiming to understand the relationships between bone elastic properties and compositions and microstructure, we applied state-of-the-art experimental modalities to assess these aspects of bone characteristics. All elastic coefficients (stiffness tensor of the transverse isotropic bone material), structure of the vascular pore network, collagen and mineral properties were measured in 52 specimens from the femoral diaphysis of 26 elderly donors. Statistical analyses and micromechanical modeling showed that vascular pore volume fraction and the degree of mineralization of bone are the most important determinants of cortical bone anisotropic mesoscopic elasticity. Though significant correlations were observed between collagen properties and elasticity, their effects in bone mesoscopic elasticity were minor in our data. This work also provides a unique set of data exhibiting a range of variations of compositional and microstructural cortical bone properties in the elderly and gives strong experimental evidence and basis for further development of biomechanical models for human cortical bone. STATEMENT OF SIGNIFICANCE: This study reports the relationships between microstructure, composition and the mesoscale anisotropic elastic properties of human femoral cortical bone in elderly. For the first time, we provide data covering the complete anisotropic elastic tensor, the microstructure of cortical vascular porosity, mineral and collagen characteristics obtained from the same or adjacent samples in each donor. The results revealed that cortical vascular porosity and degree of mineralization of bone are the most important determinants of bone anisotropic stiffness at the mesoscale. The presented data gives strong experimental evidence and basis for further development of biomechanical models for human cortical bone.
Collapse
|
21
|
Ruiz Wills C, Olivares AL, Tassani S, Ceresa M, Zimmer V, González Ballester MA, Del Río LM, Humbert L, Noailly J. 3D patient-specific finite element models of the proximal femur based on DXA towards the classification of fracture and non-fracture cases. Bone 2019; 121:89-99. [PMID: 30611923 DOI: 10.1016/j.bone.2019.01.001] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/30/2018] [Revised: 12/21/2018] [Accepted: 01/01/2019] [Indexed: 11/18/2022]
Abstract
Osteoporotic bone fractures reduce quality of life and drastically increase mortality. Minimally irradiating imaging techniques such as dual-energy X-ray absorptiometry (DXA) allow assessment of bone loss through the use of bone mineral density (BMD) as descriptor. Yet, the accuracy of fracture risk predictions remains limited. Recently, DXA-based 3D modelling algorithms were proposed to analyse the geometry and BMD spatial distribution of the proximal femur. This study hypothesizes that such approaches can benefit from finite element (FE)-based biomechanical analyses to improve fracture risk prediction. One hundred and eleven subjects were included in this study and stratified in two groups: (a) 62 fracture cases, and (b) 49 non-fracture controls. Side fall was simulated using a static peak load that depended on patient mass and height. Local mechanical fields were calculated based on relationships between tissue stiffness and BMD. The area under the curve (AUC) of the receiver operating characteristic method evaluated the ability of calculated biomechanical descriptors to discriminate fracture and control cases. The results showed that the major principal stress was better discriminator (AUC > 0.80) than the volumetric BMD (AUC ≤ 0.70). High discrimination capacity was achieved when the analysis was performed by bone type, zone of fracture and gender/sex (AUC of 0.91 for women, trabecular bone and trochanter area), and outcomes suggested that the trabecular bone is critical for fracture discrimination. In conclusion, 3D FE models derived from DXA scans might significantly improve the prediction of hip fracture risk; providing a new insight for clinicians to use FE simulations in clinical practice for osteoporosis management.
Collapse
Affiliation(s)
| | | | - Simone Tassani
- BCN MedTech, Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Mario Ceresa
- BCN MedTech, Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Veronika Zimmer
- BCN MedTech, Universitat Pompeu Fabra (UPF), Barcelona, Spain; School of Biomedical Engineering & Imaging Sciences, King's College London, London, UK
| | | | | | | | - Jérôme Noailly
- BCN MedTech, Universitat Pompeu Fabra (UPF), Barcelona, Spain
| |
Collapse
|
22
|
Solitro GF, Mainnemare F, Amirouche F, Mehta A. A novel technique with reduced computed tomography exposure to predict vertebral compression fracture: a finite element study based on rat vertebrae. Med Biol Eng Comput 2018; 57:795-805. [PMID: 30402789 DOI: 10.1007/s11517-018-1918-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2017] [Accepted: 10/21/2018] [Indexed: 10/27/2022]
Abstract
Vertebral compression fractures are a significant clinical issue with an annual incidence of approximately 750,000 cases in the USA alone. Mechanical properties of vertebrae are successfully evaluated through finite element (FE) models based on vertebrae CT. However, clinical drawbacks associated to radiation transmission encouraged to explore the possibility to use selected or reduced portions of the vertebra. The objective of our study was to develop a new procedure to predict vertebral compression fracture from sub-volumes. We reconstructed rat vertebras from micro-CT of thoracic and lumbar groups. Each vertebra was partitioned into three sub-volumes of different axial thickness. FE simulating compression tests were performed on each model to evaluate their failure load and stiffness. Using a power function, a high correlation was found for stiffness and strength. The sub-volume with three fifths thickness had a failure load of 180.7 ± 19.2 N for thoracic and of 209.5 ± 27.4 N for the lumbar vertebra. These values were not significantly different from the values found for the entire vertebra (p > 0.05). Based on our findings, failure loads and stiffnesses obtained with reduced CT scans can be successfully used to predict full vertebral failure. This sub-region analysis and power relationship suggests that one can limit radiation exposure to patients when bone characterization is needed. Graphical abstract Estimated mechanical properties in relation to the extent of the computed tomography reconstruction.
Collapse
Affiliation(s)
- Giovanni F Solitro
- Department of Orthopaedics, University of Illinois at Chicago, 835 S. Wolcott Avenue, Room E270, Chicago, IL, 60612, USA.,Department of Orthopaedic Surgery, Louisiana State University Health Science Center of Shreveport, 1501 Kings Hwy, Room 3-317, Shreveport, LA, 71104, USA
| | - Florian Mainnemare
- Department of Mechanical Engineering, ENS Cachan, Université Paris-Saclay, 61 Avenue du Président Wilson, 94235, Cachan, France
| | - Farid Amirouche
- Department of Orthopaedics, University of Illinois at Chicago, 835 S. Wolcott Avenue, Room E270, Chicago, IL, 60612, USA.
| | - Ankit Mehta
- Department of Neurosurgery, University of Illinois at Chicago, 912 S Wood St, Chicago, IL, USA
| |
Collapse
|
23
|
Biomechanical features of six design of the delta external fixator for treating Pilon fracture: a finite element study. Med Biol Eng Comput 2018; 56:1925-1938. [DOI: 10.1007/s11517-018-1830-3] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2017] [Accepted: 04/11/2018] [Indexed: 02/06/2023]
|
24
|
Budyn E, Gaci N, Sanders S, Bensidhoum M, Schmidt E, Cinquin B, Tauc P, Petite H. Human Stem Cell Derived Osteocytes in Bone-on-Chip. ACTA ACUST UNITED AC 2018. [DOI: 10.1557/adv.2018.278] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
|
25
|
Effect of rehabilitation exercise durations on the dynamic bone repair process by coupling polymer scaffold degradation and bone formation. Biomech Model Mechanobiol 2017; 17:763-775. [DOI: 10.1007/s10237-017-0991-6] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2017] [Accepted: 11/24/2017] [Indexed: 10/18/2022]
|
26
|
Szlazak K, Vass V, Hasslinger P, Jaroszewicz J, Dejaco A, Idaszek J, Scheiner S, Hellmich C, Swieszkowski W. X-ray physics-based CT-to-composition conversion applied to a tissue engineering scaffold, enabling multiscale simulation of its elastic behavior. MATERIALS SCIENCE & ENGINEERING. C, MATERIALS FOR BIOLOGICAL APPLICATIONS 2017; 95:389-396. [PMID: 30573263 DOI: 10.1016/j.msec.2017.11.044] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/20/2017] [Revised: 09/05/2017] [Accepted: 11/29/2017] [Indexed: 12/13/2022]
Abstract
Nowadays, the assessment of the mechanical competence of tissue engineering scaffolds based on computer simulations is a well-accepted technology. Typically, such simulations are performed by means of the Finite Element (FE) method, with the underlying structural model being created based on micro-computed tomography (microCT). Here, this analysis modality is applied to a new, ternary composite, consisting of PHBV, i.e. poly(3-hydroxybutyrate-co-3-hydroxyvalerate), PLGA, i.e. poly(lactic-co-glycolide), as well as of TCP, i.e. tricalcium phosphate hydrate. The studied scaffold structure is made up by fibers of this new composite material, manufactured by means of the rapid prototyping method. The data collected from microCT is utilized for adequately defining the mechanical properties of the FE model. In particular, the three-dimensional field of grey values is interpreted in terms of the underlying field of attenuation coefficients, taking into account the photon energy employed in microCT imaging, eventually allowing for calculation of the three-dimensionally distributed, voxel-specific composition of the studied material. For the sake of keeping the FE simulations as efficient as possible, groups of voxels are combined into one finite element; the grey value of the latter is obtained by volume averaging. Employing a two-step micromechanical homogenization scheme, the experimentally accessible stiffness of the three constituents (PHBV, PLGA, and TCP) is then, finite element by finite element, upscaled to the composition-dependent stiffness of the composite material. The plausibility and adequacy of the FE model is demonstrated by simulating the effects of uniaxial compression on the scaffold structure, in terms of resulting stress and strain fields, highlighting the importance of the fiber junctions (as they are the mechanically most stressed regions), and that neglecting the material heterogeneity would lead to a potentially significant underestimation of stresses and strains. Finally, a comparison is made of the employed analysis modality of microCT data with a previously pursued, simplified analysis strategy, highlighting the conceptual superiority of the former, and pointing out the application limits of the latter.
Collapse
Affiliation(s)
- Karol Szlazak
- Faculty of Materials Science and Engineering, Warsaw University of Technology, Warsaw, Poland
| | - Viktoria Vass
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Vienna, Austria
| | - Patricia Hasslinger
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Vienna, Austria
| | - Jakub Jaroszewicz
- Faculty of Materials Science and Engineering, Warsaw University of Technology, Warsaw, Poland
| | - Alexander Dejaco
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Vienna, Austria
| | - Joanna Idaszek
- Faculty of Materials Science and Engineering, Warsaw University of Technology, Warsaw, Poland
| | - Stefan Scheiner
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Vienna, Austria.
| | - Christian Hellmich
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Vienna, Austria
| | - Wojciech Swieszkowski
- Faculty of Materials Science and Engineering, Warsaw University of Technology, Warsaw, Poland
| |
Collapse
|
27
|
A Lower Limb-Pelvis Finite Element Model with 3D Active Muscles. Ann Biomed Eng 2017; 46:86-96. [PMID: 29038943 DOI: 10.1007/s10439-017-1942-1] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2017] [Accepted: 10/07/2017] [Indexed: 10/18/2022]
Abstract
A lower limb-pelvis finite element (FE) model with active three-dimensional (3D) muscles was developed in this study for biomechanical analysis of human body. The model geometry was mainly reconstructed from a male volunteer close to the anthropometry of a 50th percentile Chinese male. Tissue materials and structural features were established based on the literature and new implemented experimental tests. In particular, the muscle was modeled with a combination of truss and hexahedral elements to define its passive and active properties as well as to follow the detailed anatomy structure. Both passive and active properties of the model were validated against the experiments of Post-Mortem Human Surrogate (PMHS) and volunteers, respectively. The model was then used to simulate driver's emergency braking during frontal crashes and investigate Knee-Thigh-Hip (KTH) injury mechanisms and tolerances of the human body. A significant force and bending moment variance was noted for the driver's femur due to the effects of active muscle forces during emergency braking. In summary, the present lower limb-pelvis model can be applied in various research fields to support expensive and complex physical tests or corresponding device design.
Collapse
|
28
|
Assessment of trabecular bone tissue elasticity with resonant ultrasound spectroscopy. J Mech Behav Biomed Mater 2017; 74:106-110. [DOI: 10.1016/j.jmbbm.2017.05.037] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2017] [Revised: 05/19/2017] [Accepted: 05/29/2017] [Indexed: 11/17/2022]
|
29
|
Tissue mineral density measured at the sub-millimetre scale can provide reliable statistics of elastic properties of bone matrix. Biomech Model Mechanobiol 2017; 16:1885-1910. [DOI: 10.1007/s10237-017-0926-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2016] [Accepted: 06/08/2017] [Indexed: 12/12/2022]
|
30
|
Mao R, Guo J, Luo C, Fan Y, Wen J, Wang L. Biomechanical study on surgical fixation methods for minimally invasive treatment of hallux valgus. Med Eng Phys 2017; 46:21-26. [PMID: 28527835 DOI: 10.1016/j.medengphy.2017.04.010] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2017] [Revised: 04/18/2017] [Accepted: 04/26/2017] [Indexed: 02/03/2023]
Abstract
Hallux valgus (HV) was one of the most frequent female foot deformities. The aim of this study was to evaluate mechanical responses and stabilities of the Kirschner, bandage and fiberglass fixations after the distal metatarsal osteotomy in HV treatment. Surface traction of different forefoot regions in bandage fixation and the biomechanical behavior of fiberglass bandage material were measured by a pressure sensor device and a mechanical testing, respectively. A three-dimensional foot finite element (FE) model was developed to simulate the three fixation methods (Kirschner, bandage and fiberglass fixations) in weight bearing. The model included 28 bones, sesamoids, ligaments, plantar fascia, cartilages and soft tissue. The peak Von Mises stress (MS) and compression stress (CS) of the distal fragment were predicted from the three fixation methods: Kirschner fixation (MS=6.71MPa, CS=1.232MPa); Bandage fixation (MS=14.90MPa, CS=9.642MPa); Fiberglass fixation (MS=15.83MPa, CS=19.70MPa). Compared with the Kirschner and bandage fixation, the fiberglass fixation reduced the relative movement of osteotomy fragments and obtained the maximum CS. We concluded that fiberglass fixation in HV treatment was helpful to the bone healing of distal fragment. The findings were expected to guide further therapeutic planning of HV patient.
Collapse
Affiliation(s)
- Rui Mao
- Key Laboratory for Biomechanics and Mechanobiology of Ministry of Education, School of Biological Science and Medical Engineering, Beihang University, 100191 Beijing, China
| | - Junchao Guo
- Key Laboratory of Human Motion Analysis and Rehabilitation Technology of the Ministry of Civil Affairs, National Research Centre for Rehabilitation Technical Aids, 100176 Beijing, China
| | - Chenyu Luo
- Key Laboratory for Biomechanics and Mechanobiology of Ministry of Education, School of Biological Science and Medical Engineering, Beihang University, 100191 Beijing, China
| | - Yubo Fan
- Key Laboratory for Biomechanics and Mechanobiology of Ministry of Education, School of Biological Science and Medical Engineering, Beihang University, 100191 Beijing, China; Key Laboratory of Human Motion Analysis and Rehabilitation Technology of the Ministry of Civil Affairs, National Research Centre for Rehabilitation Technical Aids, 100176 Beijing, China
| | - Jianmin Wen
- Wangjing Hospital, China Academy of Chinese Medical Sciences, Beijing, China.
| | - Lizhen Wang
- Key Laboratory for Biomechanics and Mechanobiology of Ministry of Education, School of Biological Science and Medical Engineering, Beihang University, 100191 Beijing, China.
| |
Collapse
|
31
|
Effect of including damage at the tissue level in the nonlinear homogenisation of trabecular bone. Biomech Model Mechanobiol 2017; 16:1681-1695. [PMID: 28500359 PMCID: PMC5599493 DOI: 10.1007/s10237-017-0913-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2016] [Accepted: 04/21/2017] [Indexed: 02/05/2023]
Abstract
Being able to predict bone fracture or implant stability needs a proper constitutive model of trabecular bone at the macroscale in multiaxial, non-monotonic loading modes. Its macroscopic damage behaviour has been investigated experimentally in the past, mostly with the restriction of uniaxial cyclic loading experiments for different samples, which does not allow for the investigation of several load cases in the same sample as damage in one direction may affect the behaviour in other directions. Homogenised finite element models of whole bones have the potential to assess complicated scenarios and thus improve clinical predictions. The aim of this study is to use a homogenisation-based multiscale procedure to upscale the damage behaviour of bone from an assumed solid phase constitutive law and investigate its multiaxial behaviour for the first time. Twelve cubic specimens were each submitted to nine proportional strain histories by using a parallel code developed in-house. Evolution of post-elastic properties for trabecular bone was assessed for a small range of macroscopic plastic strains in these nine load cases. Damage evolution was found to be non-isotropic, and both damage and hardening were found to depend on the loading mode (tensile, compression or shear); both were characterised by linear laws with relatively high coefficients of determination. It is expected that the knowledge of the macroscopic behaviour of trabecular bone gained in this study will help in creating more precise continuum FE models of whole bones that improve clinical predictions.
Collapse
|
32
|
Musy SN, Maquer G, Panyasantisuk J, Wandel J, Zysset PK. Not only stiffness, but also yield strength of the trabecular structure determined by non-linear µFE is best predicted by bone volume fraction and fabric tensor. J Mech Behav Biomed Mater 2017; 65:808-813. [DOI: 10.1016/j.jmbbm.2016.10.004] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2016] [Revised: 09/20/2016] [Accepted: 10/13/2016] [Indexed: 12/11/2022]
|
33
|
Levrero-Florencio F, Manda K, Margetts L, Pankaj P. Nonlinear homogenisation of trabecular bone: Effect of solid phase constitutive model. Proc Inst Mech Eng H 2016; 231:405-414. [PMID: 28427317 DOI: 10.1177/0954411916676220] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Micro-finite element models have been extensively employed to evaluate the elastic properties of trabecular bone and, to a limited extent, its yield behaviour. The macroscopic stiffness tensor and yield surface are of special interest since they are essential in the prediction of bone strength and stability of implants at the whole bone level. While macroscopic elastic properties are now well understood, yield and post-yield properties are not. The aim of this study is to shed some light on what the effect of the solid phase yield criterion is on the macroscopic yield of trabecular bone for samples with different microstructure. Three samples with very different density were subjected to a large set of apparent load cases (which is important since physiological loading is complex and can have multiple components in stress or strain space) with two different solid phase yield criteria: Drucker-Prager and eccentric-ellipsoid. The study found that these two criteria led to small differences in the macroscopic yield strains for most load cases except for those that were compression-dominated; in these load cases, the yield strains for the Drucker-Prager criterion were significantly higher. Higher density samples resulted in higher differences between the two criteria. This work provides a comprehensive assessment of the effect of two different solid phase yield criteria on the macroscopic yield strains of trabecular bone, for a wide range of load cases, and for samples with different morphology.
Collapse
Affiliation(s)
| | - Krishnagoud Manda
- 1 Institute for Bioengineering, School of Engineering, The University of Edinburgh, Edinburgh, UK
| | - Lee Margetts
- 2 School of Mechanical, Aerospace and Civil Engineering The University of Manchester Sackville Street Manchester, UK
| | - Pankaj Pankaj
- 1 Institute for Bioengineering, School of Engineering, The University of Edinburgh, Edinburgh, UK
| |
Collapse
|
34
|
Blanchard R, Morin C, Malandrino A, Vella A, Sant Z, Hellmich C. Patient-specific fracture risk assessment of vertebrae: A multiscale approach coupling X-ray physics and continuum micromechanics. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2016; 32:e02760. [PMID: 26666734 DOI: 10.1002/cnm.2760] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/12/2015] [Accepted: 10/14/2015] [Indexed: 06/05/2023]
Abstract
While in clinical settings, bone mineral density measured by computed tomography (CT) remains the key indicator for bone fracture risk, there is an ongoing quest for more engineering mechanics-based approaches for safety analyses of the skeleton. This calls for determination of suitable material properties from respective CT data, where the traditional approach consists of regression analyses between attenuation-related grey values and mechanical properties. We here present a physics-oriented approach, considering that elasticity and strength of bone tissue originate from the material microstructure and the mechanical properties of its elementary components. Firstly, we reconstruct the linear relation between the clinically accessible grey values making up a CT, and the X-ray attenuation coefficients quantifying the intensity losses from which the image is actually reconstructed. Therefore, we combine X-ray attenuation averaging at different length scales and over different tissues, with recently identified 'universal' composition characteristics of the latter. This gives access to both the normally non-disclosed X-ray energy employed in the CT-device and to in vivo patient-specific and location-specific bone composition variables, such as voxel-specific mass density, as well as collagen and mineral contents. The latter feed an experimentally validated multiscale elastoplastic model based on the hierarchical organization of bone. Corresponding elasticity maps across the organ enter a finite element simulation of a typical load case, and the resulting stress states are increased in a proportional fashion, so as to check the safety against ultimate material failure. In the young patient investigated, even normal physiological loading is probable to already imply plastic events associated with the hydrated mineral crystals in the bone ultrastructure, while the safety factor against failure is still as high as five. Copyright © 2016 John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
- Romane Blanchard
- TU Wien-Vienna University of Technology, Institute for Mechanics of Materials and Structures, Karlsplatz 13/202, Vienna 1040, Austria
| | - Claire Morin
- CIS-EMSE, CNRS:UMR 5307, LGF, Ecole Nationale Supérieure des Mines, Saint-Etienne, F-42023, France
| | - Andrea Malandrino
- Institute for Bioengineering of Catalonia, C/Baldiri Reixac 10-12, Barcelona 08028, Spain
| | - Alain Vella
- Mechanical Engineering Department, University of Malta, Tal Qroqq, Msida MSD, 2080, Malta
| | - Zdenka Sant
- Mechanical Engineering Department, University of Malta, Tal Qroqq, Msida MSD, 2080, Malta
| | - Christian Hellmich
- TU Wien-Vienna University of Technology, Institute for Mechanics of Materials and Structures, Karlsplatz 13/202, Vienna 1040, Austria
| |
Collapse
|
35
|
Le Pense S, Chen Y. Contribution of fluid in bone extravascular matrix to strain-rate dependent stiffening of bone tissue - A poroelastic study. J Mech Behav Biomed Mater 2016; 65:90-101. [PMID: 27569757 DOI: 10.1016/j.jmbbm.2016.08.016] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2016] [Revised: 07/22/2016] [Accepted: 08/06/2016] [Indexed: 11/17/2022]
Abstract
Osteoporotic fractures represent an increasing cost to society, and its diagnosis methods based on bone density still lack accuracy in identifying risk of fracture. This is why a better understanding of mechanical behavior of bone tissue is of importance, especially when it comes to relating experimental observations to realistic physiological fall loading conditions. This study aims at exploring the stiffening effect of pore fluid in bone extravascular matrix subject to high strain rate loading that is more realistic to simulate a physiological fall. A computational approach is used, where bone tissue microstructure extracted from micro-CT images is modeled using finite elements. The solid phase of bone tissue is modeled as a poroelastic material, a porous matrix filled with fluid. When the extravascular matrix experiences certain volumetric deformation, the fluid in pores presents load carrying capacity, which consequently varies the apparent stiffness of bone tissue. It is shown that effects of fluid stiffening in bone can be significant, depending on the chosen material properties, the amount of volumetric strain in tissue and the loading rate with respect to hydraulic conductivity and drainage conditions. It is also shown that such stiffening effect is influenced by bone microstructure, and is more significant in cortical bone than in trabecular bone.
Collapse
Affiliation(s)
- Solenn Le Pense
- Institute of Mechanical, Process and Energy Engineering, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh EH14 4AS, United Kingdom
| | - Yuhang Chen
- Institute of Mechanical, Process and Energy Engineering, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh EH14 4AS, United Kingdom.
| |
Collapse
|
36
|
Evaluating the macroscopic yield behaviour of trabecular bone using a nonlinear homogenisation approach. J Mech Behav Biomed Mater 2016; 61:384-396. [DOI: 10.1016/j.jmbbm.2016.04.008] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2015] [Revised: 03/28/2016] [Accepted: 04/06/2016] [Indexed: 02/07/2023]
|
37
|
Huang H, Xiang C, Zeng C, Ouyang H, Wong KKL, Huang W. Patient-specific geometrical modeling of orthopedic structures with high efficiency and accuracy for finite element modeling and 3D printing. AUSTRALASIAN PHYSICAL & ENGINEERING SCIENCES IN MEDICINE 2015; 38:743-53. [PMID: 26577713 DOI: 10.1007/s13246-015-0402-1] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/17/2014] [Accepted: 11/05/2015] [Indexed: 11/30/2022]
Abstract
We improved the geometrical modeling procedure for fast and accurate reconstruction of orthopedic structures. This procedure consists of medical image segmentation, three-dimensional geometrical reconstruction, and assignment of material properties. The patient-specific orthopedic structures reconstructed by this improved procedure can be used in the virtual surgical planning, 3D printing of real orthopedic structures and finite element analysis. A conventional modeling consists of: image segmentation, geometrical reconstruction, mesh generation, and assignment of material properties. The present study modified the conventional method to enhance software operating procedures. Patient's CT images of different bones were acquired and subsequently reconstructed to give models. The reconstruction procedures were three-dimensional image segmentation, modification of the edge length and quantity of meshes, and the assignment of material properties according to the intensity of gravy value. We compared the performance of our procedures to the conventional procedures modeling in terms of software operating time, success rate and mesh quality. Our proposed framework has the following improvements in the geometrical modeling: (1) processing time: (femur: 87.16 ± 5.90 %; pelvis: 80.16 ± 7.67 %; thoracic vertebra: 17.81 ± 4.36 %; P < 0.05); (2) least volume reduction (femur: 0.26 ± 0.06 %; pelvis: 0.70 ± 0.47, thoracic vertebra: 3.70 ± 1.75 %; P < 0.01) and (3) mesh quality in terms of aspect ratio (femur: 8.00 ± 7.38 %; pelvis: 17.70 ± 9.82 %; thoracic vertebra: 13.93 ± 9.79 %; P < 0.05) and maximum angle (femur: 4.90 ± 5.28 %; pelvis: 17.20 ± 19.29 %; thoracic vertebra: 3.86 ± 3.82 %; P < 0.05). Our proposed patient-specific geometrical modeling requires less operating time and workload, but the orthopedic structures were generated at a higher rate of success as compared with the conventional method. It is expected to benefit the surgical planning of orthopedic structures with less operating time and high accuracy of modeling.
Collapse
Affiliation(s)
- Huajun Huang
- Department of Orthopedics, The Third Affiliated Hospital of Southern Medical University (Academy of Orthopedics · Guangdong Province), Guangzhou, 510630, China.
| | - Chunling Xiang
- Department of Orthopedics and Traumatology, Nanfang Hospital, Southern Medical University, Guangzhou, 510515, China
| | - Canjun Zeng
- Department of Orthopedics, The Third Affiliated Hospital of Southern Medical University (Academy of Orthopedics · Guangdong Province), Guangzhou, 510630, China.,Department of Anatomy, Guangdong Provincial Key Laboratory of Medical Biomechanics, School of Basic Medical Sciences, Southern Medical University, North-1838, Guangzhou, 510515, China
| | - Hanbin Ouyang
- Department of Anatomy, Guangdong Provincial Key Laboratory of Medical Biomechanics, School of Basic Medical Sciences, Southern Medical University, North-1838, Guangzhou, 510515, China
| | - Kelvin Kian Loong Wong
- Engineering Computational Biology, School of Computer Science and Software Engineering, The University of Western Australia, 35 Stirling Highway, Crawley, WA, 6000, Australia
| | - Wenhua Huang
- Department of Anatomy, Guangdong Provincial Key Laboratory of Medical Biomechanics, School of Basic Medical Sciences, Southern Medical University, North-1838, Guangzhou, 510515, China.
| |
Collapse
|
38
|
Geraldes DM, Modenese L, Phillips ATM. Consideration of multiple load cases is critical in modelling orthotropic bone adaptation in the femur. Biomech Model Mechanobiol 2015; 15:1029-42. [PMID: 26578078 PMCID: PMC5021760 DOI: 10.1007/s10237-015-0740-7] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2015] [Accepted: 10/19/2015] [Indexed: 11/02/2022]
Abstract
Functional adaptation of the femur has been investigated in several studies by embedding bone remodelling algorithms in finite element (FE) models, with simplifications often made to the representation of bone's material symmetry and mechanical environment. An orthotropic strain-driven adaptation algorithm is proposed in order to predict the femur's volumetric material property distribution and directionality of its internal structures within a continuum. The algorithm was applied to a FE model of the femur, with muscles, ligaments and joints included explicitly. Multiple load cases representing distinct frames of two activities of daily living (walking and stair climbing) were considered. It is hypothesised that low shear moduli occur in areas of bone that are simply loaded and high shear moduli in areas subjected to complex loading conditions. In addition, it is investigated whether material properties of different femoral regions are stimulated by different activities. The loading and boundary conditions were considered to provide a physiological mechanical environment. The resulting volumetric material property distribution and directionalities agreed with ex vivo imaging data for the whole femur. Regions where non-orthogonal trabecular crossing has been documented coincided with higher values of predicted shear moduli. The topological influence of the different activities modelled was analysed. The influence of stair climbing on the properties of the femoral neck region is highlighted. It is recommended that multiple load cases should be considered when modelling bone adaptation. The orthotropic model of the complete femur is released with this study.
Collapse
Affiliation(s)
- Diogo M Geraldes
- Structural Biomechanics, Department of Civil and Environmental Engineering, Skempton Building, Imperial College London, London, UK. .,Biomechanics Group, Department of Mechanical Engineering, City and Guilds Building, Imperial College London, London, UK.
| | - Luca Modenese
- Structural Biomechanics, Department of Civil and Environmental Engineering, Skempton Building, Imperial College London, London, UK.,Department of Mechanical Engineering, Sir Frederick Mappin Building, Mappin Street, The University of Sheffield, Sheffield, UK.,INSIGNEO Institute for In Silico Medicine, The Pam Liversidge Building, The University of Sheffield, Sheffield, UK
| | - Andrew T M Phillips
- Structural Biomechanics, Department of Civil and Environmental Engineering, Skempton Building, Imperial College London, London, UK
| |
Collapse
|
39
|
Stochastic multiscale modelling of cortical bone elasticity based on high-resolution imaging. Biomech Model Mechanobiol 2015. [DOI: 10.1007/s10237-015-0695-8] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
40
|
Lu Y, Engelke K, Glueer CC, Morlock MM, Huber G. The effect of in situ/in vitro three-dimensional quantitative computed tomography image voxel size on the finite element model of human vertebral cancellous bone. Proc Inst Mech Eng H 2014; 228:1208-13. [DOI: 10.1177/0954411914558654] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Quantitative computed tomography–based finite element modeling technique is a promising clinical tool for the prediction of bone strength. However, quantitative computed tomography–based finite element models were created from image datasets with different image voxel sizes. The aim of this study was to investigate whether there is an influence of image voxel size on the finite element models. In all 12 thoracolumbar vertebrae were scanned prior to autopsy (in situ) using two different quantitative computed tomography scan protocols, which resulted in image datasets with two different voxel sizes (0.29 × 0.29 × 1.3 mm3 vs 0.18 × 0.18 × 0.6 mm3). Eight of them were scanned after autopsy (in vitro) and the datasets were reconstructed with two voxel sizes (0.32 × 0.32 × 0.6 mm3 vs. 0.18 × 0.18 × 0.3 mm3). Finite element models with cuboid volume of interest extracted from the vertebral cancellous part were created and inhomogeneous bilinear bone properties were defined. Axial compression was simulated. No effect of voxel size was detected on the apparent bone mineral density for both the in situ and in vitro cases. However, the apparent modulus and yield strength showed significant differences in the two voxel size group pairs (in situ and in vitro). In conclusion, the image voxel size may have to be considered when the finite element voxel modeling technique is used in clinical applications.
Collapse
Affiliation(s)
- Yongtao Lu
- Institute of Biomechanics, Hamburg University of Technology (TUHH), Hamburg, Germany
- INSIGNEO Institute for in silico Medicine, Department of Mechanical Engineering, The University of Sheffield, Sheffield, UK
| | - Klaus Engelke
- Institute of Medical Physics, University of Erlangen-Nuremberg, Erlangen, Germany
| | - Claus-C Glueer
- Section Biomedical Imaging, Department of Radiology and Neuroradiology, Christian-Albrechts Universität zu Kiel, Kiel, Germany
| | - Michael M Morlock
- Institute of Biomechanics, Hamburg University of Technology (TUHH), Hamburg, Germany
| | - Gerd Huber
- Institute of Biomechanics, Hamburg University of Technology (TUHH), Hamburg, Germany
| |
Collapse
|
41
|
Abdalrahman T, Scheiner S, Hellmich C. Is trabecular bone permeability governed by molecular ordering-induced fluid viscosity gain? Arguments from re-evaluation of experimental data in the framework of homogenization theory. J Theor Biol 2014; 365:433-44. [PMID: 25452137 DOI: 10.1016/j.jtbi.2014.10.011] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2014] [Revised: 07/28/2014] [Accepted: 10/09/2014] [Indexed: 12/29/2022]
Abstract
It is generally agreed on that trabecular bone permeability, a physiologically important quantity, is governed by the material׳s (vascular or intertrabecular) porosity as well as by the viscosity of the pore-filling fluids. Still, there is less agreement on how these two key factors govern bone permeability. In order to shed more light onto this somewhat open issue, we here develop a random homogenization scheme for upscaling Poiseuille flow in the vascular porosity, up to Darcy-type permeability of the overall porous medium "trabecular bone". The underlying representative volume element of the macroscopic bone material contains two types of phases: a spherical, impermeable extracellular bone matrix phase interacts with interpenetrating cylindrical pore channel phases that are oriented in all different space directions. This type of interaction is modeled by means of a self-consistent homogenization scheme. While the permeability of the bone matrix equals to zero, the permeability of the pore phase is found through expressing the classical Hagen-Poiseuille law for laminar flow in the format of a "micro-Darcy law". The upscaling scheme contains pore size and porosity as geometrical input variables; however, they can be related to each other, based on well-known relations between porosity and specific bone surface. As two key results, validated through comprehensive experimental data, it appears (i) that the famous Kozeny-Carman constant (which relates bone permeability to the cube of the porosity, the square of the specific surface, as well as to the bone fluid viscosity) needs to be replaced by an again porosity-dependent rational function, and (ii) that the overall bone permeability is strongly affected by the pore fluid viscosity, which, in case of polarized fluids, is strongly increased due to the presence of electrically charged pore walls.
Collapse
Affiliation(s)
- T Abdalrahman
- Institute for Mechanics of Materials and Structures, Vienna University of Technology (TU Wien), 1040 Vienna, Austria.
| | - S Scheiner
- Institute for Mechanics of Materials and Structures, Vienna University of Technology (TU Wien), 1040 Vienna, Austria.
| | - C Hellmich
- Institute for Mechanics of Materials and Structures, Vienna University of Technology (TU Wien), 1040 Vienna, Austria.
| |
Collapse
|
42
|
Bone remodelling in humans is load-driven but not lazy. Nat Commun 2014; 5:4855. [DOI: 10.1038/ncomms5855] [Citation(s) in RCA: 162] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2014] [Accepted: 07/30/2014] [Indexed: 01/20/2023] Open
|
43
|
Geraldes DM, Phillips ATM. A comparative study of orthotropic and isotropic bone adaptation in the femur. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2014; 30:873-889. [PMID: 24753477 PMCID: PMC4272570 DOI: 10.1002/cnm.2633] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/10/2013] [Revised: 01/27/2014] [Accepted: 02/11/2014] [Indexed: 06/03/2023]
Abstract
Functional adaptation of the femur has been studied extensively by embedding remodelling algorithms in finite element models, with bone commonly assumed to have isotropic material properties for computational efficiency. However, isotropy is insufficient in predicting the directionality of bone's observed microstructure. A novel iterative orthotropic 3D adaptation algorithm is proposed and applied to a finite element model of the whole femur. Bone was modelled as an optimised strain-driven adaptive continuum with local orthotropic symmetry. Each element's material orientations were aligned with the local principal stress directions and their corresponding directional Young's moduli updated proportionally to the associated strain stimuli. The converged predicted density distributions for a coronal section of the whole femur were qualitatively and quantitatively compared with the results obtained by the commonly used isotropic approach to bone adaptation and with ex vivo imaging data. The orthotropic assumption was shown to improve the prediction of bone density distribution when compared with the more commonly used isotropic approach, whilst producing lower comparative mass, structurally optimised models. It was also shown that the orthotropic approach can provide additional directional information on the material properties distributions for the whole femur, an advantage over isotropic bone adaptation. Orthotropic bone models can help in improving research areas in biomechanics where local structure and mechanical properties are of key importance, such as fracture prediction and implant assessment.
Collapse
Affiliation(s)
- Diogo M Geraldes
- Structural Biomechanics, Department of Civil and Environmental Engineering, Skempton Building, Imperial College London, London SW7 2AZ, UK
| | | |
Collapse
|
44
|
Malandrino A, Lacroix D, Hellmich C, Ito K, Ferguson SJ, Noailly J. The role of endplate poromechanical properties on the nutrient availability in the intervertebral disc. Osteoarthritis Cartilage 2014; 22:1053-60. [PMID: 24857972 DOI: 10.1016/j.joca.2014.05.005] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/15/2013] [Revised: 04/15/2014] [Accepted: 05/07/2014] [Indexed: 02/02/2023]
Abstract
OBJECTIVE To investigate the relevance of the human vertebral endplate poromechanics on the fluid and metabolic transport from and to the intervertebral disc (IVD) based on educated estimations of the poromechanical parameter values of the bony endplate (BEP). METHODS 50 micro-models of different BEP samples were generated from μCTs of lumbar vertebrae and allowed direct determination of porosity values. Permeability values were calculated by using the micro-models, through the simulation of permeation via computational fluid dynamics. These educated ranges of porosity and permeability values were used as inputs for mechano-transport simulations to assess their effect on both the distributions of metabolites within an IVD model and the poromechanical calculations within the cartilaginous part of the endplate i.e., the cartilage endplate (CEP). RESULTS BEP effective permeability was highly correlated to local variations of porosity (R(2) ≈ 0.88). Universal patterns between bone volume fraction and permeability arose from these results and from other experimental data in the literature. These variations in BEP permeability and porosity had negligible effects on the distributions of metabolites within the disc. In the CEP, the variability of the poromechanical properties of the BEP did not affect the predicted consolidation but induced higher fluid velocities. CONCLUSIONS The present paper provides the first sets of thoroughly identified BEP parameter values that can be further used in patient-specific poromechanical studies. Representing BEP structural changes through variations in poromechanical properties did not affect the diffusion of metabolites. However, attention might be paid to alterations in fluid velocities and cell mechano-sensing within the CEP.
Collapse
Affiliation(s)
- A Malandrino
- Biomechanics and Mechanobiology, Institute for Bioengineering of Catalonia, Barcelona, Spain
| | - D Lacroix
- INSIGNEO Institute for in silico Medicine, Department of Mechanical Engineering, University of Sheffield, Sheffield, UK
| | - C Hellmich
- Institute for Mechanics of Materials and Structures, Vienna University of Technology, Vienna, Austria
| | - K Ito
- Orthopedic Biomechanics, Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - S J Ferguson
- Institute for Biomechanics, ETH, Zurich, Switzerland
| | - J Noailly
- Biomechanics and Mechanobiology, Institute for Bioengineering of Catalonia, Barcelona, Spain.
| |
Collapse
|
45
|
Ascenzi MG, Lutz A, Du X, Klimecky L, Kawas N, Hourany T, Jahng J, Chin J, Tintut Y, Nackenhors U, Keyak J. Hyperlipidemia affects multiscale structure and strength of murine femur. J Biomech 2014; 47:2436-43. [PMID: 24795172 DOI: 10.1016/j.jbiomech.2014.04.006] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2013] [Revised: 04/07/2014] [Accepted: 04/07/2014] [Indexed: 12/13/2022]
Abstract
To improve bone strength prediction beyond limitations of assessment founded solely on the bone mineral component, we investigated the effect of hyperlipidemia, present in more than 40% of osteoporotic patients, on multiscale structure of murine bone. Our overarching purpose is to estimate bone strength accurately, to facilitate mitigating fracture morbidity and mortality in patients. Because (i) orientation of collagen type I affects, independently of degree of mineralization, cortical bone׳s micro-structural strength; and, (ii) hyperlipidemia affects collagen orientation and μCT volumetric tissue mineral density (vTMD) in murine cortical bone, we have constructed the first multiscale finite element (mFE), mouse-specific femoral model to study the effect of collagen orientation and vTMD on strength in Ldlr(-/-), a mouse model of hyperlipidemia, and its control wild type, on either high fat diet or normal diet. Each µCT scan-based mFE model included either element-specific elastic orthotropic properties calculated from collagen orientation and vTMD (collagen-density model) by experimentally validated formulation, or usual element-specific elastic isotropic material properties dependent on vTMD-only (density-only model). We found that collagen orientation, assessed by circularly polarized light and confocal microscopies, and vTMD, differed among groups and that microindentation results strongly correlate with elastic modulus of collagen-density models (r(2)=0.85, p=10(-5)). Collagen-density models yielded (1) larger strains, and therefore lower strength, in simulations of 3-point bending and physiological loading; and (2) higher correlation between mFE-predicted strength and 3-point bending experimental strength, than density-only models. This novel method supports ongoing translational research to achieve the as yet elusive goal of accurate bone strength prediction.
Collapse
Affiliation(s)
- Maria-Grazia Ascenzi
- Department of Orthopaedic Surgery, University of California, Los Angeles, CA 90095, USA.
| | - Andre Lutz
- Continental Tire Company, Hannover, Germany.
| | - Xia Du
- Department of Orthopaedic Surgery, University of California, Los Angeles, CA 90095, USA.
| | | | - Neal Kawas
- Department of Molecular Oncology, John Wayne Cancer Institute, Santa Monica, CA 90404, USA.
| | - Talia Hourany
- Department of Orthopaedic Surgery, University of California, Los Angeles, CA 90095, USA.
| | - Joelle Jahng
- Department of Orthopaedic Surgery, University of California, Los Angeles, CA 90095, USA.
| | - Jesse Chin
- Department of Orthopaedic Surgery, University of California, Los Angeles, CA 90095, USA.
| | - Yin Tintut
- Department of Medicine, University of California, Los Angeles, CA 90095, USA.
| | - Udo Nackenhors
- Institute of Mechanics and Computational Mechanics, Leibniz University Hannover, 30167 Hannover, Germany.
| | - Joyce Keyak
- Department of Radiological Sciences, Medical Sciences I, Bldg 811, Room B140, University of California, Irvine, CA 92697, USA.
| |
Collapse
|