Research Papers

Dependence of Anisotropy of Human Lumbar Vertebral Trabecular Bone on Quantitative Computed Tomography-Based Apparent Density

[+] Author and Article Information
Ameet K. Aiyangar

Laboratory for Mechanical
Systems Engineering (304),
EMPA (Swiss Federal Laboratories for Materials
Science and Technology) Dübendorf,
Überlandstrasse 129,
CH-8600 Dübendorf,
Zürich 8600, Switzerland
Department of Mechanical Engineering,
University of Wisconsin-Madison,
Madison, WI 53706
e-mail: ameetaiyangar@gmail.com

Juan Vivanco

Department of Mechanical Engineering,
University of Wisconsin-Madison,
Madison, WI 53706
Facultad de Ingeniería y Ciencias,
Universidad Adolfo Ibáñez,
Viña del Mar, Chile

Anthony G. Au

VibeDx Diagnostic Corp.,
Edmonton, AB T5J4P6, Canada

Paul A. Anderson

Department of Orthopedics and Rehabilitation,
University of Wisconsin-Madison,
Madison, WI 53705

Everett L. Smith

Department of Population Health Sciences,
University of Wisconsin-Madison,
Madison, WI 53706

Heidi-Lynn Ploeg

Department of Mechanical Engineering,
University of Wisconsin-Madison,
Madison, WI 53706

1Corresponding author.

Manuscript received September 25, 2013; final manuscript received May 2, 2014; accepted manuscript posted May 14, 2014; published online July 3, 2014. Assoc. Editor: Brian D. Stemper.

J Biomech Eng 136(9), 091003 (Jul 03, 2014) (10 pages) Paper No: BIO-13-1448; doi: 10.1115/1.4027663 History: Received September 25, 2013; Revised May 02, 2014; Accepted May 14, 2014

Most studies investigating human lumbar vertebral trabecular bone (HVTB) mechanical property–density relationships have presented results for the superior–inferior (SI), or “on-axis” direction. Equivalent, directly measured data from mechanical testing in the transverse (TR) direction are sparse and quantitative computed tomography (QCT) density-dependent variations in the anisotropy ratio of HVTB have not been adequately studied. The current study aimed to investigate the dependence of HVTB mechanical anisotropy ratio on QCT density by quantifying the empirical relationships between QCT-based apparent density of HVTB and its apparent compressive mechanical properties— elastic modulus (Eapp), yield strength (σy), and yield strain (εy)—in the SI and TR directions for future clinical QCT-based continuum finite element modeling of HVTB. A total of 51 cylindrical cores (33 axial and 18 transverse) were extracted from four L1 human lumbar cadaveric vertebrae. Intact vertebrae were scanned in a clinical resolution computed tomography (CT) scanner prior to specimen extraction to obtain QCT density, ρCT. Additionally, physically measured apparent density, computed as ash weight over wet, bulk volume, ρapp, showed significant correlation with ρCTCT = 1.0568 × ρapp, r = 0.86]. Specimens were compression tested at room temperature using the Zetos bone loading and bioreactor system. Apparent elastic modulus (Eapp) and yield strength (σy) were linearly related to the ρCT in the axial direction [ESI = 1493.8 × (ρCT), r = 0.77, p < 0.01; σY,SI = 6.9 × (ρCT) − 0.13, r = 0.76, p < 0.01] while a power-law relation provided the best fit in the transverse direction [ETR = 3349.1 × (ρCT)1.94, r = 0.89, p < 0.01; σY,TR = 18.81 × (ρCT)1.83, r = 0.83, p < 0.01]. No significant correlation was found between εy and ρCT in either direction. Eapp and σy in the axial direction were larger compared to the transverse direction by a factor of 3.2 and 2.3, respectively, on average. Furthermore, the degree of anisotropy decreased with increasing density. Comparatively, εy exhibited only a mild, but statistically significant anisotropy: transverse strains were larger than those in the axial direction by 30%, on average. Ability to map apparent mechanical properties in the transverse direction, in addition to the axial direction, from CT-based densitometric measures allows incorporation of transverse properties in finite element models based on clinical CT data, partially offsetting the inability of continuum models to accurately represent trabecular architectural variations.

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


Crawford, R. P., Cann, C. E., and Keaveny, T. M., 2003, “Finite Element Models Predict in Vitro Vertebral Body Compressive Strength Better Than Quantitative Computed Tomography,” Bone, 33(4), pp. 744–750. [CrossRef] [PubMed]
Mosekilde, L., Mosekilde, L., and Danielsen, C. C., 1987, “Biomechanical Competence of Vertebral Trabecular Bone in Relation to Ash Density and Age in Normal Individuals,” Bone, 8(2), pp. 79–85. [CrossRef] [PubMed]
Nicholson, P. H. F., Cheng, X. G., Lowet, G., Boonen, S., Davie, M. W. J., Dequeker, J., and Van derPerre, G., 1997, “Structural and Material Mechanical Properties of Human Vertebral Cancellous Bone,” Med. Eng. Phys., 19(8), pp. 729–737. [CrossRef] [PubMed]
Banse, X., Devogelaer, J. P., Munting, E., Delloye, C., Cornu, O., and Grynpas, M., 2001, “Inhomogeneity of Human Vertebral Cancellous Bone: Systematic Density and Structure Patterns Inside the Vertebral Body,” Bone, 28(5), pp. 563–571. [CrossRef] [PubMed]
Hulme, P. A., Boyd, S. K., and Ferguson, S. J., 2007, “Regional Variation in Vertebral Bone Morphology and Its Contribution to Vertebral Fracture Strength,” Bone, 41(6), pp. 946–957. [CrossRef] [PubMed]
Wegrzyn, J., Roux, J. P., Arlot, M. E., Boutroy, S., Vilayphiou, N., Guyen, O., Delmas, P. D., Chapurlat, R., and Bouxsein, M. L., 2010, “Role of Trabecular Microarchitecture and Its Heterogeneity Parameters in the Mechanical Behavior of Ex Vivo Human L3 Vertebrae,” J. Bone Miner. Res., 25(11), pp. 2324–2331. [CrossRef] [PubMed]
Nazarian, A., Stauber, M., Zurakowski, D., Snyder, B. D., and Muller, R., 2006, “The Interaction of Microstructure and Volume Fraction in Predicting Failure in Cancellous Bone,” Bone, 39(6), pp. 1196–1202. [CrossRef] [PubMed]
Matsuura, M., Eckstein, F., Lochmuller, E. M., and Zysset, P. K., 2008, “The Role of Fabric in the Quasi-Static Compressive Mechanical Properties of Human Trabecular Bone From Various Anatomical Locations,” Biomech. Model Mechanobiol., 7(1), pp. 27–42. [CrossRef] [PubMed]
Chevalier, Y., Charlebois, M., Pahra, D., Varga, P., Heini, P., Schneider, E., and Zysset, P., 2008, “A Patient-Specific Finite Element Methodology to Predict Damage Accumulation in Vertebral Bodies Under Axial Compression, Sagittal Flexion and Combined Loads,” Comput. Methods Biomech. Biomed. Eng., 11(5), pp. 477–487. [CrossRef]
Chevalier, Y., and Zysset, P. K., 2012, “A Patient-Specific Computer Tomography-Based Finite Element Methodology to Calculate the Six Dimensional Stiffness Matrix of Human Vertebral Bodies,” ASME J. Biomech. Eng., 134(5), p. 051006. [CrossRef]
Mosekilde, L., Ebbeson, E. N., Tornvig, L., and Thomson, J. S., 2000, “Trabecular Bone Structure and Strength -Remodelling and Repair,” J. Musculoskeletal Neuronal Interact., 1, pp. 25–30.
Keller, T. S., 1994, “Predicting the Compressive Mechanical Behavior of Bone,” J. Biomech., 27(9), pp. 1159–1168. [CrossRef] [PubMed]
Keaveny, T. M., Pinilla, T. P., Crawford, R. P., Kopperdahl, D. L., and Lou, A., 1997, “Systematic and Random Errors in Compression Testing of Trabecular Bone,” J. Orthop. Res., 15(1), pp. 101–110. [CrossRef] [PubMed]
Kopperdahl, D. L., and Keaveny, T. M., 1998, “Yield Strain Behavior of Trabecular Bone,” J. Biomech., 31(7), pp. 601–608. [CrossRef] [PubMed]
Morgan, E. F., and Keaveny, T. M., 2001, “Dependence of Yield Strain of Human Trabecular Bone on Anatomic Site,” J. Biomech., 34(5), pp. 569–577. [CrossRef] [PubMed]
Ashman, R. B., Corin, J. D., and Turner, C. H., 1987, “Elastic Properties of Cancellous Bone: Measurement by an Ultrasonic Technique,” J. Biomech., 20(10), pp. 979–983, 985–986. [CrossRef] [PubMed]
Turner, C. H., Anne, V., and Pidaparti, R. M. V., 1997, “A Uniform Strain Criterion for Trabecular Bone Adaptation: Do Continuum-Level Strain Gradients Drive Adaptation?,” J. Biomech., 30(6), pp. 555–563. [CrossRef] [PubMed]
Currey, J. D., 1984, The Mechanical Adaptations of Bones, Princeton University Press Princeton, NJ.
Silva, M. J., Keaveny, T. M., and Hayes, W. C., 1998, “Computed Tomography-Based Finite Element Analysis Predicts Failure Loads and Fracture Patterns for Vertebral Sections,” J. Orthop. Res., 16(3), pp. 300–308. [CrossRef] [PubMed]
Kopperdahl, D. L., Roberts, A. D., and Keaveny, T. M., 1999, “Localized Damage in Vertebral Bone is Most Detrimental in Regions of High Strain Energy Density,” ASME J. Biomech. Eng., 121(6), pp. 622–628. [CrossRef]
Wachtel, E. F., and Keaveny, T. M., 1997, “Dependence of Trabecular Damage on Mechanical Strain,” J. Orthop. Res., 15(5), pp. 781–787. [CrossRef] [PubMed]
Keaveny, T. M., Morgan, E. F., Niebur, G. L., and Yeh, O. C., 2001, “Biomechanics of Trabecular Bone,” Annu. Rev. Biomed. Eng., 3(1), pp. 307–333. [CrossRef] [PubMed]
Cowin, S. C., 1986, “Wolff's Law of Trabecular Architecture at Remodeling Equilibrium,” ASME J. Biomech. Eng., 108(1), pp. 83–88. [CrossRef]
Chang, W. C. W., Christensen, T. M., Pinilla, T. P., and Keaveny, T. M., 1999, “Uniaxial Yield Strains for Bovine Trabecular Bone Are Isotropic and Asymmetric,” J. Orthop. Res., 17(4), pp. 582–585. [CrossRef] [PubMed]
Turner, C. H., 1989, “Yield Behavior of Bovine Cancellous Bone,” ASME J. Biomech. Eng., 111(3), pp. 256–260. [CrossRef]
Bevill, G., Farhamand, F., and Keaveny, T. M., 2009, “Heterogeneity of Yield Strain in Low-Density Versus High-Density Human Trabecular Bone,” J. Biomech., 42(13), pp. 2165–2170. [CrossRef] [PubMed]
Jones, D. B., Broeckmann,E., Pohl, T., and Smith, E. L., 2003, “Development of a Mechanical Testing and Loading System for Trabecular Bone Studies for Long Term Culture,” Eur. Cells Mater., 15, pp. 48–60.
Smith, E. L., and Jones, D., 2001, “Combined Perfusion and Mechanical Loading System for Explanted Bone,” U.S. Patent No. US 6,171,812 B1.
Davies, C. M., Jones, D. B., Stoddart, M. J., Koller, K., Smith, E. L., Archer, C. W., and Richards, R. G., 2006, “Mechanically Loaded Ex Vivo Bone Culture System ‘Zetos’: Systems and Culture Preparation,” Eur. Cells Mater., 11, pp. 57–75.
Garcia-Rodriguez, S., Smith, E. L., and Ploeg, H.-L., 2008, “A Calibration Procedure for a Bone Loading System,” J. Med. Devices, 2(1), pp. 1–6. [CrossRef]
Endres, S., Kratz, M., Wunsch, S., and Jones, D. B., 2009, “Zetos: A Culture Loading System for Trabecular Bone. Investigation of Different Loading Signal Intensities on Bovine Bone Cylinders,” J. Musculoskeletal Neuronal Interact., 9(3), pp. 173–183.
Vivanco, J., Garcia, S., Ploeg, H. L., Alvarez, G., Cullen, D., and Smith, E. L., 2013, “Apparent Elastic Modulus of Ex Vivo Trabecular Bovine Bone Increases With Dynamic Loading,” Proc. Inst. Mech. Eng., Part H, 227(8), pp. 904–912. [CrossRef]
Dumas, V., Perrier, A., Malaval, L., Laroche, N., Guignandon, A., Vico, L., and Rattner, A., 2009, “The Effect of Dual Frequency Cyclic Compression on Matrix Deposition by Osteoblast-Like Cells Grown in 3D Scaffolds and on Modulation of VEGF Variant Expression,” Biomaterials, 30(19), pp. 3279–3288. [CrossRef] [PubMed]
Mann, V., Huber, C., Kogianni, G., Jones, D., and Noble, B., 2006, “The Influence of Mechanical Stimulation on Osteocyte Apoptosis and Bone Viability in Human Trabecular Bone,” J. Musculoskeletal Neuronal Interact., 6(4), pp. 408–417.
Harrigan, T. P., Jasty, M., Mann, R. W., and Harris, W. H., 1988, “Limitations of the Continuum Assumption in Cancellous Bone,” J. Biomechanics, 21(4), pp. 269–275. [CrossRef]
Morgan, E. F., Yeh, O. C., Chang, W. C., and Keaveny, T. M., 2001, “Nonlinear Behavior of Trabecular Bone at Small Strains,” ASME J. Biomech. Eng., 123(1), pp. 1–9. [CrossRef]
Zysset, P. K., 2003, “A Review of Morphology-Elasticity Relationships in Human Trabecular Bone: Theories and Experiments,” J. Biomech., 36(10), pp. 1469–1485. [CrossRef] [PubMed]
Rincon-Kohli, L., and Zysset, P. K., 2009, “Multi-Axial Mechanical Properties of Human Trabecular Bone,” Biomech. Model Mechanobiol., 8(3), pp. 195–208. [CrossRef] [PubMed]
Fields, A. J., and Keaveny, T. M., 2012, “Trabecular Architecture and Vertebral Fragility in Osteoporosis,” Curr. Osteoporosis Rep., 10(2), pp. 132–140. [CrossRef]
Rho, J. Y., Ashman, R. B., and Turner, C. H., 1993, “Young's Modulus of Trabecular and Cortical Bone Material: Ultrasonic and Microtensile Measurements,” J. Biomech., 26(2), pp. 111–119. [CrossRef] [PubMed]
Keaveny, T. M., Wachtel, E. F., Zadesky, S. P., and Arramon, Y. P., 1999, “Application of the Tsai–Wu Quadratic Multiaxial Failure Criterion to Bovine Trabecular Bone,” ASME J. Biomech. Eng., 121(1), pp. 99–107. [CrossRef]
Gibson, L. J., 1985, “The Mechanical Behaviour of Cancellous Bone,” J. Biomech., 18(5), pp. 317–328. [CrossRef] [PubMed]
Snyder, B. D., Piazza, S., Edwards, W. T., and Hayes, W. C., 1993, “Role of Trabecular Morphology in the Etiology of Age-Related Vertebral Fractures,” Calcif. Tissue Int., 53(0), pp. S14–S22. [CrossRef] [PubMed]
Odgaard, A., and Linde, F., 1991, “The Underestimation of Young's Modulus in Compressive Testing of Cancellous Bone Specimens,” J. Biomech., 24(8), pp. 691–698. [CrossRef] [PubMed]
Linde, F., Hvid, I., and Madsen, F., 1992, “The Effect of Specimen Geometry on the Mechanical Behaviour of Trabecular Bone Specimens,” J. Biomech., 25(4), pp. 359–368. [CrossRef] [PubMed]
Lievers, W. B., Waldman, S. D., and Pilkey, A. K., 2010, “Minimizing Specimen Length in Elastic Testing of End-Constrained Cancellous Bone,” J. Mech. Behav. Biomed. Mater., 3(1), pp. 22–30. [CrossRef] [PubMed]
Ün, K., Bevill, G., and Keaveny, T. M., 2006, “The Effects of Side-Artifacts on the Elastic Modulus of Trabecular Bone,” J. Biomech., 39(11), pp. 1955–1963. [CrossRef] [PubMed]
Bevill, G., Easley, S. K., and Keaveny, T. M., 2007, “Side-Artifact Errors in Yield Strength and Elastic Modulus for Human Trabecular Bone and Their Dependence on Bone Volume Fraction and Anatomic Site,” J. Biomech., 40(15), pp. 3381–3388. [CrossRef] [PubMed]
Lievers, W. B., Petryshyn, A. C., Poljsak, A. S., Waldman, S. D., and Pilkey, A. K., 2010, “Specimen Diameter and “Side Artifacts” in Cancellous Bone Evaluated Using End-Constrained Elastic Tension,” Bone, 47(2), pp. 371–377. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

Sequence of steps for obtaining bone cores for compression testing. (a) Isolated human vertebra with posterior elements removed. (b) Removal of remaining adhered soft tissue from bone. (c) Sagittal slices from one half of vertebra. (d) Transverse slices from second half of vertebra. (e) Coring with a diamond tipped coring bit to obtain cylindrical bone cores. (f) Precision milling with a six-fluted end mill to obtain plano-parallel surfaces on the cores. (g) Sample bone core (center) with sapphire cylinders. Cores were color coded with india ink to separate cores from different slices.

Grahic Jump Location
Fig. 2

Sequence of steps for demarcating each vertebral core based on CT data. (a) CT scan image of cored vertebra. (b1, b2 and c1, c2) 3D bone model of each half of cored vertebra after segmentation in MIMICS 14.0. (d) 3D bone model of intact vertebra before coring (e) Merged bone models of two cored vertebral halves. (f) Coregistration of intact vertebral model and merged, cored vertebral model. (g) Boolean operation to subtract cored vertebral model from intact vertebral model to obtain segmented, 3D models of individual cores.

Grahic Jump Location
Fig. 3

(a) Schematic showing cross section of the Zetos loading unit. The load cell, PZA, and bone chamber are arranged in a series within the system. (b) A closer look at the bone chamber, where the bone specimen is placed between the two sapphire cylinders. Loading end of PZA is comprised of a convex surface to account for nonparallel surface and prevent moment transfer onto the specimen.

Grahic Jump Location
Fig. 4

Representative stress–strain curve from compression testing of vertebral trabecular bone cores. Shaded grey portion of curve (0% < ε < 0.2%) identifies the points used to calculate the apparent elastic modulus. Dashed line indicates use of the 0.2% offset technique for obtaining yield stress and strain.

Grahic Jump Location
Fig. 5

Measured apparent density as a ratio of ash weight over bulk, wet volume (ρapp) versus computed tomography derived density (ρCT)

Grahic Jump Location
Fig. 6

Variation of elastic modulus (E) with CT-derived apparent density (ρCT). SI = superior–inferior or axial direction. TR = transverse direction. The fitted linear regressions for the axial apparent elastic modulus with the physically measured apparent density (ρapp) has been included to illustrate the excellent agreement between correlations of (E) with (ρapp) and (ρCT).

Grahic Jump Location
Fig. 7

Variation of apparent yield strength (σy) with respect to CT-derived apparent density (ρCT). SI = superior–inferior or axial direction. TR = transverse direction.

Grahic Jump Location
Fig. 8

Variation of apparent yield strain (εy) with respect to to CT-derived apparent density (ρCT). SI = superior–inferior or axial direction. TR = transverse direction.

Grahic Jump Location
Fig. 9

Ratio of axial (SI)–to–transverse (TR) elastic moduli (ESI/ETR) and yield stress (σSI/σTR) for the range of density tested. Elastic moduli and yield stress computed from regression equations (3)–(6) in Table 2



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