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.

[1]  A. Freemont,et al.  End-Plate Displacement During Compression of Lumbar Vertebra-Disc-Vertebra Segments and the Mechanism of Failure , 1993, Spine.

[2]  Julia C. O'Neill,et al.  Introduction to the Design and Analysis of Experiments , 2001, Technometrics.

[3]  G A Dumas,et al.  Influence of material properties on the mechanical behaviour of the L5-S1 intervertebral disc in compression: a nonlinear finite element study. , 1991, Journal of biomedical engineering.

[4]  V K Goel,et al.  Effect of Disc Degeneration at One Level on the Adjacent Level in Axial Mode , 1991, Spine.

[5]  Madhan Shridhar Phadke,et al.  Quality Engineering Using Robust Design , 1989 .

[6]  Richard M Aspden,et al.  Statistical methods in finite element analysis. , 2002, Journal of biomechanics.

[7]  G B Andersson,et al.  A model to study the disc degeneration process. , 1994, Spine.

[8]  A M Mohsen,et al.  Patient-specific spine models. Part 1: Finite element analysis of the lumbar intervertebral disc—a material sensitivity study , 2002, Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine.

[9]  Richard O. Lynch,et al.  Introduction to Design and Analysis of Experiments , 1999, Technometrics.

[10]  J P Laible,et al.  A Poroelastic-Swelling Finite Element Model With Application to the Intervertebral Disc , 1993, Spine.

[11]  N. Broom,et al.  Biomechanics of load-bearing of the intervertebral disc: an experimental and finite element model. , 1997, Medical engineering & physics.

[12]  D W Hukins,et al.  Replacing the nucleus pulposus of the intervertebral disc. , 2001, Clinical biomechanics.

[13]  P. Brinckmann,et al.  Interlaminar Shear Stresses and Laminae Separation in a Disc: Finite Element Analysis of the L3‐L4 Motion Segment Subjected to Axial Compressive Loads , 1995, Spine.

[14]  A B Schultz,et al.  Finite element stress analysis of an intervertebral disc. , 1974, Journal of biomechanics.

[15]  J. Meakin,et al.  Replacing the nucleus pulposus of the intervertebral disk: prediction of suitable properties of a replacement material using finite element analysis , 2001, Journal of materials science. Materials in medicine.