Finite-Element Analysis for Lumbar Interbody Fusion Under Axial Loading

Patterns of Vertebral Bone Marrow Edema in the Normal Healing Process of Lumbar Interbody Fusion: Baseline Data for Diagnosis of Pathological Events

STUDY DESIGN

Retrospective investigation using a prospectively collected database.

OBJECTIVE

To examine the appearance and characteristics of vertebral bone marrow edema (BME) in the normal healing of lumbar interbody fusion.

SUMMARY OF BACKGROUND DATA

Although BME in pathological spinal conditions has been well-documented, the patterns and characteristics of BME in the normal healing process of spinal fusion remains unexplored.

MATERIALS AND METHODS

We reviewed imaging from 225 patients with normal healing following posterior lumbar interbody fusion or transforaminal lumbar interbody fusion. BME was identified on magnetic resonance imaging at the third postoperative week and categorized with respect to its appearance, including assessment of area and extension within the relevant vertebrae.

RESULTS

Three hundred eighty-nine of the 450 instrumented vertebrae (86.4%) displayed evidence BME. All instances of BME were associated with the area of contact with the endplate. The average extent of BME was 32.7±1.0%. BME within normal healing following interbody fusion could be categorized into four types: no edema (13.6%), anterior corner (36.6%), around-the-cage focal (48.0%), and diffuse (1.8%). Anterior corner BME was significantly associated with instances of single cage placement than in dual cages (42.6% vs. 24.7%, P =0.0002). Single cages had a significantly higher rate of BME than dual cages (92.0% vs. 75.3%, P <0.0001). The extent of BME was significantly greater in the single cage cohort (36.9% vs. 24.2% in dual cages, P <0.0001).

CONCLUSIONS

This serves as the first study demonstrating the patterns of BME associated with normal healing following lumbar interbody fusion procedures. Anterior corner BME and around-the-cage focal BME were the most common patterns encountered, with diffuse BME a relatively rare pattern. These findings might contribute to the better differentiation of postoperative pathological events from normal healing following lumbar interbody fusion.

LEVEL OF EVIDENCE

4.

Selection of suitable pedicle screw for degenerated cortical and cancellous bone of human lumbar spine: A finite element study

Pedicular arthrodesis is the traditional procedure in terms of increase in the biomechanical stability with higher fixation rate. The current work aims to identify the effect of three spinal pedicle screws considering cortical and cancellous degeneracy condition. Lumbar section L2-L3 is utilized and various load and moment conditions were applied to depict the various biomechanical parameters for selection of suitable screw. Three dimensional model is considered in finite element analysis to identify the various responses of pedicle screw at bone screw juncture. Computed tomography (CT) images of a healthy male were considered to generate the finite element vertebral model. Generated intact model was further utilized to develop the other implanted models of degenerated cortical and cancellous bone models. The three fused instrumented models with different cortical and cancellous degeneracy conditions were analyzed in finite element analysis. The results were obtained as stress pattern at bone screw boundary and intervertebral disc stress. FE simulated results represents significant changes in the von Mises stress due to various load and moment conditions on degenerated bones during different body movement conditions. Results have shown that among all pedicle screws, the 6.0 mm diameter screw reflects very less stress values at the juncture. Multiple results on biomechanical aspects obtained during the FE study can be considered to design a new stabilization device and may be helpful to plan surgery of critical sections.

Lumbar model generator: a tool for the automated generation of a parametric scalable model of the lumbar spine

Low back pain is a major cause of disability and requires the development of new devices to treat pathologies and improve prognosis following surgery. Understanding the effects of new devices on the biomechanics of the spine is crucial in the development of new effective and functional devices. The aim of this study was to develop a preliminary parametric, scalable and anatomically accurate finite-element model of the lumbar spine allowing for the evaluation of the performance of spinal devices. The principal anatomical surfaces of the lumbar spine were first identified, and then accurately fitted from a previous model supplied by S14 Implants (Bordeaux, France). Finally, the reconstructed model was defined according to 17 parameters which are used to scale the model according to patient dimensions. The developed model, available as a toolbox named the lumbar model generator, enables generating a population of models using subject-specific dimensions obtained from data scans or averaged dimensions evaluated from the correlation analysis. This toolbox allows patient-specific assessment, taking into account individual morphological variation. The models have applications in the design process of new devices, evaluating the biomechanics of the spine and helping clinicians when deciding on treatment strategies.

Simplifying the human lumbar spine (L3/L4) material in order to create an elemental structure for the future modeling

The human lumbar spine incorporates the best joints in nature due to its optimal static and dynamic behavior against the internal and external loads. Developing an elemental structure based on this joint requires simplification in terms of the materials employed by keeping the mechanical and anatomical behaviors of the human lumbar spine. In the present study, the finite element (FE) of two motion segments of the human lumbar spine (L3/L4) was developed based on the CT scan data as the base for vertebrae geometry, verified geometry properties for another part of two motion segments, and combination of materials and loads obtained from the validated resources. Then, simplification occurred in four continuous steps such as omitting the annual fibers of annual matrix, representing the material of the annual matrix to the nucleus, demonstrating the material of annual matrix to the endplates too, and omitting the trabecular part of vertebrae. The present study aimed to propose the method for developing the basic structure of the human lumbar spine by simplifying its materials in the above-mentioned steps, analyzing the biomechanical effects of these four steps in terms of their internal and external responses, and validating the data obtained from the FE method. The validated simplified way introduced in this study can be used for future research by making implants, prosthesis, and modeling based on the human lumbar spine in other fields such as industrial design, building structures, or joints, which results in making the model easier, cheaper, and more effective.

Hybrid Rigid-Deformable Model for Prediction of Neighboring Intervertebral Disk Loads During Flexion Movement After Lumbar Interbody Fusion at L3-4 Level

Knowledge of spinal loads in neighboring disks after interbody fusion plays an important role in the clinical decision of this treatment as well as in the elucidation of its effect. However, controversial findings are still noted in the literature. Moreover, there are no existing models for efficient prediction of intervertebral disk stresses within annulus fibrosus (AF) and nucleus pulposus (NP) regions. In this present study, a new hybrid rigid-deformable modeling workflow was established to quantify the mechanical stress behaviors within AF and NP regions of the L1-2, L2-3, and L4-5 disks after interbody fusion at L3-4 level. The changes in spinal loads were compared with results of the intact model without interbody fusion. The fusion outcomes revealed maximal stress changes (10%) in AF region of L1-2 disk and in NP region of L2-3 disk. The minimal stress change (1%) is noted at the NP region of the L1-2 disk. The validation of simulation outcomes of fused and intact lumbar spine models against those of other computational models and in vivo measurements showed good agreements. Thus, this present study may be used as a novel design guideline for a specific implant and surgical scenario of the lumbar spine disorders.

Biomechanical comparison of posterior lumbar interbody fusion and transforaminal lumbar interbody fusion by finite element analysis

BACKGROUND

The transforaminal lumbar interbody fusion (TLIF) procedure may reduce many of the risks and limitations associated with posterior lumbar interbody fusion (PLIF). However, little is known about the biomechanical difference between PLIF and TLIF.

OBJECTIVE

To determine the biomechanical difference between PLIF and TLIF by finite-element analysis.

METHODS

Three validated finite-element models of L3-5 lumbar segment were created (intact model, PLIF model, and TLIF model). To analyze the biomechanics of these models, flexion, extension, rotation, and lateral bending moments of 7.5 N-m with a compressive preload of 400 N were imposed on the superior surfaces of the L3 vertebral body.

RESULTS

The range of motion at the L4-5 level of the PLIF and TLIF models decreased for all loading cases, compared with the intact model. Differences in the range of motion between PLIF and TLIF were not significant at less than 1 degree for all loading cases. The stress of the cage was found to be high in the PLIF model at the cage-endplate interface under all loading conditions. The stress exerted on the pedicle screw was greater in TLIF than PLIF. Particularly in flexion loading, the stress experienced by the pedicle screw in the TLIF model was 70.7% greater than that in the PLIF model.

CONCLUSION

The TLIF procedure increases the approximate biomechanical stability and reduces stress at the cage-endplate interface, except for a slight increase in screw stress. Clinically, the TLIF procedure may reduce many of the risks and limitations associated with PLIF and offer a useful alternative to the PLIF procedure.

Biomechanical evaluation of three surgical scenarios of posterior lumbar interbody fusion by finite element analysis

Background

For the treatment of low back pain, the following three scenarios of posterior lumbar interbody fusion (PLIF) were usually used, i.e., PLIF procedure with autogenous iliac bone (PAIB model), PLIF with cages made of PEEK (PCP model) or titanium (Ti) (PCT model) materiel. But the benefits or adverse effects among the three surgical scenarios were still not fully understood.

Method

Finite element analysis (FEA), as an efficient tool for the analysis of lumbar diseases, was used to establish a three-dimensional nonlinear L1-S1 FE model (intact model) with the ligaments of solid elements. Then it was modified to simulate the three scenarios of PLIF. 10 Nm moments with 400 N preload were applied to the upper L1 vertebral body under the loading conditions of extension, flexion, lateral bending and torsion, respectively.

Results

Different mechanical parameters were calculated to evaluate the differences among the three surgical models. The lowest stresses on the bone grafts and the greatest stresses on endplate were found in the PCT model. The PCP model obtained considerable stresses on the bone grafts and less stresses on ligaments. But the changes of stresses on the adjacent discs and endplate were minimal in the PAIB model.

Conclusions

The PCT model was inferior to the other two models. Both the PCP and PAIB models had their own relative merits. The findings provide theoretical basis for the choice of a suitable surgical scenario for different patients.

Optimised loads for the simulation of axial rotation in the lumbar spine

Simplified loading modes (pure moment, compressive force) are usually applied in the in vitro studies to simulate flexion-extension, lateral bending and axial rotation of the spine. The load magnitudes for axial rotation vary strongly in the literature. Therefore, the results of current investigations, e.g. intervertebral rotations, are hardly comparable and may involve unrealistic values. Thus, the question 'which in vitro applicable loading mode is the most realistic' remains open. A validated finite element model of the lumbar spine was employed in two sensitivity studies to estimate the ranges of results due to published load assumptions and to determine the input parameters (e.g. torsional moment), which mostly affect the spinal load and kinematics during axial rotation. In a subsequent optimisation study, the in vitro applicable loading mode was determined, which delivers results that fit best with available in vivo measurements. The calculated results varied widely for loads used in the literature with potential high deviations from in vivo measured values. The intradiscal pressure is mainly affected by the magnitude of the compressive force, while the torsional moment influences mainly the intervertebral rotations and facet joint forces. The best agreement with results measured in vivo were found for a compressive follower force of 720N and a pure moment of 5.5Nm applied to the unconstrained vertebra L1. The results reveal that in many studies the assumed loads do not realistically simulate axial rotation. The in vitro applicable simplified loads cannot perfectly mimic the in vivo situation. However, the optimised values lead to the best agreement with in vivo measured values. Their consequent application would lead to a better comparability of different investigations.

Load- and displacement-controlled finite element analyses on fusion and non-fusion spinal implants

This study used finite element (FE) analysis with the load-controlled method (LCM) and the displacement-controlled method (DCM) to examine motion differences at the implant level and adjacent levels between fusion and non-fusion implants. A validated three-dimensional intact (INT) L1—L5 FE model was used. At the L3—L4 level, the INT model was modified to surgery models, including the artificial disc replacement (ADR) of ProDisc II, and the anterior lumbar interbody fusion (ALIF) cage with pedicle screw fixation. The LCM imposed 10 N m moments of four physiological motions and a 150 N preload at the top of L1. The DCM process was in accordance with the hybrid testing protocol. The average percentage changes in the range of motion (ROM) for whole non-operated levels were used to predict adjacent level effects (ALE%). At the implant level, the ALIF model showed similar stability with both control methods. The ADR model using the LCM had a higher ROM than the model using the DCM, especially in extension and torsion. At the adjacent levels, the ALIF model increased ALE% (at least 17 per cent) using the DCM compared with the LCM. The ADR model had an ALE% close to that of the INT model, using the LCM (average within 6 per cent), while the ALE% decreased when using the DCM. The study suggests that both control methods can be adopted to predict the fusion model at the implant level, and similar stabilization characteristics can be found. The LCM will emphasize the effects of the non-fusion implants. The DCM was more clinically relevant in evaluating the fusion model at the adjacent levels. In conclusion, both the LCM and the DCM should be considered in numerical simulations to obtain more realistic data in spinal implant biomechanics.

Finite element application in implant research for treatment of lumbar degenerative disc disease

Surgical treatment for disc degeneration can be roughly grouped as fusion, disc replacement and dynamic stabilization. The clinical efficacy and biomechanical features of the implants used for disc degenerations can be evaluated through short- or long-term follow up observation, in vitro and in vivo experiments and computational simulations. Finite element models are already making an important contribution to our understanding of the spine and its components. Models are being used to reveal the biomechanical function of the spine and its behavior when healthy, diseased or damaged. They are also providing support in the design and application of spinal instrumentation. The article reviewed the most recent studies in the application of FE models that address the issue of implant research for treatment of low back pain. The published studies were grouped and reviewed thoroughly based on the function of implants investigated. The considerations of the finite element analysis in these studies were further discussed.

Biomechanics of the lumbar spine after dynamic stabilization

Target of the study was to predict the biomechanics of the instrumented and adjacent levels due to the insertion of the DIAM spinal stabilization system (Medtronic Ltd). For this purpose, a 3-dimensional finite element model of the intact L3/S1 segment was developed and subjected to different loading conditions (flexion, extension, lateral bending, axial rotation). The model was then instrumented at the L4/L5 level and the same loading conditions were reapplied. Within the assumptions of our model, the simulation results suggested that the implant caused a reduction in range of motion of the instrumented level by 17% in flexion and by 43% in extension, whereas at the adjacent levels, no significant changes were predicted. Numerical results in terms of intradiscal pressure, relative to the intact condition, predicted that the intervertebral disc at the instrumented level was unloaded by 27% in flexion, by 51% in extension, and by 6% in axial rotation, while no variations in pressure were caused by the device in lateral bending. At the adjacent levels, a change of relative intradiscal pressure was predicted in extension, both at the L3/L4 level, which resulted unloaded by 26% and at the L5/S1 level, unloaded by 8%. Furthermore, a reduction in terms of principal compressive stress in the annulus fibrosus of the L4/L5 instrumented level was predicted, as compared with the intact condition. These numerical predictions have to be regarded as a theoretical representation of the behavior of the spine, because any finite element model represents only a simplification of the real structure.

Influences of denucleation on contact force of facet joints under whole body vibration

To investigate the influence of the injured disc, frequency, load and damping on the facet contact forces of the low lumbar spine on the condition of whole body vibration, a detailed 3-D nonlinear finite element model was created based on the actual geometrical data of embalmed vertebrae of lumbar spine. The denucleation and facetectomy, together with removal of the capsular ligaments was employed to mimic the injury conditions of lumbar spine after surgery. The compression cyclic force was assumed to mimic the dynamic loads of transport vehicles. The results show that the high frequency vibration might increase both of the value and the vibration amplitude of facet contact forces of the lumbar spine under whole body vibration. The nucleus removal may increase significantly the facet contact forces. Although damping can decrease the vibration amplitude of facet contact forces for intact models, it has less effect on the vibration amplitude of facet contact force for the denucleated models. The denucleation of intervertebral discs is more harmful to the facet articulation on the condition of whole body vibration.

Biomechanical analysis of cages for posterior lumbar interbody fusion

Interbody fusions using intervertebral cages have become increasingly common in spinal surgery. Computational simulations were conducted in order to compare different cage designs in terms of their biomechanical interaction with the spinal structures. Differences in cage design and surgical technique may significantly affect the biomechanics of the fused spine segment, but little knowledge is available on this topic. In the present study, four 3D finite element models were developed, reproducing the human L4-L5 spinal unit in intact condition and after implantation of three different cage models. The intact model consisted of two vertebral bodies and relevant laminae, facet joints, main ligaments and disc. The instrumented models reproduced the post-operative conditions resulting after implant of the different cages. The three considered devices were hollow threaded titanium cages, the BAK (Zimmer Centerpulse, Warsaw, IN, USA), the Interfix and the Interfix Fly (both by Medtronic Sofamor Danek, Memphis, TN, USA). Simulations were run imposing various loading conditions, under a constant compressive preload. A great increase in the stiffness induced on the spinal segment by all cages was observed in all the considered loading cases. Stress distributions on the bony surface were evaluated and discussed. The differences observed between the biomechanics of the instrumented models were associated with the geometrical and surgical features of the devices.

Effects of fusion-bone stiffness on the mechanical behavior of the lumbar spine after vertebral body replacement

BACKGROUND

Implants for vertebral body replacement are often inserted together with an additional stabilizing implant, e.g. an internal fixation device. During implantation bone grafts or milled bone is normally added to the anterior implant. Little is known about the stiffening effect of this fusion-bone mass on the mechanical behavior of the corresponding bone region, including the load distribution between the different parts.

METHODS

A three-dimensional finite element model of the lumbar spine was created with a vertebral body replacement at L3, a paired internal fixation device between L2 and L4, and left anterolateral fusion bone. The elastic modulus of fusion bone was varied in discrete steps between 0 MPa and 10,000 MPa. The model was loaded to simulate standing, 20 degrees flexion, 15 degrees extension and 6 degrees axial rotation in the lumbar spine.

FINDINGS

The elastic modulus of fusion bone has a considerable effect on the compressive force on vertebral body replacement and fusion bone for all loading cases studied. For extension, it also affects intersegmental rotation, the force in the erector spinae muscle, the compressive force on the internal fixator and intradiscal pressure in the adjacent discs. The elastic modulus most strongly affects the different parameters at values between 0 MPa and 500 MPa.

INTERPRETATION

Adding bone mass during vertebral body replacement reduces the loads on the ventral implant for all loading cases studied but extension when the fusion-bone stiffens. This protects the implant from fatigue. The load on the fusion bone increases with increasing elastic modulus. Thus bone grafts should be used whenever possible.