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.     

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.     

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.     

SV波斜入射下混凝土重力坝易损性分析

易损性分析是评估不同强度地震作用下混凝土重力坝各级破坏概率的有效方法。目前重力坝易损性分析通常假定地震波为垂直入射,然而在近断层区域,地震波往往是倾斜入射的,地震波斜入射对重力坝地震响应有显著影响。从太平洋地震工程研究中心数据库选取16条地震动记录,采用黏弹性人工边界结合等效节点荷载实现SV波斜入射波动输入。采用增量动力分析方法对地震动峰值加速度进行调幅,以印度Koyna混凝土重力坝为研究对象,以坝顶相对位移为抗震性能指标,建立SV波斜入射下重力坝不同震损等级的易损性曲线。结果表明,与垂直入射相比,相同震损等级和相同地震动强度下,斜入射时重力坝破坏概率减小;当PGA接近重力坝实际遭受的地震动强度时,入射角为15°和30°时破坏概率与垂直入射相比最大减小率分别为27.3%和68.2%;各地震强度下,15°和30°斜入射相对于垂直入射的破坏概率差异值最大分别达36.6%、83.9%。因此,混凝土重力坝抗震性能分析应考虑地震波斜入射的影响。研究结果也可为近断层区域混凝土重力坝安全风险评估提供参考。     

A continuum damage concrete model for earthquake analysis of concrete gravity dam–reservoir systems

In this study, the earthquake damage response of the concrete gravity dams is investigated with considering the effects of dam–reservoir interaction. A continuum damage model which is a second-order tensor and includes the strain softening behavior is selected for the concrete material. The mesh-dependent hardening technique is adopted such that the fracture energy dissipated is not affected by the finite element mesh size. The dynamic equilibrium equations of motion are solved by using the improved form of the HHT-α time integration algorithm. Two dimensional seismic analysis of Koyna gravity dam is performed by using the 1967 Koyna earthquake records. The effects of damage on the earthquake response of concrete gravity dams are discussed. Comparison of the Westergaard and Lagrangian dam–reservoir interaction solutions is made. The effects of viscous damping ratio on the damage response of the dam are also studied.     

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.     

A comparative study of the different procedures for seismic cracking analysis of concrete dams

Different procedures are compared for the three-dimensional seismic cracking analysis of gravity and arch dams during strong earthquakes. The fracture procedures include the extended finite element method with cohesive constitutive relations, crack band finite element method with plastic-damage relations, and the finite element Drucker−Prager elasto-plastic model. These procedures are used to analyze the nonlinear dynamic response of Koyna dam to the 1967 Koyna earthquake and the seismic cracking of the Dagangshan arch dam subjected to design earthquake. The cracking process and profiles of the two dams using the three different procedures are compared. The applicability and the suitability of the three procedures for seismic cracking analysis of gravity and arch dams are discussed.     

Effects of near-fault and far-fault ground motions on nonlinear dynamic response and seismic damage of concrete gravity dams

As the forward directivity and fling effect characteristics of the near-fault ground motions, seismic response of structures in the near field of a rupturing fault can be significantly different from those observed in the far field. The unique characteristics of the near-fault ground motions can cause considerable damage during an earthquake. This paper presents results of a study aimed at evaluating the near-fault and far-fault ground motion effects on nonlinear dynamic response and seismic damage of concrete gravity dams including dam-reservoir-foundation interaction. For this purpose, 10 as-recorded earthquake records which display ground motions with an apparent velocity pulse are selected to represent the near-fault ground motion characteristics. The earthquake ground motions recorded at the same site from other events that the epicenter far away from the site are employed as the far-fault ground motions. The Koyna gravity dam, which is selected as a numerical application, is subjected to a set of as-recorded near-fault and far-fault strong ground motion records. The Concrete Damaged Plasticity (CDP) model including the strain hardening or softening behavior is employed in nonlinear analysis. Nonlinear dynamic response and seismic damage analyses of the selected concrete dam subjected to both near-fault and far-fault ground motions are performed. Both local and global damage indices are established as the response parameters. The results obtained from the analyses of the dam subjected to each fault effect are compared with each other. It is seen from the analysis results that the near-fault ground motions, which have significant influence on the dynamic response of dam–reservoir–foundation systems, have the potential to cause more severe damage to the dam body than far-fault ground motions.     

基于非平稳随机地震动模型的重力坝抗震可靠度分析

基于正交展开的非平稳随机地震动模型,并考虑混凝土材料的非线性和坝体与库水之间的流固耦合,对印度Koyna重力坝进行有限元分析,得到坝顶水平位移和坝颈拉应力,结合概率密度演化方法和等价极值事件的思想,获得丰富的概率信息。这为坝体结构的随机地震反应分析和可靠度研究提供新的途径。     

A plastic-damage concrete model for earthquake analysis of dams

A new plastic-damage constitutive model for cyclic loading of concrete has been developed for the earthquake analysis of concrete dams. The rate-independent model consistently includes the effects of strain softening, represented by separate damage variables for tension and compression. A simple scalar degradation model simulates the effects of damage on the elastic stiffness and the recovery of stiffness after cracks close. To simulate large crack opening displacements, the evolution of inelastic strain is stopped beyond a critical value for the tensile damage variable. Subsequent deformation can be recovered upon crack closing. The rate-independent plastic-damage model forms the backbone model for a rate-dependent viscoplastic extension. The rate-dependent regularization is necessary to obtain a unique and mesh objective numerical solution. Damping is represented as a linear viscoelastic behaviour proportional to the elastic stiffness including the degradation damage. The plastic-damage constitutive model is used to evaluate the response of Koyna dam in the 1967 Koyna earthquake. The analysis shows two localized cracks forming and then joining at the change in geometry of the upper part of the dam. The upper portion of the dam vibrates essentially as rigid-body rocking motion after the upper cracks form, but the dam remains stable. The vertical component of ground motion influences the post-cracking response. © 1998 John Wiley & Sons, Ltd.     

Seismic failure modeling of concrete dams considering heterogeneity of concrete

Study on the failure process of high concrete dams subjected to strong earthquakes is crucial to reasonable evaluation of their seismic safety. Numerical simulation in this aspect involves dynamic failure analysis of big bulk concrete dam subjected to cyclic loading. The Rock Failure Process Analysis (RFPA) proposed by C.A. Tang, with successful applications to failure modeling of rock and concrete specimens mainly subjected to static loading, is extended for this purpose. For using the proposed model, no knowledge on the cracking route needs to be known beforehand, and no remeshing is required. Simulation of the whole process of elastic deformation, initiation and propagation of microcracks, severe damage and ultimate failure of concrete dams in earthquakes with a unified model is enabled. The model is verified through a shaking table test of an arch dam. Finally a practical gravity dam is employed as a numerical example. Considering the uncertainty in ground motion input and concrete material, typical failure process and failure modes of gravity dam are presented. Several small cracks may occur due to tension particularly at dam neck, dam faces and dam heel, and a few of them evolve into dominant ones. Relatively smaller earthquake may cause damage to the dam neck while a bigger one may bring on cracks at lower parts of the dams. Cracking at the dam bottom may incline to a direction almost perpendicular to the downstream face after propagating horizontally for a certain distance when the shaking is strong enough.     

Seismic failure modeling of concrete dams considering heterogeneity of concrete

Study on the failure process of high concrete dams subjected to strong earthquakes is crucial to reasonable evaluation of their seismic safety. Numerical simulation in this aspect involves dynamic failure analysis of big bulk concrete dam subjected to cyclic loading. The Rock Failure Process Analysis (RFPA) proposed by C.A. Tang, with successful applications to failure modeling of rock and concrete specimens mainly subjected to static loading, is extended for this purpose. For using the proposed model, no knowledge on the cracking route needs to be known beforehand, and no remeshing is required. Simulation of the whole process of elastic deformation, initiation and propagation of microcracks, severe damage and ultimate failure of concrete dams in earthquakes with a unified model is enabled. The model is verified through a shaking table test of an arch dam. Finally a practical gravity dam is employed as a numerical example. Considering the uncertainty in ground motion input and concrete material, typical failure process and failure modes of gravity dam are presented. Several small cracks may occur due to tension particularly at dam neck, dam faces and dam heel, and a few of them evolve into dominant ones. Relatively smaller earthquake may cause damage to the dam neck while a bigger one may bring on cracks at lower parts of the dams. Cracking at the dam bottom may incline to a direction almost perpendicular to the downstream face after propagating horizontally for a certain distance when the shaking is strong enough.     

基于强震监测记录的冶勒大坝动力反应特性初步分析

冶勒沥青混凝土心墙堆石坝最大坝高为124.5m,坝址区地震烈度高,地质条件复杂。大坝上布设了9台强震仪组成的强震监测台阵,自2005年12月开始投入运行至2014年12月底共获得有效记录57次,其中包括汶川地震、攀枝花地震和芦山地震记录。对典型强震记录进行时域分析和频谱分析,初步总结了冶勒大坝的动力反应特性。研究结果表明,大坝坝顶部位的动力反应主要以放大为主,并且随着底部PGA的减小,坝体放大倍数明显增大;由于地震动本身震源特性、传播途径的差异,大坝在不同地震中动力反应的频谱特性明显不同,但多次记录均显示土石坝体有较为明显的滤波作用。     

Damage evaluation of concrete gravity dams under mainshock–aftershock seismic sequences

A large mainshock may trigger numerous aftershocks within a short period, and large aftershocks have the potential to cause additional cumulative damage to structures. This paper investigates the effects and potential of aftershocks on the accumulated damage of concrete gravity dams. For that purpose, 30 as-recorded mainshock–aftershock seismic sequences are considered in this study, and a typical two-dimensional gravity dam model subjected to the selected as-recorded seismic sequences is modeled. 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 the seismic damage process of Koyna dam under mainshock–aftershock seismic sequences. According to the characteristics of the cracking damage development, the local and global damage indices are both established to study the influence of strong aftershocks on the cumulative damage of concrete gravity dams. From the results of this investigation, it is found that the as-recorded sequences of ground motions have a significant effect on the accumulated damage and on the design of concrete gravity dams.     

Effects of reservoir bottom absorption on earthquake response of concrete gravity dams

A procedure is presented to analyse the response of concrete gravity dams due to horizontal and vertical earthquake ground motion components considering dam-water interaction and partial absorption of hydrodynamic pressure waves at the reservoir bottom into the foundation medium. The effects of reservoir bottom absorption on the hydrodynamic force on a rigid dam are examined first. The harmonic response of an idealized dam cross-section is presented for a wide range of parameters characterizing the properties of the dam, the impounded water and the foundation medium. Based on these frequency response functions the effects of dam-water interaction and of reservoir bottom absorption in the response of dams due to horizontal and vertical components of ground motion are investigated.     

Earthquake response of arch dams to spatially varying ground motion

The response of two arch dams to spatially varying ground motions recorded during earthquakes is computed by a recently developed linear analysis procedure, which includes dam–water–foundation rock interaction effects and recognizes the semi‐unbounded extent of the rock and impounded water domains. By comparing the computed and recorded responses, several issues that arise in analysis of arch dams are investigated. It is also demonstrated that spatial variations in ground motion, typically ignored in engineering practice, can have profound influence on the earthquake‐induced stresses in the dam. This influence obviously depends on the degree to which ground motion varies spatially along the dam–rock interface. Thus, for the same dam, this influence could differ from one earthquake to the next, depending on the epicenter location and the focal depth of the earthquake relative to the dam site. Copyright © 2009 John Wiley & Sons, Ltd.     

APPLICATION OF BOUNDARY ELEMENT ANALYSIS FOR MULTIPLE SEISMIC CRACKING IN CONCRETE GRAVITY DAMS

The previously developed two-dimensional boundary element procedure for analysing the propagation of a single discrete crack is extended to simultaneous multiple cracking in concrete gravity dams. A brief discussion of the generalized methodology is presented and the validity of the extended procedure is verified by performing a fracture analysis of the Fongman dam and comparing the predicted rupture process with the available experimental results. The fracture response of the Koyna dam is then studied extensively under the Koyna earthquake. Both single and multiple cracking models are employed to investigate the fracture process as well as final rupture in the dam. Similar final damage involving complete separation of the crest block of the dam is predicted, irrespective of whether single or multiple crack propagation models are employed. In relation to the phenomenon of hydrodynamic uplift pressure within propagating cracks, openings of the crack on the upstream face of the dam are examined in particular. The results indicate that this phenomenon is not expected to be significant during the crack development phase, and hence unlikely to affect the final rupture characteristics of dams undergoing strong earthquake excitation.     

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.     

Dynamics of submerged intake towers including interaction with dam and foundation

The dynamics of a coupled concrete gravity dam-intake tower–reservoir water–foundation rock system is numerically studied considering two hollow slender towers submerged in reservoir of gravity dam. The system is investigated in the frequency-domain using frequency response functions of the dam and the towers, and in the time-domain using time-history seismic analysis under a real earthquake ground motion. The analyzes are separately conducted under horizontal and vertical ground motions. The coupled system is three-dimensionally modeled using finite elements by Eulerian–Lagrangian approach. It is shown that presence of the dam significantly influences the dynamic response of the towers under both horizontal and vertical excitations; however the dam is not affected by the towers. When the dam is present in the model, the water contained inside the towers has different effects if the foundation is rigid, but it alleviates the towers motion if the foundation is flexible. It is concluded that the effects of foundation interaction are of much importance in the response of tall slender towers when they are located near concrete gravity dams.     

Seismic risk ranking for large dams in South Africa

Earthquake motion is one of the extreme loads acting on large dams. Dam owners and regulators must therefore ensure that dams are safely operated and present minimal risk to the public in case of extreme loads such as floods and earthquakes. Owners of many dams or officials in charge of dam safety programs may consider comparative assessment of the seismic risk associated with their dams and establish priorities for detailed evaluation. South Africa has in excess of 100 large state-owned dams and the characteristics of these dams have been used to perform a basic seismic hazard assessment and rank the vulnerability of these dams from the lowest to highest. One of the most decisive factors that contributes to the risk of a dam is the wall type; with gravity and earthfill dams being the most vulnerable to earthquake motion. Another aspect that needs further investigation is the downstream hazard potential which, if known to a better degree of accuracy, can provide more reasonable estimates of the risk factors.