[Finite element modeling and simulation study of solid-liquid biphase fiber-reinforced lumbar intervertebral disc]

The lumbar intervertebral disc exhibits a complex physiological structure with interactions between various segments, and its components are extremely complex. The material properties of different components in the lumbar intervertebral disc, especially the water content (undergoing dynamic change as influenced by age, degeneration, mechanical loading, and proteoglycan content) - critically determine its mechanical properties. When the lumbar intervertebral disc is under continuous pressure, water seeps out, and after the pressure is removed, water re-infiltrates. This dynamic fluid exchange process directly affects the mechanical properties of the lumbar intervertebral disc, while previous isotropic modeling methods have been unable to accurately reflect such solid-liquid phase behaviors. To explore the load-bearing mechanism of the lumbar intervertebral disc and establish a more realistic mechanical model of the lumbar intervertebral disc, this study developed a solid-liquid biphasic, fiber-reinforced finite element model. This model was used to simulate the four movements of the human lumbar spine in daily life, namely flexion, extension, axial rotation, and lateral bending. The fluid pressure, effective solid stress, and liquid pressure-bearing ratio of the annulus fibrosus and nucleus pulposus of different lumbar intervertebral discs were compared and analyzed under the movements. Under all the movements, the fluid pressure distribution was closer to the nucleus pulposus, while the effective solid stress distribution was more concentrated in the outer annulus fibrosus. In terms of fluid pressure, the maximum fluid pressure of the lumbar intervertebral disc during lateral bending was 1.95 MPa, significantly higher than the maximum fluid pressure under other movements. Meanwhile, the maximum effective solid stress of the lumbar intervertebral disc during flexion was 2.43 MPa, markedly higher than the maximum effective solid stress under other movements. Overall, the liquid pressure-bearing ratio under axial rotation was smaller than that under other movements. Based on the solid-liquid biphasic modeling method, this study more accurately revealed the dominant role of the liquid phase in the daily load-bearing process of the lumbar intervertebral disc and the solid-phase mechanical mechanism of the annulus fibrosus load-bearing, and more effectively predicted the solid-liquid phase co-load-bearing mechanism of the lumbar intervertebral disc in daily life.