土石坝极限抗震能力的上限有限元法

基于极限分析上限定理,考虑堆石料内摩擦角较大以及抗剪强度具有非线性的特性,提出一个用于求解土石坝极限抗震能力的有限元计算方法。该方法通过功能平衡条件、位移协调条件、边界条件、屈服条件以及所求极限荷载,形成标准的二阶锥规划数学模型,并用内点法进行优化迭代求解,得到土石坝坡极限抗震能力的上限值。运用该方法对一个典型面板堆石坝进行抗震极限分析。结果表明:按规范设计的土石坝具有较强的抗震能力,且竖向地震力对大坝的极限抗震能力存在一定影响;此外上限有限元法具有较高的计算精度及工程实用价值。     

考虑堆石料空间变异性的土石坝坝坡地震稳定性随机有限元分析

本文探讨了筑坝堆石料的空间变异性对土石坝坝坡动力稳定性的影响。以新疆某在建高面板堆石坝为例,在蒙特卡洛法的框架下,采用基于局部平均细分法的随机有限元法模拟考虑筑坝堆石料空间变异性时土石坝的地震响应及坝坡滑移情况,通过对比随机有限元法和常规确定性有限元法的计算结果,提出:在地震动作用下,考虑筑坝材料空间变异性时,坝坡滑动体的数量、规模以及滑移量和滑动历时都有不同程度的增大,因而坝坡整体危险程度显著高于不考虑材料空间变异性的情况。坝坡各项动力安全性指标对筑坝材料空间变异性非常敏感;因而,考虑筑坝材料空间变异性时,各项安全性指标的离散性较大。     

Seismic failure probability of concrete slab on CFR dams with welded and friction contacts by response surface method

In this study, failure probability of the concrete slab on concrete-faced rockfill (CFR) dams with welded and friction contact is investigated under earthquake effects by reliability analysis. For this purpose, Torul CFR dam is selected as an example and numerical solutions are performed by considering combination of reliability analysis–finite element method. 1992 Erzincan earthquake acceleration record is used in the finite element analysis considering deconvolved-base rock input model. In this model, the ground motion to be applied to the foundation base rock is obtained by deconvolution of the free-field surface record. In the materially nonlinear analysis, Drucker–Prager model is used for concrete slab and multi-linear kinematic hardening model is utilized for rockfill. Geometrically nonlinearity is also taken into account. Viscous boundary conditions are defined in the finite element model for both foundation soil and reservoir water. The hydrodynamic pressure of the reservoir water is considered using 2D fluid finite elements based on the Lagrangian approach. Both welded contact and friction contact based on the Coulomb’s friction law are defined in the structural connections. Improved Rackwitz–Fiessler method is used with response surface method in the reliability analysis. The tensile and compression strengths of the concrete slab are utilized in the implicit limit state functions considering various thicknesses. The probability of failure of the most critical points in the concrete slab is obtained. According to this study, the probabilities of failure obtained from the CFR dam including friction contact are lower. When the welded contact is considered in joints, the probability of failure of the concrete slab is 1 due to tensile stress limit state and compression stress limit state only if concrete slab is linear. The most critical probability of failure of the concrete slab appears in the case that the concrete slab is linear and rockfill is materially nonlinear. The probability of failure of the concrete slab decreases if the nonlinearity of the concrete is considered. Also, hydrodynamic pressure decreases the reliability of the concrete slab.     

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.     

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.     

Linear and nonlinear response of concrete slab on CFR dam during earthquake

The joint between concrete slab and rockfill is designed as welded contact in the classical modeling of concrete-faced rockfill (CFR) dams and earthquake response of the CFR dams is determined by this method. In this study, linear and nonlinear response of Torul CFR Dam including interface element between concrete slab and rockfill were investigated for the duration of strong seismic excitation. The finite element analyses were performed by employing both cases, empty and full reservoir, to research the effect of the reservoir water on the earthquake response of the dam. The reservoir water was modeled with fluid finite elements by the Lagrangian approach. The Drucker-Prager model was used in nonlinear analyses for concrete slab, rockfill and soil materials. According to finite element analyses, displacement and stress components were increased by hydrodynamic pressure. The nonlinear response of the concrete slab was monitored about the peak ground acceleration (pga). This study reveals that the size of sliding zone increases with increasing acceleration amplitudes.     

地震波斜入射对高面板坝地震反应的影响

为研究地震波斜入射对高面板坝地震反应的影响,根据地震波动入射理论,采用FORTRAN进行波动荷载的编程计算,并与大型通用有限元软件ADINA相结合,实现基于黏弹性人工边界的地震波斜入射,研究P波和SV波分别以不同角度入射对高面板堆石坝地震反应的影响。结果表明,地震波斜入射时大坝地震动反应与垂直入射时明显不同,常规垂直入射的结果偏于不安全,因此在高面板坝地震反应分析和抗震设计中应考虑地震波斜入射的影响。     

土石坝地震稳定性分析方法简介

土石坝地震稳定性分析方法是随着科技的进步，特别是计算机技术的快速发展以及对大量土石坝震害现象进行调查研究的不断深入而发展起来的。本文详细介绍了几种典型分析方法的原理、步骤以及优缺点，并对这些分析方法进行了评述，最后总结了近年来土石坝地层稳定性分析方法发展概况，指出了今后的发展趋势。     

高土石坝震害与抗震措施评述

目前,我国西部地区正建或拟建一大批高土石坝,这些高坝都位于地震高烈度区,做好抗震设计非常重要。本文总结了地震荷载作用下土石坝常见的几种破坏形式,包括地震永久变形、滑坡失稳、液化、防渗体破坏、次生破坏;针对这些破坏形式,总结了地震高烈度区高土石坝的抗震措施,包括预留超高、坡脚压重、坝顶加固、抗液化措施、防渗体处理等;提出了现阶段高土石坝抗震的计算方法和抗震措施存在的缺点和不足,并对我国今后的高土石坝抗震稳定研究提出了建议。     

Seismic stability analysis of fractured rock slopes by yield design theory

This paper deals with the problem of stability of fractured rock slope located in seismic area. The rock mass is crossed by two sets of fractures which are considered to be planar, parallel and persistent. The effects of both horizontal seismic coefficient and strength characteristics of fractures are addressed. The analysis is based upon the kinematic approach of the yield design theory and the pseudo-static method as well. The fundamental inequality of the kinematic approach is invoked and the failure of the fractured rock slope is considered through simple translational mechanism involving a one-partsliding block. Rigorous upper bounds of the so-called stability factor for the structure under study of given slope angle, strength properties of constituent materials and seismic loads are obtained. The results are presented in the form of stability charts relating the estimated upper bound solutions to the friction angle of fractures.The used procedure highlights the destabilizing effects of the seismic loadings and is capable of conducting parametric studies. Furthermore, the results obtained far from various sets of data show good performance and further research work is planed to extend the analysis to include a large number of sub-problems.     

Seismic cracking and energy dissipation in concrete gravity dams

A finite element method for seismic fracture analysis of concrete gravity dams is presented. The proposed smeared crack analysis model is based on the non-linear fracture behaviour of concrete. The following features have been considered in the development of the model: (i) the strain softening of concrete due to microcracking; (ii) the rotation of the fracture band with the progressive evolution of microcrack damage in finite elements; (iii) the conservation of fracture energy; (iv) the strain-rate sensitivity of concrete fracture parameters; (v) the softening initiation criterion under biaxial loading conditions; (vi) the closing-reopening of cracks under cyclic loading conditions. The seismic fracture and energy response of dams and the significance of viscous damping models to take account of non-cracking structural energy dissipation mechanisms are discussed. The influences of global or local degradation of the material fracture resistance on the seismic cracking response of concrete dams were also studied. Two-dimensional seismic response analyses of Koyna Dam were performed to demonstrate the application of the proposed non-linear fracture mechanics model.     

Three-dimensional nonlinear seismic analysis of concrete faced rockfill dams subjected to scattered P,SV, and SH waves considering the dam–foundation interaction effects

In this study, the nonlinear seismic analysis of a typical three-dimensional concrete faced rockfill dam is reported. Three components of the Loma Prieta (Gilroy 1 station) earthquake acceleration time history are used as input excitation. The dam under study is considered as if it were located in a prismatic canyon with a trapezoidal cross-section. A nonlinear model for the rockfill material is used, and contact elements with Coulomb friction law are utilized at the slab–rockfill interface. Vertical joints in the face slab are also considered in the finite element model. A substructure method, in which the unbounded soil is modelled by the scaled boundary finite element method (SBFEM), is used to obtain the scattered motion and interaction forces along the canyon. The dam is subjected to spatially variable P, SV, and SH waves, and the effect of dam–foundation interaction and the reservoir water effects are considered. The results are compared with the non-scattered input motion analysis. Results of the analyses indicate that due to applying the scattered motion to the canyon the response of the dam and concrete face slab significantly increases. The reservoir water pressure affects the tensile stresses induced in the face slab by reducing the uplift movement of the concrete panels.Large horizontal axial forces are induced in the face slab due to out-of-phase and out-of-plane motions of the abutments. Although the normal movements of vertical joints are reduced due to the reservoir water confinement, the opening movements are still significant, and the local failure of construction joints is inevitable.     

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.     

汶川大地震中德阳市水库土坝的震害特征

汶川特大地震给四川省德阳市水库土坝造成了不同程度的震损。依据德阳市境内66座高危以上险情水库土坝震害考察资料和现场震害调查结果,对震损水库土坝进行了分类统计,研究其震害特征及震损原因,结果表明：小型水库均质土坝在此次地震中震损数量最多、险情严重,且震害类型丰富;典型震害包括裂缝、滑坡、沉陷、渗漏以及泄水建筑物与其它附属设施损坏。以德阳市3座典型震损水库土坝为例,分析了水库土坝的基本情况、震前运行状况及震损情况,进而提出了有效的震后抢险措施。     

Seismic dynamic response by approximate methods

Current seismic safety evaluation for earth dams relies on approximate methods of analysis for prediction of non-linear, transient, dynamic response. One of these approximate methods uses a Strain Reduction Factor (SRF), and has been widely applied in a variety of one-dimensional and two-dimensional soil structural analyses. A second method considered, Equivalent Temporal Damping (ETD), has been previously applied to several seismic dynamic analyses of earth dams. The relative accuracy of the two methods is assessed by comparing them with incremental plasticity, nearly exact numerical solutions. For a simple shear element subjected to deterministic and non-stationary random input accelerations, a serious overprediction of the maximum peak shear response occurs by the SRF method, whereas the ETD results agree very closely with the incremental plasticity solutions. The SRF and ETD methods are also applied to seismic dynamic response analysis of a Vertical Soil Column system, and the same trends as established in the previous case are observed. It is concluded that, in lieu of combined incremental plasticity and finite difference or finite element numerical solutions, the ETD approach is the more accurate of the methods tested for seismic dynamic response analysis of earth dams.     

Longitudinal vibrations of tall concrete faced rockfill dams in narrow canyons

Longitudinal vibrations of concrete-faced rockfill dams may cause significant compressive stresses and joint openings in the slab panels. The behavior of such dams subjected to longitudinal and vertical vibrations is investigated, based on numerical simulation of the staged construction, reservoir impoundment and seismic shaking. The static analysis uses a hyperbolic model for the rockfill, whereas the seismic analysis uses a nonlinear hysteretic model which accounts for the dynamic properties. A damage plasticity model is used for the reinforced concrete and frictional contact behavior is considered at the base and vertical sides of the concrete panels. The seismic analysis takes into account the flexibility of the canyon rock and potential dynamic rockfill settlements. An existing 150 m–high dam is used to investigate the effect of longitudinal vibrations on the compressive stresses near the slab-to-slab vertical interfaces and the opening of the joints. The effect of dynamic settlements is examined and comparisons are made to the response from upstream/downstream and combined vibrations. The effectiveness of introducing 5 cm-wide cuts at selected vertical joints to reduce slab compression in existing CFRDs is demonstrated. The presented results offer useful insight into the effect of longitudinal vibrations on the seismic behavior of CFRDs.     

Seismic rehabilitation and analysis of Chaohe earth dam

Stability of earth dams during earthquakes has been a major concern for geotechnical engineers in seismic active regions. Liquefaction induced slope failure occurred at the upstream slope of a major earth dam in the suburb of Beijing, China, during the 1976 Tangshan Earthquake. The gravelly soil with loose initial condition liquefied under relatively small ground vibration. In recent years, a major seismic rehabilitation project was carried out on a similar earth dam nearby using dumped quarry stone. Seismic stability analysis was carried out using model test, finite element simulation, and pseudo-static slope stability program after taking into account the influence of excess pore pressure.     

Nonlinear seismic response of tall concrete-faced rockfill dams in narrow canyons

The seismic behavior of tall concrete face rockfill dams in narrow canyons is investigated, based on numerical simulation of the staged construction, creep settlements, reservoir impoundment and seismic shaking of the dam. The study takes into account the flexibility of the canyon rock, the hydrodynamic effects and potential dynamic rockfill settlements. The static analysis uses a hyperbolic model for the rockfill, whereas the dynamic analysis uses a nonlinear hysteretic model, which accounts for the initial dynamic stiffness and produces hysteresis loops in agreement with the experimental data regarding the shear modulus and damping ratio. A damage plasticity model is used for the reinforced concrete, whereas frictional contact behavior is considered at the base and vertical walls of the concrete slab panels. An existing 150-m-high dam is used to investigate some key issues on the seismic behavior of such dams subjected to upstream−downstream and vertical excitation. Emphasis is placed on the evaluation of the tensile stresses within the slab panels, the compressive stresses at the slab-to-slab vertical interfaces and the opening of the joints. Moreover, the effect of potential dynamic settlements on both the slab stresses and joint openings is investigated. Recommendations for increasing the dam safety and reducing the water leakage through the dam body are given.     

基于推覆分析的钢框架结构地震倒塌判别准则研究

结构地震倒塌判别准则是工程结构强震分析的关键问题。在层损伤模型的基础上,建立了基于推覆分析的建筑结构整体损伤模型,并以国内某2层2跨平面钢框架结构拟静力试验为背景,应用有限元程序ABAQUS对平面钢框架进行了强震倒塌数值模拟。分析了钢框架结构的倒塌破坏过程,基于建议地震倒塌判别准则研究了钢框架结构的损伤演化规律。结果表明:钢框架结构在强震作用下的损伤发展顺序与塑性发展顺序一致;基于推覆分析的结构整体损伤模型能较好的体现强震作用下钢框架结构的损伤演化规律,且在上下界处收敛;强震作用下,钢框架结构的初始损伤主要由结构的残余侧移引起,而后期损伤主要由结构的承载力和刚度退化引起。     

Coupled dynamic elastic–plastic analysis of earth structures

A fully coupled finite element code based on mixture theory is developed. Prévost's multi-surface constitutive model is tailored to three-dimensional loads and used to predict effective stresses. A new viscous boundary is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb the two dilatational waves and the shear wave.Two soil deposits and two dams, with different slopes, composed by loose and dense sands have been subjected to the Pacoima accelerogram. Results show how the liquefaction propagates in the soil deposits and earth dams. The importance of the coupling between dilatancy–contractancy and filtration is highlighted by a parametric investigation. Phenomena such as liquefaction and cyclic mobility are reproduced, indicating the robustness of the constitutive model and finite element simulations. As an outcome of the parametric analysis, the seismic stability of dams cannot be improved by decreasing the upstream or downstream slopes.