Analyse des méthodes par éléments finis et méthodes sans maillage pour la déformation de corps mous en simulation chirurgicale

Ce projet de recherche s'intéresse au problème de la simulation chirurgicale qui constitue un des grands défis du domaine de l'animation en temps réel de corps virtuellement déformables. Le processus de simulation doit représenter le comportement d'un organe déformable qui réagit aux manipulations de l'utilisateur et de son outil haptique, et ce en temps réel. La modélisation repose sur la résolution d'un système d'équations aux dérivées partielles complexe pour lequel la méthode des éléments finis est généralement favorisée. Cependant, cette méthode nécessite une discrétisation du modèle déformable en une suite d'éléments géométriques connectés entre eux, un processus fastidieux et mal adapté pour des simulations qui nécessitent des découpes ou autres manipulations entrainant un changement de la topologie initiale de l'objet simulé. Ce travail cherche à confirmer l'hypothèse que des méthodes de discrétisation basées sur des particules, et donc sans maillage d'éléments, sont suffisamment rapides, précises et stables pour pallier le problème de l'approche traditionnelle des éléments finis. Pour y arriver, deux méthodes de discrétisation sans maillage sont construites. Les deux utilisent une formulation faible des équations d'équilibre de la théorie de l'élasticité linéaire en mécanique des milieux continus. Le premier modèle, la méthode basée sur les points (MBP), utilise une approximation du déplacement d'une particule basée sur un voisinage compact bâti autour de celle-ci à l'aide d'une fonction noyau. Pour intégrer l'équation d'équilibre sur le domaine non déformé, le modèle utilise le volume d'une particule décrit comme le rapport entre sa masse et la densité de son voisinage. Deux fonctions d'approximation du déplacement sont formulées : la fonction de forme “Smooth Particle Hydrodynamics” (SPH) et la fonction de forme “Moving least square” (MLS). Le deuxième modèle est basé sur la méthode Meshless Total Lagrangian Dymamics (MTLED) qui utilise également une approximation du déplacement d'une particule basée sur un voisinage compact, mais cette fois en utilisant un maillage d'éléments pour l'intégration de l'équation d'équilibre sur le domaine non déformé. Un algorithme intégrant ces deux modèles est proposé. Pour valider l'approche sans maillage, un modèle de référence basé sur la méthode des éléments finis avec des éléments tétraédriques et hexaédriques linéaires est décrit et également implémenté. La performance des méthodes sans maillage et de la méthode des éléments finis est finalement comparée. La comparaison est effectuée sur différents scénarios d'étirement et de fléchissement. Une analyse de la résolution statique du système d'équilibre et de la matrice de rigidité est également abordée. Les techniques d'intégration dans le temps et de gestion des forces de contact sont aussi testées. Les résultats obtenus amènent à l'invalidation de l'hypothèse de départ. Pour arriver à atteindre du temps réel, seules les méthodes d'intégration explicites dans le temps peuvent être utilisées sur des objets de taille considérable, c'est-à-dire comportant quelques milliers de particules. D'ailleurs, même avec une méthode d'intégration explicite dans le temps, un nombre très restreint de voisins est nécessaire, entrainant ainsi une mauvaise propagation des forces. De plus, si le nombre de particules du volume est diminué en augmentant la taille du voisinage pour ainsi accélérer le temps de calcul, la faible consistance des fonctions d'approximation apporte de graves conséquences sur la précision et la stabilité de l'animation. Pour ces raisons, l'utilisation des deux approches sans maillage étudiées lors d'une modélisation complète du comportement d'un organe dans une simulation en temps réel n'apporte aucun gain par rapport à la méthode des éléments finis.