Cette thèse a pour objectif d'étudier la dynamique spatio-temporelle des atomes refroidis par laser dans un piège magnéto-optique (PMO). Il a été montré qu'un nuage d'atomes froids dans le régime de diffusion multiple peut présenter un comportement instable sans modulation externe du système. Cependant, ces instabilités n'ont pas encore été modélisées de façon satisfaisante. Une nouvelle configuration du PMO a été mise en oeuvre pour tenter d'étudier des instabilités dans une seule direction. Ce PMO, qualifié d'anisotrope, n'utilise pas des lasers de mêmes fréquences dans chaque direction de l'espace. Il met en évidence les forts couplages existants entre les directions du piège, si bien qu'il n'est pas possible de l'utiliser pour réduire le nombre de dimensions dans lesquelles les instabilités s'établissent. Toutefois, cette étude constitue un premier pas vers une meilleure description tridimensionnelle du piège. Elle nous a notamment permis de mesurer la probabilité pour qu'un photon diffusé soit réabsorbé à l'intérieur du nuage. Cette quantité est caractéristique du PMO mais elle n'avait jusqu'à là jamais été mesurée. Nous avons également établi un modèle spatio-temporel unidimensionnel du PMO. Il est constitué d'un système d'équations non-linéaires couplées reliant la densité atomique et les intensités des faisceaux lasers. Ce système contient notamment une équation de Vlasov-Fokker-Planck, rencontrée dans de nombreux domaines de la physique. Des simulations numériques ont été effectuées dans un cas simple. D'un point de vue expérimental, l'utilisation d'une caméra rapide nous a permis de mettre en évidence la structure spatiale d'instabilités de type stochastique. / The aim of this thesis is to study the spatio-temporal dynamics of laser cooled atoms in a magneto-optical trap (MOT). Recent works have shown that in the multiple scattering regime, an atomic cloud can have an unstable behavior without external modulation of the system. Nevertheless, these instabilities have not yet been modeled in a satisfactory way. A new configuration of the MOT has been built up as a possible way to study instabilities in only one direction. This trap, called anisotropic MOT, is not made of laser beams with the same laser frequencies along each direction of space. It exhibits the strong couplings between the directions of the trap, with the result that it cannot be used to reduce the number of directions in which instabilities grow up. However, this study can be considered as a new step to a better 3D description of the MOT physics. In particular, it gives us a way to measure the probability that a scattered photon is reabsorbed inside the atomic cloud. This quantity is a characteristic of the MOT but it has never been measured so far. We also develop a 1D spatio-temporal model of the MOT. It consists in a set of coupled nonlinear equations linking the atomic density and the laser intensities. This set contains a Vlasov-Fokker-Planck equation which is used to model a lot of systems in various fields and not only in physics. Numerical simulations have been done in a simple case. In the experiment, the use of a fast video camera allows us to observe the spatial structure of one type of instabilities, the so-called stochastic instabilities.

La magnétorésistance anisotrope (AMR) des matériaux ferromagnétiques est largement utilisée comme le phénomène de base pour la mesure ou la détection de champ magnétique. En raison de la relation entre la configuration en domaines magnétiques et la résistivité macroscopique, l'application d'un champ magnétique externe modifie la résistivité des matériaux ferromagnétiques. Bien que cet effet soit largement utilisé dans des applications industrielles, certains aspects fondamentaux du comportement AMR sont encore assez mal compris. Par exemple, le rôle de la texture cristallographique dans le comportement effectif n'est pas décrit avec précision par les outils classiques de modélisation. En raison de ce lien direct entre la microstructure en domaines et l'effet AMR, les modèles de description de l'effet AMR reposent généralement sur des calculs micromagnétiques. Pour ces calculs, le nombre de degrés de liberté et d'interactions peuvent se multiplier rapidement si on recherche à décrire un comportement macroscopique (cas des polycristaux par exemple).La thèse porte sur la modélisation numérique de l'effet de magnétorésistance anisotrope des matériaux ferromagnétiques. Ce nouvel outil de modélisation 3D peut remédier à cet inconvénient majeur des approches micromagnétiques. Un modèle permettant de décrire les effets de couplage magnéto-élastique en utilisant une approche micro-macro est disponible au laboratoire GeePs. Sur la base des mêmes principes de la modélisation micro-macro, un outil de simulation de l'effet AMR en fonction de la contrainte mécanique et de la texture cristallographique des matériaux a été développé.La stratégie de modélisation est la suivante:Trois échelles de description du comportement sont introduites: le Volume Elémentaire Représentatif (VER) polycristallin (échelle macro), le monocristal ou grain, et enfin le domaine magnétique (échelle micro).Une première étape dite de localisation permet de déterminer le chargement magnéto-mécanique (champ magnétique et contrainte mécanique) à l'échelle d'un grain en fonction du chargement extérieur appliqué. L'introduction de variables internes et des lois d'évolution correspondantes permet de décrire de façon statistique l'évolution de la microstructure en domaines magnétiques sous l'influence de ce chargement local. Toujours à cette échelle, l'utilisation du modèle phénoménologique de Doring permet, pour chaque domaine, de calculer la résistivité en fonction de l'orientation relative entre aimantation locale et courant électrique. Une fois cette résistivité locale connue, une étape dite d'homogénéisation s'appuyant sur le modèle de Bruggeman permet de déterminer la résistivité macroscopique du VER polycristallin. Il est ainsi possible de prédire la variation de la résistivité entre un état initial désaimanté et un état sous chargement magnéto-mécanique quelconque.Les résultats obtenus par cette démarche ont été comparés avec succès à des résultats expérimentaux extraits de la littérature portant sur des polycristaux de Nickel, de Fer pur ou encore de Permalloy.Ensuite des simulations reproduisant les conditions de fonctionnement des capteurs AMR ont été effectuées. Ces simulations permettent de conclure qu'il est possible d'améliorer la sensibilité des capteurs AMR en générant une contrainte résiduelle biaxiale. / The anisotropic magnetoresistance (AMR) of ferromagnetic materials is widely used as the basic phenomenon for measuring or detecting magnetic field. Owing to the relationship between magnetic domain configuration and macroscopic resistivity, the application of an external magnetic field changes the resistivity of ferromagnetic materials. Although this effect is widely used in industrial applications, some basic aspects of AMR behavior are still unsufficiently understood. For example, the role of crystallographic texture is not accurately described by conventional modeling tools. As a consequence of the direct relationship between microstructure and AMR, models for AMR effect are generally based on micromagnetic calculations. For these calculations, the number of degrees of freedom and interactions can grow exponentially when investigating macroscopic behavior (case of polycrystals for example).The thesis deals with the numerical modeling of AMR effect in ferromagnetic materials. This new 3D modeling tool can overcome this major drawback of micromagnetic approaches. A model to describe the effects of magneto-elastic coupling using a micro-macro approach is available at the laboratory GeePs. Based on the same principles of micro-macro modeling, an AMR effect simulation tool has been developed including the effect of mechanical stress and the role of crystallographic texture of materials.The modeling strategy is as follows:Three scales of description of the behavior are introduced: the Representative Volume Element (RVE) of polycrystals (macro scale), the single crystal or grain, and finally the magnetic domain (micro scale).A first step, named localization, determines the magneto-mechanical loading (magnetic field and mechanical stress) within a grain depending on the external applied load. The introduction of internal variables and corresponding evolution laws allow describing in a statistical way the evolution of the magnetic domain microstructure under the influence of the local load. Also at this scale, the use of the phenomenological Doring model allows for each area, to calculate the resistivity as a function of the relative orientation between local magnetization and electric current. Once this local resistivity is known, a so-called homogenization step based on the Bruggeman model is used to determine the macroscopic resistivity of the RVE. It is thus possible to predict the variation in resistivity between an initial demagnetized state and a state under any magneto-mechanical loading.The results obtained by this approach were successfully compared to experimental results from literature on polycrystalline nickel, pure iron or Permalloy.Then simulations reproducing AMR sensors operating conditions were carried out. These simulations lead to the conclusion that it is possible to improve the sensitivity of AMR sensors by introducing an appropriate biaxial residual stress.

L'objectif de ce travail est la mise en oeuvre de la mécanique de la rupture et d'un modèle d'endommagement afin de vérifier l'intégrité du réservoir cryotechnique d'ARIANE 5 en présence de défauts. Ces derniers sont généralement assimilés à des fissures surfaciques semi-elliptiques. L'approche globale basée sur l'analyse FAD fournit une prédiction trop conservative lorsque cette méthode utilise le facteur d'intensité K et la ténacité KIc. Ce conservatisme peut toutefois être partiellement levé à condition d'introduire la notion d'intégrale J et d'utiliser une ténacité dérivée de la courbe JR - ∆a afin d'autoriser une légère avancée de la fissure dans le cas d'un matériau ductile. Cette courbe JR - ∆a s'obtient à partir d'essais sur des éprouvettes fissurées de différentes géométries (CT, SENT et SCT). L'approche locale de la rupture constitue cependant la seule alternative possible afin de statuer sur l'acceptabilité d'un défaut suffisamment important pour générer une plasticité étendue. Cette approche nécessite alors de décrire précisément les champs de contrainte et déformation en pointe de fissure. Par conséquent il est nécessaire de tenir compte du caractère anisotrope de la plasticité, lié au procédé de laminage, lors de la modélisation du comportement élasto-plastique del'alliage d'aluminium 2219 T87. Dans cette étude un critère phénoménologique est utilisé pour modéliser la plasticité dont l'anisotropie est décrite par différents coefficients. L'optimisation de ces paramètres repose sur la comparaison entre les résultats issus de simulations numériques et ceux obtenus par des essais sur des éprouvettes lisses et entaillées. L'approche local de la rupture repose également sur la compréhension des micromécanismes de rupture. L'alliage d'aluminium 2219 T87 étudié dans le cadre de cette thèse présenteune rupture ductile dont la phase de germination s'effectue sur les inclusions Al2Cu et constitue le mécanisme prépondérant de la rupture. L'endommagement est modélisé par le modèle de Gurson modifié par Needleman et Tvergaard (GTN). L'anisotopie de rupture est simulée par deux lois de germination différentes selon le sens de laminage et le sens transverse. / The purpose of this work is to use fracture mechanics and damage model to assess the structural integrity of the ARIANE 5 tank containing defects. These flaws are usually assumed to be semi-elliptical surface cracks.This situation can be analysed using the FAD (Failure Assessment Diagram) approach. This method is based on the stress intensity factor K and a material fracture toughness (Kmat) which can be the plane strain fracture toughness (KIC) for a conservative prediction. A more realistic prediction is achieved by employing the integral J and a material parameter derived from crack resistance curve JR - ∆a . This curve is obtained by using cracked test samples with different geometries (CT, SENT and SCT). The local approach to fractureis an alternative method to assess failure in particular in the case of large scale plasticity. In this case, the previous method is not able to give accurate predictions. Analyses based on local approach concepts requirean accurate evaluation of strain and stress fields near the crack tip. It is therefore important to take plastic anisotropy into account to model the plastic behavior of the rolled sheets. In this work this is done using amacroscopic phenomenological model. The parameters of this model are determined by minimizing of the difference between the simulation results and the test measurements for smooth and notched tensile test samples. The local approach to fracture is also based on the understanding of the failure micromechanismsof materials. Ductile fracture of metallic materials can be described as a three stage process. The first is void initiation at inclusions, the second is void growth and the third is void coalescence. Nucleation of defectsat Al2Cu inclusions is the main damage mechanism. In this work, the proposed model for ductile fracture is based on the Gurson model modified by Needleman and Tvergaard (GTN). This model is extended to account for fracture anisotropic by considering two nucleation gaussian laws respectively for rolling andtransverse directions.

Cette thèse a porté sur la construction de densités d'énergie de déformation permettant de décrire le comportement non linéaire de matériaux anisotropes tels que les tissus biologiques souples (ligaments, tendons, parois artérielles etc.) ou les caoutchoucs renforcés par des fibres. Les densités que nous avons proposées ont été élaborées en se basant sur la théorie mathématique des polynômes invariants et notamment sur le théorème de Noether et l'opérateur de Reynolds. Notre travail a concerné deux types de matériaux anisotropes, le premier avec une seule famille de fibre et le second avec quatre familles. Le concept de polyconvexité a également été étudié car il est notoire qu'il joue un rôle important pour s'assurer de l'existence de solutions. Dans le cas d'un matériau comportant une seule famille de fibre, nous avons démontré qu'il était impossible qu'une densité polynomiale de degré quelconque puisse prédire des essais de cisaillement avec un chargement parallèle puis perpendiculaire à la direction des fibres. Une densité polynomiale linéaire combinée avec une fonction puissance a permis de contourner cet obstacle. Dans le cas d'un matériau comportant quatre familles de fibre, une densité polynomiale a permis de prédire correctement des résultats d'essai en traction bi-axiale extraits de la littérature. Les deux densités proposées ont été implémentées avec la méthode des éléments finis et en langage C++ dans le code de calcul universitaire FER. Pour se faire, une formulation lagrangienne totale a été adoptée. L'implémentation a été validée par des comparaisons avec des solutions analytiques de référence que nous avons exhibée dans le cas de chargements simples conduisant à des déformations homogènes. Des exemples tridimensionnels plus complexes, impliquant des déformations non-homogènes, ont également été étudiés. / This thesis has focused on the construction of strain energy densities for describing the non-linear behavior of anisotropic materials such as biological soft tissues (ligaments, tendons, arterial walls, etc.) or fiber-reinforced rubbers. The densities we have proposed have been developed with the mathematical theory of invariant polynomials, particularly the Noether theorem and the Reynolds operator. Our work involved two types of anisotropic materials, the first with a single fiber family and the second with a four-fiber family. The concept of polyconvexity has also been studied because it is well known that it plays an important role for ensuring the existence of solutions. In the case of a single fiber family, we have demonstrated that it is impossible for a polynomial density of any degree to predict shear tests with a loading parallel and then perpendicular to the direction of the fibers. A linear polynomial density combined with a power-law function allowed to overcome this problem. In the case of a material made of a four-fiber family, a polynomial density allowed to correctly predict bi-axial tensile test data extracted from the literature. The two proposed densities were implemented in C++ language in the university finite element software FER by adopting a total Lagrangian formulation. This implementation has been validated by comparisons with reference analytical solutions exhibited in the case of simple loads leading to homogeneous deformations. More complex three-dimensional examples, involving non-homogeneous deformations, have also been studied.

Cette thèse concerne la détermination des paramètres R et L de matériels bobinés à hautes fréquences par la méthode des éléments finis 3D. Les travaux se sont axés sur les transformateurs de puissance.<br />Le travail réalisé a consisté à développer une méthode permettant de s'affranchir du calcul des courants de Foucault par éléments finis. Cela a abouti au concept de perméabilité magnétique complexe. Par l'utilisation d'une formulation en potentiel scalaire magnétique et de la perméabilité complexe il a été possible de mettre en place une méthode de détermination des parmètres caractéristiques d'un transformateur de puissance. Paramètres sous forme matricielle qui seront utilisés pour des simulations de réseaux électriques.

Ce mémoire traite du rehaussement d'images sismiques tridimensionnelles, constituées d'un empilement de couches géologiques, elles-mêmes interrompues par des failles sismiques. L'objectif, en apparence contradictoire, est de lisser les couches géologiques, tout en préservant certaines discontinuités : les failles. Les méthodes de filtrage adaptatif – couramment employées dans le domaine de l'imagerie sismique – ainsi que la diffusion anisotrope sont bien adaptés à cette problématique car elles permettent une prise en compte de l'orientation des couches, ainsi que de critères indiquant la possible présence de failles. Reprenant les points clefs des approches précédentes, nous proposons une équation d'évolution fondée sur la définition de trois zones de comportement différentié à l'intérieur du voisinage de chaque point. La première contient les points appartenant à la même couche géologique que le point central, la deuxième ceux correspondant aux couches voisines et la troisième les points incertains. La définition de ces zones est réalisée en chaque point à l'aide de l'orientation estimée des couches géologiques. Une évolution de cette méthode consiste à sélectionner, dans le voisinage défini précédemment, la sous-partie qui présente le moins de chances d'être traversée par une faille. Cette nouvelle approche permet non seulement de lisser l'image, mais également de renforcer la visibilité des failles. Les performances des méthodes proposées sont comparées à celles des approches classiques à l'aide d'un estimateur objectif de gain de qualité, employé sur des images de synthèse. Une comparaison visuelle de résultats obtenus sur des données réelles est également réalisée.

Aujourd'hui, la propagation de la lumière dans les<br />nano-matériaux diélectriques complexes est un sujet de recherche<br />riche et fascinant, tant pour ses implications fondamentales que<br />pour son impact technologique. Dans cette thèse, nous étudions les<br />effets d'interférence de la lumière dans les systèmes photoniques<br />quasi-ordonnés.<br /><br />Dans des milieux diélectriques aléatoires, on peut décrire le<br />mouvement des photons mutiplement diffusés par une marche<br />aléatoire de diffuseur en diffuseur: la plupart des effets<br />d'interférence se moyennent alors à zéro, mais certains survivent<br />quand même au désordre et induisent des phénomènes non-triviaux.<br />Dans des milieux qui diffusent très fortement, la lumière pourrait<br />même devenir localisée et aucun transport ne serait possible. Dans<br />les milieux ordonnés, la périodicité conduit à des lois de<br />dispersion inhabituelles où les effect collectifs d'interférence<br />dominent: le transport est fortement dépendant de la fréquence, il<br />peut être sensiblement augmenté (interférences constructives) ou<br />complètement inhibé (interférences destructives).<br /><br />Notre compréhension de la propagation des ondes lumineuses dans<br />les milieux ordonnés et désordonnés augmente rapidement, mais le<br />comportement dans le régime intermédiaire entre les deux extrèmes<br />-- ordre parfait et désordre complet -- est mal compris. Les<br />systèmes quasi-ordonnés brisent la symétrie de rotation ou de<br />translation et présentent des formes nouvelles et<br />non-conventionnelles de transport de la lumière. Les milieux<br />aléatoires qui diffusent anisotropiquement et les cristaux<br />liquides nématiques, les formes spéciales de cristaux photoniques<br />et les quasi-cristaux photoniques de Fibonacci sont des exemples<br />de systèmes quasi-ordonnés que nous étudions dans cette thèse.<br /><br /><br />Que se passe-t-il si une direction préférentielle de diffusion ou<br />un axe préférentiel de polarisabilité est présent dans un milieu<br />aléatoire ?<br />Comment la propagation de la lumière est elle modifiée dans une<br />structure périodique si un potential optique est superposé à la<br />structure cristalline ?<br /><br />Dans cette thèse nous essayerons de répondre à ces questions, avec<br />des arguments théoriques, des simulations numériques et des<br />résultats expérimentaux.

Ce travail concerne la simulation d'écoulements viscoélastiques compressibles appliquée à l'injection de polymères. La compressibilité est intégrée dans Rem3D en supposant que la densité du matériau suit une loi d'évolution du type loi de Tait. La conservation de la masse est écrite comme une équation en vitesse, pression et température, à travers des coefficients de compressibilité isotherme et de dilatation isobare. Le système obtenu est désigné "Stokes compressible" et sa résolution numérique est faite par la méthode des éléments finis mixtes. Le système obtenu est non-linéaire et non-symétrique. Le couplage thermique et l'extension à des problèmes avec surface libre sont aussi consid ér és. Le modèle viscoélastique choisi est le modèle Pom-Pom, issu de la dynamique moléculaire. L'extra-contrainte est fonction des propriétés microscopiques du matériau, comme l'orientation moléculaire et son étirement. L'élasticité est vue comme une perturbation dans le problème mécanique, et une méthode de stabilisation du type DEVSS est utilisée. L'orientation et l'étirement sont déterminés par la résolution de deux équations d'évolution via une méthode espace-temps Galerkin discontinu. Finalement, la thermoviscoélasticité est abordée brièvement. Dans le contexte de l'injection de polymères, REM3D couvre aujourd'hui toutes les phases du proc éd é. Néanmoins, la solidification et la transition liquide-solide sont approximées par un comportement du type liquide de très haute viscosité. L'introduction de la compressibilité permet de compenser le retrait du matériau par un apport supplémentaire de matière. D'un autre côté, la prise en compte d'un comportement viscoélastique détecte d'éventuelles anisotropies des propriétés de la pièce injectée. Les diverses comparaisons des résultats obtenus avec la littérature et l'expérience montre une bonne concordance, validant les modèles implémentés.

L'objet de ce travail est de présenter une partie des travaux entrepris au laboratoire CELIA (CEA, CNRS, Université Bordeaux I) dans le domaine de la modélisation numérique des écoulements fortement compressibles. Cette activité au sein de l'équipe Interaction-Fusion par Confinement Inertiel-Astrophysique, a eu pour objectif principal la mise au point et le développement de schémas numériques robustes dédiés à la simulation numérique des plasmas à haute densité d'énergie appliquée à la production d'énergie par fusion. Ces travaux se sont concrétisés par l'écriture du code CHIC (Code d'Hydrodynamique et d'Implosion du CELIA), logiciel permettant de concevoir et de restituer des expériences dans le domaine de la Fusion par Confinement Inertiel (FCI). Le modèle théorique numérique décrivant l'implosion d'une cible laser est un système d'équations aux dérivées partielles au centre duquel on trouve les équations d'Euler écrites dans le formalisme lagrangien, couplées à des équations de diffusion non linéaires modélisant le transport de l'énergie par les électrons et les photons. Dans cet exposé, après un bref rappel du contexte physique, nous décrirons les deux méthodes originales qui constituent l'ossature numérique du code CHIC. Il s'agit de deux schémas numériques d'ordre élevé du type volumes finis dédiés respectivement à la résolution des équations de l'hydrodynamique lagrangienne et à la résolution d'équations de diffusion anisotrope sur des maillages bi-dimensionnels non-structurés. Le premier schéma, dénommé EUCCLHYD (Explicit Unstructured Lagrangian HYDrodynamics), permet de résoudre les équations de la dynamique des gaz sur un maillage mobile qui se déplace à la vitesse du fluide. Il est obtenu à partir d'un formalisme général basé sur le concept de forces de sous-mailles. Dans ce cadre, les flux numériques sont exprimés en fonction des forces de sous-mailles et de la vitesse des noeuds. Leur détermination repose sur les trois principes fondamentaux suivants : compatibilité géométrique entre le mouvement des noeuds et la variation de volume des mailles (loi de conservation géométrique), compatibilité avec le second principe de la thermodynamique et conservation de l'énergie totale et de la quantité de mouvement. L'extension de ce schéma à l'ordre deux est mise en place à l'aide d'une méthode basée sur la résolution d'un problème de Riemann généralisé dans l'approximation acoustique. Le second schéma, appelé CCLAD (Cell-Centered LAgrangian Diffusion), concerne la résolution de l'équation de la chaleur anisotrope non-linéaire. La discrétisation correspondante s'appuye sur une formulation variationnelle locale au niveau des sous-mailles qui permet de construire une approximation multi-points du flux de chaleur. Cette discrétisation d'ordre élevé rend possible la résolution des équations de la diffusion anisotrope avec une précision satisfaisante sur des maillages lagrangiens fortement déformés. La précision et la robustesse de ces méthodes numériques sont démontrées sur des cas-tests représentatifs.

Le procédé de forgeage à froid permet d'obtenir des pièces de géométrie complexe et présentant une très bonne précision dimensionnelle. La forme finale est obtenue en une ou plusieurs passes, et la matière est soumise à des chargements à la fois multiaxiaux et non monotones accompagnés de forts niveaux de déformation plastique qui peuvent l'amener à s'endommager : rupture de la pièce, fissures en surface ou à cœur, ... Il est alors tout à fait primordial de connaître et mieux comprendre les mécanismes qui peuvent produire ces défauts. Le sujet de thèse concerne donc la modélisation numérique de l'endommagement en forge à froid. Dans un premier temps, le contexte bibliographique concernant les mécanismes d'endommagement et leurs modélisations ont été étudiés. Un modèle d'endommagement, basé sur le modèle de Lemaitre, a alors été implémenté dans le logiciel éléments finis Forge2005®. Dans un deuxième temps, une étude expérimentale concernant la caractérisation de deux nuances en forge à froid a été réalisée. Différents essais de traction ont été menés afin d'étudier notamment l'influence de l'orientation des inclusions sur l'endommagement. Des essais de compression/traction ont également été réalisés dans le but d'examiner l'influence éventuelle de la pré-compression sur la matière. Enfin des essais in-situ (MEB et tomographie X) ont été menés pour analyser l'endommagement au cœur de la matière et mieux en comprendre les mécanismes. Sur la base de ces observations expérimentales, plusieurs améliorations ont été apportées au modèle de Lemaitre afin de le rendre plus adapté aux chargements complexes. Parmi ces améliorations, nous pouvons citer par exemple la prise en compte de l'endommagement en compression, l'introduction d'une limite en triaxialité des contraintes négative, l'introduction d'un nouveau potentiel d'endommagement et la prise en compte de l'orientation des inclusions. Pour finir, des simulations numériques ont été réalisées par le biais du logiciel Forge2005® afin de tester le modèle d'endommagement implémenté. Des cas de validations portant sur des essais mécaniques simples (traction, TELE) ainsi que sur des cas industriels ont alors été simulés et comparés aux résultats expérimentaux.