La prédiction du clivage des aciers ferritiques a été largement étudiée à l'aide de l'approchelocale de la rupture, et des modèles macroscopiques identifiés phénoménologiquement comme celui de Beremin. Cette prédiction reste cependant difficile dans le domaine de transition ductilefragile. Cela a conduit à des études micromécaniques par les approches polycristallines afin de décrire l'évolution de la contrainte de clivage en fonction de la température pour les aciers bainitiques. Dans cette étude, on utilise une approche polycristalline, puis on développe un modèle macroscopique de prédiction du clivage d'un acier de microstructure plus simple, un acier sphéroïdisé. De nombreux résultats de la littérature indiquent que leur rupture est due à la microfissuration des carbures. Le comportement mécanique et la rupture d'un acier sphéroïdisé C35R sont obtenus par des essais de traction simple et de ténacité dans une gamme de température entre -196 et 20°C. Les analyses microstructurales sont effectuées pour déterminer la distribution de taille des grains et des carbures. Des modélisations simplifiées de la microstructure de l'acier sont proposées. L'aspect polycristallin du matériau est pris en compte. Une étude paramétrique concernant la distribution des contraintes dans les carbures est réalisée. On montre que les paramètres du modèle polycristallin n'influencent pas cette distribution à la condition de représenter le même comportement global en traction. La prédiction de la rupture par clivage est basée sur une approche probabiliste, considérant la dispersion des contraintes dans les carbures due à l'hétérogénéité des champs mécaniques issue de la modélisation polycristalline. La probabilité élémentaire de rupture des carbures est ainsi obtenue. Différents modèles de rupture macroscopiques sont alors développés, basés sur des critères en germination et propagation des microfissures. Ils sont appliqués à une éprouvette SENT afin de prédire la ténacité dans le bas de la transition ductile-fragile. La comparaison avec les résultats expérimentaux montre que l'on doit prendre en compte l'évolution de la densité volumique des microfissures avec le chargement, l'extension des microfissures de taille variable ou leur émoussement. L'importance des différents critères dépend de la température.

Si l'amorçage des fissures de fatigue est aujourd'hui bien compris dans le cas de monocristaux de métaux purs, ce phénomène s'avère plus complexe à cerner et à prédire dans le cas d'alliages métalliques polycristallins tels que l'acier AISI 316LN.D'un point de vue expérimental, notre étude s'est concentrée sur la caractérisation du comportement mécanique et l'étude, à différentes échelles, des phénomènes liés à l'amorçage des fissures de fatigue oligocyclique dans l'acier 316LN. Pour des niveaux de déformation appliquée de Δε/2 = 0,3 et 0,5%, l'adoucissement cyclique observé au cours des essais coïncide avec l'organisation des dislocations sous forme de bandes. Ces bandes sont liées à l'activation des systèmes de plus haut facteur de Schmid. Elles portent la majeure partie de la déformation et provoquent en surface des intrusions et extrusions favorisant l'apparition et la coalescence de fissures.D'un point de vue modélisation, nous avons proposé un nouveau modèle de plasticité cristalline intégrant des dislocations géométriquement nécessaires (GND) directement calculées à partir du second gradient de la rotation élastique. Implémenté dans les codes d'éléments finis AbaqusTM et Cast3mTM, ce modèle s'inspire des travaux sur le monocristal en transformations finies de Peirce et al. (1983) et de Teodosiu et al. (1993). Adapté au cas des polycristaux par Hoc (1999) et Erieau (2003), il a été enrichi par l’introduction GND selon la théorie proposée par Acharya et Bassani (2000). Les simulations réalisées sur des différents types d'agrégats (2D extrudé et 3D) montrent que la prise en compte de GND permet :- de reproduire les effets de taille de grains au niveau macroscopique et local,- de décrire plus finement les champs de contraintes calculés.Ces simulations ont permis de mettre en évidence l'influence des matrices d'élasticité et d'écrouissage sur les valeurs et l'évolution des contraintes macroscopique effective et cinématique moyenne et le rôle important des conditions aux limites lors des calculs d'agrégats. / If fatigue crack initiation is currently quite well understood for pure single crystals, its comprehension and prediction in cases of polycrystal alloys such as AISI 316LN stainless steel remain complicated.Experimentally our study focuses on the characterisation of the mechanical behaviour and on the study at different scales of the phenomenon leading to low cycle fatigue crack initiation in 316LN stainless steel. For straining amplitudes of Δε/2 = 0,3 and 0,5%, the cyclic softening observed during testing has been related to the organisation of dislocations in band structures. These bands, formed due to the activation of slip systems having the greatest Schmid's factor, carry the most part of the deformation. Their emergence at free surfaces leads to the formation of intrusions and extrusions which help cracks initiate and spread.Numerically we worked on the mesoscopic scale, proposing a new model of crystalline plasticity. This model integrates geometrically necessary dislocations (GND) directly computed from the lattice curvature. Implemented in the finite element code AbaqusTM and Cast3mTM, it is based on single crystal finite deformations laws proposed by Peirce et al. (1983) and Teodosiu et al. (1993). Extended for polycrystals by Hoc (2001) and Erieau (2003), it has been improved by the introduction of GND (Acharya and Bassani, 2000). The simulations performed on different types of aggregates (2D/3D) have shown that taking GND into account enables:- the prediction of the grain size effect on a macroscopic and on a local scale,- a finer computation of local stress field.The influence of the elasticity and interaction matrices on the values and the evolution of the isotropic and kinematic mean stresses has been shown. The importance of boundary conditions on computed mechanical fields could also be pointed out.

Modéliser la déformation plastique des polycristaux est un objectif majeur de la science des matériaux. Tous les modèles actuels comportent une partie phénoménologique n´ecessitant un ajustement de paramètres sur des résultats expérimentaux. Cette thèse vise à mettre en place un nouveau modèle, justifié physiquement, sans paramètre ajustable et adapté à la modélisation du polycristal CFC en sollicitation monotone. Afin de mesurer les champs mécaniques à l’échelle du micromètre, des mesures de microdiffraction Laue ont été réalisées sur un tricristal de cuivre à une faible déformation plastique. Ces mesures nous renseignent sur les mécanismes plastiques intervenant très près des joints de grains et définissent des états de référence pour les simulations. On montre principalement que la déformation plastique s’accompagne d’un stockage de dislocations géométriquement nécessaires (GND) aux joints de grains, en relation avec l’apparition de contraintes internes à longue distance. Des simulations de Dynamique des Dislocations dans des bicristaux ont ´et´e réalisées afin de caractériser les phénomènes physiques mis en œuvre. Ces simulations confirment que l’interaction dislocations - joints de grains s’accompagne d’un stockage de GND sous la forme de microstructures tridimensionnelles très inhomogènes. Les propriétés mécaniques induites par ce phénomène complexe peuvent être quantifiées par des lois continues élaborées à partir de l’approximation théorique d’un empilement unidimensionnel. Les lois de comportement ainsi définies ont ensuite été incorporées dans une modélisation micromécanique de plasticité cristalline, jusqu’ici dédiée au monocristal CFC. Le modèle ainsi construit a maintenant la capacité de prédire les mesures réalisées sur le tricristal de cuivre. Nous avons ainsi mis en place un modèle physiquement justifié et adapté `a la modélisation du polycristal CFC en sollicitation monotone. / The modeling of strain hardening in polycrystals is a difficult and still standing task. Current existing models are partly phenomenological, as they always consider constitutive parameters adjusted on the experiment. The aim of the present study is to design a physically based model for the basic problem of monotonic deformation in the FCC polycrystal. Laue microdiffraction is used to measure the mechanical fields in the vicinity of grain boundaries in a copper tricrystal compress at 0.2%. These measurements aim to characterize the plastic phenomena involved and to provide experimental data as bench results for the simulations. Evidences of geometrically necessary dislocations (GND) storage close to the grain boundaries are given in relation with the apparition of longrange internal stresses. Dislocations Dynamics simulations are used to study the plastic strain close to a grain boundary in Cu bicrystals. We show that close to the boundaries plastic strain is associated to the storage of heterogeneous GNDs in complex 3D microstructures. The mechanical properties associate to such microstructure can be described with continuous laws based on a theoretical approximation assuming a 1D pile-up. The corresponding constitutive laws are then introduced in a crystal plasticity model initially devoted to FCC single crystal plasticity and solved with Finite Elements simulations. The new model we propose as now the capacity to reproduce or predict the experimental results we first obtained in the Cu tricrystal. In conclusion, a physically justified model is proposed to predict plastic deformation for the FCC polycrystal in monotonic deformation.

La liaison aube-disque d'un turboréacteur est soumise en opération à un chargement complexe composé d'une forte pression et d'un mouvement relatif oscillatoire, qui génère un phénomène de fretting–fatigue. L'apparition de celui-ci réduit la résistance en fatigue des structures et engendre des dégradations qui, sous différentes formes (usure, rupture), peuvent conduire à une perte de fonctionnalité. Afin d'assurer la fiabilité des turboréacteurs, il est donc utile d'améliorer l'estimation de la durée de vie sous chargement de fretting.Cette thèse est dédiée à la fois à l'étude de l'amorçage de la première fissure et de sa micropropagation dans le domaine dit "des fissures courtes", c'est-à-dire celles qui sont encore directement influencées par la microstructure locale. Pour remplir la première exigence, on développe un modèle numérique mettant en oeuvre un critère multiaxial de fatigue. Il permet de prédire la position de l'amorçage et le nombre de cycles correspondant, en prenant en compte d'éventuels traitements de surface. Dans un second temps, une étude numérique de la direction et de la cinétique de propagation des fissures est mise en place avec prise en compte de l'influence de la microstructure du matériau, à l'aide d'un modèle de plasticité cristalline qui prend en compte les différentes familles de systèmes de glissement. La vitesse de fissuration diminue à l'approche des joints de grains, ce qui est qualitativement conformes aux observations expérimentales. / The blade-disc fixings of an aircraft engine are subjected to complex loadings composed by a high pressure and oscillatory relative motions during operation. This corresponds to a loading type called fretting. The appearance of fretting phenomenon reduces the fatigue resistance of structuresand generates damage which, in various forms (wear, rupture), can lead to a loss of functionality. In order to ensure the reliability of turbo-engines, the improvement of the models for estimating life prediction is then needed.This thesis is devoted to the study of crack initiation and crack micro-propagation, in the so called "short crack domain", where they are still in interaction with the microstructure. At first, a multiaxial fatigue model is proposed. It predicts the location and the number of cycles to initiation, taken into account the possible surface treatments. Then, a numerical study of the direction and the kinetics of crack propagation is performed with consideration to the influence of the material microstructure. A crystal plasticity model is used, that takes into account all the relevant families of crystallographic slip systems. The crack plane is assumed to be one of the crystallographic plane where the plastic slip is the most intense. The crack propagation rate decreases when approaching grain boundaries, that is qualitatively consistent with the experimental observations.

Il est communément admis aujourd'hui que le contact entre deux surfaces se compose en réalité d'une multitude de contact ponctuels entre des aspérités. Cette considération amène à une surface de contact réelle significativement différente de l'aire de contact parfaite supposée dans la théorie de Hertz.De même, elle implique également la présence d'un espace libre entre les deux surfaces en contact. Dans cette situation, l'objectif principal de ces travaux de thèse est de développer des approches numériques permettant l'analyse du contact mécanique entre surfaces rugueuses dans le but de qualifier/quantifier l'étanchéité de ce contact rugueux.Deux approches différentes sont étudiées. La première consiste à analyser le contact mécanique entre une surface rugueuse et un plan rigide au moyen de la méthode des éléments finis et d'un nouveau modèle numérique. La seconde concerne l'estimation de la transmissivité d'un contact rugueux en considérant des simulations de l'écoulement d'un fluide au sein du champ des ouvertures présent entre les deux surfaces en contact.La comparaison de ces estimations numériques avec les résultats expérimentaux révèlent des écarts importants. Dans le but de comprendre ces écarts, l'influence du modèle de comportement matériau dans de telles simulations est étudiée. La plasticité cristalline, mais également l'élévation de la température dégagée par déformation plastique seront considérés. La question de la représentativité de notre problème vis-à-vis de l'approche fluide sera également discutée. / It is now widely accepted that the contact between two surfaces is in fact a one-to-one contact between many asperities, depending on the roughness of the contact pair. This represents a strong deviation from the perfect contact assumed in engineering approaches, with the real contact area being significantly smaller than the apparent contact area.In addition, such an approach also implies the presence of a free space between surfaces in contact. Thus, the purpose of our study is the development of numerical tools to analyse the mechanical contact between rough surfaces and the tightness of such a contact.Two different approches are studied. The first one is devoted to the observation of the mechanical contact between a rough surface and a rigid plane by means of the finite element method and a new numerical model. The second approach is related to the estimation of the contact transmissivity by considering simulations of the fluid flow in the resulting free space between surfaces in contact.The comparison between experimental and numerical estimations of the resulting flow rate shows significant deviations.In order to understand these deviations, the influence of material models on contact simulations is studied. In particular, the effect of using crystal plasticity and the effect of including plasticity-driven temperature changes are considered. The question of problem representation in fluid flow approaches is also addressed.

Une méthode de corrélation d'images numériques adaptée aux mesures cinématiques dans les polycristaux a été développée dans ce travail, dans le but d'identifier les paramètres de lois de plasticité cristalline. Des mesures de champs 2D sont réalisées à la surface d'un polycristal d'acier austénitique 316LN, à partir d'une série d'images acquises au Microscope Electronique à Balayage (MEB) au cours d'un essai de traction monotone in-situ, et pour différentes tailles moyennes de grains. Pour permettre la corrélation d'images, un mouchetis adapté à l'échelle microscopique est au préalable déposé à la surface de l'éprouvette par une technique de microlithographie. Les distorsions spatiales liées à l'imagerie MEB et à la technique de marquage sont quantifiées. La connaissance de la microstructure en surface par diffraction des électrons rétrodiffusés permet de réaliser les mesures cinématiques en utilisant un maillage éléments finis non-structurés s'appuyant sur les joints de grains ou de macles. Ce même maillage est ensuite utilisé pour la simulation de chaque essai de traction sur la microstructure expérimentale, avec comme conditions aux limites les déplacements nodaux mesurés au cours du temps sur les bords du domaine. Deux lois de plasticité cristalline d'approche locale sont alors considérées pour simuler les hétérogénéités des déformations observées : la loi de Méric-Cailletaud et la loi DD-CFC développée à EDF R&D et davantage basée sur la physique. Des comparaisons entre la mesure et la simulation sont menées, en termes de déplacements, de déformations, mais aussi de systèmes de glissement activés. Finalement, une méthode inverse d'identification des paramètres de loi est proposée, se basant sur le recalage à la fois des champs de déplacements locaux et du comportement homogène du matériau. Les paramètres de l'écrouissage isotrope de la loi de Méric-Cailletaud sont ainsi identifiés pour différentes tailles moyennes de grains. Il est également montré que certain coefficients de la matrice d'interaction entre systèmes de glissement peuvent être estimés. / A digital image correlation procedure adapted to kinematic measurements in polycrystals has been developed in this work to identify parameters of crystal plasticity laws. 2D kinematic measurements are performed on the surface of 316LN austenitic steel polycrystals from a sequence of images acquired using a Scanning Electron Microscope (SEM) during in-situ tensile tests for various mean grain sizes. To enable digital image correlation, a speckle adapted to the microscopic scale is deposited onto the specimen surface by a microlithography process. Spatial distortions resulting from both patterning and SEM imaging techniques are quantified. The knowledge of the microstructure at the surface by electron backscattered diffraction allows for kinematic measurements to be performed using an unstructured finite element mesh taking as support the grain or twin boundaries. This same mesh is then used for the simulation of each tensile test on the experimental microstructure with the measured nodal displacements prescribed as boundary conditions with their time evolution. Two local crystal plasticity laws are considered to simulate the observed strain heterogeneities, namely, the Méric-Cailletaud model and the DD-CFC law developed at EDF R&D. Comparisons between measurements and simulations are performed in terms of displacements, strains but also activated slip systems. Last, an inverse identification method is proposed for the identification of the sought constitutive parameters based on both the local displacement fields and the material homogenized behavior. The parameters associated with isotropic hardening of Méric-Cailletaud law are thus identified for various mean grain sizes. It is also shown that some of the interaction parameters of slip systems can be estimated.

Les observations de microscopie électronique in situ effectuées par Daniel Caillard (CEMES, Toulouse) au cours de la déformation de cristaux de symétrie cubique centrée ont montré que les dislocations vis effectuaient des sauts de plusieurs distances inter-atomiques alors que la théorie standard de Peierls prédit des sauts de une seul distance inter-atomique. Nous avons étudié par simulation atomique le glissement d'une dislocation vis dans un cristal de fer pure. Nous montrons que la propagation de décrochement le long de la dislocation induit un échauffement local qui favorise la nucléation de décrochements supplémentaires. L'accumulation de ces décrochements permet à la dislocation de parcourir plusieurs distances inter-atomiques. Ces simulations nous permettent de proposer une théorie pour l'explication des observations de D. Caillard. / In situ straining tests in high purity α-Fe thin-foils at low temperatures have demonstrated that crystalline defects, called dislocations, have a jerky type of motion made of intermittent long jumps of several nanometers. Such an observation is in conflict with the standard Peierls mechanism for plastic deformation in bcc crystals, where the screw dislocation jumps are limited by inter-reticular distances, i.e. of a few Angstroms. Employing atomic-scale simulations, we show that although the short jumps are initially more favorable, their realization requires the propagation of a kinked profile along the dislocation line which yields coherent atomic vibrations acting as traveling thermal spikes. Such local heat bursts favor the thermally assisted nucleation of new kinks in the wake of primary ones. The accumulation of new kinks leads to long dislocation jumps like those observed experimentally. Our study constitutes an important step toward predictive atomic-scale theory for materials deformation.

L’identification de certaines lois de comportement demeure délicate surtout pour des comportements non linéaires. Plusieurs stratégies d’identification ont été développées en s’appuyant principalement sur le dialogue essai-calcul. La partie expérimentale consiste à mesurer une ou plusieurs quantités physiques qui caractérisent le comportement du matériau testé et la partie théorique permet de calculer ces mêmes quantités. Pour remonter aux propriétés matérielles recherchées, un accord entre les quantités mesurées et calculées doit être visé. Dans cette thèse, on s’intéresse au zinc qui un matériau polycristallin. Celui-ci est modélisé à travers des lois de plasticité cristalline qui permettent de prédire son comportement mécanique sous chargements complexes ainsi que l’évolution de sa microstructure. Les lois de plasticité cristalline implémentées dans le logiciel commercial Abaqus par Huang et Marin que nous avons adapté au cas du zinc ont été utilisées. L’étape d’identification de ces lois non linéaires est basée sur une approche inverse qui tient compte des hétérogénéités des déformations à l’échelle des grains. Pour cela, nous avons utilisé l’essai de nanoindentation qui est un essai hétérogène permettant d’extraire plus d’informations que les essais homogènes. A cet effet, et pour caractériser la plasticité cristalline du zinc, nous avons utilisé les courbes charge - profondeur de pénétration et les profils de déplacements mesurés sur l’empreinte résiduelle sur des grains de différentes orientations cristallographiques mesurées par EBSD. La confrontation des résultats expérimentaux et montrent la bonne adéquation entre les données expérimentales et les modèles identifiés. Les résultats obtenus ont ainsi permis de caractériser les mécanismes de déformation des monocristaux de zinc. / The identification of certain behavior laws of remains delicate especially for nonlinear behaviors. Several identification strategies have been developed based mainly on the experiment-calculation dialogue. The experimental part consists in measuring one or more physical quantities which characterize the behavior of the tested material and the theoretical part makes it possible to calculate these same quantities. To go back to the material properties sought, an agreement between the quantities measured and calculated must be aimed at. In this thesis, we are interested in zinc, which is a polycrystalline material. It is modeled through crystal plasticity laws that predict its mechanical behavior under complex loadings and the evolution of its microstructure. The crystal plasticity laws implemented in the commercial software Abaqus by Huang and Marin that we adapted to the case of zinc were used. The step of identifying these non-linear laws is based on an inverse approach that takes into account the heterogeneities of the deformations at the grain scale. For this, we used the nanoindentation test which is a heterogeneous test to extract more information than the homogeneous tests. For this purpose, and to characterize the crystal plasticity of zinc, we used the load - penetration depth curves and the displacement profiles measured on the residual imprint on grains of different crystallographic orientations measured by EBSD. The confrontation of the experimental and numerical results show the good agreement between the experimental data and the identified models. The results obtained made it possible to characterize the deformation mechanisms of zinc single crystals.

Des simulations de dynamique de dislocation (DD) sont utilisées pour l’étude de l'effet Hall-Petch (HP) et des contraintes internes à long-portée induites par les hétérogénéités de déformation dans les matériaux polycristallins.L'effet HP est reproduit avec succès grâce à des simulations de DD réalisées sur de simples agrégats polycristallins périodiques composés de 1 ou de 4 grains. De plus, l'influence de la forme des grains a été explorée en simulant des grains avec différents rapports d'aspect. Une loi généralisée de HP est proposée pour quantifier l'influence de la morphologie du grain en définissant une taille de grain effective. La valeur moyenne de la constante HP $K$ calculée avec différentes orientations cristallines à faible déformation est proche des valeurs expérimentales.Les dislocations stockées pendant la déformation sont principalement localisées à proximité des joints de grain et peuvent être traitées comme une distribution surfacique de dislocations. Nous avons utilisé des simulations DD pour calculer les contraintes associées aux parois de dislocations de différentes hauteurs, longueurs densités et caractères. Dans tous les cas, la contrainte est proportionnelle à la densité surfacique de dislocations géométriquement nécessaires (GNDs) et sa variation est capturée par un ensemble d'équations empiriques simples. Une prévision de contraintes à long-portée dans les grains est réalisée en sommant les contributions des GNDs accumulées de part et d’autre des joints de grains.L'augmentation de la contrainte interne liée au stockage de GNDs est linéaire avec la déformation plastique et est indépendante de la taille des grains. L'effet de taille observé dans les simulations de DD est attribué au seuil de déformation plastique, contrôlé par deux mécanismes concurrents : la contrainte critique de multiplication des sources et la contrainte critique de franchissement de la forêt. En raison de la localisation de la déformation dans les matériaux à gros grains, le modèle d’empilement des dislocations doit être utilisé pour prédire la contrainte critique dans ce cas. En superposant cette propriété aux analyses que nous avons fait à partir de simulations de DD dans le cas d'une déformation homogène, l'effet HP est justifié pour une large gamme de tailles de grains. / Dislocation Dynamics (DD) simulations are used to investigate the Hall-Petch (HP) effect and back stresses induced by grain boundaries (GB) in polycrystalline materials.The HP effect is successfully reproduced with DD simulations in simple periodic polycrystalline aggregates composed of 1 or 4 grains. In addition, the influence of grain shape was explored by simulating grains with different aspect ratios. A generalized HP law is proposed to quantify the influence of the grain morphology by defining an effective grain size. The average value of the HP constant K calculated with different crystal orientations at low strain is close to the experimental values.The dislocations stored during deformation are mainly located at GB and can be dealt with as a surface distribution of Geometrically Necessary Dislocations (GNDs). We used DD simulations to compute the back stresses induced by finite dislocation walls of different height, width, density and character. In all cases, back stresses are found proportional to the surface density and their spatial variations can be captured using a set of simple empirical equations. The back stress calculation inside grains is achieved by adding the contributions of GNDs accumulated at each GB facet.These back stresses are found to increase linearly with plastic strain and are independent of the grain size. The observed size effect in DD simulations is attributed to the threshold of plastic deformation, controlled by two competing mechanisms: the activation of dislocation sources and forest strengthening. Due to strain localization in coarse-grained materials, the pile-up model is used to predict the critical stress. By superposing such property to the analysis we made from DD simulations in the case of homogeneous deformation, the HP effect is justified for a wide range of grain sizes.

Les composants monocristallins fonctionnant à des températures élevées sont soumis à des conditions de chargement thermo-mécanique sévères. La géométrie et le comportement de ces composants sont très complexes. Un défi majeur est de développer des modèles mathématiques afin de prévoir l'initiation et la propagation de fissures en présence de contraintes importantes et de forts gradients de température. Dans ce cas, le comportement élastoviscoplastique fortement anisotrope du matériau étudié (superalliage à base Ni) doit être pris en compte. Le modèle correspondant doit être en mesure de rendre compte de la croissance anisotrope des fissures et de leur bifurcation dans des champs de contrainte complexes. De plus, le modèle doit être capable de prédire non seulement le taux de croissance des fissures mais aussi les chemins de fissuration. La mécanique de l'endommagement anisotrope est un cadre théorique bien adapté au développement de modèles de croissance de fissures dans les monocristaux. Au cours d'études précédentes, une loi de comportement couplant plasticité cristalline et endommagement cyclique a été développée, démontrant l'intérêt de cette approche, mais aussi ses limites, notamment du fait de la dépendance au maillage des résultats. Le développement récent de modèles non-locaux dans le cadre de la mécanique des milieux continus pourrait ainsi aider à surmonter ces difficultés. Une grande base expérimentale existe concernant l'initiation et la propagation de fissures dans les superalliages monocristallins à base de nickel. Les simulations thermomécaniques par éléments finis des aubes de turbine fournissent des informations détaillées sur la distribution des contraintes et des déformations plastiques, en particulier près de singularités géométriques comme les trous et les fentes de refroidissement. Tout d'abord, sur la base de la théorie de la plasticité cristalline qui établit un lien solide entre les contraintes et les déformations plastiques, un modèle découplé en mécanique de l'endommagement basé sur l'historique des calculs par éléments finis sera présenté. Ensuite, un modèle d'endommagement incrémental basé sur les milieux généralisés sera proposé et enfin, les prédictions du modèle pour l'initiation et la croissance de micro-fissures en résolvant le problème de dépendance au maillage seront discutés.