1
|
McPhee S, Kershaw LE, Daniel CR, Peña Fernández M, Cillán-García E, Taylor SE, Wolfram U. QCT-based computational bone strength assessment updated with MRI-derived 'hidden' microporosity. J Mech Behav Biomed Mater 2023; 147:106094. [PMID: 37741181 DOI: 10.1016/j.jmbbm.2023.106094] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Revised: 08/10/2023] [Accepted: 08/26/2023] [Indexed: 09/25/2023]
Abstract
Microdamage accumulated through sustained periods of cyclic loading or single overloading events contributes to bone fragility through a reduction in stiffness and strength. Monitoring microdamage in vivo remains unattainable by clinical imaging modalities. As such, there are no established computational methods for clinical fracture risk assessment that account for microdamage that exists in vivo at any specific timepoint. We propose a method that combines multiple clinical imaging modalities to identify an indicative surrogate, which we term 'hidden porosity', that incorporates pre-existing bone microdamage in vivo. To do so, we use the third metacarpal bone of the equine athlete as an exemplary model for fatigue induced microdamage, which coalesces in the subchondral bone. N = 10 metacarpals were scanned by clinical quantitative computed tomography (QCT) and magnetic resonance imaging (MRI). We used a patch-based similarity method to quantify the signal intensity of a fluid sensitive MRI sequence in bone regions where microdamage coalesces. The method generated MRI-derived pseudoCT images which were then used to determine a pre-existing damage (Dpex) variable to quantify the proposed surrogate and which we incorporate into a nonlinear constitutive model for bone tissue. The minimum, median, and maximum detected Dpex of 0.059, 0.209, and 0.353 reduced material stiffness by 5.9%, 20.9%, and 35.3% as well as yield stress by 5.9%, 20.3%, and 35.3%. Limb-specific voxel-based finite element meshes were equipped with the updated material model. Lateral and medial condyles of each metacarpal were loaded to simulate physiological joint loading during gallop. The degree of detected Dpex correlated with a relative reduction in both condylar stiffness (p = 0.001, R2 > 0.74) and strength (p < 0.001, R2 > 0.80). Our results illustrate the complementary value of looking beyond clinical CT, which neglects the inclusion of microdamage due to partial volume effects. As we use clinically available imaging techniques, our results may aid research beyond the equine model on fracture risk assessment in human diseases such as osteoarthritis, bone cancer, or osteoporosis.
Collapse
Affiliation(s)
- Samuel McPhee
- School of Engineering and Physical Sciences, Institute of Mechanical, Process and Energy Engineering, Heriot-Watt University, Edinburgh, UK
| | - Lucy E Kershaw
- Centre for Cardiovascular Sciences and Edinburgh Imaging, The University of Edinburgh, Edinburgh, UK
| | - Carola R Daniel
- Royal (Dick) School of Veterinary Studies and the Roslin Institute, The University of Edinburgh, Edinburgh, UK
| | - Marta Peña Fernández
- School of Engineering and Physical Sciences, Institute of Mechanical, Process and Energy Engineering, Heriot-Watt University, Edinburgh, UK
| | | | - Sarah E Taylor
- Royal (Dick) School of Veterinary Studies and the Roslin Institute, The University of Edinburgh, Edinburgh, UK
| | - Uwe Wolfram
- School of Engineering and Physical Sciences, Institute of Mechanical, Process and Energy Engineering, Heriot-Watt University, Edinburgh, UK.
| |
Collapse
|
2
|
Bachmann S, Pahr DH, Synek A. Hip joint load prediction using inverse bone remodeling with homogenized FE models: Comparison to micro-FE and influence of material modeling strategy. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2023; 236:107549. [PMID: 37084528 DOI: 10.1016/j.cmpb.2023.107549] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2022] [Revised: 03/23/2023] [Accepted: 04/12/2023] [Indexed: 05/03/2023]
Abstract
BACKGROUND AND OBJECTIVE Measuring physiological loading conditions in vivo can be challenging, as methods are invasive or pose a high modeling effort. However, the physiological loading of bones is also imprinted in the bone microstructure due to bone (re)modeling. This information can be retrieved by inverse bone remodeling (IBR). Recently, an IBR method based on micro-finite-element (µFE) modeling was translated to homogenized-FE (hFE) to decrease computational effort and tested on the distal radius. However, this bone has a relatively simple geometry and homogeneous microstructure. Therefore, the objective of this study was to assess the agreement of hFE-based IBR with µFE-based IBR to predict hip joint loading from the head of the femur; a bone with more complex loading as well as more heterogeneous microstructure. METHODS hFE-based IBR was applied to a set of 19 femoral heads using four different material mapping laws. One model with a single homogeneous material for both trabecular and cortical volume and three models with a separated cortex and either homogeneous, density-dependent inhomogeneous, or density and fabric-dependent orthotropic material. Three different evaluation regions (full bone, trabecular bone only, head region only) were defined, in which IBR was applied. µFE models were created for the same bones, and the agreement of the predicted hip joint loading history obtained from hFE and µFE models was evaluated. The loading history was discretized using four unit load cases. RESULTS The computational time for FE solving was decreased on average from 500 h to under 1 min (CPU time) when using hFE models instead of µFE models. Using more information in the material model in the hFE models led to a better prediction of hip joint loading history. Inhomogeneous and inhomogeneous orthotropic models gave the best agreement to µFE-based IBR (RMSE% <14%). The evaluation region only played a minor role. CONCLUSIONS hFE-based IBR was able to reconstruct the dominant joint loading of the femoral head in agreement with µFE-based IBR and required considerably lower computational effort. Results indicate that cortical and trabecular bone should be modeled separately and at least density-dependent inhomogeneous material properties should be used with hFE models of the femoral head to predict joint loading.
Collapse
Affiliation(s)
- Sebastian Bachmann
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Gumpendorfer Straße 7, Vienna 1060, Austria.
| | - Dieter H Pahr
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Gumpendorfer Straße 7, Vienna 1060, Austria; Division Biomechanics, Karl Landsteiner University of Health Sciences, Dr. Karl-Dorrek-Straße 30, Krems 3500, Austria
| | - Alexander Synek
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Gumpendorfer Straße 7, Vienna 1060, Austria
| |
Collapse
|
3
|
Prada DM, Galvis AF, Miller J, Foster JM, Zavaglia C. Multiscale stiffness characterisation of both healthy and osteoporotic bone tissue using subject-specific data. J Mech Behav Biomed Mater 2022; 135:105431. [DOI: 10.1016/j.jmbbm.2022.105431] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Revised: 08/22/2022] [Accepted: 08/26/2022] [Indexed: 10/31/2022]
|
4
|
Pisano AA, Fuschi P. Limit analysis of human proximal femur. J Mech Behav Biomed Mater 2021; 124:104844. [PMID: 34601433 DOI: 10.1016/j.jmbbm.2021.104844] [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: 03/16/2021] [Revised: 09/03/2021] [Accepted: 09/15/2021] [Indexed: 10/20/2022]
Abstract
A limit analysis numerical approach oriented to predict the peak/collapse load of human proximal femur, under two different loading conditions, is presented. A yield criterion of Tsai-Hu-type, expressed in principal stress space, is used to model the orthotropic bone tissues. A simplified human femur 3D model is envisaged to carry on numerical simulation of in-vitro tests borrowed from the relevant literature and to reproduce their findings. A critical discussion, together with possible future developments, is presented.
Collapse
Affiliation(s)
- A A Pisano
- University Mediterranea of Reggio Calabria, Via dell'Universitá 25, I-89124 Reggio Calabria, Italy.
| | - P Fuschi
- University Mediterranea of Reggio Calabria, Via dell'Universitá 25, I-89124 Reggio Calabria, Italy
| |
Collapse
|
5
|
Irarrázaval S, Ramos-Grez JA, Pérez LI, Besa P, Ibáñez A. Finite element modeling of multiple density materials of bone specimens for biomechanical behavior evaluation. SN APPLIED SCIENCES 2021. [DOI: 10.1007/s42452-021-04760-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
Abstract
AbstractThe finite elements method allied with the computerized axial tomography (CT) is a mathematical modeling technique that allows constructing computational models for bone specimens from CT data. The objective of this work was to compare the experimental biomechanical behavior by three-point bending tests of porcine femur specimens with different types of computational models generated through the finite elements’ method and a multiple density materials assignation scheme. Using five femur specimens, 25 scenarios were created with differing quantities of materials. This latter was applied to computational models and in bone specimens subjected to failure. Among the three main highlights found, first, the results evidenced high precision in predicting experimental reaction force versus displacement in the models with larger number of assigned materials, with maximal results being an R2 of 0.99 and a minimum root-mean-square error of 3.29%. Secondly, measured and computed elastic stiffness values follow same trend with regard to specimen mass, and the latter underestimates stiffness values a 6% in average. Third and final highlight, this model can precisely and non-invasively assess bone tissue mechanical resistance based on subject-specific CT data, particularly if specimen deformation values at fracture are considered as part of the assessment procedure.
Collapse
|
6
|
Femoral strength can be predicted from 2D projections using a 3D statistical deformation and texture model with finite element analysis. Med Eng Phys 2021; 93:72-82. [PMID: 34154777 DOI: 10.1016/j.medengphy.2021.05.012] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Revised: 04/11/2021] [Accepted: 05/18/2021] [Indexed: 11/22/2022]
Abstract
Ultimate force of the proximal human femur can be predicted using Finite Element Analysis (FEA), but the models rely on 3D computed tomography images. Landmark-based statistical appearance models (SAM) and B-Spline transformation-based statistical deformation models (SDM) have been used to estimate 3D images from 2D projections, which facilitates model generation and reduces the radiation dose. However, there is no literature on the accuracy of SDM-based FEA models of bones with respect to experimental results. In this study, a methodology for an enhanced SDM with textural information is presented. The statistical deformation and texture models (SDTMs) are based on a set of 37 quantitative CT (QCT) images. They were used to estimate 3D images from two or one projections of the set in a leave-one-out setup. These estimations where then used to create FEA models. The ultimate force predicted by FEA models estimated from two or one projection using the SDTMs were compared to the experimental ultimate force from a previous study on the same femora and to the results of standard QCT-based FEA models. High correlations between predictions and experimental measurements were found for FEA models reconstructed from 2D projections with R2=0.835 when based on two projections and R2=0.724 when using one projection. The correlations were comparable to those reached with standard QCT-based FE-models with the experimental results (R2=0.795). This study shows the high potential of SDTM-based 3D image reconstruction and FEA modelling from 2D projections to predict femoral ultimate force.
Collapse
|
7
|
Zhang L, Wang L, Fu R, Wang J, Yang D, Liu Y, Zhang W, Liang W, Yang R, Yang H, Cheng X. In Vivo
Assessment of Age‐ and Loading Configuration‐Related Changes in Multiscale Mechanical Behavior of the Human Proximal Femur Using MRI‐Based Finite Element Analysis. J Magn Reson Imaging 2020; 53:905-912. [PMID: 33075178 DOI: 10.1002/jmri.27403] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2020] [Revised: 10/03/2020] [Accepted: 10/05/2020] [Indexed: 01/08/2023] Open
Affiliation(s)
- Lingyun Zhang
- Department of Biomedical Engineering, Faculty of Environment and Life Science Beijing University of Technology Beijing China
| | - Ling Wang
- Department of Radiology Beijing Jishuitan Hospital Beijing China
| | - Ruisen Fu
- Department of Biomedical Engineering, Faculty of Environment and Life Science Beijing University of Technology Beijing China
| | - Jianing Wang
- Department of Biomedical Engineering, Faculty of Environment and Life Science Beijing University of Technology Beijing China
| | - Dongyue Yang
- Department of Biomedical Engineering, Faculty of Environment and Life Science Beijing University of Technology Beijing China
| | - Yandong Liu
- Department of Radiology Beijing Jishuitan Hospital Beijing China
| | - Wei Zhang
- Department of Radiology Beijing Jishuitan Hospital Beijing China
| | - Wei Liang
- Department of Radiology Beijing Jishuitan Hospital Beijing China
| | - Ruopei Yang
- Department of Radiology Beijing Jishuitan Hospital Beijing China
| | - Haisheng Yang
- Department of Biomedical Engineering, Faculty of Environment and Life Science Beijing University of Technology Beijing China
| | - Xiaoguang Cheng
- Department of Radiology Beijing Jishuitan Hospital Beijing China
| |
Collapse
|
8
|
Bokam P, Caille L, Germaneau A, Rigoard P, Vendeuvre T, Valle V. Identification of fractures in cancellous bone using µFE models and Heaviside-Digital volume correlation. Comput Methods Biomech Biomed Engin 2020. [DOI: 10.1080/10255842.2020.1811507] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Prasanth Bokam
- Institut Pprime. UPR 3346 CNRS, Université de Poitiers - ISAE-ENSMA, France
| | - Laetitia Caille
- Institut Pprime. UPR 3346 CNRS, Université de Poitiers - ISAE-ENSMA, France
| | - Arnaud Germaneau
- Institut Pprime. UPR 3346 CNRS, Université de Poitiers - ISAE-ENSMA, France
| | - Philippe Rigoard
- Institut Pprime. UPR 3346 CNRS, Université de Poitiers - ISAE-ENSMA, France
- Spine & Neuromodulation Function Unit. PRISMATICS Lab CHU. Poitiers, France
| | - Tanguy Vendeuvre
- Institut Pprime. UPR 3346 CNRS, Université de Poitiers - ISAE-ENSMA, France
- Spine & Neuromodulation Function Unit. PRISMATICS Lab CHU. Poitiers, France
| | - Valery Valle
- Institut Pprime. UPR 3346 CNRS, Université de Poitiers - ISAE-ENSMA, France
| |
Collapse
|
9
|
Falcinelli C, Whyne C. Image-based finite-element modeling of the human femur. Comput Methods Biomech Biomed Engin 2020; 23:1138-1161. [PMID: 32657148 DOI: 10.1080/10255842.2020.1789863] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Fracture is considered a critical clinical endpoint in skeletal pathologies including osteoporosis and bone metastases. However, current clinical guidelines are limited with respect to identifying cases at high risk of fracture, as they do not account for many mechanical determinants that contribute to bone fracture. Improving fracture risk assessment is an important area of research with clear clinical relevance. Patient-specific numerical musculoskeletal models generated from diagnostic images are widely used in biomechanics research and may provide the foundation for clinical tools used to quantify fracture risk. However, prior to clinical translation, in vitro validation of predictions generated from such numerical models is necessary. Despite adopting radically different models, in vitro validation of image-based finite element (FE) models of the proximal femur (predicting strains and failure loads) have shown very similar, encouraging levels of accuracy. The accuracy of such in vitro models has motivated their application to clinical studies of osteoporotic and metastatic fractures. Such models have demonstrated promising but heterogeneous results, which may be explained by the lack of a uniform strategy with respect to FE modeling of the human femur. This review aims to critically discuss the state of the art of image-based femoral FE modeling strategies, highlighting principal features and differences among current approaches. Quantitative results are also reported with respect to the level of accuracy achieved from in vitro evaluations and clinical applications and are used to motivate the adoption of a standardized approach/workflow for image-based FE modeling of the femur.
Collapse
Affiliation(s)
- Cristina Falcinelli
- Orthopaedic Biomechanics Laboratory, Sunnybrook Research Institute, Toronto, Canada
| | - Cari Whyne
- Orthopaedic Biomechanics Laboratory, Sunnybrook Research Institute, Toronto, Canada
| |
Collapse
|
10
|
Ovesy M, Aeschlimann M, Zysset PK. Explicit finite element analysis can predict the mechanical response of conical implant press-fit in homogenized trabecular bone. J Biomech 2020; 107:109844. [DOI: 10.1016/j.jbiomech.2020.109844] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2019] [Revised: 04/30/2020] [Accepted: 05/08/2020] [Indexed: 02/09/2023]
|
11
|
Sas A, Tanck E, Sermon A, van Lenthe GH. Finite element models for fracture prevention in patients with metastatic bone disease. A literature review. Bone Rep 2020; 12:100286. [PMID: 32551337 PMCID: PMC7292864 DOI: 10.1016/j.bonr.2020.100286] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/17/2020] [Revised: 05/04/2020] [Accepted: 05/25/2020] [Indexed: 12/17/2022] Open
Abstract
Patients with bone metastases have an increased risk to sustain a pathological fracture as lytic metastatic lesions damage and weaken the bone. In order to prevent fractures, prophylactic treatment is advised for patients with a high fracture risk. Mechanical stabilization of the femur can be provided through femoroplasty, a minimally invasive procedure where bone cement is injected into the lesion, or through internal fixation with intra- or extramedullary implants. Clinicians face the task of determining whether or not prophylactic treatment is required and which treatment would be the most optimal. Finite element (FE) models are promising tools that could support this decision process. The aim of this paper is to provide an overview of the state-of-the-art in FE modeling for the treatment decision of metastatic bone lesions in the femur. First, we will summarize the clinical and mechanical results of femoroplasty as a prophylactic treatment method. Secondly, current FE models for fracture risk assessment of metastatic femurs will be reviewed and the remaining challenges for clinical implementation will be discussed. Thirdly, we will elaborate on the simulation of femoroplasty in FE models and discuss future opportunities. Femoroplasty has already proven to effectively relieve pain and improve functionality, but there remains uncertainty whether it provides sufficient mechanical strengthening to prevent pathological fractures. FE models could help to select appropriate candidates for whom femoroplasty provides sufficient increase in strength and to further improve the mechanical benefit by optimizing the locations for cement augmentation.
Collapse
Affiliation(s)
- Amelie Sas
- Biomechanics Section, KU Leuven, Leuven, Belgium
| | - Esther Tanck
- Orthopaedic Research Laboratory, Radboud Institute for Health Sciences, Radboud University Medical Center, Nijmegen, the Netherlands
| | - An Sermon
- Department of Traumatology, University Hospitals Gasthuisberg, Leuven, Belgium and Department of Development and Regeneration, KU Leuven, Leuven, Belgium
| | | |
Collapse
|
12
|
Varga P, Willie BM, Stephan C, Kozloff KM, Zysset PK. Finite element analysis of bone strength in osteogenesis imperfecta. Bone 2020; 133:115250. [PMID: 31981754 PMCID: PMC7383936 DOI: 10.1016/j.bone.2020.115250] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/19/2019] [Revised: 01/17/2020] [Accepted: 01/19/2020] [Indexed: 12/11/2022]
Abstract
As a dedicated experimentalist, John Currey praised the high potential of finite element (FE) analysis but also recognized its critical limitations. The application of the FE methodology to bone tissue is reviewed in the light of his enthusiastic and colorful statements. In the past decades, FE analysis contributed substantially to the understanding of structure-function properties in the hierarchical organization of bone and to the simulation of bone adaptation. The systematic experimental validation of FE analysis of bone strength in anatomical locations at risk of fracture led to its application in clinical studies to evaluate efficacy of antiresorptive or anabolic treatment of bone fragility. Beyond the successful analyses of healthy or osteoporotic bone, FE analysis becomes increasingly involved in the investigation of other fragility-related bone diseases. The case of osteogenesis imperfecta (OI) is exposed, the multiscale alterations of the bone tissue and the effect of treatment summarized. A few FE analyses attempting to answer open questions in OI are then reported. An original study is finally presented that explored the structural properties of the Brtl/+ murine model of OI type IV subjected to sclerostin neutralizing antibody treatment using microFE analysis. The use of identical material properties in the four-point bending FE simulations of the femora reproduced not only the experimental values but also the statistical comparisons examining the effect of disease and treatment. Further efforts are needed to build upon the extraordinary legacy of John Currey and clarify the impact of different bone diseases on the hierarchical mechanical properties of bone.
Collapse
Affiliation(s)
- Peter Varga
- AO Research Institute Davos, Davos, Switzerland.
| | - Bettina M Willie
- Research Centre, Shriners Hospital for Children-Canada, Montreal, Canada; Department of Pediatric Surgery, McGill University, Montreal, Canada
| | - Chris Stephan
- Department of Orthopaedic Surgery, University of Michigan, Ann Arbor, USA
| | - Kenneth M Kozloff
- Department of Orthopaedic Surgery, University of Michigan, Ann Arbor, USA
| | - Philippe K Zysset
- ARTORG Centre for Biomedical Engineering Research, University of Bern, Bern, Switzerland
| |
Collapse
|
13
|
Wu Y, Morgan EF. Effect of fabric on the accuracy of computed tomography-based finite element analyses of the vertebra. Biomech Model Mechanobiol 2020; 19:505-517. [PMID: 31506861 PMCID: PMC7062572 DOI: 10.1007/s10237-019-01225-2] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2018] [Accepted: 08/31/2019] [Indexed: 10/26/2022]
Abstract
Quantitative computed tomography (QCT)-based finite element (FE) models of the vertebra are widely used in studying spine biomechanics and mechanobiology, but their accuracy has not been fully established. Although the models typically assign material properties based only on local bone mineral density (BMD), the mechanical behavior of trabecular bone also depends on fabric. The goal of this study was to determine the effect of incorporating measurements of fabric on the accuracy of FE predictions of vertebral deformation. Accuracy was assessed by using displacement fields measured via digital volume correlation-applied to time-lapse microcomputed tomography (μCT)-as the gold standard. Two QCT-based FE models were generated from human L1 vertebrae (n = 11): the entire vertebral body and a cuboid-shaped portion of the trabecular centrum [dimensions: (20-30) × (15-20) × (15-20) mm3]. For axial compression boundary conditions, there was no difference (p = 0.40) in the accuracy of the FE-computed displacements for models using material properties based on local values of BMD versus those using material properties based on local values of fabric and volume fraction. However, when using BMD-based material properties, errors were higher for the vertebral-body models (8.4-50.1%) than cuboid models (1.5-19.6%), suggesting that these properties are inaccurate in the peripheral regions of the centrum. Errors also increased when assuming that the cuboid region experienced uniaxial loading during axial compression of the vertebra. These findings indicate that a BMD-based constitutive model is not sufficient for the peripheral region of the vertebral body when seeking accurate QCT-based FE modeling of the vertebra.
Collapse
Affiliation(s)
- Yuanqiao Wu
- Department of Mechanical Engineering, Boston University, 110 Cummington Mall, Boston, MA, 02215, USA.
| | - Elise F Morgan
- Department of Mechanical Engineering, Boston University, 110 Cummington Mall, Boston, MA, 02215, USA
| |
Collapse
|
14
|
Bouxsein ML, Zysset P, Glüer CC, McClung M, Biver E, Pierroz DD, Ferrari SL. Perspectives on the non-invasive evaluation of femoral strength in the assessment of hip fracture risk. Osteoporos Int 2020; 31:393-408. [PMID: 31900541 DOI: 10.1007/s00198-019-05195-0] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/23/2019] [Accepted: 10/04/2019] [Indexed: 10/25/2022]
Abstract
UNLABELLED We reviewed the experimental and clinical evidence that hip bone strength estimated by BMD and/or finite element analysis (FEA) reflects the actual strength of the proximal femur and is associated with hip fracture risk and its changes upon treatment. INTRODUCTION The risk of hip fractures increases exponentially with age due to a progressive loss of bone mass, deterioration of bone structure, and increased incidence of falls. Areal bone mineral density (aBMD), measured by dual-energy X-ray absorptiometry (DXA), is the most used surrogate marker of bone strength. However, age-related declines in bone strength exceed those of aBMD, and the majority of fractures occur in those who are not identified as osteoporotic by BMD testing. With hip fracture incidence increasing worldwide, the development of accurate methods to estimate bone strength in vivo would be very useful to predict the risk of hip fracture and to monitor the effects of osteoporosis therapies. METHODS We reviewed experimental and clinical evidence regarding the association between aBMD and/orCT-finite element analysis (FEA) estimated femoral strength and hip fracture risk as well as their changes with treatment. RESULTS Femoral aBMD and bone strength estimates by CT-FEA explain a large proportion of femoral strength ex vivo and predict hip fracture risk in vivo. Changes in femoral aBMD are strongly associated with anti-fracture efficacy of osteoporosis treatments, though comparable data for FEA are currently not available. CONCLUSIONS Hip aBMD and estimated femoral strength are good predictors of fracture risk and could potentially be used as surrogate endpoints for fracture in clinical trials. Further improvements of FEA may be achieved by incorporating trabecular orientations, enhanced cortical modeling, effects of aging on bone tissue ductility, and multiple sideway fall loading conditions.
Collapse
Affiliation(s)
- M L Bouxsein
- Center for Advanced Orthopedic Studies, Beth Israel Deaconess Medical Center, and Department of Orthopedic Surgery, Harvard Medical School, Boston, MA, USA
| | - P Zysset
- ARTORG Center for Biomedical Engineering Research, University of Bern, Bern, Switzerland
| | - C C Glüer
- Section of Biomedical Imaging, Department of Radiology and Neuroradiology, University Medical Center of Schleswig-Holstein, Campus Kiel, Kiel, Germany
| | - M McClung
- Oregon Osteoporosis Center, Portland, OR, USA
- Mary MacKillop Institute for Health Research, Australian Catholic University, Melbourne, VIC, Australia
| | - E Biver
- Division of Bone Disease, Department of Internal Medicine Specialties, Faculty of Medicine, Geneva University Hospital, Geneva, Switzerland
| | - D D Pierroz
- International Osteoporosis Foundation (IOF), Nyon, Switzerland
| | - S L Ferrari
- Division of Bone Disease, Department of Internal Medicine Specialties, Faculty of Medicine, Geneva University Hospital, Geneva, Switzerland.
| | | |
Collapse
|
15
|
Jiang H, Robinson DL, McDonald M, Lee PVS, Kontulainen SA, Johnston JD, Yates CJ, Wark JD. Predicting experimentally-derived failure load at the distal radius using finite element modelling based on peripheral quantitative computed tomography cross-sections (pQCT-FE): A validation study. Bone 2019; 129:115051. [PMID: 31472298 DOI: 10.1016/j.bone.2019.115051] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/12/2019] [Revised: 08/06/2019] [Accepted: 08/27/2019] [Indexed: 01/08/2023]
Abstract
Dual energy X-ray absorptiometry, the current clinical criterion method for osteoporosis diagnosis, has limitations in identifying individuals with increased fracture risk, especially at the distal radius. Peripheral quantitative computed tomography (pQCT) can provide volumetric bone density data, as well as information on bone geometry, which makes it possible to establish finite element (FE) models of the distal radius from which bone strength and stiffness can be calculated. In this study, we compared experimental mechanical failure load data of the forearm with pQCT- based FE (pQCT-FE) modelling properties. Sixteen cadaveric forearm specimens were experimentally loaded until failure. Estimated stiffness and strength variables of compression, shear, bending and torsion were calculated from pQCT-FE modelling of single cross-sections of 0.2 × 0.2 × 2.4 mm of the radius pQCT image. A moderate-to-strong coefficient of determination (r2) was observed between experimental failure load and pQCT-FE variables. The highest r2 was observed for bending stiffness (r2 = 0.83). This study validates the use of pQCT-FE in the assessment of distal radius bone strength for future studies.
Collapse
Affiliation(s)
- Hongyuan Jiang
- Department of Medicine, Royal Melbourne Hospital, University of Melbourne, Melbourne, Australia
| | - Dale L Robinson
- Department of Biomedical Engineering, University of Melbourne, Melbourne, Australia
| | - Matthew McDonald
- Department of Mechanical Engineering, University of Saskatchewan, Saskatoon, Canada
| | - Peter V S Lee
- Department of Biomedical Engineering, University of Melbourne, Melbourne, Australia
| | | | - James D Johnston
- Department of Mechanical Engineering, University of Saskatchewan, Saskatoon, Canada
| | - Christopher J Yates
- Department of Medicine, Royal Melbourne Hospital, University of Melbourne, Melbourne, Australia; Department of Diabetes and Endocrinology, Royal Melbourne Hospital, Melbourne, Australia; Bone and Mineral Medicine, Royal Melbourne Hospital, Melbourne, Australia
| | - John D Wark
- Department of Medicine, Royal Melbourne Hospital, University of Melbourne, Melbourne, Australia; Department of Diabetes and Endocrinology, Royal Melbourne Hospital, Melbourne, Australia; Bone and Mineral Medicine, Royal Melbourne Hospital, Melbourne, Australia.
| |
Collapse
|
16
|
Benca E, Synek A, Amini M, Kainberger F, Hirtler L, Windhager R, Mayr W, Pahr DH. QCT-based finite element prediction of pathologic fractures in proximal femora with metastatic lesions. Sci Rep 2019; 9:10305. [PMID: 31311994 PMCID: PMC6635505 DOI: 10.1038/s41598-019-46739-y] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2018] [Accepted: 07/04/2019] [Indexed: 11/17/2022] Open
Abstract
Predicting pathologic fractures in femora with metastatic lesions remains a clinical challenge. Currently used guidelines are inaccurate, especially to predict non-impeding fractures. This study evaluated the ability of a nonlinear quantitative computed tomography (QCT)-based homogenized voxel finite element (hvFE) model to predict patient-specific pathologic fractures. The hvFE model was generated highly automated from QCT images of human femora. The femora were previously loaded in a one-legged stance setup in order to assess stiffness, failure load, and fracture location. One femur of each pair was tested in its intact state, while the contralateral femur included a simulated lesion on either the superolateral- or the inferomedial femoral neck. The hvFE model predictions of the stiffness (0.47 < R2 < 0.94), failure load (0.77 < R2 < 0.98), and exact fracture location (68%) were in good agreement with the experimental data. However, the model underestimated the failure load by a factor of two. The hvFE models predicted significant differences in stiffness and failure load for femora with superolateral- and inferomedial lesions. In contrast, standard clinical guidelines predicted identical fracture risk for both lesion sites. This study showed that the subject-specific QCT-based hvFE model could predict the effect of metastatic lesions on the biomechanical behaviour of the proximal femur with moderate computational time and high level of automation and could support treatment strategy in patients with metastatic bone disease.
Collapse
Affiliation(s)
- Emir Benca
- Department of Orthopedics and Trauma Surgery, Medical University of Vienna, 1090, Vienna, Austria.
| | - Alexander Synek
- Institute for Lightweight Design and Structural Biomechanics, TU Wien, 1060, Vienna, Austria
| | - Morteza Amini
- Institute for Lightweight Design and Structural Biomechanics, TU Wien, 1060, Vienna, Austria
| | - Franz Kainberger
- Department of Radiology, Medical University of Vienna, 1090, Vienna, Austria
| | - Lena Hirtler
- Division of Anatomy, Center for Anatomy and Cell Biology, Medical University of Vienna, 1090, Vienna, Austria
| | - Reinhard Windhager
- Department of Orthopedics and Trauma Surgery, Medical University of Vienna, 1090, Vienna, Austria
| | - Winfried Mayr
- Center for Medical Physics and Biomedical Engineering, Medical University of Vienna, 1090, Vienna, Austria
| | - Dieter H Pahr
- Institute for Lightweight Design and Structural Biomechanics, TU Wien, 1060, Vienna, Austria
| |
Collapse
|
17
|
Väänänen SP, Grassi L, Venäläinen MS, Matikka H, Zheng Y, Jurvelin JS, Isaksson H. Automated segmentation of cortical and trabecular bone to generate finite element models for femoral bone mechanics. Med Eng Phys 2019; 70:19-28. [PMID: 31280927 DOI: 10.1016/j.medengphy.2019.06.015] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2018] [Revised: 05/16/2019] [Accepted: 06/23/2019] [Indexed: 02/02/2023]
Abstract
Finite element (FE) models based on quantitative computed tomography (CT) images are better predictors of bone strength than conventional areal bone mineral density measurements. However, FE models require manual segmentation of the femur, which is not clinically applicable. This study developed a method for automated FE analyses from clinical CT images. Clinical in-vivo CT images of 13 elderly female subjects were collected to evaluate the method. Secondly, proximal cadaver femurs were harvested and imaged with clinical CT (N = 17). Of these femurs, 14 were imaged with µCT and three had earlier been tested experimentally in stance-loading, while collecting surface deformations with digital image correlation. Femurs were segmented from clinical CT images using an automated method, based on the segmentation tool Stradwin. The method automatically distinguishes trabecular and cortical bone, corrects partial volume effect and generates input for FE analysis. The manual and automatic segmentations agreed within about one voxel for in-vivo subjects (0.99 ± 0.23 mm) and cadaver femurs (0.21 ± 0.07 mm). The strains from the FE predictions closely matched with the experimentally measured strains (R2 = 0.89). The method can automatically generate meshes suitable for FE analysis. The method may bring us one step closer to enable clinical usage of patient-specific FE analyses.
Collapse
Affiliation(s)
- Sami P Väänänen
- Department of Applied Physics, University of Eastern Finland, POB 1627, FIN-70211 Kuopio, Finland; Department of Clinical Radiology, Diagnostic Imaging Center, Kuopio University Hospital, POB 100, 70029 Kuopio, Finland; Department of Orthopaedics, Traumatology and Hand Surgery, Kuopio University Hospital, POB 100, FIN-70029 Kuopio, Finland; Department of Medical Physics, Central Finland Central Hospital, Keskussairaalantie 19, FIN-40620 Jyväskylä, Finland.
| | - Lorenzo Grassi
- Department of Biomedical Engineering, Lund University, BMC D13, 221 84 Lund, Sweden.
| | - Mikko S Venäläinen
- Department of Applied Physics, University of Eastern Finland, POB 1627, FIN-70211 Kuopio, Finland; Turku Bioscience Centre, University of Turku and Åbo Akademi University, Tykistökatu 6, FIN-20520 Turku, Finland.
| | - Hanna Matikka
- Department of Clinical Radiology, Diagnostic Imaging Center, Kuopio University Hospital, POB 100, 70029 Kuopio, Finland.
| | - Yi Zheng
- Department of Physics, Technical University of Denmark, Fysikvej, building 311, 2800 Kgs. Lyngby, Denmark.
| | - Jukka S Jurvelin
- Department of Applied Physics, University of Eastern Finland, POB 1627, FIN-70211 Kuopio, Finland.
| | - Hanna Isaksson
- Department of Biomedical Engineering, Lund University, BMC D13, 221 84 Lund, Sweden.
| |
Collapse
|
18
|
Kluess D, Soodmand E, Lorenz A, Pahr D, Schwarze M, Cichon R, Varady PA, Herrmann S, Buchmeier B, Schröder C, Lehner S, Kebbach M. A round-robin finite element analysis of human femur mechanics between seven participating laboratories with experimental validation. Comput Methods Biomech Biomed Engin 2019; 22:1020-1031. [PMID: 31084272 DOI: 10.1080/10255842.2019.1615481] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Finite element analysis is a common tool that has been used for the past few decades to predict the mechanical behavior of bone. However, to our knowledge, there are no round-robin finite element analyses of long human bones with more than two participating biomechanics laboratories published yet, where the results of the experimental tests were not known in advance. We prepared a fresh-frozen human femur for a compression test in a universal testing machine measuring the strains at 10 bone locations as well as the deformation of the bone in terms of the displacement of the loading point at a load of 2 kN. The computed tomography data of the bone with a calibration phantom as well as the orientation of the bone in the testing machine with the according boundary conditions were delivered to seven participating laboratories. These were asked to perform a finite element analysis simulating the experimental setup and deliver their results to the coordinator without knowing the experimental results. Resultantly, four laboratories had deviations from the experimentally measured strains of less than 40%, and three laboratories had deviations of their numerically determined values compared to the experimental data of more than 120%. These deviations are thought to be based on different material laws and material data, as well as different material mapping methods. Investigations will be conducted to clarify and assess the reasons for the large deviations in the numerical data. It was shown that the precision of finite element models of the human femur is not yet as developed as desired by the biomechanics community.
Collapse
Affiliation(s)
- Daniel Kluess
- a Department of Orthopaedics , Rostock University Medical Center , Rostock , Germany
| | - Ehsan Soodmand
- a Department of Orthopaedics , Rostock University Medical Center , Rostock , Germany
| | - Andrea Lorenz
- b TU Wien Institute of Lightweight Design and Structural Biomechanics , Vienna , Austria
| | - Dieter Pahr
- b TU Wien Institute of Lightweight Design and Structural Biomechanics , Vienna , Austria
| | - Michael Schwarze
- c Department of Orthopaedic Surgery , Hannover Medical School , Hannover , Germany
| | - Robert Cichon
- d Chair of Mechanics and Robotics , University of Duisburg-Essen , Duisburg , Germany
| | - Patrick A Varady
- e BG Unfallklinik Murnau Institute for Biomechanics , Murnau am Staffelsee , Germany
| | - Sven Herrmann
- e BG Unfallklinik Murnau Institute for Biomechanics , Murnau am Staffelsee , Germany
| | | | - Christian Schröder
- g Orthopädie & Traumatologie/Orthopedics & Traumatology , TÜV SÜD Product Service GmbH , München/Munich , Germany
| | - Stefan Lehner
- h Faculty Mechanical Engineering and Mechatronics , Deggendorf Institute of Technology , Deggendorf , Germany
| | - Maeruan Kebbach
- a Department of Orthopaedics , Rostock University Medical Center , Rostock , Germany
| |
Collapse
|
19
|
Chandran V, Maquer G, Gerig T, Zysset P, Reyes M. Supervised learning for bone shape and cortical thickness estimation from CT images for finite element analysis. Med Image Anal 2019; 52:42-55. [DOI: 10.1016/j.media.2018.11.001] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2018] [Revised: 08/24/2018] [Accepted: 11/02/2018] [Indexed: 01/31/2023]
|
20
|
Panyasantisuk J, Dall'Ara E, Pretterklieber M, Pahr DH, Zysset PK. Mapping anisotropy improves QCT-based finite element estimation of hip strength in pooled stance and side-fall load configurations. Med Eng Phys 2018; 59:36-42. [PMID: 30131112 DOI: 10.1016/j.medengphy.2018.06.004] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2017] [Revised: 03/26/2018] [Accepted: 06/24/2018] [Indexed: 02/05/2023]
Abstract
Hip fractures are one of the most severe consequences of osteoporosis. Compared to the clinical standard of DXA-based aBMD at the femoral neck, QCT-based FEA delivers a better surrogate of femoral strength and gains acceptance for the calculation of hip fracture risk when a CT reconstruction is available. Isotropic, homogenised voxel-based, finite element (hvFE) models are widely used to estimate femoral strength in cross-sectional and longitudinal clinical studies. However, fabric anisotropy is a classical feature of the architecture of the proximal femur and the second determinant of the homogenised mechanical properties of trabecular bone. Due to the limited resolution, fabric anisotropy cannot be derived from clinical CT reconstructions. Alternatively, fabric anisotropy can be extracted from HR-pQCT images of cadaveric femora. In this study, fabric anisotropy from HR-pQCT images was mapped onto QCT-based hvFE models of 71 human proximal femora for which both HR-pQCT and QCT images were available. Stiffness and ultimate load computed from anisotropic hvFE models were compared with previous biomechanical tests in both stance and side-fall configurations. The influence of using the femur-specific versus a mean fabric distribution on the hvFE predictions was assessed. Femur-specific and mean fabric enhance the prediction of experimental ultimate force for the pooled, i.e. stance and side-fall, (isotropic: r2=0.81, femur-specific fabric: r2=0.88, mean fabric: r2=0.86,p<0.001) but not for the individual configurations. Fabric anisotropy significantly improves bone strength prediction for the pooled configurations, and mapped fabric provides a comparable prediction to true fabric. The mapping of fabric anisotropy is therefore expected to help generate more accurate QCT-based hvFE models of the proximal femur for personalised or multiple load configurations.
Collapse
Affiliation(s)
- J Panyasantisuk
- Institute for Surgical Technology and Biomechanics, University of Bern, Switzerland
| | - E Dall'Ara
- Department of Oncology and Metabolism and INSIGNEO, Institute for in silico Medicine, University of Sheffield, United Kingdom
| | | | - D H Pahr
- Institute for Lightweight Design and Structural Biomechanics, Vienna University of Technology, Austria; Department for Anatomy and Biomechanics, Karl Landsteiner Private University for Health Sciences, Austria
| | - P K Zysset
- Institute for Surgical Technology and Biomechanics, University of Bern, Switzerland.
| |
Collapse
|
21
|
Enns-Bray WS, Ferguson SJ, Helgason B. Strain rate dependency of bovine trabecular bone under impact loading at sideways fall velocity. J Biomech 2018; 75:46-52. [DOI: 10.1016/j.jbiomech.2018.04.042] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2018] [Revised: 04/23/2018] [Accepted: 04/23/2018] [Indexed: 11/16/2022]
|
22
|
Tsegai ZJ, Skinner MM, Pahr DH, Hublin J, Kivell TL. Systemic patterns of trabecular bone across the human and chimpanzee skeleton. J Anat 2018; 232:641-656. [PMID: 29344941 PMCID: PMC5835784 DOI: 10.1111/joa.12776] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/07/2017] [Indexed: 12/18/2022] Open
Abstract
Aspects of trabecular bone architecture are thought to reflect regional loading of the skeleton, and thus differ between primate taxa with different locomotor and postural modes. However, there are several systemic factors that affect bone structure that could contribute to, or be the primary factor determining, interspecific differences in bone structure. These systemic factors include differences in genetic regulation, sensitivity to loading, hormone levels, diet, and activity levels. Improved understanding of inter-/intraspecific variability, and variability across the skeleton of an individual, is required to interpret properly potential functional signals present within trabecular structure. Using a whole-region method of analysis, we investigated trabecular structure throughout the skeleton of humans and chimpanzees. Trabecular bone volume fraction (BV/TV), degree of anisotropy (DA) and trabecular thickness (Tb.Th) were quantified from high resolution micro-computed tomographic scans of the humeral and femoral head, third metacarpal and third metatarsal head, distal tibia, talus and first thoracic vertebra. We found that BV/TV is, in most anatomical sites, significantly higher in chimpanzees than in humans, suggesting a systemic difference in trabecular structure unrelated to local loading regime. Differences in BV/TV between the forelimb and hindlimb did not clearly reflect differences in locomotor loading in the study taxa. There were no clear systemic differences between the taxa in DA and, as such, this parameter might reflect function and relate to differences in joint loading. This systemic approach reveals both the pattern of variability across the skeleton and between taxa, and helps identify those features of trabecular structure that may relate to joint function.
Collapse
Affiliation(s)
- Zewdi J. Tsegai
- Department of Human EvolutionMax Planck Institute for Evolutionary AnthropologyLeipzigGermany
| | - Matthew M. Skinner
- Department of Human EvolutionMax Planck Institute for Evolutionary AnthropologyLeipzigGermany
- Skeletal Biology Research CentreSchool of Anthropology and ConservationUniversity of KentCanterburyUK
| | - Dieter H. Pahr
- Institute of Lightweight Design and Structural BiomechanicsVienna University of TechnologyViennaAustria
| | - Jean‐Jacques Hublin
- Department of Human EvolutionMax Planck Institute for Evolutionary AnthropologyLeipzigGermany
| | - Tracy L. Kivell
- Department of Human EvolutionMax Planck Institute for Evolutionary AnthropologyLeipzigGermany
- Skeletal Biology Research CentreSchool of Anthropology and ConservationUniversity of KentCanterburyUK
| |
Collapse
|
23
|
Enns-Bray W, Bahaloo H, Fleps I, Ariza O, Gilchrist S, Widmer R, Guy P, Pálsson H, Ferguson S, Cripton P, Helgason B. Material mapping strategy to improve the predicted response of the proximal femur to a sideways fall impact. J Mech Behav Biomed Mater 2018; 78:196-205. [DOI: 10.1016/j.jmbbm.2017.10.033] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2016] [Revised: 09/25/2017] [Accepted: 10/26/2017] [Indexed: 11/29/2022]
|
24
|
Chandran V, Reyes M, Zysset P. A novel registration-based methodology for prediction of trabecular bone fabric from clinical QCT: A comprehensive analysis. PLoS One 2017; 12:e0187874. [PMID: 29176881 PMCID: PMC5703488 DOI: 10.1371/journal.pone.0187874] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2017] [Accepted: 10/29/2017] [Indexed: 11/23/2022] Open
Abstract
Osteoporosis leads to hip fractures in aging populations and is diagnosed by modern medical imaging techniques such as quantitative computed tomography (QCT). Hip fracture sites involve trabecular bone, whose strength is determined by volume fraction and orientation, known as fabric. However, bone fabric cannot be reliably assessed in clinical QCT images of proximal femur. Accordingly, we propose a novel registration-based estimation of bone fabric designed to preserve tensor properties of bone fabric and to map bone fabric by a global and local decomposition of the gradient of a non-rigid image registration transformation. Furthermore, no comprehensive analysis on the critical components of this methodology has been previously conducted. Hence, the aim of this work was to identify the best registration-based strategy to assign bone fabric to the QCT image of a patient’s proximal femur. The normalized correlation coefficient and curvature-based regularization were used for image-based registration and the Frobenius norm of the stretch tensor of the local gradient was selected to quantify the distance among the proximal femora in the population. Based on this distance, closest, farthest and mean femora with a distinction of sex were chosen as alternative atlases to evaluate their influence on bone fabric prediction. Second, we analyzed different tensor mapping schemes for bone fabric prediction: identity, rotation-only, rotation and stretch tensor. Third, we investigated the use of a population average fabric atlas. A leave one out (LOO) evaluation study was performed with a dual QCT and HR-pQCT database of 36 pairs of human femora. The quality of the fabric prediction was assessed with three metrics, the tensor norm (TN) error, the degree of anisotropy (DA) error and the angular deviation of the principal tensor direction (PTD). The closest femur atlas (CTP) with a full rotation (CR) for fabric mapping delivered the best results with a TN error of 7.3 ± 0.9%, a DA error of 6.6 ± 1.3% and a PTD error of 25 ± 2°. The closest to the population mean femur atlas (MTP) using the same mapping scheme yielded only slightly higher errors than CTP for substantially less computing efforts. The population average fabric atlas yielded substantially higher errors than the MTP with the CR mapping scheme. Accounting for sex did not bring any significant improvements. The identified fabric mapping methodology will be exploited in patient-specific QCT-based finite element analysis of the proximal femur to improve the prediction of hip fracture risk.
Collapse
Affiliation(s)
- Vimal Chandran
- Institute of Surgical Technology and Biomechanics, University of Bern, Bern, Switzerland
- * E-mail:
| | - Mauricio Reyes
- Institute of Surgical Technology and Biomechanics, University of Bern, Bern, Switzerland
| | - Philippe Zysset
- Institute of Surgical Technology and Biomechanics, University of Bern, Bern, Switzerland
| |
Collapse
|
25
|
Ramos-Infante SJ, Pérez MA. In vitro and in silico characterization of open-cell structures of trabecular bone. Comput Methods Biomech Biomed Engin 2017; 20:1562-1570. [DOI: 10.1080/10255842.2017.1390086] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- S. J. Ramos-Infante
- M2BE-Multiscale in Mechanical and Biological Engineering, Instituto de Investigación en Ingeniería de Aragón (I3A), Universidad de Zaragoza Campus Río Ebro, Zaragoza, Spain
| | - M. A. Pérez
- M2BE-Multiscale in Mechanical and Biological Engineering, Instituto de Investigación en Ingeniería de Aragón (I3A), Universidad de Zaragoza Campus Río Ebro, Zaragoza, Spain
| |
Collapse
|
26
|
Fatigue failure of plated osteoporotic proximal humerus fractures is predicted by the strain around the proximal screws. J Mech Behav Biomed Mater 2017; 75:68-74. [DOI: 10.1016/j.jmbbm.2017.07.004] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2017] [Revised: 06/19/2017] [Accepted: 07/03/2017] [Indexed: 01/20/2023]
|
27
|
Taghizadeh E, Chandran V, Reyes M, Zysset P, Büchler P. Statistical analysis of the inter-individual variations of the bone shape, volume fraction and fabric and their correlations in the proximal femur. Bone 2017; 103:252-261. [PMID: 28732775 DOI: 10.1016/j.bone.2017.07.012] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/20/2017] [Revised: 06/22/2017] [Accepted: 07/11/2017] [Indexed: 10/19/2022]
Abstract
Including structural information of trabecular bone improves the prediction of bone strength and fracture risk. However, this information is available in clinical CT scans, only for peripheral bones. We hypothesized that a correlation exists between the shape of the bone, its volume fraction (BV/TV) and fabric, which could be characterized using statistical modeling. High-resolution peripheral computed tomography (HR-pQCT) images of 73 proximal femurs were used to build a combined statistical model of shape, BV/TV and fabric. The model was based on correspondence established by image registration and by morphing of a finite element mesh describing the spatial distribution of the bone properties. Results showed no correlation between the distribution of bone shape, BV/TV and fabric. Only the first mode of variation associated with density and orientation showed a strong relationship (R2>0.8). In addition, the model showed that the anisotropic information of the proximal femur does not vary significantly in a population of healthy, osteoporotic and osteopenic samples. In our dataset, the average anisotropy of the population was able to provide a close approximation of the patient-specific anisotropy. These results were confirmed by homogenized finite element (hFE) analyses, which showed that the biomechanical behavior of the proximal femur was not significantly different when the average anisotropic information of the population was used instead of patient-specific fabric extracted from HR-pQCT. Based on these findings, it can be assumed that the fabric information of the proximal femur follows a similar structure in an elderly population of healthy, osteopenic and osteoporotic proximal femurs.
Collapse
Affiliation(s)
- Elham Taghizadeh
- Institute for Surgical Technology and Biomechanics (ISTB), University of Bern, Switzerland
| | - Vimal Chandran
- Institute for Surgical Technology and Biomechanics (ISTB), University of Bern, Switzerland
| | - Mauricio Reyes
- Institute for Surgical Technology and Biomechanics (ISTB), University of Bern, Switzerland
| | - Philippe Zysset
- Institute for Surgical Technology and Biomechanics (ISTB), University of Bern, Switzerland
| | - Philippe Büchler
- Institute for Surgical Technology and Biomechanics (ISTB), University of Bern, Switzerland.
| |
Collapse
|
28
|
Lopes VM, Neto MA, Amaro AM, Roseiro LM, Paulino M. FE and experimental study on how the cortex material properties of synthetic femurs affect strain levels. Med Eng Phys 2017. [DOI: 10.1016/j.medengphy.2017.06.001] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
29
|
Maquer G, Bürki A, Nuss K, Zysset PK, Tannast M. Head-Neck Osteoplasty has Minor Effect on the Strength of an Ovine Cam-FAI Model: In Vitro and Finite Element Analyses. Clin Orthop Relat Res 2016; 474:2633-2640. [PMID: 27535284 PMCID: PMC5085938 DOI: 10.1007/s11999-016-5024-8] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/21/2016] [Accepted: 08/03/2016] [Indexed: 01/31/2023]
Abstract
BACKGROUND Osteochondroplasty of the head-neck region is performed on patients with cam femoroacetabular impingement (FAI) without fully understanding its repercussion on the integrity of the femur. Cam-type FAI can be surgically and reproducibly induced in the ovine femur, which makes it suitable for studying corrective surgery in a consistent way. Finite element models built on quantitative CT (QCT) are computer tools that can be used to predict femoral strength and evaluate the mechanical effect of surgical correction. QUESTIONS/PURPOSES We asked: (1) What is the effect of a resection of the superolateral aspect of the ovine femoral head-neck junction on failure load? (2) How does the failure load after osteochondroplasty compare with reported forces from activities of daily living in sheep? (3) How do failure loads and failure locations from the computer simulations compare with the experiments? METHODS Osteochondroplasties (3, 6, 9 mm) were performed on one side of 18 ovine femoral pairs with the contralateral intact side as a control. The 36 femurs were scanned via QCT from which specimen-specific computer models were built. Destructive compression tests then were conducted experimentally using a servohydraulic testing system and numerically via the computer models. Safety factors were calculated as the ratio of the maximal force measured in vivo by telemeterized hip implants during the sheep's walking and running activities to the failure load. The simulated failure loads and failure locations from the computer models were compared with the experimental results. RESULTS Failure loads were reduced by 5% (95% CI, 2%-8%) for the 3-mm group (p = 0.0089), 10% (95% CI, 6%-14%) for the 6-mm group (p = 0.0015), and 19% (95% CI, 13%-26%) for the 9-mm group (p = 0.0097) compared with the controls. Yet, the weakest specimen still supported more than 2.4 times the peak load during running. Strong correspondence was found between the simulated and experimental failure loads (R2 = 0.83; p < 0.001) and failure locations. CONCLUSIONS The resistance of ovine femurs to fracture decreased with deeper resections. However, under in vitro testing conditions, the effect on femoral strength remains small even after 9 mm correction, suggesting that femoral head-neck osteochondroplasty could be done safely on the ovine femur. QCT-based finite element models were able to predict weakening of the femur resulting from the osteochondroplasty. CLINICAL RELEVANCE The ovine femur provides a seemingly safe platform for scientific evaluation of FAI. It also appears that computer models based on preoperative CT scans may have the potential to provide patient-specific guidelines for preventing overcorrection of cam FAI.
Collapse
Affiliation(s)
- Ghislain Maquer
- Institute for Surgical Technology and Biomechanics, University of Bern, Stauffacherstrasse 78, 3014, Bern, Switzerland.
| | - Alexander Bürki
- Institute for Surgical Technology and Biomechanics, University of Bern, Stauffacherstrasse 78, 3014, Bern, Switzerland
| | - Katja Nuss
- Musculoskeletal Research Unit, Vetsuisse Faculty, University of Zurich, Zürich, Switzerland
| | - Philippe K Zysset
- Institute for Surgical Technology and Biomechanics, University of Bern, Stauffacherstrasse 78, 3014, Bern, Switzerland
| | - Moritz Tannast
- Musculoskeletal Research Unit, Vetsuisse Faculty, University of Zurich, Zürich, Switzerland
- Department of Orthopaedic Surgery, Inselspital, University of Bern, Bern, Switzerland
| |
Collapse
|
30
|
Knowles NK, Reeves JM, Ferreira LM. Quantitative Computed Tomography (QCT) derived Bone Mineral Density (BMD) in finite element studies: a review of the literature. J Exp Orthop 2016; 3:36. [PMID: 27943224 PMCID: PMC5234499 DOI: 10.1186/s40634-016-0072-2] [Citation(s) in RCA: 60] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/30/2016] [Accepted: 11/30/2016] [Indexed: 01/23/2023] Open
Abstract
BACKGROUND Finite element modeling of human bone provides a powerful tool to evaluate a wide variety of outcomes in a highly repeatable and parametric manner. These models are most often derived from computed tomography data, with mechanical properties related to bone mineral density (BMD) from the x-ray energy attenuation provided from this data. To increase accuracy, many researchers report the use of quantitative computed tomography (QCT), in which a calibration phantom is used during image acquisition to improve the estimation of BMD. Since model accuracy is dependent on the methods used in the calculation of BMD and density-mechanical property relationships, it is important to use relationships developed for the same anatomical location and using the same scanner settings, as these may impact model accuracy. The purpose of this literature review is to report the relationships used in the conversion of QCT equivalent density measures to ash, apparent, and/or tissue densities in recent finite element (FE) studies used in common density-modulus relationships. For studies reporting experimental validation, the validation metrics and results are presented. RESULTS Of the studies reviewed, 29% reported the use of a dipotassium phosphate (K2HPO4) phantom, 47% a hydroxyapatite (HA) phantom, 13% did not report phantom type, 7% reported use of both K2HPO4 and HA phantoms, and 4% alternate phantom types. Scanner type and/or settings were omitted or partially reported in 31% of studies. The majority of studies used densitometric and/or density-modulus relationships derived from different anatomical locations scanned in different scanners with different scanner settings. The methods used to derive various densitometric relationships are reported and recommendations are provided toward the standardization of reporting metrics. CONCLUSIONS This review assessed the current state of QCT-based FE modeling with use of clinical scanners. It was found that previously developed densitometric relationships vary by anatomical location, scanner type and settings. Reporting of all parameters used when referring to previously developed relationships, or in the development of new relationships, may increase the accuracy and repeatability of future FE models.
Collapse
Affiliation(s)
- Nikolas K. Knowles
- Graduate Program in Biomedical Engineering, The University of Western Ontario, 1151 Richmond St, London, ON Canada
- Roth|McFarlane Hand and Upper Limb Centre, Surgical Mechatronics
Laboratory, St. Josephs Health Care, 268 Grosvenor St, London, ON Canada
- Collaborative Training Program in Musculoskeletal Health Research, and
Bone and Joint Institute, The University of Western Ontario, 1151 Richmond St, London, ON Canada
| | - Jacob M. Reeves
- Roth|McFarlane Hand and Upper Limb Centre, Surgical Mechatronics
Laboratory, St. Josephs Health Care, 268 Grosvenor St, London, ON Canada
- Collaborative Training Program in Musculoskeletal Health Research, and
Bone and Joint Institute, The University of Western Ontario, 1151 Richmond St, London, ON Canada
- Department of Mechanical and Materials Engineering, The University of Western Ontario, 1151 Richmond St, London, ON Canada
| | - Louis M. Ferreira
- Graduate Program in Biomedical Engineering, The University of Western Ontario, 1151 Richmond St, London, ON Canada
- Roth|McFarlane Hand and Upper Limb Centre, Surgical Mechatronics
Laboratory, St. Josephs Health Care, 268 Grosvenor St, London, ON Canada
- Collaborative Training Program in Musculoskeletal Health Research, and
Bone and Joint Institute, The University of Western Ontario, 1151 Richmond St, London, ON Canada
| |
Collapse
|
31
|
Abstract
Beyond bone mineral density (BMD), bone quality designates the mechanical integrity of bone tissue. In vivo images based on X-ray attenuation, such as CT reconstructions, provide size, shape, and local BMD distribution and may be exploited as input for finite element analysis (FEA) to assess bone fragility. Further key input parameters of FEA are the material properties of bone tissue. This review discusses the main determinants of bone mechanical properties and emphasizes the added value, as well as the important assumptions underlying finite element analysis. Bone tissue is a sophisticated, multiscale composite material that undergoes remodeling but exhibits a rather narrow band of tissue mineralization. Mechanically, bone tissue behaves elastically under physiologic loads and yields by cracking beyond critical strain levels. Through adequate cell-orchestrated modeling, trabecular bone tunes its mechanical properties by volume fraction and fabric. With proper calibration, these mechanical properties may be incorporated in quantitative CT-based finite element analysis that has been validated extensively with ex vivo experiments and has been applied increasingly in clinical trials to assess treatment efficacy against osteoporosis.
Collapse
Affiliation(s)
- Dieter H Pahr
- Institute of Lightweight Design and Structural Biomechanics, Vienna University of Technology, Vienna, Austria
| | - Philippe K Zysset
- Institute for Surgical Technology and Biomechanics, University of Bern, Bern, Switzerland.
| |
Collapse
|
32
|
Morphology based anisotropic finite element models of the proximal femur validated with experimental data. Med Eng Phys 2016; 38:1339-1347. [DOI: 10.1016/j.medengphy.2016.08.010] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2016] [Revised: 08/05/2016] [Accepted: 08/30/2016] [Indexed: 11/21/2022]
|
33
|
Dall’Ara E, Eastell R, Viceconti M, Pahr D, Yang L. Experimental validation of DXA-based finite element models for prediction of femoral strength. J Mech Behav Biomed Mater 2016; 63:17-25. [DOI: 10.1016/j.jmbbm.2016.06.004] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2015] [Revised: 05/11/2016] [Accepted: 06/02/2016] [Indexed: 11/26/2022]
|
34
|
Modeling the Mechanical Consequences of Age-Related Trabecular Bone Loss by XFEM Simulation. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2016; 2016:3495152. [PMID: 27403206 PMCID: PMC4925952 DOI: 10.1155/2016/3495152] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/18/2016] [Accepted: 05/23/2016] [Indexed: 11/18/2022]
Abstract
The elderly are more likely to suffer from fracture because of age-related trabecular bone loss. Different bone loss locations and patterns have different effects on bone mechanical properties. Extended finite element method (XFEM) can simulate fracture process and was suited to investigate the effects of bone loss on trabecular bone. Age-related bone loss is indicated by trabecular thinning and loss and may occur at low-strain locations or other random sites. Accordingly, several ideal normal and aged trabecular bone models were created based on different bone loss locations and patterns; then, fracture processes from crack initiation to complete failure of these models were observed by XFEM; finally, the effects of different locations and patterns on trabecular bone were compared. Results indicated that bone loss occurring at low-strain locations was more detrimental to trabecular bone than that occurring at other random sites; meanwhile, the decrease in bone strength caused by trabecular loss was higher than that caused by trabecular thinning, and the effects of vertical trabecular loss on mechanical properties were more severe than horizontal trabecular loss. This study provided a numerical method to simulate trabecular bone fracture and distinguished different effects of the possible occurrence of bone loss locations and patterns on trabecular bone.
Collapse
|
35
|
The influence of the modulus-density relationship and the material mapping method on the simulated mechanical response of the proximal femur in side-ways fall loading configuration. Med Eng Phys 2016; 38:679-689. [PMID: 27185044 DOI: 10.1016/j.medengphy.2016.03.006] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2015] [Revised: 03/04/2016] [Accepted: 03/19/2016] [Indexed: 11/23/2022]
Abstract
Contributing to slow advance of finite element (FE) simulations for hip fracture risk prediction, into clinical practice, could be a lack of consensus in the biomechanics community on how to map properties to the models. Thus, the aim of the present study was first, to systematically quantify the influence of the modulus-density relationship (E-ρ) and the material mapping method (MMM) on the predicted mechanical response of the proximal femur in a side-ways fall (SWF) loading configuration and second, to perform a model-to-model comparison of the predicted mechanical response within the femoral neck for all the specimens tested in the present study, using three different modelling techniques that have yielded good validation outcome in terms of surface strain prediction and whole bone response according to the literature. We found the outcome to be highly dependent on both the E-ρ relationship and the MMM. In addition, we found that the three modelling techniques that have resulted in good validation outcome in the literature yielded different principal strain prediction both on the surface as well as internally in the femoral neck region of the specimens modelled in the present study. We conclude that there exists a need to carry out a more comprehensive validation study for the SWF loading mode to identify which combination of MMMs and E-ρ relationship leads to the best match for whole bone and local mechanical response. The MMMs tested in the present study have been made publicly available at https://simtk.org/home/mitk-gem.
Collapse
|
36
|
Nonlinear quasi-static finite element simulations predict in vitro strength of human proximal femora assessed in a dynamic sideways fall setup. J Mech Behav Biomed Mater 2016; 57:116-27. [DOI: 10.1016/j.jmbbm.2015.11.026] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2015] [Revised: 11/23/2015] [Accepted: 11/28/2015] [Indexed: 11/20/2022]
|
37
|
Fan R, Gong H, Zhang R, Gao J, Jia Z, Hu Y. Quantification of Age-Related Tissue-Level Failure Strains of Rat Femoral Cortical Bones Using an Approach Combining Macrocompressive Test and Microfinite Element Analysis. J Biomech Eng 2016; 138:041006. [DOI: 10.1115/1.4032798] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2015] [Indexed: 12/15/2022]
Abstract
Bone mechanical properties vary with age; meanwhile, a close relationship exists among bone mechanical properties at different levels. Therefore, conducting multilevel analyses for bone structures with different ages are necessary to elucidate the effects of aging on bone mechanical properties at different levels. In this study, an approach that combined microfinite element (micro-FE) analysis and macrocompressive test was established to simulate the failure of male rat femoral cortical bone. Micro-FE analyses were primarily performed for rat cortical bones with different ages to simulate their failure processes under compressive load. Tissue-level failure strains in tension and compression of these cortical bones were then back-calculated by fitting the experimental stress–strain curves. Thus, tissue-level failure strains of rat femoral cortical bones with different ages were quantified. The tissue-level failure strain exhibited a biphasic behavior with age: in the period of skeletal maturity (1–7 months of age), the failure strain gradually increased; when the rat exceeded 7 months of age, the failure strain sharply decreased. In the period of skeletal maturity, both the macro- and tissue-levels mechanical properties showed a large promotion. In the period of skeletal aging (9–15 months of age), the tissue-level mechanical properties sharply deteriorated; however, the macromechanical properties only slightly deteriorated. The age-related changes in tissue-level failure strain were revealed through the analysis of male rat femoral cortical bones with different ages, which provided a theoretical basis to understand the relationship between rat cortical bone mechanical properties at macro- and tissue-levels and decrease of bone strength with age.
Collapse
Affiliation(s)
- Ruoxun Fan
- State Key Laboratory of Automotive Simulation and Control, Jilin University, Changchun 130025, China
- Department of Engineering Mechanics, Jilin University, Nanling Campus, Changchun 130025, China e-mail:
| | - He Gong
- Professor State Key Laboratory of Automotive Simulation and Control, Jilin University, Changchun 130025, China
- Department of Engineering Mechanics, Jilin University, Nanling Campus, Changchun 130025, China e-mail:
| | - Rui Zhang
- Department of Engineering Mechanics, Jilin University, Nanling Campus, Changchun 130025, China
- Key Laboratory for Biomechanics and Mechanobiology of Ministry of Education, School of Biological Science and Medical Engineering, Beihang University, Beijing 10000, China e-mail:
| | - Jiazi Gao
- Department of Engineering Mechanics, Jilin University, Nanling Campus, Changchun 130025, China e-mail:
| | - Zhengbin Jia
- Department of Engineering Mechanics, Jilin University, Nanling Campus, Changchun 130025, China e-mail:
| | - Yanjuan Hu
- School of Mechatronic Engineering, Changchun University of Technology, Changchun 130025, China e-mail:
| |
Collapse
|
38
|
Taghizadeh E, Reyes M, Zysset P, Latypova A, Terrier A, Büchler P. Biomechanical Role of Bone Anisotropy Estimated on Clinical CT Scans by Image Registration. Ann Biomed Eng 2016; 44:2505-2517. [PMID: 26790866 DOI: 10.1007/s10439-016-1551-4] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2015] [Accepted: 01/13/2016] [Indexed: 11/27/2022]
Abstract
Image-based modeling is a popular approach to perform patient-specific biomechanical simulations. Accurate modeling is critical for orthopedic application to evaluate implant design and surgical planning. It has been shown that bone strength can be estimated from the bone mineral density (BMD) and trabecular bone architecture. However, these findings cannot be directly and fully transferred to patient-specific modeling since only BMD can be derived from clinical CT. Therefore, the objective of this study was to propose a method to predict the trabecular bone structure using a µCT atlas and an image registration technique. The approach has been evaluated on femurs and patellae under physiological loading. The displacement and ultimate force for femurs loaded in stance position were predicted with an error of 2.5% and 3.7%, respectively, while predictions obtained with an isotropic material resulted in errors of 7.3% and 6.9%. Similar results were obtained for the patella, where the strain predicted using the registration approach resulted in an improved mean squared error compared to the isotropic model. We conclude that the registration of anisotropic information from of a single template bone enables more accurate patient-specific simulations from clinical image datasets than isotropic model.
Collapse
Affiliation(s)
- Elham Taghizadeh
- Institute for Surgical Technology & Biomechanics, University of Bern, Stauffacherstrasse 78, 3014, Bern, Switzerland
| | - Mauricio Reyes
- Institute for Surgical Technology & Biomechanics, University of Bern, Stauffacherstrasse 78, 3014, Bern, Switzerland
| | - Philippe Zysset
- Institute for Surgical Technology & Biomechanics, University of Bern, Stauffacherstrasse 78, 3014, Bern, Switzerland
| | - Adeliya Latypova
- Laboratory of Biomechanical Orthopedics, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Alexandre Terrier
- Laboratory of Biomechanical Orthopedics, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Philippe Büchler
- Institute for Surgical Technology & Biomechanics, University of Bern, Stauffacherstrasse 78, 3014, Bern, Switzerland.
| |
Collapse
|
39
|
Zysset P, Pahr D, Engelke K, Genant HK, McClung MR, Kendler DL, Recknor C, Kinzl M, Schwiedrzik J, Museyko O, Wang A, Libanati C. Comparison of proximal femur and vertebral body strength improvements in the FREEDOM trial using an alternative finite element methodology. Bone 2015; 81:122-130. [PMID: 26141837 DOI: 10.1016/j.bone.2015.06.025] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/09/2015] [Revised: 06/23/2015] [Accepted: 06/29/2015] [Indexed: 01/15/2023]
Abstract
Denosumab reduced the incidence of new fractures in postmenopausal women with osteoporosis by 68% at the spine and 40% at the hip over 36 months compared with placebo in the FREEDOM study. This efficacy was supported by improvements from baseline in vertebral (18.2%) strength in axial compression and femoral (8.6%) strength in sideways fall configuration at 36 months, estimated in Newtons by an established voxel-based finite element (FE) methodology. Since FE analyses rely on the choice of meshes, material properties, and boundary conditions, the aim of this study was to independently confirm and compare the effects of denosumab on vertebral and femoral strength during the FREEDOM trial using an alternative smooth FE methodology. Unlike the previous FE study, effects on femoral strength in physiological stance configuration were also examined. QCT data for the proximal femur and two lumbar vertebrae were analyzed by smooth FE methodology at baseline, 12, 24, and 36 months for 51 treated (denosumab) and 47 control (placebo) subjects. QCT images were segmented and converted into smooth FE models to compute bone strength. L1 and L2 vertebral bodies were virtually loaded in axial compression and the proximal femora in both fall and stance configurations. Denosumab increased vertebral body strength by 10.8%, 14.0%, and 17.4% from baseline at 12, 24, and 36 months, respectively (p<0.0001). Denosumab also increased femoral strength in the fall configuration by 4.3%, 5.1%, and 7.2% from baseline at 12, 24, and 36 months, respectively (p<0.0001). Similar improvements were observed in the stance configuration with increases of 4.2%, 5.2%, and 5.2% from baseline (p≤0.0007). Differences between the increasing strengths with denosumab and the decreasing strengths with placebo were significant starting at 12 months (vertebral and femoral fall) or 24 months (femoral stance). Using an alternative smooth FE methodology, we confirmed the significant improvements in vertebral body and proximal femur strength previously observed with denosumab. Estimated increases in strength with denosumab and decreases with placebo were highly consistent between both FE techniques.
Collapse
Affiliation(s)
| | - Dieter Pahr
- Vienna University of Technology, Vienna, Austria
| | - Klaus Engelke
- University of Erlangen, Erlangen, Germany and Synarc Germany, Hamburg, Germany
| | | | | | | | | | | | | | - Oleg Museyko
- University of Erlangen-Nuremberg, Erlangen-Nuremberg, Germany
| | | | | |
Collapse
|
40
|
The influence of bone density and anisotropy in finite element models of distal radius fracture osteosynthesis: Evaluations and comparison to experiments. J Biomech 2015; 48:4116-4123. [PMID: 26542787 DOI: 10.1016/j.jbiomech.2015.10.012] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2015] [Revised: 10/07/2015] [Accepted: 10/10/2015] [Indexed: 11/23/2022]
Abstract
Continuum-level finite element (FE) models can be used to analyze and improve osteosynthesis procedures for distal radius fractures (DRF) from a biomechanical point of view. However, previous models oversimplified the bone material and lacked thorough experimental validation. The goal of this study was to assess the influence of local bone density and anisotropy in FE models of DRF osteosynthesis for predictions of axial stiffness, implant plate stresses, and screw loads. Experiments and FE analysis were conducted in 25 fresh frozen cadaveric radii with DRFs treated by volar locking plate osteosynthesis. Specimen specific geometries were captured using clinical quantitative CT (QCT) scans of the prepared samples. Local bone material properties were computed based on high resolution CT (HR-pQCT) scans of the intact radii. The axial stiffness and individual screw loads were evaluated in FE models, with (1) orthotropic inhomogeneous (OrthoInhom), (2) isotropic inhomogeneous (IsoInhom), and (3) isotropic homogeneous (IsoHom) bone material and compared to the experimental axial stiffness and screw-plate interface failures. FE simulated and experimental axial stiffness correlated significantly (p<0.0001) for all three model types. The coefficient of determination was similar for OrthoInhom (R(2)=0.807) and IsoInhom (R(2)=0.816) models but considerably lower for IsoHom models (R(2)=0.500). The peak screw loads were in qualitative agreement with experimental screw-plate interface failure. Individual loads and implant plate stresses of IsoHom models differed significantly (p<0.05) from OrthoInhom and IsoInhom models. In conclusion, including local bone density in FE models of DRF osteosynthesis is essential whereas local bone anisotropy hardly effects the models׳ predictive abilities.
Collapse
|