To provide mechanistic insight into the etiology of osteoporotic wedge fractures, we investigated the spatial distribution of tissue at the highest risk of initial failure within the human vertebral body for both forward flexion and uniform compression loading conditions. Micro-CT-based linear elastic finite element analysis was used to virtually load 22 human T9 vertebral bodies in either 5° of forward flexion or uniform compression; we also ran analyses replacing the simulated compliant disc (E=8 MPa) with stiff polymethylmethacrylate (PMMA, E=2500 MPa). As expected, we found that, compared to uniform compression, forward flexion increased the overall endplate axial load on the anterior half of the vertebra and shifted the spatial distribution of high-risk tissue within the vertebra towards the anterior aspect of the vertebral body. However, despite that shift, the high-risk tissue remained primarily within the central regions of the trabecular bone and endplates, and forward flexion only slightly altered the ratio of cortical-to-trabecular load sharing at the mid-vertebral level (mean±SD for n=22: 41.3±7.4% compression; 44.1±8.2% forward flexion). When the compliant disc was replaced with PMMA, the anterior shift of high-risk tissue was much more severe. We conclude that, for a compliant disc, a moderate degree of forward flexion does not appreciably alter the spatial distribution of stress within the vertebral body.

The mechanism of vertebral wedge fractures remains unclear and may relate to typical variations in the mechanical behavior of the intervertebral disc. To gain insight, we tested 16 individual whole discs (between levels T8 and L5) from nine cadavers (mean±SD: 66±16 years), loaded in compression at different rates (0.05-20.0% strain/s), to measure a homogenized "effective" linear elastic modulus of the entire disc. The measured effective modulus, and average disc height, were then input and varied parametrically in micro-CT-based finite element models (60-μm element size, up to 80 million elements each) of six T9 human vertebrae that were virtually loaded to 3° of moderate forward-flexion via a homogenized disc. Across all specimens and loading rates, the measured effective modulus of the disc ranged from 5.8 to 42.7MPa and was significantly higher for higher rates of loading (p<0.002); average disc height ranged from 2.9 to 9.3mm. The parametric finite element analysis indicated that, as disc modulus increased and disc height decreased across these ranges, the vertebral bone stresses increased but their spatial distribution was largely unchanged: most of the highest stresses occurred in the central trabecular bone and endplates, and not anteriorly. Taken together with the literature, our findings suggest that the effective modulus of the human intervertebral disc should rarely exceed 100MPa and that typical variations in disc effective modulus (and less so, height) minimally influence the spatial distribution but can appreciably influence the magnitude of stress within the vertebral body.