Lumbar spine finite element model for healthy subjects: development and validation

Finite element (FE) method is a proven powerful and efficient tool to study the biomechanics of the human lumbar spine. However, due to the large inter-subject variability of geometries and material properties in human lumbar spines, concerns existed on the accuracy and predictive power of one single deterministic FE model with one set of spinal geometry and material properties. It was confirmed that the combined predictions (median or mean value) of several distinct FE models can be used as an improved prediction of behavior of human lumbar spine under identical loading and boundary conditions. In light of this improved prediction, five FE models (L1-L5 spinal levels) of the human lumbar spine were developed based on five healthy living subjects with identical modeling method. The five models were extensively validated through experimental and computational results in the literature. Mesh convergence and material sensitivity analysis were also conducted. We have shown that the results from the five FE models developed in this paper were consistent with the experimental data and simulation results from the existing literature. The validated modeling method introduced in this study can be used in modeling dysfunctional lumber spines such as disc degeneration and scoliosis in future work.     

Twente spine model: A complete and coherent dataset for musculo-skeletal modeling of the lumbar region of the human spine

Musculo-skeletal modeling can greatly help in understanding normal and pathological functioning of the spine. For such models to produce reliable muscle and joint force estimations, an adequate set of musculo-skeletal data is necessary. In this study, we present a complete and coherent dataset for the lumbar spine, based on medical images and dissection measurements from one embalmed human cadaver. We divided muscles into muscle-tendon elements, digitized their attachments at the bones and measured morphological parameters. In total, we measured 11 muscles from one body side, using 96 elements. For every muscle element, we measured three-dimensional coordinates of its attachments, fiber length, tendon length, sarcomere length, optimal fiber length, pennation angle, mass, and physiological cross-sectional area together with the geometry of the lumbar spine. Results were consistent with other anatomical studies and included new data for the serratus posterior inferior muscle. The dataset presented in this paper enables a complete and coherent musculo-skeletal model for the lumbar spine and will improve the current state-of-the art in predicting spinal loading.     

Statistical shape modeling characterizes three-dimensional shape and alignment variability in the lumbar spine

The mechanics of the lumbar spine are heavily dependent on the underlying anatomy. Anatomical measures are used to assess the progression of pathologies related to low back pain and to screen patients for surgical treatment options. To describe anatomical norms and pathological differences for the population, statistical shape modeling, which uses full three-dimensional representations of bone morphology and relative alignment, can capture intersubject variability and enable comparative evaluations of subject to population. Accordingly, the objective of this study was to develop a comprehensive set of three-dimensional statistical models to characterize anatomical variability in the lumbar spine, by specifically describing the shape of individual vertebrae, and shape and alignment of the entire lumbar spine (L1-S1), with a focus on the L4-L5 and L5-S1 functional spinal units (FSU). Using CT scans for a cohort of 52 patients, lumbar spine geometries were registered to a template to establish correspondence and a principal component analysis identified the primary modes of variation. Scaling was the most prevalent mode of variation for all models. Subsequent modes of the statistical shape models of the individual bones characterized shape variation within the processes. Subsequent modes of variation for the FSU and entire spine models described alignment changes associated with disc height and lordosis. Quantification of anatomical variation in the spine with statistical models can inform implant design and sizing, assist clinicians in diagnosing pathologies, screen patients for treatment options, and support pre-operative planning.     

Three-dimensional static modeling of the lumbar spine

This paper presents three-dimensional static modeling of the human lumbar spine to be used in the formation of anatomically-correct movement patterns for a fully cable-actuated robotic lumbar spine which can mimic in vivo human lumbar spine movements to provide better hands-on training for medical students. The mathematical model incorporates five lumbar vertebrae between the first lumbar vertebra and the sacrum, with dimensions of an average adult human spine. The vertebrae are connected to each other by elastic elements, torsional springs and a spherical joint located at the inferoposterior corner in the mid-sagittal plane of the vertebral body. Elastic elements represent the ligaments that surround the facet joints and the torsional springs represent the collective effect of intervertebral disc which plays a major role in balancing torsional load during upper body motion and the remaining ligaments that support the spinal column. The elastic elements and torsional springs are considered to be nonlinear. The nonlinear stiffness constants for six motion types were solved using a multiobjective optimization technique. The quantitative comparison between the angles of rotations predicted by the proposed model and in the experimental data confirmed that the model yields angles of rotation close to the experimental data. The main contribution is that the new model can be used for all motions while the experimental data was only obtained at discrete measurement points.     

An improved multi-joint EMG-assisted optimization approach to estimate joint and muscle forces in a musculoskeletal model of the lumbar spine

Muscle force partitioning methods and musculoskeletal system simplifications are key modeling issues that can alter outcomes, and thus change conclusions and recommendations addressed to health and safety professionals. A critical modeling concern is the use of single-joint equilibrium to estimate muscle forces and joint loads in a multi-joint system, an unjustified simplification made by most lumbar spine biomechanical models. In the context of common occupational tasks, an EMG-assisted optimization method (EMGAO) is modified in this study to simultaneously account for the equilibrium at all lumbar joints (M-EMGAO). The results of this improved approach were compared to those of its conventional single-joint equivalent (S-EMGAO) counterpart, the latter method being applied to the same lumbar joints but one at a time. Despite identical geometrical configurations and passive contributions used in both models, computed outcomes clearly differed between single- and multi-joint methods, especially at larger trunk flexed postures and during asymmetric lifting. Moreover, muscle forces predicted by L5-S1 single-joint analyses do not maintain mechanical equilibrium at other spine joints crossed by the same muscles. Assuming that the central nervous system does not attempt to balance the external moments one joint at a time and that a given muscle cannot exert different forces at different joints, the proposed multi-joint method represents a substantial improvement over its single-joint counterpart. This improved approach, hence, resolves trunk muscle forces with biological integrity but without compromising mechanical equilibrium at the lumbar joints.     

A growth-based model for the prediction of fiber angle distribution in the intervertebral disc annulus fibrosus

Biomechanics and Modeling in Mechanobiology - There is a growing interest in the development of patient-specific finite element models of the human lumbar spine for both the assessment of injury...     

Biomechanical comparison of the effects of anterior,posterior and transforaminal lumbar interbody fusion on vibration characteristics of the human lumbar spine

Previous studies have compared the effects of different interbody fusion approaches on biomechanical responses of the lumbar spine to static loadings. However, very few have dealt with the whole body vibration (WBV) condition that is typically present in vehicles. This study was designed to determine the biomechanical differences among anterior, posterior and transforaminal lumbar interbody fusion (ALIF, PLIF and TLIF) under vertical WBV. A previously developed and validated finite element (FE) model of the intact L1–sacrum human lumbar spine was modified to simulate ALIF, PLIF and TLIF with bilateral pedicle screw fixation at L4–L5. Comparative studies on dynamic responses to the axial cyclic loading in these developed models were conducted. The results showed that at the fused L4–L5 level, dynamic responses of the von-Mises stress in L4 inferior and L5 superior endplates for the ALIF, PLIF and TLIF models were increased compared with the intact model. The endplate stresses in the TLIF model were lower than in the ALIF and PLIF models, but the TLIF generated greater stresses in the screws and rods compared with the ALIF and PLIF. At other levels, a decrease in dynamic responses of the disc bulge, annulus stress and intradiscal pressure was observed in all the fusion models compared with the intact one, but there was no obvious difference in these dynamic responses among the ALIF, PLIF and TLIF models. These findings might be useful in understanding vibration characteristics of the whole lumbar spine after different types of fusion surgery.     

Effect of the intra-abdominal pressure and the center of segmental body mass on the lumbar spine mechanics - a computational parametric study

Determination of physiological loads in human lumbar spine is critical for understanding the mechanisms of lumbar diseases and for designing surgical treatments. Computational models have been used widely to estimate the physiological loads of the spine during simulated functional activities. However, various assumptions on physiological factors such as the intra-abdominal pressure (IAP), centers of mass (COMs) of the upper body and lumbar segments, and vertebral centers of rotation (CORs) have been made in modeling techniques. Systematic knowledge of how these assumptions will affect the predicted spinal biomechanics is important for improving the simulation accuracy. In this paper, we developed a 3D subject-specific numerical model of the lumbosacral spine including T12 and 90 muscles. The effects of the IAP magnitude and COMs locations on the COR of each motion segment and on the joint/muscle forces were investigated using a global convergence optimization procedure when the subject was in a weight bearing standing position. The data indicated that the line connecting the CORs showed a smaller curvature than the lordosis of the lumbar spine in standing posture when the IAP was 0?kPa and the COMs were 10?mm anterior to the geometric center of the T12 vertebra. Increasing the IAP from 0 kPa to 10 kPa shifted the location of CORs toward the posterior direction (from 1.4?±?8.9 mm anterior to intervertebral disc (IVD) centers to 40.5?±?3.1 mm posterior to the IVD centers) and reduced the average joint force (from 0.78?±?0.11 Body weight (BW) to 0.31?±?0.07 BW) and overall muscle force (from 349.3?±?57.7 N to 221.5?±?84.2 N). Anterior movement of the COMs from -30 mm to 70 mm relative to the geometric center of T12 vertebra caused an anterior shift of the CORs (from 25.1?±?8.3 mm posterior to IVD centers to 7.8?±?6.2 mm anterior to IVD centers) and increases of average joint forces (from 0.78?±?0.1 BW to 0.93?±?0.1 BW) and muscle force (from 348.9?±?47.7 N to 452.9?±?58.6 N). Therefore, it is important to consider the IAP and correct COMs in order to accurately simulate human spine biomechanics. The method and results of this study could be useful for designing prevention strategies of spinal injuries and recurrences, and for enhancing rehabilitation efficiency.     

A model-based approach for estimation of changes in lumbar segmental kinematics associated with alterations in trunk muscle forces

The kinematics information from imaging, if combined with optimization-based biomechanical models, may provide a unique platform for personalized assessment of trunk muscle forces (TMFs). Such a method, however, is feasible only if differences in lumbar spine kinematics due to differences in TMFs can be captured by the current imaging techniques. A finite element model of the spine within an optimization procedure was used to estimate segmental kinematics of lumbar spine associated with five different sets of TMFs. Each set of TMFs was associated with a hypothetical trunk neuromuscular strategy that optimized one aspect of lower back biomechanics. For each set of TMFs, the segmental kinematics of lumbar spine was estimated for a single static trunk flexed posture involving, respectively, 40° and 10° of thoracic and pelvic rotations. Minimum changes in the angular and translational deformations of a motion segment with alterations in TMFs ranged from 0° to 0.7° and 0 mm to 0.04 mm, respectively. Maximum changes in the angular and translational deformations of a motion segment with alterations in TMFs ranged from 2.4° to 7.6° and 0.11 mm to 0.39 mm, respectively. The differences in kinematics of lumbar segments between each combination of two sets of TMFs in 97% of cases for angular deformation and 55% of cases for translational deformation were within the reported accuracy of current imaging techniques. Therefore, it might be possible to use image-based kinematics of lumbar segments along with computational modeling for personalized assessment of TMFs.     

A three-dimensional mathematical model of the thoracolumbar fascia and an estimate of its biomechanical effect

The thoracolumbar fascia (TLF) provides a means of attachment to the lumbar spine for several muscles including the transverse abdominis, and parts of the latissimus dorsi and internal oblique muscles. Previous biomechanical models of the lumbar spine either tend to omit the TLF on the assumption that its contribution would be negligible or incorporate only part of the TLF. Here, a three-dimensional model of the posterior and middle layers of the TLF is presented to enable its action to be included in future three-dimensional models of the spine. It is used illustratively to estimate the biomechanical influence of this structure on the lumbar spine. The formulation of the model allows the lines of action of the fibres comprising the fascia to be calculated for any posture whilst ensuring that anatomical constraints are satisfied. Application of the model suggests that the TLF produces moments primarily in flexion and extension. The simulated results demonstrate that the abdominal muscles, acting via the TLF, are capable of contributing extension moments comparable to those produced by other smaller muscles associated with the lumbar spine.     

Load-Relaxation Properties of the Human Trunk in Response to Prolonged Flexion: Measuring and Modeling the Effect of Flexion Angle

Experimental studies suggest that prolonged trunk flexion reduces passive support of the spine. To understand alterations of the synergy between active and passive tissues following such loadings, several studies have assessed the time-dependent behavior of passive tissues including those within spinal motion segments and muscles. Yet, there remain limitations regarding load-relaxation of the lumbar spine in response to flexion exposures and the influence of different flexion angles. Ten healthy participants were exposed for 16 min to each of five magnitudes of lumbar flexion specified relative to individual flexion-relaxation angles (i.e., 30, 40, 60, 80, and 100%), during which lumbar flexion angle and trunk moment were recorded. Outcome measures were initial trunk moment, moment drop, parameters of four viscoelastic models (i.e., Standard Linear Solid model, the Prony Series, Schapery''s Theory, and the Modified Superposition Method), and changes in neutral zone and viscoelastic state following exposure. There were significant effects of flexion angle on initial moment, moment drop, changes in normalized neutral zone, and some parameters of the Standard Linear Solid model. Initial moment, moment drop, and changes in normalized neutral zone increased exponentially with flexion angle. Kelvin-solid models produced better predictions of temporal behaviors. Observed responses to trunk flexion suggest nonlinearity in viscoelastic properties, and which likely reflected viscoelastic behaviors of spinal (lumbar) motion segments. Flexion-induced changes in viscous properties and neutral zone imply an increase in internal loads and perhaps increased risk of low back disorders. Kelvin-solid models, especially the Prony Series model appeared to be more effective at modeling load-relaxation of the trunk.     

Quantitative comparison of ligament formulation and pre-strain in finite element analysis of the human lumbar spine

Data has been published that quantifies the nonlinear, anisotropic material behaviour and pre-strain behaviour of the anterior longitudinal, supraspinous (SSL), and interspinous ligaments of the human lumbar spine. Additionally, data has been published on localized material properties of the SSL. These results have been incrementally incorporated into a previously validated finite element model of the human lumbar spine. Results suggest that the effects of increased ligament model fidelity on bone strain energy were moderate and the effects on disc pressure were slight, and do not justify a change in modelling strategy for most clinical applications. There were significant effects on the ligament stresses of the ligaments that were directly modified, suggesting that these phenomena should be included in FE models where ligament stresses are the desired metric.     

The viscoelastic standard nonlinear solid model: predicting the response of the lumbar intervertebral disk to low-frequency vibrations

Due to the mathematical complexity of current musculoskeletal spine models, there is a need for computationally efficient models of the intervertebral disk (IVD). The aim of this study is to develop a mathematical model that will adequately describe the motion of the IVD under axial cyclic loading as well as maintain computational efficiency for use in future musculoskeletal spine models. Several studies have successfully modeled the creep characteristics of the IVD using the three-parameter viscoelastic standard linear solid (SLS) model. However, when the SLS model is subjected to cyclic loading, it underestimates the load relaxation, the cyclic modulus, and the hysteresis of the human lumbar IVD. A viscoelastic standard nonlinear solid (SNS) model was used to predict the response of the human lumbar IVD subjected to low-frequency vibration. Nonlinear behavior of the SNS model was simulated by a strain-dependent elastic modulus on the SLS model. Parameters of the SNS model were estimated from experimental load deformation and stress-relaxation curves obtained from the literature. The SNS model was able to predict the cyclic modulus of the IVD at frequencies of 0.01 Hz, 0.1 Hz, and 1 Hz. Furthermore, the SNS model was able to quantitatively predict the load relaxation at a frequency of 0.01 Hz. However, model performance was unsatisfactory when predicting load relaxation and hysteresis at higher frequencies (0.1 Hz and 1 Hz). The SLS model of the lumbar IVD may require strain-dependent elastic and viscous behavior to represent the dynamic response to compressive strain.     

Functional implications of variation in lumbar vertebral count among hominins

As early as the 1970s, Robinson defined lumbar vertebrae according to their zygapophyseal orientation. He identified six lumbar elements in fossil Sts 14 Australopithecus africanus, one more than is commonly present in modern humans. It is now generally inferred that the modal number of lumbar vertebrae for australopiths and early Homo was six, from which the mode of five in later Homo is derived. The two central questions this study investigates are (1) to what extent do differences in human lumbar vertebral count affect lordotic shape and lumbar function, and (2) what does lumbar number variation imply about lumbar spine function in early hominins? To address these questions, I first outline a biomechanical model of lumbar number effect on lordotic function. I then identify relevant morphological differences in the human modal and extra-modal variants, which I use to test the model. These tests permit evaluation of the human L6 variant as a model for reconstructing early hominin modal number and spine function. Application of the biomechanical model in reconstructing australopith/early Homo lumbar spines highlights shared principles of Euler column strength and sagittal spine flexibility among early and modern hominins. Within modern humans, the extra-modal L6 variant has an extended series of three cranially positioned kyphotic vertebrae and strongly oblique zygapophyseal facets at the last lumbar level. Although they share the same radius and length of lumbar curvature, the L6 variant differs functionally from the L5 mode in its expanded range of sagittal flexion/extension and enhanced resistance to shear. Given the modal number of six lumbar vertebrae in australopiths and early Homo, lumbar spine mobility and strength would have been key properties of vertebral function in early bipeds whose upper and lower body segments were coupled by close approximation of the thorax and iliac crests.     

Computational study on the effect of loading alteration caused by disc degeneration on the trabecular architecture in human lumbar spine

A thorough understanding of age-related phenomena on the trabecular architecture in the human lumbar spine can help the diagnosis and prognosis of age-related architectural changes, and provide an insight into the corresponding clinical assessments. In this paper we considered the different loading conditions of the young and old lumbar spines mainly caused by disc degeneration and studied the effect of loading alteration on trabecular architecture in lumbar spines. A two-dimensional μFE models with a 40 μm pixel resolution were built to represent the full trabecular architecture in the human lumbar spine, and a topology optimization with the aid of finite element method was conducted to numerically investigate the trabecular morphological changes. Topology optimization iteratively distributes material in a design domain producing optimal layout or configuration, and it has been widely and successfully used for the study of bone remodeling. As a result of adaptive response of bone remodeling due to different loading conditions, we obtained two distinctively different trabecular architectures for the young and old lumbar spines, and we observed a strong correlation between our numerical results and the actual trabecular architecture in the literature. The proposed numerical framework and results demonstrated the potential use of the topology optimization-based numerical tool for putative treatments in advance of actual clinical procedures for the patients.     

Biomechanical design of load simulation in multiple spinal segments]

Multisegmental biomechanical studies on the lumbar spine are steadily increasing in importance. Only in this way can we acquire knowledge about the physiological behaviour of the entire lumbar spine. Furthermore, these studies allow us to analyse in vitro the biomechanics of manipulated lumbar spines after various surgical operations on the spine. A load simulator was developed to investigate multisegmental lumbar spine mobility, and its function was investigated in an initial study on 19 fresh--frozen specimens of human lumbar spine. After x-ray examination and determination of the bone mineral density, the specimens were loaded up to 10 Nm in the automatic electromechanical loading system under flexion/extension, lateral bending and axial rotation. An ultrasound-based motion analysis system was used to measure the displacements of the vertebrae involved.     

The effect of unstable loading versus unstable support conditions on spine rotational stiffness and spine stability during repetitive lifting

Lumbar spine stability has been extensively researched due to its necessity to facilitate load-bearing human movements and prevent structural injury. The nature of certain human movement tasks are such that they are not equivalent in levels of task-stability (i.e. the stability of the external environment). The goal of the current study was to compare the effects of dynamic lift instability, administered through both the load and base of support, on the dynamic stability (maximal Lyapunov exponents) and stiffness (EMG-driven model) of the lumbar spine during repeated sagittal lifts. Fifteen healthy males performed 23 repetitive lifts with varying conditions of instability at the loading and support interfaces. An increase in spine rotational stiffness occurred during unstable support scenarios resulting in an observed increase in mean and maximum Euclidean norm spine rotational stiffness (p=0.0011). Significant stiffening effects were observed in unstable support conditions about all lumbar spine axes with the exception of lateral bend. Relative to a stable control lifting trial, the addition of both an unstable load as well as an unstable support did not result in a significant change in the local dynamic stability of the lumbar spine (p=0.5592). The results suggest that local dynamic stability of the lumbar spine represents a conserved measure actively controlled, at least in part, by trunk muscle stiffening effects. It is evident therefore that local dynamic stability of the lumbar spine can be modulated effectively within a young-healthy population; however this may not be the case in a patient population.     

Effects of the rib cage on thoracic spine flexibility.

Besides protecting the internal organs of the thorax, the rib cage is the site of numerous muscle attachments. It also decreases the overall flexibility of the thoracic spine. This study developed finite element (FE) models of the thoracic spine with and without the rib cage, and the effects of the rib cage on thoracic spine flexibility were determined. The numerical models were validated by comparing the maximum rotation of the models for several loading cases with experimental data in the literature. After adapting the material properties for the discs and ligaments, the calculated maximum rotations differed from the measured median values by less than 1 degrees without the rib cage and by less than 2.5 degrees with it. The rib cage decreased the mean flexibility of the thoracic spine by 23% to 47%, depending on the loading plane. Assuming the ribs to be rigid beams required a corresponding reduction of ligament stiffnesses in order to achieve the same agreement of the maximum rotations with the measured median values. Interconnecting the FE thoracic spine model plus rib cage with the existing detailed FE lumbar spine model improves the simulation of force directions of muscles attached to the rib cage or thoracolumbar spine. In addition, such a model is suitable for determining the effects of lumbar spine implants on spinal balance.