Highlights
- •Unified least-squares-based framework for fast “nodal” and “element” strategies of mapping scalar CT voxel densities onto finite element mesh.
- •Approximation of density by finite element method using least squares projection.
- •The discontinuous variant of the finite element method is more suitable for density approximation than the continuous variant in terms of accuracy and efficiency.
- •The discontinuous zero-order finite element method preserves the density spectrum better than the continuous one.
Abstract
Background
The spatially varying mechanical properties in finite element models of bone are most
often derived from bone density data obtained via quantitative computed tomography.
The key step is to accurately and efficiently map the density given in voxels to the
finite element mesh.
Methods
The density projection is first formulated in least-squares terms and then discretized
using a continuous and discontinuous variant of the finite element method. Both discretization
variants are compared with the nodal and element approaches known from the literature.
Findings
In terms of accuracy in the L2 norm, energy distance and efficiency, the discontinuous
zero-order variant appears to be the most advantageous. The proposed variant sufficiently
preserves the spectrum of density at the edges, while keeping computational cost low.
Interpretation
The continuous finite element method is analogous to the nodal formulation in the
literature, while the discontinuous finite element method is analogous to the element
formulation. The two variants differ in terms of implementation, computational cost
and ability to preserve the density spectrum. These differences cannot be described
and measured by known indirect methods from the literature.
Keywords
To read this article in full you will need to make a payment
Purchase one-time access:
Academic & Personal: 24 hour online accessCorporate R&D Professionals: 24 hour online accessOne-time access price info
- For academic or personal research use, select 'Academic and Personal'
- For corporate R&D use, select 'Corporate R&D Professionals'
Subscribe:
Subscribe to Clinical BiomechanicsAlready a print subscriber? Claim online access
Already an online subscriber? Sign in
Register: Create an account
Institutional Access: Sign in to ScienceDirect
References
- Opportunistic application of phantom-less calibration methods for fracture risk prediction using QCT/FEA.Eur. Radiol. 2021; 31: 9428-9435https://doi.org/10.1007/s00330-021-08071-w
- Advancements in osteoporosis imaging, screening, and study of disease etiology.Curr. Osteoporosis Rep. 2021; 19: 532-541https://doi.org/10.1007/s11914-021-00699-3
- 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-689https://doi.org/10.1016/j.medengphy.2016.03.006
- CT-based internal density calibration for opportunistic skeletal assessment using abdominal CT scans.Med. Eng. Phys. 2020; 78 (ISSN 1350-4533): 55-63https://doi.org/10.1016/j.medengphy.2020.01.009
- Empirical relationships between bone density and ultimate strength: a literature review.J. Mech. Behavior Biomed. Mat. 2020; 110 (ISSN 1751-6161): 103866https://doi.org/10.1016/j.jmbbm.2020.103866
- Risk assessment of vertebral compressive fracture using bone mass index and strength predicted by computed tomography image based finite element analysis.Clin. Biomech. 2021; 85 (ISSN 0268-0033): 105365https://doi.org/10.1016/j.clinbiomech.2021.105365
- Prediction of mechanical properties of trabecular bone in patients with type 2 diabetes using damage based finite element method.J. Biomech. 2021; 123 (ISSN 0021-9290): 110495https://doi.org/10.1016/j.jbiomech.2021.110495
- Finite element analysis informed variable selection for femoral fracture risk prediction.J. Mech. Behavior Biomed. Mat. 2021; 118 (ISSN 1751-6161): 104434https://doi.org/10.1016/j.jmbbm.2021.104434
- Finite element analysis of fracture fixation.Curr. Osteoporosis Rep. 2021; 19: 403-416https://doi.org/10.1007/s11914-021-00690-y
- Biomechanics of subtrochanteric fracture fixation using short cephalomedullary nails: a finite element analysis.PLOS One. 2021; 16: 1-15https://doi.org/10.1371/journal.pone.0253862
- An improved method for the automatic mapping of computed tomography numbers onto finite element models.Med. Eng. Phys. 2004; 26 (ISSN 1350-4533): 61-69https://doi.org/10.1016/S1350-4533(03)00138-3
- A modified method for assigning material properties to FE models of bones, medical engineering.Physics. 2008; 30 (ISSN 1350-4533): 444-453https://doi.org/10.1016/j.medengphy.2007.05.006
- A new approach for assigning bone material properties from CT images into finite element models.J. Biomech. 2010; 43 (ISSN 0021-9290): 1011-1015https://doi.org/10.1016/j.jbiomech.2009.10.040
- Comparisons of node-based and element-based approaches of assigning bone material properties onto subject-specific finite element models.Med. Eng. Phys. 2015; 37 (ISSN 1350-4533): 808-812https://doi.org/10.1016/j.medengphy.2015.05.006
- An open source software tool to assign the material properties of bone for ABAQUS finite element simulations.J. Biomech. 2016; 49 (ISSN 0021-9290): 3116-3121https://doi.org/10.1016/j.jbiomech.2016.07.037
- An exclusion approach for addressing partial volume artifacts with quantitative computed tomography-based finite element modeling of the proximal tibia.Med. Eng. Phys. 2020; 76 (ISSN 1350-4533): 95-100https://doi.org/10.1016/j.medengphy.2019.10.013
- Heterogeneous material mapping methods for patient-specific finite element models of pelvic trabecular bone: a convergence study.Med. Eng. Phys. 2021; 96 (ISSN 1350-4533): 1-12https://doi.org/10.1016/j.medengphy.2021.07.012
- Multidimensional binary search trees used for associative searching.CACM. 1975; 18: 509-517https://doi.org/10.1145/361002.361007
- Automated Solution of Differential Equations by the Finite Element Method: The FEniC book. 84. Springer Science & Business Media, Heidelberg2012https://doi.org/10.1007/978-3-642-23099-8
- Gmsh, urlprefix.2021
- Efficient unstructured mesh generation by means of delaunay triangulation and Bowyer-Watson algorithm.J. Comput. Phys. 1993; 106 (ISSN 0021-9991): 125-138https://doi.org/10.1006/jcph.1993.1097
- Energy statistics: a class of statistics based on distances.J. Stat. Plan. Inference. 2013; 143: 1249-1272https://doi.org/10.1016/j.jspi.2013.03.018
- SciPy 1.0 contributors, scipy 1.0: fundamental algorithms for scientific computing in Python.Nat. Methods. 2020; 17: 261-272https://doi.org/10.1038/s41592-019-0686-2
- The Hounsfield value for cortical bone geometry in the proximal humerus-an in vitro study.Skeletal Radiol. 2012; 41: 557-568https://doi.org/10.1007/s00256-011-1255-7
- Determination of the Hounsfield value for CT-based design of custom femoral stems.J. Bone Joint Surgery. Brit. Volume. 1999; 81-B: 143-147https://doi.org/10.1302/0301-620x.81b1.0810143
- On the H1-stability of the L2-projection onto finite element spaces.Numer. Math. 2014; 126 (ISSN 0029-599X 0945-3245): 361-381https://doi.org/10.1007/s00211-013-0562-4
- FEBio: finite elements for biomechanics.J. Biomech. Eng. 2001; 134 (011005 (10 pages))https://doi.org/10.1115/1.4005694
- Visualizing with VTK: a tutorial.IEEE Comput. Graph. Appl. 2000; 20 (ISSN 0272-1716 1558-1756): 20-27https://doi.org/10.1109/38.865875
- An overview of the HDF5 technology suite and its applications.in: Proceedings of the EDBT/ICDT 2011 Workshop on Array Databases, AD ‘11, Association for Computing Machinery, New York, NY USA2011: 36-47https://doi.org/10.1145/1966895.1966900 (2011. ISBN 9781450306140)
- A computational framework for crack propagation in spatially heterogeneous materials.Philosoph. Trans. Royal Soc. A. 2021; 379: 20200291https://doi.org/10.1098/rsta.2020.0291
- QCT-FE modeling of the proximal tibia: effect of mapping strategy on convergence time and model accuracy.Med. Eng. Phys. 2021; 88 (ISSN 1350-4533): 41-46https://doi.org/10.1016/j.medengphy.2020.12.006
- Periodic table of the finite elements.Siam News. 2014; 47: 212
Article Info
Publication History
Published online: June 14, 2022
Accepted:
June 7,
2022
Received in revised form:
March 23,
2022
Received:
November 19,
2021
Identification
Copyright
© 2022 Elsevier Ltd. All rights reserved.