Au cours des deux dernières décennies, la simulation computationnelle – en particulier l'analyse des éléments finis (AFE) – est devenue un outil indispensable pour les chercheurs qui cherchent à prédire comment les disques intervertébraux réagissent aux conditions de charge complexes et réelles. En offrant une fenêtre non invasive sur les contraintes, les souches et les flux de fluides internes qui se produisent dans le disque, ces modèles aident à expliquer les mécanismes de blessure, guident la conception des implants et perfectionnent les protocoles chirurgicaux et de réadaptation.

Anatomie et fonction du disque intervertébral

Le disque intervertébral est une articulation fibrocartilagineuse spécialisée qui se trouve entre les corps vertébraux adjacents dans la colonne vertébrale. Ses rôles principaux sont de supporter et de distribuer les charges axiales, de permettre un mouvement limité (flexion, extension, flexion latérale et rotation) et d'absorber les chocs lors des activités quotidiennes telles que la marche, le levage et le saut.

Nucleus Pulposus

Dans un disque jeune et sain, il contient environ 70 à 90 % d'eau en poids, avec des protéoglycans (principalement aggrécan), des fibres de collagène de type II et de l'élastine. Cette forte teneur en eau fait du noyau un fluide essentiellement incompressible qui pressurise sous compression, lui permettant de résister aux forces verticales et de répartir radialement la pression contre l'annulaire environnant. Le noyau présente également un comportement viscoélastique : sous charge prolongée il perd du liquide (creep) et, au déchargement, il réhydrate (recovery).

Annulus Fibrosus

L'annulaire fibrosus forme l'anneau externe dur du disque. Il est composé de 15 à 25 lamelles concentriques, chacune composée principalement de fibres collagènes de type I orientées à des angles alternés (environ ±30° à l'horizontale). Cette structure lamellaire donne à l'annulaire une résistance exceptionnelle aux contraintes ténsiles et torsionnelles. L'annulaire a un double rôle : il contient le noyau pressurisé et assure la stabilité structurelle au segment spinal. L'annulaire interne se transforme progressivement vers une composition plus fibrocartilagineuse semblable au noyau, tandis que l'annulaire externe est plus dense et fibreux. Les déchirures ou fissures de l'annulaire peuvent conduire à l'hernie du disque, où la matière nucléaire pénètre vers l'extérieur et peut empiéter sur les racines nerveuses.

Plaques d'extrémité cartiagineuses

Les plaques cartilagineuses sont de fines couches de cartilage hyaline qui couvrent les surfaces supérieures et inférieures du disque, le séparant des corps vertébraux. Elles servent d'interface perméable pour l'échange de nutriments par diffusion de l'apport sanguin vertébral, un processus essentiel pour la santé du disque parce que le disque mature est avasculaire. Les plaques d'extrémité jouent également un rôle mécanique: elles distribuent uniformément les charges de compression sur le disque et aident à empêcher l'extrusion du noyau dans les vertébraux. La calcification ou les dommages aux plaques d'extrémité peuvent perturber le transport nutritionnel et accélérer la dégénérescence du disque.

Types de conditions de chargement sur le disque spinal

Pour simuler avec précision le comportement des disques, les chercheurs considèrent plusieurs modes de chargement fondamentaux, souvent combinés pour recréer des conditions réalistes in vivo. Ci-dessous est une discussion élargie des types de chargement primaire.

Compression

La compression axiale résulte d'activités telles que la position debout, l'assise et le levage. Pendant la compression, le noyau pulposus pressurise, générant une pression hydrostatique qui pousse vers l'extérieur contre les fibres annulaires. L'annulaire subit alors une tension dans la direction circonférentielle (boucle). Dans un disque sain, la distribution de la pression est relativement uniforme, mais la dégénérescence modifie cette relation.

Torsion (Rotation axiale)

La torsion se produit lorsque la colonne vertébrale est tordue, par exemple lors d'une oscillation de golf ou en tournant le torse. Comme les fibres de collagène du disque sont déjà sous tension par la pression du noyau, des charges torsionnelles supplémentaires font que les lamelles annulaires sont étirées d'un côté et comprimées de l'autre. L'angle de fibre (=30°) est optimisé pour résister à ces moments torsionnels. Cependant, une torsion excessive ou répétitive est un facteur de risque connu pour les larmes annulaires : les simulations indiquent que les contraintes maximales se concentrent dans la région postolatérale de l'annulaire, qui est également le site le plus commun pour l'hernie du disque.

Flexion et extension

Pendant la flexion, la partie antérieure des compresses de disque pendant que la région postérieure est placée en tension. Ce motif de contrainte asymétrique peut pousser le noyau pupossus vers l'arrière. Au fil du temps, la charge répétitive de flexion contribue à la dégénérescence du disque et à la défaillance annulaire postérieure. Inversement, pendant l'extension, l'annulaire postérieur est comprimé et l'annulaire antérieur est tendu. Les cycles de flexion-extension impliquent souvent une combinaison de moment de flexion et de compression axiale, ce qui les rend plus complexes à simuler.

Chargement du karité

Les forces de cisaillement se produisent lors de mouvements comme la lune avant ou le levage avec un dos arrondi; le cisaillement latéral se produit lors de la flexion latérale. La rigidité du cisaillement du disque est considérablement plus faible que sa rigidité compressive, ce qui le rend vulnérable aux blessures sous des charges de cisaillement élevées, surtout lorsqu'il est combiné à la flexion de la colonne vertébrale. De nombreuses blessures à la main résultent d'une combinaison de compression, de flexion et de cisaillement.

Chargement combiné

Dans Vivo, la charge est presque toujours multiaxiale. Par exemple, le levage d'une boîte du plancher implique la compression (du poids de l'objet et du corps supérieur), la flexion (tendance vers l'avant) et le cisaillement (inclinaison vers l'avant).Les simulations avancées intègrent ces conditions combinées en appliquant simultanément des vecteurs de force et de moment.

Techniques de simulation: Analyse des éléments finis en détail

L'analyse des éléments finis (FEA) reste la méthode de calcul la plus utilisée pour étudier la mécanique des disques. La technique discrétise la géométrie des disques en milliers ou des millions de petits éléments, chacun des propriétés spécifiques des matériaux assignés. En résolvant des ensembles d'équations différentielles partielles, FEA prédit le déplacement, la contrainte, la tension et la pression fluide dans tout le modèle sous des charges appliquées.

Géométrie et génération de mesh

Les modèles précoces ont utilisé des formes symétriques idéales (par exemple, des disques ellipsoïdaux), mais les modèles modernes tirent souvent des géométries spécifiques au patient des balayages par résonance magnétique (IRM) ou par tomographie (CT). Le noyau et l'annulaire sont généralement représentés en volumes séparés, et l'annulaire est subdivisé en plusieurs lamelles pour capturer l'orientation des fibres. La densité de mesh doit être suffisante pour résoudre les gradients de stress, surtout à proximité des plaques d'extrémité et dans l'annulaire postolatéral, où la défaillance commence souvent. Certains modèles comprennent également les vertèbres adjacentes (os corticaux, os annulés et plaques d'extrémité) et les ligaments (antérieurs longitudinales, postérieurs longitudinales, interspineux, etc.). Un modèle classique d'éléments finis de disques lombaires peut contenir de 100 000 à 500 000 éléments, selon la complexité.

Propriétés du matériau

Le disque présente un comportement complexe, non linéaire et anisotropique. Le noyau est souvent modélisé comme un matériau hyperélastique presque incompressible (par exemple, utilisant une fonction d'énergie de contrainte Mooney-Rivlin ou Ogden) combiné à une poroélasticité fluide. L'annulaire est modélisé comme un composite renforcé par des fibres : la substance du sol est hyperélastique, et les fibres de collagène sont représentées comme des ressorts à tension seulement ou comme une couche de renforcement intégrée de type rebar. Les fibres sont dotées d'une couche initiale (undulation) qui se déploie comme une augmentation de la déformation, ce qui entraîne une réponse non linéaire de contrainte-souche.

Les paramètres matériels sont généralement dérivés d'études expérimentales sur les tissus cadavériques humains. Les références clés sont : Shirazi-Adl et al. (1986), Schmidt et al. (2006) et O'Connell et al. (2012). Ces travaux fournissent des valeurs de base pour la rigidité de la fibre de collagène, la pression osmotique protéoglycane et les coefficients de perméabilité.

Conditions de limite et de chargement

Les conditions de bordure doivent reproduire l'environnement physiologique. Habituellement, la plaque inférieure de la vertèbre inférieure est entièrement fixe, tandis que les charges et/ou les déplacements sont appliqués à la plaque supérieure de la vertèbre supérieure. Le chargement peut être appliqué comme force (commandé par la charge) ou comme déplacement (commandé par la force). Pour simuler le mouvement physiologique, les charges des suiveurs, qui sont approximatives de la force de compression générée par les muscles du tronc, sont souvent superposées.

Validation et vérification

La validation consiste généralement à comparer la pression prévue du disque, la pression intradiscale, la plage segmentaire de mouvement et les profils de défaillance par rapport aux données in vitro publiées. Par exemple, un modèle validé devrait montrer que la pression intradiscale dans le disque L4-L5 de la lombaire inférieure à 1000 N est d'environ 0,5–0,7 MPa, ce qui est conforme à Wilke et al. (1992). La vérification garantit que la solution numérique est exacte et indépendante du maillage.

Principales constatations tirées des simulations

Des décennies de recherche FEA ont produit un riche corpus de connaissances sur la mécanique des disques. Ci-dessous sont les résultats clés élargis qui vont au-delà des points de balle de l'article original.

Chargement compressif : interaction en phase liquide et en phase solide

Sous compression pure, le noyau pressurise, provoquant une bourrage radiale de l'annulaire. Cette bourre est maximale dans le plan médian transverse du disque. Les fibres annulaires se redressent progressivement et s'étirent, portant une grande partie de la tension du cerceau. Les modèles poroélastiques révèlent que sous compression constante soutenue, le fluide exude du noyau à travers les plaques d'extrémité et l'annulaire, entraînant une réduction progressive de la hauteur du disque (creep). Après enlèvement de charge, le fluide imbibe en arrière, rétablissant la hauteur. Ces effets dépendants du temps sont critiques pour comprendre la dégénérescence du disque : les disques dégénérés ont une teneur en protéoglyque inférieure, une pression réduite et une tendance plus grande à la bourre diffuse.

Torsion: Risque de déformation et de déchirure de la fibre annulaire

La charge torsionnelle produit les plus fortes contraintes dans l'annulaire fibrosus, en particulier dans la région postolatérale. Les fibres du côté opposé à la direction de rotation (le « côté de la tension ») deviennent très tendues, tandis que les fibres du même côté (le « côté de la compression ») peuvent se boucler. La torsion répétée peut entraîner une délamination circonférentielle – une séparation entre les lamelles adjacentes. Les modèles d'éléments Finite qui intègrent des critères de dommages progressifs ont prédit qu'un seul événement de torsion de > 20° de rotation dans un disque dégénéré peut déclencher des larmes. Ces résultats s'alignent sur les observations cliniques selon lesquelles la charge torsion est un facteur de risque majeur pour l'hernie des disques chez les athlètes et les travailleurs manuels.

Flexion et extension : Stress asymétrique et risque de surtrusion

Pendant la flexion, l'annulaire antérieur est comprimé alors que l'annulaire postérieur est tendu. Le noyau migre postérieurement, et si l'annulaire postérieur est affaibli ou dégénéré, cela peut conduire à un gonflement ou une hernie de disque postérieur. Les simulations ont montré que le stress maximal dans l'annulaire postérieur sous flexion-compression combinée peut dépasser le seuil de défaillance tissulaire à 50-70% de la gamme de mouvements du segment de mouvement dans un disque dégénéré. En extension, l'annulaire postérieur est comprimé et l'annulaire antérieur est tendu; cela peut exacerber les larmes annulaires antérieures mais est moins souvent associé à l'hernie. De nombreuses études FEA modernes intègrent également les articulations facettes, qui portent jusqu'à 20% de la charge compressive dans l'extension et aident à limiter le mouvement torsionnel.

Couche : Vulnérabilité de l'interface Nucleus-Annulus

La charge de cisaillement, surtout lorsqu'elle est combinée à la compression axiale, induit des déformations de cisaillement à l'interface entre le noyau et l'annulaire interne. Cette région est un plan naturel de faiblesse et la délamination peut commencer là. Des études paramétriques ont montré que le rétrécissement de la hauteur du disque – une caractéristique de dégénérescence – élève encore les contraintes de cisaillement.

Demandes et répercussions cliniques

L'objectif ultime de la simulation du comportement mécanique des disques est d'améliorer les résultats des patients. Ici, nous mettons en évidence les applications clés dans la conception des appareils, la planification chirurgicale et la réadaptation.

Conception totale de remplacement de disque

L'avènement d'implants de remplacement de disques (TDR) (par exemple Charité, ProDisc, M6) a incité l'utilisation de FEA pour évaluer l'influence de ces dispositifs sur le transfert de charge vers les vertèbres adjacentes et les éléments postérieurs. Des simulations ont été utilisées pour optimiser les surfaces portantes de l'implant, les matériaux (chrome de cobalt, polyéthylène ultra-haut-moléculaire) et les géométries de la quille/stem. Un TDR bien conçu devrait rétablir le mouvement segmentaire et la pression intradiscale, réduisant ainsi le risque de maladie de segment adjacent. Par exemple, des études FEA ont montré que la surdimensionnement de l'implant ou son placement trop a posteriori peut causer une surcharge et une hypermobilité faceteuses.

Réadaptation et ergonomie

En simulant les différences entre un disque sain et un disque dégénéré sous divers schémas d'activité, les chercheurs peuvent identifier les mouvements qui imposent le plus de risque.Cette connaissance informe les protocoles de réhabilitation : par exemple, des exercices qui produisent des charges de torsion élevées ou une flexion profonde peuvent être contre-indiqués pour les patients ayant une dégénérescence précoce du disque.

Dégénérescence progressive et intégration biologique

Les recherches actuelles se tournent vers des modèles multi-échelles et multiphysiques qui intègrent non seulement la charge mécanique mais aussi la biologie : renouvellement de matrices médiée par les cellules, diffusion des nutriments et gonflement osmotique.Ces modèles « métabiologiques » peuvent simuler comment une modification de la charge – telle que celle causée par un mode de vie sédentaire ou l'obésité – peut déclencher une cascade de réponses cellulaires conduisant à une dégénérescence des disques au fil des ans ou des décennies. Ces modèles promettent d'identifier des profils de charge « dangereux » bien avant l'apparition de symptômes cliniques, ce qui permet une intervention précoce.

Orientations et défis futurs

Malgré des progrès impressionnants, la simulation du disque épinière n'est pas encore un outil clinique courant. Plusieurs défis subsistent. Premièrement, l'obtention de propriétés matérielles spécifiques au patient reste difficile – la plupart des modèles reposent sur des données cadavériques agrégées qui ne reflètent peut-être pas l'état du tissu exact de l'individu. Deuxièmement, l'intégration de changements biologiques dépendant du temps (p. ex., déshydratation du disque avec l'âge) dans les modèles prédictifs est une frontière de recherche permanente.

Conclusion

Le comportement mécanique du disque épinière dans différentes conditions de charge est un domaine riche et cliniquement vital de la recherche biomécanique. Grâce à des simulations détaillées de FEA, fondées sur l'anatomie des disques, les propriétés matérielles et les schémas physiologiques de charge, les chercheurs ont acquis une connaissance approfondie de la compression, de la torsion, de l'extension de la flexion, du cisaillement et de la charge combinée.

Pour plus de précisions, voir Nerurkar et al. (2010) sur la mécanique annulaire du fibrosus, Vergroesen et al. (2015) sur la dégénérescence et la nutrition des disques, et Schmidt et al. (2014) pour une revue complète de la modélisation des éléments finis de la colonne vertébrale lombaire.