An Elastic Finite Element Analysis of the C2-C7 Human Cervical Spine Segment

can a result and of the cord. Traumatic accidents such as motor vehicle, sports, falls, overloading, etc., are the most frequent causes of cervical injuries. No traumatic cervical injury causes are arthritis, cancer, compression fractures from osteoporosis, and inflammation of the spinal cord [2]. Therefore, it is necessary to fully understand the biomechanical response of the cervical spine to external stimuli (injury and dysfunction). Biomechanical finite element models facilitate the understanding of the mechanisms of dysfunction, spinal injury, and spinal responses to various clinical problems and treatments. In earlier cervical spine models developed to predict only spinal motion or displacement response to loading, the vertebrae are simulated as simple rigid bodies, while other connective tissues are represented as beams, truss or spring elements. Actually, anatomically precise 3D-FE models are being widely adopted to evaluate stresses, strains, displacements and rotations in complex vertebrae, facet joints, ligaments, and soft tissues. the provide on ABSTRACT A three-dimensional linear elastic FE analysis of six adjacent cervical vertebrae (C2-C7) has been developed. The model is composed of vertebrae and intervertebral discs. The objective of this research is to address the new advances in the generation of finite element meshes of the human cervical spine. A procedure to create a finite element model based on available biomedical closed volumes obtained from CT scans was introduced. The 3D biomedical parts were imported and meshed using the ABAQUS FE code. Continuous materials properties were assigned to the C2-C7 structure separating the vertebrae and discs in two different material groups. Encastered boundary conditions were used to fix translation and rotation in every node located on the inferior surface of the C7 vertebra. To assess the overall behavior of the C2-C7 spine segment in four different load modes (compression, lateral bending, flexion and extension), the model was statically loaded in four different surfaces created in the upward-facing long bony odontoid process (or dens) from the C2 vertebra. It is observed that higher strain values develop in the flexion load mode, and higher stress values in the extension load mode. It appears that the stress concentrators (created in irregular shapes and related coarse meshes) and the type of load mode are responsible for the high stresses observed. Also, it was found that at a constant pressure distribution, the spine response to the compression load mode is much lower in two orders of magnitude than the spine response to the flexion, bending and extension load modes. The study was carried out to illustrate the goodness of the finite element method on spine biomechanics research.

A three-dimensional linear elastic FE analysis of six adjacent cervical vertebrae (C2-C7) has been developed. The model is composed of vertebrae and intervertebral discs. The objective of this research is to address the new advances in the generation of finite element meshes of the human cervical spine. A procedure to create a finite element model based on available biomedical closed volumes obtained from CT scans was introduced. The 3D biomedical parts were imported and meshed using the ABAQUS FE code. Continuous materials properties were assigned to the C2-C7 structure separating the vertebrae and discs in two different material groups. Encastered boundary conditions were used to fix translation and rotation in every node located on the inferior surface of the C7 vertebra. To assess the overall behavior of the C2-C7 spine segment in four different load modes (compression, lateral bending, flexion and extension), the model was statically loaded in four different surfaces created in the upward-facing long bony odontoid process (or dens) from the C2 vertebra. It is observed that higher strain values develop in the flexion load mode, and higher stress values in the extension load mode. It appears that the stress concentrators (created in irregular shapes and related coarse meshes) and the type of load mode are responsible for the high stresses observed. Also, it was found that at a constant pressure distribution, the spine response to the compression load mode is much lower in two orders of magnitude than the spine response to the flexion, bending and extension load modes. The study was carried out to illustrate the goodness of the finite element method on spine biomechanics research. 1a. Each of these seven vertebrae in the neck has a predetermined function, creating very specific problems when they are not aligned.
Located at the top of the neck and base of the skull, the C1 vertebra is the first cervical vertebrae in the entire spine. Called atlas and shaped like a ring, this vertebra supports the head and allows it to rotate. It is held in place by muscles, and is secured by a thick and strong ligament called the transverse ligament. This bone is named the Holy Grail of the human spine, and misalignment at this joint can cause troubles for the entire body [3].
The first cervical vertebrae rotates around the dens (the superior protruding long bony element from the C2 vertebra), allowing the head to turn in different directions. The C2 vertebra, also called the Axis, is a hinge that makes the C1 vertebra pivot with the consequent turn of the head. When the C2 vertebra becomes subluxated (a misalignment in the axis vertebra), it can result in pain, difficulties to turn the head as well as creating specific secondary health issues related with sinuses, eyesight, hearing, head and even the tongue [4]. The joint Atlas (C1)-Axis (C2) is called the atlanto-axial joint. Unlike other vertebral joints, this special joint does not have an intervertebral disc. It is the most mobile section of the entire spine, and it is responsible for the 50% of the extension and flexion of the neck, and 50% of the rotation of the head [5]. In the present study, encastered boundary conditions were applied in the inferior surface of the C7 vertebra. Due to the length of its spinous process, this vertebra located at the base of the neck sometimes is called the prominent vertebra ( Figure 1). Numerous commonly open source and proprietary-software programs are available for use with DICOM data sets. For instance, to generate the FE closed volume geometry of the bone structure, segmentation mask during image processing of medical CT data is currently conducted by using specialist software such as 3D Slicer [6] or Mimics [7], and the specific meshing procedure with software such as OsiriX [8] or MeshLab [9]. 3D Slicer is an open source software platform for image informatics, medical image processing, and three-dimensional visualization. On the other hand, for processing and editing 3D triangular meshes, Mesh Lab is a powerful software that provides a set of tools for inspecting, editing, healing, cleaning, texturing, converting and rendering meshes. Features for preparing .STL models and for 3D printing processing raw data produced by 3D digitization tools/devices are also offered. STL is a file format native to the stereolithography CAD software widely used for rapid prototyping, 3D printing, and computer-aided manufacturing (CAM).

The Model
In this work, to construct the 3D-FE model, the cervical closed volume it was directly downloaded from the embodi3d web page [10]. The downloaded file was imported into the FE code ABAQUS [11] for creating an isotropic static finite element model.  The vertebrae and discs are generally modeled as a single linear elastic material. That is, the components of the vertebrae (cortical, cancellous, and posterior) and the components of the disc (nucleus and base of the annulus) are considered to form single materials but not composite materials. This assumption is clinically relevant to replicate disease or injury. As a first approach, in the present study, the entire vertebrae and discs were considered be isotropic single parts (single materials), so, they were not broken down into posterior elements, radial cortical shell, cancellous core, superior and inferior endplates, nucleus and annulus ground. During activities of daily living, it is estimated that the human cervical spine supports a compressive vertical load in the range of 120 DOI: 10.26717/BJSTR.2021.37.005937 29044 to 1200 N [12]. However, ex vivo experiments show it buckles at approximately 10 N. As differences between the estimated in vivo loads and the ex vivo load-carrying capacity exist, in this work, a pressure distribution of 10 MPa was statically applied in four different surfaces of the odontoid process, or dens (See Figure 3).
The application of a pressure load avoid the concentration of stress produced from point loads.  Figure 4. It can be seen that the highest strain values reported for the flexion and extension modes in Figure 3 are located in the intervertebral disc shown. Also, it is observed that at a constant pressure distribution, the spine response to the compression load mode is much lower in two orders of magnitude than the spine response to the flexion, bending and extension load modes.

Deformations in the compression mode do not result in strains that
would cause the material to exceed the proportional limit; that is, the material remains within the elastic limits of Hooke's law.
It appears that the stress concentrators (created in irregular shapes and related coarse meshes) and the type of load mode are responsible for the high stresses observed. Therefore, to test the sensibility of the FE mesh size in stress results, additional work has to be done by using fine mesh grids as those shown in Figure 5.    d.
It appears that the stress concentrators (created in irregular shapes and related coarse meshes) and the type of load mode are responsible for the high stresses observed.
e. It was found that at a constant pressure distribution, the spine response to the compression load mode is much lower in two orders of magnitude than the spine response to the flexion, bending and extension load modes. f.
The study was carried out to illustrate the goodness of the finite element method on spine biomechanics research.

Data Availability
The data used to support the findings of this study are included within the article.

Competing Interests
None declared.

Funding
None.

Employer
Universidad San Francisco de Quito, Ecuador.

Ethical Approval
Not required.