进水口整体稳定及基底应力计算
进水口整体稳定及基底应力计算
1. 计算总说明 1.1 计算目的及要求
某水电站进水口根据电站枢纽布置、地形、地质条件设为岸塔式进水口,“镶嵌”在L 型基础中,塔背有基岩对其起支撑作用,靠自重和岸坡岩体支撑维持稳定,加之该进水口置于土质地基上,因建基面不允许出现拉应力,因此可不进行抗倾覆稳定计算。通过对进水口整体抗浮稳定与基底应力计算,以复核其是否满足规范要求。 1.2 基本资料
进水口纵横剖面结构尺寸见附图。
3
水容重: 10k N /m
钢筋混凝土容重: 25kN /m 3
基础与混凝土之间f 值为:f '=0.4,c '=0.08Mpa 。 地基承受能力: [σ]=0. 4 M 2P a 校核洪水位: 1886. 1m 0 9设计洪水位: 1884. 0m 6 9正常蓄水位: 1900. 0m 0 0死水位: 1896. 0m 0 0多年平均最大风速: 11. 0m 0s /
k m 吹程: 0. 7
拦污栅及对应启闭机重: 52t 、8t 事故检修门及对应启闭机重: 38t 、15t 计算采用规范及参考书:
a. 《水电站进水口设计规范》DL/T5398-2007
b. 《水工建筑物荷载设计规范》DL5077-1997 c. 《水工建筑物荷载设计规范》 d. 《水电站厂房设计规范》SL266-2001 e .《公路桥涵设计通用规范》JTJ021-89
f. 《水闸设计规范》SL265-2001;
g. 《土力学》河海大学出版社出版, 卢廷浩主编; h. 某水电站首部枢纽平面布置图及主要建筑物剖面图; i. 《孔隙比与压力关系曲线表》; 1.3 计算原理及假定
岸塔式进水口根据自身的结构特点,只要基底应力在允许应力或允许抗力范围内,塔体就不会发生整体失稳,根据规范要求及进水口的布置情况,对塔体做以下假定:
1)将塔体视为刚体,在荷载作用下,岩体变形产生抗力; 2)抗力或反力强度值采用材料力学法和刚体极限平衡法求得。 3)计算地基沉降量时,不考虑周边结构对进水口塔体的影响。 此算稿按荷载组合{基本组合、特殊组合}计算。 工况1:正常运行期,正常蓄水位1900.000m ; 工况2:正常运行期,设计洪水位1884.069m ; 工况3:非常运行期,校核洪水位1886.109m ; 工况4:检修期,正常蓄水位1900.000m ; 工况5:施工完建未挡水期。 工况6:正常蓄水位+Ⅶ度地震
工况1、2为基本组合,工况3、4、5、6为特殊组合。 2. 计算过程 2.1 截面计算
取进水口底板进行计算,进水口底板顺水流方向为x 轴,垂直水流方向为y 轴,截面尺寸如图2.1-1所示。
y
x
图2.1-1 建基面示意图
基底面积:
I x =
112112⨯bh =
3
112112
⨯20⨯18.00=9720m
34
4
I y =
⨯b h =
3
⨯20⨯18.00=12000m
3
荷载作用正方向取为,垂直荷载作用:铅直向下为正;水平荷载作用:沿水流方向为正;弯矩作用:逆时针方向为正。 2.2 力矩计算
下面计算各截面形心对x 、y 轴的力矩,因为进水口左、右侧对称布置,所以只需计算进水口上下游基底应力,即可只算截面计算点对y 轴的力矩即可。
力矩计算简图见附图。 2.2.1 底板
重量为:
G 1=γQ γc V a =1.05⨯25⨯20.0⨯3.0⨯18.0=28350.0kN
l x 1=0 M
y 1
=G 1l x 1=0kN m
2.2.2 左、右边墙
G 2=γQ γc V =1.05⨯51.94⨯21.0⨯25⨯2=57263.85kN
l x 2=-1.19m
截面计算点至形心轴Y 的距离产生的弯矩
M
y
=G 2l x 2=57263.85⨯(-1.19) =-68143.98kN m
2.2.3 拦污栅中墩
G 3=γQ γc V =1.05⨯13.77⨯21.0⨯25=7590.71kN
l x 3=5.84m
截面计算点至形心轴Y 的距离产生的弯矩
M
y
=G 3l x 3=7590.71⨯5.84=44329.76kN m
2.2.4 胸墙
胸墙指的是进水口拦污栅之后的那部分实体混凝土,由于喇叭口上侧混凝土是一不规则的结构,在取其截面计算点时分两部分计算。
G 4.1=γQ γc V a =1.05⨯34.31⨯12.5⨯25=11257.97kN G 4.2=γQ γc V b =1.05⨯25⨯19.57⨯8.92=4582.32kN
l x 4.1=3.25m
、l x 4.2=-0.31m
截面计算点至形心轴Y 的距离产生的弯矩
M y M
y 4a
=G 4.1l x 4.1+G 4.2l x 4.2=11257.97⨯3.25+4582.32⨯(-0.31) =35167.88kN m
2.2.5 通气孔四周实体混凝土
由于通气孔四周混凝土在计算边墙时已经计算了一部分,因此在此只计算流道上方混凝土(包括扣除部分)。
G 5.1=γQ γc V =1.05⨯25⨯11.94⨯15.8=4952.12kN 扣除通气孔体积:
G 5.2=γQ γc V =1.05⨯25⨯2⨯15.8=829.5kN
l x 5.1=-8.28m m l x 5. 2=-8. 5
截面计算点至形心轴Y 的距离产生的弯矩
M
y 5
=G 5.1l x 5.1-G 5.2l x 5.2=4952.12⨯(-8.28) -829.5⨯(-8.5) =-33952.76kN m
2.2.6 进水口板梁
顶部板梁
G 6.1=γQ γc V =1.05⨯25⨯16.03⨯0.25⨯2=210.39kN G 6.2=γQ γc V =1.05⨯25⨯37.4⨯0.25=245.44kN G 6.3=γQ γc V =1.05⨯25⨯3.57⨯0.8=74.97kN G 6.4=γQ γc V =1.05⨯25⨯0.25⨯17.85=117.14kN G 6.5=γQ γc V =1.05⨯25⨯2.08⨯0.8=43.68kN
l x 6.1=8.73m
、l x 6.2=0.64m 、l x 6.3=-1.25m 、l x 6.4=-2.7m 、l x 6.5=-4.4m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 6
=
∑G
6i x 6i
l =210.39⨯8.73+245.44⨯0.64+74.97⨯(-1.25) +117.14⨯(-2.7) +43.68⨯(-4.4)
3 m =1391. k 6N
2.2.7 进水口启闭机排架
将启闭机排架分为两部分计算,即拦污栅部分和闸门部分。
G 7.1=γQ γc V =1.05⨯25⨯79.74=2093.18kN G 7.2=γQ γc V =1.05⨯25⨯55.66=1461.18kN
(拦污栅排架) (事故检修门排架)
l x 7.1=5.85m
、l x 7.2=-5.5m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 7
=
∑G l
7x 7
=2093.18⨯5.85+1461.18⨯(-5.5) =4209.16kN m
2.2.8 拦污栅及启闭机
G 8=1.05⨯(52+8) ⨯10=630.00kN l x 8=4.8m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 8
=G 8l x 8=630.00⨯4.8=3024.00kN m
2.2.9 事故检修门及启闭机
G 9=1.05⨯(38+15) ⨯10=556.50kN
l x 9=-5.8m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 9
=G 9l x 9=556.5⨯(-5.8) =-3227.70kN m
2.2.10 人群荷载
按《公路桥涵设计通用规范》,人群荷载一般取3kN /m 2。
G 10=1.05⨯3⨯18.00⨯20.00=1134.00kN l x 10=0
截面计算点至形心轴Y 的距离产生的弯矩
M
y 10
=G 10l x 10=0kN m
2.2.11 静水压力
根据岸塔式进水口的特点,各种工况时,其上下游水位相同,故上下游静水压力相互抵消,产生的弯矩和也为0,可不予计算。 2.2.12 扬压力
根据《水电站厂房设计规范》,扬压力分浮托力和渗透压力,进水口上下游水位相等,故渗透压力为0。 浮托力计算: 工况一:
U =1.0⨯10⨯(1900.00-1877.00) ⨯18.00⨯20.00=82800.00kN
(↑)
浮托力沿底板成对称分布,故产生的弯矩为0。 工况二:
U =1.0⨯10⨯(1884.069-1877.00) ⨯18.00⨯20.00=25448.40kN
(↑)
浮托力沿底板成对称分布,故产生的弯矩为0。 工况三:
U =1.0⨯10⨯(1886.109-1877.00) ⨯18.00⨯20.00=32792.40kN
(↑)
浮托力沿底板成对称分布,故产生的弯矩为0。 工况四:
U =1.0⨯10⨯(1900.00-1877.00) ⨯18.00⨯20.00=82800.00kN
(↑)
浮托力沿底板成对称分布,故产生的弯矩为0。 工况五:施工完建未挡水期,此时浮托力为0。 2.2.13 风荷载
垂直于建筑物表面上的风荷载标准值按下式计算:
w k =βz μz μs w o
式中:
w k ——风荷载标准值,kN /m
2
βz ——z 高程的风振系数,取1.0
: μs ——风荷载体型系数,查《水工建筑物荷载设计规范》
μs =-(0.48+0.03H /B ) μz ——风压高度变化系数。
w o ——基本风压,查全国风压基本分布图可确定为0.3kN /m
2
风荷载的作用分项系数采用1.3。
w k =βz μz μs w o
工况一:
风压高度变化系数,按《水工建筑物荷载设计规范》DL5077—1997表12.1.6查得:μz =0.68。
μs =-(0.48+0.03H /B ) =-[0.48+0.03⨯(1902.00-1900.00) /18.00]=-0.483
2
=1. 0⨯0. 6⨯8-(0. 4⨯83) =0-. 3k N 0. 0m w k =βz μz μw 986s o
/
风压力:
w =1.3⨯(-0.0986) ⨯18.00⨯(1902.00-1900.00) =-4.61kN l x 13=
23
⨯(1902.00-1900.00) +(1900.00-1881.00) =20.33m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 13
=P 13l x 13=-4.61⨯20.33=-93.83kN m
工况二:
风压高度变化系数,按《水工建筑物荷载设计规范》DL5077—1997表12.1.6
查得:μz =1.204
μs =-(0.48+0.03H /B ) =-[0.48+0.03⨯(1902.00-1884.069) /18.00]=-0.5099
2
w k =βz μz μw 1842=1. 0⨯1. 20⨯4-(0. 50⨯99) =-0. 3k N 0. m s o
/
风压力:
w =1.3⨯(-0.1842) ⨯18.00⨯(1902.00-1884.069) =-77.275kN l x 13=
23
⨯(1902.00-1884.069) +(1884.069-1881.00) =15.02m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 14
=P 13l x 13=-77.275⨯15.02=-1160.91kN m
工况三:
风压高度变化系数,按《水工建筑物荷载设计规范》DL5077—1997表12.1.6 查得:μz =1.16
μs =-(0.48+0.03H /B ) =-[0.48+0.03⨯(1902.00-1886.069) /18.00]=-0.5065
2
=1. 0⨯1. 1⨯6-(0. 50⨯65) =0-. 3k N 0. m 1 w k =βz μz μw 763s o
/
风压力:
w =1.3⨯(-0.1763) ⨯18.00⨯(1902.00-1886.109) =-65.541kN l x 13=
23
⨯(1902.00-1886.109) +(1886.109-1881.00) =15.70m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 13
=P 13l x 13=-65.541⨯15.7=-1029.19kN m
工况四:
风压高度变化系数,按《水工建筑物荷载设计规范》DL5077—1997表12.1.6 查得:μz =0.68
μs =-(0.48+0.03H /B ) =-[0.48+0.03⨯(1902.00-1900.00) /18.00]=-0.483
2=1. 0⨯0. 6⨯8-(0. 4⨯83) =0-. 3k N 0. 0m w k =βz μz μw 986s o
/
风压力:
w =1.3⨯(-0.0986) ⨯18.00⨯(1902.00-1900.00) =-4.61kN l x 13=
23
⨯(1902.00-1900.00) +(1900.00-1881.00) =20.33m
截面计算点至形心轴Y 的距离产生的弯矩
M
=P 13l x 13=-4.61⨯20.33=-93.83kN m
y 13
工况五:
风压高度变化系数,按《水工建筑物荷载设计规范》DL5077—1997表12.1.6 查得:μz =1.267
μs =-(0.48+0.03H /B ) =-[0.48+0.03⨯(1902.00-1881.00) /18.00]=-0.515
2
w k =βz μz μw 958=1. 0⨯1. 26⨯7-(0. 5⨯15) =0-. 3k N 0. m 1s o
/
风压力:
w =1.3⨯(-0.1958) ⨯18.00⨯(1902.00-1881.00) =-96.192kN l x 13=
23
⨯(1902.00-1881.00) =14.00m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 13
=P 13l x 13=-96.192⨯14.00=-1346.69kN m
2.2.14 浪压力 官厅公式:
gh v
2
=0.0076v 0
-1/12
(
gD v (
20
)
1/3
gL m v
20
=0.331v 0
-1/2.15
gD v
20
)
1/3.75
式中: h ——当h ——当
L m
gD v 0
2
=20~250时,为累计频率5%的波高;
gD v 0
2
=250~1000
时,为累计频率10%的波高;
——平均波长,m;
D ——吹程,m ;
v 0
——多年平均最大风速,m /s ,按《水工建筑物荷载设计规范》DL5077
—1997附录G 表G1查得,正常蓄水位时取1.25v 0(风速系数由内插法求得),设计洪水位时取0.924v 0,校核洪水位时取0.949v 0。 临界水深的计算:
L m 4
2
H cr =
ln
L m +2πh 1%L m -2πh 1%2πH L m
h z =
πh 1%
L m
cot
式中:
H cr ——使波浪破碎的临界水深。 h 1%
——累积频率为1%时的波高,由《水工建筑物荷载设计规范》DL5077-1997
附录G 表G2用插入法可以求得。
h z
——波浪中心线至计算水位的高度,m 。
浪压力标准值计算: 当H ≥H cr ,H ≤
P w k =
12
L m 2
时,单位长度上的浪压力标准值计算公式为:
[(h 1%+h z )(γw H +P lf ) +H P lf ]
L m 2
当H ≥H cr ,H ≥
P w k =
14
时,单位长度上的浪压力标准值计算公式为:
γW L m (h 1%+h z )
式中:
P w k
——单位长度迎水面上的浪压力标准值
h z
——波浪中心线至计算水位的高度,m 。
γW ——水的重度, kN /m 3
P lf
——建筑物底面处的剩余浪压力强度值,kN /m 2,按下式计算:
2πH L m
P lf =γw h 1%sec h
浪压力的作用分项系数采用1.2。 工况一:
gh v
20
=0.0076v 0
-1/12
(
gD v
20
)
1/3
=0.0076⨯(1.25⨯11.00)
-1/12
⨯(
9.81⨯7001.25⨯11.00
2
2
)
1/3
h =0.0076⨯(1.25⨯11.00)
-1/12
⨯(
9.81⨯7001.25⨯11.00
2
2
)
1/3
⨯(1.25⨯11.00) ÷9.81=0.39m
2
gL m v
20
-1/2.15
=0.331v 0
(
gD v
20
)
1/3.75
=0.331⨯(1.25⨯11.00)
9.81⨯7001.25⨯11.00
2
-1/2.15
⨯(
9.81⨯7001.25⨯11.00
2
22
)
1/3.75
L m =0.331⨯(1.25⨯11.00)
临界水深H cr
-1/2.15
⨯(
) 2
1/3.75
⨯(1.25⨯11.00) ÷9.81=4.91m
的计算
2
由
gD v
2
=
9.81⨯700(1.25⨯11.00)
=36.32
,比值在
gD v
20
=20~250
之间,为累计频率5%
的波高; 令
h 5%h m =
=1.95,则h m =
0.20
1900.00-1877.00
=0.0087
h 5%1.95
=0.391.95
=0.20,
h m H
m
=0.0087
根据
h m H
m
, 查《水工建筑物荷载设计规范》DL5077-1997附录G 表
h 1%h m
=2.406
G2用插入法可以求得
h 1%=2.406h m =2.406⨯0.20=0.48m
4.91+2⨯3.14⨯0.484.91-2⨯3.14⨯0.48
=0.56m
H cr =
L m 4π
2
ln
L m +2πh 1%L m -2πh 1%2πH L m
=
=
4.914⨯3.14
2
ln
h z =
πh 1%
L m
cot
3.14⨯0.48
4.91
cot
2⨯3.14⨯(1900.00-1877.00)
4.91
=0.15m
H =1900.00-1877.00=23m >H cr =0.56m H =23m >
14
L m 2=4.912
=2.46m 14
,单位长度上的浪压力标准值计算公式为:
P w k =γW L m (h 1%+h z ) =
⨯10⨯4.91⨯(0.48+0.15) =7.73kN /m
浪压力:
P w k =1.2⨯7.73⨯18.00=166.9kN
(→)
计算点至形心的距离
l x 14=h z +1900.00-1877.00=0.15+1900.00-1877.00=23.15m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 14
=p w k l x 14=-166.9⨯23.15=-3863.41kN m
工况二:
gh v 0
2
=0.0076v 0
-1/12
(
gD v 0
2
)
1/3
=0.0076⨯(0.924⨯11.00)
-1/12
⨯(
9.81⨯7000.924⨯11.00
2
) 2
1/3
h =0.0076⨯(0.924⨯11.00)
-1/12
⨯(
9.81⨯7000.924⨯11.00
2
2
)
1/3
⨯(0.924⨯11.00) ÷9.81=0.267m
-
2
gL m v 0
2
=0.331v 0
-1/2. 15
(
gD v 0
2
)
1/3. 75
=0.331⨯(0.924⨯11.00)
9.81⨯7001.204⨯11.00
2
2
⨯(
1/2. 15
9.81⨯700
2
0.924⨯11.00
2
) 2
1/3. 75
L m =0.331⨯(0.924⨯11.00)
临界水深H cr
-1/2.15
⨯()
1/3.75
⨯(0.924⨯11.00) ÷9.81=3.63m
的计算
2
由
gD v
2
=
9.81⨯700(0.924⨯11.00)
=66.47,比值在
gD v
20
=20~250
之间,为累计频率5%
的波高; 令
h 5%h m =
=1.95,则h m =
0.137
1884.069-1877.00
=0.019,
h 5%1.95
=0.2671.95
=0.137
,
h m H
m
=0.019
根据
h m H
m
查《水工建筑物荷载设计规范》DL5077-1997附录G 表
3.63+2⨯3.14⨯0.333.63-2⨯3.14⨯0.33
=0.37m
G2用插入法可以求得
h 1%h m
=2.39
h 1%=2.39h m =2.39⨯0.137=0.33m
H cr =
L m 4π
2
ln
L m +2πh 1%L m -2πh 1%2πH L m
=
=
3.634⨯3.14
2
ln
h z =
πh 1%
L m
cot
3.14⨯0.33
3.63
cot
2⨯3.14⨯(1884.069-1877.00)
3.63
=0.09m
H =1884.069-1877.00=7.069m >H cr =0.37m H =7.069m >
14
L m 2=3.632
=1.83m 14
,单位长度上的浪压力标准值计算公式为:
P w k =γW L m (h 1%+h z ) =
⨯10⨯3.63⨯(0.33+0.09) =3.81kN /m
浪压力:
P w k =1.2⨯3.81⨯18.00=82.36kN
(→)
计算点至形心的距离
l x 14=h z +1884.069-1877.00=0.09+1884.069-1877.00=7.16m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 14
=p w k l x 14=-82.36⨯7.16=-589.86kN m
工况三:
gh v
20
=0.0076v 0
-1/12
(
gD v
20
)
1/3
=0.0076⨯(0.949⨯11.00)
9.81⨯7001.16⨯11.00
2
2
-1/12
⨯(
9.81⨯7000.949⨯11.00
2
2
) 2
1/3
h =0.0076⨯(0.949⨯11.00)
-1/12
⨯()
1/3
⨯(0.949⨯11.00) ÷9.81=0.276m
-
gL m v 0
2
=0.331v 0
-1/2. 15
(
gD v 0
2
)
1/3. 75
=0.331⨯(0.949⨯11.00)
9.81⨯7000.949⨯11.00
2
2
⨯(
1/2. 15
9.81⨯700
2
0.949⨯11.00
2
) 2
1/3. 75
L m =0.331⨯(0.949⨯11.00)
临界水深H cr
-1/2.15
⨯()
1/3.75
⨯(0.949⨯11.00) ÷9.81=3.73m
的计算
2
由
gD v
2
=
9.81⨯700(0.949⨯11.00)
=63.02,比值在
gD v
20
=20~250
之间,为累计频率5%
的波高; 令
h 5%h m =
=1.95,则h m =
0.14
1886.109-1877.00
=0.016,
h 5%1.95
=0.2761.95
=0.14,
h m H
m
=0.016
根据
h m H
m
查《水工建筑物荷载设计规范》DL5077-1997附录G 表
=2.395
G2用插入法可以求得
h 1%h m
h 1%=2.395h m =2.395⨯0.14=0.34m
3.73+2⨯3.14⨯0.343.73-2⨯3.14⨯0.34
=0.39m
H cr =
L m 4π
2
ln
L m +2πh 1%L m -2πh 1%2πH L m
=
=
3.734⨯3.14
2
ln
h z =
πh 1%
L m
cot
3.14⨯0.34
3.73
cot
2⨯3.14⨯(1886.109-1877.00)
3.73
=0.1m
H =1886.109-1877.00=9.109m >H cr =0.39m H =9.109m >
L m 2=3.732
=1.87m
,单位长度上的浪压力标准值计算公式为:
14
14
P w k =γW L m (h 1%+h z ) =
⨯10⨯3.73⨯(0.34+0.1) =4.07kN /m
浪压力:
P w k =1.2⨯4.07⨯18.00=87.87kN
(→)
计算点至形心的距离
l x 14=h z +1886.109-1877.00=0.1+1886.109-1877.00=9.21m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 14
=p w k l x 14=-87.87⨯9.21=-808.92kN m
工况四:
该工况时的浪压力等于工况一时的浪压力,截面计算点至形心的轴Y 的距离产生的弯矩也与工况一的相等。 浪压力:
P w k =166.9kN
(→)
截面计算点至形心轴Y 的距离产生的弯矩
M
y 14
=-3863.41kN m
工况五:
施工完建未挡水期,浪压力为0,截面计算点至形心轴Y 的距离产生的弯矩也为0。
2.2.15 胸墙及检修门槽中的内水压力 工况一:
G 15.1=1.05⨯10⨯44.66⨯(1900.00-1888.954) =5179.80kN G 15.2=1.05⨯10⨯18.58⨯(1900.00-1888.454) =2252.51kN G 15.3=1.05⨯10⨯7.92⨯(1900.00-1887.200) =1064.45kN
计算点至形心的距离
l x 15.1=0.3m 、l x 15.2=-3.36m
、l x 15.3=-5.7m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 15
=p w k l x 15=5179.80⨯0.3+2252.51⨯(-3.36) +1064.45⨯(-5.7) =-12081.84kN m
工况二:
设计洪水位时,胸墙和门槽内无水,故内水压力为0。 工况三:
校核洪水位时,胸墙和门槽内无水,故内水压力为0。 工况四:
该工况上游水位与工况一时相同,故内水压力和其对形心轴产生的弯矩均相同。
G 15.1=5179.80kN G 15.2=2252.51kN G 15.3=1064.45kN
计算点至形心的距离
l x 15.1=0.3m 、l x 15.2=-3.36m
、l x 15.3=-5.7m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 15
=p w k l x 15=-12081.84kN m
工况五:
施工完建未挡水期,此时进水口中的内水压力为0。 2.2.16 进水口中的内水压力 G 16=γQ γc (V 1+V 2+V 3+V 4) 式中:
γQ —可变作用的分项系数;
V 1—拦污栅槽(包括栏污栅槽前) 中的水体积; V 2—喇叭口中的水体积; V 3—检修闸门槽中的水体积; V 4—边墩两侧的水体积。
工况一:
G 16.1=γQ γc V 1=1.05⨯10⨯81.06⨯(1900.00-1881.00) =16171.47kN G 16.2=γQ γc V 2=1.05⨯10⨯95.21⨯5.75=5748.42kN G 16.3=γQ γc V 3=1.05⨯10⨯29.76⨯5.2=1624.90kN
G 16.4=γQ γc V 4=1.05⨯10⨯25.41⨯(1900.00-1881.00) ⨯2=10138.99kN
截面计算点至形心轴Y 的距离为
l x 16.1=6.84m
、l x 16.2=0.78m 、l x 16.3=-7.21m 、l x 16.4=-6.17m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 16
=γQ γc (V 1l x 16.1+V 2l x 16.2+V 3l x 16.3+V 4l x 16.4)
=1.05⨯10⨯[16171.47⨯6.84+5748.42⨯0.78+1624.90⨯(-7.21) +10138.99⨯(-6.17) ]
=40823.56kN m
工况二:
G 16.1=γQ γc V 1=1.05⨯10⨯81.06⨯(1884.069-1881.00) =2612.12kN G 16.2=γQ γc V 2=1.05⨯10⨯95.21⨯5.75=5748.42kN G 16.3=γQ γc V 3=1.05⨯10⨯29.76⨯5.2=1624.90kN
G 16.4=γQ γc V 4=1.05⨯10⨯25.41⨯(1884.069-1881.00) ⨯2=1637.65kN
截面计算点至形心轴Y 的距离为:
l x 16.1=6.84m
、l x 16.2=0.78m 、l x 16.3=-7.21m 、l x 16.4=-6.17m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 16
=γQ γc (V 1l x 16.1+V 2l x 16.2+V 3l x 16.3+V 4l x 16.4)
=1.05⨯10⨯[2612.12⨯6.84+5748.42⨯0.78+1624.90⨯(-7.21) +1637.65⨯(-6.17) ]
=510.74kN m
工况三:
G 16.1=γQ γc V 1=1.05⨯10⨯81.06⨯(1886.109-1881.00) =4348.42kN G 16.2=γQ γc V 2=1.05⨯10⨯95.21⨯5.75=5748.42kN G 16.3=γQ γc V 3=1.05⨯10⨯29.76⨯5.2=1624.90kN
G 16.4=γQ γc V 4=1.05⨯10⨯25.41⨯(1886.109-1881.00) ⨯2=2726.21kN
截面计算点至形心轴Y 的距离为:
l x 16.1=6.84m
、l x 16.2=0.78m 、l x 16.3=-7.21m 、l x 16.4=-6.17m
截面计算点至形心轴Y 的距离产生的弯矩
M
y 16
=γQ γc (V 1l x 16.1+V 2l x 16.2+V 3l x 16.3+V 4l x 16.4)
=1.05⨯10⨯[4348.42⨯6.84+5748.42⨯0.78+1624.90⨯(-7.21) +2726.21⨯(-6.17) ]
=5670.63kN m
工况四:
该工况上游水位与工况一的相同,其内水压力和对形心Y 轴产生的弯矩相同。
G 16.1=16171.47kN G 16.2=5748.42kN G 16.3=1624.90kN
G 16.4=10138.99kN
截面计算点至形心轴Y 的距离产生的弯矩
M
y 16
=40823.56kN m
工况五:
该工况为施工完建未挡水期,此时进水口中的内水压力为0。 2.2.18 地震惯性力
在此计算两个方向的地震惯性力,即顺水流方向和逆水流方向。 地震惯性力按下式计算:
F i =a h ξG Ei a i /g 式中:
F i ——作用在质点i a h
的水平方向地震惯性力代表值;
——水平向设计地震加速度代表值,根据四川省地震局工程地震研究院
(四川赛思特科技有限责任公司)《东义河某水电站坝址区设计地震反应谱参数表》(2009年11月),取a h =1.29m /s 2;
ξ——地震作用的效应折减系数,根据DL5073-2000《水工建筑物抗震设计
规范》,取ξ=0.25;
G E i ——集中在质点i
的重力作用标准值;
a i ——质点
i 的动态分布系数,根据DL 5073—2000《水工建筑物抗震设计
规范》表10.1.5来确定。
根据2.2.1~2.2.9计算结果, 查据DL 5073—2000《水工建筑物抗震设计规范》表10.1.5可得:
G Ei a i =28350⨯1.0+57263.85⨯1.38+7590.71⨯1.38+11257.97⨯2.03
+4582.32⨯1+630⨯1.0+556.5⨯1.0+4952.12⨯1.78-829.5⨯1.78=153810.04kN F i =a h ξG E i a i /g =
1.29⨯0.25⨯153810.04
9.81
=5056.45kN
地震惯性力顺水流方向力矩计算:
M =
-1.29⨯0.25
9.81
⨯(28350⨯2.61+57263.85⨯15.45+7590.71⨯15.45+11257.97
⨯19.70+4582.32⨯11.78+630⨯9.25+556.5⨯7.6+4952.12⨯18.05-829.5⨯18.05) =-47215.45kN m
地震惯性力逆水流方向力矩计算:
M =
1.29⨯0.25
9.81
⨯(28350⨯2.61+57263.85⨯15.45+7590.71⨯15.45+11257.97
⨯19.70+4582.32⨯11.78+630⨯9.25+556.5⨯7.6+4952.12⨯18.05-829.5⨯18.05) =47215.45kN m
地震动水压力按拟静力法和动力法两种方法计算,计算公式分别为:
F T =0.5a h ξρw ηw AH 0(
a 2H 0
)
-0.2
拟静力法:
其作用点位于水面以下0.42H0处 式中:
F T
——单位高度塔面动水压力合力代表值;
——水平向设计地震加速度根据四川省地震局工程地震研究院(四川赛
a h
思特科技有限责任公司)《古城水电站工程地震安全性评价复核表》(2008年8月),取a h =1.29m /s 2;
ξ——地震作用的效应折减系数,取ξ=0.25; ρw
——水体质量密度标准值,ρw =9.81kN /m 3;
——塔体形状系数,根据《水工建筑物抗震设计规范》(DL5073-2000)
ηw
表10.1.6-2,取ηw =1.048;
A
——塔体沿高度平均截面与水体交线包络面积,A =360.00m 2;
a ——塔体垂直地震作用方向的迎水面最大宽度沿高度的平均值,
a =18.0m
;
。
H 0——水深,H 0=21.00m F T =0.5a h ξρw ηw A H 0(
a 2H 0
)
-0.2
=0.5⨯1.29⨯0.25⨯9.81⨯1.048⨯360.00⨯21.00⨯(=14847.29kN
18.02⨯21
)
-0.2
地震动水压力顺水流方向力矩计算:
M =-14847.29⨯0.42⨯21.00=-130953.09kN m
地震动水压力逆水流方向力矩计算:
M =14847.29⨯0.42⨯21.00=130953.09kN m
2.3 稳定、地基应力及变形计算 2.3.1 抗剪断强度和抗剪强度计算
根据地质室提供的资料混凝土与基础的抗剪(断)指标建议值为: f /=0.4 c '=0.08M P a K '=
f '∑W +c 'A
∑
P
K f =
∑V ∑U
K 0
M ∑=
∑M
s 0
式中:K '——抗剪断强度计算的抗滑稳定安全系数; K f ——抗浮稳定安全系数; K 0——抗倾覆稳定安全系数;
f ' , c '——滑动面的抗剪断的摩擦系数及粘结力,kPa ; A ——基础面受压部分的计算截面积,m 2;
∑W ——全部荷载对滑动面的法向分值,包括扬压力,kN ; ∑V ——建基面上垂直力总和,包括扬压力(不含设备总重)kN ; ∑P ——全部荷载对滑动面的切向分值,kN ; ∑U ——基建面上扬压力总和,kN 。
∑M ∑M
s
——基建面上稳定力矩总和,kN m ; ——基建面上倾覆力矩总和,kN m 。
工况一:
∑W =77980.37kN
∑V
=159593. k 8N 7=82800.00kN
∑P =171. 5k 2N ∑U
f '∑W +'c A 0. 4⨯
=
∑P
77980+. 3⨯7171. 52
K '=求。 K f =工况二:
80360
=349. 7>8
,满足规范要3. 0
∑V
∑U
=
159593. 87
=1. 93>1. ,满足规范要求。10
82800. 00
∑W =104774.52kN
∑V
=129036. k 4N 2
∑P =159. 6k 4N ∑U
f '∑W +'c A 0. 4⨯
=
∑P
10477+4. 5⨯2159. 64
=25448. k 4N 0
3. 0,满足规范要
K '=求。 K f =工况三:
80360
=442. 9>4
∑V
∑U
=
129036. 42
=5. 07>1. ,满足规范要求。10
25448. 40
∑W =100255.39kN
∑V ∑U
=131861. k 2N 9=32792. k 4N 0 80360=449. 1>4
3. 0,满足规范要
∑P =153.41kN
f '∑W +'c A 0. 4⨯
= K '=
∑P
10025+5. 3⨯9153. 41
求。 K f =工况四:
该工况的受力情况与工况一相同,故其抗剪断系数和抗滑稳定系数也与工况一时的相等。 K '=349. 7>8工况五:
2K f =1.36>1.05,满足规范要求。 ,满足规范要求。
∑V
∑U
=
131861. 29
=4. 02>1. ,满足规范要求。10
32792. 40
∑W =118599.83kN
∑V ∑U
=117413. k 3N 3=0
80360=792. 5>8
∑P =96.19kN
K '=求。 K f =工况六:
f '∑W +'c A 0. 4⨯
=
∑P
11859+9. 8⨯3
96. 19
,满足规范要2. 5
∑V
∑U
=+∞1>. 0,满足规范要求。5
(1)顺水流方向:
∑W =283105.847kN ∑U =141588.000kN
∑P =20075.25kN
一时的相等。 K '=
K
∑V =418958.847kN
该工况的受力情况与工况一相同,故其抗剪断系数和抗滑稳定系数也与工况
f '∑W +'c A 0. 4⨯
=
∑P
77980+. 3⨯720075. 25
80360
=2. 9>9
2,满足规范要求。. 5
f
=1.93>1.05,满足规范要求。
(2)逆水流方向:
∑W =283105.847kN ∑U =141588.000kN
∑P =19732.22kN ∑V
一时的相等。
f '∑W +'c A 0. 4⨯
= K '=
∑P
K
f
=418958.847kN
该工况的受力情况与工况一相同,故其抗剪断系数和抗滑稳定系数也与工况
77980+. 3⨯719732. 22
80360
=3. 0>4
2,满足规范要求。. 5
=1.93>1.05,满足规范要求。
2.3.2 应力计算
σ=
∑V
A
±
∑M
I x
x
y
±
∑M
I y
x
x
σ——基建面上计算点垂直应力; ∑V ——基建面上垂直力总和;
∑M 、∑M
x
y
——分别为基建面上垂直力对形心轴X 、Y 轴的力矩总和;
X 、Y 轴的距离;
x 、y ——分别为基建面上垂直力对形心轴
I x 、I y ——分别为基建面对形心轴
X 、Y 轴的惯性矩;
A ——基建面面积。
该进水口呈对称布置,在任何工况下∑M x =0,所以σ工况一:
∑V =159593.87kN ∑M y =7582.47kN m
σ=
159593.87
360
±
7582.47⨯1012000
=
222.93kPa (建基面上游侧)210.29kPa (建基面下游侧)
V ∑=
A
M ∑±
I y
y
x
工况二:
∑V =129036.42kN ∑M y =-18442.03kN m
σ=
129036.42
360
±
-18442.03⨯10
12000
=
275.67kPa (建基面上游侧)306.41kPa (建基面下游侧)
工况三:
∑V =131861.29kN ∑M y =-13369.49kN m
σ=
131861.29
360
±
-13369.49⨯10
12000
=
267.35kPa (建基面上游侧)289.63kPa (建基面下游侧)
工况四:
该工况的受力情况与工况一相同,故基底应力也与工况一时相等。 σ=
222.93kPa (建基面上游侧)
210.29kPa (建基面下游侧)
工况五:
∑V =117413.33kN ∑M y =-18548.7kN m
σ=
117413.33
360
±
-18548.7⨯10
12000
=
313.99kPa (建基面上游侧)344.9kPa (建基面下游侧)
工况六:
(1)顺水地震时:
∑V =159593.87kN ∑M y =-170586.07kN m
σ=
159593.87
360
±
-170586.07⨯10
12000
=
74.46kPa (建基面上游侧)358.77kPa (建基面下游侧)
(2)逆水地震时:
∑V =159593.87kN ∑M y =185751.02kN m
σ=
159593.87
360
±
185751.02⨯10
12000
=
371.40kPa (建基面上游侧)61.82kPa (建基面下游侧)
2.3.3 稳定及应力分析
经计算,可以看出该进水口的抗滑稳定系数和抗浮稳定系数均大于所要求的安全系数,满足规范要求;各工况的基底应力未出现拉应力,并且都小于进水口段岩体承载力(420kP a ),因此可以认定该塔为稳定结构。 2.3.4 变形计算
由于进水口地基为非岩基基础,其地基最终沉降量利用压缩曲线计算:
n
S ∞=
∑
i =1
e 1i -e 2i 1+e 1i
i
式中
S ∞——地基最终沉降量; i ——土层号;
n ——地基压缩层计算深度范围内的土层数;
e 1i 、e 2i ——基础底面以下第i 层土在平均自重应力作用下的孔隙比和在平均自重应力、平均附加应力共同作用下的孔隙比; h i ——第i 层土的厚度。 1)基本假定
a )地基土压缩时不发生侧向变形。
b )按基础底面中心点下的附加应力计算土层的分层压缩量。 c )最终沉降量取计算深度范围内各土层分层压缩量之和。
2)地基压缩层的计算深度应按计算层面处附加应力和自重应力之比等于0.2
的条件确定。分层厚度取4.0m 。
a) 右上角点基地应力分解
根据2.3节中计算的应力将基底压力分为竖向均匀分布压力Pn 和三角形分布压力Pt 两部分, 见表2.1。由表中应力值可知,工况六(2)中下游侧应力小于建基面以上部分土层自重应力,故该工况左下角点的沉降量计算时,应由建基面向上进行计算,故该工况单独计算。
表2.1 右上角点基底应力分解表
222.93
371.40
74.46
图2.3-1 各工况应力分解简图
b )建基面以下各分层土自重应力计算
根据互提资料联系单数据覆盖层土容重为20.0kN/m3,水的容重为10.0kN/m3。以开挖面为参照面计算。计算简图如下:
图2.3-2 分层计算简图
σs 0=γd =(20.0-10.0)×4.0=40.0kP a σs 1=σs 0+γH 1=(20.0-10.0)×4.0+40=80.0kP a
σs 2=σs 1+γH 2=(20.0-10.0)×4.0+80.0=120.0kP a σs 3=σs 2+γH 3=(20.0-10.0)×4.0+120.0=160.0kP a
σs 4=σs 3+γH 4=(20.0-10.0)×4.0+160.0=200.0kP a σs 5=σs 4+γH 5=(20.0-10.0)×4.0+200.0=240.0kP a σs 6=σs 5+γH 6=(20.0-10.0)×1.4+240.0=254.0kP a
c )右上角点附加应力计算
将基底压力分为均匀分布和三角形分布两部分。计算示意图见图2.3-1。该基础为矩形,属空间问题,故应用“角点法”求解。长度L =20m , 宽度B =18m , 故m =
l b =2018
=1.111
表2.2 矩形均布荷载作用时右上侧角点下附加应力计算表
表2.3 三角形均布荷载作用时右上侧角点下附加应力计算表
表2.5 右上侧角点下各相邻层面应力差计算
上述表中
z ——基底以下计算深度
K 4——均布荷载时角点应力系数,查SL265-2001K 5——三角形分布荷载时角点应力系数,查
附录J 。
SL265-2001附录J 。
根据上表计算结果,基底以下24.0m 处附加应力和自重应力之比小于0.20,采用该处为计算深度的下限。
根据孔隙比与压力的关系曲线,采用TK2试样进行沉降计算。由公式
a v =1000
e 1-e 2p 2-p 1
,得
e 1i -e 2i =a v (p2i -p 1i ) ÷1000=0.001a v ∆p i
沉降计算结果见表2.6。
表2.6 右上侧角点基础沉降计算表
△p i ——第i 层 上、下层面应力值(自重应力与附加应力之和)差值 e 1i ——由第i 层自重应力平均值的孔隙比
e 2i ——由第i 层自重应力平均值与附加应力平均值之和的孔隙比 h i ——第i 分层土的厚度 △s i ——第i 分层土的沉降值 Σ△s i ——第i 分层土的累计沉降值
按上述相同方法,下游侧角点地基沉降计算如下:
表2.7 左下角点基底应力分解表
表2.8 矩形均布荷载作用时左下侧角点下附加应力计算表
表2.9 三角形均布荷载作用时左下侧角点下附加应力计算表
表2.10 左下侧角点下附加应力计算汇总
表2.11 左下侧角点下各相邻层面应力差计算
表2.12 左下侧角点基础沉降计算表
计算结果表明,上游侧在工况六(1)(顺水地震工况)下地基沉降最大为59.65mm ,下游侧在工况六(2)(逆水地震工况)下地基沉降最大为51.4mm ,沉
降差最大为工况六(1)为23.36mm (以开挖面为参考坐标)。 3. 结论
通过以上稳定与基底应力计算分析(计算成果见表3-1),可以认定该进水塔的稳定性能满足要求,压应力均小于基础的允许承载力(350kP a ),该进水塔的安全性可以保证。地基的沉降量也能满足要求,但进水口基础上、下游侧间地基面不均匀,可根据地质建议对地基进行处理(固结灌浆、布置振冲碎石桩等),以减少不均匀沉降值。
进水口地基计算最大沉降量59.65mm ,沉降差23.36mm (以开挖面为参照),参考《水闸设计规范》(SL265-2001)允许值,最大沉降量小于最大允许值(15cm ),沉降差小于最大允许值(5cm ),均满足规范要求。
表3-1稳定、应力计算成果表
某水电站 进水口整体稳定、基底应力及变形计算 第0页共32页
4.附图
荷 载 计 算 简 图
某水电站 进水口整体稳定、基底应力及变形计算 第0页共32页