Parametric convergence sensitivity and validation of a finite element model of the human lumbar spine

The primary objective of this study was to generate a finite element model of the human lumbar spine (L1–L5), verify mesh convergence for each tissue constituent and perform an extensive validation using both kinematic/kinetic and stress/strain data. Mesh refinement was accomplished via convergence of strain energy density (SED) predictions for each spinal tissue. The converged model was validated based on range of motion, intradiscal pressure, facet force transmission, anterolateral cortical bone strain and anterior longitudinal ligament deformation predictions. Changes in mesh resolution had the biggest impact on SED predictions under axial rotation loading. Nonlinearity of the moment-rotation curves was accurately simulated and the model predictions on the aforementioned parameters were in good agreement with experimental data. The validated and converged model will be utilised to study the effects of degeneration on the lumbar spine biomechanics, as well as to investigate the mechanical underpinning of the contemporary treatment strategies.

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

[2]  Benjamin J. Ellis,et al.  Verification, validation and sensitivity studies in computational biomechanics , 2007, Computer methods in biomechanics and biomedical engineering.

[3]  Zvi Hashin,et al.  Continuum Theory of the Mechanics of Fibre-Reinforced Composites , 1984 .

[4]  A. M. Ahmed,et al.  Stress analysis of the lumbar disc-body unit in compression. A three-dimensional nonlinear finite element study. , 1984, Spine.

[5]  C. Whyne,et al.  Parametric finite element analysis of vertebral bodies affected by tumors. , 2001, Journal of biomechanics.

[6]  R. T. Hart,et al.  Contact analysis of a posterior cervical spine plate using a three-dimensional canine finite element model. , 1999, Journal of biomechanical engineering.

[7]  Neil R Crawford,et al.  The use of surface strain data and a neural networks solution method to determine lumbar facet joint loads during in vitro spine testing. , 2008, Journal of biomechanics.

[8]  Thomas R. Oxland,et al.  The Effect of Nucleotomy on Lumbar Spine Mechanics in Compression and Shear Loading , 2001, Spine.

[9]  L W Marks,et al.  The use of strain energy as a convergence criterion in the finite element modelling of bone and the effect of model geometry on stress convergence. , 1993, Journal of biomedical engineering.

[10]  L. Setton,et al.  Anisotropic and inhomogeneous tensile behavior of the human anulus fibrosus: experimental measurement and material model predictions. , 2001, Journal of biomechanical engineering.

[11]  Michael A. Adams,et al.  'Stress' distributions inside intervertebral discs , 1996 .

[12]  Narayan Yoganandan,et al.  Moment-rotation responses of the human lumbosacral spinal column. , 2007, Journal of biomechanics.

[13]  Thomas R Oxland,et al.  Accuracy and repeatability of a new method for measuring facet loads in the lumbar spine. , 2006, Journal of biomechanics.

[14]  P. Brinckmann,et al.  Interlaminar Shear Stresses and Laminae Separation in a Disc: Finite Element Analysis of the L3‐L4 Motion Segment Subjected to Axial Compressive Loads , 1995, Spine.

[15]  Tony M Keaveny,et al.  Quantitative computed tomography-based finite element models of the human lumbar vertebral body: effect of element size on stiffness, damage, and fracture strength predictions. , 2003, Journal of biomechanical engineering.

[16]  Gerhard A. Holzapfel,et al.  An Anisotropic Model for Annulus Tissue and Enhanced Finite Element Analyses of Intact Lumbar Disc Bodies , 2001 .

[17]  Thomas R Oxland,et al.  The Effect of Dynamic Posterior Stabilization on Facet Joint Contact Forces: An In Vitro Investigation , 2008, Spine.

[18]  G. Bergmann,et al.  Estimation of muscle forces in the lumbar spine during upper-body inclination. , 2001, Clinical biomechanics.

[19]  M. Fröhlich,et al.  Multi-segment FEA of the human lumbar spine including the heterogeneity of the annulus fibrosus , 2004 .

[20]  L. Claes,et al.  New in vivo measurements of pressures in the intervertebral disc in daily life. , 1999, Spine.

[21]  Lutz Claes,et al.  Application of a new calibration method for a three-dimensional finite element model of a human lumbar annulus fibrosus. , 2006, Clinical biomechanics.

[22]  M. Panjabi,et al.  How Does Posture Affect Coupling in the Lumbar Spine? , 1989, Spine.

[23]  Gregory P Grieve Fcsp DipTP Mmacp.  Clinical Anatomy of the Lumbar Spine and Sacrum , 1997 .

[24]  Gunnar B J Andersson,et al.  Modeling changes in intervertebral disc mechanics with degeneration. , 2006, The Journal of bone and joint surgery. American volume.

[25]  Alan Morris,et al.  A Practical Guide to Reliable Finite Element Modelling: Morris/A Practical Guide to Reliable Finite Element Modelling , 2008 .

[26]  S. Klisch,et al.  Application of a fiber-reinforced continuum theory to multiple deformations of the annulus fibrosus. , 1999, Journal of biomechanics.

[27]  Christian M Puttlitz,et al.  The micromechanical role of the annulus fibrosus components under physiological loading of the lumbar spine. , 2010, Journal of biomechanical engineering.

[28]  Lutz Claes,et al.  Stepwise reduction of functional spinal structures increase vertebral translation and intradiscal pressure. , 2007, Journal of biomechanics.

[29]  J. Lotz,et al.  Radial tensile properties of the lumbar annulus fibrosus are site and degeneration dependent , 1997, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[30]  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.

[31]  Alison C Jones,et al.  Finite element analysis of the spine: towards a framework of verification, validation and sensitivity analysis. , 2008, Medical engineering & physics.

[32]  Y K Liu,et al.  A three-dimensional nonlinear finite element model of lumbar intervertebral joint in torsion. , 1987, Journal of biomechanical engineering.

[33]  Jeffrey C Lotz,et al.  Theoretical model and experimental results for the nonlinear elastic behavior of human annulus fibrosus , 2004, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[34]  V. Goel,et al.  Load-Sharing Between Anterior and Posterior Elements in a Lumbar Motion Segment Implanted With an Artificial Disc , 2001, Spine.

[35]  Josep A Planell,et al.  Finite Element Study of a Novel Intervertebral Disc Substitute , 2005, Spine.

[36]  M M Panjabi,et al.  Disc Degeneration Affects the Multidirectional Flexibility of the Lumbar Spine , 1994, Spine.

[37]  A Rohlmann,et al.  Realistic loading conditions for upper body bending. , 2009, Journal of biomechanics.

[38]  Gunnar B J Andersson,et al.  Inclusion of regional poroelastic material properties better predicts biomechanical behavior of lumbar discs subjected to dynamic loading. , 2007, Journal of biomechanics.

[39]  M M Panjabi,et al.  Mechanical behavior of the human lumbar and lumbosacral spine as shown by three-dimensional load-displacement curves. , 1994, The Journal of bone and joint surgery. American volume.

[40]  Lutz Claes,et al.  Stepwise reduction of functional spinal structures increase range of motion and change lordosis angle. , 2007, Journal of biomechanics.

[41]  T. Keaveny,et al.  Finite Element Modeling of the Human Thoracolumbar Spine , 2003, Spine.

[42]  Manohar M. Panjabi,et al.  Clinical Biomechanics of the Spine , 1978 .

[43]  C. Puttlitz,et al.  Lower cervical spine facet cartilage thickness mapping. , 2008, Osteoarthritis and cartilage.

[44]  William C. Hutton,et al.  Can Variations in Intervertebral Disc Height Affect the Mechanical Function of the Disc? , 1996, Spine.

[45]  H. Serhan,et al.  Biomechanics of the posterior lumbar articulating elements. , 2007, Neurosurgical focus.