复杂层状岩基上重力坝极限抗震能力评估方法初探_张社荣

第46卷 第3期 2013年3月天津大学学报(自然科学与工程技术版)Journal of Tianjin University (Science and Technology )V ol.46 No.3Mar. 2013收稿日期：2012-01-09；修回日期：2012-05-23.基金项目：国家自然科学基金创新研究群体科学基金资助项目(51021004)． 作者简介：张社荣（1960— ），男，博士，教授，tjuzsr@. 通讯作者：王 超，tjudam@.复杂层状岩基上重力坝极限抗震能力评估方法初探张社荣，王 超，孙 博，王高辉(天津大学水利工程仿真与安全国家重点实验室，天津 300072)摘 要：针对复杂层状岩基上典型工程的特点，从结构系统功能的潜在失效模式出发，利用数值仿真分析方法，结合基于功能的地震破坏等级评价模型、基于断裂力学的坝体开裂行为分析模型和基于变形体突变的系统极限承载力分析评价模型，建议从坝体混凝土损伤破坏等级、开裂破坏模式、坝基岩体极限承载能力和系统稳定性突变等方面综合评价大型重力坝的极限抗震能力．实例分析结果表明，阿海重力坝溢流坝段极限抗震能力为0.550g ～0.600g .基于数值试验的极限抗震能力综合评估方法，可以全面考虑各种关键性能对极限抗震能力的影响，避免了单一收敛性或系统突变评价准则的缺陷，在一定程度上弥补了振动台模型试验的不足． 关键词：重力坝；极限抗震能力；开裂模式；层状岩基；阿海重力坝中图分类号：TV641.1 文献标志码：A 文章编号：0493-2137(2013)03-0202-09Evaluation Method of Maximum Aseismic Capability of GravityDam on Complex Layered Rock FoundationZhang Sherong ，Wang Chao ，Sun Bo ，Wang Gaohui(State Key Laboratory of Hydraulic Engineering Simulation and Safety ，Tianjin University ，Tianjin 300072，China ) Abstract ：The performance failure of the gravity dam under ultimate seismic load is on important characteristic whichmeans its aseismic capability has reached the limit ，but there is no unified standard reference to make the evaluation.Based on the characteristics of typical engineering and the potential failure modes ，with performance-based seismic damage assessment model ，cracking behavior analysis model as well as the foundation stabilization safety and limit load carrying capacity analysis model ，the maximum aseismic capability was studied and evaluated by focusing the emphasis on damage levels ，cracking extending mode ，foundation load carrying capacity and displacement jump of typical position for the gravity dams on complex layered rock foundation. Example results show that the maximum aseismic capability of the overfall section of Ahai gravity dam is evaluated as 0.550g —0.600g . The proposed method has been successfully applied to the maximum aseismic capability evaluation by considering the comprehensive influ-ence of the key performance failure ．It av oids the shortage of the single conv ergence or abrupt system changeov er criteria ，and also overcomes the shortage of the experimental research to some extent.Keywords ：grav ity dam ；maximum aseismic capability ；cracking mode ；layered rock foundation ；Ahai grav itydam在重力坝全生命周期的运行中，大坝性能可能受到诸如超标洪水、地质条件恶化、建筑材料老化乃至远远超过设计地震动水平的极端地震荷载的威胁．特别是5.12汶川地震以来，重力坝极限抗震能力受到极大关注，逐渐成为抗震安全性评价中的热点问题．然而，由于世界上重力坝震害典型实例较少[1-2]，人们对极端地震作用下重力坝的破坏机理和失效模式的认识并不充分．研究混凝土重力坝极限抗震能力及其在地震过程中的破坏模式主要有室内模型试验、数值试验及理论分析等方法．范书立等[3]采用仿真混凝土材料制作龙开口水电站溢流坝段的模型，考虑坝体-库水的动2013年3月张社荣等：复杂层状岩基上重力坝极限抗震能力评估方法初探·203·力相互作用，研究了溢流坝段的坝体开裂形式和抗震能力．但是，由于影响混凝土重力坝动力响应和破坏形态的因素很多，诸如地基和混凝土的各向异性、不确定性等，模型试验无法全部考虑，因此还需结合数值计算结果才能做出全面的评价．此方面，昇李等[4]采用超地震荷载法研究强震作用下重力坝的破坏模式；周伟等[5]利用三维非线性有限元法，研究了坝基的渐进破坏模式和坝体坝基系统的极限承载能力；Ghanaat[6]根据线弹性时程分析得到的应力需求能力比、超应力累积持时及超应力分布范围等抗震性能参数，提出了混凝土坝抗震损伤程度的定量分析方法和评价标准；沈怀至等[7]提出了基于功能的混凝土重力坝地震破坏定量评价模型；Pekau等[8]使用离散单元法(discrete element method，DEM)研究混凝土重力坝脱离体的极限抗震能力．但是，对于重大工程典型实例，基于数值试验的极限抗震能力的分析评价方法及其适用性还有待深入探讨．阿海水电站工程区地震地质背景复杂，坝址区设计烈度接近Ⅸ度，其100a超越概率2%的基岩水平地震动峰值加速度(peak ground acceleration，PGA)达0.344g．大坝建基岩体主要为板岩砂岩互层结构，复杂层状岩体地基上的大坝极限抗震能力是本工程的关键技术问题．笔者以阿海水电站工程作为分析实例，结合基于功能的地震破坏等级评价模型、基于断裂力学的坝体非线性行为分析模型和基于收敛性判据的系统极限抗震稳定性分析评价模型，采用地震动超载时域分析方法，建议从坝体混凝土损伤破坏等级、开裂破坏模式、坝基岩体塑性区破坏扩展规律、坝体坝基系统稳定性突变等方面综合评价其极限抗震能力．1 评价模型及研究方法采用数值方法分析重力坝的极限抗震能力，目前还没有统一的评价标准可以参照，通常采用非线性计算收敛性准则．但是，所采用的数学模型、破坏准则、单元类型、边界条件等因素均可能对数值计算结果及其收敛性产生不同程度的影响，且目前缺乏有效的方法来消除这些因素的影响；且从基于性能的抗震设防标准及理论出发，往往是坝体某个关键性能(如大坝蓄水功能)破坏便会影响到坝体正常运行．因此，单一的评价准则已经不能满足大坝极限抗震能力评估的要求．极限地震荷载下的系统功能失效是重力坝达到极限抗震能力的重要表征．高混凝土坝破坏机理关键[9]在于2个方面，一是非连续(节理、断层、软岩)坝基岩体在蓄水或强地震下发生变形、滑动；二是大坝混凝土在蓄水与强地震作用下发生强度断裂．极端荷载条件下，其潜在失效模式主要体现在以下3个方面[5]：超拉应力开裂；倾倒突变失稳破坏；沿坝体坝基薄弱层面滑动．它们是评判大坝是否达到极限抗震能力的关键因素．采用数值分析方法时，就需要将结构系统的屈服破坏准则、塑性区贯通、系统状态突变等判据综合考虑．1.1 基于性能的地震破坏评价模型基于坝体安全评价理念和弹塑性损伤理论[7]，首先采用线弹性分析方法，结合Ghanaat[6]提出的破坏等级划分标准，初步定量评估大坝的损伤破坏程度．坝体线弹性分析中的需求能力比(demand to ca-pacity ratio，DCR)最大允许值以地震荷载下名义抗拉强度即DCR＝2为分界点，结合超应力累积持时，将地震荷载下坝体破坏程度分为低-中等与严重损伤2个区域．然后，采用能够模拟混凝土非弹性现象且能够严格定义破坏面演化准则的连续式塑性损伤力学模型，进行坝体坝基系统的非线性分析，研究其塑性损伤行为规律．其损伤本构方程为,klkl pl,kl el,,kl kl pl,kl:():()1ijij ijσεεεε=−=−−DDd(1)式中：ijσ为有效应力；,klijD为损伤刚度矩阵；d为损伤张量，其中包含拉压损伤因子；el,,klijD为初始未损伤的弹性刚度矩阵；klε为总应变；pl,klε为塑性应变[10].1.2 基于断裂力学的坝体开裂破坏模型采用William-Warnke 五参数准则[11]判断混凝土的屈服和破坏．在屈服以前，采用等效单轴模型模拟混凝土的本构关系，在混凝土开裂后，采用无网格方法的弥散裂缝模型[12]模拟其开裂行为．弥散裂缝模型通过非弹性开裂应变等效模拟裂缝，通过调整材料软化本构关系，无需改变单元形式或重新划分单元网格；引入垂直于裂缝表面方向上的一个缺陷平面来表示在某个积分点上出现开裂．利用William-Warnke 五参数屈服准则可以确定混凝土的最终强度曲面，保证了破坏面的连续性．William-Warnke五参数(即单轴抗拉强度tf、单轴抗压强度cf、二轴等压强度cb f、拉子午线上的三轴强度1f以及压力子午线上的三轴强度2f)准则，2个子午面之间的椭圆极坐标关系为22221212(,){2()cos(2)mr r r r r r rσηη=−+−⋅·204· 天津大学学报(自然科学与工程技术版) 第46卷 第3期12222221112[4()cos54]}/r r r r r η−+−222221214()cos (2)r r r r η⎡⎤−+−⎣⎦(2)式中1r 和2r 分别为受拉子午面和受压子午面上的极坐标值，且有212312cos (2)/()ησσσσσ=−−−+}1/2222331()()σσσσ⎤−+−⎦(3)当单元达到开裂状态时，单元在最大主应力垂直的方向形成无数平行的微裂纹，单元发生损伤，即对单元的本构关系进行修正，并在本构模型中引进剪力因子考虑开裂的软化特性．能自动确定裂缝起裂、扩展的完整信息，具有较高的计算效率．基于此，在动力计算中，不必随着裂缝出现和发展而改变有限元网格划分，可以有效表达动力工况下混凝土开裂行为．混凝土拉裂后破坏面的法向刚度和剪切刚度降低，如果裂缝重新闭合，又可重新传递压应力，残留抗剪系数增大．混凝土压碎后，单元刚度变得很小，对总体刚度的贡献可以忽略．沿裂缝面建立正交坐标系，以应力应变全量的形式描述材料开裂后的力学行为(平面应力问题)，经过改进[13]，混凝土材料发生开裂后应力应变关系[14]为 22n n s s 22ns ns 0110110EEEE G μνμνμνμσενμσενμνμτγβ⎡⎤⎢⎥−−⎢⎥⎧⎫⎧⎫⎪⎪⎪⎪⎢⎥=⎨⎬⎨⎬⎢⎥−−⎪⎪⎪⎪⎢⎥⎩⎭⎩⎭⎢⎥⎢⎥⎣⎦(4)式中：E 为弹性模量；ν为泊松比；σ和τ分别为正应力和剪应力；εγ和分别为正应变和剪应变；下标n s 和分别表示垂直和平行裂缝的方向；G 为剪切模量；β为剪切抗力系数，01β≤≤；μ为随开裂应变变化的法向弹模折减系数．在有限元程序实现过程中，裂缝张开剪力传递系数取为0.5，裂缝闭合剪力传递系数取为0.9． 1.3 坝基动力极限承载能力评价模型1.3.1 基于DP 屈服准则的坝基破坏分析模型对于岩体材料，采用弹塑性模型，屈服准则采用DP (Drucker-Prager )屈服准则，即10f Ik α==(5)式中α和k 为材料参数，可以通过拟合摩尔-库伦准则得到．鉴于地应力场作用下的岩基为3向受压区，采用压力子午线下的参数，在π平面上选用库伦六边形外接圆形式的DP 准则，应用关联流动法则可以推导出k α⎧=⎪⎪⎨⎪=⎪⎩(6)式中：c 为黏聚力；ϕ为内摩擦角． 对于地基岩体，采用上述DP 屈服准则，考虑塑性屈服区扩展、联结，最后贯通，诱发其丧失承载能力．此时，便可以认为地基结构达到其承载极限，功能失效．失稳表达式为 ()0i i j m n A A φφφφ<= ∪∪∪(7)式中；φ为屈服面；i A 为局部屈服区φ∗组成的机构． 1.3.2 基于变形体突变的极限稳定性分析模型坝基渐进破坏过程的分析让人易于理解，但不同的本构模型和破坏准则以及人为主观因素的影响会导致判别结果的精度[5]，同时，考虑到地震荷载的往复特性和瞬态特性，加之材料非线性导致的应力转移现象，瞬时动力失稳尚不能说明坝体系统达到其极限抗震能力，因此，进一步从变形稳定的角度评判坝基是否达到其极限承载状态．在非线性有限元计算过程中，如果施加在结构上的荷载临近极限荷载时，坝体坝基系统的某些特征点位移会发生突变，可据此判别系统是否达到其极限承载力[5]．因此，本文以系统稳定状态突变导致非线性计算不收敛为判据，评估重力坝达到系统稳定状态极限时的抗震能力． 1.4 极限抗震能力综合研究方法综合上述系统破坏判据及对破坏状态的界定，本文从数值仿真试验的角度，建议从重力坝系统的损伤破坏等级、坝体混凝土开裂破坏模式、坝基岩体破坏扩展规律、系统稳定性突变等方面综合评价重力坝的极限抗震能力．(1) 指标E1：基于Ghanaat 坝体安全评价理念和弹塑性损伤理论，分别从线弹性和弹塑性损伤分析的角度，定量反映重力坝系统的损伤破坏等级，初步评估大坝的极限抗震能力．(2) 指标E2：从坝体坝基系统破坏模式角度评价大坝极限抗震能力．重力坝的开裂破坏是一个逐渐累积、渐进破坏的过程，当开裂破坏使大坝出现渗漏或出现坝块脱落时，大坝将丧失正常挡水功能，大坝达到极限抗震能力．同时考虑混凝土和坝基岩体材料的非线性，通过计算不同等级强震作用下大坝的损伤、开裂破坏，将坝体贯穿性裂缝的出现以及坝基2013年3月张社荣等：复杂层状岩基上重力坝极限抗震能力评估方法初探·205·损伤屈服贯通性破坏作为系统动力失效、达到其极限抗震能力的表征．(3) 指标E3：从坝基动力极限承载能力角度评价大坝极限抗震能力．大坝沿坝基面及深层薄弱层面的稳定性是整个工程抗震设计的关键和安全评价的重要方面，在地震往复荷载作用下，坝体地基系统存在失稳的可能性．随着地震动强度的增大，当某些特征点位移发生突变时，表征坝体-坝基系统稳定性状态突变，非线性计算不收敛，系统达到其极限抗震能力．最后，综合评价大坝极限抗震能力：基于上述分析结果，从系统功能失效风险的角度出发，综合坝体开裂破坏模式、坝基动力极限承载能力等方面，定义大坝的极限抗震能力为各个指标的最小值．2 评价模型的验证Koyna重力坝[15]作为少数几个在强震中破坏且有比较完整记录的重力坝之一，一直是混凝土坝动力分析的经典研究对象．Koyna混凝土重力坝坝高103m，底宽70.1m，顶宽14.78m(见图1)．1967年12月11日，该坝坝址区域遭受一次6.5级强烈地震作用，在91.7m高的库前水位、0.474g水平向与0.312g的竖向加速度峰值地震作用下，坝体发生开裂，并在裂缝处发生渗漏．已有许多学者对该坝的地震破坏过程进行了模拟分析．本文对Koyna重力坝的极限抗震能力进行分析．图1Koyna大坝示意Fig.1 Sketch of Koyna dam计算荷载考虑重力、静水压力、动水压力和Koyna地面实测地震波激励作用．坝体混凝土弹性模量[15]为31.0GPa，泊松比为0.2，密度为2500kg/m3，静态抗拉强度取 2.90MPa；动水压力根据Westergaard公式考虑；阻尼比为0.05．采用动力超载法，分析Koyna大坝的极限抗震能力．首先，按照指标E1对Koyna大坝进行极限抗震能力评价．由图2可知，当输入基岩水平峰值加速度低于0.400g时，大坝超应力累积持时均处于轻微破损区范围．当输入基岩峰值加速度大于等于0.450g 时，DCR＝1的超应力累积持时大于0.3s，且DCR＝2的超应力累积持时大于0s，大坝处于严重损伤破坏区范围．因此，基于性能的抗震评价模型，初步认为Koyna重力坝的极限抗震能力为0.450g左右．图2基于功能的Koyna大坝抗震破坏等级评价Fig.2 Performance-based damage evaluation of Koyna dam 进一步根据指标E2，分析坝体的潜在破坏模式．基于ANSYS弥散裂缝模型，图3给出了0.474g 水平向实测地震波作用下的Koyna大坝最终破坏模式与振动台模型试验[10]的结果．数值模拟得出的最终破坏形态与模型试验基本一致，也与Koyna坝实际地震开裂位置吻合，说明基于弥散裂缝模型的破坏过程分析可以有效表达动力工况下混凝土开裂行为，能够很好地模拟地震荷载作用下坝体的强震破坏模式，评价大坝的极限抗震能力．由于缺乏Koyna大坝坝基动力破坏等相关实测资料，在不考虑坝基极限承载力条件下，认为Koyna坝极限抗震能力为0.450g 左右，这与实际震害情况基本一致．（a）数值计算结果 （b）振动台模型试验结果图3数值计算破坏形态与模型试验结果比较Fig.3 Comparison of numerical failure mode andthe model test另外，笔者在文献[16]中，也从坝体开裂破坏、地震位移响应突变、塑性区贯通等方面对混凝土重力拱坝的极限抗震能力进行了初步探讨，这也佐证了文中·206· 天津大学学报(自然科学与工程技术版) 第46卷 第3期评价模型对整体大坝和相似工程极限抗震能力评价的适用性．3 工程实例分析阿海水电站属一等大(1)型工程．拦河坝为混凝土重力坝，最大坝高为138m ，水库正常蓄水位高程为1 504m ．建基岩体主要为板岩砂岩互层结构，承载能力较好；工程区地震地质背景复杂，坝址区设计烈度接近Ⅸ度，其100a 超越概率2%的基岩水平地震动峰值加速度达0.344g ．复杂层状岩基上的大坝抗震能力和抗震安全是本工程的关键技术问题． 3.1 数值模型及计算参数4号溢流坝段结构体形突异，深层滑移模式复杂，是工程薄弱环节，需进行极限抗震能力分析．数值模型(见图4)反映了坝段的主要结构特点、坝基岩体及主要地质构造分布，考虑板岩对大坝及地基应力的不利影响，通过适当的概化，对板岩单独划分单元，尽可能反映实际地质条件．坝体混凝土采用弹塑性模型，失效面采用William-Warnke 五参数强度准则；坝基岩体屈服面采用Mohr-Coulomb 屈服面的外接圆锥．计算时大坝水位为正常蓄水位，动水压力根据Westergaard 公式考虑；Rayleigh 阻尼因数根据线弹性分析得到的前两阶频率计算．模型材料参数见表1．建基面(AB /BC /CD 面)的内摩擦角和黏聚力分别为38.66°和0.61MPa ，深层AC 面的内摩擦角和黏聚力分别为41.67°和0.80MPa ．（a ）溢流坝段剖面及其滑动面示意 （b ）三维离散模型图4 溢流坝段(含滑动面)及三维离散模型Fig.4 Sketch of the overfall dam and its 3D discrete model 表1 坝体及地基材料力学参数Tab.1 Mechanical parameters of dam and foundation material材 料静态弹性 模量/MPa动态弹性 模量/MPa泊松比μ密度ρ/ (kg ·m -3)静态抗压强度/MPa动态抗压强度/MPa断裂能密度/ (N ·m ·m -3)内摩擦 角ϕ/(°)黏聚力c ′/MPa C15(碾压) 2.20×104 2.860×104 0.167 2,450 13.5 17.6 205 — — C20(碾压) 2.55×104 3.315×104 0.167 2,450 19.6 25.5 257 — — 坝体C30(常态) 3.00×1043.900×1040.1672,450 20.0 26.0 315 — — Ⅲ1 1.0×104 1.0×104 0.250 2,700 — —— 49.00 1.20Ⅲ2 8.0×1038.0×1030.270 2,700 —— — 46.40 1.00Ⅳ1 5.0×103 5.0×103 0.300 2,600 — — — 35.00 0.60 反倾节理AE 8.0×103 8.0×103 0.270 2,600 — — — 40.70 0.76 地基 岩体层间节理ED 8.0×103 8.0×103 0.2702,600 ———21.80 0.50注：混凝土动态弹性模量及抗压强度均较静态值提高30%，动态抗拉强度为动态抗压强度的10%；基岩的材料参数暂取静态值．3.2 地震动输入及研究方案坝区场地类别为Ⅰ类，100a 超越概率2%的基岩地震动水平向峰值加速度为0.344g ，100a 超越概率1%的基岩地震动水平向峰值加速度为0.415g ．基于2013年3月张社荣等：复杂层状岩基上重力坝极限抗震能力评估方法初探·207·有质量地基采用黏弹性人工边界输入，输入地震动按Koyna地震波(见图5)调幅至不同强震等级，同时考虑水平向和竖向地震，竖向地震输入加速度峰值取为水平向的2/3．为研究大坝的极限抗震能力，在设计地震(PGA 为0.344g)和校核地震(PGA为0.415g)的地震响应基础上，采用地震动超载法，又分别计算了基岩水平峰值加速度分别为0.450g、0.500g、0.550g、0.600g等多个条件下的大坝动力响应．依据前述不同方面的指标，对大坝极限抗震能力进行综合评价．图5Koyna波加速度时程归一化曲线Fig.5 Normalized curves of Koyna wave groundacceleration3.3 从结构损伤破坏等级分析极限抗震能力通过对各方案下坝体特征点动应力时程曲线对比，中墩与溢流堰顶的连接区域动力响应较大，是结构动力响应的控制性区域，该部位的破坏也将直接导致大坝蓄水功能的失效．因此将该部位的需求能力比和超应力累积持时作为控制指标．该区域混凝土动态名义抗拉强度达4.789MPa(DCR＝2的分界点)，图6给出了PGA为0.450g和0.500g时该部位的竖向正应力时程，进而得出需求能力比和超应力累积持时指标．根据文献[7]的计算方法，基于线弹性模型的不同等级地震作用下的大坝损伤破坏等级如图7所示．当地震动峰值加速度增大至0.450g时，DCR＝1的超应力累积持时达到了0.29s(小于0.30s)，DCR＝2超应力累积持时为0s，仍处于中低等损伤破坏区．继续增大至0.500g时，DCR＝1的超应力累积持时大于0.30s，且DCR＝2的超应力累积持时大于0s，坝体已处于严重损伤破坏范围内．进一步，基于ABAQUS中的塑性损伤模型，对坝体动力损伤情况进行分析．选取大坝典型部位的单元，绘制结构动力非线性损伤演化情况，如图8所示．由此可见，大坝首先在下游溢流导墙与闸墩交接（a）PGA为0.450g（b）PGA为0.500g图6竖向最大拉应力Fig.6Maximum vertical tensile stress图7基于功能的阿海大坝抗震破坏等级评价Fig.7Performance-based seismic damage evaluation of Ahaigravity dam图8PGA＝0.500g时坝体关键部位的损伤演化过程Fig.8 Damage evolution vs time when PGA＝0.500g处(E847)损伤开裂，随后是坝踵单元(E2347)，接着是溢流面下游折坡(E10206)处，最后迎水面1469m 高程处上游坝面(E490)也逐渐损伤破坏直至出现裂缝．同时，在地震作用的前4s内，坝体损伤较小且发展缓慢，地震峰值时刻后迅速发展，局部开裂已经完成．此后损伤变量接近稳定并维持在一定程度．·208·天津大学学报(自然科学与工程技术版)第46卷 第3期将中墩与溢流堰顶的连接区破损视为导致大坝蓄水功能失效的关键区域，结合基于性能的抗震评价模型，大坝的极限抗震能力指标E1可初步评价为0.550g，此时大坝已处于严重损伤破坏，需进行大坝非线性分析，从坝体-坝基的潜在失效模式，进一步评估大坝-地基系统的极限抗震能力．3.4 从坝体开裂破坏模式分析极限抗震能力将混凝土作为连续介质，基于ANSYS软件中的弥散裂缝模型对坝体在不同强度地震动作用下的裂缝分布及开裂状态进行分析．如图9所示，当PGA＝0.415g时，坝踵回填混凝土基本裂穿，闸墩与溢流堰顶连接处出现水平张开裂缝，坝段主体结构受地震影响不大；PGA＝0.450g的情况较PGA＝0.415g时并无明显恶化，当地震峰值增大至0.500g时，闸墩两侧与溢流堰顶连接处出现几乎贯通堰顶的水平向闭合裂缝，建基面也出现多处裂缝，此时中墩结构破坏导致挡水功能失效的概率明显增大；进一步，当PGA＝0.550g时，溢流堰顶已经出现了贯通结构上下游方向的贯穿性裂缝．继续增大至0.600g时，结构非线性计算不收敛．另外，数值仿真得到的坝体破坏模式与文献[3]中类似工程溢流坝段物理模型试验所得的大坝动力破坏形态相似，这也佐证了数值分析中坝体开裂破坏模型的适用性和正确性．（a）PGA＝0.415g（b）PGA＝0.450g （c）PGA＝0.500g（d）PGA＝0.550g图9不同强度地震作用下坝体裂缝分布规律Fig.9 Crack distribution under different earthquakes从坝体开裂破坏模式指标E2看，大坝的极限抗震能力为0.550g～0.600g．3.5 从坝基极限承载能力分析极限抗震能力3.5.1坝基渐进破坏规律分析坝基地质特点决定地基塑性变形发展规律，地基塑性屈服可能会造成系统突变失稳，进而影响坝体系统的极限抗震能力．基于ABAQUS分析平台，坝基岩体采用外接圆锥DP屈服准则，不同地震加速度下坝基塑性区分布如图10所示．PGA＝0.415g时，坝基塑性区出现在坝踵区域且沿反倾节理方向延伸；当PGA＝0.450g时，坝基塑性区沿坝基面向下游穿过帷幕范围，同时沿反倾节理进一步延伸，另外，坝趾处层间节理区域也开始出现塑性区；PGA＝0.500g时，坝踵塑性区向深部扩展，坝趾处塑性区进一步沿层面节理延伸；PGA增大至0.550g时，坝踵处塑性区范围较PGA＝0.500g时并无显著增大，但塑性变形量值却显著增大．地震峰值加速度至0.600g，非线性计算不收敛．（a）PGA＝0.415g（b）PGA＝0.450g（c）PGA＝0.500g（d）PGA＝0.550g图10不同强度地震作用下坝基塑性区分布(Avg：100%，单位：m)Fig.10 Plastic zone distribution under different earthquakes(Avg：100%，unit：m)3.5.2大坝抗滑稳定安全性分析同时考虑坝体坝基的材料非线性，按照现行抗震规范的承载能力极限状态表达式，结合基于有限元的刚体极限平衡法、矢量和积分的概念，建立动力抗滑稳定安全系数的计算公式为。