进水口整体稳定及基底应力计算

进水口整体稳定及基底应力计算计算总说明1.1计算目的及要求某水电站进水口根据电站枢纽布置、地形、地质条件设为岸塔式进水口，“镶嵌”在L型基础中，塔背有基岩对其起支撑作用，靠自重和岸坡岩体支撑维持稳定，加之该进水口置于土质地基上，因建基面不允许出现拉应力，因此可不进行抗倾覆稳定计算。通过对进水口整体抗浮稳定与基底应力计算，以复核其是否满足规范要求。基本资料进水口纵横剖面结构尺寸见附图。水容重： 10kN/m3钢筋混凝土容重： 25kN/m3基础与混凝土之间f值为：广=0.4，C=008Mpa。地基承受能力： b]=0.42MPa校核洪水位： 1886.109m设计洪水位： 1884.069m正常蓄水位： 1900.000m死水位： 1896.000m多年平均最大风速： 11.00m/s吹程： 0.7km拦污栅及对应启闭机重：52t、8t事故检修门及对应启闭机重： 38t、15t计算采用规范及参考书：《水电站进水口设计规范》DL/T5398-2007《水工建筑物荷载设计规范》DL5077-1997《水工建筑物荷载设计规范》《水电站厂房设计规范》SL266-2001《公路桥涵设计通用规范》JTJ021-89《水闸设计规范》SL265-2001;《土力学》河海大学出版社出版，卢廷浩主编；某水电站首部枢纽平面布置图及主要建筑物剖面图；《孔隙比与压力关系曲线表》；1.3计算原理及假定岸塔式进水口根据自身的结构特点，只要基底应力在允许应力或允许抗力范围内，塔体就不会发生整体失稳，根据规范要求及进水口的布置情况，对塔体做以下假定：1）将塔体视为刚体，在荷载作用下，岩体变形产生抗力；2） 抗力或反力强度值采用材料力学法和刚体极限平衡法求得。3） 计算地基沉降量时，不考虑周边结构对进水口塔体的影响。此算稿按荷载组合｛基本组合、特殊组合｝计算。工况1:正常运行期，正常蓄水位1900.000m；工况2:正常运行期，设计洪水位1884.069m；工况3:非常运行期，校核洪水位1886.109m；工况4:检修期，正常蓄水位1900.000m；工况5:施工完建未挡水期。工况6:正常蓄水位+VII度地震工况1、2为基本组合，工况3、4、5、6为特殊组合。计算过程2.1截面计算取进水口底板进行计算，进水口底板顺水流方向为x轴，垂直水流方向为y轴，截面尺寸如图2.1-1所示。20.0图2.1-1建基面示意图基底面积：I=—xbh3=—x20x18.003=9720m4x12 12I=—xb3h=—x203x18.00=12000m4y12 12荷载作用正方向取为，垂直荷载作用：铅直向下为正；水平荷载作用：沿水流方向为正；弯矩作用：逆时针方向为正。力矩计算下面计算各截面形心对x、y轴的力矩，因为进水口左、右侧对称布置，所以只需计算进水口上下游基底应力，即可只算截面计算点对y轴的力矩即可。力矩计算简图见附图。2.2.1底板重量为：G=YYV=1.05x25x20.0x3.0x18.0=28350.0kN1Qcal=0x1M=Gl=0kNgmy1 1x1左、右边墙G=YYV=1.05x51.94x21.0x25x2=57263.85kN2Qcl=-1.19mx2截面计算点至形心轴Y的距离产生的弯矩M二Gl二57263.85x(-1.19)=—68143.98kNgmy2x2拦污栅中墩G二丫yV=1.05x13.77x21.0x25=7590.71kN3Qcl=5.84mx3截面计算点至形心轴Y的距离产生的弯矩M=Gl=7590.71x5.84=44329.76kNgmy3x32.2.4胸墙胸墙指的是进水口拦污栅之后的那部分实体混凝土，由于喇叭口上侧混凝土是一不规则的结构，在取其截面计算点时分两部分计算。G=yyV=1.05x34.31x12.5x25=11257.97kNQcaG=yyV=1.05x25x19.57x8.92=4582.32kNQcbl =3.25m、l =-0.31mx4.1 x4.2截面计算点至形心轴Y的距离产生的弯矩MM=Gl+Gl=11257.97x3.25+4582.32x(-0.31)=35167.88kNgmyy4a 4.1x4.1 4.2x4.22.2.5通气孔四周实体混凝土由于通气孔四周混凝土在计算边墙时已经计算了一部分，因此在此只计算流道上方混凝土(包括扣除部分)。G=yyV=1.05x25x11.94x15.8=4952.12kNQc扣除通气孔体积：G=yyV=1.05x25x2x15.8=829.5kNQcl =-8.28m l =-8.5mx5.1 x5.2截面计算点至形心轴Y的距离产生的弯矩M=Gl-Gl=4952.12x(-8.28)-829.5x(-8.5)=-33952.76kNgmy5 5.1x5.1 5.2x5.22.2.6进水口板梁顶部板梁G二丫yV二1.05x25x16.03x0.25x2二210.39kNQcG=yyV二1.05x25x37.4x0.25二245.44kNQcG=yyV=1.05x25x3.57x0.8二74.97kNTOC\o"1-5"\h\zQcG二yyV二1.05x25x0.25x17.85二117.14kNQcG=yyV=1.05x25x2.08x0.8二43.68kNQcl =8.73m、l=0.64m、l=-1.25m、l =-2.7m、l =-4.4mx6.1 x6.2 x6.3 x6.4 x6.5截面计算点至形心轴Y的距离产生的弯矩M=YGl=210.39x8.73+245.44x0.64+74.97x(—1.25)+117.14x(—2.7)+43.68x(-4.4)y6 6ix6i二1391.63kNgm2.2.7进水口启闭机排架将启闭机排架分为两部分计算，即拦污栅部分和闸门部分。G二yyV二1.05x25x79.74二2093.18kN(拦污栅排架)TOC\o"1-5"\h\zQcG二yyV二1.05x25x55.66二1461.18kN(事故检修门排架)Qcl =5.85m、l =一5・5mx7.1 x7.2截面计算点至形心轴Y的距离产生的弯矩M=YGl=2093.18x5.85+1461.18x(-5.5)=4209.16kNgmy7 7x7拦污栅及启闭机G二1.05x(52+8)x10二630.00kN8l=4.8mx8截面计算点至形心轴Y的距离产生的弯矩M二Gl二630.00x4.8二3024.00kNgmy8 8x8事故检修门及启闭机G二1.05x(38+15)x10二556.50kN9l=-5.8mx9截面计算点至形心轴Y的距离产生的弯矩M二Gl二556.5x(-5.8)=—3227.70kNgmy9 9x92.2.10人群荷载按《公路桥涵设计通用规范》人群荷载一般取3kN/m2。G=1.05x3x18.00x20.00=1134.00kN10l=0x10截面计算点至形心轴Y的距离产生的弯矩M=Gl=0kNgmy10 10x102.2.11静水压力根据岸塔式进水口的特点，各种工况时，其上下游水位相同，故上下游静水压力相互抵消，产生的弯矩和也为0，可不予计算。2.2.12扬压力根据《水电站厂房设计规范》，扬压力分浮托力和渗透压力，进水口上下游水位相等，故渗透压力为0。浮托力计算：工况一：U=1.0x10x(1900.00-1877.00)x18.00x20.00=82800.00kN(T)浮托力沿底板成对称分布，故产生的弯矩为0。工况二：U=1.0x10x(1884.069-1877.00)x18.00x20.00=25448.40kN(T)浮托力沿底板成对称分布，故产生的弯矩为0。工况三：U=1.0x10x(1886.109-1877.00)x18.00x20.00=32792.40kN(T)浮托力沿底板成对称分布，故产生的弯矩为0。工况四：U=1.0X10X(1900.00-1877.00)x18.00x20.00=82800.00kN(T)浮托力沿底板成对称分布，故产生的弯矩为0。工况五：施工完建未挡水期，此时浮托力为0。2.2.13风荷载垂直于建筑物表面上的风荷载标准值按下式计算：w=Bppwkzzso式中：w——风荷载标准值，kN/m2kB——z高程的风振系数，取1.0zp——风荷载体型系数，查《水工建筑物荷载设计规范》：sp=-(0.48+0.03H/B)sp——风压高度变化系数。zw——基本风压，查全国风压基本分布图可确定为0.3kN/m2o风荷载的作用分项系数采用1.3。w=Bppwkzzso工况一：风压高度变化系数，按《水工建筑物荷载设计规范》DL5077—1997表12.1.6查得：p=0.68ozp=-(0.48+0.03H/B)=-[0.48+0.03x(1902.00-1900.00)/18.00]=-0.483sw=Bppw=1.0x0.68x(-0.483)x0.3=-0.0986kN/m2kzzso风压力：w=1.3x(-0.0986)x18.00x(1902.00-1900.00)=-4.61kNl=-x(1902.00-1900.00)+(1900.00-1881.00)=20.33mx13 3截面计算点至形心轴Y的距离产生的弯矩M=Pl=-4.61x20.33=-93.83kNgmy13 13x13工况二：风压高度变化系数，按《水工建筑物荷载设计规范》DL5077—1997表12.1.6查得：r二1.204zr二-(0.48+0.03H/B)二-[0.48+0.03x(1902.00-1884.069)/18.00]二—0.5099sw=Brrw=1.0x1.204x(-0.5099)x0.3=-0.1842kN/m2kzzso风压力：w=1.3x(-0.1842)x18.00x(1902.00—1884.069)=-77.275kNl =-x(1902.00-1884.069)+(1884.069-1881.00)=15.02mx13 3截面计算点至形心轴Y的距离产生的弯矩M=Pl=—77.275x15.02=—1160.91kNgmy14 13x13工况三：风压高度变化系数，按《水工建筑物荷载设计规范》DL5077—1997表12.1.6查得：r=1.16zr=-(0.48+0.03H/B)=-[0.48+0.03x(1902.00—1886.069)/18.00]=-0.5065sw=Brrw=1.0x1.16x(-0.5065)x0.3=-0.1763kN/m2kzzso风压力：w=1.3x(-0.1763)x18.00x(1902.00-1886.109)=-65.541kNl =-x(1902.00-1886.109)+(1886.109-1881.00)=15.70mx13 3截面计算点至形心轴Y的距离产生的弯矩M=Pl=-65.541x15.7=-1029.19kNgmy13 13x13工况四：风压高度变化系数，按《水工建筑物荷载设计规范》DL5077—1997表12.1.6查得：r=0.68zr=-(0.48+0.03H/B)=-[0.48+0.03x(1902.00-1900.00)/18.00]=-0.483sw=Brrw=1.0x0.68x(-0.483)x0.3=-0.0986kN/m2kzzso风压力：w=1.3x(-0.0986)x18.00x(1902.00-1900.00)=-4.61kNl =-x(1902.00-1900.00)+(1900.00-1881.00)=20.33mx13 3截面计算点至形心轴Y的距离产生的弯矩M=Pl=-4.61x20.33二—93.83kNgmy13 13x13工况五：风压高度变化系数，按《水工建筑物荷载设计规范》DL5077—1997表12.1.6查得：r=1.267zr=-(0.48+0.03H/B)=-[0.48+0.03x(1902.00—1881.00)/18.00]=-0.515sw=Brrw=1.0x1.267x(-0.515)x0.3=-0.1958kN/m2kzzso风压力：w=1.3x(-0.1958)x18.00x(1902.00-1881.00)=-96.192kN113=2x(1902.00-1881.00)=14.00m截面计算点至形心轴Y的距离产生的弯矩M=Pl=-96.192x14.00=-1346.69kNgmy13 13x132.2.14浪压力官厅公式：gh=0.0076v-1/12(gD)1/3v2 0 v200■gLm=0.331V-1/2.15(gD)1/3.75v2 0 v200式中：h——当gD=20〜250时，为累计频率5%的波高；v20h——当gD=250〜1000时，为累计频率10%的波高；v20L 平均波长，m;mD 吹程，m;v——多年平均最大风速，m/s，按《水工建筑物荷载设计规范》DL50770—1997附录G表G1查得，正常蓄水位时取1.25v(风速系数由内插法求得)，0设计洪水位时取0.924v，校核洪水位时取0.949v。00临界水深的计算：TTL、L+2兀hH—―mlnm 1%cr4L—2兀hm 1%7兀h2 2兀Hh— +%cot-zLLmm式中：H——使波浪破碎的临界水深。crh——累积频率为1%时的波高，由《水工建筑物荷载设计规范》DL5077-19971%附录G表G2用插入法可以求得。h 波浪中心线至计算水位的高度，m。z浪压力标准值计算：当H>H，H<时，单位长度上的浪压力标准值计算公式为：cr2P—丄口+h)(丫H+P)+HP]wk2 1%zwlf lf当H>H，H>时，单位长度上的浪压力标准值计算公式为：cr2P——YL(h+h)wk4Wm1%z式中：P——单位长度迎水面上的浪压力标准值wkh 波浪中心线至计算水位的高度，m。zY——水的重度,kN/m3WP——建筑物底面处的剩余浪压力强度值，kN/m2，按下式计算：lfP—Yhsechlfw1% Lm浪压力的作用分项系数采用1.2工况一：也—0.0076v-i/i2(竺)i/3—0.0076x(1.25x11.00)-1/12x(9.81x700)1/3v2 0 v2 1.252x11.00200h—0.0076x(1.25x11.00)-1/12x(9.81x700)1/3x(1.25x11.00)2-9.81—0.39m1.252x11.002

v20二0.331V-1/2.15(竺)1/3.75二0.331X(1.25x11.00)-1/2.15X( 9,81x700 )v200 v2 1.252x11.002=0.331X(1.25x11.00)-1/2.15X=0.331X(1.25x11.00)-1/2.15X(9.81x7001.252x11.002)1/3.75x(1.25x11.00)2一9.81=4.91m临界水深H的计算cr由gD=9.81x700二36.32，比值在iD二20一250之间，为累计频率5%v2 (1.25x11.00)2 v200的波高；令齐二1^5，则hm=崗==0"m0.201900.00-1877.00二0.0087h根据二0.0087,查《水工建筑物荷载设计规范》DL5077-1997附录G表HmhG2用插入法可以求得t%二2.406hmh=2.406h=2.406x0.20=0.48m1% m”LL+2兀h 4.91|4.91+2x3.14x0.48门H=—mIn—m 1%= In =0.56mcr 4兀 L-2兀h 4x3.144.91-2x3.14x0.48m 1%7 兀h2 2兀H 3.14x0.482 2x3.14x(1900.00-1877.00)h= +%cot = cot =0.15mzLL 4.91 4.91mmH=1900.00-1877.00=23m>H=0.56mcrL4.91H=23m>—m= =2.46m，单位长度上的浪压力标准值计算公式为:22P=1yL(h+h)=1x10x4.91x(0.48+0.15)=7.73kN/mwk4Wm1%z4浪压力:P=1.2x7.73x18.00=166.9kN(T)wk计算点至形心的距离l=h+1900.00-1877.00=0.15+1900.00-1877.00=23.15mx14 z截面计算点至形心轴Y的距离产生的弯矩M=pl=-166.9x23.15=-3863.41kNgmy14 wkx14

工况二：譽=°.0076v-1/12(¥)1/3=0・0076x(0・924x11・00)一1/12x(。麗：70.』/300h=0-00766°924x11-00)-1/12x(0.9241x70.O02)1/3x(0-924x11W)29^=0267m号=。・3叫屮乎”75=。・331x（0・924川.00）-1/，15x（。置::霊”7500L—0.331x(0.924x11.00)-1/2.15x(9.81x700)1加5x(0.924x11.00)2-9.81—3.63mm 1.2042x11.002临界水深H的计算9.81x9.81x700由巴v20^（0.924x11.00）2=氏47，比值在畔=20"由巴v200的波高；人h令^%

hm—人h令^%

hm—1.95，则hm=毎=需=O'137,h —Hm0.1371884.069-1877.00=0.019h根据—0.019,查《水工建筑物荷载设计规范》DL5077-1997附录G表HmhG2用插入法可以求得—2.39hmh1%=2.39h—2.39x0.137=h1%mHcrL L+2兀h——mlnm 1%4HcrL L+2兀h——mlnm 1%4兀 L一2兀hm3.63|3.63+2x3.14x0.33 …— In —0.37m4x3.14 3.63-2x3.14x0.331%Pwk兀h2 2兀H+%cot-LLmm3」4x°332cot2x3・14x(1884.°69一1877.°0)—0.09m3.631884.069-1877.00—7.069m>H—0.37mcr3.637.069m> —363—1.83m，单位长度上的浪压力标准值计算公式为:22=—YL(h+h)—x10x3.63x(0.33+0.09)—3.81kN/m44Wm1%z浪压力：P二1.2x3.81x18.00二82.36kN(T)wk计算点至形心的距离l二h+1884.069-1877.00二0.09+1884.069-1877.00二7.16mx14 z截面计算点至形心轴Y的距离产生的弯矩M二pl =—82.36x7.16=—589.86kNgmy14 wkx14工况三：吕=0・0076v—1/12(乎)1/3=。・0076x(。・949x^00)-1/12x(。麗：需02)1/300h=0.0076x(0.949x11.00)-1/12x(9.81x700)1/3x(0.949x11.00)2-9.81=0.276m1.162x11.002gLm=0.331v-1/2.15(gD)1/3.75=0.331x(0.949x11.00)-1/2.15x( 9.81x700)1/375v2 0 v2 0.9492x11.00200L—0.331x(0.949xll.00)-1/2」5x( 9.81x700 )1/3.75x(0.949x11.00)2-9.81=3.73mm 0.9492x11.002，临界水深H的计算cr由gD= 9.81x700 =63.02，比值在gD=20一250之间，为累计频率5%v2 (0.949x11.00)2 v200的波高；h h0.276令一5%=1.95，贝Uh=—5%= =0.14,h m1.95 1.95mh m—h m—Hm0.141886.109-1877.00=0.016h根据=0.016,查《水工建筑物荷载设计规范》DL5077-1997附录G表HmhG2用插入法可以求得t%=2.395hmh=2.395h=2.395x0.14=0.34m1% m” L L+2兀h3.73, 3.73+2x3.14x0.34 门H=—mIn—m 1%= In =0.39mcr 4兀 L-2兀h 4x3.14 3.73-2x3.14x0.34m 1%’兀h2 2兀H 3.14x0.342 2x3.14x(1886.109-1877.00)「h= 1%cot= cot =0.1mzLL 3.73 3.73mmH=1886.109-1877.00=9.109m>H=0.39mcrL 3.73H=9.109m>f= =1.87m，单位长度上的浪压力标准值计算公式为:22P=1YL(h+h)=1xioX3.73x(0.34+0.1)=4.07kN/mwk4Wm1%z4浪压力：P二1.2x4.07X18.00二87.87kN(T)wk计算点至形心的距离l二h+1886.109-1877.00二0.1+1886.109-1877.00二9.21mx14 z截面计算点至形心轴Y的距离产生的弯矩M=pl—-87.87x9.21=—808.92kNgmy14 wkx14工况四：该工况时的浪压力等于工况一时的浪压力，截面计算点至形心的轴Y的距离产生的弯矩也与工况一的相等。浪压力：P—166.9kN(T)wk截面计算点至形心轴Y的距离产生的弯矩M —-3863.41kNgmy14工况五：施工完建未挡水期，浪压力为0,截面计算点至形心轴丫的距离产生的弯矩也为0。2.2.15胸墙及检修门槽中的内水压力工况一：G—1.05x10x44.66x(1900.00-1888.954)—5179.80kN15.1G—1.05x10x18.58x(1900.00-1888.454)—2252.51kN15.2G—1.05x10x7.92x(1900.00-1887.200)—1064.45kN15.3计算点至形心的距离l —0.3m、l —-3.36m、l —-5.7mx15.1 x15.2 x15.3截面计算点至形心轴丫的距离产生的弯矩M—pl—5179.80x0.3+2252.51x(—3.36)+1064.45x(—5.7)——12081.84kNgmy15 wkx15工况二：设计洪水位时，胸墙和门槽内无水，故内水压力为0。工况三：校核洪水位时，胸墙和门槽内无水，故内水压力为0。工况四:该工况上游水位与工况一时相同，故内水压力和其对形心轴产生的弯矩均相同。G二5179.80kN15.1G=2252.51kN15.2G=1064.45kN15.3计算点至形心的距离l=0.3m、l=-3.36m、l =-5.7mx15.1 x15.2 x15.3截面计算点至形心轴Y的距离产生的弯矩M=pl=-12081.84kNgmy15 wkx15工况五：施工完建未挡水期，此时进水口中的内水压力为0。2.2.16进水口中的内水压力G=yy(V+V+V+V)16Qc1 2 3 4式中：Y厂可变作用的分项系数；Q—拦污栅槽(包括栏污栅槽前)中的水体积；1—喇叭口中的水体积；2—检修闸门槽中的水体积；3—边墩两侧的水体积。4工况一：G=yyV=1.05x10x81.06x(1900.00—1881.00)=16171.47kNQc1G=yyV=1.05x10x95.21x5.75=5748.42kNQc2G=yyV=1.05x10x29.76x5.2=1624.90kNQc3G二丫丫V二1.05x10x25.41x(1900.00-1881.00)x2二10138.99kNQc4截面计算点至形心轴Y的距离为l =6.84m、l=0.78m、l =-7.21m、l =-6.17mx16.1 x16.2 x16.3 x16.4截面计算点至形心轴Y的距离产生的弯矩M=yy(Vl +Vl+Vl+Vl)y16 Qc1x16.1 2x16.2 3x16.3 4x16.4=1.05x10x[16171.47x6.84+5748.42x0.78+1624.90x(-7.21)+10138.99x(—6.17)]二40823.56kNgm工况二：G二yyV二1.05x10x81.06x(1884.069-1881.00)二2612.12kNQc1G=yyV=1.05x10x95.21x5.75二5748.42kNQc2G=yyV=1.05x10x29.76x5.2二1624.90kNQc3G二yyV二1.05x10x25.41x(1884.069-1881.00)x2二1637.65kNQc4截面计算点至形心轴Y的距离为：l =6.84m、l=0.78m、l =-7.21m、l =-6.17mx16.1 x16.2 x16.3 x16.4截面计算点至形心轴Y的距离产生的弯矩M二yy(Vl +Vl+Vl+Vl)y16 Qc1x16.1 2x16.2 3x16.3 4x16.4=1.05x10x[2612.12x6.84+5748.42x0.78+1624.90x(-7.21)+1637.65x(—6.17)]二510.74kNgm工况三：G二yyV二1.05x10x81.06x(1886.109-1881.00)二4348.42kN16.1Qc1G二yyV二1.05x10x95.21x5.75二5748.42kN16.2Qc2G二yyV二1.05x10x29.76x5.2二1624.90kNQc3G二yyV二1.05x10x25.41x(1886.109-1881.00)x2二2726.21kNQc4截面计算点至形心轴Y的距离为：l =6.84m、l=0.78m、l =-7.21m、l =-6.17mx16.1 x16.2 x16.3 x16.4截面计算点至形心轴Y的距离产生的弯矩M二yy(Vl +Vl+Vl+Vl)y16 Qc1x16.1 2x16.2 3x16.3 4x16.4=1.05x10x[4348.42x6.84+5748.42x0.78+1624.90x(-7.21)+2726.21x(-6.17)]二5670.63kNgm工况四：该工况上游水位与工况一的相同，其内水压力和对形心Y轴产生的弯矩相同。G二16171.47kN16.1G=5748.42kN16.2G=1624.90kN16.3G =10138.99kN16.4截面计算点至形心轴Y的距离产生的弯矩M=40823.56kNgmy16工况五：该工况为施工完建未挡水期，此时进水口中的内水压力为0。2.2.18地震惯性力在此计算两个方向的地震惯性力，即顺水流方向和逆水流方向。地震惯性力按下式计算：F=agGa/gihEii式中：F——作用在质点i的水平方向地震惯性力代表值；ia——水平向设计地震加速度代表值，根据四川省地震局工程地震研究院h（四川赛思特科技有限责任公司）《东义河某水电站坝址区设计地震反应谱参数表》（2009年11月），取a=1.29m/s2；hg——地震作用的效应折减系数，根据DL5073-2000《水工建筑物抗震设计规范》，取g=0.25；G——集中在质点i的重力作用标准值；Eia——质点i的动态分布系数，根据DL5073—2000《水工建筑物抗震设计i规范》表10.1.5来确定。根据2.2.1〜2.2.9计算结果，查据DL5073—2000《水工建筑物抗震设计规范》表10.1.5可得:

Ga二28350x1.0+57263.85x1.38+7590.71x1.38+11257.97x2.03Eii+4582.32x1+630x1.0+556.5x1.0+4952.12x1.78-829.5x1.78二153810.04kN=agGa/=agGa/ghEii=5056.45kN911地震惯性力顺水流方向力矩计算:M=-x x(28350x2.61+57263.85x15.45+7590.71x15.45+11257.979.81x19.70+4582.32x11.78+630x9.25+556.5x7.6+4952.12x18.05-829.5x18.05)=-47215.45kNgm地震惯性力逆水流方向力矩计算:x0.25M=二—x(28350x2.61+57263.85x15.45+7590.71x15.45+11257.979.81x19.70+4582.32x11.78+630x9.25+556.5x7.6+4952.12x18.05-829.5x18.05)=47215.45kNgm地震动水压力按拟静力法和动力法两种方法计算，计算公式分别为：aF=0.5agpnAH(——)-0.2拟静力法：Thww02H0其作用点位于水面以下0.42H0处式中：F——单位高度塔面动水压力合力代表值；Ta——水平向设计地震加速度根据四川省地震局工程地震研究院（四川赛h思特科技有限责任公司）《古城水电站工程地震安全性评价复核表》（2008年8月），取a=1.29m/s2;hg——地震作用的效应折减系数，取g=0.25;水体质量密度标准值，p=9.81kN/m3；w塔体形状系数，根据《水工建筑物抗震设计规范》（DL5073-2000）表10.1.6-2，取n=1.048；wA——塔体沿高度平均截面与水体交线包络面积，A=360.00m2；a——塔体垂直地震作用方向的迎水面最大宽度沿高度的平均值

a二18.0m;H 水深，H=21.00m。00F=0.5agpnAH(丄)-0.2Thww02H018.0=0.5x1.29x0.25x9.81x1.048x360.00x21.00x( )-0.22x21=14847.29kN地震动水压力顺水流方向力矩计算:M二-14847.29x0.42x21.00二-130953.09kNgm地震动水压力逆水流方向力矩计算:M二14847.29x0.42x21.00二130953.09kNgm稳定、地基应力及变形计算2.3.1抗剪断强度和抗剪强度计算根据地质室提供的资料混凝土与基础的抗剪断）指标建议值为：根据地质室提供的资料混凝土与基础的抗剪断）指标建议值为：ff/二0.4f'工W+c'AK=工Pc'=0.08MPa工VK=—f乙U工MK=s0乙M0式中：K'——抗剪断强度计算的抗滑稳定安全系数；抗浮稳定安全系数；K——抗倾覆稳定安全系数；0f',c' 滑动面的抗剪断的摩擦系数及粘结力，kPa;A 基础面受压部分的计算截面积，m2；工W——全部荷载对滑动面的法向分值，包括扬压力，kN；工V——建基面上垂直力总和，包括扬压力（不含设备总重）kN;工P——全部荷载对滑动面的切向分值，kN;工U——基建面上扬压力总和，kN。工M 基建面上稳定力矩总和，kNgm;s工M0 基建面上倾覆力矩总和，kNgm。

工况一：工W=77980.37kN工V=159593.87kN求。工况二：工P=171.52kN工U工况一：工W=77980.37kN工V=159593.87kN求。工况二：工P=171.52kN工U=82800.00kNK'=广gFA=°4x77980・37780X360=349.78>3.0,满足规范要171.52工V15959387KP=莎0000=1.93>1.10，满足规范要求。工W=104774.52kN工V=129036.42kN工P=159.64kN求。工况三：159.64工V12903642h二玉吋=5-07>x-10，满足规范要求。工W=100255.39kN工V=131861.29kN工P=153.41kN工U求。工况三：159.64工V12903642h二玉吋=5-07>x-10，满足规范要求。工W=100255.39kN工V=131861.29kN工P=153.41kN工U=32792.40kN求。求。153.41K<=广¥TA=°4X10025]39「80X360=449.14>3.0,满足规范要153.41工V13186129KfP=远丽40=4.02>1.10，满足规范要求。工况四：该工况的受力情况与工况一相同，故其抗剪断系数和抗滑稳定系数也与工况一时的相等。K'二349.78>2.5，满足规范要求。二1.36>1.05，满足规范要求。工况五：

工W=118599.83kN工V工W=118599.83kN工V=117413.33kN求。96.19K'=广¥FA=04小8599.83+80X360=792.58>2.5，满足规范要求。96.19工VK=Zu十>1.05，满足规范要求。工况六：1）顺水流方向：工W=283105.847kN工U=141588.000kN工P=20075.25kN工V=418958.847kN该工况的受力情况与工况一相同，故其抗剪断系数和抗滑稳定系数也与工况一时的相等。f工W+CA0.4x77980.37+80x360K=Z-20075.25=2.99>2.5，满足规范要求。K二1.93>1.05，满足规范要求。f2）逆水流方向：工W=283105.847kN工U=141588.000kN工P=19732.22kN工V=418958.847kN该工况的受力情况与工况一相同，故其抗剪断系数和抗滑稳定系数也与工况一时的相等。f工W+c'A0.4X77980.37+80x360 屮卄岳亠K'= = =3.04>2.5，满足规范要求。乙P 19732.22K二1.93>1.05，满足规范要求。f2.3.2应力计算工V工My工MxQ— + x—+ X—A~I -Ixy

b 基建面上计算点垂直应力；工V——基建面上垂直力总和；工M、工M——分别为基建面上垂直力对形心轴X、丫轴的力矩总和;xyx、y——分别为基建面上垂直力对形心轴X、丫轴的距离；I、I——分别为基建面对形心轴X、Y轴的惯性矩；xyAMxA 基建面面积。AMx该进水口呈对称布置，在任何工况下工M=0,所以bx工况一：工V=159593.87kN 工M=7582.47kNgmy7582.47x10 222.93kP（建基面上游侧）b= + =360 12000 210.29kPa（建基面下游侧）工况二：工V=129036.42kN 工M=-18442.03kNgmy129036.42-18442.03x10 275.67kPa（建基面上游侧）b= 土 =360 12000 306.41kPa（建基面下游侧）工况三：工V=131861.29kN 工M=-13369.49kNgmy-13369.49x10267.35kPa（建基面上游侧）b= + =360 12000 289.63kPa（建基面下游侧）工况四：该工况的受力情况与工况一相同，故基底应力也与工况一时相等222.93kPa（建基面上游侧）b=210.29kPa（建基面下游侧）工况五：工V=117413.33kN 工M=-18548.7kNgmy117413.33-18548.7x10 313.99kPa（建基面上游侧）b= + =360 12000 344.9kPa（建基面下游侧）工况六：（1）顺水地震时：工V=159593.87kN 工M=—170586.07kNgmy—170586.07x10 74.46kP（建基面上游侧）G二 土 二360 12000 358.77kPa（建基面下游侧）（2）逆水地震时：工V=159593.87kN 工M=185751.02kNgmy185751.02x10 371.40kPa（建基面上游侧）G= + =360 12000 61.82kPa（建基面下游侧）2.3.3稳定及应力分析经计算，可以看出该进水口的抗滑稳定系数和抗浮稳定系数均大于所要求的安全系数，满足规范要求；各工况的基底应力未出现拉应力，并且都小于进水口段岩体承载力（420kPa），因此可以认定该塔为稳定结构。2.3.4变形计算由于进水口地基为非岩基基础，其地基最终沉降量利用压缩曲线计算：S=£e1i—e2ihs 1+e 1i=1 1i式中S——地基最终沉降量；OOi——土层号；n——地基压缩层计算深度范围内的土层数；e、e——基础底面以下第i层土在平均自重应力作用下的孔隙比和在平均1i2i自重应力、平均附加应力共同作用下的孔隙比；h——第i层土的厚度。i1）基本假定a）地基土压缩时不发生侧向变形。b）按基础底面中心点下的附加应力计算土层的分层压缩量。c）最终沉降量取计算深度范围内各土层分层压缩量之和。2）地基压缩层的计算深度应按计算层面处附加应力和自重应力之比等于0.2

的条件确定。分层厚度取4.0m。a）右上角点基地应力分解根据2.3节中计算的应力将基底压力分为竖向均匀分布压力Pn和三角形分布压力Pt两部分,见表2.1。由表中应力值可知，工况六（2）中下游侧应力小于建基面以上部分土层自重应力，故该工况左下角点的沉降量计算时，应由建基面向上进行计算，故该工况单独计算。表2.1右上角点基底应力分解表工况四工况四工况工况一工况二工况三工况四工况五工况六（1）工况六（2）上游侧应力（kpa）222.93275.67267.35222.93313.9974.46371.40下游侧应力（kpa）210.29306.41289.63210.29344.90358.7761.82竖直均布压力Pn（kpa）222.93275.67267.35222.93313.9974.46371.40O—Omin so三角形分布压力Pt（kpa）12.6430.7422.2812.6430.91284.31309.59p-pmax min图2.3-1各工况应力分解简图b）建基面以下各分层土自重应力计算根据互提资料联系单数据覆盖层土容重为20.0kN/m3，水的容重为10.0kN/m3。以开挖面为参照面计算。计算简图如下：4.0m4.0m4.0m4.0m4.0m4.0m1图2.3-2分层计算简图c=丫d=(20.0—10.0)X4.0=40.0kPas0c=c+丫H=(20.0—10.0)X4.0+40=80.0kPaTOC\o"1-5"\h\zs1 s0 1c=c+yH=(20.0—10.0)X4.0+80.0=120.0kPas2 s1 2c=c+yH=(20.0—10.0)X4.0+120.0=160.0kPas3 s2 3c=c+yH=(20.0—10.0)X4.0+160.0=200.0kPas4 s3 4c=c+yH=(20.0—10.0)X4.0+200.0=240.0kPas5 s4 5c=c+yH=(20.0—10.0)X1.4+240.0=254.0kPas6 s5 6右上角点附加应力计算将基底压力分为均匀分布和三角形分布两部分。计算示意图见图2.3—1。该基础为矩形，属空间问题，故应用“角点法”求解。长度L=20m,宽度B=18m,l20故m=-= =1.111b18表2.2矩形均布荷载作用时右上侧角点下附加应力计算表Z(m)n=z/Bm=L/B1K4az1=KPn(kpa)4工况一工况二工况三工况四工况五工况六(1)工况六(2)00.0001.110.25055.7368.9266.8455.7378.5018.6192.8540.2221.110.24855.3668.4566.3955.3677.9718.4992.2380.4441.110.23953.1765.7563.7653.1774.8917.7688.58120.6671.110.21948.8760.4358.6048.8768.8316.3281.41160.8891.110.19443.3253.5751.9543.3261.0114.4772.17201.1111.110.16837.5346.4145.0137.5352.8612.5462.5324.01.3331.110.14432.1539.7538.5532.1545.2810.7453.5625.41.4111.110.13630.4137.6136.4730.4142.8410.1650.67

表2.3三角形均布荷载作用时右上侧角点下附加应力计算表Z(m)n=z/Bm=L/B1K5ez2=Kfn(kpa)工况一工况二工况三工况四工况五工况六（1）工况六（2）00.0001.110.0000.000.000.000.000.000.000.0040.2221.110.0330.421.030.750.421.039.5110.3680.4441.110.0570.721.761.280.721.7716.3117.76120.6671.110.0690.872.121.530.872.1319.5721.31160.8891.110.0710.892.171.570.892.1820.0721.85201.1111.110.0670.842.051.490.842.0618.9820.6724.01.3331.110.0610.761.861.350.761.8717.2118.7325.41.4111.110.0580.731.791.300.731.8016.5318.00表2.4右上侧角点下附加应力计算汇总项目层厚(m)az(kpa)0.2acz(kPa)工况一工况二工况三工况四工况五工况六（1）工况六（2）az1-az2az1+az2az1+az2oz1-0z2az1+az2az1+az2oz1-0z2基础以下深度Z=0m055.7368.9266.8455.7378.5018.6192.850Z=4m454.9369.4867.1354.9379.0028.0081.878Z=8m452.4567.5165.0452.4576.6634.0770.8216Z=12m448.0062.5460.1448.0070.9535.8960.1124Z=16m442.4355.7453.5242.4363.1934.5350.3232Z=20m436.6948.4646.5036.6954.9231.5241.8640Z=24m431.3841.6139.9031.3847.1527.9434.8248表2.5右上侧角点下各相邻层面应力差计算项目层厚（m）△P(kPa)i工况一工况二工况三工况四工况五工况六（1）工况六（2）基础以下深度Z=0m055.7368.9266.8455.7378.5018.6192.85Z=4m494.93109.48107.1394.93119.0068.00121.87Z=8m4132.45147.51145.04132.45156.66114.07150.82Z=12m4168.00182.54180.14168.00190.95155.89180.11Z=16m4202.43215.74213.52202.43223.19194.53210.32Z=20m4236.69248.46246.50236.69254.92231.52241.86Z=24m4271.38281.61279.90271.38287.15267.94274.82上述表中z 基底以下计算深度均布荷载时角点应力系数，查L265-2001附录J。三角形分布荷载时角点应力系数，查SL265-2001附录J。根据上表计算结果，基底以下24.0m处附加应力和自重应力之比小于0.20,采用该处为计算深度的下限。根据孔隙比与压力的关系曲线，采用TK2试样进行沉降计算。由公式a二1000ei-e2，得v P-P21e一e=a(p一p)十1000=0.001aAp1i2iv2i1i vi沉降计算结果见表2.6。表2.6右上侧角点基础沉降计算表项目hi(mm)av(MPa-i)eiiVs=eii-e2ihi1+e i1i工况一工况二工况三工况四工况五工况六(1)工况六(2)基础以下深度Z=4m40000.0800.359.2929.6159.5519.2929.60111.7076.878Z=8m40000.0800.348.9589.0829.0538.9588.99311.0016.915Z=12m40000.0800.348.4908.3668.3818.4908.1899.9876.993Z=16m40000.0800.338.2847.9868.0328.2847.7579.2987.269Z=20m40000.0800.338.2437.8747.9348.2437.6358.8987.589Z=24m40000.0800.338.3477.9768.0378.3477.7538.7647.931^△s(mm)i51.61450.89950.98951.61449.92859.6543.574表中△p 第i层上、下层面应力值(自重应力与附加应力之和)差值ie——由第i层自重应力平均值的孔隙比1ie——由第i层自重应力平均值与附加应力平均值之和的孔隙比2ih——第i分层土的厚度i△s——第i分层土的沉降值i^△s——第i分层土的累计沉降值i按上述相同方法，下游侧角点地基沉降计算如下：

表2.7左下角点基底应力分解表工况工况一工况二工况三工况四工况五工况六(1)工况六(2)上游侧应力(kpa)222.93275.67267.35222.93313.9974.46371.40下游侧应力(kpa)210.29306.41289.63210.29344.90358.7761.82竖直均布压力Pn(kpa)210.29306.41289.63210.29344.90358.7761.82O—Omin so三角形分布压力Pt(kpa)12.6430.7422.2812.6430.91284.31309.59p-pmax min表2.8矩形均布荷载作用时左下侧角点下附加应力计算表Z(m)n=z/Bm=L/B1K4az1=KPn(kpa)4工况一工况二工况三工况四工况五工况六(1)工况六(2)00.0001.110.25052.5776.6072.4152.5786.2389.6915.4540.2221.110.24852.2276.0971.9252.2285.6489.0915.3580.4441.110.23950.1673.0869.0850.1682.2685.5714.74120.6671.110.21946.1067.1763.4946.1075.6078.6413.55160.8891.110.19440.8659.5456.2840.8667.0269.7112.01201.1111.110.16835.4051.5948.7635.4058.0760.4010.4124.01.3331.110.14430.3244.1841.7630.3249.7351.738.9125.41.4111.110.13628.6941.8039.5128.6947.0648.958.43表2.9三角形均布荷载作用时左下侧角点下附加应力计算表Z(m)n=z/Bm=L/B1K5az2=Kfn(kpa)工况一工况二工况三工况四工况五工况六(1)工况六(2)00.0001.110.0000.000.000.000.000.000.000.0040.2221.110.0330.421.030.750.421.039.5110.3680.4441.110.0570.721.761.280.721.7716.3117.76120.6671.110.0690.872.1