Technical Brief

Elastic Anisotropy of Trabecular Bone in the Elderly Human Vertebra

[+] Author and Article Information
Ginu U. Unnikrishnan, John A. Gallagher, Amira I. Hussein

Orthopaedic and Developmental Biomechanics Laboratory,
Department of Mechanical Engineering,
Boston University,
Boston, MA 02215

Glenn D. Barest

Department of Radiology,
Boston University,
Boston, MA 02118

Elise F. Morgan

Orthopaedic and Developmental Biomechanics Laboratory,
Department of Mechanical Engineering,
Boston University,
Boston, MA 02215;
Department of Biomedical Engineering,
Boston University,
Boston, MA 02215;
Orthopaedic Surgery,
Boston University,
Boston, MA 02118
e-mail: efmorgan@bu.edu

Manuscript received October 9, 2014; final manuscript received July 31, 2015; published online October 1, 2015. Assoc. Editor: Joel D Stitzel.

J Biomech Eng 137(11), 114503 (Oct 01, 2015) (6 pages) Paper No: BIO-14-1504; doi: 10.1115/1.4031415 History: Received October 09, 2014; Revised July 31, 2015

Knowledge of the nature of the elastic symmetry of trabecular bone is fundamental to the study of bone adaptation and failure. Previous studies have classified human vertebral trabecular bone as orthotropic or transversely isotropic but have typically obtained samples from only selected regions of the centrum. In this study, the elastic symmetry of human vertebral trabecular bone was characterized using microfinite element (μFE) analyses performed on 1019 cubic regions of side length equal to 5 mm, obtained via thorough sampling of the centrums of 18 human L1 vertebrae (age = 81.17 ± 7.7 yr; eight males and ten females). An optimization procedure was used to find the closest orthotropic representation of the resulting stiffness tensor for each cube. The orthotropic elastic constants and orientation of the principal elastic axes were then recorded for each cube and were compared to the constants predicted from Cowin's fabric-based constitutive model (Cowin, 1985, “The Relationship Between the Elasticity Tensor and the Fabric Tensor,” Mech. Mater., 4(2), pp. 137–147.) and the orientation of the principal axes of the fabric tensor, respectively. Deviations from orthotropy were quantified by the “orthotropic error” (van Rietbergen et al., 1996, “Direct Mechanics Assessment of Elastic Symmetries and Properties of Trabecular Bone Architecture,” J. Biomech., 29(12), pp. 1653–1657), and deviations from transverse isotropy were determined by statistical comparison of the secondary and tertiary elastic moduli. The orthotropic error was greater than 50% for nearly half of the cubes, and the secondary and tertiary moduli differed from one another (p < 0.0001). Both the orthotropic error and the difference between secondary and tertiary moduli decreased with increasing bone volume fraction (BV/TV; p ≤ 0.007). Considering only the cubes with an orthotropic error less than 50%, only moderate correlations were observed between the fabric-based and the μFE-computed elastic moduli (R2 ≥ 0.337; p < 0.0001). These results indicate that when using a criterion of 5 mm for a representative volume element (RVE), transverse isotropy or orthotropy cannot be assumed for elderly human vertebral trabecular bone. Particularly at low values of BV/TV, this criterion does not ensure applicability of theories of continuous media. In light of the very sparse and inhomogeneous microstructure found in the specimens analyzed in this study, further work is needed to establish guidelines for selecting a RVE within the aged vertebral centrum.

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Fig. 1

Schematic representation of the procedure adopted to find the best orthotropic representation of each cube of trabecular bone: All δij terms in COPT (d) are set to zero to obtain CORT (e); no other changes are made when obtaining CORT from COPT. (a) μCT image of vertebra, (b) micro-FE analyses of the cube, (c) trabecular fabric-Eigenvectors of the cube, (d) nearly orthotropic stiffness matrix after optimization, and (e) orthotropic stiffness matrix after assigning nonorthotropic terms of zero.

Grahic Jump Location
Fig. 2

Histogram of orthotropic error (a) and correlations between orthotropic error and BV/TV (b) and between E2/E1 and BV/TV (c)

Grahic Jump Location
Fig. 3

Correlations of the Young's moduli ((a)–(c)), shear moduli ((d)–(f)), and Poisson's ratio ((g)–(i)) from the μFE model versus the morphology-based model: The correlation coefficients were nearly unchanged when the circled data point was excluded




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In