丰满水电站坝体电梯井有限元分析计算研究

2023-03-23杨苗苗

东北水利水电 2023年3期

谭欣，杨苗苗

（中水东北勘测设计研究有限责任公司，吉林长春 130021）

1 工程概况

丰满水电站全面治理（重建）工程是按恢复电站原任务和功能，在原大坝下游120.00 m 处新建一座大坝并利用原有三期电站。治理工程实施后，水电站以发电为主，兼顾防洪、城市及工业供水、灌溉、生态环境保护，并具有旅游、水产养殖等效益。电站枢纽工程主要建筑物由碾压混凝土重力坝、坝身泄洪系统、左岸泄洪兼导流洞、坝后式引水发电系统、GIS 开关站、过鱼设施及原有三期电站组成。大坝设计地震采用基准期100 年内超越概率2%，相应地震水平加速度为131 gal。

本文计算选取的电梯井位于左岸8 号挡水坝段，电梯井底高程取194.00 m，顶高程269.50 m，总高75.50 m，高程194.00～240.00 m 部分电梯井位于坝体内部，高程240.00 m 以上电梯井位于坝体外部下游侧。电梯井断面尺寸为6.00 m×6.40 m。

2 电梯井有限元计算

2.1 抗震计算方法

采用阵型分解反应谱法对电梯井进行抗震计算，采用NB 35047-2015《水电工程水工建筑物抗震设计规范》标准设计反应谱。首先进行结构静力计算，采用振型分解反应谱法分别对顺水流方向、垂直水流方向及竖直方向地震荷载下的结构响应进行计算，然后采用SRSS 法[1]对各方向结构响应进行组合作为总的地震荷载响应。

2.2 混凝土配筋方法

大体积混凝土部分依据DL/T 5057-2009《水工混凝土结构设计规范》附录D（非杆件体系钢筋混凝土结构的线弹性应力图形法配筋计算原则）进行配筋验算。整理结构混凝土部分各截面内力，然后按结构力学方法配筋。

2.3 计算参数

1）地震动参数

丰满电站场址100年超越概率2%和100年超越概率1%的地震水平加速度分别为131 gal和172 gal，相应的地震动反应谱特征周期均为0.20 s。竖向设计地震加速度的代表值取水平向设计地震加速度代表值的2/3，动力放大系数最大值采用2.0。

2）材料参数

坝基岩体及坝体混凝土材料参数见表1。

表1 材料参数表

2.4 计算过程

基于弹性实体单元理论，采用通用的有限元计算软件ANSYS 进行计算，建立大坝、基岩及附属电梯井的三维实体模型，采用Solid45 单元模拟大坝混凝土、基岩、电梯井混凝土。计算工况采用地震工况。基岩侧面施加法向约束，坝体碾压混凝土横河向施加侧向约束，底面基岩全约束。荷载主要有自重、上下游水压力、扬压力、浪压力、地震荷载[2]。考虑坝体和库水的相互作用，地震时坝上游面动水压力折算为单位地震加速度相应的坝面附加质量[3]。有限元计算模型坝基在上下游及深度方向上各延伸150.00 m，横水流方向计算边界取8 号坝段部分。根据有限元计算经验，综合考虑计算精度与效率，在此次计算重点关注的电梯井附近区域网格大小采用0.50 m，远离电梯井部分的坝体混凝土采用1.50 m，基岩部分采用10.00 m。

3 计算结果及分析

电梯井各高程应力变形及内力计算结果不同，高程250.00～269.50 m 网格划分比较规则，本文选取250.00 m 高程计算结果进行分析，通过提取薄壁结构部分应力并积分最终得到各个剖面的内力。沿电梯井周边各剖面布置见图1。

图1 高程250.00 m 计算截面及路径示意图

根据提取的应力沿路径积分最终得出的各截面内力汇总见表2。由表2 可知，250.00 m 高程电梯井周边轴力、弯矩、剪力最终计算结果均较小，最大弯矩值仅为12.92 kN·m，结构配筋只需按规范规定的最小配筋率计算即可满足设计标准。通过采用有限元程序计算电梯井周围内力及配筋，节省了工程材料用量，为工程设计提供了有效参考。