Comparison of eight published static finite element models of the intact lumbar spine: Predictive power of models improves when combined together
Introduction
Accurate and clinically relevant modeling of complex biological systems such as the human lumbar spine remains challenging, yet promising, with the potential to substantially enhance the quality of patient care. Due to its ability to represent intricate systems with material nonlinearities, irregular loading, and geometrical and material domains, the finite element (FE) method has been recognized as an important computational tool in various biomedical fields (Zhang and Teo, 2008) and has been widely adopted for describing spinal biomechanics (Schmidt et al., 2013). In comparison to in vitro or in vivo approaches, computational methods are advantageous in offering cost efficient and powerful response solutions while at the same time effectively dealing with the ethical concerns related to the use of live animals in experiments. Moreover, use of computational models may greatly diminish the need for experimental investigations that utilize post mortem human and animal specimens. For example, finite element models provide improved insight into the functional mechanisms of the spine by assessing the isolated effect of various parameters independently – a feature that has been invaluable with respect to the design/optimization of spinal implants (Fagan et al., 2002a, Schmidt et al., 2013, Zhang and Teo, 2008).
Despite the proven success of computational studies in other disciplines, the FE method׳s role in clinical spine research has sometimes been questioned (Viceconti et al., 2005). The uncertainty and high variability of tissue material properties, the anatomical complexity of spinal structures (Panjabi et al., 1992, Panjabi et al., 1993), and the unknown loading (Rohlmann et al., 2009, Wilke et al., 1998) and boundary conditions, particularly in vivo, has cast doubt on the accuracy and reliability of FE model predictions. The inherent geometric and material property differences among individuals and alterations in these parameters due to age, sex and degeneration may limit the widespread applicability of the reported results. To gain confidence in and to enhance the predictive quality of FE models, recommendations have been made on how to develop suitable models in order to address research questions within an adequate degree of predictive accuracy (Anderson et al., 2007, Jones and Wilcox, 2008, Oreskes et al., 1994, Roache, 1998, Viceconti, 2011, Viceconti et al., 2005). These standards comprise three main steps: code verification, sensitivity analyses of uncertain model input parameters, and task-specific validations of the model.
The verification of the code poses the least concern as the vast majority of computational studies nowadays employ extensively verified, commercially available FE software. The analysis of the sensitivity to alterations in geometrical (Dupont et al., 2002, Meijer et al., 2011, Natarajan and Andersson, 1999, Niemeyer et al., 2012, Noailly et al., 2007, Robin et al., 1994), material (Fagan et al., 2002b, Lee and Teo, 2005, Rao and Dumas, 1991, Shirazi-Adl, 1994a, Zander et al., 2004) or loading parameters (Dreischarf et al., 2011, Dreischarf et al., 2012, Rohlmann et al., 2009); however, demands more time and effort and has hence only occasionally been carried out. It has been shown that the range of motion (RoM) of a lumbar motion segment is strongly affected by the disc height (Meijer et al., 2011, Natarajan and Andersson, 1999, Niemeyer et al., 2012, Robin et al., 1994) and material properties (e.g. ligament properties (Zander et al., 2004). Furthermore, appropriate loading conditions (Dreischarf et al., 2011, Dreischarf et al., 2012) are necessary to realistically simulate relevant tasks under maximal voluntary motion measured in vivo (Pearcy, 1985, Pearcy et al., 1984, Pearcy and Tibrewal, 1984, Wilke et al., 2001).
The term ‘validation’ merits attention as it remains controversial. Validation is commonly used to indicate that model predictions are consistent with observations. However, it is intractable to completely validate numerical models because it is not possible to account for the multiplicity of their inherent degrees of freedom in an experiment (Oreskes et al., 1994). It is, however, generally accepted that greater number and diversity of corroborating observations between a model and experimental data increases the probability that the model predictions are not flawed (Oreskes et al., 1994, Viceconti et al., 2005). To increase the confidence in a model, the number of free independent parameters employed to construct the model should remain low to decrease the risk of non-uniqueness. Detailed experimental data on the lumbar spine that would allow for a thorough validation of model predictions remain, however, limited. For example, measurements are often only performed at a single level. Model validation is therefore often performed by comparing the calculated results with the limited data that is available from in vitro studies (Moramarco et al., 2010, Zander et al., 2009). However, experimental setups, specimens, loading and boundary conditions substantially differ among various studies (Brinckmann and Grootenboer, 1991, Kettler et al., 2011, Rohlmann et al., 2001b, Wilke et al., 1994), and these differences are often neglected with regard to the resulting data. Furthermore, the validation of numerical models should preferably include as many relevant outputs as possible (Woldtvedt et al., 2011), as some may be more sensitive to model assumptions than others under specific loading conditions. Moreover, for clinically relevant parameters such as the facet joint forces (FJF), which have considerable dependence on loading and geometry, almost no in vivo data exist (Wilson et al., 2006).
Well-established FE models should incorporate the aforementioned three steps to meet the conditions for a meaningful numerical study. Despite these requirements, most FE studies account for only one spinal geometry with one set of material properties and are validated with very few available experimental data. This raises questions with regard to the reliability/comparability of their predictions under various conditions, on the range of results of these numerical predictions, and on their agreement with in vitro values. Concerns also exist when attempting to validate predictions with in vivo data under complex combined loading modes (e.g. compression and bending). To address these issues, one may compare the salient predictions of peer-reviewed models obtained under nearly identical loading and boundary conditions. For this purpose and due to the importance and complexity of the lumbar spine, this novel multicenter study was undertaken to compare the results of eight well-established FE models of the lumbar spine that have been developed, validated and applied for many years in different research centers around the globe. Tasks simulated consist of pure and combined bending, torsion and compression loads in order to better compare model predictions with each other and with the published in vitro and in vivo data. The objective is to evaluate the predictive power of individual estimations versus the median of all estimations. It is hypothesized that the median predictions of FE models when combined could more closely approximate the experimental data than the predictions of individual models.
Section snippets
Inclusion criteria
Ten different research groups, working in the field of spinal FE modeling were invited to participate in the present study. Only validated models of the lumbar spine (L1-5) that were previously published in peer reviewed journals were considered. A model was considered to be validated when its predictions compared favorably with available measurements under simple loading conditions. From ten groups, eight agreed to participate, one declined due to lack of resources and one did not respond to
Participating groups
Seven of eight groups completed all calculations for the first part of this study. Due to resource limitations, one group only presented results under pure moments and not pure compression. Six of the eight groups participated in the second part as two groups did not participate due to resource limitations. One of the six participating groups was not able to deliver results for the load case upper body flexion, due to convergence problems.
Part 1 – pure moments and pure compression
Under pure moments, the median total L1-5 rotations of
Discussion
Over the last few decades, the finite element method has been used to investigate the biomechanical behavior of the lumbar spine. These FE models are usually based on only one specific or one idealized average subject with unique mechanical and geometrical characteristics. Thus, with a few exceptions (Little and Adam, 2013, Niemeyer et al., 2012), the effect of inter-subject variability in geometry has mostly not been accounted for in modeling efforts. In addition, the crucial role of
Conflict of interest statement
There are no conflicts of interests.
Acknowledgments
This study was partly supported by the Deutsche Forschungsgemeinschaft (DFG PU 510/2-1) and by the Bundesinstitut für Sportwissenschaft (MiSpEx – Network).
References (76)
- et al.
A comparative study of two trunk biomechanical models under symmetric and asymmetric loadings
J. Biomech.
(2010) - et al.
Computation of trunk equilibrium and stability in free flexion-extension movements at different velocities
J. Biomech.
(2008) - et al.
in vitro axial preload application during spine flexibility testing: towards reduced apparatus-related artefacts
J. Biomech.
(2000) - et al.
Euler stability of the human ligamentous lumbar spine. Part II: Experiment
Clin. Biomech.
(1992) - et al.
Optimised loads for the simulation of axial rotation in the lumbar spine
J. Biomech.
(2011) - et al.
A non-optimized follower load path may cause considerable intervertebral rotations
J. Biomech.
(2010) - et al.
Anatomic Facet Replacement System (AFRS) Restoration of Lumbar Segment Mechanics to Intact: A Finite Element Study and in vitro Cadaver Investigation
SAS J.
(2007) - et al.
Moment-rotation responses of the human lumbosacral spinal column
J. Biomech.
(2007) - et al.
Stepwise reduction of functional spinal structures increase vertebral translation and intradiscal pressure
J. Biomech.
(2007) - et al.
A new laser scanning technique for imaging intervertebral disc displacement and its application to modeling nucleotomy
Clin. Biomech. (Bristol, Avon)
(2008)
Stepwise reduction of functional spinal structures increase range of motion and change lordosis angle
J. Biomech.
Finite element analysis of the spine: towards a framework of verification, validation and sensitivity analysis
Med. Eng. Phys.
An accurate validation of a computational model of a human lumbosacral segment
J. Biomech.
Geometry strongly influences the response of numerical models of the lumbar spine--a probabilistic finite element analysis
J. Biomech.
How does the geometry affect the internal biomechanics of a lumbar spine bi-segment finite element model? Consequences on the validation process
J. Biomech.
Effects of degenerated intervertebral discs on intersegmental rotations, intradiscal pressures, and facet joint forces of the whole lumbar spine
Comput. Biol. Med.
Influence of material properties on the mechanical behaviour of the L5-S1 intervertebral disc in compression: a nonlinear finite element study
J. Biomed. Eng.
Realistic loading conditions for upper body bending
J. Biomech.
The use of surface strain data and a neural networks solution method to determine lumbar facet joint loads during in vitro spine testing
J. Biomech.
What have we learned from finite element model studies of lumbar intervertebral discs in the past four decades?
J. Biomech.
Nonlinear stress analysis of the whole lumbar spine in torsion--mechanics of facet articulation
J. Biomech.
Load-bearing and stress analysis of the human spine under a novel wrapping compression loading
Clin. Biomech.
Extr. Clin. Relev. Data Finite Elem. Simul.
Intradiscal pressure together with anthropometric data – a data set for the validation of models
Clin. Biomech. (Bristol, Avon)
Accuracy and repeatability of a new method for measuring facet loads in the lumbar spine
J Biomech
Influence of ligament stiffness on the mechanical behavior of a functional spinal unit
J. Biomech.
Influence of different artificial disc kinematics on spine biomechanics
Clin. Biomech.
Finite element application in implant research for treatment of lumbar degenerative disc disease
Med. Eng. Phys.
Verification, validation and sensitivity studies in computational biomechanics
Comput. Methods Biomech. Biomed. Eng.
Evaluating and managing acute low back pain in the primary care setting
J. Gen. Intern. Med.
Parametric convergence sensitivity and validation of a finite element model of the human lumbar spine
Comput. Methods Biomech. Biomed. Eng.
Development and Validation of a Three Dimensional High Resolution Nonlinear Finite Element Model of an L3/L4 Funcional Spinal Unit
Change of disc height, radial disc bulge, and intradiscal pressure from discectomy. An in vitro investigation on human lumbar discs
Spine
Stress analysis of the disc adjacent to interbody fusion in lumbar spine
Med. Eng. Phys.
Optimised in vitro applicable loads for the simulation of lateral bending in the lumbar spine
Med. Eng. Phys.
The role of disc, facets and fibres in degenerative process: a sensitivity study
Studies in Health Technology and Informatics
Finite element analysis in spine research. Proceedings of the Institution of Mechanical Engineers. Part H
J. Eng. Med.
Cited by (319)
Assessment of a fully-parametric thoraco-lumbar spine model generator with articulated ribcage
2024, Journal of BiomechanicsFinite element investigation of the influence of a new transpedicular vertebral implant positioning on biomechanical responses of the spine segment
2024, Medical Engineering and Physics