Anisotropic computational modelling of bony structures from CT data:​ an almost automatic procedure

Toniolo, Ilaria<sup>1,2</sup>; Salmaso, Claudia<sup>2</sup>; Bruno, Giovanni<sup>3</sup>; De Stefani, Alberto<sup>3</sup>; Stefanini, Cesare<sup>1</sup>; Gracco, Antonio Luigi Tiberio<sup>3</sup>; Carniel, Emanuele Luigi<sup>2,4</sup>

Affiliations

¹The BioRobotics Institute, Scuola Superiore Sant’Anna, Pisa, Italy

²Department of Industrial Engineering, University of Padova, Italy

³Department of Neurosciences, University of Padova, Italy

⁴Centre for Mechanics of Biological Materials, University of Padova, Italy

Abstract and keywords

Background and objective: The use of modelling techniques that combine CT data and bone tissue micromechanics is spreading in computational biomechanics. Finite Element models show great potential in surgical planning of intervention and in prediction of stress and strain fields through a non-invasive method. The main challenge pertains to the reliable characterization of bone mechanical behaviour. An almost automatic procedure is here defined, which provides computational models of bony structures considering the actual anisotropy of bone tissue response. The innovative aspect resides on the automatic detection of the directions of anisotropy as the eigenvectors of a three-dimensional distribution matrix of HU values.

Methods: The procedure combines CT data and micromechanics modelling techniques. Regarding a specific location, the procedure reports both the orthotropic elastic constants, by the analysis of the local HU value, and the anisotropic material directions, by the analysis of the HU values distribution around the specific location.

Results: The procedure returns the distribution of bone tissue orthotropic elasticity tensor. The procedure proves to correctly respect the differentiation between cortical and trabecular bone. Principal directions show to be consistent with experimental data from ultrasound measurements. Regarding the material mapping from voxel to FE model, the developed strategies show to be reliable, leading to marginal errors (lower than 10%) for most of CT voxels (more than 90%). The computational analyses of typical structural loading conditions lead to strain values that are comparable with results from strain gauges experimentations. The development and the exploitation of FE models of different bony structures allow assessing the reliability of the procedure for cortical bone.

Conclusions: The results highlight the potentialities of the procedure in providing accurate patientspecific biomechanical models of bony structures starting from CT data. The accuracy and the automatism of the procedure are important factors for the development of real time clinical tools. The main limitations of this work remain the not fully automatism and the reliability assessment, which is based mainly on cortical bone regions only

Keywords: Bone mechanics; Orthotropic elasticity; CT data; Micromechanical modelling

Links

DOI: 10.1016/j.cmpb.2020.105319

PubMed: 31951872

WoS: 000533562800014

History

Received: 2019-07-31

Revised: 2019-11-27

Accepted: 2020-01-05

Cited Works (12)

Year	Entry
1975	Reilly DT, Burstein AH. The elastic and ultimate properties of compact bone tissue. J Biomech. 1975;8(6):393-405.
2009	Espinoza Orías AA, Deuerling JM, Landrigan MD, Renaud JE, Roeder RK. Anatomic variation in the elastic anisotropy of cortical bone tissue in the human femur. J Mech Behav Biomed Mater. July 2009;2(3):255-263.
2007	Taddei F, Schileo E, Helgason B, Cristofolini L, Viceconti M. The material mapping strategy influences the accuracy of ct-based finite element models of bones: an evaluation against experimental measurements. Med Eng Phys. 2007;29(9):973-979.
2006	Taddei F, Cristofolini L, Martelli S, Gill HS, Viceconti M. Subject-specific finite element models of long bones: an in vitro evaluation of the overall accuracy. J Biomech. 2006;39(13):2457-2467.
2004	Taddei F, Pancanti A, Viceconti M. An improved method for the automatic mapping of computed tomography numbers onto finite element models. Med Eng Phys. 2004;26(1):61-69.
2011	Trabelsi N, Yosibash Z. Patient-specific finite-element analyses of the proximal femur with orthotropic material properties validated by experiments. J Biomech Eng. June 2011;133(6):061001.
2001	Morgan EF, Keaveny TM. Dependence of yield strain of human trabecular bone on anatomic site. J Biomech. 2001;34(5):569-577.
2016	Maggiano IS, Maggiano CM, Clement JG, Thomas CDL, Carter Y, Cooper DML. Three‐dimensional reconstruction of Haversian systems in human cortical bone using synchrotron radiation‐based micro‐CT: morphology and quantification of branching and transverse connections across age. J Anat. May 2016;228(5):719-732.
1984	Ashman RB, Cowin SC, Van Buskirk WC, Rice JC. A continuous wave technique for the measurement of the elastic properties of cortical bone. J Biomech. 1984;17(5):349-361.
2011	Rathnayaka K, Sahama T, Schuetz MA, Schmutz B. Effects of CT image segmentation methods on the accuracy of long bone 3D reconstructions. Med Eng Phys. March 2011;33(2):226-233.
2009	Vasilić B, Rajapakse CS, Wehrli FW. Classification of trabeculae into three‐dimensional rodlike and platelike structures via local inertial anisotropy. Med Phys. July 2009;36(7):3280-3291.
2004	Viceconti M, Davinelli M, Taddei F, Cappello A. Automatic generation of accurate subject-specific bone finite element models to be used in clinical studies. J Biomech. 2004;37(10):1597-1605.

Cited By (1)

Year	Entry
2021	Sabet FA, Koric S, Idkaidek A, Jasiuk I. High-performance computing comparison of implicit and explicit nonlinear finite element simulations of trabecular bone. Comput Methods Prog Biomed. March 2021;200:105870.