水电站框格重力坝应力及稳定分析
重力坝的稳定分析
2.3
提高坝体抗滑稳定的工程措施
利用水重 将坝基开挖成倾向上游的斜面 在坝踵下设齿墙 抽水措施 加固地基 横缝灌浆 预加应力措施
重力坝的应力分析
【目的】 检验大坝在施工期和运用期是否满
足强度要求； 根据应力分布情况进行坝体混凝土
标号分区； 为研究坝体某些部位的应力集中和
配筋等提供依据。
E1
E1>E2 A处易出现拉应力
A E2
u y
0
影响坝体应力的各种因素
4、纵缝
n=0 时， 无影响
n>0 时， 坝踵合成应力
A
A处易出现拉应力
u y
0
n<0 时， 坝踵合成应力 对强度有利
影响坝体应力的各种因素
5、考虑分期施工
A
A
σ1
σ3
σ1
σ3
σ2
ቤተ መጻሕፍቲ ባይዱσy呈由原来直线分布变为折线分布，在坝
施工期
d 1
2kg / cm2 (拉)
影响坝体应力的各种因素
1、地基变形模量（E）
E 两端
AE E
沿坝底压应力σy，中间
坝踵A处易出现拉应力
u y
0
影响坝体应力的各种因素
2、地基变弹性模量
A E1 E3 E1
E3>E1 时，
A处易出现拉应力
u y
0
影响坝体应力的各种因素
3、坝体异弹性模量
扬 压 力
d
(
d y
p
pud )m2
u x
(p
puu ) u m1
重力坝的稳定性分析
缺点：丌能考虑岩体受力后所产生变形的影响， 极限状态不允许工作状态也有较大出入。
7.2 重力坝的稳定性分析
7.2.2 重力坝的稳定分析
有限单元法
可计算地基受力后的应力场 和位移场，并可模拟地基中软弱 结构面局部化效应及多场耦合作 用效应等，研究地基破坏发展全 过程。
根据重力坝设计规范的相关规定，承载能力极限状态： 坝体断面、结构及坝基岩体进行强度（应力）和抗滑稳定 性计算，必要时进行抗浮、抗倾验算；对需抗震设防的坝 及结构，尚需按抗震设计规范进行抗震验算。
稳定问题的类型：
抗滑稳定：坝体沿抗剪能力丌足的薄弱面产生滑动； 抗浮稳定：坝体在上、下游水荷载作用下，产生向上 浮起破坏形式； 抗倾稳定：上游坝踵以下岩体受拉产生斜裂缝及坝趾 以下岩体受压发生压碎而产生倾倒滑移破坏形式。
优点：可以考虑复杂地基的 局部化效应及材料的非线性本构 关系，模拟地基及坝体变形不破 坏的全过程等；
缺点：对有限元计算结果应 用及稳定判据的应用上尚需进一
7.2 重力坝的稳定性分析
7.2.2 重力坝的稳定分析
地质力等荷载，能形象地显示滑移破 坏的过程。
7.2 重力坝的稳定性分析
7.2.1 重力坝的工作原理
岩基上重力坝的基本剖面呈三角形，上游面通常是垂 直的或稍倾向上游的三角形断面。
工作原理：主要依靠坝体重量，在坝体和地基的接触 面上产生抗滑力来抵抗库水的推力，以达到稳定的要求。
Q
P
H
G
U
7.2 重力坝的稳定性分析
7.2.2 重力坝的稳定分析
稳定问题的提出：
优点：直观模拟坝体不地基 体系中的主要影响因素及变形不 破坏全过程；
《重力坝稳定分析》课件
重力坝失稳的后果
洪水泛滥：可能导致下游地区洪 水泛滥，造成人员伤亡和财产损 失
经济影响：可能导致下游地区经 济活动受到严重影响，影响当添加标题
生态破坏：可能导致下游生态环 境破坏，影响生物多样性和生态 系统平衡
社会影响：可能导致下游地区社 会稳定受到威胁，影响社会和谐 与稳定
01
提高重力坝稳定性的措施
优化坝体结构设计
优化坝体形状：采用流线型 设计，降低水流阻力
采用高强度混凝土：提高坝 体强度，降低变形
增加坝体厚度：提高坝体稳 定性，降低变形风险
采用预应力技术：提高坝体 抗拉强度，降低变形风险
加强材料强度和稳定性
选用高强度混凝土：提高混凝土 的抗压强度和抗拉强度
采用预应力混凝土：通过预应力 钢筋提高混凝土的抗拉强度和抗 剪强度
添加标题
添加标题
添加标题
添加标题
扬压力：水流通过重力坝时产生 的压力，影响坝体的稳定性
库水压力和扬压力的变化：库水 压力和扬压力的变化会影响坝体 的稳定性，需要及时监测和调整
温度变化和冰压力
温度变化：温度变化会 导致坝体材料性能变化， 影响坝体稳定性
冰压力：冰压力会导致 坝体受力不均，影响坝 体稳定性
公司
重力坝稳定分析
单击此处添加副标题
汇报人：
目录
单击添加目录项标题
01
重力坝的稳定性概述
02
重力坝稳定分析的方法
03
影响重力坝稳定的因素
04
重力坝稳定分析的实践应用
05
提高重力坝稳定性的措施
06
01
添加章节标题
01
重力坝的稳定性概述
重力坝的工作原理
重力坝的稳定性主要取决于 坝体的自重和地基的承载力
重力坝的稳定及应力分析
2. 公式：
K'
f ' ( W U ) c ' A
P
3.抗剪断参数的选定
对于大型工程，在设计阶段， f ′,c′应由野外及室内试验 成果决定。在规划阶段，可以参考规范给定的数值选用：
4.安全系数[K′] 设计规范规定： 不分等级,基本荷载组合：采用3.0; 特殊荷载组合：（1）采用2.5；（2）采 用不小于2.3。
地基的接触面、坝体折坡处或坝体断面
削弱的部位（如廊道、泄水管道等部 位）。
1） 基本假定
i.
坝体混凝土为均质、连续、各向同性 的弹性材料； 不考虑两侧坝体的影响，各坝段独立 工作； 假定坝体水平截面上的正应力σy按直 线分布，不考虑廊道等对坝体应力的 影响。
ii.
iii.
2） 边缘应力的计算
一般情况下，坝体的最大应力和 最小应力都出现在坝面，所以应该 首先校核坝体边缘应力是否满足强
坝 踵 坝 踵 坝 趾Fra bibliotek硬 库 满
软
Ec—— Er——
基坝 岩体
2、地基变形弹模对坝体 应力的影响 3、坝体异弹模对坝体应 力的影响 4、纵缝对坝体应力的影 响 5、分期施工对坝体应力 的影响(见下图) 6、坝踵断裂对坝体应力 的影响
坝体主应力分布示意图
影响坝体应力的主要因素有：
1)
地基变形对坝体应力的影响；
2 2
2u Pu
2 d Pd
3）内部应力的计算
1 、坝内水平截面上的正应力 σy 假 定和σy在水平截面上直线分布。 2、坝体内剪应力τ。 3、坝内水平正应力σx。 4、坝内主应力σ1和σ2。 5、考虑扬压力时的计算方法。
考虑扬压力作用时的应力计算
汉中市洋县酉水电站重力坝的强度和稳定性分析
的摩擦系数 ； ｔ ｇ 表示滑动表面存在 的真正的摩擦系数 。 为
保 证 安 全 ，有 要 求规 定 ，ｔ ｇ ０ ≤ｔ ｇ 。 但 上式 考 虑 软 基 在 地 基 上 滑 动 的情 况 ，而 忽 略 了对 安全 系数 的考 虑 ，对 此 ，需 引进 一 个 安 全 系数 Ｋ：
０期 第 １ ３卷 第 １
２ ０１ ３住
中 国
水
运
Ｖ ｏｌ １ ３ ０ｃ ｔｏ ｂｅｒ
Ｎ ｏ．１ Ｏ ２ ０１ Fra bibliotek １ ０月
Ｏ ｈ ｉ １ ３ ａ Ｗａｔ ｅｒ Ｔ ｒ ａｎｓ ｐ ｏｒ ｔ
汉 中市洋县酉水 电站重力坝的强度和稳定性分析
苏炳煌
（ 安 溪 县水 利 局 ，福 建 安 溪 ３ ６ ２ ４ ０ ０）
３ ． 沿 重 力 坝 基 面 的 稳 定性
Ｋ ：＿ 厂 ’ ∑ ＋ ｃ ’ Ａ ∑ｐ
离散性 ， 数值不 固定 ， 其稳定程度远低于 ｔ ｇ
和ｔ ｇ 的安 全 系 数 不 同 ，则 抗 剪 断 公 式 为 ：
随着重力坝 的高度增加，∑ Ｐ和￡ ｗ 和坝高平 方之 间呈正
比 相 关 ，部 分凝 聚 力 和 坝 高是 呈 正 相 关 的 。如 此 一 来 ，大 坝 越 高 ，相 应 的 安 全储 备 就越 低 。为 此 ，需 对 抗 剪 强 度 重 新 考 虑 ，得 出剪 摩 公 式 如 下 ：
ｆ 和 ｃ 代 表 混 凝 土 和 岩 基 间 的抗 剪 断 系 数 ，经 多 次 试 验 计 算 ，大 多数 情 况 下 ，ｆ＝ １ ． ０ — １ ． ３，Ｃ ＝０ ． ３ — １ ． ５ ＭＰ ａ 。
适 用恶 劣 条 件 和 不 同形 状 的 河谷 中 ；
第3节重力坝的稳定分析
第3节重力坝的稳定分析重力坝是一种常见的水利工程结构，广泛应用于水电站、灌溉渠道和排水系统等领域。
稳定性是设计和构建重力坝时必须考虑的重要因素之一。
本文将介绍重力坝的稳定性分析方法，以帮助读者更好地理解和应用在实际工程中。
一、稳定性分析的基本原理重力坝的稳定性分析是指通过力学的方法来评估坝体在受到水流、地震和土体压力等外力作用下的稳定性。
其基本原理是根据力的平衡和破坏准则对坝体进行分析。
稳定性分析的结果直接关系到坝体是否能够保持安全稳定，因此是设计中至关重要的环节。
稳定性分析通常包括静力分析和动力分析两个方面。
静力分析主要考虑坝体受到静水压力的作用，以及坝体自重和地震力等因素。
动力分析则关注坝体在地震和水流等动力荷载作用下的响应和变形。
二、静力分析方法1. 基本假设静力分析方法的基本假设是坝体在静水压力下处于静力平衡状态。
在分析中，可以假设坝体为刚体，计算坝顶的受力和坝底的抗力，以确定坝体的稳定性。
2. 受力计算在静力分析中，需要计算坝体所受的静水压力。
静水压力由上游水体的水位、坝体几何形状和水的密度决定。
通过计算坝顶的受力和坝底的抗力，可以确定坝体的受力情况。
3. 稳定性评估稳定性评估主要考虑坝体所受力矩和抗力矩之间的平衡关系。
如果受力矩大于抗力矩，坝体就会发生倾覆或滑移等破坏形式。
因此，需要通过计算力矩的大小来评估坝体的稳定性。
三、动力分析方法1. 地震力分析地震是重力坝的重要设计参数之一，也是动力分析的关键内容之一。
地震力分析需要考虑地震的频率和幅值，以及坝体的响应特性。
通过建立地震波模型和采用动力计算方法，可以得到坝体在地震作用下的响应和变形情况。
2. 水动力分析水动力分析主要考虑坝体在水流作用下的稳定性。
水流作用会对坝体施加水平力和垂直力，产生坝体的变形和振动。
通过建立水动力模型和采用数值计算方法，可以预测坝体的响应和变形情况，进而评估坝体的稳定性。
四、稳定性分析的实际应用稳定性分析方法在实际工程中有广泛应用。
重力坝应力稳定分析
程序中各符号的意义：YT——坝顶高程；YF——坝底高程；TD——坝顶宽度；YU——上游起坡高程；NN——上游坡比；YD——下游起坡高程；ND——下游坡比；Y8——正常工况时的上游蓄水位；Y9——正常工况时的下游水位；VZ——正常工况时的风速；FD——水库吹程；Y3——校核工况时的上游水位；Y4——校核工况时的下游水位；VF——校核工况时的风速；RD——坝体容重；YN——淤沙高程；RN——淤沙的浮容重；UT——扬压力折减点至坝踵距离；UH——折减系数；F——抗剪摩擦系数；FF——抗剪断摩擦系数；FC——抗剪断凝聚力；RC——砼强度标准值；YC——重要性系数；YG——廊道形心坐标（高程）；XG——廊道形心坐标（距上游距离）。
YU——上游起坡高程；NU——上游坡比；ND——下游坡比；RB——反弧半径；BY——鼻坎高程；BA——挑角；BT——鼻坎宽；DT——胸墙宽；DZ——胸墙高；DO——堰顶直线段长；LB——坝段宽；YT——坝顶高程；YF——坝基高程；YO——堰顶高程；RD——坝体容重；HS——定型设计水头；UT——扬压力折减点至坝踵距离；UH——折减系数；F——抗剪摩擦系数；FF——抗剪断摩擦系数；FC——抗剪断凝聚力；RC——砼强度标准值；YC——重要性系数；YN——淤沙高程；RN——淤沙的浮容重；NF——内摩擦角；VM——平均风速；FD——水库吹程；Y8——正常工况时的上游水位；Y9——正常工况时的下游水位；Y3——校核工况时的上游水位；Y4——校核工况时的下游水位；Q2——通过鼻坎的单宽流量；VQ——通过鼻坎的平均流速；DG——整个闸墩的重量；DX——闸墩重心与坝轴线距离；DL——闸墩厚度；GG——闸门重量；GX——闸门重心至坝轴线的距离；YG——廊道形心坐标（高程）；XG——廊道形心坐标（距上游距离）；ZU——上游第二起坡高程；N1——上游第二起坡的坡比；AH——上游堰面曲线宽度；BH——上游堰面曲线高度；KW——下游堰面曲线方程系数；NW——下游堰面曲线方程指数；。
重力坝稳定和应力计算
坝体强度承载能力极限状态计算及坝体稳定承载能力极限状态计算（一）、基本资料坝顶高程：1107.0 m校核洪水位（P = 0.5 %）上游：1105.67 m下游：1095.18 m 正常蓄水位上游：1105.5 m下游：1094.89 m死水位：1100.0 m混凝土容重：24 KN/m3坝前淤沙高程：1098.3 m泥沙浮容重：5 KN/m3混凝土与基岩间抗剪断参数值：f `= 0.5c `= 0.2 Mpa坝基基岩承载力：［f］= 400 Kpa坝基垫层混凝土：C15坝体混凝土：C1050年一遇最大风速：v 0 = 19.44 m/s多年平均最大风速为：v 0 `= 12.9 m/s吹程D = 1000 m（二）、坝体断面1、非溢流坝段标准剖面(1)荷载作用的标准值计算（以单宽计算）A 、正常蓄水位情况（上游水位1105.5m ，下游水位1094.89m ） ① 竖向力（自重）W 1 = 24×5×17 = 2040 KN W 2 = 24×10.75×8.6 /2 = 1109.4 KNW 3 = 9.81×（1094.5-1090）2×0.8 /2 = 79.46 KN ∑W = 3228.86 KNW 1作用点至O 点的力臂为： (13.6-5) /2 = 4.3 m W 2作用点至O 点的力臂为：m 067.16.83226.13=⨯- W 3作用点至O 点的力臂为：m 6.58.0)10905.1094(3126.13=⨯-⨯-竖向力对O点的弯矩（顺时针为“－”，逆时针为“+”）：M OW1 = 2040×4.3 = 8772 KN·mM OW2 = －1109.4×1.067 = －1183.7 KN·mM OW3 = －79.46×5.6 = －445 KN·m∑M OW = 7143.3 KN·m②静水压力（水平力）P1 = γH12 /2 = 9.81×(1105.5－1090)2 /2= －1178.4 KNP2 =γH22 /2 =9.81×(1094.89-1090)2 /2 = 117.3KN∑P = －1061.1 KNP1作用点至O点的力臂为：(1105.5－1090)/3 = 5.167mP2作用点至O点的力臂为：(1094.89－1090)/3 = 1.63m静水压力对O点的弯矩（顺时针为“－”，逆时针为“+”）：M OP1 = 1178.4×5.167 = －6089 KN·mM OP2 = 117.3×1.63 = 191.2 KN·m∑M OP = －5897.8 KN·m③扬压力扬压力示意图请见下页附图：H1 = 1105.5－1090 = 15.5 mH2 = 1094.89－1090 = 4.89 m(H1 －H1) = 15.5－4.89 = 10.61 m计算扬压力如下：U1 = 9.81×13.6×4.89 = 652.4 KNU2 = 9.81 ×13.6×10.61 /2 = 707.8 KN∑U = 1360.2 KNU1作用点至O点的力臂为：0 mU2作用点至O点的力臂为：13.6 / 2－13.6 / 3 = 2.267m 竖向力对O点的弯矩（顺时针为“－”，逆时针为“+”）：M OU1 = 0 KN·mM OU2 = －707.8×2.267 = －1604.6 KN·m∑M OU = －1604.6 KN·m④浪压力（直墙式）浪压力计算简图如下：由确定坝顶超高计算时已知如下数据：单位：m使波浪破碎的临界水深计算如下：%1%122ln 4h L h L L H m m m cr πππ-+=将数据代入上式中得到： 013.183.02644.783.02644.7ln 4644.7=-+=πππcr H 由判定条件可知，本计算符合⑴H ≥H cr 和H ≥L m /2，单位长度上的浪压力标准值按下式计算：)(41%1Z m W Wkh h L P +=γ 式中：γw ──水的重度 = 9.81 KN/m 3其余计算参数已有计算结果。
第四节重力坝的应力分析
W 6 M
a
B
B2
b
12 M B2
坝体内部应力计算图
（三）内部应力的计算
2、坝体内剪应力τ。呈抛物线分布
a1 b1x c1x2
3、坝内水平正应力σx。呈三次分布
x a2 b2 x c2 x2 d2 x3 接 近 直
线，对中小型工程可近似假定σx呈直线分布
运用期坝踵垂直应力不应出现拉应力（计扬压力）， 可按下式计算：
WR M RTR 0
AR
JR
WR B
6M R B2
0
规范要求，坝踵和坝体上游面的垂直应力的核算应按作用的标准值分别 计算作用的短期组合和长期组合。
2)短期组合下游坝面的垂直拉应力核算
•施工期属短暂状况，坝体下游面的垂直拉应力应不大于
x a3 b3x
4、坝内主应力σ1和σ2。(见课本P56式(3-48))
(四)考虑扬压力时的计算方法。
1、边缘应力 1)计算W和M计入扬压力 2)计算u 、xu、d、xd时
pu用 pu-puu代入 pd 用pd-pdd代入 2、求解坝内应力 可先不计扬压力计算、x、y，然后再 叠加由扬压力引起的应力
1、水平截面上的正应力σyu、σyd。 2、剪应力τu和τd。 3、水平正应力σxu和σxd。 4、主应力σ1u，σ2u和σ1d，σ2d。
1、水平截面上的正应力σyu、σyd。
根据偏心受压公式，坝
体上下游边缘垂直正应
力：
yu
W B
6M B2
yd
W B
6M B2
2、剪应力
4、边缘主应力
由于两个主应力面互相正交，由微分体
第一章 重力坝(4 应力分析)
边缘应力计算
？计算截面？荷载与应力的正方向规定？ 计算截面？荷载与应力的正方向规定？
1）水平截面上的正应力σyu、σyd。 水平截面上的正应力σ
按偏心受压公式计算
σ yu
σ yd
∑ W + 6∑ M = 2
B B
∑ W − 6∑ M = 2
B B
边缘应力计算 2）剪应力τu和τd。 ）剪应力
取上游坝面的微分体
由上游坝面微分体，根据平衡条件Σ 由上游坝面微分体，根据平衡条件ΣＦy = 0
σ1ucos2φudx=σyud x -Pusin2φudx
σ1u = σ yu cos 2 φu − Pu tgφu 2
σ 1u = ( 1 + n 2 )σ yu − pu n 2
σ 1d = ( 1 + m )σ yd − pd m
6.0 69.3 66.0 46.2 60.0
4.0 0.0 8.5
σxu
φu dy
根据平衡条件Σ 根据平衡条件ΣＦy=0
τ u dy + σ yu dx = p uφuds
dx τ u = (p u − σ yu ) dy
τu σyu
dx
τ u = (p u − σ yu )n
τ d = (σ yd − p d )m
边缘应力计算 水平正应力σ 3）水平正应力σxu和σxd 取上游坝面的微分体 σxu
重力坝稳定计算
某重力坝为三级建筑物，建在山区峡谷地区， 某重力坝为三级建筑物，建在山区峡谷地区，坝 顶五交通要求，上游设计洪水位为66.0m 66.0m， 顶五交通要求，上游设计洪水位为66.0m，相应是的下 游水位为4m 坝址处基岩面高程0.0m 坝基为凝灰岩， 4m， 0.0m， 游水位为4m，坝址处基岩面高程0.0m，坝基为凝灰岩， =0.55， 坝底与基岩之间的摩擦系数 f =0.55，扬压力折减系数 0.34。坝体断面见图，多年平均最大风速为16m/s 16m/s， 为0.34。坝体断面见图，多年平均最大风速为16m/s， 吹程D=3km 计算中暂不计泥沙压力， D=3km， 吹程D=3km，计算中暂不计泥沙压力，坝体材料为细骨 料混凝土，容重为24KN/m 根据上述资料要求， 料混凝土，容重为24KN/m3，根据上述资料要求，进行 坝体抗滑稳定计算，验算其安全系数能否满足要求， 坝体抗滑稳定计算，验算其安全系数能否满足要求， 若不满足，指出改进措施。 若不满足，指出改进措施。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第35卷第1期2010年2月广西大学学报：自然科学版Journal of Guangxi University ：Nat Sci Ed Vol.35No.1Feb.2010收稿日期：2009-09-25；修订日期：2009-10-31基金项目：广西自然科学基金资助项目（桂科自0991063）通讯联系人：张仲卿（1940-），男，四川内江人，广西大学教授，博士研究生导师；E-mail ：zhangzq@ 。
文章编号：1001-7445（2010）01-0157-05水电站框格重力坝应力及稳定分析张益智，李金友，唐若晗，张仲卿，赵小莲（广西大学土木建筑工程学院广西南宁530004）摘要：为研究广西拉浪水电站扩建工程破坝开孔对坝体稳定性的影响，运用了规范规定的材料力学法并结合三维有限元软件，对4个施工步骤共7种工况下进水口及整个坝体应力及稳定性进行分析，结果表明：材料力学法计算的孔口底部截面最大压应力为0.67MPa ，最小压应力为0.34MPa ，坝基截面最大压应力为0.63MPa ，最小压应力为0.06MPa ，均满足《混凝土重力坝设计规范SL319-2005》中坝体上游面的垂直应力不出现拉应力且不大于混凝土允许压应力的要求。
进水口以上坝体抗滑稳定安全系数K 为1.24，Kᶄ为5.28，全坝体抗滑稳定安全系数K 为1.10，Kᶄ为3.10，均大于《混凝土重力坝设计规范SL319-2005》规定的K 为1.0与Kᶄ为2.5，满足抗滑稳定要求。
在孔口周围、压力管上弯段及下弯段存在局部拉应力区，建议在局部拉应力区合理配筋。
关键词：水电站扩建工程；重力坝；进水口；应力状况；稳定分析中图分类号：TV311文献标识码：AStress and stability analysis of sash gravitydam for a hydropower stationZHANG Yi-zhi ，LI Jin-you ，TANG Ruo-han ，ZHANG Zhong-qing ，ZHAO Xiao-lian（College of Civil Engineering and Architecture ，Guangxi University ，Nanning 530004，China ）Abstract ：In order to study the stability of the gravity dam for Lalang hydropower station in Guangxi after an intake hole digged in the extension project ，material mechanics method and 3D finite ele-ment method are adopted to analyze stress and stability under four construction steps in seven condi-tions.The result indicates that maximum compressive stress of bottom section for intake hole is 0.67MPa and the minimum one is 0.34MPa through material mechanics method.And these two stresses of the whole dam are 0.63MPa and 0.06MPa.All of these stresses are less than the ones allowed in criterion （SL319-2005）.The stability safety coefficients K and K ᶄof the dam body above the in-take hole are 1.24and 5.28.The ones of the whole dam are 1.10and 3.10.All of these coeffi-cients are larger than the ones allowed in criterion （SL319-2005）.There are some tensile stressfields around the intake hole ；upper curved section and lower curved section of the pressure tube.It is suggested that these areas should be reinforced with steel bars.Key words ：extension project for hydropower station ；gravity dam ；intake ；stress condition ；stability analysis1工程概况为充分利用水利资源发电，以缓解当地的电力供需矛盾，改善当地的投资环境，促进地方经济发展，851广西大学学报：自然科学版第35卷提高人民的物质和文化生活水平，经过有关部门批复对广西拉浪水电站进行扩建。
扩建电站为坝后式电站，装机容量为1.9ˑ104kW，多年平均年发电量3.8ˑ107kW·h，采用原正常蓄水位177.00m，死水位174.00m，设计洪水位177.15m，校核洪水位179.58m，坝高32.00m，坝顶高程181.40m。
扩建工程主要建筑物由引水建筑物、门机排架、发电厂房、尾水渠和主变压器场等组成。
2扩建工程施工方案及特点为了满足下游工厂及居民的用水需求，扩建工程采用不放库水的施工方案。
发电引水建筑物布置于右岸重力坝段中部，由坝式进水口和坝内式钢筋混凝土矩形压力管组成，中心线垂直于坝轴线。
其中坝式进水口结构施工采用15m（长）ˑ8m（宽）的临时闸门挡住坝前库水，拆除原坝段（沿坝轴线18m宽）框格填碴重力坝隔墩，新建两孔宽高4.50m，宽4.75m的进水口，中间用1m宽的中墩隔开。
这种采用临时闸门挡水的破坝施工方案在当地实属少见，而且施工时在18m宽坝段上开12m宽的进水口，孔口左侧坝体宽度仅为2.6m，右侧坝体宽度仅为3.4m，坝体自重明显减小，而坝前水位并没有降低，仍为正常水压，进水口两侧薄坝体要承受十多米高的水压力，所以施工时进水口处坝体应力状况会恶化，而且坝体稳定性会降低。
根据国家电力监管委员会大坝安全监察中心安全监察处2009年“关于补充拉浪水电站扩建工程安全评价所需资料的函”、“拉浪水电站扩建工程对大坝安全影响评价专家组会议纪要”，决定对破坝施工时的进水口及坝体的应力、稳定性进行重点研究［1-4］。
针对进水口破坝施工的特殊性，进水口施工拟采用4个步骤进行，具体如下：①进行中墩开挖和浇筑。
中墩采用一次浇注完成。
这样形成施工工作面，有利于下一步的工作面展开。
中墩开挖在原重力坝中部空腔体内进行，左侧开挖面与空腔左侧内侧面相平，右侧预留50cm的距离做为施工工作空间。
②进行左侧结构开挖和浇筑。
左侧结构的开挖要求在中墩钢筋砼结构强度达到70%设计强度后，方可进行。
根据左侧结构尺寸要求先开挖左下游半截结构，开挖深度3.3m，待该半截钢筋砼结构达70%设计强度并对顶板上部进行回填灌浆合格后，再进行挡水面板开挖和砼结构浇筑及回填灌浆。
③进行右侧结构开挖和浇筑，右侧结构开挖和浇筑的要求同左侧。
④进水口结构顶部灌浆及质量检测，进水口框型结构钢筋砼达70%设计强度后，方可进行顶部回填灌浆。
全部回填灌浆检查合格后，方可进行临时挡水闸门拆除。
3坝体及进水口应力稳定分析3.1材料力学法及计算结果分析重力坝坝基截面及进水口底部截面的垂直应力计算公式采用《混凝土重力坝设计规范SL319-2005》中规定的垂直应力计算公式，坝体抗滑稳定计算采《混凝土重力坝设计规范SL319-2005》中规定的抗剪断强度公式或抗剪强度公式。
计算工况为最危险的施工工况，即采用临时闸门挡住库水，破坝后中墩未浇筑，孔口开挖尺寸沿坝轴线为12m宽，高度为6.5m，上游水位为校核洪水位179.58m，下游无水工况，并考虑扬压力作用，坝基及进水口截面垂直正应力与安全系数计算结果见表1及表2。
表1截面垂直应力计算结果Tab.1Calculation result of section perpendicularity stress截面位置∑W/kN A/m2M/（kN·m）J/m4x/mσy min/MPaσy max/MPa 坝基截面 1.30ˑ105379.00 3.79ˑ10514011.010.530.060.63进水口底部截面 3.55ˑ10470.56 2.32ˑ104813.2 5.880.340.67第1期张益智等：水电站框格重力坝应力及稳定分析表2抗滑稳定安全系数K 及Kᶄ计算结果Tab.2Calculation result of slide resistance coefficient K and Kᶄ部位∑W /kN A /m 2∑P /kPa f f ᶄcᶄ/kPa K Kᶄ全坝体 1.30ˑ105379.008.25ˑ1040.70.8400 1.10 3.10进水口以上坝体3.55ˑ10470.562.01ˑ1040.71.010001.245.283.2有限元法由于进水口施工工况以及坝体与压力管作用力的复杂，用规范规定的材料力学法计算只能进行坝基截面的简化应力计算，难以反映坝体各部分真实的应力情况，故采用大型通用有限元软件ANSYS 对坝体及进水口结构应力进行分析［1-10］。
图1有限元计算模型Fig.1Finite element calculation model 3.2.1计算模型及材料参数有限元计算模型取进水口施工坝段18m 进行模拟，坝段坐标原点取在进水口上弯段处，沿坝长指向左岸为X 正方向，竖直向上为Y 正方向，顺河向为Z 正方向，坝基上游取1.5倍坝高，下游取2倍坝高，压力管末端距水轮机轴线6m 。
有限元分析采用3－D 结构实体单元SOLID45建立三维实体结构模型，并采用接触单元TARGE170及CONTA174模拟3条伸缩缝。
整个模型共划分68122个单元，76645个节点，其中坝体单元42602个，节点49188个，模型见图1。
计算中取旧坝体C10混凝土弹性模量E =17.5GPa ，泊松比ν=0.167，密度为2400kg /m 3，新浇C20混凝土弹性模量E =25.5GPa ，泊松比ν=0.167，密度为2400kg /m 3，基岩弹性模量E =0.75GPa ，泊松比ν=0.250，密度为2650kg /m 3。