Stochastic finite element analysis of biological systems: comparison of a simple intervertebral disc model with experimental results

Statistical methods allow the effects of uncertainty to be incorporated into finite element models. This has potential benefits for the analysis of biological systems where natural variability can give rise to substantial uncertainty in both material and geometrical properties. In this study, a simple model of the intervertebral disc under compression was created and analysed as both a deterministic and a stochastic system. Factorial analysis was used to determine the important parameters to be included in the stochastic analysis. The predictions from the model were compared to experimental results from 21 sheep discs. The size and shape of the distribution of the axial deformations predicted by the model was consistent with the experimental results given that the number of model solutions far exceeded the number of experimental results. Stochastic models could be valuable in determining the range and most likely value of stress in a tissue or implant.