Sensitivity of Intervertebral Disc Finite Element Models to Internal Geometric and Non-geometric Parameters

Finite element models are useful for investigating internal intervertebral disc (IVD) behaviours without using disruptive experimental techniques. Simplified geometries are commonly used to reduce computational time or because internal geometries cannot be acquired from CT scans. This study aimed to (1) investigate the effect of altered geometries both at endplates and the nucleus-anulus boundary on model response, and (2) to investigate model sensitivity to material and geometric inputs, and different modelling approaches (graduated or consistent fibre bundle angles and glued or cohesive inter-lamellar contact). Six models were developed from 9.4 T MRIs of bovine IVDs. Models had two variations of endplate geometry (a simple curved profile from the centre of the disc to the periphery, and precise geometry segmented from MRIs), and three variations of NP-AF boundary (linear, curved, and segmented). Models were subjected to axial compressive loading (to 0.86 mm at a strain rate of 0.1/s) and the effect on stiffness and strain distributions, and the sensitivity to modelling approaches was investigated. The model with the most complex geometry (segmented endplates, curved NP-AF boundary) was 3.1 times stiffer than the model with the simplest geometry (curved endplates, linear NP-AF boundary), although this difference may be exaggerated since segmenting the endplates in the complex geometry models resulted in a shorter average disc height. Peak strains were close to the endplates at locations of high curvature in the segmented endplate models which were not captured in the curved endplate models. Differences were also seen in sensitivity to material properties, graduated fibre angles, cohesive rather than glued inter-lamellar contact, and NP:AF ratios. These results show that FE modellers must take care to ensure geometries are realistic so that load is distributed and passes through IVDs accurately.