Statistical factorial analysis on the poroelastic material properties sensitivity of the lumbar intervertebral disc under compression, flexion and axial rotation

A statistical factorial analysis approach was conducted on a poroelastic finite element model of a lumbar intervertebral disc to analyse the influence of six material parameters (permeabilities of annulus, nucleus, trabecular vertebral bone, cartilage endplate and Young's moduli of annulus and nucleus) on the displacement, fluid pore pressure and velocity fields. Three different loading modes were investigated: compression, flexion and axial rotation. Parameters were varied considering low and high levels in agreement with values found in the literature for both healthy and degenerated lumbar discs. Results indicated that annulus stiffness and cartilage endplate permeability have a strong effect on the overall fluid- and solid-phase responses in all loading conditions studied. Nucleus stiffness showed its main relevance in compression while annulus permeability influenced mainly the annular pressure field. This study confirms the permeability's central role in biphasic modelling and highlights for the lumbar disc which experiments of material property characterization should be performed. Moreover, such sensitivity study gives important guidelines in poroelastic material modelling and finite element disc validation.