Application of a new calibration method for a three-dimensional finite element model of a human lumbar annulus fibrosus

Biomechanical evaluation of a spherical lumbar interbody device at varying levels of subsidence

Background

Ulf Fernström implanted stainless steel ball bearings following discectomy, or for painful disc disease, and termed this procedure disc arthroplasty. Today, spherical interbody spacers are clinically available, but there is a paucity of associated biomechanical testing. The primary objective of the current study was to evaluate the biomechanics of a spherical interbody implant. It was hypothesized that implantation of a spherical interbody implant, with combined subsidence into the vertebral bodies, would result in similar ranges of motion (RoM) and facet contact forces (FCFs) when compared with an intact condition. A secondary objective of this study was to determine the effect of using a polyetheretherketone (PEEK) versus a cobalt chrome (CoCr) implant on vertebral body strains. We hypothesized that the material selection would have a negligible effect on vertebral body strains since both materials have elastic moduli substantially greater than the annulus.

Methods

A finite element model of L3-L4 was created and validated by use of ROM, disc pressure, and bony strain from previously published data. Virtual implantation of a spherical interbody device was performed with 0, 2, and 4 mm of subsidence. The model was exercised in compression, flexion, extension, axial rotation, and lateral bending. The ROM, vertebral body effective (von Mises) strain, and FCFs were reported.

Results

Implantation of a PEEK implant resulted in slightly lower strain maxima when compared with a CoCr implant. For both materials, the peak strain experienced by the underlying bone was reduced with increasing subsidence. All levels of subsidence resulted in ROM and FCFs similar to the intact model.

Conclusions

The results suggest that a simple spherical implant design is able to maintain segmental ROM and provide minimal differences in FCFs. Large areas of von Mises strain maxima were generated in the bone adjacent to the implant regardless of whether the implant was PEEK or CoCr.

Does anterior lumbar interbody fusion promote adjacent degeneration in degenerative disc disease? A finite element study

BACKGROUND

The increase in the number of anterior lumbar interbody fusions being performed carries with it the potential for the long-term complication of adjacent segmental degeneration. While its exact mechanism remains uncertain, adjacent segment degeneration has become much more widespread. Using a nonlinear, three-dimensional finite element model to analyze and compare the biomechanical influence of anterior lumbar interbody fusion and lumbar disc degeneration on the superior adjacent intervertebral disc, we attempt to determine if anterior lumbar interbody fusion aggravates adjacent segment degeneration.

METHODS

A normal three-dimensional non-linear finite element model of L3-5 has been developed. Three different grades of disc degeneration models (mild, moderate, severe) and one anterior lumbar interbody fusion model were developed by changing either the geometry or associated material properties of the L4-5 segment. The 800 N pre-compressive loading plus 10 Nm moments simulating flexion, extension, lateral bending and axial rotation in five steps was imposed on the L3 superior endplate of all models. The intradiscal pressure, intersegmental rotation range and Tresca stress of the annulus fibrosus in the L3-4 segment were investigated.

RESULTS

The intradiscal pressure, intersegmental rotation range and Tresca stress of the L3-4 segment in the fusion model are higher than in the normal model and different degeneration models under all motion directions. The intradiscal pressures in the three degenerative models are higher than in the normal model in flexion, extension and lateral bending, whereas in axial rotation, the value of the mild degeneration model is lower. The intersegmental rotation ranges in the three degenerative models are higher than in the normal model in flexion and extension. The values for the mild degeneration model in lateral bending and all the degeneration models in axial rotation are lower than in the normal model. The Tresca stresses are higher in the three degenerative models than in the normal model.

CONCLUSION

Anterior lumbar interbody fusion has more adverse biomechanical influence than disc degeneration on the adjacent upper disc and may aggravate the adjacent upper segmental degeneration.

The micromechanical role of the annulus fibrosus components under physiological loading of the lumbar spine

To date, studies that have investigated the kinematics of spinal motion segments have largely focused on the contributions that the spinal ligaments play in the resultant motion patterns. However, the specific roles played by intervertebral disk components, in particular the annulus fibrosus, with respect to global motion is not well understood in spite of the relatively large literature base with respect to the local ex vivo mechanical properties of the tissue. The primary objective of this study was to implement the nonlinear and orthotropic mechanical behavior of the annulus fibrosus in a finite element model of an L4/L5 functional spinal unit in the form of a strain energy potential where the individual mechanical contributions of the ground substance and fibers were explicitly defined. The model was validated biomechanically under pure moment loading to ensure that the individual role of each soft tissue structure during load bearing was consistent throughout the physiologically relevant loading range. The fibrous network of the annulus was found to play critical roles in limiting the magnitude of the neutral zone and determining the stiffness of the elastic zone. Under flexion, lateral bending, and axial rotation, the collagen fibers were observed to bear the majority of the load applied to the annulus fibrosus, especially in radially peripheral regions where disk bulging occurred. For the first time, our data explicitly demonstrate that the exact fiber recruitment sequence is critically important for establishing the range of motion and neutral zone magnitudes of lumbar spinal motion segments.

FEM Simulation of Non-Progressive Growth from Asymmetric Loading and Vicious Cycle Theory: Scoliosis Study Proof of Concept

Scoliosis affects about 1-3% of the adolescent population, with 80% of cases being idiopathic. There is currently a lack of understanding regarding the biomechanics of scoliosis, current treatment methods can be further improved with a greater understanding of scoliosis growth patterns. The objective of this study is to develop a finite element model that can respond to loads in a similar fashion as current spine biomechanics models and apply it to scoliosis growth. Using CT images of a non-scoliotic individual, a finite element model of the L3-L4 vertebra was created. By applying asymmetric loading in accordance to the 'vicious cycle' theory and through the use of a growth modulation equation it is possible to determine the amount of growth each region of the vertebra will undergo; therefore predict scoliosis growth over a period of time. This study seeks to demonstrate how improved anatomy can expand researchers current knowledge of scoliosis.

Biomechanical modeling of the lateral decubitus posture during corrective scoliosis surgery

BACKGROUND

Patient prone positioning in scoliosis surgeries modifies the spinal curves prior to instrumentation. However, the biomechanical effects of the lateral decubitus posture, used in anterior approaches and minimally invasive techniques, have not yet been investigated. The objectives were to develop and validate a finite element model simulating the spinal changes resulting from this positioning.

METHODS

The 3D pre-op reconstructed geometries of six adolescent patients with idiopathic scoliosis were used to develop personalized finite element models of the spine, which integrated a three-step method simulating the lateral posture. Clinical indices were measured on pre- and intra-operative radiographs to validate the finite element model.

FINDINGS

The major Cobb angle and apical vertebral translation were reduced by 44% and 37% respectively between the pre- and intra-op postures. Using appropriately oriented gravity forces and boundary conditions, the finite element model simulations represented adequately these changes, with average differences of 4 degrees for the major Cobb angle and 4mm for the apical vertebral translation with the radiographic values.

INTERPRETATION

Lateral decubitus positioning significantly reduces the spinal deformities prior to instrumentation, as demonstrated by the finite element model. This study is a first step in the development of a modeling tool for the optimal adjustments of intra-operative positioning, which remains to be further investigated with complementary clinical studies.

Restoration of compressive loading properties of lumbar discs with a nucleus implant-a finite element analysis study

BACKGROUND CONTEXT

Discectomy is a common procedure for treating sciatica. However, both the operation and preceding herniated disc alter the biomechanical properties of the spinal segment. The disc mechanics are also altered in patients with chronic contained herniation. The biomechanical properties of the disc can potentially be restored with an elastomeric nucleus replacement implanted via minimally invasive surgery.

PURPOSE

The purpose of this study was to determine whether the compressive characteristics of the intervertebral disc after a nucleotomy can be restored with an elastomeric nucleus replacement.

STUDY DESIGN

A finite element model of the L4-L5 intervertebral disc was created to investigate the effect of the implantation of an elastomeric nucleus replacement on the biomechanical properties of the disc under axial loading.

METHOD

A L4-L5 physiologic intervertebral disc model was constructed and then modified to contain a range by volume of nucleotomies and nucleus replacements. The material properties of the nucleus replacement were based on experimental data for an elastomeric implant. The compressive stiffness, radial annular bulge, and stress distribution of the nucleotomy and nucleus replacement models were investigated under displacement-controlled loading.

RESULTS

Removal of nucleus pulposus from the physiologic disc reduced the force necessary to compress the disc 2 mm by 50%, altered the von Mises stress distribution, and reduced the outward radial annular bulge. Replacing the natural nucleus pulposus of the physiologic disc with an artificial nucleus reduced the force required to compress the disc 2 mm by 10%, indicating a restoration of disc compressive stiffness. The von Mises stress distribution and annular bulge observed in the disc with an artificial nucleus were similar to that observed in the physiologic disc.

CONCLUSION

This study demonstrates that despite having different material properties, a nucleus replacement implant can restore the axial compressive mechanical properties of a disc after a discectomy. The implant carries compressive load and transfers the load into annular hoop stress.

Response analysis of the lumbar spine during regular daily activities--a finite element analysis

A non-linear poroelastic finite element model of the lumbar spine was developed to investigate spinal response during daily dynamic physiological activities. Swelling was simulated by imposing a boundary pore pressure of 0.25 MPa at all external surfaces. Partial saturation of the disc was introduced to circumvent the negative pressures otherwise computed upon unloading. The loading conditions represented a pre-conditioning full day followed by another day of loading: 8h rest under a constant compressive load of 350 N, followed by 16 h loading phase under constant or cyclic compressive load varying in between 1000 and 1600 N. In addition, the effect of one or two short resting periods in the latter loading phase was studied. The model yielded fairly good agreement with in-vivo and in-vitro measurements. Taking the partial saturation of the disc into account, no negative pore pressures were generated during unloading and recovery phase. Recovery phase was faster than the loading period with equilibrium reached in only approximately 3h. With time and during the day, the axial displacement, fluid loss, axial stress and disc radial strain increased whereas the pore pressure and disc collagen fiber strains decreased. The fluid pressurization and collagen fiber stiffening were noticeable early in the morning, which gave way to greater compression stresses and radial strains in the annulus bulk as time went by. The rest periods dampened foregoing differences between the early morning and late in the afternoon periods. The forgoing diurnal variations have profound effects on lumbar spine biomechanics and risk of injury.

Effect of multilevel lumbar disc arthroplasty on spine kinematics and facet joint loads in flexion and extension: a finite element analysis

Total disc arthroplasty (TDA) has been successfully used for monosegmental treatment in the last few years. However, multi-level TDA led to controversial clinical results. We hypothesise that: (1) the more artificial discs are implanted, the stronger the increases in spinal mobility and facet joint forces in flexion and extension; (2) deviations from the optimal implant position lead to strong instabilities. A three-dimensional finite element model of the intact L1-L5 human lumbar spine was created. Additionally, models of the L1-L5 region implanted with multiple Charité discs ranging from two to four levels were created. The models took into account the possible misalignments in the antero-posterior direction of the artificial discs. All these models were exposed to an axial compression preload of 500 N and pure moments of 7.5 Nm in flexion and extension. For central implant positions and the loading case extension, a motion increase of 51% for two implants up to 91% for four implants and a facet force increase of 24% for two implants up to 38% for four implants compared to the intact spine were calculated. In flexion, a motion decrease of 5% for two implants up to 8% for four implants was predicted. Posteriorly placed implants led to a better representation of the intact spine motion. However, lift-off phenomena between the core and the implant endplates were observed in some extension simulations in which the artificial discs were anteriorly or posteriorly implanted. The more artificial discs are implanted, the stronger the motion increase in flexion and extension was predicted with respect to the intact condition. Deviations from the optimal implant position lead to unfavourable kinematics, to high facet joint forces and even to lift-off phenomena. Therefore, multilevel TDA should, if at all, only be performed in appropriate patients with good muscular conditions and by surgeons who can ensure optimal implant positions.

The effect of different design concepts in lumbar total disc arthroplasty on the range of motion, facet joint forces and instantaneous center of rotation of a L4-5 segment

Although both unconstrained and constrained core lumbar artificial disc designs are in clinical use, the effect of their design on the range of motion, center of rotations, and facet joint forces is not well understood. It is assumed that the constrained configuration causes a fixed center of rotation with high facet forces, while the unconstrained configuration leads to a moving center of rotation with lower loaded facets. The authors disagree with both assumptions and hypothesized that the two different designs do not lead to substantial differences in the results. For the different implant designs, a three-dimensional finite element model was created and subsequently inserted into a validated model of a L4-5 lumbar spinal segment. The unconstrained design was represented by two implants, the Charité disc and a newly developed disc prosthesis: Slide-Disc. The constrained design was obtained by a modification of the Slide-Disc whereby the inner core was rigidly connected to the lower metallic endplate. The models were exposed to an axial compression preload of 1,000 N. Pure unconstrained moments of 7.5 Nm were subsequently applied to the three anatomical main planes. Except for extension, the models predicted only small and moderate inter-implant differences. The calculated values were close to those of the intact segment. For extension, a large difference of about 45% was calculated between both Slide-Disc designs and the Charité disc. The models predicted higher facet forces for the implants with an unconstrained core compared to an implant with a constrained core. All implants caused a moving center of rotation. Except for axial rotation, the unconstrained and constrained configurations mimicked the intact situation. In axial rotation, only the Slide- Disc with mobile core reproduced the intact behavior. Results partially support our hypothesis and imply that different implant designs do not lead to strong differences in the range of motion and the location of center of rotations. In contrast, facet forces appeared to be strongly dependent on the implant design. However, due to the great variability in facet forces reported in the literature, together with our results, we could speculate that these forces may be more dependent on the individual spine geometry rather than a specific implant design.

A new approach for assigning bone material properties from CT images into finite element models

Generation of subject-specific finite element (FE) models from computed tomography (CT) datasets is of significance for application of the FE analysis to bone structures. A great challenge that remains is the automatic assignment of bone material properties from CT Hounsfield Units into finite element models. This paper proposes a new assignment approach, in which material properties are directly assigned to each integration point. Instead of modifying the dataset of FE models, the proposed approach divides the assignment procedure into two steps: generating the data file of the image intensity of a bone in a MATLAB program and reading the file into ABAQUS via user subroutines. Its accuracy has been validated by assigning the density of a bone phantom into a FE model. The proposed approach has been applied to the FE model of a sheep tibia and its applicability tested on a variety of element types. The proposed assignment approach is simple and illustrative. It can be easily modified to fit users' situations.

Effect of nucleus replacement device properties on lumbar spine mechanics

STUDY DESIGN

A validated nonlinear three-dimensional finite element model of a single lumbar motion segment (L3-L4) was used to evaluate a range of moduli for ideally conforming nucleus replacement devices.

OBJECTIVE

The objective of the current study was to determine the biomechanical effects of nucleus replacement technology for a variety of implant moduli. We hypothesized that there would be an optimal modulus for a nucleus replacement that would provide loading in the surrounding bone and anulus similar to the intact state.

SUMMARY OF BACKGROUND DATA

Nucleus pulposus replacements are interventional therapies that restore stiffness and height to mildly degenerated intervertebral discs. Currently a wide variety of nucleus replacement technologies with a large range of mechanical properties are undergoing preclinical testing.

METHODS

A finite element model of L3-L4 was created and validated using range of motion, disc pressure, and bony strains from previously published data. The intact model was altered by changing the mechanical properties of the nucleus pulposus to represent a wide range of nucleus replacement technologies (E = 0.1, 1, 4, and 100 MPa). All of the models were exercised in compression, flexion, extension, lateral bending, and axial rotation. Vertebral body strain, peak anulus fibrosus shear strain, initial bone remodeling stimulus, range of motion, and center of rotation were analyzed.

RESULTS

A nucleus replacement modulus of 1 and 4 MPa resulted in vertebral body strains similar to the intact model. The softest device indicated increased loading in the AF and bone resorption adjacent to the implant. Areas of strain maxima and bone formation were observed adjacent to the implant for the stiffest device.

CONCLUSION

The current study predicted an optimal nucleus replacement of 1 to 4 MPa. An overly stiff implant could result in subsidence, which would preclude the benefit of disc height increase or restoration. Conversely, an overly soft implant could accelerate a degenerative cascade in the anulus.

Dependency of disc degeneration on shear and tensile strains between annular fiber layers for complex loads

BACKGROUND

One of the first signs of disc degeneration is the formation of circumferential tears within the annulus fibrosus. It is assumed that high shear and tensile strains between the lamellae mainly cause the initiation of these failures. However, it is not known which load application and which degree of disc degeneration could lead to the highest strains and therefore, might induce the formation of tears. Therefore, the aim of this finite element (FE) study was, to find load combinations that would yield highest shear and tensile strains in differently degenerated discs.

MATERIALS AND METHODS

A three-dimensional FE-model of a motion segment L4-5 was utilized in different degrees of disc degeneration (healthy, mild, moderate, and severe). The degenerated models consider the reduction of disc height, endplate curvatures, the osteophyte formation, the increase of nucleus compressibility, and the decrease of fiber and ligament stiffness. An axial compression load of 500 N together with moments of 7.5 Nm in single and combined load directions were simulated.

RESULTS

High strains for the healthy and degenerated discs were predicted for load combinations, particularly for the combination of lateral bending plus flexion or extension. The maximum strains were located in the postero-lateral region of the disc. In comparison to the healthy disc, the maximum strains increased slightly for the mildly and moderately degenerated disc. Strains decreased strongly for the severely degenerated disc. With progressive degeneration, the size of the region of maximum strains diminished and the location transferred from the inner annulus to the adjacent bony endplates.

CONCLUSIONS

The results could be a possible explanation for the initiation of circumferential tears. The mildly degenerated disc model, which represents early stages of life, suggests that circumferential tears could primarily occur at these stages, especially for the load combinations of lateral bending plus axial rotation and lateral bending plus flexion.

Prospective design delineation and subsequent in vitro evaluation of a new posterior dynamic stabilization system

STUDY DESIGN

Finite element and in vitro study.

OBJECTIVE

Finite element calculations to delineate a dynamic fixator and confirmation with an in vitro experiment.

SUMMARY AND BACKGROUND DATA

In the last few years, there was a paradigm shift from rigid to dynamic fixation of spinal segments. However, some so-called dynamic implants like the Dynesys performed still stiffer than anticipated. The aim of this study was to optimize a dynamic stabilization system.

METHODS

The development steps of this implant design can be summarized in a development loop. First, a finite element model of an intact human L4-L5 segment was used to delineate implant stiffness parameters for the implant, in consideration of clinical concerns and safety aspects. These data were used in a second step, leading to the final implant design. This development process was completed with an appropriate in vitro experiment. The optimal axial and bending stiffness were computed to reduce the spinal motion by 30%. For the validation process, in vitro tests were performed on 6 human lumbar spinal segments L2-L3 with a median age of 52. The model and the specimens were loaded with pure unconstrained moments of 7.5 Nm in flexion, extension, lateral bending, and axial rotation.

RESULTS

This study demonstrated the advantages of employing a finite element model for the implant parameter delineation. It was possible to prospectively outline the needed stiffness parameters for a desired spinal range of motion achievement.

CONCLUSION

In summary, FEM may accelerate the development and the realization of a new implant design.

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.

Interaction between finite helical axes and facet joint forces under combined loading

STUDY DESIGN

Finite element study.

OBJECTIVE

To investigate the interaction between the finite helical axis and facet joint loads under combined loading.

SUMMARY OF BACKGROUND DATA

Finite helical axes (FHA) in a functional spinal unit can indicate mechanical disorders and are relevant for the development of new arthoplasty techniques. The facet joints protect the intervertebral discs from excessive movements. The relationship between the FHAs and facet joint forces is not well-understood, because previous studies have separated both, spinal motion and facet forces.

METHODS

A finite element model of a lumbar spinal segment L4-L5 was used to simulate axial compression load of 500 N together with moments starting from 0 to 7.5 Nm in single anatomic main planes. Load combinations of 7.5 Nm were generated by changing the load direction in steps of 15 degrees between each pair of the 3 anatomic mainplanes.

RESULTS

For single axes loading, the FHAs were found to be in the center of the disc under small moments, independently from load directions. The facet joints were only slightly loaded. Higher moments increased the forces in facet joints up to 105 N in axial rotation, followed by extension (50 N) and lateral bending (36 N). Combined moments did not essentially increase the facet forces compared with the same moment applied in an anatomic main direction. High facet forces might have directed the FHAs to migrate posteriorly, especially for axial rotation. This situation resulted in FHAs outside the disc toward the compressed facet joint.

CONCLUSION

For clinical practice, patients immediately after the operation, or patients with facet joint arthritis should reduce or avoid axial rotation alone or in combination with other load applications, especially axial rotation plus lateral bending or flexion.

Which axial and bending stiffnesses of posterior implants are required to design a flexible lumbar stabilization system?

Dynamic stabilization devices have been introduced to clinics as an alternative to rigid fixation. The stiffness of these devices varies widely, whereas the optimal stiffness, achieving a predefined stabilization of the spine, is unknown. This study was focused on the determination of stiffness values for posterior stabilization devices achieving a flexible, semi-flexible or rigid connection between two vertebrae. An extensively validated finite element model of a lumbar spinal segment L4-5 with an implanted posterior fixation device was used in this study. The model was exposed to pure moments of 7.5 and 20Nm around the three principal anatomical directions, simulating flexion, extension, lateral bending and axial rotation. In parametrical studies, the influence of the axial and bending fixator stiffness on the spinal range of motion was investigated. In order to examine the validity of the computed results, an in-vitro study was carried out. In this, the influence of two posterior stabilization devices (DSS and rigidly internal fixator) on the segmental stabilization was investigated. The finite element (FE)-model predicted that each load direction caused a pairing of stiffness relations between axial and bending stiffness. In flexion and extension, however, the bending stiffness had a neglectable effect on the segmental stabilization, compared to the axial stiffness. In contrast, lateral bending and axial rotation were influenced by both stiffness parameters. Except in axial rotation, the model predictions were in a good agreement with the determined in-vitro data. In axial rotation, the FE-model predicted a stiffer segmental behavior than it was determined in the in-vitro study. It is usually expected that high stiffness values are required for a posterior stabilization device to stiffen a spinal segment. We found that already small stiffness values were sufficient to cause a stiffening. Using these data, it may possible to develop implants for certain clinical indications.

Finite element analysis of the spine: towards a framework of verification, validation and sensitivity analysis

A number of papers have recently emphasised the importance of verification, validation and sensitivity testing in computational studies within the field of biomechanical engineering. This review examines the methods used in the development of spinal finite element models with a view to a standardised framework of verification, validation and sensitivity analysis. The scope of this paper is restricted to models of the vertebra, the intervertebral disc and short spinal segments. In the case of single vertebral models, specimen-specific methods have been developed, which allow direct validation against experimental tests. The focus of intervertebral disc modelling has been on representing the complex material properties and further sensitivity testing is required to fully understand the relative roles of these input parameters. In order to construct complex multi-component short segment models, many geometric and material parameters are required, some of which are yet to be fully characterised. There are also major challenges in terms of short segment model validation. Throughout the review, areas of good practise are highlighted and recommendations for future development are proposed, taking a step towards more robust spinal modelling procedures, promoting acceptance from the wider biomechanics community.

Total disc replacement positioning affects facet contact forces and vertebral body strains

STUDY DESIGN

A validated nonlinear three-dimensional finite element (FE) model of a single lumbar motion segment (L3-L4) was used to evaluate the effects of total disc replacement (TDR). The model was implanted with a fixed-bearing TDR (ProDisc-L) at 2 surgically relevant positions and exercised about the 3 anatomic axes. Facet forces, range of motion (RoM), and vertebral body strains were evaluated.

OBJECTIVE

The objective of the current study was to evaluate how TDR implantation and positioning affects facet joint forces and vertebral body strains. We hypothesized that facet contact forces (FCFs) would increase with TDR to compensate for the loss of periprosthetic load-bearing structures, and that vertebral body strains would increase in the region around the metallic footplates.

SUMMARY OF BACKGROUND DATA

TDR has the potential to replace fusion as the gold standard for the treatment of painful degenerative disc disease. However, complications after TDR include index level facet arthrosis and implant subsidence. Alterations in facet and vertebral body loading after TDR and their dependence on implant positioning are not fully understood.

METHODS

An FEM of L3-L4 was created and validated using RoM, disc pressure, and bony strains from previously published data. A TDR was incorporated into the L3-L4 spine model. All models were subjected to a compressive follower load of 500 N and moments of 7.5 Nm about the 3 anatomic axes.

RESULTS

Overall RoM and FCFs tended to increase with TDR. FCFs increased by an order of magnitude during flexion. Posterior placement of the device resulted in an unloading of the facets during extension. Areas of strain maxima were observed in the anterior portion of the vertebral body during flexion after TDR. The area of initial bone resorption signal under the metal footplate was greater when the device was anteriorly placed.

CONCLUSION

The current study predicted a decrease in segmental rotational stiffness resulting from TDR. This resulted from the removal of load bearing soft tissue structures, and caused increased loading in the facets. Additionally, vertebral body strains were generally higher after TDR, and tended to increase with decreased rotational stiffness. Posterior placement of the device provided a more physiologic load transfer to the vertebral body.

A new laser scanning technique for imaging intervertebral disc displacement and its application to modeling nucleotomy

BACKGROUND

Nucleotomy is a standard procedure for treating disc prolapse. It can reduce intervertebral disc height, flattening and displacing the disc, which could lead to a painful narrowing of the foramina due to nerve root compression. The purpose of this study was to investigate the disc displacement of a complete spinal segment with and without nucleotomy. We hypothesized that a nucleotomy under a certain load combination might amplify disc displacement.

METHODS

A laser scanner was developed for recording three-dimensional disc displacement of six loaded L4-5 specimens for three conditions: intact, disc with vertebral bodies and subsequent nucleotomy. Specimens were exposed to pure moments of 7.5 N m in the three principal anatomical directions. Disc displacement was obtained at maximal deflection. A finite element model was validated and subsequently utilized to determine disc displacement. The task of the finite element model was to provide supplemental data for the posterolateral region, which could not be measured from intact specimens.

FINDINGS

Disc displacement measurements of intact specimens were limited to the anterior part of discs, whereas the finite element model was able to provide the missing data of the dorsal disc region. The simulation of load combinations showed that the highest disc displacement was 1.9 mm at the lateral or posterolateral region. The nucleotomy increased the disc displacement up to 2.1mm, whereas the displacement zenith migrated posterolaterally.

INTERPRETATION

These results could be a possible explanation for disadvantages of nucleotomy as a treatment. With the methodology presented here, we would be able to assess the performance of nucleus implants by determining the disc displacement map. This could also give us appropriate information of the annular deformation, which is needed for the development of motion preserving implants.

The relation between the instantaneous center of rotation and facet joint forces - A finite element analysis

BACKGROUND

The instantaneous center of rotation in a functional spinal unit is an indicator for mechanical disorders and is relevant for the development of motion preserving techniques. In addition to the intervertebral disc, the facet joints also play a major role for load transmission through the spine, providing stability to it. The relationship between the rotation center and facet joint forces is not fully understood, since previous studies have separated both; spinal motion and facet forces.

METHODS

A finite element model of a L4-5 lumbar spinal segment was exposed to an axial compression preload of 500 N. Pure unconstrained moments of 7.5 Nm were additionally applied in the three anatomical main planes. The instantaneous center of rotation and the facet joint forces were investigated.

FINDINGS

For small moments, the center of rotation was found to be almost in the center of the disc, no matter what motion direction. With an increasing flexion moment, the center of rotation moved anteriorly. The facet joints remained unloaded in flexion. With proceeding extension movement, the center of rotation moved posteriorly. The facet forces increased up to 50 N. In lateral bending, with increasing moment the center of rotation migrated posteriorly in the ipsilateral side of the disc. The forces in the facet joints rose to 36 N. In axial rotation, the center of rotation migrated towards the compressed facet joint with increasing moment. Axial rotation yielded the maximum facet forces with 105 N.

INTERPRETATION

The determination of the rotation center is highly sensible against measurement resolution obtained during in vivo and in vitro studies. This finite element method can be used to complement the knowledge of the rotation center location from former experimental findings.

The risk of disc prolapses with complex loading in different degrees of disc degeneration - a finite element analysis

BACKGROUND

Disc prolapses can result from various complex load situations and degenerative changes in the intervertebral disc. The aim of this finite element study was to find load combinations that would lead to the highest internal stresses in a healthy and in degenerated discs.

METHODS

A three-dimensional finite element model of a lumbar spinal segment L4-L5 in different grades of disc degeneration (healthy, mild, moderate, and severe) were generated, in which the disc height reduction, the formation of osteophytes and the increasing of nucleus' compressibility were considered. The intradiscal pressure in the nucleus, the fiber strains, and the shear strains between the annulus and the adjacent endplates under pure and complex loads were investigated.

RESULTS

In all grades of disc degeneration the intradiscal pressure was found to be highest in flexion. The shear and fiber strains predicted a strong increase under lateral bending+flexion for the healthy disc and under axial rotation and lateral bending+axial rotation for all degenerated discs, mostly located in the postero-lateral annulus. Compared to the healthy disc, the mildly degenerated disc indicated an increase of the intradiscal pressure and of the fiber strains, both of 25% in axial rotation. The shear strains showed an increase of 27% in axial rotation+flexion. As from the moderately degenerated disc all measurement parameters strongly decreased.

INTERPRETATION

The results support how specifically changes associated with disc degeneration might contribute to risk of prolapse. Thus, the highest risk of prolapses can be found for healthy and mildly degenerated discs.

Parametric equations to represent the profile of the human intervertebral disc in the transverse plane

Computational and finite element models of the spine are used to investigate spine and disc mechanics. Subject specific data for the transverse profile of the disc could improve the geometric accuracy of these models. The current study aimed to develop a mathematical algorithm to describe the profile of the disc components, using subject-specific data points. Using data points measured from pictures of human intervertebral discs sectioned in the transverse plane, parametric formulae were derived that mapped the outer profile of the anulus and nucleus. The computed anulus and nucleus profile were a similar shape to the discs from which they were derived. The computed total disc area was similar to the experimental data. The nucleus:disc area ratios were sensitive to the data points defined for each disc. The developed formulae can be easily implemented to provide patient specific data for the disc profile in computational models of the spine.

Creep associated changes in intervertebral disc bulging obtained with a laser scanning device

BACKGROUND

Lumbar disc bulging has been determined with different methods in the past. Reported methods of bulging assessment were limited to a direct physical contact, were two-dimensional and were time-consuming. Assessing the three-dimensional contour of a biological object under load would imply that the tissue would creep and therefore changes its contour. For that purpose, we were interested how fast the contour has to be assessed and how creeping would counteract on the intradiscal pressure and disc height.

METHODS

For that purpose, a laser based three-dimensional contour scanner was developed. This scanner was especially designed to be mounted in a spine tester. For 15 min a static compression of 500 N was applied to seven human lumbar segments having all ligaments, facets and arches removed. Disc height, intradiscal pressure and disc contour were time dependently measured.

FINDINGS

Load application reduced the disc height by 1.14 mm. The further decrease showed a typical creep behavior whereas the intradiscal pressure slightly but significantly decreased from 0.49 to 0.48 MPa. Cross-sectional disc contours showed that bulging was largest anterolateral followed by the anterior region. The creeping also increased the disc circumference. This effect varied region dependently having a maximum of 0.1 mm posterolateral.

INTERPRETATION

Results suggest that geometries of biological tissues should be obtained within one minute avoiding superimposing creep effects. This new method might be used to evaluate disc injuries, degeneration and disc treatments. Measuring disc contours under different loads and conditions yields the outer annular strain distribution. This is a prerequisite for the development of cell seeded and tissue engineered implants.

Application of a calibration method provides more realistic results for a finite element model of a lumbar spinal segment

BACKGROUND

An important step in finite element modeling is the process of validation to derive clinical relevant data. It can be assumed that defect states of a finite element model, which have not been validated before, may predict wrong results. The purpose of this study was to show the differences in accuracy between a calibrated and a non-calibrated finite element model of a lumbar spinal segment for different clinical defects.

METHODS

For this study, two geometrically identical finite element models were used. An in vitro experiment was designed, deriving data for the calibration. Frequently used material properties were obtained from the literature and transferred into the non-calibrated model. Both models were validated on three clinical defects: bilateral hemifacetectomy, nucleotomy and interspinous defects, whereas in vitro range of motion data served as control points. Predictability and accuracy of the calibrated and non-calibrated finite element model were evaluated and compared.

FINDINGS

Both finite element models could mimic the intact situation with a good agreement. In the defects, the calibrated model predicted motion behavior with excellent agreement, whereas the non-calibrated model diverged greatly.

INTERPRETATION

Investigating the biomechanical performance of implants and load distribution of different spinal structures by numerical analysis requires not only good agreement with the intact segment, but also with the defect states, which are initiated prior to implant insertion. Because of more realistic results the calibration method may be recommended, however, it is more time consuming.

Intradiscal pressure, shear strain, and fiber strain in the intervertebral disc under combined loading

STUDY DESIGN

Finite element study.

OBJECTIVE

To investigate intradiscal pressure, shear strain between anulus and adjacent endplates, and fiber strain in the anulus under pure and combined moments.

SUMMARY OF BACKGROUND DATA

Concerning anulus failures such as fissures and disc prolapses, the mechanical response of the intervertebral disc during combined load situations is still not well understood.

METHODS

A 3-dimensional, nonlinear finite element model of a lumbar spinal segment L4-L5 was used. Pure unconstraint moments of 7.5 Nm in all anatomic planes with and without an axial preload of 500 N were applied to the upper vertebral body. The load direction was incrementally changed with an angle of 15 degrees between the 3 anatomic planes to realize not only moments in the principle motion planes but also moment combinations.

RESULTS

Intradiscal pressure was highest in flexion and lowest in lateral bending. Load combinations did not increase the pressure. A combination of lateral bending plus flexion or lateral bending plus extension strongly increased the maximum shear strains. Lateral bending plus axial rotation yielded the highest increase in fiber strains, followed by axial rotation plus flexion or axial rotation plus extension. The highest shear and fiber strains were both located posterolaterally. An additional axial preload tended to increase the pressure, the shear, and fiber strains essentially for all load scenarios.

CONCLUSIONS

Combined moments seem to lead to higher stresses in the disc, especially posterolaterally. This region might be more susceptible to disc failure and prolapses. These results may help clinicians better understand the mechanical causes of disc prolapses and may also be valuable in developing preventive clinical strategies and postoperative treatments.

Nonlinear finite element analysis of anular lesions in the L4/5 intervertebral disc

Degenerate intervertebral discs exhibit both material and structural changes. Structural defects (lesions) develop in the anulus fibrosus with age. While degeneration has been simulated in numerous previous studies, the effects of structural lesions on disc mechanics are not well known. In this study, a finite element model (FEM) of the L4/5 intervertebral disc was developed in order to study the effects of anular lesions and loss of hydrostatic pressure in the nucleus pulposus on the disc mechanics. Models were developed to simulate both healthy and degenerate discs. Degeneration was simulated with either rim, radial or circumferential anular lesions and by equating nucleus pressure to zero. The anulus fibrosus ground substance was represented as a nonlinear incompressible material using a second-order polynomial, hyperelastic strain energy equation. Hyperelastic material parameters were derived from experimentation on sheep discs. Endplates were assumed to be rigid, and annulus lamellae were assumed to be vertical in the unloaded state. Loading conditions corresponding to physiological ranges of rotational motion were applied to the models and peak rotation moments compared between models. Loss of nucleus pulposus pressure had a much greater effect on the disc mechanics than the presence of anular lesions. This indicated that the development of anular lesions alone (prior to degeneration of the nucleus) has minimal effect on disc mechanics, but that disc stiffness is significantly reduced by the loss of hydrostatic pressure in the nucleus. With the degeneration of the nucleus, the outer innervated anulus or surrounding osteo-ligamentous anatomy may therefore experience increased strains.

Stepwise reduction of functional spinal structures increase vertebral translation and intradiscal pressure

To date, there are only a few studies that provide data to efficiently calibrate finite element models for the spine due to its complexity. In a recent study, we quantified the range of motion rotation and the lordosis angle. This paper provides complementary results regarding two more parameters, intradiscal pressure and vertebral translation. All parameters were obtained as a function of stepwise anatomical reduction, loading direction and magnitude. Eight lumbar spinal segments (L4-5) with a median age of 52 years (38-59 years) and no signs of disc degeneration were used for the in vitro testing. A miniaturized pressure probe was implanted into the nucleus. An ultrasound-based motion-tracking system was employed to record spatial movements of several landmarks on the specimens. The center of L4, the anterior, posterior, left and right point of the lower endplate of L4 were digitized as landmarks and its translation was determined. Specimens were loaded with pure moments (1-10Nm) in the three principal anatomical planes at a loading rate of 1.0 degrees /s. Anatomy was stepwise reduced by cutting different ligaments, facet capsules and joints and removing nucleus. Translation analysis showed that the L4 center point had its largest displacement in sagittal direction and almost none vertically. Removal of the supra- and interspinous, flaval ligaments showed a slight increase and further removal of structures, a higher increase of translation. Axial rotation also was accompanied with L4 to elevate when torsion was applied. This effect was found to be larger with progressing defects. Nucleotomy exhibited the most unstable situation for specimens. Results of the intradiscal pressure indicated a large increase after removing the facet capsules and joints. Furthermore, it was found that intradiscal pressure correlated well with data of range of motion for rotation. Predicting and simulating clinical defects, surgical intervention or treatment methods requires a well performed calibration based on in vitro data, whereas it is important to adapt all including structures with as many known parameters as possible. Results provided by these studies may be used as a database for researchers aiming to calibrate or validate finite element models of L4-5 segments.