Deux approches ont été utilisées pour étudier le comportement des sols compressibles lors de l'élargissement d'un remblai : étude sur modèles réduits centrifugés et analyse numérique en éléments finis. Pour l'étude en centrifugeuse, la construction des remblais pendant la rotation du modèle a été effectuée grâce à une trémie développée pour cette recherche. Sept conteneurs d'argile Speswhite reconstituée ont été préparés et équipés de capteurs de tassements de surface, de pressions interstitielles et de déplacements horizontaux. Sept essais ont été réalisés, les premiers pour mettre au point les procédures d'essai et les quatre derniers pour comparer deux variantes d'élargissement d'un remblai. L'analyse des résultats des mesures et des observations faites sur les modèles a montré l'importance des conditions de préparation des massifs d'argile et d'un suivi détaillé de toute l'histoire des sols depuis leur mise en place dans le conteneur jusqu'à la fin de l'essai. Le comportement des argiles pendant la centrifugation dépend d'une combinaison de consolidation et de fluage qui dépend de leur état de surconsolidation initiale et complique l'interprétation des mesures effectuées sur le modèle. La modélisation numérique a été effectuée au moyen de CESAR-LCPC. Elle a porté sur les modèles testés en centrifugeuse et sur le comportement d'un remblai autoroutier élargi vingt-cinq ans après sa construction. Le calcul en éléments finis, de type élastoplastique couplé, nécessitait des données géotechniques qui n'étaient pas toutes disponibles et les comparaisons effectuées avec les mesures sont restées pour cette raison incomplètes, bien que les phénomènes soient reproduits qualitativement de façon acceptable. Ces recherches montrent qu'il est indispensable de mieux définir les propriétés mécaniques des sols pour déboucher sur des résultats d'application pratique.

Cette thèse examine la réponse de différents types de lois de comportement en chargement et en déchargement et modélise numériquement les problèmes du creusement dans un massif du sol. Deux méthodes de modélisation d'excavation par éléments finis, l'une proposée par Ghaboussi et Pecknold (1984) et Brown et Booker (1985) et l'autre par Borja et al. (1989) ont été programmées dans le code CESAR-LCPC. Ces deux méthodes ont été ensuite comparées pour la simulation par étapes de l'excavation de massifs de sol élastique et élastoplastique. Nous avons également développé la notion de chargement-déchargement-rechargement multiaxial basée sur le signe du travail incrémental. L'application de ce principe (qui englobe le cas unidirectionnel) à une excavation pour différents types de lois de comportement, nous permet de connaître l'étendue des zones en chargement et/ou en déchargement autour de l'ouvrage. Un modèle élastique non-linéaire (Wong et Duncan), un modèle élastoplastique parfait à critère de Mohr-Coulomb et un modèle élastoplastique avec écrouissage (Nova) ont été étudiés plus en détail dans cette étude. Ces modèles ont été d'abord appliqués à la simulation d'essais triaxiaux, puis à un problème d'excavation. Ce travail s'achève par la modélisation d'un ouvrage réel : une excavation réalisée devant un rideau de palplanche butonné dans un milieu sableux (Université de Karlsruhe). Dans cette partie, nous avons comparé avec les observations sur site, les résultats des modélisations en éléments finis avec le modèle élastoplastique à critère de Mohr-Coulomb et avec le modèle de Nova.

Cette thèse étudie les écoulements hydrauliques autour des tunnels creusés en terrains aquifères et analyse les déformations induites par la consolidation du massif en présence de sols compressibles. Le mémoire comporte plusieurs synthèses bibliographiques portant sur les problèmes liés aux écoulements hydrauliques vers un ouvrage souterrain, sur les méthodes de calculs existantes pour l'étude de ces écoulements, sur la réponse hydromécanique des massifs de sols au creusement de tunnel (observations de chantiers), sur les techniques proposées pour simuler la construction d'un tunnel par la méthode des éléments finis et sur les approches de calculs développées pour le traitement de la consolidation d'un massif en présence d'une nappe à surface libre. Plusieurs études numériques ont été réalisées à l'aide du progiciel CESAR-LCPC. La méthode des éléments finis a été utilisée pour analyser les écoulements autour de tunnels en conditions bidimensionnelles et tridimensionnelles, en régimes permanent et transitoire. Différentes configurations sont étudiées, considérant un massif indéformable, pour étudier le rôle de l’anisotropie de perméabilité des massifs, leur stratification, et le rôle du revêtement de la paroi du tunnel et de différentes dispositions proposées pour atténuer l'effet des écoulements au front en cours de construction. La réponse d'un massif de sols compressibles au creusement d'un tunnel a été analysée par des calculs hydromécaniques couplés et un modèle de comportement élastoplastique adapté aux sols anisotropes MELANIE. Les résultats de calcul sont comparés à ceux obtenus par des analyses drainée et non drainée. L'approche non drainée a, d'autre part, été adoptée pour étudier la réponse des massifs relativement raides à la méthode de creusement par prédécoupage et pour analyser le comportement de la prévoûte au cours des différentes phases de la construction. La thèse présente, de plus, une contribution au développement du code de calcul CESAR-LCPC. La programmation du modèle MELANIE dans CESAR-LCPC a été adaptée pour traiter la variation de la contrainte de préconsolidation et des modules d'élasticité avec la profondeur, et pour actualiser le coefficient de perméabilité du sol en fonction de la variation de l'indice des vides. Une approche est, par ailleurs, proposée pour le traitement de la consolidation en présence d'un niveau piézométrique variable en utilisant les généralisations de Richards (1931) et de Bishop (I960) pour le traitement de l'écoulement et du couplage au-dessus de la surface libre. Cette approche a été programmée dans un nouveau module s'intégrant dans le code CESAR-LCPC et a été appliquée à l'étude de la consolidation d'un massif sous le chargement d'une fondation et suite à la réalisation d'un tunnel. Ce travail s'achève par la modélisation d'un ouvrage instrumenté : le tunnel de Grimsby, creusé dans un massif argileux compressible. Les calculs ont été réalisés avec différentes lois de comportement (MELANIE, élasticité isotrope et anisotrope et élastoplasticité isotrope), et avec différentes hypothèses pour le traitement du couplage au-dessus du niveau piézométrique ; les résultats ont été comparés aux tassements de surface mesurés sur une période de onze ans après la construction.

Over the last decades, there has been considerable attention in the accurate mathematical modeling and numerical simulations of free surface wave propagation in near-shore environments. A physical correct description of the large scale phenomena, which take place in the shallow water region, must account for strong nonlinear and dispersive effects, along with the interaction with complex topographies. First, a study on the behavior in nonlinear regime of different Boussinesq-type models is proposed, showing the advantage of using fully-nonlinear models with respect to weakly-nonlinear and weakly dispersive models (commonly employed). Secondly, a new flexible strategy for solving the fully-nonlinear and weakly-dispersive Green-Naghdi equations is presented, which allows to enhance an existing shallow water code by simply adding an algebraic term to the momentum balance and is particularly adapted for the use of hybrid techniques for wave breaking. Moreover, the first discretization of the Green-Naghdi equations on unstructured meshes is proposed via hybrid finite volume/ finite element schemes. Finally, the models and the methods developed in the thesis are deployed to study the physical problem of bore formation in convergent alluvial estuary, providing the first characterization of natural estuaries in terms of bore inception. / Ces dernières décennies, une attention particulière a été portée sur la modélisation mathématique et la simulation numérique de la propagation de vagues en environnements côtiers. Une description physiquement correcte des phénomènes à grande échelle, qui apparaissent dans les régions d'eau peu profonde, doit prendre en compte de forts effets non-linéaires et dispersifs, ainsi que l'interaction avec des bathymétries complexes. Dans un premier temps, une étude du comportement en régime non linéaire de différents modèles de type Boussinesq est proposée, démontrant l'avantage d'utiliser des modèles fortement non-linéaires par rapport à des modèles faiblement non-linéaires et faiblement dispersifs (couramment utilisés). Ensuite, une nouvelle approche flexible pour résoudre les équations fortement non-linéaires et faiblement dispersives de Green-Naghdi est présentée. Cette stratégie permet d'améliorer un code "shallow water" existant par le simple ajout d'un terme algébrique dans l'équation du moment et est particulièrement adapté à l'utilisation de techniques hybrides pour le déferlement des vagues. De plus, la première discrétisation des équations de Green-Naghdi sur maillage non structuré est proposée via des schémas hybrides Volume Fini/Élément Fini. Finalement, les modèles et méthodes développés dans la thèse sont appliqués à l'étude du problème physique de la formation du mascaret dans des estuaires convergents et alluviaux. Cela a amené à la première caractérisation d'estuaire naturel en terme d'apparition de mascaret.

Les argiles raides sont généralement rencontrées à grande profondeur (plus de 300 mètres). Ces matériaux se caractérisent par leur teneur en eau relativement faible. De nombreux projets industriels justifient l'intérêt récent porté à ces matériaux. La thèse traite en particulier d'argiles raides susceptibles de répondre à des sollicitations par des déformations élastiques, plastiques et viscoplastiques. C'est la prise en compte simultanée d'une rupture engendrant des déformations plastiques et des déformations différées de types viscoplastiques qui font l'originalité de cette thèse. Cette nouvelle approche constitue une amélioration indispensable dans la description du comportement mécanique d'argiles raides. Ce travail est présenté en trois parties : la première partie est consacrée à l'étude expérimentale et à la modélisation du comportement mécanique des argiles raides. Dans cette partie, nous décrivons les matériaux considérés, et les essais réalisés. Nous discutons les résultats et élaborons un modèle viscoplastique avec rupture. La deuxième partie, « Application du modèle VPR au calcul d'un tunnel de section circulaire », est consacrée à l'élaboration d'une solution originale semi analytique et d'un algorithme implanté dans GEOMEC91. La troisième partie, « Application du modèle VPR au calcul 2d axisymétrique des tunnels soutenus et étude paramétrique », montre très clairement l'influence du modèle VPR sur la convergence du tunnel au moment de la pose du soutènement, et par conséquent sur les sollicitations de ce dernier. Les résultats des calculs montrent une forte influence de la cohésion à court terme sur la convergence du tunnel.

La détermination de la forme d'un obstacle élastique immergé dans un milieu fluide à partir de mesures du champ d'onde diffracté est un problème d'un vif intérêt dans de nombreux domaines tels que le sonar, l'exploration géophysique et l'imagerie médicale. A cause de son caractère non-linéaire et mal posé, ce problème inverse de l'obstacle (IOP) est très difficile à résoudre, particulièrement d'un point de vue numérique. De plus, son étude requiert la compréhension de la théorie du problème de diffraction direct (DP) associé, et la maîtrise des méthodes de résolution correspondantes. Le travail accompli ici se rapporte à l'analyse mathématique et numérique du DP élasto-acoustique et de l'IOP. En particulier, nous avons développé un code de simulation numérique performant pour la propagation des ondes associée à ce type de milieux, basé sur une méthode de type DG qui emploie des éléments finis d'ordre supérieur et des éléments courbes à l'interface afin de mieux représenter l'interaction fluide-structure, et nous l'appliquons à la reconstruction d'objets par la mise en oeuvre d'une méthode de Newton régularisée.

La détermination de la forme d'un obstacle élastique immergé dans un milieu fluide à partir de mesures du champ d'onde diffracté est un problème d'un vif intérêt dans de nombreux domaines tels que le sonar, l'exploration géophysique et l'imagerie médicale. A cause de son caractère non-linéaire et mal posé, ce problème inverse de l'obstacle (IOP) est très difficile à résoudre, particulièrement d'un point de vue numérique. De plus, son étude requiert la compréhension de la théorie du problème de diffraction direct (DP) associé, et la maîtrise des méthodes de résolution correspondantes. Le travail accompli ici se rapporte à l'analyse mathématique et numérique du DP élasto-acoustique et de l'IOP. En particulier, nous avons développé un code de simulation numérique performant pour la propagation des ondes associée à ce type de milieux, basé sur une méthode de type DG qui emploie des éléments finis d'ordre supérieur et des éléments courbes à l'interface afin de mieux représenter l'interaction fluide-structure, et nous l'appliquons à la reconstruction d'objets par la mise en oeuvre d'une méthode de Newton régularisée. / The determination of the shape of an elastic obstacle immersed in water from some measurements of the scattered field is an important problem in many technologies such as sonar, geophysical exploration, and medical imaging. This inverse obstacle problem (IOP) is very difficult to solve, especially from a numerical viewpoint, because of its nonlinear and ill-posed character. Moreover, its investigation requires the understanding of the theory for the associated direct scattering problem (DP), and the mastery of the corresponding numerical solution methods. The work accomplished here pertains to the mathematical and numerical analysis of the elasto-acoustic DP and of the IOP. More specifically, we have developed an efficient numerical simulation code for wave propagation associated to this type of media, based on a DG-type method using higher-order finite elements and curved edges at the interface to better represent the fluid-structure interaction, and we apply it to the reconstruction of objects with the implementation of a regularized Newton method.