Brevet FR3009406A1 — publié le 20150206

La présente invention concerne un procédé de détermination d'un tracé de rayon d'une onde ultrasonore traversant une structure à interfaces diffractantes. Elle concerne également un programme d'ordinateur correspondant. Une application de l'invention est la simulation, par méthode algorithmique numérique, d'une inspection ultrasonore de structure réalisée par au moins un couple émetteur/récepteur d'ondes ultrasonores. Un domaine industriel particulièrement intéressé par un tel procédé est notamment le contrôle non destructif, c'est-à-dire l'inspection d'une structure sans avoir à l'altérer ou la détruire, notamment pour des applications nucléaires (l'inspection d'une cuve de réacteur), aéronautiques (le contrôle de l'intégrité de la structure d'un avion), de transport ou de métallurgie. L'objectif final est de décrire le plus complètement possible, sous forme de tracés de rayons simulés, la propagation bidimensionnelle ou tridimensionnelle d'ondes ultrasonores transitant entre un point source (i.e. l'émetteur ultrasonore) et un point d'observation, entre lesquels se trouve une structure présentant une ou plusieurs interfaces diffractantes, ces interfaces diffractantes pouvant engendrer des hétérogénéités de propagation ultrasonore et inclure des défauts de structure (par exemple des fissures, inclusions, délaminages, décollements,....). Il est notamment pertinent de prendre en compte tous les effets de diffraction d'ondes conformément, d'une part, au principe de Huygens-Fresnel selon lequel chaque point d'une interface atteint par une onde se comporte comme une source de diffraction secondaire et, d'autre part, à la Théorie Géométrique de la Diffraction reposant, dans sa forme élémentaire, sur le principe de Fermat selon lequel une onde se propage d'un point à un autre sur des trajectoires, de types tracés de rayons, telles que la durée du parcours soit localement stationnaire, c'est-à-dire en pratique minimale. La simulation d'une inspection ultrasonore de structure par contrôle non destructif peut avantageusement être combinée à une inspection expérimentale de la structure selon cette méthode afin de pouvoir interpréter correctement les résultats obtenus expérimentalement. D'une façon générale, il existe de nombreuses solutions algorithmiques permettant de modéliser par simulation la propagation d'une onde à l'aide de rayons. Ces solutions ont été développées essentiellement dans le domaine de la géophysique et des ondes sismiques, mais certaines ont pu être adaptées à l'échelle du contrôle non destructif par émission/réception d'ondes ultrasonores. Une première famille de solutions est basée sur le principe du lancer de rayons (de l'anglais « ray shooting »). Cette technique est particulièrement adaptée pour calculer un champ émis par une source dans une zone définie de l'espace. Le principe est de créer un premier rayon incident issu de la source, dont les caractéristiques dépendent de celles de la source (dimensions géométriques, fréquence de l'onde, longueur du paquet émis), puis de modifier les caractéristiques du rayon en fonction de l'environnement dans lequel il se propage, et des obstacles qu'il rencontre. Le défaut de cette famille de solutions est principalement qu'elle n'est pas adaptée pour déterminer efficacement et rapidement les rayons possibles entre deux points prédéterminés de l'espace puisque l'analyse du champ émis se fait à partir du point d'émission et de proche en proche sans limitation de la simulation à un point de réception prédéterminé. Un exemple de méthode développée selon cette première famille de solutions dans le domaine des ondes sismiques est donné dans l'article de Moser, intitulé « Shortest path calculation of seismic rays », publié dans Geophysics, vol. 56, n °1, pages 59-67, Janvier 1991. Dans ce document, la structure observée à partir d'un point d'émission d'onde sismique est discrétisée par maillage, chaque point de ce maillage se comportant alors comme un point de diffraction secondaire. La propagation de l'onde sismique est ensuite construite en respectant les principes de la Théorie Géométrique de la Diffraction par la constitution de proche en proche d'un arbre des temps de propagation minimaux. A chaque étape de la construction de l'arbre, une optimisation de chaque branche est réalisée par une implémentation particulière de l'algorithme de Dijkstra. Une fois l'arbre complètement construit, il est possible de déterminer un tracé de rayon entre le point d'émission que constitue la racine de l'arbre et un point de réception quelconque de l'arbre, en parcourant l'arbre selon un chemin unique. Mais la construction de l'arbre et longue et coûteuse en calculs. Une amélioration de cette méthode est proposée dans l'article de Zhang et al, intitulé « A global minimum traveltime ray tracing algorithm of wavefront expanding with interface points as secondary sources », publié dans Chinese Journal of Geophysics, vol. 49, n°4, pages 1046-1053, 2006. Dais ce document, seuls sont conservés comme points de diffraction secondaire ceux qui se situent aux interfaces entre milieux dont les modes de propagation diffèrent. Cela donne des tracés de rayons lissés dans un arbre des temps de propagation minimaux simplifié. Un autre exemple de méthode améliorée basée sur le principe enseigné par Moser est donné dans l'article de Zhao et al, intitulé « Minimum travel time tree algorithm for seismic ray tracing: improvement in efficiency », publié dans Journal of Geophysics and engineering, vol. 1, pages 245-251, septembre 2004. Dans ce document, il est prévu que la structure observée puisse être constituée de milieux hétérogènes dont les modes volumiques de propagation diffèrent, certains modes permettant une vitesse continument variable selon différents axes d'un repère cartésien. Les rayons continument variables apportent un véritable avantage dans la modélisation des milieux hétérogènes, mais là encore, la construction d'un arbre des temps de propagation minimaux est nécessaire, avec optimisation de chaque branche par application d'une implémentation particulière de l'algorithme de Dijkstra, ce qui est long et fastidieux. Un autre exemple de méthode développée selon cette première famille de solutions dans le domaine du contrôle non destructif par ultrasons est celle décrite dans l'article de Gengembre et al, intitulé « Pencil method in elastodynamics: application to ultrasonic field computation », publié dans Ultrasonics, vol. 38, pages 495-499, novembre 2000. Une onde y est représentée sous la forme de tubes de rayons dont la section variable évolue dans le temps afin de respecter le principe de conservation de l'énergie au cours de la propagation de l'onde. La propagation d'une onde est donc modélisée en calculant la trajectoire et la déformation des tubes de rayons au fur et à mesure de la propagation. Mais comme dans les deux exemples précédents, la méthode n'est pas adaptée pour déterminer efficacement et rapidement les rayons possibles entre deux points prédéterminés de l'espace. En outre, cette méthode fondée sur l'acoustique géométrique ne prend en compte que les phénomènes de réflexions/transmissions aux interfaces mais pas la diffraction engendrée par des irrégularités des interfaces ou par des défauts de structure. Une deuxième famille de solutions est basée sur le principe du tracé de rayon entre deux points prédéterminés (de l'anglais « two points ray tracing »). Cette approche algorithmique n'a pas pour but d'effectuer un calcul de champ, mais de connaître, à partir d'une onde émise en un point de l'espace, ce qui est reçu en un autre point de l'espace. Cette approche est particulièrement bien adaptée à la problématique d'un calcul d'écho, qui consiste à obtenir le signal effectivement reçu sur un capteur récepteur (un ensemble de points de « réception »), le signal pouvant être le résultat de la sommation de plusieurs ondes, donc de plusieurs rayons. L'invention s'applique plus particulièrement à cette deuxième famille de solutions en portant sur un procédé de détermination d'un tracé de rayon d'une onde ultrasonore entre un point d'émission et un point de réception par simulation d'une inspection ultrasonore de structure à au moins une interface diffractante disposée entre ces points d'émission et de réception, comportant les étapes suivantes : discrétisation de la structure en une pluralité de points de diffraction secondaire, et sauvegarde de ces points de diffraction secondaire en mémoire, construction, par un processeur ayant accès à la mémoire en lecture/écriture, d'un graphe orienté, dont les noeuds comportent le point d'émission, le point de réception et l'ensemble des points de diffraction secondaire, et dont les liens orientés sont associés à des poids représentant des temps de propagation entre les noeuds qu'ils relient, et sélection par le processeur, parmi l'ensemble des tracés de rayons possibles parcourant le graphe orienté du point d'émission au point de réception, du tracé de rayon minimisant la somme des poids des liens parcourus. Un tel procédé permet de tracer des rayons entre des points d'émission et de réception connus par avance et procède par optimisation d'un chemin parcouru dans un graphe orienté entre ces points d'émission et de réception connus par avance, cette optimisation consistant en une minimisation de temps de propagation de l'onde. Il se distingue en cela des procédés de lancer de rayons précédemment évoqués. En effet, alors que certains des procédés selon la première famille de solutions procèdent par optimisation de proche en proche de la construction d'un arbre des temps de propagation minimaux qui est ensuite parcouru de façon déterministe pour le lancer de rayons, cette deuxième famille de solution procède à l'inverse d'abord par la construction d'un graphe orienté non optimisé, puis par l'optimisation d'un parcours dans ce graphe orienté. Un exemple de méthode développée selon cette deuxième famille de solutions dans le domaine des ondes sismiques est évoqué dans la demande de brevet EP 0 541 451 A1, dans lequel l'optimisation du parcours dans un graphe est réalisée par l'algorithme de Dijkstra. Mais dans ce document, cette optimisation est exécutée pour rechercher un rayon optimal dans le voisinage d'un rayon existant, sans que la détermination du rayon existant ne soit clairement détaillée. Notamment, la gestion d'une structure à zones hétérogènes de propagation n'est pas abordée. En outre, la méthode détaillée dans ce document est adaptée au domaine des ondes sismiques mais pas au domaine des inspections par ultrasons. Il peut ainsi être souhaité de prévoir un procédé de détermination d'un tracé de rayon qui permette de s'affranchir d'au moins une partie des problèmes et contraintes précités. Il est donc proposé un procédé de détermination d'un tracé de rayon d'une onde ultrasonore entre un point d'émission et un point de réception par simulation d'une inspection ultrasonore de structure à au moins une interface diffractante disposée entre ces points d'émission et de réception, comportant les étapes suivantes : discrétisation de la structure en une pluralité de points de diffraction secondaire, et sauvegarde de ces points de diffraction secondaire en mémoire, construction, par un processeur ayant accès à la mémoire en lecture/écriture, d'un graphe orienté dont les noeuds comportent le point d'émission, le point de réception et l'ensemble des points de diffraction secondaire, et dont les liens orientés sont associés à des poids représentant des temps de propagation entre les noeuds qu'ils relient, sélection par le processeur, parmi l'ensemble des tracés de rayons possibles parcourant le graphe orienté du point d'émission au point de réception, du tracé de rayon minimisant la somme des poids des liens parcourus, selon lequel : plusieurs valeurs de vitesses de propagation associées à plusieurs modes de propagation incluant au moins un mode de propagation de surface le long de l'interface diffractante sont conservées en mémoire, et lors de la construction du graphe orienté, chaque lien orienté retenu est associé par le processeur à une vitesse de propagation parmi celles conservées en mémoire en fonction de la localisation de ce lien par rapport à l'interface diffractante, pour le calcul de son poids. Ainsi, les hétérogénéités de propagation engendrées par la présence d'une ou plusieurs interfaces diffractantes sont prises en compte par l'association à chaque lien du graphe construit d'une vitesse de propagation caractéristique d'un mode de propagation ultrasonore particulier. En prévoyant de prendre en compte non seulement les modes de propagations volumiques mais également les modes de propagations surfaciques, une meilleure simulation des effets des interfaces diffractantes est en outre assurée. De façon optionnelle, l'ensemble des points de diffraction secondaire conservés en mémoire est constitué de points de ladite au moins une interface diffractante et le processeur est conçu pour construire le graphe en excluant tout lien orienté traversant cette interface diffractante. De façon optionnelle également, ladite au moins une interface diffractante comporte au moins l'un des éléments de l'ensemble constitué d'une surface de pièce disposée dans la structure observée entre le point d'émission et le point de réception, d'une interface entre zones hétérogènes d'une même pièce et d'un défaut de structure. De façon optionnelle également, lorsque ladite au moins une interface diffractante comporte un défaut de structure modélisé par au moins une facette polygonale plane présentant des extrémités, les seuls points de diffraction secondaire retenus sur ce défaut de structure sont situés aux extrémités de chaque facette le modélisant. De façon optionnelle également : - chaque zone de propagation homogène considérée comme isotrope dans la structure est associée à au moins un mode de propagation lui-même associé à une vitesse de propagation caractéristique de ce mode de propagation, et - chaque zone de propagation homogène considérée comme anisotrope dans la structure est associée à au moins un mode de propagation lui-même associé à une pluralité de vitesses de propagation en fonction d'une direction et d'un sens de propagation, cette pluralité de vitesses de propagation étant conservée en mémoire sous la forme d'une surface des lenteurs. De façon optionnelle également, la construction du graphe orienté comporte les étapes suivantes : - construction d'un graphe orienté intermédiaire dont les noeuds comportent uniquement les points de diffraction secondaire et dont les liens orientés sont associés à des poids représentant des distances ou temps de propagation entre les noeuds qu'ils relient, puis - construction du graphe orienté en complétant le graphe orienté intermédiaire par l'ajout d'une pluralité de points de passage ou groupes de points de passage imposés incluant le point d'émission et le point de réception. Ces points de passage ou groupes de points de passage imposés, lorsqu'ils sont autres que le point d'émission et le point de réception, correspondent à des contraintes d'interaction entre le tracé de rayon et les surfaces diffractantes ou les défauts d'une pièce d'étude. De façon optionnelle également, la construction du graphe orienté intermédiaire comporte les étapes suivantes : classement des points de diffraction secondaire dans l'ordre croissant de leur éloignement au point d'émission, calcul des coefficients d'une matrice M modélisant ce graphe orienté intermédiaire, chaque coefficient de cette matrice M portant la valeur d'une distance ou d'un temps de propagation d'onde ultrasonore entre les ième et j-ème points de diffraction secondaire dans l'ordre du classement, et dans lequel les règles suivantes sont appliquées : tout coefficient pour lequel j La configuration générale d'une installation 10 d'inspection ultrasonore de structure 12 à au moins une interface diffractante est illustrée sur la figure 1. Cette installation 10, apte à mettre en oeuvre un procédé de détermination d'un tracé de rayon d'une onde ultrasonore entre des points d'émission et de réception de la structure 12, comporte les éléments suivants : un capteur plan 14 émetteur d'ondes ultrasonores, par exemple un capteur plan de contact émettant des ondes ultrasonores longitudinales de fréquence égale à environ 5 MHz, un sabot d'émission 16, par exemple en polyméthacrylate de méthyle ou PMMA (de l'anglais « PolyMethyle MethAcrylate), sur lequel est disposé le capteur plan émetteur 14, une pièce d'étude 18, par exemple formée d'un bloc parrallélépipédique rectangle d'acier austénitique de géométrie d'extrusion bidimensionnelle comportant un affouillement 20 aux bordures arrondies, sur laquelle est disposé le sabot d'émission 16, un sabot de réception 22, par exemple en PPMA, par exemple identique au sabot d'émission 16, un capteur plan 24 récepteur d'ondes ultrasonores, disposé sur le sabot de réception 22, par exemple un capteur plan de contact apte à recevoir des ondes ultrasonores longitudinales de fréquence égale à environ 5 MHz, et un système informatique 26 comportant une unité centrale 28 et des périphériques de saisie/affichage. L'unité centrale 28 du système informatique 26 est reliée aux capteurs plans émetteur 14 et récepteur 24 pour engendrer/émettre des signaux ultrasonores à l'aide du capteur plan émetteur 14 et recevoir/traiter des signaux ultrasonores à l'aide du capteur plan récepteur 24. Elle est en outre programmée pour la mise en oeuvre d'un procédé de détermination d'un tracé de rayon d'une onde ultrasonore selon l'invention. Les sabots d'émission et réception 16, 22 sont conformés pour imposer un angle optimal d'émission/réception des ondes ultrasonores par exemple longitudinales transitant dans la pièce d'étude 18, par exemple un angle proche de 60°. Ils orientent pour cela les capteurs plans par rapport à la surface de la pièce d'étude 18 selon un angle prédéterminé. Une onde ultrasonore émise par le capteur plan émetteur 14 est ainsi orientée de façon optimale par réfraction entre le sabot d'émission 16 en PMMA et la pièce d'étude 18 en acier austénitique. De même, une onde ultrasonore est reçue selon une orientation optimale par le capteur plan récepteur 24 après réfraction entre le sabot d'émission 16 en PMMA et la pièce d'étude 18 en acier austénitique. L'inspection de la structure 12 se fait alors selon une technique de mesure de temps de vol d'onde diffractée dite TOFD (de l'anglais « Time Of Flight Diffraction »), selon laquelle deux transducteurs (i.e. les capteurs plans émetteur 14 et récepteur 24) placés face à face à une distance prédéterminée se déplacent le long de la surface d'une pièce inspectée, l'un émettant des impulsions ultrasonores, l'autre recevant les ondes qui se sont alors propagées dans la pièce. La structure 12, inspectée à l'aide du système informatique 26 et des capteurs plans émetteur 14 et récepteur 24, comporte l'ensemble constitué de la pièce d'étude 18 et d'au moins une partie des sabots d'émission et réception 16, 22. Elle comporte donc au moins une interface diffractante, à la surface de la pièce d'étude 18, entre cette pièce et chaque sabot d'émission et réception 16, 22, disposée nécessairement sur le trajet des ondes ultrasonores émises et reçues par les capteurs plans émetteur 14 et récepteur 24. La pièce d'étude 18 présente en outre, dans l'exemple illustré par la figure 1, deux défauts de structure D1 et D2 se présentant sous la forme de segments verticaux situés sous l'affouillement 20 dans la coupe plane de la figure 1. Ces défauts de structures D1 et D2 constituent également des interfaces diffractantes susceptibles de se trouver sur le trajet d'ondes ultrasonores émises ou reçues par les capteurs plans émetteur 14 et récepteur 24. Enfin, la pièce d'étude 18 pourrait être composée de plusieurs zones hétérogènes à modes de propagation d'ondes différents, leurs interfaces se comportant alors également comme des interfaces diffractantes susceptibles de se trouver sur le trajet d'ondes ultrasonores émises ou reçues par les capteurs plans émetteur 14 et récepteur 24. Un procédé selon l'invention a pour but de déterminer de façon efficace et précise un tracé T de rayon d'onde ultrasonore entre un point d'émission imposé PE situé dans la structure 12, par exemple plus précisément dans le sabot d'émission 16 à proximité du capteur plan émetteur 14, et un point de réception imposé PR situé également dans la structure 12, par exemple dans un plan de coupe vertical de la pièce d'étude 18 dont une section S est représentée sur la figure 1, ce plan de coupe étant situé en aval des défauts de structure D1 et D2 dans le sens de transmission des ondes émises à partir du point d'émission PE. Pour la mise en oeuvre d'un tel procédé, l'unité centrale 28 du système informatique 26, telle que représentée schématiquement sur la figure 2, peut par exemple être constituée d'un dispositif informatique tel qu'un ordinateur classique comportant un processeur 30 associé à une ou plusieurs mémoires ou zones de mémoire 32A, 32B, 32C (par exemple une mémoire RAM fonctionnellement divisée en plusieurs zones) pour le stockage de fichiers de données et d'instructions que le processeur 30 exécute sous forme de programmes d'ordinateurs. Par exemple, une première zone de mémoire 32A est affectée au stockage de données d'entrée utilisées par le processeur 30. Une deuxième zone de mémoire 32B est affectée au stockage de programmes d'ordinateurs exécutables par le processeur 30. Et une troisième zone de mémoire 32C est affectée au stockage de données de sortie fournies par le processeur 30. L'unité centrale 28 telle qu'illustrée sur la figure 2 comporte ainsi fonctionnellement deux programmes d'ordinateurs ou deux fonctions d'un même programme d'ordinateur 34 et 36 stocké(e)s en mémoire 32B. On notera en effet que les programmes d'ordinateurs 34 et 36 sont présentés comme distincts, mais cette distinction est purement fonctionnelle. Ils pourraient tout aussi bien être regroupés selon toutes les combinaisons possibles en un ou plusieurs logiciels. Leurs fonctions pourraient aussi être au moins en partie micro programmées ou micro câblées dans des circuits intégrés dédiés. Ainsi, en variante, le dispositif informatique mettant en oeuvre l'unité centrale 28 pourrait être remplacé par un dispositif électronique composé uniquement de circuits numériques (sans programme d'ordinateur) pour la réalisation des mêmes actions. Les données d'entrée stockées dans la première zone de mémoire 32A comportent : des données 38 relatives à la géométrie de la pièce d'étude 18 sans tenir compte de ses éventuels défauts, des données 40 relatives à la géométrie des défauts, ici D1 et D2, de la pièce d'étude 18, des données 42 relatives à des paramètres physico-techniques liés aux modes de propagation possibles dans la structure 12 et à un pas de discrétisation souhaité pour la mise en oeuvre d'un procédé de détermination de tracé de rayon selon l'invention, des données 44 relatives à des points de passage ou groupes de points de passage imposés pour un tracé de rayon, ces points de passage ou groupes de points de passage incluant au moins le point d'émission PE et le point de réception PR. Les données 38 comportent par exemple une description des surfaces délimitant la pièce d'étude 18, ainsi que des interfaces internes dans le cas d'une pièce à plusieurs zones hétérogènes. Cela peut prendre la forme d'un fichier de CAO (« conception Assistée par Ordinateur »), par exemple au format DXF (de l'anglais « Drawing eXchange Format ») compatible avec le logiciel AutoCAD (marque déposée). Dans le cas simple de la pièce d'étude 18 à géométrie bidimensionnelle de la figure 1, la pièce est complètement définie dans un repère orthonormé cartésien bidimensionnel (x, y) ayant par exemple pour unité le millimètre. Le fichier de CAO décrit alors les interfaces de la pièce d'étude 18 sous la forme d'un ensemble de segments et d'arcs de cercle, l'affouillement 20 creusé dans la forme rectangulaire de la pièce 18 étant en particulier représenté par deux quarts de cercles reliés par un segment. L'exemple donné ici en deux dimensions est très facilement généralisable à une description tridimensionnelle. Les données 40 comportent une description spécifique des défauts de la pièce d'étude 18. Ces données 40 sont distinguées des données 38 parce que, dans une variante de réalisation optionnelle mais préférée de l'invention qui sera détaillée ci-dessous avec les données 42, elles ne sont pas traitées de la même façon. Mais elles peuvent également prendre la forme d'un fichier de CAO similaire au précédent. D'une façon générale, chaque défaut de structure peut être représenté sous la forme d'une facette polygonale plane ou de plusieurs facettes polygonales planes reliées entre elles. Pour un défaut, ces facettes planes décrivent sa surface, ouverte ou fermée. Mises bout-à-bout, elles permettent de reproduire un défaut de surface complexe. En deux dimensions, ces facettes planes sont représentées par des segments, disposés bout-à-bout lorsqu'il y en a plusieurs pour un même défaut. Dans l'exemple simple de la figure 1 où les défauts D1 et D2 sont plans et verticaux, chaque défaut D1 ou D2 est ainsi constitué d'une unique facette et représenté par un segment vertical à deux extrémités. Les ondes ultrasonores peuvent se propager sur les facettes ou être diffractées sur les bords des facettes. Les données 38 et 40 constituent une description complète des interfaces diffractantes de la structure 12. Les données 42 comportent tout d'abord une description des différents modes de propagation, volumiques et surfaciques, possibles dans la structure 12 à l'aide de paramètres physiques. Au moins un mode de propagation est défini pour chaque milieu et chaque interface diffractante présents dans la structure 12. En particulier, dans l'exemple de la figure 1, un premier mode de propagation volumique peut être défini pour le milieu en PMMA des sabots d'émission et réception 16, 22. Un deuxième mode de propagation volumique peut être défini pour le milieu en acier austénitique de la pièce d'étude 18 sachant qu'elle ne présente pas de zones hétérogènes. Un premier mode de propagation surfacique peut être défini pour les surfaces délimitant la pièce d'étude 18. Enfin, conformément à la variante de réalisation préférée mentionnée précédemment, selon laquelle les données 40 ne sont pas traitées de la même façon que les données 38, un deuxième mode de propagation surfacique, indépendant du précédent, peut être défini pour les défauts de structures. Si l'on considère que les milieux de propagation sont isotropes, chaque mode de propagation peut être défini par une unique vitesse de propagation qui le caractérise. Ainsi, dans l'exemple de la figure 1, un quadruplet de vitesses de propagation définit à lui seul les quatre modes de propagation volumique et surfacique possibles. Si l'on considère au contraire que les milieux de propagation sont anisotropes, chaque mode de propagation peut être défini par une unique surface des lenteurs qui le caractérise, cette surface des lenteurs définissant une vitesse de propagation par direction et sens de propagation dans le milieu considéré. Ainsi, dans l'exemple de la figure 1, un quadruplet de surfaces des lenteurs définit à lui seul les quatre modes de propagation volumique et surfacique possibles. Grâce aux paramètres physiques précités, il est possible de prendre en compte des conversions de modes de propagation au cours d'un tracé de rayon, comme par exemple le passage d'une onde de volume transversale ou longitudinale en une onde de surface de Rayleigh et réciproquement. Les données 42 comportent en outre une définition d'un pas de discrétisation souhaité des interfaces diffractantes définies par les données 38 et 40, ce pas de discrétisation constituant un paramètre technique pour la mise en oeuvre d'un procédé de détermination de tracé de rayon selon l'invention. Il représente la distance entre deux points de discrétisation consécutifs d'une interface diffractante, par exemple choisie à 0,5 mm. Conformément au fonctionnement d'un tel procédé qui sera détaillé en référence à la figure 3, plus le pas de discrétisation est fin, plus la détermination du tracé de rayon est précise. Conformément à la variante de réalisation préférée mentionnée précédemment, selon laquelle les données 40 ne sont pas traitées de la même façon que les données 38, le pas de discrétisation ne s'applique qu'aux interfaces définies par les données 38, les seuls points de discrétisation retenus sur les défauts de structure étant les extrémités des facettes qui les modélisent. La taille des facettes de chaque défaut peut en revanche être définie librement et indépendamment du pas de discrétisation s'appliquant aux interfaces définies par les données 38. Ainsi, pour les défauts D1 et D2 représentés par deux segments verticaux sur la figure 1, seules les extrémités des deux segments sont retenues comme points de discrétisation, indépendamment du pas de discrétisation défini. En effet, il est considéré qu'aucun tracé de rayon ne peut traverser les facettes d'un défaut de structure interne à la pièce d'étude 18 alors qu'il peut traverser les surfaces qui la délimitent selon le pas de discrétisation défini pour ces surfaces. Un tracé de rayon peut en revanche longer une facette de défaut. Par exemple, une onde ultrasonore peut tout à fait se propager autour des facettes d'un défaut de type inclusion formant une structure fermée. Les données 44 fournissent les coordonnées des points de passage ou groupes de points de passage imposés pour un tracé de rayon. Comme indiqué précédemment, ces points de passage ou groupes de points de passage correspondent pour la plupart à des contraintes d'interaction entre le tracé de rayon et les surfaces diffractantes ou les défauts de la pièce d'étude 18. Ils sont classés dans l'ordre de passage, du premier (le point d'émission PE) au dernier (le point de réception PR). Ils sont tous situés dans la structure 12, mais indifféremment à l'intérieur, en surface ou à l'extérieur de la pièce d'étude 18. Les données 44 définissent donc au minimum deux points de passage, PE et PR, mais peuvent en définir d'autres si l'on souhaite déterminer des tracés de rayons en imposant des contraintes précises de passages intermédiaires entre ces deux points. Par ailleurs, il peut être prévu de définir un quadruplet de vitesses de propagation ou de surfaces des lenteurs par paire de points de passage ou groupes de points de passage successifs, comme cela sera détaillé plus loin. Dans ce cas, il faut en tenir compte dans la description des différents modes de propagation par les données 42. Plusieurs modes de propagation peuvent être prévus pour chaque milieu et chaque interface diffractante de la structure 12 en fonction de la paire de points de passage ou groupes de points de passage successifs concernée. A partir des données d'entrée 38, 40 et 42 précitées, le premier programme d'ordinateur 34 stocké dans la deuxième zone de mémoire 32B est conçu pour discrétiser les interfaces diffractantes de la structure 12 et pour construire un graphe orienté intermédiaire dont les noeuds sont les points de discrétisation des surfaces diffractantes, ces points étant alors considérés comme constituant l'ensemble des points de diffraction secondaire de la structure 12. Les liens orientés de ce graphe intermédiaire, construit selon une méthode qui sera détaillée en référence à la figure 3 pour un premier mode de réalisation ou selon une méthode qui sera détaillée en référence à la figure 9 pour un deuxième mode de réalisation, sont associés à des poids qui, selon le mode de réalisation choisi, représentent des temps de propagation ou des distances entre les noeuds qu'ils relient. Le graphe orienté intermédiaire est sauvegardé sous la forme d'un fichier 46 dans la troisième zone de mémoire 32C. A partir des données d'entrée 44 et du graphe orienté intermédiaire sauvegardé dans le fichier 46, le deuxième programme d'ordinateur 36 stocké dans la deuxième zone de mémoire 32B est conçu pour compléter le graphe orienté intermédiaire à l'aide des points de passage ou groupes de points de passage imposés de manière à construire un graphe orienté final et pour fournir une détermination de tracé de rayon par optimisation d'un parcours dans le graphe orienté final. Le graphe orienté final est constitué de noeuds comportant le point d'émission PE, le point de réception PR, les autres points de passage ou groupes de points de passage imposés et l'ensemble des points de diffraction secondaire. Ses liens sont associés à des poids représentant des temps de propagation entre les noeuds qu'ils relient. Il est sauvegardé, par mise à jour du graphe orienté intermédiaire, dans le fichier 46. La détermination de tracé de rayon peut comporter plusieurs données de sortie. Des données de sortie 48, sauvegardées dans la troisième zone de mémoire 32C, définissent géométriquement le tracé de rayon déterminé par optimisation de parcours dans le graphe orienté final. Des données de sortie 50, sauvegardées dans la troisième zone de mémoire 32C, définissent le temps de vol d'une onde ultrasonore associée à ce tracé de rayon. Enfin, des données de sortie 52, sauvegardées dans la troisième zone de mémoire 32C, définissent la succession des modes de propagation le long de ce tracé de rayon. On note donc que selon l'invention, il est supposé qu'un tracé de rayon est constitué d'une série de trajets élémentaires rectilignes entre une pluralité de points comprenant le point d'émission, le point de réception, les autres points de passage imposés et une partie des points de diffraction secondaire dus aux interfaces diffractantes de la structure 12. Cette hypothèse est réaliste et permet de fournir une bonne simulation par optimisation de parcours dans un graphe orienté. Conformément à un premier mode de réalisation de l'invention, le procédé de la figure 3, mis en oeuvre par exécution des programmes d'ordinateurs 34 et 36 à l'aide du processeur 30, comporte une première étape 100 d'acquisition des données d'entrée 38, 40 et 42. Cette première étape d'acquisition 100 est réalisée par exécution du premier programme d'ordinateur 34. Au cours d'une étape de discrétisation 102 suivante, à l'aide des données 38 relatives aux interfaces de diffraction secondaire de la structure 12 n'incluant pas les défauts de la pièce d'étude 18 et à l'aide du pas de discrétisation défini dans les données d'entrée 42, des premiers points de diffraction secondaire à retenir pour la construction du graphe orienté intermédiaire précité sont définis. Cette discrétisation permet de s'affranchir d'une définition analytique de la pièce d'étude 18 en termes de surfaces externes et d'interfaces internes entre zones hétérogènes éventuelles. Il est ainsi possible de discrétiser une pièce disposant de surfaces et interfaces irrégulières quelconques. Au cours d'une autre étape de discrétisation 104, à l'aide des données 40 relatives aux défauts de structure D1 et D2 et conformément à la variante de réalisation préférée mentionnée précédemment, selon laquelle les données 40 ne sont pas traitées de la même façon que les données 38, seules les extrémités des facettes de défauts de structure, ici les quatre extrémités des deux segments modélisant les défauts D1 et D2, sont ajoutées en tant que points de diffraction secondaire supplémentaires à retenir pour la construction du graphe orienté intermédiaire. A l'issue des deux étapes 102 et 104, réalisées par exécution du premier programme d'ordinateur 34, N points de diffraction secondaire sont ainsi définis et retenus pour la construction du graphe orienté intermédiaire. Ils sont classés dans l'ordre croissant de leur éloignement au point d'émission PE et stockés en mémoire 32C. Ces N points sont illustrés sur la figure 4, ceux correspondant à la surface externe de la pièce d'étude 18 étant représentés par des croix obliques « x » et ceux correspondant aux défauts D1 et D2 étant représentés par des croix droites « + ». Les interfaces internes que constituent les défauts sont quant à elles représentées en traits interrompus mixtes. Au cours d'une étape de construction 106 suivante, réalisée par exécution du premier programme d'ordinateur 34, le graphe orienté intermédiaire est construit par calcul d'une matrice M de taille NxN le modélisant. Un coefficient quelconque de cette matrice a pour valeur un temps de vol, ou temps de propagation d'onde ultrasonore, entre le i-ème point P, et le j-ème point P, de diffraction secondaire. Ce temps de vol est aisément calculable à partir de la distance P,P, entre les deux points considérés et de la vitesse de propagation de l'onde ultrasonore dans le milieu du volume ou de l'interface diffractante (surface ou défaut de la pièce d'étude 18) où se trouve le trajet élémentaire [P,, PJ] et dans la direction de ce trajet élémentaire. La vitesse de propagation à appliquer est déterminée à l'aide des paramètres physiques définis dans les données d'entrée 42. Il convient en outre d'appliquer un certain nombre de règles définissant certains trajets élémentaires comme impossibles donc de temps de vol infini noté 09. Etant donné qu'un temps de vol est en réalité nécessairement positif ou nul, une valeur négative, par exemple « -1 », peut être affectée aux temps de vol infinis 09. Ces règles permettent de faire le lien entre les contraintes physiques portant sur un tracé de rayon et la résolution théorique du problème mathématique d'optimisation de parcours dans un graphe orienté. Elles permettent en outre d'accélérer le calcul de la matrice M. Ces règles comportent par exemple : une première règle selon laquelle une onde ultrasonore émise depuis le point d'émission PE ne peut pas se rétro-propager, une deuxième règle selon laquelle un trajet élémentaire traversant une interface diffractante n'est pas valide parce qu'il ne respecte pas le principe de Huygens selon lequel une interface est interprétée comme un ensemble de sources secondaires, toute onde traversant cette interface devant alors être considérée comme suivant un premier trajet élémentaire aboutissant en un point de diffraction secondaire de cette interface puis un deuxième trajet élémentaire débutant en ce même point, - une troisième règle selon laquelle un trajet élémentaire traversant un défaut de structure n'est pas valide, pour les mêmes raisons que celles de la deuxième règle. La première règle a pour conséquence que v (i, j), j La deuxième règle a pour conséquence que pour tout couple de points de diffraction (P,, Pi), si le trajet élémentaire [P,, Pi] traverse une interface diffractante, par exemple la surface externe de la pièce d'étude 18 ou toute interface interne entre zones hétérogènes de cette pièce, alors mj,i = 09 (ou -1). Une infraction à cette deuxième règle est illustrée sur la figure 5 par des trajets élémentaires référencés t2 relatifs à des traversées de surface externe de la pièce d'étude 18. Ces trajets élémentaires t2 sont donc associés à un temps de vol infini. La troisième règle est par défaut traitée de la même manière que la deuxième, sauf que l'interface diffractante traversée considérée est celle constituée par une facette d'un défaut. Une infraction à cette troisième règle est illustrée sur la figure 5 par un trajet élémentaire référencé t3 relatif à une traversée du défaut D1. Ce trajet élémentaire t3 est donc associé à un temps de vol infini. Le calcul de la matrice M à l'étape 106 se fait donc de la façon suivante : 1. v (i, j), j On notera qu'il existe toujours une unique vitesse de propagation à associer au trajet élémentaire [P,, PJ] dans le cas 3.b de l'étape 106 puisque, aucune intersection n'étant détectée entre le trajet élémentaire [Pi, Pj] et l'une quelconque des interfaces diffractantes, celui-ci est nécessairement intégralement situé dans une même zone homogène de la structure 12, qu'elle soit volumique ou surfacique, isotrope ou anisotrope. Dans l'exemple de la figure 1, il s'agit donc juste de localiser le trajet élémentaire [Pi, Pj] à l'intérieur, à l'extérieur, sur une facette d'un défaut ou sur la surface externe de la pièce d'étude 18 pour connaître le mode de propagation associé puis déterminer la vitesse de propagation correspondante grâce aux paramètres physiques des données d'entrée 42. A titre d'exemple, les trajets élémentaires illustrés sur la figure 5 en traits pleins sont des trajets pour lesquels il est possible d'associer un temps de vol non infini. Ceux référencés t4 sont localisés dans le volume de la pièce d'étude 18 et celui référencé t5 est localisé sur sa surface externe. Au cours d'une étape 108, réalisée par exécution du premier programme d'ordinateur 34, une matrice V de même taille NxN que la matrice M est calculée. Chaque coefficient v,,, de cette matrice V définit la nature du trajet élémentaire [Pi, Pj] lorsque ce trajet élémentaire est associé à un temps de vol non infini. Plus précisément, l'un quelconque de ces coefficients porte une ou plusieurs valeurs définissant la localisation du trajet élémentaire (à l'intérieur, sur une facette de défaut, en surface ou à l'extérieur de la pièce), la nature de la propagation associée (volumique ou surfacique) et le mode de propagation associé (longitudinal, transversal, quasi longitudinal, quasi transversal, Rayleigh, ...). Cette matrice V, qui permet alors de donner connaissance de la nature d'une onde en tout point d'un parcours constitué de trajets élémentaires, est utilisée le cas échéant pour l'application d'un modèle d'amplitude d'onde ultrasonore le long de ce parcours. On notera que l'étape 108 peut être indifféremment exécutée avant, pendant ou après l'étape 106. Après exécution des étapes 100 à 108, le premier programme d'ordinateur 34 fournit les matrices calculées M et V et les enregistre dans le fichier 46. Au cours d'une étape de mise à jour 110, réalisée par exécution du deuxième programme d'ordinateur 36, L points de passage définis dans les données d'entrée 44, parmi lesquels le point d'émission PE et le point de réception PR, sont ajoutés au graphe orienté intermédiaire modélisé par la matrice M s'ils n'y sont pas déjà présents. En effet, les points de passage, pour qu'ils aient un sens physique, doivent pour la plupart appartenir à une surface ou un défaut de la pièce d'étude 18 puisque ce sont les seuls éléments diffractants de la structure 12 et donc les seuls susceptibles de créer un rayon élémentaire. Il est donc probable que ces points de passage fassent déjà partie de la matrice M, mais pas nécessairement. Par exemple, pour une plus grande liberté d'implémentation d'un procédé selon l'invention et une meilleure précision de la simulation, la position des points de passage peut ne pas correspondre aux points de discrétisation des interfaces diffractantes pris en considération dans la matrice M. Il est effectivement important que la discrétisation des contraintes de passage soit indépendante de la discrétisation des interfaces diffractantes pour une meilleure flexibilité. Au cours de l'étape 110, il est donc nécessaire de déterminer quels sont les K points de passage, parmi les L définis dans les données d'entrée 44, qui doivent être ajoutés à la matrice M, avec 0 « K « L. On calcule ainsi une matrice Mf de taille (N+K)x(N+K) dans laquelle les points de passage non encore présents dans la matrice M sont insérés en fonction de leur éloignement par rapport au point d'émission PE. K lignes et K colonnes sont ainsi ajoutées à la matrice M pour former la matrice Mf. Les nouveaux coefficients ajoutés sont calculés selon les mêmes règles que celles appliquées à l'étape 106. Cela revient à ajouter de nouveaux trajets élémentaires valides ayant comme origine ou destination l'un des nouveaux points de passage ajoutés. La matrice Mf modélise le graphe orienté final et est stockée dans le fichier 46. Au cours de cette même étape 110, on calcule également une matrice Vf de taille (N+K)x(N+K) obtenue à partir de la matrice V. De même, les points de passage non encore présents dans la matrice V sont insérés dans la matrice Vf en fonction de leur éloignement par rapport au point d'émission PE. K lignes et K colonnes sont ainsi ajoutées à la matrice V pour former la matrice Vf dont les nouveaux coefficients ajoutés sont calculés comme à l'étape 108 et qui est elle aussi stockée dans le fichier 46. L'ajout a posteriori des points de passage à l'étape 110 dans les matrices M et V déjà calculées aux étapes 106 et 108 produit un gain considérable en temps de calcul. En effet, si plusieurs tracés de rayons doivent être déterminés dans une même configuration de structure (pièce d'étude, défauts et nature de la propagation dans chaque zone homogène de propagation), seuls les points de passage imposés changent d'un tracé à l'autre. Les calculs faits aux étapes 106 et 108 n'ont ainsi pas besoin d'être répétés puisqu'ils sont indépendants de ces points de passage. L'étape 110 est suivie d'une étape 112 de parcours optimisé du graphe orienté final modélisé par la matrice Mf. Ce parcours optimisé est un parcours de moindre coût entre le point d'émission PE et le point de réception PR, en passant par tous les autres points de passage imposés, les coûts étant représentés par les coefficients de la matrice Mf, c'est-à-dire les poids des liens orientés du graphe qu'elle modélise. C'est un problème bien connu de la théorie des graphes. Conformément à ce premier mode de réalisation, l'optimisation de parcours est répétée indépendamment L-1 fois entre chaque paire de points de passage successifs par application d'une méthode d'optimisation de parcours de graphe orienté telle que l'algorithme de Dijkstra. Cet algorithme bien connu ne sera pas détaillé. Il permet de sélectionner le tracé de rayon le plus court en temps de vol entre chaque point de passage, c'est-à-dire le tracé de rayon réel respectant le principe de Fermat généralisé. Il est en outre compatible avec le choix de modéliser les temps de vol infinis par une valeur négative telle que « -1 » parce qu'il est conçu pour refuser les trajets élémentaires à temps de vol négatif. Au cours d'une dernière étape 114, le résultat de l'optimisation 112 est traité pour fournir les données de sortie 48, 50 et 52. Tout d'abord, les L-1 parcours optimisés à l'étape 112 sont concaténés pour fournir une définition géométrique du tracé de rayon complet déterminé par l'algorithme de Dijkstra. Cette définition géométrique du tracé est paramétrée à l'aide des données de sortie 48. Le temps de vol ultrasonore associé, défini par les données de sortie 50, est obtenu en sommant les L-1 temps de vol optimisés à l'étape 112. Enfin, grâce à la matrice Vf, la succession des natures de trajets élémentaires, sélectionnés à l'étape 112 pour former le tracé de rayon complet, est enregistrée, par exemple sous forme d'une liste de valeurs formant les données de sortie 52. La connaissance de la succession des natures de trajets élémentaires permet d'appliquer un modèle d'amplitude et de déphasage d'ondes ultrasonores le long du tracé de rayon, par exemple un modèle issu de la Théorie de l'Elastodynamique Géométrique ou de la Théorie Géométrique de la Diffraction, ce qui permet d'en déduire les variations d'amplitude et de phase d'une onde ultrasonore portée par le tracé de rayon déterminé. Un exemple de résultat des étapes 112 et 114 est par exemple illustré sur la figure 6, dans un cas simple de tracé de rayon déterminé entre le point d'émission PE et le point de réception PR sans aucun autre point de passage imposé. Le tracé de rayon obtenu, à l'aide d'une seule application de l'algorithme de Dijkstra entre PE et PR, porte la référence T. Il est associé à un temps de vol égal à la somme des trajets élémentaires sélectionnés par l'algorithme de Dijkstra. Il est également associé à une liste des natures de ces trajets élémentaires. Dans l'exemple illustré par la figure 6, entre le point d'émission PE et un point P1 de la surface externe de la pièce d'étude 18, la nature du trajet est définie par un mode de propagation volumique longitudinal dans le matériau PMMA du sabot d'émission 16. Entre le point P1 et un point P2 d'extrémité du défaut D1, la nature du trajet est définie par un mode de propagation volumique longitudinal dans l'acier austénitique de la pièce d'étude 18. Entre le point P2 et un point P3 de la surface externe de la pièce d'étude 18 situé dans un arrondi de l'affouillement 20, la nature du trajet est définie par un mode de propagation volumique longitudinal dans l'acier austénitique de la pièce d'étude 18. Entre le point P3 et un autre point P4 de la surface externe de la pièce d'étude 18 situé dans le même arrondi de l'affouillement 20, la nature du trajet est définie par un mode de propagation surfacique longitudinal le long de la surface externe de la pièce d'étude 18. Enfin, entre le point P4 et le point de réception PR, la nature du trajet est définie par un mode de propagation volumique longitudinal dans l'acier austénitique de la pièce d'étude 18. Le procédé décrit précédemment peut être répété entre les étapes 110 et 114 pour une pluralité de points de réception PR situés le long de la section S mentionnée en référence à la figure 1. On obtient ainsi par exemple la pluralité de tracés de rayons illustrée sur la figure 7. Chacun de ces tracés est associé à un temps de vol et à une liste de natures des trajets élémentaires qui le composent. On remarque que ces tracés de rayons prennent en compte des diffractions en surface de la pièce d'étude 18 au niveau de l'affouillement 20 et sur les deux défauts D1, D2. Le premier mode de réalisation a été détaillé avec un exemple de structure bidimensionnelle, isotrope ou anisotrope. Mais l'invention est aisément généralisable à des exemples de structures tridimensionnelles, isotropes ou anisotropes. Il suffit d'adapter les étapes de discrétisation 102, 104 pour réaliser ces discrétisations selon un maillage tridimensionnel. Ensuite, les étapes de calculs des matrices M et V obéissent aux mêmes règles, les points de diffraction secondaire obtenus par discrétisation pouvant toujours être ordonnés selon leur distance au point d'émission PE. Le premier mode de réalisation permet de déterminer un tracé de rayon entre le point d'émission PE et le point de réception PR, y compris en imposant des contraintes de points de passage entre ces deux points. En revanche, la figure 8 illustre une contrainte imposée au tracé de rayon T à déterminer qui ne peut pas être traitée simplement par le premier mode de réalisation. Il s'agit d'imposer à ce trajet T qu'il passe par un unique point d'un ensemble de points prédéterminés, par exemple un ensemble de points de discrétisation d'une surface contre laquelle une onde ultrasonore longitudinale est susceptible de se réfléchir. Dans l'exemple illustré par la figure 8, cet ensemble de points constitue la surface externe supérieure R de la pièce d'étude 18, dans sa partie située entre l'affouillement 20 et la section S (qui a été décalée vers la droite de la figure pour illustrer clairement la nouvelle contrainte). Le principe est donc de déterminer le tracé de rayon optimisant le temps de vol entre PE et PR sous la contrainte d'une réflexion d'onde longitudinale contre la surface R. Un contournement possible de cette nouvelle contrainte complexe serait d'imposer un point de passage sur la surface R, ce que le premier mode de réalisation sait résoudre. Cependant, la position de ce point de passage ne peut pas être choisie au hasard et doit obéir aux lois de la physique des ondes pour que le tracé de rayon déterminé représente la propagation effective de l'onde et ait ainsi un sens physique : plus précisément, il doit respecter la loi de Snell-Descartes, qui dit que l'angle d'incidence et l'angle de réflexion sont égaux au point de réflexion sur la surface R. Il faudrait donc calculer a priori ce point de réflexion, puis le définir comme point de passage, ce qui ferait perdre tout son intérêt au procédé de détermination de tracé selon le premier mode de réalisation de l'invention. Une autre solution consiste à proposer un deuxième mode de réalisation de l'invention, tel que celui illustré sur la figure 9, capable de prendre en compte cette nouvelle contrainte de passage parmi un unique point d'un ensemble prédéterminé de points. Grâce à ce deuxième mode de réalisation, il est en outre possible, comme cela sera détaillé ci-dessous, de prendre en compte plusieurs quadruplets de vitesses de propagation ou de surfaces des lenteurs le long d'un tracé de rayon à déterminer. La première étape 200 de ce deuxième mode de réalisation d'un procédé de détermination de tracé de rayon selon l'invention diffère de l'étape 100 du premier mode de réalisation en ce que seules les données d'entrée 38, 40 et le pas de discrétisation défini dans les données 42 sont acquises. Les paramètres physiques des modes de propagation ne sont pas exploités à ce stade du procédé. Cette première étape d'acquisition 200 est réalisée par exécution du premier programme d'ordinateur 34. Les étapes de discrétisation 202 et 204 suivantes sont identiques aux étapes 102 et 104 et sont également réalisées par exécution du premier programme d'ordinateur 34. Au cours de l'étape de construction 206 suivante, réalisée par exécution du premier programme d'ordinateur 34, le graphe orienté intermédiaire est construit par calcul de la matrice M le modélisant comme à l'étape 106. En revanche, dans ce deuxième mode de réalisation, un coefficient quelconque de la matrice M a pour valeur, non plus un temps de vol, mais simplement la distance P,P, entre le i-ème point P, et le j-ème point P, de diffraction secondaire puisque les paramètres physiques définis dans les données d'entrée 42 ne sont pas encore pris en compte. Par ailleurs, les mêmes règles que celles appliquées à l'étape 106, définissant certains trajets élémentaires comme impossibles donc de distances infinies, sont appliquées dans cette étape 206 également. Le calcul de la matrice M à l'étape 206 se fait donc de la façon suivante : 1. v (i, j), j Après exécution des étapes 200 à 208, le premier programme d'ordinateur 34 fournit les matrices calculées M et V et les enregistre dans le fichier 46. Au cours d'une étape de mise à jour 210A, similaire à l'étape 110 et réalisée par exécution du deuxième programme d'ordinateur 36, plusieurs groupes disjoints de points de passage définis dans les données d'entrée 44, parmi lesquels le point d'émission PE et le point de réception PR, sont ajoutés au graphe orienté intermédiaire modélisé par la matrice M. La notion de groupe de points de passage selon ce deuxième mode de réalisation généralise la notion de point de passage selon le premier mode de réalisation en ce qu'un groupe de points de passage est défini comme étant un ensemble de points de diffraction secondaire voisins, définissant une courbe (en 2D) ou une surface (en 3D) discrétisée, dans lequel un unique point non prédéterminé doit être sélectionné comme point de passage imposé pour le tracé de rayon à déterminer. Un point de passage imposé, tel que les points PE et PR, est ainsi un cas particulier de groupe de points de passage à un seul élément. On note L le nombre total de points de diffraction secondaire contenus dans l'union de ces groupes disjoints. Ces points de diffraction secondaire sont ajoutés au graphe orienté intermédiaire modélisé par la matrice M s'ils n'y sont pas déjà présents. Comme dans le cas particulier des points de passage évoqué dans le premier mode de réalisation, il est probable qu'une partie des points des groupes de points de passage fassent déjà partie de la matrice M, mais pas nécessairement. Par exemple, une surface de réflexion imposée définie comme un groupe de points de passage peut être discrétisée selon un pas de discrétisation différent de celui des interfaces diffractantes définies par les données 38, même si elle constitue une partie de ces interfaces diffractantes. Si ce pas de discrétisation est deux fois plus petit par exemple, alors un point sur deux de ce groupe sera déjà présent dans la matrice M et les autres points seront à ajouter. Au cours de l'étape 210A, il est donc nécessaire de déterminer quels sont les K points, parmi les L définis dans les données d'entrée 44, qui doivent être ajoutés à la matrice M, avec 0 « K « L. On calcule ainsi la matrice Mf de taille (N+K)x(N+K) dans laquelle les points des groupes de points de passage non encore présents dans la matrice M sont insérés en fonction de leur éloignement par rapport au point d'émission PE. K lignes et K colonnes sont ainsi ajoutées à la matrice M pour former la matrice Mf. Les nouveaux coefficients ajoutés sont calculés selon les mêmes règles que celles appliquées à l'étape 206. La matrice Mf résultante est stockée dans le fichier 46. Comme pour la matrice M, la matrice Mf est, dans ce deuxième mode de réalisation, une matrice dont les coefficients représentent des distances entre points. Elle ne modélise pas encore le graphe orienté final à partir duquel sera déterminé le tracé de rayon. Au cours de cette même étape 210A, on calcule également la matrice Vf de taille (N+K)x(N+K) obtenue à partir de la matrice V. De même, les points des groupes de points de passage non encore présents dans la matrice V sont insérés dans la matrice Vf en fonction de leur éloignement par rapport au point d'émission PE. K lignes et K colonnes sont ainsi ajoutées à la matrice V pour former la matrice Vf dont les nouveaux coefficients ajoutés sont calculés comme à l'étape 208 et qui est elle aussi stockée dans le fichier 46. Ensuite, au cours d'une étape 210B, qui n'a pas d'équivalent dans le premier mode de réalisation et qui est réalisée par exécution du deuxième programme d'ordinateur 36, une matrice Mf' est calculée à partir des matrices Mf, Vf calculées à l'étape 210A et à l'aide des paramètres physiques définis dans les données d'entrée 42. Il est possible, au cours de cette étape, de prendre en compte plusieurs quadruplets de vitesses de propagation ou de surfaces des lenteurs, notamment un quadruplet de vitesses de propagation ou de surfaces des lenteurs par paire de groupes de points de passage successifs. Un coefficient mt,,, quelconque de cette matrice Mf' a pour valeur un temps de vol entre le i-ème point P, et le j-ème point P, de diffraction secondaire. La matrice Mf' modélise donc le graphe orienté final à partir duquel sera déterminé le tracé de rayon. Ce temps de vol est aisément calculable à partir de la distance P,P, entre les deux points considérés, donnée par le coefficient mf,,, de la matrice Mf, et à partir de la vitesse de propagation de l'onde ultrasonore dans le milieu du volume, de la facette de défaut ou de la surface où se trouve le trajet élémentaire [P,, PJ] et dans la direction de ce trajet élémentaire. La vitesse de propagation à appliquer est déterminée à l'aide des paramètres physiques définis dans les données d'entrée 42, éventuellement de la position du trajet élémentaire [P,, PJ] par rapport aux groupes de points de passage, et du coefficient vf,,, de la matrice Vf qui définit la nature d'une onde ultrasonore le long du trajet élémentaire [P,, Il convient en outre d'appliquer un certain nombre de règles supplémentaires définissant certains trajets élémentaires comme impossibles donc de temps de vol infini noté 09, compte tenu des groupes disjoints de points de passage définis dans les données d'entrée 44. Pour un groupe G de points de passage {P,, PJ}, j dans lequel on impose un unique point de passage à sélectionner parmi le ou les points P,, PJ, ces règles comportent : une première règle supplémentaire selon laquelle un trajet élémentaire ne peut pas relier deux points de diffraction situés de part et d'autre du groupe G dans la liste ordonnée des points de diffraction secondaire, et une deuxième règle supplémentaire selon laquelle un trajet élémentaire ne peut pas relier deux points de diffraction différents du groupe G. La première règle supplémentaire a pour conséquence que v (k,1), k j nek,1 = 09 (ou -1). La deuxième règle supplémentaire a pour conséquence que v (k,1), i k j, i 1 j et k # 1 = mfkj = 00 (ou -1). Le calcul de la matrice Mf' à l'étape 210B se fait donc de la façon suivante : 1. v (k,1), mfk,1 = 00 = mf'k,l = 00 (ou -1), 2. v (k,1), k = 1 = rem = 0, 3. v (k, 1), k j nek,1 = cx), 4. V (k, 1), i k j, 1 j et k 1 rem = 09, et 5. pour tout autre coefficient mf', , = , où Vprop est la vitesse de vprop propagation associée au milieu volumique ou surfacique dans lequel se trouve le trajet élémentaire [P,, PJ] et éventuellement à la direction et au sens de [P,, Pl cette vitesse de propagation étant déterminable à partir de la valeur du coefficient vf,,, de la matrice Vf et des paramètres physiques des données d'entrée 42 (i.e. quadruplets de vitesses de propagation ou de surfaces des lenteurs). L'étape 210B est suivie d'une étape 212 de parcours optimisé du graphe orienté final modélisé par la matrice Mf'. Comme dans le premier mode de réalisation, ce parcours optimisé est un parcours de moindre coût entre le point d'émission PE et le point de réception PR, en passant par tous les autres points de passage ou groupes de points de passage imposés, les coûts étant représentés par les coefficients de la matrice Mf', c'est-à-dire les poids des liens orientés du graphe qu'elle modélise. C'est de nouveau un problème bien connu de la théorie des graphes. Mais cette fois-ci, conformément à ce deuxième mode de réalisation, l'optimisation de parcours est réalisée une seule fois entre les points PE et PR puisque les coefficients de la matrice Mf' intègrent les contraintes de passages définies par les points de passage ou groupes de points de passage imposés (étapes 210B.3 et 210B.4 précitées). Elle est de nouveau réalisée par application de l'algorithme de Dijkstra. Au cours d'une dernière étape 214, le résultat de l'optimisation 212 est traité pour fournir les données de sortie 48, 50 et 52. Cette étape, identique à l'étape 114 à l'exception de la concaténation de plusieurs parcours optimisés qui n'est plus utile dans ce deuxième mode de réalisation, n'a pas besoin d'être détaillée. La figure 10 illustre un exemple de configuration de la matrice Mf', de taille n=N+K, dans le cas particulier où un seul groupe G de points de passage imposés est défini dans les données d'entrée 44 en dehors des points PE et PR. Ce groupe G impose au tracé de rayon de passer par l'un parmi trois points notés Pi, Pi+15 Pi+2. La matrice Mf' de la figure 10 comporte une première zone ZA1, strictement inférieure à sa diagonale, dans laquelle tous les coefficients sont à la valeur infinie, par application de la règle de l'étape 206.1 à la matrice Mf. Sa diagonale comporte des coefficients tous nuls, par application de la règle de l'étape 2106.2. Elle comporte une deuxième zone ZA2 dans laquelle tous les coefficients sont à la valeur infinie, par application de la règle de l'étape 210B.3. Elle comporte une troisième zone ZA3 dans laquelle tous les coefficients sont à la valeur infinie, par application de la règle de l'étape 210B.4. Enfin, par application de la règle de l'étape 210B.5, les autres coefficients sont calculés à l'aide de la matrice Vf et d'un quadruplet {T} de vitesses de propagation (cas de milieux isotropes) ou de surfaces des lenteurs (cas de milieux anisotropes). La figure 11 illustre un autre exemple de configuration de la matrice Mf', de taille n=N+K, dans le cas particulier où deux groupes G1, G2 de points de passage imposés sont définis dans les données d'entrée 44 en dehors des points PE et PR. Le groupe G1 impose au tracé de rayon de passer par l'un parmi trois points notés P,, Pi+15 Pi+2. Le groupe G2 impose au tracé de rayon de passer par l'un parmi deux points notés Pm, Pm+1. En outre, selon une variante de réalisation possible de l'invention évoquée précédemment, plusieurs quadruplets de vitesses de propagation ou de surfaces des lenteurs sont définis et applicables entre les groupes de points de passage imposés. Par exemple, dans l'exemple de la figure 11, un premier quadruplet {T,} s'applique entre le point d'émission et le groupe G1, un deuxième quadruplet {T2} s'applique entre le groupe G1 et le groupe G2, et le premier quadruplet {T,} s'applique de nouveau entre le groupe G2 et le point de réception PR. On peut ainsi définir des variations de modes de propagation entre les différentes contraintes imposées par les données d'entrée 44. Notamment, cela rend possible le fait d'autoriser des conversions de mode à l'occasion de réflexions contre certaines surfaces définies en tant que groupes de points de passage, par exemple une conversion de mode volumique longitudinal en mode volumique transversal ou inversement à l'occasion d'une réflexion. La matrice Mf' de la figure 11 comporte une première zone ZB1, strictement inférieure à sa diagonale, dans laquelle tous les coefficients sont à la valeur infinie, par application de la règle de l'étape 206.1 à la matrice Mf. Sa diagonale comporte des coefficients tous nuls, par application de la règle de l'étape 210B.2. Elle comporte une deuxième zone ZB2 dans laquelle tous les coefficients sont à la valeur infinie, par application de la règle de l'étape 210B.3 au groupe G1. Elle comporte une troisième zone ZB3 dans laquelle tous les coefficients sont à la valeur infinie, par application de la règle de l'étape 210B.3 au groupe G2. Elle comporte une quatrième zone ZB4 dans laquelle tous les coefficients sont à la valeur infinie, par application de la règle de l'étape 210B.4 au groupe G1. Elle comporte une cinquième zone ZB5 dans laquelle tous les coefficients sont à la valeur infinie, par application de la règle de l'étape 210B.4 au groupe G2. Enfin, par application de la règle de l'étape 210B.5, les autres coefficients sont calculés à l'aide de la matrice Vf et des quadruplets {Ti}, {T2} de vitesses de propagation (cas de milieux isotropes) ou de surfaces des lenteurs (cas de milieux anisotropes). La figure 12 illustre un autre exemple de configuration de la matrice Mf', de taille n=N+K, dans le cas particulier où un seul groupe G à unique point de passage imposé est défini dans les données d'entrée 44 en dehors des points PE et PR. Ce groupe G impose au tracé de rayon de passer par un unique point noté P. La matrice Mf' de la figure 12 comporte une première zone ZC1, strictement inférieure à sa diagonale, dans laquelle tous les coefficients sont à la valeur infinie, par application de la règle de l'étape 206.1 à la matrice Mf. Sa diagonale comporte des coefficients tous nuls, par application de la règle de l'étape 2106.2. Elle comporte une deuxième zone ZC2 dans laquelle tous les coefficients sont à la valeur infinie, par application de la règle de l'étape 210B.3. La règle de l'étape 210B.4 ne s'applique pas parce que le groupe G ne comporte qu'un seul point. Enfin, par application de la règle de l'étape 210B.5, les autres coefficients sont calculés à l'aide de la matrice Vf et d'un quadruplet {T,} de vitesses de propagation (cas de milieux isotropes) ou de surfaces des lenteurs (cas de milieux anisotropes). Il convient de remarquer que la matrice Mf', qui est triangulaire supérieure, peut se décomposer dans cette configuration en deux sous-matrices triangulaires supérieures. Cela signifie qu'il est possible de traiter cette configuration indifféremment selon le premier ou le deuxième mode de réalisation. Selon le deuxième mode de réalisation, une seule optimisation de trajet est réalisée sur toute la matrice Mf'. Selon le premier mode de réalisation, deux optimisations de trajets sont réalisées indépendamment sur les deux sous-matrices, entre le point d'émission PE et le point de passage imposé P, d'une part, puis entre P, et le point de réception PR d'autre part. La figure 13 illustre un autre exemple de configuration de la matrice Mf', de taille n=N+K, dans le cas particulier où aucune contrainte de point de passage ou groupe de points de passage n'est imposée. Dans ce cas très simple, la matrice Mf' de la figure 13 comporte une première zone ZD1, strictement inférieure à sa diagonale, dans laquelle tous les coefficients sont à la valeur infinie, par application de la règle de l'étape 206.1 à la matrice Mf. Sa diagonale comporte des coefficients tous nuls, par application de la règle de l'étape 210B.2. Enfin, par application de la règle de l'étape 210B.5, les autres coefficients sont calculés à l'aide de la matrice Vf et d'un quadruplet {T} de vitesses de propagation (cas de milieux isotropes) ou de surfaces des lenteurs (cas de milieux anisotropes). Dans cette configuration particulière, les premier et deuxième modes de réalisation sont équivalents. Il apparaît clairement qu'un procédé de détermination d'un tracé de rayon tel que l'un de ceux décrits précédemment permet de simuler avec précision la propagation d'une onde ultrasonore le long d'un tracé de rayon en prenant en compte, lors de la détermination de ce tracé, les modifications de son mode de propagation, de sa polarisation ainsi éventuellement que de son amplitude et de sa phase. Selon ce procédé avantageux, les vitesses de propagation spécifiques à chaque mode de propagation possible sont prédéterminées, connues et prises en compte en tant que données d'entrée, que le milieu soit isotrope (une valeur de vitesse par mode de propagation) ou anisotrope (une surface des lenteurs par mode de propagation). A chaque zone homogène de la structure inspectée est associé un mode de propagation de sorte que les phénomènes de conversions de modes aux interfaces puissent être simulés. En incluant également la possibilité de modes de propagations surfaciques, une meilleure simulation des effets des interfaces diffractantes est en outre assurée. En contrôle non destructif, on observe par exemple souvent que lorsqu'une onde volumique atteint une surface avec des incidences particulières, cela engendre la création d'ondes se propageant le long de la surface comme l'onde de Rayleigh ou l'onde de tête. En traitant en outre spécifiquement les phénomènes de diffraction sur défauts internes à la structure inspectée, un procédé tel que l'un de ceux décrits précédemment affine encore la simulation. En étant capable enfin de modéliser des contraintes de type points de passage ou groupes de points de passage imposés, et en permettant des conversions de modes de propagations avant et après ces points ou groupes de points de passage, un procédé selon tel que l'un de ceux décrits précédemment permet de modéliser un grand nombre de contraintes. Un procédé selon l'invention permet ainsi de sélectionner l'un ou plusieurs des effets précités dans une structure comportant une pièce d'étude donnée et de calculer de façon générique le trajet représentant le type de propagation résultant de ces effets afin d'obtenir le temps de vol et éventuellement la phase, la forme du front et l'amplitude des ondes ultrasonores en des points d'observation situés à l'intérieur ou à l'extérieur de la pièce. En variant les effets et contraintes précités, une grande variété de tracés de rayons peut être déterminée. Cette capacité à simuler une grande variété de tracés de rayons est très utile à la compréhension de la physique d'un contrôle non destructif et donc au diagnostic, car généralement de nombreuses ondes ultrasonores se propagent dans les pièces inspectées et sont issues d'interactions très variées avec le milieu de propagation. En particulier, le deuxième mode de réalisation permet, par son formalisme plus élaboré que le premier, de déterminer des tracés de rayons avec ou sans contraintes, sans que la présence de contraintes ne ralentisse son exécution. En outre, il ne nécessite qu'une seule exécution de l'algorithme de Dijkstra ou autre algorithme équivalent d'optimisation de trajet dans un graphe orienté. Ainsi, même lorsque le premier mode de réalisation est applicable, le deuxième mode de réalisation devient clairement avantageux dès lors que le tracé de rayon à déterminer inclut au moins deux points de passage imposés en plus du point d'émission PE et du point de réception PR. On notera par ailleurs que l'invention n'est pas limitée aux modes de réalisation décrits précédemment. Il apparaîtra en effet à l'homme de l'art que diverses modifications peuvent être apportées aux modes de réalisation décrits ci-dessus, à la lumière de l'enseignement qui vient de lui être divulgué. Dans les revendications qui suivent, les termes utilisés ne doivent pas être interprétés comme limitant les revendications aux modes de réalisation exposés dans la présente description, mais doivent être interprétés pour y inclure tous les équivalents que les revendications visent à couvrir du fait de leur formulation et dont la prévision est à la portée de l'homme de l'art en appliquant ses connaissances générales à la mise en oeuvre de l'enseignement qui vient de lui être divulgué.25 REVENDICATIONS1. Procédé de détermination d'un tracé (T) de rayon d'une onde ultrasonore entre un point d'émission (PE) et un point de réception (PR) par simulation d'une inspection ultrasonore de structure (12) à au moins une interface diffractante disposée entre ces points d'émission et de réception, comportant les étapes suivantes : discrétisation (102, 104 ; 202, 204) de la structure (12) en une pluralité de points de diffraction secondaire, et sauvegarde de ces points de diffraction secondaire en mémoire (32C), construction (106, 108, 110 ; 206, 208, 210A, 210B), par un processeur (30) ayant accès à la mémoire (32A, 32B, 32C) en lecture/écriture, d'un graphe orienté dont les noeuds comportent le point d'émission (PE), le point de réception (PR) et l'ensemble des points de diffraction secondaire, et dont les liens orientés sont associés à des poids représentant des temps de propagation entre les noeuds qu'ils relient, sélection (112 ; 212) par le processeur (30), parmi l'ensemble des tracés de rayons possibles parcourant le graphe orienté du point d'émission (PE) au point de réception (PR), du tracé de rayon minimisant la somme des poids des liens parcourus, caractérisé en ce que : plusieurs valeurs de vitesses de propagation (42) associées à plusieurs modes de propagation incluant au moins un mode de propagation de surface le long de l'interface diffractante sont conservées en mémoire (32A), et lors de la construction (106, 108, 110 ; 206, 208, 210A, 210B) du graphe orienté, chaque lien orienté retenu est associé (106, 110 ; 210B) par le processeur (30) à une vitesse de propagation parmi celles conservées en mémoire (32A) en fonction de la localisation de ce lien par rapport à l'interface diffractante, pour le calcul de son poids. 2. Procédé de détermination d'un tracé (T) de rayon d'une onde ultrasonore selon la revendication 1, dans lequel l'ensemble des points de diffraction secondaire conservés en mémoire (32C) est constitué de points de ladite au moins une interfacediffractante et le processeur (30) est conçu pour construire le graphe en excluant tout lien orienté traversant cette interface diffractante. 3. Procédé de détermination d'un tracé (T) de rayon d'une onde ultrasonore selon la revendication 1 ou 2, dans lequel ladite au moins une interface diffractante comporte au moins l'un des éléments de l'ensemble constitué d'une surface de pièce (18) disposée dans la structure observée (12) entre le point d'émission (PE) et le point de réception (PR), d'une interface entre zones hétérogènes d'une même pièce (18) et d'un défaut de structure (D1, D2). 4. Procédé de détermination d'un tracé (T) de rayon d'une onde ultrasonore selon la revendication 3, dans lequel, lorsque ladite au moins une interface diffractante comporte un défaut de structure (D1, D2) modélisé par au moins une facette polygonale plane présentant des extrémités, les seuls points de diffraction secondaire retenus sur ce défaut de structure (D1, D2) sont situés aux extrémités de chaque facette le modélisant. 5. Procédé de détermination d'un tracé (T) de rayon d'une onde ultrasonore selon l'une quelconque des revendications 1 à 4, dans lequel : chaque zone de propagation homogène considérée comme isotrope dans la structure (12) est associée à au moins un mode de propagation lui-même associé à une vitesse de propagation caractéristique de ce mode de propagation, et chaque zone de propagation homogène considérée comme anisotrope dans la structure (12) est associée à au moins un mode de propagation lui-même associé à une pluralité de vitesses de propagation en fonction d'une direction et d'un sens de propagation, cette pluralité de vitesses de propagation étant conservée en mémoire sous la forme d'une surface des lenteurs. 6. Procédé de détermination d'un tracé (T) de rayon d'une onde ultrasonore selon l'une quelconque des revendications 1 à 5, dans lequel la construction (106, 108, 110 ; 206, 208, 210A, 210B) du graphe orienté comporte les étapes suivantes : construction (106 ; 206) d'un graphe orienté intermédiaire dont les noeuds comportent uniquement les points de diffraction secondaire et dont les liens orientés sont associés à des poids représentant des distances ou temps de propagation entre les noeuds qu'ils relient, puisconstruction (110 ; 210A, 210B) du graphe orienté en complétant le graphe orienté intermédiaire par l'ajout d'une pluralité de points de passage ou groupes de points de passage imposés incluant le point d'émission (PE) et le point de réception (PR). 7. Procédé de détermination d'un tracé (T) de rayon d'une onde ultrasonore selon la revendication 6, dans lequel la construction (106 ; 206) du graphe orienté intermédiaire comporte les étapes suivantes : classement des points de diffraction secondaire dans l'ordre croissant de leur éloignement au point d'émission (PE), calcul des coefficients d'une matrice M modélisant ce graphe orienté intermédiaire, chaque coefficient de cette matrice M portant la valeur d'une distance ou d'un temps de propagation d'onde ultrasonore entre les i-ème et j-ème points de diffraction secondaire dans l'ordre du classement, et dans lequel les règles suivantes sont appliquées : tout coefficient pour lequel j 8. Procédé de détermination d'un tracé (T) de rayon d'une onde ultrasonore selon l'une quelconque des revendications 1 à 7, entre le point d'émission (PE) et le point de réception (PR) via au moins un groupe de points intermédiaires parmi lesquels un unique point de passage imposé à déterminer, dans lequel la construction (210A, 210B) du graphe orienté en complétant le graphe orienté intermédiaire comporte les étapes suivantes : classement des points de diffraction secondaire, incluant le groupe de points intermédiaires, dans l'ordre croissant de leur éloignement au point d'émission (PE), calcul des coefficients d'une matrice Mf' modélisant le graphe orienté, chaque coefficient mf',,, de cette matrice Mf' portant la valeur d'un temps de propagation d'onde ultrasonore entre les i-ème et j-ème points de diffraction secondaire dans l'ordre du classement, et dans lequel les règles supplémentaires suivantes sont appliquées :tout coefficient pour lequel les i-ème et j-ème points de diffraction secondaire sont situés de part et d'autre du groupe de points intermédiaires est d'une valeur indicatrice d'un temps de propagation infini, et tout coefficient pour lequel les i-ème et j-ème points de diffraction secondaire sont différents et dans le groupe de points intermédiaires est d'une valeur indicatrice d'un temps de propagation infini. 9. Procédé de détermination d'un tracé (T) de rayon d'une onde ultrasonore selon l'une quelconque des revendications 1 à 7, entre le point d'émission (PE) et le point de réception (PR) via au moins un autre point de passage imposé, dans lequel les étapes de construction (106) du graphe intermédiaire et de sélection (112) du tracé de rayon optimal sont exécutées indépendamment pour chaque paire de points de passage imposés successifs entre le point d'émission (PE) et le point de réception (PR). 10. Procédé de détermination d'un tracé (T) de rayon d'une onde ultrasonore selon l'une quelconque des revendications 1 à 9, entre le point d'émission (PE) et le point de réception (PR) via au moins un autre point de passage imposé ou groupe de points intermédiaires parmi lesquels un unique point de passage imposé à déterminer, dans lequel un p-uplet de modes de propagation est défini pour chaque paire de points de passage imposés ou groupes de points intermédiaires successifs, p étant le nombre de zones volumiques ou surfaciques de propagation homogènes dans la structure (12). 11. Programme d'ordinateur (34, 36) téléchargeable depuis un réseau de communication et/ou enregistré sur un support (32B) lisible par ordinateur et/ou exécutable par un processeur (30), caractérisé en ce qu'il comprend des instructions pour l'exécution des étapes d'un procédé de détermination d'un tracé (T) de rayon d'une onde ultrasonore selon l'une quelconque des revendications 1 à 10, lorsque ledit programme (34, 36) est exécuté sur un ordinateur (28).