Role of posterior elements in the disc bulging of a degenerated cervical spine

Finite element analysis of two-level discontinuous cervical hybrid revision surgery strategy to reduce biomechanical responses of adjacent segments

Background

Hybrid surgery (HS) combined cervical disc arthroplasty (CDA) with anterior cervical discectomy and fusion (ACDF) is emerging, but its biomechanical effects as a revision surgery (RS) on adjacent segments were unclear.

Objectives

This finite element (FE) study aimed to investigate the biomechanical characteristics of HS to treat two-level discontinuous ASD in ACDF RS.

Methods

A C2-T1 intact FE model was established and modified to a primary C5/6 ACDF model and five RS models. These RS models' segments C4/5 and C6/7 were revised using cage plus plate (C), zero-profile devices (P), and Bryan disc (D), respectively, generating C-C-C, P-C-P, D-C-P, P-C-D, and D-C-D models. In the intact and C5/6 ACDF models, a 1.0 Nm moment was used to produce the range of motion (ROM). A displacement load was applied to all RS models, to achieve a total ROM match that of the primary C5/6 ACDF model.

Results

In the P-C-P model, biomechanical responses including ROM, Intradiscal pressure (IDP), Facet joint force (FJF), and Maximum von Mises stresses of discs at segments C3/4 and C7/T1 were slightly lower than the C-C-C model. The biomechanical response parameters at segments C3/4 and C7/T1 of P-C-D, D-C-P, and D-C-D were smaller than those in C-C-C and P-C-P models. D-C-D had the most significant effect on reducing all biomechanical responses among all RS models in segments C3/4 and C7/T1. Moreover, the disc stress cloud maps showed that the maximum von Mises stress of the C3/4 disc was higher than that of C7/T1.

Conclusions

D-C-D, P-C-D, and D-C-P are good RS choices for reducing the biomechanical responses, and D-C-D was the best choice. P-C-P can be the best recommendation when it does not meet the CDA indications. This study provided a biomechanical reference for hybrid surgical decision-making in the ACDF RS for preventing ASD recurrence.

Development of a spinopelvic complex finite element model for quantitative analysis of the biomechanical response of patients with degenerative spondylolisthesis

Research on degenerative spondylolisthesis (DS) has focused primarily on the biomechanical responses of pathological segments, with few studies involving muscle modelling in simulated analysis, leading to an emphasis on the back muscles in physical therapy, neglecting the ventral muscles. The purpose of this study was to quantitatively analyse the biomechanical response of the spinopelvic complex and surrounding muscle groups in DS patients using integrative modelling. The findings may aid in the development of more comprehensive rehabilitation strategies for DS patients. Two new finite element spinopelvic complex models with detailed muscles for normal spine and DS spine (L4 forwards slippage) modelling were established and validated at multiple levels. Then, the spinopelvic position parameters including peak stress of the lumbar isthmic-cortical bone, intervertebral discs, and facet joints; peak strain of the ligaments; peak force of the muscles; and percentage difference in the range of motion were analysed and compared under flexion-extension (F-E), lateral bending (LB), and axial rotation (AR) loading conditions between the two models. Compared with the normal spine model, the DS spine model exhibited greater stress and strain in adjacent biological tissues. Stress at the L4/5 disc and facet joints under AR and LB conditions was approximately 6.6 times greater in the DS spine model than in the normal model, the posterior longitudinal ligament peak strain in the normal model was 1/10 of that in the DS model, and more high-stress areas were found in the DS model, with stress notably transferring forwards. Additionally, compared with the normal spine model, the DS model exhibited greater muscle tensile forces in the lumbosacral muscle groups during F-E and LB motions. The psoas muscle in the DS model was subjected to 23.2% greater tensile force than that in the normal model. These findings indicated that L4 anterior slippage and changes in lumbosacral-pelvic alignment affect the biomechanical response of muscles. In summary, the present work demonstrated a certain level of accuracy and validity of our models as well as the differences between the models. Alterations in spondylolisthesis and the accompanying overall imbalance in the spinopelvic complex result in increased loading response levels of the functional spinal units in DS patients, creating a vicious cycle that exacerbates the imbalance in the lumbosacral region. Therefore, clinicians are encouraged to propose specific exercises for the ventral muscles, such as the psoas group, to address spinopelvic imbalance and halt the progression of DS.

Biomechanical study between percutaneous vertebroplasty combined with cement pedicle plasty improves vertebral biomechanical stability: A finite element analysis

OBJECTIVE

To investigate the biomechanical effects of percutaneous vertebroplasty combined with cement pedicle plasty (PVCPP) on the unstable osteoporotic vertebral fractures (OVFs) through finite element (FE) analysis. The study compares the biomechanical stability of finite element models between percutaneous vertebroplasty (PVP) and percutaneous vertebroplasty combined with cement pedicle plasty.

METHODS

Two patients with unstable OVFs underwent computed tomography (CT) examination at the thoracolumbar vertebral body levels, respectively. The CT images were reconstructed into three-dimensional finite element models to simulate stress conditions across six dimensions and to evaluate the vertebral von Mises stress before and after bone cement reinforcement.

RESULTS

The study found that stress distribution differed between groups mainly at the pedicle base. In the surgical vertebral bodies, the maximum stress in the PVP group decreased during flexion and left bending, while it increased in other states. In the PVCPP group, all maximum stresses decreased. In the inferior vertebral bodies, the maximum stress in the PVP group generally increased, while it decreased in the PVCPP group. In the superior vertebral bodies, postoperatively, the maximum stress in the PVP group generally increased, while it almost remained unchanged in the PVCPP group. PVP group had higher cement stress and displacement.

CONCLUSION

PVCPP is an effective treatment method for patients with unstable OVFs. It can quickly relieve pain and enhance the stability of the three columns, thereby reducing the risk of some complications.

The effect of large channel-based foraminoplasty on lumbar biomechanics in percutaneous endoscopic discectomy: a finite element analysis

BACKGROUND

This study aimed to evaluate the effect of foraminoplasty using large-channel endoscopy during TESSYS on the biomechanics of the lumbar spine.

METHODS

A complete lumbar spine model, M1, was built using 3D finite elements, and models M2 and M3 were constructed to simulate the intraoperative removal of the superior articular process of L5 using a trephine saw with diameters of 5 mm and 8.5 mm, respectively, and applying normal physiological loads on the different models to simulate six working conditions-anterior flexion, posterior extension, left-right lateral bending, and left-right rotation-to investigate the displacement and facet joint stress change of the surgical segment, and the disc stress change of the surgical and adjacent segments.

RESULTS

Compared with the M1 model, the M2 and M3 models showed decreased stress at the L4-5 left FJ and a significant increase in stress at the right FJ in forward flexion. In the M2 and M3 models, the L4-5 FJ stresses were significantly greater in left lateral bending or left rotation than in right lateral bending or right rotation. The right FJ stress in M3 was greater during left rotation than that in M2, and that in M2 was greater than that in M1. The L4-5disc stress in the M3 model was greater during posterior extension than that in the M1 and M2 models. The L4-5disc stress in the M3 model was greater in the right rotation than in the M2 model, and that in the M2 model was greater than that in the M1 model.

CONCLUSION

Foraminoplasty using large-channel endoscopy could increase the stress on the FJ and disc of the surgical segment, which suggested unnecessary and excessive resection should be avoided in PTED to minimize biomechanical disruption.

The impact of disc degeneration on the dynamic characteristics of the lumbar spine: a finite element analysis

The dynamics of disc degeneration was analyzed to determine the effect of disc degeneration at the L4-L5 segment on the dynamic characteristics of the total lumbar spine. A three-dimensional nonlinear finite element model of the L1-S1 normal lumbar spine was constructed and validated. This normal model was then modified to construct two degeneration models with different degrees of degeneration (mild, moderate) at the L4-L5 level. Modal analysis, harmonic response analysis, and transient dynamics analysis were performed on the total lumbar spine when experiencing following compressive loading (500 N). As the degree of disc degeneration increased, the vibration patterns corresponding to the first three orders of the model's intrinsic frequency were basically unchanged, with the first order being in the left-right lateral bending direction, the second order being in the forward-flexion and backward-extension direction, and the third order being in the axial stretching direction. The nucleus pulposus pressure peaks corresponding to the first order intrinsic frequency for the harmonic response analysis are all on the right side of the model, with sizes of 0.053 MPa, 0.061 MPa, and 0.036 MPa, respectively; the nucleus pulposus pressure peaks corresponding to the second order intrinsic frequency are all at the rear of the model, with sizes of 0.13 MPa, 0.087 MPa, and 0.11 MPa, respectively; and the nucleus pulposus pressure peaks corresponding to the third order intrinsic frequency are all at the front of the model, with sizes of 0.19 MPa, 0.22 MPa, and 0.22 MPa, respectively. The results of the transient analysis indicated that over time, the response curves of the healthy model, the mild model, and the moderate model all exhibited cyclic response characteristics. Intervertebral disc degeneration did not adversely affect the vibration characteristics of the entire lumbar spine system. Intervertebral disc degeneration significantly altered the dynamics of the degenerative segments and their neighboring normal segments. The process of disc degeneration gradually shifted the load from the nucleus pulposus to the annulus fibrosus when the entire lumbar spine was subjected to the same vibratory environment.

Finite element modeling and analysis of effect of preexisting cervical degenerative disease on the spinal cord during flexion and extension

Recent studies have emphasized the importance of dynamic activity in the development of myelopathy. However, current knowledge of how degenerative factors affect the spinal cord during motion is still limited. This study aimed to investigate the effect of various types of preexisting herniated cervical disc and the ligamentum flavum ossification on the spinal cord during cervical flexion and extension. A detailed dynamic fluid-structure interaction finite element model of the cervical spine with the spinal cord was developed and validated. The changes of von Mises stress and maximum principal strain within the spinal cord in the period of normal, hyperflexion, and hyperextension were investigated, considering various types and grades of disc herniation and ossification of the ligamentum flavum. The flexion and extension of the cervical spine with spinal canal encroachment induced high stress and strain inside the spinal cord, and this effect was also amplified by increased canal encroachments and cervical hypermobility. The spinal cord might evade lateral encroachment, leading to a reduction in the maximum stress and principal strain within the spinal cord in local-type herniation. Although the impact was limited in the case of diffuse type, the maximum stress tended to appear in the white matter near the encroachment site while compression from both ventral and dorsal was essential to make maximum stress appear in the grey matter. The existence of canal encroachment can reduce the safe range for spinal cord activities, and hypermobility activities may induce spinal cord injury. Besides, the ligamentum flavum plays an important role in the development of central canal syndrome.Significance. This model will enable researchers to have a better understanding of the influence of cervical degenerative diseases on the spinal cord during extension and flexion.

Analyzing isolated degeneration of lumbar facet joints: implications for degenerative instability and lumbar biomechanics using finite element analysis

The facet joint contributes to lumbar spine stability as it supports the weight of body along with the intervertebral discs. However, most studies on the causes of degenerative lumbar diseases focus on the intervertebral discs and often overlook the facet joints. This study aimed to investigate the impact of facet joint degeneration on the degenerative changes and diseases of the lumbar spine. A finite element model of the lumbar spine (L1-S1) was fabricated and validated to study the biomechanical characteristics of the facet joints. To simulate degeneration of the facet joint, the model was divided into four grades based on the number of degenerative segments (L4-L5 or L4-S1) and the contact condition between the facet joint surfaces. Finite element analysis was performed on four spine motions: flexion, extension, lateral bending, and axial torsion, by applying a pure moment to the upper surface of L1. Important parameters that could be used to confirm the effect of facet joint degeneration on the lumbar spine were calculated, including the range of motion (ROM) of the lumbar segments, maximum von Mises stress on the intervertebral discs, and reaction force at the facet joint. Facet joint degeneration affected the biomechanical characteristics of the lumbar spine depending on the movements of the spine. When analyzed by dividing it into degenerative onset and onset-adjacent segments, lumbar ROM and the maximum von Mises stress of the intervertebral discs decreased as the degree of degeneration increased in the degenerative onset segments. The reaction force at the facet joint decreased with flexion and increased with lateral bending and axial torsion. In contrast, lumbar ROM of the onset-adjacent segments remained almost unchanged despite severe degeneration of the facet joint, and the maximum von Mises stress of the intervertebral discs increased with flexion and extension but decreased with lateral bending and axial torsion. Additionally, the facet joint reaction force increased with extension, lateral bending, and axial rotation. This analysis, which combined the ROM of the lumbar segment, maximum von Mises stress on the intervertebral disc, and facet joint reaction force, confirmed the biomechanical changes in the lumbar spine due to the degeneration of isolated facet joints under the load of spinal motion. In the degenerative onset segment, spinal instability decreased, whereas in the onset-adjacent segment, a greater load was applied than in the intact state. When conducting biomechanical studies on the lumbar spine, considering facet joint degeneration is important since it can lead to degenerative spinal diseases, including adjacent segment diseases.

Stability simulation analysis of targeted puncture in L4/5 intervertebral space for PELD surgery

Introduction: The application prospects of percutaneous endoscopic lumbar discectomy (PELD) as a minimally invasive spinal surgery method in the treatment of lumbar disc herniation are extensive. This study aims to find the optimal entry angle for the trephine at the L4/5 intervertebral space, which causes less lumbar damage and has greater postoperative stability. To achieve this, we conduct a three-dimensional simulated analysis of the degree of damage caused by targeted puncture-based trephine osteotomy on the lumbar spine. Methods: We gathered clinical CT data from patients to construct a lumbar model. This model was used to simulate and analyze the variations in trephine osteotomy volume resulting from targeted punctures at the L4/5 interspace. Furthermore, according to these variations in osteotomy volume, we created Finite Element Analysis (FEA) models specifically for the trephine osteotomy procedure. We then applied mechanical loads to conduct range of motion and von Mises stress analyses on the lumbar motion unit. Results: In percutaneous endoscopic interlaminar discectomy, the smallest osteotomy volume occurred with a 20° entry angle, close to the base of the spinous process. The volume increased at 30° and reached its largest at 40°. In percutaneous transforaminal endoscopic discectomy, the largest osteotomy volume was observed with a 50° entry angle, passing through the facet joints, with smaller volumes at 60° and the smallest at 70°. In FEA, M6 exhibited the most notable biomechanical decline, particularly during posterior extension and right rotation. M2 and M3 showed significant differences primarily in rotation, whereas the differences between M3 and M4 were most evident in posterior extension and right rotation. M5 displayed their highest stress levels primarily in posterior extension, with significant variations observed in right rotation alongside M4. Conclusion: The appropriate selection of entry sites can reduce lumbar damage and increase stability. We suggest employing targeted punctures at a 30° angle for PEID and at a 60° angle for PTED at the L4/5 intervertebral space. Additionally, reducing the degree of facet joint damage is crucial to enhance postoperative stability in lumbar vertebral motion units.

Biomechanical analysis of sandwich vertebrae in osteoporotic patients: finite element analysis

Objective

The aim of this study was to investigate the biomechanical stress of sandwich vertebrae (SVs) and common adjacent vertebrae in different degrees of spinal mobility in daily life.

Materials and methods

A finite element model of the spinal segment of T10-L2 was developed and validated. Simultaneously, T11 and L1 fractures were simulated, and a 6-ml bone cement was constructed in their center. Under the condition of applying a 500-N axial load to the upper surface of T10 and immobilizing the lower surface of L2, moments were applied to the upper surface of T10, T11, T12, L1, and L2 and divided into five groups: M-T10, M-T11, M-T12, M-L1, and M-L2. The maximum von Mises stress of T10, T12, and L2 in different groups was calculated and analyzed.

Results

The maximum von Mises stress of T10 in the M-T10 group was 30.68 MPa, 36.13 MPa, 34.27 MPa, 33.43 MPa, 26.86 MPa, and 27.70 MPa greater than the maximum stress value of T10 in the other groups in six directions of load flexion, extension, left and right lateral bending, and left and right rotation, respectively. The T12 stress value in the M-T12 group was 29.62 MPa, 32.63 MPa, 30.03 MPa, 31.25 MPa, 26.38 MPa, and 26.25 MPa greater than the T12 stress value in the other groups in six directions. The maximum stress of L2 in M-T12 in the M-L2 group was 25.48 MPa, 36.38 MPa, 31.99 MPa, 31.07 MPa, 30.36 MPa, and 32.07 MPa, which was greater than the stress value of L2 in the other groups. When the load is on which vertebral body, it is subjected to the greatest stress.

Conclusion

We found that SVs did not always experience the highest stress. The most stressed vertebrae vary with the degree of curvature of the spine. Patients should be encouraged to avoid the same spinal curvature posture for a long time in life and work or to wear a spinal brace for protection after surgery, which can avoid long-term overload on a specific spine and disrupt its blood supply, resulting in more severe loss of spinal quality and increasing the possibility of fractures.

Biomechanical analysis of the cervical spine segment as a method for studying the functional and dynamic anatomy of the human neck

BACKGROUND

Traditionally, dynamic and functional anatomy, in particular the dynamic anatomy of the neck, is studied on cadaveric material. However, the development of in vivo visualization technologies and in silico modeling has made it possible to expand these possibilities. Despite significant progress in the study of dynamic and functional anatomy of the neck by means of in silico methods, the issues of validating the developed models and taking into account the pronounced nonlinearity of soft tissues as well as local anisotropy remain open. The aim of this study was to develop a virtual dynamic anatomical model of the human neck and reproduce the dynamic processes in the cervical spine from this model using the finite element method.

MATERIALS AND METHODS

Reverse engineering was used to generate a dynamic anatomical model of the neck from CT data (both male, 24 and 22 years old). Two segments of the cervical spine (C3-C5, C2-T1) were isolated from the resulting model for finite element analysis. Finite element mesh generation and contact interactions were performed using the HyperMesh software (Altair Engineering Inc, Troy, Michigan, USA). The anisotropic hyperelastic Holzapfel-Gasser-Ogden model was used to describe the material behavior of the fibrous rings of the disc. Material modeling and finite element analysis were performed using Abaqus CAE 6.14 software (Simulia, Johnston, Rhode Island, USA).

RESULTS

A technique for creating a virtual dynamic anatomical model of the neck was elaborated and implemented. The model includes 79 major anatomical structures of the neck segmented from radiological data. A finite element analysis of the cervical spine was performed. The results of finite element analysis of the C3-C5 segment under axial load were compared with in vitro data. The proposed model shows nonlinear deformation of the disc under static loading; the model predicted displacement values agree well with the experimental ones. The displacement of the С3-С5 central vertebra with an axial load of 800 N reaches a value of 0.65 mm. For the segment C2-T1, data on intradiscal pressure, stress plots and displacements during flexion were obtained. The maximum stress value of 10.036 MPa is observed in the C3-C4 disc.

CONCLUSION

Simulation results using the proposed methodology are in good agreement with experimental data. The generated biomechanical models allow describing dynamic phenomena in the cervical spine and obtaining a wide range of quantitative properties of anatomical objects, which are otherwise inaccessible to classical methods for studying dynamic and functional anatomy.

Rapid increase in intradiscal pressure in porcine cervical spine units negatively impacts annulus fibrosus strength

Vertebral endplate fracture is generally accepted to occur as a result of excessive, typically rapid, pressure developed within the intervertebral disc. Clinical evidence of disc disorders later in life following endplate fracture suggests that the disc is also impacted by the fracture event. There is an abundance of evidence to indicate that loss of disc pressure at the time of fracture alters annulus fibrosus loading which could impact disc health long-term. However, little is known regarding the isolated impact of the rapid disc pressurization during the fracture event. This study sought to determine if rapid disc pressurization impacts mechanical properties of the annulus. Intervertebral discs of porcine cervical functional spine units (FSU) were pressurized via rapid injection of hydraulic fluid (standard motor oil). FSUs were subsequently grouped according to whether the pressurization resulted in endplate fracture ('fracture'; n = 9) or not ('non-fracture'; n = 13). Annular tensile properties and lamellar adhesion strength were quantified for both fractured and non-fractured samples. A third group was also examined in which FSUs did not undergo the pressurization protocol ('control'; n = 10). Annular tensile properties were not found to be statistically different between the three conditions (p > 0.05). Lamellar adhesion strength, in contrast, was found to be, on average, 53% lower in specimens which were pressurized compared to control, non-pressurized samples (p = 0.03). Further, no statistical differences in adhesion strength was found between fractured and non-fractured pressurized samples (p = 0.990) indicating that rapid pressurization , and not the presence of an endplate fracture, impacted the mechanical properties of the annulus. Decreased lamellar adhesion strength increases the potential of delamination, which can lead to disc disorders.