Effect of ground-motion asynchronism on the equivalent acceleration of earth dams

Under seismic loads the deformability of an earth dam may induce several effects, including ground-motion amplification and asynchronism between different points of the dam embankment. The paper analyses the asynchronous effects occurring in two existing earth dams, representing well-documented case histories: the El Infiernillo Dam (Mexico) and the Camastra Dam (Italy). Asynchronous effects are analysed by theoretical predictions of the dam seismic response by adopting an advanced dynamic approach, which takes into account the main features that dam soils exhibit under cyclic loading conditions. For different potentially unstable masses within the dam embankment, equivalent accelerations were computed as the ratio between the resultant of the inertial forces and the weight of the volume V associated to the unstable mass. With the exception of very cortical sliding surfaces – not significant for dam stability – in most of the analysed cases the equivalent seismic coefficients do not exceed the peak acceleration at the dam base.     

Amplification of flood discharge caused by the cascading failure of landslide dams

The cascading failure of multiple landslide dams can trigger a larger peak flood discharge than that caused by a single dam failure.Therefore,for an accurate numerical simulation,it is essential to elucidate the primary factors affecting the peak discharge of the flood caused by a cascading failure,which is the purpose of the current study.First,flume experiments were done on the cascading failure of two landslide dams under different upstream dam heights,downstream dam heights,and initial downstream reservoir water volumes.Then,the experimental results were reproduced using a numerical simulation model representing landslide dam erosion resulting from overtopping flow.Finally,the factors influencing the peak flood discharge caused by the cascading failure were analyzed using the numerical simulation model.Experimental results indicated that the inflow discharge into the downstream dam at the time when the downstream dam height began to rapidly erode was the main factor responsible for a cascading failure generating a larger peak flood discharge than that generated by a single dam failure.Furthermore,the results of a sensitivity analysis suggested that the upstream and downstream dam heights,initial water volume in the reservoir of the downstream dam,upstream and downstream dam crest lengths,and distance between two dams were among the most important factors in predicting the flood discharge caused by the cascading failure of multiple landslide dams.     

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.     

Modeling of dam break wave propagation in a partially ice-covered channel

During the last four decades, several numerical formulations and specialized software have been developed in response to studies about dam break (DB) wave propagation and its hydraulic and environmental impacts on downstream hydraulic structures and valleys. These methods cannot, however, be used to predict wave propagation within partially covered channels or reservoirs located upstream of hydraulic structures. In fact, such problems require the modelling of the complex transition from a free surface flow into a pressurized one. Because rivers or channels partially covered with ice sheets are typical examples commonly met in winter in northern climates, it is vitally important to assess ice-cover effects on the DB wave propagation and develop appropriate tools to predict resulting hydrodynamic loads on hydraulic structures downstream. This paper proposes an original numerical formulation to model wave propagation and hydrodynamic pressure in partially covered channels. The proposed formulation uses one-dimensional St. Venant equations to simulate open-water flow and water hammer equations to simulate pressure flow within the partially covered channel. To illustrate the use of the hydrodynamic pressures obtained, a case study is presented where a channel cover and a dam located downstream are modelled using finite elements to investigate their dynamic structural response.     

Circular tunnels in sand: dynamic response and efficiency of seismic analysis methods at extreme lining flexibilities

The paper discusses the seismic response of circular tunnels in dry sand and investigates the efficiency of current seismic analysis methods at extreme lining flexibilities. Initially, a dynamic centrifuge test on a flexible circular model tunnel, embedded in dry sand, is analyzed by means of rigorous full dynamic analysis of the coupled soil–tunnel system, applying various non-linear soil and soil–tunnel interface models. The numerical results are compared to the experimental ones, aiming to better understand the recorded response and calibrate the numerical models. Then a series of numerical analyses are conducted using the validated numerical model, in order to investigate the effect of the tunnel lining rigidity on the dynamic response of the soil–tunnel system. In parallel, the accuracy of currently used simplified analysis methods is evaluated, by comparing their predictions with the results of the a priori more accurate and well validated numerical models. The comparative analyses allow us to highlight and discuss several crucial aspects of the soil-tunnel system seismic response, including (1) the post-earthquake residual values of the lining forces, which are amplified with the increase of the flexibility of the tunnel and (2) the importance of the soil-tunnel interface conditions. It is finally concluded that simplified analysis methods may provide a reasonable framework for the analysis at a preliminary stage, under certain conditions.     

Modeling the impact of dam removal on channel evolution and sediment delivery in a multiple dam setting

Dam removal can generate geomorphic disturbances, including channel bed and bank erosion and associated abrupt/pulsed release and downstream transfer of reservoir sediment, but the type and rate of geomorphic response often are hard to predict. The situation gets even more complex in systems which have been impacted by multiple dams and a long and complex engineering history. In previous studies one-dimensional (1-D) models were used to predict aspects of post-removal channel change. However, these models do not consider two-dimensional (2-D) effects of dam removal such as bank erosion processes and lateral migration. In the current study the impacts of multiple dams and their removal on channel evolution and sediment delivery were modeled by using a 2-D landscape evolution model (CAESAR-Lisflood) focusing on the following aspects: patterns, rates, and processes of geomorphic change and associated sediment delivery on annual to decadal timescales. The current modeling study revealed that geomorphic response to dam removal (i.e., channel evolution and associated rates of sediment delivery) in multiple dam settings is variable and complex in space and time. Complexity in geomorphic system response is related to differences in dam size, the proximity of upstream dams, related buffering effects and associated rates of upstream sediment supply, and emerging feedback processes as well as to the presence of channel stabilization measures. Modeled types and rates of geomorphic adjustment, using the 2-D landscape evolution model CAESAR-Lisflood, are similar to those reported in previous studies. Moreover, the use of a 2-D method showed some advantages compared to 1-D models, generating spatially varying patterns of erosion and deposition before and after dam removal that provide morphologies that are more readily comparable to field data as well as features like the lateral re-working of past reservoir deposits which further enables the maintenance of sediment delivery downstream.     

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.     

Channel adjustments of the lower Trinity River,Texas, downstream of Livingston Dam

Channel cross‐sectional changes since construction of Livingston Dam and Lake Livingston in 1968 were studied in the lower Trinity River, Texas, to test theoretical models of channel adjustment, and to determine controls on the spatial extent of channel response. High and average flows were not significantly modified by the dam, but sediment transport is greatly reduced. The study is treated as an opportunistic experiment to examine the effects of a reduction in sediment supply when discharge regime is unchanged. Channel scour is evident for about 60 km downstream, and the general phenomena of incision, widening, coarsening of channel sediment and a decrease in channel slope are successfully predicted, in a qualitative sense, by standard models of channel response. However, there is no consistent channel response within this reach, as various qualitatively different combinations of increases, decreases or no change in width, depth, slope and roughness occur. These multiple modes of adjustment are predicted by the unstable hydraulic geometry model. Between about 60 km and the Trinity delta 175 km downstream of the dam, no morphological response to the dam is observed. Rather than a diminution of the dam's effects on fluvial processes, this is due to a fundamental change in controls of the fluvial system. The downstream end of the scour zone corresponds to the upstream extent of channel response to Holocene sea level rise. Beyond 60 km downstream, the Trinity River is characterized by extensive sediment storage and reduced conveyance capacity, so that even after dam construction sediment supply still exceeds transport capacity. The channel bed of much of this reach is near or below sea level, so that sea level rise and backwater effects from the estuary are more important controls on the fluvial system than upstream inputs. Copyright © 2005 John Wiley & Sons, Ltd.     

Significance of small strain damping and dilation parameters in numerical modeling of free-field lateral spreading centrifuge tests

Lateral spreads are complex dynamic phenomena that are challenging to represent numerically. In this paper numerical models are developed and calibrated using the displacement, acceleration, and pore water pressure time histories recorded in a free-field lateral spreading centrifuge test. The calibrated numerical model then is used to predict another free-field lateral spreading centrifuge test using the same soil profile but different input acceleration time history. The computed response shows good agreement with the centrifuge test measurements. This paper demonstrates that even in a large strain problem, such as lateral spreading, small strain damping plays an important role in numerical simulation results; it also shows the need to have pressure dependent dilation parameters in the employed soil constitutive model implemented in order to simultaneously reproduce measurements of pore water pressure, acceleration and lateral displacement.     

Study of damage to a high concrete dam subjected to underwater shock waves

To study the effect of a strong underwater shock wave on a concrete dam, this research aims to improve hammer impact methods in model tests. Six 1:200 scale models were designed and tested under distributed impact loads. A device was deployed for a direct measurement of the impact force at the upstream face of the dams. The model dam bases were anchored to prevent displacement. The experimental results indicate that the top part of the concrete dam is a weak zone, and the impact failure initiates with a fracture on the top of the dam. The peak value of impact stress increases when the second crack appears in the concrete dam from the upstream face to the downstream face. And, the level of the second crack in the dam body is lower as the peak value of impact stress increases. In this study, dynamic analysis was conducted by calculating the results to verify the effectiveness of a device to directly measure the impact force. This method may be used to approximately forecast the damage of concrete dam and may also be useful in other engineering applications.     

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.     

A simplified model for lateral response of large diameter caisson foundations—Linear elastic formulation

The transient response of large embedded foundation elements of length-to-diameter aspect ratio D/B=2–6 is characterized by a complex stress distribution at the pier–soil interface that cannot be adequately represented by means of existing models for shallow foundations or flexible piles. On the other hand, while three-dimensional (3D) numerical solutions are feasible, they are infrequently employed in practice due to their associated cost and effort. Prompted by the scarcity of simplified models for design in current practice, we here develop an analytical model that accounts for the multitude of soil resistance mechanisms mobilized at their base and circumference, while retaining the advantages of simplified methodologies for the design of non-critical facilities. The characteristics of soil resistance mechanisms and corresponding complex spring functions are developed on the basis of finite element simulations, by equating the stiffness matrix terms and/or overall numerically computed response to the analytical expressions derived by means of the proposed Winkler model. Sensitivity analyses are performed for the optimization of the truncated numerical domain size, the optimal finite element size and the far-field dynamic boundary conditions to avoid spurious wave reflections. Numerical simulations of the transient system response to vertically propagating shear waves are next successfully compared to the analytically predicted response. Finally, the applicability of the method is assessed for soil profiles with depth-varying properties. The formulation of frequency-dependent complex spring functions including material damping is also described, while extension of the methodology to account for nonlinear soil behavior and soil–foundation interface separation is described in the conclusion and is being currently investigated.     

Elasto-plastic earthquake shear-response of one-dimensional earth dam models

A simplified analysis procedure for the non-linear hysteretic earthquake-response of earth dams is presented. The dam is modelled as a one-dimensional hysteretic shear-wedge subjected to base excitation. The hysteretic stress-strain behaviour of the dam materials is modelled by using elasto-plastic constitutive equations based on multi-surface kinematic plasticity theory. The method is based on a Galerkin formulation of the equations of motion in which the solution is expanded using eigenmodes of the linearized problem defined over the spatial domain occupied by the dam. The technique is applied to analyse the non-linear dynamic response of an earth dam subject to two very different input ground motions. The following investigations are presented: (i) comparison between the results obtained using two soil models depicting different nonlinear properties, (ii) comparison between the results of the one-mode and the multi-mode solution expansions, (iii) comparison with the results obtained through an elaborate finite element representation of the dam, and finally, (iv) comparison with the results obtained through the Makdisi-Seed¹¹ iterative procedure for earth dam analysis. The comparisons show that the proposed technique can be used to determine adequately the transient earthquake response of long earth dams. Furthermore, the efficiency and low computational cost make the technique very attractive; it can easily and systematically be extended to two- and three-dimensional calculations of earth dam response.     

Frequency analysis of annual maximum suspended sediment concentrations in Abiod wadi,Biskra (Algeria)

Sediment transport processes in the Mediterranean semi‐arid areas have interested a great number of hydrologists and statisticians. Frequency analysis (FA) procedures are commonly applied for several hydrological events such as floods and droughts. However, in general, FA is not widely applied to treat suspended sediment concentration (SSC), especially in semi‐arid regions. In the present study, an FA was performed on SSC data from 1979 to 1991 observed upstream from the Foum El Gherza dam, which is located at Biskra in the South–East of Algeria. This dam is problematic in terms of silting. Therefore, this study is mainly motivated by providing an FA model and risk evaluation in order to assist dam managers to better evaluate the potential of silting resulting from SSC transport to the reservoir. Probability distributions commonly used in hydrology were tested to SSC data recorded at the M'chounech station on Abiod wadi located upstream the dam. All the FA steps were considered; including classical techniques (e.g. goodness‐of‐fit tests, statistical criteria) as well as recently developed tools (tail distribution classification). Their application led to the selection of the lognormal distribution (LN2) to fit the considered data, and hence an accurate risk assessment could be obtained. Because of the presence of a slight trend, non‐stationary models are also considered in the present study. Copyright © 2013 John Wiley & Sons, Ltd.     

A new macromodel for the assessment of the seismic response of infilled RC frames

Numerous research studies have proved that numerical models aiming at an accurate evaluation of the seismic response of RC framed buildings cannot ignore the inelastic behaviour of infills and the interaction between infill and frame elements. To limit the high computational burden of refined non-linear finite element models, in the latest decades, many researchers have developed simplified infill models by means of single or multiple strut-elements. These models are low time-consuming and thus adequate for static and dynamic analyses of multi-storey structures. However, their simulation of the seismic response is sometimes unsatisfying, particularly in the presence of infill walls with regular or (particularly) irregular distributions of openings. This paper presents a new 2D plane macro-element, which provides a refined simulation of the non-linear cyclic response of infilled framed structures at the expense of a limited computational cost. The macro-element consists of an articulated quadrilateral panel, a single 1D diagonal link, and eight 2D links and is able to model the shear and flexural behaviour of the infill and the non-linear flexural/sliding interaction between infill and surrounding frame. The proposed macro-element has been implemented into the open source software OpenSees and used to simulate the response of single-storey, single-span RC infilled frame prototypes tested by other authors. The above prototypes are selected as made of different masonry units and characterised by full or open geometric configuration.     

Numerical investigation of the response of the Yele rockfill dam during the 2008 Wenchuan earthquake

In this paper the seismic response of a well-documented Chinese rockfill dam, Yele dam, is simulated and investigated employing the dynamic hydro-mechanically (HM) coupled finite element (FE) method. The objective of the study is to firstly validate the numerical model for static and dynamic analyses of rockfill dams against the unique monitoring data on the Yele dam recorded before and during the Wenchuan earthquake. The initial stress state of the dynamic analysis is reproduced by simulating the geological history of the dam foundation, the dam construction and the reservoir impounding. Subsequently, the predicted seismic response of the Yele dam is analysed, in terms of the deformed shape, crest settlements and acceleration distribution pattern, in order to understand its seismic behaviour, assess its seismic safety and provide indication for the application of any potential reinforcement measures. The results show that the predicted seismic deformation of the Yele dam is in agreement with field observations that suggested that the dam operated safely during the Wenchuan earthquake. Finally, parametric studies are conducted to explore the impact of two factors on the seismic response of rockfill dams, i.e. the permeability of materials comprising the dam body and the vertical ground motion.     

Effect of channel deepening on tidal flow and sediment transport: part I—sandy channels

Natural tidal channels often need deepening for navigation purposes (to facilitate larger vessels). Deepening often leads to tidal amplification, salinity intrusion, and increasing sand and mud import. These effects can be modelled and studied by using detailed 3D models. Reliable simplified models for a first quick evaluation are however lacking. This paper presents a simplified model for sand transport in prismatic and converging tidal channels. The simplified model is a local model neglecting horizontal sand transport gradients. The latter can be included by coupling (as post-processing) the simplified model to a 2DH or 3D flow model. Basic sand transport processes in stratified tidal flow are studied based on the typical example of the tidal Rotterdam Waterway in The Netherlands. The objective is to gain quantitative understanding of the effects of channel deepening on tidal penetration, salinity intrusion, tidal asymmetry, residual density-driven flow, and the net tide-integrated sand transport. We firstly study the most relevant tidal parameters at the mouth and along the channel with simple linear tidal models and numerical 2DH and 3D tidal models. We then present a simplified model describing the transport of sand (TSAND) in tidal channels. The TSAND model can be used to compute the variation of the depth-integrated suspended sand transport and total sand transport (incl. bed-load transport) over the tidal cycle. The model can either be used in stand-alone mode or with computed near-bed velocities from a 3D hydrodynamic model as input data.     

土层平均剪切波速确定及其可靠性验证——基于钻井台阵强震动记录

震害资料显示,场地条件对地震动特性以及工程结构破坏程度影响显著。为减少因场地效应而造成的经济损失和社会影响,在进行场地地震反应分析时,需最大限度地减小因场地土层模型参数的不确定性引起的地震动评估偏差,为工程结构地震反应分析选取并生成适当的地震动输入。随着强震动观测技术的逐渐发展,大量可靠的钻井台阵记录为地震过程中场地观测点的动力反应提供了直接数据。以美国加州地区La Cienega钻井台阵强震动观测数据为基础,利用互相关函数,对不同强度地震作用下场地土层的平均剪切波速进行分析,并在此基础上,以Cyclic 1D为模拟平台,建立一维自由场地地震反应有限元分析模型。分析结果表明：通过钻井台阵地震动观测数据识别,得到场地平均剪切波速,能够反映该场地的动力特性,数值模拟计算结果和台阵地震动记录基本吻合,可为数值模型参数选取提供依据。     

Stochastic dynamic response of dam–reservoir–foundation systems to spatially varying earthquake ground motions

In this paper, stochastic dynamic responses of dam–reservoir–foundation systems subjected to spatially varying earthquake ground motions are investigated using the displacement-based fluid finite elements. For this purpose, variable-number-node two-dimensional (2D) fluid finite elements based on the Lagrangian approach is programmed in FORTRAN language and incorporated into a computer program SVEM, which is used for stochastic dynamic analysis of solid systems subjected to spatially varying earthquake ground motion. The spatially varying earthquake ground motion model includes incoherence, wave-passage and site-response effects. The incoherence effect is examined by considering the Harichandran and Vanmarcke coherency model. The effect of the wave passage is investigated by using various wave velocities. Homogeneous medium and firm soil types are selected for considering the site-response effect where the foundation supports are constructed. The Sar?yar concrete gravity dam, constructed in Turkey is selected for numerical example. The ground motion is described by filtered white noise and applied to each support point of the 2D finite element model of the dam–reservoir–foundation system. The record of Kocaeli earthquake in 1999 is used in the analyses. Displacements, stresses and hydrodynamic pressures occurring on the upstream face of the dam are calculated for four cases. It is concluded that spatially varying earthquake ground motions have important effects on the stochastic dynamic response of dam–reservoir–foundation systems.     

Nonlinear dynamic analysis of Meloland Road Overpass using three-dimensional continuum modeling approach

This paper presents a three-dimensional (3D) continuum nonlinear analysis of the Meloland Road Overpass (MRO) near El Centro, California. The modeling methodology and the computational tools are discussed in detail. The performance of the computational model is evaluated by comparing the computed responses with the responses recorded at the bridge site during the 1979 Imperial Valley and 2010 El Mayor-Cucapah earthquakes. Amongst the recorded earthquake events at the bridge site, these two events caused the strongest shaking. The comparison shows that the 3D model is potentially an effective tool for detailed analysis of a full bridge system including foundation soils, pile foundations, embankments, supporting columns, and the bridge structure itself in a unified system without relying on any ancillary models such as Winkler springs. Additional response parameters such as displacements, rockings, and bending moments are also evaluated although none of these was measured during the seismic events.