Regional variations of mechanical responses of IVD to 7 different motions: An in vivo study combined with FEA and DFIS

The abnormal mechanical behaviour of a lumbar intervertebral disc (IVD) is commonly recognized as a direct indicator of intervertebral disc degeneration (IDD). However, current methods cannot evaluate the patient-specific mechanical performance of an IVD in vivo during movement. This study establishes a patient-specific (PS) model that combines the kinematics parameters of the lumbar spine obtained with a dual fluoroscopic imaging system (DFIS) and a finite-element (FE) method for the first time to reveal the mechanical behaviours of IVDs in vivo under seven motions. Three healthy participants were recruited for this study. CT images were obtained to create finite-element models of L3-L5 spine segments. Meanwhile, participants were required to take specific positions including upright standing, flexion, extension, left and right lateral bending, as well as left and right axial torsion in the DFIS. The in vivo kinematic parameters, calculated by registering the CT images with images obtained with DFIS, were considered as loading conditions in FE simulations. Significant differences of von Mises stresses and principal strains were found between PS model and GN model which employing a generalized moment as loading conditions, former resulting in up to 76.74 % lower strain than the GN model. Also, considerable differences were observed for five anatomical regions of the IVD (L3-L5). Under all motions, the stress in the centre region (nucleus pulposus) was the lowest, while the stress in the posterior region was the highest in extension motion. Therefore, activities such as stretching with an extension, should be avoided by patients with a herniated disc, in which the posterior region was the herniation site. The PS model combining in vivo kinematics and FE simulations shows the potential in the design and assessment of patient-specific implants.