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

Finite element models of intervertebral disc: recent advances and prospects

OBJECTIVES

The incidence rate of intervertebral disc degeneration (IVDD) is increasing year by year, which brings great harm to our health. The change of biomechanical factors is an important reason for IVDD. Therefore, more and more studies use finite element (FE) models to analyze the biomechanics of spine.

METHODS

In this review, literatures which reported the FE model of intervertebral disc (IVD) were reviewed. We summarized the types and constructional methods of the FE models and analyzed the applications of some representative FE models.

RESULTS

The most widely used model was the nonlinear model which considers the behavior of porous elastic materials. As more advanced methods, More and more models which involve penetration parameters were used to simulate the biological behavior and biomechanical properties of IVD.

CONCLUSIONS

Personalized modeling should be carried out in order to better provide accurate basis for the diagnosis and treatment of the disease. In addition, microstructure, cell behavior and complex load should be considered in the process of model construction to build a more realistic model.

Effect of vehicular vibrations on L-4 lumbar vertebrae - A finite element study

Lower Back Pain (LBP) is a global health issue, with increasing prevalence, partly attributed to vehicular vibrations experienced by motorcyclists. The L4 lumbar vertebra is responsible for greater mobility and flexibility of the body, but also is the most crucial body element affected by vehicular vibrations. Anthropometric properties, types of speed humps, and vehicle types are the critical variables that impact bone health during riding, need to be studied. To understand the potential zones of injury, computational simulation can be performed under the influence of vehicle vibrations while crossing different types of speed humps at varying speeds. In the present study, finite element method (FEM) is used to evaluate stress and deformation in the bone. The L4 cortical bone is modelled by considering the CT-Scan data and assumed to be homogeneous and isotropic material. Vibration data is collected using two vehicle types (Type I and Type II) on four different humps (Trapezoidal, Bitumen Semi-circular, Rubber Semi-circular, and Rumble strip). The bone's dynamic behavior is studied using FEM simulation, which involved static structural, modal and transient dynamic analyses. The findings from static analysis indicate that the most concentrated stress is located in the lower pedicle region and is an expected commonplace for injuries because of vibrations. In transient dynamic analysis, Type I vehicle showed a 25 % higher stress than Type II.

In Vitro Biomechanical Experiment on the Effect of Unilateral Partial Facetectomy Performed by Percutaneous Endoscopy on the Stability of Lumbar Spine

OBJECTIVES

This study's purpose is to investigate the lumbar biomechanical effects of unilateral partial facetectomy (UPF) of different facet joint (FJ) portions under percutaneous endoscopy.

METHODS

Forty fresh calf spine models were used to simulate UPF under a physiological load performed through three commonly used needle insertion points (IPs): (1) The apex of the superior FJ (as the first IP); (2) The midpoint of the ventral side of the superior FJ (as the second IP); (3) The lowest point of the ventral side of the superior FJ (as the third IP). The range of motion (ROM) and the L4/5 intradiscal maximum pressure (IMP) were measured and analyzed under a physiological load in all models during flexion, extension, left-right lateral flexion, and left-right axial rotation.

RESULTS

When UPF was performed through the second IP, the ROM of the lumbar spine and the L4/5 IMP in the calf spine models were not statistically different from the intact calf spine model.

CONCLUSIONS

UPF through the second IP resulted in a minimal impact on the biomechanics of the lumbar spine. Thus, it might be considered the most appropriate IP for UPF.

Lumbar Total Disc Replacement Leads to Increased Subsequent Facet Injections Compared to Anterolateral Lumbar Interbody Fusions

Study DesignRetrospective Matched Cohort.ObjectiveDespite known consequences to the facet joints following lumbar total disc replacement (TDR), there is limited data on facet injection usage for persistent postoperative pain. This study uses real-world data to compare the usage of therapeutic lumbar facet injections as a measure of symptomatic facet arthrosis following single-level, stand-alone TDR vs anterolateral lumbar interbody fusion (ALIF/LLIF).MethodsThe PearlDiver database was queried for patients (2010-2021) with lumbar degenerative disc disease who received either a single-level, stand-alone TDR or ALIF/LLIF. All patients were followed for ≥2 years and excluded if they had a history of facet injections or spinal trauma, fracture, infection, or neoplasm. The two cohorts were matched 1:1 based on age, sex, insurance, year of operation, and medical comorbidities. The primary outcome was the use of therapeutic lumbar facet injections at 1-, 2-, and 5-year follow-up. Secondary outcomes included subsequent lumbar surgeries and surgical complications.ResultsAfter 1:1 matching, each cohort had 1203 patients. Lumbar facet injections occurred significantly more frequently in the TDR group at 1-year (6.07% vs 1.66%, P < .0001), 2-year (8.40% vs 3.74%%, P < .0001), and 5-year (11.47% vs 6.40%, P < .0001) follow-up. 5-year injection-free probability curves demonstrated an 87.1% injection-free rate for TDR vs 92.9% for ALIF/LLIF. There was no clinical difference in the incidence of subsequent lumbar surgeries or complications.ConclusionCompared with ALIF/LLIF, patients who underwent TDR received significantly more facet injections, suggesting a greater progression of symptomatic facet arthrosis. TDR was not protective against reoperations compared to ALIF/LLIF.

Sensitivity of Lumbar Total Joint Replacement Contact Stresses Under Misalignment Conditions-Finite Element Analysis of a Spine Wear Simulator

A novel total joint replacement (TJR) that treats lumbar spine degeneration was previously assessed under Mode I and Mode IV conditions. In this study, we relied on these previous wear tests to establish a relationship between finite element model (FEM)-based bearing stresses and in vitro wear penetration maps. Our modeling effort addressed the following question of interest: Under reasonably worst-case misaligned conditions, do the lumbar total joint replacement (L-TJR) polyethylene stresses and strains remain below values associated with Mode IV impingement wear tests? The FEM was first formally verified and validated using the risk-informed credibility assessment framework established by ASME V&V 40 and FDA guidance. Then, based on criteria for unreasonable misuse outlined in the surgical technique guide, a parametric analysis of reasonably worst-case misalignment using the validated L-TJR FEM was performed. Reasonable misalignment was created by altering device positioning from the baseline condition in three scenarios: Axial Plane Convergence (20-40°), Axial Plane A-P Offset (0-4 mm), and Coronal Plane Tilt (±20°). We found that, for the scenarios considered, the contact pressures, von Mises stresses, and effective strains of the L-TJR-bearing surfaces remained consistent with Mode I (clean) conditions. Specifically, the mechanical response values fell below those determined under Mode IV (worst-case) boundary conditions, which provided the upper-bound benchmarks for the study (peak contact pressure 83.3 MPa, peak von Mises stress 32.2 MPa, and peak effective strain 42%). The L-TJR was judged to be insensitive to axial and coronal misalignment under the in vitro boundary conditions imposed by the study.

Enhancing Biomechanical Simulations Based on a Posteriori Error Estimates: The Potential of Dual-Weighted Residual-Driven Adaptive Mesh Refinement

The finite-element method (FEM) is a well-established procedure for computing approximate solutions to deterministic engineering problems described by partial differential equations. FEM produces discrete approximations of the solution with a discretisation error that can be quantified with a posteriori error estimates. The practical relevance of error estimates for biomechanics problems, especially for soft tissue where the response is governed by large strains, is rarely addressed. In this contribution, we propose an implementation of a posteriori error estimates targeting a user-defined quantity of interest, using the dual-weighted residual (DWR) technique tailored to biomechanics. The proposed method considers a general setting that encompasses three-dimensional geometries and model nonlinearities, which appear in hyperelastic soft tissues. We take advantage of the automatic differentiation capabilities embedded in modern finite-element software, which allows the error estimates to be computed generically for a large class of models and constitutive laws. First, we validate our methodology using experimental measurements from silicone samples and then illustrate its applicability for patient-specific computations of pressure ulcers on a human heel.

Combined effect of artificial cervical disc replacement and facet tropism on the index-level facet joints: a finite element study

BACKGROUND

Artificial Cervical Disc Replacement (ACDR) is an effective treatment for cervical degenerative disc diseases. However, clinical information regarding the facet joint alterations after ACDR was limited. Facet tropism is common in the sub-axial cervical spine. Our previous research indicated that facet tropism could lead to increased pressure on the cervical facet joints. This study aimed to assess the impact of facet tropism on the facet contact force and facet capsule stress after ACDR.

METHODS

A C2-T1 cervical finite element model was constructed from computed tomography (CT) scans of a 28-year-old male volunteer. Symmetrical, moderate asymmetrical (7 degrees tropism), and severe asymmetrical (14 degrees tropism) models were created at the C5/C6 level by altering the facet orientation at the C5-C6 level. The C5/C6 ACDR was simulated in the intact, moderate asymmetrical and severe asymmetrical models. A 75-N follower load with 1.0-Nm moments was applied to the top of C2 vertebra in the models to simulate flexion, extension, lateral bending, and axial rotation with the T1 vertebra fixed. The range of motions (ROMs) under all moments, facet contact forces (FCFs) and facet capsule strains were tested.

RESULTS

In the asymmetrical model, the right FCFs considerably increased under flexion, extension, right bending, left rotation, especially under right bending the right sided FCF of the severe asymmetrical model was about 5.44 times of the neutral position, and 3.14 times of the symmetrical model. and concentrated on the cephalad part of the facets. The facet capsule stresses on both sides remarkably increased under extension, lateral bending and right rotation. In the moderate and severe asymmetrical models, the capsule strain was greater on both sides of each position than in the symmetric model.

CONCLUSIONS

The face tropism increased facet contact force and facet capsule strain after ACDR, especially under extension, lateral bending, and rotation, and also could result in abnormal stress distribution on the facet joint surface and facet joint capsule. The results suggest that face tropism might be a risk factor for post-operative facet joint degeneration progression after ACDR. Facet tropism may be noteworthy when ACDR is considered as a surgical option.

Comparative finite element analysis of posterior short segment fixation constructs with or without intermediate screws in the fractured vertebrae for the treatment of type a thoracolumbar fracture

Six-screw short-segment posterior fixation for thoracolumbar fractures, which involves intermediate screws at the fractured vertebrae has been proposed to reduce the rates of kyphosis recurrence and implant failure. Yet, little is known about the mechanisms and biomechanical responses by which intermediate screws at the fracture vertebrae enhance fixation strength. The objective of this study was to investigate the biomechanical properties that are associated with the augmentation of intermediate screws in relation to the severity of type A thoracolumbar fracture using finite element analysis. Short-segment stabilization models with or without augmentation screws at fractured vertebrae were established based on finite element model of moderate compressive fractures, severe compressive fractures and burst fractures. The spinal stiffness, stresses at the implanted hardware, and axial displacement of the bony defect were measured and compared under mechanical loading conditions. All six-screw stabilization showed a decreased range of motion in extension, lateral bending, and axial rotation compared to the traditional four-screw fixation models. Burst thoracolumbar fracture benefited more from augmentation of intermediate screws at the fracture vertebrae. The stress of the rod in six-screw models increased while decreased that of pedicle screws. Our results suggested that patients with more unstable fractures might achieve greater benefits from augmentation of intermediate screws at the fracture vertebrae. Augmentation of intermediate screws at the fracture vertebrae is recommended for patients with higher wedge-shaped or burst fractures to reduce the risk of hardware failure and postoperative re-collapse of injured vertebrae.

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.

Validation and Estimation of Obesity-Induced Intervertebral Disc Degeneration through Subject-Specific Finite Element Modelling of Functional Spinal Units

(1) Background: Intervertebral disc degeneration has been linked to obesity; its potential mechanical effects on the intervertebral disc remain unknown. This study aimed to develop and validate a patient-specific model of L3-L4 vertebrae and then use the model to estimate the impact of increasing body weight on disc degeneration. (2) Methods: A three-dimensional model of the functional spinal unit of L3-L4 vertebrae and its components were developed and validated. Validation was achieved by comparing the range of motions (RoM) and intradiscal pressures with the previous literature. Subsequently, the validated model was loaded according to the body mass index and estimated stress, deformation, and RoM to assess disc degeneration. (3) Results: During validation, L3-L4 RoM and intradiscal pressures: flexion 5.17° and 1.04 MPa, extension 1.54° and 0.22 MPa, lateral bending 3.36° and 0.54 MPa, axial rotation 1.14° and 0.52 MPa, respectively. When investigating the impact of weight on disc degeneration, escalating from normal weight to obesity reveals an increased RoM, by 3.44% during flexion, 22.7% during extension, 29.71% during lateral bending, and 33.2% during axial rotation, respectively. Also, stress and disc deformation elevated with increasing weight across all RoM. (4) Conclusions: The predicted mechanical responses of the developed model closely matched the validation dataset. The validated model predicts disc degeneration under increased weight and could lay the foundation for future recommendations aimed at identifying predictors of lower back pain due to disc degeneration.

[Effect of fiber volume and material property and nucleus pulposus area on intervertebral disc mechanical behavior]

The material properties and volume proportion of the fibers as well as the cross-sectional area proportion of nucleus pulposus vary greatly in different studies. The effect of these factors on the mechanical behavior of intervertebral discs (IVDs) are uncertain. The IVDs finite element models with different parameters were created to investigate the pressure, height, rotation, stress, and strain of the IVDs under loads: pure compression, rotation after compression or axial moment after compression. The results showed that the material properties of fibers had great impact on the mechanical behavior of IVDs, especially on the rotation angle. When the fiber volume ratio was small, its changes had a significant impact on the rotation angle of the IVDs. The area proportions of nucleus pulposus had relatively little effect on the mechanical behavior of IVDs. The IVDs rotation should be observed when validating the model. By adjusting the elastic modulus or volume ratio of fibers within a reasonable range, a model that could simulate the mechanical behavior of normal IVDs could be obtained. It was reasonable to make the area proportion of nucleus pulposus within 25%-50% for the IVDs finite element model. This study provides guidance and reference for finite element modeling of the IVDs and the investigation of the IVDs degeneration mechanism.

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

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

Cadaveric biomechanical studies of ADDISC total lumbar disc prosthesis

BACKGROUND

Most total disc replacements provide excessive mobility and not reproduce spinal kinematics, inducing zygapophyseal joint arthritic changes and chronic back pain. In cadaveric lumbosacral spines, we studied if a new lumbar disc prosthesis kinematics mimics the intact intervertebral disc.

METHODS

In eight cold preserved cadaveric lumbosacral spines, we registered the movement ranges in flexion, extension, right and left lateral bending, and rotation in the intact status, post-discectomy, and after our prosthesis implantation, comparing them for each specimen.

FINDINGS

Comparing the intact lumbosacral spine with the L4-L5 prosthesis implanted specimens, we saw statistically significant differences in lateral bending and right rotation but not in the full range of rotation. Analyzing segments, we also noticed statistically significant differences at L4-L5 in flexion-extension and rotation. On the other hand, the L4-L5 discectomy, compared to the baseline spine condition, showed a statistically significant mobility increase in flexion, extension, lateral bending, and axial rotation, with an abnormal instantaneous center of rotation, which destabilizes the segment partly due to anterior annulus surgical removal. Disc prosthesis implantation reversed these changes in instantaneous center of rotation, but the prosthesis failed to restore the initial range of motion due to the destabilization of the ligaments in the operated disc.

INTERPRETATION

The ADDISC total disc replacement reproduces the intact disc kinematics and Instantaneous Center of Rotation, but the prosthesis fails to restore the initial range of motion due to ligament destabilization. More studies will be necessary to define a technique that restores the damaged ligaments when implanting the prosthesis.

Effects of open and minimally invasive Transforaminal Lumbar Interbody Fusion (TLIF) surgical techniques on mechanical behaviour of fused L3-L4 FSU: A comparative finite element study

For predicting the biomechanical effects of the fusion procedure, finite element (FE) analysis is widely used as a preclinical tool. Although several FE studies examined the efficacies of various fusion surgical techniques, comparative studies on Open and minimally invasive (MIS) transforaminal lumbar interbody fusion (TLIF) procedures incorporating a follower coordinate system have not been investigated yet. The current FE study evaluates the ranges of motion (ROM) and load distributions of Open-TLIF and MIS-TLIF implanted models, under physiological loading such as compression, flexion, extension and lateral bending. The most noteworthy finding from the investigation is that both the fusion procedures significantly reduced the ROMs of the implanted segment (L3-L4) and full model (L1-L5) by at least 89 % and 44 %, respectively, compared to the intact model. For all loading situations, over 95 % of the implanted models' cancellous bone volume was subjected to von Mises strains ranging from 0.0003 to 0.005. The maximum von Mises strain was observed to be localized on a small amount of cancellous bone volume (<5 %). The likelihood of adjacent segment degeneration is higher in the case of MIS-TLIF due to the higher stress (22-53 MPa) and strain (0.018-0.087) noticed on the upper facet of the L3 vertebra.

Bibliometric and Visualization Analysis of Biomechanical Research on Lumbar Intervertebral Disc

Background

Biomechanical research on the lumbar intervertebral disc (IVD) provides valuable information for the diagnosis, treatment, and prevention of related diseases, and has received increasing attention. Using bibliometric methods and visualization techniques, this study investigates for the first time the research status and development trends in this field, with the aim of providing guidance and support for subsequent research.

Methods

The Science Citation Index Expanded (SCI-Expanded) within the Web of Science Core Collection (WoSCC) database was used as the data source to select literature published from 2003 to 2022 related to biomechanical research on lumbar IVD. VOSviewer 1.6.19 and CiteSpace 6.2.R2 visualization software, as well as the online analysis platform of literature metrology, were utilized to generate scientific knowledge maps for visual display and data analysis.

Results

The United States is the most productive country in this field, with the Ulm University making the largest contribution. Wilke HJ is both the most prolific author and one of the highly cited authors, while Adams MA is the most cited author. Spine, J Biomech, Eur Spine J, Spine J, and Clin Biomech are not only the journals with the highest number of publications, but also highly cited journals. The main research topics in this field include constructing and validating three-dimensional (3D) finite element model (FEM) of lumbar spine, measuring intradiscal pressure, exploring the biomechanical effects and related risk factors of lumbar disc degeneration, studying the mechanical responses to different torque load combinations, and classifying lumbar disc degeneration based on magnetic resonance images (MRI), which are also the hot research themes in recent years.

Conclusion

This study systematically reviews the knowledge system and development trends in the field of biomechanics of lumbar IVD, providing valuable references for further research.

Study of tractor side tilt operation on intervertebral disc injury between L4 and L5 in drivers

The tilting of the cab seat when the tractor is in deep ploughing operation changes the sitting position of the driver, which may accelerate lumbar spine injury. This paper adopts the musculoskeletal model and the finite element model of the lumbar L4-L5 segment to predict the maximum Von-Mises stress and maximum strain of the driver's lumbar L4-L5 segment intervertebral disc. In this study, we used 3D motion capture to obtain the driver's spine position spatial data when the tractor tilted at different angles. A tractor-driver musculoskeletal model and a finite element model of the lumbar spine L4-L5 segments were created in AnyBody™ and Abaqus, respectively. The tractor-driver musculoskeletal model was used to calculate the load of the driver's lumbar spine L4-L5 segment at different angles of tractor tilt, which was used as the load condition of the finite element model of the lumbar spine L4-L5 segment, and then the influence of tractor tilt angle and vibration on the driver's lumbar spine L4-L5 disc was studied. The results show that the maximum Von-Mises stress and maximum strain of the driver's lumbar L4-L5 intervertebral disc will increase due to the tilt. The maximum Von-Mises stress occurs in the annulus II, and the maximum strain occurs in the upper end plate of the intervertebral disc. With the occurrence of tilt, the position of the maximum Von-Mises stress changes, which can lead to disc injury to the driver, and vibration may exacerbate this injury.

Stability and Instrumentation Stresses Among Sacropelvic Fixation Techniques With Novel Porous Fusion/Fixation Implants: A Finite Element Study

BACKGROUND

Sacropelvic fixation is frequently combined with thoracolumbar instrumentation for correcting spinal deformities. This study aimed to characterize sacropelvic fixation techniques using novel porous fusion/fixation implants (PFFI).

METHODS

Three T10-pelvis finite element models were created: (1) pedicle screws and rods in T10-S1, PFFI bilaterally in S2 alar-iliac (S2AI) trajectory; (2) fixation in T10-S1, PFFI bilaterally in S2AI trajectory, triangular implants bilaterally above the PFFI in a sacro-alar-iliac trajectory (PFFI-IFSAI); and (3) fixation in T10-S1, PFFI bilaterally in S2AI trajectory, PFFI in sacro-alar-iliac trajectory stacked cephalad to those in S2AI position (2-PFFI). Models were loaded with pure moments of 7.5 Nm in flexion-extension, lateral bending, and axial rotation. Outputs were compared against 2 baseline models: (1) pedicle screws and rods in T10-S1 (PED), and (2) pedicle screws and rods in T10-S1, and S2AI screws.

RESULTS

PFFI and S2AI resulted in similar L5-S1 motion; adding another PFFI per side (2-PFFI) further reduced this motion. Sacroiliac joint (SIJ) motion was also similar between PFFI and S2AI; PFFI-IFSAI and 2-PFFI demonstrated a further reduction in SIJ motion. Additionally, PFFI reduced max stresses on S1 pedicle screws and on implants in the S2AI position.

CONCLUSION

The study shows that supplementing a long construct with PFFI increases the stability of the L5-S1 and SIJ and reduces stresses on the S1 pedicle screws and implants in the S2AI position.

CLINICAL RELEVANCE

The findings suggest a reduced risk of pseudarthrosis at L5-S1 and screw breakage. Clinical studies may be performed to demonstrate applicability to patient outcomes.

LEVEL OF EVIDENCE

Not applicable (basic science study).

Comparison of biomechanical parameters of two Chinese cervical spine rotation manipulations based on motion capture and finite element analysis

Objective: The purpose of this study was to obtain the stress-strain of the cervical spine structure during the simulated manipulation of the oblique pulling manipulation and the cervical rotation-traction manipulation in order to compare the mechanical mechanism of the two manipulations. Methods: A motion capture system was used to record the key kinematic parameters of operating the two manipulations. At the same time, a three-dimensional finite element model of the C0-T1 full healthy cervical spine was established, and the key kinematic parameters were loaded onto the finite element model in steps to analyze and simulate the detailed process of the operation of the two manipulations. Results: A detailed finite element model of the whole cervical spine including spinal nerve roots was established, and the validity of this 3D finite element model was verified. During the stepwise simulation of the two cervical spine rotation manipulations to the right, the disc (including the annulus fibrosus and nucleus pulposus) and facet joints stresses and displacements were greater in the oblique pulling manipulation group than in the cervical rotation-traction manipulation group, while the spinal cord and nerve root stresses were greater in the cervical rotation-traction manipulation group than in the oblique pulling manipulation group. The spinal cord and nerve root stresses in the cervical rotation-traction manipulation group were mainly concentrated in the C4/5 and C5/6 segments. Conclusion: The oblique pulling manipulation may be more appropriate for the treatment of cervical spondylotic radiculopathy, while cervical rotation-traction manipulation is more appropriate for the treatment of cervical spondylosis of cervical type. Clinicians should select cervical rotation manipulations for different types of cervical spondylosis according to the patient's symptoms and needs.

Biomechanical effect of endplate defects on the intermediate vertebral bone in consecutive two-level anterior cervical discectomy and fusion: a finite element analysis

BACKGROUND

Intermediate vertebral collapse is a newly discovered complication of consecutive two-level anterior cervical discectomy and fusion (ACDF). There have been no analytical studies related to the effects of endplate defects on the biomechanics of the intermediate vertebral bone after ACDF. This study aimed to compare the effects of endplate defects on the intermediate vertebral bone biomechanics in the zero-profile (ZP) and cage-and-plate (CP) methods of consecutive 2-level ACDF and to determine whether collapse of the intermediate vertebra is more likely to occur using ZP.

METHODS

A three-dimensional finite element (FE) model of the intact cervical spine (C2-T1) was constructed and validated. The intact FE model was then modified to build ACDF models and imitate the situation of endplate injury, establishing two groups of models (ZP, IM-ZP and CP, IM-ZP). We simulated cervical motion, such as flexion, extension, lateral bending and axial rotation, and compared the range of motion (ROM), upper and lower endplate stress, fusion fixation device stress, C5 vertebral body stress, intervertebral disc internal pressure (intradiscal pressure, or IDP) and the ROM of adjacent segments in the models.

RESULTS

There was no significant difference between the IM-CP model and the CP model in the ROM of the surgical segment, upper and lower endplate stress, fusion fixation device stress, C5 vertebral body stress, IDP, or ROM of the adjacent segments. Compared with the CP model, the endplate stress of the ZP model is significantly higher in the flexion, extension, lateral bending and axial rotation conditions. Compared with the ZP model, endplate stress, screw stress, C5 vertebral stress and IDP in IM-ZP were significantly increased under flexion, extension, lateral bending and axial rotation conditions.

CONCLUSIONS

Compared to consecutive 2-level ACDF using CP, collapse of the intermediate vertebra is more likely to occur using ZP due to its mechanical characteristics. Intraoperative endplate defects of the anterior lower margin of the middle vertebra are a risk factor leading to collapse of the middle vertebra after consecutive 2-level ACDF using ZP.

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.

ADDISC lumbar disc prosthesis: Analytical and FEA testing of novel implants

The intact intervertebral disc is a six-freedom degree elastic deformation structure with shock absorption. "Ball-and-socket" TDR do not reproduce these properties inducing zygapophyseal joint overload. Elastomeric TDRs reproduce better normal disc kinematics, but repeated core deformation causes its degeneration. We aimed to create a new TDR (ADDISC) reproducing healthy disc features. We designed TDR, analyzed (Finite Element Analysis), and measured every 500,000 cycles for 10 million cycles of the flexion-extension, lateral bending, and axial rotation cyclic compression bench-testing. In the inlay case, we weighted it and measured its deformation. ADDISC has two semi-spherical articular surfaces, one rotation centre for flexion, another for extension, the third for lateral bending, and a polycarbonate urethane inlay providing shock absorption. The first contact is between PCU and metal surfaces. There is no metal-metal contact up to 2000 N, and CoCr28Mo6 absorbs the load. After 10 million cycles at 1.2-2.0 kN loads, wear 140.96 mg (35.50 mm³), but no implant failures. Our TDR has a physiological motion range due to its articular surfaces' shape and the PCU inlay bumpers, minimizing the facet joint overload. ADDISC mimics healthy disc biomechanics and Instantaneous Rotation Center, absorbs shock, reduces wear, and has excellent long-term endurance.

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.

Computational lumbar spine models: A literature review

BACKGROUND

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

METHODS

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

FINDINGS

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

INTERPRETATION

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

Effect of sagittal screw angle and distance of screw apex to superior endplate on adjacent segment disease after posterolateral lumbar fusion: a retrospective study

Background

Numerous complications of lumbar fusion surgery have been reported, with adjacent segment disease (ASD) being one of the most important. Few studies describe the effect of sagittal, horizontal screw angles and distance of pedicle screw apex to superior endplate on the incidence of ASD in lumbar spine. The purpose of this retrospective study is to evaluate the hypothesis that unsatisfactory pedicle screw insertion positions would increase the likelihood of ASD.

Methods

Outpatients with lumbar spinal stenosis underwent posterolateral lumbar fusion at L4-S1 with a least 2-year follow-up were studied. ASD at L3–L4 was defined as a condition in which intervertebral disk narrowing, posterior vertebral opening, and vertebral slippage progress at the last follow-up in comparison with the postoperative. Independent t test was performed to compare data between two groups; Spearman analysis was performed to analyze the relationship between two continuous variables. Multivariate binary logistic models were performed to identify the independent risk factors of ASD. The receiver operating characteristic (ROC) curve was performed to measure model discrimination and Hosmer–Lemeshow (H–L) test was used to measure calibration. ROC curve evaluated the discrimination ability of sagittal screw angle and distance in predicting incidence of ASD.

Results

Patients in ASD group exhibit higher incidence of osteoporosis, higher Visual analogue scale (VAS), Oswestry disability index (ODI), bigger sagittal screw angle, shorter distance of pedicle screw apex to superior endplate than those in non-ASD group (p < 0.05). VAS, ODI at the last follow-up were positively correlated with Pfirrmann grade of L3–4 disk and sagittal screw angle, while negatively correlated with distance of screw apex to superior endplate (p < 0.05). Multivariate binary logistic model indicated that follow-up time (odds ratio [OR] 1.637, 95% confidence interval [CI] 1.186–2.260), distance of screw apex to superior endplate (OR 0.150, 95% CI 0.067–0.336), sagittal screw angle (OR 2.404, 95% CI 1.608–3.594) were statistically significant. The models showed great discrimination and calibration. The area under the curve of ASD identified by sagittal angle and distance was 0.895 and the cut-off values were 5.500° and 6.250 mm, respectively.

Conclusion

Sagittal screw angle and distance of screw apex to superior endplate were significantly associated with the risk of ASD.

The application of finite element analysis to determine the optimal UIV of growing-rod treatment in early-onset scoliosis

Objectives: To analyze the stress distribution in the proximal vertebral body and soft tissue of dual growing-rod (GR) with different upper instrumented vertebra (UIV) to determine the optimal UIV.

Methods: A ten-year-old male EOS case treated with GR was selected. Based on spiral computed tomography (CT) scanning performed in 0.6 mm thick slices, a finite element model (FEM) of the preoperative state (M0, the original spine state) of the patient was created. Subsequently, four models with different UIV fixations were numerically analyzed by FEM, including M1 (UIV = T1, i.e., the upper-end vertebrae (UEV) of the upper thoracic curve), M2 (UIV = T2), M3 (UIV = T3) and M4 (UIV = T4, i.e., the lower end vertebrae (LEV) of the upper thoracic curve). Displacement and maximum stress in the proximal vertebral body and soft tissue were measured and compared among the five models.

Results: The spine model was fixed with the sacrum, and the gravity conditions were imposed on each vertebral body according to the research of Clin and Pearsall. The results are as follows:M4 model has the largest overall displacement, while M1 has the least displacement among the four models. Except M2, the maximum normalized stress of UIV increases with the downward movement of UIV. M1 has the lowerest annulus fibrosus stress and highest joint capsule stress, which is characterized by the vertebrae backward leaning, while M4 is the opposite. The supraspinous ligament stress of M3 and M4 is significantly higher than that of M1 and M2. This suggests that UIV downshift increases the tendency of the proximal vertebral bodies to bend forward, thereby increasing the tension of the posterior ligaments (PL).

Conclusion: The UIV of the GR is recommended to be close to the UEV of the upper thoracic curve, which can reduce the stress of the proximal PL, thereby reducing the occurrence of proximal junctional kyphosis (PJK).

Patient-Specific Finite Element Modeling of the Whole Lumbar Spine Using Clinical Routine Multi-Detector Computed Tomography (MDCT) Data-A Pilot Study

(1) Background: To study the feasibility of developing finite element (FE) models of the whole lumbar spine using clinical routine multi-detector computed tomography (MDCT) scans to predict failure load (FL) and range of motion (ROM) parameters. (2) Methods: MDCT scans of 12 subjects (6 healthy controls (HC), mean age ± standard deviation (SD): 62.16 ± 10.24 years, and 6 osteoporotic patients (OP), mean age ± SD: 65.83 ± 11.19 years) were included in the current study. Comprehensive FE models of the lumbar spine (5 vertebrae + 4 intervertebral discs (IVDs) + ligaments) were generated (L1−L5) and simulated. The coefficients of correlation (ρ) were calculated to investigate the relationship between FE-based FL and ROM parameters and bone mineral density (BMD) values of L1−L3 derived from MDCT (BMDQCT-L1-3). Finally, Mann−Whitney U tests were performed to analyze differences in FL and ROM parameters between HC and OP cohorts. (3) Results: Mean FE-based FL value of the HC cohort was significantly higher than that of the OP cohort (1471.50 ± 275.69 N (HC) vs. 763.33 ± 166.70 N (OP), p < 0.01). A strong correlation of 0.8 (p < 0.01) was observed between FE-based FL and BMDQCT-L1-L3 values. However, no significant differences were observed between ROM parameters of HC and OP cohorts (p = 0.69 for flexion; p = 0.69 for extension; p = 0.47 for lateral bending; p = 0.13 for twisting). In addition, no statistically significant correlations were observed between ROM parameters and BMDQCT- L1-3. (4) Conclusions: Clinical routine MDCT data can be used for patient-specific FE modeling of the whole lumbar spine. ROM parameters do not seem to be significantly altered between HC and OP. In contrast, FE-derived FL may help identify patients with increased osteoporotic fracture risk in the future.

Biomechanical effects of an oblique lumbar interbody fusion combined with posterior augmentation: a finite element analysis

Background

Oblique lateral interbody fusion (OLIF) is widely used to treat lumbar degenerative disc disease. This study aimed to evaluate the biomechanical stability of OLIF, OLIF including posterior pedicle screw and rod (PSR), and OLIF including cortical screw and rod (CSR) instrumentation through finite element analysis.

Methods

A complete L2-L5 finite element model of the lumbar spine was constructed. Surgical models of OLIF, such as stand-alone, OLIF combined with PSR, and OLIF combined with CSR were created in the L3-L4 surgical segments. Range of motion (ROM), end plate stress, and internal fixation peak stress were compared between different models under the same loading conditions.

Results

Compared to the intact model, ROM was reduced in the OLIF model under all loading conditions. The surgical models in order of increasing ROM were PSR, CSR, and stand-alone; however, the difference in ROM between BPS and CSR was less than 0.4° and was not significant under any loading conditions. The stand-alone model had the highest stress on the superior L4 vertebral body endplate under all loading conditions, whereas the end plate stress was relatively low in the BPS and CSR models. The CSR model had the highest internal fixation stress, concentrated primarily at the end of the screw.

Conclusions

OLIF alone significantly reduces ROM but does not provide sufficient stability. Addition of posterior PSR or CSR internal fixation instrumentation to OLIF surgery can significantly improve biomechanical stability of the segment undergoing surgery.

Development of a Finite Element Model of the Pediatric Thoracic and Lumbar Spine, Ribcage, and Pelvis with Orthotropic Region-Specific Vertebral Growth

Finite element (FE) modeling of the spine has increasingly been applied in orthopedic precision-medicine approaches. Previously published FE models of the pediatric spine growth have made simplifications in geometry of anatomical structures, material properties, and representation of vertebral growth. To address those limitations, a comprehensive FE model of a pediatric (10-year-old) osteo-ligamentous thoracic and lumbar spine (T1-L5 with intervertebral discs (IVDs) and ligaments), ribcage, and pelvis with age- and level-specific ligament properties and orthotropic region-specific vertebral growth was developed and validated. Range of motion (ROM) measures, namely lateral bending, flexion-extension, and axial rotation, of the current 10 YO FE model were generally within reported ranges of scaled in vitro adult ROM data. Changes in T1-L5 spine height, as well as kyphosis (T2-T12) and lordosis (L1-L5) angles in the current FE model for two years of growth (from ages 10 to 12 years) were within ranges reported from corresponding pediatric clinical data. The use of such comprehensive pediatric FE models can provide clinically relevant insights into normative and pathological biomechanical responses of the spine, and also contribute to the development and optimization of clinical interventions for spine deformities.

Optimization of a lumbar interspinous fixation device for the lumbar spine with degenerative disc disease

Interspinous spacer devices used in interspinous fixation surgery remove soft tissues in the lumbar spine, such as ligaments and muscles and may cause degenerative diseases in adjacent segments its stiffness is higher than that of the lumbar spine. Therefore, this study aimed to structurally and kinematically optimize a lumbar interspinous fixation device (LIFD) using a full lumbar finite element model that allows for minimally invasive surgery, after which the normal behavior of the lumbar spine is not affected. The proposed healthy and degenerative lumbar spine models reflect the physiological characteristics of the lumbar spine in the human body. The optimum number of spring turns and spring wire diameter in the LIFD were selected as 3 mm and 2 turns, respectively—from a dynamic range of motion (ROM) perspective rather than a structural maximum stress perspective—by applying a 7.5 N∙m extension moment and 500 N follower load to the LIFD-inserted lumbar spine model. As the spring wire diameter in the LIFD increased, the maximum stress generated in the LIFD increased, and the ROM decreased. Further, as the number of spring turns decreased, both the maximum stress and ROM of the LIFD increased. When the optimized LIFD was inserted into a degenerative lumbar spine model with a degenerative disc, the facet joint force of the L3-L4 lumbar segment was reduced by 56%–98% in extension, lateral bending, and axial rotation. These results suggest that the optimized device can strengthen the stability of the lumbar spine that has undergone interspinous fixation surgery and reduce the risk of degenerative diseases at the adjacent lumbar segments.

A comparative finite element analysis of artificial intervertebral disc replacement and pedicle screw fixation of the lumbar spine

Titanium alloy-based Pedicle screw-rod fusion is a very common technique to provide higher fusion regularity than other methods. In recent times, Carbon-fibre-reinforced (CFR)-PEEK rod is used to better reduce the fusion rate. Alternatively, total disc replacement (TDR) is also very common for the non-fusion treatment method for degenerative disc disease (DDD). This study aims to investigate flexibility (ROM), stability, stress condition in implant, implant adjacent bone of the implanted lumbar spine during different physiological movements and loading environments. The finite element (FE) intact model of the lumbar spine (L2-L5) with two-level pedicle screw-rod fusion at L3-L4-L5 and two-level artificial disc replacement was developed. CFR-PEEK was taken for rod for semi-rigid fusion. UHMWPE was taken as core part of the artificial disc. The FE models were simulated under 6, 8 and 10 Nm moments in left right lateral bending, flexion and extension movements. The total ROM was reduced for two-level pedicle screw fixation and increased for the artificial disc replacement model during flexion extension compared to the intact spine. The total ROM was reduced by around 54% and 25% for two-level fixation and increased by 30% and 19.5% for artificial disc replacement spine, under flexion-extension and left-right lateral bending respectively. For screw fixation, the ROM increased by 15% and 18% reduced by 4.5% and 14% for disc replacement at the adjacent segments for flexion-extension and left-right lateral bending.

Low-Modulus PMMA Has the Potential to Reduce Stresses on Endplates after Cement Discoplasty

Cement discoplasty has been developed to treat patients with advanced intervertebral disc degeneration. In discoplasty, poly(methylmethacrylate) (PMMA) bone cement is injected into the disc, leading to reduced pain and certain spinal alignment correction. Standard PMMA-cements have much higher elastic modulus than the surrounding vertebral bone, which may lead to a propensity for adjacent fractures. A PMMA-cement with lower modulus might be biomechanically beneficial. In this study, PMMA-cements with lower modulus were obtained using previously established methods. A commercial PMMA-cement (V-steady^®, G21 srl) was used as control, and as base cement. The low-modulus PMMA-cements were modified by 12 vol% (LA12), 16 vol% (LA16) and 20 vol% (LA20) linoleic acid (LA). After storage in 37 °C PBS from 24 h up to 8 weeks, specimens were tested in compression to obtain the material properties. A lower E-modulus was obtained with increasing amount of LA. However, with storage time, the E-modulus increased. Standard and low-modulus PMMA discoplasty were compared in a previously developed and validated computational lumbar spine model. All discoplasty models showed the same trend, namely a substantial reduction in range of motion (ROM), compared to the healthy model. The V-steady model had the largest ROM-reduction (77%), and the LA20 model had the smallest (45%). The average stress at the endplate was higher for all discoplasty models than for the healthy model, but the stresses were reduced for cements with higher amounts of LA. The study indicates that low-modulus PMMA is promising for discoplasty from a mechanical viewpoint. However, validation experiments are needed, and the clinical setting needs to be further considered.

Investigation of Reaction Forces in the Thoracolumbar Fascia during Different Activities: A Mechanistic Numerical Study

Spinal instability remains a complex phenomenon to study while the cause of low back pain continues to challenge researchers. The role of fascia in biomechanics adds to the complexity of spine biomechanics but offers a new window from which to investigate our spines. Specifically, the thoracolumbar fascia may have an important role in spine biomechanics, and thus the purpose of this study was to access the mechanical influence of the thoracolumbar fascia on spine biomechanics during different simulated activities. A numerical finite element model of the lumbar spine inclusive of the intra-abdominal and intra-muscular regions as well as the thoracolumbar fascia was constructed and validated. Four different loading scenarios were simulated while deformation, stress, pressure, and reaction forces between the thoracolumbar fascia and spine were measured. Model validation was accomplished through comparison to in vivo and ex vivo published studies. Force transmission between the thoracolumbar fascia and the spine increased 40% comparing kyphotic and squatting lifting patterns. Further, the importance of reciprocating paraspinal and intra-abdominal pressures was demonstrated. It was also found that tension in the thoracolumbar fascia remains even in a simulated prone position. This numerical analysis allowed for an objective interpretation of the loads conveyed through the thoracolumbar fascia in different positional or lifting scenarios. Based on validation studies, it would appear to be a viable experimental platform from which insight can be derived. The loads in the thoracolumbar fascia vary considerably based on simulated tasks and are linked to the pressures in the paraspinal and intra-abdominal regions.

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

Background

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

Methods

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

Results

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

Conclusions

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

Sensitivity of Intervertebral Disc Finite Element Models to Internal Geometric and Non-geometric Parameters

Finite element models are useful for investigating internal intervertebral disc (IVD) behaviours without using disruptive experimental techniques. Simplified geometries are commonly used to reduce computational time or because internal geometries cannot be acquired from CT scans. This study aimed to (1) investigate the effect of altered geometries both at endplates and the nucleus-anulus boundary on model response, and (2) to investigate model sensitivity to material and geometric inputs, and different modelling approaches (graduated or consistent fibre bundle angles and glued or cohesive inter-lamellar contact). Six models were developed from 9.4 T MRIs of bovine IVDs. Models had two variations of endplate geometry (a simple curved profile from the centre of the disc to the periphery, and precise geometry segmented from MRIs), and three variations of NP-AF boundary (linear, curved, and segmented). Models were subjected to axial compressive loading (to 0.86 mm at a strain rate of 0.1/s) and the effect on stiffness and strain distributions, and the sensitivity to modelling approaches was investigated. The model with the most complex geometry (segmented endplates, curved NP-AF boundary) was 3.1 times stiffer than the model with the simplest geometry (curved endplates, linear NP-AF boundary), although this difference may be exaggerated since segmenting the endplates in the complex geometry models resulted in a shorter average disc height. Peak strains were close to the endplates at locations of high curvature in the segmented endplate models which were not captured in the curved endplate models. Differences were also seen in sensitivity to material properties, graduated fibre angles, cohesive rather than glued inter-lamellar contact, and NP:AF ratios. These results show that FE modellers must take care to ensure geometries are realistic so that load is distributed and passes through IVDs accurately.

Prediction of biomechanical responses of human lumbar discs - a stochastic finite element model analysis

BACKGROUND

Accurate biomechanical investigation of human intervertebral discs (IVDs) is difficult because of their complicated structural and material features.

AIM

To investigate probabilistic distributions of the biomechanical responses of the IVD by considering varying nonlinear structural and material properties using a stochastic finite element (FE) model.

METHODS

A FE model of a L3-4 disc was reconstructed, including the nucleus pulposus (NP), annular matrix and fibers. A Monte Carlo method was used to randomly generate 500 sets of the nonlinear material properties and fiber orientations of the disc that were implemented into the FE model. The FE model was analyzed under seven loading conditions: a 500 N compressive force, a 7.5Nm moment simulating flexion, extension, left-right lateral bending, and left-right axial rotation, respectively. The distributions of the ranges of motion (ROMs), intradiscal pressures (IDP), fiber stresses and matrix strains of the disc were analyzed.

RESULTS

Under the compressive load, the displacement varied between 0.29 mm and 0.76 mm. Under the 7.5Nm moment, the ROMs varied between 3.0° and 6.0° in primary rotations. The IDPs varied within 0.3 MPa under all the loading conditions. The maximal fiber stress (3.22 ± 0.64 MPa) and matrix strain (0.27 ± 0.12%) were observed under the flexion and extension moments, respectively.

CONCLUSION

The IVD biomechanics could be dramatically affected by the structural and material parameters used to construct the FE model. The stochastic FE model that includes the probabilistic distributions of the structural and material parameters provides a useful approach to analyze the statistical ranges of the biomechanical responses of the IVDs.

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.

Numerical analysis of the mechanical behaviour of intact and implanted lumbar functional spinal units: Effects of loading and boundary conditions

The objective of this study was to develop an improved finite element (FE) model of a lumbar functional spinal unit (FSU) and to subsequently analyse the deviations in load transfer owing to implantation. The effects of loading and boundary conditions on load transfer in intact and implanted FSUs and its relationship with the potential risk of vertebral fracture were investigated. The FE models of L1-L5 and L3-L4 FSUs, intact and implanted, were developed using patient-specific CT-scan dataset and segmentation of cortical and cancellous bone regions. The effect of submodelling technique, as compared to artificial boundary conditions, on the elastic behaviour of lumbar spine was examined. Applied forces and moments, corresponding to physiologic movements, were used as loading conditions. Results indicated that the loading and boundary conditions considerably affect stress-strain distributions within a FSU. This study, based on an improved FE model of a vertebra, highlights the importance of using the submodelling technique to adequately evaluate the mechanical behaviour of a FSU. In the intact FSU, strains of 200-400 µε were observed in the cancellous bone of vertebral body and pedicles. High equivalent stresses of 10-25 MPa and 1-5 MPa were generated around the pars interarticularis for cortical and cancellous regions, respectively. Implantation caused reductions of 85%-92% in the range of motion for all movements. Insertion of the intervertebral cage resulted in major deviations in load transfer across a FSU for all movements. The cancellous bone around cage experienced pronounced increase in stresses of 10-15 MPa, which indicated potential risk of failure initiation in the vertebra.

Effect of single and multilevel artificial inter-vertebral disc replacement in lumbar spine: A finite element study

Degenerative disc disease (DDD) in lumbar spine is one of the major musculoskeletal disorders that cause low back pain (LBP). The intervertebral disc structure and dynamics of the lumbar spine are significantly affected by lumbar DDD, leading to a reduced range of motion (ROM), muscle weakness and gradual degradation. Spinal fusion and inter-vertebral disc replacement prostheses are two major surgical methods used for treating lumbar DDD. The aim of this present study is to examine biomechanical impacts of single level (L3-L4 and L4-L5) and multi level (L3-L4-L5) inter-vertebral disc replacement in lumbar spine (L2-L5) and to compare the performance with intact spine. Finite element (FE) analysis has been used to compare the mobility and stress distribution of all the models for four physiological movements, namely flexion, extension, left and right lateral bending under 6, 8 and 10 Nm moments. Spinal fusion implants completely restrict the motion of the implanted segment and increase disc stress at the adjacent levels. In contrast to that, the results single level ADR models showed closer ROM and disc stress to natural model. At the spinal segments adjacent to the implantation, single level ADR shows lower chance of disc degeneration. However, significantly increased ROM was observed in case of double level ADR.

Biomechanical effects of laminectomies in the human lumbar spine: a finite element study

BACKGROUND CONTEXT

Previous studies have analyzed the effect of laminectomy on intervertebral disc (IVD), facet-joint-forces (FJF), and range of motion (ROM), while only two have specifically reported stresses at the pars interarticularis (PI) with posterior element resection. These studies have been performed utilizing a single subject, questioning their applications to a broader population.

PURPOSE

We investigate the effect of graded PI resection in a three-dimensional manner on PI stress to provide surgical guidelines for avoidance of iatrogenic instability following lumbar laminectomy. Additionally, quantified FJF and IVD stresses can provide further insight into the development of adjacent segment disease.

STUDY DESIGN

Biomechanical finite element (FE) method investigation of the lumbar spine.

METHODS

FE models of the lumbar spine of three subjects were created using the open-source finite element software, FEBio. Single-level laminectomy, two-level laminectomy, and ventral-to-dorsal PI resection simulations were performed with varying degrees of PI resection from 0% to 75% of the native PI. These models were taken through cardinal ROM under standard loading conditions and PI stresses, FJF, IVD stresses, and ROM were analyzed.

RESULTS

The three types of laminectomy simulated in this study showed a nonlinear increase in PI stress with increased bone resection. Axial rotation generated the most stress at the PI followed by flexion, extension and lateral bending. At 75% bone resection all three types of laminectomy produced PI stresses that were near the ultimate strength of human cortical bone during axial rotation. FJF decreased with increased bone resection for the three laminectomies simulated. This was most notable in axial rotation followed by extension and lateral bending. IVD stresses varied greatly between the nonsurgical models and likewise the effect of laminectomy on IVD stresses varied between subjects. ROM was mostly unaffected by the laminectomies performed in this study.

CONCLUSIONS

Regarding the risk of iatrogenic spondylolisthesis, the combined results are sufficient evidence to suggest surgeons should be particularly cautious when PI resection exceeds 50% bone resection for all laminectomies included in this study. Lastly, the effects seen in FJF and IVD stresses indicate the degree to which the remainder of the spine must experience compensatory biomechanical changes as a result of the surgical intervention.

Development and validation of a timely and representative finite element human spine model for biomechanical simulations

Numerous spine Finite Element (FE) models have been developed to assess spinal tolerances, spinal loadings and low back pain-related issues. However, justified simplifications, in terms of tissue decomposition and inclusion, for such a complex system may overlook crucial information. Thus, the purpose of this research was to develop and validate a comprehensive and representative spine FE model inclusive of an accurate representation of all major torso elements. A comprehensive model comprised of 273 tissues was developed via a novel FE meshing method to enhance computational feasibility. A comprehensive set of indirect validation tests were carried out to validate every aspect of the model. Under an increasing angular displacement of 24°-41°, the lumbar spine recorded an increasing moment from 5.5 to 9.3 Nm with an increase in IVD pressures from 0.41 to 0.66 MPa. Under forward flexion, vertical vertebral displacements simulated a 6% and 13% maximum discrepancy for intra-abdominal and intramuscular pressure results, all closely resembling previously documented in silico measured values. The developed state-of-the-art model includes most physiological tissues known to contribute to spinal loadings. Given the simulation's accuracy, confirmed by its validation tests, the developed model may serve as a reliable spinal assessment tool.

Biomechanical analysis of lumbar spine with interbody fusion surgery and U-shaped lumbar interspinous spacers

Previous research indicates whole-body vibration may lead to low back pain. The aim of this study is assessing the dynamic characteristics of a lumbar spine with Coflex and Coflex-F (commercial implants used as lumbar interspinous spacers) and effect of lumbar interbody fusion surgery. A transient dynamic analysis is performed on three numerical lumbar spine models under the loading condition of a vertical sinusoidal force of ±40 N with a compressive follower preload of 400 N. Also, Coflex-F model with and without interbody fusion surgery is analyzed under the same loading condition. The results show that the maximum value and vibration amplitude of von Mises stress in annulus ground substance (AGS) and intradiscal pressure (IDP) at implanted segment decrease from healthy model to Coflex model, and Coflex-F model. By contrast, for adjacent segments the maximum value of implanted models are larger than that of healthy model. The maximum value of endplates with and without cage are 2.44 MPa and 1.73 MPa (L4 inferior endplate), 1.94 MPa and 1.42 MPa (L5 superior endplate), respectively. The vibration amplitude of Coflex-F model with fusion surgery is smaller than that without fusion surgery. Coflex and Coflex-F not only protect implanted segment but also have a negative effect on adjacent segments. Inserting cage for Coflex-F model can absorb vibration energy at adjacent segments.

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

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

Simulation of L-4 lumbar spine model of motorist exposed to vibration from speed hump

BACKGROUND

The motorcycle is often used in recurring travel between locations, dense traffic, poor conditioned roads and thus the repetitive loading on the musculoskeletal system of the rider leads to risk factors associated with musculoskeletal disorders. This study was aimed to examine the biomechanical response of the L-4 lumbar spine while riding the motorcycle on the speed hump at 20 km/h.

METHODS

Three-dimensional (3D) model of the L-4 lumbar spine was reconstructed based on the CT scan data obtained from the subjects. Material properties of the L-4 lumbar spine were assumed to be isotropic and homogenous. Mesh convergence and sensitivity analyses were performed and validated before simulation. Static and dynamic analyses were accomplished using quasi-static and steady-state dynamic analyses.

RESULTS

Static analysis results show that the highest stress concentrations were found around the pedicle and spinal canal. It is an expected commonplace for injuries because of loading. The dynamic simulation results showed the major resonance of the L-4 lumbar spine model is about 8-40 Hz. The stress, displacement, velocity, and acceleration value declines beyond 40 Hz as the frequency increases.

CONCLUSIONS

The simulation specifies the symmetric and unsymmetrical distributions of vibration magnitude regions of the lumbar spine. This study provides the modelling of the lumbar spine (L-4) and validated the effect of overloading failure as well as identified the biomechanical behaviour.

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.

The use of triangular implants to enhance sacropelvic fixation: a finite element investigation

BACKGROUND CONTEXT

Long thoracolumbar fixation and fusion have become a consolidated treatment for severe spinal disorders. Concomitant sacropelvic fixation with S2 alar-iliac (S2AI) screws is frequently performed to limit instrumentation failure and pseudarthrosis at the lumbosacral junction.

PURPOSE

This study explored the use of triangular titanium implants in different configurations in which the implants supplemented standard sacropelvic fixation with S2AI screws in order to further increase the stability of S2AI fixation.

STUDY DESIGN

Finite element study.

METHODS

Four T10-pelvis instrumented models were built: pedicle screws and rods in T10-S1 (PED); pedicle screws and rods in T10-S1, and bilateral S2 alar-iliac screws (S2AI); pedicle screws and rods in T10-S1, bilateral S2AI screws, and triangular implants inserted bilaterally in a sacral alar-iliac trajectory (Tri-SAI); pedicle screws and rods in T10-S1, bilateral S2AI screws and two bilateral triangular titanium implants inserted in a lateral trajectory (Tri-Lat). The models were tested under pure moments of 7.5 Nm in flexion-extension, lateral bending and axial rotation.

RESULTS

SIJ motion was reduced by 50% to 66% after S2AI fixation; the addition of triangular titanium implants in either a SAI or a lateral trajectory further reduced it. S2AI, Tri-SAI, and Tri-Lat resulted in significantly lower stresses in S1 pedicle screws when compared to PED. Triangular implants had a protective effect on the maximal stresses in S2AI screws, especially when placed in the SAI trajectory. Sacropelvic fixation did not have any protective effect on the posterior rods.

CONCLUSIONS

Supplementing S2AI screws with triangular implants had a protective effect on the S2AI screws themselves, as well as the S1 pedicle screws, in the tested model.

CLINICAL SIGNIFICANCE

Triangular implants can substantially reduce the residual flexibility of the SIJ with respect to S2AI fixation alone, suggesting a possible role in patients needing reinforced fixation. In vivo investigation is needed to determine if these in vitro effects translate into clinically important differences.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

ICR in human cadaveric specimens: An essential parameter to consider in a new lumbar disc prosthesis design

STUDY DESIGN

Biomechanical study in cadaveric specimens.

BACKGROUND

The commercially available lumbar disc prostheses do not reproduce the intact disc's Instantaneous centre of Rotation (ICR), thus inducing an overload on adjacent anatomical structures, promoting secondary degeneration.

AIM

To examine biomechanical testing of cadaveric lumbar spine specimens in order to evaluate and define the ICR of intact lumbar discs.

MATERIAL AND METHODS

Twelve cold preserved fresh human cadaveric lumbosacral spine specimens were subjected to computerized tomography (CT), magnetic resonance imaging (MRI) and biomechanical testing. Kinematic studies were performed to analyse range of movements in order to determine ICR.

RESULTS

Flexoextension and lateral bending tests showed a positive linear correlation between the angle rotated and the displacement of the ICR in different axes.

DISCUSSION

ICR has not been taken into account in any of the available literature regarding lumbar disc prosthesis. Considering our results, neither the actual ball-and-socket nor the withdrawn elastomeric nucleus models fit the biomechanics of the lumbar spine, which could at least in part explain the failure rates of the implants in terms of postoperative failed back syndrome (low back pain). It is reasonable to consider then that an implant should also adapt the equations of the movement of the intact ICR of the joint to the post-surgical ICR.

CONCLUSIONS

This is the first cadaveric study on the ICR of the human lumbar spine. We have shown that it is feasible to calculate and consider this parameter in order to design future prosthesis with improved clinical and biomechanical characteristics.