The finite element method is widely used for studying the intervertebral disc at the organ level due to its ability to model complex geometries. An indispensable requirement for proper modelling of the intervertebral disc is a reliable porohyperelastic framework that captures the elaborate underlying mechanics. The increased complexity of such models requires significant computational power that is available within high-performance computing systems. 

The objective of this study is to present such a framework, validated both against literature and experiments, aiming to enable intervertebral disc research to benefit from state-of-the-art computational resources.