Development, calibration and validation of impact-specific cervical spine models: A novel approach using hybrid multibody and finite-element methods

BACKGROUND AND OBJECTIVE

Spinal cord injuries can have a severe impact on athletes' or patients' lives. High axial impact scenarios like tackling and scrummaging can cause hyperflexion and buckling of the cervical spine, which is often connected with bilateral facet dislocation. Typically, finite-element (FE) or musculoskeletal models are applied to investigate these scenarios, however, they have the drawbacks of high computational cost and lack of soft tissue information, respectively. Moreover, material properties of the involved tissues are commonly tested in quasi-static conditions, which do not accurately capture the mechanical behavior during impact scenarios. Thus, the aim of this study was to develop, calibrate and validate an approach for the creation of impact-specific hybrid, rigid body - finite-element spine models for high-dynamic axial impact scenarios.

METHODS

Five porcine cervical spine models were used to replicate in-vitro experiments to calibrate stiffness and damping parameters of the intervertebral joints by matching the kinematics of the in-vitro with the in-silico experiments. Afterwards, a five-fold cross-validation was conducted. Additionally, the von Mises stress of the lumped FE-discs was investigated during impact.

RESULTS

The results of the calibration and validation of our hybrid approach agree well with the in-vitro experiments. The stress maps of the lumped FE-discs showed that the highest stress of the most superior lumped disc was located anterior while the remaining lumped discs had their maximum in the posterior portion.

CONCLUSION

Our hybrid method demonstrated the importance of impact-specific modeling. Overall, our hybrid modeling approach enhances the possibilities of identifying spine injury mechanisms by facilitating dynamic, impact-specific computational models.

Quantifying the Importance of Active Muscle Repositioning a Finite Element Neck Model in Flexion Using Kinematic, Kinetic, and Tissue-Level Responses

PURPOSE

Non-neutral neck positions are important initial conditions in impact scenarios, associated with a higher incidence of injury. Repositioning in finite element (FE) neck models is often achieved by applying external boundary conditions (BCs) to the head while constraining the first thoracic vertebra (T1). However, in vivo, neck muscles contract to achieve a desired head and neck position generating initial loads and deformations in the tissues. In the present study, a new muscle-based repositioning method was compared to traditional applied BCs using a contemporary FE neck model for forward head flexion of 30°.

METHODS

For the BC method, an external moment (2.6 Nm) was applied to the head with T1 fixed, while for the muscle-based method, the flexors and extensors were co-contracted under gravity loading to achieve the target flexion.

RESULTS

The kinematic response from muscle contraction was within 10% of the in vivo experimental data, while the BC method differed by 18%. The intervertebral disc forces from muscle contraction were agreeable with the literature (167 N compression, 12 N shear), while the BC methodology underpredicted the disc forces owing to the lack of spine compression. Correspondingly, the strains in the annulus fibrosus increased by an average of 60% across all levels due to muscle contraction compared to BC method.

CONCLUSION

The muscle repositioning method enhanced the kinetic response and subsequently led to differences in tissue-level responses compared to the conventional BC method. The improved kinematics and kinetics quantify the importance of repositioning FE neck models using active muscles to achieve non-neutral neck positions.

Smoothed particle hydrodynamics implementation to enhance vertebral fracture finite element model in a cervical spine segment under compression

Spinal cord injuries (SCIs) can arise from compression loading when a vertebra fractures and bone fragments are pushed into the spinal canal. Experimental studies have demonstrated the importance of both fracture initiation and post-fracture response in the investigation of vertebral fractures and spinal canal occlusion resulting from compression. Finite element models, such as the Global Human Body Models Consortium (GHBMC) model, focused on predicting the initiation location of fractures using element erosion to model hard tissue fracture. However, the element erosion method resulted in a loss of material and structural support during compression, which limited the ability of the model to predict the post-fracture response. The current study aimed to improve the post-fracture response by combining strain-based element erosion with smoothed particle hydrodynamics (SPH) to preserve the volume of the trabecular bone during compression fracture. The proposed implementation was evaluated using a model comprising two functional spinal units (FSUs) (C5-C6-C7) extracted from the GHBMC 50th percentile male model, and loaded under central compression. The original and enhanced models were compared to experimental force-displacement data and measured occlusion of the spinal canal. The enhanced model with SPH improved the shape and magnitude of the force-displacement response to be in good agreement with the experimental data. In contrast to the original model, the enhanced SPH model demonstrated occlusion on the same order of magnitude as reported in the experiments. The SPH implementation improved the post-fracture response by representing the damaged material post-fracture, providing structural support throughout compression loading and material flow leading to occlusion.

Effect of muscle pre-tension and pre-impact neck posture on the kinematic response of the cervical spine in simulated low-speed rear impacts

Computational human body models (HBMs) can identify potential injury pathways not easily accessible through experimental studies, such as whiplash induced injuries. However, previous computational studies investigating neck response to simulated impact conditions have neglected the effect of pre-impact neck posture and muscle pre-tension on the intervertebral kinematics and tissue-level response. The purpose of the present study was addressing this knowledge gap using a detailed neck model subjected to simulated low-acceleration rear impact conditions, towards improved intervertebral kinematics and soft tissue response for injury assessment. An improved muscle path implementation in the model enabled the modeling of muscle pre-tension using experimental muscle pre-stretch data determined from previous cadaver studies. Cadaveric neck impact tests and human volunteer tests with the corresponding cervical spine posture were simulated using a detailed neck model with the reported boundary conditions and no muscle activation. Computed intervertebral kinematics of the model with pre-tension achieved, for the first time, the S-shape behavior of the neck observed in low severity rear impacts of both cadaver and volunteer studies. The maximum first principal strain in the muscles for the model with pre-tension was 27% higher than that without pre-tension. Although, the pre-impact neck posture was updated to match the average posture reported in the experimental tests, the change in posture was generally small with only small changes in vertebral kinematics and muscle strain. This study provides a method to incorporate muscle pre-tension in HBM and quantifies the importance of pre-tension in calculating tissue-level distractions.

Enhancing the Biofidelity of an Upper Cervical Spine Finite Element Model within the Physiologic Range of Motion and Its Effect On the Full Ligamentous Neck Model Response

Contemporary finite element neck models are developed in a neutral posture; however, evaluation of injury risk for out-of-position impacts requires neck model repositioning to non-neutral postures, with much of the motion occurring in the upper cervical spine (UCS). Current neck models demonstrate a limitation in predicting the intervertebral motions within the UCS within the range of motion (ROM), while recent studies have highlighted the importance of including the tissue strains resulting from repositioning FE neck models to predict injury risk. In the current study, the ligamentous cervical spine from a contemporary neck model (GHBMC M50 v4.5) was evaluated in flexion, extension and axial rotation by applying moments from 0 to 1.5 Nm in 0.5 Nm increments, as reported in experimental studies and corresponding to the physiologic loading of the UCS. Enhancements to the UCS model were identified, including the C0-C1 joint-space and alar ligament orientation. Following geometric enhancements, an analysis was undertaken to determine the UCS ligament laxities, using a sensitivity study followed by an optimization study. The ligament laxities were optimized to UCS-level experimental data from the literature. The mean percent difference between UCS model response and experimental data improved from 55% to 23% with enhancements. The enhanced UCS model was integrated with a ligamentous cervical spine (LS) model and assessed with independent experimental data. The mean percent difference between the LS model and the experimental data improved from 46% to 35% with the integration of the enhanced UCS model.

Biomechanical Study of Cervical Endplate Removal on Subsidence and Migration in Multilevel Anterior Cervical Discectomy and Fusion

Study Design

This study compares four cervical endplate removal procedures, validated by finite element models.

Purpose

To characterize the effect of biomechanical strength and increased contact area on the maximum von Mises stress, migration, and subsidence between the cancellous bone, endplate, and implanted cage.

Overview of Literature

Anterior cervical discectomy and fusion (ACDF) has been widely used for treating patients with degenerative spondylosis. However, no direct correlations have been drawn that incorporate the impact of the contact area between the cage and the vertebra/endplate.

Methods

Model 1 (M1) was an intact C2C6 model with a 0.5 mm endplate. In model 2 (M2), a cage was implanted after removal of the C4-C5 and C5-C6 discs with preservation of the osseous endplate. In model 3 (M3), 1 mm of the osseous endplate was removed at the upper endplate. Model 4 (M4) resembles M3, except that 3 mm of the osseous endplate was removed.

Results

The range of motion (ROM) at C2C6 in the M2-M4 models was reduced by at least 9º compared to the M1 model. The von Mises stress results in the C2C3 and C3C4 interbody discs were significantly smaller in the M1 model and slightly increased in the M2-M3 and M3-M4 models. Migration and subsidence decreased from the M2-M3 model, whereas further endplate removal increased the migration and subsidence as shown in the transition from M3 to M4.

Conclusions

The M3 model had the least subsidence and migration. The ROM was higher in the M3 model than the M2 and M4 models. Endplate preparation created small stress differences in the healthy intervertebral discs above the ACDF site. A 1 mm embedding depth created the best balance of mechanical strength and contact area, resulting in the most favorable stability of the construct.

Finite element analysis of the effect of anterior dynamic plating on two-level anterior cervical discectomy fusion biomechanics

BACKGROUND

Limitations of anterior cervical discectomy and fusion (ACDF) relate to mechanical failure of the construct after recurring subsidence and migration. This study aims to evaluate the effect of the maximum rotation of variable angle screws on the range of motion (ROM), cage migration, and subsidence.

METHODS

Five finite element (FE) models were developed from a C2-C7 cervical spine model. The first model was an intact C2-C7 spine model, and the second model was an altered C2-C7 model with C4-C6 cage insertion and a 2-level static plate. The other three models were altered C2-C7 models with the same C4-C6 cage insertion and a 2-level dynamic plate.

RESULTS

ROM of C4-C6 in the static plate model was reduced by about 14º from the intact model, while only reduced by about 9^o in dynamic plate models. The maximum migration and subsidence at the cage-endplate interface in the dynamic plate models were lower than that in the static plate model under all moments. The von-Mises stress of the C3-C4 and C6-C7 discs in the dynamic plate models was lower than that in the static plate model.

CONCLUSION

Results indicate dynamic plating has promising potential (higher ROM and lower von Mises stress of discs) for stabilization in multilevel ACDF than static plate, though both dynamic plate and static plate has lower ROM than the intact model. Lower screw rotational angle has superior biomechanical performance (lower migration and subsidence) to higher rotational angle in multilevel applications regardless of loading.

Sex, Age and Stature Affects Neck Biomechanical Responses in Frontal and Rear Impacts Assessed Using Finite Element Head and Neck Models

The increased incidence of injury demonstrated in epidemiological data for the elderly population, and females compared to males, has not been fully understood in the context of the biomechanical response to impact. A contributing factor to these differences in injury risk could be the variation in geometry between young and aged persons and between males and females. In this study, a new methodology, coupling a CAD and a repositioning software, was developed to reposture an existing Finite element neck while retaining a high level of mesh quality. A 5th percentile female aged neck model (F05_75YO) and a 50th percentile male aged neck model (M50_75YO) were developed from existing young (F05_26YO and M50_26YO) neck models (Global Human Body Models Consortium v5.1). The aged neck models included an increased cervical lordosis and an increase in the facet joint angles, as reported in the literature. The young and the aged models were simulated in frontal (2, 8, and 15 g) and rear (3, 7, and 10 g) impacts. The responses were compared using head and relative facet joint kinematics, and nominal intervertebral disc shear strain. In general, the aged models predicted higher tissue deformations, although the head kinematics were similar for all models. In the frontal impact, only the M50_75YO model predicted hard tissue failure, attributed to the combined effect of the more anteriorly located head with age, when compared to the M50_26YO, and greater neck length relative to the female models. In the rear impacts, the F05_75YO model predicted higher relative facet joint shear compared to the F05_26YO, and higher relative facet joint rotation and nominal intervertebral disc strain compared to the M50_75YO. When comparing the male models, the relative facet joint kinematics predicted by the M50_26YO and M50_75YO were similar. The contrast in response between the male and female models in the rear impacts was attributed to the higher lordosis and facet angle in females compared to males. Epidemiological data reported that females were more likely to sustain Whiplash Associated Disorders in rear impacts compared to males, and that injury risk increases with age, in agreement with the findings in the present study. This study demonstrated that, although the increased lordosis and facet angle did not affect the head kinematics, changes at the tissue level were considerable (e.g., 26% higher relative facet shear in the female neck compared to the male, for rear impact) and relatable to the epidemiology. Future work will investigate tissue damage and failure through the incorporation of aged material properties and muscle activation.

A BIOMECHANICAL STUDY OF THE EFFECTS OF FLEXION ANGLE ON THE INDUCTION MECHANISM OF CERVICAL SPONDYLOSIS

Lesions in facet joints such as bone hyperplasia and degenerative changes in the intervertebral discs, can compress nerve roots and the spinal cord, leading to cervical spondylosis (CS). Lesions in these parts of the spine are commonly related to abnormal loads caused by bad posture of the cervical spine. This study aimed to understand the potential mechanical effects of load amplitude on cervical spine motion to provide a theoretical basis for the biomechanical causes of CS, and to provide a reference for preventing of the condition. In this study, a finite element model of the normal human cervical spine (C1-C7) was established and validated using an infrared motion capture system to analyze the effects of flexion angle on the stresses experienced by intervertebral discs, the anterior edge of the vertebral body, the pedicle, uncinate and facet joints. Our analysis indicated that the intervertebral disc load increased by at least 70% during the 20∘ to 45∘ flexion of the neck with 121% load increase in the vertebrae. In the intervertebral discs, the stress was largest at C4-C5, and the stress was moderate at C5-C6. These results are consistent with clinical CS prone site research. According to Wolff’s law, when bones are placed under large stresses, hyperplasia can result to allow adaptation to large loads. Increased cervical spine flexion angles caused the proliferation of bone in the above-mentioned parts of the spine and can accelerate accelerating the appearance of CS.

Importance of the cervical capsular joint cartilage geometry on head and facet joint kinematics assessed in a Finite element neck model

Finite element human neck models (NMs) aim to predict neck response and injury at the tissue level; however, contemporary models are most often assessed using global response such as head kinematics. Additionally, many NMs are developed from subject-specific imaging with limited soft tissue resolution in small structures such as the facet joints in the neck. Such details may be critical to enable NMs to predict tissue-level response. In the present study, the capsular joint cartilage (CJC) geometry in a contemporary NM was enhanced (M50-CJC) based on literature data. The M50-CJC was validated at the segment and full neck levels and assessed using relative facet joint kinematics (FJK), capsular ligament (CL) and intervertebral disc (IVD) strains, a relative vertebral rotation assessment (IV-NIC) and head kinematics in frontal and rear impact. The validation ratings at the segment level increased from 0.60 to 0.64, with improvements for modes of deformation associated with the facet joints, while no difference was noted at the head kinematic level. The improved CJC led to increased FJK rotation (188%) and IVD strain (152.2%,) attributed to the reduced facet joint gap. Further enhancements of the capsular joint representation or a link between the FJK and CL injury risk are recommended. Enhancements at the tissue level demonstrated a large effect on the IVD strain, but were not apparent in global metrics such as head kinematics. This study demonstrated that a biofidelic and detailed geometrical representation of the CJC contributes significantly to the predicted joint response, which is critical to investigate neck injury risk at the tissue level.

Vestibulocollic and Cervicocollic Muscle Reflexes in a Finite Element Neck Model During Multidirectional Impacts

Active neck musculature plays an important role in the response of the head and neck during impact and can affect the risk of injury. Finite element Human Body Models (HBM) have been proposed with open and closed-loop controllers for activation of muscle forces; however, controllers are often calibrated to specific experimental loading cases, without considering the intrinsic role of physiologic muscle reflex mechanisms under different loading conditions. This study aimed to develop a single closed-loop controller for neck muscle activation in a contemporary male HBM based on known reflex mechanisms and assess how this approach compared to current open-loop controllers across a range of impact directions and severities. Controller parameters were optimized using volunteer data and independently assessed across twelve impact conditions. The kinematics from the closed-loop controller simulations showed good average CORA rating to the experimental data (0.699) for the impacts following the ISO/TR9790 standard. Compared to previously optimized open-loop activation strategy, the average difference was less than 9%. The incorporation of the reflex mechanisms using a closed-loop controller can provide robust performance for a range of impact directions and severities, which is critical to improving HBM response under a larger spectrum of automotive impact simulations.

Cage-screw and anterior plating combination reduces the risk of micromotion and subsidence in multilevel anterior cervical discectomy and fusion-a finite element study

BACKGROUND CONTEXT

Anterior cervical discectomy and fusion (ACDF) is widely used to treat patients with spinal disorders, where the cage is a critical component to achieve satisfactory fusion results. However, it is still not clear whether a cage with screws or without screws will be the best choice for long-term fusion as the micromotion (sliding distance) and subsidence (penetration) of the cage still take place repeatedly.

PURPOSE

This study aims to examine the effect of cage-screws on the biomechanical characteristics of the human spine, implanted cage, and associate hardware by comparing the micromotion and subsidence.

STUDY DESIGN

A finite element (FE) analysis study.

METHODS

A FE model of a C3-C5 cervical spine with ACDF was developed. The spinal segment was modeled with the removal of the anterior longitudinal ligament (ALL), posterior longitudinal ligament (PLL), and discectomy was then implanted with a cage-screw system. Three models were analyzed: the first was the original spine (S1 model), the second, S2, was implanted with cages and anterior plating, and the third, S3, was implanted with a cage-screw system in addition to the anterior plate. All investigations were under 1 N•m in flexion, extension, lateral bending, and axial rotation situations.

RESULTS

Finite element analysis (FEA) demonstrated that range of motion (ROM) at C3-C4 in the S2 model was significantly reduced more than that in the S3 model, while the ROM at both C4-C5 in the S3 model was reduced more than that in the S2 model in all simulations. The ROM at C3-C5 in the S1 model was reduced by over 5° in the S2 and S3 models in all loading conditions. The micromotion and subsidence at all contacts of C3-C5 in the S3 model were lower than that in the S2 model in all flexion, extension, bending, and axial simulations. The subsidence and micromotion could be seen in the barrier area of the S2 model, while they occurred near the edge of the screw in the S3 model.

CONCLUSIONS

These results showed that the cage-screw and anterior plating combination has promising potential to reduce the risk of micromotion and subsidence of implanted cages in two or more level ACDFs.

CLINICAL SIGNIFICANCE

The use of double segmental fixation with cage-screw anterior plating combination constructs may increase the stiffness of the construct and reduce the incidence of clinical and radiographic pseudarthrosis following multilevel ACDF, which in turn, could decrease the need for revision surgeries or supplemental posterior fixation.

The importance of intervertebral disc material model on the prediction of mechanical function of the cervical spine

BACKGROUND

Linear elastic, hyperelastic, and multiphasic material constitutive models are frequently used for spinal intervertebral disc simulations. While the characteristics of each model are known, their effect on spine mechanical response requires a careful investigation. The use of advanced material models may not be applicable when material constants are not available, model convergence is unlikely, and computational time is a concern. On the other hand, poor estimations of tissue's mechanical response are likely if the spine model is oversimplified. In this study, discrepancies in load response introduced by material models will be investigated.

METHODS

Three fiber-reinforced C2-C3 disc models were developed with linear elastic, hyperelastic, and biphasic behaviors. Three different loading modes were investigated: compression, flexion and extension in quasi-static and dynamic conditions. The deformed disc height, disc fluid pressure, range of motion, and stresses were compared.

RESULTS

Results indicated that the intervertebral disc material model has a strong effect on load-sharing and disc height change when compression and flexion were applied. The predicted mechanical response of three models under extension had less discrepancy than its counterparts under flexion and compression. The fluid-solid interaction showed more relevance in dynamic than quasi-static loading conditions. The fiber-reinforced linear elastic and hyperelastic material models underestimated the load-sharing of the intervertebral disc annular collagen fibers.

CONCLUSION

This study confirmed the central role of the disc fluid pressure in spinal load-sharing and highlighted loading conditions where linear elastic and hyperelastic models predicted energy distribution different than that of the biphasic model.

Brain response of a computational head model for prescribed skull kinematics and simulated football helmet impact boundary conditions

Computational human body models (HBM) present a novel approach to predict brain response in football impact scenarios, with prescribed kinematic boundary conditions for the HBM skull typically used at present. However, computational optimization of helmets requires simulation of the coupled helmet and HBM model; which is much more complex and has not been assessed in the context of brain deformation and existing simplified approaches. In the current study, two boundary conditions and the resulting brain deformations were compared using a HBM head model: (1) a prescribed skull kinematics (PK) boundary condition using measured head kinematics from experimental impacts; and (2) a novel detailed simulation of a HBM head and neck, helmet and linear impactor (HBM-S). While lateral and rear impacts exhibited similar levels of maximum principal strain (MPS) in the brain tissue using both boundary conditions, differences were noted in the frontal orientation (at 9.3 m/s, MPS was 0.39 for PK, 0.54 for HBM-S). Importantly, both PK and HBM-S boundary conditions produced a similar distribution of MPS throughout the brain for each impact orientation considered. Within the corpus callosum and thalamus, high MPS was associated with lateral impacts and lower values with frontal and rear impacts. The good correspondence of both boundary conditions is encouraging for future optimization of helmet designs. A limitation of the PK approach is the need for experimental head kinematics data, while the HBM-S can predict brain response for varying impact conditions and helmet configurations, with potential as a tool to improve helmet protection performance.

Variations in User Implementation of the CORA Rating Metric

The CORA rating metric is frequently used in the field of injury biomechanics to compare the similarity of response time histories. However, subjectivity exists within the CORA metric in the form of user-customizable parameters that give the metric the flexibility to be used for a variety of applications. How these parameters are customized is not always reported in the literature, and it is unknown how these customizations affect the CORA scores. Therefore, the purpose of this study was to evaluate how variations in the CORA parameters affect the resulting similarity scores. A literature review was conducted to determine how the CORA parameters are commonly customized within the literature. Then, CORA scores for two datasets were calculated using the most common parameter customizations and the default parameters. Differences between the CORA scores using customized and default parameters were statistically significant for all customizations. Furthermore, most customizations produced score increases relative to the default settings. The use of standard deviation corridors and exclusion of the corridor component were found to produce the largest score differences. The observed differences demonstrated the need for researchers to exercise transparency when using customized parameters in CORA analyses.

A Novel Modelling Methodology Which Predicts the Structural Behaviour of Vertebral Bodies under Axial Impact Loading: A Finite Element and DIC Study

Cervical spine injuries (CSIs) arising from collisions are uncommon in contact sports, such as rugby union, but their consequences can be devastating. Several FE modelling approaches are available in the literature, but a fully calibrated and validated FE modelling framework for cervical spines under compressive dynamic-impact loading is still lacking and material properties are not adequately calibrated for such events. This study aimed to develop and validate a methodology for specimen-specific FE modelling of vertebral bodies under impact loading. Thirty-five (n = 35) individual vertebral bodies (VBs) were dissected from porcine spine segments, potted in bone cement and μCT scanned. A speckle pattern was applied to the anterior faces of the bones to allow digital image correlation (DIC), which monitored the surface displacements. Twenty-seven (n = 27) VBs were quasi-statically compressively tested to a load up to 10 kN from the cranial side. Specimen-specific FE models were developed for fourteen (n = 14) of the samples in this group. The material properties were optimised based on the experimental load-displacement data using a specimen-specific factor (kGSstatic) to calibrate a density to Young’s modulus relationship. The average calibration factor arising from this group was calculated (K¯GSstatic) and applied to a control group of thirteen (n = 13) samples. The resulting VB stiffnesses was compared to experimental findings. The final eight (n = 8) VBs were subjected to an impact load applied via a falling mass of 7.4kg at a velocity of 3.1ms⁻¹. Surface displacements and strains were acquired from the anterior VB surface via DIC, and the impact load was monitored with two load cells. Specimen-specific FE models were created for this dynamic group and material properties were assigned again based on the density–Young’s modulus relationship previously validated for static experiments, supplemented with an additional factor (KGSdynamic). The optimised conversion factor for quasi-static loading, K¯GSstatic, had an average of 0.033. Using this factor, the validation models presented an average numerical stiffness value 3.72% greater than the experimental one. From the dynamic loading experiments, the value for KGSdynamic was found to be 0.14, 4.2 times greater than K¯GSstatic. The average numerical stiffness was 2.3% greater than in the experiments. Almost all models presented similar stiffness variations and regions of maximum displacement to those observed via DIC. The developed FE modelling methodology allowed the creation of models which predicted both static and dynamic behaviour of VBs. Deformation patterns on the VB surfaces were acquired from the FE models and compared to DIC data, achieving high agreement. This methodology is now validated to be fully applied to create whole cervical spine models to simulate axial impact scenarios replicating rugby collision events.

Numerical investigation of the relative effect of disc bulging and ligamentum flavum hypertrophy on the mechanism of central cord syndrome

BACKGROUND

The pathogenesis of the central cord syndrome is still unclear. While there is a consensus on hyperextension as the main traumatic mechanism leading to this condition, there is yet to be consensus in studies regarding the pathological features of the spine (intervertebral disc bulging or ligamentum flavum hypertrophy) that could contribute to clinical manifestations.

METHODS

A comprehensive finite element model of the cervical spine segment and spinal cord was used to simulate high-speed hyperextension. Four stenotic cases were modelled to study the effect of ligamentum flavum hypertrophy and intervertebral disc bulging on the von Mises stress and strain.

FINDINGS

During hyperextension, the downward displacement of the ligamentum flavum and a reduction of the spinal canal diameter (up to 17%) led to a dynamic compression of the cord. Ligamentum flavum hypertrophy was associated with stress and strain (peak of 0.011 Mpa and 0.24, respectively) in the lateral corticospinal tracts, which is consistent with the histologic pattern of the central cord syndrome. Linear intervertebral disc bulging alone led to a higher stress in the anterior and posterior funiculi (peak 0.029 Mpa). Combined with hypertrophic ligamentum flavum, it further increased the stress and strain in the corticospinal tracts and in the posterior horn (peak of 0.023 Mpa and 0.35, respectively).

INTERPRETATION

The stenotic typology and geometry greatly influence stress and strain distribution resulting from hyperextension. Ligamentum flavum hypertrophy is a main feature leading to central cord syndrome.

Head and Neck Response of an Active Human Body Model and Finite Element Anthropometric Test Device During a Linear Impactor Helmet Test

It has been proposed that neck muscle activation may play a role in head response resulting from impacts in American Football. The importance of neck stiffness and active musculature in the standard linear impactor helmet test was assessed using a detailed head and neck finite element (FE) model from a current human body model (HBM) compared to a validated hybrid III head and neck FE model. The models were assessed for bare-head and helmeted impacts at three speeds (5.5, 7.4, and 9.3 m/s) and three impact orientations. The HBM head and neck was assessed without muscle activation and with a high level of muscle activation representing a braced condition. The HBM and hybrid III had an average cross-correlation rating of 0.89 for acceleration in the primary impact direction, indicating excellent correspondence regardless of muscle activation. Differences were identified in the axial head acceleration, attributed to axial neck stiffness (correlation rating of 0.45), but these differences did not have a large effect on the overall head response using existing head response metrics (head injury criteria, brain injury criteria, and head impact power). Although responses that develop over longer durations following the impact differed slightly, such as the moment at the base of the neck, this occurred later in time, and therefore, did not considerably affect the short-term head kinematics in the primary impact direction. Though muscle activation did not play a strong role in the head response for the test configurations considered, muscle activation may play a role in longer duration events.

Comparison of Two Intervertebral Disc Failure Models in a Numerical C4-C5 Trauma Model

The intervertebral disc (IVD) is essential for the mobility and stability of the spine. During flexion-distraction injuries, which are frequent at the cervical spine level, the IVD is often disrupted. Finite element studies have been done to investigate injury mechanisms and patterns at the cervical spine. However, they rarely include IVD failure model. The aim of this paper was to implement and compare two types of IVD failure models and their impact on hyperflexion and hyperflexion-compression injuries simulations. The failure models were tested on a detailed C4-C5 finite elements model. The first failure model consisted in a maximal strain model applied to the elements of the annulus and nucleus. The second failure model consisted in the implementation of a rupture plane in the middle of the IVD with a tied interface created between the two sections. This interface is defined by threshold stress values of detachment in traction and shearing. The two failure models were tested in flexion only and in flexion-compression. The model without inclusion of an IVD failure model was also tested. Loads at failure and injury patterns were reported. Both failure models produce failure loads that were consistent with experimental data. Injury patterns observed were in agreement with experimental and numerical studies. However, in flexion-compression, the rupture plane model simulation reached important energy error due to high deformations in the IVD elements. Also, without inclusion of an IVD failure model, energy error forced the end of the simulation in flexion-compression. Therefore, inclusion of IVD failure model is important since it leads to realistic results, but the maximal strain failure model is recommended.

Shear stiffness in the lower cervical spine: Effect of sequential posterior element injury

The aim of this study was to determine the effect of the posterior ligaments and facet joints on the shear stiffness of lower cervical functional spinal units in anterior, posterior, and lateral shear. Five functional spinal units were loaded in anterior, posterior, and right lateral shear up to 100 N using a custom-designed apparatus in a materials testing machine. Specimens were tested in three conditions: intact, with the posterior ligaments severed, and with the facet joints removed. There was a significant decrease in anterior stiffness in the 20-100 N load range from 186 (range: 98-327) N/mm in the intact condition to 105 (range: 78-142) N/mm in the disc-only condition (p = 0.03). Posterior stiffness between these condition decreased significantly from 134 (range: 92-182) N/mm to 119 (range: 83-181) N/mm (p = 0.03). There was no significant effect of posterior ligament removal on shear stiffness. No significant differences were found in the lateral direction or in the 0-20 N range for any direction. Under a 100-N shear load, the facet joints played a significant role in the stiffness of the cervical spine in the anterior-posterior direction, but not in the lateral direction.

Objective Evaluation of Whole Body Kinematics in a Simulated, Restrained Frontal Impact

The use of human body models as an additional data point in the evaluation of human-machine interaction requires quantitative validation. In this study a validation of the Global Human Body Models Consortium (GHBMC) average male occupant model (M50-O v. 4.5) in a restrained frontal sled test environment is presented. For vehicle passengers, frontal crash remains the most common mode, and the most common source of fatalities. A total of 55-time history traces of reaction loads and kinematics from the model were evaluated against corresponding PMHS data (n = 5). Further, the model's sensitivity to the belt path was studied by replicating two documented PMHS cases with prominent lateral and medial belt paths respectively. Results were quantitatively evaluated using open source CORA software. A tradeoff was observed; better correlation scores were achieved on gross measures (e.g. reaction loads), whereas better corridor scores were achieved on localized measures (rib deflections), indicating that subject specificity may dominate the comparison at localized anatomical regions. On an overall basis, the CORA scores were 0.68, 0.66 and 0.60 for force, body kinematics and chest wall kinematics. Belt force responses received the highest grouped CORA score of 0.85. Head and sternum kinematics earning a 0.8 and 0.7 score respectively. The model demonstrated high sensitivity to belt path, resulting in a 20-point increase in CORA score when the belt was routed closer to analogous location of data collection. The human model demonstrated overall reasonable biofidelity and sensitivity to countermeasures in frontal crash kinematics.

On the importance of retaining stresses and strains in repositioning computational biomechanical models of the cervical spine

Human body models are created in a specific posture and often repositioned and analyzed without retaining stresses that result from repositioning. For example, repositioning a human neck model within the physiological range of motion to a head-turned posture prior to an impact results in initial stresses within the tissues distracted from their neutral position. The aim of this study was to investigate the effect of repositioning on the subsequent kinetics, kinematics, and failure modes, of a lower cervical spine motion segment, to support future research at the full neck level. Repositioning was investigated for 3 modes (tension, flexion, and extension) and 3 load cases. The model was repositioned and loaded to failure in one continuous load history (case 1), or repositioned then restarted with retained stresses and loaded to failure (case 2). In case 3, the model was repositioned and then restarted in a stress-free state, representing current repositioning methods. Not retaining the repositioning stresses and strains resulted in different kinetics, kinematics, or failure modes, depending on the mode of loading. For the motion segment model, the differences were associated with the intervertebral disc fiber reorientation and load distribution, because the disc underwent the largest deformation during repositioning. This study demonstrated that repositioning led to altered response and tissue failure, which is critical for computational models intended to predict injury at the tissue level. It is recommended that stresses and strains be included and retained for subsequent analysis when repositioning a human computational neck model.