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.     

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.     

A unified approach for long-term behavior and seismic response of AAR-affected concrete dams

Concrete dams suffering from alkali-aggregate reaction (AAR) exhibit swelling and deterioration of concrete or even cracking over a long period. The deterioration of concrete may significantly affect the dynamic behavior of the structures, and it is necessary to estimate seismic safety of the deteriorated dams subjected to strong earthquakes. A unified approach is presented in this paper for long-term behavior and seismic response analysis of AAR-affected concrete dams by combining AAR kinetics, effects of creep and plastic-damage model in the finite element method. The proposed method is applied to a gravity dam and an arch dam. The long-term behavior of the AAR-affected dams is first predicted in terms of anisotropic swelling, spatially non-uniform deterioration of concrete, and cracking initiation and propagation with the development of AAR. The seismic response of the deteriorated dams is subsequently analyzed based on the state of the structures at the end of the long-term analysis. The AAR-induced expansion displacements obtained from the proposed method are in good agreement with the measured ones in the long-term operation. The simulated cracking patterns in the dams caused by the continuing AAR are also similar to the field observation. The results from the seismic analysis show that AAR-induced deterioration of concrete and cracking may lead to more severe damage cracking in the dams during earthquake. The dynamic displacements are also increased compared with the dams that are not suffering from AAR. The seismic safety of the AAR-affected concrete dams is significantly reduced because of the AAR-induced deterioration of concrete and cracking.     

Numerical simulation of reinforcement strengthening for high‐arch dams to resist strong earthquakes

This paper focuses on analyzing the nonlinear seismic response of high‐arch dams with cantilever reinforcement strengthening. A modified embedded‐steel model is presented to evaluate the effects of the strengthening measure on alleviating the extension and opening of cracks under strong earthquakes. By stiffening reinforced steel, this model can easily consider the steel–concrete interaction for lightly reinforced concrete (RC) members without the need of dividing them into RC and plain concrete zones. The new tensile constitutive relations of reinforced steel are derived from the load–deformation relationship of RC members in direct tension. This model has been implemented in the finite element code and its applicability is verified by two numerical simulations for RC tests. Subsequently, numerical analyses for a 210‐m high‐arch dam (Dagangshan arch dam) are conducted with and without the presence of cantilever reinforcement. Numerical results show that reinforcement strengthening can reduce the nonlinear response of the arch dam, e.g. joint opening and crest displacement, and limit the extension and opening width of concrete cracks. Copyright © 2008 John Wiley & Sons, Ltd.     

Analysis of seismic disaster failure mechanism and dam-break simulation of high arch dam

Based on a Chinese national high arch dam located in a meizoseismal region, a nonlinear numerical analysis model of the damage and failure process of a dam-foundation system is established by employing a 3-D deformable distinct element code(3DEC) and its re-development functions. The proposed analysis model considers the dam-foundation-reservoir coupling effect, infl uence of nonlinear contact in the opening and closing of the dam seam surface and abutment rock joints during strong earthquakes, and radiation damping of far fi eld energy dissipation according to the actual workability state of an arch dam. A safety assessment method and safety evaluation criteria is developed to better understand the arch dam system disaster process from local damage to ultimate failure. The dynamic characteristics, disaster mechanism, limit bearing capacity and the entire failure process of a high arch dam under a strong earthquake are then analyzed. Further, the seismic safety of the arch dam is evaluated according to the proposed evaluation criteria and safety assessment method. As a result, some useful conclusions are obtained for some aspects of the disaster mechanism and failure process of an arch dam. The analysis method and conclusions may be useful in engineering practice.     

Earthquake response of concrete arch dams: a plastic–damage approach

There are several alternatives to evaluate seismic damage‐cracking behavior of concrete arch dams, among which damage theory is the most popular. A more recent option introduced for this purpose is plastic–damage (PD) approach. In this study, a special finite element program coded in 3‐D space is developed on the basis of a well‐established PD model successfully applied to gravity dams in 2‐D plane stress state. The model originally proposed by Lee and Fenves in 1998 relies on isotropic damaged elasticity in combination with isotropic tensile and compressive plasticity to capture inelastic behaviors of concrete in cyclic or dynamic loadings. The present implementation is based on the rate‐dependent version of the model, including large crack opening/closing possibilities. Moreover, with utilizing the Hilber–Hughes–Taylor time integration scheme, an incremental–iterative solution strategy is detailed for the coupled dam–reservoir equations while the damage–dependent damping stress is included. The program is initially validated, and then, it is employed for the main analyses of the Koyna gravity dam in a 3‐D modeling as well as a typical concrete arch dam. The former is a major verification for the further examination on the arch dam. The application of the PD model to an arch dam is more challenging because the governing stress condition is multiaxial, causing shear damage to become more important than uniaxial states dominated in gravity dams. In fact, the softening and strength loss in compression for the damaged regions under multiaxial cyclic loadings affect its seismic safety. Copyright © 2013 John Wiley & Sons, Ltd.     

Seismic safety of high concrete dams

China is a country of high seismicity with many hydropower resources. Recently,a series of high arch dams have either been completed or are being constructed in seismic regions,of which most are concrete dams. The evaluation of seismic safety often becomes a critical problem in dam design. In this paper,a brief introduction to major progress in the research on seismic aspects of large concrete dams,conducted mainly at the Institute of Water Resources and Hydropower Research(IWHR) during the past 60 years,is presented. The dam site-specific ground motion input,improved response analysis,dynamic model test verification,field experiment investigations,dynamic behavior of dam concrete,and seismic monitoring and observation are described. Methods to prevent collapse of high concrete dams under maximum credible earthquakes are discussed.     

An investigation into the seismic behaviour of dams using dynamic centrifuge modelling

Dynamic response of dams is significantly influenced by foundation stiffness and dam-foundation interaction. This in turn, significantly effects the generation of hydrodynamic pressures on upstream face of a concrete dam due to inertia of reservoir water. This paper aims at investigating the dynamic response of dams on soil foundation using dynamic centrifuge modelling technique. From a series of centrifuge tests performed on model dams with varying stiffness and foundation conditions, significant co-relation was observed between the dynamic response of dams and the hydrodynamic pressures developed on their upstream faces. The vertical bearing pressures exerted by the concrete dam during shaking were measured using miniature earth pressure cells. These reveal the dynamic changes of earth pressures and changes in rocking behaviour of the concrete dam as the earthquake loading progresses. Pore water pressures were measured below the dam and in the free-field below the reservoir. Analysis of this data provides insights into the cyclic shear stresses and strains generated below concrete dams during earthquakes. In addition, the sliding and rocking movement of the dam and its settlement into the soil below are discussed.     

Behavior of large fill dams during earthquake and earthquake damage

This paper describes the behavior of and damage to large fill dams during earthquakes. Although fill dams have been constructed since olden times, there are only records of recent earthquake damage. We analyzed the behavior of fill dams of at least 15 m in height during earthquakes, and investigated the damage to the fill dams. The fact that fill dams that are designed and constructed based on the modern design standards have never been failed but have been only slightly affected by past earthquakes in Japan and other nations suggests that such large fill dams are highly earthquake resistant. The modern conventional design methods (soil engineering control of the bank soil, method for evaluating slip safety, and evaluation method of the foundation ground) are thus fully adequate for determining earthquake resistance.     

大坝抗震安全性研究进展

论证了广大坝抗震安全性研究的实践与发展现状。目前大坝在地震作用下的应力与变形分析方法主要有拟静力法和动力响应分析法,并依据大坝混凝土的抗拉强度判断大坝的安全性;各国规范体现的抗震设防弹念和大坝材料的容许应力差别很大。坝址河谷不同高程处地震动状态不尽相同、河谷两恻同一高程处地震动也不一样。混凝土材料的强度与加载速度、应变速率有关;地震时大坝不同部位的应变速率不相同、同一部位的应变速率也随时间变化;混凝土的动态强度既与应变速率有关。也与应变历史等其它因素有关。大坝河谷地震动的输入机理和模型研究、混凝土的动态强度的变化规律探索、大坝抗震安全性评价准则的完善与创新等将有待深入。通过以上内容针对性分析,提出了大坝抗震评价的一些合理建议、方法以及进一步的研究方向。     

汶川地震中宝珠寺重力坝地震响应的三维有限元模拟

汶川地震中宝珠寺水电站遭受的地震烈度为8°(相当于水平峰值加速度0.2 g),远超过大坝的设计地震水平(0.1 g),震后大坝未见明显震损.为解释大坝在地震中的抗震现象,构建了坝址区三维模型.考虑坝体横缝非线性以及三个方向地震作用的不同组合方式,对汶川地震中大坝的动力响应进行有限元模拟.在此基础上,针对震后提高的抗震设防标准,进一步选取典型坝段,采用二维弹塑性方法对大坝进行抗震复核并分析可能的破坏模式.模拟结果表明:横河向地震分量起主导作用而顺河向地震作用相对较弱是宝珠寺重力坝在汶川地震中免于发生损坏的主要原因.坝顶混凝土发生挤压破碎缘于永久横缝在地震中高频渐开渐合行为引起的剧烈碰撞.宝珠寺重力坝对设计地震0.27 g的强震可以保持整体的安全性,对校核地震0.32 g的强震整体安全性降低,水库正常运行及抵抗余震的能力将受到影响.     

Non‐linear seismic response of arch dams with contraction joint opening and joint reinforcements

A comprehensive study of non‐linear seismic response of arch dams with contraction joint opening and joint reinforcements has been conducted. A numerical model of contraction joint reinforcements is presented for optimization control of the joint opening. The objective of this control is to reduce the joint opening and expectantly to balance the sustained loads between the horizontal and the vertical components of the dam, thus avoiding an overstress in the cantilever while retaining the release of arch tensile stresses to some extent. Several parameter studies such as critical element size and required number of joints to be simulated for convergence are also performed. As an engineering application, a 292‐m high arch dam (the Xiaowan arch dam) and the Big Tujunga dam are analysed in detail. The results demonstrate that the joint opening and the corresponding load transfer from the arch to cantilever components of the dam during strong earthquakes are substantial. It is also evident that by providing sufficient strength and reinforcement flexibility, the joint opening can be controlled to some extent. However, the stress redistribution due to reinforcement control is not sufficient to avoid the overstress in the cantilever for the Xiaowan arch dam. Thus, alternative measures are discussed. Copyright © 2000 John Wiley & Sons, Ltd.     

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.     

Experimental study of seismic overloading of large arch dam

A dynamic overloading model test has been carried out on a shaking table for an arch dam of 278 m in height to investigate its behaviours under strong earthquake. The model system included the arch dam with contraction joints, part of reservoir, partial foundation rock with topographic feature near the dam. A damping boundary consisting of viscous liquid has been used to simulate the effect of dynamic energy emission to far field, which made the dynamic interaction between dam and foundation in model arch dam system be represented properly. Three sets of different seismic waves of design level have been used as the input to investigate the difference in the responses of arch dam. Artificial waves of different levels have been used to verify the behaviours of arch dam under seismic overloading. Since the opening of joints during strong earthquake reduced the response acceleration and tensile arch stress, cantilever stress on downstream face exceeded the tensile strength first for the model dam. And the arch dam responded in a non‐linear way when input seismic load increased. Some cracks appeared near abutments, and the damage made the natural frequency of arch dam to drop obviously, but the static function did not seem to change for the model tested. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

Parallel computation of seismic analysis of high arch dam

Parallel computation programs are developed for three-dimensional meso-mechanics analysis of fully-graded dam concrete and seismic response analysis of high arch dams （ADs）, based on the Parallel Finite Element Program Generator （PFEPG）. The computational algorithms of the numerical simulation of the meso-structure of concrete specimens were studied. Taking into account damage evolution, static preload, strain rate effect, and the heterogeneity of the meso-structure of dam concrete, the fracture processes of damage evolution and configuration of the cracks can be directly simulated. In the seismic response analysis of ADs, all the following factors are involved, such as the nonlinear contact due to the opening and slipping of the contraction joints, energy dispersion of the far-field foundation, dynamic interactions of the dam-foundation- reservoir system, and the combining effects of seismic action with all static loads. The correctness, reliability and efficiency of the two parallel computational programs are verified with practical illustrations.     

Nonlinear seismic response analysis of arch dam-foundation systems- part I dam-foundation rock interaction

The arch dam–foundation rock dynamic interaction and the nonlinear opening and closing effects of contact joints on arch dam are important to the seismic response analysis of arch dams. Up to date, there is not yet a reasonable and rigorous procedure including the two factors in seismic response analysis. The methods for the analysis of arch dam–foundation rock dynamic interaction in frequency domain are not suitable to the problem with nonlinear behaviors, in this paper, so an analysis method in time domain is proposed by combining the explicit finite element method and the transmitting boundary, and the dynamic relaxation technique is adopted to obtain the initial static response for dynamic analysis. Moreover, the influence of arch dam–foundation dynamic interaction with energy dispersion on seismic response of designed Xiaowan arch dam in China is studied by comparing the results of the proposed method and the conventional method with the massless foundation, and the local material nonlinear and nonhomogeneous behaviors of foundation rock are also considered. The reservoir water effect is assumed as Westergaard added mass model in calculation. The influence of the closing–opening effects of contact joints of arch dam on the seismic response will be studied in another paper.     

Seismic damage analysis of the outlet piers of arch dams using the finite element sub-model method

This study aims to analyze seismic damage of reinforced outlet piers of arch dams by the nonlinear finite element(FE) sub-model method. First, the dam–foundation system is modeled and analyzed, in which the effects of infinite foundation, contraction joints, and nonlinear concrete are taken into account. The detailed structures of the outlet pier are then simulated with a refined FE model in the sub-model analysis. In this way the damage mechanism of the plain(unreinforced) outlet pier is analyzed, and the effects of two reinforcement measures(i.e., post-tensioned anchor cables and reinforcing bar) on the dynamic damage to the outlet pier are investigated comprehensively. Results show that the plain pier is damaged severely by strong earthquakes while implementation of post-tensioned anchor cables strengthens the pier effectively. In addition, radiation damping strongly alleviates seismic damage to the piers.     

Integrative seismic safety evaluation of a high concrete arch dam

An integrative seismic safety evaluation of an arch dam should include all sources of nonlinearities, dynamic interactions between different components and the external loads. The present paper investigates the calibration procedure and nonlinear seismic response of an existing high arch dam. The first part explains the conducted analyses for the static and thermal calibrations of the dam based on site measurements. The second part investigates the nonlinear seismic analysis of the calibrated model considering the effect of joints, cracking of mass concrete, reservoir–dam–rock interaction, hydrodynamic pressure inside the opened joints and the geometric nonlinearity. Penetration of the water inside the opened joints accelerates the damage process. The integrative seismic assessment of a case study shows that the dam will fail under the maximum credible earthquake scenario. The dam is judged to be severely damaged with extensive cracking and the joints undergo opening/sliding. A systematic procedure is proposed for seismic and post-seismic safety of dams.