The current study measures the mechanical behavior of both natural and the monobloc elastomeric disc prosthesis (CadiscTM-L) by employing a finite element method (FEM) to study the fiber-reinforced constitutive formulation provided in the literature. The three-dimensional geometry was created by computed tomography (CT) scan imaging technique. Frontal pure rotational, sagittal, and axial momentum of 7.5 N·m were applied on the top of L3 while the lower half of the L5 was fixed in all directions. This investigation was performed considering two stages: (1) intact L3–L5 lumbar spine (INT model), and (2) Cadisc implemented between L4 and L5 (IMP model). The numerical results for the INT model were validated by experimental data from the literature. Several parameters including the inter-segmental rotation, range of motion in flexion-extension, axial rotation and lateral bending were analyzed. Our numerical results show that the IMP model has a 50% reduction in the ‘range of motion’ and a 33% reduction in flexion in lateral bending compared to the INT model. These outcomes of this paper reveal the feasibility of applying a fibre-reinforced constitutive formulation to generate an accurate three-dimensional FEM model.

The investigation of dynamic response of intervertebral disc is beneficial for the development of new synthetic and engineered tissues for treating diseased or injured disc. There are limited experimental studies on comparing the effect of loading mode and rate on global response of intervertebral disc. In this study, in-vitro experiments were performed using a total of 24 porcine motion segments. The harvested specimens were assigned to prolong and 2 different cyclic loadings. Both disc deformations and water contents were measured to investigate how the mode and rate of loading affect the response of intervertebral disc. In parallel, a backward FE poroelastic model combined with in-vitro experiments were used to find the material properties of intervertebral discs. The experimental result showed that the final disc height loss under creep loading was significantly greater than cyclic groups. Increasing the frequency of cyclic loading decreased the disc height loss. The water content decreased significantly in cyclic loading from those in prolong loading. The backward FE models showed that, the elastic modulus of anulus fibrosus and nucleus pulposus were 2.43 (±0.48) MPa and 1.46 (±0.29) MPa, respectively. The hydraulic permeability was 2.08 (±0.42) ×10-16m4/Ns, and the Poisson’s ratio was 0.21 (±0.03). In conclusion, this study investigated how the loading mode and rate affect porcine intervertebral disc deformation. It is found that dynamic stiffness is greater at higher frequencies which resulted from interactions between the solid phase and fluid flow within the disc.