Effect of graded facetectomy on biomechanics of Dynesys dynamic stabilization system

Mechanical influence of facet tropism in patients with chronic discogenic pain disorder

Aims

The presence of facet tropism has been correlated with an elevated susceptibility to lumbar disc pathology. Our objective was to evaluate the impact of facet tropism on chronic lumbosacral discogenic pain through the analysis of clinical data and finite element modelling (FEM).

Methods

Retrospective analysis was conducted on clinical data, with a specific focus on the spinal units displaying facet tropism, utilizing FEM analysis for motion simulation. We studied 318 intervertebral levels in 156 patients who had undergone provocation discography. Significant predictors of clinical findings were identified by univariate and multivariate analyses. Loading conditions were applied in FEM simulations to mimic biomechanical effects on intervertebral discs, focusing on maximal displacement and intradiscal pressures, gauged through alterations in disc morphology and physical stress.

Results

A total of 144 discs were categorized as 'positive' and 174 discs as 'negative' by the results of provocation discography. The presence of defined facet tropism (OR 3.451, 95% CI 1.944 to 6.126) and higher Adams classification (OR 2.172, 95% CI 1.523 to 3.097) were important predictive parameters for discography-'positive' discs. FEM simulations showcased uneven stress distribution and significant disc displacement in tropism-affected discs, where loading exacerbated stress on facets with greater angles. During varied positions, notably increased stress and displacement were observed in discs with tropism compared to those with normal facet structure.

Conclusion

Our findings indicate that facet tropism can contribute to disc herniation and changes in intradiscal pressure, potentially exacerbating disc degeneration due to altered force distribution and increased mechanical stress.

Biomechanical Analysis of a Newly Proposed Surgical Combination (MIS Screw-Rod System for Indirect Decompression+ Interspinous Fusion System for long Term Spinal Stability) in Treatment of Lumbar Degenerative Diseases

OBJECTIVE

The aim of this study is to analyze the biomechanical stability of a newly proposed surgical combination (minimally invasive surgery of screw-rod system for indirect decompression + interspinous fusion system for long term spinal stability) in treatment of lumbar degenerative diseases.

METHODS

The three-dimensional (3D) computed tomography (CT) image data of an adult healthy male volunteer were selected. An intact model of L4/5 was further established and validated by using Mimic and 3-matic, 3D slicer, abaqus, Python. Four surgical models were constructed. The biomechanical stability among these surgical modes was compared and analyzed using finite element analysis.

RESULTS

The maximum von mises on fixation system in surgical models 2 and 3 exhibited comparable values. This finding suggested that the increase in interspinous fusion did not result in a significant elevation in maximum von mises on fixation system. Compared with the third surgical model, the fourth model, which received less average von mises experienced by the screw in contact with both cancellous and cortical bone. The findings indicated that the inclusion of facet joint fusion in surgical procedures might not be necessary to increase the average von Mises stress experienced by the screw in contact with both cancellous and cortical bone.

CONCLUSIONS

The biomechanical stability of the newly proposed surgical combination (MIS screw-rod for indirect decompression + interspinous fusion for long term spinal stability technique) was not lower than that of the other surgical combination groups, and it might not be necessary to perform facet joint fusion during the surgery.

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

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

Hybrid cortical bone trajectory and modified cortical bone trajectory techniques in transforaminal lumbar interbody fusion at L4-L5 segment: A finite element analysis

Background

The academia has increasingly acknowledged the superior biomechanical performance of the hybrid fixation technique in recent years. However, there is a lack of research on the hybrid fixation technique using BCS (Bilateral Cortical Screws) and BMCS (Bilateral Modified Cortical Screws). This study aims to investigate the biomechanical performance of the BCS and BMCS hybrid fixation technique in transforaminal lumbar interbody fusion (TLIF) at the L4-L5 segment in a complete lumbar-sacral finite element model.

Methods

Three cadaver specimens are used to construct three lumbar-sacral finite element models. The biomechanical properties of various fixation technologies (BCS-BCS, BMCS-BMCS, BMCS-BCS, and BCS-BMCS) are evaluated at the L4-5 segment with a TLIF procedure conducted, including the range of motion (ROM) of the L4-5 segment, as well as the stress experienced by the cage, screws, and rods. The testing is conducted under specific loading conditions, including a compressive load of 400 N and a torque of 7.5Nm, subjecting the model to simulate flexion, extension, lateral bending, and rotation.

Results

No significant variations are seen in the ROM at the L4-5 segment when comparing the four fixation procedures during flexion and extension. However, when it comes to lateral bending and rotation, the ROM is ordered in descending order as BCS-BCS, BCS-BMCS, BMCS-BMCS, and BMCS-BCS. The maximum stress experienced by the cage is observed to be highest within the BMCS-BCS technique during movements including flexion, extension, and lateral bending. Conversely, the BMCS-BMCS technique exhibits the highest cage stress levels during rotational movements. The stress applies to the screws and rods order the sequence of BCS-BCS, BCS-BMCS, BMCS-BCS, and BMCS-BMCS throughout all four working conditions.

Conclusion

The BMCS-BCS technique shows better biomechanical performance with less ROM and lower stress on the internal fixation system compared to other fixation techniques. BMCS-BMCS technology has similar mechanical performance to BMCS-BCS but has more contact area between screws and cortical bone, making it better for patients with severe osteoporosis.

Asymmetric Load Transmission Induces Facet Joint Subchondral Sclerosis and Hypertrophy in Patients with Idiopathic Adolescent Scoliosis: Evaluation Using Finite Element Model and Surgical Specimen

Adolescent idiopathic scoliosis (AIS) with thoracic curvature primarily progresses from the thoracolumbar region, causing abnormal twisting and rotation of the spinal column. This results in unbalanced, asymmetric loads on each vertebrae and increased demands on the thoracic facet joints to withstand rotational stress from adjacent vertebrae. However, no studies have focused on the stress distribution on the facet joints of the thoracic spine in patients with AIS. This study aimed to investigate the mechanical loading and its distribution on the thoracic facet joints of AIS patients using finite element (FE) analysis and surgical specimens. FE models of the thoracic spine were created from a total of 13 female AIS patients (Lenke type 1, n = 4; Lenke type 2, n = 4; Lenke type 3, n = 5). A load of 200 N on the T3 vertebrae and 30 N each on the bilateral superior articular processes were applied vertically to quantify the contact force on the facet joints from T3 to T11. In addition, morphological and histological analyses were performed on the inferior articular processes obtained during surgery. FE analysis demonstrated that contact forces of the facet joint progressively increased from the mid to lower thoracic spine of the concave side, reaching a maximum around the apex. More than 91% of the load was transmitted by the facet joints at the concave side, resulting in facet joint subchondral sclerosis and hypertrophy. The apical facet joint in AIS helps counteract rotational stress between vertebrae and transfers most stress through the concave side. In conclusion, this study found that asymmetric load transfer in the facet joints leads to subchondral sclerosis and hypertrophy. These findings can enhance our understanding of the stress loading on facet joints and the resulting biological changes and help clarify the mechanisms involved in scoliosis progression. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

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.

Biomechanical evaluation of the hybrid pedicle screw-cortical bone trajectory technique in transforaminal lumbar interbody fusion to adjacent segment degeneration-finite element analysis

BACKGROUND

Transforaminal lumbar interbody fusion is an effective surgical treatment of intervertebral disk herniation. However, its clinical efficacy for adjacent segment disk degeneration (ASDD) after hybrid bilateral pedicle screw - bilateral cortical screw (pedicle screw at L4 and cortical bone trajectory screw at L5) and hybrid bilateral cortical screw - bilateral pedicle screw (bilateral cortical screw at L4 and bilateral pedicle screw at L5) remains undiscovered. Therefore, the aim of this study is to evaluate the effect of the hybrid bilateral pedicle screw - bilateral cortical screw and hybrid bilateral cortical screw - bilateral pedicle screw on the adjacent segment via a 3-dimensional (3D) finite element (FE) analysis.

METHODS

Four human cadaveric lumbar spine specimens were provided by the anatomy teaching and research department of Xinjiang Medical University. Four finite element models of L1-S1 lumbar spine segment were generated. For each of these, four lumbar transforaminal lumbar interbody fusion models at L4-L5 segment with the following instruments were created: hybrid bilateral pedicle screw - bilateral cortical screw, bilateral cortical screw - bilateral cortical screw (bilateral cortical screw at both L4 and L5 segments), bilateral pedicle screw - bilateral pedicle screw (bilateral pedicle screw at both L4 and L5 segments), and hybrid bilateral cortical screw - bilateral pedicle screw. A 400-N compressive load with 7.5 Nm moments was applied for the simulation of flexion, extension, lateral bending, and rotation. The range of motion of L3-L4 and L5-S1 segments and von Mises stress of the intervertebral disc at the adjacent segment were compared.

RESULTS

Hybrid bilateral pedicle screw - bilateral cortical screw has the lowest range of motion at L3-L4 segment in flexion, extension, and lateral bending, and the highest disc stress in all motions, while the range of motion at L5-S1 segment and disc stress was lower than bilateral pedicle screw - bilateral pedicle screw in flexion, extension, and lateral bending, and higher than bilateral cortical screw - bilateral cortical screw in all motions. The range of motion of hybrid bilateral cortical screw - bilateral pedicle screw at L3-L4 segment was lower than bilateral pedicle screw - bilateral pedicle screw and higher than bilateral cortical screw - bilateral cortical screw in flexion, extension, and lateral bending, and the range of motion at L5-S1 segment was higher than bilateral pedicle screw - bilateral pedicle screw in flexion, lateral bending, and axial rotation. The disc stress at L3-L4 segment was lowest and more dispersed in all motions, and the disc stress at L5-S1 segment was higher than bilateral pedicle screw - bilateral pedicle screw in lateral bending and axial rotation, but more dispersed.

CONCLUSION

Hybrid bilateral cortical screw - bilateral pedicle screw decreases the impact on adjacent segments after spinal fusion, reduces the iatrogenic injury to the paravertebral tissues, and provides throughout decompression of the lateral recess.

Hybrid pedicle screw and modified cortical bone trajectory technique in transforaminal lumbar interbody fusion at L4-L5 segment: finite element analysis

BACKGROUND

Investigate the biomechanical properties of the hybrid fixation technique with bilateral pedicle screw (BPS) and bilateral modified cortical bone trajectory screw (BMCS) in L4-L5 transforaminal lumbar interbody fusion (TLIF).

METHODS

 Three finite element (FE) models of the L1-S1 lumbar spine were established according to the three human cadaveric lumbar specimens. BPS-BMCS (BPS at L4 and BMCS at L5), BMCS-BPS (BMCS at L4 and BPS at L5), BPS-BPS (BPS at L4 and L5), and BMCS-BMCS (BMCS at L4 and L5) were implanted into the L4-L5 segment of each FE model. The range of motion (ROM) of the L4-L5 segment, von Mises stress of the fixation, intervertebral cage, and rod were compared under a 400-N compressive load with 7.5 Nm moments in flexion, extension, bending, and rotation.

RESULTS

 BPS-BMCS technique has the lowest ROM in extension and rotation, and BMCS-BMCS technique has the lowest ROM in flexion and lateral bending. The BMCS-BMCS technique showed maximal cage stress in flexion and lateral bending, and the BPS-BPS technique in extension and rotation. Compared to the BPS-BPS and BMCS-BMCS technique, BPS-BMCS technique presented a lower risk of screw breakage and BMCS-BPS technique presented a lower risk of rod breakage.

CONCLUSION

 The results of this study support that the use of the BPS-BMCS and BMCS-BPS techniques in TLIF surgery for offering the superior stability and a lower risk of cage subsidence and instrument-related complication.

Computational lumbar spine models: A literature review

BACKGROUND

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

METHODS

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

FINDINGS

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

INTERPRETATION

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

The influence of the viscoelastic property of polycarbonate urethane as an artificial disc core material under various physiological motions at the L4-L5 level

Intervertebral disc (IVD) degeneration is one of the musculoskeletal disorders due to the Degenerative Disc Disease (DDD), that cause low back pain (LBP) and leads to a reduced range of motion. Spinal fusion and arthroplasty are the other surgical procedures that could replace the disc affected by DDD against artificial disc replacement (ADR). This study aims to analyse the biomechanical behaviour of proposed core material as Polycarbonate Urethane (PCU) in the L4-L5 lumbar segment for ADR with Ti-6Al-4V and Co-28Cr-6M as endplate materials and compare it to the performance of an ultra-high molecular weight polyethylene (UHMWPE) core. Finite element methods have been approached to measure the overall stress distribution along with other physiological motions like Flexion (FLEX), Extension (EXT), Axial rotation (AR) and Lateral bending (LB), respectively. Preload of 450 N compressive load, 8 N-m for Flex, 6 N-m for EXT, 6 N-m for AR and 4 N-m for LB are applied. It could be concluded that Ti-6Al-4V - PCU and Co-28Cr-6M - PCU is the best composition for the ADR for the L4-L5 level.

Biomechanical Comparison between Isobar and Dynamic-Transitional Optima (DTO) Hybrid Lumbar Fixators: A Lumbosacral Finite Element and Intersegmental Motion Analysis

Biomechanical performance of longitudinal component in dynamic hybrid devices was evaluated to display the load-transfer effects of Dynesys cord spacer or Isobar damper-joint dynamic stabilizer on junctional problem based on various disc degenerations. The dynamic component was adapted at the mildly degenerative L3–L4 segment, and the static component was fixed at the moderately degenerative L4–L5 segment under a displacement-controlled mode for the finite element study. Furthermore, an intersegmental motion behavior was analyzed experimentally on the synthetic model under a load-controlled mode. Isobar or DTO hybrid fixator could reduce stress/motion at transition segment, but compensation was affected at the cephalic adjacent segment more than the caudal one. Within the trade-off region (as a motion-preserving balance between the transition and adjacent segments), the stiffness-related problem was reduced mostly in flexion by a flexible Dynesys cord. In contrast, Isobar damper afforded the effect of maximal allowable displacement (more than peak axial stiffness) to reduce stress within the pedicle and at facet joint. Pedicle-screw travel at transition level was related to the extent of disc degeneration in Isobar damper-joint (more than Dynesys cord spacer) attributing to the design effect of axial displacement and angular rotation under motion. In biomechanical characteristics relevant to clinical use, longitudinal cord/damper of dynamic hybrid lumbar fixators should be designed with less interface stress occurring at the screw-vertebral junction and facet joint to decrease pedicle screw loosening/breakage under various disc degenerations.

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

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

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

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

Biomechanical comparison of vertebral augmentation and cement discoplasty for the treatment of symptomatic Schmorl's node: a finite element analysis

Percutaneous vertebral augmentation (PVA) and percutaneous cement discoplasty (PCD) are two relatively new minimally invasive surgeries for symptomatic Schmorl's reported in recent decade. However, the clinical evidence for the effectiveness of these two surgeries is insufficient. The purpose of this study was to compare the biomechanical benefits and risks of the two surgeries in order to analyze their biomechanical differences and effectiveness. We reconstructed Five lumbar finite element models via computed tomography data, including control model, PVA-ideal model, PVA-nonideal model, PCD-ideal model, and PCD-nonideal model. The stress and strain of Schmorl's nodes, bone marrow edema zone (BMEZ), affected endplate, and the overall stability of segment were analyzed and compared. The validity of our models was confirmed. As a result, the PVA-ideal model can significantly reduce the stress of Schmorl's node and the strain of BMEZ, while this effect is inappreciable in PVA-nonideal model. The PCD-ideal model significantly reduced the strain of Schmorl's nodes and BMEZ, and significantly improve segmental stability, but also resulted in a significant increase in the stress of Schmorl's nodes, BMEZ and endplates. The PCD-nonideal model not only lacks blocking effect, but also sharply increases the strain of Schmorl's nodes and BMEZ. Thus, We recommend that both PVA and PCD surgeries in ideal distribution facilitated a more stable paranodular biomechanical microenvironment. However, due to the possibility of poor biomechanical outcomes caused by the non-ideal cement distribution, the non-ideal distribution of bone cement needs to be remedied in practice.

Relationship between the elastic modulus of the cage material and the biomechanical properties of transforaminal lumbar interbody fusion: A logarithmic regression analysis based on parametric finite element simulations

BACKGROUND AND OBJECTIVE

Conventional method for evaluating the biomechanical effects of a specific elastic modulus of cage (cage-E) on spinal fusions requires establishing a "one-on-one" biomechanical model, which seems laborious and inefficient when dealing with the emergence of numerous cage materials with various cage-Es. We aim to offer a much convenient method to instantly predicting the biomechanical effects of any targeted cage-E on transforaminal lumbar interbody fusion (TLIF) by using a parametric finite element (FE) analysis to determining the regression relationship between cage-E and biomechanical properties of TLIF.

MATERIALS AND METHODS

A L4/5 FE TLIF construct was modeled. Cage-E was linearly increased from 0.1 GPa (cancellous bone) to 110 GPa (titanium alloy). The function equations for assessing the influence of cage-E on the biomechanical indexes of TLIF were established using a logarithmic regression analysis.

EXPERIMENTAL RESULTS

As cage-E increased from 0.1 GPa to 110 GPa, all the biomechanical indexes initially increased or decayed rapidly, and then slowed over time. Logarithmic regression models and functional equations were successfully established between cage-E and these indexes (P<0.0001). Their determination coefficients ranged from 0.72 to 0.99. The range of motions decreased from 0.37-1.10° to 0.20-1.07°. The mean stresses of the central and peripheral grafts reduced from 0.10-0.41 and 0.25-0.42 MPa to 0.03-0.04 and 0.19-0.27 MPa, respectively. In addition, the maximum stresses of the screw-bone interface and posterior instrumentation reduced from 11.76-25.04 and 8.91-84.68 MPa to 9.71-18.92 and 6.99-70.59 MPa, respectively. Finally, the maximum stresses of the cage and endplate increased from 0.28-1.35 MPa and 3.90-8.63 MPa to 14.86-36.16 MPa and 11.01-36.55 MPa, respectively.

CONCLUSIONS

The decrease of cage-E reduces the risks of cage subsidence, cage breakage, and pseudarthrosis, while increasing the risk of instrumentation failure. The logarithmic regression models optimally demonstrate the relationship between cage-E and biomechanical properties of TLIF. The functional equations based on these models can be adopted to predict the biomechanical effects of any targeted cage-Es on TLIF, which effectively simplifies the procedures for the biomechanical assessments of cage materials.

Biomechanical analysis of stand-alone lumbar interbody cages versus 360° constructs: an in vitro and finite element investigation

OBJECTIVE

Low fusion rates and cage subsidence are limitations of lumbar fixation with stand-alone interbody cages. Various approaches to interbody cage placement exist, yet the need for supplemental posterior fixation is not clear from clinical studies. Therefore, as prospective clinical studies are lacking, a comparison of segmental kinematics, cage properties, and load sharing on vertebral endplates is needed. This laboratory investigation evaluates the mechanical stability and biomechanical properties of various interbody fixation techniques by performing cadaveric and finite element (FE) modeling studies.

METHODS

An in vitro experiment using 7 fresh-frozen human cadavers was designed to test intact spines with 1) stand-alone lateral interbody cage constructs (lateral interbody fusion, LIF) and 2) LIF supplemented with posterior pedicle screw-rod fixation (360° constructs). FE and kinematic data were used to validate a ligamentous FE model of the lumbopelvic spine. The validated model was then used to evaluate the stability of stand-alone LIF, transforaminal lumbar interbody fusion (TLIF), and anterior lumbar interbody fusion (ALIF) cages with and without supplemental posterior fixation at the L4-5 level. The FE models of intact and instrumented cases were subjected to a 400-N compressive preload followed by an 8-Nm bending moment to simulate physiological flexion, extension, bending, and axial rotation. Segmental kinematics and load sharing at the inferior endplate were compared.

RESULTS

The FE kinematic predictions were consistent with cadaveric data. The range of motion (ROM) in LIF was significantly lower than intact spines for both stand-alone and 360° constructs. The calculated reduction in motion with respect to intact spines for stand-alone constructs ranged from 43% to 66% for TLIF, 67%-82% for LIF, and 69%-86% for ALIF in flexion, extension, lateral bending, and axial rotation. In flexion and extension, the maximum reduction in motion was 70% for ALIF versus 81% in LIF for stand-alone cases. When supplemented with posterior fixation, the corresponding reduction in ROM was 76%-87% for TLIF, 86%-91% for LIF, and 90%-92% for ALIF. The addition of posterior instrumentation resulted in a significant reduction in peak stress at the superior endplate of the inferior segment in all scenarios.

CONCLUSIONS

Stand-alone ALIF and LIF cages are most effective in providing stability in lateral bending and axial rotation and less so in flexion and extension. Supplemental posterior instrumentation improves stability for all interbody techniques. Comparative clinical data are needed to further define the indications for stand-alone cages in lumbar fusion surgery.

Long-term Outcomes after Microsurgical Decompression of Lumbar Foraminal Stenosis and Adverse Effects of Preoperative Scoliosis: A Prospective Cohort Study

Lumbar foraminal stenosis is a common disorder, with surgical treatment varying from simple decompression to interbody fusion. It is often associated with degenerative lumbar scoliosis, but the effects of scoliosis on outcomes are unclear. The objectives of this study were to clarify long-term outcomes after microsurgical decompression of lumbar foraminal stenosis through Wiltse’s approach and to determine the effects of scoliosis on these outcomes. A total of 86 consecutive patients with lumbar foraminal stenosis were prospectively followed after microsurgical decompression. They were categorized in multiple subcohorts with follow-up durations ranging from 6 months to 5 years. Outcomes were assessed using the Short Form 36 questionnaire (average physical scores and bodily pain scores). Local Cobb angle of the operative segment was measured preoperatively, and its effects on outcomes were analyzed. Average physical scores improved significantly from 33.8 (95% confidence interval [CI]: 29.1–38.5) preoperatively to 59.5 (95% CI: 54.6–64.3) at 6 months postoperatively and remained improved for 5 years. Bodily pain scores improved significantly from 23.7 (95% CI: 18.7–28.6) preoperatively to 56.3 (95% CI: 51.2–61.6) at 6 months postoperatively and remained improved for 5 years. Patients with preoperative scoliosis (local Cobb angle >10 degrees) had poorer outcomes: average physical scores were worse by 9.6 points (p = 0.07) and bodily pain scores were worse by 12.1 points (p = 0.02), compared with patients without scoliosis (local Cobb angle ≤10 degrees). Microsurgical foraminal decompression produced overall excellent outcomes in patients with lumbar foraminal stenosis. Preoperative scoliosis attenuated these beneficial effects.

3D pull-out finite element simulation of the pedicle screw-trabecular bone interface at strain rates

Biomedical experimental studies such as pull-out (PO), screw loosening experience variability mechanical properties of fresh bone, legal procedures of cadaver bone samples and time-consuming problems. Finite Element Method (FEM) could overcome experimental problems in biomechanics. However, material modelling of bone is quite difficult, which has viscoelastic and viscoplastic properties. The study presents a bone material model which is constructed at the strain rates with the Johnson-Cook (JC) material model, one of the robust constitutive material models. The JC material constants of trabecular bone are determined by the curve fitting method at strain rates for the 3D PO finite element simulation, which defines the screw-bone interface relationship. The PO simulation is performed using the Abaqus/CAE software program. Bone fracture mechanisms are simulated with dynamic/explicit solutions during the PO phenomenon. The paper exposes whether the strain rate has effects on the PO performance. Moreover, simulation reveals the relationship between pedicle screw diameter and PO performance. The results obtained that the maximum pull-out force (POF) improves as both the screw diameter and the strain rate increase. For 5.5 mm diameter pedicle screw POFs were 487, 517 and 1708 N at strain rate 0.00015, 0.015 and 0.015 s^-1, respectively. The FOFs obtained from the simulation of the other screw were 730, 802 and 2008 N at strain rates 0.00015, 0.0015 and 0.015, respectively. PO phenomenon was also simulated realistically in the finite element analysis (FEA).

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

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

Influence of structural and material property uncertainties on biomechanics of intervertebral discs - Implications for disc tissue engineering

This study investigated how variations of structural and material properties of human intervertebral discs (IVDs) affect the biomechanical responses of the IVDs under simulated physiological loading conditions using a stochastic finite element (SFE) model. An SFE method, which combined an anatomic FE model of human lumbar L3-4 segment and probabilistic analysis of its structural and material properties, was used to generate a dataset of 500 random disc samples with varying structural and material properties. The sensitivity of the biomechanical responses, including intervertebral displacements/rotations, intradiscal pressures (IDP), fiber stresses and matrix strains of annulus fibrosus (AF), were systematically quantified under various physiological loading conditions, including a 500N compression and 7.5Nm moments in the 3 primary rotations. Significant variations of the IDPs, IVD displacements/rotations, and stress/strain distributions were found using the dataset of 500 ramdom disc samples. Under all the loading conditions, the IDPs were positively correlated with the Poisson's ratio of the NP (r = 0.46 to 0.75, p = 0.004-0.001) and negatively with the Young's modulus of the annulus matrix (r = -0.48 to -0.65, p = 0.003-0.001). The primary intervertebral rotations were significantly affected by the Young's modulus of the annulus matrix (r = -0.44 to -0.71, p = 0.001-0.032) and the orientations of the annular fibers (r = -0.45 to -0.69, p = 0.001-0.029). The heterogeneity of structures and material properties of the IVD had distinct effects on the biomechanical performances of the IVD. These data could help improve our understanding of the intrinsic biomechanics of the IVD and provide references for optimal design of tissue engineered discs by controlling structural and material properties of the disc components.

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

BACKGROUND

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

AIM

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

METHODS

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

RESULTS

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

CONCLUSION

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

Preoperative facet joint arthropathy does not impact long-term clinical outcomes after lumbar-stability-preserving decompression and dynesys fixation

To evaluate the impact of the preoperative severity of facet joint arthropathy on long-term functional outcomes and spinopelvic parameters in patients undergoing lumbar-stability-preserving decompression and Dynesys fixation. In this retrospective study, 88 patients undergoing combined surgery at our hospital from 2008 to 2015 were included. The patients were divided into two groups, the less and more than mean degeneration groups, based on preoperative facet joint arthropathy of the index level(s). The clinical outcomes were the Visual Analogue Scale (VAS) score, the Oswestry Disability Index (ODI) score and spinopelvic parameters. The mean follow-up durations for the less and more than mean degeneration groups were 84.83 ± 27.58 and 92.83 ± 20.45 months, respectively. The combined surgery significantly improved VAS and ODI scores, and increased sacral slope (SS) regardless of preoperative arthropathy severity. In addition, facet joint arthropathy at adjacent levels continued to worsen after surgery in both arthropathy severity groups. Preoperative facet joint arthropathy did not influence most long-term clinical outcomes in patients undergoing lumbar-stability-preserving decompression and Dynesys fixation. This combined surgery may be suitable for patients with facet joint arthropathy regardless of disease severity.

Different fixation pattern for thoracolumbar fracture of ankylosing spondylitis: A finite element analysis

The objective of this study is to establish an ankylosing spondylitis (AS) thoracolumbar fracture finite element (FE) model and provide a proper posterior fixation choice from the biomechanical perspective. The ankylosing spondylitis T9-L5 FE model was built and the range of motion (ROM) was compared to previous studies. The L1 transverse fracture was simulated and was separately fixed by five different patterns. The pull force and yielding force of the screws, the von Mises stress of the internal fixation, and the displacement of fracture site were analyzed to evaluate the proper fixation pattern for thoracolumbar fracture of AS. ROM of AS model was obviously restricted comparing to the normal vertebral experimental data. All the fixation patterns can stabilize the fracture. At least four levels of fixation can reduce the von Mises stress of the internal fixation. Four levels fixation has a higher pull force than the six levels fixation. Skipped level fixation did not reduce the stress, pull force and yielding force. The kyphosis correction did not change the biomechanical load. At least 4 levels fixation was needed for AS thoracolumbar fracture. The cemented screws should be chosen in 4 levels fixation to increase the holding of the screws. The skipped fixation has no advantage. The kyphosis correction can be chosen after weighing the pros and cons.

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

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

Pedicle Screw System May Not Control Severe Spinal Rotational Instability

STUDY DESIGN

An in vitro biomechanical study.

OBJECTIVE

The purpose of this study is to discuss whether pedicle screw systems can control spinal rotational instability in a functional spinal unit of lumbar spine on human cadaver.

SUMMARY OF BACKGROUND DATA

Rotational experiments using deer lumbar cadaveric models showed that rotational range of motion (ROM) of the model fixed by a pedicle screw system with crosslinking after total facetectomy for both the sides was larger than that in the intact model, and stated that spinal rotational instability could not be controlled using a pedicle screw system.

METHODS

A rotation experiment using 10 functional spinal units (L3-4) of lumbar spine on human cadavers was performed by preparing the four models (intact model, damaged model, pedicle screw model, and crosslink (CL) model) in stages, then calculating and comparing rotational ROM among the four models.

RESULTS

Rotational ROM in the CL model was still larger than that of the intact model in all the samples. And, rotational ROM decreased in the order of damaged model >> pedicle screw model > CL model > intact model. Statistical analysis revealed significant differences between all models (P < 0.001).

CONCLUSIONS

Pedicle screw systems may not control severe spinal rotational instability in human lumbar cadaveric models with total facetectomy on both the sides. This may represent a major biomechanical drawback to the pedicle screw system.

LEVEL OF EVIDENCE

N/A.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

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

Evaluation of nucleus pulposus fluid velocity and pressure alteration induced by cartilage endplate sclerosis using a poro-elastic finite element analysis

The nucleus pulposus (NP) in the intervertebral disk (IVD) depends on diffusive fluid transport for nutrients through the cartilage endplate (CEP). Disruption in fluid exchange of the NP is considered a cause of IVD degeneration. Furthermore, CEP calcification and sclerosis are hypothesized to restrict fluid flow between the NP and CEP by decreasing permeability and porosity of the CEP matrix. We performed a finite element analysis of an L3-L4 lumbar functional spine unit with poro-elastic constitutive equations. The aim of the study was to predict changes in the solid and fluid parameters of the IVD and CEP under structural changes in CEP. A compressive load of 500 N was applied followed by a 10 Nm moment in extension, flexion, lateral bending, and axial rotation to the L3-L4 model with fully saturated IVD, CEP, and cancellous bone. A healthy case of L3-L4 physiology was then compared to two cases of CEP sclerosis: a calcified cartilage endplate and a fluid constricted sclerotic cartilage endplate. Predicted NP fluid velocity increased for the calcified CEP and decreased for the calcified + less permeable CEP. Decreased NP fluid velocity was prominent in the axial direction through the CEP due to a less permeable path available for fluid flux. Fluid pressure and maximum principal stress in the NP were predicted to increase in both cases of CEP sclerosis compared to the healthy case. The porous medium predictions of this analysis agree with the hypothesis that CEP sclerosis decreases fluid flow out of the NP, builds up fluid pressure in the NP, and increases the stress concentrations in the NP solid matrix.

Muscle-driven and torque-driven centrodes during modeled flexion of individual lumbar spines are disparate

Lumbar spine biomechanics during the forward-bending of the upper body (flexion) are well investigated by both in vivo and in vitro experiments. In both cases, the experimentally observed relative motion of vertebral bodies can be used to calculate the instantaneous center of rotation (ICR). The timely evolution of the ICR, the centrode, is widely utilized for validating computer models and is thought to serve as a criterion for distinguishing healthy and degenerative motion patterns. While in vivo motion can be induced by physiological active structures (muscles), in vitro spinal segments have to be driven by external torque-applying equipment such as spine testers. It is implicitly assumed that muscle-driven and torque-driven centrodes are similar. Here, however, we show that centrodes qualitatively depend on the impetus. Distinction is achieved by introducing confidence regions (ellipses) that comprise centrodes of seven individual multi-body simulation models, performing flexion with and without preload. Muscle-driven centrodes were generally directed superior–anterior and tail-shaped, while torque-driven centrodes were located in a comparably narrow region close to the center of mass of the caudal vertebrae. We thus argue that centrodes resulting from different experimental conditions ought to be compared with caution. Finally, the applicability of our method regarding the analysis of clinical syndromes and the assessment of surgical methods is discussed.

Finite element modelling of hybrid stabilization systems for the human lumbar spine

Intersomatic fusion is a very popular treatment for spinal diseases associated with intervertebral disc degeneration. The effects of three different hybrid stabilization systems on both range of motion and intradiscal pressure were investigated, as there is no consensus in the literature about the efficiency of these systems. Finite element simulations were designed to predict the variations of range of motion and intradiscal pressure from intact to implanted situations. After hybrid stabilization system implantation, L4-L5 level did not lose its motion completely, while L5-S1 had no mobility as a consequence of disc removal and fusion process. BalanC hybrid stabilization system represented higher mobility at the index level, reduced intradiscal pressure of adjacent level, but caused to increment in range of motion by 20% under axial rotation. Higher tendency by 93% to the failure was also detected under axial rotation. Dynesys hybrid stabilization system represented more restricted motion than BalanC, and negligible effects to the adjacent level. B-DYN hybrid stabilization system was the most rigid one among all three systems. It reduced intradiscal pressure and range of motion at the adjacent level except from motion under axial rotation being increased by 13%. Fracture risk of B-DYN and Dynesys Transition Optima components was low when compared with BalanC. Mobility of the adjacent level around axial direction should be taken into account in case of implantation with BalanC and B-DYN systems, as well as on the development of new designs. Having these findings in mind, it is clear that hybrid systems need to be further tested, both clinically and numerically, before being considered for common use.

Material properties of human lumbar intervertebral discs across strain rates

BACKGROUND CONTEXT

The use of finite element (FE) methods to study the biomechanics of the intervertebral disc (IVD) has increased over recent decades due to their ability to quantify internal stresses and strains throughout the tissue. Their accuracy is dependent upon realistic, strain-rate dependent material properties, which are challenging to acquire.

PURPOSE

The aim of this study was to use the inverse FE technique to characterize the material properties of human lumbar IVDs across strain rates.

STUDY DESIGN

A human cadaveric experimental study coupled with an inverse finite element study.

METHODS

To predict the structural response of the IVD accurately, the material response of the constituent structures was required. Therefore, compressive experiments were conducted on 16 lumbar IVDs (39±19 years) to obtain the structural response. An FE model of each of these experiments was developed and then run through an inverse FE algorithm to obtain subject-specific constituent material properties, such that the structural response was accurate.

RESULTS

Experimentally, a log-linear relationship between IVD stiffness and strain rate was observed. The material properties obtained through the subject-specific inverse FE optimization of the annulus fibrosus (AF) fiber and AF fiber ground matrix allowed a good match between the experimental and FE response. This resulted in a Young modulus of AF fibers (-MPa) to strain rate (ε˙, /s) relationship of YMAF=31.5ln(ε˙)+435.5, and the C₁₀ parameter of the Neo-Hookean material model of the AF ground matrix was found to be strain-rate independent with an average value of 0.68 MPa.

CONCLUSIONS

These material properties can be used to improve the accuracy, and therefore predictive ability of FE models of the spine that are used in a wide range of research areas and clinical applications.

CLINICAL SIGNIFICANCE

Finite element models can be used for many applications including investigating low back pain, spinal deformities, injury biomechanics, implant design, design of protective systems, and degenerative disc disease. The accurate material properties obtained in this study will improve the predictive ability, and therefore clinical significance of these models.

A novel narrow surface cage for full endoscopic oblique lateral lumbar interbody fusion: A finite element study

OBJECTIVE

To evaluate the mechanical response of a narrow surface cage that we designed for full endoscopic oblique lateral lumbar interbody fusion (FE-OL-LIF).

METHODS

A finite element (FE) model of lumbar spine was developed and validated. To simulate the FE-OL-LIF, the functional spinal unit (FSU) L4-L5 was assembled with a narrow surface polyetheretherketone (PEEK) cage, two rods and four screws. 500N load combined with 7.5Nm moment was applied to the surgical models. Effect of the cage width on the stress was studied.

RESULTS

Range of motion (ROM) in the surgical models significantly decreased by 88% in flexion, 91% extension, 85% in right and left lateral bending, 75% in right and left axial rotation as compared to the intact model. Width of the cage slightly decreased the ROM in all loading scenarios. Flexion produced the highest stress in the cages and endplates. In all loading cases, the maximum stresses of cages and endplates were both lower than their yield stress.

CONCLUSIONS

In engineering analysis, the novel narrow-surface cage had a strength to support spine activities. 9 mm width cage was recommended in FE-OL-LIF. This study provided engineering evidence and technical advice to improve the design of minimally invasive cage. Fatigue test and cadaver trial shall be improved.

Finite Element Based-Analysis for Pre and Post Lumbar Fusion of Adult Degenerative Scoliosis Patients

STUDY DESIGN

Pre-post cohort finite elements (FE).

OBJECTIVES

To investigate the effect of adjacent load transfer pre and post fusion surgery of lumbar scoliotic spines using FE models.

SUMMARY OF BACKGROUND DATA

Adult degenerative scoliosis (ADS) results from age-related changes, leading to segmental instability, deformity, and stenosis. FE study is capable of capturing the biomechanical parameters internal to the bones and connective soft tissues of the spine, which is difficult to measure by experimental approaches. Literature that describes the underlying mechanisms responsible for spinal fusion in scoliosis patients is limited, and FE study with larger subject sample size should be conducted.

METHODS

Twenty three-dimensional nonlinear FE models of the lumbosacral spine were created from pre (Cobb angle: 28.1° ± 10.5°) and post scoliosis surgery in vivo CT scans. During surgery, pedicle screws and rods were implanted at lumbar and sacral levels. A compressive load and six different moments (flexion, extension, right lateral bending, left lateral bending, right axial rotation, left axial rotation) were applied to the top level of each model. Outcome measures were range of motion (RoM), intradiscal pressure (IDP), and facet joint forces (FJF). Spinal fusion did alter the mechanical function of the scoliotic spine.

RESULTS

Scoliotic spine presented abnormal and asymmetrical kinetic and kinematic behavior. RoM: At the adjacent level, spinal fusion surgery produced a statically significantly increased left and right later bending intersegmental rotation (p < .006) in comparison to presurgical scoliosis models. At the fused level, spinal fusion surgery produced a statically significantly reduced intersegmental rotation in all the loading conditions (p = .001) in comparison to presurgical scoliosis models. IDP: At the fused level, spinal fusion surgery produced a much lower IDP in all of the loading conditions (p = .001). FJF: At the adjacent level, spinal fusion surgery produced a considerably larger left lateral rotation FJF (p = .001) in comparison to presurgical scoliosis models. At the fused level, spinal fusion surgery produced considerably lower FJF in all the loading conditions (p = .001) in comparison to presurgical scoliosis models.

CONCLUSIONS

This study was the first to investigate the effect of adjacent load transfer before and after fusion surgery using in vivo CT scans of 10 scoliotic spines. A posterior fusion has only a minor effect on mechanical behavior and a large effect on pressure and forces at the adjacent level. As expected, a large effect in the kinematics and kinetics was found at the fused level.

LEVEL OF EVIDENCE

Level 3.

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

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

Estrogen Deficiency Exacerbates Intervertebral Disc Degeneration Induced by Spinal Instability in Rats

STUDY DESIGN

An experimental animal study of osteoporosis (OP) and intervertebral disc degeneration (IDD).

OBJECTIVE

The aim of this study was to clarify the effects of estrogen deficiency and supplement on cervical IDD induced by bilateral facetectomy in rats.

SUMMARY OF BACKGROUND DATA

The relationship between IDD and OP is still controversy with the wide prevalence in aged people.

METHODS

Seventy-two Sprague-Dawley female rats were randomly divided into ovariectomy (OVX) group, facet joints resection of C4-6 (FR), FR-OVX group, estrogen replacement therapy (ERT, based on the FR-OVX group) group, and sham group. Specimens of C4-6 segment were harvested at 12 and 24 weeks. The microstructures of C5 vertebrae, vertebral endplate lesions and calcification, and IDD of C5/6 disc were evaluated by micro-computed tomography (micro-CT) and histology. The protein and gene levels of aggrecan, Col2α1, matrix metalloprotease (MMP)-3, and MMP-13 in the C5/6 and C4/5 discs were measured.

RESULTS

Microstructures of C5 vertebral body were weakened significantly after ovariectomy, while restored effectively with estradiol supplementation. The facetectomy led to significant IDD, and the IDD was aggravated when combined with OVX. The IDD of the ERT group was alleviated effectively and similar to that of the FR group in intervertebral disc height, vertebral endplate lesions and calcification, and disc degeneration scores. In addition, the estrogen supplement maintained the extracellular matrix by decreasing MMP-3 and MMP-13, and increasing aggrecan and Col2α1 expression.

CONCLUSION

The present study demonstrated that estrogen deficiency exacerbated IDD induced by spinal instability, while estrogen supplementation alleviated the progression of disc degeneration related to osteoporosis.

LEVEL OF EVIDENCE

N/A.

Validation of an In Vivo Medical Image-Based Young Human Lumbar Spine Finite Element Model

Mathematical models of the human spine can be used to investigate spinal biomechanics without the difficulties, limitations, and ethical concerns associated with physical experimentation. Validation of such models is necessary to ensure that the modeled system behavior accurately represents the physics of the actual system. The goal of this work was to validate a medical image-based nonlinear lumbosacral spine finite element model of a healthy 20-yr-old female subject under physiological moments. Range of motion (ROM), facet joint forces (FJF), and intradiscal pressure (IDP) were compared with experimental values and validated finite element models from the literature. The finite element model presented in this work was in good agreement with published experimental studies and finite element models under pure moments. For applied moments of 7.5 N·m, the ROM in flexion–extension, axial rotation, and lateral bending were 39 deg, 16 deg, and 28 deg, respectively. Excellent agreement was observed between the finite element model and experimental data for IDP under pure compressive loading. The predicted FJFs were lower than those of the experimental results and validated finite element models for extension and torsion, likely due to the nondegenerate properties chosen for the intervertebral disks and morphology of the young female spine. This work is the first to validate a computational lumbar spine model of a young female subject. This model will serve as a valuable tool for predicting orthopedic spinal injuries, studying the effect of intervertebral disk replacements using advanced biomaterials, and investigating soft tissue degeneration.

Clinical and Radiological Outcomes of Foraminal Decompression Using Unilateral Biportal Endoscopic Spine Surgery for Lumbar Foraminal Stenosis

Background

Since open Wiltse approach allows limited visualization for foraminal stenosis leading to an incomplete decompression, we report the short-term clinical and radiological results of unilateral biportal endoscopic foraminal decompression using 0° or 30° endoscopy with better visualization.

Methods

We examined 31 patients that underwent surgery for neurological symptoms due to lumbar foraminal stenosis which was refractory to 6 weeks of conservative treatment. All 31 patients underwent unilateral biportal endoscopic far-lateral decompression (UBEFLD). One portal was used for viewing purpose, and the other was for surgical instruments. Unilateral foraminotomy was performed under guidance of 0° or 30° endoscopy. Clinical outcomes were analyzed using the modified Macnab criteria, Oswestry disability index, and visual analogue scale. Plain radiographs obtained preoperatively and 1 year postoperatively were compared to analyze the intervertebral angle (IVA), dynamic IVA, percentage of slip, dynamic percentage of slip (gap between the percentage of slip on flexion and extension views), slip angle, disc height index (DHI), and foraminal height index (FHI).

Results

The IVA significantly increased from 6.24° ± 4.27° to 6.96° ± 3.58° at 1 year postoperatively (p = 0.306). The dynamic IVA slightly decreased from 6.27° ± 3.12° to 6.04° ± 2.41°, but the difference was not statistically significant (p = 0.375). The percentage of slip was 3.41% ± 5.24% preoperatively and 6.01% ± 1.43% at 1-year follow-up (p = 0.227), showing no significant difference. The preoperative dynamic percentage of slip was 2.90% ± 3.37%; at 1 year postoperatively, it was 3.13% ± 4.11% (p = 0.720), showing no significant difference. The DHI changed from 34.78% ± 9.54% preoperatively to 35.05% ± 8.83% postoperatively, which was not statistically significant (p = 0.837). In addition, the FHI slightly decreased from 55.15% ± 9.45% preoperatively to 54.56% ± 9.86% postoperatively, but the results were not statistically significant (p = 0.705).

Conclusions

UBEFLD using endoscopy showed a satisfactory clinical outcome after 1-year follow-up and did not induce postoperative segmental spinal instability. It could be a feasible alternative to conventional open decompression or fusion surgery for lumbar foraminal stenosis.

Effect of Asymmetric Tension on Biomechanics and Metabolism of Vertebral Epiphyseal Plate in a Rodent Model of Scoliosis

OBJECTIVE

To investigate the effect of asymmetric tension on idiopathic scoliosis (IS) and to understand its pathogenic mechanism.

METHODS

The rodent model of scoliosis was established using Sprague-Dawley rats with left rib-tethering from T₆ to T₁₂ , tail and shoulder amputation, and high-cage feeding. Vertebrae epiphyseal cartilage plates were harvested from the convex and concave sides. To analyze differences on the convex and concave sides, finite element analysis was carried out to determine the mechanical stress. Protein expression on epiphyseal cartilage was evaluated by western blot. Micro-CT was taken to evaluate the bone quality of vertebral on both sides.

RESULTS

Scoliosis curves presented in X-ray radiographs of the rats. Finite element analysis was carried out on the axial and transverse tension of the spine. Stresses of the convex side were -170.14, -373.18, and -3832.32 MPa (X, Y, and Z axis, respectively), while the concave side showed stresses of 361.99, 605.55, and 3661.95 MPa. Collagen type II, collagen type X, Sox 9, RunX2, VEGF, and aggrecan were expressed significantly more on the convex side (P < 0.05). There was asymmetric expression of protein on the epiphyseal cartilage plate at molecular level. Compared with the convex side, the concave side had significantly lower value in the BV/TV and Tb.N, but higher value in the Tb.Sp (P < 0.05). There was asymmetry of bone quality in micro-architecture.

CONCLUSIONS

In this study, asymmetric tension contributed to asymmetry in protein expression and bone quality on vertebral epiphyseal plates, ultimately resulting in asymmetry of anatomy. In addition, asymmetry of anatomy aggravated asymmetric tension. It is the first study to show that there is an asymmetrical vicious circle in IS.

Extraforaminal approach of biportal endoscopic spinal surgery: a new endoscopic technique for transforaminal decompression and discectomy

This study was performed to describe the extraforaminal approach of biportal endoscopic spinal surgery (BESS) as a new endoscopic technique for transforaminal decompression and discectomy and to demonstrate the clinical outcomes of this new procedure for the first time. Twenty-one patients (27 segments) who underwent the extraforaminal approach of BESS between March 2015 and April 2016 were enrolled according to the inclusion and exclusion criteria. The operative time (minutes/level) and complications after the procedure were recorded. The visual analog scale (VAS) score was checked to assess the degree of radicular leg pain preoperatively and at the time of the last follow-up. The modified Macnab criteria were used to examine the clinical outcomes at the time of the last follow-up. The mean duration of the follow-up period was 14.8 months (minimum duration 12 months). The mean operative time was 96.7 minutes for one level. The mean VAS score for radicular leg pain dropped from a preoperative score of 7.5 ± 0.9 to a final follow-up score of 2.5 ± 1.2 (p < 0.001). The final outcome according to the modified Macnab criteria was excellent in 5 patients (23.8%), good in 12 (57.2%), fair in 4 (19.0%), and poor in 0. Therefore, excellent or good results (a satisfied outcome) were obtained in 80.9% of the patients. Complications were limited to one dural tear (4.8%). The authors found that the extraforaminal approach of BESS was a feasible and advantageous endoscopic technique for the treatment of foraminal lesions, including stenosis and disc herniation. They suggest that this technique represents a useful, alternative, minimally invasive method that can be used to treat lumbar foraminal stenosis and disc herniation.

Biomechanical testing of a PEEK-based dynamic instrumentation device in a lumbar spine model

BACKGROUND

The purpose of this study was to investigate the range-of-motion after posterior polyetheretherketone-based rod stabilisation combined with a dynamic silicone hinge in order to compare it with titanium rigid stabilisation.

METHODS

Five human cadaveric lumbar spines with four vertebra each (L2 to L5) were tested in a temperature adjustable spine-testing set-up in four trials: (1) native measurement; (2) kinematics after rigid monosegmental titanium rod instrumentation with anterior intervertebral bracing of the segment L4/5; (3) kinematics after hybrid posterior polyetheretherketone rod instrumentation combined with a silicone hinge within the adjacent level (L3/4) and (4) kinematics after additional decompression with laminectomy of L4 and bilateral resection of the inferior articular processes (L3). During all steps, the specimens were loaded quasi-statically with 1°/s with pure moment up to 7.5Nm in flexion/extension, lateral bending and axial rotation.

FINDINGS

In comparison to the native cadaveric spine, both the titanium device and polyetheretherketone-based device reduce the range-of-motion within the level L4/5 significantly (flexion/extension: reduction of 77%, p<0.001; lateral bending: reduction of 62%, p<0.001; axial rotation: reduction of 71%, p<0.001). There was a clear stabilisation effect after hybrid-instrumentation within the level L3/4, especially in flexion/extension (64%, p<0.001) and lateral bending (62%, p<0.001) but without any effect on the axial rotation. Any temperature dependency has not been observed.

INTERPRETATION

Surprisingly, the hybrid device compensates for laminectomy L4 and destabilising procedure within the level L3/4 in comparison to other implants. Further studies must be performed to show its effectiveness regarding the adjacent segment instability.

Effect of Graded Facetectomy on Lumbar Biomechanics

Facetectomy is an important intervention for spinal stenosis but may lead to spinal instability. Biomechanical knowledge for facetectomy can be beneficial when deciding whether fusion is necessary. Therefore, the aim of this study was to investigate the biomechanical effect of different grades of facetectomy. A three-dimensional nonlinear finite element model of L3–L5 was constructed. The mobility of the model and the intradiscal pressure (IDP) of L4-L5 for standing were inside the data from the literature. The effect of graded facetectomy on intervertebral rotation, IDP, facet joint forces, and maximum von Mises equivalent stresses in the annuli was analyzed under flexion, extension, left/right lateral bending, and left/right axial rotation. Compared with the intact model, under extension, unilateral facetectomy increased the range of intervertebral rotation (IVR) by 11.7% and IDP by 10.7%, while the bilateral facetectomy increased IVR by 40.7% and IDP by 23.6%. Under axial rotation, the unilateral facetectomy and the bilateral facetectomy increased the IVR by 101.3% and 354.3%, respectively, when turned to the right and by 1.1% and 265.3%, respectively, when turned to the left. The results conclude that, after unilateral and bilateral facetectomy, care must be taken when placing the spine into extension and axial rotation posture from the biomechanical point of view.

Material properties of bovine intervertebral discs across strain rates

The intervertebral disc (IVD) is a complex structure responsible for distributing compressive loading to adjacent vertebrae and allowing the vertebral column to bend and twist. To study the mechanical behaviour of individual components of the IVD, it is common for specimens to be dissected away from their surrounding tissues for mechanical testing. However, disrupting the continuity of the IVD to obtain material properties of each component separately may result in erroneous values. In this study, an inverse finite element (FE) modelling optimisation algorithm has been used to obtain material properties of the IVD across strain rates, therefore bypassing the need to harvest individual samples of each component. Uniaxial compression was applied to ten fresh-frozen bovine intervertebral discs at strain rates of 10^-3-1/s. The experimental data were fed into the inverse FE optimisation algorithm and each experiment was simulated using the subject specific FE model of the respective specimen. A sensitivity analysis revealed that the IVD's response was most dependent upon the Young's modulus (YM) of the fibre bundles and therefore this was chosen to be the parameter to optimise. Based on the obtained YM values for each test corresponding to a different strain rate (ε̇), the following relationship was derived:YM=35.5lnε̇+527.5. These properties can be used in finite element models of the IVD that aim to simulate spinal biomechanics across loading rates.