某水库除险加固中双曲拱坝应力及稳定计算分析

吴泽华，刘鹏程

（1.上海千年城市规划工程设计股份有限公司淮海设计院，江苏宿迁 223800；2.徐州市南水北调工程管理中心，江苏徐州 221200）

1 工程概况

某水库大坝为细石混凝土砌块石双曲拱坝，坝顶高程151.80 m，防浪墙顶高程152.70 m，最大坝高42.0 m，坝顶宽3.0 m（含防浪墙），坝顶弧长191.0 m。冲砂孔埋在坝内，内径为0.6 m，采用钢管衬套，放水孔中心高程为115.00 m，底高程为114.70 m，坝后接闸阀。溢洪道位于大坝中部，为钢筋混凝土实用堰，堰顶高程150.00 m，溢洪道净宽73.0 m，挑流消能，鼻坎高程147.14 m，反弧半径R=4.0 m，挑射角15°。引水隧洞布置在大坝左岸，采用圆形有压洞设计，内径1.8 m，长1.78 km，进洞口底高程为122.50 m，出洞口底高程为98.00 m，隧洞距离出口50.0 m设置一直径为4.0 m圆筒式调压井。水库正常水位150.00 m，下游无水，设计洪水位152.05 m，校核洪水位152.62 m。下文通过拱坝应力计算、重力墩稳定计算和坝肩稳定计算，进一步研究坝体结构安全，为水库除险加固设计施工提供指导。

2 拱坝应力计算

拱坝应力采用北京水利规划设计院开发的GADAP25拱坝软件进行计算，其力学原理基础是拱坝试荷载法，计算方法采用反力参数法[1]。

2.1 计算荷载组合

计算荷载组合分为4种工况[2]：1）在正常蓄水位情况下均匀水砂压力+温度变化（温降）；2）在设计洪水位情况下均匀水砂压力+温度变化（温升）；3）在校核洪水位情况下均匀水砂压力+温度变化（温升）；4）在死水位情况下均匀水砂压力+温度变化（温升）+地震力。前2种工况为基本组合，后2种为特殊组合。

2.2 计算参数

1）拱坝坝体

根据地质报告资料，拱坝砌体饱和容重为2.3 t/m3，弹性模量取Er=10 GPa（10×105t/m2），砌体温度线性膨胀系数采用0.000 007/℃。

2）坝基岩体

变形模量弹模E0=23 GPa，泊松比μ=0.25。

3）温度变化荷载

按SL 25-2006《砌石坝设计规范》中的简化计算法计算。

库水位以上：

库水位以下：

当坝体厚度L＜10.0 m时：ρ1=e-0.000671L3；

式中：Axi——最热与最冷月多年平均气温差值之半；Asho——库表面水温年变幅；y——水深。根据气象情况，同时参考同类水库及初设情况取值，Axi和Asho均取8℃。

该工程拱坝根据不同的荷载组合，温度计算参数Tm2及Td2，这里仅列出2种荷载工况，见表1。

表1 拱坝温度荷载电算参数

4）地震荷载

根据国家地震局、建设部发布的《中国地震烈度区划图（1990）》，云霄县地震基本烈度为Ⅶ度，相应加速度为0.10g，计算时先按45°增量从0变化至360°下每个方向的地震应力，再进行组合，找出其中最大的应力值。

2.3 计算结果及分析

根据上述计算方法及计算参数，此工程拱坝应力计算成果如表2，3所示，应力成果以压为正、拉为负[3]。

从表2，3的结果可以看出：在校核工况组合下，出现的拉应力最大值为1.923 MPa，在地震荷载组合下，出现的拉应力最大值为0.887 MPa；在校核工况组合下，出现的压应力最大值为3.692 MPa，在地震荷载组合下，出现的压应力最大值为1.950 MPa。

表2 拱坝左岸上下游面应力计算成果MPa

顶拱端上下游坝面主应力值见表4。

表4 顶拱端上下游坝面主应力值成果表MPa

根据以上计算结果可知，大坝压应力和拉应力分布规律仅在112.00 m高程个别应力点略大于规范容许值，大坝应力计算基本满足要求。

表3 拱坝右岸上下游面应力计算成果MPa

3 重力墩稳定验算

3.1 基本情况

拱坝右岸重力墩基础位于强风化花岗岩上，墩底高程142.30 m，墩顶高程152.80 m，墩顶长×宽为12.20 m×13.50 m，墩底长×宽为15.35 m×13.50 m。重力墩砌体容重为2.35 t/m3，根据地质勘查报告，夯填粘土的容重为1.6 t/m3，粘聚力c=2.2 kN/m2，内摩擦角φ=15°，砌体与基础的摩擦系数取0.50。拱坝重力墩稳定计算荷载组合分基本组合和特殊组合[4]：基本组合，正常蓄水位为自重+静水压力+扬压力+泥砂压力+浪压力+土压力+拱端推力，设计洪水位为自重+静水压力+扬压力+泥砂压力+浪压力+土压力+动水压力+拱端推力；特殊组合，校核蓄水位为自重+静水压力+扬压力+泥砂压力+浪压力+土压力+拱端推力，正常蓄水位为自重+静水压力+扬压力+泥砂压力+浪压力+土压力+地震荷载+拱端推力。

3.2 计算结果及分析

3.2.1 拱端推力荷载计算

拱的推力按照GADAP25程序计算，拱端推力荷载计算成果见表5。

表5 拱端推力荷载计算成果表

3.2.2 地震荷载计算

地震荷载采用拟静力法计算，按照SL 203-97《水工建筑物抗震设计规范》相关规定计算。沿建筑物高度作用于质点i的水平向地震惯性力代表值按下式计算：

式中：Fi为作用于质点i的水平向地震惯性力代表值；αh为与地震设计烈度相对应的地震加速度峰值，该工程αh=7度为0.1g，8度为0.2g；ε为地震作用的效应折减系数，拱坝一般取0.1；GEi为集中于质点i的重力作用标准值；αi为质点i的动态分布系数坝顶取3，坝基取1，按坝高线性递减。地震荷载计算表见表6。

表6 地震荷载计算表

3.2.3 其他荷载计算

当拱坝设有重力墩时，重力墩的应力和稳定分析，应符合SL 319-2005《混凝土重力坝设计规范》和SL 25-2006的规定，自重、扬压力等荷载计算参照规范中相关公式进行计算，计算结果见表7。

表7 其他荷载计算成果表

经计算，上下游边缘计算结果见表8。

表8 上下游边缘应力σz计算成果表Pa

从表8可以看出，在各种荷载组合作用下，均未产生拉应力，故上下游边缘应力满足要求。

3.2.4 抗滑稳定验算

重力墩抗滑稳定验算采用刚体极限平衡法进行计算，计算公式：

式中变量符号的含义参见林明财在硕士论文《拱坝重力墩抗滑稳定可靠度分析》一文[5]，不再详述。

经计算可知，在基本荷载组合一作用下抗滑稳定安全系数为2.4，基本荷载组合二作用下抗滑稳定安全系数为1.65，均大于1.05；在特殊荷载组合一作用下抗滑稳定安全系数为1.58，特殊荷载组合二作用下抗滑稳定安全系数为2.3，均大于1.00。故该重力墩的抗滑稳定满足要求。

4 结语

大坝应力及稳定计算直接关系到大坝安全性，是大坝设计的重要内容之一，尤其是在拱坝除险加固中判断大坝结构是否安全及是否能够继续使用起着决定性作用，若因计算不准确而判断大坝已不能正常使用，需要重新建设或者采取大量措施而造成投资增加；若因计算不准确而判断大坝结构安全，后期使用中可能会发生严重安全事故造成人身伤害及经济损失。实例研究表明，采用GADAP25拱坝软件开展大坝应力计算，能够准确预估不同工况下大坝应力分布情况，采用水利水电工程PC1500程序集进行稳定计算，结果符合《混凝土拱坝设计规范》相关要求，计算结果较为保守，满足设计需要。