Numerical simulation of seismic damage and cracking of concrete slabs of high concrete face rockfill dams

Based on the damage constitutive model for concrete, the Weibull distribution function was used to characterize the random distribution of the mechanical properties of materials by finely subdividing concrete slab elements, and a concrete random mesoscopic damage model was established. The seismic response of a 100-m high concrete face rockfill dam(CFRD), subjected to ground motion with different intensities, was simulated with the three-dimensional finite element method(FEM), with emphasis on exploration of damage and the cracking process of concrete slabs during earthquakes as well as analysis of dynamic damage and cracking characteristics during strong earthquakes. The calculated results show that the number of damaged and cracking elements on concrete slabs grows with the duration of earthquakes. With increasing earthquake intensity, the damaged zone and cracking zone on concrete slabs grow wider. During a 7.0-magnitude earthquake, the stress level of concrete slabs is low for the CFRD, and there is almost no damage or slight damage to the slabs. While during a 9.0-magnitude strong earthquake, the percentages of damaged elements and macrocracking elements continuously ascend with the duration of the earthquake, peaking at approximately 26% and 5% at the end of the earthquake, respectively. The concrete random mesoscopic damage model can depict the entire process of sprouting, growing, connecting, and expanding of cracks on a concrete slab during earthquakes.

Seismic responses of high concrete face rockfill dams:A case study

Seismic responses of the Zipingpu concrete face rockfill dam were analyzed using the finite element method. The dynamic behavior of rockfill materials was modeled with a viscoelastic model and an empirical permanent strain model. The relevant parameters were obtained either by back analysis using the field observations or by reference to parameters of similar rockfill materials. The acceleration responses of the dam,the distribution of earthquake-induced settlement, and the gap propagation under the concrete slabs caused by the settlement of the dam were analyzed and compared with site investigations or relevant studies. The mechanism of failure of horizontal construction joints was also analyzed based on numerical results and site observations. Numerical results show that the input accelerations were considerably amplified near the top of the dam, and the strong shaking resulted in considerable settlement of the rockfill materials, with a maximum value exceeding 90 cm at the crest.As a result of the settlement of rockfill materials, the third-stage concrete slabs were separated from the cushion layer. The rotation of the cantilever slabs about the contacting regions, under the combined action of gravity and seismic inertial forces, led to the failure of the construction joints and tensile cracks appeared above the construction joints. The effectiveness and limitations of the so-called equivalent linear method are also discussed.

Dynamic response of concrete face rockfill dam affected by polarity reversal of near-fault earthquake

In China,an increasing number of high concrete face rockfill dams(CFRDs)are located in high intensity earthquake zones,some of which are close to the seismic fault line.Recordings suggest that near-fault ground motions are characterized by large one-sided velocity pulses.The conventional dynamic analysis of dams,however,neglects the features of strong ground movements.In this study,under different ground motion levels some numerical dynamic studies considering the one-sided broadband pulses of near-fault earthquakes are presented for CFRDs based on a generalized plasticity model for rockfill materials.The results indicate that the displacements of dam crest corresponding to positive and reverse input of near-fault ground motion make a significant difference,while the displacements of the dam crest under artificial seismic waves are similar.Furthermore,using the horizontal and vertical components as simultaneous excitations near the faults,the displacements of the dam crest before and after reversing the motion produce a larger difference than that using a single component.More importantly,the difference of horizontal displacements of the dam crest caused by polarity reversal of near-fault ground motions increases with the increase of earthquake intensity.Due to the randomness and uncertainties of earthquakes,using a stochastic near-field motion input as excitation without considering the polarity(i.e.,positive vs reversed waveform),does not necessarily obtain a conservative result.

Fuzzy Earthwork Dynamic Allocation and Optimization for Construction of High Concrete Face Rockfill Dam

Due to the complexity of earthwork allocation system for the construction of high concrete face rockfill dam,traditional allocation and planning are not able to function properly in the construction process with strong randomness.In this paper,the working mechanism of earthwork dynamic allocation system is analyzed comprehensively and a solution to fuzzy earthwork dynamic allocation is proposed on the basis of uncertain factors in the earthwork allocation of a hydropower project.Under the premise of actual situation and the experience of the construction site,an all-coefficient-fuzzy linear programming mathematical model with fuzzy parameters and constraints for earthwork allocation is established according to the structure unit weighted ranking criteria.In this way,the deficiency of certain allocation model can be overcome.The application results indicate that the proposed method is more rational compared with traditional earthwork allocation.

Seepage simulation of high concrete-faced rockfill dams based on generalized equivalent continuum model

This research focused on the three-dimensional(3 D) seepage field simulation of a high concrete-faced rockfill dam(CFRD) under complex hydraulic conditions. A generalized equivalent continuum model of fractured rock mass was used for equivalent continuous seepage field analysis based on the improved node virtual flow method. Using a high CFRD as an example, the generalized equivalent continuum range was determined, and a finite element model was established based on the terrain and geological conditions, as well as structural face characteristics of the dam area. The equivalent seepage coefficients of different material zones or positions in the dam foundation were calculated with the Snow model or inverse analysis. Then, the 3 D seepage field in the dam area was calculated under the normal water storage conditions, and the corresponding water head distribution, seepage flow, seepage gradient, and seepage characteristics in the dam area were analyzed. The results show that the generalized equivalent continuum model can effectively simulate overall seepage patterns of the CFRD under complex hydraulic conditions and provide a reference for seepage analysis of similar CFRDs.

A simplified physically-based breach model for a high concrete-faced rockfill dam:A case study

A simplified physically-based model was developed to simulate the breaching process of the Gouhou concrete-faced rockfill dam （CFRD）, which is the only breach case of a high CFRD in the world. Considering the dam height, a hydraulic method was chosen to simulate the initial scour position on the downstream slope, with the steepening of the downstream slope taken into account; a headcut erosion formula was adopted to simulate the backward erosion as well. The moment equilibrium method was utilized to calculate the ultimate length of a concrete slab under its self-weight and water loads. The calculated results of the Gouhou CFRD breach case show that the proposed model provides reasonable peak breach flow, final breach width, and failure time, with relative errors less than 15% as compared with the measured data. Sensitivity studies show that the outputs of the proposed model are more or less sensitive to different parameters. Three typical parametric models were compared with the proposed model, and the comparison demonstrates that the proposed physically-based breach model performs better and provides more detailed results than the parametric models.

Reservoir water effects on earthquake performance evaluation of Torul Concrete-Faced Rockfill Dam

This study presents earthquake performance analysis of the Torul Concrete-Faced Rockfill （CFR） Dam with two-dimensional dam-soil and dam-soil-reservoir finite element models. The Lagrangian approach was used with fluid elements to model impounded water. The interface elements were used to simulate the slippage between the concrete face slab and the rockfill. The horizontal component of the 1992 Erzincan earthquake, with a peak ground acceleration of 0.515g, was considered in time-history analysis. The Drucker-Prager model was preferred in nonlinear analysis of the concrete slab, rockfill and foundation soil. The maximum principal stresses and the maximum displacements in two opposite directions were compared by the height of the concrete slab according to linear time-history analysis to reveal the effect of reservoir water. The changes of critical displacements and principal stresses with time are also shown in this paper. According to linear and nonlinear time-history analysis, the effect of the reservoir water on the earthquake performance of the Torul CFR Dam was investigated and the possible damage situation was examined. The results show that the hydrodynamic pressure of reservoir water leads to an increase in the maximum displacements and principal stresses of the dam and reduces the earthquake performance of the dam. Although the linear time-history analysis demonstrates that the earthquake causes a momentous damage to the concrete slab of the Torul CFR Dam, the nonlinear time-history analysis shows that no evident damage occurs in either reservoir case.

Seismic response of concrete gravity dam reinforced with FRP sheets on dam surface

This paper aims at exploring the effects of anti-seismic reinforcement with the fiber-reinforced polymer （FRP） material bonded to the dam surface in dam engineering. Time-history analysis was performed to simulate the seismic failure process of a gravity dam that was assumed to be reinforced at the locations of slope discontinuity at the downstream surface, part of the upstream face, and the dam heel. A damage model considering the influence of concrete heterogeneity was used to model the nonlinearity of concrete. A bond-slip model was applied to the interface between FRP and concrete, and the reinforcement mechanism was analyzed through the bond stress and the stress in FRP. The results of the crack pattern, displacement, and acceleration of the reinforced dam were compared with those of the original one. It is shown that FRP, as a reinforcement material, postpones the occurrence of cracks and slows the crack propagation, and that cracks emanating from the upstream surface and downstream surface are not connected, meaning that the reinforced dam can retain water-impounding function when subjected to the earthquake. Anti-seismic reinforcement with FRP is therefore beneficial to improving the seismic resistant capability of concrete dams.

超深覆盖层上面板坝防渗墙连接板接缝变形机制及优化分析

针对超深软弱覆盖层上的面板坝工程止水缝安全问题,从变形和应力2个方面开展优化研究,揭示了防渗墙连接板趾板面板防渗体系的变形模式,定位了薄弱位置,提出了优化方案,并量化了工程应用效果.结果表明,防渗墙与趾板间沉降差集中在防渗墙与连接板之间的止水缝处,设置多块连接板无法起到逐步过渡沉降差的作用.防渗墙下游侧设置侧向支撑、制作倾斜形止水缝、趾板下设支撑墙均可减小防渗墙连接板间止水缝沉陷变形.相比垂直止水缝和多块连接板的传统方案,所提的防渗墙顶部拓宽连接板倾斜形止水缝以及防渗墙顶部拓宽倾斜形止水缝趾板下设支撑墙的深厚覆盖层上面板坝防渗体系综合优化方案可使防渗墙与连接板之间止水缝沉陷变形分别减小41.5%和53.8%.

250 m级特高面板堆石坝混凝土面板顺坡向应力控制措施研究

面板应力变形控制是建设250 m级特高混凝土面板堆石坝的关键难题.大量工程实测案例和研究计算表明,特高面板坝压应力常处于较高水平,对于坝轴向应力来说,采用柔性缝降低压应力措施已被工程界普遍接受,然而,目前对于如何降低面板顺坡向压应力尚缺少有效的解决方法.本文以拟建的240 m级GS面板坝为研究对象,从特殊增模区布置、面板施工顺序、面板低高程柔性横缝、高压力区高强度混凝土材料等4个方面提出改善面板顺坡向应力的工程措施.文中通过三维有限元精细模拟方法分析了提出的顺坡向应力改善措施的效果,文中提出的措施可为250 m面板堆石坝面板应力控制提供借鉴.

堆石料本构模型对混凝土面板坝应力变形计算结果的影响研究

堆石料本构模型选择是影响混凝土面板坝应力变形计算结果的主要因素。以某典型面板坝为例,采用邓肯E-B非线性弹性模型和“南水”双屈服面弹塑性模型,对坝体填筑过程和蓄水过程进行了三维有限元模拟;研究了两种模型计算的坝体和面板的位移和应力分布差异;分析了两种模型计算结果呈现差异的原因。两种模型计算结果最显著的差异体现在两个方面:(1)“南水”模型计算的坝体沉降、面板挠度等位移指标小于E-B模型,其原因是“南水”模型计算的坝体小主应力和变形模量更高。(2)E-B模型计算结果显示蓄水后面板底部顺坡向受拉;而“南水”模型计算的面板顺坡向应力为全断面受压。前者是由E-B模型的各向同性弹性本质决定的,垫层料因水压力作用顺坡向膨胀,产生沿坡面向上的位移;后者是其F1屈服面持续扩张产生塑性体积收缩导致的,塑性体缩抵消了顺坡向弹性膨胀,使垫层料产生沿坡面向下的位移。

高纬度高寒区面板堆石坝防渗面板修复方案及经验

JLBLK水电站大坝是建于高纬度高寒区的百米级全断面堆石面板坝,最大坝高140.3 m。通过降水处理,现场发现18^(#)、19^(#)面板挤压破损程度严重且有加剧趋势,25^(#)、26^(#)、27^(#)面板断裂。经与国内外同类工程比较,水库渗漏量总体偏大,变形量偏大,且尚未收敛,三维流变反演成果表明后期尚有17%的沉降量。为保证修复效果满足后期变形和渗透稳定需要,针对18^(#)、19^(#)面板水上破损区采用高于原设计强度混凝土分小区进行补强,并结合工程特点和设计方案提出适合本工程的修复施工方案,修复后面板各项指标均满足要求。

水循环荷载作用下高面板堆石坝长期变形特性研究

高面板堆石坝的长期变形是面板坝建设中的关键技术问题之一,已有研究大多将长期变形归结为堆石体蠕变的结果,但某些高面板堆石坝的变形原型监测资料显示,仅从堆石体蠕变角度难以解释高面板堆石坝长期变形的机理。通过水布垭面板堆石坝17a的变形监测资料分析和室内低频循环荷载作用下堆石料变形特性试验,探讨了实际运行中库水位周期性变化形成的循环荷载对高面板堆石坝长期变形特性的影响。结果表明:(1)坝体填筑和初次蓄水引起堆石体的瞬时变形和蠕变,堆石体的后期变形主要是水循环荷载长期作用的结果。(2)水循环荷载作用下,堆石体的变形表现为两种形式,其一为不可恢复的永久变形(亦称为残余变形),其二为可恢复的弹性变形。(3)室内低频循环荷载作用试验得到的堆石料变形及其增长规律很好地验证了水循环荷载与堆石坝后期变形的相互关系以及后期变形的力学意义。研究成果为高面板堆石坝长期变形特性研究提供了新的思路,并为具有往复大消落特点的高面板堆石坝变形及安全控制提供借鉴。

强震区抽蓄电站沥青混凝土面板堆石坝抗震稳定性分析

我国西北强震区某拟建抽水蓄能电站上水库采用沥青混凝土面板堆石坝,最大坝高161 m。为分析该坝抗震稳定性,在大坝非线性静力分析基础上,采用改进的等效粘弹性模型和有限元动力时程法等理论,对大坝开展三维有限元地震反应和坝坡抗震稳定性计算。结果表明,大坝在100年1%超越概率的校核地震(地震加速度为0.48g)作用下,坝体及面板的动力反应分布规律合理。其中,坝体顺河向、坝轴向、竖向加速度极值分别为13.93、13.91、13.80 m/s^(2)。考虑静动叠加后,在沥青混凝土面板反弧段的拉应变极值达0.95%,但小于改性沥青混凝土抗拉强度的一般允许值。因此,大坝即使遭遇0.48g的校核地震作用,亦能保持较好的整体安全性,不会出现重大抗震安全事故。

P波空间斜入射下沥青混凝土面板堆石坝响应特性及面板破坏评价

地震波入射方位和斜入射角度对于沥青混凝土面板堆石坝地震响应结果有着显著影响。本文从地震动传播机制出发,基于波场叠加原理,推导了P波以任意入射方位角和斜入射角空间斜入射下弹性半空间自由场计算公式,建立了P波空间斜入射波动输入模型。以不同应变率下的沥青混凝土单轴动态拉伸试验结果为基础,提出了考虑抗拉强度随应变率实时变化的沥青混凝土面板破坏评价方法。考虑16种不同入射方位角和斜入射角下的P波空间斜入射工况,从面板应力和加速度两方面分析了入射方位角和斜入射角对某实际沥青混凝土面板地震响应特性的影响,并进行了沥青混凝土面板的抗拉破坏评价,分析了静态评价方法与动态评价方法的差异。结果表明:所建立的P波空间斜入射波动输入模型能够准确模拟半空间自由场,波动输入数值解与解析解符合良好;与垂直输入射相比,面板顺水流向、垂直水流向加速度峰值最大分别增加135.9%和92.7%,竖向加速度峰值最大减少68.3%;空间斜输入导致面板动应力增大,主拉应力和主压应力最大分别增加3.6倍和2.7倍,忽略入射方位角和斜入射角可能会严重低估面板的地震响应。传统将静态强度提高30%作为动态强度的面板抗拉破坏判别方法过于严格,本文提出的考虑应变率效应的沥青混凝土面板破坏评价方法较为符合实际。

大倾角坝基上混凝土面板堆石坝滑移变形特征研究

抽水蓄能电站是实现“双碳”战略目标的重要举措,混凝土面板堆石坝作为抽水蓄能电站上水库大坝的优选坝型之一,往往坝基顺河向坡度较陡,坝体容易产生顺坡向滑移变形,坝坡抗滑稳定性较差,危及抽水蓄能电站上水库大坝运行安全。因此,建立大倾角坝基上混凝土面板堆石坝数值模型,研究大倾角坝基上混凝土面板堆石坝滑移变形的规律,揭示大倾角坝基上混凝土面板堆石坝产生滑移变形的机制,提出减小其滑移变形的物理力学和几何结构设计措施。研究结果表明,较大下滑力和较小抗滑力是大倾角坝基上混凝土面板堆石坝容易产生滑移变形的根本原因;坝基倾角越大,重力沿坡面分力越大,坝体滑移变形越大,坝坡安全系数越小(当坝基倾角大于15°时,坝坡安全系数小于规范要求的最小值1.5);增设压坡能够有效减小坝体滑移变形,显著提高坝坡抗滑稳定性,压坡平台高度越高、宽度越宽和黏聚力、内摩擦角越大,坝体滑移变形越小,坝坡安全系数越大,推荐压坡平台高度和宽度取1/2倍主坝最大坝高,压坡密实度(黏聚力和内摩擦角)与主坝堆石密实度相当。

本构模型对面板堆石坝动力特性影响及敏感性分析

为研究堆石料本构模型及模型参数对面板堆石坝动力响应特性的影响,基于有限元方法和动力时程分析,分别采用邓肯-张E-B模型和南水模型对某沥青混凝土堆石坝进行动力计算,比较了两种本构模型下坝体的静动位移、动应力和地震永久变形的分布规律和数值大小,在此基础上针对堆石料参数设计了正交试验,采用极差和方差分析方法研究了参数敏感性。结果表明,两种本构模型下面板坝的动位移分布规律相似但数值大小不同,动应力大小和分布规律均不同,参数K和Rd分别是邓肯-张E-B模型和南水模型中对坝体动位移最敏感的参数。

基于伪随机场的面板坝运行期库岸渗漏成因分析

四川某混凝土面板堆石坝坝址区地形条件复杂,坝后量水堰处渗水量相比同类工程观测数据偏大,为此,基于饱和—非饱和渗流原理,利用自编程准三维渗流分析程序建立渗流分析模型,在伪随机场探测结果的基础上分别对导流洞封堵区、左岸防渗帷幕、止水及右岸渗漏通道进行敏感性分析。结果表明,导流洞封堵区渗水及左岸趾板下死水位以下的上层防渗帷幕失效是造成左岸渗流量增大的重要原因;左岸面板周边缝止水的不同失效性仅对面板附近渗压造成局部影响,表层止水对测点渗压影响更大;右岸存在渗漏通道时,右岸测点拟合效果更好,右岸山体渗流量增大,但量水堰处流量并未同等程度增加。

近断层地震动组合斜入射下沥青混凝土面板堆石坝响应研究

西部地区水利工程多建于复杂地震断裂带附近,因此近场地震动的脉冲效应对水工建筑物所产生的影响不可忽略。近断层地震动往往由P波、SV波等共同组成,到达地表后明显带有一定角度,垂直入射假定不再适用,因此有必要开展近断层脉冲型地震动组合斜入射对沥青混凝土面板堆石坝地震响应影响的研究。基于地表二维地震动反演确定P波、SV波斜入射时程,构建近断层地震动组合斜输入模型,选取脉冲型和非脉冲型各5条地震动记录,研究近断层地震动组合斜入射对沥青混凝土面板加速度、应力和坝体永久变形的影响规律。结果表明:与地震动垂直入射相比,面板水平向加速度随P波角度的增大而增大,随SV波角度的增大而减小,竖向加速度及坝体永久变形的规律则相反;面板拉应力随着P波、SV波角度的增大而增大;近断层地震动的脉冲特性对堆石坝的响应有显著放大作用,其中面板拉应力最大增加了116%。

考虑施工质量与填筑进度的面板坝精细有限元分析

针对混凝土面板堆石坝施工质量空间上的差异性导致的坝料物理力学参数具有空间差异性的问题,提出了考虑实际施工质量与填筑进度(面貌)的精细有限元分析方法。首先基于试坑试验数据,采用空间插值算法对大坝压实质量进行空间估计;然后建立各分区坝料本构模型参数与压实质量之间的定量关系,实现对任意单元模型参数的空间估计与赋值;最后提出考虑填筑进度(面貌)的精细建模方法,并结合工程实例对混凝土面板堆石坝的应力变形等性态进行了精细模拟与分析。结果表明:相比于基于设计参数模拟得到的沉降量平均误差10.46 cm,采用考虑实际施工质量与填筑进度的精细分析方法模拟得到的沉降量平均误差仅为2.41 cm,其结果与实测值具有更好的一致性。研究成果为混凝土面板堆石坝的结构性态精准分析提供了新的途径。