Ce travail se situe à l’intersection des trois disciplines : méthodes numériques, techniques expérimentales et du machine learning, a pour but de proposer une famille de techniques d’identification par analyse inverse des lois de comportement en mécanique. Dans le domaine d’identification des matériaux, l’indentation instrumentée est particulièrement attractive, car elle permet de procéder à des essais non-destructifs sur l’échantillon ou sur une structure en service. L’essai d’indentation, similaire à un test de dureté, consiste à enfoncer la pointe de l’indenteur à une faible profondeur dans la matière tout en enregistrant le déplacement en fonction de la force appliquée. L’identification des propriétés élastoplastiques des matériaux est basée alors sur l’exploitation de la courbe force-déplacement (courbe P-h). Toutefois, le problème inverse est souvent mal posé et des problèmes d’unicité mènent à la notion de paires de "matériaux mystiques" produisant, dans des conditions d’essai donnés, des courbes P-h identiques, malgré des propriétés différentes. L’idée de notre travail est de compléter la procédure d’identification en faisant appel à des dispositifs expérimentaux récents, notamment à la microscopie laser, permettant de mesurer la carte 3D de l’empreinte résiduelle obtenue après le retrait de l’indenteur. Pour aborder la question de la richesse d’information de l’empreinte par rapport à la courbe P-h seule, nous proposons de construire, dans un espace affine réduit, la variété des formes d’empreinte admissibles au sens d’une loi de comportement et du modèle d’éléments finis de l’essai. La mesure de la dimension intrinsèque nous indique alors le nombre maximal de paramètres potentiellement identifiables. Cela nous permet de proposer et de valider numériquement des nouveaux procédés expérimentaux, plus représentatifs, à partir des données synthétiques, ainsi que des algorithmes d’identification associés. La prise en compte de l’erreur de modèle et de l’erreur de mesure, nous mène ensuite à proposer un ensemble d’algorithmes de projection d’empreintes expérimentales, réalisées en collaboration avec l’INSA de Rennes sur la variété synthétique. Nous abordons alors le problème d’identification des propriétés d’écrouissage de plusieurs matériaux de complexité croissante et départageons des "jumeaux mystiques" par des essais de multi-indentation, basés sur l’exploitation de l’empreinte seule ou en complément de la courbe P-h. / The thesis lies at the intersection of three disciplines : numerical methods, experimental techniques, and machine learning. The primary aim of this work is to develop a group of algorithms for characterization by inverse analysis of a material’s constitutive law. In the field of material characterization, indentation test is especially attractive since it is considered non-destructive, and may be performed even on a structure in service. The test, similar to a hardness test, consists in penetrating an indenter into the surface of the material. The force exerted on the indenter is recorded against the penetration depth over a series of time instants, leading to a force-displacement (P-h) curve, which is the most frequently used source of information for the identification of material properties. However, the inverse problem based solely on this curve tends to be ill-posed, leading to nonunique identification solution, i.e., the "mystical material pair", for whom the corresponding force-displacement curves are almost identical despite the very different material properties. The basic idea is then to complete the identification process with innovative experimental measurements, such as laser microscope, which allows measuring the 3D residual imprint after the withdrawal of the indenter. To address the advantage of this measurement over P-h curve, we propose to construct, within a reduced affine space, a manifold of shapes admissible to the postulated constitutive law, experimental and simulation setups, based on synthetic data. The intrinsic dimensionality of the manifold limits the number of identifiable parameters allowing to validate numerically experimental procedures. Considering both the model and measurement errors, we develop a series of local manifold learning algorithms to solve the inverse problem iteratively for experimental results obtained in cooperation with INSA de Rennes. This approach allows us to characterize diverse metallic materials of increasing complexity, based on actual experimental measurements. For example, for the Hollomon’s law, the mystical pair is alleviated in using a single imprint, while for the Voce law, a multi-depth experimental protocol is proposed to differentiate mystical siblings.

Cette thèse porte sur le développement d'un traducteur ultrasonore multi-élément capacitif à couplage air (MEACUT) et son utilisation dans le domaine du contrôle non destructif (CND) de matériaux. Un modèle est employé pour simuler numériquement ce traducteur, et pour optimiser sa conception. Un prototype est ensuite fabriqué, puis caractérisé expérimentalement pour quantifier ses performances. Son originalité réside dans le fait qu'il possède une large bande passante en fréquence, tout en offrant la possibilité d'une focalisation variable. Ce prototype est alors employé pour la détection d'un endommagement causé par impact, dans une plaque composite. Il est clairement constaté que la résolution spatiale du procédé d'inspection employé (C-scan) est fortement améliorée grâce aux performances techniques du MEACUT. Enfin, un modèle hybride 3D est développé pour simuler, rapidement et intégralement, ce procédé de CND. Le très bon accord obtenu entre prédictions numériques et mesures expérimentales laisse présager que cet outil de simulation pourra servir à mettre au point d'autres expérimentations de CND, qui pourront à leur tour exploiter les performances du MEACUT. / This thesis focuses on the development of a Multi-Element Air-coupled Capacitive Ultrasonic Transducer (MEACUT) and its use in the field of non-destructive testing (NDT) of materials. A numerical model is used to simulate the translator, and to optimize its design. A prototype is then built and experimentally characterized to quantify its performance. Its originality lies in the fact that it has a broad frequency bandwidth while offering the possibility of a dynamic focusing. This prototype is then used for the detection of damage caused by impact, in a composite plate. It is clearly found that the spatial resolution of the inspection process employed (C-scan) is greatly improved thanks to the technical performance MEACUT. Finally, a 3D hybrid model is developed to simulate quickly and fully, the process of NDT. The good agreement obtained between numerical predictions and experimental measurements suggests that this simulation tool can be used to develop other NDT experiments, which may in turn exploit the performance of MEACUT.

Le principe du contrôle par flux de fuite magnétique (Magnetic Flux Leakage MFL) consiste à aimanter la pièce à contrôler par un champ magnétique et à détecter à l'aide d'un capteur magnétique les fuites des lignes du champ qui résultent de la présence d'un défaut dans la pièce. Dans le but d'améliorer les performances d'un dispositif de détection, le CEA et la société Vallourec collaborent pour développer un modèle numérique dédié au contrôle virtuel des défauts longitudinaux dans les tubes ferromagnétiques. Le dispositif expérimental comprend un circuit magnétique tournant à une vitesse constante autour du tube qui défile. Dans le cadre de cette thèse, on débute le problème de la modélisation sans tenir compte des effets de la vitesse de rotation, il s'agit donc de résoudre un problème d'électromagnétisme en régime magnétostatique.Pour résoudre ce problème, on propose de comparer une approche semi-analytique basée sur le formalisme des équations intégrales (EI) et une approche purement numérique utilisant les éléments finis (EF).Dans la première partie de cette thèse, après avoir établi le formalisme théorique par EI, un premier modèle considérant des matériaux ferromagnétiques à perméabilité magnétique constante (régime linéaire) a été mis en œuvre en 2D. Ce modèle a été appliqué pour un exemple de système extrait de la littérature et validé numériquement par une comparaison des résultats EI/EF. Pour une meilleure détection, il est opportun de saturer magnétiquement la pièce. Le matériau ferromagnétique est alors caractérisé par une courbe B(H) non-linéaire. Par conséquent, la deuxième partie de la thèse a été consacrée à la mise en œuvre du modèle en régime non linéaire qui tient compte de cette caractéristique.Différentes méthodes de discrétisation ont été étudiées afin de réduire le nombre d'inconnues et le temps de calcul. L'originalité de la thèse réside dans l'utilisation des fonctions d'interpolation d'ordre élevé (polynôme de Legendre) pour une discrétisation des équations intégrales par une approche de type Galerkin. Les premiers essais de validation numérique de ce modèle ont été effectués sur un système MFL simplifié. Des premiers essais de validation expérimentale pour des données obtenues par EF ont été effectués en deux phases : La première a consisté à vérifier le distribution du champ magnétique pour un tube sain et en régime magnétostatique. La deuxième phase a consisté à calculer la réponse d'un défaut dans le tube ferromagnétique en tenant en compte les effets éventuels de la rotation du circuit magnétique par rapport au tube. / The principle of the non destructive testing by magnetic flux leakage (MFL) is to magnetize the part to be inspected by a magnetic field and to detect a flaw thanks to magnetic leakage field lines due to the strong decreasing of the magnetic permeability in the flawed region. In order to improve the performance of detection, the CEA and the Vallourec society collaborate to develop a numerical model dedicated to the virtual NDT of longitudinal defects in ferromagnetic tubes. The experimental system includes a magnetic circuit rotating at a constant speed around the tube to be inspected. The modeling task is started without considering the effects of the rotational speed, so the magnetostatic regime is considered to solve the modeling problem. In the framework of this thesis, we propose to compare a semi-analytical approach based on the formalism of integral equations method (IEM) and a purely numerical approach using finite element method (FEM).In the first part of this thesis, the theoretical formalism was established. A first simple discretization scheme is been implemented in the linear regime considering a constant magnetic permeability. This first numerical model has been validated for a simplified MFL configuration extracted and modified from the literature.For better detection, it is wishable to magnetically saturate the piece under-test. The ferromagnetic material is then characterized by a B(H) curve. Therefore, the second part of the thesis was devoted to the implementation of the model in the non-linear regime that takes into account this non-linear characteristic. Different discretization schemes have been studied in order to reduce the number of unknowns and the computational time. The originality of the thesis lies in the use of basis function of high order (Legendre polynomials) associated to a Galerkin approach for the discretization of integral equations. The first numerical result has been validated on a simplified MFL system. The first results of the experimental validation based on simulated data obtain by FEM have been performed in two steps. The first one consists to verify the distribution of the magnetic field for a ferromagnetic tube without any defect and in the magnetostatic regime. The objective of the second one was to compute the response of the flaw and to evaluate the effects of the rotational speed of the magnetic circuit around the tube.

Les travaux de cette thèse portent sur la détection et le diagnostic des défauts naissants dans les systèmes d’ingénierie et industriels, par des approches statistiques non-paramétriques. Un défaut naissant est censé provoquer comme tout défaut un changement anormal dans les mesures des variables du système. Ce changement est imperceptible mais aussi imprévisible dû à l’important rapport signal-sur défaut, et le faible rapport défaut-sur-bruit caractérisant le défaut naissant. La détection et l’identification d’un changement général nécessite une approche globale qui prend en compte la totalité de la signature des défauts. Dans ce cadre, la divergence de Kullback-Leibler est proposée comme indicateur général de défauts, sensible aux petites variations anormales cachées dans les variations du bruit. Une approche d’analyse spectrale globale est également proposée pour le diagnostic de défauts ayant une signature fréquentielle. L’application de l’approche statistique globale est illustrée sur deux études différentes. La première concerne la détection et la caractérisation, par courants de Foucault, des fissures dans les structures conductrices. La deuxième application concerne le diagnostic des défauts de roulements dans les machines électriques tournantes. En outre, ce travail traite le problème d’estimation de l’amplitude des défauts naissants. Une analyse théorique menée dans le cadre d’une modélisation par analyse en composantes principales, conduit à un modèle analytique de la divergence ne dépendant que des paramètres du défaut. / This phD dissertation deals with the detection and diagnosis of incipient faults in engineering and industrial systems by non-parametric statistical approaches. An incipient fault is supposed to provoke an abnormal change in the measurements of the system variables. However, this change is imperceptible and also unpredictable due to the large signal-to-fault ratio and the low fault-to-noise ratio characterizing the incipient fault. The detection and identification of a global change require a ’global’ approach that takes into account the total faults signature. In this context, the Kullback-Leibler divergence is considered to be a ’global’ fault indicator, which is recommended sensitive to abnormal small variations hidden in noise. A ’global’ spectral analysis approach is also proposed for the diagnosis of faults with a frequency signature. The ’global’ statistical approach is proved on two application studies. The first one concerns the detection and characterization of minor cracks in conductive structures. The second application concerns the diagnosis of bearing faults in electrical rotating machines. In addition, the fault estimation problem is addressed in this work. A theoretical study is conducted to obtain an analytical model of the KL divergence, from which an estimate of the amplitude of the incipient fault is derived.

Induction hardening is a complex process that requires regular verification of the case depth and the microstructure. Today this verification is done by destructive testing where a cross section is analysed. The case depth is measured by several Vickers hardness impressions and the microstructure is analysed by a light optic microscope. This master thesis was performed in collaboration with Volvo GTO in Skövde, with the target to find out if an ultrasonic backscattering method is usable for analysing the case depths on induction hardened crankshafts. This was done by verifying the results provided by ultrasonic testing with the results provided by destructive testing. The measurements were done on several crankshafts with different case depths and steel composition. The result, for crankshafts with normal case depths, shows a good correlation between the destructive testing (DT) and the ultrasonic testing (UT). The mean value for the discrepancy was 0.39 mm with a standard deviation of 0.20 mm. The reproducibility of the method has a standard deviation of 0.074 mm.

Este trabalho descreve o desenvolvimento de um algoritmo para avaliar o grau de corrosão em paredes de oleodutos, através de ultrassom. O algoritmo desenvolvido será implementado em um circuito digital, fazendo parte de um sistema que poderá ser embarcado em pigs ultrassônicos para inspeção de dutos à procura de corrosão. O modelo matemático conhecido como filtro casado é usado como base do algoritmo, além de um método de detecção de picos para localização dos sinais de eco de ultrassom, provenientes das reflexões do sinal nas paredes dos dutos. Um dos problemas é a definição do sinal (e do filtro) que será utilizado, mas por meio de sinais sintetizados (signal design) uma grande gama de alternativas foi analisada e apenas algumas foram selecionadas para testes em experimentos. Por último os melhores sinais foram utilizados na inspeção de um corpo de prova feito à partir de um duto corroído, e os resultados são analisados. A detecção da primeira parede do oleoduto apresenta resultados promissores em comparação com a medição tridimensional das mesmas paredes. O sistema foi capaz de reconstruir o perfil da parede interna do duto com elevada verossimilhança, mesmo em regiões severamente corroídas. A detecção da segunda parede apresenta resultados satisfatórios para baixos níveis de corrosão, porém alguns problemas foram encontrados para níveis mais altos de corrosão, prejudicando a medição direta da espessura da parede nestas condições. / This work describes the development of an algorithm to evaluate the degree of corrosion in pipeline walls, using ultrasound. The developed algorithm will be implemented in a dedicated digital circuit that can be embedded in ultrasonic pigs for pipeline corrosion inspection. The algorithm is based on the mathematic model known as matched filter, besides a peak detection method for the location of ultrasound echoes, originated from the signal reflections in the pipeline walls. One problem is the determination of the signal (and filter) that should be used, but using a signal design approach a large number of possibilities could be analyzed and only a few was tested in experiments. Finally, the best signals were used to inspect a test specimen made of a corroded pipeline section, and the results were analyzed. The first wall detection shows promising results when compared to the results of a three-dimensional inspection of the same walls. It was possible to reconstruct the profile of the internal wall with good likelihood, even at badly corroded areas. The second wall detection shows some satisfactory results on low corrosion levels, but show some problems with high corrosion levels, compromising the direct measurement of the wall of the pipeline.

The diploma thesis describes non-destructive methods for diagnostics of reinforcement. It is focused on the method of electromagnetic indicators. In the first part the methods are described and compared. The profometer PM-630 is also described. The second part describes experimental measurements with profometer on reinforced concrete structures. The results are described and interpreted. The effects, which have influence on measurement, are described. The recommendations for practical measurement with profometer are described. The experimental models were created for calibration of profometer. The part of this thesis is also the determination of load bearing capacity of reinforced concrete bridge in Nosislav.

La caractérisation non destructive de l’adhérence des structures du type couche sur substrat est un enjeu industriel et académique important. Ce type d’échantillon est en effet utilisé pour de nombreuses applications et sa durée de vie dépend en grande partie de la qualité d’adhérence des films au substrat. Celle-ci modifie sensiblement le comportement dispersif des ondes acoustiques de surface se propageant dans de ce type de structure. Pour générer et détecter ces ondes, un dispositif Ultrasons-Laser a été privilégié. Dans un premier temps, nous avons cherché à contourner les difficultés d’interprétation rencontrées habituellement dans le contrôle de l’adhérence par ondes acoustiques de surface. Les variations d’épaisseur de la couche peuvent en effet avoir une influence sur la dispersion des ondes comparable à celle due à l’adhérence. Pour ce faire, des films polymères dont l’épaisseur est quasi-constante sont employés et apposés sur un substrat en aluminium. Ces films possèdent en plus la propriété d’être transparents. Cela permet de focaliser l’impulsion laser générant les ondes acoustiques à travers le film, directement à la surface du substrat et de placer ainsi la source acoustique à l’interface film-substrat. L’influence de la position de la source sur le comportement dispersif des ondes acoustiques de surface et par conséquent sur le contrôle de la qualité d’adhérence est alors étudiée expérimentalement ainsi qu’au travers de simulations par éléments finis. Finalement, une caractérisation de l’adhérence de différents échantillons est effectuée grâce aux courbes de dispersion obtenues à l’aide de la méthode Matrix-Pencil appliquée aux résultats expérimentaux. En utilisant un algorithme d’inversion, les raideurs d’interface caractéristiques de l’adhérence des échantillons analysés sont estimées. / The non-destructive characterization of the adhesion of layer-on-substrate structures is an important issue in industrial and academic domains. This type of sample is indeed used for many applications and its lifetime depends mainly on the adhesion of the film to the substrate. This one changes significantly the dispersive behavior of the surface acoustic waves. To generate and detect these waves, a Laser-Ultrasonics setup has been used. First, we are looking to bypass the interpretation difficulties usually encountered in the control of adhesion by surface acoustic waves. Indeed, the layer thickness variations influence the dispersion of the waves in a similar way to the adhesion. Consequently, the polymer films used have a quasi-constant thickness and they are deposited directly on an aluminum substrate. In addition, these films are also transparent. It allows us to generate directly the acoustic waves on the substrate surface, at the interface between the film and the substrate, by focusing the laser pulse through the film. In this way, the influence of the source location on the dispersive behavior of the surface acoustic waves and thus on the adhesion quality control may be studied experimentally and by using finite element simulations. Finally, a characterization of the adhesion of several samples is performed using the dispersion curves obtained applying the Matrix-Pencil method to the experimental results. An inversion algorithm allows us to estimate the interfacial stiffnesses corresponding to the adhesion of the samples.

Cette thèse a pour objectif de développer des méthodes ultrasonores pour l’END de l’adhésion de joints collés. Pour aborder ce problème, les assemblages sont réalisés avec des substrats en aluminium (isotrope, élastique) et un adhésif de type époxy(isotrope, viscoélastique). Selon la géométrie des assemblages, deux méthodes sont proposées pour obtenir une information quantitative sur le niveau d’adhésion. La première est adaptée à des échantillons de type tri-couches avec recouvrement total.Elle consiste à analyser la transmission d’ondes ultrasonores planes de volume à travers l’assemblage immergé dans l'eau.Les conséquences d'un mauvais traitement de surface des substrats sur la mesure des modules de viscoélasticité du joint adhésif sont étudiées. Il a été montré que des interphases dégradées provoquent une anisotropie apparente des modules mesurés pour le joint de colle. Cette anisotropie a été quantifiée à l'aide de deux paramètres β1 et β2 dont les valeurs permettent de révéler la qualité des interphases. Ensuite, les modules élastiques (ou raideurs kL et kT) des interphases ont été estimés, en supposant les propriétés de l'adhésif connues. On montre que leurs valeurs sont maximales lorsque l'adhésion est nominale, et qu’elles diminuent franchement lorsque l'adhésion est dégradée, mais dans des proportions différentes. Des mesures de la contrainte à rupture, réalisées sur des échantillons préparés dans les mêmes conditions, corroborent la chute des modules des interphases. La seconde méthode est adaptée à des échantillons collés avec un recouvrement partiel. Elle es tbasée sur la mesure du coefficient de transmission d'ondes de Lamb se propageant d'un substrat à l'autre, à travers la zone de recouvrement. Une étude de sensibilité numérique (par éléments finis) des coefficients de transmission des modes de Lamb a montré que les propriétés mécaniques des interphases (modélisées par des raideurs surfaciques) peuvent être évaluées si les autres caractéristiques de l’assemblage sont connues. Des mesures expérimentales de ces coefficients de transmission ont ensuite été réalisées avec deux échantillons. L’un d’eux possède des interphases à adhésion nominale et l’autre des interphases dégradées. Une confrontation entre les résultats des mesures obtenus pour les différents modes et les simulations numériques permet de déterminer les valeurs des raideurs d’interfaces pour chaque échantillon. Là encore, il est observé qu’une mauvaise adhésion se traduit par des valeurs faibles des raideurs d’interfaces, qui peuvent être quantifiées, cette fois,grâce aux ondes ultrasonores guidées. / The aim of this thesis is to develop NDT ultrasonic methods for the evaluation of the adhesion at interfaces in bonded joints.To address this problem, the assemblies are made of aluminum substrates (isotropic, elastic) and epoxy adhesive (isotropic,viscoelastic). Depending on the geometry of the studied assemblies, two methods are proposed to obtain quantitativeinformation on the adhesion level. The first method is suitable for three-layered plate-like samples. It consists in analyzingthe transmission of bulk ultrasonic plane waves through the assembly immersed in water. The consequences of a inapropriatesurface treatment of the substrates on the measurement of viscoelastic moduli of the joint are studied. It has been shown thatdegraded interphases cause an apparent anisotropy of the measured elastic moduli of the joint. This anisotropy was quantifiedusing two parameters β1 and β2 whose values can reveal the quality of the interphases. Then, the elastic moduli (or stiffnesseskL and kT) of the interphases were estimated, assuming that the adhesive layer properties are known. It was shown that theirvalues are higher when the adhesion is nominal, and are strongly decreased when the adhesion is degraded. Measurements ofthe mechanical strength, on samples prepared under the same conditions, confirmed the drops in these elastic moduli. Thesecond method is more suitable for lap joint samples. It is based on the Lamb wave transmission coefficient measured fromone substrate to the other, across the overlap zone. A numerical sensitivity analysis (finite element model) of the transmissioncoefficients of Lamb waves showed that the mechanical properties of the interphases (modeled by interfacial stiffnesses) canbe evaluated if the other characteristics of the assembly are well known. Experimental measurements of the transmissioncoefficients were then performed on two samples. One of them has interphases with nominal adhesion and the other hasdegraded interphases. A comparison between the measured results obtained for the different modes and numericalsimulations was used to determine the values of the interfacial stiffnesses for each sample. Again, it was observed that pooradhesion leads to low values of the interfacial stiffnesses of the interphases, that can be quantified using guided ultrasonicwaves.

La tomographie en contraste de diffraction est une nouvelle technique non destructive d'imagerie synchrotron qui caractérise la microstructure et l'orientation des grains dans les matériaux polycristallins en trois dimensions (3D). En la combinant avec la tomographe par contraste de phase. Il est pour la première fois possible d'observer in situ la propagation 3D des fissures de fatigue courtes au sein d'un ensemble de grains entièrement caractérisé (orientation et forme). L'approche combinée, appelée « tomographie tri-dimensionnelle par rayons X des fissures courtes et de la microstructure »(T3DXFM), a été développée sur l’alliage de titane métastable "Beta21S". Une grande partie de ce travail porte sur le développement de la méthodologie T3DXFM. Dans le jeu de données combinées, chaque point de la surface de rupture 3D peut être associé à une structure de données multidimensionnelle contenant des variables décrivant l'orientation des grains, l'orientation locale de la surface de rupture ainsi que l'histoire de la propagation. La méthode utilise un maillage de surface composé de triangles qui décrit la fissure (en d'autres termes: la surface de rupture) dans l'état de propagation mesuré au dernier cycle de fatigue réalisé. Les orientations des grains, les différents fronts de la fissure, les vitesses de croissance locales ainsi que les joints de grains peuvent être visualisés en attribuant des couleurs à ce maillage. Des outils d'extraction des figures de pôle ont été créés et mis en œuvre. Un algorithme a été développé qui est capable de mesurer la vitesse de propagation locale 30 d'une fissure contenant des branchements. / X-Ray Diffraction Contrast Tomography (DCT) is a recently developed, non-destructive synchrotron imaging technique which characterizes microstructure and grain orientation in polycrystalline materials in three dimensions (3D). By combining it with propagation based phase contrast tomography (PCT) it is for the first lime possible to observe in situ the 3D propagation behavior of short fatigue cracks (SFCs) within a set of fully characterized grains (orientation and shape). The combined approach, termed 3D X-ray Tomography of short cracks and Microstructure (3DXTSM), has been developed on the metastable beta titanium alloy "Beta21S". A large part of this work deals with the development of the 3DXTSM methodology. In the combined dataset, each point on the 3D fracture surface can be associated with a multidimensional data structure containing variables describing the grain orientation, the local fracture surface normal and the propagation history. The method uses a surface mesh composed of triangles that describes the crack (in other words: the fracture surface) in the last propagation state measured. Grain orientations, crack fronts, local growth rates and grain boundaries can be visualized by assigning colors to this mesh. The data structure can be interrogated in a number of different ways. Tools for extracting pole figures and pole density distribution functions have been implemented. An algorithm was developed that is capable of measuring the 3D local growth rate of a crack containing branches. The accuracy of the grain boundaries as reconstructed with OCT was evaluated and the elastic constants of Beta21S were determined.