Academic literature on the topic 'Barrages – Modèles mathématiques'

Résumé Les besoins en eau, en perpétuelle croissance, nécessitent une mobilisation des eaux de surface. La construction de barrages, menacés par un envasement précoce, nécessite une exploitation rationnelle, moyennant des outils d’aide à la quantification et à la prédiction des dépôts de sédiments. Ainsi, le comblement des retenues peut être simulé pour différentes périodes d’exploitation. La maîtrise du processus transport-déposition des sédiments, constitue un centre d’intérêt vers lequel convergent plusieurs approches : prédiction, modélisation stochastique et modélisations mathématique et physique. De multiples interactions existent entre les matériaux solides formant le lit, ceux transportés à proximité du fond et ceux se trouvant en suspension, se traduisant par divers modes de transport. Cette contribution présente un modèle numérique qui se prête au calcul par ordinateur dont le comportement morphologique du lit d’un cours d’eau peut être facilement simulé. Les étapes de son élaboration sont décrites avec détail, le code de calcul ainsi produit est mis en valeur sur un site test de validation. Un tronçon de 17 km est étudié, entre un barrage en exploitation et une station hydrométrique, sise à l’amont de celui-ci avec une bathymétrie à l’appui. Le transport solide dans le cours d’eau naturel est représenté par un système d’équations unidimensionnelles décrivant un mélange Eau-Sédiment et traduisant les lois de conservation. La méthode utilisée aux différences finies est appliquée, avec un schéma implicite. Elle est du premier ordre en temps et du deuxième ordre en espace. Les équations des écoulements à surface libre, en régime non permanent et graduellement varié, sont utilisées qui, associées à l’équation de continuité solide, forment le système de Saint-Venant-Exner. Une formule de charriage est utilisée pour exprimer le débit solide. Les résultats obtenus expliquent la contribution du transport solide par charriage dans le comblement de la retenue d’un barrage étudié en zone semi-aride.

Le barrage d’Arzal est un ouvrage clé pour la gestion des niveaux sur le bassin de la Vilaine aval et une réserve d’eau douce alimentant la plus importante unité de production d’eau potable de Bretagne. Cette réserve est multi-usage et peut engendrer des conflits complexes à gérer, notamment en période de sécheresse. Le barrage est équipé d’une écluse, entraînant des intrusions d’eau de mer lors des éclusées, et d’un système de siphons, permettant de capter la majeure partie des eaux salées parasites et de les rediriger vers l’estuaire. Cependant, en étiage, une fraction d’eau salée peut remonter dans le plan d’eau et dégrader la qualité de l’eau à la prise d’eau de l’usine. Dans un contexte de changement climatique, les périodes de sécheresse ainsi que les conflits d’usage pourraient devenir plus intenses. L’amélioration des connaissances des phénomènes complexes de diffusion des eaux salées en amont du barrage permettrait de piloter au mieux l’ouvrage, en validant ou identifiant des pratiques de gestion. L’étude de solutions de limitation transitoires « simples » pourrait également permettre d’optimiser cette gestion. Le projet transversal avec quatre parties prenantes (EPTB Vilaine, imaGeau, Saur et Inria), réalisé de 2016 à 2018, a permis de répondre à ces objectifs. Les résultats de cette étude montrent l’importance des outils de mesure de salinité en continu innovants en temps réel pour l’aide à la décision de restrictions de certains usages si la ressource est menacée, via un retour d’expérience sur trois ans de mesures. Il montre également l’utilisation et la valorisation de la donnée, en mode prédictif, via la construction d’un modèle mathématique 3D de diffusion de la salinité, intégrant l’ensemble des ouvrages. La modélisation a également permis de simuler et d’étudier les performances de solutions de limitation des phénomènes de remontées d’eaux salées afin de choisir la mieux adaptée.

Les zones méditerranéennes et semi - arides sont caractérisées par des régimes pluviométriques irréguliers dans le temps et dans l'espace. L'érosion hydrique des bassins versants et la sédimentation des cours d'eau constituent des phénomènes complexes très répandus qui posent de nombreux problèmes aux ingénieurs chargés des projets de mise en valeur. En raison de son ampleur, cette érosion constitue une contrainte majeure au développement de l'agriculture et à la promotion des activités rurales en Algérie. Une approche statistique de quantification des apports solides est développée dans cette étude. Les données utilisées sont recueillies au niveau de cinq sous-bassins du bassin versant de l'Oued Mina drainés par les stations hydrométriques, en amont du barrage de Sidi M'hamed Ben Aouda. La période d'observation s'étale sur 22 ans (1973 à 1995). Il existe une grande variabilité inter-annuelle et intrannuelle des apports liquides et solides. En plus de cette variation temporelle, il a été remarqué une variabilité spatiale des apports solides et liquides. Les valeurs maximales du transport solide sont enregistrées au début d'automne et à la fin du printemps. Le transport solide d'automne reste le plus élevé, dépassant significativement les autres saisons. Cette variabilité trouve son explication dans la pauvreté du couvert végétal durant cette saison et la nature agressive des pluies d'automne. La charge spécifique du bassin versant de l'oued Haddad est de l'ordre de 212 t/km2/an. Le bassin versant de l'oued Taht a une charge spécifique inférieure à celle du bassin de l'oued Haddad (Ds=191 t/km2/an). Les bassins versants de l'oued el Abd présentent des charges spécifiques de l'ordre de 117 t/km2/an à Ain El Hamara et de 65 t/km2/an à Takhmaret. La charge spécifique de l'oued Mina au barrage de Sidi M'hamed Ben Aouda est de 160 t/km2/an. Des modèles mathématiques liant l'apport solide au coefficient d'écoulement ont été proposés pour les cinq sous-bassins drainés par les stations hydrométriques.

Un modèle numérique et non linéaire pour les problèmes de dynamique des sols est présenté. La formulation variationnelle du problème ainsi que l'approximation numérique des équations sont détaillées. L'influence des accélérations du fluide et du solide sur la propagation des ondes est également discutée. Pour prendre en compte le milieu semi-infini, une frontière absorbante basée sur l'approximation paraxiale est développée pour les milieux poreux satures. Les résultats des applications numériques réalisées a l'aide du code de calcul Gefdyn sont confrontes aux mesures réelles lors de benchmark internationaux sur le calcul des barrages et des effets de sites qui garantissent le niveau des recherches effectuées

Les barrages hydro-électriques construits dans les climats froids subissent des poussées des glaces importantes. Ces charges doivent être prises en compte pour l’évaluation de la sécurité de ces structures. Bien que les barrages aient été construits et exploités dans différents pays pendant de nombreuses années, il n’y a pas de consensus sur les critères de charge de glace. Les valeurs de conception varient considérablement – plus de 3 fois d’une autorité internationale à l’autre. La nature complexe des poussées des glaces implique de nombreux facteurs et incertitudes qui empêchent encore la pleine compréhension des processus en jeu, malgré les avancées récentes. Cette recherche doctorale vise à donner une nouvelle approche aux calculs des poussées des glaces sur les barrages et, ainsi, faire progresser l’objectif d’harmonisation des critères de conception pour les poussées des glaces. L’approche méthodologique de cette recherche consiste à compiler les modèles existants de comportement rhéologique de la glace et les inclure dans la simulation numérique des principaux événements historiques de poussées des glaces mesurés au Canada. La construction d’un modèle numérique commence par un choix de processus et d’hypothèses clés, ainsi que par des modèles rhéologiques à inclure. Par conséquent, une revue compréhensive de la littérature sur le comportement rhéologique de la glace a été réalisée. Cette thèse décrit un modèle de fluage 3D du comportement rhéologique pour deux types de glace (glace columnaire et glace de neige). Le modèle construit peut être implémenté dans un logiciel d’éléments finis pour simuler le comportement de la glace sous les charges statiques. Au cours du projet, le logiciel ANSYS a été utilisé pour construire un modèle d’éléments finis pour des simulations d’événements mesurés. Le modèle d’éléments finis a été validé et étalonné par rapport à des expériences en laboratoire provenant de la littérature. Par la suite, il a été utilisé dans la simulation d’événements de poussées des glaces sur des barrages au Canada. Seule expansion thermique de la glace a été simulée par le modèle d’éléments finis. Les fluctuations de niveau d’eau ont été indirectement incluses dans le modèle par la simulation de la glace confinée latéralement. L’influence de différents facteurs sur la poussée des glaces simulées a été étudiée. Les conseils et les lignes directrices pour le choix et l’implémentation du modèle rhéologique de fluage dans un logiciel d’éléments finis commercial sont également fournis. La thèse est basée sur des résultats d’études qui ont été publiés dans deux articles de colloques et dans un article de journal. La thèse inclut l’introduction des poussées des glaces sur les barrages, les objectifs du projet, les revues de littérature sur les processus influents des poussées des glaces, sur les propriétés rhéologiques des deux types de glace fréquents sur un réservoir et la modélisation du comportement de la glace, ainsi que la présentation des données de terrain, la description du modèle, les études non encore publiées, la conclusion et les recommandations. La thèse discute du lien entre les expériences en laboratoire sur la glace et les mesures in situ des poussées des glaces sur les barrages. L’étude documente le fait que le fluage élastique retardé de la glace a un effet limité sur les poussées des glaces résultantes pour les grandeurs de grains de 3 mm à 20 mm, et son calcul peut être omis dans des évaluations approximatives des poussées des glaces avec cettes grandeurs de grains moyennes. Il est également montré que la variabilité du coefficient de dilatation thermique de la glace avec la température n’est pas importante pour les températures étudiées. L’importance des conditions initiales et du comportement de la glace en tension pour les poussées des glaces résultantes est également montrée. La thèse met l’accent sur l’importance de la considération correcte du confinement latéral et du type de glace pour les calculs de poussées des glaces sur les barrages. Il a été montré que la hausse du niveau d’eau et le confinement latéral de la glace peuvent produire très grandes forces de poussées, si la glace est principalement columnaire. De ce fait, il devient important d’évaluer la possibilité d’avoir des hivers secs ayant peu d’accumulation de neige sur la glace (les conditions favorables à la formation de la glace columnaire) pour la sélection de la charge de conception.
Hydroelectric dams are subjected to ice loads in cold climates. Withstanding such loads is an important safety requirement for the structures. Although dams have been built and operated by different countries for many years, there is no consensus on a design ice load criteria. Design values vary considerably – by more than 3 times from one international authority to another. The complex nature of ice loads involves many factors and uncertainties that still prevent the full understanding of ice field processes in spite of much progress done in this direction. This doctoral research aims to give new insights into reservoir ice loads, and, thus, to move forward the objective of harmonization of ice load design criteria. The main approach of this research is a compilation of existing models of ice material behaviour and their inclusion in the numerical simulation of key historical ice load events measured in Canada. A construction of a numerical model begins with a choice of key processes and key assumptions as well as material models to include. Therefore, a comprehensive literature review on the topic is required. During the project, diverse literature on ice material behaviour was reviewed. This thesis describes a 3D creep material model for rheological behaviour of two types of ice, namely columnar (S2) ice and snow (T1) ice. The constructed model can be implemented into software to simulate behaviour of ice under the static loads. In this project, the commercial software ANSYS was used for this purpose to build a finite element model (FEM) for simulations of field events. After the FEM was built, it was validated and calibrated against laboratory experiments taken from the literature. Subsequently, the resulting FEM was used for simulation of key Canadian field events. Only thermal expansion of ice was simulated by the FEM with indirect accounting of the impact of water level fluctuations through the modelling of ice cover under lateral confinement. The influence of different factors on the simulated ice load was studied, namely, the influence (i) of selected material constants, (ii) of lateral confinement, (iii) of grain size and of delayed elastic creep, (iv) of variable or constant coefficient of thermal dilatation, (v) of ice type, (vi) of initial stresses as well as (vii) the influence of ice behaviour in tension. The advice and guidelines for the choice and implementation of a creep material model into commercial finite element software are also given. The thesis is based on studies conducted for two conference papers and one journal paper. The thesis includes introduction to the ice loads on dams, objectives of the project, literature reviews on governing ice loading processes, on material properties of two types of ice of frequent occurrence on a reservoir and on modelling of ice behaviour as well as field data presentation, model description, yet unpublished studies, conclusion and recommendations. The thesis discusses the linkage between laboratory experiments on ice and in situ measurements of reservoir ice loads on dams through the numerical modelling of field events. It documents the fact that delayed elastic creep has a limited effect on the resulting ice load for grain sizes of 3 mm – 20 mm, and its calculation can be omitted in approximate assessments of field ice loads for ice cover with these average grain sizes. It is also shown that variability of coefficient of thermal dilatation of ice with temperature is not important for the temperature range studied. The importance of initial conditions and ice behaviour in tension for the resulting ice load is also shown. The thesis emphasizes the importance of correct accounting of lateral confinement and ice type for ice load calculations. It is shown that when water level rises and the ice is confined laterally, if the reservoir consists primarily of columnar ice, the loads can be very substantial. As such, it is recommended that engineers pay particular attention to the possibility of having cold dry winters with little snow accumulations on an ice cover (favorable conditions for columnar ice formation) when choosing their design criteria.

Cette thèse vise à modéliser numériquement des barrages en enrochements pendant les phases de construction et de mise en eau. Ceci est fait en se basant sur des essais en laboratoire et l’analyse des données de l'instrumentation insitu. Des essais triaxiaux en compression monotone à déformation/contrainte contrôlée sont effectués sur des assemblages compactés de particules de granit. Les échantillons utilisés en laboratoire ont initialement été préparés à différentes densités sèches et à différentes granulométries. Une loi de comportement élastoplastique parfait, en considérant l’élasticité non linéaire est implémentée dans le logiciel FLAC. Ceci est fait pour modéliser le comportement des particules de roche pendant les essais. Le nouveau modèle est basé sur le comportement hyperbolique, en considérant le critère de rupture de Mohr-Coulomb avec une règle d'écoulement non associé pour considérer la dilatance à cause du cisaillement. La technique de relaxation des contraintes de Nobari-Duncan est également implémentée dans FLAC avec les procédures modifiées et les nouveaux algorithmes pour reproduire les déformations causées par le mouillage. À l'aide de neuf paramètres, le modèle réussi à reproduire le comportement observé pendant les essais triaxiaux et à caractériser l'assemblage des particules de roche. Les effets de temps et de saturation sont également discutés en détail. Finalement, le modèle est utilisé pour la modélisation des étapes de la construction et de la mise en eau du barrage en enrochement LG4 construit au Québec, Canada. Les résultats de la modélisation sont comparés avec les données de l'instrumentation. Aussi, les efforts antérieurs de modélisation de ce barrage et les avantages du nouveau modèle sont discutés.
This thesis contributes to the field of numerical modeling of rockfill dams during construction and impoundment based on laboratory research and field instrumentation. Data is presented from monotonic triaxial compression tests, in both stress and strain-controlled conditions, on compacted assemblages of granite rock particles prepared in the laboratory using different initial dry densities and grain size distribution. A nonlinear elastic-perfectly plastic constitutive model is implemented into the commercial software, FLAC, to simulate the behavior of the rock particles during the tests. The new implemented model is the traditional hyperbolic model coupled with the Mohr-Coulomb plastic failure criteria, and with a non-associated flow rule to consider shear dilatancy. The Nobari-Duncan stress relaxation technique with modified procedures and algorithms is also implemented in FLAC to reproduce the collapse deformations due to wetting. Using nine parameters, the model is successful to capture the observed behavior in the triaxial tests and to characterize the assemblage of rock particles. The effects of time and saturation are also discussed in detail. Finally, the model is used to simulate the construction and impoundment stages of LG4 rockfill dam in Quebec, Canada. The results of the modeling of both stages are compared with the instrumentation data and previous efforts of the modeling of this dam and the advantages of the new model are discussed.

La géométrie et la flexibilité d'un canyon sont les paramètres qui affectent grandement la valeur des périodes naturelles dans les barrages en terre. Le canyon entourant des barrages peut être considéré comme un domaine illimité. Pour prendre en compte ces deux effets, le canyon a été modélisé par SBFEM et le barrage en terre, à géométrie limitée, par FEM. La technique hybride SBFEM-FEM pour l'analyse tridimensionnelle dynamique de l'interaction sol-barrage a été validée avec les résultats disponibles dans la littérature. Comme la matrice de rigidité dynamique du domaine non borné est complexe et dépendante de la fréquence, la méthode classique de superposition de modes n'est pas simple pour le système d'interaction sol-structure. Ainsi, pour obtenir la fréquence propre fondamentale, le barrage a été excité en direction amont-aval. Les périodes naturelles du barrage de terre pour des canyons de formes géométriques et de coefficient de impédance différents ont été obtenues. Ils se sont avérés avoir des effets significatifs sur la période naturelle. Les résultats ont été comparés aux données enregistrées réelles. Il a été constaté que les graphiques proposés dans cette étude peuvent être utilisés par des concepteurs de barrages pour l'estimation des périodes naturelles des barrages en terre dans des canyons de formes et de propriétés matérielles différentes. Plusieurs fonctions d'amplification correspondant à différentes conditions de canyon ont été obtenues en appliquant un déplacement uniforme à la limite du canyon. Une étude approfondie a été réalisée pour examiner les effets de la géométrie et de la flexibilité du canyon sur la réponse en régime permanent du barrage. Ces deux effets ont influencé de manière importante les fonctions d'amplification. Alors que la flexibilité du canyon affecte de manière significative la valeur de la fonction d'amplification maximale, cette valeur ne change pas pour les barrages en terre dans lesquels les canyons ont des formes différentes et la même longueur. De plus, la réponse latérale du barrage de terre dans le domaine temporel a été calculée pour analyser les effets susmentionnés lors d'un tremblement de terre réel. Les fonctions d'amplification proposées ont été utilisées pour comparer les spectres de réponse enregistrés du barrage d'El Infiernillo lors des tremblements de terre de 1966 avec la fonction d'amplification calculée. Un accord raisonnable a été observé entre eux. La méthode linéaire équivalente (EQL) a été implémentée dans le FEM. La technique FEMSBFEM a été étendue pour prendre en compte l'effet du comportement non linéaire des barrages en terre. Il a été observé que le comportement non linéaire affecte grandement la fréquence naturelle, la fonction d'amplification et l'accélération de crête maximale du barrage de terre situé dans les canyons. Les effets de la géométrie et de la flexibilité du canyon sur le comportement non linéaire ont été examinés, et on a vu qu'en augmentant la flexibilité du canyon, l'effet de la non-linéarité était diminué. Le barrage d'El Infiernillo a été modélisé par FEM-SBFEM non linéaire 3D, et une comparaison de la fonction d'amplification de crête obtenue par la méthode proposée avec les données enregistrées montre la précision du FEM-SBFEM non linéaire.
The canyon surrounding a dam can be assumed as an unbounded domain, and the geometry and flexibility of a canyon are parameters that greatly affect the values of natural periods in earth dams. In this thesis, in order to take into account these two effects, canyons are modeled by SBFEM, and earth dams, which have limited geometries, are modeled by FEM. The hybrid FEM-SBFEM technique used for the dynamic three-dimensional analysis of soil-earth dam interactions is validated with results available in the literature. Because the dynamic-stiffness matrix of the unbounded domain is complex and frequency-dependent, the classical mode-superposition method is not straightforward for a soil-structure interaction system, and thus, to obtain their fundamental natural frequencies, the modeled dams were excited in the upstream-downstream direction. The natural periods of earth dams in canyons with different geometries shapes and impedance ratios are obtained, and are found to have significant effects on the dams’ natural periods. The results are compared with actual recorded data, and it is found that the graphs put forward in this study may be used by practical engineers for the estimation of natural periods of earth dams in canyons with different shapes and material properties. Several amplification functions corresponding to different canyon conditions are obtained by applying a uniform displacement at the canyons’ boundaries. A comprehensive study is performed to examine the effects of canyon geometry and flexibility on the steady-state responses of the dams, and it is found that these two effects significantly influence the amplification functions. While the flexibility of the canyon does affect the maximum amplification function value, this value does not change for earth dams in canyons that have different shapes but the same length. In addition, the lateral responses of earth dams in the time domain are computed in order to analyze the aforementioned effects under an actual earthquake. The proposed amplification functions are used to compare the recorded response spectra of the El Infiernillo dam under the two 1966 earthquakes with the calculated amplification function, and a reasonable agreement is observed between them. The equivalent linear method (EQL) is implemented into the FEM, and the FEM-SBFEM technique is extended in order to take into consideration the effect of earth dams’ nonlinear behavior. It is observed that such nonlinear behavior greatly affects the natural frequency, the amplification function, and peak crest acceleration of earth dams located in canyons. The effects of canyon geometry and flexibility on the nonlinear behavior are examined, and it is found that by increasing canyon flexibility, the effect of nonlinearity is decreased. The El Infiernillo dam is modeled by the 3D nonlinear FEM-SBFEM, and comparison of the crest amplification function obtained by the proposed method with the recorded data shows the accuracy of the nonlinear FEM-SBFEM.

Les enrochements sont des matériaux granulaires constitués de blocs rocheux dont la taille peut atteindre plusieurs dizaines de centimètres. Les barrages en enrochements, constitués par ces matériaux grossiers, présentent des déformations relativement importantes au cours du temps et peuvent également tasser au moment de leur remplissage ou lors d'une crue accidentelle. Ces déformations semblent liées à des ruptures de blocs à l'intérieur de l'ouvrage, mais ne sont malheureusement pas connues après la construction.

La prévision de ces tassements est donc particulièrement importante pour s'assurer de l'intégrité structurelle de l'ouvrage à long terme.

Le comportement des matériaux granulaires étant fortement lié à la nature discrète du milieu, un modèle discret est proposé afin de prendre en compte les particularités des enrochements avec des paramètres locaux au sens physique clairs. Les phénomènes de rupture sont les principales caractéristiques des enrochements. La démarche retenue dans ce travail s'est attachée à développer un modèle numérique discret capable de prendre en compte la fissuration progressive et différée des blocs rocheux. Chaque bloc est représenté par un assemblage de particules, initialement liées par une cohésion qui peut diminuer progressivement au cours du chargement. Un modèle d'endommagement interfacial est proposé pour décrire cette décroissance progressive au cours du temps. Le modèle comporte deux surfaces seuil (endommagement et de rupture). Un temps caractéristique d'endommagement règle la vitesse d'endommagement. L'effet de l'eau est introduit par couplage avec les paramètres du modèle de cohésion, par la diminution des forces de pesanteur et du coefficient de frottement local. La modélisation adoptée est de type « Non Smooth Contact Dynamics », où les grains et particules sont supposés rigides. Le modèle est implémenté sur la plateforme logicielle LMGC90.

Des simulations numériques à l'échelle des contacts, à l'échelle du bloc et à l'échelle de l'enrochement sont effectuées. Des premières simulations consistent à simuler un ensemble de blocs rocheux incassables placés dans une colonne d'enrochements progressivement remplie en eau. Les simulations révèlent que la prise en compte du déjaugeage et d'une diminution réaliste du coefficient de frottement local ne génèrent que des déplacements faibles dans le milieu. Ces simulations permettent de vérifier que le phénomène essentiel à l'origine des tassements dans un barrage en enrochement est lié à la rupture des blocs rocheux.

Des simulations de compression de blocs cassables (constitués de particules) et de compression œdométrique d'un ensemble de ces blocs sont réalisées. Les essais montrent que le temps caractéristique d'endommagement permet de mettre à l'échelle le comportement du matériau. Les résultats indiquent que l'histoire des évènements de rupture et des déformations sont bien corrélées. Les simulations numériques d'essais de compression œdométrique d'enrochements modèles de barrage mettent en évidence des réponses assimilables au fluage ou à la relaxation. Ces essais reproduisent les phénomènes de base avec très peu de paramètres.

Le travail de thèse concerne le développement d'un modèle numérique pour la détermination des propriétés hydromécaniques équivalentes pour les roches fracturées. Le rapport de thèse est organisé en 7 chapitres. Le premier chapitre décrit les objectifs de la thèse et présente brièvement une synthèse des travaux réalisés sur les roches fracturées. L'accent est mis sur l'Élément Composite et la détermination des propriétés hydrauliques et mécaniques équivalentes des roches fracturées. Le second chapitre présente les outils utilisés dans le travail de thèse, à savoir la méthode des éléments finis, l'Elément Composite, le couplage écoulement/contraintes, le Réseau de Neurones Artificiels et la méthode de Monte-Carlo. Le troisième chapitre présente le concept de la méthode numérique proposée pour obtenir les propriétés hydrauliques équivalentes des roches fracturées. Elle décrit aussi les données nécessaires à ('application de cette méthode, à savoir celles sur les fractures dominantes. Les chapitres 4 et 5 présentent la détermination des propriétés hydrauliques et mécaniques équivalentes des roches sans prendre en compte le couplage contraintes/écoulement. Le 6ème chapitre présente la détermination des propriétés hydrauliques équivalentes en considérant le couplage contraintes/écoulement. Ce chapitre comporte également la construction d'une base de données sur l'influence du couplage contraintes/perméabilité sur les propriétés hydrauliques équivalentes. Le dernier chapitre décrit un procédé basé sur une combinaison des approches implicite et explicite pour l'analyse d'un barrage poids. Les données de la méthode implicite sont déterminées en utilisant la méthode numérique proposée dans ce travail
This thesis concerns the use of the "Numerical Test System" for the determination of the equivalent hydraulic and mechanical properties of fractured rock masses. The thesis is composed of 7 chapters. The first chapter describes the objective of the thesis and briefly presents the previous studies conducted on fractured rock masses, with a particular focus on the Composite Element Method as weil as the estimation of the equivalent hydraulic and mechanical properties of fractured rocks. The second chapter presents the tools used in the thesis development, namely the finite element method, the composite element method, the coupling theory, the artificial neural network and the Monte Carlo Method. The third chapter presents the concept of the "Numerical Test System (NTS)" for determination of the equivalent hydraulic properties offractured rock. This method requires field data on dominating fractures. Chapters 4 and 5 present the determination of the equivalent hydraulic and mechanical properties and the relevant REYs neglecting the seepage/stress coupling. The 6th chapter presents the determination of the equivalent hydraulic properties considering seepage/stress coupling. This Chapter also establishes a stress/permeability database for describing the relationship between stresses and the equivalent hydraulic property. The last chapter describes a procedure based on a combination of the implicit and explicit approaches for the analysis of a grravity dam. The eqivalent properties of the implicit method are determined using the NTS