Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Bien qu'assez répandus dans le milieu naturel et dans diverses activités industrielles, les écoulements granulaires denses possèdent des caractéristiques particulières qui semblent dépendre fortement des géométries observées. Lorsque le fond qui supporte la coulée peut être érodé, les mécanismes déterminant la position de l'interface liquide-solide ou la vitesse à laquelle la pente est dévalée ne sont pas encore bien compris. Ceci provient du manque de connaissance de la rhéologie des écoulements granulaires et constitue le coeur de cette thèse.<br /><br />Ce travail s'articule principalement en trois parties. Premièrement, dans le but de fournir une description détaillée des écoulements sur tas, une étude expérimentale systématique des écoulements stationnaires a été menée dans un canal. Nous avons ainsi pu mettre en évidence l'importance des parois latérales du système dans le cadre des écoulements à la surface d'un tas : les propriétés de l'écoulement sont complètement contrôlé par la largeur du canal. Ce constat nous a ensuite conduit à la proposition d'une modélisation fondée sur une rhéologie récemment développée.<br />Cette rhéologie a été ensuite étendue à trois dimensions et appliquée à notre configuration. Des simulations numériques permettent de reproduire correctement les propriétés des écoulements stationnaires. Nous avons enfin étudié un phénomène transitoire particulier afin de comprendre la manière dont l'épaisseur coulante évolue. Malgré la présence de certaines limitations, ce travail montre qu'une description locale des écoulements granulaire est appropriée et que la rhéologie proposée ouvre la voie à des modélisations plus complètes des écoulements granulaires denses.

Le présent travail à caractère théorique et numérique étudie le comportement mécanique des interfaces sol granulaire-structure ainsi que son application aux pieux isolés sous chargement statique. A partir du modèle de comportement élastoplastique d'interface bidimensionnel "MEPI 2D", un modèle tridimensionnel dit "MEPI 3D" est développé afin de mettre en évidence l'effet des contraintes de confinement tridimensionnelles régnant à l'interface au cours du cisaillement. Ces modèles considèrent un écrouissage déviatorique, prennent en compte l'état caractéristique (notion de contractance-dilatance) et intègrent l'état critique aux grands déplacements tangentiels. Les prévisions du comportement de l'interface à l'aide de MEPI 3D sont satisfaisantes. Parallèlement à cette contribution rhéologique, une méthodologie générale de modélisation numérique par éléments finis d'un essai de chargement de pieu est proposée. Les modèles MEPI 2D et MEPI 3D sont implantés dans le code de calcul aux éléments finis CESAR-LCPC. L'analyse numérique par CESAR-LCPC d'un pieu modèle en chambre d'étalonnage (Cermes, type 1), de pieux modèles sur site réel (Imperial College Pile, type 2) et d'un pieu en vraie grandeur (caisson du LCPC, type 3) est menée considérant les modélisations axisymétriques (MEPI 2D, MEPI 3D) et tridimensionnelles (MEPI 3D). En comparant les prévisions numériques aux données expérimentales, il est montré que le comportement des pieux est correctement reproduit. Par ailleurs, la comparaison de l'influence des modèles MEPI 2D et MEPI 3D est abordée.

Notre principale contribution à l'étude du matériau gel a été la mise au point d'un dispositif original de mesure directe de la pression osmotique dans les gels. Des expériences effectuées sur des gels cylindriques d'acide polyacrylique partiellement ionisés ont montré que la pression osmotique varie linéairement avec la concentration en contre ions en dessous du seuil de condensation. Ce résultat est conforme aux théories existantes. Le comportement, sous compression oedométrique, d'empilements lâches de billes de gel gonflées dans une solution de salinité donnée est non linéaire, et la caractéristique force-déformation peut être décrite par morceaux par une loi de puissance. Celle-ci présente trois régimes : un régime de consolidation décrit par un exposant de l'ordre de 1.6, comparable à l'exposant donné par la loi de Hertz. Dans le régime consolidé, le nombre de contacts augmente avec la déformation. Le comportement est alors gouverné par un exposant de l'ordre de 2.5, plus grand que l'exposant de Hertz. Enfin, apparaît un comportement fortement non linéaire décrit par un exposant de l'ordre de 3.5 : tous les contacts sont établis et le comportement est dominé par la contribution des particules fortement déformées. Lors des cycles de compression, l'hystérésis est principalement due à la dissipation visqueuse dans les billes sans expulsion de solvant. Sous compression isotrope générée au sein de l'empilement, contraint à occuper un volume constant, par variation de la taille des particules obtenue en faisant varier la salinité de la solution, le comportement de l'empilement est encore non linéaire avec le paramètre de contrôle de l'expérience. Le comportement observé s'apparente à un régime consolidé, décrit par un exposant supérieur à celui obtenu dans le même régime en compression oedométrique.

Ce travail est destiné à mettre au point une méthode de dimensionnement à l'état limite ultime des structures en béton armé fondée sur le calcul à la rupture. La méthode propose une analyse bidimensionnelle, en contrainte plane, dans laquelle les armatures sont représentées comme des milieux continus unidimensionnels plongés dans le béton considéré comme un milieu continu bidimensionnel (modélisation mixte). Appliquant la méthode à la prise en compte de l'effort tranchant dans le dimensionnement des poutres en flexion, on propose de modéliser le renforcement transversal des poutres par la procédure d'homogénéisation en calcul à la rupture, tout en conservant une modélisation discrète des armatures longitudinales. Ayant choisi des critères de résistance simplifiés pour le béton et les armatures, on montre ensuite les limites d'une approche unidimensionnelle du dimensionnement des poutres par le calcul à la rupture. On développe donc une approche numérique fondée sur la programmation linéaire dans le cadre de la méthode des éléments finis. Cette dernière approche permet ainsi d'évaluer l'influence de l'effort tranchant sur la résistance des poutres et de déterminer le domaine de validité de l'approche unidimensionnelle proposée précédemment. La modélisation et l'approche numérique sont finalement validées par comparaison avec des résultats expérimentaux.