两种库水附加质量模型的重力坝动力响应研究_黄耀英.docVIP

——————————————————————————————————————————————— 两种库水附加质量模型的重力坝动力响应研究_黄耀英 第40卷第7期 ?? 2009年4月 ????文章编号:1001-4179(2009)07-0064-03 ?? 人??民??长??江 ?? Yangtze??River ?? Vo.l40,No.7Apr.,2009 两种库水附加质量模型的重力坝动力响应研究 黄耀英??孙大伟??田??斌 (三峡大学土木水电学院,湖北宜昌443002) 摘要:对两种不可压缩性库水附加质量模型的混凝土重力坝动力响应进行了研究,得到结论:??相对有限元模 型计算结果比较,韦斯特伽德模型进行折半处理后计算得到的大坝自振频率较直接使用韦斯特伽德模型计算的精度高;??当水体区域向上游取3倍坝高可以得到比较稳定的自振频率;??由于大坝是弹性体,上游面坝体结点的等效结点荷载不仅受该结点动水附加质量本身的影响,同时受到该结点周边结点的动水附加质量的影响。 关??键??词:库水附加质量;地震响应;混凝土重力坝;抗震分析中图分类号:TV642.3??????文献标识码:A ????坝水耦合作用是混凝土大坝抗震分析中的一个重要课题。自从Westergaard提出垂直刚性坝面的库水动水压力理论以来, 这一问题已经历了70多年的深入研究,但仍有比较多的问题不甚清楚。目前工程上常假定库水为不可压缩流体,此时,库水的动力作用就相当于在大坝系统的质量矩阵中加上一个附加质量矩阵,关于附加质量矩阵的计算,工程上常用的有两种方法,即按有限元计算的方法和按广义韦斯特伽德公式计算的方法。陈厚群等[4]通过白山拱坝模型实测以及三向电拟试验求得振型动水压力,得到结论:流固耦合的有限单元数学模型给出的结果与试验结果很符合,而广义韦斯特伽德公式给出的结果偏大,建议对广义韦斯特伽德公式进行折半修正。李瓒等[5]对二滩拱坝采用有限元模型和韦斯特伽德模型进行动力分析计算的结果也表明韦斯特伽德模型夸大了动水压力的影响。李瓒等[5]认为应取3倍以上坝高的水库域。文献[6,7]采用有限元模型计算动水附加质量时,对上游水体取3倍坝高进行有限单元离散。但对比分析上述两种计算方法对大坝自振频率的影响以及对大坝地震响应的影响,有关文献报导的仍显不足。因此,本文对上述两种不可压缩性库水附加质量模型的混凝土重力坝动力响应进行研究。 [1~7] 在大坝系统的质量矩阵中加上一个附加质量矩阵[Mp],此时,有限元动力控制方程为: [M*]{y??}+[C]{y??}+[K]{y}=-[M*]{??yg}(1)式中{y}、{y??}、{??y}分别为相对位移、速度和加速度;{??yg}为地面加速度;[M*]=[M]+[Mp];[M]、[K]、[C]分别为整体质量矩阵、劲度矩阵和阻尼矩阵;[Mp]为附加质量矩阵。 关于附加质量矩阵的计算,工程上常用的有两种方法,即按有限元模型计算的方法和按广义韦斯特伽德模型计算的方法。当采用广义韦斯特伽德模型计算动水压力时,作用在上游坝面上的动水压力及附加质量矩阵为: Pi=??i??yni7??i=??iii8 [Mpi]=??iAi(??)T??i i [Mp]= (2)(3)(4)(5) ?? [Mpi] 式中Pi为上游坝面的动水压力,以压为正;??为库水的质量密度;Hi为包含节点i的铅直截面上的水深;Zi为包含节点i的铅直截面上从坝基算起到节点i的高度;??yni为在节点i的坝面法向加速度;Ai为与节点i有关的辅助面积;??i为在节点i的坝面法线方向的余弦向量;[Mpi]为节点i的库水附加质量,可直接将其组合到大坝系统的质量矩阵中。 当采用有限元模型计算动水的附加质量时,需要对水体进行有限单元离散,此时,附加质量矩阵为: [Mp]=??[S]T[H]-1[S](6)式中??为水的密度;[H]为动水压力劲度矩阵;[S]为影响矩 阵。 本文采用振型叠加法或纽马克法(NewMark)求解动力控制 1??基本原理 当大坝受到地震运动激励时,坝基随着地震产生刚性位移,同时大坝还发生弹性位移,因而产生库水对坝体的动水压力 (包括刚体位移加速度和弹性位移加速度产生的动水压力),而动水压力又影响着坝的弹性位移加速度,所以它们相互影响着,这就需要把坝体和水体作为一个整体进行动力反应分析。 当假设库水为不可压缩流体时,库水的动力作用就相当于收稿日期:2008-