Nonlinear earthquake analysis of high arch dam–water–foundation rock systems

A nonlinear finite element model for earthquake response analysis of arch dam–water–foundation rock systems is proposed in this paper. The model includes dynamic dam–water and dam–foundation rock interactions, the opening of contraction joints, the radiation damping of semi‐unbounded foundation rock, the compressibility of impounded water, and the upstream energy propagating along the semi‐unbounded reservoir. Meanwhile, a new equivalent force scheme is suggested to achieve free‐field input in the model. The effects of the earthquake input mechanism, joint opening, water compressibility, and radiation damping on the earthquake response of the Ertan arch dam (240 m high) in China are investigated using the proposed model. The results show that these factors significantly affect the earthquake response of the Ertan arch dam. Such factors should therefore be considered in the earthquake response analysis and earthquake safety evaluation of high arch dams. Copyright © 2011 John Wiley & Sons, Ltd.     

Experimental and finite element studies of the forced vibration response of morrow point dam

Measured accelerations and water pressures obtained during a recent forced vibration test on a large thin arch dam at high water are compared to predictions from a finite element model for which water compressibility is both included and neglected. The numerical model is calibrated using the antisymmetric response data because they are only slightly affected by water compressibility; good agreement is obtained. In the effort to reproduce the symmetric response data, for which water compressibility plays a strong role, the calibrated model shows better correlation when water compressibility is included, but the agreement is still inadequate. A successful isolation of the fundamental water resonance from the experimental data shows significantly different features from those of the numerical water model, indicating possible inaccuracy in the assumed geometry and/or boundary conditions for the reservoir. Some other results at low water level are also included.     

Natural frequencies of arch dam reservoir system—by a mapping finite element method

A simple mapping finite element method is used to calculate the coupled natural frequencies and mode shapes of realistic arch dam reservoir systems in which the dam is circular cylindrical with non-uniform cross-section. This method, in which both the dam and the reservoir domains are mapped into geometrically simpler shapes using cylindrical-polar transformations, is found to give accurate results, achieved simply and economically. Results of analysis show that hydrodynamic interaction has a substantial effect on the coupled natural frequencies and mode shapes; also that the effect of water compressibility in the type of dams considered can be ignored without significant loss of accuracy. A simple method is also presented for predicting the water compressibility effect before undertaking detailed response analysis.     

Numerical simulation of earthquake excited dam-reservoirs with irregular geometries using an immersed boundary method

In this work, a ghost-cell immersed boundary method is proposed for the hydrodynamic response of earthquake excited dam-reservoirs. The numerical method employs a second order accurate two-step projection algorithm including compressibility effects in pressure field due to earthquake. The effects of reservoir bottom absorption are treated by introducing damping terms into the momentum equations. Hydrodynamic response of earthquake excited dam with a sloping face is simulated to demonstrate the accuracy of the present numerical method. Numerical results compared with previous numerical and analytical solutions show that the present immersed boundary method can accurately compute the hydrodynamic forces on inclined and curved dam faces including the effects of water compressibility and reservoir bottom absorption for the possibility of resonance. The proposed numerical method was shown to have significant advantages in computational time and memory usage for the hydrodynamic simulation of large dam-reservoirs with arbitrary geometries. Hydrodynamic forces on a double curvature arch dam subjected to real earthquake induced ground motion are also simulated to demonstrate the capability of the method.     

Study of the earthquake response of pine flat dam

The earthquake response of Pine Flat Dam is examined by a study of time history responses computed for a large set of earthquake ground acceleration records whose time axes have been systematically varied. Linear elastic behaviour is assumed. Topics considered include an investigation of the importance of the presence of water, water compressibility and the vertical component of ground motion; an evaluation of the accuracy of the lumped, added mass representation of the water; and a determination of the intensity of earthquake required to initiate non-linear behaviour in both the dam and water.     

Seismic response analysis of arch dam-water-rock foundation systems

The effect of water compressibility on the seismic responses of arch dams is not well understood. In this paper, a numerical model is developed with rigorous representation of the dynamic interaction between arch dam-water-rock foundation. The model is applied to the seismic response a nalysis of an arch dam with a height of 292m designed to aseismic intensity of IX. It is shown that consideration of the water compressibility clearly decreases the stress responses at key positions of the dam, while the added mass model gives a conservative estimate.     

Coupled hydrodynamic response of concrete gravity dams using finite and infinite elements

This paper presents the application of the finite element method for analysing the two-dimensional response of reservoir-dam systems subjected to horizontal ground motion. The interaction between the dam and the reservoir as well as the compressibility of water has been taken into account. The complete system has been considered to be composed of two substructures, namely the reservoir and the dam. To take into account the large extent of the reservoir, it has been idealized using specially developed infinite elements coupled with standard finite elements while the dam is represented using finite elements alone. Structural damping of the dam and radiation damping in the fluid phase have been accounted for in the analysis. It is concluded that the effect of radiation damping is considerable at high frequencies of excitation. The coupled response of the system is significantly large at and near the fundamental natural frequency of the system in comparison to the uncoupled responses. The method is computationally quite economical, capable of taking into account the arbitrary geometry of the system and is recommended for practical application. Further applications and extensions of the approach to three dimensional analyses are possible.     

Earthquake analysis of gravity dams including hydrodynamic interaction

A general procedure for analysis of the response of gravity dams, including hydrodynamic interaction and compressibility of water, 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 a dam, and the material behaviour is assumed to be linearly elastic The complete system is considered as composed of two substructures—the dam, represented as a finite element system, and the reservoir, as a continuum of infinite length in the upstream direction governed by the wave equation. The structural displacements of the dam (including effects of water) are expressed as a linear combination of the modes of vibration of the dam with the reservoir empty. The effectiveness of this analytical formulation lies in its being able to produce excellent results by considering only the first few modes. The complex frequency response for the modal displacements are obtained first. The responses to arbitrary ground motion are subsequently obtained with the aid of the Fast Fourier Transform algorithm An example analysis is presented to illustrate results obtained from this method. It is concluded that the method is very effective and efficient and is capable of producing results to any desired degree of accuracy by including the necessary number of modes of vibration of the dam.     

Seismic safety analysis of tall concrete dams,investigation and insights on critical challenges

This paper discusses critical and potentially controversial issues related to the seismic safety of tall concrete dams. These include the seismic input at a dam site, the effective treatment of the damage-rupture process, and the consideration of compressibility of reservoir water for hydrodynamic pressure. Major challenges to currently popular but questionable treatments of these critical problems are presented. Insights and additional research on these critical challenges are emphasized and explained based on prior published works of the author. More reasonable alternatives to dealing with these potentially controversial problems are provided in light of engineering practice in China. First, the design seismic input at depth as deconvoluted from an arbitrarily selected recorded accelerogram at a control point of an artificially developed free-field surface with the elevation of the dam crest is difficult for engineering projects to accept as appropriate. It may be more reasonable to use the design seismic incident motions as half of the ground surface motions from seismic safety analyses obtained from deterministic or probability approaches conducted by seismologists according to approved standards or guidelines. Second, since seismic damage to the dam must be estimated separately following uniaxial tensile and compressive experimental damage evolution rules, a simplified and realistic nonlinear elastic model is proposed as an alternative to the plastic-damage coupling model, which is very complex and includes assumptions based on a number of uncertainties. Finally, the effect of the reflection coefficient for compressibility of reservoir water on hydrodynamic pressures is very sensitive. The notion that the applied unified reflection coefficient at the reservoir bottom could be frequency-dependent and exhibit a significant variability in space as confirmed by field tests is questionable. To neglect the compressibility of reservoir water it may be closer to engineering practice at present.     

A solution procedure for the earthquake analysis of arch dam-reservoir systems with compressible water

Two processes using the Newmark implicit integration scheme are presented for the analysis of the earthquake response of a three-dimensional model for arch dam-reservoir systems including the effect of compressibility of the water. The solid structure and fluid regions are modelled separately, and the forcing functions at the interface are due to the hydrodynamic pressures from the reservoir acting on the upstream face of the dam wall, and the accelerations from the dam wall acting in turn on the reservoir. For the purposes of an initial investigation, elastic properties are assumed for the material of the dam, whilst in the reservoir radiation damping at the upstream boundary has been included, but bottom absorption has not. The excitation is provided by means of a combisweep which is fashioned so that its continuously varying frequencies pass through the fundamental frequencies of both the arch dam-reservoir system and the reservoir alone. Consequently the response is highly resonant, thus providing a severe test for the numerical procedures. From the numerical results obtained for an example problem it is concluded that both schemes provide an acceptable solution to the problem posed, and the possibility of enhancement to cater for more complex situations is discussed.     

Numerical simulation of reservoir sediment and effects on hydro‐dynamic response of arch dams

Based on the dynamic theory for saturated porous media by Biot (Journal of the Acoustical Society of America 1956; 28 : 168–178), a numerical model is presented to analyse the reflection behaviours of reservoir sediment and compared with those from the visco‐elastic model. It is concluded that the two models give very similar results of reflection coefficient α within the frequency range of interest. Then, using the two models, the change of the reflection coefficients α with various sedimentation parameters and excitation frequencies are studied in detail. The results are further used in the analysis of response functions of hydro‐dynamic pressures on, and structural displacements of the Xiang Hong Dian arch dam, for which some results from a field vibration test are available. It appears that effects of water compressibility with sediment reflection on hydro‐dynamic pressures and structural response are not significant for this specific case. Copyright © 2001 John Wiley & Sons, Ltd.     

Three‐dimensional models of reservoir sediment and effects on the seismic response of arch dams

The important effects of bottom sediments on the seismic response of arch dams are studied in this paper. To do so, a three‐dimensional boundary element model is used. It includes the water reservoir as a compressible fluid, the dam and unbounded foundation rock as viscoelastic solids, and the bottom sediment as a two‐phase poroelastic domain with dynamic behaviour described by Biot's equations. Dynamic interaction among all those regions, local topography and travelling wave effects are taken into account. The results obtained show the important influence of sediment compressibility and permeability on the seismic response. The former is associated with a general change of the system response whereas the permeability has a significant influence on damping at resonance peaks. The analysis is carried out in the frequency domain considering time harmonic excitation due to P and S plane waves. The time‐domain results obtained by using the Fourier transform for a given earthquake accelerogram are also shown. The possibility of using simplified models to represent the bottom sediment effects is discussed in the paper. Two alternative models for porous sediment are tested. Simplified models are shown to be able to reproduce the effects of porous sediments except for very high permeability values. Copyright © 2004 John Wiley & Sons, Ltd.     

Hydrodynamic pressures and response of gravity dams to vertical earthquake component

Hydrodynamic pressures and structural response of concrete gravity dams, including dam-reservoir interaction, due to the vertical component of earthquake ground motions are investigated. The response of the dam is approximated by the deformations in the fundamental mode of vibration, and the effects of deformability of bed rock on hydrodynamic pressures are recognized in the analysis. Expressions for the complex frequency response functions for the dam displacement, dam acceleration and lateral hydrodynamic force are derived. These results along with the Fast Fourier Transform algorithm are utilized to compute the time-history of responses of dams of 100, 300 and 600 ft height, with full reservoir, for different values of elastic modulus of mass concrete: 3.0, 3.5, 4.0, 4.5 and 5.0 million psi, to the vertical component of El Centro, 1940, and Taft, 1952, ground motions. It is concluded that the hydrodynamic forces caused by vertical ground motion are affected substantially by damreservoir interaction and depend strongly on the modulus of elasticity of the dam. The dam response to the vertical component of ground motion is compared with that due to the horizontal component. It is concluded that because the vertical component of ground motion causes significant hydrodynamic forces in the horizontal direction on a vertical upstream face, responses to the vertical component of ground motion are of special importance in analysis of concrete gravity dams subjected to earthquakes.     

库水模型及库水深度对拱坝动力响应的影响

基于接触非线性有限元模型,以锦屏一级拱坝为例,库水分别采用附加质量模型、可压缩流体有限元模型、不可压缩流体有限元模型计算了正常蓄水位及运行低水位时坝体的动力响应,结果表明:库水模型对拱坝动力响应有较大影响,随库水深度的增大,各模型计算结果差异增大;相比于流体可压缩模型,采用不可压缩流体模型所得动力响应普遍偏大;运行低水位工况,由于静水压力减小导致拱效应减弱,从而降低了拱坝的整体性,因此运行低水位工况各缝开度普遍高于正常蓄水位工况,且其拉应力范围较大,因此,运行低水位工况将对抗震设计起控制作用。     

An efficient approach for frequency-domain and time-domain hydrodynamic analysis of dam–reservoir systems

An efficient procedure is developed for the hydrodynamic analysis of dam–reservoir systems. The governing equations of hydrodynamic pressure in the frequency as well as time domain are derived in the framework of the scaled boundary finite element method. The water compressibility and absorption of reservoir sediments can be conveniently taken into consideration. By extending the reservoir to infinity with uniform cross-section, only the dam–reservoir interface needs to be discretized to model the fluid domain, and the hydrodynamic pressure in the stream direction is solved analytically. Several numerical examples including a gravity dam with an inclined upstream face and an arch dam with a reservoir of arbitrary cross-section are provided to demonstrate the computational efficiency and accuracy of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Hydrodynamic and foundation interaction effects in frequency response functions for concrete gravity dams

The linear response of idealized dam cross-sections to harmonic horizontal or vertical ground motion is presented for a range of the important system parameters characterizing the properties of the dam, foundation rock and impounded water. Based on these frequency response functions, the separate effects of interaction between the dam and water and interaction between the dam and foundation, and the combined effects of the two sources of interaction, on dynamic response of dams are investigated.     

Seismic response of intake towers including dam–tower interaction

The seismic response of the intake–outlet towers has been widely analyzed in recent years. The usual models consider the hydrodynamic effects produced by the surrounding water and the interior water, characterizing the dynamic response of the tower–water–foundation–soil system. As a result of these works, simplified added mass models have been developed. However, in all previous models, the surrounding water is assumed to be of uniform depth and to have infinite extension. Consequently, the considered added mass is associated with only the pressures created by the displacements of the tower itself. For a real system, the intake tower is usually located in proximity to the dam and the dam pressures may influence the equivalent added mass. The objective of this paper is to investigate how the response of the tower is affected by the presence of the dam. A coupled three‐dimensional boundary element‐finite element model in the frequency domain is employed to analyze the tower–dam–reservoir interaction problem. In all cases, the system response is assumed to be linear, and the effect of the internal fluid and the soil–structure interaction effects are not considered. The results suggest that unexpected resonance amplifications can occur due to changes in the added mass for the tower as a result of the tower–dam–reservoir interaction. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

Effects of reservoir bottom absorption and dam-water-foundation rock interaction on frequency response functions for concrete gravity dams

The linear response of an idealized concrete gravity dam monolith to harmonic horizontal or vertical ground motion is presented for a range of the important system parameters that characterize the properties of the dam, foundation rock, impounded water and reservoir bottom materials. Based on these frequency response functions, the effects of alluvium and sediments at the reservoir bottom on the response of the dam, including its interaction with the impounded water and foundation rock, are investigated. It is shown that the partial absorption of hydrodynamic pressure waves by the reservoir bottom materials has an important effect on the dynamic response of concrete gravity dams.