Hostname: page-component-586b7cd67f-t8hqh Total loading time: 0 Render date: 2024-11-22T14:42:20.997Z Has data issue: false hasContentIssue false

Micromechanics-based conversion of CT data into anisotropic elasticity tensors, applied to FE simulations of a mandible

Published online by Cambridge University Press:  15 March 2011

Christian Hellmich
Affiliation:
Vienna University of Technology (TU Wien), Institute for Mechanics of Materials and Structures, A-1040 Vienna, Austria
Cornelia Kober
Affiliation:
Hamburg University for Applied Sciences, Faculty of Life Sciences, 21033 Hamburg, Germany
Bodo Erdmann
Affiliation:
Zuse Institute, 14195 Berlin-Dahlem, Germany
Get access

Abstract

Computer Tomographic (CT) image data have become a standard basis for structural analyses of bony organs. In this context, regression functions between stiffness components and Hounsfields units (HU) from Computer Tomography, related to X-ray attenuation coefficients, are widely used for the definition of the (actually inhomogeneous and anisotropic) material behavior inside the organ. Herein, we suggest to derive the functional dependence of the fully orthotropic stiffness tensors on the Hounsfield units from the physical information contained in the X-ray attenuation coefficients: (i) Based on voxel average rules for the X-ray attenuation coefficients, we assign to each voxel the volume fraction occupied by water (marrow) and that occupied by solid bone matrix. (ii) By means of a continuum micromechanics representation for bone, which is based on voxel-invariant (species and whole bone-specific) stiffness properties of solid bone matrix and of water, we convert the aforementioned volume fractions into voxel-specific orthotropic stiffness tensor components. The micromechanics model, in combination with the average rule for X-ray attenuation coefficients, predicts a quasi-linear relationship between axial Young's modulus and HU, and highly nonlinear relationships for both circumferential and radial Young's moduli as well as for the shear moduli in all principal material directions. Corresponding whole-organ Finite Element analyses of a partially dentulous human mandible characterized by atrophy of the alveolar ridge show that volumetric strain concentrations/peaks within the organ are decreased when considering material anisotropy, and increased when considering material inhomogeneity.

Type
Research Article
Copyright
Copyright © Materials Research Society 2009

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

REFERENCES

1. Boccaccio, A., Lamberti, L., Pappalettere, C., Carano, A., and Cozzani, M., J. Biomech. 39, 2907 (2006)Google Scholar
2. Lotz, J.C., Cheal, E. J., and Hayes, W. C., J. Biomech. Eng. 113, 353 (1991)Google Scholar
3. Couteau, B., Hobatho, M.-C., Darmana, R., Brignola, J.-C., and Arlaud, J.-Y.. J. Biomech. 31, 383 (1998)Google Scholar
4. Crawley, E.O., E. O., Evans, W. D., and Owen, G. M., Phys. Med. Biol. 33, 1113 (1988)Google Scholar
5. Fritsch, A. and Hellmich, Ch., J. Theor. Biol. 244, 597 (2007)Google Scholar
6. Hellmich, Ch., Kober, C., Erdmann, B., Ann. Biomed. Eng. 36, 108 (2008)Google Scholar
7. Kober, C., C., Erdmann, B., Hellmich, Ch., Sader, R., and Zeilhofer, H.-F., Comput. Methods Biomech. Biomed. Eng. 9, 91 (2006)Google Scholar
8.NLM. The visible human project. National Library of Medicine 4 (2005), http://www.nlm.nih.gov/research/visible.Google Scholar