基于断裂力学的重力坝有限元非线性动力分析

针对国外位于强震区的某重力坝,利用基于Hillerborg断裂力学理论的混凝土非线性本构模型,精细地建立了三维地基与坝体相互作用的有限元计算模型。根据坝址的地震动参数,拟合了水平向、竖直向的两条人工地震波,同时施加到重力坝体上,进行了非线性动力时程法的分析计算,揭示了重力坝复杂的地震响应变化规律。数值计算成果满足规范设计要求,为重力坝结构抗震设计提供了一定的科学参考。     

基于Marc二次开发实现粘弹性人工边界和振型分解反应谱法

新版《水电工程水工建筑物抗震设计规范》（NB35047-2015）中规定:"抗震设防类别为甲类的混凝土坝应考虑远域地基的辐射阻尼效应"。针对通用商业有限元软件Marc在振型分解反应谱法中的不足,且缺少时程分析法中模拟地基辐射阻尼效应的粘弹性人工边界,采用Fortran语言编制相应的独立程序及二次开发程序,以便在Mrac软件中精确实现新规范要求下的混凝土坝抗震安全评价。数值算例和工程实例分析结果验证了基于Marc二次开发实施思路和自编程序的正确性;重力坝抗震薄弱部位主要为坝体断面突变处,考虑无限地基的辐射阻尼效应后,坝体地震动力响应明显下降。     

高拱坝体系整体抗震安全评价方法研究

采用非线性地震波动反应分析方法，将坝体、地基、库水的强震反应本质上作为满足体系中接触面边界约束条件的波传播问题，在时域内以显式有限元方法求解。在不同概率水平的地震波作用下进行高拱坝与地基体系的有限元时程分析，并建议了以坝体位移反应发生突变为依据的拱坝体系整体失稳判断准则，从而建立了进行高拱坝体系整体抗震安全评价的方法，并以小湾拱坝为算例，对其抗震安全度作出初步评价。     

基于IDA的金安桥混凝土重力坝潜在失效模式研究

为了研究混凝土重力坝在地震动荷载作用下的潜在失效模式,以金安桥碾压混凝土重力坝5号非溢流坝段为例,运用粘弹性边界法和流固耦合法建立了反映重力坝在地震动作用下动力响应特征的坝体-地基-库水抗震分析模型。基于增量动力分析（IDA）法:绘制了以相对位移转角为x轴（损伤指标,DM）和峰值地面加速度为y轴（强度指标,IM）的IDA曲线簇;分析了金安桥大坝在极端荷载作用下的潜在失效模式和其在不同峰值地面加速度下重力坝的损伤破坏过程。结果表明:金安桥大坝在地震动荷载作用下,可能发生功能失效的地方多出现在坝体折坡处、碾压分区交界处、坝踵与坝基交界处、廊道顶等应力集中处。因此,加强对这些区域的抗震防护有利于提高大坝整体的抗震水平。     

关于坝体震害预测的研究

本文在“嘉峪关市地震危险性分析”工作的基础上，以大草滩水库土石坝为例，研究了坝体震害预测的方法。用动力非线性有限元法计算并分析了坝体的非线性动力特性，给出了５０年超越概率为１０％和２％两种风险水平下，坝体表面可能出现的拉裂区域和坝体内可能出现的塑性破坏区域，并给出了有代表性的一些地点的加速度时程和动剪应力时程。本文还用拟静力法对坝体的整体稳定性进行了分析计算，给出上述两种风险水平下，在可能遇到的各种情况下，坝体的整体安全可靠性能。     

强震作用下深厚砂质覆盖层土坝有效应力动力分析

利用基于Biot的饱和多孔介质理论和砂土多重机构模型的动力分析有限元程序FLIP,对遭受M6.7地震的国外某深厚砂质覆盖层土坝进行有效应力动力分析,研究坝体和地基的动力反应特性及其超静孔隙水压力的分布规律。通过对坝体加速度和永久变形的计算结果与现场实测数据的比较分析,证明两者之间存在一定差异,但计算结果基本上反映坝体加速度与永久变形的实际分布特征,从而说明采用的数值计算方法和本构模型具有一定精度。根据计算结果可以得出:坝体无液化发生;坝底上游浅层地基可能会发生局部液化,但范围较小,可以不进行加固处理;坝趾附近浅层地基可能会发生较大范围的液化,因此须采取相应的抗液化加固措施。     

考虑坝体材料非线性的混凝土拱坝地震反应分析

本文基于可信概率水准的破坏性强震作用,针对小湾高拱坝进行了考虑坝体材料非线性的拱坝地震反应分析。在分析模型中,同时考虑了无约束域地震能量辐射效应和近域地基材料非均匀性的影响。为了实现非线性条件下的静、动力组合分析,利用显式有限元结合修正的黏弹性人工边界的开放系统时域静、动力统一分析方法进行了求解,对在自重作用下的初始静力解计算采用了动力松弛技术。     

基于高阶双渐近透射边界的重力坝-层状地基动力相互作用分析

时域高阶双渐近透射边界能够同时模拟层状介质中行波和快衰波的传播,具有很高的计算精度和计算效率.本文将高阶双渐近透射边界推广应用到多层层状地基系统弹性波传播问题的模拟,采用广义特征值分解分析该透射边界的数值稳定性,通过移谱法消除导致数值不稳定的虚假模态.将高阶双渐近透射边界以超单元的形式直接嵌入到近场有限元方程,建立了有限元-高阶双渐近透射边界时域耦合分析模型,并将其应用到重力坝-层状地基动力相互作用分析.数值算例分析结果表明,该时域耦合分析模型具有很高的精度和计算效率,适用于实际重力坝工程的地震响应分析.     

基于能量耗散的重力坝地震响应时频特征分析

为了从能量角度研究重力坝地震响应的时频特征,在重力坝非线性动力分析基础上,探讨了重力坝地震过程中振动能量的时域耗散机制;采用小波频域多层次分解技术研究了其动响应的分频段能量特征,得到了坝体典型位置动响应分频段振动能量随高程的变化规律.通过分析发现:结构地震能量耗散为时域上不可逆的增加,坝体损伤集中出现在地震过程的某个时间段,地震动峰值后坝体损伤状态基本稳定;小波分解可以较全面地描述结构动响应能量的分频段特征.对于本文算例,在坝踵和上游折坡附近,重力坝地震响应的振动能量以4～8 Hz频段为主,这与输入的地震信号分频段特征一致;而坝顶附近则以1 ～4 Hz的振动能量为主,高频能量分量的比重随高程增大而逐渐减小.     

深厚覆盖层上某核电护岸结构地震响应及稳定性分析

以某深厚覆盖层上核电护岸工程为研究对象,介绍了护岸抗震稳定分析的原理和分析方法。采用等价线性法对代表性断面进行有限元动力时程分析,得到SL1和SL2地震动作用下护岸结构的加速度响应、安全系数、液化区分布状况以及震后残余变形。通过分析地基不同土层交界面加速度响应变化来研究深厚覆盖层对加速度峰值的动力放大效应影响。     

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.     

Seismic stability assessment of an arch dam-foundation system

A seismic stability assessment of arch dam-foundation systems is presented using a comprehensive approach,in which the main factors that significantly influence the seismic response of an arch dam-foundation system are considered.A large scale finite element model with over 1 million degrees of freedom is constructed for the Baihetan arch dam(289 m high),which is under construction in the Southwest of China.In particular,the complicated geological conditions with faults intersecting interlayer shear weakness zones at the dam base and the dam abutment resisting force body is modeled in the analysis.Three performance indices are adopted to assess the seismic stability of the arch dam.The results demonstrate that the opening of the joints of the Baihetan arch dam is small and the water stop installed between the joints would not be torn during a design earthquake.The yielding formed in the interface between the dam and foundation does not reach the grouting curtain that would remain in an elastic state after an earthquake.The yielding zones occurring on the upper portion of the dam faces extend 1/8 thickness of block section into the dam body and thus cantilever blocks need not be concerned with sliding stability.The faults and interlayer shear weakness zones in the near field foundation exhibit severe yielding,and a potential sliding surface is penetrated.Although the factor of safety against sliding of the surface fluctuates with a decreased trend during an earthquake,the minimum instantaneous value reaches 1.02 and is still larger than 1.0.Therefore,a conclusion is drawn that the Baihetan arch dam-foundation system will remain stable under the design earthquake.     

AN INVESTIGATION INTO THE BASE SLIDING RESPONSE OF RIGID CONCRETE GRAVITY DAMS TO DYNAMIC LOADING

A series of dynamic slip tests on a concrete gravity dam model was conducted on a shaking table. The aim of the experiments was to investigate the dynamically induced sliding and overturning characteristics of a typical low height gravity dam monolith cracked at its base. Tests indicated that downstream sliding is the main instability that could be expected during an earthquake. Dynamic, finite element analyses of the experimental model, using a Lagrangian contact surface algorithm, were also performed. A comparison of the experimental and analytical responses indicated that the seismically induced slip can be predicted reasonably by such a contact surface algorithm implemented in a standard finite element package. A comparison of observed displacements with Newmark's sliding block displacements indicated that a conservative estimate of seismic induced slip of a gravity dam could be obtained by using Newmark's sliding block concept, generally adopted for earth dams and embankments.     

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.     

Seismic stability of concrete gravity dams

Linear finite element analyses are commonly used to simulate the behaviour of gravity dam—foundation systems. However, the foundation is generally unable to develop any significant tensile stresses. Therefore any tension occurring in the vicinity of the dam—foundation interface is largely fictitious. Moreover, the traditional overturning and sliding stability criteria have little meaning in the context of the oscillatory response of dams during earthquakes. In this study, time domain analyses using non-linear contact elements located at the dam—foundation interface have been used to determine the dynamic sliding and uplifting response of gravity dam monoliths considering various elastic foundation properties. The magnitudes of the relative interface displacements, of the percentage of base not in contact (PBNC) and of the compressive stresses at the heel or toe of the dam have been used to monitor the seismic stability. The numerical results have shown that the non-linear behaviour of the dam—foundation interface reduces the seismic response of the system, indicating the possibility of more rational and economical designs. The PBNC was identified as the critical seismic stability response parameter for all analyses except for very flexible foundation conditions where the maximum values of relative interface displacements need to be considered.     

Seismic assessment of concrete gravity dams using capacity estimation and damage indexes

A new concept to determine state of the damage in concrete gravity dams is introduced. The Pine Flat concrete gravity dam has been selected for the purpose of the analysis and its structural capacity, assuming no sliding plane and rigid foundation, has been estimated using the two well‐known methods: nonlinear static pushover (SPO) and incremental dynamic analysis (IDA). With the use of these two methods, performance and various limit states of the dam have been determined, and three damage indexes have been proposed on the basis of the comparison of seismic demands and the dam's capacity. It is concluded that the SPO and IDA can be effectively used to develop indexes for seismic performance evaluation and damage assessment of concrete gravity dams. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

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

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

挡土墙地震被动土压力的拟动力分析

对地震土压力的研究是地震区挡土墙安全设计的一项重要课题。地震条件下,目前的研究主要是给出了土压力的近似拟静力解析解。本文采用可考虑动力荷载下的周期和纵波及横波效应的拟动力方法,对挡土墙后的地震被动土压力进行分析。在挡土墙后平面滑裂面假设的基础上,考虑了水平和垂直向地震加速度、纵波速度、横波速度、挡土墙摩擦角、填土内摩擦角、填土坡角对地震被动土压力的影响。与Mononobe-Okabe理论的拟静力法不同的是,用本方法得出了沿墙身地震被动土压力是非线性变化的结果,这更符合地震条件下土压力的变化规律。     

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

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