Cette thèse a pour but le développement d'un modèle constitutif hydromécanique prenant en compte le couplage de l'endommagement et de la plasticité dans les géomatériaux non-saturés. Le couplage de ce modèle avec les transferts hydrauliques est rendu possible par son implémentation dans le code aux éléments finis $Theta$-Stock. Un algorithme local a été développé spécifiquement dans ce but. Des problèmes hydromécaniques complètement couplés, tels que la création de la zone endommagée par excavation autour d'un tunnel ou la création d'un endommagement de dessiccation et d'humidification ont été simulés. Une contrainte doublement effective incorporant les effets de la succion et de l'endommagement (considéré isotrope) a été définie en s'appuyant sur des bases thermodynamiques. Cette approche a l'avantage de considérer une unique variable de contrainte étant thermodynamiquement conjuguée aux déformations élastiques. Une formulation hyperélastique dépendante de la pression de confinement est utilisée pour décrire le comportement dans le domaine élastique. L'évolution des rigidités élastiques avec l'endommagement est comparée en considérant deux hypothèses : le principe des déformations équivalentes et le principe de l'énergie élastique équivalente. L'hypothèse d'équivalence des déformations permet d'introduire la contrainte doublement effective dans les équations de plasticité et ainsi de coupler plasticité et endommagement. Les équations de plasticité sont dérivées du modèle de Barcelone (Alonso extit{et al.}, 1990), lui même basé sur le modèle de Cam-Clay. Deux critères distincts d'endommagement et de plasticité sont définis qui peuvent être activés aussi bien indépendamment que simultanément. Les surfaces de charges étant exprimées en fonction de la contrainte effective et de la succion, elles évoluent automatiquement en fonction de la succion et de l'endommagement dans l'espace des contraintes totales. Cela permet de représenter la transition d'un comportement ductile vers un comportement fragile lorsque le matériau est asséché. Un algorithme local explicite a été développé pour gérer le couplage des deux phénomènes dissipatifs. L'implémentation du modèle constitutif dans le code aux éléments finis $Theta$-Stock permet l'étude de problèmes hydromécaniques complètement couplés, les lois de transfert dépendant elles aussi de l'état de saturation du matériau. Ce modèle a d'abord été appliqué à la simulation du développement de micro-fissures lors du séchage d'une éprouvette. L'apparition de l'endommagement est expliquée principalement par le gradient de pression très important créé à la surface de l'échantillon lorsqu'on applique une forte variation de succion. Enfin, un problème à échelle réelle est simulé. L'excavation d'un tunnel, la désaturation du sol environnant dû à la ventilation, ainsi que la création de la zone endommagée par excavation sont étudiés. L'étendue des zones l'endommagement et de déformations plastiques autour du tunnel est étudiée

Lors de séismes forts, la propagation des ondes sismiques dans les sols met en jeu des non linéarités de comportement qui se manifestent différemment selon le niveau de sollicitation. En effet, pour de faibles déformations (généralement <10^{-6}), une loi de comportement linéaire (i.e. module et amortissement indépendants du niveau de sollicitation) permet de reproduire les observations expérimentales sur site. En revanche, pour des déformations plus élevées, une loi de comportement non linéaire hystérétique est nécessaire pour décrire l'évolution de la rigidité et des dissipations énergétiques au cours de la sollicitation sismique. De plus, comme les séismes forts sont caractérisés par des amplitudes et des durées plus importantes, le rôle de la pression interstitielle ne peut pas être négligé pour les sols saturés sous fortes sollicitations (mobilité cyclique et liquéfaction). Ces phénomènes peuvent conduire à l'annulation des contraintes effectives et devenir cause de dommages sévères pour les structures et les ouvrages. L'analyse proposée élargit l'applicabilité des modèles de calcul actuels pour une analyse plus fine du risque sismique. En partant d'une formulation aux éléments finis décrivant la propagation des ondes sismiques suivant la direction verticale en prenant en compte le chargement en 3D (l'approche " 1D-3C" : une direction-trois composantes) dans les sols nonlinéaires secs, de nouvelles stratégies pour la prise en compte du rôle de l'eau sont développées. Le modèle de comportement est basé sur la relation entre la pression interstitielle et le travail de la contrainte de cisaillement. Ce modèle décrivant l'évolution de la pression interstitielle considère l'état de contrainte tridimensionnel du matériau. Le modèle est validé par comparaison avec des résultats expérimentaux. L'approche " 1D-3C " a été utilisée pour modéliser la réponse des sols pour 4 séismes réels : le séisme de Superstition Hills en 1987 aux États-Unis (M_{w}=6.7); le séisme de Tohoku en 2011 au Japon (M_{w}=9.1 ); le séisme de Kushiro en 1993 au Japon (M_{w}=7.8) et le séisme d'Emilie Romagne en 2012 en Italie (M_{w}=5.9). Pour les trois premiers séismes, des enregistrements en profondeur et en surface sont disponibles. L'étude de ces trois premiers cas rend possible la validation du modèle par comparaison des données enregistrées et calculées. Le modèle peut donc être considéré comme un outil fiable pour la prédiction de la réponse sismique des sols saturés

Le terme gridshell désigne une coque discrète qui est obtenue par déformation élastique d'une grille bidirectionnelle continue plane sans rigidité en cisaillement puis rigidifiée par une troisième direction de barres. Ainsi défini, un gridshell a un potentiel structural intéressant et peut répondre à des exigences architecturales complexes. La recherche de forme de ces structures a été menée à travers l'histoire principalement par deux méthodes, la méthode du filet inversé et la relaxation dynamique. Ces deux méthodes permettent d'obtenir une forme approchée de celle proposée par l'architecte, dérivant d'une grille à plat et de conditions aux limites partiellement ou complètement imposées. Dans le cadre de cette thèse, nous nous sommes intéressés à générer un gridshell sur une surface à forme et contours imposés. Un outil numérique se basant sur la méthode du compas a été développé. Il permet de mailler un réseau de Tchebychev sur une surface connaissant son équation cartésienne. Un autre outil permettant le maillage se basant sur un calcul en éléments finis explicite a été mis en œuvre. La particularité de cette technique est de pouvoir tenir en compte des propriétés mécaniques de la structure et de simuler le comportement du gridshell. Des applications des deux méthodes sur des formes architecturalement intéressantes ont permis de voir les limitations de la possibilité de mailler une forme avec un réseau de Tchebychev. La méthode du compas a ensuite été couplée à des algorithmes métaheuristiques types génétiques. L'algorithme résultant permet d'optimiser un gridshell en minimisant la courbure dans les barres et donc les contraintes dans la structure introduites lors de la mise en forme. Il a été mis en œuvre et testé pour plusieurs surfaces / Gridshells are often defined as structures that have the shape and rigidity of a double curvature shell but they consist of a grid and not a continuous surface. They are obtained by elastic deformation of a two-way grid initially flat. The deformed grid is then rigidified using a third direction of bars. Thus, a gridshell has an interesting structural potential and can respond to complex architectural requirements. Two methods have been used through out history for the form finding of gridshells, the inversion method and the dynamic relaxation method. Both techniques lead to a deformed grid which is a result of calculations. The form obtained is closed to the one proposed by the architect. A numerical tool based on the compass method is developed in this thesis. It allows mapping aTchebychev net on an imposed form and imposed boundary conditions.Another tool based on an explicit dynamic finite element calculationis proposed. The particularity of this technique is to be able to take into account the mechanical properties of the structure and to simulate the gridshell behavior. Applications of both methods on differents forms show the limitations of mapping a Tchebychev net on an imposed form. The compass method has been coupled with geneticalgorithms. The algorithm optimizes the gridshell by minimizing the curvature in bars in order not to break the bars during the construction. It has been implemented and tested on several surfaces

L'industrie de la pêche est confrontée à plusieurs contraintes comme l'augmentation du prix du carburant et la diminution des stocks de poisson. Pour améliorer la rentabilité, nous avons fait l'optimisation des chaluts de pêche (pélagique et de fond) afin d'augmenter leurs efficacités énergétiques. Deux facteurs principaux ont été utilisés : la traînée du chalut et sa surface balayée. Notre travail porte sur l'optimisation du plan du chalut par modification du nombre de mailles des pièces du filet et par changement de la longueur des câbles. Nous avons développé une méthode de recherche successive appelée SOT et nous l'avons comparée à une méthode de recherche aléatoire dans le cas d'un chalut pélagique avec un gain net de 54%. La SOT a été utilisée pour optimiser le plan des pièces du filet d'un chalut de fond. Une amélioration du gain énergétique de 38% a été obtenue, ainsi qu'une réduction de 45 jours de pêche. Par souci de risque de baisse de capturabilité, nous avons développé une fonction objectif qui dépend de la distribution de poisson. Trois fonctions de même densité surfacique ont été testées. L'amélioration de l'efficacité énergétique calculée par la méthode SOT sur un chalut de fond est respectivement de 52%, 16% et 32% selon les trois distributions. Avec la méthode SRT (méthode de raffinement séquentiel), nous obtenons des réductions de 56%, 30% et 39% pour les trois distributions. Finalement nous avons adapté la méthode SOT pour optimiser les longueurs des câbles du chalut. Nous obtenons dans un cas une amélioration de l'efficacité énergétique de 49%. L'utilité de ce travail devrait être étendue hors de la pêche comme l'écologie en général. / Nowadays, fishing industry is facing several constraints such as increased fuel prices and declining fish stocks. In this work, we perform optimization of fishing trawls (pelagic and bottom) in order to improve their energy efficiencies. Two main factors were used: the trawl drag and its mouth swept area. Our work focuses on optimizing trawl design by altering mesh numbers and changing cable length. We have developed a method called SOT and compared it to a random search method in the pelagic trawl case with a net energy gain of 54%. The SOT was used in design optimization of bottom trawl panels. Improved energy gain of 38% was obtained and a reduction of 45 total number of fishing days. In order to bypass catchability reduction, we developed an objective function which depends on fish distribution. Three functions with same bottom surface density were tested. Improving energy efficiency calculated by the SOT method on a bottom trawl is respectively 52%, 16% and 32% for the three distributions. With the SRT method (sequential refinement), we obtain reductions of 56%, 30% and 39% for the three aforementioned distributions. Finally we adapted the SOT method to optimize trawl cable lengths. In one case, we obtain energy efficiency improved by 49%. The usefulness of this work should be extended beyond fishery to ecology in general.

Ce travail de thèse vise à une meilleure compréhension et prise en compte des mécanismes d'endommagement ductile au cours des procédés de mise en forme à froid. Le tréfilage, le laminage de fils et le laminage à pas de pèlerin sont pris comme exemples d'application. Une grande attention est portée également à la méthodologie d'identification des paramètres des modèles d'endommagement implémentés. Trois approches de l'endommagement ductile ont été étudiées : des critères de rupture phénoménologique, des modèles phénoménologiques couplés et des modèles micromécaniques. Ces modèles ont été implémentés dans Forge®, et les algorithmes ont dû être adaptés à sa formulation mixte vitesse-pression et à son élément fini (P1+/P1). Parallèlement aux développements numériques, différents essais mécaniques ont été effectués sur trois matériaux différents (acier haut carbone, acier inoxydable et alliage de zirconium) pour identifier les paramètres des modèles d'écrouissage et d'endommagement. Des essais de traction in situ sous micro-tomographie aux rayons X ont été exploités pour l'identification des mécanismes d'endommagement ductile (germination, croissance et coalescence), ainsi que l'identification du modèle micromécanique. Enfin, la dernière partie de ce travail a consisté à réaliser des études comparatives de ces différents modèles d'endommagement pour les différents procédés de mise en forme et différents matériaux. Concernant le tréfilage et le laminage de l'acier inoxydable, un bon accord entre les simulations numériques et les résultats expérimentaux a été obtenu. Pour l'acier à haute teneur en carbone (perlitique), le modèle micro-mécanique GTN a donné le meilleur résultat à la fois qualitativement et quantitativement pour le tréfilage ultime. Par ailleurs, la comparaison entre les différents modèles sur différents procédés (laminage du fil avec l'acier haut carbone, laminage à pas de pèlerin avec l'alliage de zirconium) a aidé d'une part à mettre en évidence le rôle important du troisième invariant de la contrainte déviatorique sur la localisation de l'endommagement dans les procédés de mise en forme où le cisaillement est dominant. D'autre part, elle montre que le processus d'identification lui-même devrait se fonder sur la mesure de la microstructure afin de garantir un résultat précis lors de l'application aux procédés. / The present PhD thesis aims at a better understanding and modeling of ductile damage mechanisms during cold forming processes, with wire drawing, rolling and cold pilgering as examples. In addition, special attention is paid to implemented damage models parameters identification methodology. All three approaches of ductile damage were investigated: uncoupled phenomenological fracture criteria; coupled phenomenological models; micromechanical model. These models have been implemented in Forge®, which required adaptation of algorithms to its mixed velocity-pressure formulation and to its finite element (P1+/P1). Parallel to the numerical work, various mechanical tests on three different materials (high carbon steel, stainless steel and zirconium alloy) were carried out for work hardening, and damage models parameters identification. In situ X-ray micro-tomography tensile tests have also been exploited for the identification of ductile damage mechanisms (nucleation, growth and coalescence) as well as the identification of micromechanical model. Finally, we carried out comparative studies of these models on our three abovementioned forming processes and materials. Regarding wire drawing and rolling of stainless steel, good agreement between numerical simulations and experimental results was found. For high carbon perlitic steel ultimate wire drawing, the GTN micro-mechanical model has given the best result, both qualitatively and quantitatively. Moreover, the comparison of the different models on different processes (wire rolling on high carbon steel, cold pilgering on zirconium alloy) highlights on the one hand the important role of the third deviatoric stress invariant in damage localization for shear-dominated forming processes. It shows on the other hand that the identification process itself should be based on microstructure measurements to provide accurate results in forming application.

L’objectif principal de la thèse est de répondre à des questions liées aux étapes clé d’un processus de l’adaptation de modèles de plaques. Comme l’adaptativité dépend des estimateurs d’erreurs fiables, une part importante du rapport est dédiée au développement des méthodes numériques pour les estimateurs d’erreurs aussi bien dues à la discrétisation qu’au choix du modèle. Une comparaison des estimateurs d’erreurs de discrétisation d’un point de vue pratique est présentée. Une attention particulière est prêtée a la méthode de résiduels équilibrés (en anglais, "equilibrated residual method"), laquelle est potentiellement applicable aux estimations des deux types d’erreurs, de discrétisation et de modèle.Il faut souligner que, contrairement aux estimateurs d’erreurs de discrétisation, les estimateurs d’erreur de modèle sont plus difficiles à élaborer. Le concept de l’adaptativité de modèles pour les plaques est implémenté sur la base de la méthode de résiduels équilibrés et de la famille hiérarchique des éléments finis de plaques. Les éléments finis dérivés dans le cadre de la thèse, comprennent aussi bien les éléments de plaques minces et que les éléments de plaques épaisses. Ces derniers sont formulés en s’appuyant sur une théorie nouvelle de plaque, intégrant aussi les effets d’étirement le long de l’épaisseur. Les erreurs de modèle sont estimées via des calculs élément par élément. Les erreurs de discrétisation et de modèle sont estimées d’une manière indépendante, ce qui rend l’approche très robuste et facile à utiliser. Les méthodes développées sont appliquées sur plusieurs exemples numériques. Les travaux réalisés dans le cadre de la thèse représentent plusieurs contributions qui visent l’objectif final de la modélisation adaptative, ou une procédure complètement automatique permettrait de faire un choix optimal du modèle de plaques pour chaque élément de la structure. / The primary goal of the thesis is to provide some answers to the questions related to the key steps in the process of adaptive modeling of plates. Since the adaptivity depends on reliable error estimates, a large part of the thesis is related to the derivation of computational procedures for discretization error estimates as well as model error estimates. A practical comparison of some of the established discretization error estimates is made. Special attention is paid to what is called equilibrated residuum method, which has a potential to be used both for discretization error and model error estimates. It should be emphasized that the model error estimates are quite hard to obtain, in contrast to the discretization error estimates. The concept of model adaptivity for plates is in this work implemented on the basis of equilibrated residuum method and hierarchic family of plate finite element models.The finite elements used in the thesis range from thin plate finite elements to thick plate finite elements. The latter are based on a newly derived higher order plate theory, which includes through the thickness stretching. The model error is estimated by local element-wise computations. As all the finite elements, representing the chosen plate mathematical models, are re-derived in order to share the same interpolation bases, the difference between the local computations can be attributed mainly to the model error. This choice of finite elements enables effective computation of the model error estimate and improves the robustness of the adaptive modeling. Thus the discretization error can be computed by an independent procedure.Many numerical examples are provided as an illustration of performance of the derived plate elements, the derived discretization error procedures and the derived modeling error procedure. Since the basic goal of modeling in engineering is to produce an effective model, which will produce the most accurate results with the minimum input data, the need for the adaptive modeling will always be present. In this view, the present work is a contribution to the final goal of the finite element modeling of plate structures: a fully automatic adaptive procedure for the construction of an optimal computational model (an optimal finite element mesh and an optimal choice of a plate model for each element of the mesh) for a given plate structure.

Le travail présenté dans cette thèse propose des nouvelles méthodes de diagnostic des défauts dans les machines à courant alternatif. Les procédures de diagnostic décrites sont fiables, originales, peu coûteuses et simples à mettre en œuvre. Elles présentent l’avantage d’être non invasives et viennent éliminer le principal inconvénient présenté par les autres méthodes de diagnostic déjà connues qui ont recours à la comparaison avec un état supposé sain et connu. L’analyse porte principalement sur le champ magnétique de dispersion, plus particulièrement ses composantes axiales et radiales qui présentent des sensibilités différentes aux défauts. À cet effet, les phénomènes dus aux défauts de court-circuit entre spires dans le bobinage statorique, l’excentricité et la rupture de barres rotoriques, sont étudiés. Il est abordé en particulier les effets de ces défauts dans le cas de la machine asynchrone et celui de la machine synchrone sur le champ axial et sur le champ radial. Pour chaque défaut, des signatures spécifiques sont mises en évidence et expliquées par la modélisation analytique et la modélisation numérique qui comprend les méthodes par circuits électriques couplés et méthodes par éléments finis. Des défauts ont été créés au rotor et au stator sur différentes machines, afin de valider expérimentalement les méthodes de diagnostic proposées. / The work described in this thesis proposes new procedures to diagnose faults in AC machines. The diagnostic procedures described are reliable, original, inexpensive and simple to implement. They have the advantage of being noninvasive and just get rid from the main drawback presented by other diagnostic methods based on a comparison with a healthy state assumed to be known. The analysis focuses on the magnetic field dispersion in the vicinity of the machine, especially its radial an axial distribution which presents different sensitivity according to various faults. To this end, the phenomena due to inter-turn short-circuit faults in the stator winding, rotor eccentricity and broken rotor bars, are studied in the case of an induction machine and a synchronous machine. For each fault, specific signatures are identified and justified by analytical modeling and numerical method of analysis including coupled electric circuit and finite element methods. Defects have been created in the rotor and stator on different machines in order to validate experimentally the suggested diagnostic procedures.

L’opération de cellules solaires sous un flux de lumière concentrée permet d’accroître leur efficacité. Les densités élevées de courant généré limitent toutefois ce gain et font chuter l’efficacité par effet Joule. L’aire métallisée des contacts en face avant des cellules solaires doit alors être augmentée avec le facteur de concentration. Cela provoque des pertes par ombrage qui deviennent très importantes dans le cas des contacts conventionnels en grille. Dans ce mémoire, une nouvelle architecture de contact se basant sur la technologie des vias traversants est étudiée et comparée aux contacts conventionnels. Celle-ci pourrait permettre une réduction simultanée des pertes par effet Joule et par ombrage pour un faible coût. Un modèle de simulations combinant l’évaluation de résistance série par éléments finis 3D et le calcul des performances avec un circuit équivalent 1D a été développé. Cela a permis d’optimiser des cellules solaires InGaP/(In)GaAs/Ge avec les différentes géométries de contacts sous lumière hautement concentrée. Un gain d’efficacité allant jusqu’à 3%abs et une augmentation de la densité de puissance générée de l’ordre de 20% ont ainsi été obtenus avec les cellules à vias traversants. Un modèle de circuit distribué 3D a également été développé pour simuler l’illumination non-uniforme des cellules solaires. Il a alors été montré que la dégradation des performances associée à ce phénomène est égale ou inférieure dans les cellules à vias traversants par rapport aux cellules conventionnelles. Enfin, des simulations par éléments finis ont été réalisées pour estimer l’impact de la qualité des parois des vias sur les recombinaisons non-radiatives. Il a alors pu être déterminé que ce mécanisme de pertes ne devrait pas constituer un frein au développement de la technologie des vias traversants.

Cette thèse porte sur la modélisation des ondes électromagnétiques en milieu complexe et à haute fréquence. Ceci reste un défi, de part les besoins en temps de calcul et en place mémoire que cela suscite. L'idée de la méthode Finite Element Tearing and Interconnect (FETI) est de diviser le domaine en sous-domaines distincts dans lesquels des problèmes électromagnétiques peuvent être résolus de manière indépendante. Une modification de la méthode classique, appelée FETI-DPEM2-full, est ici proposée afin d'améliorer le processus numérique. Les champs diffractés par différents objets ont ainsi été simulés et comparés avec succès à des champs expérimentaux. A partir de ces champs, les techniques d'imagerie micro-onde permettent de déterminer les paramètres électromagnétiques des objets illuminés. Ce problème inverse peut être réécrit sous la forme d'un problème d'optimisation faisant appel à la résolution de deux problèmes directs à chaque itération. Une méthode de type quasi-Newton a donc été couplée efficacement avec la méthode FETI-DPEM2-full. Des reconstructions obtenues à partir de champs mesurés sont présentées et valident la chaîne méthodologique proposée. / This PhD work is devoted to the electromagnetic modeling of threedimensional large-scale wave propagation problems, which is very challenging in terms of memory and computation time. The main idea of the Finite Element Tearing and Interconnect (FETI) method is to divide the domain into non-overlapping subdomains where each electromagnetic problem can be independently solved. A modification of the classical FETI method, called the FETI-DPEM2-full method, is proposed herein to improve the numerical process. The fields scattered by various structures have thus been computed and successfully compared with measured fields. From these experiments, quantitative microwave imaging algorithms attempt to estimate the physical parameters of the illuminated target. This inverse problem is recasted into a minimization problem where two forward problems are required at each iteration step. A quasi-Newton algorithm has thus been efficiently coupled with the FETI-DPEM2-full method. Reconstructions of various targets from measured scattered fields have been successfully performed, which validates the effectiveness of the proposed methodology.

Dans l'étude des éléments combustibles des réacteurs à eau pressurisée, cette thèse s'inscrit dans la compréhension et la modélisation du comportement viscoplastique du dioxyde d'uranium (UO2) à l'échelle du polycristal. Lors de fonctionnement de type incidentel du réacteur, le combustible subit une forte élévation de la température avec un gradient thermique de la pastille engendrant des déformations viscoplastiques contrôlées par des mouvements de dislocations. D'abord, un modèle de plasticité cristalline a été développé de manière à décrire l’anisotropie viscoplastique du matériau en fonction de la température et de la vitesse de sollicitation. Des simulations par éléments finis (EF) sur monocristaux ont permis d’identifier que les trois modes de glissement généralement observés dans l'UO2 sont importants pour décrire le comportement anisotrope du matériau. Dans un second temps, les coefficients de la matrice d'interactions entre dislocations ont été déterminés spécifiquement pour l’UO2 afin d’améliorer la modélisation des polycristaux. En effet, en calculant par EF les dislocations géométriquement nécessaires, qui sont responsables d’une forte augmentation de la densité de dislocations stockées dans les polycristaux, les interactions entre dislocations permettent de simuler l’effet dé taille de grain et l’écrouissage des pastilles. Finalement, le modèle, adapté pour les polycristaux, a été validé par comparaison avec les essais expérimentaux sur pastille et par comparaison du comportement intra-granulaire simulé avec des mesures EBSD. Grâce à cette dernière comparaison, il est possible de remonter indirectement aux hétérogénéités de déformation dans les grains / This thesis is part of the study of fuel elements of pressurized water reactors and, more specifically, focus on the understanding and modelling of the viscoplastic behavior of uranium dioxide (UO$_2$) at polycrystalline scale. During the incidental operation of the reactor, the fuel undergoes a strong increase of temperature and thermal gradient between the center and the periphery of the pellet leading to viscoplastic strains due to dislocation movement mechanisms. First, a crystal plasticity model was developed in order to describe the viscoplastic anisotropy of the material considering the temperature and the loading rate. Finite element (FE) simulations on single crystals enabled to highlight that the three slip modes generally observed in UO$_2$ are crucial to describe the anisotropic behavior of the material. Secondly, coefficients of the interaction matrix have been identified specifically for UO$_2$ in order to improve the polycrystal modelling. Indeed, by calculating geometrically necessary dislocations (GNDs), which are responsible of the great increase of the stored dislocation density in polycrystals, the interactions between dislocations enable to simulate de grain size sensitivity and hardening of the fuel pellet. Finally, the model adapted for polycrystals, have been validated by comparing FE simulations with pellet compression tests and by comparing the simulated intra-granular behavior with EBSD measurements. Thanks to the latter comparison, it is possible to indirectly compare the strain heterogeneities in the grains