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.

We introduce a mathematical framework for the inverse medium problem arising commonly in geotechnical site characterization and geophysical probing applications, when stress waves are used to probe the material composition of the interrogated medium. Specifically, we attempt to recover the spatial distribution of Lame's parameters ( and μ) of an elastic semi-infinite arbitrarily heterogeneous medium, using surface measurements of the medium's response to prescribed dynamic excitations. The focus is on characterizing near-surface deposits, and to this end, we develop a method that is implemented directly in the time-domain, is driven by the full waveform response collected at receivers on the surface, while the domain of interest is truncated using Perfectly-Matched-Layers (PMLs) to limit the originally semi-infinite extent of the physical domain. There are two key issues associated with the problem at hand: (a) the forward problem, namely the numerical simulation of the wave motion in the domain of interest; and (b) the framework and strategies for tackling the inverse problem. To address the forward problem, it is necessary that the domain of interest be truncated, and the resulting finite domain be forced to mimic the physics of the original problem: to this end, we introduce unsplit-field PMLs, and develop and implement two new formulations, one fully-mixed and one hybrid (mixed coupled with a non-mixed approach) that model wave motion within the, now PML-truncated, domain. To address the inverse problem, we adopt a partial-differential-equation-constrained optimization framework that results in the usual triplet of an initial-and-boundary-value forward problem, a final-and-boundary-value adjoint problem, and a time-independent boundary-value control problem. This triplet of boundary-value-problems is used to guide the optimizer to the target profile of the spatially distributed Lame parameters. Given the multiplicity of solutions, we assist the optimizer, by deploying regularization schemes, continuation schemes (regularization factor and source-frequency content), as well as a physics-driven simple procedure to bias the search directions. We report numerical examples attesting to the quality, stability, and efficiency of the forward wave modeling. We also report moderate success with numerical experiments targeting inversion of both smooth and sharp profiles in two dimensions. / text

This thesis investigates the numerical simulation of three-dimensional, mechanical deformation problems in the context of large deformations. The main focus lies on the prediction of non-linearly elastic, incompressible material. Based on the equilibrium of forces, we present the weak formulation of the large deformation problem. The discrete version can be derived by using linearisation techniques and an adaptive mixed finite element method. This problem turns out to be a saddle point problem that can, among other methods, be solved via the Bramble-Pasciak conjugate gradient method or the minimal residual algorithm. With some modifications the resulting simulation can be improved but we also address remaining limitations. Some numerical examples show the capability of the final FEM software. In addition, we briefly discuss the special case of linear elasticity with small deformations. Here we directly derive a linear weak formulation with a saddle point structure and apply the adaptive mixed finite element method. It is shown that the presented findings can also be used to treat the nearly incompressible case.

This thesis investigates the numerical simulation of three-dimensional, mechanical deformation problems in the context of large deformations. The main focus lies on the prediction of non-linearly elastic, incompressible material. Based on the equilibrium of forces, we present the weak formulation of the large deformation problem. The discrete version can be derived by using linearisation techniques and an adaptive mixed finite element method. This problem turns out to be a saddle point problem that can, among other methods, be solved via the Bramble-Pasciak conjugate gradient method or the minimal residual algorithm. With some modifications the resulting simulation can be improved but we also address remaining limitations. Some numerical examples show the capability of the final FEM software. In addition, we briefly discuss the special case of linear elasticity with small deformations. Here we directly derive a linear weak formulation with a saddle point structure and apply the adaptive mixed finite element method. It is shown that the presented findings can also be used to treat the nearly incompressible case.

No description available.