La compréhension des phénomènes thermiques et mécaniques mis en jeu lors de la mise en forme des matériaux est généralement réalisée avec l’aide de simulations numériques. Ces simulations montrent leurs limites pour les procédés qui conduisent à de très grandes déformations de la matière. Dans ce cas, de très fortes distorsions du maillage se produisent pendant le calcul, entrainant une augmentation de l’erreur, voire l’arrêt prématuré de la simulation. Cette étude porte sur le développement d’une stratégie de remaillage adaptative afin d’éviter les distorsions des éléments pendant les simulations en grandes transformations. La méthode proposée a été intégrée dans un environnement de calcul utilisant le solveur ABAQUS/Explicit, un mailleur 3D et un algorithme de transfert de champ.La méthode h-adaptative en combinaison avec un critère de contrôle basé sur l’endommagement et un estimateur d’erreur de type Zienkiewicz-Zhu Z2 (SPR-amélioré) ont été implantés. Le maillage initial est remplacé par un nouveau maillage avec le niveau de qualité désiré par l’utilisateur, tout en minimisant le nombre des degrés de liberté. Cette technique s’est montrée robuste et entièrement automatique pour déterminer la taille optimale des nouveaux éléments. Une fois le nouveau maillage généré, toutes les variables doivent être soigneusement transférées. Plusieurs techniques de transfert sont décrites et comparées. Des améliorations permettant d’augmenter leurs efficacités en termes de diffusion de l’information et de stabilité numérique ont été proposées. Une attention particulière est portée à la restauration de l'équilibre mécanique local du système. Les différentes techniques développées ont permis de modéliser différents procédés entrainant de grandes déformations élastoplastiques avec endommagement. Dans toutes les applications testées, il a été montré une amélioration de la précision et de la qualité des résultats numériques obtenus. Pour des opérations d’usinage, des mesures de champs cinématiques à travers la technique de corrélation d’images ont été réalisées afin de déterminer les champs de déformation en pointe d’outil. Ces mesures ont servi à la validation de la simulation numérique à l’échelle locale. La comparaison des champs cinématiques expérimentaux avec ceux issus du calcul éléments finis met en évidence la robustesse du processus d’adaptation du maillage proposée pour retranscrire les phénomènes locaux observés expérimentalement. En effet, la reproduction de l’écoulement de la matière sur les bords et la géométrie du copeau sont en très bonne corrélation avec les résultats expérimentaux. Ce développement a permis de proposer une description nouvelle du processus de formation des bandes de cisaillement. / In this work, a fully automated adaptive remeshing strategy, based on a tetrahedral element to simulate various 3D metal forming processes, was proposed. The aim of this work is to solve problems associated with the severe mesh distortion that occurs during the computation and which may be incompatible with the evolution of the physical behavior of the FE solution. Indeed, the quality of the mesh conditions affects the accuracy of the calculations. The proposed strategy is integrated in a computational platform which integrates a finite element solver (Abaqus/Explicit), 3D mesh generation and a field transfer algorithm.The base idea is to use the h-adaptive methodology in the combination with a damage-criterion error and Zienkiewicz-Zhu Z2 type error estimator (SPR-improved) to locally control the mesh modification-as-needed. Once a new mesh is generated, all history-dependent variables need to be carefully transferred between subsequent meshes. Therefore, different transfer techniques are described and compared. An important part of this work concerns the presentation of the proposed modification of the field transfer operator and a special attention is given to restore the local mechanical equilibrium of the system. During the large elasto-plastic deformation simulation with damage, the necessary steps for remeshing the mechanical structure are presented. The several types of applications are also given. For all studied applications, the above strategy can improve the accuracy and quality of numerical results. It also has benefits to decide how refined a mesh needs to be to reach a particular level of accuracy, or how coarse the mesh can be without unacceptably impacting solution accuracy.For the machining processes, kinematic field measurements using Digital image Correlation were performed to validate the numerical simulation at the local level. The comparison of the experimental kinematic fields and those resulting from the FE calculation highlights the robustness of the proposed mesh adaptation process which can transcribe the experimental local phenomena. Also, the reproduction of the material flow at the edges and the chip are correlated with the experimental results accurately. Finally, the physical study of the numerical results can be allowed to propose an innovative description of ASB formation.

les travaux développés dans cette thèse portent sur l'optimisation du compromis débit-distorsion pour des maillages triangulaires quantifiés par un codeur géométrique. De nombreux travaux considèrent que la géométrie est la composante la plus importante d'un maillage polygonal. Les codeurs issus de ces travaux appelés codeurs géométriques, incluent une étape de remaillage et une transformée en ondelettes qui permet une analyse multirésolution efficace. L'objectif de nos travaux est de proposer un codeur géométrique incluant une allocation binaire qui optimise la qualité visuelle du maillage reconstruit en fonction du débit. Le but de l'allocation est de minimiser la distance surface-surface entre le maillage d'entrée et le maillage quantifié en fonction d'un débit de consigne. Cette distance qui traduit la différence géométrique entre deux maillages demande un lourd processus d'un point de vue calculatoire. Utiliser une approximation de cette mesure comme critère de distorsion est donc préférable si l'on souhaite une allocation rapide. Nous avons alors montré que sous certaines hypothèses cette mesure pouvait être estimée par une somme pondérée des erreurs de quantification des coefficients d'ondelettes. De plus, les particularités statistiques des coefficients d'ondelettes<br />géométriques permettent l'utilisation de modèles théoriques pour le<br />débit et la distorsion des sous-bandes de coefficients. Finalement,<br />nous proposons un codeur géométrique incluant une allocation rapide et performante qui optimise la quantification des coefficients pour que la qualité visuelle de l'objet reconstruit soit maximisée sous la contrainte d'un débit total fixé. Expérimentalement, l'algorithme<br />proposé donne de meilleurs résultats que toutes les méthodes de l'état de l'art.

Ce travail de thèse concerne l'analyse structurelle des maillages triangulaires surfaciques, ainsi que leur traitement en vue de l'amélioration de leur qualité (remaillage) ou de leur simplification. Dans la littérature, le repositionnement des sommets d'un maillage est soit traité de manière locale, soit de manière globale mais sans un contrôle local de l'erreur géométrique introduite, i.e. les solutions actuelles ne sont pas globales ou introduisent de l'erreur géométrique non-contrôlée. Les techniques d'approximation de maillage les plus prometteuses se basent sur une décomposition en primitives géométriques simples (plans, cylindres, sphères etc.), mais elles n'arrivent généralement pas à trouver la décomposition optimale, celle qui optimise à la fois l'erreur géométrique de l'approximation par les primitives choisies, et le nombre et le type de ces primitives simples. Pour traiter les défauts des approches de remaillage existantes, nous proposons une méthode basée sur un modèle global, à savoir une modélisation graphique probabiliste, intégrant des contraintes souples basées sur la géométrie (l'erreur de l'approximation), la qualité du maillage et le nombre de sommets du maillage. De même, pour améliorer la décomposition en primitives simples, une modélisation graphique probabiliste a été choisie. Les modèles graphiques de cette thèse sont des champs aléatoires de Markov, ces derniers permettant de trouver une configuration optimale à l'aide de la minimisation globale d'une fonction objectif. Nous avons proposé trois contributions dans cette thèse autour des maillages triangulaires 2-variétés : (i) une méthode d'extraction statistiquement robuste des arêtes caractéristiques applicable aux objets mécaniques, (ii) un algorithme de segmentation en régions approximables par des primitives géométriques simples qui est robuste à la présence de données aberrantes et au bruit dans la position des sommets, (iii) et finalement un algorithme d'optimisation de maillages qui cherche le meilleur compromis entre l'amélioration de la qualité des triangles, la qualité de la valence des sommets, le nombre de sommets et la fidélité géométrique à la surface initiale.

L'objectif de cette thèse était de proposer, puis de valider expérimentalement, un ensemble de méthodes algorithmiques permettant le remaillage parallèle de maillages distribués, en s'appuyant sur une méthode séquentielle de remaillage préexistante. Cet objectif a été atteint par étapes : définition de structures de données et de schémas de communication adaptés aux maillages distribués, permettant le déplacement à moindre coût des interfaces entre sous-domaines sur les processeurs d'une architecture à mémoire distribuée ; utilisation d'algorithmes de répartition dynamique de la charge adaptés aux techniques parallèles de remaillage ; conception d'algorithmes parallèles permettant de scinder le problème global de remaillage parallèle en plusieurs sous-tâches séquentielles, susceptibles de s'exécuter concurremment sur les processeurs de la machine parallèle. Ces contributions ont été mises en oeuvre au sein de la bibliothèque parallèle PaMPA, en s'appuyant sur les briques logicielles MMG3D (remaillage séquentiel de maillages tétraédriques) et PT-Scotch (repartitionnement parallèle de graphes). La bibliothèque PaMPA offre ainsi les fonctionnalités suivantes : communication transparente entre processeurs voisins des valeurs portées par les noeuds, les éléments, etc. ;remaillage, selon des critères fournis par l'utilisateur, de portions du maillage distribué, en offrant une qualité constante, que les éléments à remailler soient portés par un unique processeur ou bien répartis sur plusieurs d'entre eux ; répartition et redistribution de la charge des maillages pour préserver l'efficacité des simulations après remaillage.

Cette thèse est consacrée à l'optimisation de la topologie et de la forme de systèmesmultiphysiques motivés par des applications de l'industrie aéronautique. Nouscalculons les dérivées de forme de fonctions de coût arbitraires pour un modèlefluide, thermique et mécanique faiblement couplé. Nous développons ensuite unalgorithme de type gradient adapté à la résolution de problèmes d'optimisation deformes sous contraintes qui ne requiert par de réglage de paramètres nonphysiques. Nous introduisons ensuite une méthode variationnelle qui permet decalculer des intégrales le long de rayons sur un maillage par la résolution d'unproblème variationnel qui ne requiert pas la détermination explicite de ces lignessur la discrétisation spatiale. Cette méthode nous a ainsi permis d'imposer unecontrainte de non-mélange de phases pour une application à l'optimisationd'échangeurs de chaleur bi-tubes. Tous ces ingrédients ont été employés pour traiterune variété de cas tests d'optimisation de formes pour des systèmes multi-physiques2-d ou 3-d. Nous avons considéré des problèmes à une seule, deux ou bien troisphysiques couplées en 2-d, et des problèmes de tailles relativement élevées en 3-dpour la mécanique, la conduction thermique, l'optimisation de profils aérodynamiques,et de la forme de systèmes en interaction fluide-structure. Un dernier chapitred'ouverture est consacré à l'étude de modèles homogénéisées d'ordres élevés pour lessystèmes elliptiques perforés. Ces équations d'ordres élevés englobent les troisrégimes homogénéisés classiques associés à divers rapports d'échelles pour la tailledes obstacles. Elles pourraient permettre, dans de futurs travaux, de développer denouvelles méthodes d'optimisation pour les systèmes fluides caractérisés par desmotifs multi-échelles, ainsi que couramment rencontré dans la conception deséchangeurs thermiques industriels. / This work is devoted to shape and topology optimization of multiphysics systemsmotivated by aeronautic industrial applications. Shape derivatives of arbitraryobjective functionals are computed for a weakly coupled thermal fluid-structuremodel. A novel gradient flow type algorithm is then developed for solving genericconstrained shape optimization problems without the need for tuning non-physicalmetaparameters. Motivated by the need for enforcing non-mixing constraints in thedesign of liquid-liquid heat exchangers, a variational method is developed in orderto simplify the numerical evaluation of geometric constraints: it allows to computeline integrals on a mesh by solving a variational problem without requiring theexplicit knowledge of these lines on the spatial discretization. All theseingredients allowed us to implement a variety of 2-d and 3-d multiphysics shapeoptimization test cases: from single, double or three physics problems in 2-d, tomoderately large-scale 3-d test cases for structural design, thermal conduction,aerodynamic design and a fluid-structure interacting system. A final opening chapterderives high order homogenized equations for perforated elliptic systems. These highorder equations encompass the three classical regimes of homogenized modelsassociated with different obstacle's size scalings. They could allow, in futureworks, to develop new topology optimization methods for fluid systems characterizedby multi-scale patterns as commonly encountered in industrial heat exchanger designs.

Ce travail de thèse concerne l'analyse structurelle des maillages triangulaires surfaciques, ainsi que leur traitement en vue de l'amélioration de leur qualité (remaillage) ou de leur simplification. Dans la littérature, le repositionnement des sommets d'un maillage est soit traité de manière locale, soit de manière globale mais sans un contrôle local de l'erreur géométrique introduite, i.e. les solutions actuelles ne sont pas globales ou introduisent de l'erreur géométrique non-contrôlée. Les techniques d'approximation de maillage les plus prometteuses se basent sur une décomposition en primitives géométriques simples (plans, cylindres, sphères etc.), mais elles n'arrivent généralement pas à trouver la décomposition optimale, celle qui optimise à la fois l'erreur géométrique de l'approximation par les primitives choisies, et le nombre et le type de ces primitives simples. Pour traiter les défauts des approches de remaillage existantes, nous proposons une méthode basée sur un modèle global, à savoir une modélisation graphique probabiliste, intégrant des contraintes souples basées sur la géométrie (l'erreur de l'approximation), la qualité du maillage et le nombre de sommets du maillage. De même, pour améliorer la décomposition en primitives simples, une modélisation graphique probabiliste a été choisie. Les modèles graphiques de cette thèse sont des champs aléatoires de Markov, ces derniers permettant de trouver une configuration optimale à l'aide de la minimisation globale d'une fonction objectif. Nous avons proposé trois contributions dans cette thèse autour des maillages triangulaires 2-variétés : (i) une méthode d'extraction statistiquement robuste des arêtes caractéristiques applicable aux objets mécaniques, (ii) un algorithme de segmentation en régions approximables par des primitives géométriques simples qui est robuste à la présence de données aberrantes et au bruit dans la position des sommets, (iii) et finalement un algorithme d'optimisation de maillages qui cherche le meilleur compromis entre l'amélioration de la qualité des triangles, la qualité de la valence des sommets, le nombre de sommets et la fidélité géométrique à la surface initiale. / The work in this thesis concerns structural analysis of 2-manifold triangular meshes, and their processing towards quality enhancement (remeshing) or simplification. In existing work, the repositioning of mesh vertices necessary for remeshing is either done locally or globally, but in the latter case without local control on the introduced geometrical error. Therefore, current results are either not globally optimal or introduce unwanted geometrical error. Other promising remeshing and approximation techniques are based on a decomposition into simple geometrical primitives (planes, cylinders, spheres etc.), but they generally fail to find the best decomposition, i.e. the one which jointly optimizes the residual geometrical error as well as the number and type of selected simple primitives. To tackle the weaknesses of existing remeshing approaches, we propose a method based on a global model, namely a probabilistic graphical model integrating soft constraints based on geometry (approximation error), mesh quality and the number of mesh vertices. In the same manner, for segmentation purposes and in order to improve algorithms delivering decompositions into simple primitives, a probabilistic graphical modeling has been chosen. The graphical models used in this work are Markov Random Fields, which allow to find an optimal configuration by a global minimization of an objective function. We have proposed three contributions in this thesis about 2-manifold triangular meshes : (i) a statistically robust method for feature edge extraction for mechanical objects, (ii) an algorithm for the segmentation into regions which are approximated by simple primitives, which is robust to outliers and to the presence of noise in the vertex positions, (iii) and lastly an algorithm for mesh optimization which jointly optimizes triangle quality, the quality of vertex valences, the number of vertices, as well as the geometrical fidelity to the initial surface.