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

Role of intra-lamellar collagen and hyaluronan nanostructures in annulus fibrosus on lumbar spine biomechanics: insights from molecular mechanics-finite element-based multiscale analyses

Annulus fibrosus' (AF) ability to transmit multi-directional spinal motion is contributed by a combination of chemical interactions among biomolecular constituents-collagen type I (COL-I), collagen type II (COL-II), and proteoglycans (aggrecan and hyaluronan)-and mechanical interactions at multiple length scales. However, the mechanistic role of such interactions on spinal motion is unclear. The present work employs a molecular mechanics-finite element (FE) multiscale approach to investigate the mechanistic role of molecular-scale collagen and hyaluronan nanostructures in AF, on spinal motion. For this, an FE model of the lumbar segment is developed wherein a multiscale model of AF collagen fiber, developed from COL-I, COL-II, and hyaluronan using the molecular dynamics-cohesive finite element multiscale method, is incorporated. Analyses show AF collagen fibers primarily contribute to axial rotation (AR) motion, owing to angle-ply orientation. Maximum fiber strain values of 2.45% in AR, observed at the outer annulus, are 25% lower than the reported values. This indicates native collagen fibers are softer, attributed to the softer non-fibrillar matrix and higher interfibrillar sliding. Additionally, elastic zone stiffness of 8.61 Nm/° is observed to be 20% higher than the reported range, suggesting native AF lamellae exhibit lower stiffness, resulting from inter-collagen fiber bundle sliding. The presented study has further implications towards the hierarchy-driven designing of AF-substitute materials.

Biomechanical changes of oblique lumbar interbody fusion with different fixation techniques in degenerative spondylolisthesis lumbar spine: a finite element analysis

OBJECTIVE

There is a dearth of comprehensive research on the stability of the spinal biomechanical structure when combining Oblique Lumbar Interbody Fusion (OLIF) with internal fixation methods. Hence, we have devised this experiment to meticulously examine and analyze the biomechanical changes that arise from combining OLIF surgery with different internal fixation techniques in patients diagnosed with degenerative lumbar spondylolisthesis.

METHODS

Seven validated finite element models were reconstructed based on computed tomography scan images of the L3-L5 segment. These models included the intact model, a stand-alone (S-A) OLIF model, a lateral screw rod (LSR) OLIF model, a bilateral pedicle screw (BPS) OLIF model, an unilateral pedicle screw (UPS) OLIF model, a bilateral CBT (BCBT) OLIF model, and an unilateral CBT(UCBT) OLIF model. The range of motion (ROM), as well as stress levels in the cage, L4 lower endplate, L5 upper endplate, and fixation constructs were assessed across these different model configurations.

RESULTS

S-A model had the highest average ROM of six motion modes, followed by LSR, UPS, UCBT, BPS and BCBT. The BCBT model had a relatively lower cage stress than the others. The maximum peak von Mises stress of the fixation constructs was found in the LSR model. The maximum peak von Mises stress of L4 lower endplate was found in the S-A model. The peak von Mises stress on the L4 lower endplate of the rest surgical models showed no significant difference. The maximum peak von Mises stress of the L5 upper endplate was found in the S-A model. The minimum peak von Mises stress of the L5 upper endplate was found in the BCBT model. No significant difference was found for the peak von Mises stress of the L5 upper endplate among LSR, BPS, UPS and UCBT models.

CONCLUSION

Among the six different fixation techniques, BCBT exhibited superior biomechanical stability and minimal stress on the cage-endplate interface. It was followed by BPS, UCBT, UPS, and LSR in terms of effectiveness. Conversely, S-A OLIF demonstrated the least stability and resulted in increased stress on both the cage and endplates. Combining OLIF with BCBT fixation technique enhanced biomechanical stability compared to BPS and presented as a less invasive alternative treatment for patients with degenerative lumbar spondylolisthesis.

Anterior percutaneous full-endoscopic transcorporeal decompression for cervical disc herniation: a finite element analysis and long-term follow-up study

OBJECTIVE

The purpose of this study was to investigate the long-term consequences on the cervical spine after Anterior transcorporeal percutaneous endoscopy cervical discectomy (ATc-PECD) from the biomechanical standpoint.

METHODS

A three-dimensional model of the normal cervical spine C2-T1 was established using finite element method. Subsequently, a disc degeneration model and degeneration with surgery model were constructed on the basis of the normal model. The same loading conditions were applied to simulate flexion, extension, lateral bending and axial rotation of the cervical spine. We calculated the cervical range of motion (ROM), intradiscal pressure, and intravertebral body pressure under different motions for observing changes in cervical spine biomechanics after surgery. At the same time, we combined the results of a long-term follow-up of the ATc-PECD, and used imaging methods to measure vertebral and disc height and cervical mobility, the Japanese Orthopaedic Association (JOA) score and visual analog scale (VAS) score were used to assess pain relief and neurological functional recovery.

RESULTS

The long-term follow-up results revealed that preoperative JOA score, neck VAS score, hand VAS score, IDH, VBH, and ROM for patients were 9.49 ± 2.16, 6.34 ± 1.68, 5.14 ± 1.48, 5.95 ± 0.22 mm, 15.41 ± 1.68 mm, and 52.46 ± 9.36° respectively. It changed to 15.71 ± 1.13 (P < 0.05), 1.02 ± 0.82 (P < 0.05), 0.77 ± 0.76 (P < 0.05), 4.73 ± 0.26 mm (P < 0.05), 13.67 ± 1.48 mm (P < 0.05), and 59.26 ± 6.72° (P < 0.05), respectively, at 6 years postoperatively. Finite element analysis showed that after establishing the cervical spondylosis model, the overall motion range for flexion, extension, lateral bending, and rotation decreased by 3.298°, 0.753°, 3.852°, and 1.131° respectively. Conversely, after establishing the bone tunnel model, the motion range for these actions increased by 0.843°, 0.65°, 0.278°, and 0.488° respectively, consistent with the follow-up results. Moreover, analysis of segmental motion changes revealed that the increased cervical spine mobility was primarily contributed by the surgical model segments. Additionally, the finite element model demonstrated that bone tunneling could lead to increased stress within the vertebral bodies and intervertebral discs of the surgical segments.

CONCLUSIONS

Long-term follow-up studies have shown that ATc-PECD has good clinical efficacy and that ATc-PECD can be used as a complementary method for CDH treatment. The FEM demonstrated that ATc-PECD can lead to increased internal stresses in the vertebral body and intervertebral discs of the operated segments, which is directly related to cervical spine degeneration after ATc-PECD.

Traction Spurs in the Lumbar Spine: A Historical Overview and Future Perspectives

Numerous studies have explored the connection between lumbar osteophytes, their pathophysiology, and instability since Macnab's 1971 report on traction spurs as an indicator of lumbar instability. This study provides a narrative historical overview of traction spurs, a classic finding that suggests lumbar instability. It summarizes the causes of anterior lumbar vertebral osteophytes, the relationship between traction spurs and lumbar spinal instability, and the clinical significance of traction spurs. Vertebral osteophytes are grouped into two categories, namely, traction spurs or claw spurs, which represent different stages of the same pathological process. Traction spurs are indicative of instability and occur in the early stage of disc degeneration, characterized by temporary dysfunction or instability. Traction spur formation following fusion surgery can predict union or nonunion, and it serves as an indicator of preoperative and postoperative segmental instability. The relationship between traction spurs and radiographic instability, as well as their association with imaging findings such as CT and MRI, has been clarified. Additionally, finite element analysis and mechanical testing have been used to investigate the significance of traction spurs. However, further research is needed to verify that traction spurs are an accurate indicator of pre- and postoperative lumbar instability.

Effect of compressive and tensile forces on glucose concentration and cell viability within the intervertebral disc: A finite element study

Understanding the role of mechanical force on tissue nutrient transport is essential, as sustained force may affect nutrient levels within the disc and initiate disc degeneration. This study aims to evaluate the time-dependent effects of different compressive force amplitudes as well as tensile force on glucose concentration and cell viability within the disc. Based on the mechano-electrochemical mixture theory, a multiphasic finite element model of the lumbar intervertebral disc was developed. The minimum glucose concentration and minimum cell density in both normal and degenerated discs were predicted for different compressive force amplitudes, tensile force, and corresponding creep time. Under high compressive force, the minimum glucose concentration exhibited an increasing and then decreasing trend with creep time in the normal disc, whereas that of the degenerated disc increased, then decreased, and finally increased again. At steady state, a higher compressive force was accompanied by a lower glucose concentration distribution. In the degenerated disc, the minimum cell density was negatively correlated with creep time, with a greater range of affected tissue under a higher compressive force. For tensile force, the minimum glucose concentration of the degenerated disc raised over time. This study highlighted the importance of creep time, force magnitude, and force type in affecting nutrient concentration and cell viability. Sustained weight-bearing activities could deteriorate the nutrient environment of the degenerated disc, while tensile force might have a nonnegligible role in effectively improving nutrient levels within the degenerated disc.

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.

Biomechanical comparison of polyetheretherketone rods and titanium alloy rods in transforaminal lumbar interbody fusion: a finite element analysis

BACKGROUND

Whether polyetheretherketone (PEEK) rods have potential as an alternative to titanium alloy (Ti) rods in transforaminal lumbar interbody fusion (TLIF) remains unclear, especially in cases with insufficient anterior support due to the absence of a cage. The purpose of this study was to investigate biomechanical differences between PEEK rods and Ti rods in TLIF with and without a cage.

METHODS

An intact L1-L5 lumbar finite element model was constructed and validated. Accordingly, four TLIF models were developed: (1) Ti rods with a cage; (2) PEEK rods with a cage; (3) Ti rods without a cage; and (4) PEEK rods without a cage. The biomechanical properties were then compared among the four TLIF constructs.

RESULTS

With or without a cage, no obvious differences were found in the effect of PEEK rods and Ti rods on the range of motion, adjacent disc stress, and adjacent facet joint force. Compared to Ti rods, PEEK rods increase the average bone graft strain (270.8-6055.2 µE vs. 319.0-8751.6 µE). Moreover, PEEK rods reduced the stresses on the screw-rod system (23.1-96.0 MPa vs. 7.2-48.4 MPa) but increased the stresses on the cage (4.6-35.2 MPa vs. 5.6-40.9 MPa) and endplates (5.7-32.5 MPa vs. 6.6-37.6 MPa).

CONCLUSIONS

Regardless of whether a cage was used for TLIF, PEEK rods theoretically have the potential to serve as an alternative to Ti rods because they may provide certain stability, increase the bone graft strain, and reduce the posterior instrumentation stress, which might promote bony fusion and decrease instrumentation failure.

Stepwise reduction of bony density in patients induces a higher risk of annular tears by deteriorating the local biomechanical environment

BACKGROUND CONTEXT

The relationship between osteoporosis and intervertebral disc degeneration (IDD) remains unclear. Considering that annular tear is the primary phenotype of IDD in the lumbar spine, the deteriorating local biomechanical environment may be the main trigger for annular tears.

PURPOSE

To investigate whether poor bone mineral density (BMD) in the vertebral bodies may increase the risk of annular tears via the degradation of the local biomechanical environment.

STUDY DESIGN

This study was a retrospective investigation with relevant numerical mechanical simulations.

PATIENT SAMPLE

A total of 64 patients with low back pain (LBP) and the most severe IDD in the L4-L5 motion segment were enrolled.

OUTCOME MEASURES

Annulus integration status was assessed using diffusion tensor fibre tractography (DTT). Hounsfield unit (HU) values of adjacent vertebral bodies were employed to determine BMD. Numerical simulations were conducted to compute stress values in the annulus of models with different BMDs and body positions.

METHODS

The clinical data of the 64 patients with low back pain were collected retrospectively. The BMD of the vertebral bodies was measured using the HU values, and the annulus integration status was determined according to DTT. The data of the patients with and without annular tears were compared, and regression analysis was used to identify the independent risk factors for annular tears. Furthermore, finite element models of the L4-L5 motion segment were constructed and validated, followed by estimating the maximum stress on the post and postlateral interfaces between the superior and inferior bony endplates (BEPs) and the annulus.

RESULTS

Patients with lower HU values in their vertebral bodies had significantly higher incidence rates of annular tears, with decreased HU values being an independent risk factor for annular tears. Moreover, increased stress on the BEP-annulus interfaces was associated with a stepwise reduction of bony density (ie, elastic modulus) in the numerical models.

CONCLUSIONS

The stepwise reduction of bony density in patients results in a higher risk of annular tears by deteriorating the local biomechanical environment. Thus, osteoporosis should be considered to be a potential risk factor for IDD biomechanically.

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

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.

IVD fibrosis and disc collapse comprehensively aggravate vertebral body disuse osteoporosis and zygapophyseal joint osteoarthritis by posteriorly shifting the load transmission pattern

BACKGROUND

Disuse is a typical phenotype of osteoporosis, but the underlying mechanism has yet to be identified in elderly patients. Disc collapse and intervertebral disc (IVD) fibrosis are two main pathological changes in IVD degeneration (IDD) progression, given that these changes affect load transmission patterns, which may lead to disuse osteoporosis of vertebral bodies and zygapophyseal joint (ZJ) osteoarthritis (ZJOA) biomechanically.

METHODS

Clinical data from 59 patients were collected retrospectively. Patient vertebral bony density, ZJOA grade, and disc collapse status were judged via CT. The IVD fibrosis grade was determined based on the FA measurements. Regression analyses identified potential independent risk factors for osteoporosis and ZJOA. L4-L5 numerical models with and without disc collapse and IVD fibrosis were constructed; stress distributions on the bony endplate (BEP) and zygapophyseal joint (ZJ) cartilages were computed in models with and without disc collapse and IVD fibrosis.

RESULTS

A significantly lower disc height ratio and significantly greater FA were recorded in patients with ZJOA. A significant correlation was observed between lower HU values and two parameters related to IDD progression. These factors were also proven to be independent risk factors for both osteoporosis and ZJOA. Correspondingly, compared to the intact model without IDD. Lower stress on vertebral bodies and greater stress on ZJOA can be simultaneously recorded in models of disc collapse and IVD fibrosis.

CONCLUSIONS

IVD fibrosis and disc collapse simultaneously aggravate vertebral body disuse osteoporosis and ZJOA by posteriorly shifting the load transmission pattern.

Biomedical analysis of four fixation systems in treatment of type II traumatic spondylolisthesis of the axis: a finite element analysis

This study aimed to compare the properties and safety of self-designed plates in type II traumatic spondylolisthesis of the axis with those of traditional devices via finite element (FE) analysis. We constructed a hangman's fracture FE model from the occipital bone (C0) level to the C3 level. Then, FE models were constructed for the following four fixation systems: an anterior cervical L-shaped plate with four vertebral screws (4-ACLP), or six screws (6-ACLP), an anterior cervical orion plate (ACOP), and a posterior fixation system. A preloaded compressive force of 50 N and a moment of 1.5 N·m were applied to each model under six working conditions. The mobility of the C2/3 segment decreased significantly in four fixation models. In the Mises stress cloud diagram, 4-ACLP showed a better stress distribution in both the bone graft and fixation system than 6-ACLP and ACOP. The resultant force of 4-ACLP was lower but higher than ACOP in axial force. Additionally, the cage in the 4-ACLP configuration experienced the highest stress in the six working conditions. Hence, this novel self-designed plate has the potential to mitigate the operational difficulties, provide sufficient stability, reduce the risk of plate or screw fractures, and improve bone fusion.

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.

Intraoperative capsule protection can reduce the potential risk of adjacent segment degeneration acceleration biomechanically: an in silico study

BACKGROUND

The capsule of the zygapophyseal joint plays an important role in motion segmental stability maintenance. Iatrogenic capsule injury is a common phenomenon in posterior approach lumbar interbody fusion operations, but whether this procedure will cause a higher risk of adjacent segment degeneration acceleration biomechanically has yet to be identified.

METHODS

Posterior lumbar interbody fusion (PLIF) with different grades of iatrogenic capsule injury was simulated in our calibrated and validated numerical model. By adjusting the cross-sectional area of the capsule, different grades of capsule injury were simulated. The stress distribution on the cranial motion segment was computed under different loading conditions to judge the potential risk of adjacent segment degeneration acceleration.

RESULTS

Compared to the PLIF model with an intact capsule, a stepwise increase in the stress value on the cranial motion segment can be observed with a step decrease in capsule cross-sectional areas. Moreover, compared to the difference between models with intact and slightly injured capsules, the difference in stress values was more evident between models with slight and severe iatrogenic capsule injury.

CONCLUSION

Intraoperative capsule protection can reduce the potential risk of adjacent segment degeneration acceleration biomechanically, and iatrogenic capsule damage on the cranial motion segment should be reduced to optimize patients' long-term prognosis.

Biomechanical evaluation of different oblique lumbar interbody fusion constructs: a finite element analysis

BACKGROUND

Finite element analysis (FEA) was performed to investigate the biomechanical differences between different adjunct fixation methods for oblique lumbar interbody fusion (OLIF) and to further analyze its effect on adjacent segmental degeneration.

METHODS

We built a single-segment (Si-segment) finite element model (FEM) for L4-5 and a double-segment (Do-segment) FEM for L3-5. Each complete FEM was supplemented and modified, and both developed two surgical models of OLIF with assisted internal fixation. They were OLIF with posterior bilateral percutaneous pedicle screw (TINA system) fixation (OLIF + BPS) and OLIF with lateral plate system (OLIF + LPS). The range of motion (ROM) and displacement of the vertebral body, cage stress, adjacent segment disc stress, and spinal ligament tension were recorded for the four models during flexion/extension, right/left bending, and right/left rotation by applying follower load.

RESULTS

For the BPS and LPS systems in the six postures of flexion, extension, right/left bending, and right/left rotation, the ROM of L4 in the Si-segment FEM were 0.32°/1.83°, 0.33°/1.34°, 0.23°/0.47°, 0.24°/0.45°, 0.33°/0.79°, and 0.34°/0.62°; the ROM of L4 in the Do-segment FEM were 0.39°/2.00°, 0.37°/1.38°, 0.23°/0.47°, 0.21°/0.44°, 0.33°/0.57°, and 0.31°/0.62°, and the ROM of L3 in the Do-segment FEM were 6.03°/7.31°, 2.52°/3.50°, 4.21°/4.38°, 4.21°/4.42°, 2.09°/2.32°, and 2.07°/2.43°. BPS system had less vertebral displacement, less cage maximum stress, and less spinal ligament tension in Si/Do-segment FEM relative to the LPS system. BPS system had a smaller upper adjacent vertebral ROM, greater intervertebral disc stress in terms of left and right bending as well as left and right rotation compared to the LPS system in the L3-4 of the Do-segment FEM. There was little biomechanical difference between the same fixation system in the Si/Do-segment FEM.

CONCLUSIONS

Our finite element analysis showed that compared to OLIF + LPS, OLIF + BPS (TINA) is more effective in reducing interbody stress and spinal ligament tension, and it better maintains the stability of the target segment and provides a better fusion environment to resist cage subsidence. However, OLIF + BPS (TINA) may be more likely to cause adjacent segment degeneration than OLIF + LPS.

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

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

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.

Preliminary exploration of the biomechanical properties of three novel cervical porous fusion cages using a finite element study

BACKGROUND

Porous cages are considered a promising alternative to high-density cages because their interconnectivity favours bony ingrowth and appropriate stiffness tuning reduces stress shielding and the risk of cage subsidence.

METHODS

This study proposes three approaches that combine macroscopic topology optimization and micropore design to establish three new types of porous cages by integrating lattices (gyroid, Schwarz, body-centred cubic) with the optimized cage frame. Using these three porous cages along with traditional high-density cages, four ACDF surgical models were developed to compare the mechanical properties of facet articular cartilage, discs, cortical bone, and cages under specific loads.

RESULTS

The facet joints in the porous cage groups had lower contact forces than those in the high-density cage group. The intervertebral discs in all models experienced maximum stress at the C5/6 segment. The stress distribution on the cortical bone surface was more uniform in the porous cage groups, leading to increased average stress values. The gyroid, Schwarz, and BCC cage groups showed higher average stress on the C5 cortical bone. The average stress on the surface of porous cages was higher than that on the surface of high-density cages, with the greatest difference observed under the lateral bending condition. The BCC cage demonstrated favourable mechanical stability.

CONCLUSION

The new porous cervical cages satifies requirements of low rigidity and serve as a favourable biological scaffold for bone ingrowth. This study provides valuable insights for the development of next-generation orthopaedic medical devices.

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 response of lumbar intervertebral disc in daily sitting postures: a poroelastic finite element analysis

This study aims to establish and validate a poroelastic L4-L5 finite element model to evaluate the effect of different sitting postures and their durations on the mechanical responses of the disc. During the sustained loading conditions, the height loss, fluid loss and von-Mises stress gradually increased, but the intradiscal pressure decreased. The varying rates of aforementioned parameters were more significant at the initial loading stage and less so at the end. The predicted values in the flexed sitting posture were significantly greater than other postures. The extended sitting posture caused an obvious von-Mises stress concentration in the posterior region of the inter-lamellar matrix. From the biomechanical perspective, prolonged sitting may pose a high risk of lumbar disc degeneration, and therefore adjusting the posture properly in the early stage of sitting time may be useful to mitigate that. Additionally, upright sitting is a safer posture, while flexed sitting posture is more harmful.

Topology optimization and dynamic characteristic evaluation of W-shaped interspinous process device

For reducing the adjacent segment degeneration of the lumbar spine, the interspinous process device as a kind of flexible non-fusion device was designed to overcome the deficiencies associated with rigid fusion devices. However, it was not clear how the interspinous process device influenced the human spine system, especially the lumbar spine under a vibration environment. This study was designed to evaluate the effect of StenoFix under the vibration condition and also to optimize the structure of the device to obtain better biomechanical performance. A finite element model of the intact lumbar spine was developed and validated. The surgical finite element model was constructed by implanting the interspinous process device StenoFix. Using topology optimization, a new device StenoFix-new was redesigned. The results showed that the interspinous process device decreased vibration amplitudes of annulus stress and intradiscal pressure under vibration at the surgical level. The redesigned StenoFix-new with the smaller stiffness exhibited a better dynamic flexibility performance than StenoFix. In addition, the range of motions of StenoFix-new was closer to the intact model than StenoFix at the surgical level. These results might encourage the designers to give more consideration to the dynamic characteristics of the human spine on the premise of ensuring the safety and strength of implanted devices.

Biomechanical Comparison of Different Surgical Approaches for the Treatment of Adjacent Segment Diseases after Primary Transforaminal Lumbar Interbody Fusion: A Finite Element Analysis

BACKGROUND AND OBJECTIVE

Adjacent segment disease (ASD) is a well-known complication after interbody fusion. Revision surgery is necessary for symptomatic ASD to further decompress and fix the affected segment. However, no optimal construct is accepted as a standard in treating ASD. The purpose of this study was to compare the biomechanical effects of different surgical approaches for the treatment of ASD after primary transforaminal lumbar interbody fusion (TLIF).

METHODS

A finite element model of the L1-S1 was conducted based on computed tomography scan images. The primary surgery model was developed with a single-level TLIF at L4-L5 segment. The revision surgical models were developed with anterior lumbar interbody fusion (ALIF), lateral lumbar interbody fusion (LLIF), or TLIF at L3-L4 segment. The range of motion (ROM), intradiscal pressure (IDP), and the stress in cages were compared to investigate the biomechanical influences of different surgical approaches.

RESULTS

The results indicated that all the three surgical approaches can stabilize the spinal segment by reducing the ROM at revision level. The ROM and IDP at adjacent segments of revision model of TLIF was greater than those of other revision models. While revision surgery with ALIF and LLIF had similar effects on the ROM and IDP of adjacent segments. Compared among all the surgical models, cage stress in revision model of TLIF was the maximum in extension and axial rotation.

CONCLUSION

The IDP at adjacent segments and stress in cages of revision model of TLIF was greater than those of ALIF and LLIF. This may be that direct extension of the surgical segment in the same direction results in stress concentration.

Annulus Calibration Increases the Computational Accuracy of the Lumbar Finite Element Model

STUDY DESIGN

Mechanical simulations.

OBJECTIVE

Inadequate calibration of annuli negatively affects the computational accuracy of finite element (FE) models. Specifically, the definition of annulus average radius (AR) does not have uniformity standards. Differences between the elastic moduli in the different layers and parts of the annulus were not fully calibrated when a linear elastic material is used to define its material properties. This study aims to optimize the computational accuracy of the FE model by calibrating the annulus.

METHODS

We calibrated the annulus AR and elastic modulus in our anterior-constructed lumbar model by eliminating the difference between the computed range of motion and that measured by in vitro studies under a flexion-extension loading condition. Multi-indicator validation was performed by comparing the computed indicators with those measured in in vitro studies. The computation time required for the different models has also been recorded to evaluate the computational efficiency.

RESULTS

The difference between computed and measured ROMs was less than 1% when the annulus AR and elastic modulus were calibrated. In the model validation process, all the indicators computed by the calibrated FE model were within ±1 standard deviation of the average values obtained from in vitro studies. The maximum difference between the computed and measured values was less than 10% under nearly all loading conditions. There is no apparent variation tendency for the computational time associated with different models.

CONCLUSION

The FE model with calibrated annulus AR and regional elastic modulus has higher computational accuracy and can be used in subsequent mechanical studies.

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

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

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.

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.

The influence of over-distraction on biomechanical response of cervical spine post anterior interbody fusion: a comprehensive finite element study

Introduction: Anterior cervical discectomy and fusion (ACDF) has been considered as the gold standard surgical treatment for cervical degenerative pathologies. Some surgeons tend to use larger-sized interbody cages during ACDF to restore the index intervertebral disc height, hence, this study evaluated the effect of larger-sized interbody cages on the cervical spine with ACDF under both static and cyclic loading. Method: Twenty pre-operative personalized poro-hyperelastic finite element (FE) models were developed. ACDF post-operative models were then constructed and four clinical scenarios (i.e., 1) No-distraction; 2) 1 mm distraction; 3) 2 mm distraction; and 4) 3 mm distraction) were predicted for each patient. The biomechanical responses at adjacent spinal levels were studied subject to static and cyclic loading. Non-parametric Friedman statistical comparative tests were performed and the p values less than 0.05 were reflected as significant. Results: The calculated intersegmental range of motion (ROM) and intradiscal pressure (IDP) from 20 pre-operative FE models were within the overall ranges compared to the available data from literature. Under static loading, greater ROM, IDP, facet joint force (FJF) values were detected post ACDF, as compared with pre-op. Over-distraction induced significantly higher IDP and FJF in both upper and lower adjacent levels in extension. Higher annulus fibrosus stress and strain values, and increased disc height and fluid loss at the adjacent levels were observed in ACDF group which significantly increased for over-distraction groups. Discussion: it was concluded that using larger-sized interbody cages (the height of ≥2 mm of the index disc height) can result in remarkable variations in biomechanical responses of adjacent levels, which may indicate as risk factor for adjacent segment disease. The results of this comprehensive FE investigation using personalized modeling technique highlight the importance of selecting the appropriate height of interbody cage in ACDF surgery.

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 different semi-rigid junctional fixation techniques using finite element analysis

BACKGROUND

Proximal junctional failure is a common complication attributed to the rigidity of long pedicle screw fixation constructs used for surgical correction of adult spinal deformity. Semi-rigid junctional fixation achieves a gradual transition in range of motion at the ends of spinal instrumentation, which could lead to reduced junctional stresses, and ultimately reduce the incidence of proximal junctional failure. This study investigates the biomechanical effect of different semi-rigid junctional fixation techniques in a T8-L3 finite element spine segment model.

METHODS

First, degeneration of the intervertebral disc was successfully implemented by altering the height. Second, transverse process hooks, one- and two-level clamped tapes, and one- and two-level knotted tapes instrumented proximally to three-level pedicle screw fixation were validated against ex vivo range of motion data of a previous study. Finally, the posterior ligament complex forces and nucleus pulposus stresses were quantified.

FINDINGS

Simulated range of motions demonstrated the fidelity of the general model and modelling of semi-rigid junctional fixation techniques. All semi-rigid junctional fixation techniques reduced the posterior ligament complex forces at the junctional zone compared to pedicle screw fixation. Transverse process hooks and knotted tapes reduced nucleus pulposus stresses, whereas clamped tapes increased nucleus pulposus stresses at the junctional zone.

INTERPRETATION

The relationship between the range of motion transition and the reductions in posterior ligament complex and nucleus pulposus stresses was complex and dependent on the fixation techniques. Clinical trials are required to compare the effectiveness of semi-rigid junctional fixation techniques in terms of reducing proximal junctional failure incidence rates.

Finite element analysis of optimized novel additively manufactured non-articulating prostheses for cervical total disc replacement

Ball-and-socket designs of cervical total disc replacement (TDR) have been popular in recent years despite the disadvantages of polyethylene wear, heterotrophic ossification, increased facet contact force, and implant subsidence. In this study, a non-articulating, additively manufactured hybrid TDR with an ultra-high molecular weight polyethylene core and polycarbonate urethane (PCU) fiber jacket, was designed to mimic the motion of normal discs. A finite element (FE) study was conducted to optimize the lattice structure and assess the biomechanical performance of this new generation TDR with an intact disc and a commercial ball-and-socket Baguera^®C TDR (Spineart SA, Geneva, Switzerland) on an intact C5-6 cervical spinal model. The lattice structure of the PCU fiber was constructed using the Tesseract or the Cross structures from the IntraLattice model in the Rhino software (McNeel North America, Seattle, WA) to create the hybrid I and hybrid II groups, respectively. The circumferential area of the PCU fiber was divided into three regions (anterior, lateral and posterior), and the cellular structures were adjusted. Optimal cellular distributions and structures were A2L5P2 in the hybrid I and A2L7P3 in the hybrid II groups. All but one of the maximum von Mises stresses were within the yield strength of the PCU material. The range of motions, facet joint stress, C6 vertebral superior endplate stress and path of instantaneous center of rotation of the hybrid I and II groups were closer to those of the intact group than those of the Baguera^®C group under 100 N follower load and pure moment of 1.5 Nm in four different planar motions. Restoration of normal cervical spinal kinematics and prevention of implant subsidence could be observed from the FE analysis results. Superior stress distribution in the PCU fiber and core in the hybrid II group revealed that the Cross lattice structure of a PCU fiber jacket could be a choice for a next-generation TDR. This promising outcome suggests the feasibility of implanting an additively manufactured multi-material artificial disc that allows for better physiological motion than the current ball-and-socket design.

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

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

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.

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.

A biomechanical comparison of posterior fixation approaches in lumbar fusion using computed tomography based lumbosacral spine modelling

Extreme lateral interbody fusion (XLIF) may be performed with a standalone interbody cage, or with the addition of unilateral or bilateral pedicle screws; however, decisions regarding supplemental fixation are predominantly based on clinical indicators. This study examines the impact of posterior supplemental fixation on facet micromotions, cage loads and load-patterns at adjacent levels in a L4-L5 XLIF at early and late fusion stages. CT data from an asymptomatic subject were segmented into anatomical regions and digitally stitched into a surface mesh of the lumbosacral spine (L1-S1). The interbody cage and posterior instrumentation (unilateral and bilateral) were inserted at L4-L5. The volumetric mesh was imported into finite element software for pre-processing, running nonlinear static solves and post-processing. Loads and micromotions at the index-level facets reduced commensurately with the extent of posterior fixation accompanying the XLIF, while load-pattern changes observed at adjacent facets may be anatomically dependent. In flexion at partial fusion, compressive stress on the cage reduced by 54% and 72% in unilateral and bilateral models respectively; in extension the reductions were 58% and 75% compared to standalone XLIF. A similar pattern was observed at full fusion. Unilateral fixation provided similar stability compared to bilateral, however there was a reduction in cage stress-risers with the bilateral instrumentation. No changes were found at adjacent discs. Posterior supplemental fixation alters biomechanics at the index and adjacent levels in a manner that warrants consideration alongside clinical information. Unilateral instrumentation is a more efficient option where the stability requirements and subsidence risk are not excessive.

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.

Will the adjustment of insertional pedicle screw positions affect the risk of adjacent segment diseases biomechanically? An in-silico study

Background

The fixation-induced biomechanical deterioration will increase the risk of adjacent segment diseases (ASD) after lumbar interbody fusion with Bilateral pedicle screw (BPS) fixation. The accurate adjustment of insertional pedicle screw positions is possible, and published studies have reported its mechanical effects. However, no studies clarified that adjusting insertional screw positions would affect the postoperative biomechanical environment and the risk of ASD. The objective of this study was to identify this issue and provide theoretical references for the optimization of insertional pedicle screw position selections.

Methods

The oblique lumbar interbody fusion fixed by BPS with different insertional positions has been simulated in the L4-L5 segment of our previously constructed and validated lumbosacral model. Biomechanical indicators related to ASD have been computed and recorded under flexion, extension, bending, and axial rotation loading conditions.

Results

The change of screw insertional positions has more apparent biomechanical effects on the cranial than the caudal segment. Positive collections can be observed between the reduction of the fixation length and the alleviation of motility compensation and stress concentration on facet cartilages. By contrast, no pronounced tendency of stress distribution on the intervertebral discs can be observed with the change of screw positions.

Conclusions

Reducing the fixation stiffness by adjusting the insertional screw positions could alleviate the biomechanical deterioration and be an effective method to reduce the risk of ASD caused by BPS.

Elastic Fibers in the Intervertebral Disc: From Form to Function and toward Regeneration

Despite extensive efforts over the past 40 years, there is still a significant gap in knowledge of the characteristics of elastic fibers in the intervertebral disc (IVD). More studies are required to clarify the potential contribution of elastic fibers to the IVD (healthy and diseased) function and recommend critical areas for future investigations. On the other hand, current IVD in-vitro models are not true reflections of the complex biological IVD tissue and the role of elastic fibers has often been ignored in developing relevant tissue-engineered scaffolds and realistic computational models. This has affected the progress of IVD studies (tissue engineering solutions, biomechanics, fundamental biology) and translation into clinical practice. Motivated by the current gap, the current review paper presents a comprehensive study (from the early 1980s to 2022) that explores the current understanding of structural (multi-scale hierarchy), biological (development and aging, elastin content, and cell-fiber interaction), and biomechanical properties of the IVD elastic fibers, and provides new insights into future investigations in this domain.

Effect of degenerative factors on cervical spinal cord during flexion and extension: a dynamic finite element analysis

Spinal cord injury (SCI) is a global problem that brings a heavy burden to both patients and society. Recent investigations indicated degenerative disease is taking an increasing part in SCI with the growth of the aging population. However, little insight has been gained about the effect of cervical degenerative disease on the spinal cord during dynamic activities. In this work, a dynamic fluid-structure interaction model was developed and validated to investigate the effect of anterior and posterior encroachment caused by degenerative disease on the spinal cord during normal extension and flexion. Maximum von-Mises stress and maximum principal strain were observed at the end of extension and flexion. The abnormal stress distribution caused by degenerative factors was concentrated in the descending tracts of the spinal cord. Our finding indicates that the excessive motion of the cervical spine could potentially exacerbate spinal cord injury and enlarge injury areas. Stress and strain remained low compared to extension during moderate flexion. This suggests that patients with cervical degenerative disease should avoid frequent or excessive flexion and extension which could result in motor function impairment, whereas moderate flexion is safe. Besides, encroachment caused by degenerative factors that are not significant in static imaging could also cause cord compression during normal activities.

Biomechanical Evaluation of Spinal Column after Percutaneous Cement Discoplasty: A Finite Element Analysis

Objective

To compare the biomechanical properties of percutaneous cement discoplasty (PCD) in the spinal column between different implant‐endplate friction.

Methods

A validated L3‐Scarumfinite element (FE) model was modified for simulation. In the PCD model, the L4/5 level was modified based on model 1 (M1) and model 2 (M2). In M1, the interaction between bone cement and endplate was defined as face‐to‐face contact with a friction coefficient of 0.3; in M2, the contact was defined as a Tie constraint. 7.5 N m moments of four physiological motions and axial load of 15, 100 and 400 N preload were imposed at the top of L3. The range of motion (ROM) and interface stress analysis of endplates, annulus fibrosus and bone cement of the operated level were calculated for comparisons among the three models.

Results

The ROM of M1 and M2 increased when compared with the intact model during flexion (FL) (17.5% vs 10.0%), extension (EX) (8.8% vs −8.8%), left bending (LB) (19.0% vs −17.2%) and left axial rotation (LR) (34.6% vs −3.8%). The stress of annulus fibrosus in M1 and M2 decreased in FL (−48.4% vs −57.5%), EX (−25.7% vs −14.7%), LB (−47.5% vs −52.4%), LR (−61.4% vs −68.7%) and axis loading of 100 N (−41.5% vs −15.3%), and 400 N (−27.9% vs −27.3%). The stress of upper endplate of M1 and M2 increased in FL (24.6% vs 24.7%), LB (82.2% vs 89.5%), LR (119% vs 62.4%) and axis loading of 100 N (64.6% vs 45.5%), and 400 N (58.2% vs 24.3%), but was similar in EX (2.9% vs 0.3%). The stress of lower endplate of M1 and M2 increased in FL (170.9% vs 175.0%), EX (180.8% vs 207.7%), LB (302.6% vs 274.7%), LR (332.4% vs 132.8%) and axis loading of 100 N (350.7% vs 168.6%), and 400 N (165.2% vs 106.7%).

Conclusion

Percutaneous cement discoplasty procedure could make effect on the mobility or stiffness. The fusion of bone cement and endplate might have more biomechanical advantages, including of the decreasing rate of implant subsidence and dislocation, and the increase spine stability.

From Mechanobiology to Mechanical Repair Strategies: A Bibliometric Analysis of Biomechanical Studies of Intervertebral Discs

Neck pain and low back pain are major challenges in public health, and intervertebral disc (IVD) biomechanics is an important multidisciplinary field. To date, no bibliometric literature review of the relevant literature has been performed, so we explored the emerging trends, landmark studies, and major contributors to IVD biomechanics research. We searched the Web of Science core collection (1900–2022) using keywords mainly composed of “biomechanics” and “intervertebral disc” to conduct a bibliometric analysis of original papers and their references, focusing on citations, authors, journals, and countries/regions. A co-citation analysis and clustering of the references were also completed. A total of 3189 records met the inclusion criteria. In the co-citation network, cluster #0, labeled as “annulus fibrosus tissue engineering”, and cluster #1, labeled as “micromechanical environment”, were the biggest clusters. References by MacLean et al and Holzapfel et al were positioned exactly between them and had high betweenness centrality. There existed a research topic evolution between mechanobiology and mechanical repair strategies of IVDs, and the latter had been identified as an emerging trend in IVD biomechanics. Numerous landmark studies had contributed to several fields, including mechanical testing of normal and pathological IVDs, mechanical evaluation of new repair strategies and development of finite element model. Adams MA was the author most cited by IVD biomechanics papers. Spine, the European Spine Journal, and the Journal of Biomechanics were the three journals where the most original articles and their references have been published. The United States has contributed most to the literature (n = 1277 papers); however, the research output of China is increasing. In conclusion, the present study suggests that IVD repair is an emerging trend in IVD biomechanics.

A finite element analysis on different bone cement forms and injection volumes injected into lumbar vertebral body in percutaneous kyphoplasty

Background

To investigate the stress changes between different bone cement forms and injection volumes in adjacent vertebrae after percutaneous kyphoplasty (PKP) by establishing a three-dimensional finite element model of osteoporosis.

Methods

A male healthy volunteer was selected. CT of scans L1 to L3 vertebrae were imported into Mimics 21.0 software.The vertebral model of osteoporosiswas established based on previous literature reference. The models were divided into three groups: unilateral, bilateral integration and bilateral separation groups, with each group injecting 2 ml, 4,ml and 6 ml of bone cement, respectively. In all models, a vertical compressive load of 500 N, anterior flexion/posterior extension, left/right bending, and left/right rotation were applied with a moment of 7.5 N/m, of which 85% was applied to the anterior mid-column and 15% to the posterior column. The stress changes between adjacent vertebrae under different conditions were calculated.

Results

After percutaneous kyphoplasty was applied to the L2 vertebral body, some differences can be found between the effects of different cement injection volumes and cement morphology on adjacent structures. There was no major difference between the groups when the bone cement injection volume was 2 ml. When the amount of bone cement injected was 4 ml, the bone cement morphology of the bilateral integration group (BIG) produced less stress between adjacent vertebral bodies. The minimum stress was 14.95 MPa in the L3 vertebral body in posterior extension. Whereas the stress levels on adjacent intervertebral structures, BIG shaped bone cement shows some superiority. In addition, the adjacent vertebrae and intervertebral structures are subjected to less stress during left and right rotation.

Conclusions

The present finite element study suggested that bilateral integration bone cement is a suitable form of cement injection, and when the injection volume is 4 ml, reduces stress on adjacent segments by approximately 15% while maintaining the stability of the injected vertebral body.

Effect of the In Situ Screw Implantation Region and Angle on the Stability of Lateral Lumbar Interbody Fusion: A Finite Element Study

Objective

To investigate the effect of the in situ screw implantation region and angle on the stability of lateral lumbar interbody fusion (LLIF) from a biomechanical perspective.

Methods

A validated L2‐4 finite element (FE) model was modified for simulation. The L3‐4 fused segment undergoing LLIF surgery was modeled. The area between the superior and inferior edges and the anterior and posterior edges of the vertebral body (VB) is divided into four zones by three parallel lines in coronal and horizontal planes. In situ screw implantation methods with different angles based on the three parallel lines in coronal plane were applied in Models A, B, and C (A: parallel to inferior line; B: from inferior line to midline; C: from inferior line to superior line). In addition, four implantation methods with different regions based on the three parallel lines in horizontal plane were simulated as types 1–2, 1–3, 2–2, and 2–3 (1–2: from anterior line to midline; 1–3: from anterior line to posterior line; 2–2: parallel to midline; 2–3: from midline to posterior line). L3‐4 ROM, interbody cage stress, screw‐bone interface stress, and L4 superior endplate stress were tracked and calculated for comparisons among these models.

Results

The L3‐4 ROM of Models A, B, and C decreased with the extent ranging from 47.9% (flexion‐extension) to 62.4% (lateral bending) with no significant differences under any loading condition. Types 2–2 and 2–3 had 45% restriction, while types 1–2 and 1–3 had 51% restriction in ROM under flexion‐extension conditions. Under lateral bending, types 2–2 and 2–3 had 70.6% restriction, while types 1–2 and 1–3 had 61.2% restriction in ROM. Under axial rotation, types 2–2 and 2–3 had 65.2% restriction, while types 1–2 and 1–3 had 59.3% restriction in ROM. The stress of the cage in types 2–2 and 2–3 was approximately 20% lower than that in types 1–2 and 1–3 under all loading conditions in all models. The peak stresses at the screw‐bone interface in types 2–2 and 2–3 were much lower (approximately 35%) than those in types 1–2 and 1–3 under lateral bending, while no significant differences were observed under flexion‐extension and axial rotation. The peak stress on the L4 superior endplate was approximately 30 MPa and was not significantly different in all models under any loading condition.

Conclusions

Different regions of entry‐exit screws induced multiple screw trajectories and influenced the stability and mechanical responses. However, different implantation angles did not. Considering the difficulty of implantation, the ipsilateral‐contralateral trajectory in the lateral middle region of the VB can be optimal for in situ screw implantation in LLIF surgery.

Finite element modeling of temporal bone graft changes in XLIF: Quantifying biomechanical effects at adjacent levels

Extreme lateral interbody fusion allows for the insertion of a large-footprint interbody cage while maintaining the presence of natural stabilizing ligaments and the facets. It is unclear how the load-distribution mechanisms through these structures alter with temporal changes in the bone graft. The aim of this study was to examine the effects of temporal bone graft changes on load distribution among the cage, graft, and surrounding spinal structures using finite element analysis. Thoracolumbosacral spine computed tomography data from an asymptomatic male subject were segmented into anatomical regions of interest and digitally stitched to generate a surface mesh of the lumbar spine (L1-S1). The interbody cage was inserted into the L4-L5 region during surface meshing. A volumetric mesh was generated and imported into finite element software for pre-processing, running nonlinear static solves, and post-processing. Temporal stiffening was simulated in the graft region with unbonded (Soft Callus, Temporal Stages 1-3, Solid Graft) and bonded (Partial Fusion, Full Fusion) contact. In flexion and extension, cage stress reduced by 20% from the soft callus to solid graft state. Force on the graft was directly related to its stiffness, and load-share between the cage and graft improved with increasing graft stiffness, regardless of whether contact was fused with the endplates. Fused contact between the cage-graft complex and the adjacent endplates shifted load-distribution pathways from the ligaments and facets to the implant, however, these changes did not extend to adjacent levels. These results suggest that once complete fusion is achieved, the existing load paths are seemingly diminished.

Single-level cervical disc arthroplasty in the spine with reversible kyphosis: A finite element study

Background

Our previous studies found the single-level cervical disc arthroplasty (CDA) might be a feasible treatment for the patients with reversible kyphosis (RK). Theoretically, the change of cervical alignment from lordosis to RK comes with the biomechanical alteration of prostheses and cervical spine. However, the biomechanical data of CDA in the spine with RK have not been reported. This study aimed at establishing finite element (FE) models to (1) explore the effects of RK on the biomechanics of artificial cervical disc; (2) investigate the biomechanical differences of single-level anterior cervical discectomy and fusion (ACDF) and CDA in the cervical spine with RK.

Methods

The FE models of the cervical spine with lordosis and RK were constructed, then three single-level surgical models were developed: (1) RK + ACDF; (2) RK + CDA; (3) lordosis + CDA. A 73.6-N follower load combined with 1 N·m moment was applied at the C2 vertebra to produce cervical motion.

Results

At the surgical level, "lordosis + CDA" had the greatest ROM (except for flexion) while "RK + ACDF" had the minimum ROM. However, at adjacent levels, the ROM of "RK + ACDF" increased by 4.05% to 38.04% in comparison to "RK + CDA." "RK + ACDF" had the greatest prosthesis interface stress, while the maximum prosthesis interface stress of "RK + CDA" was at least 2.15 times higher than "lordosis + CDA." Similarly, "RK + ACDF" had the greatest intradiscal pressure (IDP) at adjacent levels, while the IDP of "RK + CDA" was 1.6 to 6.7 times higher than "lordosis + CDA." At the surgical level, "RK + CDA" had the greatest facet joint stress (except for extension), which was 1.9 to 11.2 times higher than "lordosis + CDA." At the adjacent levels, "RK + CDA" had the greatest facet joint stress (except for extension), followed by "RK + ACDF" and "lordosis + CDA" in descending order.

Conclusions

RK significantly changed the biomechanics of CDA, which is demonstrated by the decreased ROM and the significantly increased prosthesis interface stress, IDP, and facet joint stress in the "RK + CDA" model. Compared with ACDF, CDA overall exhibited a better biomechanical performance in the cervical spine with RK, with the increased ROM of surgical level and facet joint stress and the decreased ROM of adjacent levels, prosthesis interface stress, and IDP.

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.

Computational model predicts risk of spinal screw loosening in patients

PURPOSE

Pedicle screw loosening is a frequent complication in lumbar spine fixation, most commonly among patients with poor bone quality. Determining patients at high risk for insufficient implant stability would allow clinicians to adapt the treatment accordingly. The aim of this study was to develop a computational model for quantitative and reliable assessment of the risk of screw loosening.

METHODS

A cohort of patient vertebrae with diagnosed screw loosening was juxtaposed to a control group with stable fusion. Imaging data from the two cohorts were used to generate patient-specific biomechanical models of lumbar instrumented vertebral bodies. Single-level finite element models loading the screw in axial or caudo-cranial direction were generated. Further, multi-level models incorporating individualized joint loading were created.

RESULTS

The simulation results indicate that there is no association between screw pull-out strength and the manifestation of implant loosening (p = 0.8). For patient models incorporating multiple instrumented vertebrae, CT-values and stress in the bone were significantly different between loose screws and non-loose screws (p = 0.017 and p = 0.029, for CT-values and stress, respectively). However, very high distinction (p = 0.001) and predictability (R²_Pseudo = 0.358, AUC = 0.85) were achieved when considering the relationship between local bone strength and the predicted stress (loading factor). Screws surrounded by bone with a loading factor higher than 25% were likely to be loose, while the chances of screw loosening were close to 0 with a loading factor below 15%.

CONCLUSION

The use of a biomechanics-based score for risk assessment of implant fixation failure might represent a paradigm shift in addressing screw loosening after spondylodesis surgery.

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).