Geometric determinants of the mechanical behavior of image-based finite element models of the intervertebral disc

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² = 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.