This DPhil project is set in the context of single crystal elasticity-plasticity finite element modelling. Its core objective was to develop and implement a methodology for predicting the evolution of texture in single and dual-phase material systems. This core objective has been successfully achieved. Modelling texture evolution entails essentially modelling large deformations (as accurately as possible) and taking account of the deformation mechanisms that cause texture to change. The most important deformation mechanisms are slip and twinning. Slip has been modelled in this project and care has been taken to explore conditions where it is the dominant deformation mechanism for the materials studied. Modelling slip demands that one also models dislocations since slip is assumed to occur by the movement of dislocations. In this project a model for geometrically necessary dislocations has been developed and validated against experimental measurements. A texture homogenisation technique which relies on interpretation of EBSD data in order to allocate orientation frequencies based on representative area fractions has been developed. This has been coupled with a polycrystal plasticity RVE framework allowing for arbitrarily sized RVEs and corresponding allocation of crystallographic orientation. This has enabled input of experimentally measured initial textures into the CPFE model allowing for comparison of predictions against measured post-deformation textures, with good agreement obtained. The effect of texture on polycrystal physical properties has also been studied. It has been confirmed that texture indeed has a significant role in determining the average physical properties of a polycrystal. The thesis contributes to the following areas of micro-mechanics materials research: (i) 3D small deformation crystal plasticity finite element (CPFE) modelling, (ii) geometrically necessary dislocation modelling, (iii) 3D large deformation CPFE modelling, (iv) texture homogenisation methods, (v) single and dual phase texture evolution modelling, (vi) prediction of polycrystal physical properties, (vii) systematic calibration of the power law for slip based on experimental data, and (viii) texture analysis software development (pole figures and Kearns factors).

Les ultrasons laser font partie des techniques expérimentales d'inspection en développement pour des problématiques de contrôle non destructif toujours plus exigeantes. Ces méthodes sans contact permettent dans notre cas de générer et détecter des ondes de surface se propageant sur des distances de quelques dizaines de microns.L'objectif de cette thèse est de mettre au point et tester une méthode ultrasons lasers pour l'étude d'échantillons d'aciers polycristallins. En effet ces matériaux à la base d'un ensemble de pièces manufacturées sont composés d'une multitudes de cristaux anisotropes de tailles et d'orientations variables. Parmi les paramètres physiques d'intérêt l'évaluation des constantes élastiques ainsi que la détermination de texturations sont des facteurs essentiels dans toute application CND et leur simulation numérique. Dans un premier temps, les variations de vitesses de propagation dues à l'anisotropie et aux l'orientation différentes des cristaux sont mesurées grâce au montage expérimental en surface d'échantillons d'acier austénitique. Les résultats sont exploités sous la forme de distributions statistiques de vitesses de propagation et un protocole d'inversion est proposé pour l'évaluation des constantes élastiques. Dans une seconde partie, un échantillon d'acier moulé présentant de larges macrograins est analysé pour la mise en évidence de sous-structures au sein d'un milieu biphasique.Finalement un échantillon de soudure possédant une texturation typique avec grains allongés est étudié et la présence d'une orientation préférentielle est soulignée. / Laser ultrasonics are among the most interesting experimental technics to answer the need of industrials in NDT applications. These non contact methods allow to generate and detect elastic surface wave propagating on typical distances of several tens of microns.The main objective of this thesis is to test laser ultrasonics methods for polycristalline steel study. These materials are composed of several anisotropic crystals with differents sizes and orientations. Among the essential parameters to describe such material, the determination of elastic parameters and texturations are key to good understanting and for the modelisation of material in NDT studies.In a first part, velocities variations due to elastic anisotropy and crystalline orientation of grains is assessed by measuring surface acoustic wave velocity. The experimental results are used in the form of statistical distribution of velocity combined with a numerical model to evaluate elastic parameters.Stainless casted steel is a complex material presenting large biphasic grains. They are here studied to show the presence of sub-structures.Finally, a weld sample presenting typical structuration of ellongated grains is studied to show the presence of preferential orientation.

Dans les réacteurs nucléaires à eau sous-pression, le combustible est constitué de pastilles d’oxyde d’uranium (UO2), dont le comportement ne peut être simulé qu'à travers une modélisation multi-échelles et multi-physiques, tenant compte à la fois de la thermo-mécanique et de la physico-chimie relative aux produits de fission. L’évolution récente des modèles et des moyens de calcul a permis de développer les simulations à l’échelle de la microstructure et d’accroitre les possibilités de couplage. Ce travail concerne le développement d'un modèle de comportement thermo-mécanique de l’UO2 à l’échelle du polycristal. Le comportement du VER est analysé en termes de réponse effective et de phénomènes de localisation. Nous nous intéressons notamment aux valeurs de pression hydrostatique, qui pilotent le transport des produits de fission. La robustesse des résultats obtenus en fonction du choix du maillage éléments finis est étudiée. Une série de calculs est présentée afin de trouver un compromis satisfaisant en termes de discrétisation pour une estimation correcte des contraintes locales. Une première étude propose de retrouver des mesures expérimentales de dé cohésion intergranulaire sur le combustible en introduisant des modèles de zones cohésives dans le VER. Afin d'analyser le comportement micromécanique de l’UO2 en irradiation, un chargement de type rampe de puissance est appliqué au polycristal. L’analyse des distributions locales de contraintes donne lieu à une discussion sur l’effet de l’incompatibilité intergranulaire sur le transport des produits de fission. / In Pressurized Water Reactors (PWR), fuel is generally composed of uranium dioxide pellets piled up within a zirconium tubular cladding. Modeling of the fuel behavior in nominal and accidental conditions requires multi-scale models in order to take into account both the thermo-mechanical behavior of the pellets and the effects of fission gases. Recent development of micromecanical tools of simulation has improved coupling possibilities. Our study proposes to set a full micromechanical model for uranium dioxide dealing with both mechanics and fission products transport at the scale of a polycristalline aggregate. Both the effective behavior of the RVE and stress localization effects are studied. Hydrostatic pressure, which directly controls the diffusion of fission gases, is given a particular focus. The numerical robustness of results is also debated in terms of mesh refinement. A series of simulations leads to a satisfying compromise between accuracy and calculation time. A study compares experimental measurement of intergranular crack opening and simulation results obtained using cohesive models. The micromecanical behavior of uranium dioxide during irradiation is analysed by submitting the polycristalline RVE to transient irradiation. The local stress distribution leads to a debate on the consequences of intergranular strain incompatibility on fission gases diffusion.

Im Legierungssytem Ni-Mn-Ga tritt bei Zusammensetzungen nahe der stöchiometrischen Zusammensetzung Ni2MnGa der magnetische Formgedächtniseffekt auf. Darunter versteht man die Dehnung durch Bewegung von Zwillingsgrenzen im Magnetfeld. Einkristalle aus Ni-Mn-Ga mit einer tetragonalen 5M-Martensitstruktur zeigen magnetisch erzeugbare Dehnungen von bis zu 6 %. Diese großen Dehnungen verbunden mit der schnellen Schaltfrequenz von Magnetfeldern machen den Effekt interessant für technische Anwendungen z. B. als Aktoren. Derartige Einkristalle sind schwierig und teuer herzustellen, weshalb für technische Anwendung Polykristalle von Interesse sind. Diese lassen sich im Allgemeinen leichter und preiswerter herstellen. Um den magnetischen Formgedächtniseffekt in Polykristalle einzustellen, werden grobkörnige, texturierte Proben mittels des Verfahrens der gerichteten Erstarrung hergestellt. Die Gefügeuntersuchungen erfolgen mit metallographischen Schliffen und die Kornorientierungen werden mit der EBSD-Technik bestimmt. Um das Gefüge zu vergröbern, werden Glühungen nach einer aufgebrachten Warmverformung untersucht. Zur Verringerung der für die Bewegung der Zwillingsgrenzen notwendigen Spannung (Zwillingsspannung) werden die Proben im Druckversuch mechanisch trainiert. Dabei kann die Zwillingsspannung teilweise unter die magnetisch erzeugbare Spannung auf die Zwillingsgrenzen (Magnetospannung) abgesenkt werden. Eine weitere Absenkung der Zwillingsspannung wird durch eine plattenförmige Probengeometrie mit Dicken im Bereich der Korndurchmesser erreicht. An derartigen Proben wird magnetisch rückstellbare freie Dehnung durch Zwillingsgrenzenbewegung erzielt.

Dans ce travail est proposée une méthodologie générale de micromécanique expérimentale multi-échelle des polycristaux. Elle a été appliquée dans le cas d'un polycristal de sel gemme, qui en plus d'avoir des applications industrielles de stockage d'énergie et de déchets, constitue un matériau modèle de micromécanique présentant une déformation plastique aussi bien à l'ambiante qu'à haute température. La déformation ductile à l'échelle de la microstructure opère par la plasticité cristalline intra-granulaire traditionnelle, mais aussi des mécanismes de déformation inter-granulaires, tels que le glissement aux joints de grains. Nous avons dans un premier temps quantifié précisément la part de chacun de ces mécanismes locaux dans la déformation macroscopique du sel en se basant sur la technique de corrélation d'images numériques (CIN), obtenues au cours d'un essai de compression uni-axiale in-situ dans la chambre d'un microscope électronique à balayage (MEB). Afin d'augmenter la précision de cette quantification, des motifs spéciaux gravés aux interfaces des grains par micro-lithograhie ont été proposés. Ensuite, les observations surfaciques (par MEB) ont été étendues au cœur du matériau grâce à la micro-tomographie à rayons X et à la technique de corrélation d'images volumiques (CIV). Pour ce faire, des particules micrométriques de cuivre (3 % en volume) ont été dispersées dans le matériau lors de son élaboration, afin d'avoir un marquage local volumique adapté pour la CIV. Différentes microstructures (en termes de taille moyenne de grain) ont été considérées. De nouvelles procédures de CIV ont permis d'accéder à la répartition tridimensionnelle de la déformation ductile à l'échelle de la microstructure polycristalline avec une précision inferieure à la taille moyenne de grain. Les mécanismes de déformation observés à cœur d'échantillon sous chargement uni-axial sont cohérents avec ceux identifiés par les observations surfaciques. L'importance des mécanismes inter-granulaires dans la déformation ductile et dans l'endommagement diffus du sel a été confirmée. Une caractérisation tridimensionnelle de la microstructure par DCT (Diffraction Contrast Tomography) a été effectuée et comparée à des mesures surfaciques d'orientation cristalline par EBSD (Electron BackScattered Diffraction). Enfin, la comparaison des champs de déformation surfacique et volumique obtenus sur les mêmes échantillons a permis de retrouver les mêmes organisations et développements des localisations de déformation ductile en surface et en volume, et de les relier aux conditions de chargement et à la microstructure / We develop in this study new experimental methodologies for the multi-scale experimental investigation of the micromechanics of polycrystalline materials. These methodologies are applied to synthetic halite (NaCl), which is a convenient model polycristal due to its viscoplastic behavior at both ambient and high temperatures (350°C). In addition, halite is used for industrial applications such as underground energy and waste storage. The ductile deformation at the scale of the microstructure operates not only through conventional intra-granular plasticity, but also through inter-granular deformation mechanisms, such as grain-boundary sliding (GBS). First, we precisely quantify the relative contribution of each of these local mechanisms to the macroscopic deformation of halite. For this purpose, we apply digital image correlation (DIC) technique to high resolution images obtained during uniaxial compression tests in the chamber of a scanning electron microscope (SEM). The DIC algorithms have been modified to account for the discontinuous kinematics at grain boundries. We also propose a method to improve accuracy of GBS quantification, which consists in creating specific artificial patterns across grain-boundaries by electron beam lithography. The results show that GBS is present from the beginning of plastic deformation of the polycrystal. The 2D observations (using SEM) are complemented by 3D volume investigations using X-ray computed microtomography and Digital Volume Correlation (DVC) techniques. In order to obtain local volume markers differing in contrast (density) from NaCl and adapted to DVC, micrometric copper particles (3 % in volume) are dispersed into the material during its elaboration. Various microstructures (in terms of average grain size) are considered. New DVC protocols allow us to obtain the three-dimensional distribution of ductile deformation at the scale of the polycrystalline microstructure, with a spatial resolution finer than the average grain size. 3D and 2D local mechanical fields are compared on the same samples submitted to uniaxial compression. The strain patterns and the deformation mechanisms observed in depth of the sample are consistent with those identified by 2D observations. The results show the same organization and development of strain localization bands in relation with the loading conditions and microstructure, both at the surface and in volume. The importance of inter-granular mechanisms for the plastic deformation and diffuse damage of halite is also confirmed in 3D. Finally, in view of a further numerical model of the plasticity of the polycrystal, the three-dimensional polycrystalline microstructure is characterized by diffraction contrast tomography and compared to 2D measurements obtained by electron BackScattered diffraction

La texture marquée des tôles et des gaines en alliages de zirconium se traduit par une forte anisotropie du comportement mécanique. L'objectif de l'étude est de caractériser et de modéliser le comportement anisotrope de tôles en alliage de Zircaloy-4 recristallisé. La caractérisation de l'anisotropie du comportement est réalisée au travers d'essais mécaniques conventionnels (chargements en traction et en cisaillement) sur tôles en utilisant la méthode de corrélation d'images numériques. Dans un premier temps, un modèle a été identifié à partir de cette base expérimentale sur tôle. La loi est validée par des calculs éléments finis d'essais de traction sur éprouvettes plates entaillées. Dans un second temps, la transférabilité du modèle de la tôle vers le tube a été étudiée. Pour les chargements uniaxiés, la transférabilité est avérée. En revanche, pour les chargements biaxiés la transférabilité est moins bonne. Une réidentification des paramètres gérant l'anisotropie du comportement en intégrant à la base d'identification un essai équibiaxié sur tube a été nécessaire. / Zirconium alloy sheet and clad are strongly textured materials, resulting in sharp anisotropic mecanical behavior. The purpose of this work is to characterize and model the anisotropic behavior of recrystallized Zircaloy-alloy sheets. Anisotropy is investigated by usual mechanical tests (tensile and shear loadings) performed on sheet material using digital image correlation measurments. A model is identified from this experimental database obtained on sheet material. Finite element simulations of tensile notched tests are used to validate the law. Then, the model transferability from the sheet to the rod is studied. The transferability is suitable for uniaxial loading. On the contrary, the transferability is not fully adequate for biaxial loadings. Therefore, a new identification of parameters dealing with anisotropy from enriched database with an equibiaxial rod test is necessary.

La propagation des fissures courtes de fatigue dans un matériau polycristallin dépend fortement de la microstructure. Bien que de nombreuses études de caractérisation et de modélisation existent sur le sujet, la prédiction du chemin et de la vitesse de propagation de ce type de fissure n'est pas encore possible aujourd'hui.Afin de bien comprendre les mécanismes de propagation, la caractérisation in-situ d'un échantillon par la tomographie aux rayons X a été réalisée à l'ESRF en combinant deux techniques de caractérisation. La tomographie par Contraste de Diffraction (DCT) qui est une méthode non destructive permettant de caractériser en 3D la morphologie et l'orientation des grains constitutifs de la microstructure, à l'état non-déformé, et la tomographie par Contraste de Phase (PCT) qui permet d'obtenir la forme de fissure à divers étapes de la vie de l'éprouvette. Grâce à ces informations, il est possible de simuler la propagation de fissure en utilisant un maillage réaliste reconstruit à partir des images tomographiques. Dans ce travail, une étude de l'anisotropie de comportement élastique est effectuée dans un maillage microstructural 3D reconstruit à partir des images tomographiques. Cette étude permet de comparer les tenseurs de déformation élastique moyennés à chaque grain avec les mesures expérimentale. Ensuite, une nouvelle méthodologie est proposée pour simuler la propagation de fissure. Issue d'une simulation en plasticité cristalline, la direction et la vitesse de la propagation de fissure est déterminée par un post-traitement, ce qui permet de propager la fissure par remaillage. Cette méthode est appliquée dans un premier temps à un monocristal pré-fissuré pour prédire le trajet de fissuration en fonction des systèmes de glissement activés. L'ensemble de la démarche est enfin appliqué au polycristal complet imagé par tomographie. Le rôle du joint de grains et la vitesse de propagation sont également analysés. En comparant les résultats de simulation avec les mesures expérimentales, le critère de la propagation de fissure est discuté. / The short fatigue crack propagation in polycrystal materials depends strongly on microstructure. Although numerous studies of characterisation and of simulation, the prediction of the short fatigue crack propagation remains a challenge.In order to understand the mechanisms of short fatigue crack propagation, an in-situ characterisation by X-ray tomography was carried out at ESRF, using two techniques of tomography. Diffraction Contrast Tomography (DCT) that is a non-destructive method can be used to obtain 3D morphology and grain orientations in an undeformed state of polycrystal materials. Couple with Phase Contrast Tomography (PCT), it allows to characterise the short fatigue crack propagation at different loading stages. Access to this information, it is possible to simulate the short fatigue crack propagation using a 3D reel microstructural mesh reconstructed from the tomographic images.In this work, the elastic anisotropic behaviour in a 3D microstructural mesh is performed. The elastic strain tensors averaged in grains are also compared to the experimental measurements. Then, a new numerical approach is proposed to simulate crack propagation. From a crystal plasticity FE simulation, the crack growth direction is determined by a post processing. Next, the crack is propagated through remeshing. This approach is firstly applied to the single crystals, then to the polycrystal mesh reconstructed from the tomographic images. The grain boundary effects and the crack growth rate are also analysed. By comparing between simulation and experimental crack, the damage indicator is discussed at the end.

Dans ce travail est proposée une méthodologie générale de micromécanique expérimentale multi-échelle des polycristaux. Elle a été appliquée dans le cas d'un polycristal de sel gemme, qui en plus d'avoir des applications industrielles de stockage d'énergie et de déchets, constitue un matériau modèle de micromécanique présentant une déformation plastique aussi bien à l'ambiante qu'à haute température. La déformation ductile à l'échelle de la microstructure opère par la plasticité cristalline intra-granulaire traditionnelle, mais aussi des mécanismes de déformation inter-granulaires, tels que le glissement aux joints de grains. Nous avons dans un premier temps quantifié précisément la part de chacun de ces mécanismes locaux dans la déformation macroscopique du sel en se basant sur la technique de corrélation d'images numériques (CIN), obtenues au cours d'un essai de compression uni-axiale in-situ dans la chambre d'un microscope électronique à balayage (MEB). Afin d'augmenter la précision de cette quantification, des motifs spéciaux gravés aux interfaces des grains par micro-lithograhie ont été proposés. Ensuite, les observations surfaciques (par MEB) ont été étendues au cœur du matériau grâce à la micro-tomographie à rayons X et à la technique de corrélation d'images volumiques (CIV). Pour ce faire, des particules micrométriques de cuivre (3 % en volume) ont été dispersées dans le matériau lors de son élaboration, afin d'avoir un marquage local volumique adapté pour la CIV. Différentes microstructures (en termes de taille moyenne de grain) ont été considérées. De nouvelles procédures de CIV ont permis d'accéder à la répartition tridimensionnelle de la déformation ductile à l'échelle de la microstructure polycristalline avec une précision inferieure à la taille moyenne de grain. Les mécanismes de déformation observés à cœur d'échantillon sous chargement uni-axial sont cohérents avec ceux identifiés par les observations surfaciques. L'importance des mécanismes inter-granulaires dans la déformation ductile et dans l'endommagement diffus du sel a été confirmée. Une caractérisation tridimensionnelle de la microstructure par DCT (Diffraction Contrast Tomography) a été effectuée et comparée à des mesures surfaciques d'orientation cristalline par EBSD (Electron BackScattered Diffraction). Enfin, la comparaison des champs de déformation surfacique et volumique obtenus sur les mêmes échantillons a permis de retrouver les mêmes organisations et développements des localisations de déformation ductile en surface et en volume, et de les relier aux conditions de chargement et à la microstructure / We develop in this study new experimental methodologies for the multi-scale experimental investigation of the micromechanics of polycrystalline materials. These methodologies are applied to synthetic halite (NaCl), which is a convenient model polycristal due to its viscoplastic behavior at both ambient and high temperatures (350°C). In addition, halite is used for industrial applications such as underground energy and waste storage. The ductile deformation at the scale of the microstructure operates not only through conventional intra-granular plasticity, but also through inter-granular deformation mechanisms, such as grain-boundary sliding (GBS). First, we precisely quantify the relative contribution of each of these local mechanisms to the macroscopic deformation of halite. For this purpose, we apply digital image correlation (DIC) technique to high resolution images obtained during uniaxial compression tests in the chamber of a scanning electron microscope (SEM). The DIC algorithms have been modified to account for the discontinuous kinematics at grain boundries. We also propose a method to improve accuracy of GBS quantification, which consists in creating specific artificial patterns across grain-boundaries by electron beam lithography. The results show that GBS is present from the beginning of plastic deformation of the polycrystal. The 2D observations (using SEM) are complemented by 3D volume investigations using X-ray computed microtomography and Digital Volume Correlation (DVC) techniques. In order to obtain local volume markers differing in contrast (density) from NaCl and adapted to DVC, micrometric copper particles (3 % in volume) are dispersed into the material during its elaboration. Various microstructures (in terms of average grain size) are considered. New DVC protocols allow us to obtain the three-dimensional distribution of ductile deformation at the scale of the polycrystalline microstructure, with a spatial resolution finer than the average grain size. 3D and 2D local mechanical fields are compared on the same samples submitted to uniaxial compression. The strain patterns and the deformation mechanisms observed in depth of the sample are consistent with those identified by 2D observations. The results show the same organization and development of strain localization bands in relation with the loading conditions and microstructure, both at the surface and in volume. The importance of inter-granular mechanisms for the plastic deformation and diffuse damage of halite is also confirmed in 3D. Finally, in view of a further numerical model of the plasticity of the polycrystal, the three-dimensional polycrystalline microstructure is characterized by diffraction contrast tomography and compared to 2D measurements obtained by electron BackScattered diffraction

Dans ce mémoire, des méthodes spectrales basées sur la transformée de Fourier rapide ("fast Fourier transform" en anglais notée "FFT") sont développées pour résoudre les équations de champs et d’évolution des densités de dislocations polarisées ou géométriquement nécessaires dans la théorie de la mécanique des champs de dislocations ("Field Dislocations Mechanics" en anglais et notée "FDM") et de son extension phénoménologique et mésoscopique ("Phenomenological Mesoscopic Field Dislocations Mechanics" en anglais et notée "PMFDM"). Dans un premier temps, une approche spectrale a été développée pour résoudre les équations élasto-statiques de la FDM pour la détermination des champs mécaniques locaux provenant des densités de dislocations polarisées et des hétérogénéités élastiques présentes dans les matériaux de microstructure supposée périodique et au comportement élastique linéaire. Les champs élastiques sont calculés de façon précise et sans oscillation numérique même lorsque les densités de dislocations sont concentrées sur un seul pixel (pour les problèmes à deux dimensions) ou sur un seul voxel (pour les problèmes à trois dimensions). Ces résultats sont obtenus grâce à l’application de formules de différenciation spatiale pour les dérivées premières et secondes dans l’espace de Fourier basées sur des schémas à différences finies combinées à la transformée de Fourier discrète. Les résultats obtenus portent sur la détermination précise des champs élastiques des dislocations individuelles de types vis et coin, et des champs élastiques d’interaction entre des inclusions de géométries variées et différentes distributions de densités de dislocations telles que les dipôles ou les boucles de dislocations dans un matériau composite biphasé et des microstructures tridimensionnelles. Dans un second temps, une approche spectrale a été développée pour résoudre de façon rapide et stable l’équation d’évolution spatio-temporelle des densités de dislocations dans la théorie FDM. Cette équation aux dérivées partielles, de nature hyperbolique, requiert une méthode spectrale avec des filtres passe-bas afin de contrôler à la fois les fortes oscillations inhérentes aux approches FFT et les instabilités numériques liées à la nature hyperbolique de l’équation de transport. La validation de cette approche a été effectuée par des comparaisons avec les solutions exactes et les méthodes éléments finis dans le cadre de la simulation des phénomènes physiques d’annihilation ou d’extension/annihilation de boucles de dislocations. En dernier lieu, une technique numérique pour la résolution des équations de la PMFDM est développée dans le cadre d’une formulation FFT pour un comportement élasto-visco-plastique avec la prise en compte de la contribution des dislocations géométriquement nécessaires et statistiquement stockées ainsi que des conditions de saut de la distorsion plastique aux interfaces de type joint de grains ou joint de phases. Cette technique est par la suite appliquée à la simulation de la déformation plastique de structures modèles telles que des microstructures périodiques à canaux et des polycristaux métalliques / Fast Fourier transform (FFT)-based methods are developed to solve both the elasto-static equations of the Field Dislocation Mechanics (FDM) theory and the dislocation density transport equation of polarized or geometrically necessary dislocation (GND) densities for FDM and its mesoscopic extension, i.e. the Phenomenological Mesoscopic Field Dislocations Mechanics (PMFDM). First, a numerical spectral approach is developed to solve the elasto-static FDM equations in periodic media for the determination of local mechanical fields arising from the presence of both polarized dislocation densities and elastic heterogeneities for linear elastic materials. The elastic fields are calculated in an accurate fashion and without numerical oscillation, even when the dislocation density is restricted to a single pixel (for two-dimensional problems) or a single voxel (for three-dimensional problems). These results are obtained by applying the differentiation rules for first and second derivatives based on finite difference schemes together with the discrete Fourier transform. The results show that the calculated elastic fields with the present spectral method are accurate for different cases considering individual screw and edge dislocations, the interactions between inhomogeneities of various geometries/elastic properties and different distributions of dislocation densities (dislocation dipoles, polygonal loops in two-phase composite materials). Second, a numerical spectral approach is developed to solve in a fast, stable and accurate fashion, the hyperbolic-type dislocation density transport equation governing the spatial-temporal evolution of dislocations in the FDM theory. Low-pass spectral filters are employed to control both the high frequency oscillations inherent to the Fourier method and the fast-growing numerical instabilities resulting from the hyperbolic nature of the equation. The method is assessed with numerical comparisons with exact solutions and finite element simulations in the case of the simulation of annihilation of dislocation dipoles and the expansion/annihilation of dislocation loops. Finally, a numerical technique for solving the PMFDM equations in a crystal plasticity elasto-viscoplastic FFT formulation is proposed by taking into account both the time evolutions of GND and SSD (statistically stored dislocations) densities as well as the jump condition for plastic distortion at material discontinuity interfaces such as grain or phase boundaries. Then, this numerical technique is applied to the simulation of the plastic deformation of model microstructures like channel-type two-phase composite materials and of polycrystalline metals

<p dir="ltr">We investigate the crystal plasticity and shock response of single crystal and polycrystalline pentaerythritol tetranitrate (PETN) using mesoscale finite element simulations. The model includes the Mie-Grüneisen Equation of State and a single crystal plasticity model. Simulations with single crystals with different orientations are tested using our plasticity model under shock compression to explore shear stress and slip. Parameters regarding the Mie-Grüneisen Equation of State are also verified in various orientations from 0.50 to 1.75 km/s. A polycrystalline PETN sample with varying grain sizes and orientations are subjected to shock loading with impact velocities ranging from 0.25 to 0.75 km/s. We study how differences in shock orientation affect slip and stress in PETN at different shock strengths.</p>