0
Research Papers

Experimental Validation of Finite Element Analysis of Human Vertebral Collapse Under Large Compressive Strains

[+] Author and Article Information
Hadi S. Hosseini

Institute for Surgical Technology
and Biomechanics,
University of Bern,
Stauffacherstr. 78,
Bern CH-3014, Switzerland
e-mail: hadi.seyed@istb.unibe.ch

Allison L. Clouthier, Philippe K. Zysset

Institute for Surgical Technology
and Biomechanics,
University of Bern,
Stauffacherstr. 78,
Bern CH-3014, Switzerland

Contributed by the Bioengineering Division of ASME for publication in the JOURNAL OF BIOMECHANICAL ENGINEERING. Manuscript received July 9, 2013; final manuscript received December 29, 2013; accepted manuscript posted January 2, 2014; published online March 24, 2014. Assoc. Editor: Pasquale Vena.

J Biomech Eng 136(4), 041006 (Mar 24, 2014) (10 pages) Paper No: BIO-13-1307; doi: 10.1115/1.4026409 History: Received July 09, 2013; Revised December 29, 2013; Accepted January 02, 2014

Osteoporosis-related vertebral fractures represent a major health problem in elderly populations. Such fractures can often only be diagnosed after a substantial deformation history of the vertebral body. Therefore, it remains a challenge for clinicians to distinguish between stable and progressive potentially harmful fractures. Accordingly, novel criteria for selection of the appropriate conservative or surgical treatment are urgently needed. Computer tomography-based finite element analysis is an increasingly accepted method to predict the quasi-static vertebral strength and to follow up this small strain property longitudinally in time. A recent development in constitutive modeling allows us to simulate strain localization and densification in trabecular bone under large compressive strains without mesh dependence. The aim of this work was to validate this recently developed constitutive model of trabecular bone for the prediction of strain localization and densification in the human vertebral body subjected to large compressive deformation. A custom-made stepwise loading device mounted in a high resolution peripheral computer tomography system was used to describe the progressive collapse of 13 human vertebrae under axial compression. Continuum finite element analyses of the 13 compression tests were realized and the zones of high volumetric strain were compared with the experiments. A fair qualitative correspondence of the strain localization zone between the experiment and finite element analysis was achieved in 9 out of 13 tests and significant correlations of the volumetric strains were obtained throughout the range of applied axial compression. Interestingly, the stepwise propagating localization zones in trabecular bone converged to the buckling locations in the cortical shell. While the adopted continuum finite element approach still suffers from several limitations, these encouraging preliminary results towardsthe prediction of extended vertebral collapse may help in assessing fracture stability in future work.

FIGURES IN THIS ARTICLE
<>
Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.

References

Magerl, F., Aebi, M., Gertzbein, S., Harms, J., and Nazarian, S., 1994, “A Comprehensive Classification of Thoracic and Lumbar Injuries,” Eur. Spine J., 3, pp. 184–201. [CrossRef]
Johnell, O. and Kanis, J., 2005, “Epidemiology of Osteoporotic Fractures,” Osteoporosis Int., 16, pp. S3–S7. [CrossRef]
Johnell, O. and Kanis, J., 2006, “An Estimate of the Worldwide Prevalence and Disability Associated With Osteoporotic Fractures”. Osteoporosis Int., 17, pp. 1726–1733. [CrossRef]
Kanis, J., McCloskey, E., Johansson, H., Oden, A., Melton, L., and Khaltaev, N., 2008, “A Reference Standard for the Description of Osteoporosis,” Bone, 42, pp. 467–475. [CrossRef]
Old, J. and Calvert, M., 2004, “Vertebral Compression Fractures in the Elderly,” Am. Fam. Physician, 69, pp. 111–116.
Lydick, E., Martin, A., and Yawn, B., 1996. “Impact of Fears on Quality of Life in Patients With a Silent Disease: Osteoporosis,” Clin. Ther., 18(6), pp. 1307–1315. [CrossRef]
NIH, 2000, “NIH Consensus Statement: Osteoporosis Prevention, Diagnosis, and Therapy,” Technical Report, National Institutes of Health, Bethesda, MD.
Nazarian, A. and Müller, R., 2004, “Time-Lapsed Microstructural Imaging of Bone Failure Behavior,” J. Biomech., 37(1), pp. 55–65. [CrossRef]
Thurner, P., Wyss, P., Voide, R., Stauber, M., Stampanoni, M., Sennhauser, U., and Müller, R., 2006, “Time-Lapsed Investigation of Three-Dimensional Failure and Damage Accumulation in Trabecular Bone Using Synchrotron Light,” Bone, 39(2), pp. 289–299. [CrossRef]
Hulme, P., Ferguson, S., and Boyd, S., 2008, “Determination of Vertebral Endplate Deformation Under Load Using Micro-Computed Tomography,” J. Biomech., 41(1), pp. 78–85. [CrossRef]
Bay, B., Smith, T., Fyhrie, D., and Saad, M., 1999, “Digital Volume Correlation: Three-Dimensional Strain Mapping Using X-Ray Tomography,” Exp. Mech., 39, pp. 217–226. [CrossRef]
Verhulp, E., Rietbergen, B., and Huiskes, R., 2004, “A Three-Dimensional Digital Image Correlation Technique for Strain Measurements in Microstructures,” J. Biomech., 37(9), pp. 1313–1320. [CrossRef]
Silva, M. and Gibson, L., 1997, “Modeling the Mechanical Behavior of Vertebral Trabecular Bone: Effects of Age-Related Changes in Microstructure,” Bone, 21(2), pp. 191–199. [CrossRef]
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]
Bevill, G. and Keaveny, T. M., 2009, “Trabecular Bone Strength Predictions Using Finite Element Analysis of Micro-Scale Images at Limited Spatial Resolution,” Bone, 44(4), pp. 579–584. [CrossRef]
Charlebois, M., Jirasek, M., and Zysset, P., 2010, “A Nonlocal Constitutive Model for Trabecular Bone Softening in Compression,” Biomech. Model. Mechanobiol., 9, pp. 597–611. [CrossRef]
Dall'Ara, E., Schmidt, R., and Zysset, P., 2012, “Microindentation can Discriminate Between Damaged and Intact Human Bone Tissue,” Bone, 50(4), pp. 925–929. [CrossRef]
Zysset, P., 1994, “A Constitutive Law for Trabecular Bone,” Ph.D. thesis, Swiss Federal Institute of Technology (EPFL), Lausanne, Switzerland.
DallAra, E., Pahr, D., Varga, P., Kainberger, F., and Zysset, P., 2012, “QCT-Based Finite Element Models Predict Human Vertebral Strength in vitro Significantly Better Than Simulated DEXA,” Osteoporosis Int., 23, pp. 563–572. [CrossRef]
Hosseini, H. S., Pahr, D. H., and Zysset, P. K., 2012, “Modeling and Experimental Validation of Trabecular Bone Damage, Softening and Densification Under Large Compressive Strains,” J. Mech. Behav. Biomed. Mater., 15, pp. 93–102. [CrossRef]
Pahr, D. and Zysset, P., 2008, “Influence of Boundary Conditions on Computed Apparent Elastic Properties of Cancellous Bone,” Biomech. Model. Mechanobio., 7, pp. 463–476. [CrossRef]
Homminga, J. J., 2003, Towards a Rational Definition of Osteoporosis, Technische Universiteit, Eindhoven, Netherlands.
Adams, M. A., Pollintine, P., Tobias, J. H., Wakley, G. K., and Dolan, P., 2006, “Intervertebral Disc Degeneration can Predispose to Anterior Vertebral Fractures in the Thoracolumbar Spine,” J. Bone Miner. Res., 21(9), pp. 1409–1416. [CrossRef]
Chevalier, Y., Pahr, D., and Zysset, P. K., 2009, “The Role of Cortical Shell and Trabecular Fabric in Finite Element Analysis of the Human Vertebral Body,” ASME J. Biomech. Eng., 131(11), p. 111003. [CrossRef]
Hosseini, H. S., 2013, “Constitutive Modeling of Trabecular Bone in Large Strain Compression,” Ph.D. thesis, University of Bern, Bern.
Schwiedrzik, J., Wolfram, U., and Zysset, P., 2013, “A Generalized Anisotropic Quadric Yield Criterion and Its Application to Bone Tissue at Multiple Length Scales,” Biomech. Model. Mechanobiol. (submitted).
Poh, L. and Swaddiwudhipong, S., 2009, “Gradient-Enhanced Softening Material Models,” Int. J. Plast., 25(11), pp. 2094–2121. [CrossRef]
Poh, L. H. and Swaddiwudhipong, S., 2009, “Over-Nonlocal Gradient Enhanced Plastic-Damage Model for Concrete,” Int. J. Solids Struct., 46(2526), pp. 4369–4378. [CrossRef]
Pahr, D. H. and Zysset, P. K., 2009, “From High-Resolution CT Data to Finite Element Models: Development of an Integrated Modular Framework,” Comput. Methods Biomech. Biomed. Eng., 12(1), pp. 45–57. [CrossRef]
Chevalier, Y., Charlebois, M., Pahr, 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]
Pahr, D. H. and Zysset, P., 2009, “A Comparison of Enhanced Continuum FE With Micro FE Models of Human Vertebral Bodies,” J. Biomech., 42(4), pp. 455–462. [CrossRef]
Zysset, P. K., Goulet, R. W., and Hollister, S. J., 1998, “A Global Relationship Between Trabecular Bone Morphology and Homogenized Elastic Properties,” ASME J. Biomech. Eng., 120(5), pp. 640–646. [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]
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]
Zysset, P. K. and Rincon, L., 2006, “An Alternative Fabric-Based Yield and Failure Criterion for Trabecular Bone,” Mech. Biol. Tissue, pp. 457–470. [CrossRef]
Rincon-Kohli, L., 2003, “Identification of a Multiaxial Failure Criterion for Human Trabecular Bone,” Ph.D. thesis, Swiss Federal Institute of Technology (EPFL), Lausanne, Switzerland.
Belytschko, T., Liu, W. K., and Moran, B., 2000, Nonlinear Finite Elements for Continua and Structures. John Wiley and Sons, New York.
Marsden, J. E. and Hughes, T. J. R., 1983, Mathematical Foundations of Elasticity, Prentice- Hall, Englewood Cliffs, NJ.
Lin, L. I.-K., 1989, “A Concordance Correlation Coefficient to Evaluate Reproducibility,” Biometrics, 45(1), pp. 255–268. [CrossRef]
Eswaran, S. K., Gupta, A., and Keaveny, T. M., 2007, “Locations of Bone Tissue at High Risk of Initial Failure During Compressive Loading of the Human Vertebral Body,” Bone, 41(4), pp. 733–739. [CrossRef]

Figures

Grahic Jump Location
Fig. 2

Compression device used to apply the axial load to the vertebral body. The device is mounted on the XtremeCT machine after each loading step.

Grahic Jump Location
Fig. 1

(a) A vertebral body cleaned from the soft tissues and the posterior elements, and (b) the PMMA end-cap, upon which the shape of the endplate is printed

Grahic Jump Location
Fig. 3

(a) The μCT reconstruction of a vertebral body, and (b) the hFE mesh showing the distribution of the homogenized BV/TV for each element

Grahic Jump Location
Fig. 4

The BV/TV of the transverse slices calculated from the (a) μCT image, and (b) FE mesh in their intact state

Grahic Jump Location
Fig. 5

Three-dimensional representation of a vertebral body collapse up to very large compression. Trabecular bone densification is represented by the volumetric strain ln J. The corresponding applied compressive deformation is distinguished by the red circle on the force-displacement curves.

Grahic Jump Location
Fig. 6

Qualitative analysis of the vertebral body collapse in five specimens. Trabecular bone densification is represented by the volumetric strain ln J. The corresponding applied compressive deformation is distinguished by the red circle on the force-displacement curves.

Grahic Jump Location
Fig. 7

Qualitative analysis of the vertebral body collapse in four specimens. Trabecular bone densification is represented by the volumetric strain ln J. The corresponding applied compressive deformation is distinguished by the red circle on the force-displacement curves.

Grahic Jump Location
Fig. 8

Qualitative analysis of the vertebral body collapse in four specimens where the hFE simulations failed in the correct prediction of the localization pattern. Trabecular bone densification is represented by the volumetric strain ln J. The corresponding applied compressive deformation is distinguished by the red circle on the force-displacement curves.

Grahic Jump Location
Fig. 9

Quantitative analysis of Δ(BV/TV) for 13 vertebrae up to 68% of apparent compressive strain

Tables

Errata

Discussions

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