C4–C5 segment finite element model development,validation, and load-sharing investigation

Detailed cervical spine models are necessary to better understand cervical spine response to loading, improve our understanding of injury mechanisms, and specifically for predicting occupant response and injury in auto crash scenarios. The focus of this study was to develop a C4–C5 finite element model with accurate representations of each tissue within the segment. This model incorporates more than double the number of elements of existing models, required for accurate prediction of response. The most advanced material data available were then incorporated using appropriate nonlinear constitutive models to provide accurate predictions of response at physiological levels of loading. This tissue-scale segment model was validated against a wide variety of experimental data including different modes of loading (axial rotation, flexion, extension, lateral bending, and translation), and different load levels. In general, the predicted response of the model was within the single standard deviation response corridors for both low and high load levels. Importantly, this model demonstrates that appropriate refinement of the finite element mesh, representation at the tissue level, and sufficiently detailed material properties and constitutive models provide excellent response predictions without calibration of the model to experimental data. Load sharing between the disc, ligaments, and facet joints was investigated for various modes of loading, and the dominant load-bearing structure was found to correlate with typical anatomical injury sites for these modes of loading. The C4–C5 model forms the basis for the development of a full cervical spine model. Future studies will focus on tissue-level injury prediction and dynamic response.     

Biomechanical analysis of the anterior cervical fusion

This paper presents a biomechanical analysis of the cervical C5–C6 functional spine unit before and after the anterior cervical discectomy and fusion. The aim of this work is to study the influence of the medical procedure and its instrumentation on range of motion and stress distribution. First, a three-dimensional finite element model of the lower cervical spine is obtained from computed tomography images using a pipeline of image processing, geometric modelling and mesh generation software. Then, a finite element study of parameters' influence on motion and a stress analysis at physiological and different post-operative scenarios were made for the basic movements of the cervical spine. It was confirmed that the results were very sensitive to intervertebral disc properties. The insertion of an anterior cervical plate influenced the stress distribution at the vertebral level as well as in the bone graft. Additionally, stress values in the graft decreased when it is used together with a cage.     

Experimental flexion/extension data corridors for validation of finite element models of the young, normal cervical spine

Finite element (FE) modeling is an important tool for studying the cervical spine in normal, injured and diseased conditions. To understand the role of mechanical changes on the spine as it goes from a normal to a diseased or injured state, experimental studies are needed to establish the external response of young, normal cervical spinal segments compared to injured or degenerated cervical spinal segments under physiologic loading. It is important to differentiate injured or degenerated specimens from young, normal specimens to provide accurate experimental results necessary for the validation of FE models. This study used seven young, normal fresh adult cadaver cervical spine segments C2-T1 ranging in age from 20 to 51 years. Prior to testing, the spines were graded in three ways: specimen quality, facet degeneration and disc degeneration. Spine segments were tested in flexion/extension, and the range of loads applied to the specimens was 0.33, 0.5, 1.0, 1.5 and 2.0 Nm. These loads resulted in rotations in the direction of loading as the primary response to loading. In general, results for young, normal specimens showed greater flexibility in flexion and less flexibility in extension than results previously reported in the literature. The flexion/extension curves are asymmetric with a greater magnitude in flexion than in extension. These experimental results will be used to validate FE models of young, normal cervical spines.     

An accurate finite element model of the cervical spine under quasi-static loading

Cervical disc injury due to impact has been observed in clinical and biomechanical investigations; however, there is a lack of data that helps to elucidate the mechanisms of disc injury during these collisions. Therefore, it is necessary to understand the behavior of the cervical spine under different types of loading situations. A three dimensional finite element (FE) model for the multi-level cervical spine segment (C0-C7) was developed using computed tomography (CT) data and applied to study the internal stresses and strains of the intervertebral discs under quasi-static loading conditions. The intervertebral discs were treated as nonlinear, anisotropic and incompressible subjected to large deformations. The model accuracy was validated by comparing it with previously published experimental and numerical results for different movements. It was shown that the use of a fiber reinforced model to describe the behavior of the annulus of the discs would predict higher maximum shear strains than an isotropic one, being therefore important the use of complex constitutive models in order to be able to detect the appearance of injured zones, since those strains and stresses are supposed to be related with damage to soft tissues. Several movements were analyzed: flexion, extension and axial rotation, obtaining that the maximum shear stresses in the disc were higher for a flexo-extension movement.     

Biomechanical analysis of the anterior cervical fusion

This paper presents a biomechanical analysis of the cervical C5-C6 functional spine unit before and after the anterior cervical discectomy and fusion. The aim of this work is to study the influence of the medical procedure and its instrumentation on range of motion and stress distribution. First, a three-dimensional finite element model of the lower cervical spine is obtained from computed tomography images using a pipeline of image processing, geometric modelling and mesh generation software. Then, a finite element study of parameters' influence on motion and a stress analysis at physiological and different post-operative scenarios were made for the basic movements of the cervical spine. It was confirmed that the results were very sensitive to intervertebral disc properties. The insertion of an anterior cervical plate influenced the stress distribution at the vertebral level as well as in the bone graft. Additionally, stress values in the graft decreased when it is used together with a cage.     

Finite element analysis of moment-rotation relationships for human cervical spine

A comprehensive, geometrically accurate, nonlinear C0-C7 FE model of head and cervical spine based on the actual geometry of a human cadaver specimen was developed. The motions of each cervical vertebral level under pure moment loading of 1.0 Nm applied incrementally on the skull to simulate the movements of the head and cervical spine under flexion, tension, axial rotation and lateral bending with the inferior surface of the C7 vertebral body fully constrained were analysed. The predicted range of motion (ROM) for each motion segment were computed and compared with published experimental data. The model predicted the nonlinear moment-rotation relationship of human cervical spine. Under the same loading magnitude, the model predicted the largest rotation in extension, followed by flexion and axial rotation, and least ROM in lateral bending. The upper cervical spines are more flexible than the lower cervical levels. The motions of the two uppermost motion segments account for half (or even higher) of the whole cervical spine motion under rotational loadings. The differences in the ROMs among the lower cervical spines (C3-C7) were relatively small. The FE predicted segmental motions effectively reflect the behavior of human cervical spine and were in agreement with the experimental data. The C0-C7 FE model offers potentials for biomedical and injury studies.     

Application of a stand-alone interbody fusion cage based on a novel porous TiO2/glass ceramic--2: Biomechanical evaluation after implantation in the sheep cervical spine]

Animals are becoming more and more common as in vivo models for the human spine. Especially the sheep cervical spine is stated to be of good comparability and usefulness in the evaluation of in vivo radiological, biomechanical and histological behaviour of new bone replacement materials, implants and cages for cervical spine interbody fusion. In preceding biomechanical in vitro examinations human cervical spine specimens were tested after fusion with either a cubical stand-alone interbody fusion cage manufactured from a new porous TiO2/glass composite (Ecopore) or polymethylmethacrylate (PMMA) after discectomy. Following our first experience with the use of the new material and its influence on the primary stability after in vitro application we carried out fusions of 20 sheep cervical spines levels with either PMMA or an Ecopore-cage, and performed radiological examinations during the following 2-4 months. In this second part of the study we intended the biomechanical evaluation of the spine segments with reference to the previously determined morphological findings, like subsidence of the implants, significant increase of the kyphosis angle and degree of the bony fusion along with the interpretation of the results. 20 sheep cervical spines segments with either PMMA- or Ecopore-fusion in the levels C2/3 and C4/5 were tested, in comparison to 10 native corresponding sheep cervical spine segments. Non-destructive biomechanical testing was performed, including flexion/extension, lateral bending and axial rotation using a spine testing apparatus. Three-dimensional range of motion (ROM) was evaluated using an ultrasound measurement system. In the native spine segments C2/3 and C4/5 the ROM increased in cranio-caudal direction particulary in flexion/extension, less pronounced in lateral flexion and axial rotation (p < 0.05). The overall ROM of both tested segments was greatest in lateral flexion, reduced to 52% in flexion/extension and to 16% in axial rotation. After 2 months C2/3- and C4/5-segments with PMMA-fusion and C2/3-segments with Ecopore-interposition showed decrease of ROM in lateral flexion in comparison to the native segments, indicating increasing stiffening. However, after 4 months all operated segments, independent from level or implanted material, were stiffer than the comparable native segments. The decrease of the ROM correlated with the radiological-morphological degree of fusion. Our evaluation of the new porous TiO2/glass composite as interbody fusion cage has shown satisfactory radiological results as well as distinct biomechanical stability and fusion of the segments after 4 months in comparison to PMMA. After histological analysis of the bone-biomaterial-interface, further examinations of this biomaterial previous to an application as alternative to other customary cages in humans are necessary.     

Comparison of the anatomical morphology of cervical vertebrae between humans and macaques: related to a spinal cord injury model

Non-human primates are most suitable for generating cervical experimental models, and it is necessary to study the anatomy of the cervical spine in non-human primates when generating the models. The purpose of this study was to provide the anatomical parameters of the cervical spine and spinal cord in long-tailed macaques (Macaca fascicularis) as a basis for cervical spine-related experimental studies. Cervical spine specimens from 8 male adult subjects were scanned by micro-computed tomography, and an additional 10 live male subjects were scanned by magnetic resonance imaging. The measurements and parameters from them were compared to those of 12 male adult human subjects. Additionally, 10 live male subjects were scanned by magnetic resonance imaging, and the width and depth of the spinal cord and spinal canal and the thickness of the anterior and posterior cerebrospinal fluid were measured and compared to the relevant parameters of 10 male adult human subjects. The tendency of cervical parameters to change with segmental changes was similar between species. The vertebral body, spinal canal, and spinal cord were significantly flatter in the human subjects than in the long-tailed macaques. The cerebrospinal fluid space in the long-tailed macaques was smaller than that in the human subjects. The anatomical features of the cervical vertebrae of long-tailed macaques provide a reference for establishing a preclinical model of cervical spinal cord injury.     

Effect of muscle wrapping on model estimates of neck muscle strength

Of the computational models of the cervical spine reported in the literature, not one takes into account the changes in muscle paths due to the underlying vertebrae. Instead, all model the individual muscle paths as straight-line segments. The major aim of this study was to quantify the changes in muscle moment arm, muscle force and joint moment due to muscle wrapping in the cervical spine. Five muscles in a straight-line model of the cervical spine were wrapped around underlying vertebrae, and the results obtained from this model were compared against the original. The two models were then validated against experimental and computational data. Results show that muscle wrapping has a significant effect on muscle moment arms and therefore joint moments and should not be neglected.     

Effect of muscle wrapping on model estimates of neck muscle strength

Of the computational models of the cervical spine reported in the literature, not one takes into account the changes in muscle paths due to the underlying vertebrae. Instead, all model the individual muscle paths as straight-line segments. The major aim of this study was to quantify the changes in muscle moment arm, muscle force and joint moment due to muscle wrapping in the cervical spine. Five muscles in a straight-line model of the cervical spine were wrapped around underlying vertebrae, and the results obtained from this model were compared against the original. The two models were then validated against experimental and computational data. Results show that muscle wrapping has a significant effect on muscle moment arms and therefore joint moments and should not be neglected.     

Influence of cervical disc degeneration after posterior surgical techniques in combined flexion-extension--a nonlinear analytical study

Laminectomy and facetectomy are surgical techniques used for decompression of the cervical spinal stenosis. Recent in vitro and finite element studies have shown significant cervical spinal instability after performing these surgical techniques. However, the influence of degenerated cervical disk on the biomechanical responses of the cervical spine after these surgical techniques remains unknown. Therefore, a three-dimensional nonlinear finite element model of the human cervical spine (C2-C7) was created. Two types of disk degeneration grades were simulated. For each grade of disk degeneration, the intact as well as the two surgically altered models simulating C5 laminectomy with or without C5-C6 total facetectomies were exercised under flexion and extension. Intersegmental rotational motions, internal disk annulus, cancellous and cortical bone stresses were obtained and compared to the normal intact model. Results showed that the cervical rotational motion decreases with progressive disk degeneration. Decreases in the rotational motion due to disk degeneration were accompanied by higher cancellous and cortical bone stress. The surgically altered model showed significant increases in the rotational motions after laminectomies and facetectomies when compared to the intact model. However, the percentage increases in the rotational motions after various surgical techniques were reduced with progressive disk degeneration.     

Importance of nonlinear and multivariable flexibility coefficients in the prediction of human cervical spine motion

The flexibility matrix currently forms the basis for multibody dynamics models of cervical spine motion. While studies have aimed to determine cervical motion segment behavior, their accuracy and utility have been limited by both experimental and analytical assumptions. Flexibility terms have been primarily represented as constants despite the spine's nonlinear stiffening response. Also, nondiagonal terms, describing coupled motions, of the matrices are often omitted. Currently, no study validates the flexibility approach for predicting vertebral motions; nor have the effects of matrix approximations and simplifications been quantified. Therefore, the purpose of this study is to quantify flexibility relationships for cervical motion segments, examine the importance of nonlinear components of the flexibility matrix, and determine the extent to which multivariable relationships may alter motion prediction. To that end, using unembalmed human cervical spine motion segments, a full battery of flexibility tests were performed for a neutral orientation and also following an axial pretorque. Primary and coupled matrix components were described using linear and piecewise nonlinear incremental constants. A third matrix approach utilized multivariable incremental relationships. Measured motions were predicted using structural flexibility methods and evaluated using RMS error between predicted and measured responses. A full set of flexibility relationships describe primary and coupled motions for C3-C4 and C5-C6. A flexibility matrix using piecewise incremental responses offers improved predictions over one using linear methods (p<0.01). However, no significant improvement is obtained using nonlinear terms represented by a multivariable functional approach (p<0.2). Based on these findings, it is suggested that a multivariable approach for flexibility is more demanding experimentally and analytically while not offering improved motion prediction.     

An image-based method to automatically propagate bony landmarks: application to computational spine biomechanics

In attempts to elucidate the underlying mechanisms of spinal injuries and spinal deformities, several experimental and numerical studies have been conducted to understand the biomechanical behavior of the spine. However, numerical biomechanical studies suffer from uncertainties associated with hard- and soft-tissue anatomies. Currently, these parameters are identified manually on each mesh model prior to simulations. The determination of soft connective tissues on finite element meshes can be a tedious procedure, which limits the number of models used in the numerical studies to a few instances. In order to address these limitations, an image-based method for automatic morphing of soft connective tissues has been proposed. Results showed that the proposed method is capable to accurately determine the spatial locations of predetermined bony landmarks. The present method can be used to automatically generate patient-specific models, which may be helpful in designing studies involving a large number of instances and to understand the mechanical behavior of biomechanical structures across a given population.     

Comparison of Cervical Spine Anatomy in Calves,Pigs and Humans

Background Context

Animals are commonly used to model the human spine for in vitro and in vivo experiments. Many studies have investigated similarities and differences between animals and humans in the lumbar and thoracic vertebrae. However, a quantitative anatomic comparison of calf, pig, and human cervical spines has not been reported.

Purpose

To compare fundamental structural similarities and differences in vertebral bodies from the cervical spines of commonly used experimental animal models and humans.

Study Design

Anatomical morphometric analysis was performed on cervical vertebra specimens harvested from humans and two common large animals (i.e., calves and pigs).

Methods

Multiple morphometric parameters were directly measured from cervical spine specimens of twelve pigs, twelve calves and twelve human adult cadavers. The following anatomical parameters were measured: vertebral body width (VBW), vertebral body depth (VBD), vertebral body height (VBH), spinal canal width (SCW), spinal canal depth (SCD), pedicle width (PW), pedicle depth (PD), pedicle inclination (PI), dens width (DW), dens depth (DD), total vertebral width (TVW), and total vertebral depth (TVD).

Results

The atlantoaxial (C1–2) joint in pigs is similar to that in humans and could serve as a human substitute. The pig cervical spine is highly similar to the human cervical spine, except for two large transverse processes in the anterior regions ofC4–C6. The width and depth of the calf odontoid process were larger than those in humans. VBW and VBD of calf cervical vertebrae were larger than those in humans, but the spinal canal was smaller. Calf C7 was relatively similar to human C7, thus, it may be a good substitute.

Conclusion

Pig cervical vertebrae were more suitable human substitutions than calf cervical vertebrae, especially with respect to C1, C2, and C7. The biomechanical properties of nerve vascular anatomy and various segment functions in pig and calf cervical vertebrae must be considered when selecting an animal model for research on the spine.     

Inertial properties and loading rates affect buckling modes and injury mechanisms in the cervical spine

Cervical spine injuries continue to be a costly societal problem. Future advancements in injury prevention depend on improved physical and computational models which, in turn, are predicated on a better understanding of the responses of the neck during dynamic loading. Previous studies have shown that the tolerance of the neck is dependent on its initial position and its buckling behavior. This study uses a computational model to examine the mechanical factors influencing buckling behavior during impact to the neck. It was hypothesized that the inertial properties of the cervical spine influence the dynamics during compressive axial loading. The hypothesis was tested by performing parametric analyses of vertebral mass, mass moments of inertia, motion segment stiffness, and loading rate. Increases in vertebral mass resulted in increasingly complex kinematics and larger peak loads and impulses. Similar results were observed for increases in stiffness. Faster loading rates were associated with higher peak loads and higher-order buckling modes. The results demonstrate that mass has a great deal of influence on the buckling behavior of the neck, particularly with respect to the expression of higher-order modes. Injury types and mechanisms may be substantially altered by loading rate because inertial effects may influence whether the cervical spine fails in a compressive mode, or a bending mode.     

Biomechanical investigation of post-operative C5 palsy due to ossification of the posterior longitudinal ligament in different types of cervical spinal alignment

Post-operative C5 palsies are among the most common complications seen after cervical surgery for ossification of the posterior longitudinal ligament (OPLL). Although C5 palsy is a well-known complication of cervical spine surgery, its pathogenesis is poorly understood and depends on many other factors. In this study, a finite element model of the cervical spine and spinal cord-nerve roots complex structures was developed. The changes in stress in the cord and nerve roots, posterior shift of the spinal cord, and displacement and elongation of the nerve roots after laminectomy for cervical OPLL were analyzed for three different cervical sagittal alignments (lordosis, straight, and kyphosis). The results suggest that high stress concentrated on the nerve roots after laminectomy could be the main cause of C5 palsy because ossification of ligaments increases spinal cord shifting and root displacement. The type of sagittal alignment had no influence on changes in cord stress after laminectomy, although cases of kyphosis with a high degree of occupying ratio resulted in greater increases in nerve root stress after laminectomy. Therefore, kyphosis with a high OPLL occupying ratio could be a risk factor for poor surgical outcomes or post-operative complications and should be carefully considered for surgical treatment.     

A three-dimensional model of the human cervical spine for impact simulation

A three-dimensional analytical model of the cervical spine is described. The cervical vertebrae and the head are modeled as rigid bodies which are interconnected by deformable elements representing the intervertebral disks, facet joints, ligaments and muscles. A special pentahedral continuum element for representing the articular facets is described which effectively maintains stability of the cervical spine in both lateral and frontal plane accelerations, which is very difficult with multi-spring models of the facets. A simplified representation is used for the spine and body below the level of T1. The neck musculature is modeled by over 100 muscle elements representing 22 major muscle groups in the neck. The model has been validated for frontal and sideways impact accelerations by simulating published experimental data. Results are also presented to show the effects of the stretch reflex response on the dynamics of the head and neck under moderate acceleration.     

Animation of in vitro biomechanical tests.

Interdisciplinary communication of three-dimensional kinematic data arising from in vitro biomechanical tests is challenging. Complex kinematic representations such as the helical axes of motion (HAM) add to the challenge. The difficulty increases further when other quantities (i.e. load or tissue strain data) are combined with the kinematic data. The objectives of this study were to develop a method to graphically replay and animate in vitro biomechanical tests including HAM data. This will allow intuitive interpretation of kinematic and other data independent of the viewer's area of expertise. The value of this method was verified with a biomechanical test investigating load-sharing of the cervical spine. Three 3.0 mm aluminium spheres were glued to each of the two vertebrae from a C2-3 segment of a human cervical spine. Before the biomechanical tests, CT scans were made of the specimen (slice thickness=1.0 mm and slice spacing=1.5 mm). The specimens were subjected to right axial torsion moments (2.0 Nm). Strain rosettes mounted to the anterior surface of the C3 vertebral body and bilaterally beneath the facet joints on C3 were used to estimate the force flow through the specimen. The locations of the aluminium spheres were digitised using a space pointer and the motion analysis system. Kinematics were measured using an optoelectronic motion analysis system. HAMs were calculated to describe the specimen kinematics. The digitised aluminium sphere locations were used to match the CT and biomechanical test data (RMS errors between the CT and experimental points were less than 1.0 mm). The biomechanical tests were "replayed" by animating reconstructed CT models in accordance with the recorded experimental kinematics, using custom software. The animated test replays allowed intuitive analysis of the kinematic data in relation to the strain data. This technique improves the ability of experts from disparate backgrounds to interpret and discuss this type of biomechanical data.     

Moment arms of the human neck muscles in flexion, bending and rotation

There is a paucity of data available for the moment arms of the muscles of the human neck. The objective of the present study was to measure the moment arms of the major cervical spine muscles in vitro. Experiments were performed on five fresh-frozen human head-neck specimens using a custom-designed robotic spine testing apparatus. The testing apparatus replicated flexion-extension, lateral bending and axial rotation of each individual intervertebral joint in the cervical spine while all other joints were kept immobile. The tendon excursion method was used to measure the moment arms of 30 muscle sub-regions involving 13 major muscles of the neck about all three axes of rotation of each joint for the neutral position of the cervical spine. Significant differences in the moment arm were observed across sub-regions of individual muscles and across the intervertebral joints spanned by each muscle (p<0.05). Overall, muscle moment arms were larger in flexion-extension and lateral bending than in axial rotation, and most muscles had prominent moment arms in at least 2 out of the 3 joint motions investigated. This study emphasizes the importance of detailed representation of a muscle's architecture in prediction of its torque capacity about the individual joints of the cervical spine. The dataset produced may be useful in developing and validating computational models of the human neck.     

Finite element model of the human neck during omni-directional impacts. Part II: relation between cervical curvature and risk of injury

A detailed 3D FE model of the human neck was used to assess a possible relationship between risk of injury and cervical spine curvature for various impacts. A FE model was previously developed, representing the head and neck of a 50th percentile human with a normal lordotic curvature. The model behaviour was omni-directionally validated for various impacts using published results. For the present study, the model was deformed in order to obtain a straight and a kyphotic curvature, and for each geometry, rear-end, frontal, lateral and oblique impact were simulated. Although results showed similar kinematic patterns, significant differences were found in the distribution and peak values of ligament elongations, forces and moments along the cervical spine for the three configurations. It was concluded that the variability observed on the curvature of the human cervical spine may have a significant influence both on the behaviour and on the risk of injury of the neck during impact.