第40卷第7期2009年4月 人 民 长 江Yangtze River Vol.40,No.7Apr.,2009收稿日期:2008-08-20作者简介:黄耀英,男,三峡大学土木水电学院,讲师,博士。 文章编号:1001-4179(2009)07-0064-03两种库水附加质量模型的重力坝动力响应研究黄耀英 孙大伟 田 斌(三峡大学土木水电学院,湖北宜昌443002)摘要:对两种不可压缩性库水附加质量模型的混凝土重力坝动力响应进行了研究,得到结论:①相对有限元模型计算结果比较,韦斯特伽德模型进行折半处理后计算得到的大坝自振频率较直接使用韦斯特伽德模型计算的精度高;②当水体区域向上游取3倍坝高可以得到比较稳定的自振频率;③由于大坝是弹性体,上游面坝体结点的等效结点荷载不仅受该结点动水附加质量本身的影响,同时受到该结点周边结点的动水附加质量的影响。关 键 词:库水附加质量;地震响应;混凝土重力坝;抗震分析中图分类号:TV642.3 文献标识码:A 坝水耦合作用是混凝土大坝抗震分析中的一个重要课题。自从Westergaard提出垂直刚性坝面的库水动水压力理论以来,这一问题已经历了70多年的深入研究,但仍有比较多的问题不甚清楚。目前工程上常假定库水为不可压缩流体,此时,库水的动力作用就相当于在大坝系统的质量矩阵中加上一个附加质量矩阵,关于附加质量矩阵的计算,工程上常用的有两种方法[1~7],即按有限元计算的方法和按广义韦斯特伽德公式计算的方法。陈厚群等[4]通过白山拱坝模型实测以及三向电拟试验求得振型动水压力,得到结论:流固耦合的有限单元数学模型给出的结果与试验结果很符合,而广义韦斯特伽德公式给出的结果偏大,建议对广义韦斯特伽德公式进行折半修正。李瓒等[5]对二滩拱坝采用有限元模型和韦斯特伽德模型进行动力分析计算的结果也表明韦斯特伽德模型夸大了动水压力的影响。李瓒等[5]认为应取3倍以上坝高的水库域。文献[6,7]采用有限元模型计算动水附加质量时,对上游水体取3倍坝高进行有限单元离散。但对比分析上述两种计算方法对大坝自振频率的影响以及对大坝地震响应的影响,有关文献报导的仍显不足。因此,本文对上述两种不可压缩性库水附加质量模型的混凝土重力坝动力响应进行研究。1 基本原理当大坝受到地震运动激励时,坝基随着地震产生刚性位移,同时大坝还发生弹性位移,因而产生库水对坝体的动水压力(包括刚体位移加速度和弹性位移加速度产生的动水压力),而动水压力又影响着坝的弹性位移加速度,所以它们相互影响着,这就需要把坝体和水体作为一个整体进行动力反应分析。当假设库水为不可压缩流体时,库水的动力作用就相当于在大坝系统的质量矩阵中加上一个附加质量矩阵[Mp],此时,有限元动力控制方程为:[M]{¨y}+[C]{y}+[K]{y}=-[M]{¨yg}(1)式中{y}、{y}、{¨y}分别为相对位移、速度和加速度;{¨yg}为地面加速度;[M]=[M]+[Mp];[M]、[K]、[C]分别为整体质量矩阵、劲度矩阵和阻尼矩阵;[Mp]为附加质量矩阵。关于附加质量矩阵的计算,工程上常用的有两种方法,即按有限元模型计算的方法和按广义韦斯特伽德模型计算的方法。当采用广义韦斯特伽德模型计算动水压力时,作用在上游坝面上的动水压力及附加质量矩阵为:Pi=αi¨yni(2)αi=78ρHi(Hi-Zi槡)(3)[Mpi]=αiAi(λi)Tλi(4)[Mp]=[Mpi](5)式中Pi为上游坝面的动水压力,以压为正;ρ为库水的质量密度;Hi为包含节点i的铅直截面上的水深;Zi为包含节点i的铅直截面上从坝基算起到节点i的高度;¨yni为在节点i的坝面法向加速度;Ai为与节点i有关的辅助面积;λi为在节点i的坝面法线方向的余弦向量;[Mpi]为节点i的库水附加质量,可直接将其组合到大坝系统的质量矩阵中。当采用有限元模型计算动水的附加质量时,需要对水体进行有限单元离散,此时,附加质量矩阵为:[Mp]=ρ[S]T[H]-1[S](6)式中ρ为水的密度;[H]为动水压力劲度矩阵;[S]为影响矩阵。本文采用振型叠加法或纽马克法(NewMark)求解动力控制 第7期 黄耀英等:两种库水附加质量模型的重力坝动力响应研究方程,时间步长取为0.01s,计算表明,做线弹性分析时,振型叠加法和纽马克法的计算结果基本一致。其中,纽马克法的阻尼矩阵采用瑞利(Rayleigh)比例阻尼,即阻尼矩阵假设为质量矩阵和劲度矩阵的线性组合,其表达式为:[C]=α0[M]+α1[K](7)式中[C]、[M]和[K]分别为结构阻尼矩阵、质量矩阵和劲度矩阵;α0、α1为系数。首先求出结构的工程频率f1和f2,然后得到相应的圆频率ω1和ω2,再结合水工建筑物抗震设计规范假定的阻尼比λ1和λ2,由式(8)可计算得到系数α0和α1。例如:采用韦斯特伽德模型进行柯依那(Koyna)重力坝满库下大坝的地震响应时,本文计算的α0=1.15103、α1=0.001732。λ1=α02ω1+α1ω12λ2=α02ω2+α1ω2}2(8) 当采用振型叠加法时,则直接按规范给定前若干阶振型的阻尼比λ2(i=1,2,…,n)。2 算例分析2.1 基本资料印度柯依那(Koyna)重力坝[3],最大坝高103m,采用有限元法对该坝及水体区域进行动力响应分析,材料参数取为:动弹性模量E=31.5GPa,容重γ=26.5kN/m3,各阶振型的阻尼比均取为0.05。按平面应变问题考虑,利用衰减三角级数叠加模型人工生成地震波,基底顺河向均匀输入地震波,地震历程10s。由于本文重点不是进行大坝地震破坏分析,所以峰值加速度调至0.1g,如图1所示。图1 输入地震波2.2 计算内容(1)计算大坝在空库、满库下的动力特性。其中,满库计算时分别采用了按有限元法计算与按广义韦斯特伽德(韦氏)公式计算,同时给出其它文献的计算值进行比较。其中,按韦氏公式计算时,计算了折半处理和不折半处理下大坝满库自振频率;按有限元法计算时,对水体区域所取范围进行了敏感性分析。(2)分别进行了空库和满库下大坝的动力响应时程分析。其中,按有限元法计算动水压力时,水体区域向上游取3倍坝高。2.3 计算结果2.3.1 自振频率不同工况下大坝自振频率见表1。由表1数据可见:(1)本文计算的大坝自振频率与文献[3]计算的结果相近;(2)动水附加质量提高了大坝的整体质量,使得大坝满库时的自振频率较空库时小;(3)由于韦斯特伽德模型夸大了动水压力作用,由它计算出的库水附加质量偏大,所以韦斯特伽德模型进行折半计算后得到的大坝自振频率较直接使用韦斯特伽德模型计算的精度高(相对有限元法计算结果比较);(4)对水体区域范围进行敏感性分析表明,当水体区域向上游取3倍坝高可以得到比较稳定的自振频率结果。表1 混凝土重力坝的自振频率比较Hz工况阶 次12345空库3.2008.43911.11016.38024.780文献[3]空库3.1658.33310.98916.12925.000文献[3]满库2.7327.24610.87014.93023.260韦氏模型满库2.5276.66010.96012.72018.740韦氏模型折半满库2.8067.40311.01014.30021.690有限元模型满库水体区域取1倍坝高2.7107.22710.98015.06023.780 水体区域取2倍坝高2.7257.28410.99015.09023.810 水体区域取3倍坝高2.7267.28610.99015.09023.810 水体区域取4倍坝高2.7267.28610.99015.09023.810 水体区域取5倍坝高2.7267.28610.99015.09023.8102.3.2 地震响应不同工况下大坝的地震响应见表2。图2~5给出了空库和有限元法满库时坝顶处的顺河向位移、速度和加速度以及坝颈处垂直应力时程。表2中位移以向下游为正,应力以拉为正,位移、速度和加速度均为相对值。表2 大坝地震响应工况最大顺河向位移/cm最小顺河向位移/cm最大顺河向速度/(m·s-1)最小顺河向速度/(m·s-1)最大顺河向加速度/(m·s-2)最小顺河向加速度/(m·s-2)最大垂直向应力/MPa最小垂直向应力/MPa空库 1.181-1.0130.231-0.2364.904-4.5851.622-1.741韦氏模型满库 2.201-2.1930.369-0.3385.952-6.6363.144-3.521韦氏模型折半满库1.601-1.3430.220-0.2284.196-4.6512.111-2.464有限元模型满库 2.609-2.4990.401-0.4077.582-7.4683.681-3.829由表2数据可见:(1)由于库水附加质量的影响,大坝满库的地震响应较空库时大;例如:空库时计算的坝顶最大顺河向位移为1.181cm,韦斯特伽德模型满库计算的坝顶最大顺河向位移为2.201cm,增大86%,有限元模型满库计算的坝顶最大顺河向位移2.609cm,增大121%。(2)韦斯特伽德模型满库的地震响应较韦斯特伽德模型折半满库时大;例如:韦斯特伽德模型满库计算的坝颈最大拉应力为3.144MPa,折半满库计算的坝颈最大拉应力为2.111MPa。(3)有限元模型满库的地震响应较韦斯特伽德模型满库时大,原因为韦斯特伽德模型的库水附加质量是假设大坝为刚性获得的,而实际大坝是弹性体,上游面坝体结点的等效结点荷载不仅受该结点库水附加质量本身的影响,同时受到该结点周边结点的库水附加质量的影响。56 人 民 长 江2009年 图2 坝顶顺河向位移响应图3 坝顶顺河向速度响应图4 坝顶顺河向加速度响应图5 坝颈垂直向应力响应3 结语对两种不可压缩性库水附加质量模型的混凝土重力坝动力响应进行了研究,得到结论:(1)由于韦斯特伽德模型夸大了动水压力的作用,由它计算出的库水附加质量偏大,韦斯特伽德模型进行折半计算后得到的大坝自振频率较直接使用韦斯特伽德模型计算的精度高(相对有限元模型计算结果)。(2)对水体区域范围进行敏感性分析表明,当水体区域向上游取3倍坝高可以得到比较稳定的自振频率结果。(3)韦斯特伽德模型满库的地震响应较韦斯特伽德模型折半满库时大,而有限元模型满库的地震响应较韦斯特伽德模型满库时大。原因为韦斯特伽德模型的库水附加质量是假设大坝为刚性获得的,而实际大坝是弹性体,上游面坝体结点的等效结点荷载不仅受该结点库水附加质量本身的影响,同时受到该结点周边结点的库水附加质量的影响。(4)由于韦斯特伽德模型夸大了动水压力的作用,虽然可对该模型进行修正处理,但具有不确定性;另外,韦斯特伽德模型是假设大坝为刚性体获得的,其计算的上游面坝体结点的等效结点荷载难以合理考虑周边结点的库水附加质量的影响,所以建议采用按有限元模型计算库水压力的附加质量。参考文献:[1] 中华人民共和国水利部.水工建筑物抗震设计规范.北京:水利电力出版社,1997.[2] 华东水利学院.弹性力学问题的有限单元法.北京:水利电力出版社,1978.[3] 赵光恒.结构动力学.北京:中国水利水电出版社,1996.[4] 陈厚群,侯顺载,杨大伟.地震条件下拱坝库水相互作用的试验研究.水利学报,1989,(7):29-39.[5] 李瓒,陈兴华,郑建波等.混凝土拱坝设计.北京:中国电力出版社,2000.[6] 谢能刚.高地震区高拱坝体型的多目标动力优化设计方法的研究.南京:河海大学博士学位论文,1999.[7] 黄耀英.混凝土坝接缝动力行为研究.南京:河海大学硕士学位论文,2003.(编辑:赵凤超)66