La thèse présente une méthodologie pour construire des éléments finis de type " solide-coque " avec intégration réduite en vue des applications à la simulation de la mise en forme tel que l'emboutissage des tôles où ces éléments finis doivent présenter de bonnes aptitudes à modéliser la flexion mais également les situations de laminage de la tôle. A partir des éléments volumiques à 8 nœuds et 3 degrés de liberté par nœud (les 3 composantes du déplacement), un neuvième nœud est rajouté au centre de l'élément. Ce neuvième nœud n'est pourvu que d'un seul degré de liberté, le déplacement le long de la direction de l'épaisseur. Cette direction privilégiée a un nombre de points d'intégration supérieur ou égal à 3 mais l'intégration est réduite au centre de l'élément diminuant très sensiblement les temps CPU par rapport à une intégration complète. Un soin particulier a été pris pour contrôler tous les modes à énergie nulle dus à l'intégration réduite. Ce nœud supplémentaire permet une distribution linéaire de la déformation normale. Avec les lois de comportement complètement 3D ces nouveaux éléments solide-coque donnent des résultats similaires en flexion à ceux obtenus avec des éléments coques et état plan de contrainte. Le neuvième nœud joue le rôle d'un paramètre supplémentaire pour l'interpolation quadratique du déplacement dans la direction de l'épaisseur. Ce degré de liberté a une signification physique et un effort équivalent à une pression normale peut être prescrit. Dans les situations de pression normale et dans le cas du contact, la contrainte normale obtenue est physique ce qui n'est pas le cas de nombreux éléments solide-coque de la littérature. Le pincement ou le laminage des tôles est correctement modélisé. Pour valider ces éléments, un module d'emboutissage en U avec passage et laminage de la bande de tôle sur des rouleaux a été construit au laboratoire. La comparaison entre les efforts d'emboutissage calculés et mesurés est très bonne ainsi que la géométrie des bandes de tôle obtenue après retour élastique.

La simulation des procédés de mise en forme des composites à renforts tissés de type RTM est un enjeu majeur pour les industries de pointe mettant en œuvre ce type de matériaux. Au cours de ces procédés, la préforme tissée est souvent soumise à des déformations importantes. La connaissance et la simulation du comportement mécanique de la préforme à l'échelle macroscopique et à l'échelle mésoscopique s'avère souvent nécessaire pour optimiser la phase de conception de pièces composites formées par de tels procédés. Une analyse du comportement mésoscopique des préformes tissées de composites est d'abord proposée. Une loi de comportement hyperélastique isotrope transverse est développée, permettant de décrire le comportement mécanique de chacun des modes de déformation de la mèche : élongation dans la direction des fibres, compaction et distorsion dans le plan d'isotropie de la mèche, cisaillement le long des fibres. Une méthodologie est proposée pour identifier les paramètres de cette loi de comportement à l'aide d'essais sur la mèche et sur le tissu, et une validation par comparaison avec des essais expérimentaux est présentée. Une analyse du comportement macroscopique des renforts interlocks est ensuite proposée : une loi de comportement hyperélastique orthotrope est développée et implémentée. Cette loi, extension de la loi de comportement pour la mèche, est également basée sur une description phénoménologique des modes de déformation de la préforme. Une méthode d'identification des paramètres de cette loi de comportement est mise en œuvre, utilisant des essais expérimentaux classiques dans le contexte des renforts tissés (tension uniaxiale, compression, bias extension test, flexion). Cette seconde loi de comportement est validée par comparaison avec des essais de flexion et d'emboutissage hémisphérique.

En région montagneuse, les infrastructures et les voies de communications sont soumises à de nombreux risques naturels dont les phénomènes d'origine gravitaire. Au-delà du danger pour les usagers, les conséquences des interruptions du trafic deviennent problématiques d'un point de vue économique et il devient indispensable de sécuriser les itinéraires. La mise en place d'écrans de filets pare-pierres est une des solutions possibles pour la protection contre les éboulements rocheux. Cette thèse porte sur l'étude des écrans souples ou filets métalliques de protection contre les chutes de blocs et plus précisément sur l'écran développé par l'entreprise GTS. Le filet constitutif de ces écrans se différencie par rapport aux systèmes conventionnels par un comportement orthotrope, dû à un maillage spécifique. Dans cette étude nous caractérisons le comportement de ces écrans de filets sous des chargements statiques et dynamiques de type impact par une approche couplant l'expérimentation et la modélisation numérique. L'étude procède pas à pas, les divers constituants sont évalués de façon quasi-statique, en laboratoire, et également in-situ pour reproduire les conditions réelles d'utilisation, en particulier l'aspect dynamique. Une attention particulière concernant les dissipateurs d'énergie, qui représentent l'élément centrale de ce type de structure, nous a permis de mettre au point un élément fusible robuste et fiable. Une campagne d'essais en grandeur réel sur les écrans de filets étudiés a permis de valider deux classes énergétiques (3000 kJ et 5000 kJ) selon les recommandations européennes. Les données recueillis lors des expérimentations ont permis de calibrer et valider différentes modélisations numériques de type éléments finis et éléments discrets. La pertinence de la modélisation a été évaluée au niveau des différentes échelles étudiées, échelle d'une maille, échelle d'une nappe, échelle du dissipateur d'énergie et échelle de la structure entière. Les performances et les limites des deux approches, MEF (méthode des éléments finis) et MED (méthodes des éléments discrets) ont été évaluées pour ce qui est de nos modélisations.

Ce travail de recherche porte sur le comportement des milieux poreux (triphasiques), plus particulièrement les sols non saturés sous sollicitations hydro-mécaniques. Un modèle constitutif élastoplastique couplé est développé. Ce modèle original est formulé selon les principes suivants: une loi constitutive est développée pour décrire le comportement de chaque phase (squelette solide, liquide, et gaz). Ensuite, des relations de couplage sont ajoutées entre chacune des phases. Pour le comportement du squelette solide, une loi élastoplastique non associée est adoptée, avec deux surfaces de charges, en cisaillement et en compression. La partie hydrique est décrite par une formulation qui permet de prendre en compte l'effet d'hystérésis. Ce modèle a été enrichi par une relation de couplage hydromécanique qui permet d'exprimer la pression d'entrée d'air en fonction de la porosité. Ensuite, le couplage complet se fait avec la contrainte effective de Bishop en utilisant une nouvelle définition du paramètre de succion χ grâce à laquelle, les différents phénomènes présents dans la réponse des milieux poreux sous différentes sollicitations peuvent être reproduits. Ce modèle est validé par une confrontation à des données expérimentales issues de la littérature sur différents types de sol (sable, limon,...). Le modèle est implanté dans le code aux éléments finis Cast3M. L'analyse de problèmes particuliers, tels que la mise en œuvre d'un cas test d'un sol d'assise soumis à un cycle pluvial, ainsi que l'étude de la stabilité d'une pente, permette de montrer la capacité du modèle à reproduire le comportement des milieux poreux non saturés.

Les instabilités fluides peuvent dégrader le confinement du plasma au sein des tokamaks. Étant données les échelles spatio-temporelles, on choisit les modèles fluides obtenus à partir de la dérivation des modèles cinétiques. On dérive plusieurs modèles hiérarchiques de la MagnétoHydroDynamique (MHD) et en particulier les modèles de la MHD réduite du Current Hole et de l'équilibre de Grad-Shafranov. Une des difficulté de l'ensemble de ces modèles est de respecter l'équation modélisant l'absence de monopôles magnétiques. Pour assurer cette condition en tout point du domaine, le champ magnétique est réécrit avec un potentiel vecteur. L'utilisation de potentiels fait apparaître des équations faisant intervenir des dérivées d'ordre supérieurs. La stratégie numérique développée est l'utilisation de la méthode des éléments finis avec des éléments C1 de Bell. Sur un maillage non structuré, ces éléments ont l'intérêt de présenter une base réduite définir exclusivement avec des variables aux noeuds du maillage. Les modèles de MHD réduite du Current Hole et de Grad-Shafranov ont été résolus avec ces éléments. La résolution du cas test de Grad-Shafranov avec les conditions de bords exactes a permis d'obtenir l'ordre optimale de 5. La résolution du système du Current Hole avec ces éléments, validée par l'obtention du paramètre η1/3, a permis l'observation de développement d'instabilités en dents de scies.

La technique de renforcement des sols par inclusions linéaires rigides, c'est-à-dire capables de reprendre des efforts de flexion et cisaillement, connaît depuis quelques décennies un essor important bien que la simulation et le calcul de tels ouvrages reste un problème difficile, en raison de la très forte hétérogénéité du sol renforcé combinée au grand nombre d'inclusions mis en jeu dans de tels ouvrages. Une modélisation multi phasique a récemment été développée permettant d'aborder le dimensionnement de tels ouvrages de géotechnique avec des temps de calcul considérablement réduits par rapport à ceux d'une modélisation numérique classique, où les inclusions de renforcement doivent être discrétisées individuellement. On s'intéresse dans le cadre de ce travail à la mise en oeuvre de ce modèle dans le cadre du comportement élastoplastique des différents constituants, lorsque les effets de flexion et de cisaillement des inclusions sont pris en compte. Après avoir développé un certain nombre de solutions analytiques de référence, un outil numérique basé sur la méthode des éléments finis est mis au point, puis appliqué à la simulation du comportement de différentes structures renforcées, telles que les radiers de fondation sur groupes de pieux, ou la stabilisation de pentes instables par l'incorporation de telles inclusions.

Ce travail, à l'interface entre les mathématiques appliquées et la physique, s'articule autour de la modélisation numérique des vésicules biologiques, un modéle pour les globules rouges du sang. Pour cela, le modéle de Canham et Helfrich est adopté pour décrire le comportement des vésicules. La modélisation numérique utilise la méthode des fonctions de niveau dans un cadre éléments finis. Un nouvel algorithme de résolution numérique combinant une technique de multiplicateurs de Lagrange avec une adaptation automatique de maillages garantit la conservation exacte des volumes et des surfaces. Cet algorithme permet donc de dépasser une limitation cruciale actuelle de la méthode des fonctions de niveau, à savoir les pertes de masse couramment observées dans ce type de problémes. De plus, les propriétés de convergence de la méthode des fonctions de niveau se trouvent ainsi grandement améliorées, comme l'indiquent de nombreux tests numériques. Ces tests comprennent notamment des problémes d'advection élémentaires, des mouvements par courbure moyenne ainsi que des mouvements par diffusion de surface. Concernant l'équilibre statique des vésicules, une condition générale d'équilibre d'Euler-Lagrange est obtenue à l'aide d'outils de dérivation de forme. En dynamique, le mouvement d'une vésicule sous l'action d'un écoulement de cisaillement est étudié dans le cadre des nombres de Reynolds élevés. L'effet du confinement est considéré, et les régimes classiques de chenille de char et de basculement sont retrouvés. Finalement, pour la premiére fois, l'effet des termes inertiels est étudié et on montre qu'au delà d'une valeur critique du nombre de Reynolds, la vésicule passe d'un mouvement de basculement à un mouvement de chenille de char.

Ce travail de thèse a pour objectif l'élaboration et la caractérisation de composites Alumine/Zircone obtenus par voies conventionnelles, et dédiés à un usage orthopédique. Ces composites présentent une biocompatibilité prouvée, d'excellentes propriétés mécaniques ainsi qu'une grande stabilité. Ils sont plus résistants, plus fiables que l'alumine ou la zircone seules et permettent d'envisager des composants de tailles et formes plus exigeantes mécaniquement. Actuellement ces composites semblent les plus adaptés pour la réalisation de prothèses orthopédiques mais peuvent encore être optimisés via la modification des microstructures. La première partie de ce travail a concerné l'étude de la réalisation industrielle de composants de grande taille à partir d'une poudre. Les différentes étapes de l'élaboration sont traitées : pressage des composants, frittage et usinage. Cette première partie est majoritairement consacrée à l'étude des gradients thermique dans une sphère lors du frittage. Nous montrons qu'il est possible de modéliser et de mesurer les gradients thermiques dans le matériau de manière très réaliste, ainsi que d'obtenir des ordres de grandeur des contraintes mécaniques. On pourra ainsi envisager de tester numériquement les cycles de frittage en fonction de la géométrie des pièces frittées. Nous étudions par ailleurs la possibilité de réaliser un usinage des composants après un traitement de préfrittage, qui permettrait de diminuer les coûts et simplifier l'élaboration de composants de grande taille. La seconde partie de ce mémoire a permis de montrer que différents types de microstructures, présentant des propriétés mécaniques différentes, peuvent être obtenues par simple mélange de poudre. Ceci est possible par l'utilisation d'un traitement thermique adapté, la variation du taux de zircone et grâce à l'ajout de dopants (Si, Ca et Mg) jouant sur la mobilité des joints de grains d'alumine. Lors de l'utilisation de Ca ou Mg, le taux de zircone et la température ont un effet prépondérant sur l'aspect des microstructures, permettant d'obtenir des micro/micro-composites (< 16vol% de zircone et >1500°C) et nano/nano-composites (25vol% de zircone et T < 1500°C). Seul l'ajout conjoint de Si et de Ca pour des échantillons contenant 2,5vol% de zircone permet de conduire à des micro/nano-composites avec une grande proportion de zircone intragranulaire. Les observations sur des composites avec un taux de zircone proche du taux de percolation (16vol%) permettent de mettre en évidence l'ensemble des types de renforcement observés dans la littérature, en fonction de la température et des dopants utilisés. La variété des microstructures obtenues permet de progresser dans le contrôle des microstructures des composites alumine-zircone, mais aussi d'envisager d'autres applications de ces composites en fonction des mécanismes de renforcement observés et de leurs propriétés mécaniques et structurales.

Cette thèse est consacrée au développement d'outils pour FreeFem++ destinés à faciliter la résolution des problèmes d'optimisation. Ce travail se compose de deux parties principales. La première consiste en la programmation, la validation et l'exploitation d'interfaces permettant l¿utilisation de routines d'optimisation directement dans le logiciel. La seconde comprend le développement de solutions pour le calcul automatisé des dérivées, toujours au sein de FreeFem++, en exploitant les paradigmes de la différentiation automatique. FreeFem++ est un environnement de développement intégré dédié à la résolution numérique d¿équations aux dérivées partielles en dimension 2 et 3. Son langage ergonomique permet à l'utilisateur d'exploiter aisément ses nombreux outils de création de maillages, de résolution de systèmes linéaires, ainsi que ses bibliothèques d'éléments finis, etc... Nous introduisons les nouvelles routines d'optimisation désormais accessibles depuis la bibliothèque de modules du logiciel. En particulier, le logiciel libre d'optimisation sous contraintes IPOPT, qui implémente une méthode de points intérieurs très robuste pour l¿optimisation en grande dimension. Nous appliquons avec succès ces algorithmes à une série de problèmes concrets parmi lesquels la résolution numérique de problèmes de sur- faces minimales, la simulation de condensats de Bose-Einstein, ou encore un problème de positionnement inverse en mécanique des fluides. Une version prototypique de FreeFem++ contenant les outils de différentiation automatique est présentée, après avoir exposé les principes fondamentaux de cette méthode de calcul de dérivées pour le calcul scientifique.

Lorsqu'une onde acoustique se propage dans un milieu granulaire, elle est susceptible de provoquer la mobilité des grains, aussi infime soit-elle. Inversement, la mobilité d'un grain dans une matrice fluide peut induire un champ acoustique et dans les deux cas, l'énergie acoustique peut être transférée à la fois au travers des pores et des contacts entre grains. Nous avons mis au point un modèle original permettant de considérer ces deux modes de transfert d'énergie pour simuler la propagation d'ondes acoustiques dans les milieux granulaires immergés. Dans le cas des milieux granulaires secs, l'inertie du fluide est telle que l'énergie transférée dans l'air peut être négligée et le milieu modélisé avec des algorithmes de type "dynamique moléculaire". Au contraire, dans le cas de milieux immergés, l'énergie portée par le fluide ne peut pas être négligée et nous montrons que la méthode des domaines fictifs basée sur les multiplicateurs de Lagrange distribués permet de coupler les équations de la dynamique et l'équation d'onde. Nous utilisons la méthode des éléments finis pour propager l'onde dans le fluide, les grains étant modélisés en 2D par des sphères rigides et incompressibles afin de satisfaire les hypothèses de l'algorithme de dynamique moléculaire. Les résultats du modèle sur des expériences numériques simples mais pour lesquelles existent des solutions analytiques de l'acoustique mettent en évidence la validité du nouveau modèle. Nous en donnons une illustration pour l'étude des interactions subies par un empilement réaliste de multiples grains mobiles soumis à un signal acoustique.