Research Papers

Quantitative Analysis of the Nonlinear Displacement–Load Behavior of the Lumbar Spine

[+] Author and Article Information
Andrew D. Hanlon, Daniel J. Cook, Matthew S. Yeager

Department of Neurosurgery,
Allegheny General Hospital,
420 East North Avenue,
Pittsburgh, PA 15212

Boyle C. Cheng

Department of Neurosurgery,
Allegheny General Hospital,
420 East North Avenue,
Pittsburgh, PA 15212
e-mail: boylecheng@yahoo.com

1Corresponding author.

Manuscript received July 11, 2013; final manuscript received May 22, 2014; accepted manuscript posted May 29, 2014; published online June 11, 2014. Assoc. Editor: Brian D. Stemper.

J Biomech Eng 136(8), 081009 (Jun 11, 2014) (7 pages) Paper No: BIO-13-1313; doi: 10.1115/1.4027754 History: Received July 11, 2013; Revised May 22, 2014; Accepted May 29, 2014

There is currently no universal model or fitting method to characterize the visco-elastic behavior of the lumbar spine observed in displacement versus load hysteresis loops. In this study, proposed methods for fitting these loops, along with the metrics obtained, were thoroughly analyzed. A spline fitting technique was shown to provide a consistent approximation of spinal kinetic behavior that can be differentiated and integrated. Using this tool, previously established metrics were analyzed using data from two separate studies evaluating different motion preservation technologies. Many of the metrics, however, provided no significant differences beyond range of motion analysis. Particular attention was paid to how different definitions of the neutral zone capture the high-flexibility region often seen in lumbar hysteresis loops. As a result, the maximum slope was introduced and shown to be well defined. This new parameter offers promise as a descriptive measurement of spinal instability in vitro and may have future implications in clinical diagnosis and treatment of spinal instability. In particular, it could help in assigning treatments to specific stabilizing effects in the lumbar spine.

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


Adams, M. A., and Dolan, P., 2005, “Spine Biomechanics,” J. Biomech., 38(10), pp. 1972–1983. [CrossRef] [PubMed]
Goel, V. K., Panjabi, M. M., Patwardhan, A. G., Dooris, A. P., and Serhan, H., 2006, “Test Protocols for Evaluation of Spinal Implants,” J. Bone Joint Surg. Am., 88(Suppl 2), pp. 103–109. [CrossRef] [PubMed]
Panjabi, M. M., 1988, “Biomechanical Evaluation of Spinal Fixation Devices: I. A Conceptual Framework,” Spine13(10), pp. 1129–1134. [CrossRef] [PubMed]
Panjabi, M. M., Oxland, T. R., Yamamoto, I., and Crisco, J. J., 1994, “Mechanical Behavior of the Human Lumbar and Lumbosacral Spine as Shown by Three-Dimensional Load-Displacement Curves,” J. Bone Joint Surg. Am., 76(3), pp. 413–424. [PubMed]
Chutkan, N. B., Zhou, H., Akins, J. P., and Wenger, K. H., 2008, “Effects of Facetectomy and Crosslink Augmentation on Motion Segment Flexibility in Posterior Lumbar Interbody Fusion,” Spine33(22), pp. E828–835. [CrossRef] [PubMed]
Thompson, R. E., Barker, T. M., and Pearcy, M. J., 2003, “Defining the Neutral Zone of Sheep Intervertebral Joints During Dynamic Motions: An in vitro Study,” Clin. Biomech.18(2), pp. 89–98. [CrossRef]
Gay, R. E., Ilharreborde, B., Zhao, K., Zhao, C., and An, K. N., 2006, “Sagittal Plane Motion in the Human Lumbar Spine: Comparison of the in vitro Quasistatic Neutral Zone and Dynamic Motion Parameters,” Clin. Biomech.21(9), pp. 914–919. [CrossRef]
Hongo, M., Gay, R. E., Hsu, J. T., Zhao, K. D., Ilharreborde, B., Berglund, L. J., and An, K. N., 2008, “Effect of Multiple Freeze-Thaw Cycles on Intervertebral Dynamic Motion Characteristics in the Porcine Lumbar Spine,” J. Biomech., 41(4), pp. 916–920. [CrossRef] [PubMed]
Phillips, F. M., Tzermiadianos, M. N., Voronov, L. I., Havey, R. M., Carandang, G., Renner, S. M., Rosler, D. M., Ochoa, J. A., and Patwardhan, A. G., 2009, “Effect of the Total Facet Arthroplasty System After Complete Laminectomy-Facetectomy on the Biomechanics of Implanted and Adjacent Segments,” Spine J., 9(1), pp. 96–102. [CrossRef] [PubMed]
Panjabi, M. M., 1992, “The Stabilizing System of the Spine. Part II. Neutral Zone and Instability Hypothesis,” J. Spinal Disord., 5(4), pp. 390–396. [CrossRef] [PubMed]
Wilke, H. J., Wenger, K., and Claes, L., 1998, “Testing Criteria for Spinal Implants: Recommendations for the Standardization of in vitro Stability Testing of Spinal Implants,” Eur. Spine J., 7(2), pp. 148–154. [CrossRef] [PubMed]
Busscher, I., van der Veen, A. J., van Dieen, J. H., Kingma, I., Verkerke, G. J., and Veldhuizen, A. G., 2010, “In vitro Biomechanical Characteristics of the Spine: A Comparison Between Human and Porcine Spinal Segments,” Spine, 35(2), pp. E35–42. [CrossRef] [PubMed]
Deguchi, M., Cheng, B. C., Sato, K., Matsuyama, Y., and Zdeblick, T. A., 1998, “Biomechanical Evaluation of Translaminar Facet Joint Fixation. A Comparative Study of Poly-l-Lactide Pins, Screws, and Pedicle Fixation,” Spine, 23(12), pp. 1307–1312. [CrossRef] [PubMed]
Crawford, N. R., Peles, J. D., and Dickman, C. A., 1998, “The Spinal Lax Zone and Neutral Zone: Measurement Techniques and Parameter Comparisons,” J. Spinal Disord., 11(5), pp. 416–429. [CrossRef] [PubMed]
Pflugmacher, R., Schleicher, P., Gumnior, S., Turan, O., Scholz, M., Eindorf, T., Haas, N. P., and Kandziora, F., 2004, “Biomechanical Comparison of Bioabsorbable Cervical Spine Interbody Fusion Cages,” Spine29(16), pp. 1717–1722. [CrossRef] [PubMed]
Pflugmacher, R., Schleicher, P., Schaefer, J., Scholz, M., Ludwig, K., Khodadadyan-Klostermann, C., Haas, N. P., and Kandziora, F., 2004, “Biomechanical Comparison of Expandable Cages for Vertebral Body Replacement in the Thoracolumbar Spine,” Spine (Phila Pa 1976), 29(13), pp. 1413–1419. [CrossRef] [PubMed]
Mimura, M., Panjabi, M. M., Oxland, T. R., Crisco, J. J., Yamamoto, I., and Vasavada, A., 1994, “Disc Degeneration Affects the Multidirectional Flexibility of the Lumbar Spine,” Spine19(12), pp. 1371–1380. [CrossRef] [PubMed]
Cheng, L., Yu, Y., Zhu, R., Lv, H., Jia, Y., Zeng, Z., Chen, B., and Ding, Z., 2011, “Structural Stability of Different Reconstruction Techniques Following Total Sacrectomy: A Biomechanical Study,” Clin. Biomech.26(10), pp. 977–981. [CrossRef]
Lee, M. J., Bransford, R. J., Bellabarba, C., Chapman, J. R., Cohen, A. M., Harrington, R. M., and Ching, R. P., 2010, “The Effect of Bilateral Laminotomy Versus Laminectomy on the Motion and Stiffness of the Human Lumbar Spine: A Biomechanical Comparison,” Spine35(19), pp. 1789–1793. [CrossRef] [PubMed]
Smolders, L. A., Kingma, I., Bergknut, N., van der Veen, A. J., Dhert, W. J., Hazewinkel, H. A., van Dieen, J. H., and Meij, B. P., 2012, “Biomechanical Assessment of the Effects of Decompressive Surgery in Non-Chondrodystrophic and Chondrodystrophic Canine Multisegmented Lumbar Spines,” Eur. Spine J., 21(9), pp. 1692–1699. [CrossRef] [PubMed]
Panjabi, M. M., and Timm, J. P., 2007, “Development of Stabilimax NZ From Biomechanical Principles,” SAS J., 1(1), pp. 2–7. [CrossRef]
Cook, D. J., Yeager, M. S., and Cheng, B. C., 2012, “Interpedicular Travel in the Evaluation of Spinal Implants: An Application in Posterior Dynamic Stabilization,” Spine37(11), pp. 923–931. [CrossRef] [PubMed]
Lundgren, J., 2007, “Splinefit,” Matlab Central. Available at: http://www.mathworks.com/matlabcentral/fileexchange/13812-splinefit
Kranenburg, H. J., Meij, B. P., Onis, D., van der Veen, A. J., Saralidze, K., Smolders, L. A., Huizinga, J. G., Knetsch, M. L., Luijten, P. R., Visser, F., Voorhout, G., Dhert, W. J., Hazewinkel, H. A., and Koole, L. H., 2012, “Design, Synthesis, Imaging, and Biomechanics of a Softness-Gradient Hydrogel Nucleus Pulposus Prosthesis in a Canine Lumbar Spine Model,” J. Biomed. Mater. Res. B: Appl. Biomater., 100(8), pp. 2148–2155. [CrossRef] [PubMed]
Arthur, A., Cannella, M., Keane, M., Singhatat, W., Vresilovic, E., and Marcolongo, M., 2010, “Fill of the Nucleus Cavity Affects Mechanical Stability in Compression, Bending, and Torsion of a Spine Segment, Which has Undergone Nucleus Replacement,” Spine35(11), pp. 1128–1135. [CrossRef] [PubMed]
Cannella, M., Arthur, A., Allen, S., Keane, M., Joshi, A., Vresilovic, E., and Marcolongo, M., 2008, “The Role of the Nucleus Pulposus in Neutral Zone Human Lumbar Intervertebral Disc Mechanics,” J. Biomech., 41(10), pp. 2104–2111. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

Radiographic image of an FSU from the nucleus augmentation study in the augmented condition

Grahic Jump Location
Fig. 2

Posterior view of lumbar spine specimen with the Stabilimax device at L4–L5

Grahic Jump Location
Fig. 3

A representative hysteresis loop of a functional spinal unit in flexion–extension in the third cycle of flexibility testing with a spline fit. The fit approximates the flexible region, the stiff region, and the transition between them.

Grahic Jump Location
Fig. 4

A hysteresis loop with the NZ shaded. The NZ slopes were averaged over the shaded region (NZ = neutral zone).

Grahic Jump Location
Fig. 5

A hysteresis loop with the LZ shaded. The LZ slopes were averaged over the shaded region. The zero load extrapolation lines were used to determine the bounds of the LZ as described by Crawford et al. [14] (LZ = lax zone).

Grahic Jump Location
Fig. 6

Average range of motion values with standard errors over the three conditions for the specimens which underwent a nucleus augmentation treatment (FE = flexion–extension and LB = lateral bending)

Grahic Jump Location
Fig. 7

Average maximum slope values with standard errors over the three conditions for the specimens which underwent a nucleus augmentation treatment (F = flexion, E = extension, LLB = left lateral bending, and RLB = right lateral bending)

Grahic Jump Location
Fig. 8

Average range of motion values with standard errors over the five conditions for the index level specimens in the Stabilimax study (FE = flexion–extension, LB = lateral bending, AT = axial torsion, and AC = axial compression)

Grahic Jump Location
Fig. 9

Average maximum slope values with standard errors over the five conditions for the specimens which underwent a Stabilimax treatment (F = flexion, E = extension, LLB = left lateral bending, RLB = right lateral bending, RT = right torsion, and LT = left torsion)

Grahic Jump Location
Fig. 10

A representative derivative of the hysteresis loop for a specimen in the third cycle of flexion. A clear area of high flexibility is seen between−1 and 0.75 Nm. The boundaries of Panjabi's neutral zone and Crawford's lax zone fail to capture that region.



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