基于AN和SYS的阿海水电站静动力分析.docVIP

基于ANSYS的阿海水电站混凝土碾压重力坝静动力分析 邹小红1，于跃2，张彩秀2 1河海大学 水利水电学院，江苏 南京 （210098） 2昆明勘测设计研究院，云南 昆明 （650051） E-mail：、2 计算内容及方法 2.1 计算内容 通过阿海重力坝方案进行线性静动力有限元分析，，，，；，DL5073-1997《水工建筑物抗震设计规范》中所规定的振型分解反应谱法按弹性材料进行计算[5]。 2.2 计算工况 计算考虑正常蓄水位工况、设计洪水位工况、校核洪水位工况、施工期工况和地震工况，各工况荷载组合见表1[3,4]。 表1 计算工况及荷载组合 计算工况 坝体自重 静水压力(m) 扬压力 泥沙压力 浪压力 上游水位 下游水位 厂房下游水位 基本组合 正常蓄水位 √ 1504.00 1411.866（1.5台机） 1411.866（1.5台机） √ √ √ 设计洪水位(P=0.2％ 重力坝的地震分析一般有动力法和静力法，，，，:首先是建立模型后求得模态解，；，SRSS”方式，10阶振型。坝体在地震作用下发生振动时，库水在上游坝面将产生动水压力，这相当于增加了坝体运动的惯性力，从而加大了地震荷载。所以需要考虑动水压力，， 式中：Pw(h)为地震动水压力；ah为地震加速度水平向代表值；H0为水深；阻尼比采用Rayleigh比例阻尼假定，阻尼系数ξ=0.05。 本次计算中同时计入水平向和竖向地震作用。计算参数如下：水平向设计地震加速度代表值ah＝0.344g，竖向设计地震加速度代表值av＝2*ah/3；特征周期Tg＝0.2（Ⅰ级场地）；设计反应谱最大值的代表值βmax＝2.0。动弹性模量为静弹性模量的1.3倍。根据《水工建筑物设计规范（DL5073－2000）》，总的地震作用效应为竖向地震作用效应乘以0.5的遇合系数后与水平向地震作用效应直接相加。 3 静力计算结果分析 3.1 静态位移场 坝体顺河向位移：在正常蓄水位工况、设计洪水位工况及校核洪水位工况中，坝体在上游水压力的作用下倾向下游，各坝段的最大顺河向位移均在坝顶下游侧角点，且均发生在校核洪水位工况，厂房坝段位移值最大，为5.45cm。在施工期工况中，坝体均倾向上游，各坝段的最大位移均出现在坝顶上游侧角点，左冲沙底孔坝段最大，为-1.33cm。坝体竖向位移：在正常蓄水位工况、设计洪水位工况及校核洪水位工况中，除了溢流坝段竖直向最大沉降出现在校核水位工况外，其它各坝段均在正常蓄水位工况下沉降较大，厂房坝段最大，为-3.73cm。在施工期工况中，各坝段的最大沉降均出现在坝顶上游侧角点，厂房坝段最大，为-3.80cm。 3.2 静态应力场 竖向正应力：在正常蓄水位工况、设计洪水位工况及校核洪水位工况中，各坝段坝踵均有拉应力集中现在，厂房坝段拉应力最大，为1.73MPa，拉应力宽度为2.9％，且不超过上游帷幕；溢流坝段拉应力宽度最大，为4.7％，但未超过上游帷幕。在上述三个工况中，各坝段坝趾均出现压应力集中现在，但在应力集中区域之外，坝基压应力均小于其承载力。厂房坝段竖向应力最大，发生在校核洪水位工况，为-8.46MPa，在距坝趾4m处折减为-3.05MPa。在施工期工况中，坝基均为压应力，且坝踵比坝趾大，厂房坝段坝踵处最大，为-4.77MPa。垂直岩层方向压应力：在正常蓄水位工况、设计洪水位工况、校核洪水位工况及施工期工况中，各坝段在坝趾附近压应力较集中，但在应力集中区域之外，坝基压应力均小于其承载力。厂房坝段最大，发生在校核洪水位工况，为-12.00MPa，距厂坝分缝处约4m处折减为-4.44MPa。 4 动力计算结果分析 4.1 模态分析 用于确定结构或机器部件的振动特性(固有频率和振型)。一般来说，nn个自振频率，n个主振型。但这n个主振型中一般只有前几个是起决定作用的。为了计算准确，10阶自振频率，，，，，，各坝段，，，，，，，0.07，，(2)各坝段在满库时的第一阶频率分别为2.27Hz，2.42Hz，1.76Hz和1.50Hz，，，，，(3)总体上来讲，，，，，， 参考文献 [1] 荆旭春，等.宝泉电站下水库重力坝地震动力响应分析[J].人民黄河，200，(5)：56－60. [2] 罗业辉，赵海涛, 邓仕涛.应用ANSYS软件进行碾压混凝土重力坝非线性有限元静力和动力[J].西北水电，200，()：18－22. [3] 中华人民共和国电力工业部．DL5108-1999[S].北京：中国水利水电出版社，1999．[4] 中华人民共和国电力工业部．DL5077-1997[S].北京：中国水利水电出版，1997．[5] 中华人民共和国电力工业部．DL5073-1997[S].北京：中国水利水电出版，1997．