Dans ce travail, nous nous sommes intéressés au développement et à l'analyse numérique de méthodes numériques capables de résoudre efficacement les problèmes de Helmholtz à 2D, notamment en régime moyenne et haute fréquence. La méthode que nous proposons s'inscrit dans la lignée des méthodes de type Galerkin discontinues (DG). Dans chaque élément du maillage, la solution est approchée en utilisant une superposition d'ondes planes. La continuité de la solution aux interfaces est renforcée en utilisant des multiplicateurs de Lagrange. La méthodologie proposée est une procédure en deux étapes: nous résolvons d'abord des problèmes locaux bien posés et ensuite un système global issu de la condition de continuité imposée sur les interfaces. Les plus importantes propriétés de la méthode sont: (a) les problèmes locaux obtenus sont associés à des matrices Hermitiennes et définies positives et (b) le système global, à résoudre dans la deuxième étape, est associé à une matrice Hermitienne et semi-définie positive. Les résultats numériques obtenus montrent la supériorité de la méthode proposée par rapport aux méthodes de type élément fini standard, mais aussi par rapport à d'autres méthodes de type DG, comme par exemple celle développée par Farhat et al (2003).

Le but de cette thèse est de modéliser les vibrations linéaires des liquides (ballottement) dans des environnements de faible gravité ou de microgravité. L'étude de la dynamique des engins spatiaux contenant des liquides, fait partie des applications possibles. Du fait du faible niveau de gravité, les phénomènes de tension superficielle (capillarité) à la surface libre du liquide sont prépondérants. Leur prise en compte introduit, d'une part, des forces supplémentaires proportionnelles à la courbure de la surface libre, et d'autre part une condition sur l'angle de contact entre la surface libre du liquide et la paroi interne du réservoir. La première partie de ce document traite le problème de recherche de la position d'équilibre du liquide au sein du réservoir, c'est-à-dire vise à déterminer la géométrie du ménisque formé par la surface libre du liquide à l'équilibre. La deuxième partie est consacrée à l'étude du ballottement du liquide autour de cette position d'équilibre. L'originalité de ce travail est de traiter ces deux problèmes (statique et dynamique) d'un point de vu tridimensionnel, en utilisant des formulations variationnelles évitant le calcul explicite des courbures de la surface libre. Dans les deux parties, la discrétisation a pu ainsi être faite en utilisant des Eléments Finis linéaires, et on a pu, dans la partie dynamique, aboutir à un problème spectral symétrique. Les résultats obtenus, en statique et en dynamique, ont été validés sur des cas simples où des solutions analytiques ou semi-analytiques sont disponibles.

Cette thèse porte sur la modélisation par éléments finis des procédés de mises en forme à froid, que sont le tréfilage et le laminage. Tout d'abord le comportement mécanique des aciers haut carbone a été mesuré grâce à une large campagne d'essais expérimentaux tout au long de cette gamme de mise en forme et une progressive anisotropie mécanique a été observée au cours du tréfilage. Puis, la simulation numérique du tréfilage et du laminage a été réalisée à l'aide de FORGE2005®. Le résultat principal concerne la prédiction de l'élargissement en fin de laminage qui est très imprécise avec une loi isotrope (sous estimation de la largeur de 10%). Cette sous-estimation passe à 5% avec une loi de comportement anisotrope. Ensuite, une troisième partie a porté sur l'étude microstructurale couplée à une analyse des mécanismes d'endommagement des aciers perlitiques au cours du tréfilage et du laminage. L'anisotropie mécanique provient de l'alignement des colonies de perlite au tréfilage et par l'apparition d'une orientation cristallographique préférentielle. Trois mécanismes d'endommagement ont pu être identifiés au cours du tréfilage. Lors du passage au laminage, les cinétiques de propagation de l'endommagement sont modifiées. La simulation a permis d'apporter des informations supplémentaires et de valider les observations expérimentales. Enfin, des calculs d'optimisation du tréfilage ont été effectués et ont permis d'étudier la sensibilité des fonctions objectifs (endommagement et force de tréfilage) aux paramètres d'optimisation (géométrie de filière). De plus, cette étude a mis en évidence que les solutions optimales diffèrent en fonction du choix de la fonction coût et qu'il est possible de diminuer l'endommagement sans trop augmenter la force de tréfilage et le risque de rupture.

Le travail présenté dans ce mémoire est consacré à l'application des méthodes de changement d'échelles pour la modélisation et le calcul des ouvrages construits en milieux rocheux fracturés. Après avoir présenté les deux approches usuellement adoptées, l'approche discontinue et l'approche continue, et leurs limites en termes de pertinence, dans une première partie essentiellement bibliographique (chap.1), nous avons formulé à partir des hypothèses sur les caractéristiques géométriques et mécaniques des deux constituants (matrice et joints), le comportement macroscopique en élasto plasticité du milieu homogénéisé aussi bien dans le cas périodique que dans le cas des milieux aléatoirement hétérogènes (chap.2). Les deux traits importants de ce comportement sont l'absence d'écruissage et l'absence de terme rendant compte de l'interaction entre les joints se coupant. Afin de valider le code d'éléments finis en déplacement que nous avons mis au point (chap.4), nous avons construit une solution analytique (convergence du tunnel circulaire) pour servir de référence (chap.3). Pour terminer, nous avons procédé à l'application du code numérique construit sur deux problèmes typiques de la Mécanique des Roches (chap.5).

Une étude sur la durabilité d'assemblages collés en milieu chaud et humide est proposée dans ce manuscrit. L'objectif a été de construire un outil numérique capable de prédire le comportement mécanique après vieillissement des assemblages collés en vue de leurs dimensionnements. A cette fin, une importante étude expérimentale a été menée sur un adhésif structural de type époxy chargé élastomère et des assemblages collés simple recouvrement acier/époxy/acier, destinés à des applications du secteur automobile. Les résultats de ces essais ont permis le développement d'un modèle de diffusion non linéaire physiquement motivé, capable de prédire la pénétration de l'eau dans l'adhésif pour une température donnée et les contraintes qui en résultent. Ces dernières ont été obtenues grâce à une loi de comportement thermohydromécanique couplée, développée spécialement pour l'adhésif structural. Toutes ces lois ont été intégrées dans un modèle éléments finis capable de simuler l'essai de traction/cisaillement sur l'assemblage avec ou sans vieillissement. L'impact du vieillissement sur la tenue de l'adhésif et/ou de l'interface dans l'assemblage a été discuté pour améliorer la compréhension des observations expérimentales.

La gamme des alliages à mémoire de forme offre une variété de comportements mécaniques : effet mémoire de forme, super élasticité, amortissement... Ces matériaux peuvent être exploités en vue d'améliorer les solutions techniques existantes, mais ils sont également à l'origine de nouvelles solutions et de produits innovants. Ces travaux de thèse s'attachent au développement d'un outil numérique d'aide au dimensionnement d'applications en AMF. L'accent est mis sur le couplage thermique-mécanique et l'exploitation des AMF sous forme d'actionneur. Le modèle du comportement macroscopique des AMF développé est basé sur une étude de la thermodynamique de la transformation de phase martensitique et propose une description de l'état interne du matériau basé sur un choix original des variables internes : la fraction volumique de martensite et la déformation de transformation propre de la martensite. Le modèle est validé sur une campagne de caractérisation du comportement mécanique et thermique du matériau. Le point de vue macroscopique adopté autorise l'utilisation raisonnable du modèle avec la méthode des Éléments finis. Le travail d'intégration numérique du modèle conduit à un algorithme robuste et efficace adapté au calcul de structures.

Le contrôle actif nécessite un modèle numérique représentatif de la structure réelle dont on souhaite diminuer les vibrations. Les méthodes de recalage sont les plus efficaces pour obtenir ce modèle. Les plus répandues se basent sur la minimisation d'une fonction objectif construite à partir de la solution d'Équations aux Dérivées Partielles (EDP) paramétrées. Le coût d'évaluation de cette fonction peut vite exploser lorsque les modèles sont trop complexes ou trop nombreux, ce qui arrive quand on souhaite une grande famille de structures similaires ou une structure dont le comportement varie à cause d'un vieillissement ou de phénomènes aléatoires. On parle alors de recalage multimodèle. Afin de construire une approximation de la fonction coût, nous introduisons une nouvelle méthode de résolution des EDP paramétrées, la Méthode Éléments Finis sur Algèbre Polynomiale (MÉFAP). Elle présente l'avantage d'introduire les variabilités paramétriques dans le modèle numérique sans changer la base éléments finis. Ceci est réalisé grâce à un anneau de polynômes multivariable. Nous mettons en œuvre la MÉFAP afin d'obtenir une approximation de l'erreur en relation de comportement modifiée, qui est un estimateur de la qualité d'un modèle numérique vis-à-vis de résultats expérimentaux. Nous traitons des cas de recalage simple puis du recalage multimodèle. Les exemples présentés sont représentatifs d'un ensembles de cartes électroniques. Ils comprennent des cas 1D ou 2D, piézoélectriques ou purement mécaniques, des structures virtuelles ou réelles, des modèles déterministes ou stochastiques.

Dans l'industrie pétrolière, de nombreuses campagnes sismiques visant à trouver de nouveaux champs pétroliers sont effectuées en mer. Afin de reproduire numériquement les résultats de ces campagnes d'acquisitions sismiques marine, les méthodes de différences finies ou d'éléments finis sont le plus souvent utilisées. Dans cette thèse, nous travaillons avec la méthode des éléments spectraux (Komatitsch et Tromp, 2002), qui est une méthode d'éléments finis explicite à haut degré d'interpolation polynomiale. Pour de l'acquisition sismique marine, la propagation d'ondes s'effectue dans des modèles composés de deux parties : une partie fluide (la couche d'eau épaisse homogène représentant l'océan) et une partie solide (les roches composant le fond des océans). La partie fluide étant moins intéressante que la partie solide où peuvent se trouver les réservoirs d'hydrocarbures, nous voulons une méthode numérique où la propagation d'ondes dans le fluide est la moins coûteuse possible en temps calcul. Un moyen de gagner ce temps est d'utiliser une méthode intégrant des pas de temps locaux (Diaz et Joly, 2005) où la partie fluide est traitée à l'aide d'un pas de temps plus grand que dans la partie solide. Pour cela, il faut s'assurer que l'augmentation du pas de temps sans modifier le maillage ne change pas la condition de stabilité CFL dans le fluide. Cette condition est calculée à partir de la vitesse maximale de propagation des ondes dans le domaine, et des pas de temps et d'espace. Or, les modèles utilisés présentent généralement un fort contraste des vitesses à l'interface fluide-solide en faveur du solide. C'est donc le pas de temps de ce dernier qui gouverne la condition CFL. Donc augmenter le pas de temps seulement dans le fluide ne changera pas la condition CFL et la simulation numérique pourra s'effectuer sans risquer l'explosion de la solution pour cause d'instabilité du schéma numérique. Mon travail a consisté à intégrer un tel procédé dans la méthode des éléments spectraux du code SPECFEM à 2D. Comme on a des pas de temps différents de part et d'autre de l'interface fluide-solide, les données à échanger entre les deux milieux ne sont pas coordonnées car à des temps différents. Pour palier à cette difficulté, on impose la conservation de l'énergie à notre système pour la construction des conditions de raccord. Cette manoeuvre implique néanmoins la résolution d'un système linéaire l'interface mais son coût est négligeable comparé aux calculs économisés dans la partie fluide. Afin de se rapprocher plus encore de la réalité des modèles pétroliers, j'ai travaillé sur l'implémentation de couches absorbantes de type CPML (convolutional perfectly match layer, voir Komatitsch et Martin 2007 ) qui permettent de simuler un milieu infini en absorbant toutes les ondes y pénétrant et donc évitent des réflexions parasites des bords du domaine. Ces conditions ont été implémentées pour des modèles fluide-solide à 2D.

De nos jours, l'intégration entre les différentes étapes de modélisation d'un produit mécanique est une nécessité afin de réduire le temps de conception et d'optimisation, de valider rapidement de nouveaux concepts et ainsi de gagner en coût, en délai et en qualité. Cela se traduit par une meilleure communication entre les différents outils de CAO et d'analyse. Cette tendance nécessite l'intégration des données entre les différentes tâches et l'unicité du modèle numérique. Cette intégration doit avoir le sens large du terme, c'est-à-dire de la conception vers le calcul et inversement. Notre activité de recherche tente à améliorer cette intégration à travers le développement d'un outil de reconstruction du modèle CAO à partir des résultats de simulation (maillage déformé).

Les travaux de cette thèse concernent le rayonnement acoustique d'une source périodique en temps placée dans un conduit Infini, contenant un guide en écoulement parallèle cisaille. Le phénomène est modélise a l'aide de l'équation de Galbrun, dont l'inconnue u est la perturbation de déplacement. L'objectif de cette étude est de développer une méthode éléments Nis, susceptible d'être étendue à des géométries et des écoulements plus complexes. Cette thèse fait suite a celle de Guillaume Legendre qui a établi, dans le cas d'un ecoulement uniforme, une formulation dite régularisée de l'´equation de Galbrun afin de corriger un défaut d'ellipticité. Le but de ce manuscrit est détendre cette méthode à un écoulement non uniforme. La difficulte supplémentaire vient du fait que la vorticite ψ = rot u (qui intervient dans le terme de régularisation) ne peut plus être calculée a priori car le cisaillement induit un couplage entre acoustique et hydrodynamique. En régime dissipatif, nous avons explicite ψ en fonction de u à l'aide d'une convolution (le long des lignes de courant). Si l'ecoulement est lent, cette formule de convolution (qui devient une intégrale très oscillante) peut être approchée par une formule differentielle beaucoup plus simple dont l'utilisation conduit a un modèle ”faible Mach”. Des idées similaires ont ensuite été utilisées pour résoudre le problème non dissipatif, a l'aide de couches PML. Les deux approches (exacte et ”faible Mach”) ont été validées par des tests numériques en 2D et en 3D.