Simulation of Inhomogeneous Rather Than Homogeneous Poroelastic Tissue Material Properties Within Disc Annulus and Nucleus Better Predicts Cervical Spine Response

Abstract

A finite element (FE) modeling of homogeneous and inhomogeneous poroelastic tissue material properties within disc anulus fibrosus (AF) and nucleus pulposus (NP). To test the hypothesis that simulation of inhomogeneous poroelastic tissue material properties within AF and NP quadrants, rather than homogeneous properties within regions of AF and NP without quadrants, would better predict the cervical spine biomechanics. In order to represent tissue swelling and creep deformation behavior more physiologically in FE models, disc poroelastic tissue material properties should be modeled appropriately. Past studies show an existence of inhomogeneous rather than homogeneous nature of the tissue properties in various quadrants of AF and NP, and this has been simulated in a single-segment FE lumbar model with only compression analysis. This article simulated these tissue properties in a multisegmental cervical spine and reported the results of both compression and moment loads. Two three-dimensional FE models of a C3-T1 segment were developed. Model I included homogeneous poroelastic tissue properties in AF and NP, whereas Model II included inhomogeneous poroelastic tissue properties in AF and NP quadrants. Biomechanical responses of the FE models under diurnal compression and moment loads were compared with corresponding in vivo published studies. Model II with disc quadrant-based inhomogeneous poroelastic tissue properties predicted better, mainly in flexion and extension, than the Model I with homogeneous tissue properties when compared with the corresponding in vivo results, thereby confirming the current study hypothesis. Inhomogeneous tissue properties govern segmental behavior mainly during sagittal plane motions, with a root-mean-square difference of nearly 50% across the motion segments. The current data justify the need to simulate inhomogeneous tissue properties within disc quadrants for any FE model analysis. Model II can be further used to understand the biomechanical effects of quadrant-based degenerative poroelastic tissue properties on cervical spine behavior. Future experiments are necessary to support the current study results.