Dans le contexte actuel de transition énergétique, les acteurs de la production électrique (dont EDF) s'orientent de plus en plus vers les technologies "bas carbone", permettant de développer leur mix énergétique et d'assurer une production toujours plus respectueuse de l'environnement. Parmi l'ensemble des alternatives progressivement industrialisées au cours de ces dernières décennies, l'énergie éolienne voit son développement s'accélérer. De nouveaux projets voient le jour, notamment avec le développement de parcs éoliens posés en mer le long des côtes françaises.Afin de garantir une sécurité maximale des installations vis-à-vis des défaillances mécaniques, les ingénieurs doivent mettre en place un ensemble de vérifications ayant pour but d'assurer l'intégrité structurelle pendant une durée de vie donnée. Cette étape, nécessaire à la certification et donc à l'industrialisation du projet, exige un investissement numérique important. Dans le cas particulier de l'estimation de l'endommagement à durée de vie, un ensemble complet d'analyses doit être effectué afin de simuler les réactions structurelles en différents points de la conception et ce, pour l'ensemble des conditions environnementales de chargement possibles (cas de charge). Au total, quelques milliers de simulations correspondant à l'ensemble des situations probables sont à prévoir pour la seule estimation de l'endommagement à durée de vie. De plus, la prise en compte des différents phénomènes physiques du problème (aérodynamique, hydrodynamique, mécanique élastique), ainsi que la considération des non-linéarités des réponses liées à l'asservissement de la turbine, impliquent l'utilisation de codes de résolution temporels et multi-physiques coûteux pour chacune des situations de chargement à simuler.Afin de réduire l'investissement lié à l'estimation de cet endommagement, l'utilisation de méthodes numériques de substitution (également appelées métamodèles) est une alternative prometteuse. Ces méthodologies ont montré leur efficacité dans divers domaines de l'ingénierie permettant d'approcher, par diverses hypothèses statistiques, la réponse de modèles numériques en ne considérant qu'un nombre réduit de simulations.Ce travail de thèse s'est focalisé sur le développement d'outils numériques pour le dimensionnement efficace des structures éoliennes en mer et plus particulièrement sur l'utilisation des métamodèles de Krigeage (ou méthode de régression par processus gaussien) pour l'estimation de l'endommagement mécanique. Sous certaines hypothèses, cette méthode de substitution se prête particulièrement bien à l'approximation des réponses de modèles physiques, notamment grâce à la mise en place facilitée d'architectures itératives d'enrichissement ("active learning strategy"). Deux domaines ont principalement été étudiés : l'estimation rapide des quantités d'endommagement structurel par le développement de la méthode "Adaptive Kriging for Damage Assessment" (AK-DA) et les possibilités d'utilisation de cette dernière dans le cadre d'études de la fiabilité conceptuelle au regard de la tenue à l'endommagement. L'applicabilité de ces méthodes a été illustrée via l'utilisation d'exemples numériques inspirés du contexte industriel et de ses contraintes. / The mechanical certification of wind turbine structures is required for the funding of newoffshore projects on the French coasts. In order to ensure a maximal safety level of installations,a series of structural analyzes are required by the certification bodies. Amongst all, thedamage based computations represent an important numerical effort for EDF. The presentedworks focus on the applicability and the performances of Kriging metamodels for the estimationof the lifetime cumulated damage of offshore wind turbine structures (AK-DA approach)and the damage based reliability assessment of new designs (AK-MCS/AK-DA coupling).

Le 32CrMoV13 nitruré est un matériau à gradient de propriétés utilisé pour la réalisation de bagues de roulements aéronautiques. Il présente l'avantage d'être dur en surface et ductile à cœur, facilitant ainsi la lutte contre les effets de la pollution des contacts par les lubrifiants. De plus, afin d'améliorer les performances des roulements, la tendance est à l'utilisation de céramique (nitrure de silicium Si3N4) pour les éléments roulants à la place de l'acier (M50). Ainsi, le travail de thèse proposé ici est l'étude du comportement en fatigue de contact en présence d'indentations pour deux couples de matériaux (32CrMoV13 nitruré / M50 et 32CrMoV13 nitruré / Si3N4). Ce comportement est principalement orienté en terme de durée de vie du contact indenté par l'intermédiaire d'un critère d'amorçage de fissures. Une étude des propriétés mécaniques du 32CrMoV13 nitruré, nécessaires à l'application du critère, a été réalisée. La loi d'écrouissage et l'ouverture micro-plastique cyclique de la couche nitrurée du 32CrMoV13 nitruré ont été identifiées. Puis, par l'intermédiaire de défauts géométrique de type trous cylindriques, les propriétés d'amorçage de fissures ont été étudiées lors d'essais de compression uniaxiale. Le comportement en fatigue de contact en présence d'indentations a été étudié sur machine bi-disques. Ces essais ont permis de caractériser la différence de comportement du 32CrMoV13 nitruré, et indenté, dans un contact acier / acier et acier / céramique. L'effet du glissement a notamment été mis en évidence. Ces études expérimentales ont été couplées avec des simulations numériques afin d'étalonner le modèle physique d'amorçage de fissures et d'en valider sa bonne représentativité. En conclusion un modèle de prédiction de durée de vie en fatigue de roulement a été proposé. Il est particulièrement adapté aux cas où des contraintes résiduelles sont présentes, par exemple suite à l'indentation des surfaces. Le modèle s'appuie sur la caractérisation mécanique des couches traitées thermo-chimiquement de l'acier 32CrMoV13 nitruré.

Le 32CrMoV13 nitruré est un matériau à gradient de propriétés utilisé pour la réalisation de bagues de roulements aéronautiques. Il présente l’avantage d’être dur en surface et ductile à cœur, facilitant ainsi la lutte contre les effets de la pollution des contacts par les lubrifiants. De plus, afin d’améliorer les performances des roulements, la tendance est à l’utilisation de céramique (nitrure de silicium Si3N4) pour les éléments roulants à la place de l’acier (M50). Ainsi, le travail de thèse proposé ici est l’étude du comportement en fatigue de contact en présence d’indentations pour deux couples de matériaux (32CrMoV13 nitruré / M50 et 32CrMoV13 nitruré / Si3N4). Ce comportement est principalement orienté en terme de durée de vie du contact indenté par l’intermédiaire d’un critère d’amorçage de fissures. Une étude des propriétés mécaniques du 32CrMoV13 nitruré, nécessaires à l’application du critère, a été réalisée. La loi d’écrouissage et l’ouverture micro-plastique cyclique de la couche nitrurée du 32CrMoV13 nitruré ont été identifiées. Puis, par l’intermédiaire de défauts géométrique de type trous cylindriques, les propriétés d’amorçage de fissures ont été étudiées lors d’essais de compression uniaxiale. Le comportement en fatigue de contact en présence d’indentations a été étudié sur machine bi-disques. Ces essais ont permis de caractériser la différence de comportement du 32CrMoV13 nitruré, et indenté, dans un contact acier / acier et acier / céramique. L’effet du glissement a notamment été mis en évidence. Ces études expérimentales ont été couplées avec des simulations numériques afin d’étalonner le modèle physique d’amorçage de fissures et d’en valider sa bonne représentativité. En conclusion un modèle de prédiction de durée de vie en fatigue de roulement a été proposé. Il est particulièrement adapté aux cas où des contraintes résiduelles sont présentes, par exemple suite à l’indentation des surfaces. Le modèle s’appuie sur la caractérisation mécanique des couches traitées thermo-chimiquement de l’acier 32CrMoV13 nitruré. / The nitrided 32CrMoV13 steel is a gradient properties material used for the production of aeronautical bearing rings. The principal advantages are its surface hardness and core ductility, useful against debris denting. Moreover, to improve rolling bearing performances, ceramic balls (silicon nitride Si3N4) instead of steel (M50). The present study aims to investigate the contact fatigue performance in presence of indentation for two pairs of materials (nitrided 32CrMoV13 steel / M50 and nitrided 32CrMoV13 steel / Si3N4). This performance is mainly oriented in terms of lifetime of the indented contact by a crack initiation criterion. A study of the mechanicals properties of the nitrided 32CrMoV13 steel, required for the application of the criterion, is then performed. It allows to characterize the hardening law and the cyclic micro-plastic behavior of the nitrided layer of the nitrided 32CrMoV13 steel. Crack initiation properties of the nitrided surface layer have been evaluated by compression tests on thin hollow cylindrical specimens with small cylindrical holes. The contact fatigue behavior in presence of indentations was studied on a high-speed twin-disks machine. These tests allowed to characterize the difference in behavior of nitrided 32CrMoV13 steel (indented) in a steel-to-steel contact or steel-to-ceramic contact. The major effect of sliding has been highlighted. All these experimental studies were coupled with numerical simulations to calibrate the physical model of crack initiation and to show its capabilities. In conclusion a new model of rolling contact fatigue life has been proposed. It is specifically adapted to operation in presence of residual stresses such as resulting from debris denting. The model uses mechanical properties that have been identified for thermo-chemically treated layers of nitrided 32CrMoV13 steel.

Lorsqu'on découpe un feuillard à l'aide d'un outil, ou lorsqu'on découpe un tube selon son axe, au fur et à mesure que l'on propage la fissure qui traduit la découpe il arrive que des ondulations de flambage perturbent les deux bords libres générés par la propagation de la fissure. Cette étude vise à analyser les origines de ces ondulations. Nous avons mené une campagne expérimentale, dans laquelle des tubes en acier inox avec différentes géométries (rayon/épaisseur) sont " découpés " selon une génératrice. Une instrumentation adéquate, couplant des mesures ponctuelles, à l'aide de jauges de déformation, et une méthode champ par corrélation d'image, nous a permis de correctement mettre en exergue la phénoménologie, en particulier les cinématiques induites à l'échelle géométrique de la fissure (front de fissure) ainsi qu'à l'échelle du tube, avec les longueurs d'onde de flambage observées à l'aval de la fissure. La modélisation numérique menée en non linéaire géométrique (flambage), matériau (déchirure ductile), et conditions aux limites (contact) est aussi abordée à l'aide du code de calcul Abaqus/Standard. Pour la gestion de la propagation de la fissure, deux modèles de rupture sont proposés. Le premier modèle dit zone cohésive est développé et implanté dans le code Abaqus via la subroutine UEL. Pour la deuxième modélisation, nous avons utilisé le modèle dit " d'endommagement ductile " du code Abaqus. La modélisation via des éléments massifs ou des éléments coques volumiques ainsi que l'utilisation de ces modèles de rupture permettent de corroborer les observations expérimentales. Ces travaux montrent que l'augmentation de la charge inhérente au déplacement de l'outil de " découpe ", induit une extension dans la direction circonférentielle et donc une striction dans la direction radiale amenant finalement la rupture. Lors de la rupture, un " sillage plastique " apparait, relativement large, près et parallèle aux bords de la fissure. " Confiné " par les autres parties du tube qui restent élastiques, des contraintes de compression axiale résiduelles apparaissent dans ce sillage plastique, à l'aval de la fissure, leur intensité est suffisante pour produire les ondulations des bords libres qui traduisent un flambage local. Les contraintes résiduelles liées à l'opération de découpe induisent donc le flambage.

La corrosion sous contrainte (CSC) est un phénomène synergique d'endommagement qui résulte d'un processus de corrosion (dissolution, adsorption) et d'une rupture mécanique (fissuration). Les mécanismes de couplage mécano-électrochimique en pointe de fissure nécessaires à la compréhension du phénomène sont encore mal connus puisqu'ils dépendent du système d'étude (métal/milieu agressif) et font intervenir de nombreux facteurs mécaniques et électrochimiques. Dans cette thèse, nous nous proposons d'étudier les interactions réciproques entre la dissolution et l'état de contrainte mécanique en pointe de fissure (facteurs d'intensité des contraintes) pour le cas du Zircaloy-4 en milieu aqueux halogéné. Les éprouvettes sont d'abord pré-fissurées par fatigue à l'air avec la technique du Load-Shedding, ce qui permet alors de maîtriser le facteur d'intensité des contraintes résiduel de pré-fissuration par fatigue. Ensuite, une pré-oxydation thermique est réalisée pour produire une couche de film passif sur les surfaces de l'éprouvette. Les réactions électrochimiques sont alors concentrées en pointe de fissure qui induit également une concentration des effets mécaniques. Des techniques de corrélation d'images sont développées dans le but d'identifier les facteurs d'intensité des contraintes et d'estimer la longueur de fissure en temps réel. Ainsi des essais originaux de CSC, permettant d'imposer les facteurs d'intensité en pointe de fissure, sont conduits et les résultats montrent alors clairement les effets synergiques entre l'évolution des mesures mécaniques et des courants de dissolution. D'autre part, l'existence d'un facteur d'intensité des contraintes seuil de propagation en CSC, présentant une forte dépendance à l'histoire de chargement et à la variation de l'état de contrainte locale en pointe de fissure (variation des facteurs d'intensité des contraintes), est mise en évidence. Ceci montre alors que la plasticité doit évoluer pour que la dissolution se produise. Ainsi, le temps caractéristique d'arriver de nouvelles dislocations en pointe de fissure ne doit pas dépasser la durée caractéristique des réactions de passivation. Enfin sur la base des résultats expérimentaux, un modèle de loi de propagation phénoménologique de fissure en CSC est proposé et ses paramètres sont identifiés et validés à partir de différents essais de CSC.

L'objectif de ce travail est de proposer une méthodologie simple pour évaluer l'intégrité et la durée de vie d'échangeurs thermiques soudés. Une approche à deux échelles est proposée. Une description macroscopique avec la prise en compte de la structure de l'échangeur est menée pour permettre des calculs thermoélastiques par éléments finis. La réponse mécanique de l'échangeur pour deschargements thermiques, cycliques, simples est évaluée. Notamment, les zones de concentration decontraintes sont repérées. A partir cette étude, une étude micromécanique du comportement dumatériau composant l'échangeur est menée. Le matériau considéré est un acier 316L. Soncomportement élastoplastique est identifié avec un écrouissage isotrope et cinématique. La tenuemécanique pour des chargements en fatigue oligocyclique est évaluée à l'aide d'un dispositif deflexion 4 points alternée et un critère de Manson-Coffin est identifié. Ce critère est utilisé pour évaluerle nombre de cycle admissible par l'échangeur pour une amplitude de chargement donnée. Pour cela,la déformation plastique attendue dans l'échangeur est évaluée à partir d'une équivalence en énergieaux endroits où la contrainte se concentre. Les prédictions du modèle ont été comparées de manièresatisfaisante avec les résultats expérimentaux menés sur un échangeur test, pour la réponsethermoélastique que pour l'évaluation du nombre de cycles à rupture.

Lorsqu’on découpe un feuillard à l’aide d’un outil, ou lorsqu’on découpe un tube selon son axe, au fur et à mesure que l’on propage la fissure qui traduit la découpe il arrive que des ondulations de flambage perturbent les deux bords libres générés par la propagation de la fissure. Cette étude vise à analyser les origines de ces ondulations. Nous avons mené une campagne expérimentale, dans laquelle des tubes en acier inox avec différentes géométries (rayon/épaisseur) sont « découpés » selon une génératrice. Une instrumentation adéquate, couplant des mesures ponctuelles, à l’aide de jauges de déformation, et une méthode champ par corrélation d’image, nous a permis de correctement mettre en exergue la phénoménologie, en particulier les cinématiques induites à l’échelle géométrique de la fissure (front de fissure) ainsi qu’à l’échelle du tube, avec les longueurs d’onde de flambage observées à l’aval de la fissure. La modélisation numérique menée en non linéaire géométrique (flambage), matériau (déchirure ductile), et conditions aux limites (contact) est aussi abordée à l’aide du code de calcul Abaqus/Standard. Pour la gestion de la propagation de la fissure, deux modèles de rupture sont proposés. Le premier modèle dit zone cohésive est développé et implanté dans le code Abaqus via la subroutine UEL. Pour la deuxième modélisation, nous avons utilisé le modèle dit « d’endommagement ductile » du code Abaqus. La modélisation via des éléments massifs ou des éléments coques volumiques ainsi que l’utilisation de ces modèles de rupture permettent de corroborer les observations expérimentales. Ces travaux montrent que l’augmentation de la charge inhérente au déplacement de l’outil de « découpe », induit une extension dans la direction circonférentielle et donc une striction dans la direction radiale amenant finalement la rupture. Lors de la rupture, un « sillage plastique » apparait, relativement large, près et parallèle aux bords de la fissure. « Confiné » par les autres parties du tube qui restent élastiques, des contraintes de compression axiale résiduelles apparaissent dans ce sillage plastique, à l’aval de la fissure, leur intensité est suffisante pour produire les ondulations des bords libres qui traduisent un flambage local. Les contraintes résiduelles liées à l’opération de découpe induisent donc le flambage. / When one uses a tool to cut a sheet metal, or a tube according to his axis, as one propagates the crack which translates cutting it arrives that undulations of buckling disturb the two free edges generated by the propagation of the crack. This study aims at analyzing the origins of this behavior. We conducted an experimental campaign, in which stainless steel tubes with various geometries (radius/thickness) are « cut out » according to a generator. An adequate instrumentation, coupling of specific measurements, using strain gauges, and a field method, by digital image correlation, allowed us accurately to put forward phenomenology, in particular the kinematics at the scale of the crack (ahead of crack tip) and at the level of tube, with the wavelengths of buckling observed at the downstream of the crack tip. The numerical modeling taking into account nonlinearities of material (ductile tear), geometry (buckling) and boundary conditions (contact) is also approached using the code Abaqus/Standard. For the management of the crack propagation, two rupture models are proposed. The first model called cohesive zone is developed and implemented in the Abaqus code via the user routine UEL. For the second modeling, we used the model called “ductile damage model” in the Abaqus code. Modeling via solid elements or shell continuum elements as well as the use of these rupture models make it possible to corroborate the experimental observations. These studies show that the increase of the load inherent in the displacement of the tool of « cutting » induced a circumferential extension of the tube that leads to a local necking in the radial direction bringing the rupture finally. During the failure, a “plastic wake” appears, relatively wide, close and parallel to crack lips. Constrained by other parts of the tube which remain elastic, sufficient axial residual compressive stresses produced in this plastic wake produce the undulations which represents a local buckling. The residual stresses related to the operation of cutting thus induce buckling.

L'objectif de ce travail est la caractérisation de l’endommagement thermique et mécanique dans le mortier par les ondes acoustiques non linéaires. La corrélation entre les paramètres acoustiques linéaires et non-linéaires étudiée est basée sur les essais expérimentaux et la modélisation. Pour l’endommagement thermique, des éprouvettes cylindriques ont été préparées et caractérisées par l'étude de la porosité et la saturation. Ensuite, l'acoustique linéaire et l’acoustique non linéaire (génération d'harmoniques) ont été appliquées afin de quantifier l’endommagement. Les essais acoustiques linéaires ont prouvé que les vitesses transversales, longitudinales et le module d'Young du mortier diminue en fonction de la température. Les essais acoustiques non linéaires ont montré l'augmentation du bêta fonction de la température. Pour l’endommagement mécanique, le phénomène d'autocicatrisation est suivi par la perméabilité et les essais acoustiques. Les essais de perméabilité ont montré que le débit d'air et la taille de la fissure diminue rapidement au cours du premier mois, puis lentement durant le reste du processus d'autocicatrisation. D'autre part, les tests acoustiques non linéaires ont montré que « alpha » et « bêta » diminuent durant le processus de l’autocicatrisation qui signifie que les paramètres non linéaires sont un bon indicateur pour caractériser ce phénomène. A partir des résultats expérimentaux, une corrélation polynomiale de 2ème degré a été établie entre les paramètres non linéaires et l’index d’endommagement. Les résultats de ce travail représentent un bon départ pour étudier le phénomène de l'autocicatrisation par les ondes acoustiques non linéaires / The objective of this work is the characterization of heat and mechanical damage in the mortar by the nonlinear acoustic waves. The correlation between non-linear/linear acoustic parameters and damage in mortar is studied based on experiments and modelling. For the heat damage, cylindrical specimens were prepared and were characterized by studying the porosity and saturation. Indeed, the linear acoustic (UPV) and non-linear acoustic (Higher harmonic generation) were applied to characterize the damage. The linear acoustic tests have shown that velocities and modulus of Young of the mortar decreases in function of the temperature. The non-linear acoustic tests have shown that beta increases in function of the temperature.For the mechanical damage. The self-healing phenomenon was characterized by the permeability and the acoustic tests. Indeed, the permeability tests have shown that the airflow and the crack size decreases quickly in the first month then slowly for the rest of the self-healing process. On the other hand, the non-linear acoustic tests shown that the alpha and beta decreases according to the self healing process which means that the nonlinear parameters are a good indicators to characterize the self-healing. Moreover, the analysis of the experimental results indicates that the frequency resonant technique is more efficient to characterize the defects in the mortar than the higher harmonic generation. A polynomial correlations of a 2nd degree was established between the nonlinear parameters and the index damage. The findings of this work should be most appropriate as a foundation for the study of the self healing by the nonlinear acoustic waves

L'objectif de ce travail est de proposer un modèle de comportement avec endommagement ductile pour la simulation des procédés de mise en forme de tôles minces qui peut bien représenter le comportement des matériaux sous des trajets de chargement complexes en grandes déformations plastiques. Basé sur la thermodynamique des processus irréversibles, les équations de comportement couplé à l’endommagement tiennent compte des anisotropies initiales et induites, de l’écrouissage isotrope et cinématique et de l’endommagement isotrope ductile. Les effets de fermeture des microfissures, de triaxialité des contraintes et de l’angle de Lode sont introduits pour influencer l’évolution de l’endommagement sous une large gamme de triaxialité des contraintes. La distorsion de la surface de charge est introduite via un tenseur déviateur qui gouverne la distorsion de la surface de charge. A des fins de comparaison, les courbes limites de formage sont tracées basées sur l’approche M-K.Des essais sont conduits sur trois matériaux pour les besoins d’identification et de validation des modèles proposés. L’identification utilise un couplage entre le code ABAQUS et un programme MATLAB via un script en langage Python. Après l’implémentation numérique du modèle dans ABAQUS/Explicite et une étude paramétrique systématique, plusieurs procédés de mise en forme de structures minces sont simulés. Des comparaisons expériences-calculs montrent les performances prédictives de la modélisation proposée / The objective of this work is to propose a “highly” predictive material model for sheet metal forming simulation which can well represent the sheet material behavior under complex loading paths and large plastic strains. Based on the thermodynamics of irreversible processes framework, the advanced fully coupled constitutive equations are proposed taking into account the initial and induced anisotropies, isotropic and kinematic hardening as well as the isotropic ductile damage. The microcracks closure, the stress triaxiality and the Lode angle effects are introduced to influence the damage rate under a wide range of triaxiality ratios. The distortion of the yield surface is described by replacing the usual stress deviator tensor by a ‘distorted stress’ deviator tensor, which governs the distortion of the yield surfaces. For comparisons, the FLD and FLSD models based on M-K approach are developed.A series of experiments for three materials are conducted for the identification and validation of the proposed models. For the parameters identification of the fully coupled CDM model, an inverse methodology combining MATLAB-based minimization software with ABAQUS FE code through the Python script is used. After the implementation of the model in ABAQUS/Explicit and a systematic parametric study, various sheet metal forming processes have been numerically simulated. At last, through the comparisons between experimental and numerical results including the ductile damage initiation and propagation, the high capability of the fully coupled CDM model is proved

L’objectif de ce travail est de développer une modélisation prédictive du comportement et de la rupture ductile des matériaux métalliques à anisotropies initiales et induites par l’endommagement. La thermodynamique des processus irréversibles est utilisée comme cadre pour la formulation proposée. Le modèle de comportement est élastoplastique anisotrope avec écrouissage non linéaire isotrope et cinématique en grandes déformations plastiques, avec une théorie non associée à normalité associée, basée sur des normes de contraintes quadratiques. L’endommagement ductile anisotrope est décrit par un tenseur du second ordre symétrique dont l’évolution est décrites par des relations de type Lemaitre-Desmorat. Le couplage fort comportement-endommagement est réalisé dans le cadre de l’hypothèse de l’équivalence en énergie totale où l’effet de l’endommagement sur le comportement est introduit par un tenseur « effet d’endommagement » d’ordre quatre symétrique de type Murakami. Après une caractérisation expérimentale des mécanismes physiques de l’endommagement dans l’acier AISI 316L, le modèle de comportement avec endommagement a été identifié. Une fois discrétisé et implémenté dans le code de calcul de structures ABAQUS/Explicit®, une étude paramétrique et de nombreuses simulations numériques de l’endommagement anisotrope en mise en forme de quelques structures ont été réalisées et discutées en détail / The objective of this work is to develop a predictive modeling of behavior and ductile fracture of metallic materials with initial anisotropy and induced by the ductile anisotropic damage. Thermodynamics of irreversible processes is used as a framework for the proposed formulation. The model is anisotropic elastoplastic with non-linear isotropic and kinematic hardening under large plastic strains. It is formulated in the framework of the non-associative plasticity theory with associative normality rule and based on quadratic equivalent stress. The anisotropic ductile damage is described by a symmetric second-rank tensor whose evolution is described by Lemaitre /Desmorat type relationships. The strong damage-behavior coupling is done under the assumption of total energy equivalence where the effect of the anisotropic damage is introduced by a fourth-rank symmetric damage-effect tensor of Murakami kind.After an experimental characterization of the main physical mechanisms of anisotropic damage in stainless steel AISI 316L, the behavior model with damage has been identified. Once discretized and implemented in the computer code ABAQUS / Explicit ®, a parametric study and many numerical simulations of anisotropic damage in some metal forming processes have been carried out and discussed in detail