Mechanical response of a lumbar motion segment in axial torque alone and combined with compression

The effect of polymethylmethacrylate augmentation on the primary stability of stand-alone implant construct versus posterior stabilization in oblique lumbar interbody fusion with osteoporotic bone quality- a finite element study

BACKGROUND CONTEXT

Oblique lumbar interbody fusion (OLIF) can provide an ideal minimally invasive solution for achieving spinal fusion in an older, more frail population where decreased bone quality can be a limiting factor. Stabilization can be achieved with bilateral pedicle screws (BPS), which require additional incisions and longer operative time. Alternatively, a novel self-anchoring stand-alone lateral plate system (SSA) can be used, where no additional incisions are required. Based on the relevant literature, BPS constructs provide greater primary biomechanical stability compared to lateral plate constructs, including SSA. This difference is further increased by osteoporosis. Screw augmentation in spinal fusion surgeries is commonly used; however, in the case of OLIF, it is a fairly new concept, lacking a consensus-based guideline.

PURPOSE

This comparative finite element (FE) study aimed to investigate the effect of PMMA screw augmentation on the primary stability of a stand-alone implant construct versus posterior stabilization in OLIF with osteoporotic bone quality.

STUDY DESIGN

The biomechanical effect of screw augmentation was studied inside an in-silico environment using computer-aided FE analysis.

METHODS

A previously validated and published L2-L4 FE model with normal and osteoporotic bone material properties was used. Geometries based on the OLIF implants (BPS, SSA) were created and placed inside the L3-L4 motion segment with increasing volumes (1-6 cm3) of PMMA augmentation. A follower load of 400 N and 10 Nm bending moment (in the three anatomical planes) were applied to the surgical FE models with different bone material properties. The operated L3-L4 segmental range of motion (ROM), the inserted cage's maximal caudal displacements, and L4 cranial bony endplate principal stress values were measured.

RESULTS

The nonaugmented values for the BPS construct were generally lower compared to SSA, and the difference was increased by osteoporosis. In osteoporotic bone, PMMA augmentation gradually decreased the investigated parameters and the difference between the two constructs as well. Between 3 cm3 and 4 cm3 of injected PMMA volume per screw, the difference between augmented SSA and standard BPS became comparable.

CONCLUSIONS

Based on this study, augmentation can enhance the primary stability of the constructs and decrease the difference between them. Considering leakage as a possible complication, between 3 cm3 and 4 cm3 of injected PMMA per screw can be an adequate amount for SSA augmentation. However, further in silico, and possibly in vitro and clinical testing is required to thoroughly understand the investigated biomechanical aspects.

CLINICAL SIGNIFICANCE

This study sheds light on the possible biomechanical advantage offered by augmented OLIF implants and provides a theoretical augmentation amount for the SSA construct. Based on the findings, the concept of an SSA device with PMMA augmentation capability is desirable.

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.

The effect of different fixation systems on oblique lumbar interbody fusion under vibration conditions

Despite the fact that lower back pain caused by degenerative lumbar spine pathologies seriously affects the quality of life, however, there is a paucity of research on the biomechanical properties of different auxiliary fixation systems for its primary treatment (oblique lumbar interbody fusion) under vibratory environments. In order to study the effects of different fixation systems of OLIF surgery on the vibration characteristics of the human lumbar spine under whole-body vibration (WBV), a finite element (FE) model of OLIF surgery with five different fixation systems was established by modifying a previously established model of the normal lumbar spine (L1-S1). In this study, a compressive follower load of 500 N and a sinusoidal axial vertical load of ±40 N at the frequency of 5 Hz with a duration of 0.6 s was applied. The results showed that the bilateral pedicle screw fixation model had the highest resistance to cage subsidence and maintenance of disc height under WBV. In contrast, the lateral plate fixation model exerted very high stresses on important tissues, which would be detrimental to the patient's late recovery and reduction of complications. Therefore, this study suggests that drivers and related practitioners who are often in vibrating environments should have bilateral pedicle screws for OLIF surgery, and side plates are not recommended to be used as a separate immobilization system. Additionally, the lateral plate is not recommended to be used as a separate fixation system.

The biomechanical effect of lumbopelvic distance reduction on reconstruction after total sacrectomy: a comparative finite element analysis of four techniques

BACKGROUND CONTEXT

Following total sacrectomy, lumbopelvic reconstruction is essential to restore continuity between the lumbar spine and pelvis. However, to achieve long-term clinical stability, bony fusion between the lumbar spine and the pelvic ring is crucial. Reduction of the lumbopelvic distance can promote successful bony fusion. Although many lumbopelvic reconstruction techniques (LPRTs) have been previously analyzed, the biomechanical effect of lumbopelvic distance reduction (LPDR) has not been investigated yet.

PURPOSE

To evaluate and compare the biomechanical characteristics of four different LPRTs while considering the effect of LPDR.

STUDY DESIGN/SETTING

A comparative finite element (FE) study.

METHODS

The FE models following total sacrectomy were developed to analyze four different LPRTs, with and without LPDR. The closed-loop reconstruction (CLR), the sacral-rod reconstruction (SRR), the four-rod reconstruction (FRR), and the improved compound reconstruction (ICR) techniques were analyzed in flexion, extension, lateral bending, and axial rotation. Lumbopelvic stability was assessed through the shift-down displacement and the relative sagittal rotation of L5, while implant safety was evaluated based on the stress state at the bone-implant interface and within the rods.

RESULTS

Regardless of LPDR, both the shift-down displacement and relative sagittal rotation of L5 consistently ranked the LPRTs as ICR

CONCLUSIONS

LPDR significantly improved both lumbopelvic stability and implant safety in all reconstruction techniques after total sacrectomy. LPDR reduced the shift-down displacement of L5, the relative sagittal rotation of L5, and the stress values at the bone-implant interface. Furthermore, in the ICR and SRR techniques, LPDR decreased the peak stress values within the rods. All four investigated LPRTs demonstrated suitability for lumbopelvic reconstruction, with the ICR technique exhibiting the highest lumbopelvic stiffness.

CLINICAL SIGNIFICANCE

LPDR creates a biomechanically advantageous environment following total sacrectomy; therefore, it has the potential to impact the design of custom-made 3D-printed or traditional LPRTs. However, to confirm the findings of the current FE study, long-term clinical trials are recommended.

Collapse

Key Words

• Biomechanics
• Finite element analysis
• Lumbopelvic distance reduction
• Lumbopelvic fixation
• Lumbopelvic reconstruction
• Lumbopelvic stabilization
• Total sacrectomy

Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Mate Turbucz School of PhD Studies, Semmelweis University, Üllői Str. 26, Budapest, Hungary; In Silico Biomechanics Laboratory, National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary
Agoston Jakab Pokorni School of PhD Studies, Semmelweis University, Üllői Str. 26, Budapest, Hungary; In Silico Biomechanics Laboratory, National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary
Benjamin Hajnal School of PhD Studies, Semmelweis University, Üllői Str. 26, Budapest, Hungary; In Silico Biomechanics Laboratory, National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary
Kristof Koch School of PhD Studies, Semmelweis University, Üllői Str. 26, Budapest, Hungary; National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary
Zsolt Szoverfi National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary
Peter Pal Varga National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary
Aron Lazary National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary; Department of Spine Surgery, Department of Orthopaedics, Semmelweis University, Üllői Str. 78/b, Budapest, Hungary
Peter Endre Eltes In Silico Biomechanics Laboratory, National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary; National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary; Department of Spine Surgery, Department of Orthopaedics, Semmelweis University, Üllői Str. 78/b, Budapest, Hungary.

Collapse

Simulation of cancer progression in bone by a virtual thermal flux with a case study on lumbar vertebrae with multiple myeloma

BACKGROUND

Two main problems examining the mechanism of cancer progression in the tissues using the computational models are lack of enough knowledge on the effective factors for such events in vivo environments and lack of specific parameters in the available computational models to simulate such complicated reactions.

METHODS

In this study, it was tried to simulate the progression of cancerous lesions in the bone tissues by an independent parameter from the anatomical and physiological characteristics of the tissues, so to degrade the orthotropic mechanical properties of the bone tissues, a virtual temperature was determined to be used by a well-known framework for simulation of damages in the composite materials. First, the reliability of the FE model to simulate hyperelastic response in the intervertebral discs (IVDs) and progressive failure in the bony components were verified by simulation of some In-Vitro tests, available in the literature. Then, the progression of the osteolytic damage was simulated in a clinical case with multiple myeloma in the lumbar vertebrae.

RESULTS

The FE model could simulate stress-shielding and diffusion of lesion in the posterior elements of the damaged vertebra which led to spinal stenosis. The load carrying shares associated with the anterior half and the posterior half of the examined vertebral body and the posterior elements were estimated equal to 41 %, 47 % and 12 %, respectively for the intact condition, that changed to 14 %, 16 % and 70 %, when lesion occupied one third of the vertebral body.

CONCLUSION

Correlation of the FE results with the deformation shapes, observed in the MRIs for the clinical case study, indicated appropriateness of the procedure, proposed for simulation of the progressive osteolytic damage in the vertebral segments. The future studies may follow simulation of tumor growth for various metastatic tissues using the method, established here.

Impact of osteoporosis and Cement-Augmented fusion on adjacent spinal levels Post-Fusion Surgery: Patient-Specific finite element analysis

Cement-augmentation is a technique commonly used during posterior lumbar instrumented fusion (PLIF) to reinforce compromised osteoporotic vertebral bone, minimize the risk of loosening screws, enhance stability, and improve overall surgical outcomes. In this study, we introduce a novel segmented vertebral body regional modeling approach to investigate the effects of osteoporosis and cement-augmented lumbar fusion on disc biomechanics at spinal levels adjacent to the fused vertebrae. Using our previously validated personalized-poroelastic-osteoligamentous FE model of the spine, fusion was simulated at L4-L5, and the biomechanics of adjacent levels were studied for 30 patients (non-osteoporotic patients (N = 15), osteoporotic patients (N = 15)). PLIF models, with and without cement-augmentation, were developed and compared after an 8 h-rest period (200 N), following a 16 h-cyclic compressive loading of 500-1000 N (40 and 20 min, respectively). Movement in different directions (flexion/ extension/ lateral bending/ axial rotation) was simulated using 10Nm moment before and after cyclic loading. The material mapping algorithm was validated by comparing the results of voxel-based and parametric models. The FE cement-augmented models, subject to daily activity loading, demonstrated significant differences in disc height loss and fluid loss as compared to non-cemented models. The calculated axial stress and fiber strain values were also significantly higher for these models. This work demonstrates that although osteoporosis does not significantly alter the time-dependent characteristics of adjacent IVDs post-surgery, cement-augmentation increases the risk of adjacent segment disease (ASD) incidence. A holistic understanding of the trade-offs and long-term complex interplay between structural reinforcement modalities, including cement augmentation, and altered biomechanics warrants further investigation.

Inter-Specimen Analysis of Diverse Finite Element Models of the Lumbar Spine

Over the past few decades, there has been a growing popularity in utilizing finite element analysis to study the spine. However, most current studies tend to use one specimen for their models. This research aimed to validate multiple finite element models by comparing them with data from in vivo experiments and other existing finite element studies. Additionally, this study sought to analyze the data based on the gender and age of the specimens. For this study, eight lumbar spine (L2-L5) finite element models were developed. These models were then subjected to finite element analysis to simulate the six fundamental motions. CT scans were obtained from a total of eight individuals, four males and four females, ranging in age from forty-four (44) to seventy-three (73) years old. The CT scans were preprocessed and used to construct finite element models that accurately emulated the motions of flexion, extension, lateral bending, and axial rotation. Preloads and moments were applied to the models to replicate physiological loading conditions. This study focused on analyzing various parameters such as vertebral rotation, facet forces, and intradiscal pressure in all loading directions. The obtained data were then compared with the results of other finite element analyses and in vivo experimental measurements found in the existing literature to ensure their validity. This study successfully validated the intervertebral rotation, intradiscal pressure, and facet force results by comparing them with previous research findings. Notably, this study concluded that gender did not have a significant impact on the results. However, the results did highlight the importance of age as a critical variable when modeling the lumbar spine.

The Impact of Spine Pathology on Posterior Ligamentous Complex Structure and Function

PURPOSE OF REVIEW

Spinal ligament is an important component of the spinal column in mitigating biomechanical stress. Particularly the posterior ligamentous complex, which is composed of the ligamentum flavum, interspinous, and supraspinous ligaments. However, research characterizing the biomechanics and role of ligament health in spinal pathology and clinical context are scarce. This article provides a comprehensive review of the implications of spinal pathology on the structure, function, and biomechanical properties of the posterior ligamentous complex.

RECENT FINDINGS

Current research characterizing biomechanical properties of the posterior ligamentous complex is primarily composed of cadaveric studies and finite element modeling, and more recently incorporating patient-specific anatomy into finite element models. The ultimate goal of current research is to understand the relative contributions of these ligamentous structures in healthy and pathological spine, and whether preserving ligaments may play an important role in spinal surgical techniques. At baseline, posterior ligamentous complex structures account for 30-40% of spinal stability, which is highly dependent on the intrinsic biomechanical properties of each ligament. Biomechanics vary widely with pathology and following rigid surgical fixation techniques and are generally maladaptive. Often secondary to morphological changes in the setting of spinal pathology, but morphological changes in ligament may also serve as a primary pathology. Biomechanical maladaptations of the spinal ligament adversely influence overall spinal column integrity and ultimately predispose to increased risk for surgical failure and poor clinical outcomes. Future research is needed, particularly in living subjects, to better characterize adaptations in ligaments that can provide targets for improved treatment of spinal pathology.

Comparison of anterior column reconstruction techniques after en bloc spondylectomy: a finite element study

Total en bloc spondylectomy (TES) effectively treats spinal tumors. The surgery requires a vertebral body replacement (VBR), for which several solutions were developed, whereas the biomechanical differences between these devices still need to be completely understood. This study aimed to compare a femur graft, a polyetheretherketone implant (PEEK-IMP-C), a titan mesh cage (MESH-C), and a polymethylmethacrylate replacement (PMMA-C) using a finite element model of the lumbar spine after a TES of L3. Several biomechanical parameters (rotational stiffness, segmental range of motion (ROM), and von Mises stress) were assessed to compare the VBRs. All models provided adequate initial stability by increasing the rotational stiffness and decreasing the ROM between L2 and L4. The PMMA-C had the highest stiffness for flexion-extension, lateral bending, and axial rotation (215%, 216%, and 170% of intact model), and it had the lowest segmental ROM in the instrumented segment (0.2°, 0.5°, and 0.7°, respectively). Maximum endplate stress was similar for PMMA-C and PEEK-IMP-C but lower for both compared to MESH-C across all loading directions. These results suggest that PMMA-C had similar or better primary spinal stability than other VBRs, which may be related to the larger contact surface and the potential to adapt to the patient's anatomy.

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.

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.

Effect of graded posterior element and ligament removal on annulus stress and segmental stability in lumbar spine stenosis: a finite element analysis study

The study aimed to investigate the impact of posterior element and ligament removal on the maximum von Mises stress, and maximum shear stress of the eight-layer annulus for treating stenosis at the L3-L4 and L4-L5 levels in the lumbar spine. Previous studies have indicated that laminectomy alone can result in segmental instability unless fusion is performed. However, no direct correlations have been established regarding the impact of posterior and ligament removal. To address this gap, four models were developed: Model 1 represented the intact L2-L5 model, while model 2 involved a unilateral laminotomy involving the removal of a section of the L4 inferior lamina and 50% of the ligament flavum between L4 and L5. Model 3 consisted of a complete laminectomy, which included the removal of the spinous process and lamina of L4, as well as the relevant connecting ligaments between L3-L4 and L4-L5 (ligament flavum, interspinous ligament, supraspinous ligament). In the fourth model, a complete laminectomy with 50% facetectomy was conducted. This involved the same removals as in model 3, along with a 50% removal of the inferior/superior facets of L4 and a 50% removal of the facet capsular ligaments between L3-L4 and L4-L5. The results indicated a significant change in the range of motion (ROM) at the L3-L4 and L4-L5 levels during flexion and torque situations, but no significant change during extension and bending simulation. The ROM increased by 10% from model 1 and 2 to model 3, and by 20% to model 4 during flexion simulation. The maximum shear stress and maximum von-Mises stress of the annulus and nucleus at the L3-L4 levels exhibited the greatest increase during flexion. In all eight layers of the annulus, there was an observed increase in both the maximum shear stress and maximum von-Mises stress from model 1&2 to model 3 and model 4, with the highest rate of increase noted in layers 7&8. These findings suggest that graded posterior element and ligament removal have a notable impact on stress distribution and range of motion in the lumbar spine, particularly during flexion.

Influence of Simulated State of Disc Degeneration and Axial Stiffness of Coupler in a Hybrid Performance Stabilisation System on the Biomechanics of a Spine Segment Model

Spinal fusion surgery leads to the restriction of mobility in the vertebral segments postoperatively, thereby causing stress to rise at the adjacent levels, resulting in early degeneration and a high risk of adjacent vertebral fractures. Thus, to address this issue, non-fusion surgery applies some pedicle screw-based dynamic stabilisation systems to provide stability and micromotion, thereby reducing stress in the fusion segments. Among these systems, the hybrid performance stabilisation system (HPSS) combines a rigid rod, transfer screw, and coupler design to offer a semi-rigid fixation method that preserves some mobility near the fusion site and reduces the adjacent segment compensatory effects. However, further research and confirmation are needed regarding the biomechanical effects of the dynamic coupler stiffness of the HPSS on the intrinsic degenerated adjacent segment. Therefore, this study utilised the finite element method to investigate the impact of the coupler stiffness of the HPSS on the mobility of the lumbar vertebral segments and the stress distribution in the intervertebral discs under flexion, extension, and lateral bending, as well as the clinical applicability of the HPSS on the discs with intrinsic moderate and severe degeneration at the adjacent level. The analytical results indicated that, regardless of the degree of disc degeneration, the use of a dynamic coupler stiffness of 57 N/mm in the HPSS may reduce the stress concentrations at the adjacent levels. However, for severely degenerated discs, the postoperative stress on the adjacent segments with the HPSS was still higher compared with that of the discs with moderate degeneration. We conclude that, when the discs had moderate degeneration, increasing the coupler stiffness led to a decrease in disc mobility. In the case of severe disc degeneration, the effect on disc mobility by coupler stiffness was less pronounced. Increasing the coupler stiffness ked to higher stress on intervertebral discs with moderate degeneration, while its effect on stress was less pronounced for discs with severe degeneration. It is recommended that patients with severe degeneration who undergo spinal dynamic stabilisation should remain mindful of the risk of accelerated adjacent segment degeneration.

A comparison of intervertebral ligament properties utilized in a thoracic spine functional unit through kinematic evaluation

Ligament properties in the literature are variable, yet scarce, but needed to calibrate computational models for spine clinical research applications. A comparison of ligament stiffness properties and their effect on the kinematic behavior of a thoracic functional spinal unit (FSU) is examined in this paper. Six unique ligament property sets were utilized within a volumetric T7-T8 finite element (FE) model developed using computer-aided design (CAD) spinal geometry. A 7.5 Nm moment was applied along three anatomical planes both with and without costovertebral (CV) joints present. Range of Motion (RoM) was assessed for each property set and compared to published experimental data. Intact and serial ligament removal procedures were implemented in accordance with experimental protocol. The variance in both kinematic behavior and comparability with experimental data among property sets emphasizes the role nonlinear characterization plays in determining proper kinematic behavior in spinal FE models. Additionally, a decrease in RoM variation among property sets was exhibited when the model setup incorporated the CV joint. With proper assessment of the source and size of each ligament, the material properties considered here could be expanded and justified for implementation into thoracic spine clinical studies.

Development and validation of lumbar spine finite element model

The functional biomechanics of the lumbar spine have been better understood by finite element method (FEM) simulations. However, there are still areas where the behavior of soft tissues can be better modeled or described in a different way. The purpose of this research is to develop and validate a lumbar spine section intended for biomechanical research. A FE model of the 50th percentile adult male (AM) Total Human Model for Safety (THUMS) v6.1 was used to implement the modifications. The main modifications were to apply orthotropic material properties and nonlinear stress-strain behavior for ligaments, hyperelastic material properties for annulus fibrosus and nucleus pulposus, and the specific content of collagenous fibers in the annulus fibrosus ground substance. Additionally, a separation of the nucleus pulposus from surrounding bones and tissues was implemented. The FE model was subjected to different loading modes, in which intervertebral rotations and disc pressures were calculated. Loading modes contained different forces and moments acting on the lumbar section: axial forces (compression and tension), shear forces, pure moments, and combined loading modes of axial forces and pure moments. The obtained ranges of motion from the modified numerical model agreed with experimental data for all loading modes. Moreover, intradiscal pressure validation for the modified model presented a good agreement with the data available from the literature. This study demonstrated the modifications of the THUMS v6.1 model and validated the obtained numerical results with existing literature in the sub-injurious range. By applying the proposed changes, it is possible to better model the behavior of the human lumbar section under various loads and moments.

Development, calibration and validation of a comprehensive customizable lumbar spine FE model for simulating fusion constructs

Instrumentation alters the biomechanics of the spine, and therefore prediction of all output quantities that have critical influence post-surgically is significant for engineering models to aid in clinical predictions. Geometrical morphological finite element models can bring down the development time and cost of custom intact and instrumented models and thus aids in the better inference of biomechanics of surgical instrumentation on patient-specific diseased spine segments. A comprehensive hexahedral morphological lumbosacral finite element model is developed in this work to predict the range of motions, disc pressures, and facet contact forces of the intact and instrumented spine. Facet contact forces are needed to predict the impact of fusion surgeries on adjacent facet contacts in bending, axial rotation, and extension motions. Extensive validation in major physiological loading regimes of the pure moment, pure compression, and combined loading is undertaken. In vitro, experimental corridor results from six different studies reported in the literature are compared and the generated model had statistically significant comparable values with these studies. Flexion, extension and bending moment rotation curves of all segments of the developed model were favourable and within two separately established experimental corridor windows as well as recent simulation results. Axial torque moment rotation curves were comparable to in vitro results for four out of five lumbar functional units. The facet contact force results also agreed with in vitro experimental results. The current model is also computationally efficient with respect to contemporary models since it uses significantly smaller number of elements without losing the accuracy in terms of response prediction. This model can further be used for predicting the impact of different instrumentation techniques on the lumbar vertebral column.

Biomechanical Effects of a Novel Pedicle Screw W-Type Rod Fixation for Lumbar Spondylolysis: A Finite Element Analysis

Lumbar spondylolysis involves anatomical defects of the pars interarticularis, which causes instability during motion. The instability can be addressed through instrumentation with posterolateral fusion (PLF). We developed a novel pedicle screw W-type rod fixation system and evaluated its biomechanical effects in comparison with PLF and Dynesys stabilization for lumbar spondylolysis via finite element (FE) analysis. A validated lumbar spine model was built using ANSYS 14.5 software. Five FE models were established simulating the intact L1-L5 lumbar spine (INT), bilateral pars defect (Bipars), bilateral pars defect with PLF (Bipars_PLF), Dynesys stabilization (Bipars_Dyn), and W-type rod fixation (Bipars_Wtyp). The range of motion (ROM) of the affected segment, the disc stress (DS), and the facet contact force (FCF) of the cranial segment were compared. In the Bipars model, ROM increased in extension and rotation. Compared with the INT model, Bipars_PLF and Bipars_Dyn exhibited remarkably lower ROMs for the affected segment and imposed greater DS and FCF in the cranial segment. Bipars_Wtyp preserved more ROM and generated lower stress at the cranial segment than Bipars_PLF or Bipars_Dyn. The injury model indicates that this novel pedicle screw W-type rod for spondylolysis fixation could return ROM, DS, and FCF to levels similar to preinjury.

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.

Recent Advances in Coupled MBS and FEM Models of the Spine—A Review

How back pain is related to intervertebral disc degeneration, spinal loading or sports-related overuse remains an unanswered question of biomechanics. Coupled MBS and FEM simulations can provide a holistic view of the spine by considering both the overall kinematics and kinetics of the spine and the inner stress distribution of flexible components. We reviewed studies that included MBS and FEM co-simulations of the spine. Thereby, we classified the studies into unidirectional and bidirectional co-simulation, according to their data exchange methods. Several studies have demonstrated that using unidirectional co-simulation models provides useful insights into spinal biomechanics, although synchronizing the two distinct models remains a key challenge, often requiring extensive manual intervention. The use of a bidirectional co-simulation features an iterative, automated process with a constant data exchange between integrated subsystems. It reduces manual corrections of vertebra positions or reaction forces and enables detailed modeling of dynamic load cases. Bidirectional co-simulations are thus a promising new research approach for improved spine modeling, as a main challenge in spinal biomechanics is the nonlinear deformation of the intervertebral discs. Future studies will likely include the automated implementation of patient-specific bidirectional co-simulation models using hyper- or poroelastic intervertebral disc FEM models and muscle forces examined by an optimization algorithm in MBS. Applications range from clinical diagnosis to biomechanical analysis of overload situations in sports and injury prediction.

Shorter screw lengths in dynamic Dynesys fixation have less screw loosening: From clinical investigation to finite-element analysis

BACKGROUND

The dynamic Dynesys Stabilization System preserves lumbar mobility at instrumented levels. This study investigated the effect of screw length on screw loosening (SL) after dynamic Dynesys fixation and screw displacement during lumbar motion, using clinical investigation and finite-element (FE) analysis.

METHODS

Clinical data of 50 patients with degenerative spondylolisthesis treated with decompression and Dynesys fixation in 2011 were analyzed retrospectively. Horizontal sliding displacement and vertical displacement of screw tips at L4 were analyzed postoperatively using displacement-controlled FE analysis at the L4-L5 level with screw lengths 45 (long screw), 36 (median screw), and 27 (short screw), and 6.4 mm in diameter, under flexion, extension, lateral bending, and rotation.

RESULTS

In 13 patients (13/50, 26%), 40 screws (40/266, 15%) were loose at mean follow-up of 101.3 ± 4.4 months. Radiographic SL at 35, 40, 45, and 50 mm were 7.7%, 10.7%, 12.1%, and 37.5%, respectively, regardless of the fixation level ( p = 0.009). FE analysis revealed that the long screw model with corresponding longer lever arm had maximal horizontal sliding displacement under all directions and maximal vertical displacement, except for lateral bending.

CONCLUSION

Shorter screws in Dynesys fixation may help avoid dynamic SL. Clinically, 50 mm screws showed the greatest SL and median screw screws demonstrated the least displacement biomechanically.

Anatomical parameters alter the biomechanical responses of adjacent segments following lumbar fusion surgery: Personalized poroelastic finite element modelling investigations

Introduction: While the short-term post-operative outcome of lumbar fusion is satisfying for most patients, adjacent segment disease (ASD) can be prevalent in long-term clinical observations. It might be valuable to investigate if inherent geometrical differences among patients can significantly alter the biomechanics of adjacent levels post-surgery. This study aimed to utilize a validated geometrically personalized poroelastic finite element (FE) modeling technique to evaluate the alteration of biomechanical response in adjacent segments post-fusion. Methods: Thirty patients were categorized for evaluation in this study into two distinct groups [i.e., 1) non-ASD and 2) ASD patients] based on other long-term clinical follow-up investigations. To evaluate the time-dependent responses of the models subjected to cyclic loading, a daily cyclic loading scenario was applied to the FE models. Different rotational movements in different planes were superimposed using a 10 Nm moment after daily loading to compare the rotational motions with those at the beginning of cyclic loading. The biomechanical responses of the lumbosacral FE spine models in both groups were analyzed and compared before and after daily loading. Results: The achieved comparative errors between the FE results and clinical images were on average below 20% and 25% for pre-op and post-op models, respectively, which confirms the applicability of this predictive algorithm for rough pre-planning estimations. The results showed that the disc height loss and fluid loss were increased for the adjacent discs in post-op models after 16 h of cyclic loading. In addition, significant differences in disc height loss and fluid loss were observed between the patients who were in the non-ASD and ASD groups. Similarly, the increased stress and fiber strain in the annulus fibrosus (AF) was higher in the adjacent level of post-op models. However, the calculated stress and fiber strain values were significantly higher for patients with ASD. Discussion: Evaluating the biomechanical response of pre-op and post-op modeling in the non-ASD and ASD groups showed that the inherent geometric differences among patients cause significant variations in the estimated mechanical response. In conclusion, the results of the current study highlighted the effect of geometrical parameters (which may refer to the anatomical conditions or the induced modifications regarding surgical techniques) on time-dependent responses of lumbar spine biomechanics.

Implications of range of motion requirements for the laxity of ligaments in a lumbar finite element model

Finite element models of the lumbar spine often adopt ligament properties from tensile tests without accounting for possible differences between testing and in situ initial ligament length. Such differences could result in laxities or preloads at the beginning of a simulation that would affect the ligament forces, tangent stiffness, and the posture at which they fail. In vivo and in vitro human experimental data reported laxities or preloads. However, laxities or preloads, which could also result from postural differences, are often neglected in simulation studies. This study proposes a numerical methodology to identify ranges of ligament laxities or preloads compatible with the selected tensile ligament properties, the model, and the range of motion (RoM) the model aims to simulate. The approach assumes that ligaments should remain in a safe elongation range for the complete RoM, and that each ligament should play a significant mechanical role in at least one load case. The methodology was applied to the functional spinal unit (FSU) models using the RoM from healthy subjects and ligament properties from the literature. Without laxity, some ligaments reached their elongation at failure within the RoM. Laxity ranges varied considerably (from -9.2 mm preload to 10.7 mm laxity) and flexion was the most critical load case to determine them. Their effect on the mobility response was also assessed. The effect on the mobility of a FSU was also assessed. While the proposed method cannot determine an exact laxity value, it is simple and it can be applied to any model to identify a plausible range of ligament initial length.

Comparison of the biomechanical effects of lumbar disc degeneration on normal patients and osteoporotic patients: A finite element analysis

BACKGROUND

Some older patients who suffered from both conditions (disc degeneration and osteoporosis) have higher surgical risks and longer postoperative recovery times. Understanding the relation between disc degeneration and osteoporosis is fundamental to know the mechanisms of orthopedic disorders and improve clinical treatment. However, there is a lack of finite element (FE) studies to predict the combined effects of disc degeneration and osteoporosis. So the aim of the present study is to explore the differences of biomechanical effects of lumbar disc degeneration on normal patients and osteoporotic patients.

METHODS

A normal lumbar spine finite element model (FEM) was developed based on the geometric information of a healthy male subject (age 35 years; height 178 cm; weight 65 kg). This normal lumbar spine FEM was modified to build three lumbar spine degeneration models simulating mild, moderate and severe grades of disc degeneration at the L4-L5 segment. Then the degenerative lumbar spine models for osteoporotic patients were constructed on the basis of the above-mentioned degeneration models. Firstly, the normal model (flexion: 8 Nm; extension: 6 Nm; lateral bending: 6 Nm; torsion: 4 Nm) and degenerative models (10 Nm) were calibrated under pure moment load, respectively. Secondly, under a 400 N follower load, the 7.5 Nm moments of different directions were applied on all models to simulate different motion postures. Finally, under the above loading conditions, we calculated and analyzed the range of motion (ROM), Mises stress in cortical (MSC1), Mises stress in endplate (MSE), Mises stress in cancellous (MSC2), and Mises stress in post (MSP).

RESULTS

Compared with disc degeneration patients without osteoporosis, the ROM, MSC1, and MSE of osteoporosis patients with various disc degeneration decreased in all postures, while the MSC2 and MSP increased. With increase in the degree of disc degeneration, the reduction proportions of ROM and MSE in osteoporotic patients gradually increased, while the reduction percentages in MSC1 of osteoporotic patients gradually decreased. The increase percentages of MSC2 in osteoporotic patients gradually increased. Given the progressive changes of disc degeneration, the changes in MSP in osteoporosis patients were uneven.

CONCLUSION

In summary, the effect of disc degeneration on flexibility in the two kinds of patients (osteoporosis and non-osteoporosis patients) was nearly same. By comparing the remaining biomechanical parameters (MSC1, MSE, MSC2, and MSP), we found that degenerated intervertebral discs caused changes in loading patterns of osteoporosis patients. Disc degeneration reduced the Mises stress in the cortical and endplate, which increased the Mises stress in the cancellous and post. That is to say, in order to cope with the changes in bone stresses caused by disc degeneration and osteoporosis, clinicians should be more careful in choosing the surgical option for osteoporotic patients with disc degeneration.

Development and model form assessment of an automatic subject-specific vertebra reconstruction method

BACKGROUND

Current spine models for analog bench models, surgical navigation and training platforms are conventionally based on 3D models from anatomical human body polygon database or from time-consuming manual-labelled data. This work proposed a workflow of quick and accurate subject-specific vertebra reconstruction method and quantified the reconstructed model accuracy and model form errors.

METHODS

Four different neural networks were customized for vertebra segmentation. To validate the workflow in clinical applications, an excised human lumbar vertebra was scanned via CT and reconstructed into 3D CAD models using four refined networks. A reverse engineering solution was proposed to obtain the high-precision geometry of the excised vertebra as gold standard. The 3D model evaluation metrics and a finite element analysis (FEA) method were designed to reflect the model accuracy and model form errors.

RESULTS

The automatic segmentation networks achieved the best Dice score of 94.20% in validation datasets. The accuracy of reconstructed models was quantified with the best 3D Dice index of 92.80%, 3D IoU of 86.56%, Hausdorff distance of 1.60 mm, and the heatmaps and histograms were used for error visualization. The FEA results showed the impact of different geometries and reflected partial surface accuracy of the reconstructed vertebra under biomechanical loads with the closest percentage error of 4.2710% compared to the gold standard model.

CONCLUSIONS

In this work, a workflow of automatic subject-specific vertebra reconstruction method was proposed while the errors in geometry and FEA were quantified. Such errors should be considered when leveraging subject-specific modelling towards the development and improvement of treatments.

Adjacent segments biomechanics following lumbar fusion surgery: a musculoskeletal finite element model study

PURPOSE

This study exploits a novel musculoskeletal finite element (MS-FE) spine model to evaluate the post-fusion (L4-L5) alterations in adjacent segment kinetics.

METHODS

Unlike the existing MS models with idealized representation of spinal joints, this model predicts stress/strain distributions in all passive tissues while organically coupled to a MS model. This generic (in terms of musculature and material properties) model uses population-based in vivo vertebral sagittal rotations, gravity loads, and an optimization algorithm to calculate muscle forces. Simulations represent individuals with an intact L4-L5, a preoperative severely degenerated L4-L5 (by reducing the disc height by ~ 60% and removing the nucleus incompressibility), and a postoperative fused L4-L5 segment with either a fixed or an altered lumbopelvic rhythm with respect to the intact condition (based on clinical observations). Changes in spine kinematics and back muscle cross-sectional areas (due to intraoperative injuries) are considered based on in vivo data while simulating three activities in upright/flexed postures.

RESULTS

Postoperative changes in some adjacent segment kinetics were found considerable (i.e., larger than 25%) that depended on the postoperative lumbopelvic kinematics and preoperative L4-L5 disc condition. Postoperative alterations in adjacent disc shear, facet/ligament forces, and annulus stresses/strains were greater (> 25%) than those found in intradiscal pressure and compression (< 25%). Kinetics of the lower (L5-S1) and upper (L3-L4) adjacent segments were altered to different degrees.

CONCLUSION

Alterations in segmental rotations mainly affected adjacent disc shear forces, facet/ligament forces, and annulus/collagen fibers stresses/strains. An altered lumbopelvic rhythm (increased pelvis rotation) tends to mitigate some of these surgically induced changes.

Comparison of the susceptibility to implant failure in the lateral, posterior, and transforaminal lumbar interbody fusion: A finite element analysis

OBJECTIVE

There are several techniques for lumbar interbody fusion, and implant failure following lumbar interbody fusion can be troublesome. This study aimed to compare the stress in posterior implant and peri-screw vertebral bodies among lateral lumbar interbody fusion (LLIF), posterior lumbar interbody fusion (PLIF), and transforaminal lumbar interbody fusion (TLIF) and to select the technique that is least likely to cause implant failure.

METHODS

We created an intact L3-L5 model and simulated the LLIF, PLIF, and TLIF techniques at L4-L5 using finite element methods. All models at the lower portion of L5 were fixed and imposed a preload of 400 N and a moment of 7.5 Nm on the upper portion of L3 to simulate flexion, extension, lateral bending, and axial rotation. We investigated the peak stresses and stress concentration in the posterior implant and peri-screw vertebral bodies for the LLIF, PLIF, and TLIF techniques.

RESULTS

The extension, flexion, bending, and rotation peak stresses and stress concentration in the posterior implant, and the peri-screw vertebral bodies, were the lowest in LLIF, followed by PLIF and TLIF, respectively.

CONCLUSIONS

It was found that implant failure was least likely to occur in LLIF, followed by PLIF and TLIF, respectively. Hence, surgeons should be aware of these factors when selecting an appropriate surgical technique and be careful for implant failure during postoperative follow-up.

Biomechanical study of oblique lumbar interbody fusion (OLIF) augmented with different types of instrumentation: a finite element analysis

Background

To explore the biomechanical differences in oblique lumbar interbody fusion (OLIF) augmented by different types of instrumentation.

Methods

A three-dimensional nonlinear finite element (FE) model of an intact L3-S1 lumbar spine was built and validated. The intact model was modified to develop five OLIF surgery models (Stand-alone OLIF; OLIF with lateral plate fixation [OLIF + LPF]; OLIF with unilateral pedicle screws fixation [OLIF + UPSF]; OLIF with bilateral pedicle screws fixation [OLIF + BPSF]; OLIF with translaminar facet joint fixation + unilateral pedicle screws fixation [OLIF + TFJF + UPSF]) in which the surgical segment was L4–L5. Under a follower load of 500 N, a 7.5-Nm moment was applied to all lumbar spine models to calculate the range of motion (ROM), equivalent stress peak of fixation instruments (ESPFI), equivalent stress peak of cage (ESPC), equivalent stress peak of cortical endplate (ESPCE), and equivalent stress average value of cancellous bone (ESAVCB).

Results

Compared with the intact model, the ROM of the L4–L5 segment in each OLIF surgery model decreased by > 80%. The ROM values of adjacent segments were not significantly different. The ESPFI, ESPC, and ESPCE values of the OLIF + BPSF model were smaller than those of the other OLIF surgery models. The ESAVCB value of the normal lumbar model was less than the ESAVCB values of all OLIF surgical models. In most postures, the ESPFI, ESPCE, and ESAVCB values of the OLIF + LPF model were the largest. The ESPC was higher in the Stand-alone OLIF model than in the other OLIF models. The stresses of several important components of the OLIF + UPSF and OLIF + TFJF + UPSF models were between those of the OLIF + LPF and OLIF + BPSF models.

Conclusions

Our biomechanical FE analysis indicated the greater ability of OLIF + BPSF to retain lumbar stability, resist cage subsidence, and maintain disc height. Therefore, in the augmentation of OLIF, bilateral pedicle screws fixation may be the best approach.

Biomechanical Analysis of Different Internal Fixation Combined with Different Bone Grafting for Unstable Thoracolumbar Fractures in the Elderly

This research was developed to accurately evaluate the unstable fractures of thoracolumbar before and after surgery and discuss the treatment timing and methods. Three-dimensional (3D) finite element method was adopted to construct the T₁₂-L₅ segment model of human body. The efficiency of percutaneous kyphoplasty (PKP) and percutaneous vertebroplasty (PVP), two commonly used internal fixation procedures, was retrospectively compared. A total of 150 patients with chest fracture who received PKP or PVP surgery in our hospital, and 104 patients with the same symptoms who received conservative treatment were collected and randomly rolled into PVP group (75 cases), PKP group (75 cases), and control group (104 cases). Visual analog scale (VAS) score and Oswestry disability index (ODI) of patients were collected before and after surgery and 2, 12, and 24 months after surgery. Then, the anterior and central height of the patient's cone and the kyphosis angle were calculated by X-ray. Lumbar minimally invasive fusion system and lumbar pedicle screw rod system were established by computer-aided design (CAD), and the biomechanical characteristics were analyzed. The results showed that there was no substantial difference in VAS score and ODI score between PKP and PVP (P > 0.05), but they were higher than those of the control group (P < 0.05). The anterior edge and middle height of vertebra in the two groups were higher than those in control group (P < 0.05), and the increase in PKP group was more substantial (P < 0.05). The kyphosis of the two groups was smaller than that of the control group (P < 0.05), and the decrease of the kyphosis of the PKP group was more substantial (P < 0.05). In summary, the thoracolumbar segment model established by 3D finite element method was an effective model, and it was verified on patients that both PKP and PVP could achieve relatively satisfactory efficacy. The implantation of the new internal fixation system had no obvious effect on the lumbar movement. This work provided a novel idea and method for the treatment of senile thoracolumbar unstable fracture, as well as experimental data of biomechanics for the operation of senile unstable fracture.

Biomechanical changes at the adjacent segments induced by a lordotic porous interbody fusion cage

Biomechanical changes at the adjacent segments after interbody fusion are common instigators of adjacent segment degeneration (ASD). This study aims to investigate how the presence of a lordotic porous cage affects the biomechanical performance of the adjacent segments. A finite element model (FEM) of a lumbar spine implanted with a lordotic cage at L3-L4 was validated by in-vitro testing. The stress distribution on the cage and range of motion (ROM) of L3-L4 were used to assess the stability of the implant. Three angles of cage (0° = non-restoration, 7° = normal restoration and 11° = over-restoration) were modelled with different porosities (0%, 30% and 60%) and evaluated in the motions of flexion, extension, lateral bending and rotation. The ROM, intervertebral disc pressure (IDP) and facet joint force (FJF) were used to evaluate biomechanical changes at the adjacent segments in each model. The results indicated that porous cages produced more uniform stress distribution, but cage porosity did not influence the ROM, IDP and FJF at L2-L3 and L4-L5. Increasing the cage lordotic angle acted to decrease the ROM and IDP, and increase the FJF of L4-L5, but did not alter the ROM of L2-L3. In conclusion, changes in ROM, IDP and FJF at the adjacent segments were mainly influenced by the lordotic angle of the cage and not by the porosity. A larger angle of lordotic cage was shown to reduce the ROM and IDP, and increase the FJF of the lower segment (L4-L5), but had little effect on the ROM of the upper segment (L2-L3).

A biomechanical study of ligament tethers strengthening for the prevention of proximal junctional kyphosis after posterior long-segment spinal fusion

BACKGROUND

Proximal junctional kyphosis is a known complication of posterior long-segment thoracolumbar fusion. Here, the biomechanical effectiveness of ligament tethers strengthening and vertebral body augmentation, in proximal junctional kyphosis prevention was explored using the finite element analysis.

METHODS

Based on a validated model of T1-L5 with the pedicle screw system instrumented T8-L5, strengthening models with different strategies were created to assess the range of motion in proximal vertebrae, vertebrae stress, pedicle screw stress, and pressure on intervertebral discs during extension, flexion, lateral bending, and axial rotation motions. Strengthening strategies included two- and three-level posterior ligament tethers (TE-T7-T9 and TE-T6-T9), and tethers with T7 &T8 vertebral body augmentation (TECE-T7-T9 and TECE-T6-T9).

FINDINGS

Compared to the spinal fusion model, the ligament tethers strengthening significantly reduced the flexion-extension range of motion difference among the proximal vertebrae. During the flexion-extension motion, the T8 vertebra stresses in the TE-T7-T9, TE-T6-T9, TECE-T7-T9, and TECE-T6-T9 models were distinctively reduced, the values decreased by 26.8%, 28.3%, 28.8%, and 9.6%, respectively, during flexion, and by 21.9%, 35.2%, 23%, and 18.6%, respectively, during extension. In the strengthening models, the maximum stresses on the T7/T8 intervertebral disc in the TE-T6-T9 model were reduced by 13.8% during flexion and by 14.7% during extension.

INTERPRETATION

Based on our results, the ligament strengthening configuration of the three-level posterior tethers produced a more gradual transition in range of motion, vertebrae stresses, and intervertebral discs stress between the fused and non-fused segments, especially during flexion-extension, which may significantly decrease the proximal junctional kyphosis biomechanical risk.

Biphasic Properties of PVAH (Polyvinyl Alcohol Hydrogel) Reflecting Biomechanical Behavior of the Nucleus Pulposus of the Human Intervertebral Disc

PVAH is a mixture of solid and fluid, but its mechanical behavior has usually been described using solid material models. The purpose of this study was to obtain material properties that can reflect the mechanical behavior of polyvinyl alcohol hydrogel (PVAH) using finite element analysis, a biphasic continuum model, and to optimize the composition ratio of PVAH to replace the nucleus pulposus (NP) of the human intervertebral disc. Six types of PVAH specimens (3, 5, 7, 10, 15, 20 wt%) were prepared, then unconfined compression experiments were performed to acquire their material properties using the Holmes–Mow biphasic model. With an increasing weight percentage of PVA in PVAH, the Young’s modulus increased while the permeability parameter decreased. The Young’s modulus and permeability parameter were similar to those of the NP at 15 wt% and 20 wt%. The range of motion, facet joint force, and NP pressures measured from dynamic motional analysis of the lumbar segments with the NP model also exhibited similar values to those with 15~20 wt% PVAH models. Considering the structural stability and pain of the lumbar segments, it appears that 20 wt% PVAH is most suitable for replacing the NP.

Development of a finite element lumbar spine model to predict intervertebral disc herniation risk factors

The objective of the current study was to develop a lumbar motion segment FE model that predicts disc herniation risk. The posterolateral nucleus extrusion force and disc pressure increased as the amount of flexion and magnitude of compression was increased in all loading scenarios. The nucleus extrusion force and posterior stress in the annulus both increased when exposed to a combination of compression and flexion. Results of the current study confirmed the authors hypothesis that the model would accurately predict herniation risk when exposed to a biomechanical environment known to cause herniations.

Stability Evaluation of Different Oblique Lumbar Interbody Fusion Constructs in Normal and Osteoporotic Condition - A Finite Element Based Study

Introduction: In developed countries, the age structure of the population is currently undergoing an upward shift, resulting a decrease in general bone quality and surgical durability. Over the past decade, oblique lumbar interbody fusion (OLIF) has been globally accepted as a minimally invasive surgical technique. There are several stabilization options available for OLIF cage fixation such as self-anchored stand-alone (SSA), lateral plate-screw (LPS), and bilateral pedicle screw (BPS) systems. The constructs’ stability are crucial for the immediate and long-term success of the surgery. The aim of this study is to investigate the biomechanical effect of the aforementioned constructs, using finite element analysis with different bone qualities (osteoporotic and normal).

Method: A bi-segmental (L2–L4) finite element (FE) model was created, using a CT scan of a 24-year-old healthy male. After the FE model validation, CAD geometries of the implants were inserted into the L3–L4 motion segment during a virtual surgery. For the simulations, a 150 N follower load was applied on the models, then 10 Nm of torque was used in six general directions (flexion, extension, right/left bending, and right/left rotation), with different bone material properties.

Results: The smallest segmental (L3–L4) ROM (range of motion) was observed in the BPS system, except for right bending. Osteoporosis increased ROMs in all constructs, especially in the LPS system (right bending increase: 140.26%). Osteoporosis also increased the caudal displacement of the implanted cage in all models (healthy bone: 0.06 ± 0.03 mm, osteoporosis: 0.106 ± 0.07 mm), particularly with right bending, where the displacement doubled in SSA and LPS constructs. The displacement of the screws inside the L4 vertebra increased by 59% on average (59.33 ± 21.53%) due to osteoporosis (100% in LPS, rotation). BPS-L4 screw displacements were the least affected by osteoporosis.

Conclusions: The investigated constructs provide different levels of stability to the spine depending on the quality of the bone, which can affect the outcome of the surgery. In our model, the BPS system was found to be the most stable construct in osteoporosis. The presented model, after further development, has the potential to help the surgeon in planning a particular spinal surgery by adjusting the stabilization type to the patient’s bone quality.

Biomechanics of Lumbar Spine Injury in Road Barrier Collision-Finite Element Study

Literature and field data from CIREN database have shown that lumbar spine injuries occur during car crashes. There are multiple hypotheses regarding how they occur; however, there is no biomechanical explanation for these injuries during collisions with road safety barriers (RSBs). Therefore, the objective of this study was to investigate the mechanics of vertebral fractures during car collisions with concrete RSBs. The finite element method was used for the numerical simulations. The global model of the car collision with the concrete RSB was created. The lumbar spine kinematics were extracted from the global simulation and then applied as boundary conditions to the detailed lumbar spine model. The results showed that during the collision, the occupant was elevated, and then dropped during the vehicle landing. This resulted in axial compression forces 2.6 kN with flexion bending moments 34.7 and 37.8 Nm in the L2 and L3 vertebrae. It was shown that the bending moment is the result of the longitudinal force on the eccentricity. The lumbar spine index for the L1-L5 section was 2.80, thus indicating a lumbar spine fracture. The minimum principal strain criterion of 7.4% and damage variable indicated L2 and L3 vertebrae and the inferior part of L1, as those potentially prone to fracture. This study found that lumbar spine fractures could occur as a consequence of vehicle landing during a collision with a concrete RSB mostly affecting the L1-L3 lumbar spine section. The fracture was caused by a combination of axial forces and flexion bending moments.

A biomechanical analysis of anterior cervical discectomy and fusion alone or combined cervical fixations in treating compression-extension injury with unilateral facet joint fracture: a finite element study

Background

Compression-extension injury with unilateral facet joint fracture is one of the most devastating injuries of subaxial cervical spine. However, it is not yet clear which fixation technique represents the optimal choice in surgical management. This study aims to assess the construct stability at the operative level (C4/C5 cervical spine) following anterior cervical discectomy and fusion (ACDF) alone and combined fixation techniques (posterior-anterior fixations).

Methods

A previously validated three-dimensional C2-T1 finite element model were modified to simulate surgical procedures via the anterior-only approach (ACDF) and combined cervical approach [(transarticular screw, lateral mass screw, unilateral pedicle screw, bilateral pedicle screw) + ACDF, respectively] for treating compression-extension injury with unilateral facet joint fracture at C4/C5 level. Construct stability (range of rotation, axial compression displacement and anterior shear displacement) at the operative level was comparatively analyzed.

Results

In comparison with combined fixation techniques, a wider range of motion and a higher maximum von Mises stress was found in single ACDF. There was no obvious difference in range of motion among transarticular screw and other posterior fixations in the presence of anterior fixation. In addition, the screws inserted by transarticular screw technique had high stress concentration at the middle part of the screw but much less than 500 MPa under different conditions. Furthermore, the variability of von Mises stress in the transarticular screw fixation device was significantly lower than ACDF but no obvious difference compared with other posterior fixations.

Conclusions

Of the five fixation techniques, ACDF has proven poor stability and high structural stress. Compared with lateral and pedicle screw, transarticular screw technique was not worse biomechanically and less technically demanding to acquire in clinical practice. Therefore, our study suggested that combined fixation technique (transarticular screw + ACDF) would be a reasonable treatment option to acquire an immediate stabilization in the management of compression-extension injury with unilateral facet joint fracture. However, clinical aspects must also be regarded when choosing a reconstruction method for a specific patient.

Comparative biomechanical analysis of rigid vs. flexible fixation devices for the lumbar spine: A geometrically patient-specific poroelastic finite element study

BACKGROUND AND OBJECTIVE

Lumbar spinal stenosis (LSS), or the narrowing of the spinal canal, continues to be the leading preoperative diagnosis for adults older than 65 years who undergo spine surgery. Although the treatment of LSS depends on its severity, the optimal surgical technique towards decreasing the risk of adjacent segment disease (ASD) remains elusive. This study aimed to comparatively analyze spinal biomechanics with rigid and flexible fixation devices (i.e., rigid and dynamic posterolateral fusion (PLF) and interspinous process (ISP) devices) during daily activities.

METHODS

Using a validated parametric poroelastic finite element modeling approach, 8 subject-specific pre-operative models were developed, and their validity was evaluated. Parametric FE models of the lumbar spines were then regenerated based on post-operation images for (A) rigid PLF (B) dynamic PLF (C) rigid ISP device (Coflex) and (D) flexible ISP device (DIAM) at L4-L5 level. Biomechanical responses for instrumented and adjacent intervertebral discs (IVDs) were analyzed and compared subject to static and cyclic loading.

RESULTS

The preoperative models were well comparable with previous works in literature. The postoperative results for the PLF and Coflex rigid systems, demonstrated greater ROM; higher values of stress and strain in the AF region; and increased disc height and fluid loss at the adjacent levels, as compared with the pre-op models and the post-op results of the flexible systems (i.e., dynamic PLF and DIAM). The calculated forces on the facet joint were of smaller magnitude for the ISP devices as compared to the PLF, particularly during extension.

CONCLUSIONS

This study demonstrates that the dynamic PLF construct and DIAM implants could be effective to maintain the natural poroelastic characteristics of adjacent IVDs, which could be beneficial for enhancing long-term clinical outcomes. FEM provides clinicians with an invaluable patient-specific quantitative tool for informed surgical planning and discerning follow-up management.

Influence of posterior pedicle screw fixation at L4-L5 level on biomechanics of the lumbar spine with and without fusion: a finite element method

Background

Posterior pedicle screw (PS) fixation, a common treatment method for widespread low-back pain problems, has many uncertain aspects including stress concentration levels, effects on adjacent segments, and relationships with physiological motions. A better understanding of how posterior PS fixation affects the biomechanics of the lumbar spine is needed. For this purpose, a finite element (FE) model of a lumbar spine with posterior PS fixation at the L4–L5 segment level was developed by partially removing facet joints (FJs) to imitate an actual surgical procedure. This FE study aimed to investigate the influence of the posterior PS fixation system on the biomechanics of the lumbar spine before and after fusion by determining which physiological motions have the most increase in posterior instrumentation (PI) stresses and FJ loading.

Results

It was determined that posterior PS fixation increased FJ loading by approximately 35% and 23% at the L3–L4 adjacent level with extension and lateral bending motion, respectively. This increase in FJ loading at the adjacent level could point to the possibility that adjacent segment disease has developed or progressed after posterior lumbar interbody fusion. Furthermore, analyses of peak von Mises stresses on PI showed that the maximum PI stresses of 272.1 MPa and 263.7 MPa occurred in lateral bending and flexion motion before fusion, respectively.

Conclusions

The effects of a posterior PS fixation system on the biomechanics of the lumbar spine before and after fusion were investigated for all physiological motions. This model could be used as a fundamental tool for further studies, providing a better understanding of the effects of posterior PS fixation by clearing up uncertain aspects.

Finite element study on the influence of pore size and structure on stress shielding effect of additive manufactured spinal cage

The stress shielding effect occurs when the orthopedic implant reduces the load delivered to the bone, causing inefficient stress transfer to the host bone. The usage of porous additive manufactured (AM) cages reduces the stress shielding effect and promotes bone ingrowth also. The purpose of this work is to study the stress and deformation on porous hybrid spinal cages under different loading conditions using Finite Element Analysis (FEA). The spinal cages consisting of three porous structures with pore sizes ranging from 0.4 to 0.6 mm were investigated for stress shielding and fatigue strength. The results showed a significant reduction in stress shielding for the studied designs and conclude that the pore size has a greater significant effect on stress shielding than the porous structure in spinal cages.

Biomechanical effects of lumbar fusion surgery on adjacent segments using musculoskeletal models of the intact, degenerated and fused spine

Adjacent segment disorders are prevalent in patients following a spinal fusion surgery. Postoperative alterations in the adjacent segment biomechanics play a role in the etiology of these conditions. While experimental approaches fail to directly quantify spinal loads, previous modeling studies have numerous shortcomings when simulating the complex structures of the spine and the pre/postoperative mechanobiology of the patient. The biomechanical effects of the L4–L5 fusion surgery on muscle forces and adjacent segment kinetics (compression, shear, and moment) were investigated using a validated musculoskeletal model. The model was driven by in vivo kinematics for both preoperative (intact or severely degenerated L4–L5) and postoperative conditions while accounting for muscle atrophies. Results indicated marked changes in the kinetics of adjacent L3–L4 and L5–S1 segments (e.g., by up to 115% and 73% in shear loads and passive moments, respectively) that depended on the preoperative L4–L5 disc condition, postoperative lumbopelvic kinematics and, to a lesser extent, postoperative changes in the L4–L5 segmental lordosis and muscle injuries. Upper adjacent segment was more affected post-fusion than the lower one. While these findings identify risk factors for adjacent segment disorders, they indicate that surgical and postoperative rehabilitation interventions should focus on the preservation/restoration of patient’s normal segmental kinematics.

Stress analysis of intervertebral disc during occupational activities

BACKGROUND AND OBJECTIVE

Manual material handling activities cause large compression of the intervertebral disc of the lumbar spine. Intradiscal pressure (IDP) has generally been employed to predict the risk of low back injury. As an alternative to in vivo measurements, either motion analysis or finite element (FE) analysis has been used to estimate IDP. The purpose of this study is to propose a new biomechanical method that integrates FE analysis with motion analysis, in order to estimate the stresses and deformations of the intervertebral disc of the lumbar spine during occupational activities.

METHODS

In the proposed method, motion analysis is performed first by using motion capture data, and the results are employed as input data to FE analysis at specific times of interest during motion. In this method, an in-house interface program is used to scale an initial reference FE model to the subject of experiment, and transformed to the corresponding posture at a specific time during motion. The muscle forces and GRF obtained from motion analysis are applied to FE analysis as boundary and loading conditions. For a total of eighteen occupational activities, the IDP, shear stress, and strain of the L4-L5 segment are estimated.

RESULTS

Under each in vivo activity, the predicted IDP was in overall agreement with the available in vivo data. For lifting activities according to lift origin position, the maximum IDP occurred in the far-knee position immediately after lifting. As the lift origin position moved away from the spine, the stresses and strains in the disc increased.

CONCLUSIONS

This new proposed method is expected to allow the estimation of the stresses and deformations in the intervertebral disc during various occupational activities.

Biomechanical modeling of spinal ligaments: finite element analysis of L4-L5 spinal segment

The complex mechanical structure of spine is usually simplified in finite element (FE) modes. In this study, different 3D models of L4-L5 spinal segment distinguished by their ligament modelling were developed (1D truss, 2D shell and 3D space truss elements). All models could be considered validated with respect to range of motion and intradiscal pressure, although their ligament stresses/forces were substantially different. The models with 2D shell and 3D space truss ligaments showed the stress distribution and identified the potential failure/injury locations in ligaments. The model with 3D space truss ligaments showed the stress/force direction (representing collagen fiber directions).

Development of a flexible instrumented lumbar spine finite element model and comparison with in-vitro experiments

Surgical corrections of degenerative lumbar scoliosis and sagittal malalignment are associated with significant complications, such as rod fractures and pseudarthrosis, particularly in the lumbosacral junction. Finite element studies can provide relevant insights to improve performance of spinal implants. The aim of the present study was to present the development of non-instrumented and instrumented Finite Element Models (FEMs) of the lumbopelvic spine and to compare numerical results with experimental data available in the literature. The lumbo-pelvic spine FEM was based on a CT-scan from an asymptomatic volunteer representing the 50th percentile male. In a first step a calibration of mechanical properties was performed in order to obtain a quantitative agreement between numerical results and experimental data for defect stages of spinal segments. Then, FEM results were compared in terms of range of motion and strains in rods to in-vitro experimental data from the literature for flexible non-instrumented and instrumented lumbar spines. Numerical results from the calibration process were consistent with experimental data, especially in flexion. A positive agreement was obtained between FEM and experimental results for the lumbar and sacroiliac segments. Instrumented FEMs predicted the same trends as experimental in-vitro studies. The instrumentation configuration consisting of double rods and an interbody cage at L5-S1 maximally reduced range of motion and strains in main rods and thus had the lowest risk of pseudarthrosis and rod fracture. The developed FEMs were found to be consistent with published experimental results; therefore they can be used for further post-operative complication investigations.

Exploring lumbar and lower limb kinematics and kinetics for evidence that lifting technique is associated with LBP

Purpose

To investigate if lumbar and lower limb kinematics or kinetics are different between groups with and without a history of LBP during lifting. Secondly, to investigate relationships between biomechanical variables and pain ramp during repeated lifting.

Methods

21 LBP and 20 noLBP participants completed a 100-lift task, where lumbar and lower limb kinematics and kinetics were measured during lifting, with a simultaneous report of LBP intensity every 10 lifts. Lifts were performed in a laboratory setting, limiting ecological validity.

Results

The LBP group used a different lifting technique to the noLBP group at the beginning of the task (slower and more squat-like). Kinetic differences at the beginning included less peak lumbar external anterior shear force and greater peak knee power demonstrated by the LBP group. However, at the end of the task, both groups lifted with a much more similar technique that could be classified as more stoop-like and faster. Peak knee power remained greater in the LBP group throughout and was the only kinetic difference between groups at the end of the lifting task. While both groups lifted using a more comparable technique at the end, the LBP group still demonstrated a tendency to perform a slower and more squat-like lift throughout the task. Only one of 21 variables (pelvic tilt at box lift-off), was associated with pain ramp in the LBP group. Conclusions: Workers with a history of LBP, lift with a style that is slower and more squat-like than workers without any history of LBP. Common assumptions that LBP is associated with lumbar kinematics or kinetics such as greater lumbar flexion or greater forces were not observed in this study, raising questions about the current paradigm around ‘safe lifting’.

Biomechanical comparison of four types of instrumentation constructs for revision surgery in lumbar adjacent segment disease: A finite element study

BACKGROUND

Different constructs are applied in revision surgery (RS) for adjacent segment disease (ASD) aiming to further decompress and fixate the affected segment(s) in two ways: replacing or preserving the primary implants. This study aimed to compare the biomechanical properties of four constructs with different configurations.

METHODS

An T12-L5 finite element (FE) model was constructed and validated. Primary surgery was performed at L4-L5 and instrumented from L3 to L5. Thereafter, RS was undertook by decompressing L2-L3 and fixated with implant-replacing construct A, or implant-preserving construct B, C or D. Range of motion (ROM) and intervertebral disc pressure (IDP) were compared. Maximum von Mises stress on the rods between Construct A and B was evaluated.

RESULTS

An obvious reduction of ROM was observed when the FE model was instrumented with four constructs respectively. The overall changing characteristics of ROM were approximately identical among four constructs. The changing characteristic of IDP among four constructs was similar. The degree of IDP reduction of Construct B was comparable to Construct A, while that of Construct C was comparable to Construct D. Maximum von Mises stress on the rods between Construct A and B indicated that no stress concentration was recorded at the locking part of the connector rod.

CONCLUSIONS

The biomechanics of implant-preserving constructs were comparable to the traditional implant-replacing construct. The location of side-by-side connector could not affect the stability of Construct C and D. Construct B might be an optimal choice in RS for less dissection, less complication and more convenience in manipulation.

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.

Intervertebral Disc Mechanics With Nucleotomy: Differences Between Simple and Dual Loading

Painful herniated discs are treated surgically by removing extruded nucleus pulposus (NP) material (nucleotomy). NP removal through enzymatic digestion is also commonly performed to initiate degenerative changes to study potential biological repair strategies. Experimental and computational studies have shown a decrease in disc stiffness with nucleotomy under single loading modalities, such as compression-only or bending-only loading. However, studies that apply more physiologically relevant loading conditions, such as compression in combination with bending or torsion, have shown contradicting results. We used a previously validated bone-disc-bone finite element model (Control) to create a Nucleotomy model to evaluate the effect of dual loading conditions (compression with torsion or bending) on intradiscal deformations. While disc joint stiffness decreased with nucleotomy under single loading conditions, as commonly reported in the literature, dual loading resulted in an increase in bending stiffness. More specifically, dual loading resulted in a 40% increase in bending stiffness under flexion and extension and a 25% increase in stiffness under lateral bending. The increase in bending stiffness was due to an increase and shift in compressive stress, where peak stresses migrated from the NP-annulus interface to the outer annulus. In contrast, the decrease in torsional stiffness was due to greater fiber reorientation during compression. In general, large radial strains were observed with nucleotomy, suggesting an increased risk for delamination or degenerative remodeling. In conclusion, the effect of nucleotomy on disc mechanics depends on the type and complexity of applied loads.

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

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

Biomechanical influence of the surgical approaches, implant length and density in stabilizing ankylosing spondylitis cervical spine fracture

Ankylosing spondylitis cervical spine fractures (ASCFs) are particularly unstable and need special consideration when selecting appropriate internal fixation technology. However, there is a lack of related biomechanical studies. This study aimed to investigate the biomechanical influence of the pattern, length, and density of instrumentation for the treatment of ASCF. Posterior, anterior, and various combined fixation approaches were constructed using the finite element model (FEM) to mimic the surgical treatment of ASCFs at C5/6. The rate of motion change (RMC) at the fractured level and the internal stress distribution (ISD) were observed. The results showed that longer segments of fixation and combined fixation approaches provided better stability and lowered the maximal stress. The RMC decreased more significantly when the length increased from 1 to 3 levels (302% decrease under flexion, 134% decrease under extension) than from 3 to 5 levels (22% decrease under flexion, 23% decrease under extension). Longer fixation seems to be more stable with the anterior/posterior approach alone, but 3-level posterior fixation may be the most cost-effective option. It is recommended to perform surgery with combined approaches, which provide the best stability. Long skipped-screwing posterior fixation is an alternative technique for use in ASCF patients.

Validation of a finite element model of the thoracolumbar spine to study instrumentation level variations in early onset scoliosis correction

Growth-guidance constructs are an alternative to growing rods for the surgical treatment of early onset scoliosis (EOS). Constructs containing ultra-high molecular weight polyethylene (UHMWPE) sublaminar tape have been proposed as an improvement to the traditional Luque trolley. Ideally, a certain minimum number of levels is instrumented, thus offering the best balance between providing adequate spinal fixation and minimizing surgical exposure and spinal mobility reduction. The objective of the current study was to validate a parametric FE model of the thoracolumbar spine including its ability to predict the biomechanical effects of varying the number of levels instrumented with UHMWPE sublaminar tape in a growth-guidance construct for EOS correction. In a first step, the material properties of the L4-L5 segment in the model were calibrated relative to literature data. Next, whole thoracolumbar spine behavior was verified relative to literature data as well. Subsequently, rods, screws, and sublaminar tape were implemented in the model and a simulation of a previously performed in vitro experiment, in which the range of motion (ROM) of porcine spine segments was measured for different tape configurations, was performed. Good agreement between in vitro and FE-results was found for the changes in ROM before and after instrumentation. Good agreement for changes in ROM was obtained when varying the number of instrumented levels as well, indicating that the model can be a useful tool to evaluate the effects of construct composition variations. The present study was limited by the fact that only normal spine curvatures were analyzed and the fact that results of porcine spine experiments were compared to results of human FE models. Nevertheless, the good agreement in results, even at a detailed level, supports the idea that the model can ultimately be used as a pre-operative planning tool to evaluate different construct designs. The FE model of the thoracolumbar spine was successfully validated and was able to capture the biomechanical effect of construct component variations.