Estimating Facet Joint Apposition with Specimen-Specific Computer Models of Subaxial Cervical Spine Kinematics

Computational models of experimental data can provide a noninvasive method to estimate spinal facet joint biomechanics. Existing models typically consider each vertebra as one rigid-body and assume uniform facet cartilage thickness. However, facet deflection occurs during motion, and cervical facet cartilage is nonuniform. Multi rigid-body computational models were used to investigate the effect of specimen-specific cartilage profiles on facet contact area estimates. Twelve C6/C7 segments underwent non-destructive intervertebral motions. Kinematics and facet deflections were measured. Three-dimensional models of the vertebra and cartilage thickness estimates were obtained from pre-test CT data. Motion-capture data was applied to two model types (2RB: C6, C7 vertebrae each one rigid body; 3RB: left and right C6 posterior elements, and C7 vertebrae, each one rigid body) and maximum facet mesh penetration was compared. Constant thickness cartilage (CTC) and spatially-varying thickness cartilage (SVTC) profiles were applied to the facet surfaces of the 3RB model. Cartilage apposition area (CAA) was compared. Linear mixed-effects models were used for all quantitative comparisons. The 3RB model significantly reduced penetrating mesh elements by accounting for facet deflections (p = 0.001). The CTC profile resulted in incongruent facet articulation, whereas realistic congruence was observed for the SVTC profile. The SVTC profile demonstrated significantly larger CAA than the CTC model (p < 0.001).