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.     

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.     

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.     

Effect of reservoir bottom on earthquake response of concrete dams

The absorption of hydrodynamic pressure waves at the reservoir bottom has dominant effects on the structural response of the dam when subjected to ground motion. In the present study, a model is proposed for the absorption effects of the reservoir bottom in the earthquake analysis of dams. The model utilizes the wave reflection coefficient approach and is based on the solution of the wave equation in a sediment layer of viscoelastic material with a constant thickness overlying an elastic, semi-infinite foundation. Numerical studies were conducted to evaluate the effect of the sediment layer thickness and material properties as well as the effect of reflection of waves from the underlying rock. It is shown that the current approach of assuming the wave reflection coefficient at the reservoir bottom based on the characteristics of the sediment material and excluding the effect of the reflected waves from the underlying rock, may significantly underestimate the seismic response of the dam.     

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.     

The effects of reservoir geometry on the seismic response of gravity dams

Conventional seismic analysis of gravity dams assumes that the behaviour of the dam–water–soil system can be represented using a 2‐D model since dam vertical contraction joints between blocks allow them to vibrate independently from each other. The 2‐D model assumes the reservoir to be infinite and of constant width, which is not true for certain types of reservoirs. In this paper, a boundary element method (BEM) model in the frequency domain is used to investigate the influence of the reservoir geometry on the hydrodynamic dam response. Important conceptual conclusions about the dam–reservoir system behaviour are obtained using this model. The results show that the reservoir shape influences the seismic response of the dam, making it necessary to account for 3‐D effects in order to obtain accurate results. In particular, the 3‐D pressure and displacement responses can be substantially larger than those computed with the 2‐D model. Copyright © 2007 John Wiley & Sons, Ltd.     

预应力混凝土结构非线性地震反应分析

本文针对预应力混凝土与普通混凝土不同的受力特点,提出了适合于分析预应力混凝土结构抗震性能的非线性有限元模型。模型首先将混凝土结构离散为杆单元,然后对各杆单元按分层组合原理分成许多混凝土层、预应力钢筋层和普通钢筋层,计算混凝土分层单元、预应力钢筋和普通钢筋的应力和应变,最后,根据钢筋与混凝土之间的粘结－滑移关系高速钢筋变形,根据混凝土弹性模量调整结构及杆单元变形,通过对普通混凝土构件和三个预应力混凝     

大开间小型混凝土砌块10层模型房屋抗震性能试验研究(Ⅰ)

我国抗震设计规范（GB50011－2001）规定在6、7、8度区，混凝土小砌块结构分别可以建七、六、五层。由于混凝土小砌块结构的最大优势在于10－20层（与混凝土框架或框架剪力墙结构比）。本文按1／4比例制作了10层混凝土小砌块结构模型，按7度设防要求实施构造措施，通过振动台试验研究模型结构的抗震性能。结果表明，模型结构完全能够满足在7度区“小震不坏，大震不倒”的要求。圈梁、构造柱以及水平拉结筋构成的约束体系抗震作用是明显的。试验利用砌块的非注芯孔灌注铁砂来模拟墙体出平面动力效应。在本模型的构造措施下，平出面反应不对结构破坏起控制作用。     

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.     

土石坝地震响应的非线性剪切条模型与对比分析

对一维剪切条计算模型进行改进,提出了土石坝非线性地震反应的简化计算方法。首先将坝体沿坝高离散为一系列的具有不同剪切刚度与阻尼比等参数特性的层状体系,建立了各层的振动控制方程及其边值条件,进而采用数学物理方程方法进行了求解,确定了体系的振动特性,并根据振型叠加原理和Duhamel积分确定了坝体地震反应的线弹性解。采用等价线性化方法考虑坝料的动力非线性性质,通过对线弹性地震响应的反复迭代计算,使得各层土的模量和阻尼比与其相应的剪应变水平相协调,确定出与非线性坝体系统相等效的线性解答,并将所得到的地震响应作为非线性地震响应的近似解。最后,以均质坝和心墙坝作为算例进行了具体的数值计算,将所得结果与有限元数值解进行对比分析,论证了所提方法的适用性和合理性。     

桥梁非线性地震反应分析若干问题研究现状

评述了钢筋混凝土桥墩非线性地震反应分析模型与损伤评估、桩－土－桥梁结构动力相互作用试验与理论研究和伸缩缝处的碰撞效应等问题的研究现状，指出了进一步的发展趋势。     

An alternative approach for computing seismic response with accidental eccentricity

Accidental eccentricity is a non-standard assumption for seismic design of tall buildings. Taking it into consideration requires reanalysis of seismic resistance, which requires either time consuming computation of natural vibration of eccentric structures or finding a static displacement solution by applying an approximated equivalent torsional moment for each eccentric case. This study proposes an alternative modal response spectrum analysis(MRSA) approach to calculate seismic responses with accidental eccentricity. The proposed approach, called the Rayleigh Ritz Projection-MRSA(RRP-MRSA), is developed based on MRSA and two strategies:(a) a RRP method to obtain a fast calculation of approximate modes of eccentric structures; and(b) an approach to assemble mass matrices of eccentric structures. The efficiency of RRP-MRSA is tested via engineering examples and compared with the standard MRSA(ST-MRSA) and one approximate method, i.e., the equivalent torsional moment hybrid MRSA(ETM-MRSA). Numerical results show that RRP-MRSA not only achieves almost the same precision as ST-MRSA, and is much better than ETM-MRSA, but is also more economical. Thus, RRP-MRSA can be in place of current accidental eccentricity computations in seismic design.     

Numerical analysis of permafrost effects on the seismic site response

Some of the damage to the infrastructure observed in past earthquakes occurred in Alaska could be related to the existence of permafrost. However, only limited research has been carried out so far to investigate the effects of permafrost on the seismic site response. Permafrost with relatively high shear wave velocity (1000–1500 m/s) extensively exists in the interior of Alaska and causes anomaly in the shear wave velocity profile that may alter the site response. In current design practices, permafrost has been treated as bedrock and its potential effects on site response are ignored. A systematic investigation was conducted to understand the effects of permafrost on the ground motion characteristics using one-dimensional equivalent linear analysis for the MCE, AASHTO and IBC Design Earthquake level hazards. The average surface displacement, velocity and acceleration response spectra for a typical permafrost site were obtained and the worst case scenario was identified. The results show that the presence of permafrost can significantly alter the ground motion characteristics and it may not be conservative to ignore the effects of permafrost in the seismic design of civil structures.     

基于地震记录的混凝土拱坝模态识别结果评价研究

利用混凝土拱坝地震记录识别的模态参数，可以揭示结构在地震过程中实际动力特性的变化情况，为结构地震反应分析和震后损伤评估提供重要信息。首先对模态识别常用方法的基本原理进行介绍；然后利用龙羊峡拱坝两次地震观测数据，采用"输入—输出"型和"仅考虑输出"型两类方法对大坝模态参数进行识别；最后将获得的模态识别结果与工程经验值和其他学者的研究成果进行对比，以分析评价识别结果的合理性。相关研究成果可为基于地震记录的混凝土拱坝模态参数识别的工程应用提供经验借鉴。     

A probabilistic approach for estimating the seismic response of elasto‐plastic SDOF systems

The seismic response of elasto‐plastic structures to both recorded and generated accelerograms is characterized by a large scattering of the results, even for accelerograms with similar peak ground acceleration values and frequency content. According to current code recommendations a design value of the seismic response of an elasto‐plastic structure can be computed as the mean of the responses to a certain number of spectrum‐fitting generated accelerograms. A more effective probabilistic approach is presented herein. It allows the analyst to calculate a design value of the seismic response characterized by a predefined non‐exceedance probability using a limited number of generated accelerograms. The results of the performed analyses are presented in diagrams that can be used for structural design applications. The applicability of the proposed method is demonstrated in the case of an elasto‐plastic structural system and the results are compared with those obtained applying current code recommendations. Copyright © 2005 John Wiley & Sons, Ltd.     

Assessment of seismic design response factors of concrete wall buildings

To verify the seismic design response factors of high-rise buildings, five reference structures, varying in height from 20- to 60-stories, were selected and designed according to modern design codes to represent a wide range of concrete wall structures. Verified fiber-based analytical models for inelastic simulation were developed, considering the geometric nonlinearity and material inelasticity of the structural members. The ground motion uncertainty was accounted for by employing 20 earthquake records representing two seismic scenarios, consistent with the latest understanding of the tectonic setting and seismicity of the selected reference region (UAE). A large number of Inelastic Pushover Analyses (IPAs) and Incremental Dynamic Collapse Analyses (IDCAs) were deployed for the reference structures to estimate the seismic design response factors. It is concluded that the factors adopted by the design code are adequately conservative. The results of this systematic assessment of seismic design response factors apply to a wide variety of contemporary concrete wall buildings with various characteristics.     

Evaluation of the seismic response of a homogeneous earth dam

The response of an earth dam to seismic loading is studied through displacement-based analyses and finite element, effective stress dynamic analyses. Displacement-based analyses are carried out using both empirical relationships and the decoupled approach in which the deformable response of the soil is accounted for through ground response analyses, and the resulting accelerograms are used in the sliding block analysis. The FE analyses are carried out using a constitutive model capable to reproduce soil non-linearity, calibrated against laboratory measurements of the stiffness at small strains. The influence of the assumed input motion and bedrock depth on the seismic response of the dam is also studied.     

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.     

基于超级单元的钢-混凝土混合结构动力特性及地震反应研究

本文提出了基于超级单元的钢框架－混凝土剪力墙结构体系动力特性及地震反应分析的一种方法，该方法考虑了混合结构体系中钢框架部分的受力和变形特点，即节点柔性（半刚性）、梁柱效应（p-δ效应）及P－△效应（以下合称P－Delta效应），推导了考虑三种因素情况下，钢框架的层间侧移刚度及钢框架－混凝土剪力墙结构体系的超级单元刚度矩阵，采用刚度法求解钢框架－混凝土混合结构周期和振型。该方法物理概念清楚，计算简便。     

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.