Effect of disc degeneration on the mechanical behavior of the human lumbar spine: a probabilistic finite element study

Assessment of a fully-parametric thoraco-lumbar spine model generator with articulated ribcage

The present paper describes a novel user-friendly fully-parametric thoraco-lumbar spine CAD model generator including the ribcage, based on 22 independent parameters (1 posterior vertebral body height per vertebra + 4 sagittal alignment parameters, namely pelvic incidence, sacral slope, L1-L5 lumbar lordosis, and T1-T12 thoracic kyphosis). Reliable third-order polynomial regression equations were implemented in Solidworks to analytically calculate 56 morphological dependent parameters and to automatically generate the spine CAD model based on primitive geometrical features. A standard spine CAD model, representing the case-study of an average healthy adult, was then created and positively assessed in terms of spinal anatomy, ribcage morphology, and sagittal profile. The immediate translation from CAD to FEM for relevant biomechanical analyses was successfully demonstrated, first, importing the CAD model into Abaqus, and then, iteratively calibrating the constitutive parameters of one lumbar and three thoracic FSUs, with particular interest on the hyperelastic material properties of the IVD, and the spinal and costo-vertebral ligaments. The credibility of the resulting lumbo-sacral and thoracic spine FEM with/without ribcage were assessed and validated throughout comparison with extensive in vitro and in vivo data both in terms of kinematics (range of motion) and dynamics (intradiscal pressure) either collected under pure bending moments and complex loading conditions (bending moments + axial compressive force).

Annulus Calibration Increases the Computational Accuracy of the Lumbar Finite Element Model

STUDY DESIGN

Mechanical simulations.

OBJECTIVE

Inadequate calibration of annuli negatively affects the computational accuracy of finite element (FE) models. Specifically, the definition of annulus average radius (AR) does not have uniformity standards. Differences between the elastic moduli in the different layers and parts of the annulus were not fully calibrated when a linear elastic material is used to define its material properties. This study aims to optimize the computational accuracy of the FE model by calibrating the annulus.

METHODS

We calibrated the annulus AR and elastic modulus in our anterior-constructed lumbar model by eliminating the difference between the computed range of motion and that measured by in vitro studies under a flexion-extension loading condition. Multi-indicator validation was performed by comparing the computed indicators with those measured in in vitro studies. The computation time required for the different models has also been recorded to evaluate the computational efficiency.

RESULTS

The difference between computed and measured ROMs was less than 1% when the annulus AR and elastic modulus were calibrated. In the model validation process, all the indicators computed by the calibrated FE model were within ±1 standard deviation of the average values obtained from in vitro studies. The maximum difference between the computed and measured values was less than 10% under nearly all loading conditions. There is no apparent variation tendency for the computational time associated with different models.

CONCLUSION

The FE model with calibrated annulus AR and regional elastic modulus has higher computational accuracy and can be used in subsequent mechanical studies.

A swelling-based biphasic analysis on the quasi-static biomechanical behaviors of healthy and degenerative intervertebral discs

BACKGROUND AND OBJECTIVE

The degeneration of intervertebral discs is significantly dependent of the changes in tissue composition ratio and tissue structure. Up to the present, the effects of degeneration on the quasi-static biomechanical responses of discs have not been well understood. The goal of this study is to quantitatively analyze the quasi-static responses of healthy and degenerative discs.

METHODS

Four biphasic swelling-based finite element models are developed and quantitatively validated. Four quasi-static test protocols, including the free-swelling, slow-ramp, creep and stress-relaxation, are implemented. The double Voigt and double Maxwell models are further used to extract the immediate (or residual), short-term and long-term responses of these tests.

RESULTS

Simulation results show that both the swelling-induced pressure in the nucleus pulposus and the initial modulus decrease with degeneration. In the free-swelling test of discs possessing healthy cartilage endplates, simulation results show that over 80% of the total strain is contributed by the short-term response. The long-term response is dominant for discs with degenerated permeability in cartilage endplates. For the creep test, over 50% of the deformation is contributed by the long-term response. In the stress-relaxation test, the long-term stress contribution occupies approximately 31% of total response and is independent of degeneration. Both the residual and short-term responses vary monotonically with degeneration. In addition, both the glycosaminoglycan content and permeability affect the engineering equilibrium time constants of the rheologic models, in which the determining factor is the permeability.

CONCLUSIONS

The content of glycosaminoglycan in intervertebral soft tissues and the permeability of cartilage endplates are two critical factors that affect the fluid-dependent viscoelastic responses of intervertebral discs. The component proportions of the fluid-dependent viscoelastic responses depend also strongly on test protocols. In the slow-ramp test, the glycosaminoglycan content is responsible for the changes of the initial modulus. Since existing computational models simulate disc degenerations only by altering disc height, boundary conditions and material stiffness, the current work highlights the significance of biochemical composition and cartilage endplates permeability in the biomechanical behaviors of degenerated discs.

Biomechanical evaluation of a novel intervertebral disc repair technique for large box-shaped ruptures

Objective: The purpose of this study was to analyze the feasibility of repairing a ruptured intervertebral disc using a patch secured to the inner surface of the annulus fibrosus (AF). Different material properties and geometries for the patch were evaluated. Methods: Using finite element analysis, this study created a large box-shaped rupture in the posterior-lateral region of the AF and then repaired it with a circular and square inner patch. The elastic modulus of the patches ranged from 1 to 50 MPa to determine the effect on the nucleus pulposus (NP) pressure, vertical displacement, disc bulge, AF stress, segmental range of motion (ROM), patch stress, and suture stress. The results were compared against the intact spine to determine the most suitable shape and properties for the repair patch. Results: The intervertebral height and ROM of the repaired lumbar spine was similar to the intact spine and was independent of the patch material properties and geometry. The patches with a modulus of 2-3 MPa resulted in an NP pressure and AF stresses closest to the healthy disc, and produced minimal contact pressure on the cleft surfaces and minimal stress on the suture and patch of all models. Circular patches caused lower NP pressure, AF stress and patch stress than the square patch, but also caused greater stress on the suture. Conclusion: A circular patch with an elastic modulus of 2-3 MPa secured to the inner region of the ruptured annulus fibrosus was able to immediately close the rupture and maintain an NP pressure and AF stress similar to the intact intervertebral disc. This patch had the lowest risk of complications and produced the greatest restorative effect of all patches simulated in this study.

Biomechanical comparison of pedicle screw fixation strength among three different screw trajectories using single vertebrae and one-level functional spinal unit

Three key factors are responsible for the biomechanical performance of pedicle screw fixation: screw mechanical characteristics, bone quality and insertion techniques. To the best of the authors' knowledge, no study has directly compared the biomechanical performance among three trajectories, i.e., the traditional trajectory (TT), modified trajectory (MT) and cortical bone trajectory (CBT), in a porcine model. This study compared the pullout strength and insertion torque of three trajectory methods in single vertebrae, the pullout strength and fixation stiffness including flexion, extension, and lateral bending in a one-level instrumented functional spinal unit (FSU) that mimics the in vivo configuration were clarified. A total of 18 single vertebrae and 18 FSUs were randomly assigned into three screw insertion methods (n = 6 in each trajectory group). In the TT group, the screw converged from its entry point, passed completely inside the pedicle, was parallel to the superior endplate, was located in the superior third of the vertebral body and reached to at least the anterior third of the vertebral body. In the MT group, the convergent angle was similar to that of the TT method but directed caudally to the anterior inferior margin of the vertebral body. The results of insertion torque and pullout strength in single vertebrae were analyzed; in addition, the stiffness and pullout strength in the one-level FSU were also investigated. This study demonstrated that, in single vertebrae, the insertion torque was significantly higher in CBT groups than in TT and MT groups (p < 0.05). The maximal pullout strength was significantly higher in MT groups than in TT and CBT groups (p < 0.05). There was no significant difference in stiffness in the three motions among all groups. The maximal pullout strength in FSUs of MT and CBT groups were significantly higher than the TT groups (p < 0.05). We concluded that either MT or CBT provides better biomechanical performance than TT in single vertebrae or FSUs. The lack of significance of stiffness in FSUs among three methods suggested that MT or CBT could be a reasonable alternative to TT if the traditional trajectory was not feasible.

Computational lumbar spine models: A literature review

BACKGROUND

Computational spine models of various types have been employed to understand spine function, assess the risk that different activities pose to the spine, and evaluate techniques to prevent injury. The areas in which these models are applied has expanded greatly, potentially beyond the appropriate scope of each, given their capabilities. A comprehensive understanding of the components of these models provides insight into their current capabilities and limitations.

METHODS

The objective of this review was to provide a critical assessment of the different characteristics of model elements employed across the spectrum of lumbar spine modeling and in newer combined methodologies to help better evaluate existing studies and delineate areas for future research and refinement.

FINDINGS

A total of 155 studies met selection criteria and were included in this review. Most current studies use either highly detailed Finite Element models or simpler Musculoskeletal models driven with in vivo data. Many models feature significant geometric or loading simplifications that limit their realism and validity. Frequently, studies only create a single model and thus can't account for the impact of subject variability. The lack of model representation for certain subject cohorts leaves significant gaps in spine knowledge. Combining features from both types of modeling could result in more accurate and predictive models.

INTERPRETATION

Development of integrated models combining elements from different model types in a framework that enables the evaluation of larger populations of subjects could address existing voids and enable more realistic representation of the biomechanics of the lumbar spine.

The Biomechanical Effects of Different Bag-Carrying Styles on Lumbar Spine and Paraspinal Muscles: A Combined Musculoskeletal and Finite Element Study

OBJECTIVES

Bags such as handbags, shoulder bags, and backpacks are commonly used. However, it is difficult to assess the biomechanical effects of bag-carrying styles on the lumbar spine and paraspinal muscles using traditional methods. This study aimed to evaluate the biomechanical effects of bag-carrying styles on the lumbar spine.

METHODS

We developed a hybrid model that combined a finite element (FE) model of the lumbar spine and musculoskeletal models of three bag-carrying styles. The image data was collected from a 26-years-old, 176 cm and 70 kg volunteer. OpenSim and ABAQUS were used to do the musculoskeletal analysis and finite analysis. Paraspinal muscle force, intervertebral compressive force (ICF), and intervertebral shear force (ISF) on L1 were calculated and loaded into the FE model to assess the stress distribution on the lumbar spine.

RESULTS

Different paraspinal muscle activation occurred in the three bag-carrying models. The increase in the ICF generated by all three bags was greater than the bags' weights. The handbag produced greater muscle force, ICF, ISF, and peak stress on the nucleus pulposus than the backpack and shoulder bag of the same weight. Peak stress on the intervertebral discs in the backpack model and the L1-L4 segments of the shoulder bag model increased linearly with bag weight, and increased exponentially with bag weight in the handbag model.

CONCLUSION

Unbalanced bag-carrying styles (shoulder bags and handbags) led to greater muscle force, which generated greater ICF, ISF, and peak stress on the lumbar spine. The backpack produced the least burden on the lumbar spine and paraspinal muscles. Heavy handbags should be used carefully in daily life.

The role of nerve fibers and their neurotransmitters in regulating intervertebral disc degeneration

Intervertebral disc degeneration (IVDD) has been the major contributor to chronic lower back pain (LBP). Abnormal apoptosis, senescence, and pyroptosis of IVD cells, extracellular matrix (ECM) degradation, and infiltration of immune cells are the major molecular alternations during IVDD. Changes at tissue level frequently occur at advanced IVD tissue. Ectopic ingrowth of nerves within inner annulus fibrosus (AF) and nucleus pulposus (NP) tissue has been considered as the primary cause for LBP. Innervation at IVD tissue mainly included sensory and sympathetic nerves, and many markers for these two types of nerves have been detected since 1940. In fact, in osteoarthritis (OA), beyond pain transmission, the direct regulation of neuropeptides on functions of chondrocytes have attracted researchers' great attention recently. Many physical and pathological similarities between joint and IVD have shed us the light on the neurogenic mechanism involved in IVDD. Here, an overview of the advances in the nervous system within IVD tissue will be performed, with a discussion on in the role of nerve fibers and their neurotransmitters in regulating IVDD. We hope this review can attract more research interest to address neuromodulation and IVDD itself, which will enhance our understanding of the contribution of neuromodulation to the structural changes within IVD tissue and inflammatory responses and will help identify novel therapeutic targets and enable the effective treatment of IVDD disease.

Computational Modeling Intervertebral Disc Pathophysiology: A Review

Lower back pain is a medical condition of epidemic proportion, and the degeneration of the intervertebral disc has been identified as a major contributor. The etiology of intervertebral disc (IVD) degeneration is multifactorial, depending on age, cell-mediated molecular degradation processes and genetics, which is accelerated by traumatic or gradual mechanical factors. The complexity of such intertwined biochemical and mechanical processes leading to degeneration makes it difficult to quantitatively identify cause–effect relationships through experiments. Computational modeling of the IVD is a powerful investigative tool since it offers the opportunity to vary, observe and isolate the effects of a wide range of phenomena involved in the degenerative process of discs. This review aims at discussing the main findings of finite element models of IVD pathophysiology with a special focus on the different factors contributing to physical changes typical of degenerative phenomena. Models presented are subdivided into those addressing role of nutritional supply, progressive biochemical alterations stemming from an imbalance between anabolic and catabolic processes, aging and those considering mechanical factors as the primary source that induces morphological change within the disc. Limitations of the current models, as well as opportunities for future computational modeling work are also discussed.

Glycine-Serine-Threonine Metabolic Axis Delays Intervertebral Disc Degeneration through Antioxidant Effects: An Imaging and Metabonomics Study

Although intervertebral disc degeneration (IDD) can be described as different stages of change through biological methods, this long and complex process cannot be defined in stages by single or simple combination of biological techniques. Under the background of the development of nuclear magnetic resonance (NMR) technology and the emerging metabonomics, we based on animal models and expanded to the study of clinical human degeneration models. The characteristics of different stages of IDD were analyzed by omics. Omics imaging combined with histology, cytology, and proteomics was used for screening of the intervertebral disc (IVD) of research subjects. Furthermore, mass spectrometry nontargeted metabolomics was used to explore profile of metabolites at different stages of the IDD process, to determine differential metabolic pathways and metabolites. NMR spectroscopy was used to qualitatively and quantitatively identify markers of degeneration. NMR was combined with mass spectrometry metabolomics to explore metabolic pathways. Metabolic pathways were determined through protein molecular biology and histocytology of the different groups. Distinguishing advantages of magnetic resonance spectroscopy (MRS) for analysis of metabolites and effective reflection of structural integrity and water molecule metabolism through diffusion tensor imaging (DTI) were further used to verify the macrometabolism profile during degeneration. A corresponding model of in vitro metabolomics and in vivo omics imaging was established. The findings of this study show that a series of metabolic pathways associated with the glycine-serine-threonine (Gly-Ser-Thr) metabolic axis affects carbohydrate patterns and energy utilization efficiency and ultimately delays disc degeneration through antioxidant effects.

TELD with limited foraminoplasty has potential biomechanical advantages over TELD with large annuloplasty: an in-silico study

Background

Facetectomy, an important procedure in the in–out and out–in techniques of transforaminal endoscopic lumbar discectomy (TELD), is related to the deterioration of the postoperative biomechanical environment and poor prognosis. Facetectomy may be avoided in TELD with large annuloplasty, but iatrogenic injury of the annulus and a high grade of nucleotomy have been reported as risk factors influencing poor prognosis. These risk factors may be alleviated in TELD with limited foraminoplasty, and the grade of facetectomy in this surgery can be reduced by using an endoscopic dynamic drill.

Methods

An intact lumbo-sacral finite element (FE) model and the corresponding model with adjacent segment degeneration were constructed and validated to evaluate the risk of biomechanical deterioration and related postoperative complications of TELD with large annuloplasty and TELD with limited foraminoplasty. Changes in various biomechanical indicators were then computed to evaluate the risk of postoperative complications in the surgical segment.

Results

Compared with the intact FE models, the model of TELD with limited foraminoplasty demonstrated slight biomechanical deterioration, whereas the model of TELD with large annuloplasty revealed obvious biomechanical deterioration. Degenerative changes in adjacent segments magnified, rather than altered, the overall trends of biomechanical change.

Conclusions

TELD with limited foraminoplasty presents potential biomechanical advantages over TELD with large annuloplasty. Iatrogenic injury of the annulus and a high grade of nucleotomy are risk factors for postoperative biomechanical deterioration and complications of the surgical segment.

Biomechanical Investigation Between Rigid and Semirigid Posterolateral Fixation During Daily Activities: Geometrically Parametric Poroelastic Finite Element Analyses

While spinal fusion using rigid rods remains the gold standard treatment modality for various lumbar degenerative conditions, its adverse effects, including accelerated adjacent segment disease (ASD), are well known. In order to better understand the performance of semirigid constructs using polyetheretherketone (PEEK) in fixation surgeries, the objective of this study was to analyze the biomechanical performance of PEEK versus Ti rods using a geometrically patient-specific poroelastic finite element (FE) analyses. Ten subject-specific preoperative models were developed, and the validity of the models was evaluated with previous studies. Furthermore, FE models of those lumbar spines were regenerated based on postoperation images for posterolateral fixation at the L4–L5 level. Biomechanical responses for instrumented and adjacent intervertebral discs (IVDs) were analyzed and compared subjected to static and cyclic loading. The preoperative model results were well comparable with previous FE studies. The PEEK construct demonstrated a slightly increased range of motion (ROM) at the instrumented level, but decreased ROM at adjacent levels, as compared with the Ti. However, no significant changes were detected during axial rotation. During cyclic loading, disc height loss, fluid loss, axial stress, and collagen fiber strain in the adjacent IVDs were higher for the Ti construct when compared with the intact and PEEK models. Increased ROM, experienced stress in AF, and fiber strain at adjacent levels were observed for the Ti rod group compared with the intact and PEEK rod group, which can indicate the risk of ASD for rigid fixation. Similar to the aforementioned pattern, disc height loss and fluid loss were significantly higher at adjacent levels in the Ti rod group after cycling loading which alter the fluid–solid interaction of the adjacent IVDs. This phenomenon debilitates the damping quality, which results in disc disability in absorbing stress. Such finding may suggest the advantage of using a semirigid fixation system to decrease the chance of ASD.

Automated Pipeline to Generate Anatomically Accurate Patient-Specific Biomechanical Models of Healthy and Pathological FSUs

State-of-the-art preoperative biomechanical analysis for the planning of spinal surgery not only requires the generation of three-dimensional patient-specific models but also the accurate biomechanical representation of vertebral joints. The benefits offered by computational models suitable for such purposes are still outweighed by the time and effort required for their generation, thus compromising their applicability in a clinical environment. In this work, we aim to ease the integration of computerized methods into patient-specific planning of spinal surgery. We present the first pipeline combining deep learning and finite element methods that allows a completely automated model generation of functional spine units (FSUs) of the lumbar spine for patient-specific FE simulations (FEBio). The pipeline consists of three steps: (a) multiclass segmentation of cropped 3D CT images containing lumbar vertebrae using the DenseVNet network, (b) automatic landmark-based mesh fitting of statistical shape models onto 3D semantic segmented meshes of the vertebral models, and (c) automatic generation of patient-specific FE models of lumbar segments for the simulation of flexion-extension, lateral bending, and axial rotation movements. The automatic segmentation of FSUs was evaluated against the gold standard (manual segmentation) using 10-fold cross-validation. The obtained Dice coefficient was 93.7% on average, with a mean surface distance of 0.88 mm and a mean Hausdorff distance of 11.16 mm (N = 150). Automatic generation of finite element models to simulate the range of motion (ROM) was successfully performed for five healthy and five pathological FSUs. The results of the simulations were evaluated against the literature and showed comparable ROMs in both healthy and pathological cases, including the alteration of ROM typically observed in severely degenerated FSUs. The major intent of this work is to automate the creation of anatomically accurate patient-specific models by a single pipeline allowing functional modeling of spinal motion in healthy and pathological FSUs. Our approach reduces manual efforts to a minimum and the execution of the entire pipeline including simulations takes approximately 2 h. The automation, time-efficiency and robustness level of the pipeline represents a first step toward its clinical integration.

Biomechanical modelling of the facet joints: a review of methods and validation processes in finite element analysis

There is an increased interest in studying the biomechanics of the facet joints. For in silico studies, it is therefore important to understand the level of reliability of models for outputs of interest related to the facet joints. In this work, a systematic review of finite element models of multi-level spinal section with facet joints output of interest was performed. The review focused on the methodology used to model the facet joints and its associated validation. From the 110 papers analysed, 18 presented some validation of the facet joints outputs. Validation was done by comparing outputs to literature data, either computational or experimental values; with the major drawback that, when comparing to computational values, the baseline data was rarely validated. Analysis of the modelling methodology showed that there seems to be a compromise made between accuracy of the geometry and nonlinearity of the cartilage behaviour in compression. Most models either used a soft contact representation of the cartilage layer at the joint or included a cartilage layer which was linear elastic. Most concerning, soft contact models usually did not contain much information on the pressure-overclosure law. This review shows that to increase the reliability of in silico model of the spine for facet joints outputs, more needs to be done regarding the description of the methods used to model the facet joints, and the validation for specific outputs of interest needs to be more thorough, with recommendation to systematically share input and output data of validation studies.

The biomechanical effects of graded upper articular process arthroplasty on lumbar spine: A finite element study

BACKGROUND

Superior articular process arthroplasty is important for intervertebral foramen microscopy but may lead to spinal instability. Currently, there has been no relevant study in relation to the biomechanical analysis of superior articular process arthroplasty. Hence, this study is intended to verify biomechanical effects after unilateral S1 superior articular process arthroplasty.

METHODS

Eight finite element (FE) models of lumbosacral vertebrae (L4-S) were constructed, and the superior articular process formation was simulated with the help of Geomagic studio. Then, the models were imported into Nastran software after optimization. Normal load and appropriate torque were applied to simulate forward flexion, back extension, lateral flexion and lateral rotation. In the end, changes of lumbar range of motion (ROM) and structural stress were compared with those of normal model.

RESULTS

Compared with the normal model, formed from ventral to dorsal (Longitudinal), the larger motion of lumbar spine and the greater larger stress of articular process showed statistical significance (P < 0.05) in most of directions when the forming range was greater than 3/5. Formed from the apex to the base (transverse), the larger motion of lumbar spine and the greater stress of articular process showed statistical significance (P < 0.05) in most of directions when the forming range was great than 1/5.

CONCLUSION

When conducting unilateral S1 articular process arthroplasty from ventral to dorsal, the forming range is recommended to be less than 3/5 of the superior articular process. Notably, it is not advisable to form from the apex to the base.

The therapeutic effect of TBK1 in intervertebral disc degeneration via coordinating selective autophagy and autophagic functions

Introduction

While its innate immune function has been known, recent works of literature have focused on the role of Tank binding kinase 1 (TBK1) in regulating autophagy and it is unknown whether TBK1 protects against intervertebral disc degeneration (IVDD) through affecting autophagy.

Objectives

Here, we aim to explore whether TBK1 is implicated in the pathogenesis of IVDD, and investigated the potential mechanism.

Methods

Western blotting and immunohistochemistry were used to detect the TBK1 expression in human and rat NP tissue. After TBK1 overexpression in NP cells with lentivirus transfection, autophagic flux, apoptosis and senescence percentage were assessed. Si-RNA , a utophagy inhibitors and protein phosphatase inhibitors were applied to study the mechanism of autophagy regulation. In vivo study, we further evaluated the therapeutic action of lentivirus-TBK1(Lv-TBK1)injection in a rodent IVDD model.

Results

The TBK1 level was reduced in rat and human NP tissue. TBK1 overexpression protected against apoptosis and premature senescence. These functions of TBK1 were abolished by chloroquine-medicated autophagy inhibition.P-TBK1, an activation form of TBK, is involved in selective autophagy through directly phosphorylating P62 at Ser 403, and the activation of TBK1 is also dependent on Parkin manner. TBK1 also activated NPCs autophagy to relieve puncture injury in vivo.

Conclusion

We demonstrated that TBK1 overexpression attenuated senescence and apoptosis and promoted NPCs survival via upregulating autophagy. TBK1 represents a promising avenue for IVDD treatment.

The response surface method-genetic algorithm for identification of the lumbar intervertebral disc material parameters

Long-term compressive load on the lumbar intervertebral disc (IVD) might lead to lumbar IVD herniation. Exploring the material parameters of normal and degenerative enucleated IVDs is the basis for studying their mechanical behavior. According to the inverse analysis principle of the parameter estimation, an optimization method was proposed to identify the parameters of the porous material of the lumbar IVD based on finite element inverse analysis. The poroelastic finite element models were established in line with the compression creep experiment. The material parameters were combined by Box-Behnken design (BBD), and the response surface (RS) models were constructed using a quadratic polynomial with cross terms and optimized by genetic algorithm (GA). The results showed that the simulation result of the best material parameter combination had a good agreement with the experiment. Compared with the normal lumbar IVD, the elastic modulus and permeability decreased, and Poisson's ratio increased for the enucleated disk, resulting in a significant difference in mechanical properties. The algorithm used in this study can reduce the parameter identification error compared with only the RS method, and decrease the number of finite element simulations compared with only the GA.

Development and Validation of a Computationally Efficient Finite Element Model of the Human Lumbar Spine: Application to Disc Degeneration

INTRODUCTION

This study develops and validates an accurate, computationally efficient, 3-dimensional finite element model (FEM) of the human lumbar spine. Advantages of this simplified model are shown by its application to a disc degeneration study that we demonstrate is completed in one-sixth the time required when using more complicated computed tomography (CT) scan-based models.

METHODS

An osseoligamentous FEM of the L1-L5 spine is developed using simple shapes based on average anatomical dimensions of key features of the spine rather than CT scan images. Pure moments of 7.5 Nm and a compressive follower load of 1000 N are individually applied to the L1 vertebra. Validation is achieved by comparing rotations and intradiscal pressures to other widely accepted FEMs and in vitro studies. Then degenerative disc properties are modeled and rotations calculated. Required computation times are compared between the model presented in this paper and other models developed using CT scans.

RESULTS

For the validation study, parameter values for a healthy spine were used with the loading conditions described above. Total L1-L5 rotations for flexion, extension, lateral bending, and axial rotation under pure moment loading were calculated as 20.3°, 10.7°, 19.7°, and 10.3°, respectively, and under a compressive follower load, maximum intradiscal pressures were calculated as 0.68 MPa. These values compare favorably with the data used for validation. When studying the effects of disc degeneration, the affected segment is shown to experience decreases in rotations during flexion, extension, and lateral bending (24%-56%), while rotations are shown to increase during axial rotation (14%-40%). Adjacent levels realize relatively minor changes in rotation (1%-6%). This parametric study required 17.5 hours of computation time compared to more than 4 days required if utilizing typical published CT scan-based models, illustrating one of the primary advantages of the model presented in this article.

CONCLUSIONS

The FEM presented in this article produces a biomechanical response comparable to widely accepted, complex, CT scan-based models and in vitro studies while requiring much shorter computation times. This makes the model ideal for conducting parametric studies of spinal pathologies and spinal correction techniques.