Abstract
Clinical investigations show that the cervical spine presents wide inter-individual variability, where its motion patterns and load sharing strongly depend on the anatomy. The magnitude and scope of cervical diseases, including disc degeneration, stenosis, and spondylolisthesis, constitute serious health and socioeconomic challenges that continue to increase along with the world”s growing aging population. Although complex exact finite element (FE) modeling is feasible and reliable for biomechanical studies, its clinical application has been limited as it is time-consuming and constrained to the input geometry, typically based on one or few subjects. The objective of this study was twofold: first to develop a validated parametric subject-specific FE model that automatically updates the geometry of the lower cervical spine based on different individuals; and second to investigate the motion patterns and biomechanics associated with typical cervical spine diseases. Six healthy volunteers participated in this study upon informed consent. 26 parameters were identified and measured for each vertebra in the lower cervical spine from Lateral and AP radiographs in neutral, flexion and extension viewpoints in the standing position. The lower cervical FE model was developed including the typical vertebrae (C3-C7), intervertebral discs, facet joints, and ligaments using ANSYS (PA, USA). In order to validate the FE model, the bottom surface of C7 was fixed, and a 73.6N preload together with a 1.8 N.m pure moment were input into the model in both flexion and extension. The results were compared to experimental studies from literature. Disc degeneration disease (DDD) was used as an example, where the geometry of C5-C6 disc was changed in the model to simulate 3 different grades of disc degeneration (mimicking grades 1 to 3), and the resulting biomechanical responses were evaluated. The average ranges of motion (ROM) were found to be 4.84 (±0.73) degrees and 5.36 (±0.68) degrees for flexion and extension for C5-C6 functional unit, respectively, in alignment with literature. The total ROM of the model with disc generation grades 2 and 3 was found to have decreased significantly as compared to the intact model. In contrast, the axial stresses on the degenerated discs were significantly higher than the intact discs for all 3 degeneration grades. Our preliminary results show that this novel validated subject-specific FE model provides a potential valuable tool for noninvasive time and cost effective analyses of cervical spine biomechanical (kinematic and kinetic) changes associated with various diseases. The model also provides an opportunity for clinicians to use quantitative data towards subject-specific informed therapy and surgical planning. Ongoing and future work includes expanding the studied population to investigate individuals with different cervical spine afflictions.
