- 1Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI, United States
- 2School of Mechanical Engineering, Vellore Institute of Technology, Chennai, India
Introduction: Degenerative cervical myelopathy (DCM) is a leading cause of spinal cord dysfunction and often requires surgical decompression and instrumentation. Anterior cervical discectomy and fusion (ACDF) is the most commonly performed procedure, but the biomechanical response of instrumented cervical spines under rear-impact loading remains unclear.
Methods: A validated finite element model of the head–neck complex with active musculature, spinal cord, and three instrumented constructs was used to simulate a low-speed rear impact corresponding to a 2.6 m/s change in velocity. The following surgical conditions were modeled at C4–C7: ACDF, posterior cervical laminectomy and fusion (PCLF), and laminoplasty (LP). Spinal cord stress and strain, ligament strain (including anterior longitudinal and capsular ligaments), and implant stress were evaluated at the index, superior, and inferior segments.
Results: Compared to ACDF, both PCLF and LP reduced spinal cord stress and strain at the index level, with PCLF showing the greatest reduction. LP preserved motion but increased facet capsular ligament strain at the index and adjacent levels. Anterior longitudinal ligament strains remained below 36%, within the reported 40–45% failure range, while PCLF screws exceeded the 900 MPa yield strength of Ti–6Al–4V at all levels except C4, indicating a reduced hardware safety margin. ACDF and LP implants remained within safe stress limits.
Conclusion: Under rear-impact loading, PCLF reduced spinal cord stress and strain relative to ACDF but increased hardware stresses, whereas LP preserved motion at the cost of higher facet capsular ligament strain. These findings suggest that PCLF may be preferred when minimizing spinal cord loading is a priority, while LP may be selected when motion preservation is critical, with awareness of construct-specific trade-offs in ligament and implant loading.
Introduction
Degenerative cervical myelopathy (DCM) is the leading cause of spinal cord dysfunction in adults, resulting from chronic cervical spinal cord compression due to degenerative changes (Smith et al., 2021). This condition significantly impacts public health, with its progressive nature reducing quality of life. Epidemiological data indicate an increasing prevalence, with reporting that 62% of individuals over 40 exhibit degenerative cervical changes, underscoring the need for effective management. Traffic accidents further exacerbate this burden, contributing to over 1.25 million global deaths and 20–50 million non-fatal injuries annually, including a substantial number of cervical spine traumas (World Health Organization, 2023). Neck injuries account for 66% of occupant trauma in automotive collisions, demonstrating the intersection between DCM and post-operative vulnerabilities (Guo et al., 2024).
Surgical decompression with instrumentation remains the primary intervention for DCM, aiming to restore mechanical stability and alleviate symptoms. ACDF is the most frequently performed procedure, achieving fusion rates exceeding 95%. In the United States, approximately 150,000 cervical spine surgeries are performed annually using the ACDF, compared to 30,000 using the posterior approach (Epstein, 2019). The number of anterior cervical surgeries is projected to increase by 13% between 2020 and 2040, while posterior cervical procedures are expected to rise by 19% in response to an aging population and increasing case complexity (Mohamad Radzi et al., 2021). Despite the efficacy of these interventions, rigid instrumentation alters spinal biomechanics, potentially increasing vulnerability under dynamic loading conditions such as automotive impacts.
Rear impact collisions account for 57.7% of global car accident cases (Razzak et al., 2022). Whiplash injuries, although classified as minor trauma, can result in capsular ligament damage, anterior longitudinal ligament tears, disc fiber disruption, and direct spinal cord trauma (Panjabi et al., 1998). Guo et al. (2024) report that 50% of neck injury come from automotive crashes, with frontal impacts alone causing severe spinal fractures in 5,592 ± 1,170 cases annually (Guo et al., 2024; Parenteau and Viano, 2014). Patients with prior cervical instrumentation may experience catastrophic secondary injuries, yet the biomechanical response of instrumented spine to rear impacts remains unclear. White et al. (2016) compared arthroplasty and fusion in frontal impacts, omitting spinal cord injury (White et al., 2016). Huang et al. (2018) analyzed anterior longitudinal ligament strain post-fusion during rear impact but did not assess posterior approaches or spinal cord injury (Huang et al., 2018). FE model have been enhanced the understanding of cervical spine biomechanics under automotive impact; however, most studies evaluated non-instrumented spines, the limited instrumented spine models did not include the spinal cord. Guo et al. (2024) evaluated ACDF under frontal impacts, increased vertebral destruction with additional fusion levels but did not investigate rear impact conditions or alternative surgical techniques such as PCLF and LP (Guo et al., 2024). The absence of comprehensive analyses examining instrumented spines, including spinal cord, under rear impact loading limits understanding of post-operative trauma risks and stability of cervical spine. Under dynamic rear-impact loading, muscle activation can reduce capsular ligament (CL) strain from 28% to 13% in 7 g simulations (Fice et al., 2011). Facet joints are a primary source of chronic post-whiplash neck pain in about 60% of cases, and cadaver studies have reported CL strains exceeding a sub-failure threshold of 35% at 5 g and 10.5 g (Panjabi et al., 1998; Panjabi et al., 2004; Lord et al., 1996). For the anterior longitudinal ligament (ALL), failure-related strain thresholds of 40%–45% have been reported (Henderson et al., 2005). Instrumented constructs are evaluated by comparing implant stresses to the 900 MPa yield strength of Ti-6Al-4V (Kulkarni et al., 2024; Tadepalli and Grosland, 2025).
This study addresses these gaps by utilizing a validated FE model of the head-neck, incorporating active musculature, disc, ligaments, and the spinal cord, to compare biomechanical outcomes of ACDF, PCLF, and LP following rear impact trauma. A low-speed rear impact (2.6 m/s velocity change) is simulated to quantify spinal cord stress, ligament strain, and instrumentation stress across these surgical interventions. The findings provide critical data on the stability and injury risks of surgically altered cervical spines, informing surgical decision-making and occupant protection strategies in prevalent rear impact trauma.
Materials and methods
We used our validated 3D head-neck FE model (John et al., 2019). This model is developed using a multi-block hexahedral meshing technique (John JD. et al., 2017). It encompasses the osteoligamentous cervical spine column from C1 to T1 (Figure 1A), integrating vertebral bodies, posterior elements, intervertebral disks, facet joints with their cartilage, and ligaments. The subaxial cervical spinal column model also featured the spinal cord, active musculature, and a rigid head (Vedantam et al., 2023a). The vertebral components are detailed with layers representing cortical and cancellous bone, topped and bottomed by superior and inferior endplates. For intervertebral disks, we distinguished between annulus and nucleus pulposus, each modeled with specific elements to capture their unique material properties (Jobin JD. et al., 2017). The annulus fibrosus represented by quadratic shell elements, reflecting its nonlinear behavior, while the nucleus pulposus is modeled with solid hexahedral elements. This model is comprehensive, featuring all critical ligaments, the spinal cord, and cervical muscles using the Hill-type muscle model. This includes representation of 23 pairs of cervical spine muscles, with placement informed by FE study (John et al., 2019). We also include muscle activation through the active Hill-type muscle, contributing to a dynamic understanding of muscular interactions within the cervical spine (Harinathan et al., 2023). The acceleration time-history (2.6 m/s over 120 ms) shown in Figure 2. Muscle activation for flexors muscle group is activated from 74 ms to 120 ms after impact, and extensor muscle group activation is set to 70% of the flexor level. The spinal cord-spinal canal interface as frictionless surface-to-surface contact, representing the lubricated interface created by cerebrospinal fluid. Table 1 Contains material properties of the cervical spine, spinal cord and surgical construct components.
Figure 1. (A) FE model of head and neck (B) ACDF at C4-C7 (C) PCLF at C4-C6 and fusion C3-C7 (D) LP at C4-C6.
The ACDF FE modelling mimics the surgical insertion of a graft between the C4-C5, C5-C6, and C6-C7 vertebral bodies, aligning with the procedural norms. Utilizing CATIA V6 software (Dassault systems, USA). Figure 1B shows, 3D models of the anterior cervical plate and its screws, which have variable angles. These screws, standardized to a length of 18 mm and a diameter of 3 mm, were integrated into our FE model (Vedantam et al., 2023a). To ensure the graft’s stability, constraints mimicking the surgical attachment to the adjacent vertebral bodies were applied. The material properties of the surgical instruments were consistent with current research findings. Contact was simulated by using a ‘tied’ constraint for the screw-vertebra and an automatic surface-to-surface contact for the screw-plate and plate-vertebra. The ACDF’s biomechanical impact was demonstrated by modeling a 1 mm ventral CSF column at the decompressed sites, reflecting the surgical outcome.
In the PCLF simulation, excision was performed on the C5-C6 segment’s spinous process, lamina, and ligamentum flavum, extending to adjacent segments. Using CATIA V6, (Dassault systems, USA). We simulated the fusion from C4 to C7 (Figure 1C), incorporating a custom-modeled titanium rod with a variable curvature to align with the lateral mass screws. Contact was simulated by using a ‘tied’ constraint for the screw-vertebra, screw-rod and an automatic surface-to-surface contact for the rod-vertebra (Vedantam et al., 2023b). Post-surgical biomechanical changes were represented by a dorsal shift of the spinal cord by 2 mm at the treated segments, mirroring the expected post-decompression displacement (Ashana et al., 2021).
For the open-door LP model, a full-thickness trough was created on one side of the lamina at C4–C6, and a partial-thickness trough was created on the contralateral side to function as a hinge (Figure 1D). Following the lamina and ligament modifications, the laminae were rotated towards the hinge, simulating the surgical procedure (Vedantam et al., 2023b). The implementation of a titanium plate, meshed using hexahedral elements through ANSA software (BETA CAE systems, USA), secured the LP, Contact was simulated by using a ‘tied’ constraint for the screw-vertebra and an automatic surface-to-surface contact for the plate-vertebra (Vedantam et al., 2023b).
Our model underwent a rear impact simulation, where linear acceleration was applied to the T1 vertebrae, inducing a velocity change of 2.6 m/s (5.82 mph) in the Gx direction (Figure 2), aligning with real-world data derived from cadaver experiments in our facility (Stemper et al., 2004). The model incorporated a 4.06 kg head mass, represented as a point mass. Head positioning matched the experimental setup, with the Frankfort plane horizontal and the occipital condyles directly above T1. Passive and active muscle structures were included in the model, the applied acceleration pulse spanning 120 ms to capture the dynamic response. The FE analyses were conducted using the LS DYNA R10.0.0 solver, MATLAB for postprocessing simulations.
We quantitatively assessed key parameters to understand the biomechanics of an instrumented spine subjected to a low-speed rear impact. Specifically, we measured spinal cord stress and strain, range of motion, capsular ligament strain, anterior longitudinal ligament strain and instrument stress. Interventions included both anterior (ACDF) and posterior (PCLF, LP) approaches, and these metrics were compared with the ACDF The index level, which averages the operated levels (C4-C5, C5-C6, C6-C7); the superior level (C3-C4); and the inferior level (C6-C7).
Results
Intact head-neck model was validated using range of motions derived from human cadaver experiments, focusing on segments ranging from C2–C3 to C6–C7. This validation process confirmed that the FE model response to loading aligned with experimental data (John et al., 2019; Purushothaman and Yoganandan, 2022). In the comparative analysis of biomechanical outcomes following ACDF, PCLF, and LP, we assessed the range of motion (ROM) at three critical spinal levels: superior, index, and inferior. The index level represents the average outcomes across the three decompressed regions.
At the superior level, PCLF showed a slight ROM reduction (−8%) compared to ACDF’s significant decrease (−59%). LP exhibited a lesser ROM reduction. The index level observed a substantial increase in ROM for both PCLF (70%) and LP (302%) relative to ACDF, indicating increased mobility at the decompressed levels. At the inferior level, ROM increases were noted for PCLF (88%) and LP (35%) compared to ACDF (Figure 3). These findings suggest that PCLF and LP increases ROM, particularly at the index levels, with LP offering a substantial increases in ROM. PCLF and LP allow for increased ROM at the index levels, with LP showing a profound increase in ROM, thereby highlighting its potential for preserving motion post-surgery.
Superior Level: PCLF and LP showed a reduction in stress by −71% and −39% respectively compared to ACDF. Index Level: At the decompressed levels, there was a notable decrease in stress with PCLF (−81%) and LP (−51%) when compared with ACDF (Figure 4). Inferior Level: Both PCLF and LP exhibited lower stress reductions at −41% and −29% respectively. Figure 5 illustrates the spatial distribution of stress across the spinal cord in the mid-sagittal plane. When a peak scale value is set to 25 kPa, ACDF reaches peak stress in 90 ms, while LP reaches at 120 ms, and PCLF stays within the scale limit through out the event.
Superior Level: Strain reductions were −70% for PCLF and −14% for LP. Index Level: PCLF and LP demonstrated decreases in strains by −57% and −38%, respectively. Inferior Level: the reductions were −14% for PCLF and −27% for LP (Figure 6). The data indicate that both PCLF and LP are effective in reducing spinal cord stress and strain compared to ACDF, with PCLF showing more pronounced reductions in von Mises stress, particularly at the index and superior levels. Figure 7 illustrates the spatial distribution of the strain across the spinal cord in the mid-sagittal plane. The peak scale value is set to 10%. ACDF reaches peak strain at 90 ms, especially at adjacent levels, while LP reaches at 120 ms, and PCLF stays within the scale limit.
Superior Level: CL strain decreased by 71% in PCLF but increased by 103% in LP. At the index level, CL strain increased to 4 times in PCLF and 6 times in LP. At the inferior adjacent level, both procedures resulted in 2 times increase in CL strain (Figure 8).
Superior Level: ALL strain slightly increased by 34% in PCLF and decreased by 1% in LP. Inferior Level: Substantial increases were seen, with PCLF showing a 151% increase and LP a 56% increase (Figure 9).
Instrument stress was evaluated in terms of peak von Mises stress. For ACDF, screw stresses ranged from 145 MPa at C6 to 398 MPa at C7. PCLF lateral mass screws displayed stresses from 210 MPa at C4 to 3,330 MPa at C3. In LP, lateral mass screw stresses were lower, from 99 MPa at C4 to 171 MPa at C6, and lamina screws showed 125 MPa at C5 to 251 MPa at C4. For plates/rods, the plate exhibited a peak stress of 587 MPa, the PCLF rod showed a much lower stress of 57 MPa, and the LP plates peaked at 527 MPa at C4.
Discussion
In this study, we conducted simulations on three different cervical surgical interventions: ACDF, PCLF, and LP, simulated under the rear impact by applying a linear acceleration corresponding to a change in velocity of 2.6 m/s at T1 vertebrae, comparing their biomechanical outcomes with ACDF. The intact head-neck model was validated using range of motion data derived from human cadaver experiments for segments C2–C3 to C7–T1, ensuring the FE model’s response under load accurately reflected experimental data (Panjabi et al., 2004; Stemper et al., 2004). The rationale behind comparing between anterior approaches and posterior approaches, this comparison allows for an understanding of how different approaches alter spinal cord and column biomechanics under the rear impact.
For both fusion techniques, ACDF and PCLF, there was a substantial reduction in ROM at the index level. In contrast, LP, as a motion-preserving intervention, showed a higher ROM, as shown in cadaver studies (Huang et al., 2023). The ROM at the inferior level was higher than at the superior adjacent levels, consistent with previous in vivo and in vitro studies indicating hyperextension in the lower cervical spine during the initial phase of rear impact (John et al., 2019). The timing of the maximum S-curve, defined as the moment when C2–C3 reaches peak flexion, was observed at 40-45 ms across all surgical interventions, followed by an extension phase (John et al., 2019).
In rear impact scenarios characterized by extension movements, during extension the posterior elements of the spine are predominantly engaged. PCLF uses lateral mass screws and rods to effectively immobilize the spine, thus mitigating forces associated with such impacts (Graham et al., 1996; Lindsey and Miclau, 2000; Kiapour et al., 2022). In contrast, ACDF involves anterior fusion, which does not completely restrict extension movements, potentially increasing spinal cord stress and strain. While the range of motion is directly proportional to spinal cord stress and strain (Breig et al., 1966; Henderson et al., 2005). However, our study reveals that LP, despite allowing a higher range of motion than ACDF, results in minimal spinal cord stress and strain. This outcome is attributed to the posterior approach’s capacity to expand the fluid volume behind the spinal cord, acting as a buffer against impact forces (Ashana et al., 2021). These points demonstrate the biomechanical benefits of posterior surgical interventions in maintaining spinal integrity and reducing injury risk during rear impacts.
Muscle activation played a role, as evidenced by a study from Fice et al., which showed a reduction in CL strain from 28% to 13% with active muscles in a 7 g rear impact scenario (Fice et al., 2011). Study highlights the potential of predict soft tissue deformation during impacts, which could enhance our understanding of pain response mechanisms (Corrales and Cronin, 2021). Given that cervical facet joints are implicated as a primary source of chronic neck pain post-whiplash in 60% of cases (Lord et al., 1996). Previous cadaver studies have documented capsular ligament (CL) strains surpassing the 35% threshold, considered indicative of sub-traumatic failure, at load levels of 5 g and 10.5 g (Panjabi et al., 1998; Panjabi et al., 2004). In contrast, our investigation observed an average CL strain of 32.4% in PCLF procedures, which remains below this threshold. However, LP (LP) procedures at the index levels yielded a significantly higher average CL strain of 44.7%, thereby exceeding the sub-traumatic failure threshold of 35%. For the ALL, which is crucial in extension movements, the study by Yoganandan et al. 2000 suggests a strain threshold range of 40%-45% (Yoganandan et al., 2000). The average ALL strains were below this failure threshold, at under 36%, indicating a lower risk of damage. sOverall, posterior constructs reduced spinal cord stress and strain compared with ACDF during rear-impact loading. By increasing the posterior canal diameter, PCLF and LP permit dorsal spinal cord translation during extension-dominated motion, which reduces ventral contact and peak cord responses at the index levels. However, open-door laminoplasty produced the highest capsular ligament strains because the unilateral trough and hinge disrupt the posterior tension band and alter facet kinematics, leading to greater posterior translation, extension, and relative motion across the facets ligaments.
Instrument failure was evaluated by comparing the observed stress values to the yield strength of titanium alloy (Ti-6Al-4V), which is 900 MPa (Kulkarni et al., 2024; Tadepalli et al., 2011). In PCLF, lateral mass screws exhibited stress levels that surpassed the yield strength at all levels except C4, indicating a pronounced risk of failure for these screws under dynamic loads. However, the rod used in PCLF remained well within the safety limit. For ACDF and LP, the instruments exhibited stress levels below the yield threshold.
This study is limited by the use of low-speed rear impact simulations and lack of white and gray matter-specific tissue tension due to the nonavailability of human spinal cord white and gray matter tissue properties. Instrumentation was modeled as fully fused to bone using tied screw–bone and plate–bone interfaces, representing an ideal healed construct. This does not capture early postoperative states with micromotion which could increase motion and alter implant stresses and should be explored in future sensitivity analyses. Although we validated our FE model with experimental data, the interpretation of our results must consider the limitations of FE models for injury simulation. Our findings are based on a single male cervical spine model with neutral lordosis and average tissue properties, which limits generalizability. Variations in spinal curvature, vertebral geometry, muscle morphology, and sex-related anatomy may alter the quantitative response. Future work should use subject-specific models or parametric variations (e.g., different lordosis angles, bone density, and muscle cross-sectional area) to capture inter-individual variability. The outputs for the inferior adjacent segment may be unreliable due to the boundary condition effect, specifically the restriction of six degrees of freedom at the T1 vertebra. This study is limited to a single impact scenario (change in velocity 2.6 m/s) and assumes fully fused bone–implant interfaces. Future work should explore a range of impact speeds and directions, as well as partial bonding or screw loosening at the bone–implant interface (e.g., frictional or degraded contacts), to assess sensitivity to these factors and better represent early postoperative behavior. PCLF lateral mass screws experienced peak von Mises stresses above the assumed Ti-6Al-4V yield strength of 900 MPa. Screws were modeled as linear elastic and fully bonded to bone, so post-yield plasticity and fatigue were not represented; stresses above 900 MPa therefore mark regions where yielding would be expected rather than true elastic behavior beyond yield. This reduced safety margin for the PCLF hardware under extension-dominated rear-impact loading should be interpreted as a conservative indication of potential screw yielding or fracture, and future studies should incorporate elasto-plastic screw models and fatigue assessment to better quantify hardware-related risk. The neck position was neutral at the time of rear impact; however, neck rotation, flexion, or extension postures may alter spinal cord responses. True validation of the spinal cord responses is not feasible in vivo, and in vitro, cadaver studies will not be able to measure spinal cord stress or accurately recreate muscle activation. Global neck injury criteria commonly used in rear-impact assessment (e.g., NIC) were not computed, because the present analysis focused on segmental kinematics and tissue-level responses of the spinal cord and capsular ligaments. Despite these limitations, this is the first study to measure spinal cord and column responses in instrumented spine during rear impact.
Conclusion
This study evaluated the biomechanical response of ACDF, PCLF, and LP under rear-impact loading using a validated finite element model. Compared to ACDF, posterior approaches reduced spinal cord stress and strain at the index level, with PCLF showing the greatest reduction despite similar reductions in segmental ROM for ACDF and PCLF. LP preserved the highest range of motion but resulted in elevated capsular ligament strain at the index level (44.7%) exceeding the sub-traumatic failure threshold of 35%, whereas PCLF remained below this threshold (32.4%); anterior longitudinal ligament strains remained below 36%, within the reported 40%–45% failure. Despite using a low-speed impact and long fusion, this is the first study to quantify spinal cord and implant responses in instrumented spines during rear impact. Clinically, these findings suggest that PCLF may be preferred when minimizing spinal cord loading under rear-impact conditions is a priority in rear impact scenario. These findings provide critical insights for surgical planning in trauma susceptible patients with prior cervical instrumentation.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
BH: Data curation, Formal Analysis, Investigation, Methodology, Software, Validation, Writing – original draft, Writing – review and editing. DJ: Conceptualization, Data curation, Resources, Supervision, Writing – review and editing. NY: Conceptualization, Methodology, Project administration, Resources, Writing – review and editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. Authors acknowledge the support of ANRF-SERB-CRG grant (File no CRG/2023/004727) for the work.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Ashana, A. O., Ajiboye, R. M., Sheppard, W. L., Ishmael, C. R., Cohen, J. Y., Beckett, J. S., et al. (2021). Spinal cord drift following laminoplasty versus laminectomy and fusion for cervical spondylotic myelopathy. Int. J. Spine Surg. 15 (2), 205–212. doi:10.14444/8028
Breig, A., Turnbull, I., and Hassler, O. (1966). Effects of mechanical stresses on the spinal cord in cervical spondylosis: a study on fresh cadaver material. J. Neurosurg. 25 (1), 45–56. doi:10.3171/jns.1966.25.1.0045
Corrales, M. A., and Cronin, D. S. (2021). Importance of the cervical capsular joint cartilage geometry on head and facet joint kinematics assessed in a finite element neck model. J. Biomechanics 123, 110528. doi:10.1016/j.jbiomech.2021.110528
Epstein, N. E. (2019). A review of complication rates for anterior cervical diskectomy and fusion (ACDF). Surg. Neurol. Int. 10, 100. doi:10.25259/SNI-191-2019
Fice, J. B., Cronin, D. S., and Panzer, M. B. (2011). Cervical spine model to predict capsular ligament response in rear impact. Ann. Biomed. Eng. 39 (8), 2152–2162. doi:10.1007/s10439-011-0315-4
Graham, A. W., Swank, M. L., Kinard, R. E., Lowery, G. L., and Dials, B. E. (1996). Posterior cervical arthrodesis and stabilization with a lateral mass plate. Clinical and computed tomographic evaluation of lateral mass screw placement and associated complications. Spine (Phila Pa 1976) 21 (3), 323–328. doi:10.1097/00007632-199602010-00014
Guo, L. X., Zhang, D. X., and Zhang, M. (2024). Destruction mechanism of anterior cervical discectomy and fusion in frontal impact. Med. Biol. Eng. Comput. 62, 3855–3873. doi:10.1007/s11517-024-03167-z
Harinathan, B., Jebaseelan, D., Yoganandan, N., and Vedantam, A. (2023). Effect of cervical stenosis and rate of impact on risk of spinal cord injury during whiplash injury. Spine (Phila Pa 1976). doi:10.1097/BRS.0000000000004759
Henderson, F. C., Geddes, J. F., Vaccaro, A. R., Woodard, E., Berry, K. J., and Benzel, E. C. (2005). Stretch-associated injury in cervical spondylotic myelopathy: new concept and review. Neurosurgery 56 (5), 1101–1113. doi:10.1227/01.NEU.0000157929.85251.7C
Huang, H., Nightingale, R. W., and Dang, A. B. C. (2018). Biomechanics of coupled motion in the cervical spine during simulated whiplash in patients with pre-existing cervical or lumbar spinal fusion: a finite element study. Bone and Jt. Res. 7 (1), 28–35. doi:10.1302/2046-3758.71.BJR-2017-0100.R1
Huang, X., Liu, D., Yang, Y., Qiu, H., Ma, Z., Lei, W., et al. (2023). A novel surgical technique for cervical laminoplasty in patients with multilevel cervical spondylotic myelopathy: a case report and literature review. Front. Surg. 10, 1078138. doi:10.3389/fsurg.2023.1078138/full
John, J. D., Arun, M. W. J., Yoganandan, N., and Kurpad, S. N. (2017a). Mapping block-based morphing for subject-specific spine finite element models.
John, J. D., Arun, M. W. J., and Kumar, S. (2017b). Narayan yoganandan. Cervical spine finite element model with anatomically accurate asymmetric intervertebral discs. Summer Biomechanics, Bioeng. Biotransport Conference (SB3C).
John, J. D., Saravana Kumar, G., and Yoganandan, N. (2019). Rear-impact neck whiplash: role of head inertial properties and spine morphological variations on segmental rotations. J. Biomechanical Eng. 141 (11), 111008. doi:10.1115/1.4043666
Kiapour, A., Massaad, E., Joukar, A., Hadzipasic, M., Shankar, G. M., Goel, V. K., et al. (2022). Biomechanical analysis of stand-alone lumbar interbody cages versus 360° constructs: an in vitro and finite element investigation. J. Neurosurg. Spine 36 (6), 928–936. doi:10.3171/2021.9.SPINE21558
Kulkarni, A. G., Kumar, P., Shetty, G. M., Roy, S., Manickam, P. S., Dhason, R., et al. (2024). Finite element analysis comparing the biomechanical parameters in multilevel posterior cervical instrumentation model involving lateral mass screw versus transpedicular screw fixation at the C7 vertebra. Asian Spine J. 18 (2), 163–173. doi:10.31616/asj.2023.0231
Lindsey, R. W., and Miclau, T. (2000). Posterior lateral mass plate fixation of the cervical spine. J. South Orthop. Assoc. 9 (1), 36–42.
Lord, S. M., Barnsley, L., Wallis, B. J., and Bogduk, N. (1996). Chronic cervical zygapophysial joint pain after whiplash: a placebo-controlled prevalence study. Spine 21 (15), 1737–1744. doi:10.1097/00007632-199608010-00005
Mohamad Radzi, M. A., Zainal Abidin, A. N. S., Shabadin, A., Azman, N. S., Mohd Faudzi, S. A., Paiman, N. F., et al. (2021). Analysis of injury severity involving of passenger vehicles’ occupants in frontal impact collisions. JSAEM 1 (3), 430–438. doi:10.56381/jsaem.v5i3.185
Panjabi, M. M., Cholewicki, J., Nibu, K., Grauer, J., and Vahldiek, M. (1998). Capsular ligament stretches during in vitro whiplash simulations. J. Spinal Disord. 11 (3), 227–232. doi:10.1097/00002517-199806000-00009
Panjabi, M. M., Ito, S., Pearson, A. M., and Ivancic, P. C. (2004). Injury mechanisms of the cervical intervertebral disc during simulated whiplash. Spine 29 (11), 1217–1225. doi:10.1097/00007632-200406010-00011
Parenteau, C. S., and Viano, D. C. (2014). Spinal fracture-dislocations and spinal cord injuries in motor vehicle crashes. Traffic Inj. Prev. 15 (7), 694–700. doi:10.1080/15389588.2013.867434
Purushothaman, Y., and Yoganandan, N. (2022). Gender differences in cervical spine motions and loads with head supported mass using finite element models. J. Eng. Sci. Med. Diagn Ther. 5 (4), 041004. doi:10.1115/1.4054856
Razzak, J. A., Bhatti, J., Wright, K., Nyirenda, M., Tahir, M. R., and Hyder, A. A. (2022). Improvement in trauma care for road traffic injuries: an assessment of the effect on mortality in low-income and middle-income countries. Lancet 400 (10348), 329–336. doi:10.1016/S0140-6736(22)00887-X
Smith, S. S., Stewart, M. E., Davies, B. M., and Kotter, M. R. N. (2021). The prevalence of asymptomatic and symptomatic spinal cord compression on magnetic resonance imaging: a systematic review and meta-analysis. Glob. Spine J. 11 (4), 597–607. doi:10.1177/2192568220934496
Stemper, B. D., Yoganandan, N., and Pintar, F. A. (2004). Response corridors of the human head-neck complex in rear impact. Annu. Proc. Assoc. Adv. Automot. Med. 48, 149–163.
Tadepalli, S. C., Gandhi, A. A., Fredericks, D. C., Grosland, N. M., and Smucker, J. (2011). Cervical laminoplasty construct stability: an experimental and finite element investigation. Iowa Orthop. J. 31, 207–214.
Tadepalli, S. C., and Grosland, N. M. (2025). Cervical laminoplasty construct stability: an experimental and finite element investigation.
Vedantam, A., Purushothaman, Y., Harinathan, B., Scripp, S., Budde, M. D., and Yoganandan, N. (2023a). Spinal cord stress after anterior cervical diskectomy and fusion: results from a patient-specific finite element model. Ann. Biomed. Eng. 51 (5), 1040–1051. doi:10.1007/s10439-022-03118-5
Vedantam, A., Harinathan, B., Warraich, A., Budde, M. D., and Yoganandan, N. (2023b). Differences in spinal cord biomechanics after laminectomy, laminoplasty, and laminectomy with fusion for degenerative cervical myelopathy. J. Neurosurg. Spine 39, 28–39. doi:10.3171/2023.3.SPINE2340
White, N. A., Danelson, K. A., Gayzik, F. S., Hsu, W., Powers, A., and Stitzel, J. D. (2016). Effects of cervical arthroplasty on neck response during a simulated rotary-wing aircraft impact. Int. J. Crashworthiness 21 (4), 323–337. doi:10.1080/13588265.2016.1178365
Keywords: cervical spine, spinal cord compression, spinal cord injury, finite element model, spinal cord strain, rear impact
Citation: Harinathan B, Jebaseelan D and Yoganandan N (2026) Effect of rear impact on the instrumented cervical spine: a finite element study. Front. Bioeng. Biotechnol. 13:1681077. doi: 10.3389/fbioe.2025.1681077
Received: 06 August 2025; Accepted: 01 December 2025;
Published: 07 January 2026.
Edited by:
Feng Zhu, Johns Hopkins University, United StatesReviewed by:
Fusako Sato, Japan Automobile Research Institute (JARI), JapanJianyin Lei, Taiyuan University of Technology, China
Copyright © 2026 Harinathan, Jebaseelan and Yoganandan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Davidson Jebaseelan, ZGF2aWRzb24uamRAdml0LmFjLmlu