Incremental dynamic analysis of concrete gravity dams including base and lift joints

The growth in computer processing power has made it possible to use time-consuming analysis methods such as incremental dynamic analysis(IDA) with higher accuracy in less time.In an IDA study,a series of earthquake records are applied to a structure at successively increasing intensity levels,which causes the structure to shift from the elastic state into the inelastic state and finally into collapse.In this way,the limit-states and capacity of a structure can be determined.In the present research,the IDA of a concrete gravity dam considering a nonlinear concrete behavior,and sliding planes within the dam body and at the dam-foundation interface,is performed.The influence of the friction angle and lift joint slope on the response parameters are investigated and the various limit-states of the dam are recognized.It is observed that by introducing a lift joint,the tensile damage can be avoided for the dam structure.The lift joint sliding is essentially independent of the base joint friction angle and the upper ligament over the inclined lift joint slides into the upstream direction in strong earthquakes.     

Two-dimensional dynamic analysis of concrete gravity and embankment dams including hydrodynamic effects

An analysis procedure in the frequency domain is developed for determining the earthquake response of two-dimensional concrete gravity and embankment dams including hydrodynamic effects; responses of the elastic dams and compressible water are assumed linear. The dam and fluid domain are treated as substructures and modelled with finite elements. The only geometric restriction is that an infinite fluid domain must maintain a constant depth beyond some point in the upstream direction. For such an infinite uniform region, a finite element discretization over the depth is combined with a continuum representation in the upstream direction. The fluid domain model approximately accounts for interaction between the fluid and underlying foundation medium through a damping boundary condition applied along the reservoir bottom, while the dam foundation is assumed rigid. Several examples are presented to demonstrate the accuracy of the fluid domain model and to illustrate dam responses obtained from the analysis procedure.     

Hydrodynamic effects in the dynamic response of concrete gravity dams

Hydrodynamic effects in the acceleration response of concrete gravity dams to harmonic ground motion are investigated. The effects include the presence of water, compressibility of water, interaction between the fluid and underlying foundation medium, shape of the fluid domain, and the extent of excitation applied to an infinite fluid domain under vertical ground motion.     

Seismic fracture analysis of concrete gravity dams

A numerical procedure for evaluation of the fracture process of gravity dams during strong earthquakes is presented. The BEM is used to discretize the dam reservoir system including the crack surfaces, and stress intensity factors at the crack tip are employed in a stage by stage procedure which simulates the crack extension. For each stage of constant crack length the mode superposition technique is applied; this is made possible by simulating the impact process of crack closing by a load pulse applied at the contact points which permits the structural stiffness to be assumed unchanged. To verify the proposed procedure, a cantilever beam model structure made of gypsum was tested on a shaking table. Good correlation with the numerical results was obtained, from which it is concluded that the procedure can be employed for evaluation of the crack propagation process in concrete structures subjected to dynamic loadings.     

Nonlinear dynamic response of concrete gravity dams: cavitation effect

A finite element method for the dynamic analysis of concrete gravity dams is presented. Displacement based formulation is used for both fluid and structural domains. During severe ground motion, the impounding fluid in the reservoir may separate from the dam and cause forming of micro bubbles. As a result, the compressibility of water is reduced. This nonlinear phenomenon of the reservoir is termed cavitation. When the direction of the ground motion is changed, the micro bubble's region of fluid collapses, and an impact will occur. By using different damping ratios in the fluid and solid domains the spurious oscillations which were caused by the impact are removed. The cavitation is confined to the upper part of the reservoir, where it has an effect of paramount importance on the tensile stresses. To illustrate the cavitation effect, the response of the non-overflow monolith of the Pine Flat dam subjected to the first 6.5 s of the May 1940 El-Centro, California earthquake, is considered. In order that the cavitation phenomenon take place more widely, maximum acceleration was scaled to give an amplitude of 1 g.     

Two‐dimensional modelling of ice cover effects for the dynamic analysis of concrete gravity dams

A two‐dimensional numerical model for determining the effects of the presence of an ice cover on the dynamic behaviour of large gravity dams is presented. Analytical predictions are compared to results obtained during a series of extensive dynamic tests on a large gravity dam. Data were obtained during summer and severe winter conditions to investigate the dynamic interactions between the dam, foundation, reservoir and the ice cover. The analysis includes ice‐reservoir interaction as well as the effects of water compressibility, flexible foundation and reservoir bottom absorption. Good agreement with the experimental findings is obtained. Copyright © 2002 John Wiley & Sons, Ltd.     

Earthquake analysis of concrete gravity dams including base sliding

A numerical method, the hybrid frequency-time domain (HFTD) procedure, is used to compute the earthquake response of concrete gravity dams, including sliding along the interface between the dam base and the foundation rock. The solution procedure accounts for the non-linear base sliding behaviour and the frequency-dependent response of the impounded water and flexible foundation rock. A Coulomb friction model represents the force-displacement relationship for sliding at the base interface. Using the solution procedure, an analysis of a typical dam (122 m high) shows that base sliding will occur during a moderate earthquake but the sliding displacement will be a tolerable amount when dam-foundation rock interaction is considered.     

Seismic stability of concrete gravity dams

Linear finite element analyses are commonly used to simulate the behaviour of gravity dam—foundation systems. However, the foundation is generally unable to develop any significant tensile stresses. Therefore any tension occurring in the vicinity of the dam—foundation interface is largely fictitious. Moreover, the traditional overturning and sliding stability criteria have little meaning in the context of the oscillatory response of dams during earthquakes. In this study, time domain analyses using non-linear contact elements located at the dam—foundation interface have been used to determine the dynamic sliding and uplifting response of gravity dam monoliths considering various elastic foundation properties. The magnitudes of the relative interface displacements, of the percentage of base not in contact (PBNC) and of the compressive stresses at the heel or toe of the dam have been used to monitor the seismic stability. The numerical results have shown that the non-linear behaviour of the dam—foundation interface reduces the seismic response of the system, indicating the possibility of more rational and economical designs. The PBNC was identified as the critical seismic stability response parameter for all analyses except for very flexible foundation conditions where the maximum values of relative interface displacements need to be considered.     

An application of damage mechanics for seismic analysis of concrete gravity dams

This paper discusses the local approach of fracture using damage mechanics concepts to evaluate the seismic response of concrete gravity dams. A constitutive model for plain concrete, subjected to tensile stresses, is presented. The mesh-dependent hardening technique is adopted such that the fracture energy dissipated is not affected by the finite element mesh size. The model is implemented in conjunction with the Hilber, Hughes Taylor alpha algorithm for time marching. Koyna dam is utilized to validate the proposed formulation. The importance of initial damage prior to the advent of an earthquake is also investigated. A 60 m concrete gravity dam is therefore selected and subjected to ground motion typical of eastern North America. Five scenarios of initial damage are presented and the results confirm the importance of accounting for the initial state for the seismic safety evaluation of an existing dam.     

Earthquake analysis of concrete gravity dams including dam-water-foundation rock interaction

A general procedure for analysis of the response of concrete gravity dams, including the dynamic effects of impounded water and flexible foundation rock, to the transverse (horizontal) and vertical components of earthquake ground motion is presented. The problem is reduced to one in two dimensions, considering the transverse vibration of a monolith of the dam. The system is analysed under the assumption of linear behaviour for the concrete, foundation rock and water. The complete system is considered as composed of three substructures—the dam, represented as a finite element system, the fluid domain, as a continuum of infinite length in the upstream direction, and the foundation rock region as a viscoelastic half-plane. The structural displacements of the dam are expressed as a linear combination of Ritz vectors, chosen as normal modes of an associated undamped dam-rock system. The effectiveness of this analytical formulation lies in its being able to produce excellent results by considering only a few Ritz vectors. The generalized displacements due to earthquake motion are computed by synthesizing their complex frequency responses using Fast Fourier Transform procedures. The stress responses are calculated from the displacements. An example analysis is presented to illustrate results obtained from this analytical procedure. Computation times for several analyses are presented to illustrate the effectiveness of the procedure.     

Seismic fragility assessment of concrete gravity dams

Many concrete gravity dams have been in service for over 50 years, and over this period important advances in the methodologies for evaluation of natural phenomena hazards have caused the design‐basis events for these dams to be revised upwards. Older existing dams may fail to meet revised safety criteria and structural rehabilitation to meet such criteria may be costly and difficult. Fragility assessment provides a tool for rational safety evaluation of existing facilities and decision‐making by using a probabilistic framework to model sources of uncertainty that may impact dam performance. This paper presents a methodology for developing fragilities of concrete gravity dams to assess their performance against seismic hazards. The methodology is illustrated using the Bluestone Dam on the New River in West Virginia, which was designed in the late 1930s. The seismic fragility assessment indicated that sliding along the dam–foundation interface is likely if the dam were to be subjected to an earthquake with a magnitude of the maximum credible earthquake (MCE) specified by the U.S. Army Corps of Engineers. Moreover, there will likely be tensile cracking at the neck of the dam at this level of seismic excitation. However, loss of control of the reservoir is unlikely. Copyright © 2003 John Wiley & Sons, Ltd.     

Seismic behaviour of cracked concrete gravity dams

A finite element model of incremental displacement constraint equations (IDCE), based on an existing node‐to‐surface concept, is implemented to deal with dynamic contact surfaces in the seismic behaviour analysis of cracked concrete gravity dams. After verification for sliding, rocking and impact, the IDCE model is applied to study the seismic responses of concrete gravity dams with different profiles and crack locations for a variety of parameters, such as coefficient of friction, water level and type of earthquake, as well as impact damping based on the concept of coefficient of restitution. It is revealed that cracked concrete gravity dams can experience not only sliding and rocking modes, but also the drifting mode in some cases of crack either at the base or at a height. Downstream sliding is normally accompanied by rocking, especially for the cases of crack at a height. Due to rocking and drifting, a cracked dam may still acquire a certain amount of residual sliding even if the effective coefficient of friction is relatively high. Copyright © 2005 John Wiley & Sons, Ltd.     

Investigation of the relationship of seismic intensity measures and the accumulation of damage on concrete gravity dams using incremental dynamic analysis

Nonlinear analysis tools are gaining prominence for the design and evaluation of concrete gravity dams. The performance limits of concrete gravity dams within the framework of performance based design are challenging to determine in comparison to those used for the assessments based on linear elastic analyses. The uncertainty in quantifying the behavior of these systems and the strong dependence of the behavior on the ground motion play an important role. The purpose of the study is to quantify the damage levels on a representative monolith using incremental dynamic analysis (IDA). For this purpose, the constitutive model utilized was calibrated first to the existing experimental results to verify the ability of the utilized cracking model to simulate the crack propagation process. Next, the relation between the damage levels on the monolith and the ground motion characteristics was investigated. The results of the conducted IDA showed that the engineering demand parameters (EDP) such as the crest displacement and acceleration showed weak correlation with the damage states. The spectral velocity and the peak ground acceleration were determined to be better predictors for the damage on the monolith. Copyright © 2015 John Wiley & Sons, Ltd.     

The effects of strong motion duration on the dynamic response and accumulated damage of concrete gravity dams

Strong motion duration is one of the challenging characteristics of ground motion, which affects the cumulative damage of structures significantly. Many researchers have conducted investigations related to the effects of strong motion duration on the response of building structures. However, the corresponding studies of concrete gravity dams are limited. In this paper, the effects of strong motion duration on the accumulated damage of concrete gravity dams are investigated. A Concrete Damaged Plasticity (CDP) model including the strain hardening or softening behavior is selected for the concrete material. This model is used to evaluate the nonlinear dynamic response and seismic damage process of Koyna dam during 1976 Koyna earthquake. Subsequently, the damage analyses of Koyna dam subjected to earthquake motions with different strong motion durations are performed. 20 as-recorded accelerograms, which are modified to match a 5% damped target spectrum, are considered in this study. Strong motion durations are obtained based on the definition of significant duration. According to the characteristics of the cracking damage development, both local and global damage indices are established as the response parameters. The results show that strong motion duration is positively correlated to the accumulated damage for events with similar response spectrum, and has significant influence on the cumulative damage of the dam. Longer duration will lead to greater accumulation damage to which aseismic design of the dam should be given attention.     

Earthquake safety assessment of concrete arch and gravity dams

Based on research studies currently being carried out at Dalian University of Technology, some important aspects for the earthquake safety assessment of concrete dams are reviewed and discussed. First, the rate-dependent behavior of concrete subjected to earthquake loading is examined, emphasizing the properties of concrete under cyclic and biaxial loading conditions. Second, a modified four-parameter Hsieh-Ting-Chen viscoplastic consistency model is developed to simulate the rate-dependent behavior of concrete. The earthquake response of a 278m high arch dam is analyzed, and the results show that the strain-rate effects become noticeable in the inelastic range. Third, a more accurate non-smooth Newton algorithm for the solution of three-dimensional frictional contact problems is developed to study the joint opening effects of arch dams during strong earthquakes. Such effects on two nearly 300m high arch dams have been studied. It was found that the canyon shape has great influence on the magnitude and distribution of the joint opening along the dam axis. Fourth, the scaled boundary finite element method presented by Song and Wolf is employed to study the dam-reservoir-foundation interaction effects of concrete dams. Particular emphases were placed on the variation of foundation stiffness and the anisotropic behavior of the foundation material on the dynamic response of concrete dams. Finally, nonlinear modeling of concrete to study the damage evolution of concrete dams during strong earthquakes is discussed. An elastic-damage mechanics approach for damage prediction of concrete gravity dams is described as an example. These findings are helpful in understanding the dynamic behavior of concrete dams and promoting the improvement of seismic safety assessment methods.     

Effects of concrete cracking on the earthquake response of gravity dams

Tensile stresses exceeding the tensile strength of concrete can develop in concrete dams subjected to earthquake ground motion. This study examines the earthquake response of gravity dams including tensile cracking of the concrete. The interaction between the dam and compressible water is included in the analysis using a numerical procedure for computing the non-linear dynamic response of fluid-structure systems. The crack band theory is used to model tensile cracking with modifications to allow for the large finite elements necessary for dam analysis. The earthquake response of a typical gravity dam monolith shows that concrete cracking is an important non-linear phenomenon. Cracking is concentrated near the base of the dam and near the discontinuities in the face slope. The extensive cracking, which develops due to ground motion typical of maximum credible earthquakes, may affect the stability of dams during and after strong earthquakes.     

Effects of sedimentary material on the response of concrete gravity dams

The transient response of a concrete gravity dam to a selected earthquake record is determined in order to evaluate the significance of the sedimentary material accumulating on the bottom of the reservoir. The sediment is modelled as a two-phase medium (fluid-filled, poroelastic solid), and the analysis is carried out accounting for all interactions between system components (dam, water, sediment, foundation).     

Seismic cracking and energy dissipation in concrete gravity dams

A finite element method for seismic fracture analysis of concrete gravity dams is presented. The proposed smeared crack analysis model is based on the non-linear fracture behaviour of concrete. The following features have been considered in the development of the model: (i) the strain softening of concrete due to microcracking; (ii) the rotation of the fracture band with the progressive evolution of microcrack damage in finite elements; (iii) the conservation of fracture energy; (iv) the strain-rate sensitivity of concrete fracture parameters; (v) the softening initiation criterion under biaxial loading conditions; (vi) the closing-reopening of cracks under cyclic loading conditions. The seismic fracture and energy response of dams and the significance of viscous damping models to take account of non-cracking structural energy dissipation mechanisms are discussed. The influences of global or local degradation of the material fracture resistance on the seismic cracking response of concrete dams were also studied. Two-dimensional seismic response analyses of Koyna Dam were performed to demonstrate the application of the proposed non-linear fracture mechanics model.     

Earthquake analysis of concrete gravity dams including reservoir bottom absorption and dam-water-foundation rock interaction

The available substructure method for the earthquake analysis of concrete gravity dams, including the dynamic effects of the impounded water and the flexible foundation rock, is extended to include the effects of alluvium and sediments invariably present at the bottom of actual reservoirs. Modelled approximately by a reservoir bottom that partially absorbs incident hydrodynamic pressure waves, these effects are incorporated into the continuum solution for the hydrodynamic pressure. The dam-water-foundation rock system is idealized as a two-dimensional system and analysed under the assumption of linear behaviour. An example earthquake analysis is presented to demonstrate the results obtained from the analytical procedure. Computation times for several cases illustrate the efficiency of the analytical procedure. In particular, the additional computation time required to include reservoir bottom absorption is shown to be very small.     

Seismic stability of concrete gravity dams strengthened by rockfill buttressing

Rockfill buttressing resting on the downstream face of masonry or concrete gravity dam is often considered as a strengthening method to improve the stability of existing dam for hydrostatic and seismic loads. Simplified methods for seismic stability analysis of composite concrete-rockfill dams are discussed. Numerical analyses are performed using a nonlinear rockfill model and nonlinear dam-rockfill interface behavior to investigate the effects of backfill on dynamic response of composite dams. A typical 35 m concrete gravity dam, strengthened by rockfill buttressing is considered. The results of analyses confirm that backfill can improve the seismic stability of gravity dams by exerting pressure on the dam in opposition to hydrostatic loads. According to numerical analyses results, the backfill pressures vary during earthquake base excitations and the inertia forces of the backfill are the main source for those variations. It is also shown that significant passive (or active) pressure cannot develop in composite dams with a finite backfill width. A simplified model is also proposed for dynamic analysis of composite dam by replacing the backfill with by a series of vertical cantilever shear beams connected to each other and to the dam by flexible links.