On the limits of finite element models created from (micro)CT datasets and used in studies of bone-implant-related biomechanical problems

Marcián, Petr; Borák, Libor; Zikmund, Tomáš; Horáčková, Ladislava; Kaiser, Jozef; Joukal, Marek; Wolff, Jan

Affiliations

¹Institute of Solid Mechanics, Mechatronics and Biomechanics, Faculty of Mechanical Engineering, Brno University of Technology, Brno, Czech Republic

²CEITEC - Central European Institute of Technology, Brno University of Technology, Czech Republic

³Department of Anatomy, Faculty of Medicine, Masaryk University, Brno, Czech Republic

⁴Department of Oral and Maxillofacial Surgery, Division for Regenerative Orofacial Medicine, University Hospital Hamburg-Eppendorf, Hamburg, Germany

⁵Fraunhofer Research Institution for Additive Manufacturing Technologies IAPT, Hamburg, Germany

Abstract and keywords

Patient-specific approach is gaining a wide popularity in computational simulations of biomechanical systems. Simulations (most often based on the finite element method) are to date routinely created using data from imaging devices such as computed tomography which makes the models seemingly very complex and sophisticated. However, using a computed tomography in finite element calculations does not necessarily enhance the quality or even credibility of the models as these depend on the quality of the input images. Low-resolution (medical-)CT datasets do not always offer detailed representation of trabecular bone in FE models and thus might lead to incorrect calculation of mechanical response to external loading. The effect of image resolution on mechanical simulations of bone-implant interaction has not been thoroughly studied yet. In this study, the effect of image resolution on the modeling procedure and resulting mechanical strains in bone was analyzed on the example of cranial implant. For this purpose, several finite element models of bone interacting with fixation-screws were generated using seven computed tomography datasets of a bone specimen but with different image resolutions (ranging from micro-CT resolution of 25 μm to medical-CT resolution of 1250 μm). The comparative analysis revealed that FE models created from images of low resolution (obtained from medical computed tomography) can produce biased results. There are two main reasons: 1. Medical computed tomography images do not allow generating models with complex trabecular architecture which leads to substituting of the intertrabecular pores with a fictitious mass; 2. Image gray value distribution can be distorted resulting in incorrect mechanical properties of the bone and thus in unrealistic or even completely fictitious mechanical strains. The biased results of calculated mechanical strains can lead to incorrect conclusion, especially when bone-implant interaction is investigated. The image resolution was observed not to significantly affect stresses in the fixation screw itself; however, selection of bone material representation might result in significantly different stresses in the screw.

Keywords: Finite element method; Computational modeling; Computed tomography; Image resolution; Bone tissue; Mechanical strain

Links

DOI: 10.1016/j.jmbbm.2021.104393

PubMed: 33647729

WoS: 000635269700003

History

Received: 2020-07-23

Revised: 2021-02-10

Accepted: 2021-02-10

Online: 2021-02-18

Cited Works (13)

Year	Entry
1999	Turner CH, Rho J, Takano Y, Tsui TY, Pharr GM. The elastic properties of trabecular and cortical bone tissues are similar: results from two microscopic measurement techniques. J Biomech. April 1999;32(4):437-441.
1991	Snyder SM, Schneider E. Estimation of mechanical properties of cortical bone by computed tomography. J Orthop Res. May 1991;9(3):422-431.
1995	van Rietbergen B, Weinans H, Huiskes R, Odgaard A. A new method to determine trabecular bone elastic properties and loading using micromechanical finite-element models. J Biomech. January 1995;28(1):69-81.
2012	Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. July 2012;9(7):671-675.
2009	Motherway JA, Verschueren P, Van der Perre G, Vander Sloten J, Gilchrist MD. The mechanical properties of cranial bone: the effect of loading rate and cranial sampling position. J Biomech. September 18, 2009;42(13):2129-2135.
1997	Rho J-Y, Tsui TY, Pharr GM. Elastic properties of human cortical and trabecular lamellar bone measured by nanoindentation. Biomaterials. 1997;18(20):1325-1330.
2009	Yoganandan N, Pintar FA, Zhang J, Baisden JL. Physical properties of the human head: mass, center of gravity and moment of inertia. J Biomech. June 19, 2009;42(9):1177-1192.
1994	Keller TS. Predicting the compressive mechanical behavior of bone. J Biomech. September 1994;27(9):1159-1168.
2017	Costa MC, Tozzi G, Cristofolini L, Danesi V, Viceconti M, Dall’Ara E. Micro finite element models of the vertebral body: validation of local displacement predictions. PLoS One. July 11, 2017;12(7):e0180151.
2004	Bayraktar HH, Morgan EF, Niebur GL, Morris GE, Wong EK, Keaveny TM. Comparison of the elastic and yield properties of human femoral trabecular and cortical bone tissue. J Biomech. January 2004;37(1):27-35.
2012	Wirth AJ, Müller R, van Lenthe GH. The discrete nature of trabecular bone microarchitecture affects implant stability. J Biomech. April 5, 2012;45(6):1060-1067.
2005	Shefelbine SJ, Simon U, Claes L, Gold A, Gabet Y, Bab I, Müller R, Augat P. Prediction of fracture callus mechanical properties using micro-CT images and voxel-based finite element analysis. Bone. March 2005;36(3):480-488.
2008	Schileo E, Taddei F, Cristofolini L, Viceconti M. Subject-specific finite element models implementing a maximum principal strain criterion are able to estimate failure risk and fracture location on human femurs tested in vitro. J Biomech. 2008;41(2):356-367.

Cited By (1)

Year	Entry
2022	Ortner L. Comparison of Peri-Implant Loading in Screw-Bone Constructs Predicted by Homogenized and Micro FE Models [Master's thesis]. Vienna University of Technology; May 11, 2022.