Cette thèse a pour objectif l"établissement de méthodes numériques pour l'optimisation multi-objectif de structures soumises à des facteurs incertains. Au cœur de ce travail, nous nous sommes focalisés sur l'adaptation du chaos polynomial pour l'évaluation non intrusive de la part de l'incertain. Pour atteindre l'objectif fixé, nous sommes confrontés à deux verrous : l'un concerne les coûts élevés de calcul d'une simulation unitaire par éléments finis, l'autre sa précision limitée. Afin de limiter la charge de calcul pour la construction du chaos polynomial, nous nous sommes concentrés sur la construction d'un chaos polynomial creux. Nous avons également développé un programme d’échantillonnage basé sur l’hypercube latin personnalisé prenant en compte la précision limitée de la simulation. Du point de vue de la modélisation nous avons proposé une approche multi-fidèle impliquant une hiérarchie de modèles allant des simulations par éléments finis complètes jusqu'aux surfaces de réponses en passant par la réduction de modèles basés sur la physique. Enfin, nous avons étudié l'optimisation multi-objectif de structures sous incertitudes. Nous avons étendu le modèle PCE des fonctions objectif à la prise en compte des variables déterministes de conception. Nous avons illustré notre travail sur des exemples d'emboutissage et sur la conception optimale des structures en treillis. / This thesis aims at Multi-Objective Optimization under Uncertainty in structural design. We investigate Polynomial Chaos Expansion (PCE) surrogates which require extensive training sets. We then face two issues: high computational costs of an individual Finite Element simulation and its limited precision. From numerical point of view and in order to limit the computational expense of the PCE construction we particularly focus on sparse PCE schemes. We also develop a custom Latin Hypercube Sampling scheme taking into account the finite precision of the simulation. From the modeling point of view,we propose a multifidelity approach involving a hierarchy of models ranging from full scale simulations through reduced order physics up to response surfaces. Finally, we investigate multiobjective optimization of structures under uncertainty. We extend the PCE model of design objectives by taking into account the design variables. We illustrate our work with examples in sheet metal forming and optimal design of truss structures.

De nombreuses techniques de résolution d'équations aux dérivées partielles sans maillage ont été développées dans la dernière décennie, proposant une alternative attrayante lorsque les éléments finis atteignent leurs limites. Notre travail se concentre sur l'étude de l'approximation diffuse, de ses applications au lissage et a la résolution des équations différentielles : les éléments diffus. Cependant, les solutions proposées s'appliquent aussi à d'autres méthodes et de nombreux résultats numériques illustrent chaque développement théorique. Dans un premier temps, nous étudions les techniques d'approximation sans maillage et les comparons à l'aide d'une méthode hybride. L'approximation myope est construite en modifiant le critère de construction des splines et montre ainsi la différence de nature entre moindres carres glissants et interpolateurs radiaux (krigeage, splines). Nous développons ensuite l'approximation diffuse Hermite dont l'application au calcul de la courbure des surfaces triangulées aboutit à une technique de haute résolution. Un algorithme de reverse engineering de modèle CAO s'appuie sur cette nouvelle technique et met en évidence son potentiel. L'intégration numérique est un ingrédient essentiel pour une méthode sans maillage. L'analyse du patch test nous conduit a la définition d'une technique robuste et assurant de bons rangs de convergence. Nous l'appliquons à l'approximation Hermite pour construire un modèle de poutre 3d pour le forage pétrolier. Un algorithme de contact roche-structure original est propose pour ce problème.

L’acier à double phase (DP) a été développé par l'industrie automobile pour le but de réduire le poids, l'amélioration de la performance de la sécurité et l'efficacité énergétique. Habituellement, l'acier DP contient des îlots de martensite dure noyée dans une matrice de ferrite doux. La synergie entre ces deux phases avec la microstructure inhomogène présente d'excellentes propriétés mécaniques. Les propriétés mécaniques (comportements de plasticité et de dégâts) d'acier DP sont principalement dérivés de sa microstructure, par exemple, la fraction de volume, la taille, la distribution et la morphologie de chaque phase constituante. Les approches micromécaniques sont largement appliquées pour prédire la plasticité et d'autres propriétés mécaniques de l'acier DP selon divers scénarios de chargement. Dans ce travail, la modélisation micromécanique de l'acier DP a été réalisée en utilisant des microstructures réelles ou artificielles. Une véritable microstructure est obtenue à partir de l'image métallographique, tandis qu'un générateur de microstructure artificielle à l'aide d'un algorithme d'affectation de phase améliorée basée sur l'optimisation de la topologie matériau est proposé d'étudier les propriétés mécaniques. Dans ce générateur artificiel, un processus d'affectation de phase est réalisé sur une mosaïque de Voronoï modifié pour obtenir une mesure représentative de l'élément de volume (VER) avec une bonne convergence. La méthode proposée comprend également une réduction appropriée décomposition orthogonale (POD) des courbes de débit (instantanés), qui sont calculés en utilisant le schéma asymptotique homogénéisation d'extension (AEH), pour identifier le contrôle des paramètres optimaux pour l'acier DP. Cette méthode numérique est vérifiée en utilisant DP590 et DP980 aciers qui indiquent un bon accord avec la contrainte d'écoulement à partir de mesures et prédiction de RVE basés sur de vraies microstructures. Les prédictions des modèles de déformation plastique, y compris des bandes de cisaillement en utilisant la microstructure artificielle ressemblent étroitement le comportement mécanique réel dans des conditions de chargement similaires. En outre, une interpolation a été adoptée pour obtenir une corrélation entre ces paramètres de contrôles basés sur l'identification des différents aciers DP. En outre, un modèle de substitution bi-niveau réduit est élaboré et présenté pour identifier les paramètres matériels du critère de rupture de Mohr-Coulomb (MMC). En utilisant cette méthode, le processus d'identification devient possible avec un nombre limité de tests Expérimentaux. La méthode combine des éléments critiques locaux associés à des modèles globaux. Le modèle de substitution de la souche de fracture construit en utilisant l'approximation diffuse et les éléments locaux, réduit le coût de calcul pour la recherche des paramètres matériels. Des simulations de fracturation sont effectuées globales pour mettre à jour la déformation à la rupture de la cible et pour calculer le déplacement de l'apparition de la panne correspondante. Des résultats probants sont obtenus par application successive de la conception de l'expérience (DOE) et l'amélioration des algorithmes de transformation de l'espace de conception. Le protocole d'identification proposée est validé avec de l'acier DP590. Robustesse de la méthode est confirmée par des valeurs initiales différentes. Ces investigations numériques fournissent nouvelle direction pour les simulations multi-échelles de la plasticité et de dégâts des comportements d'acier DP. De plus, ils contribuent efficacement à combler le fossé entre la recherche scientifique et à l'application de l'ingénierie des matériaux hétérogènes. / Dual-phase (DP) steel has been developed by automotive industry for the purpose of weight reduction, improvement in safety performance and fuel efficiency. Usually, DP steel contains hard martensite islands embedded in a soft ferrite matrix. Synergy between these two phases with the inhomogeneous microstructure exhibits excellent mechanical properties. The mechanical properties (plasticity and damage behaviors) of DP steel are mostly derived from its microstructure, e.g., volume fraction, size, distribution and morphology of each constituent phase. Micromechanical approaches are vastly applied to predict plasticity and other mechanical properties of DP steel under various loading scenarios. In this work, micromechanical modelling of DP steel has been performed using real or artificial microstructures. A real microstructure is obtained from metallographic image, while an artificial microstructure generator with an enhanced phase assignment algorithm based on material topology optimization is proposed to investigate the mechanical properties. In this artificial generator, phase assignment process is performed on a modified Voronoï tessellation to achieve the tailored representative volume element (RVE) with a good convergence. The proposed method also includes a proper orthogonal decomposition (POD) reduction of flow curves (snapshots), which are computed using the asymptotic extension homogenization (AEH) scheme, to identify the optimal controlling parameters for DP steel. This numerical method is verified using DP590 and DP980 steels that indicate a good agreement with the flow stress from measurements and RVE prediction based on real microstructures. Predictions of plastic strain patterns including shear bands using the artificial microstructure closely resemble the actual mechanical behavior under similar loading conditions. Moreover, an interpolation has been adopted to obtain a correlation between these controlling parameters based on the identification for various DP steels. Additionally, a bi-level reduced surrogate model is developed and presented to identify the material parameters of the Mohr-Coulomb (MMC) fracture criterion. Using this method, the identification process becomes feasible with a limited number of experimental tests. The method combines local critical elements associated with global models. The surrogate model of fracture strain constructed using the diffuse approximation and the local elements, reduced the computational cost for searching material parameters. Global fracture simulations are performed to update the target fracture strain and to compute the corresponding failure onset displacement. Convincing results are obtained via successive application of design of experiment (DOE) and enhanced design space transformation algorithms. The proposed identification protocol is validated with DP590 steel. Robustness of the method is confirmed with different initial values. These numerical investigations provide new direction for multiscale simulations of the plasticity and damage behaviors of DP steel. Moreover, they efficiently contribute to bridge the gap between scientific research and engineering application of heterogeneous materials.