Dans la première partie de ce travail nous considérons des problèmes hyperboliques du premier ordre linéaires où des problèmes paraboliques linéaires dégénérés en temps. En utilisant une méthode de matrice de masse singulière, nous proposons une méthode d’élément finis permettant d’avoir des estimations d’erreur en espace optimale pour l’élément fini de Lagrange P1 par exemple. Nous appliquons ces résultats au cas d’un système parabolique utilisé en electroradiologie. La seconde partie est consacrée aux polynômes Lp extrémaux à l’extérieur du cercle unité associés à une mesure de la forme générale α = βa + βs + γ, où βa est régulière, βs singulière et γ discrète. Dans un premier temps nous considérons βs = 0, et nous avons généralisé au cas Lp des résultats connus dans le cas L2. Dans le cas où βs = 0 nous montrons les mêmes résultats (formules d’optimalité) mais en utilisant d’autres fonctions de régularité. / In the first part of this work, we deal with, linear hyperbolic problems of first order or linear parabolic problems, which are degenerated with respect to the time operator. By using a singular mass matrix technique, we propose a finite element method allowing to get optimal error estimates with respect to space for the Lagrange first order finite element for example. Then our method is applied to a parabolic system degenerated with respect to time which is used in electrocardiology. The second part of this work is dedicated to extremal polynomials in Lp , outside to the unit circle associated to a measure α, with a general form given by α = βa + βs + γ. The regular part is denoted βa , the singular part βs and the discrete part γ. In a first step we take βs = 0, and we generalized to the Lp case the known results in the L2 case. When the singular part is non zero, by using different regularity functions, we get the same optimality formulae.

Les équations de Maxwell, qui gouvernent le comportement des champs électromagnétiques mésoscopiques, ne forment pas un ensemble clos d'équations. Il faut donc les compléter avec des relations constitutives, qui caractérisent un matériau. Nous avons voulu donner des relations fort générales, et les simplifier ensuite, en explicitant toutes les restrictions physiques que cela comportait, jusqu'à obtenir un système fini d'équations aux dérivées partielles non-linéaires couplées résoluble numériquement (méthode des éléments finis). Cela nous a permis d'étudier la diffraction d'un champ électromagnétique par un système, ayant la forme d'un cylindre ou d'un cristal photonique, fini ou infini, générant une deuxième et une troisième harmonique. Ces différentes composantes du champ électrique interagissent entre elles, par l'effet Kerr-optique ou la déplétion de l'onde de pompe. Nous avons aussi généralisé à l'optique non-linéaire quelques propriétés des matériaux bien connues dans le régime linéaire, comme l'anisotropie ou la perte d'énergie électrique.

Le travail présenté est une analyse critique des éléments finis poroélastiques basés sur la théorie de Biot, permettant la simulation du comportement vibroacoustique des matériaux poreux à squelette déformable (mousses polymères, en particulier). Dans un premier temps, une étude de convergence montre que les critères de maillage, habituellement utilisés pour les éléments monophasiques, s'appliquent eu égard aux longueurs d'onde de Biot, mais restent insuffisants pour des applications tridimensionnelles. Une deuxième partie est consacrée à la validation expérimentale. Une attention particulière est portée sur la caractérisation expérimentale des matériaux. La première validation concerne la mesure d'impédance en conduit d'un échantillon résonnant. Différentes conditions aux limites sont réalisées afin de tester la validité des lois de comportement isotrope ou isotrope transverse pour le squelette. La deuxième validation porte sur la vibration d'une plaque amortie par collage d'une couche poreuse. Bien que les tendances soient bonnes, la difficulté à caractériser le comportement mécanique du squelette est mise en relief. Dans la dernière partie, on s'intéresse à l'analyse du comportement des matériaux poroélastiques. Un calcul de la répartition des puissances dissipées et réactives est développé à partir de la formulation par éléments finis. Appliqué au cas d'une couche poreuse collée à une plaque, il permet de montrer l'importance de la dissipation due à la viscoélasticité du squelette. Un modèle de plaque équivalente, considérant la couche poreuse comme un milieu viscoélastique, sujet au cisaillement, a été développé. Les bons résultats obtenus permettent un allégement conséquent des calculs en ce qui concerne le comportement vibratoire de cette structure.

En ce début de 21e siècle, le bruit des transports est l'une des principales nuisances. Un important travail est réalisé par les acousticiens afin d'améliorer l'efficacité des écrans acoustiques.<br />Actuellement, le choix des protections anti-bruit se porte souvent sur des solutions triviales ou issues d'études paramétriques. Or le CSTB possède des codes de calculs numériques avancés dédiés à la simulation de la propagation acoustique en milieu extérieur. Le but de cette thèse est de coupler ces codes prévisionnels avec un outil mathématique à développer permettant une recherche systématique des solutions optimales.<br />Une étude bibliographique a permis de dégager des algorithmes applicables au cas de l'optimisation de paramètres caractérisant un écran acoustique. Ils ont été appliqués avec succès aux écrans multi-diffracteurs. Des méthodes d'optimisation multiobjectif ont en outre permis d'optimiser simultanément le prix de revient de la protection et son efficacité acoustique. Dans une optique de développement durable, un travail a été réalisé sur l'efficacité de protections acoustiques de faible hauteur en milieu urbain dans le but de créer des espaces calmes. Des mesures sur modèles réduits ont été confrontées aux simulations acoustiques et l'optimisation de la forme de ces protections a permis de dégager des solutions optimales. Enfin, les méthodes d'optimisation ont été appliquées au cas de l'optimisation de trajectoires.<br />Les principaux objectifs de la thèse ont été atteints, notamment par la création d'un outil utilisant des algorithmes d'optimisation (Nelder Mead et Algorithmes Génétiques) ainsi que les codes de calculs de simulation acoustique en milieu extérieur.

Alternatives aux méthodes de Monte-Carlo pour le traitement des problèmes de propagation d'incer- titudes dans les modèles mécaniques structuraux, les méthodes d'éléments finis stochastiques (MEFS) connaissent un succès grandissant depuis une dizaine d'années, concrétisée par de nombreux travaux de recherche internationaux. Le présent travail est une contribution à ces recherches et son but est double. D'une part, considérant la MEFS spectrale (Ghanem & Spanos, 1991), qui est actuellement très utilisée, nous en faisons une présentation détaillée afin d'en cerner les limites. Cette méthode, essentiellement valable pour les problèmes linéaires, présente l'internet de permettre l'emploi, non seulement de variables aléatoires mais également de processus ou champs stochastiques pour la modélisation probabiliste des paramètres incertains du modèle. Deux applications, l'une sur une barre en traction, l'autre sur une structure formée de poutres modélisant un assemblage bois, permettent de juger de ses possibilités. D'autre part, et c'est le coeur du travail, nous proposons une MEFS originale pour la résolution de problèmes mécaniquement non linéaires. Cette approche, de type surface de réponse, comprend deux étapes clés : une projection de la réponse mécanique non linéaire sur une base de polynomes d'Hermite et une approximation de cette réponse par B-splines cubiques interpolantes pour le calcul des coeficients de la projection. Sa pertinence vis-à-vis de l'estimation des moments de la réponse est jugée à travers quatre exemples de complexité variable : un treillis élastoplastique, une sphère élastoplastique sous pression et un cylindre en contact avec un plan, étudiée en phases élastique puis élastoplastique.

Dans le cadre de la conception de turbo-machines aéronautiques, l'amélioration du rendement est assurée, entre autres, par la diminution du jeu fonctionnel entre les parties tournantes et statiques, telles que les roues aubagées et le carter qui les entoure. Cette stratégie a des conséquences sur le comportement vibratoire des turbo-machines en favorisant l'apparition de contacts structuraux entre les sommets d'aubes et le carter, par exemple. Lorsque plusieurs points de contact co-existent, des phénomènes d'interaction modale peuvent apparaître. La simulation de ces phénomènes, potentiellement dangereux, dans un cadre indutriel est au cœur de notre étude. Du fait de la dimension du problème à traiter, des méthodes de réduction modale sont utilisées : une méthode à interfaces fixes, la méthode de Craig-Bampton, et une méthode à interfaces libres, la méthode de Craig-Martinez. L'étude porte notamment sur la performance de ces méthodes de réduction lorsqu'elles sont associées à une non linéarité de type contact. Le phénomène d'interaction modale est tout d'abord étudié sur des modèles 2D simplifiés avec l'analyse de la sensibilité des régimes d'interaction à la taille des modèles réduits, puis sur les modèles 3D industriels pour lesquels la solution vraie n'est pas accessible. La gestion du contact dans le 3D fait appel aux B-splines bicubiques surfaciques afin d'assurer de meilleures propriétés à la surface de contact et de faciliter la résolution numérique. Les résultats obtenus à l'aide des deux méthodes de réduction modale sont comparés et permettent de déterminer la méthode la plus appropriée pour l'industrialisation du code présentée en annexe de ce mémoire.

Cette thèse est composée de trois parties. Dans la première, nous considérons un système d'équations Reynolds Averaged Navier-Stokes en 3D, modélisant le couplage de deux fluides turbulents ( par exemple, océan/atmosphére). Nous proposons un schéma numérique, et nous montrons sa convergence vers l'unique solution du modèle.<br />La seconde partie est consacrée à une extension des méthodes spectrales dans des géométries complexes. Cette nouvelle méthode s'appuie sur deux idées: traitement des conditions aux limites de Dirichlet par pénalisation, en suivant la méthode de Nitsche, et une approximation de la géométrie par des pavés, en utilisant une octree (par exemple). <br />Nous donnons des erreurs de projection polyômiale et des estimations a priori. <br />Enfin, la dernière partie est consacrée au calucl scientifique où on a implémenté en C++ et validé cette méthode dans le logiciel FreeFem3d.

Pour les problèmes linéaires, différentes démarches existent pour estimer la qualité de la solution éléments finis. Elles conduisent généralement à l'estimation d'une norme de l'erreur globale. Mais le choix d'une précision globale fondée sur une norme du déplacement est souvent délicat car il n'y a pas de lien direct quantitatif avec une erreur sur des quantités mécaniques locales (appelées quantités d'intérêt) ayant un intérêt pour l'ingénieur. L'objectif de ce travail est de développer et de mettre en oeuvre dans Code_Aster une méthode d'estimation d'erreur en quantité d'intérêt pour le calcul de structures industrielles.<br /><br />Une étude bibliographique a permis de situer l'état de l'art et de faire le point sur les méthodes d'estimation d'erreur globale et locale et sur les techniques d'adaptation de maillages afin de déterminer celles qui sont les plus adaptées pour le calcul de structures en milieu industriel.<br /><br />Ensuite, une comparaison a été faite afin de déterminer la pertinence de l'utilisation de l'erreur en quantité d'intérêt dans un processus adaptatif par rapport à l'utilisation de l'erreur en norme de l'énergie, d'un raffinement uniforme ou d'un raffinement local dans la zone d'intérêt. L'étude a montré que l'utilisation de l'erreur en quantités d'intérêt est la stratégie qui permet toujours d'obtenir la meilleure précision pour un nombre d'éléments donné : cette stratégie est la plus pertinente dans un contexte de calcul de structures industrielles.<br /><br />Compte tenu des résultats précédents, une méthode d'estimation d'erreur en quantité d'intérêt a été développée et implémentée dans Code_Aster. Elle est basée sur la relation fondamentale de l'erreur en quantité d'intérêt et sur un estimateur d'erreur globale en norme de l'énergie de type résidus explicites. Cet estimateur ne permettant pas d'accéder à des bornes de l'erreur mais à une valeur approchée, il s'apparente donc à un indicateur d'erreur. Pour réaliser cette estimation, la résolution d'un problème dual est nécessaire. Les chargements de ce problème, particulier pour chaque quantité d'intérêt considérée, sont fournis pour de nombreuses quantités d'intérêt (notamment pour la contrainte de Von Mises). Le calcul final de l'erreur en quantité d'intérêt, à partir du calcul primal et du calcul dual est automatique. Une quantité d'intérêt pour l'estimation de l'erreur sur les facteurs d'intensité des contraintes a également été développée en 2D et 3D. Des exemples numériques ont montré le bon comportement de cet indicateur pour guider un processus d'adaptation de maillage pour les différentes quantités d'intérêt.<br /><br />Enfin, la stratégie de contrôle d'erreur développée a été utilisée, avec Code_Aster, dans le cadre d'études de structures industrielles. Deux études effectuées au département AMA, à EDF R&D sont présentées. La première consiste à étudier la contrainte de Von Mises dans un rotor HP, la seconde s'intéresse à la contrainte verticale dans les goujons d'un assemblage boulonné.

Les travaux de recherche menés dans le cadre de cette thèse concernent l'étude du comportement des groupes et des réseaux de micropieux. L'interaction entre les micropieux d'un groupe, ainsi que l'effet d'inclinaison des micropieux au sein d'un réseau sont traités à travers les trois parties de cette thèse.<br />La première partie est consacrée à la définition, la classification et les domaines d'application des micropieux. Les différentes méthodes usuelles de calcul des micropieux (isolés, en groupes ou en réseaux) sont également présentées.<br />La deuxième partie présente la synthèse de divers essais expérimentaux réalisés, dans le cadre du projet national FOREVER, sur des groupes et/ou des réseaux de micropieux. Ces essais sont réalisés en vraie grandeur ou en modèle réduit (centrifugeuse, chambre d'étalonnage ou cuve expérimentale).<br />La troisième partie traite de deux approches numériques différentes pour l'étude du comportement des micropieux. La méthode des fonctions de transfert de charge est utilisée à travers le programme de calcul des groupes de pieux (GOUPEG) qui est développé pour tenir compte de l'effet de l'inclinaison des micropieux. Le cas d'un réseau élémentaire de 4 micropieux et d'un groupe équivalent est étudié. Une étude paramétrique sur l'effet de l'inclinaison des micropieux au sein d'un réseau élémentaire est présentée. Enfin, la méthode des éléments finis est utilisée (le logiciel CESAR-LCPC) pour l'étude du comportement d'un réseau élémentaire de 4 micropieux (en double chevalet). Cette dernière approche fournit des résultats intéressants sur le comportement du massif du sol sous les micropieux inclinés.

Cette étude, réalisée au Centre des Matériaux de l'Ecole des Mines de Paris, porte sur les problèmes d'identification et de transférabilité des lois de fissuration des aciers utilisés à haute température. Une approche globale, fondée sur les paramètres C* et J de la mécanique non linéaire de la rupture, a été utilisée pour caractériser l'amorçage et la propagation des fissures en<br />fluage. Les nuances d'aciers étudiées sont : les aciers ferritiques 1Cr-1Mo-1/4V (chaud et froid, travaillant à 540°C et 250°C) utilisés dans les centrales thermiques et l'acier inoxydable austénitique 316 L(N) utilisé dans les centrales nucléaires. Au cours de cette thèse, une base de<br />données a été mise en place, elle regroupe plusieurs essais de fatigue, de fluage, de fatigue-fluage, et de relaxation. Sa particularité est de contenir plusieurs essais de fluage (27 essais), réalisés à différentes températures (550°C à 650°C) et sur trois différentes géométries. La pertinence du paramètre C* pour décrire la fissuration en fluage est analysée par une étude systématique des singularités de contraintes en élasto-viscoplasticité sous plusieurs modes de chargement (différents taux de triaxialité). Il a été montré que, outre le paramètre C*, un deuxième terme non singulier, noté Q*, est nécessaire pour décrire les variables locales au voisinage de la pointe de fissure. Les valeurs de ce paramètre de confinement sont toujours<br />négatives. La conséquence en est que les conditions dapplication des lois classiques de fissuration corrélant la vitesse de fissuration et le paramètre C* (da/dt - C*), seront sécurisantes pour les applications industrielles.<br />A travers cette étude, on a également montré que pour les aciers ferritiques, la période d'incubation des fissures est importante, donc une corrélation de type Ti - C* a été retenue pour prédire le temps à l'amorçage. Pour l'acier inoxydable austénitique, la phase pertinente est celle de<br />la propagation des fissures, ainsi une courbe maîtresse (da/dt - C*) a été établie pour cet acier. Pour cette identification une nouvelle méthodologie de dépouillement des essais de fluage a été mise en place. Enfin, la propagation des fissures a été simulée numériquement par la technique de relâchement des noeuds, permettant ainsi, de valider les expressions analytiques retenues pour dépouiller les essais de fluage.