Biomechanical influence of disk properties on the load transfer of healthy and degenerated disks using a poroelastic finite element model.

Spine degeneration is a pathology that will affect 80% of the population. Since the intervertebral disks play an important role in transmitting loads through the spine, the aim of this study was to evaluate the biomechanical impact of disk properties on the load carried by healthy (Thompson grade I) and degenerated (Thompson grades III and IV) disks. A three-dimensional parametric poroelastic finite element model of the L4/L5 motion segment was developed. Grade I, grade II, and grade IV disks were modeled by altering the biomechanical properties of both the annulus and nucleus. Models were validated using published creep experiments, in which a constant compressive axial stress of 0.35 MPa was applied for 4 h. Pore pressure (PP) and effective stress (S(E)) were analyzed as a function of time following loading application (1 min, 5 min, 45 min, 125 min, and 245 min) and discal region along the midsagittal profile for each disk grade. A design of experiments was further implemented to analyze the influence of six disk parameters (disk height (H), fiber proportion (%F), drained Young's modulus (E(a),E(n)), and initial permeability (k(a),k(n)) of both the annulus and nucleus) on load-sharing for disk grades I and IV. Simulations of grade I, grade III, and grade IV disks agreed well with the available published experimental data. Disk height (H) had a significant influence (p<0.05) on the PP and S(E) during the entire loading history for both healthy and degenerated disk models. Young's modulus of the annulus (E(a)) significantly affected not only S(E) in the annular region for both disk grades in the initial creep response but also S(E) in the nucleus zone for degenerated disks with further creep response. The nucleus and annulus permeabilities had a significant influence on the PP distribution for both disk grades, but this effect occurred at earlier stages of loading for degenerated than for healthy disk models. This is the first study that investigates the biomechanical influence of both geometrical and material disk properties on the load transfer of healthy and degenerated disks. Disk height is a significant parameter for both healthy and degenerated disks during the entire loading. Changes in the annulus stiffness, as well as in the annulus and nucleus permeability, control load-sharing in different ways for healthy and degenerated disks.

[1]  A. Nachemson.  Disc Pressure Measurements , 1981, Spine.

[2]  Lutz Claes,et al.  Validity and interobserver agreement of a new radiographic grading system for intervertebral disc degeneration: Part I. Lumbar spine , 2006, European Spine Journal.

[3]  G. Drouin,et al.  Rheological properties of the human lumbar spine ligaments. , 1991, Journal of biomedical engineering.

[4]  B R Simon,et al.  1985 Volvo Award in Biomechanics: Poroelastic Dynamic Structural Models of Rhesus Spinal Motion Segments , 1985, Spine.

[5]  Josep A Planell,et al.  Statistical factorial analysis on the poroelastic material properties sensitivity of the lumbar intervertebral disc under compression, flexion and axial rotation. , 2009, Journal of biomechanics.

[6]  V C Mow,et al.  Degeneration affects the anisotropic and nonlinear behaviors of human anulus fibrosus in compression. , 1998, Journal of biomechanics.

[7]  Lachlan J. Smith,et al.  The elastic fibre network of the human lumbar anulus fibrosus: architecture, mechanical function and potential role in the progression of intervertebral disc degeneration , 2009, European Spine Journal.

[8]  G. Andersson,et al.  Relationship between disc injury and manual lifting: A poroelastic finite element model study , 2008, Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine.

[9]  G A Dumas,et al.  Influence of material properties on the mechanical behaviour of the L5-S1 intervertebral disc in compression: a nonlinear finite element study. , 1991, Journal of biomedical engineering.

[10]  W C Hutton,et al.  The effect of posture on the role of the apophysial joints in resisting intervertebral compressive forces. , 1980, The Journal of bone and joint surgery. British volume.

[11]  Leonid Kalichman,et al.  Lumbar facet joint osteoarthritis: a review. , 2007, Seminars in arthritis and rheumatism.

[12]  Antonius Rohlmann,et al.  Analysis of the influence of disc degeneration on the mechanical behaviour of a lumbar motion segment using the finite element method. , 2006, Journal of biomechanics.

[13]  A Shirazi-Adl,et al.  Poroelastic creep response analysis of a lumbar motion segment in compression. , 1996, Journal of biomechanics.

[14]  Keita Ito,et al.  Fluid flow and convective transport of solutes within the intervertebral disc. , 2004, Journal of biomechanics.

[15]  J. P. Thompson,et al.  Preliminary Evaluation of a Scheme for Grading the Gross Morphology of the Human Intervertebral Disc , 1990, Spine.

[16]  Louis Y. Cheng,et al.  Dependence of mechanical behavior of the murine tail disc on regional material properties: a parametric finite element study. , 2005, Journal of biomechanical engineering.

[17]  P. Prithvi Raj,et al.  Intervertebral Disc: Anatomy‐Physiology‐Pathophysiology‐Treatment , 2008, Pain practice : the official journal of World Institute of Pain.

[18]  M. Biot.  General Theory of Three‐Dimensional Consolidation , 1941 .

[19]  G. Andersson,et al.  Intervertebral disc degeneration. Summary of an AAOS/NIH/ORS workshop, September 2005. , 2006, The Journal of bone and joint surgery. American volume.

[20]  J. Lotz,et al.  Compression-induced degeneration of the intervertebral disc: an in vivo mouse model and finite-element study. , 1998, Spine.

[21]  Pierre-Jean Arnoux,et al.  Finite element investigation of the loading rate effect on the spinal load-sharing changes under impact conditions. , 2009, Journal of biomechanics.

[22]  Lutz Claes,et al.  The risk of disc prolapses with complex loading in different degrees of disc degeneration - a finite element analysis. , 2007, Clinical biomechanics.

[23]  Ian H Parkinson,et al.  Antero-postero differences in cortical thickness and cortical porosity of T12 to L5 vertebral bodies. , 2006, Joint, bone, spine : revue du rhumatisme.

[24]  P J Reckwerdt,et al.  A proposal for a standard electronic anthropomorphic phantom for radiotherapy. , 1999, Medical physics.

[25]  P Kurowski,et al.  The Relationship of Degeneration of the Intervertebral Disc to Mechanical Loading Conditions on Lumbar Vertebrae , 1986, Spine.

[26]  L. Denaro,et al.  Intervertebral disc degeneration - an autopsy study , 2000, European Journal of Orthopaedic Surgery & Traumatology.

[27]  A. Shirazi-Adl,et al.  Finite–Element Simulation of Changes in the Fluid Content of Human Lumbar Discs: Mechanical and Clinical Implications , 1992, Spine.

[28]  B. Vernon‐roberts,et al.  Degenerative changes in the intervertebral discs of the lumbar spine and their sequelae. , 1977, Rheumatology and rehabilitation.

[29]  G B Andersson,et al.  The influence of lumbar disc height and cross-sectional area on the mechanical response of the disc to physiologic loading. , 1999, Spine.

[30]  P Brinckmann,et al.  Measurement of the Distribution of Axial Stress on the End-Plate of the Vertebral Body , 1981, Spine.

[31]  Hutton Wc,et al.  Do bending, twisting, and diurnal fluid changes in the disc affect the propensity to prolapse? A viscoelastic finite element model , 1996 .

[32]  A E Goodship,et al.  In Vivo Stress Measurement Can Predict Pain on Discography , 1996, Spine.

[33]  D. Elliott,et al.  Effects of Degeneration on the Biphasic Material Properties of Human Nucleus Pulposus in Confined Compression , 2005, Spine.

[34]  Howard S. An,et al.  The relationship between facet joint osteoarthritis and disc degeneration of the lumbar spine: an MRI study , 1999, European Spine Journal.

[35]  V K Goel,et al.  Effect of Disc Degeneration at One Level on the Adjacent Level in Axial Mode , 1991, Spine.

[36]  E C Teo,et al.  Poroelastic Analysis of Lumbar Spinal Stability in Combined Compression and Anterior Shear , 2004, Journal of spinal disorders & techniques.

[37]  M. Adams,et al.  Internal Intervertebral Disc Mechanics as Revealed by Stress Profilometry , 1992, Spine.

[38]  M. Adams,et al.  What is Intervertebral Disc Degeneration, and What Causes It? , 2006, Spine.

[39]  S. Aharinejad,et al.  Morphometric analysis of vertebrae and intervertebral discs as a basis of disc replacement. , 1990, The American journal of anatomy.

[40]  J P Laible,et al.  A Poroelastic-Swelling Finite Element Model With Application to the Intervertebral Disc , 1993, Spine.

[41]  Gunnar B J Andersson,et al.  Recent Advances in Analytical Modeling of Lumbar Disc Degeneration , 2004, Spine.

[42]  Hutton Wc,et al.  Can variations in intervertebral disc height affect the mechanical function of the disc , 1996 .

[43]  J. Paige Little,et al.  Parametric equations to represent the profile of the human intervertebral disc in the transverse plane , 2007, Medical & Biological Engineering & Computing.

[44]  V C Mow,et al.  The anisotropic hydraulic permeability of human lumbar anulus fibrosus. Influence of age, degeneration, direction, and water content. , 1999, Spine.

[45]  N Yoganandan,et al.  Finite element model of the human lower cervical spine: parametric analysis of the C4-C6 unit. , 1997, Journal of biomechanical engineering.