Biomechanical study of intervertebral disc degeneration

Degeneration and age affect the biomechanics of the intervertebral disc, by reducing its stiffness, flexibility and shock absorption capacities against daily movement and spinal load. The biomechanical characterization of intervertebral discs is achieved by conducting mechanical testing to vertebra-disc-vertebra segments and applying axial, shear, bend and torsion loads, statically or dynamically, with load magnitudes corresponding to the physiological range. However, traditional testing does not give a view of the load and deformation states of the disc components: nucleus pulposus, annulus fibrosus and endplate. Thus, the internal state of stress and strains of the disc can only be predicted by numerical methods, one of which is the finite element method. The objective of this thesis was, to study the biomechanics of degenerated intervertebral discs to load conditions in compression, bending and torsion, by using mechanical testing and a finite element model of disc degeneration, based on magnetic resonance imaging (MRI). Therefore, lumbar discs obtained from cadavers corresponding to spinal levels L2-L3 and L4-L5 with mild to severe degeneration were used. Intervertebral osteochondrosis and spondylosis deformans were identified, being the disc space collapse, the most striking feature. Next, all discs were tested to static and dynamic load conditions, the results gained corresponded to the disc stiffness (in compression, bending and torsion), stress relaxation and dynamic response. Of these, the stiffness response was used to validate the disc model. The testing results suggest that discs with advanced degeneration over discs with mild degeneration are, less rigid in compression, less stiffer under bending and torsion, showed less radial bulge, and reduce their viscoelastic and damping properties. This study shows that degeneration has an impact on the disc biomechanical properties which can jeopardize normal functionality. Development of one finite element model of disc degeneration started by choosing a MRI of a L2-L3 disc. Segmentation of vertebra bone and disc materials followed, and were based on pixel brightness and radiology fundamentals, then a finite element mesh was created to account for the disc irregular shape. The disc materials were modeled as hyperelastic and the bone materials were modeled as orthotropic and isotropic. Adjustment of material properties was based on integrity of the annulus fibrosus, giving a stiffness value matching that of a mild degeneration disc. Then, validation of the model was performed, and included a study of the distributions of stress and strain under loads of compression, bending and torsion. The results from all load simulations show that the disc undergoes large deformations. In contrast, the vertebrae are subjected to higher stress but with negligible deformations. In compression, the model predicted formation of symmetrical disc bulge which agree with the testing behavior. The nucleus pulposus showed to be the principal load carrier with negative principal stresses and strains. In bending and torsion, the annulus fibrosus showed to be the principal load carrier with large symmetrical principal strains and stresses for the former loading and large shearing for the latter. The study showed the importance of soft tissue deformation, mostly noticed in advanced degeneration. In contrast, the higher stresses in the vertebra over those of the intervertebral disc showed the relevance of bone predisposition to fracture. Such kind of studies, should contribute to the understanding of the biomechanics of the intervertebral disc.La degeneración y edad afectan la biomecánica del disco intervertebral, reduciendo la capacidad de rigidez, flexibilidad y atenuación de impactos, contra el movimiento y carga del raquis. La caracterización biomecánica del disco se realiza con ensayos mecánicos a segmentos de vértebra-disco-vértebra y aplicando cargas axiales, cortantes, flexión y torsión, estáticas ó dinámicas, con magnitudes de carga según el intervalo fisiológico. Sin embargo, las pruebas tradicionales no dan una visión de los estados de carga y deformación de los componentes del disco: núcleo pulposo, anillo fibroso y placa terminal. Por lo tanto, el estado interno de esfuerzos y deformaciones del disco, solo puede ser predicho con métodos numéricos, uno de los cuales es el método de elemento finito. El objetivo de esta tesis fue, estudiar la biomecánica de discos intervertebrales degenerados a las condiciones de carga en compresión, flexión y torsión, mediante el uso de ensayos mecánicos y de un modelo de elementos finitos de la degeneración de disco, basado en imágenes con resonancia magnética (MRI). Por lo tanto, se usaron discos lumbares L2-L3 y L4-L5 obtenidos de cadáveres, con degeneración leve a severa. Se identificó osteocondrosis intervertebral y espondilosis deformante, siendo el colapso del espacio intervertebral el aspecto más relevante. Luego, todos los discos fueron ensayados a condiciones de carga estática y dinámica, y los resultados correspondieron a la rigidez del disco (a compresión, flexión y torsión), a la relajación de tensiones y a la respuesta dinámica. De éstos, la rigidez fue usada para validar el modelo de disco. Los resultados de los ensayos sugieren que los discos con degeneración avanzada sobre aquellos con degeneración leve son, menos rigidos a compresión, menos rigidos a flexión y torsión, presentan menor protuberancia radial, y reducen sus propiedades viscoelásticas y de amortiguamiento. El estudio muestra que la degeneración impacta las propiedades biomecánicas del disco, poniendo en riesgo la funcionalidad normal. El desarollo de un modelo de elementos finitos de la degeneración de disco inició eligiendo una secuencia de resonancia magnética de un disco L2-L3. La segmentación de los materiales del disco y de las vértebras se realizó basado en intensidad de brillo del pixel y en fundamentos de radiología, y se creó una malla de elementos finitos correspondiente a la forma irregular del disco. Los materiales del disco se modelaron como hiperelásticos y los tejidos óseos se modelaron como materiales ortotrópicos e isotrópicos. El ajuste de propiedades de los materiales fue basado en la integridad del anillo fibroso, y dio una rigidez correspondiente a la de un disco con degeneración leve. Luego, se realizó la validación del modelo, e incluyó un estudio de las distribuciones de esfuerzo y deformación a las condiciones de carga en compresión, flexión y torsión. Los resultados de todas las simulaciones de carga mostraron que el disco es sometido a grandes deformaciones. En contraste, las vértebras fueron sometidas a mayores esfuerzos pero con deformaciones insignificantes. En compresión, el modelo predijo la formación de una protuberancia radial simétrica, en concordancia con la experimentación. El núcleo pulposo mostró ser el portador principal de carga, con tensiones y deformaciones principales negativas. En flexión y torsión, el anillo fibroso mostró ser el portador principal de carga, con grandes deformaciones y tensiones principales simétricas para la primera carga, y con grandes tensiones cortantes para la segunda carga. El estudio mostró la importancia de las deformaciones de los tejidos blandos, principalmente notados en la degeneración avanzada. Por el contrario, las tensiones mayores en los cuerpos vertebrales sobre aquellas del disco intervertebral mostraron la relevancia de la predisposición a las fracturas óseas. Este tipo de estudio debe contribuir a la comprensión de la biomecánica del disco intervertebral