Replacing the nucleus pulposus of the intervertebral disk: prediction of suitable properties of a replacement material using finite element analysis

An axisymmetric finite element model of a human lumbar disk was developed to investigate the properties required of an implant to replace the nucleus pulposus. In the intact disk, the nucleus was modeled as a fluid, and the annulus as an elastic solid. The Young's modulus of the annulus was determined empirically by matching model predictions to experimental results. The model was checked for sensitivity to the input parameter values and found to give reasonable behavior. The model predicted that removal of the nucleus would change the response of the annulus to compression. This prediction was consistent with experimental results, thus validating the model. Implants to fill the cavity produced by nucleus removal were modeled as elastic solids. The Poisson's ratio was fixed at 0.49, and the Young's modulus was varied from 0.5 to 100 MPa. Two sizes of implant were considered: full size (filling the cavity) and small size (smaller than the cavity). The model predicted that a full size implant would reverse the changes to annulus behavior, but a smaller implant would not. By comparing the stress distribution in the annulus, the ideal Young's modulus was predicted to be approximately 3 MPa. These predictions have implications for current nucleus implant designs.