No description available.

Le comportement de grands barrages en enrochement d'Hydro-Québec à la suite des tassements différentiels de l'enrochement lors de la mise en eau a amené l'entreprise à se questionner sur les relations contraintes-déformations de l'enrochement et ce qui les influence. Pour le barrage Romaine-2, il était d'autant plus important que l'enrochement soit rigide et ne connaisse pas de grands tassements différentiels puisque le noyau est en béton asphaltique. Il est donc plus étroit et a des pentes plus raides que les ouvrages zonés en enrochement à noyau de till qu'Hydro-Québec à l'habitude de construire. Pour évaluer les propriétés et les tassements de l'enrochement, une revue des facteurs d'influence de ces paramètres est importante et a été réalisée dans le cadre de ce mémoire. Ensuite, plus spécifiquement, l'influence de ces facteurs sur l'enrochement du barrage RO-2 a été analysée. Il a été conclu que l'enrochement utilisé pour la construction du barrage Romaine-2 est de bonne qualité, a été mis en place adéquatement et a montré de hautes masses volumiques, de faibles tassements et des modules de déformation élevés.

La contribution significative de cette thèse se situe dans le domaine de la déformation sismique de différents types de barrages en enrochement. La déformation excessive permanente à la crête d'un barrage causée par les secousses du sol menacera la sécurité du barrage. Les caractéristiques de ces barrages et la gravité des tremblements de terre (quantifiée par des mesures d'intensité, c'est-à-dire des IM) sont identifiées comme les composants les plus critiques impliquées dans la performance sismique du barrage. Dans la première partie, une étude est réalisée sur la déformation sismique (tassement de la crête) des barrages en enrochement avec masque en béton (CFRD) et sa corrélation avec les IM de mouvements de sol. De cette manière, l'importance des principales caractéristiques des charges cycliques dans une étude de cas CFRD typique modélisée numériquement est soulignée. On en conclut qu'en plus de l'amplitude et de la durée du mouvement de sol, la réponse du barrage est fortement affectée par les caractéristiques de fréquence du séisme. En rassemblant les histoires de cas de CFRD disponibles et en analysant leur tassement sismique, un nouveau IM est proposé pour décrire efficacement le potentiel d'un séisme à induire un règlement dans les CFRD. Une nouvelle relation prédictive est alors établie qui relie la IM suggérée et l'amplitude du tremblement de terre avec le tassement de crête étudié des cas CFRDs. Dans la deuxième partie, les données disponibles pour les barrages en enrochement à noyau de terre (ECRD) comprenant les valeurs de tassement rapportées et les mouvements de sol enregistrés lors des tremblements de terre sont analysées. Deux nouvelles approches sont développées en introduisant deux nouveaux IM et en les corrélant au règlement observé des cas ECRD. Les IM proposés pour les ECRD prennent en compte l'influence des aspects essentiels concernant le comportement non linéaire des barrages lors de séismes violents: dégradation de la rigidité des matériaux, augmentation de l'énergie dissipée et augmentation des périodes de vibration des barrages. Ces problèmes sont affectés par la déformation induite par le tremblement de terre, à travers laquelle la corrélation avec la déformation d'un barrage est découverte. Les analyses ont montré que l'occurrence du tassement est directement affectée par les amplitudes de déformation de cisaillement. Compte tenu de l'importance de la déformation de cisaillement, la dernière partie de cette thèse est consacrée aux caractéristiques dépendantes de la déformation des ECRD. À cette fin, des études empiriques ont été menées sur plusieurs séries chronologiques d'accélération pour les ECRD qui avaient été frappés par des tremblements de terre au Japon. Une série de graphiques et de relations a été établie pour estimer: (i) l'amplitude de la déformation de cisaillement induite; (ii) la diminution du module de cisaillement des noyaux des barrages; et (iii) augmentation des périodes fondamentales des barrages. Cette thèse fournit un assemblage approprié d'outils qui peuvent être utilisés dans la pratique de l'ingénierie, que ce soit à des fins d'analyse dynamique ou de conception. Les graphiques et les relations présentés ici sont basés sur l'examen et l'analyse de la performance réelle de nombreux barrages sous les tremblements de terre. Ils répondent non seulement aux lacunes des méthodes empiriques précédentes, mais sont également précis et efficaces. / The significant contribution of this thesis is in the area of earthquake-induced deformation of different types of rockfill dams. The permanent excess deformation at a dam's crest caused by ground-shaking will threaten the dam's safety. The characteristics of these dams and the severity of the earthquakes (quantified by intensity measures, i.e., IMs) are identified as the most critical components involved in the dam's seismic performance. In the first part, a study is performed on the seismic deformation (crest settlement) of concrete-face rockfill dams (CFRDs) and its correlation with ground motion IMs. In this way, the importance of cyclic loads' main characteristics on a numerically modelled typical CFRD case study is underlined. It is concluded that in addition to the amplitude and duration of ground motion, the dam's response is strongly affected by the frequency characteristics of the earthquake. By gathering available CFRD case histories and analyzing their seismic settlement, a new IM is proposed to efficiently describe an earthquake's potential to induce settlement in CFRDs. A new predictive relationship is then established which relates the suggested IM and earthquake magnitude with the surveyed crest settlement of the CFRD cases. In the second part, the data available for earth-core rockfill dams (ECRDs) comprising the reported settlement values and the ground motions recorded during the earthquakes are analyzed. Two novel approaches are developed by introducing two new IMs and correlating them to the observed settlement of the ECRD cases. The IMs proposed for ECRDs take into account the influence of the essential aspects regarding the nonlinear behaviour of the dams during severe earthquakes: stiffness degradation of the materials, increase in dissipated energy, and increase in the dams' periods of vibration. These issues are affected by the earthquake-induced strain, through which the correlation with a dam's deformation is discovered. The analyses showed that the occurrence of settlement is directly affected by the shear strain amplitudes. Given the importance of the shear strain, the last part of this thesis is devoted to the strain-dependent characteristics of ECRDs. To this end, empirical studies were performed on several acceleration time series for ECRDs that had been struck by earthquakes in Japan. A series of graphs and relationships were established to estimate: (i) the amplitude of the induced shear strain; (ii) the decreased shear modulus of the dams' cores; and (iii) increased fundamental periods of the dams. This thesis provides an appropriate assembly of tools that can be utilized in engineering practice, either for dynamic analysis or design purposes. The graphs and relationships presented herein are based on reviewing and analyzing the actual performance of numerous dams under earthquakes. They not only address the shortcomings of previous empirical methods, but are also accurate and efficient.

Thèse (M.Sc.Géogr.)--Université Laval, 2007. / Titre de l'écran-titre (visionné le 5 mai 2008). Bibliogr.

Le projet de recherche a été réalisé dans le but premier de mettre à niveau, d'adapter et de rendre opérationnelle la méthode de Newmark pour les types de barrages et les types de séismes habituellement rencontrés au Québec. Pour ce faire, une Banque de Signaux Sismiques pour l'Est du Canada (BDSSEC), contenant 451 accélérogrammes, a été produite de même qu'un code de calcul permettant de calculer les déplacements permanents, par double intégration numérique, selon la méthode de Newmark. À partir d'analyses statistiques, des modèles prédictifs applicables au contexte sismologique de l'Est canadien ont été développés en considérant un comportement rigide-plastique (méthode dite du "bloc rigide") et un comportement viscoélastique linéaire équivalent (méthode dite "flexible"). La comparaison entre le nouveau modèle rigide-plastique pour l'Est avec ceux existants développés à partir de signaux de l'Ouest américain a montré que les déplacements permanents évalués pour l'Est peuvent être significativement inférieurs à ceux de l'Ouest pour des magnitudes et des distances épicentrales comparables compte tenu de leur contenu fréquentiel plus fort en hautes fréquences. Le modèle prédictif flexible a montré qu'il pouvait prédire des déplacements permanents du même ordre de grandeur que Makdisi et Seed (1977) et de Rathje et Antonakos (2011). Les déplacements permanents sont amplifiés par la méthode flexible comparativement à la méthode rigide principalement lorsque le ratio entre la période prédominante de l'ouvrage et la période moyenne de la sollicitation se situe autour de 1. Le second but consistait à tester l'applicabilité de la méthode du bruit de fond ambiant (H/V) comme outil expérimental pour estimer la période prédominante des barrages en remblai et en enrochement (Tp). Les études comparatives ont montré, de façon générale, une bonne concordance entre les calculs analytiques et les mesures expérimentales tant et aussi longtemps que les hypothèses à la base des modèles sont respectées.

La gestion des interventions sur les ouvrages en béton, publics et privés, n’est pas toujours optimale. Ceci est dû à l’absence d’un protocole complet capable de guider de manière homogène les ingénieurs dans leurs choix relatifs à la priorité ainsi qu’aux types de travaux à effectuer sur leurs ouvrages. Ce manque d’assistance décisionnelle se traduit par des interventions qui ne sont pas toujours optimales en termes de coûts et de durabilité. Ce mémoire présente la méthodologie développée en collaboration avec les ingénieurs d’Hydro-Québec dans le but d’uniformiser, rationaliser et intégrer dans un outil de gestion les différentes étapes requises pour réussir des interventions optimales et de qualité. Trois activités sont essentielles afin de systématiser l’approche de réfection ou d’entretien sur les ouvrages en béton: le diagnostic des dégradations, l’évaluation de la pertinence des interventions et le choix des méthodes et produits de réfection appropriés. Après avoir fait une revue de documentation exhaustive afin de bien cerner le problème, les chapitres de ce mémoire décrivent successivement les approches retenues pour chaque activité en vue d’orienter efficacement les ingénieurs dans la gestion de leurs infrastructures. Plus précisément, un outil informatique a été développé de manière à incorporer les méthodes qui résultent de ce projet de recherche. Un logiciel de ce genre est essentiel, puisqu’en plus de simplifier l’utilisation à grande échelle des divers outils d’évaluation développés, il permet d’uniformiser les méthodologies de travail à l’intérieur d’une organisation de grande taille comme Hydro-Québec. / Hydro-Québec is amongst Canada’s largest hydroelectric companies. It is responsible for maintaining an acceptable level of performance of many concrete dams. In order to do so efficiently, from cost and durability standpoints, the different stages of the condition survey required in the evaluation of the structure must be homogeneized and integrated in a tool capable of assisting the engineer in his choices. These stages are the diagnosis of the causes of concrete distresses, the evaluation of the need of repair or maintenance actions and the selection of the appropriate materials and repair procedures. This document describes the methodologies developed and later incorporated in a computer program, which is able to guide the engineer in a systematic and rational way.

Hydro-Québec reçoit plusieurs demandes de la part de compagnies minières pour pouvoir utiliser la crête de leurs ouvrages en remblai comme route d'accès aux ressources naturelles. La problématique liée à ces demandes d'autorisation de passage est le fait que lors de la conception d'un barrage en remblai, la circulation des charges lourdes en crête n'est pas considérée. Il est donc pertinent de développer des outils décisionnels pour les gestionnaires des barrages de manière à ce qu'ils soient en mesure d'octroyer ou d'appliquer certaines restrictions sur les passages demandés. Ces critères décisionnels doivent être basés, entres autres, sur la déformation permanente subie par le noyau du barrage. Concrètement, des instruments de mesure de contraintes et de déformations ont été installés dans l'extrémité supérieure du noyau d'argile du barrage en remblai de St-Timothée. Les contraintes horizontales et verticales mesurées lors du passage de camions en crête permettent de reproduire le cheminement de contraintes en laboratoire lors des essais de déformations permanentes. L'échantillon d'argile recompacté en laboratoire est soumis à des contraintes déviatorique cycliques allant jusqu'à des valeurs très élevées par rapport à ce qui a été enregistré sur le terrain. La contrainte axiale cyclique à laquelle l'argile accumule une certaine déformation permanente suite à sa post-compaction est de 5 kPa et celle à laquelle l'instabilité du matériau est observée correspond à 55 kPa. La charge en crête de l'ouvrage qui engendre ces incréments de contraintes axiales au sein du noyau est déterminée avec les données enregistrées avec les instruments mis en place dans le barrage. Le nombre de passages nécessaire pour engendrer des déformations permanentes de 50 mm et 100 mm au sein du noyau peut être estimé en fonction de différentes charges circulant en crête du barrage. Ceci est effectué en considérant une variation de la profondeur dans le noyau sur laquelle la charge en crête contribue à sa déformation permanente.

L'érosion régressive est l'un des mécanismes initiaux de l'érosion interne qui se produit dans les fondations des barrages en remblai et des digues. Une fois que la charge hydraulique de l'écoulement souterrain atteint des conditions critiques, des particules de la couche granulaire se détachent et sont transportées vers le pied aval avec l'écoulement sous le toit du matériau cohésif. Un conduit d'érosion étiré de l'amont vers l'aval se forme dans la fondation à travers ce processus, mettant en danger la stabilité du corps supérieur du barrage. Des études antérieures ont révélé que deux mécanismes différents sont responsables de la progression du conduit d'érosion, à savoir les érosions primaires et secondaires, c'est-à-dire l'érosion au niveau de la tête du conduit et du lit du conduit, respectivement. Plusieurs critères ont été développés sur la base de l'érosion secondaire, tels que la méthode de Sellmeijer pour déterminer le gradient global critique requis pour la progression du conduit d'érosion. Cependant, la plupart des investigations récentes se sont concentrées sur les conditions hydrauliques locales, car l'application de tels critères globaux sur le terrain pourrait être difficile en raison de l'hétérogénéité du terrain et de la stratification complexe du sol. Cette étude vise principalement à modéliser expérimentalement l'érosion régressive. Pour atteindre cet objectif, une cellule d'essai rectangulaire avec sortie de type trou est conçue pour contourner les inconvénients et les limitations des dispositifs existants. La configuration améliorée est composée de plusieurs rangées de capteurs de pression le long du chemin du conduit d'érosion, permettant de suivre librement les conditions hydrauliques le long des méandres du conduit. La compaction horizontale, l'application automatique de charge, le système d'acquisition de données de haute précision et fréquence et la procédure d'essai basée sur l'instrumentation sont d'autres améliorations ajoutées à la cellule. Les essais d'érosion avec la configuration améliorée sur différents sables fins uniformes montrent que le gradient global critique (i$\sf_{G\_cr}$) dans la configuration de sortie de type trou peut être expliqué par la conductivité hydraulique (k) et la taille de grain représentative (d$_{\mathsf70}$). La comparaison des données d'essai de la présente étude avec les données de la littérature démontre en outre que i$\sf_{G\_cr}$ corrèle parfaitement avec le facteur d'échelle (F$\sf_s$) qui intègre à la fois k et d$_{\mathsf70}$. En conséquence, le critère modifié de Sellmeijer est examiné pour son inclusion d'un nombre excessif de paramètres dans la prédiction du gradient global critique. La révision de la procédure de développement de la méthode de Sellmeijer révèle que le modèle est surajusté et inclut des paramètres intercorrélés. L'analyse statistique des données expérimentales utilisées pour modifier le critère de Sellmeijer montre qu'un modèle alternatif avec un nombre moindre de paramètres existe qui peut prédire le gradient global critique encore plus précisément. Le modèle proposé, développé sur la base des données expérimentales avec une sortie de type pente, est ensuite validé avec les données d'essai obtenues à partir de la configuration améliorée de cette étude. Il est démontré que le modèle est capable de prédire avec précision les deux types de configurations de sortie. Une analyse supplémentaire sur la capacité de prédiction du critère modifié de Sellmeijer rejette l'hypothèse précédente concernant sa surestimation 2 fois supérieure dans la prédiction du gradient pour les sorties de type trou. De plus, l'historique des charges et le gradient local sont étudiés dans les tests d'érosion sur les sables fins. Il est constaté que l'historique de charge hydraulique, c'est-à-dire l'application de charge en une seule étape ou en plusieurs étapes, n'a aucune influence sur le processus d'érosion. L'analyse des gradients locaux mesurés en amont du conduit d'érosion à l'aide des capteurs de la cellule reflète la relation mutuelle entre la progression du conduit d'érosion et le gradient local à la tête du conduit. Les résultats des essais indiquent la relation linéaire entre le gradient global critique et le gradient local. Une méthode géostatistique (krigeage) est utilisée pour étudier davantage l'influence de la distance entre les capteurs sur l'ampleur du gradient local mesuré. Il est démontré que le krigeage pourrait ne pas être un outil approprié pour l'interpolation de la charge hydraulique à proximité de la tête du conduit d'érosion. Une approche différente est utilisée, démontrant que le gradient local augmente de manière exponentielle avec la diminution de la distance entre les capteurs. / Backward erosion piping is one the initiating mechanisms of internal erosion that occurs in the foundations of embankment dams and dikes. Once the hydraulic load of the underground seepage reaches the critical conditions, particles of the granular layer are detached and transported towards the downstream toe with the seepage under the roof of the cohesive material. A stretched pipe from upstream to downstream forms in the foundation through this process which endangers the stability of the upper dam body. Previous studies found out that two different mechanisms are responsible for the progression of the pipe, namely primary and secondary erosions, i.e., erosion at the pipe tip and the pipe bed, respectively. Several criteria have been developed based on the secondary erosion such as Sellmeijer's method to determine the critical global gradient required for the pipe progression. Most recent investigations, however, focused on the local hydraulic conditions since application of such global criteria in the field could be challenging due to the field heterogeneity and complex soil stratification. This study aims primarily to model backward erosion experimentally. To achieve this end, a rectangular test cell with hole-type exit is designed which improves the drawbacks and limitations of the existing setups. The improved setup is composed of several rows of pressure sensors along the pipe pathway enabling to track freely meandering pipe's hydraulic conditions. Horizontal compaction, automatic load application, high accuracy and frequency data acquisition system and instrumentation-based testing procedure are other improvements added to the cell. The erosion tests with the improved setup on different fine uniform sands show that the critical global gradient (i$\sf_{G\_cr}$) in hole-type exit configuration can be explained by the hydraulic conductivity (k) and the representative grain size ($_{\mathsf70}$). Comparing the present study's test data with the data from the literature further demonstrates that i$\sf_{G\_cr}$ correlates perfectly with the scale factor (F$\sf_s$) which incorporates both k and $_{\mathsf70}$. Consequently, Sellmeijer's modified criterion is scrutinized for its inclusion of an excessive number of parameters in predicting the critical global gradient. Revisiting the development procedure of Sellmeijer's method reveals that the model is overfitted and it includes intercorrelated parameters. Statistical analysis on the experimental data used to modify Sellmeijer's criterion shows that an alternative model with a smaller number of parameters exists that can predict critical global gradient even more accurate. The proposed model developed based on the experimental data with slope-type exit is further validated with the test data obtained from this study's improved setup. It is shown that the model is capable of accurate prediction for both types of exit configurations. Additional analysis on the prediction capability of the Sellmeijer's modified criterion dismisses the prior hypothesis concerning its 2-time overestimation in prediction of the gradient for the hole-type exits. Furthermore, the loading history and local gradient are investigated in the erosion tests on the fine sands. It is found out that hydraulic loading history, i.e., single stage or multistage load application does not have any influence on the erosion process. The analysis of the local gradients measured upstream of the pipe using the cell sensors reflects the mutual relationship between the pipe progression and the local gradient at the pipe tip. The test results indicate the linear relationship between the critical global gradient and the local gradient. A geostatistical method (kriging) is used to further study the influence of the sensors distance on the magnitude of the measured local gradient. It is demonstrated that kriging might not be an appropriate tool to use for interpolation of the pressure head in the vicinity of the pipe tip. A different approach is used which demonstrates that local gradient increases exponentially with decreasing distance between the sensors.

La problématique de surverse du noyau de barrages en remblai lors de crues causées par une défaillance des systèmes de contrôle du niveau du réservoir préoccupent les gestionnaires de ces ouvrages. Une surverse de ce type peut entraîner des vitesses d’écoulement beaucoup plus importantes que celles de conception et causer un écoulement important dans des zones où aucun écoulement n’est initialement prévu. Pour évaluer l’impact des propriétés hydrauliques des composantes du barrage, de la géométrie du barrage et des caractéristiques de la crue, un modèle numérique 2D d’éléments finis a été conçu en utilisant le logiciel COMSOL Multiphysics 5.3. Le modèle numérique a permis de réaliser une série d’analyses de sensibilité sur un barrage modélisé typique qui a montré l’importance de la perméabilité de la crête et du noyau ainsi que de la hauteur maximale atteinte par la crue sur la vitesse à laquelle arrive la surverse et sur les vitesses d’érosion atteintes. L’érosion de contact selon le critère de Brauns (1985) ne semblait pas problématique lors des analyses de sensibilité alors que l’érosion interne selon le critère de Konrad & Côté (2013) paraissait très probable. Ensuite, une étude de cas a été réalisée portant sur un barrage existant d’Hydro-Québec. En raison des perméabilités très élevées de la crête pour ce barrage, les vitesses en crête dans le cas de hautes crues peuvent causer un écoulement turbulent à l’aval du noyau, sans toutefois être problématique pour cet ouvrage au niveau de l’érosion de contact. Finalement, la convergence numérique a été grandement meilleure lors de l’utilisation de courbes de rétention basée sur le modèle de van Genuchten (1980) comparativement à celles basées sur le modèle de Brooks & Corey (1964). / The issue of core overtopping in embankment dams during floods caused by failures of water-level control systems is a concern to dam managers. Core overtopping can generate flow velocities far greater than those considered during design and can result in intense flow in zones where no flow was intended. A 2D finite-element numeric model running on the COMSOL Multiphysics 5.3 software platform was generated to evaluate the impact of soil permeability, dam geometry and flood shape on seepage and internal flow systems. This numerical model was used to conduct a series of sensitivity analyses on a conceptual standard dam. These analyses have shown the importance of crest and core permeability and maximum height reached by the flood on the rate at wich core overtopping occurs and on erosion velocities generated. Contact erosion based on the Brauns (1985) criteria did not seem to be a concern for these sensitivity analyses as opposed to internal erosion based on the Konrad & Côté (2013) criteria which seemed highly probable. Subsequently, a case-study on an existing Hydro- Québec earthfill dam was also carried out. Because of the high permeability of the crest components of this structure, significant velocities and turbulent flow could be obtained in the crest material in cases of high floods, without however being problematic in terms of contact erosion. Furthermore, numerical convergence of the model has been greatly improved by using the van Genuchten-based (1980) water retention curves instead of the Brooks & Corey-based (1964) curves.

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.