Skip to main content

Finite element analysis after rod fracture of the spinal hybrid elastic rod system

Abstract

Background

The spinal hybrid elastic (SHE) rod dynamic stabilization system can provide sufficient spine support and less adjacent segment stress. This study aimed to investigate the biomechanical effects after the internal fracture of SHE rods using finite element analysis.

Methods

A three-dimensional nonlinear finite element model was developed. The SHE rod comprises an inner nitinol stick (NS) and an outer polycarbonate urethane (PCU) shell (PS). The fracture was set at the caudal third portion of the NS, where the maximum stress occurred. The resultant intervertebral range of motion (ROM), intervertebral disc stress, facet joint contact force, screw stress, NS stress, and PCU stress were analyzed.

Results

When compared with the intact spine model, the overall trend was that the ROM, intervertebral disc stress, and facet joint force decreased in the implanted level and increased in the adjacent level. When compared with the Ns-I, the trend in the Ns-F decreased and remained nearly half effect. Except for torsion, the PS stress of the Ns-F increased because of the sharing of NS stress after the NS fracture.

Conclusions

The study concluded the biomechanical effects still afford nearly sufficient spine support and gentle adjacent segment stress after rod fracture in a worst-case scenario of the thinnest PS of the SHE rod system.

Peer Review reports

Background

Supplemental posterior spinal instrumentation is considered the gold standard for the surgical treatment of spinal degenerative diseases, deformities, and fractures [1]. There is still much debate regarding the indications and clinical outcomes of spinal arthrodesis [2]. Common fusion surgery combined decompression osteotomy limits physiological motion and increases the pressure on the adjacent spine. In recent years it has caused a significant inevitable increase in adjacent segment disease, failed back surgery syndrome and mechanical complications [3]. The spinal rod fracture is a serious instrumentation failure and often requires reoperation. It is associated with the older, greater body mass index, larger sagittal rod contour, presence of connectors and crossing thoracolumbar and lumbosacral junctions. Of note, the rod material is another major factor and the rate of rod fracture was significantly higher with cobalt chromium rods than with titanium alloy or stainless-steel rods [4].

Semi-rigid dynamic stabilization systems have been introduced to support the spine and preserve physiological functions [5]. Nevertheless, the existing dynamic systems often cause complications, including bulky implants, complex structures, difficult operation, or new materials without sufficient rigidity [6, 7]. The high rate of early complications and re-revisions made surgeons hesitant to perform non-fusion dynamic stabilizations.

The mechanical structure and material selection should be improved simultaneously for the further development of a perfect dynamic stabilization system [8, 9]. The spinal hybrid elastic (SHE) rod is a semi-rigid pedicle screw-based rod intended for use with universal clinical pedicle screws. This novel construct was created using an inner semi-rigid nitinol stick (NS) and outer flexible polycarbonate urethane (PCU) shell (PS). The SHE rod dynamic stabilization system is used to support the spine and reduce loading across the adjacent segment to limit degeneration while preserving motion.

However, all materials are subjected to fatigue failure due to cyclic loads [10]. To mimic extreme failure conditions after implantation, the rod was set to a rigid fracture. This study aimed to investigate the biomechanical effects on the conditions after the internal fracture of a SHE rod. The hypothesis is that the internal fracture of the SHE rods can still provide nearly sufficient spine support and gentle adjacent segment stress.

Materials and methods

Intact lumbar spine model (INT)

This study used a three-dimensional, nonlinear finite element spinal model validated in our previous studies, constructed using ANSYS 14.5 (ANSYS Inc., Canonsburg, PA, USA). A detailed description of the material properties of an INT has been reported [11,12,13]. The INT includes the vertebrae, intervertebral discs, endplates, posterior bony elements, and all seven ligaments. These intervertebral discs comprised a ground substance, the hyperelastic annulus fibrosus, and the incompressible nucleus pulposus, with 12 double cross-linked fiber layers embedded in the ground substance. The facet joint and the annulus fibrosus were constructed with nonlinear material properties in this model. The facet joint was treated as having a non-contact behavior, and the friction coefficient was set at 0. The initial gap between a pair of facet surfaces was kept 0.5 mm. The parameters of annulus fibrosus were selected as the Young’s modulus of 5.36 MPa and the Poisson’s ratio of 0.45. The nonlinear annulus ground substance was simulated using a hyperelastic Mooney-Rivlin solid model with material constants of C10 = 0.42 and C01 = 0.105. To verify the reliability of the range of motion (ROM) and facet contact force of the INT, they were compared with those of the earlier in vitro tests. It was confirmed that this INT had a stiffness like that of cadaveric lumbar spine studies [14]. The material properties of the spine are listed in Table 1.

Table 1 Material properties of the implants used in the finite element model

Implanted model

The implant system was bilaterally inserted into the L3-L4 level of the INT according to standard surgical procedures. A set of the SHE rod system comprises four conical screws without threads (diameter = 6.4 mm, length = 45 mm) and two SHE rods (diameter = 5.5 mm, length = 30 mm). The SHE rod comprises an inner semi-rigid NS (diameter = 5.5 mm, length = 25 mm) and an outer flexible PS. The material parameters of each element were obtained from previous research [15]. The interface between the pedicle screw and the bone was simulated using total bonded contact elements to provide stable support. No sliding or separation between these edges is allowed under tensile force and compressive force. The interface of the fracture surface of the rod was treated as standard contact to provide semi-rigid fixation for lumbar stability in the finite element model. This type of contact has specified compressive strength but does not resist tensile force. Standard contact elements were also simulated on the interfaces between the pedicle screw/PS and NS/PS. The fracture was set to the distal third portion of the inner NS of the SHE rods where the maximum stress occurred. The SHE rod systems with an intact NS (Ns-I) and fractured NS (Ns-F) were compared with the INT (Fig. 1).

Fig. 1
figure 1

Three finite element spine models were established: intact spine (INT), implanted with a spinal hybrid elastic (SHE) intact rod system (Ns-I), and implanted with a SHE fractured rod system (Ns-F)

Boundary and loading conditions

With the model constrained at the bottom of the L5 vertebra, the first step of loading was applying 150 N preload on the superior surface of the L1 vertebra. In the second load step, a higher pure unconstrained moment in 0.36 Nm increments was applied to ensure that the resultant ROM of the L1 to L5 vertebrae would match all four physiological motions. A load was applied with flexion 24°, extension 12.6°, torsion 18.8° and lateral bending 24.8°. The boundary load is the maximum load under convergence. A displacement control method was applied to predict adjacent segment effects after spinal implantation [16]. The cosine function method to determine the ROM used by Hsieh [17] was applied to assess the degree of each motion segment. All values were normalized with respect to intact. The resultant intervertebral ROM and stress of the intervertebral disc and facet joint contact forces were analyzed. Distortion energy theory was applied to the intervertebral discs. The von Mises stress of each model was obtained after applying torque in each direction of the model.

Results

Intervertebral ROM

At the implanted L3-L4 level, the ROM in both models decreased compared with the INT. The ROM of the Ns-I and Ns-F decreased to 40 and 66% during flexion, 38 and 64% during extension, 82 and 95% during torsion, and 43 and 76% during lateral bending, respectively. In contrast, the ROM in both models increased compared with the INT at the adjacent L2-L3 level. The ROM of the Ns-I and Ns-F increased to 119 and 111% during flexion, 115 and 109% during extension, 105 and 101% during torsion, and 117 and 108% during lateral bending, respectively. The normalized results are in Fig. 2. After the NS fracture, there was still nearly half the stabilizing effect of the Ns-F, except during torsion, compared with that of the Ns-I. The results with a scale are listed in Table 2.

Fig. 2
figure 2

Comparison of the range of motion (ROM) during flexion. The ROM of the fractured nitinol stick (Ns-F) were between those of the intact spine (INT) and intact nitinol stick (Ns-I)

Table 2 Comparison of the intervertebral range of motion (degree) in the finite element models

Intervertebral disc stresses

At the implanted L3-L4 level, the disc stress in both models decreased compared with the INT. The disc stress of the Ns-I and Ns-F decreased to 67 and 79% during flexion, 57 and 67% during extension, 81 and 96% during torsion, and 45 and 76% during lateral bending, respectively. In contrast, the disc stress in both models increased compared with the INT at the adjacent L2-L3 level. The disc stresses of the Ns-I and Ns-F increased to 127 and 115% during flexion, 114 and 108% during extension, 106 and 100% during torsion, and 124 and 110% during lateral bending, respectively (Fig. 3 and Fig. 4). Overall, the changes in the Ns-F were approximately half of those in the Ns-I (Table 3).

Fig. 3
figure 3

The von Mises stress distribution in the adjacent (L2-L3) and implanted (L3-L4) levels during flexion. The stresses of the fractured nitinol stick (Ns-F) are between those of the intact model (INT) and intact nitinol stick (Ns-I)

Fig. 4
figure 4

Comparison of the disc stress during flexion. The disc stresses of the fractured nitinol stick (Ns-F) were between those of the intact model (INT) and intact nitinol stick (Ns-I)

Table 3 Comparison of the intervertebral disc stresses (KPa) in the finite element models

Facet joint contact force

There was no facet force during flexion in any of the models. At the implanted L3-L4 level, the facet force in both models decreased compared with the INT. The facet forces of the Ns-I and Ns-F decreased to 67 and 89% during flexion, 57 and 67% during extension, 81 and 96% during torsion, and 45 and 76% during lateral bending, respectively. In contrast, the facet force in both models increased compared with the INT at the adjacent L2-L3 level. The facet forces of the Ns-I and Ns-F increased to 127 and 115% during flexion, 114 and 108% during extension, 106 and 100% during torsion, and 124 and 110% during lateral bending, respectively. Overall, the changes in the Ns-F were approximately half of those in the Ns-I (Table 4).

Table 4 Comparison of the facet joint contact forces (N) in the finite element models

Screw stresses

The screw stress of the Ns-F decreased to 66% during flexion, 72% during extension, 41% during torsion, and 60% during lateral bending compared with the Ns-I (Fig. 5). The screw stress of all the models was between 75.6 and 158 MPa, with the maximum stress recorded in the Ns-I during lateral bending (Table 5).

Fig. 5
figure 5

The von Mises stress distribution on the pedicle screw, outer PCU shell (PS) and inner nitinol stick (NS) during flexion. The PS stress of the fractured nitinol stick (Ns-F) was higher than that of intact nitinol stick (Ns-I). In contrast, the screw stress and NS stress were lower in Ns-F group. The outer PS stress was shielded by the inner NS when the NS was intact. Once the inner NS was broken, the screw and NS lost attachment in the Ns-F, and the PS took over the stress. Hence, the NS stress and screw stress decreased, and the PS stress increased

Table 5 Comparison of the maximum stress on implants (MPa) in the finite element models

PS stress

The PS stress of the Ns-F increased to 343% during flexion, 245% during extension, and 290% during lateral bending compared with the Ns-I (Fig. 5). In contrast, the PS stress decreased by 24% during torsion. The PS stress of all the models was between 4.7 and 19.7 MPa, with the maximum stress recorded in the Ns-F during lateral bending (Table 5).

NS stress

The NS stress of the Ns-F decreased to 73% during flexion, 54% during extension, 36% during torsion, and 32% during lateral bending compared with the Ns-I (Fig. 5). The NS stress of all the models was between 21.5 and 219 MPa, with the maximum stress recorded in the Ns-I during torsion (Table 5).

Discussion

This study aimed to investigate the biomechanical effects after the internal fractures of SHE rods. The SHE rod system is a pedicle screw-based construct that incorporates universal rigid fixation technology, combined with the benefits of flexible materials. However, there are intrinsic and extrinsic biomechanical risk factors for any fusion device after fatigue failure [18]. Intrinsic properties include the material type and diameter. Semi-rigid nitinol, the material of the inner SHE rod, is comparable to or superior to titanium rods with regard to biomechanical evidence [19]. Flexible PCU, the material of the outer SHE rod, decreases the stress under the same load by NS stress shielding. Flexion is the most valued physiological movement of the spine owing to its high frequency. The load-controlled findings were the major results, while the other loads were often regarded as minor references. Furthermore, the maximum stress occurred in the caudal third, especially during flexion. Therefore, the fracture was set in the caudal third portion of the NS.

The entire diameter of the SHE rod, including the PS and NS, is the same as that of clinical spinal rods. The universal 5.5-mm diameter can easily promote application and marketing. Increasing the diameter of the NS can increase the rigidity and improve the strength. Increasing the thickness of the PCU can enhance elasticity to provide sufficient support in the case of NS fractures. Currently, the fabrication process hardly makes the PS thinner than 0.5 mm. Our study is the worst-case test scenario for the thinnest PS and thickest NS after rod fracture. The findings confirmed that there was still approximately half of the biomechanical support after the NS fracture. Therefore, when the NS strength is sufficient to cover the physiological conditions, increasing the PS thickness can make the clinical application more reliable.

Other prior studies have proposed the use of Dynesys and PEEK rods to stabilize the lumbar spine dynamically and the protection of adjacent levels by the nonrigid material [5, 6]. The Dynesys consists of a polyethylene terephthalate cord and PCU spacer that are elastic to resist fracture. That is why complications associated with the Dynesys were infection, screw loosening and screw fracture [7]. However, compared to the bulky spacer of Dynesys, the implantation of SHE rods is easier and more intuitive. One concern regarding PEEK rods is their durability and fatigue fracture. The hard surface of PEEK rods is associated with similar clinical risks to metal rods. It has been shown to initiate scratching at the rod-screw interface [20]. That is why PEEK rods need larger diameter and should be locked in specific pedicle screws [6].

The extrinsic factors can be divided into bending design, notch sensitivity, and cyclic loading. All rods studied in the present model were straight, whereas there may be pre-bent rods in future clinical use. The NS enveloped by the insulated PS was locked by a titanium screw head and nut. This decreases the risk of notch defects and electrochemical corrosion on the surface of the SHE rods and prevents repeated metal strain and fatigue after cyclic loading [21]. Simulated physiological in vivo accelerated 10-year loading had no significant degradative effect on a PCU-nitinol implant [22]. Hence, the hybrid use of PCU-nitinol may lead to the development of durable spinal implants with better clinical compliance (Fig. 6).

Fig. 6
figure 6

The von Mises stress distribution on the inner intact nitinol stick (a) (Ns-I) and fractured nitinol stick (b) (Ns-F) during flexion. The NS stress markedly increased on the proximal third in Ns-F group. This phenomenon was compatible with the connecting region of the screw and rod

The outer PS stress was shielded by the inner NS when the NS was intact. Once the inner NS was broken, the screw and NS lost attachment in the Ns-F, and the PS took over the stress. Hence, the NS stress and screw stress decreased, and the PS stress increased. However, there was an exception while the PS stress decreased during torsion. This was due to the multidirectional torsion shared by the spinal ligaments, intervertebral disc, and facet joint [23]. Theoretically, the stress distributed on the screws and rods returns to the spine after a rod fracture. However, our results showed that screw stress and NS stress decreased but did not disappear. It also indirectly proved that the SHE rods bore partial support from the PS after rod fracture (Fig. 5). Hence, SHE rod fractures differ from common spinal rod fractures.

This study has some limitations. All spinal models were healthy and had no pathological properties, such as degenerative disc diseases, compression fractures, and spondylolisthesis. The INT without decompression was incompatible with the conventional SHE rod implantation surgery. The interface between the pedicle screw and the PS was modeled as always bonded, as were the PS and NS. This study aimed to evaluate the biomechanical effects on the rod, implants contacting the rod, and spinal tissues. Moreover, the interface between the pedicle screw and the bone was assumed to be a union. It was simulated to be bonded contact. Hence, an assessment of the stress distribution of the screw and bone was not reported. Finally, this model demonstrates highly idealized precise implantation and fabrication techniques.

In conclusion, in the worst-case scenario of the thinnest PS, the function remained nearly half after rod fracture of the SHE rod system.

Availability of data and materials

The raw data and materials used and analyzed during the current study are provided in a supplementary file. Further detailed datasets are available from the author JYH on reasonable request.

Abbreviations

SHE:

Spinal hybrid elastic

INT:

Intact lumbar spine model

NS:

Nitinol stick

PCU:

Polycarbonate urethane

PS:

PCU shell

Ns-I:

Intact nitinol stick

Ns-F:

Fractured nitinol stick

ROM:

Range of motion

References

  1. Levin JM, Tanenbaum JE, Steinmetz MP, Mroz TE, Overley SC. Posterolateral fusion (PLF) versus transforaminal lumbar interbody fusion (TLIF) for spondylolisthesis: a systematic review and meta-analysis. Spine J. 2018;18(6):1088–98.

    Article  Google Scholar 

  2. Yavin D, Casha S, Wiebe S, Feasby TE, Clark C, Isaacs A, et al. Lumbar fusion for degenerative disease: a systematic review and meta-analysis. Neurosurgery. 2017;80(5):701–15.

    Article  Google Scholar 

  3. Zhong ZM, Deviren V, Tay B, Burch S, Berven SH. Adjacent segment disease after instrumented fusion for adult lumbar spondylolisthesis: incidence and risk factors. Clin Neurol Neurosurg. 2017;156:29–34.

    Article  Google Scholar 

  4. Smith JS, Shaffrey E, Klineberg E, Shaffrey CI, Lafage V, Schwab FJ, et al. Prospective multicenter assessment of risk factors for rod fracture following surgery for adult spinal deformity. J Neurosurg Spine. 2014;21(6):994–1003.

    Article  Google Scholar 

  5. Ozer AF, Cevik OM, Erbulut DU, Yaman O, Senturk S, Oktenoglu T, et al. A novel modular dynamic stabilization system for the treatment of degenerative spinal pathologies. Turk Neurosurg. 2019;29(1):115–20.

    PubMed  Google Scholar 

  6. Benezech J, Garlenq B, Larroque G. Flexible stabilisation of the degenerative lumbar spine using PEEK rods. Adv Orthop. 2016;2016:7369409.

    Article  Google Scholar 

  7. Pham MH, Mehta VA, Patel NN, Jakoi AM, Hsieh PC, Liu JC, et al. Complications associated with the Dynesys dynamic stabilization system: a comprehensive review of the literature. Neurosurg Focus. 2016;40(1):E2.

    Article  Google Scholar 

  8. Kok D, Firkins PJ, Wapstra FH, Veldhuizen AG. A new lumbar posterior fixation system, the memory metal spinal system: an in-vitro mechanical evaluation. BMC Musculoskelet Disord. 2013;14:269.

    Article  Google Scholar 

  9. Kojio K, Furukawa M, Motokucho S, Shimada M, Sakai M. Structure-mechanical property relationships for polycarbonate urethane elastomers with novel soft segments. Macromolecules. 2009;42(21):8322–7.

    Article  CAS  Google Scholar 

  10. Smith JS, Shaffrey CI, Ames CP, Demakakos J, Fu KM, Keshavarzi S, et al. Assessment of symptomatic rod fracture after posterior instrumented fusion for adult spinal deformity. Neurosurgery. 2012;71:862–7.

    Article  Google Scholar 

  11. Chen SH, Zhong ZC, Chen CS, Chen WJ, Hung C. Biomechanical comparison between lumbar disc arthroplasty and fusion. Med Eng Phys. 2009;31(2):244–53.

    Article  Google Scholar 

  12. Liu CL, Zhong ZC, Shih SL, Hung CH, Lee YE, Chen CS. Influence of Dynesys system screw profile on adjacent segment and screw. J Spinal Disord Tech. 2010;23(6):410–7.

    Article  Google Scholar 

  13. Liu CL, Zhong ZC, Hsu HW, Shih SL, Wang ST, Hung CH, et al. Effect of the cord pretension of the Dynesys dynamic stabilisation system on the biomechanics of the lumbar spine: a finite element analysis. Eur Spine J. 2011;20(11):1850–8.

    Article  Google Scholar 

  14. Yamamoto I, Panjabi MM, Crisco T, Oxland T. Three-dimensional movement of the whole lumbar spine and lumbosacral joint. Spine. 1989;14:1256–60.

    Article  CAS  Google Scholar 

  15. Chevalier V, Arbab-Chirani R, Arbab-Chirani S, Calloch S. An improved model of 3-dimensional finite element analysis of mechanical behavior of endodontic instruments. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2010;109(3):111–21.

    Article  Google Scholar 

  16. Zhong ZC, Chen SH, Hung CH. Load- and displacement-controlled finite element analyses on fusion and non-fusion spinal implants. Proc Inst Mech Eng H. 2009;223(2):143–57.

    Article  Google Scholar 

  17. Hsieh YY, Kuo YJ, Chen CH, Wu LC, Chiang CJ, Lin CL. Biomechanical assessment of vertebroplasty combined with cement-augmented screw fixation for lumbar burst fractures: a finite element analysis. Appl Sci. 2020;10:2133.

    Article  CAS  Google Scholar 

  18. Barton C, Noshchenko A, Patel V, Cain C, Kleck C, Burger E. Risk factors for rod fracture after posterior correction of adult spinal deformity with osteotomy: a retrospective case-series. Scoliosis. 2015;10:30.

    Article  Google Scholar 

  19. Massey PA, Hoge S, Nelson BG, Ogden AL, Mody MG, Myers M, et al. Nitinol memory rods versus titanium rods: a biomechanical comparison of posterior spinal instrumentation in a synthetic Corpectomy model. Global Spine J. 2021;11(3):277–82.

    Article  Google Scholar 

  20. Kurtz SM, Lanman TH, Higgs G, Macdonald DW, Berven SH, Isaza JE, et al. Retrieval analysis of PEEK rods for posterior fusion and motion preservation. Eur Spine J. 2013;22(12):2752–9.

    Article  Google Scholar 

  21. Stambough JL, Genaidy AM, Huston RL, Serhan H, El-khatib F, Sabri EH. Biomechanical assessment of titanium and stainless steel posterior spinal constructs: effects of absolute/relative loading and frequency on fatigue life and determination of failure modes. J Spinal Disord. 1997;10:473–81.

    CAS  PubMed  Google Scholar 

  22. Desai M, Bakhshi R, Zhou X, Odlyha M, You Z, Seifalian AM, et al. A sutureless aortic stent-graft based on a nitinol scaffold bonded to a compliant nanocomposite polymer is durable for 10 years in a simulated in vitro model. J Endovasc Ther. 2012;19(3):415–27.

    Article  Google Scholar 

  23. Bezci SE, Klineberg EO, O'Connell GD. Effects of axial compression and rotation angle on torsional mechanical properties of bovine caudal discs. J Mech Behav Biomed Mater. 2018;77:353–9.

    Article  Google Scholar 

Download references

Acknowledgements

This study was approved by the National Taiwan University Hospital and its Jinshan branch. The authors wish to thank all the members who participated in this study.

Funding

The authors, their immediate families, and any research foundations to which they are affiliated have not received any financial payments or other benefits from any commercial entity related to the subject of this article.

Author information

Authors and Affiliations

Authors

Contributions

JYH participated in the conception and design of the study, acquisition, analysis, and interpretation of the data, and the drafting and final preparation of the manuscript. YYH participated in the conception and design of the study, acquisition, analysis, and interpretation of the data, and the drafting and final preparation of the manuscript. CSS and SMC participated in the conception and design of the study, acquisition, analysis, and interpretation of the data, and the final preparation and revision of the manuscript. JHW and PQC participated in the conception and design of the study, acquisition, analysis, and interpretation of the data, and the final preparation of the manuscript. All authors have read and approved the final manuscript, study conception and design, data analysis and interpretation, and manuscript writing.

Corresponding author

Correspondence to Yi-You Huang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supp.

Table 1. The biomechanical analysis of the three finite element models in flexion.

Additional file 2: Supp.

Table 2The biomechanical analysis of the three finite element models in extension.

Additional file 3: Supp.

Table 3. The biomechanical analysis of the three finite element models in torsion.

Additional file 4: Supp.

Table 4. The biomechanical analysis of the three finite element models in lateral bending.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hsieh, JY., Chen, CS., Chuang, SM. et al. Finite element analysis after rod fracture of the spinal hybrid elastic rod system. BMC Musculoskelet Disord 23, 816 (2022). https://doi.org/10.1186/s12891-022-05768-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12891-022-05768-x

Keywords