Application of the polynomial chaos expansion to approximate the homogenised response of the intervertebral disc