The purpose of this study was to investigate the relationship between cervical, thoracic, and lumbar spinal alignments in one automotive occupant seated posture. An image dataset of the spinal column in the automotive seated posture, previously acquired by an upright open magnetic resonance imaging (MRI) system, was re-analyzed in this study. Spinal alignments were presented by the geometrical centers of the vertebral bodies extracted from the image data. Cervical, thoracic, and lumbar spinal alignments were analyzed separately with multidimensional scaling (MDS). Based on distribution maps of cervical, thoracic, and lumbar spinal alignments created by MDS, representative spinal alignment patterns of the cervical, thoracic, and lumbar spines and the relationship between cervical, thoracic, and lumbar spinal alignments were investigated. As a result, this study found a correlation between cervical and thoracic spinal alignments in an automotive occupant seated posture. According to representative spinal alignment patterns illustrated by the distribution map of spinal alignments, subjects who had kyphotic cervical spinal alignment tended to have less kyphotic thoracic spinal alignment, while subjects who had lordotic cervical spinal alignment tended to have more kyphotic thoracic spinal alignment. For lumbar spinal alignments, no prominent relationship was found between cervical and thoracic spinal alignment in the seated condition of this study.

In studying traumatic brain injury (TBI), it has been long hypothesized that the head is more vulnerable to injury from impacts in certain directions or locations, as the relationship between impact force and the resulting neurological outcome is complex and can vary significantly between individual cases. Many studies have identified head angular acceleration to be the putative cause of brain trauma, but it is not well understood how impact location can affect the resulting head kinematics and tissue strain. Here, we identify the susceptibility of the head to accelerations and brain strain from normal forces at contact points across the surface of the skull and jaw using a three-dimensional, 20-degree-of-freedom rigid-body head and cervical spine model. We find that head angular acceleration and brain tissue strain resulting from an input force can vary by orders of magnitude based on impact location on the skull, with the mandible as the most vulnerable region. Conversely, head linear acceleration is not sensitive to contact location. Using these analyses, we present an optimization scheme to distribute helmet padding thickness to minimize angular acceleration, resulting in a reduction of angular acceleration by an estimated 25% at the most vulnerable contact point compared to uniform thickness padding. This work gives intuition behind the relationship between input force and resulting brain injury risk, and presents a framework for developing and evaluating novel head protection gear.

Repeated loading of ligamentous tissues during repetitive occupational and physical tasks even within physiological ranges of motion has been implicated in the development of pain and joint instability. The pathophysiological mechanisms of pain after repetitive joint loading are not understood. Within the cervical spine, excessive stretch of the facet joint and its capsular ligament has been implicated in the development of pain. Although a single facet joint distraction (FJD) at magnitudes simulating physiologic strains is insufficient to induce pain, it is unknown whether repeated stretching of the facet joint and ligament may produce pain. This study evaluated if repeated loading of the facet at physiologic nonpainful strains alters the capsular ligament's mechanical response and induces pain. Male rats underwent either two subthreshold facet joint distractions (STFJDs) or sham surgeries each separated by 2 days. Pain was measured before the procedure and for 7 days; capsular mechanics were measured during each distraction and under tension at tissue failure. Spinal glial activation was also assessed to probe potential pathophysiologic mechanisms responsible for pain. Capsular displacement significantly increased (p = 0.019) and capsular stiffness decreased (p = 0.008) during the second distraction compared to the first. Pain was also induced after the second distraction and was sustained at day 7 (p < 0.048). Repeated loading weakened the capsular ligament with lower vertebral displacement (p = 0.041) and peak force (p = 0.014) at tissue rupture. Spinal glial activation was also induced after repeated loading. Together, these mechanical, physiological, and neurological findings demonstrate that repeated loading of the facet joint even within physiologic ranges of motion can be sufficient to induce pain, spinal inflammation, and alter capsular mechanics similar to a more injurious loading exposure.

Many clinical case series have reported the predisposing factors for C5 palsy and have presented comparisons of the two types of laminoplasty. However, there have been no biomechanical studies focusing on cervical spinal cord and nerve root following laminoplasty. The purpose of this study is to investigate biomechanical changes in the spinal cord and nerve roots following the two most common types of laminoplasty, open-door and double-door laminoplasty, for cervical ossification of the posterior longitudinal ligament (OPLL). A finite element (FE) model of the cervical spine and spinal cord with nerve root complex structures was developed. Stress changes in the spinal cord and nerve roots, posterior shift of the spinal cord, and displacement of the cervical nerve roots were analyzed with two types of cervical laminoplasty models for variations in the degree of canal occupying ratio and shape of the OPLL. The shape and degree of spinal cord compression caused by the OPLL had more influence on the changes in stress, posterior shift of the spinal cord, and displacement of the nerve root than the type of laminoplasty. The lateral-type OPLL resulted in imbalanced stress on the nerve roots and the highest nerve root displacement. Type of laminoplasty and shape and degree of spinal cord compression caused by OPLL were found to influence the changes in stress and posterior displacement of the cervical spinal cord and nerve roots. Lateral-type OPLL might contribute to the development of C5 palsy due to the imbalanced stress and tension on the nerve roots after laminoplasty.

Advanced computational human body models (HBM) enabling enhanced safety require verification and validation at different levels or scales. Specifically, the motion segments, which are the building blocks of a detailed neck model, must be validated with representative experimental data to have confidence in segment and, ultimately, full neck model response. In this study, we introduce detailed finite element motion segment models and assess the models for quasi-static and dynamic loading scenarios. Finite element segment models at all levels in the lower human cervical spine were developed from scans of a 26-yr old male subject. Material properties were derived from the in vitro experimental data. The segment models were simulated in quasi-static loading in flexion, extension, lateral bending and axial rotation, and at dynamic rates in flexion and extension in comparison to previous experimental studies and new dynamic experimental data introduced in this study. Single-valued experimental data did not provide adequate information to assess the model biofidelity, while application of traditional corridor methods highlighted that data sets with higher variability could lead to an incorrect conclusion of improved model biofidelity. Data sets with continuous or multiple moment–rotation measurements enabled the use of cross-correlation for an objective assessment of the model and highlighted the importance of assessing all motion segments of the lower cervical spine to evaluate the model biofidelity. The presented new segment models of the lower cervical spine, assessed for range of motion and dynamic/traumatic loading scenarios, provide a foundation to construct a biofidelic model of the spine and neck, which can be used to understand and mitigate injury for improved human safety in the future.

A finite element model was used to compare the biomechanical properties of a novel anterior transpedicular screw artificial vertebral body system (AVBS) with a conventional anterior screw plate system (ASPS) for fixation in the lower cervical spine. A model of the intact cervical spine (C3–C7) was established. AVBS or ASPS constructs were implanted between C4 and C6. The models were loaded in three-dimensional (3D) motion. The Von Mises stress distribution in the internal fixators was evaluated, as well as the range of motion (ROM) and facet joint force. The models were generated and analyzed by mimics , geomagic studio , and ansys software. The intact model of the lower cervical spine consisted of 286,382 elements. The model was validated against previously reported cadaveric experimental data. In the ASPS model, stress was concentrated at the connection between the screw and plate and the connection between the titanium mesh and adjacent vertebral body. In the AVBS model, stress was evenly distributed. Compared to the intact cervical spine model, the ROM of the whole specimen after fixation with both constructs is decreased by approximately 3 deg. ROM of adjacent segments is increased by approximately 5 deg. Facet joint force of the ASPS and AVBS models was higher than those of the intact cervical spine model, especially in extension and lateral bending. AVBS fixation represents a novel reconstruction approach for the lower cervical spine. AVBS provides better stability and lower risk for internal fixator failure compared with traditional ASPS fixation.

While abnormal loading is widely believed to cause cervical spine disc diseases, in vivo cervical disc deformation during dynamic neck motion has not been well delineated. This study investigated the range of cervical disc deformation during an in vivo functional flexion–extension of the neck. Ten asymptomatic human subjects were tested using a combined dual fluoroscopic imaging system (DFIS) and magnetic resonance imaging (MRI)-based three-dimensional (3D) modeling technique. Overall disc deformation was determined using the changes of the space geometry between upper and lower endplates of each intervertebral segment (C3/4, C4/5, C5/6, and C6/7). Five points (anterior, center, posterior, left, and right) of each disc were analyzed to examine the disc deformation distributions. The data indicated that between the functional maximum flexion and extension of the neck, the anterior points of the discs experienced large changes of distraction/compression deformation and shear deformation. The higher level discs experienced higher ranges of disc deformation. No significant difference was found in deformation ranges at posterior points of all the discs. The data indicated that the range of disc deformation is disc level dependent and the anterior region experienced larger changes of deformation than the center and posterior regions, except for the C6/7 disc. The data obtained from this study could serve as baseline knowledge for the understanding of the cervical spine disc biomechanics and for investigation of the biomechanical etiology of disc diseases. These data could also provide insights for development of motion preservation surgeries for cervical spine.

Dynamic cervical spine loading can produce facet capsule injury. Despite a large proportion of neck pain being attributable to the C2/C3 facet capsule, potential mechanisms are not understood. This study replicated low-speed frontal and rear-end traffic collisions in occiput-C3 human cadaveric cervical spine specimens and used kinematic and full-field strain analyses to assess injury. Specimens were loaded quasi-statically in flexion and extension before and after dynamic rotation of C3 at 100 deg/s. Global kinematics in the sagittal plane were tracked at 1 kHz, and C2/C3 facet capsule full-field strains were measured. Dynamic loading did not alter the kinematics from those during quasi-static (QS) loading, but maximum principal strain (MPS) and shear strain (SS) were significantly higher (p = 0.028) in dynamic flexion than for the same quasi-static conditions. The full-field strain analysis demonstrated that capsule strain was inhomogeneous, and that the peak MPS generally occurred in the anterior aspect and along the line of the C2/C3 facet joint. The strain magnitude in dynamic flexion continued to rise after the rotation of C3 had stopped, with a peak MPS of 12.52 ± 4.59% and a maximum SS of 5.34 ± 1.60%. The peak MPS in loading representative of rear-end collisions approached magnitudes previously shown to induce pain in vivo, whereas strain analysis using linear approaches across the facet joint was lower and may underestimate injury risk compared to full-field analysis. The time at which peak MPS occurred suggests that the deceleration following a collision is critical in relation to the production of injurious strains within the facet capsule.

Clinically, spinal cord injuries (SCIs) are radiographically evaluated and diagnosed from plain radiographs, computed tomography (CT), and magnetic resonance imaging. However, it is difficult to conclude that radiographic evaluation of SCI can directly explain the fundamental mechanism of spinal cord damage. The von-Mises stress and maximum principal strain are directly associated with neurological damage in the spinal cord from a biomechanical viewpoint. In this study, the von-Mises stress and maximum principal strain in the spinal cord as well as the cord cross-sectional area (CSA) were analyzed under various magnitudes for contusion, dislocation, and distraction SCI mechanisms, using a finite-element (FE) model of the cervical spine with spinal cord including white matter, gray matter, dura mater with nerve roots, and cerebrospinal fluid (CSF). A regression analysis was performed to find correlation between peak von-Mises stress/peak maximum principal strain at the cross section of the highest reduction in CSA and corresponding reduction in CSA of the cord. Dislocation and contusion showed greater peak stress and strain values in the cord than distraction. The substantial increases in von-Mises stress as well as CSA reduction similar to or more than 30% were produced at a 60% contusion and a 60% dislocation, while the maximum principal strain was gradually increased as injury severity elevated. In addition, the CSA reduction had a strong correlation with peak von-Mises stress/peak maximum principal strain for the three injury mechanisms, which might be fundamental information in elucidating the relationship between radiographic and mechanical parameters related to SCI.

Mathematical cervical spine models allow for studying of impact loading that can cause whiplash associated disorders (WAD). However, existing models only cover the male anthropometry, despite the female population being at a higher risk of sustaining WAD in automotive rear-end impacts. The aim of this study is to develop and validate a ligamentous cervical spine intended for biomechanical research on the effect of automotive impacts. A female model has the potential to aid the design of better protection systems as well as improve understanding of injury mechanisms causing WAD. A finite element (FE) mesh was created from surface data of the cervical vertebrae of a 26-year old female (stature 167 cm, weight 59 kg). Soft tissues were generated from the skeletal geometry and anatomical literature descriptions. Ligaments were modeled with nonlinear elastic orthotropic membrane elements, intervertebral disks as composites of nonlinear elastic bulk elements, and orthotropic anulus fibrosus fiber layers, while cortical and trabecular bones were modeled as isotropic plastic–elastic. The model has geometrical features representative of the female cervical spine—the largest average difference compared with published anthropometric female data was the vertebral body depth being 3.4% shorter for the model. The majority the cervical segments compare well with respect to biomechanical data at physiological loads, with the best match for flexion–extension loads and less biofidelity for axial rotation. An average female FE ligamentous cervical spine model was developed and validated with respect to physiological loading. In flexion–extension simulations with the developed female model and an existing average male cervical spine model, a greater range of motion (ROM) was found in the female model.

Numerical finite element (FE) models of the neck have been developed to simulate occupant response and predict injury during motor vehicle collisions. However, there is a paucity of data on the response of young cervical spine segments under dynamic loading in flexion and extension, which is essential for the development or validation of tissue-level FE models. This limitation was identified during the development and validation of the FE model used in this study. The purpose of this study was to measure the high rotation rate loading response of human cervical spine segments in flexion and extension, and to investigate a new tissue-level FE model of the cervical spine with the experimental data to address a limitation in available data. Four test samples at each segment level from C2–C3 to C7–T1 were dissected from eight donors and were tested to 10 deg of rotation at 1 and 500 deg/s in flexion and extension using a custom built test apparatus. There was strong evidence (p < 0.05) of increased stiffness at the higher rotation rate above 4 deg of rotation in flexion and at 8 deg and 10 deg of rotation in extension. Cross-correlation software, Cora, was used to evaluate the fit between the experimental data and model predictions. The average rating was 0.771, which is considered to demonstrate a good correlation to the experimental data.

A majority of the middle-aged population exhibit cervical spondylosis that may require decompression and fusion of the affected level. Minimally invasive cervical fusion is an attractive option for decreasing operative time, morbidity, and mortality rates. A novel interfacet joint spacer (DTRAX facet screw system, Providence Medical) promises minimally invasive deployment resulting in decompression of the neuroforamen and interfacet fusion. The present study investigates the effectiveness of the device in minimizing intervertebral motion to promote fusion, decompression of the nerve root during bending activity, and performance of the implant to adhere to anatomy during repeated bending loads. We observed flexion, extension, lateral bending, and axial rotation resonant overshoot mode (ROM) in cadaver models of c-spine treated with the interfacet joint spacer (FJ spacer) as stand-alone and supplementing anterior plating. The FJ spacer was deployed bilaterally at single levels. Specimens were placed at the limit of ROM in flexion, extension, axial bending, and lateral bending. 3D images of the foramen were taken and postprocessed to quantify changes in foraminal area. Stand-alone spacer specimens were subjected to 30,000 cycles at 2 Hz of nonsimultaneous flexion-extension and lateral bending under compressive load and X-ray imaged at regular cycle intervals for quantitative measurements of device loosening. The stand-alone FJ spacer increased specimen stiffness in all directions except extension. 86% of all deployments resulted in some level of foraminal distraction. The rate of effective distraction was maintained in flexed, extended, and axially rotated postures. Two specimens demonstrated no detectable implant loosening (<0.25 mm). Three showed unilateral subclinical loosening (0.4 mm maximum), and one had subclinical loosening bilaterally (0.5 mm maximum). Results of our study are comparable to previous investigations into the stiffness of other stand-alone minimally invasive technologies. The FJ spacer system effectively increased stiffness of the affected level comparable to predicate systems. Results of this study indicate the FJ spacer increases foraminal area in the cervical spine, and decompression is maintained during bending activities. Clinical studies will be necessary to determine whether the magnitude of decompression observed in this cadaveric study will effectively treat cervical radiculopathy; however, results of this study, taken in context of successful decompression treatments in the lumbar spine, are promising for the continued development of this product. Results of this biomechanical study are encouraging for the continued investigation of this device in animal and clinical trials, as they suggest the device is well fixated and mechanically competent.

Cervical spine ligaments have an important role in providing spinal cord stability and restricting excessive movements. Therefore, it is of great importance to study the mechanical properties and model the response of these ligaments. The aim of this study is to characterize the aging effects on the failure properties and model the damage of three cervical spine ligaments: the anterior and the posterior longitudinal ligament and the ligamentum flavum. A total of 46 samples of human cadaveric ligaments removed within 24–48 h after death have been tested. Uniaxial tension tests along the fiber direction were performed in physiological conditions. The results showed that aging decreased the failure properties of all three ligaments (failure load, failure elongation). Furthermore, the reported nonlinear response of cervical ligaments has been modeled with a combination of the previously reported hyperelastic and damage model. The model predicted a nonlinear response and damage region. The model fittings are in agreement with the experimental data and the quality of agreement is represented with the values of the coefficient of determination close to 1.

Current understanding of the biomechanics of cervical spine injuries in head-first impact is based on decades of epidemiology, mathematical models, and in vitro experimental studies. Recent mathematical modeling suggests that muscle activation and muscle forces influence injury risk and mechanics in head-first impact. It is also known that muscle forces are central to the overall physiologic stability of the cervical spine. Despite this knowledge, the vast majority of in vitro head-first impact models do not incorporate musculature. We hypothesize that the simulation of the stabilizing mechanisms of musculature during head-first osteoligamentous cervical spine experiments will influence the resulting kinematics and injury mechanisms. Therefore, the objective of this study was to document differences in the kinematics, kinetics, and injuries of ex vivo osteoligamentous human cervical spine and surrogate head complexes that were instrumented with simulated musculature relative to specimens that were not instrumented with musculature. We simulated a head-first impact (3 m/s impact speed) using cervical spines and surrogate head specimens (n = 12). Six spines were instrumented with a follower load to simulate in vivo compressive muscle forces, while six were not. The principal finding was that the axial coupling of the cervical column between the head and the base of the cervical spine (T1) was increased in specimens with follower load. Increased axial coupling was indicated by a significantly reduced time between head impact and peak neck reaction force (p = 0.004) (and time to injury (p = 0.009)) in complexes with follower load relative to complexes without follower load. Kinematic reconstruction of vertebral motions indicated that all specimens experienced hyperextension and the spectrum of injuries in all specimens were consistent with a primary hyperextension injury mechanism. These preliminary results suggest that simulating follower load that may be similar to in vivo muscle forces results in significantly different impact kinetics than in similar biomechanical tests where musculature is not simulated.

The effects of degeneration and surgery on cervical spine mechanics are commonly evaluated through in vitro testing and finite element models derived from these tests. The objectives of the current study were to estimate the load applied to the C2 vertebra during in vivo functional flexion-extension and to evaluate the effects of anterior cervical arthrodesis on spine kinetics. Spine and head kinematics from 16 subjects (six arthrodesis patients and ten asymptomatic controls) were determined during functional flexion-extension using dynamic stereo X-ray and conventional reflective markers. Subject-specific inverse dynamics models, including three flexor muscles and four extensor muscles attached to the skull, estimated the force applied to C2. Total force applied to C2 was not significantly different between arthrodesis and control groups at any 10 deg increment of head flexion-extension (all p values ≥ 0.937). Forces applied to C2 were smallest in the neutral position, increased slowly with flexion, and increased rapidly with extension. Muscle moment arms changed significantly during flexion-extension, and were dependent upon the direction of head motion. The results suggest that in vitro protocols and finite element models that apply constant loads to C2 do not accurately represent in vivo cervical spine kinetics.

Knowledge of the distributions of cervical-spine curvature is needed for computational studies of cervical-spine injury in motor-vehicle crashes. Many methods of specifying spinal curvature have been proposed, but they often involve qualitative assessment or a large number of parameters. The objective of this study was to develop a quantitative method of characterizing cervical-spine curvature using a small number of parameters. 180 sagittal X-rays of subjects seated in automotive posture with their necks in neutral, flexed, and extended postures were collected in the early 1970s. Subjects were selected to represent a range of statures and ages for each gender. X-rays were reanalyzed using advanced technology and statistical methods. Coordinates of the posterior margins of the vertebral bodies and dens were digitized. Bézier splines were fit through the coordinates of these points. The interior control points that define the spline curvature were parameterized as a vector angle and length. By defining the length as a function of the angle, cervical-spine curvature was defined with just two parameters: superior and inferior Bézier angles. A classification scheme was derived to sort each curvature by magnitude and type of curvature (lordosis versus S-shaped versus kyphosis; inferior or superior location). Cervical-spine curvature in an automotive seated posture varies with gender and age but not stature. Average values of superior and inferior Bézier angles for cervical spines in flexion, neutral, and extension automotive postures are presented for each gender and age group. Use of Bézier splines fit through posterior margins offers a quantitative method of characterizing cervical-spine curvature using two parameters: superior and inferior Bézier angles.

To understand the role seating plays in the support of posture and spinal articulation, it is necessary to study the interface between a human and the seat. However, a method to quantify lumbar curvature in commercially available unmodified seats does not currently exist. This work sought to determine if the lumbar curvature for normal ranges of seated posture could be documented by using body landmarks located on the anterior portion of the body. The development of such a methodology will allow researchers to evaluate spinal articulation of a seated subject while in standard, commercially available seats and chairs. Anterior measurements of boney landmarks were used to quantify the relative positions of the ribcage and pelvis while simultaneous posterior measurements were made of lumbar curvature. The relationship between the anterior and the posterior measures was compared. The predictive capacity of this approach was evaluated by determining linear and second-order regressions for each of the four postures across all subjects and conducting a leave-one-out cross validation. The relationships between the anterior and posterior measures were approximated by linear and second-order polynomial regressions (r 2 = 0.829, 0.935 respectively) across all postures. The quantitative analysis showed that openness had a significant relationship with lumbar curvature, and a first-order regression was superior to a second-order regression. Average standard errors in the prediction were 5.9° for the maximum kyphotic posture, 9.9° for the comfortable posture, 12.8° for the straight and tall, and 22.2° for the maximum lordotic posture. These results show predictions of lumbar curvature are possible in seated postures by using a motion capture system and anterior measures. This method of lumbar curvature prediction shows potential for use in the assessment of seated spinal curvatures and the corresponding design of seating to accommodate those curvatures; however, additional inputs will be necessary to better predict the postures as lordosis is increased.

The facet joint contributes to the normal biomechanical function of the spine by transmitting loads and limiting motions via articular contact. However, little is known about the contact pressure response for this joint. Such information can provide a quantitative measure of the facet joint’s local environment. The objective of this study was to measure facet pressure during physiologic bending in the cervical spine, using a joint capsule-sparing technique. Flexion and extension bending moments were applied to six human cadaveric cervical spines. Global motions (C2-T1) were defined using infra-red cameras to track markers on each vertebra. Contact pressure in the C5-C6 facet was also measured using a tip-mounted pressure transducer inserted into the joint space through a hole in the postero-inferior region of the C5 lateral mass. Facet contact pressure increased by 67.6 ± 26.9 kPa under a 2.4 Nm extension moment and decreased by 10.3 ± 9.7 kPa under a 2.7 Nm flexion moment. The mean rotation of the overall cervical specimen motion segments was 9.6 ± 0.8° and was 1.6 ± 0.7° for the C5-C6 joint, respectively, for extension. The change in pressure during extension was linearly related to both the change in moment (51.4 ± 42.6 kPa/Nm) and the change in C5-C6 angle (18.0 ± 108.9 kPa/deg). Contact pressure in the inferior region of the cervical facet joint increases during extension as the articular surfaces come in contact, and decreases in flexion as the joint opens, similar to reports in the lumbar spine despite the difference in facet orientation in those spinal regions. Joint contact pressure is linearly related to both sagittal moment and spinal rotation. Cartilage degeneration and the presence of meniscoids may account for the variation in the pressure profiles measured during physiologic sagittal bending. This study shows that cervical facet contact pressure can be directly measured with minimal disruption to the joint and is the first to provide local pressure values for the cervical joint in a cadaveric model.

With motion-sparing disk replacement implants gaining popularity as an alternative to anterior cervical discectomy and fusion (ACDF) for the treatment of certain spinal degenerative disorders, recent laboratory investigations have studied the effects of disk replacement and implant design on spinal kinematics and kinetics. Particularly relevant to cervical disk replacement implant design are any postoperative changes in solid stresses or contact conditions in the articular cartilage of the posterior facets, which are hypothesized to lead to adjacent-level degeneration. Such changes are commonly investigated using finite element methods, but significant simplification of the articular geometry is generally employed. The impact of such geometric representations has not been thoroughly investigated. In order to assess the effects of different models of cartilage geometry on load transfer and contact pressures in the lower cervical spine, a finite element model was generated using cadaver-based computed tomography imagery. Mesh resolution was varied in order to establish model convergence, and cadaveric testing was undertaken to validate model predictions. The validated model was altered to include four different geometric representations of the articular cartilage. Model predictions indicate that the two most common representations of articular cartilage geometry result in significant reductions in the predictive accuracy of the models. The two anatomically based geometric models exhibited less computational artifact, and relatively minor differences between them indicate that contact condition predictions of spatially varying thickness models are robust to anatomic variations in cartilage thickness and articular curvature. The results of this work indicate that finite element modeling efforts in the lower cervical spine should include anatomically based and spatially varying articular cartilage thickness models. Failure to do so may result in loss of fidelity of model predictions relevant to investigations of physiological import.

While the individual tissues of the vertebral joint demonstrate viscoelastic properties, the global viscoelastic properties of the lumbar vertebral joint are not well established. This study investigated how changes in displacement rate influenced the mechanical response of the porcine cervical spine (a surrogate or model for the human lumbar spine) exposed to acute anterior shear failure loading. Thirty porcine cervical spine specimens (15 C3-C4 and 15 C5-C6) were placed under a 1600 N compressive load and subsequently loaded in anterior shear to failure at one of three randomly assigned displacement rates (1 mm/s, 4 mm/s, or 16 mm/s). Ultimate anterior shear force, ultimate displacement, average stiffness, and energy stored until failure were calculated. Load rate in the elastic region was also calculated to compare the load rates used in this study to those used in previous studies. Changes in displacement rate affected the C3-C4 and C5-C6 specimens differently. C5-C6 specimens tested at 16 mm/s had an ultimate force that was 28% and 23% higher than at 1 ( p = 0.0215 ) and 4 mm/s ( p = 0.0461 ) , respectively. The average stiffness to failure of the C5-C6 specimens tested at 16 mm/s was 52% higher than at 4 mm/s ( p = 0.0289 ) . No such differences were found for the C3-C4 specimens. An increase in the anterior shear displacement rate did not necessarily demonstrate viscoelasticity of the vertebral joint. Specimen intervertebral levels were affected differently by changes in anterior shear displacement rate, which may have been a result of anatomical and postural differences between the two levels. Future studies should further investigate the effect of displacement rate on the spine and the inconsistencies between different specimen levels.