Dawn M. Elliott

The intervertebral disc is an important structure for load transfer through the spine. Its injury and degeneration have been linked to pain and spinal fractures. Disc injury and spine fractures are associated with high stresses; however, these stresses cannot be measured, necessitating the use of finite element (FE) models. These models should include the disc's complex structure, as changes in disc geometry have been linked to altered mechanical behavior. However, image‐based models using disc‐specific structures have yet to be established. This study describes a multiphasic FE modeling approach for noninvasive estimates of subject‐specific intervertebral disc mechanical behavior based on medical imaging. The models ( n = 22) were used to study the influence of disc geometry on the predicted global mechanical response (moments and forces), internal local disc stresses, and tractions at the interface between the disc and the bone. Disc geometry was found to have a strong influence on the predicted moments and forces on the disc ( R 2 = 0.69–0.93), while assumptions regarding the side curvature (bulge) of the disc had only a minor effect. Strong variability in the predicted internal disc stresses and tractions was observed between the models (mean absolute differences of 5.1%–27.7%). Disc height had a systematic influence on the internal disc stresses and tractions at the disc‐to‐bone interface. The influence of disc geometry on mechanics highlights the importance of disc‐specific modeling to estimate disc injury risk, loading on the adjacent vertebral bodies, and the mechanical environment present in disc tissues.

Background Cadaveric intervertebral discs are often studied for a variety of research questions, and outcomes are interpreted in the in vivo context. Unfortunately, the cadaveric disc does not inherently represent the LIVE condition, such that the disc structure (geometry), composition (T2 relaxation time), and mechanical function (opening pressure, OP) measured in the cadaver do not necessarily represent the in vivo disc. Methods We conducted serial evaluations in the Yucatan minipig of disc geometry, T2 relaxation time, and OP to quantify the changes that occur with progressive dissection and used axial loading to restore the in vivo condition. Results We found no difference in any parameter from LIVE to TORSO; thus, within 2 h of sacrifice, the TORSO disc can represent the LIVE condition. With serial dissection and sample preparation the disc height increased (SEGMENT height 18% higher than TORSO), OP decreased (POTTED was 67% lower than TORSO), and T2 time was unchanged. With axial loading, an imposed stress of 0.20–0.33 MPa returned the disc to in vivo, LIVE disc geometry and OP, although T2 time was decreased. There was a linear correlation between applied stress and OP, and this was conserved across multiple studies and species. Conclusion To restore the LIVE disc state in human studies or other animal models, we recommend measuring the OP/stress relationship and using this relationship to select the applied stress necessary to recover the in vivo condition.