Topping-off surgery vs posterior lumbar interbody fusion for degenerative lumbar disease: a finite element analysis

Abstract

Background Adjacent segment disease (ASD) is a common complication after posterior lumbar interbody fusion (PLIF). Recently, a topping-off surgery (non-fusion with Coflex) has been developed to reduce the risk of ASD, yet whether and how the topping-off surgery can relieve ASD remains unclear. The purpose of this study was to explore the biomechanical effect of PLIF and Coflex on the adjacent segments via finite element (FE) analysis and discuss the efficacy of Coflex in preventing ASD.

Methods A FE model of L3-5 segments was generated from the CT of a healthy volunteer via three commercially available software. Coflex and PLIF devices were modeled and implanted with the segment model in the FE software. In the FE model, a pre-compressive load of 500 N, which was equal to two thirds of the human body mass, was applied to the top surface of the L3. In addition, four types of moments (anteflexion, rear protraction, bending and axial rotation) set as 10 Nm were successively applied to the FE model combined with this pre-compressive load. Then the range of motion (ROM), the torsional rigidity and the maximum von Mises equivalent stress on L3-L4 disc and the implant were analyzed.

Results Both Coflex and PLIF reduced ROM. However, no significant difference was found in the maximum von Mises equivalent stress of adjacent segment disc between the two devices. Interestingly enough, both systems increased the torsional rigidity at the adjacent lumbar segment, and PLIF had a more significant increase. The Coflex implant had a larger maximum von Mises equivalent stress.

Conclusions Both Coflex and PLIF reduced ROM at L3-L4, and thus improved the lumbar stability. Under the same load, both devices had almost the same maximum von Mises equivalent stress as the normal model on the adjacent intervertebral disc. But it is worthy to notice the torsional rigidity of PLIF was higher than that of Coflex. It indicated that the lumbar treated with PLIF would undertake a larger load to reach ROM of Coflex. Therefore, we presumed that ADS was related to higher torsional rigidity.

Background

Lumbar spinal stenosis (LSS), a common senile disease, has an incidence of 9.3% in the elderly [1]. LSS is a major cause of surgery in people over 65[2]. Clinical observations show the onset age of LSS is becoming younger on account of the sedentary lifestyle. LSS is clinically manifested by numbness and radiative pain in the buttocks and lower limbs [3]. Some patients may have lower back pain. These symptoms exert negative impacts on the patient’s physical functions. Some bedridden LSS patients even develop pneumonia and deep vein thrombosis (DVT).

Posterior lumbar interbody fusion (PLIF), as the “gold standard” in surgical practice, can significantly relieve the symptoms of nerve root compression. Although the primary segment disc compression had a low recurrence rate after PLIF, during the long-term follow-up, many patients are found to develop LSS at an adjacent segment and recurrent relative nerve root compression. Some authors have shown PLIF surgery may increase the incidence of adjacent segment disease (ASD) [3][4]. Currently, a new “topping-off” technique using Coflex, a U-shaped elastic non-fusion interspinous instrument placed between two adjacent spinous processes, is often applied in clinical practice. Some literature has demonstrated Coflex implantation can reduce the risk of ASD [3,4].

Various studies have proved that Colfex is safe and effective for LSS. However, the effect of the Coflex on the adjacent segment is still not clear. Lee N [4] believes ASD after PLIF is caused by the compensatory increase in the mobility of adjacent lumbar segment and the added pressure on the intervertebral disc and intervertebral joints. We presumed that torsional rigidity might be related to ASD but few studies were concerned with the torsional rigidity. To figure out the specific mechanism of how the Coflex reduced the rate of ASD and the relationship between torsional rigidity and ASD, we conducted this study to explore the biomechanical effect of PLIF and Coflex on the adjacent segments via finite element analysis and discussed the efficacy of Coflex in preventing ASD.

Methods

1.1 Establishment of basic model

The geometrical specifications of the spine were obtained from 64 spiral CT images of a 28-year-old male without spine injury or radiographic evidence of degeneration. He received a SOMATOM SENSATION 64 spiral CT (Siemens, Munich, Germany) scan for health examination at our hospital. The CT images were used with his consent. The CT scans of the L3-L5 lumbar spine with Coflex or PLIF were obtained at 1-mm intervals. The threshold segmentation was used in MIMICS19 based on the CT data to obtain the masks of intervertebral disc, nucleus pulposus, vertebral body and the implant. Using the masks, the 3-D solid model was constructed as shown in Figure 1.1 (1). The contour of the model was refined and the CT image distortion was decreased as much as possible. Then the STL file was output with the quality of the mesh maintained. In GEOMAGIC12, the STL vertebral model was de-noised and smoothed with all the details maintained. All vertebral models were offset by 1 mm to make a smaller vertebral cancellous bone model. The STP files of all models were output and imported into UG software, and then the final models were obtained by Boolean subtraction [5–9].

Threshold segmentation was used in MIMICS19 based on postoperative CT to obtain the Coflex implant masks. The preliminary solid model was established by 3-D model construction using the masks. The output STL file was imported into the UG software, and the Coflex model was designed based on the preliminary model using the UG software, as shown in Figure 1.1(2).

The nail-rod model and the cage model were designed and assembled using the UG software. The cage was individually designed to make the upper and lower surfaces of the cage completely match the adjacent vertebral body surfaces. The Arbeitsgemeinschaftfür Osteosynthesefragen (AO) spinal internal fixation standard was used as reference. The screws were implanted in the center of the pedicle and a 20 ° cohesion angle was used in the L4-L5 pedicle screws, as shown in Figure 1.1 (3) (4).

The normal lumbar spine, Coflex and K-ROD assembly models were imported into ANSYS 17. Then ANSYS analysis files were generated for each model as shown in Figure 1.1(5). Table 1.1(1) shows the material properties of the bones, intervertebral disc, implants, and various ligaments are shown in. In ANSYS, the simulation of ligaments was achieved by adding a spring unit (tension only). Table 1.1(2) shows the properties of the spring unit. [10]

1.2 Pre-compressive load

A 500N pre-compressive load equal to two thirds of the human body mass (50kg, 500N) [5] was applied on the upper surface of L3 lumbar vertebra. An additional moment of 10N*m (10000N*mm) was exerted in each direction as described in the literature [5–12] to simulate lumbar anteflexion, rear protraction and rotation (Figure 1.2).

1.3 Simplified restriction

A six-degree-of-freedom rigidly fixed constraint was applied to the lower face of L5 lumbar according to the previous literature [11,12] to verify that the lower face of L5 vertebral would not produce displacement and rotation with added moment.

1.4 Bind

All ligaments were simplified into spring units as shown in Table 1.2(2). The BOND connection was used between the vertebral body and the intervertebral disc. The facet joint surface was subjected to frictional simulation and the friction coefficient was set as 0.2 [4–9,11–13]. “No Separation” was used to connect the surfaces of Coflex and those of the spinous processes. BOND connection was applied in all the K-ROD components, cage and vertebral body connections.

1.5 Meshing

To ensure the comparability of the models and prompt calculation, and also to avoid the calculation error caused by meshing, a 10-node tetrahedral mesh was used in mesh 7. A more refined mesh of 1.5mm was used in L3/4 intervertebral disc and the inner implant. The octahedral mesh was used in the cancellous bone. The same meshing method was applied in the normal lumbar, Coflex and PLIF implant models. The mesh of 3mm was used in all the other parts. The meshing was performed using ANSYS software (ANSYS WORKBENCH 15.0). The units used are shown in Table 1.1(1). No warnings and failures were reported. Table 1.5 shows the number of nodes and elements in the three models.

2 The verification method

2.1 Verification of L4, L5 segment axial compression (displacement-load curve)

After establishing the L4/L5 lumbar model, a pre-compressive load ranging from 500 to 2000 N was applied on the upper surface of L4. The axial displacement of the lumbar vertebral model was calculated and compared with the references [14–18], as shown in Figure 2.1

2.2Measurement of lumbar range of motion and torsional rigidity

The origin of the local coordinate system was defined as the center of endplate in the distal vertebral body. Perpendicular to the endplate in the sagittal position was the X-axis. Parallel to the endplate was the Z-axis (right) and the Y-axis (left).

With the distal vertebral body fixed, two nodes on the upper surface of the adjacent lumbar were selected. The coordinates of the two nodes before pre-compressive load application were (X1, Y1, Z1), (X2, Y2, Z2). The coordinates of the two nodes after pre-compressive load application were (X3, Y3, Z3), (X4, Y4, Z4). The lines connecting the two points before and after pre-compressive load application formed the angle θ (i.e. the maximum angular displacement). [11,19–21] The formula was presented as follows:

[Due to technical limitations, the formula could not be displayed here. Please see the supplementary files section to see the formula.

Torsional rigidity was the moment required to produce a unit torsion of the lumbar and was measured via the following formula: Torsional rigidity(k) = M/θ, (M: the exerted moment; θ:the range of motion).

2.3 Verification of maximum von Mises equivalent stress on L3/L4 intervertebral disc

The uppermost L3 is not subject to any constraints. An evenly distributed 500 N pre-compressive load was exerted on the endplate of the L3 vertebral body; the motor moment was 10 Nm. [8,19–20]. The maximum von Mises equivalent stress on the L3/L4 intervertebral disc was calculated by FE analysis and compared with the references.

2.4 Evaluation of the sensitivity of the elements

The low-medium, medium, medium-high, optimized and high densities of five representative meshes were performed to determine the number of elements, as shown in Fig 2.4. An analysis was performed using the five densities, and the results were converged from the medium- to high mesh densities. Table 2.2.4 lists the number of elements, number of nodes, computational time, max model displacement and equivalent stress on the adjacent disc for each density.

Results

3 Verification of experimental results

3.1 Verification of L4, L5 segment axial compression (displacement-load curve)

Compared with the results from the literature [1418], ours was approximated to that of Virgin’s [18], and in-between of those from Virgin’s and Markolf’s. The axial compression stiffness of our model was of an intermediate level among the results in literature. The curve was obviously nonlinear, which was in accordance with the literature (Figure 3.1).

3.2 Measurement of lumbar ROM and torsional rigidity

The ROM and torsional rigidity of this study is shown in Table 3.2(1). The ROM and torsional rigidity of the models were in good accordance with those reported in literature, as shown in Table 3.2(2).

3.3 Verification of maximum von Mises equivalent stress on L3/L4 intervertebral disc

The maximum von Mises equivalent stress on the L3-L4 intervertebral disc in this model was in accordance with those from the other literature, as shown in Table 3.3.

3.4 Evaluation of the sensitivity of the elements

For high and medium-high mesh, the convergence and accuracy were increased; however, they were not practiced on account of the increased computational time. Therefore, the optimized mesh with a high degree of convergence was selected for the analysis as shown in Table 3.4 and Fig. 3.4.

4. L3-L4 lumbar ROM and torsional rigidity

In the normal model, the lumbar ROM of the L3-L4 was higher than the others. The Coflex and PLIF models exhibited significant differences in anteflexion, rear protraction and right lateral bending. The Coflex model had a significantly smaller L3-L4 axial rotation angle than the normal model and the rotation in upright position had the smallest ROM. The right axial rotation showed the largest ROM. However, the ROMs of upright, rear protraction and left lateral bending showed no significant change compared with the PLIF model. In the PLIF model, the up-right position had the smallest L3-L4 ROM and the biggest right axial rotation. The ROM of PLIF showed significant reductions, but the left lateral bending did not change significantly, compared with the normal and Coflex models, as shown in Table 4(1).

In terms of the L3-L4 lumbar spine torsional rigidity, the PLIF model was significantly higher than the normal and Coflex models, and the Coflex model was higher than the normal. In the normal model, the rear protraction had the smallest stiffness and the left lateral bending had the largest stiffness. In the Coflex model, the stiffness of L3-L4 right axial rotation was the smallest and the right lateral bending rotation was the largest; compared with the normal, the stiffness of anteflexion, rear protraction, left lateral bending and left axial rotation was significantly increased. The average L3-L4 lumbar torsional rigidity reached 3.708858 N*m/°, which was significantly larger than the normal model. In the PLIF model, compared with Coflex, L3-L4 lumbar torsional rigidity was significantly increased in all the positions except in rear protraction and the right lateral bending. All the data are presented in Table 4(2).

5Maximum von Mises equivalent stress on L3/L4 intervertebral disc

The maximum von Mises equivalent stress on L3-L4 intervertebral discs had no significant difference among the three models. In the normal model, the intervertebral disc had the largest maximum von Mises equivalent stress when in the right lateral bending position. And the smallest stress appeared when in the upright position. The maximum von Mises equivalent stresses in all positions were smaller than the other two models (Figure 5(1)). In the Coflex group, the maximum von Mises equivalent stress of the intervertebral disc was the largest in right lateral bending and the smallest in the upright position (Figure 5(2)). The maximum von Mises equivalent stress of L3-L4 intervertebral disc in the PLIF model did not differ significantly from that of the Coflex model (Figure 5(3)). All the data are presented in Table 5.

6 Maximum von Mises equivalent stress of implant

In this study, we performed the finite element analysis of Coflex, cage, and nail rods in each model to measure the maximum von Mises equivalent stress of the implant. It showed that Coflex had the smallest maximum von Mises equivalent stress when in the upright position, and the smallest maximum von Mises equivalent stress appeared in the right axial rotation. The maximum von Mises equivalent stress of the cage was the smallest compared with the other two parts, and the smallest stress appeared in the upright position and the largest in right axial rotation. The smallest stress on the nail was in anteflexion, unlike the other two implants. The largest stress was at the right lateral bending, as shown Figure 6(1)-(3). All data are presented in Table 6.

Discussion

PLIF surgery, the “the golden standard” in treating severe LSS, can significantly relieve LSS symptoms (e.g. radical pain and palsy), but it can also accelerate the degeneration of the adjacent segment. Paul Park [21] reported a high incidence of ASD after PLIF surgery. According to two other studies, the incidence even reached 100% [22]. Shinya Okuda [22] reported a case of repeated ASD after PLIF in which the patient had four times of operation at different adjacent lumbar segments. Although the patient’s symptoms resolved after the operation, the time between surgeries was increasingly shorter and the sacral slope decreased.

The “topping-off” surgical method—using either a hybrid stabilization device (HSD) or an inter-spinous process device (IPD)—is being widely used in treating ASD. Khoueir[23] et al. have classified posterior dynamic stabilization devices into three categories: (1) HSD with pedicle screw/rod instrument; (2) IPD such as Wallis and Coflex; (3) total facet replacement system. The device placed at the default adjacent segment after PLIF can effectively reduce the ROM and loading force on the adjacent segment. The systematic review by Po-Hsin Chou [24] showed that the fusion-alone group had a higher incidence of radiographic ASD (52.6%) and symptomatic ASD (11.6%) than revision surgery group (8.1%). Besides, the HSD and fusion groups had a higher incidence of radiographic ASD at the supra-adjacent segment (10.5% and 24.7%, respectively) than the IPD (1%). In a review of 91 cases, Lu K [25] evaluated and compared the incidence of ASD in the PLIF and the “topping-off” (an IPD was put at the adjacent segment proximal to the PLIF instrument) groups, showing that radiographic ASD occurred in 20 cases (48%) in the PLIF group and 3 (6%) in the “topping-off” group. The PLIF group had 9 symptomatic cases while only 3 were found in the “topping-off” group. He concluded that the “topping-off” device reduced the risk of ASD. Both PLIF and Wallis relieved the LSS symptoms [25]. In addition, Wallis significantly reduced the adjacent segment degeneration.

Coflex is an inter-spinous fixation device for the “topping-off” technology. It is a elastic U-shaped structure implanted after possessing the inter-spinous and supra-spinous ligaments. Qu SD [26]believed that Coflex relieved the nerve root compression by propping up the spinous process, maintaining the lordosis of the implanted segment, and reducing the ligamentum flavum invasion of the spinal canal. The distraction force also enlarged the intervertebral foramen and reduced the load on the intervertebral disc and the facet joint in flexion. The meta-analysis by Ai-min Li [27] showed that Coflex was more effective than PLIF in terms of decompression, the Oswestry Disability Index (ODI), length of hospital stay (LOS) and blood loss. Wei yuan [28] held that Coflex surgery had significantly less blood loss, shorter hospital stay and operative time than PLIF (p < 0.001), and it also had a lower reoperation rate for ASD than PLIF, yet with no statistical difference (11.1% vs. 4.8%, p = 0.277). He also thought Coflex was not suitable for the patients with lumbar disc herniation.

To investigate the mechanism of Coflex in reducing ASD incidence, mechanical analysis of the supra-lumbar segment is needed. Finite element analysis can simulate real physical systems (geometry and load cases) and perform the measurement using a mathematical approximation program. With simple and mutual elements, a finite unknown quantity can approximate a real system of infinitely unknown quantity.

In this study, no significant difference of intervertebral disc load force on the adjacent segment was found between Coflex and PLIF under a 500N pre-compressive load. The increased adjacent disc pressure of Coflex ranged from –5%–10%, while the pressure of PLIF ranged from –10% to 2%, indicating the higher load force on the adjacent disc was not a primary risk factor of ASD.

Coflex significantly reduced the anteflexion and rear protraction ROM of the adjacent segment, and the lateral bending was also decreased. The angle of axial rotation had almost no change. By comparison, PLIF reduced the ROM in all motions. The decreased ROM could promote lumbar stability in both Coflex and PLIF models. Two factors might contribute to the difference in Coflex activity between our result and that from the literature: first, ours was an idealized model and the micro-motion of the Coflex facet was smaller than the actual one; second, the intervertebral facet joint in our model was intact, unlike those in the literature, which were partially resected.

The increased torsional rigidity can impose a heavier burden on the adjacent segment. PLIF increases the torsional rigidity, which means the patient may have a higher load burden on the adjacent disc and facet joint in rotational activity. The higher load burden would then lead to degeneration of the adjacent intervertebral segments and the nucleus pulposus would be prone to protrusion and nerve root compression. As a result, the patients who has received PLIF would have a higher incidence of ADS and suffer from nerve root compression. Hence, higher torsional rigidity could be a primary risk factor of ASD.

Rigorous conclusions cannot be made based on the results of this study because of the small sample size and the number of working conditions for each finite element model. A larger sample size and more working conditions should be included in the further study. The maximum von Mises equivalent stress on the intervertebral disc needs to be statistically analyzed and compared. A comprehensive comparison of the effects of Coflex and PLIF on adjacent segments is also necessary in the future study.

Conclusion

Both Coflex and PLIF reduced lumbar ROM and therefore provided stability at the surgical and adjacent segments. Under the same pressure, both devices had almost the same maximum von Mises stress as the normal model on the adjacent intervertebral disc. But it is worthy to notice that the torsional rigidity of both PLIF and Coflex models were higher than that of the normal model, and PLIF produced an even higher burden on the adjacent segment than Coflex. It indicates that the lumbar of patient treated with PLIF would undertake a larger load to reach ROM of Coflex. Therefore, we presumed that ADS was related to higher torsional rigidity.

Declaration

Availability of data and material:

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Funding:

This research was supported by Medical Health Science and Technology Project of Hangzhou City (Grant No.:  20170533B38).
Ethics approval and consent to participate:

This article contains a study with human participants with his consent and no any animals performed by any of the authors.

Competing Interests

The authors declare no competing interests.

Contributions

YPF conceived and designed the study and wrote this manuscript. YPF searched database and reviewed studies. YPF performed data analysis and prepare the initial draft of the manuscript. LLZ and YPF participated in interpretation of data, helped in drafting the manuscript and critically reviewed the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank Nanjing medical University library system for all of the help and resources and our deepest gratitude goes to the anonymous reviewers and editors for their careful work and thoughtful suggestions that have helped improve this paper substantially.

Authors' information

Not applicable

Abbreviations

Lumbar spinal stenosis (LSS), deep vein thrombosis (DVT), Posterior lumbar interbody fusion (PLIF), adjacent segment disease (ASD), finite element (FE), Arbeitsgemeinschaftfür Osteosynthesefragen (AO), hybrid stabilization device (HSD), inter-spinous process device (IPD), Oswestry Disability Index (ODI), length of hospital stay (LOS)

References

[1] Ishimoto Y, Yoshimura N, Muraki S, et al. (2012) Prevalence of symptomatic lumbar spinal stenosis and its association with physical performance in a population-based cohort in Japan: the Wakayama Spine Study[J]. Osteoarthritis & Cartilage, 20(10).

[2] Deyo RA, Mirza SK, Martin BI, et al. (2010) Trends, major medical complications, and charges associated with surgery for lumbar spinal stenosis in older adults. JAMA; 303:1259–65.

[3] Zaina F, Tomkins-Lane C, Carragee E, et al. (2016) Surgical versus non-surgical treatment for lumbar spinal stenosis[J]. Cochrane Database of Systematic Reviews,, 1(14):CD010264.

[4] Lee, N, Shin, D. A, Kim, K. N. et al. (2016) Paradoxical Radiographic Changes of Coflex Interspinous Device with Minimum 2-Year Follow-Up in Lumbar Spinal Stenosis[J]. World Neurosurgery. 85:177–184.

[5] CHEN Zhi ming， MA Hua song， ZHAO Ji e， et al (2010). Three—dimensional finite element analyses of unilateral pedicles crews fixation in lumbar spine. Chinese Journal of Spine and SpinalCord. 20(8):684688

[6] YANG Ming-jie ZENG Cheng LI Li-jun, et al. (2018) Establishment of 3D-finite element model for analysis of biomechanical stability of extraforaminal or transforaminal lumbar interbody fusion [J]. Journal of Tongji University(Medical Science), 39(3):41–47.

[7] YANG Ming-jie ZENG Cheng LI Li-jun, et al. (2018) Biomechanical analyses of extraforaminal lumbar interbody fusion [J]. Journal of Tongji University (Medical Science), 39(4):51–55.

[8] YAN Jia-zhi WU Zhi-hong WANG Xue-song, et al. (2009) Finite element analysis on stress change of lumbar spine [J]. NATIONAL MEDICAL JOURNAL OF CHINA, 89(17):1162–1165.

[9] Hao jian, Piao zhe, Li jihai, et al. (2012) Establishment of a normal human lumbar three-dimensional finite element model based on CT image and reverse engineering methods [J]. JOURNAL OF CLINICAL REHABILITATIVE TISSUE ENGINEERING RESEARCH, 16(4):593–596.

[10] Putzer M, Auer S, Malpica W, et al. (2016) A numerical study to determine the effect of ligament stiffness on kinematics of the lumbar spine during flexion[J]. BMC Musculoskeletal Disorders, 17(1):95.

[11] Shin JK, Lim BY, Goh TS, et al. (2018) Effect of the screw type (S2-alar-iliac and iliac), screw length, and screw head angle on the risk of screw and adjacent bone failures after a spinopelvic fixation technique: A finite element analysis. PLoS One. Aug 16;13(8).

[12] Zhao Y, Li J, Wang D, Liu Y, Tan J, Zhang S. (2012) Comparison of stability of two kinds of sacro-iliac screws in the fixation of bilateral sacral fractures in a finite element model. Injury; 43(4):490–4.

[13] Driscoll M, Aubin C-E, Moreau A, Parent S. (2011) Biomechanical comparison of fusionless growth modulation corrective techniques in pediatric scoliosis. Medical & biological engineering & computing; 49(12):1437–45.

[14] Schultz AB,Warwich DN,Berkson MH,et al. (1979) Mechanical properties of human lumbar spine motion segments-Part I: responses in flexion,extension,lateral bending,and torsion[J]. J Biomech Eng, 101:46.

[15] Andersson GB,Schultz AB. (1979) Effects of fluid injection on mechanical properties of intervertebral discs[J]. J Biomech,12:453.

[16] Tencer AF,Ahmed AM,Burke DL. (1982) Some static mechanical properties of the lumbar intervertebral joint, intact and injured[J]. J Biomech Eng,104(3):193.

[17] Wilke HJ,Neef P,Caimi M,et al. (1999) New in vivo measurements of pressures in the intervertebral disc in daily life[J]. Spine (Phila Pa 1976), 24(8):755.

[18] Virgin WJ. (1951) Experimental investigations into the physical properties of the intervertebral disc[J]. J Bone Joint Surg Br, 33-B(4):607.

[19] Huang yufeng et al, (2015) Establishment of normal lumbosacral vertebral three-dimensional finite element [J]. Orthopaedic Biomechanics Materials and Clinical Study.15(5):16.

[20] Vadapalli S，Sairyo K，Goel VK，et al. (2006) Biomechanical rationale for using polyetheretherketone （PEEK） spacers for lumbar interbody fusion：a finite element study[J]. Spine. 31（26）：E992-E998.

[21] Park P, Garton H J, Gala V C, et al. (2004) Adjacent Segment Disease after Lumbar or Lumbosacral Fusion: Review of the Literature[J]. Spine. 29(17):1938–1944.

[22] Okuda S, Oda T, Yamasaki R, et al. (2014) Repeated adjacent-segment degeneration after posterior lumbar interbody fusion[J]. Journal of Neurosurgery: Spine, 20(5):538–541.

[23] P. Khoueir, K. A. Kim, and M. Y. Wang. (2007) “Classification of posterior dynamic stabilization devices,” Neurosurgical Focus, vol.22,articleE3,

[24] Po-Hsin C, Hsi-Hsien L, An H S, et al. (2017) Could the Topping-Off Technique Be the Preventive Strategy against Adjacent Segment Disease after Pedicle Screw-Based Fusion in Lumbar Degenerative Diseases? A Systematic Review[J]. BioMed Research International, 2017:1–13.

[25] Lu K, Liliang P C, Wang H K, et al. (2015) Reduction in adjacent-segment degeneration after multilevel posterior lumbar interbody fusion with proximal DIAM implantation[J]. Journal of neurosurgery. Spine, 23(2):1–7.

[26] Qu SD, Hai Y, Su QJ, Qu SP. (2015) Finite element analysis of the refined interspinous dynamic system based on Coflex. Zhongguo Zuzhi Gongcheng Yanjiu.19(22):3571–3578

[27] Li A M, Li X, Yang Z. (2017) Decompression and coflex interlaminar stabilisation compared with conventional surgical procedures for lumbar spinal stenosis: A systematic review and meta-analysis[J]. International Journal of Surgery, 40(Complete):60–67.

[28] Yuan W, Su Q J, Liu T, et al. (2017) Evaluation of Coflex interspinous stabilization following decompression compared with decompression and posterior lumbar interbody fusion for the treatment of lumbar degenerative disease: A minimum 5-year follow-up study[J]. Journal of Clinical Neuroscience, 35:24–29.

Tables

Table 1.1（1）Material properties of the tissues and implants.

*：Intervertebral fusion cage was used to supplement the k-rod internal fixation system.

Table 1.1（2）Ligament stiffness matrix in N/mm with the according ranges of the intervals.

The abbreviations are as follows: ALL anterior longitudinal, PLL posterior longitudinal, ISL intraspinous, SSL supraspinous, LF flavum, IT intertransverse

Table 1.5 Number of nodes and elements in the models.

Table 3.2(1) ROM and torsional rigidity in different motions of this study

Table 3.2(2) Comparison of ROM and torsional rigidity  (N.m/°)

Table 3.3 vonmises stress in the intervertebral disc of L3-L4 comparison with other literatures (Mpa)

Table 3.4 Analysis of the state and influence on the interpretation by densities of the mesh.

Table 4(1) Finite element analysis of L3-L4 lumbar ROM (°)

Table 4(2) L3-L4 lumbar torsional rigidity (N*m/°)

Table 5 Maximum von Mises equivalent stress (Mpa) of L3/L4 intervertebral disc

Table 6 Maximum von Mises equivalent stress (Mpa) of implants