An accurate validation of a computational model of a human lumbosacral segment

Comparative FEM study on intervertebral disc modeling: Holzapfel-Gasser-Ogden vs. structural rebars

Introduction: Numerical modeling of the intervertebral disc (IVD) is challenging due to its complex and heterogeneous structure, requiring careful selection of constitutive models and material properties. A critical aspect of such modeling is the representation of annulus fibers, which significantly impact IVD biomechanics. This study presents a comparative analysis of different methods for fiber reinforcement in the annulus fibrosus of a finite element (FE) model of the human IVD. Methods: We utilized a reconstructed L4-L5 IVD geometry to compare three fiber modeling approaches: the anisotropic Holzapfel-Gasser-Ogden (HGO) model (HGO fiber model) and two sets of structural rebar elements with linear-elastic (linear rebar model) and hyperelastic (nonlinear rebar model) material definitions, respectively. Prior to calibration, we conducted a sensitivity analysis to identify the most important model parameters to be calibrated and improve the efficiency of the calibration. Calibration was performed using a genetic algorithm and in vitro range of motion (RoM) data from a published study with eight specimens tested under four loading scenarios. For validation, intradiscal pressure (IDP) measurements from the same study were used, along with additional RoM data from a separate publication involving five specimens subjected to four different loading conditions. Results: The sensitivity analysis revealed that most parameters, except for the Poisson ratio of the annulus fibers and C01 from the nucleus, significantly affected the RoM and IDP outcomes. Upon calibration, the HGO fiber model demonstrated the highest accuracy (R2 = 0.95), followed by the linear (R2 = 0.89) and nonlinear rebar models (R2 = 0.87). During the validation phase, the HGO fiber model maintained its high accuracy (RoM R2 = 0.85; IDP R2 = 0.87), while the linear and nonlinear rebar models had lower validation scores (RoM R2 = 0.71 and 0.69; IDP R2 = 0.86 and 0.8, respectively). Discussion: The results of the study demonstrate a successful calibration process that established good agreement with experimental data. Based on our findings, the HGO fiber model appears to be a more suitable option for accurate IVD FE modeling considering its higher fidelity in simulation results and computational efficiency.

Study on the biodynamic characteristics and internal vibration behaviors of a seated human body under biomechanical characteristics

To provide reference and theoretical guidance for establishing human body dynamics models and studying biomechanical vibration behavior, this study aimed to develop and verify a computational model of a three-dimensional seated human body with detailed anatomical structure under complex biomechanical characteristics to investigate dynamic characteristics and internal vibration behaviors of the human body. Fifty modes of a seated human body were extracted by modal method. The intervertebral disc and head motions under uniaxial white noise excitation (between 0 and 20 Hz at 1.0, 0.5 and 0.5 m/s2 r.m.s. for vertical, fore-aft and lateral direction, respectively) were computed by random response analysis method. It was found that there were many modes of the seated human body in the low-frequency range, and the modes that had a great impact on seated human vibration were mainly distributed below 13 Hz. The responses of different positions of the spine varied greatly under the fore-aft and lateral excitation, but the maximum stress was distributed in the lumbar under different excitations, which could explain why drivers were prone to lower back pain after prolonged driving. Moreover, there was a large vibration coupling between the vertical and fore-aft direction of an upright seated human body, while the vibration couplings between the lateral and other directions were very small. Overall, the study could provide new insights into not only the overall dynamic characteristics of the human body, but also the internal local motion and biomechanical characteristics under different excitations.

Modeling the Effect of Annulus Fibrosus Stiffness on the Stressed State of a Vertebral L1 Body and Nucleus Pulposus

The investigation examines the transference of stiffness from intervertebral discs (IVDs) to the lumbar body of the L1 vertebra and the interactions among adjacent tissues. A computational model of the vertebra was developed, considering parameters such as cortical bone thickness, trabecular bone elasticity, and the nonlinear response of the nucleus pulposus to external loading. A nonlinear dynamic analysis was performed, revealing certain trends: a heightened stiffness of the annulus fibrosus correlates with a significant reduction in the vertebral body's ability to withstand external loading. At a supplied displacement of 6 mm, the vertebra with a degenerative disc reached its yielding point, whereas the vertebrae with a healthy annulus fibrosus exhibited a strength capacity exceeding 20%. The obtained findings and proposed methodology are potentially useful for biomedical engineers and clinical specialists in evaluating the condition of the annulus fibrosus and predicting its influence on the bone components of the spinal system.

Muscle-driven forward dynamic active hybrid model of the lumbosacral spine: combined FEM and multibody simulation

Most spine models belong to either the musculoskeletal multibody (MB) or finite element (FE) method. Recently, coupling of MB and FE models has increasingly been used to combine advantages of both methods. Active hybrid FE-MB models, still rarely used in spine research, avoid the interface and convergence problems associated with model coupling. They provide the inherent ability to account for the full interplay of passive and active mechanisms for spinal stability. In this paper, we developed and validated a novel muscle-driven forward dynamic active hybrid FE-MB model of the lumbosacral spine (LSS) in ArtiSynth to simultaneously calculate muscle activation patterns, vertebral movements, and internal mechanical loads. The model consisted of the rigid vertebrae L1-S1 interconnected with hyperelastic fiber-reinforced FE intervertebral discs, ligaments, facet joints, and force actuators representing the muscles. Morphological muscle data were implemented via a semi-automated registration procedure. Four auxiliary bodies were utilized to describe non-linear muscle paths by wrapping and attaching the anterior abdominal muscles. This included an abdominal plate whose kinematics was optimized using motion capture data from upper body movements. Intra-abdominal pressure was calculated from the forces of the abdominal muscles compressing the abdominal cavity. For the muscle-driven approach, forward dynamics assisted data tracking was used to predict muscle activation patterns that generate spinal postures and balance the spine without prescribing accurate spinal kinematics. During calibration, the maximum specific muscle tension and spinal rhythms resulting from the model dynamics were evaluated. To validate the model, load cases were simulated from -10° extension to +30° flexion with weights up to 20 kg in both hands. The biomechanical model responses were compared with in vivo literature data of intradiscal pressures, intra-abdominal pressures, and muscle activities. The results demonstrated high agreement with this data and highlight the advantages of active hybrid modeling for the LSS. Overall, this new self-contained tool provides a robust and efficient estimation of LSS biomechanical responses under in vivo similar loads, for example, to improve pain treatment by spinal stabilization therapies.

Biomechanical stability of oblique lateral interbody fusion combined with four types of internal fixations: finite element analysis

Objective: Using finite element analysis to identify the optimal internal fixation method for oblique lateral lumbar interbody fusion (OLIF), providing guidance for clinical practice. Methods: A finite element model of the L4 - L5 segment was created. Five types of internal fixations were simulated in the generated L4-L5 finite element (FE) model. Then, six loading scenarios, i.e., flexion, extension, left-leaning, right-leaning, rotate left, and rotate right, were simulated in the FE models with different types of fixations. The biomechanical stability of the spinal segment after different fixations was investigated. Results: Regarding the range of motion (ROM) of the fused segment, OLIF + Bilateral Pedicle Screws (BPS) has a maximum ROM of 1.82° during backward bending and the smallest ROM in all directions of motion compared with other models. In terms of the von Mises stress distribution on the cage, the average stress on every motion direction of OLIF + BPS is about 17.08MPa, and of OLIF + Unilateral Vertebral Screw - Pedicle Screw (UVS-PS) is about 19.29 MPa. As for the von Mises stress distribution on the internal fixation, OLIF + BPS has the maximum internal fixator stress in left rotation (31.85 MPa) and OLIF + Unilateral Pedicle Screw (UPS) has the maximum internal fixator stress in posterior extension (76.59 MPa). The data of these two models were smaller than those of other models. Conclusion: OLIF + BPS provides the greatest biomechanical stability, OLIF + UPS has adequate biomechanical stability, OLIF + UVS-PS is inferior to OLIF + UPS synthetically, and OLIF + Double row vertical screw (DRVS) and Individual OLIF (IO) do not present significant obvious advantages.

Influence of foot excitation and shin posture on the vibration behavior of the entire spine inside a seated human body

Due to ethical issues and simplification of traditional biomechanical models, experimental methods and traditional computer methods were difficult to quantify the effects of foot excitation and shin posture on vibration behavior of the entire spine inside a seated human body under vertical whole-body vibration. This study developed and verified different three-dimensional (3D) finite element (FE) models of seated human body with detailed anatomical structure under the biomechanical characteristics to predict vibration behavior of the entire spine inside a seated human body with different foot excitation (with and without vibration) and shin posture (vertical and tilt posture). Random response analysis was performed to study the transmissibility of the entire spine to seat under vertical white noise excitation between 0 and 20 Hz at 0.5 m/s2 r.m.s. The results showed that although the foot excitation could reduce the fore-aft transmissibility in the cervical spine (23% reduction), it could significantly increase that in the lumbar spine (52% increase), which resulted in complex alternating stresses at lumbar spine and made the lumbar spine more vulnerable to injury in long-term vibration environment. Moreover, the shin tilt posture made the maximum fore-aft transmissibility in the lumbar spine move to the upper lumbar spine. The study provided new insights into the influence of foot excitation and shin posture on the vibration behavior of the entire spine inside a seated human body. Foot excitation exposed the lumbar spine to complex alternating stresses and made it more vulnerable to injury in long-term whole body vibration.

Toward a mechanically biocompatible intervertebral disc: Engineering of combined biomimetic annulus fibrosus and nucleus pulposus analogs

Intervertebral disc (IVD) degeneration and accompanying lower back pain impose global medical and societal challenges, affecting over 600 million people worldwide. The IVD complex fibrocartilaginous structure is responsible for the spine biomechanical function. The nucleus pulposus (NP), composed of swellable glycosaminoglycan (GAG), transfers compressive loads to the surrounding fiber-reinforced annulus fibrosus (AF) lamellae, which stretches under tension. Together, these substructures allow the IVD to withstand extremely high and complex loads. Key to mimic the complete disc must consider the properties of its substructures. This study presents three novel substructures-a biomimetic silk-reinforced composite lamella for the AF, a GAG analog for the NP, and a novel biomimetic combined AF-NP construct. The biomimetic AF demonstrates nonlinear, hyperelastic, and anisotropic behavior similar to the native human AF, while the NP analog demonstrates mechanical behavior similar to the human NP. The synergized biomimetic AF-NP demonstrates similar behavior to the unconfined NP, with significantly increased deformations indicating improved performance. Validation of the AF-NP construct mechanics using a finite element model yields results compatible with native human IVD under various physiological loadings. The ability of our AF-NP construct to mimic the native IVD offers a revolutionary concept for the potential development of a fully functional IVD.

Biomechanical and clinical studies on lumbar spine fusion surgery: a review

Low back pain is associated with degenerative disc diseases of the spine. Surgical treatment includes fusion and non-fusion types. The gold standard is fusion surgery, wherein the affected vertebral segment is fused. The common complication of fusion surgery is adjacent segment degeneration (ASD). The ASD often leads to revision surgery, calling for a further fusion of adjacent segments. The existing designs of nonfusion type implants are associated with clinical problems such as subsidence, difficulty in implantation, and the requirement of revision surgeries. Various surgical approaches have been adopted by the surgeons to insert the spinal implants into the affected segment. Over the years, extensive biomechanical investigations have been reported on various surgical approaches and prostheses to predict the outcomes of lumbar spine implantations. Computer models have been proven to be very effective in identifying the best prosthesis and surgical procedure. The objective of the study was to review the literature on biomechanical studies for the treatment of lumbar spinal degenerative diseases. A critical review of the clinical and biomechanical studies on fusion spine surgeries was undertaken. The important modeling parameters, challenges, and limitations of the current studies were identified, showing the future research directions.

Contribution of Shape Features to Intradiscal Pressure and Facets Contact Pressure in L4/L5 FSUs: An In-Silico Study

Finite element models (FEMs) of the spine commonly use a limited number of simplified geometries. Nevertheless, the geometric features of the spine are important in determining its FEM outcomes. The link between a spinal segment's shape and its biomechanical response has been studied, but the co-variances of the shape features have been omitted. We used a principal component (PCA)-based statistical shape modelling (SSM) approach to investigate the contribution of shape features to the intradiscal pressure (IDP) and the facets contact pressure (FCP) in a cohort of synthetic L4/L5 functional spinal units under axial compression. We quantified the uncertainty in the FEM results, and the contribution of individual shape modes to these results. This parameterisation approach is able to capture the variability in the correlated anatomical features in a real population and sample plausible synthetic geometries. The first shape mode ([Formula: see text]) explained 22.6% of the shape variation in the subject-specific cohort used to train the SSM, and had the largest correlation with, and contribution to IDP (17%) and FCP (11%). The largest geometric variation in ([Formula: see text]) was in the annulus-nucleus ratio.

A novel approach for tetrahedral-element-based finite element simulations of anisotropic hyperelastic intervertebral disc behavior

Intervertebral discs are microstructurally complex spinal tissues that add greatly to the flexibility and mechanical strength of the human spine. Attempting to provide an adjustable basis for capturing a wide range of mechanical characteristics and to better address known challenges of numerical modeling of the disc, we present a robust finite-element-based model formulation for spinal segments in a hyperelastic framework using tetrahedral elements. We evaluate the model stability and accuracy using numerical simulations, with particular attention to the degenerated intervertebral discs and their likely skewed and narrowed geometry. To this end, 1) annulus fibrosus is modeled as a fiber-reinforced Mooney-Rivlin type solid for numerical analysis. 2) An adaptive state-variable dependent explicit time step is proposed and utilized here as a computationally efficient alternative to theoretical estimates. 3) Tetrahedral-element-based FE models for spinal segments under various loading conditions are evaluated for their use in robust numerical simulations. For flexion, extension, lateral bending, and axial rotation load cases, numerical simulations reveal that a suitable framework based on tetrahedral elements can provide greater stability and flexibility concerning geometrical meshing over commonly employed hexahedral-element-based ones for representation and study of spinal segments in various stages of degeneration.

Stress analysis of the thoracolumbar junction in the process of backward fall: An experimental study and finite element analysis

The aim of the present study was to evaluate the biomechanical mechanism of injuries of the thoracolumbar junction by the methods of a backward fall simulation experiment and finite element (FE) analysis (FEA). In the backward fall simulation experiment, one volunteer was selected to obtain the contact force data of the sacrococcygeal region during a fall. Utilizing the fall data, the FEA simulation of the backward fall process was given to the trunk FE model to obtain the stress status of local bone structures of the thoracolumbar junction during the fall process. In the fall simulation test, the sacrococcygeal region of the volunteer landed first; the total impact time was 1.14±0.58 sec, and the impact force was up to 4,056±263 N. The stress of thoracic (T)11 was as high as 42 MPa, that of the posterior margin and the junction of T11 was as high as 70.67 MPa, and that of the inferior articular process and the superior articular process was as high as 128 MPa. The average stress of T12 and the anterior margin of lumbar 1 was 25 MPa, and that of the endplate was as high as 21.7 MPa, which was mostly distributed in the back of the endplate and the surrounding cortex. According to the data obtained from the fall experiment as the loading condition of the FE model, the backward fall process can be simulated to improve the accuracy of FEA results. In the process of backward fall, the front edge of the vertebral body and the root of vertebral arch in the thoracolumbar junction are stress concentration areas, which have a greater risk of injury.

A finite element model of the human lower thorax to pelvis spinal segment: Validation and modal analysis

BACKGROUND

Several finite element (FE) models have been developed to study the effects of vibration on human lumbar spine. However, the authors know of no published results so far that have proposed computed tomography-based FE models of whole lumbar spine including the pelvis to conduct dynamic analysis.

OBJECTIVE

To create and validate a three-dimensional ligamentous FE model of the human lower thorax to pelvis spinal segment (T12-Pelvis) and provide a detailed simulation environment to investigate the dynamic characteristics of the lumbar spine under whole body vibration (WBV).

METHODS

The T12-Pelvis model was generated based on volume reconstruction from computed tomography scans and validated against the published experimental data. FE modal analysis was implemented to predict dynamic characteristics associated with the first-order vertical resonant frequency and vibration mode of the model with upper body mass of 40 kg under WBV.

RESULTS

It was found that the current FE model was validated and corresponded closely with the published data. The obtained results from the modal analysis indicated that the first-order vertical resonant frequency of the T12-Pelvis model was 6.702 Hz, and the lumbar spine mainly performed vertical motion with a small anteroposterior motion. It was also found that shifting the upper body mass centroid onwards or rearwards from the normal upright sitting posture reduced the vertical resonant frequency.

CONCLUSIONS

These findings may be helpful to better understand vibration response of the human spine, and provide important information to minimize injury and discomfort for these WBV-exposed occupational groups.

A microstructure-based model for a full lamellar-interlamellar displacement and shear strain mapping inside human intervertebral disc core

The determinant role of the annulus fibrosus interlamellar zones in the intervertebral disc transversal and volumetric responses and hence on their corresponding three-dimensional conducts have been only revealed and appreciated recently. Their consideration in disc modeling strategies has been proven to be essential for the reproduction of correct local strain and displacement fields inside the disc especially in the unconstrained directions of the disc. In addition, these zones are known to be the starting areas of annulus fibrosus circumferential tears and disc delamination failure mode, which is often judged as one of the most dangerous disc failure modes that could evolve with time leading to disc hernia. For this latter reason, the main goal of the current contribution is to incorporate physically for the first time, the interlamellar zones, at the scale of a complete human lumbar intervertebral disc, in order to allow a correct local vision and replication of the different lamellar-interlamellar interactions and an identification of the interlamellar critical zones. By means of a fully tridimensional chemo-viscoelastic constitutive model, which we implemented into a finite element code, the physical, mechanical and chemical contribution of the interlamellar zones is added to the disc. The chemical-induced volumetric response is accounted by the model for both the interlamellar zones and the lamellae using experimentally-based fluid kinetics. Computational simulations are performed and critically discussed upon different simple and complex physiological movements. The disc core and the interlamellar zones are numerically accessed, allowing the observation of the displacement and shear strain fields that are compared to direct MRI experiments from the literature. Important conclusions about the correct lamellar-interlamellar-nucleus interactions are provided thanks to the developed model. The critical interlamellar spots with the highest delamination potentials are defined, analyzed and related to the local kinetics and microstructure.

Influence of nucleotomy on the load sharing in the spinal facet joint under the loading scenarios of different human postures

One-in-five people suffer from chronic low back pain (LBP). The incidence of this disease has doubled since 1950s and affects not only the elderly, but also the young population. However, the mechanism of LBP is still unknown. A possible location where the LBP may develop is the facet joint and it has been revealed that the intervertebral disc (IVD) nucleotomy may be a trigger for LBP. The aim of the present study was to investigate the influence of IVD nucleotomy on the load sharing in the spinal facet joint under the loading scenarios of different postures. Finite element (FE) models of the intact and nucleotomised L4 - L5 spinal segments were generated from the clinical CT images. Seven human postures, including upright, 5° extension, 5° flexion, ± 6° lateral bending and ± 2° axial rotation, were simulated. The resultant forces in the fact joint were compared between the intact and the nucleotomised cases. It was revealed that the IVD nucleotomy significantly increased the forces in the facet joints under the loading scenarios of upright, 5° extension and 5° flexion. The IVD nucleotomy increased the force in the ipsilateral facet joint but decreased the force on the contralateral side under the loading scenarios of ± 2° axial rotation. However, the IVD nucleotomy made little influence on the resultant forces in both facet joints in the postures of ± 6° lateral bending. In conclusion, the IVD nucleotomy can cause an increase in the overall force in the facet joint, and thus may serve as a possible explanation for the LBP and a main contributing factor for the pain complaints.

Calibration and validation of a novel hybrid model of the lumbosacral spine in ArtiSynth-The passive structures

In computational biomechanics, two separate types of models have been used predominantly to enhance the understanding of the mechanisms of action of the lumbosacral spine (LSS): Finite element (FE) and musculoskeletal multibody (MB) models. To combine advantages of both models, hybrid FE-MB models are an increasingly used alternative. The aim of this paper is to develop, calibrate, and validate a novel passive hybrid FE-MB open-access simulation model of a ligamentous LSS using ArtiSynth. Based on anatomical data from the Male Visible Human Project, the LSS model is constructed from the L1-S1 rigid vertebrae interconnected with hyperelastic fiber-reinforced FE intervertebral discs, ligaments, and facet joints. A mesh convergence study, sensitivity analyses, and systematic calibration were conducted with the hybrid functional spinal unit (FSU) L4/5. The predicted mechanical responses of the FSU L4/5, the lumbar spine (L1-L5), and the LSS were validated against literature data from in vivo and in vitro measurements and in silico models. Spinal mechanical responses considered when loaded with pure moments and combined loading modes were total and intervertebral range of motions, instantaneous axes and centers of rotation, facet joint contact forces, intradiscal pressures, disc bulges, and stiffnesses. Undesirable correlations with the FE mesh were minimized, the number of crisscrossed collagen fiber rings was reduced to five, and the individual influences of specific anatomical structures were adjusted to in vitro range of motions. Including intervertebral motion couplings for axial rotation and nonlinear stiffening under increasing axial compression, the predicted kinematic and structural mechanics responses were consistent with the comparative data. The results demonstrate that the hybrid simulation model is robust and efficient in reproducing valid mechanical responses to provide a starting point for upcoming optimizations and extensions, such as with active skeletal muscles.

Development of a multiscale model of the human lumbar spine for investigation of tissue loads in people with and without a transtibial amputation during sit-to-stand

Quantification of lumbar spine load transfer is important for understanding low back pain, especially among persons with a lower limb amputation. Computational modeling provides a helpful solution for obtaining estimates of in vivo loads. A multiscale model was constructed by combining musculoskeletal and finite element (FE) models of the lumbar spine to determine tissue loading during daily activities. Three-dimensional kinematic and ground reaction force data were collected from participants with ([Formula: see text]) and without ([Formula: see text]) a unilateral transtibial amputation (TTA) during 5 sit-to-stand trials. We estimated tissue-level load transfer from the multiscale model by controlling the FE model with intervertebral kinematics and muscle forces predicted by the musculoskeletal model. Annulus fibrosis stress, intradiscal pressure (IDP), and facet contact forces were calculated using the FE model. Differences in whole-body kinematics, muscle forces, and tissue-level loads were found between participant groups. Notably, participants with TTA had greater axial rotation toward their intact limb ([Formula: see text]), greater abdominal muscle activity ([Formula: see text]), and greater overall tissue loading throughout sit-to-stand ([Formula: see text]) compared to able-bodied participants. Both normalized (to upright standing) and absolute estimates of L4-L5 IDP were close to in vivo values reported in the literature. The multiscale model can be used to estimate the distribution of loads within different lumbar spine tissue structures and can be adapted for use with different activities, populations, and spinal geometries.

Influence of different postures under vertical impact load on thoracolumbar burst fracture

Clinical studies have extensively shown that burst fractures can cause severe and long-term neurological deficits. However, the mechanism of burst fracture is not clear, and the influence of different spinal postures on burst fracture is still unidentified. The study aimed at investigating the influence of different postures under vertical impact load on thoracolumbar burst fracture. A detailed nonlinear finite element model of T12-L2 segment was developed to investigate these problems. In this work, a rigid ball was used to vertically impact the finite element spinal segment, which emulated the process of burst fracture as in experimental condition. During the process, amounting to 9 different postures (normal, flexion, extension, right/left lateral bending of 8°, right/left axial rotation of 4° and 8°) were studied. Totally five failure modes were observed. Six different parameters, including vertebral height, vertebral bulge, interpedicular widening, vertebral kyphotic angle, posterior vertebral body angle, and joint facet contact force, were observed to evaluate the corresponding severity of burst fracture. Burst fracture in extension was the severest, and the loss of vertebral height in flexion was the most. The different postures in these simulations changed the morphology of intervertebral disc and facet joints force, resulting in different types of fracture.

Design of Patient Specific Spinal Implant (Pedicle Screw Fixation) using FE Analysis and Soft Computing Techniques

BACKGROUND

This work uses genetic algorithm (GA) for optimum design of patient specific spinal implants (pedicle screw) with varying implant diameter and bone condition. The optimum pedicle screw fixation in terms of implant diameter is on the basis of minimum strain difference from intact (natural) to implantation at peri-prosthetic bone for the considered six different peri-implant positions.

METHODS

This design problem is expressed as an optimization problem using the desirability function, where the data generated by finite element analysis is converted into an artificial neural network (ANN) model. The finite element model is generated from CT scan data. Thereafter all the ANN predictions of the microstrain in six positions are converted to unitless desirability value varying between 0 and 1, which is then combined to form the composite desirability. Maximization of the composite desirability is done using GA where composite desirability should be made to go up as close as possible to 1. If the composite desirability is 1, then all 'strain difference values in 6 positions' are 0.

RESULTS

The optimum solutions obtained can easily be used for making patient-specific spinal implants.

Effect of the intervertebral disc on vertebral bone strength prediction: a finite-element study

BACKGROUND CONTEXT

Osteoporotic vertebral fractures (OVFs) are a prevalent skeletal condition in the elderly but the mechanism behind these fractures remain unclear due to the complex biomechanical interplay between spinal segments such as the vertebra and intervertebral discs (IVDs).

PURPOSE

To investigate the biomechanical influence of IVDs by (1) comparing finite element (FE)-predicted failure load with experimentally measured failure load of functional spinal units (FSUs) and (2) comparing this correlation with those of FE-predicted failure load and bone mineral density (BMD) of the single central vertebra with experimentally measured failure load.

STUDY DESIGN

A computational biomechanical analysis.

PATIENT SAMPLE

Ten thoracic FSUs consisting of a central vertebra, the adjacent IVDs, and the upper and lower halves of the adjacent vertebrae were harvested from formalin-fixed human donors (4 males, 6 females; mean age of 82±9 years).

OUTCOME MEASURES

The outcome measures included the prediction of vertebral strength and determination of BMD in FSUs and the single central vertebra and the correlation of both measures with experimentally measured vertebral strength of the FSUs.

METHODS

The FSUs underwent clinical multidetector computed tomography (MDCT) (spatial resolution: 250×250×600 μm³). BMD was determined for the FSUs from the MDCT images of the central vertebrae. FE-predicted failure load was calculated in the single central vertebra of the FSUs alone and the entire FSUs. Experimentally measured failure load of the FSUs was determined in a uniaxial biomechanical test.

RESULTS

BMD of the central vertebrae correlated significantly with experimentally measured failure load (R²=0.66, p<.02), whereas FE-predicted failure load of the central vertebra showed no significant correlation with experimentally measured failure load (p=.07). However, FE-predicted failure load of FSUs best predicted experimentally measured failure load of FSUs (R²=0.93, p<.0001).

CONCLUSIONS

This study demonstrated that routine clinical MDCT images can be an accurate and feasible tool for prediction of OVFs using patient-specific FE analysis of FSU models.

CLINICAL SIGNIFICANCE

Improved management of OVFs is essential amidst current clinical challenges. Implementation of a vertebral strength assessment tool could result in more accurate prediction of osteoporotic fracture risk and aid clinicians with better targeted early treatment strategies.

Numerical simulation of lateral and transforaminal lumbar interbody fusion, two minimally invasive surgical approaches

The present study aims to compare spinal stability after two different minimally invasive techniques, the lateral lumbar interbody fusion (LLIF) and the transforaminal lumbar interbody fusion (TLIF) approaches. Two nonlinear three-dimensional finite element (FE) models of the L4-L5 functional spinal unit (FSU) were subjected to the loads that usually act on the lumbar spine. Findings show that the LLIF approach yields better results for torsion load case, due to the larger surface area of the implant. For extension, flexion and lateral bending loads, the TLIF approach presents smaller displacements probably due to the anterior placement of the cage and to the smaller damaged area of the annulus fibrosus.

Minimally invasive transforaminal and anterior lumbar interbody fusion surgery at level L5-S1

This paper presents a finite element analysis to investigate the biomechanical changes caused by transforaminal (TLIF) and anterior lumbar interbody fusion (ALIF) at the L5-S1 level, applying two different implants: T_PAL (TLIF) and SynFix (ALIF). The main objective is to determine which one is more stable for patients. Numerical simulations of segmental motion show that, in the early postoperative phase, displacements and rotation angles obtained in ALIF are greater than the corresponding ones obtained in TLIF, as well as the principal stress values on the ligaments. So, TLIF performed with T_PAL is more stable than ALIF, especially during the recovery phase.

Comparison of Biomechanical Performance of Five Different Treatment Approaches for Fixing Posterior Pelvic Ring Injury

Background

A large number of pelvic injuries are seriously unstable, with mortality rates reaching 19%. Approximately 60% of pelvic injuries are related to the posterior pelvic ring. However, the selection of a fixation method for a posterior pelvic ring injury remains a challenging problem for orthopedic surgeons. The aim of the present study is to investigate the biomechanical performance of five different fixation approaches for posterior pelvic ring injury and thus provide guidance on the choice of treatment approach in a clinical setting.

Methods

A finite element (FE) model, including the L3-L5 lumbar vertebrae, sacrum, and full pelvis, was created from CT images of a healthy adult. Tile B and Tile C types of pelvic fractures were created in the model. Five different fixation methods for fixing the posterior ring injury (PRI) were simulated: TA1 (conservative treatment), TA2 (S1 screw fixation), TA3 (S1 + S2 screw fixation), TA4 (plate fixation), and TA5 (modified triangular osteosynthesis). Based on the fixation status (fixed or nonfixed) of the anterior ring and the fixation method for PRI, 20 different FE models were created. An upright standing loading scenario was simulated, and the resultant displacements at the sacroiliac joint were compared between different models.

Results

When TA5 was applied, the resultant displacements at the sacroiliac joint were the smallest (1.5 mm, 1.6 mm, 1.6 mm, and 1.7 mm) for all the injury cases. The displacements induced by TA3 and TA2 were similar to those induced by TA5. TA4 led to larger displacements at the sacroiliac joint (2.3 mm, 2.4 mm, 4.8 mm, and 4.9 mm), and TA1 was the worst case (3.1 mm, 3.2 mm, 6.3 mm, and 6.5 mm).

Conclusions

The best internal fixation method for PRI is the triangular osteosynthesis approach (TA5), followed by S1 + S2 screw fixation (TA3), S1 screw fixation (TA2), and plate fixation (TA4).

A patient specific computational biomechanical model for the entire lumbosacral spinal unit with imposed spondylolysis

BACKGROUND

A biomechanical model of the lumbosacral spinal unit between L1-S1 was developed to investigate the behavior of normal and select pathological states. Our aims were to generate predictive structural models for mechanical deformation including critical stresses in the spine components and to investigate the probability of subsequent lumbar spine fractures in the presence of unilateral spondylolysis.

METHODS

A non-linear three-dimensional finite element pathology-free model of the L1-S1 lumbosacral unit was generated using patient-specific computerized tomography scans and calibrated by comparing it to experimental data of a range of motion modes consisting of flexion, extension, left and right lateral bending, and left and right axial rotation. Unilateral and bilateral pars defects were created on the isthmus of L5 to simulate spondylolysis.

FINDINGS

Results showed that under flexion, left lateral bending and right axial rotation, stresses were higher on the contralateral L5 pars-interarticularis, whereas, no significant changes occurred on the left-right isthmus of the L2-L4 and S1. Significant changes in the range of motion compared to the pathology-free model were observed in bilateral spondylolysis not only adjacent to the pars defect area but also in other lumbar spine levels.

INTERPRETATION

The proposed pathology-free lumbosacral unit model showed good correlation with experimental tests for all loading cases. In unilateral spondylolysis, a subsequent pars defect was observed within the same vertebra. The overall modeling approach can be used to study different pathological states.

Evaluation of the capability of the simulated dual energy X-ray absorptiometry-based two-dimensional finite element models for predicting vertebral failure loads

Prediction of the vertebral failure load is of great importance for the prevention and early treatment of bone fracture. However, an efficient and effective method for accurately predicting the failure load of vertebral bones is still lacking. The aim of the present study was to evaluate the capability of the simulated dual energy X-ray absorptiometry (DXA)-based finite element (FE) model for predicting vertebral failure loads. Thirteen dissected spinal segments (T11/T12/L1) were scanned using a HR-pQCT scanner and then were mechanically tested until failure. The subject-specific three-dimensional (3D) and two-dimensional (2D) FE models of T12 were generated from the HR-pQCT scanner and the simulated DXA images, respectively. Additionally, the areal bone mineral density (aBMD) and areal bone mineral content (aBMC) of T12 were calculated. The failure loads predicted by the simulated DXA-based 2D FE models were more moderately correlated with the experimental failure loads (R² = 0.66) than the aBMC (R² = 0.61) and aBMD (R² = 0.56). The 2D FE models were slightly outperformed by the HR-pQCT-based 3D FE models (R² = 0.71). The present study demonstrated that the simulated DXA-based 2D FE model has better capability for predicting the vertebral failure loads than the densitometric measurements but is outperformed by the 3D FE model. The 2D FE model is more suitable for clinical use due to the low radiation dose and low cost, but it remains to be validated by further in vitro and in vivo studies.

Is there any advantage of using stand-alone cages? A numerical approach

Background

Segment fusion using interbody cages supplemented with pedicle screw fixation is the most common surgery for the treatment of low back pain. However, there is still much controversy regarding the use of cages in a stand-alone fashion. The goal of this work is to numerically compare the influence that each surgery has on lumbar biomechanics.

Methods

A non-linear FE model of the whole lumbar spine was developed to compare between two types of cages (OLYS and NEOLIF) with and without supplementary fixation. The motion of the whole spine was analysed and the biomechanical environment of the adjacent segments to the operated one was studied. Moreover, the risk of subsidence of the cages was qualitatively evaluated.

Results

A great ROM reduction occurred when supplementary fixation was used. This stiffening increased the stresses at the adjacent levels. It might be hypothesised that the overloading of these segments could be related with the clinically observed adjacent disc degeneration. Meanwhile, the stand-alone cages allowed for a wider movement, and therefore, the influence of the surgery on adjacent discs was much lower. Regarding the risk of subsidence, the contact pressure magnitude was similar for both intervertebral cage designs and near the value of the maximum tolerable pressure of the endplates.

Conclusions

A minimally invasive posterior insertion of an intervertebral cage (OLYS or NEOLIF) was compared using a stand-alone design or adding supplementary fixation. The outcomes of these two techniques were compared, and although stand-alone cage may diminish the risk of disease progression to the adjacent discs, the spinal movement in this case could compromise the vertebral fusion and might present a higher risk of cage subsidence.

Electronic supplementary material

The online version of this article (10.1186/s12938-019-0684-8) contains supplementary material, which is available to authorized users.

The Mechanical Role of the Radial Fiber Network Within the Annulus Fibrosus of the Lumbar Intervertebral Disc: A Finite Elements Study

The annulus fibrosus (AF) of the intervertebral disc (IVD) consists of a set of concentric layers composed of a primary circumferential collagen fibers arranged in an alternating oblique orientation. Moreover, there exists an additional secondary set of radial translamellar collagen fibers which connects the concentric layers, creating an interconnected fiber network. The aim of this study was to investigate the mechanical role of the radial fiber network. Toward that goal, a three-dimensional (3D) finite element model of the L3–L4 spinal segment was generated and calibrated to axial compression and pure moment loading. The AF model explicitly recognizes the two heterogeneous networks of fibers. The presence of radial fibers demonstrated a pronounced effect on the local disc responses under lateral bending, flexion, and extension modes. In these modes, the radial fibers were in a tensile state in the disc region that subjected to compression. In addition, the circumferential fibers, on the opposite side of the IVD, were also under tension. The local stress in the matrix was decreased in up to 9% in the radial fibers presence. This implies an active fiber network acting collectively to reduce the stresses and strains in the AF lamellae. Moreover, a reduction of 26.6% in the matrix sideways expansion was seen in the presence of the radial fibers near the neutral bending axis of the disc. The proposed biomechanical model provided a new insight into the mechanical role of the radial collagen fibers in the AF structure. This model can assist in the design of future IVD substitutes.

Stand-alone lumbar cage subsidence: A biomechanical sensitivity study of cage design and placement

BACKGROUND AND OBJECTIVE

Spinal degeneration and instability are commonly treated with interbody fusion cages either alone or supplemented with posterior instrumentation with the aim to immobilise the segment and restore intervertebral height. The purpose of this work is to establish a tool which may help to understand the effects of intervertebral cage design and placement on the biomechanical response of a patient-specific model to help reducing post-surgical complications such as subsidence and segment instability.

METHODS

A 3D lumbar functional spinal unit (FSU) finite element model was created and a parametric model of an interbody cage was designed and introduced in the FSU. A Drucker-Prager Cap plasticity formulation was used to predict plastic strains and bone failure in the vertebrae. The effect of varying cage size, cross-sectional area, apparent stiffness and positioning was evaluated under 500 N preload followed by 7.5 Nm multidirectional rotation and the results were compared with the intact model.

RESULTS

The most influential cage parameters on the FSU were size, curvature congruence with the endplates and cage placement. Segmental stiffness was higher when increasing the cross-sectional cage area in all loading directions and when the cage was anteriorly placed in all directions but extension. In general, the facet joint forces were reduced by increasing segmental stiffness. However, these forces were higher than in the intact model in most of the cases due to the displacement of the instantaneous centre of rotation. The highest plastic deformations took place at the caudal vertebra under flexion and increased for cages with greater stiffness. Thus, wider cages and a more anteriorly placement would increase the volume of failed bone and, therefore, the risk of subsidence.

CONCLUSIONS

Cage geometry plays a crucial role in the success of lumbar surgery. General considerations such as larger cages may be applied as a guideline, but parameters such as curvature or cage placement should be determined for each specific patient. This model provides a proof-of-concept of a tool for the preoperative evaluation of lumbar surgical outcomes.

Effect of two-level pedicle-screw fixation with different rod materials on lumbar spine: A finite element study

BACKGROUND

Pedicle-screw-rod fixation system is very popular surgical remedy for degenerative disc disease. It is important to observe load vs. spinal motion characteristic for better understanding of clinical problems and treatment of spinal instability associated with low-back pain.

OBJECTIVE

The objective of this study is to understand the effect [range of motion (ROM) and intervertebral foramen height] of pedicle-screw fixation with three rod materials on lumbar spine under three physiological loading conditions.

METHOD

A three-dimensional finite element (FE) model of lumbar to sacrum (L1-S) vertebrae with pedicle-screw-rod fixation at L3-L5 level is developed. Three rod materials [titanium alloy (Ti6Al4V), ultra-high molecular weight poly ethylene (UHMWPE) and poly-ether-ether-ketone (PEEK)] are used for two-level fixation and the FE models are simulated for axial rotation, lateral bending and flexion-extension under ±10 Nm moment and 500 N compressive load and compared with the intact (natural) model.

RESULT & DISCUSSION

For axial rotation, lateral bending and flexion, ROM increased 2.8, 4.5 and 1.83 times respectively for UHMWPE, and 3.7, 7.2 and 2.15 times respectively for PEEK in comparison to Ti6Al4V. As ROM is 49, 29 and 31% of the intact model during axial rotation, lateral bending and flexion respectively, PEEK rod produced better motion flexibility than Ti6Al4V and UHMWPE rod. Foramen height increased insignificantly by 2.21% for the PEEK rod with respect to the intact spine during flexion. For the PEEK rod, maximum stress of 40 MPa during axial rotation is much below the yield stress of 98 MPa.

CONCLUSION

Ti6Al4V pedicle-screw-rod fixation system highly restricted the ROM of the spine, which is improved by using UHMWPE and PEEK, having lower stiffness. The foramen height did not vary significantly for any implant materials. In terms of ROM and maximum stress, PEEK rod may be considered for a better implant design to get better ROM and thus mobility.

Human body modeling method to simulate the biodynamic characteristics of spine in vivo with different sitting postures

The aim of this study is to model the computational model of seated whole human body including skeleton, muscle, viscera, ligament, intervertebral disc, and skin to predict effect of the factors (sitting postures, muscle and skin, buttocks, viscera, arms, gravity, and boundary conditions) on the biodynamic characteristics of spine. Two finite element models of seated whole body and a large number of finite element models of different ligamentous motion segments were developed and validated. Static, modal, and transient dynamic analyses were performed. The predicted vertical resonant frequency of seated body model was in the range of vertical natural frequency of 4 to 7 Hz. Muscle, buttocks, viscera, and the boundary conditions of buttocks have influence on the vertical resonant frequency of spine. Muscle played a very important role in biodynamic response of spine. Compared with the vertical posture, the posture of lean forward or backward led to an increase in stress on anterior or lateral posterior of lumbar intervertebral discs. This indicated that keeping correct posture could reduce the injury of vibration on lumbar intervertebral disc under whole-body vibration. The driving posture not only reduced the load of spine but also increased the resonant frequency of spine.

The Effects of Physiological Biomechanical Loading on Intradiscal Pressure and Annulus Stress in Lumbar Spine: A Finite Element Analysis

The present study was conducted to examine the effects of body weight on intradiscal pressure (IDP) and annulus stress of intervertebral discs at lumbar spine. Three-dimensional finite element model of osseoligamentous lumbar spine was developed subjected to follower load of 500 N, 800 N, and 1200 N which represent the loads for individuals who are normal and overweight with the pure moments at 7.5 Nm in flexion and extension motions. It was observed that the maximum IDP was 1.26 MPa at L1-L2 vertebral segment. However, the highest increment of IDP was found at L4-L5 segment where the IDP was increased to 30% in flexion and it was more severe at extension motion reaching to 80%. Furthermore, the maximum annulus stress also occurred at the L1-L2 segment with 3.9 MPa in extension motion. However, the highest increment was also found at L4-L5 where the annulus stress increased to 17% in extension motion. Based on these results, the increase of physiological loading could be an important factor to the increment of intradiscal pressure and annulus fibrosis stress at all intervertebral discs at the lumbar spine which may lead to early intervertebral disc damage.

Mechanical assessment of the effects of metastatic lytic defect on the structural response of human thoracolumbar spine

To investigate the effects of a clinical lytic defect on the structural response of human thoracolumbar functional spinal unit. A novel CT-compatible mechanical test system was used to image the deformation of a T12-L1 motion segment and measure the change in strain response under compressive loads ranging from 50 to 750 N. A lytic lesion (LM) with cortex involvement (33% by volume) was introduced to the upper vertebral body and the CT experiments were repeated. Finite element models, established from the CT volumes, were used to investigate the defect's effects on the structural response and the state of principal and shear stresses within the affected and adjacent vertebrae. The lytic lesion resulted in severe loss of the vertebral structural competence, resulting in significant, non-linear, and asymmetric increase in the experimentally measured strains and computed stresses within both vertebrae (p < 0.01). At the cortex, the tensile strains were significantly increased, while compressive strains significantly decreased, (p < 0.05). Both the vertebral bone and cortex regions adjacent to the defect showed significant increase in computed compressive, tensile, and shear stresses (p < 0.01). Changes in stress and strain distribution within the affected and adjacent vertebral bone and the experimentally observed bulging and buckling of the vertebral cortices suggested that initiation of catastrophic vertebral failure may occur under load magnitudes encountered in daily living. Although the effect of LM on the global deformation of the spine was well-predicted, our results show that FE predictions of local strain changes must be carefully assessed for clinical relevance. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1808-1819, 2016.

Influence of different fusion techniques in lumbar spine over the adjacent segments: A 3D finite element study

The most conventional technique to treat the intervertebral disc degeneration consists on fusing the affected segment with a posterior screw fixation and sometimes with the insertion of a cage in the intersomatic space. However, this kind of surgeries had controversial results in the adjacent discs. The aim of this work was to prove the stabilization of the spine and the decompression of the disc and to analyze the influence over the adjacent segments. With this purpose, four different models were built and simulated under different loading conditions. The stabilization of the spine was ensured by the screw fixation which reduced dramatically the relative motion in the affected segment. On the other hand, the pore pressure showed a high fall in the operated models proving the decompression of the neural structures. In the adjacent segments, the ROM increased up to 50% in the upper disc and 70% in the lower one. The pore pressure and principal stresses also increased after both surgeries. The observed results suggested that the fusion procedure could trigger a cascade degeneration effect over the adjacent discs, while it is also seen that cage insertion helps to maintain disc height in a better way than screw fixation only.

On the Use of Biaxial Properties in Modeling Annulus as a Holzapfel-Gasser-Ogden Material

Besides the biology, stresses and strains within the tissue greatly influence the location of damage initiation and mode of failure in an intervertebral disk. Finite element models of a functional spinal unit (FSU) that incorporate reasonably accurate geometry and appropriate material properties are suitable to investigate such issues. Different material models and techniques have been used to model the anisotropic annulus fibrosus, but the abilities of these models to predict damage initiation in the annulus and to explain clinically observed phenomena are unclear. In this study, a hyperelastic anisotropic material model for the annulus with two different sets of material constants, experimentally determined using uniaxial and biaxial loading conditions, were incorporated in a 3D finite element model of a ligamentous FSU. The purpose of the study was to highlight the biomechanical differences (e.g., intradiscal pressure, motion, forces, stresses, strains, etc.) due to the dissimilarity between the two sets of material properties (uniaxial and biaxial). Based on the analyses, the biaxial constants simulations resulted in better agreements with the in vitro and in vivo data, and thus are more suitable for future damage analysis and failure prediction of the annulus under complex multiaxial loading conditions.

Finite element analysis predicts experimental failure patterns in vertebral bodies loaded via intervertebral discs up to large deformation

Vertebral compression fractures are becoming increasingly common. Patient-specific nonlinear finite element (FE) models have shown promise in predicting yield strength and damage pattern but have not been experimentally validated for clinically relevant vertebral fractures, which involve loading through intervertebral discs with varying degrees of degeneration up to large compressive strains. Therefore, stepwise axial compression was applied in vitro on segments and performed in silico on their FE equivalents using a nonlocal damage-plastic model including densification at large compression for bone and a time-independent hyperelastic model for the disc. The ability of the nonlinear FE models to predict the failure pattern in large compression was evaluated for three boundary conditions: healthy and degenerated intervertebral discs and embedded endplates. Bone compaction and fracture patterns were predicted using the local volume change as an indicator and the best correspondence was obtained for the healthy intervertebral discs. These preliminary results show that nonlinear finite element models enable prediction of bone localisation and compaction. To the best of our knowledge, this is the first study to predict the collapse of osteoporotic vertebral bodies up to large compression using realistic loading via the intervertebral discs.

Finite element based nonlinear normalization of human lumbar intervertebral disc stiffness to account for its morphology

Disc degeneration, usually associated with low back pain and changes of intervertebral stiffness, represents a major health issue. As the intervertebral disc (IVD) morphology influences its stiffness, the link between mechanical properties and degenerative grade is partially lost without an efficient normalization of the stiffness with respect to the morphology. Moreover, although the behavior of soft tissues is highly nonlinear, only linear normalization protocols have been defined so far for the disc stiffness. Thus, the aim of this work is to propose a nonlinear normalization based on finite elements (FE) simulations and evaluate its impact on the stiffness of human anatomical specimens of lumbar IVD. First, a parameter study involving simulations of biomechanical tests (compression, flexion/extension, bilateral torsion and bending) on 20 FE models of IVDs with various dimensions was carried out to evaluate the effect of the disc's geometry on its compliance and establish stiffness/morphology relations necessary to the nonlinear normalization. The computed stiffness was then normalized by height (H), cross-sectional area (CSA), polar moment of inertia (J) or moments of inertia (Ixx, Iyy) to quantify the effect of both linear and nonlinear normalizations. In the second part of the study, T1-weighted MRI images were acquired to determine H, CSA, J, Ixx and Iyy of 14 human lumbar IVDs. Based on the measured morphology and pre-established relation with stiffness, linear and nonlinear normalization routines were then applied to the compliance of the specimens for each quasi-static biomechanical test. The variability of the stiffness prior to and after normalization was assessed via coefficient of variation (CV). The FE study confirmed that larger and thinner IVDs were stiffer while the normalization strongly attenuated the effect of the disc geometry on its stiffness. Yet, notwithstanding the results of the FE study, the experimental stiffness showed consistently higher CV after normalization. Assuming that geometry and material properties affect the mechanical response, they can also compensate for one another. Therefore, the larger CV after normalization can be interpreted as a strong variability of the material properties, previously hidden by the geometry's own influence. In conclusion, a new normalization protocol for the intervertebral disc stiffness in compression, flexion, extension, bilateral torsion and bending was proposed, with the possible use of MRI and FE to acquire the discs' anatomy and determine the nonlinear relations between stiffness and morphology. Such protocol may be useful to relate the disc's mechanical properties to its degree of degeneration.

Biomechanical analysis of two-step traction therapy in the lumbar spine

Traction therapy is one of the most common conservative treatments for low back pain. However, the effects of traction therapy on lumbar spine biomechanics are not well known. We investigated biomechanical effects of two-step traction therapy, which consists of global axial traction and local decompression, on the lumbar spine using a validated three-dimensional finite element model of the lumbar spine. One-third of body weight was applied on the center of the L1 vertebra toward the superior direction for the first axial traction. Anterior translation of the L4 vertebra was considered as the second local decompression. The lordosis angle between the superior planes of the L1 vertebra and sacrum was 44.6° at baseline, 35.2° with global axial traction, and 46.4° with local decompression. The fibers of annulus fibrosus in the posterior region, and intertransverse and posterior longitudinal ligaments experienced stress primarily during global axial traction, these stresses decreased during local decompression. A combination of global axial traction and local decompression would be helpful for reducing tensile stress on the fibers of the annulus fibrosus and ligaments, and intradiscal pressure in traction therapy. This study could be used to develop a safer and more effective type of traction therapy.

Comparison of eight published static finite element models of the intact lumbar spine: predictive power of models improves when combined together

Finite element (FE) model studies have made important contributions to our understanding of functional biomechanics of the lumbar spine. However, if a model is used to answer clinical and biomechanical questions over a certain population, their inherently large inter-subject variability has to be considered. Current FE model studies, however, generally account only for a single distinct spinal geometry with one set of material properties. This raises questions concerning their predictive power, their range of results and on their agreement with in vitro and in vivo values. Eight well-established FE models of the lumbar spine (L1-5) of different research centers around the globe were subjected to pure and combined loading modes and compared to in vitro and in vivo measurements for intervertebral rotations, disc pressures and facet joint forces. Under pure moment loading, the predicted L1-5 rotations of almost all models fell within the reported in vitro ranges, and their median values differed on average by only 2° for flexion-extension, 1° for lateral bending and 5° for axial rotation. Predicted median facet joint forces and disc pressures were also in good agreement with published median in vitro values. However, the ranges of predictions were larger and exceeded those reported in vitro, especially for the facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with measured in vivo values. In light of high inter-subject variability, the generalization of results of a single model to a population remains a concern. This study demonstrated that the pooled median of individual model results, similar to a probabilistic approach, can be used as an improved predictive tool in order to estimate the response of the lumbar spine.