A l'opposé du dimensionnement de fondations d'ouvrages sur les sols, régi par des règles de calcul validées, le dimensionnement de fondations au rocher est mal maîtrisé. Jusqu'à présent, ceci n'a pas posé de problème majeur, car le dimensionnement de fondations se trouvait plutôt limité par la résistance du béton que par celle du massif rocheux. Mais la construction d'ouvrages d'art de plus en plus majestueux nécessite aujourd'hui une meilleure connaissance de ce domaine de la mécanique des roches. La première partie recense les méthodes les mieux adaptées pour déterminer la déformabilité et la résistance d'un massif rocheux. Elle présente les différentes méthodes de dimensionnement de fondations superficielle et semi-profonde au rocher. Enfin, elle analyse les textes réglementaires français et étrangers. Il existe peu de méthodes de dimensionnement de fondations qui tiennent compte du caractère discontinu des massifs rocheux ; le cas de fondations soumises à des efforts latéraux et à des moments renversants n'est quasiment pas traité. La deuxième partie consiste en l'analyse numérique du comportement de fondations superficielle et semi-profonde sur massif rocheux à l'aide du code de calcul par éléments distincts UDEC. Les résultats montrent que les discontinuités du massif rocheux ont un rôle primordial sur le comportement sous effort incliné de fondations semi-profondes. Selon leur pendage, les discontinuités peuvent diminuer considérablement la charge limite admissible. Quelques modèles analytiques simples à peu de blocs sont ensuite proposés, afin d'estimer - pour un dimensionnement préliminaire - la charge limite de rupture d'un puits marocain sous effort latéral. Les résultats numériques obtenus sont bidimensionnels; les résultats tridimensionnels peuvent être estimés à l'aide de coefficients de transfert 2D/3D, qu'il serait intéressant de vérifier à l'aide d'un code de calcul par éléments distincts tridimensionnel. Enfin, il serait nécessaire de valider les résultats de modélisation numérique sur un site réel.

Les phénomènes d'érosion sont une cause majeure dans la rupture des ouvrages hydrauliques tels que les digues ou les barrages. Le territoire compte 10000 km de digues, d'où un enjeu important en terme de risque d'inondation. Une approche plus générale est essentielle à une compréhension fine des phénomènes. La modélisation des mécanismes de l'érosion doit permettre l'élaboration de solutions pour sécuriser les ouvrages et les rendre moins sensibles à ces phénomènes. Ce travail est consacré à la modélisation d'une interface entre un milieu de sol et un milieu d'écoulement. Le mouvement de cette interface est régi par l'érosion du sol sous l'effet de l'écoulement. L'écoulement est calculé par les équations de Stokes munies d'un terme de pénalisation permettant la représentation unifiée des comportements des deux milieux. Dans une première partie, nous aborderons la problématique de la représentation de l'érosion. La seconde partie traite plus précisément de la représentation unifiée des comportements des milieux, notamment grâce à l'utilisation des méthodes de pénalisation. La troisième partie constitue une revue des principales méthodes permettant le suivit d'une interface. Une description précise de la méthode Level Set est donnée. Dans la quatrième partie, nous abordons la résolution numérique des équations d'écoulement. La cinquième partie rassemble l'ensemble des tests de la méthode. Enfin, la dernière partie est consacrée à l'application du modèle sur trois cas : le " Hole Erosion Test ", le phénomène de suffusion et l'érosion d'un lit de rivière autour d'une pile de pont.

Le cadre général de ce travail concerne la modélisation du comportement mécanique sous impact de merlons de protection contre les chutes de blocs rocheux. Les ouvrages ciblés sont à technologie cellulaire, et sont composés de gabions à enveloppe grillagée remplis de matériaux granulaires (pierres concassées ou mélange de sable et de déchiquetas de pneus). L'absence de recommandations relatives au dimensionnement de ce type de structure de protection a conduit à la mise en place du projet ANR REMPARe couplant expérimentations in situ et modélisations numériques. Dans le cadre de ce projet, une méthode numérique innovante couplant les approches discrète et continue a été développée. Localement, la méthode des éléments discrets est employée pour conserver la précision de modélisation dans les zones fortement sollicitées proches de l'impact. Le remblai technique aval, peu sollicité, a été discrétisé par une méthode continue aux différences finies. La jonction entre ces deux zones fait l'objet d'un couplage basé sur des critères de conservation d'énergie et de condition cinématique en déplacement. Un modèle de comportement pour les gabions a été proposé et calibré à partir d'essais expérimentaux réalisés à l'échelle de la cellule. Des simulations d'impacts sur un merlon de protection prototype ont été menées puis validées par comparaison avec les mesures issues des expérimentations en vraie grandeur. Afin d'améliorer la compréhension des mécanismes engendrés, l'étude paramétrique portant sur les propriétés des matériaux et les caractéristiques d'impact a permis d'évaluer leur influence sur le mécanisme d'impact. / This study focuses on numerical modelling of the mechanical behaviour of rockfall protection structure submitted to impact. The concerned structures are constituted by cell assembly composed by wire mesh cages filled with granular material (crushed stone, mixture of sand and shredded tires). The lack of design norms led to the set up of the ANR REMPARE project which is based on real scale experiment and numerical modelling. In this project, an innovative numerical method, combining continuous and discrete approaches, has been developed. Locally, the discrete element method is used to maintain the accuracy of modelling in the highly stressed areas near the impacted zone. Downstream, the embankment, lowly stressed, is simulated by a finite difference method. The link between the two domains is ensured by a combined method which is based on energy conservation and respect of border kinematic conditions. A model for the geocells is proposed and calibrated from experimental tests carried out at the cell scale. Impacts on rockfall protection prototype are simulated and validated by comparison with measurements obtained on real scale experiments. Parametric studies on material's properties and characteristics of impact are carried out to evaluate their influence on the impact mechanisms.

L’état de contrainte des galeries et tunnels souterrains évolue dans le temps, tout d’abord en raison de l’avancement du front de taille lors de la construction, mais aussi en raison des changements graduels des propriétés du terrain ou des caractéristiques mécaniques des matériaux constitutifs de l’ouvrage (réduction de la cohésion, vieillissement). Ces changements sont causés par des actions combinées comme: les processus d’altération dus aux infiltrations des eaux météoriques et de réseaux concessionnaires, mais aussi par l’apparition des fissures qui induisent par ailleurs des redistributions de contraintes et de déformations, ainsi que par le fluage des matériaux. Dans cette étude, un travail préliminaire, nous permet de quantifier l’état de contrainte initial lors de la construction des ouvrages. Notre approche basée sur la modélisation éléments finis 2D permet de simuler l’état d’équilibre interne des excavations, nos résultats sont en adéquation avec ceux issus des méthodes analytiques de la littérature. Cet état initial est ensuite calibré via une modélisation 3D. Puis, une macro-modélisation à base d’éléments finis nous permet d’envisager différents scénarios de détérioration des soutènements, notamment à l’extrados des galeries de couloirs d’accès et tunnels de ligne de métro Parisien. Un des scénarios étudié prend en considération la réduction de la résistance mécanique du piédroit et de la zone interface entre le sol et la paroi, ainsi que le fluage de l’encaissant et du revêtement. En utilisant une loi qui permet de simuler la perte de cohésion, paramètre qui évolue en fonction de la déformation inélastique et en fonction du temps, le comportement mécanique des zones dégradées est ainsi reproduit par un modèle viscoplastique de type adoucissant ou à écrouissage négatif. Nous mettons bien en exergue les processus de dégradation observés sur ces galeries durant leur exploitation. Les différentes modélisations permettent de mieux évaluer l’état des différents composants de l’ouvrage selon le type de dégradation et selon le temps traduisant ainsi son vieillissement. Elles permettent également de jauger la sensibilité et la pertinence de différents paramètres. In fine, nous évaluons les pathologies de ces structures, ce qui permet de fournir des indicateurs quant à la gestion de ce patrimoine. / The stress state of underground tunnels and galleries evolves over time , firstly because of the advancement of the working face during construction , then due to gradual changes in soil properties and mechanical properties of materials component of the infrastructures (reduction of cohesion, aging). These changes are caused by combined actions such as: weathering processes due to infiltration of meteoric waters and concessionaires networks, but also by the cracks appearance that induce further redistribution of stress and strain, and the creep. A particular attention is paid to the determination of the initial state of stresses applied on the structure right after the construction. In order to approach this initial stress state which corresponds to the construction phase of underground structures, we propose a methodology based on comparing the numerical model with analytical methods taken from literature. Then we compare 2D and 3D numerical models. Thereafter, a macro- modeling based on the finite elements method allows considering different scenarios of the tunnel-linings deterioration, specifically at the extrados of the galleries and Paris metro line tunnels. One of the studied cases takes into account the mechanical strength loss of the abutment and the interface between the ground and the wall area including their creeps. Therefore, Elasto-viscoplastic with strain softening constitutive model is used in this study to reproduce the degraded zone behavior. The numerical models show that we can reproduce the underground structure behavior during degradation and assess the state of different structure components. The approach used here for the Parisian metro tunnels may be used by the other underground structures managers.

La modélisation à l’aide de l’approche RANS (Reynolds Averaged Navier Stokes) a été menée en trois dimensions, en considérant les régimes permanent et non permanent, dans le but de simuler les écoulements au niveau des jonctions, bifurcations et déversoirs d’orage. A travers ces trois exemples d’étude, plusieurs stratégies de modélisation (différentes combinaison lois de paroi/modèles de turbulence, différents algorithmes de couplage pression/vitesse, différents schémas de discrétisation spatiale, différentes conditions aux limites, différents types et tailles de maille, etc.) ont été testées et évaluées à l’aide de plusieurs indicateurs de performance (de type RMS –Root Mean Square) en s’appuyant sur les données de vitesses (vitesses moyennes in situ et champs de vitesses obtenus par PIV en laboratoire), hauteurs d’eau, débits (répartition de débits en bifurcation en laboratoire ou débits déversés in situ). Les résultats obtenus sont transposables aux autres cas de singularités et ouvrages spéciaux rencontrés en réseau d’assainissement et montrent que : i) les résultats des simulations 3D sont sensibles à la rugosité, aux conditions limites de hauteur et de vitesse ; ii) dans les trois cas d’étude, les schémas de discrétisation du second ordre et l’algorithme de couplage pression/vitesse PISO sont appropriés ; ii) la loi de paroi scalable couplée aux modèles de turbulence de type k-ε pour le cas des jonctions (avec un débit latéral inférieur ou égal au débit principal) ou des déversoirs semblent convenir, tandis que le modèle de turbulence RSM associé à la loi de paroi enhanced ou scalable permet de mieux représenter les écoulements à travers les bifurcations ou au niveau des jonctions lorsque le débit latéral est dominant. Sur la base de ces résultats, une méthodologie de modélisation plus générale définie en six étapes et fondée sur le guide proposée par Jakeman et al. (2006) a été mise au point. La méthodologie ainsi définie a été utilisée pour i) améliorer l’instrumentation du déversoir OTHU (Observatoire de Terrain en Hydrologie Urbaine) situé à Ecully, à partir de la simulation de sa courbe de fonctionnement et en fournissant les incertitudes sur les débits déversés obtenus ; ii) simuler l’implémentation des capteurs débitmétriques à l’aval d’une jonction. Elle a permis de concevoir et de dimensionner un nouveau dispositif de maîtrise des flux d’eau et de polluants déversés (technologie DSM – dispositif de surveillance et de maîtrise des flux déversés au milieu naturel). Ce dispositif a fait l’objet d’un dépôt de brevet international. Enfin, la mise en œuvre de cette méthodologie a été à l’origine de la découverte d’une nouvelle structure d’écoulement dans la branche latérale d’une bifurcation à 90°. L’analyse des résultats des simulations des écoulements mettant en évidence cette nouvelle structure a montré qu’il était possible de prédire l’apparition de cette dernière à partir du rapport d’aspect et du nombre de Froude. / The understanding of sewer flows behaviour is a key component for better urban drainage monitoring and management. However, these flows are conveyed across singularities (such as bends, drops, deviations, etc.) and special structures (combined sewer overflows –CSOs–, channels junction, dividing flow structures, etc.). These singularities and specific structures exhibit complex geometries, leading to open channel turbulent, three-dimensional and multiphase (pollutants and storm and sewer waters) flows. Using three-dimensional CFD (Computational Fluid Dynamics) platform allows a better understanding of mechanisms of contaminants transport through these structures and singularities, leading to a better sewer monitoring. In this thesis, 3D-RANS (Reynolds Averaged Navier Stokes) modelling approach under steady-state conditions is used in order to study flows within CSOs, junctions and bifurcations. Through these three structures, several modelling strategies (wall law/turbulence model combination, velocity/pressure coupling algorithm, spatial discretisation schemes, boundary conditions, computational mesh –shape and size of cells–, etc.) are tested and evaluated thanks to performance indicators (such as RMS –Root Mean Square– indicators) based on velocities (in situ mean velocities and PIV velocity fields obtained in laboratory), water depths and discharge (discharge repartition for bifurcation in laboratory or in situ overflow discharge, for CSOs). Results deriving from these tests are transposable to other singularities or special structures encountered in sewer network and suggest that: i) simulated CFD results are sensitive to the roughness coefficient; ii) for the three studied structures, second-order discretisation schemes and SIMPLE or PISO velocity/pressure coupling algorithm are appropriate; iii) scalable wall function associated to the group of k-ε turbulence model for junction flows (with a lateral inflow lower or equal to the main inflow) or for CSOs is appropriate, whereas RSM turbulence model associated to enhanced wall function allows a better representation of bifurcation flows or junctions flows when the lateral inflow is greater than the main inflow. Based on these results and on Jakeman et al. (2006) guideline, a six steps-methodology focused on the using of RANS approach modelling has been proposed. This six steps-methodology is used in order to i) enhance the monitoring of an OTHU (Observatoire de Terrain en Hydrologie Urbaine) CSO located at Ecully accounting for uncertainties on overflow discharge values ;ii) simulate the performance of flowmeters downstream of a junction, defining the best location for these sensors. This methodology is used to design the new overflow discharge measurement device. This device is an international patent. Finally, the application of the methodology led to point out a new flow structure, occurring in the downstream lateral branch of a 90° bifurcation.

Les instabilités de versants sont en partie responsables de l'apparition de désordres ou de la réactivation de désordres antérieurs dans les structures de tunnels anciens de faible profondeur. Il est cependant nécessaire de bien comprendre ce qui résulte de ces contextes d'interaction et ce qui peut être dû au seul vieillissement rhéologique de ces structures anciennes. Les rétro-analyses réalisées ont conduit ainsi à considérer différentes situations de tunnels entrant dans un versant : tunnel entrant perpendiculairement dans le versant (tête de tunnel), tunnel oblique par rapport au versant, ou tunnel longeant le versant sous une faible couverture. Les pathologies identifiées dans la structure de ces ouvrages varient en fonction de ces situations. Dans cette thèse, différents mécanismes d'instabilités responsables de dommages récurrents ont été étudiés dans le cas d'un tunnel longeant un versant. Pour cela, une partie de la recherche s'est intéressée aux facteurs d'influence à travers une série d'approches numériques paramétriques. Certaines sont relatives à l'altitude du tunnel dans le versant et à sa distance par rapport à la surface libre du versant. Puis sont mis en évidence des facteurs d'influence relatifs aux caractéristiques géomorphologiques du versant et à son évolution naturelle ainsi que des facteurs relatifs à la méthode de construction de l'ouvrage. Dans une seconde partie, des rétro-analyses d'étude de cas ont permis de caractériser différents contextes d'interactions et les pathologies des tunnels anciens concernés. Des simulations numériques ont été réalisées afin d'illustrer le comportement de ces ouvrages et d'expliquer les pathologies identifiées. Enfin dans une troisième partie, une méthodologie d'analyse des différents contextes d'interaction entre tunnels anciens et versants instables a été établie, conduisant à redéfinir le concept de zone d'influence géotechnique dans ce type de situation. / Tunnels and unstable slopesABSTRACT: The slope instabilities are partly responsible for the appearance of disorders or reactivation of previous disorders in old shallow tunnel structures. However, it is necessary to understand what results from these interaction contexts and what is due solely to rheological aging of these old structures. Retro-analyzes have led to the consideration of different situations of tunnels entering a hillslope: tunnel entering a slope (tunnel head), tunnel oblique to the slope or shallow tunnel parallel to the slope. Identified pathologies of structures vary depending on these situations. In this thesis, different mechanisms responsible for damage to structures were studied in the case of a shallow tunnel parallel to a slope. For this, a part of the research has focused on the influence factors through a series of parametric numerical approaches. Some are related to the elevation of the tunnel in the slope and its distance from the free surface of the slope. Then are highlighted influencing factors related to geomorphological characteristics of the hillslope and its natural evolution as well as factors related to the construction method of the tunnel. In the second part, case study back-analyzes have characterized different contexts of interaction and pathologies of old tunnels concerned. Numerical simulations were performed to illustrate the behaviour of these structures and explain tunnel lining damages. Finally in the third part, a methodology for analyzing different contexts of interaction between old tunnels and unstable slopes has been established, leading to redefine the concept of geotechnical influenced zone in this situation.

Cette thèse a pour cadre général l'étude du comportement des ouvrages de génie civil, sous le double aspect, de la loi de comportement non-linéaire des sols et des roches, et de la résolution des problèmes du creusement des tunnels par la méthode de calcul aux éléments finis. Deux catégories de lois de comportement non-linéaires ont été plus particulièrement considérées : modèles hyperboliques (modèle de Duncan et Chang, 1970 ; modèle modifié de Hardin et Drnevich, 1972) et modèles élastoplastiques avec écrouissage (modèle Cam-Clay modifié, 1968) ; dans le dernier, différentes variantes ont été proposées pour tenir compte de la non-linéarité de la partie élastique, trois variantes de type hyperélastique, considérés à partir d'une fonction énergie complémentaire et obéissent aux principes de la thermodynamique Ces modèles ont été appliqués à la simulation du creusement d'une galerie à section circulaire en milieu homogène. Les variantes élastiques non-linéaires du modèle Cam-Clay modifié donnent, dans le cas de contraintes initiales anisotropes, en piédroit des zones en compression et en clé des zones en extension à l'intérieur du massif. L'étude des chemins de contrainte à proximité de la galerie met en évidence l'influence de l'état des contraintes initiales sur les mécanismes différents de rupture induits par le déconfinement. Les modèles retenus ont également été appliqués à la construction du tube Sud du tunnel de Las Planas, creusé dans des marnes. Le comportement de la marne de Las Planas a été étudié à partir d'essais triaxiaux réalisés sous fortes contraintes (essais isotropes et de cisaillement). Les variantes élastiques non-linéaires de type hyperélastique du modèle Cam- Clay modifié et le modèle de Duncan et Chang (1970) conduisent aux résultats les plus satisfaisants pour la plupart des essais de cisaillement. En revanche, la variante élastique linéaire du modèle Cam-Clay modifié surestime systématiquement les déformations volumiques. La modélisation de creusement du tunnel Sud de Las Planas à l'aide des variantes hyperélastiques du modèle Cam-Clay modifié conduit également à des simulations plus représentatives des mesures expérimentales que les autres modèles. Les modèles hyperboliques tendent à surestimer les déplacements induits, par le creusement, dans le massif encaissant.

Les rejets urbains par temps de pluie entraînent des dégradations importantes dans les milieux récepteurs où ils sont rejetés sans traitement préalable. Parmi les campagnes de mesure françaises, effectuées depuis les années 70 à l'exutoire de bassins séparatifs ou unitaires, 17 d'entre elles ont été sélectionnées pour être intégrées dans la base de données constituée au CERGRENE avec l'aide financière de l'Agence de l'Eau Seine-Normandie. Son exploitation a pour objectifs de caractériser la pollution rejetée par temps de pluie (évaluation de masses annuellement rejetées, identification des caractéristiques des événements pluvieux les mieux reliées aux concentrations et aux masses par événement, description des événements les plus pollués, répartition de la masse en cours d'écoulement) et de fournir des éléments sur les dimensions des ouvrages de traitement. Il a fallu tout d'abord reconstituer les caractéristiques des événements pluvieux qui manquaient (volume écoulé, concentration moyenne ou masse produite). Pour cela, des équations de régression ont été établies site par site. Les ajustements obtenus ne sont pas tous de bonne qualité, à cause du nombre réduit de mesures disponibles ou du type des événements utilisés. Quand ils ne sont pas représentatifs de l'ensemble des événements de la campagne, la fiabilité des équations n'est pas assurée et leur domaine d'application ne contient pas tous les événements. Les masses annuelles ont pu être évaluées pour chaque site grâce aux reconstitutions précédentes. L'emploi des équations de régression a permis d'établir des intervalles de confiance autour des valeurs. Ces intervalles sont souvent très larges, et le nombre de sites est assez restreint. Par conséquent, la seule relation qui apparaît est celle qui est la plus évidente et qui existe entre les masses annuelles et les surfaces des sites. Concernant la répartition des masses de pollution en cours d'écoulement, leur forme est très variable d'un événement à un autre, sur un site donné. En moyenne, la pollution contenue dans la première partie de l'écoulement est légèrement supérieure à celle de la deuxième partie. Il y a néanmoins un nombre non négligeable d'événements pour lesquels la majeure partie de la masse de pollution est rejetée en fin d'écoulement. Lorsqu'on compare les sites, la pollution contenue dans les premiers flots (dans 20% de l'écoulement) varie sur les sites unitaires dans des intervalles qui sont les mêmes que sur les sites séparatifs. Pour traiter les rejets urbains de temps de pluie, des ouvrages de décantation (par stockage ou par traitement au fil de l'eau) sont reconnus être des outils efficaces. Nous avons calculé la masse interceptée par des ouvrages fictifs en faisant varier leurs dimensions (respectivement volume de stockage ou débit de traitement suivant le type d'ouvrage). Les reconstitutions utilisées conduisent à des intervalles de confiance très larges. Les résultats dépendent de la période de retour des événements pluvieux qui ont été utilisés. En première approximation, pour intercepter 80% de la masse produite sur une durée d'environ un an, il faut des volumes de l'ordre de 100 m3/ha act ou des débits de l'ordre de 10 1/s/ha act. Pour préciser ces chiffres, il faudrait disposer de longues séries d'événement de plusieurs années. Il serait alors possible en fonction de l'objectif (par exemple intercepter l'événement dont la période de retour est fixée ou intercepter la pollution rejetée pendant une année de pluviométrie moyenne) d'indiquer quelle dimension est la plus efficace.

Le clouage est une technique de renforcement des sols en place par des inclusions passives. Dans la première partie de ce rapport, on présente les résultats expérimentaux et théoriques d'une recherche approfondie sur le mécanisme d'interaction sol-renforcement rigide et la mobilisation de la résistance au cisaillement du sol lors du cisaillement direct de la masse du sol cloué. La méthode des éléments finis a été utilisée pour la simulation numérique des essais de cisaillement direct effectués au CERMES en 1980-1981 sur des éprouvettes de grandes dimensions (40 x 60 x 40 cm) de limon argileux renforcé par des inclusions perpendiculaires à la surface de cisaillement. Cette étude a permis de mettre en évidence la mobilisation d'une cohésion apparente et une diminution de l'angle de frottement interne de la masse du sol renforcé et d'étudier l'effet de certains paramètres (contrainte normale appliquée ; rigidité et nombre de renforcements) sur l’effet des inclusions, le comportement de la masse du sol renforcé et la mobilisation de sa résistance globale au cisaillement. La deuxième partie de ce rapport présente une analyse du comportement des ouvrages de soutènement en sol renforcé. On étudie, d'une part, les remblais en Terre Armée, et d'autre part, les soutènements in-situ par clouage. Si la technique de la Terre Armée a fait l'objet de plusieurs recherches et observations sur ouvrages en vraie grandeur, peu d'expérimentations ont été effectuées jusqu'à présent sur des soutènements par clouage. Ces quelques expérimentations montrent une certaine analogie entre le comportement de ces deux types d'ouvrages. Cependant, il existe des différences fondamentales qui concernent en particulier : - Les modes de constructions, - La rigidité des renforcements, - L'inclinaison des renforcements, - La technique de mise en place des renforcements, - La nature du sol. Notre travail a essentiellement pour objectif d'étudier l'influence de ces paramètres sur le comportement global de l'ouvrage sur les efforts mobilisés dans les inclusions ainsi que sur le déplacement de la paroi. On a utilisé une simulation numérique à l'aide d'un programme de calcul en éléments finis, permettant une modélisation phase par phase de ces deux techniques. Pour vérifier cette modélisation on compare les résultats de calculs, d'une part, avec les mesures sur ouvrages réels, et d'autre part, avec observations sur modèles réduits.

L'objectif de ce travail est de modéliser le comportement hydro-mécanique de massifs rocheux fissurés. Ceci concerne précisément l'amélioration et l'extension d'un modèle par éléments discrets appelé BRIG3D. Ce modèle assimile un massif rocheux fissuré à un assemblage de blocs rigides qui interagissent au niveau de leurs interfaces. La déformation de ces interfaces est reliée aux déplacements relatifs des blocs. Sous l'effet des sollicitations du massif, chaque bloc, en interaction avec d'autres, se déplace vers une position d'équilibre. Les améliorations concernent la reformulation mécanique de BRIG3D. De nouvelles descriptions concernant le mouvement d'un corps rigide, la position des interfaces et le champ des contraintes sont introduites. Ce travail de recherche concerne également l'introduction d'un modèle par éléments frontières décrivant un écoulement laminaire, permanent et plan aux interfaces des blocs. Le couplage de ces deux modèles à permis d'analyser le comportement hydro-mécanique de massifs rocheux fissurés soumis à diverses sollicitations. Une étude particulière a été consacrée à l'analyse du comportement d'une fondation de barrage.