Srirangam Kumaresan

Biomechanics of the cervical spine Part 2. Cervical spine soft tissue responses and biomechanical modeling.

OBJECTIVE The responses and contributions of the soft tissue structures of the human neck are described with a focus on mathematical modeling. Spinal ligaments, intervertebral discs, zygapophysial joints, and uncovertebral joints of the cervical spine are included. Finite element modeling approaches have been emphasized. Representative data relevant to the development and execution of the model are discussed. A brief description is given on the functional mechanical role of the soft tissue components. Geometrical characteristics such as length and cross-sectional areas, and material properties such as force-displacement and stress-strain responses, are described for all components. Modeling approaches are discussed for each soft tissue structure. The final discussion emphasizes the normal and abnormal (e.g., degenerative joint disease, iatrogenic alteration, trauma) behaviors of the cervical spine with a focus on all these soft tissue responses. A brief description is provided on the modeling of the developmental biomechanics of the pediatric spine with a focus on soft tissues. Relevance. Experimentally validated models based on accurate geometry, material property, boundary, and loading conditions are useful to delineate the clinical biomechanics of the spine. Both external and internal responses of the various spinal components, a data set not obtainable directly from experiments, can be determined using computational models. Since soft tissues control the complex structural response, an accurate simulation of their anatomic, functional, and biomechanical characteristics is necessary to understand the behavior of the cervical spine under normal and abnormal conditions such as facetectomy, discectomy, laminectomy, and fusion.

Geometric and Mechanical Properties of Human Cervical Spine Ligaments

This study characterized the geometry and mechanical properties of the cervical ligaments from C2-T1 levels. The lengths and cross-sectional areas of the anterior longitudinal ligament, posterior longitudinal ligament, joint capsules, ligamentum flavum, and interspinous ligament were determined from eight human cadavers using cryomicrotomy images. The geometry was defined based on spinal anatomy and its potential use in complex mathematical models. The biomechanical force-deflection, stiffness, energy, stress, and strain data were obtained from 25 cadavers using in situ axial tensile tests. Data were grouped into middle (C2-C5) and lower (C5-T1) cervical levels. Both the geometric length and area of cross section, and the biomechanical properties including the stiffness, stress, strain, energy, and Youngs modulus, were presented for each of the five ligaments. In both groups, joint capsules and ligamentum flavum exhibited the highest cross-sectional area (p < 0.005), while the longitudinal ligaments had the highest length measurements. Although not reaching statistical significance, for all ligaments, cross-sectional areas were higher in the C5-T1 than in the C2-C5 group; and lengths were higher in the C2-C5 than in the C5-T1 group with the exception of the flavum (Table 1 in the main text). Force-deflection characteristics (plots) are provided for all ligaments in both groups. Failure strains were higher for the ligaments of the posterior (interspinous ligament, joint capsules, and ligamentum flavum) than the anterior complex (anterior and posterior longitudinal ligaments) in both groups. In contrast, the failure stress and Youngs modulus were higher for the anterior and posterior longitudinal ligaments compared to the ligaments of the posterior complex in the two groups. However, similar tendencies in the structural responses (stiffness, energy) were not found in both groups. Researchers attempting to incorporate these data into stress-analysis models can choose the specific parameter(s) based on the complexity of the model used to study the biomechanical behavior of the human cervical spine.

Contribution of disc degeneration to osteophyte formation in the cervical spine: a biomechanical investigation

Cervical spine disorders such as spondylotic radiculopathy and myelopathy are often related to osteophyte formation. Bone remodeling experimental—analytical studies have correlated biomechanical responses such as stress and strain energy density to the formation of bony outgrowth. Using these responses of the spinal components, the present study was conducted to investigate the basis for the occurrence of disc‐related pathological conditions. An anatomically accurate and validated intact finite element model of the C4‐C5‐C6 cervical spine was used to simulate progressive disc degeneration at the C5‐C6 level. Slight degeneration included an alteration of material properties of the nucleus pulposus representing the dehydration process. Moderate degeneration included an alteration of fiber content and material properties of the anulus fibrosus representing the disintegrated nature of the anulus in addition to dehydrated nucleus. Severe degeneration included decrease in the intervertebral disc height with dehydrated nucleus and disintegrated anulus. The intact and three degenerated models were exercised under compression, and the overall force—displacement response, local segmental stiffness, anulus fiber strain, disc bulge, anulus stress, load shared by the disc and facet joints, pressure in the disc, facet and uncovertebral joints, and strain energy density and stress in the vertebral cortex were determined. The overall stiffness (C4‐C6) increased with the severity of degeneration. The segmental stiffness at the degenerated level (C5‐C6) increased with the severity of degeneration. Intervertebral disc bulge and anulus stress and strain decreased at the degenerated level. The strain energy density and stress in vertebral cortex increased adjacent to the degenerated disc. Specifically, the anterior region of the cortex responded with a higher increase in these responses. The increased strain energy density and stress in the vertebral cortex over time may induce the remodeling process according to Wolffs law, leading to the formation of osteophytes.

Finite element analysis of the cervical spine: a material property sensitivity study

OBJECTIVE The study determined the effect of variations in the material properties of the cervical spinal components on the output of the finite element analysis (external and internal responses of the cervical spine) under physiologic load vectors. DESIGN A three-dimensional (3D) anatomically accurate finite element model comprising of the C4-C5-C6 cervical spine unit including the three vertebrae, two interconnecting intervertebral discs, and the anterior and posterior ligament complex is used. BACKGROUND The effect of material property variations of spinal components on the human lumbar spine biomechanics is extensively studied. However, a similar investigation of the cervical spine is lacking. METHODS Parametric studies on the variations in the material properties of all the cervical spine components including the cortical shell, cancellous core, endplates, intervertebral disc, posterior elements and ligaments were conducted by exercising the 3D finite element model under flexion, extension, lateral bending and axial torsion loading modes. Low, basic and high material property cases for each of the six components under all the four physiologic loading modes were considered in the finite element analysis. A total of 432 results were evaluated to analyze the external angular rotation, and the internal stresses in the middle vertebral body, the superior and inferior endplates and the two intervertebral discs. RESULTS Variations in the material properties of the different cervical spinal components produced dissimilar changes in the external and internal responses. Variations in the material properties of the cancellous core, cortical shell, endplates and posterior element structures representing the hard tissues did not affect the external angular motion, and the internal stresses of the inferior and superior intervertebral discs under all four loading modes. In contrast, variations in the material properties of the intervertebral disc and ligament structures representing the soft tissues significantly altered the angular motion, and the stresses in the inferior and superior intervertebral discs of the cervical spine. CONCLUSION The material properties of the soft tissue structures have a preponderant effect on the external and internal responses of the cervical spine compared with the changes in the material properties of the hard tissue structures. RELEVANCE Bone remodeling (e.g., osteophyte) secondary to degeneration of the human cervical joints may be explained by a change in the material property of the soft tissues, coupled with an increase in stress (due to these material property variations) in the spinal components. Consequently, to accurately predict the biomedical effects of cervical spine degeneration, it is critical to accurately determine the material property of these components.

Finite-element models of the human head

A review is presented of the existing finite-element (FE) models for the biomechanics of human head injury. Finite element analysis can be an important tool in describing the injury biomechanics of the human head. Complex geometric and material properties pose challenges to FE modelling. Various assumptions and simplifications are made in model development that require experimental validation. More recent models incorporate anatomic details with higher precision. The cervical vertebral column and spinal cord are included. Model results have been more qualitative than quantitative owing to the lack of adequate experimental validation. Advances include transient stress distribution in the brain tissue, frequency responses, effects of boundary conditions, pressure release mechanism of the foramen magnum and the spinal cord, verification of rotation and cavitation theories of brain injury, and protective effects of helmets. These theoretical results provide a basic understanding of the internal biomechanical responses of the head under various dynamic loading conditions. Basic experimental research is still needed to determine more accurate material properties and injury tolerance criteria, so that FE models can fully exercise their analytical and predictive power for the study and prevention of human head injury.

Finite Element Applications in Human Cervical Spine Modeling

The authors present a comprehensive state-of-the-art and critical review of the finite element models of the human cervical spine. They also focused on the developments in model construction (geometry generation), constitutive law (material property) identification, loading and boundary condition details, and validation, the most important phase. A data base of available experimental sources is also provided, which can be used by the modeler for validating the finite element model. The potential developments in finite element modeling of the human cervical spine are discussed.

Finite element modeling of the cervical spine: role of intervertebral disc under axial and eccentric loads.

An anatomically accurate, three-dimensional, nonlinear finite element model of the human cervical spine was developed using computed tomography images and cryomicrotome sections. The detailed model included the cortical bone, cancellous core, endplate, lamina, pedicle, transverse processes and spinous processes of the vertebrae; the annulus fibrosus and nucleus pulposus of the intervertebral discs; the uncovertebral joints; the articular cartilage, the synovial fluid and synovial membrane of the facet joints; and the anterior and posterior longitudinal ligaments, interspinous ligaments, capsular ligaments and ligamentum flavum. The finite element model was validated with experimental results: force-displacement and localized strain responses of the vertebral body and lateral masses under pure compression, and varying eccentric anterior-compression and posterior-compression loading modes. This experimentally validated finite element model was used to study the biomechanics of the cervical spine intervertebral disc by quantifying the internal axial and shear forces resisted by the ventral, middle, and dorsal regions of the disc under the above axial and eccentric loading modes. Results indicated that higher axial forces (compared to shear forces) were transmitted through different regions of the disc under all loading modes. While the ventral region of the disc resisted higher variations in axial force, the dorsal region transmitted higher shear forces under all loading modes. These findings may offer an insight to better understand the biomechanical role of the human cervical spine intervertebral disc.

Finite element modeling of the C4–C6 cervical spine unit

This study was conducted to develop a detailed, three-dimensional, anatomically accurate finite element model of the human cervical spine structure using close-up computed tomography scans and to validate against experimental data. The finite element model of the three vertebra segment C4-C6 unit consisted of 9178 solid elements and 1193 thin shell elements. The force-displacement response under axial compression correlated well with experimental data. Because of the inclusion of three levels in the spinal structure, it was possible to determine the internal mechanics of the various components at each level. The applicability of the model was illustrated by adopting appropriate material properties from literature. Results indicated that, the stresses in the anterior column were higher compared to the posterior column at the inferior level, while the opposite was found to be true at the superior level. The superior and inferior endplate stresses were higher in the middle vertebral body compared to the adjacent vertebrae. In addition, the stresses in the cancellous core of the middle, unconstrained vertebral body were higher. The present three-dimensional finite element model offers an additional facet to a better understanding of the biomechanics of the human cervical spine.

Finite element modeling approaches of human cervical spine facet joint capsule

The human cervical spine facet joint capsule was modeled using four nonlinear finite element approaches: slideline, contact surface, hyperelastic, and fluid models. Slideline elements and contact surface definitions were used in the first two models to simulate the synovial fluid between the articulating cartilages. Incompressible solid elements approximated the synovial fluid in the hyperelastic model. Hydrostatic fluid elements idealized the synovial fluid in the fluid model. The finite element analysis incorporated geometric, material and contact nonlinearities. All models were subjected to compression, flexion, extension, and lateral bending. The fluid model idealization better approximates the actual facet joint anatomy and its behavior than the gap assumption in the slideline and contact surface models, and the solid element simulation in the hyperelastic model.

Finite element model of the human lower cervical spine: parametric analysis of the C4-C6 unit.

In this study, a three-dimensional finite element model of the human lower cervical spine (C4-C6) was constructed. The mathematical model was based on close-up CT scans from a young human cadaver. Cortical shell, cancellous core, endplates, and posterior elements including the lateral masses, pedicle, lamina, and transverse and spinous processes, and the intervertebral disks, were simulated. Using the material properties from literature, the 10,371-element model was exercised under an axial compressive mode of loading. The finite element model response agreed with literature. As a logical step, a parametric study was conducted by evaluating the biomechanical response secondary to changes in the elastic moduli of the intervertebral disk and the endplates. In the stress analysis, the minimum principal compressive stress was used for the cancellous core of the vertebral body and von Mises stress was used for the endplate component. The model output indicated that an increase in the elastic modulii of the disk resulted in an increase in the endplate stresses at all the three spinal levels. In addition, the inferior endplate of the middle vertebral body responded with the highest mean compressive stress followed by its superior counterpart. Furthermore, the middle vertebral body produced the highest compressive stresses compared to its counterparts. These findings appear to correlate with experimental results as well as common clinical experience wherein cervical fractures are induced due to external compressive forces. As a first step, this model will lead to more advanced simulations as additional data become available.