《[精品论文]泰安抽水蓄能电站水利枢纽--------上水库库盆、材料及导流建筑物设计计算书.doc》由会员分享，可在线阅读，更多相关《[精品论文]泰安抽水蓄能电站水利枢纽--------上水库库盆、材料及导流建筑物设计计算书.doc（43页珍藏版）》请在三一文库上搜索。
1、课件之家精心整理资料-欢迎你的欣赏 课件之家精心整理资料-欢迎你的欣赏 目录目录 目录 .- 1 - 第一章 坝体计算 - 1 - 1.1 防浪墙顶高程及坝顶高程确定 .- 2 - 1.1.1 防浪墙顶高程确定.- 2 - 1.1.2 防浪墙底高程的确定.- 4 - 1.1.3 坝顶高程的确定.- 4 - 1.1.4 坝顶面宽度的确定.- 4 - 1.2 防浪墙应力稳定计算及配筋计算 .- 5 - 1.2.1 防浪墙应力稳定计算.- 5 - 1.2.2 防浪墙配筋计算.- 7 - 1.3 面板的计算 .- 8 - 1.3.1 面板的厚度计算.- 8 - 1.3.2 面板的配筋计算.- 8 -
2、1.4 结论 .- 9 - 第二章 趾板计算 - 9 - 21 趾板剖面尺寸 - 9 - 22 趾板配筋 .- 12 - 23 结论 .- 13 - 第三章 坝体稳定和变形计算 - 13 - 31 边坡稳定计算 .- 13 - 311 计算公式- 13 - 边坡稳定计算的 FORTRAN 语言程序 - 16 - 3.2 变形计算 - 38 - 3.3 结论- 40 - 第四章 导流隧洞的计算 .- 40 - 4.1 进水口体型计算- 40 - 4.2 洞身断面尺寸- 42 - 第一章坝体计算 泰安抽水蓄能电站水利枢纽-上水库库盆及导流建筑物设计 - 2 - 1.1 防浪墙顶高程及坝顶高程确定
3、1.1.1 防浪墙顶高程确定 防浪墙顶高程由水库静水位加波浪爬高,壅高及安全超高决定，坝顶高程 计 算分别考虑正常情况和非常情况，取所得坝顶高程的较大值。即 坝顶高程=设计洪水位+h设 坝顶高程=校核洪水位+h校 其中，正常情况下坝顶高程应高于静水位 0.5m 以上，非常情况下不低于静水 位。 图 1-1 坝顶超高计算简图 由于坝顶设防浪墙,超高计算采用重力坝超高计算公式(图 1-1) （1-1）hchhh 01 2 式中：h水库静水位以上的超高（m） 2h1累计频率 1%的波浪高度（m） h0 波浪中心线高出静水位的高度（m） hc 取决于坝的级别和计算情况的安全超高（m）,由表 1-1 确
4、 定。 表 1-1 土坝坝顶安全超高值（m） 坝的级别 运用情况 123 4，5 正常 1.51.00.70.5 非常 0.70.50.40.3 本工程坝的级别为级，故安全超高正常情况下取 1.5m,非常情况下为 0.7m 由于计算风速200 最小超高值(m) 234 57 本设计坝高在 50100m,根据表 1-2,最小超高值为 3m,防浪墙顶高程比正常 蓄水位高出 3.8m,3.0m,满足地震区安全超高的要求.故防浪墙顶高程取为 414.35m. 1.1.2 防浪墙底高程的确定 根据混凝土面板堆石坝设计规范DL/T5016-1999 中的 5.2.2 规定:防浪墙 顶高出坝顶 11.2m,
5、防浪墙与面板顶部的接缝高程即防浪墙的底面高程宜高于水 库正常蓄水位。但对于本坝，坝高中等，现在的止水技术又可以承受一定的水头， 为节省坝体的填筑方量，综合分析，取防浪墙高 3.8m,底面高程为 410.55m,与正 常蓄水位齐平.防浪墙具体尺寸见图 1-2. 图 1-2 防浪墙剖面图 1.1.3 坝顶高程的确定 根据混凝土面板堆石坝设计规范DL/T5016-1999 中的 5.2.2 规定:防浪墙 顶高出坝顶 1.01.2m。 现取防浪墙高出坝顶 1.2m,则坝顶高程为：413.95-1.2=412.75m. 1.1.4 坝顶面宽度的确定 根据混凝土面板堆石坝设计规范SL228-98 中的第
6、5.1.1 条规定：坝顶宽 度应由运行布置坝顶设施和施工的要求确定，亦按照坝高不同采用 58m，100m 课件之家精心整理资料-欢迎你的欣赏 课件之家精心整理资料-欢迎你的欣赏 以上的高坝宜适当加宽，如坝顶有交通要求时，坝顶的宽度还用遵照有关规定选 用。本工程中最小坝顶宽度按公式：Bmin=0.1h 计算。本工程最大坝高约为 98.0m，所以 Bmin=0.1 98=9.8m。因考虑交通要求，现取坝顶净宽 9.5m，上游防 浪墙厚 0.5m，则坝顶总宽为 10.0m。 1.2 防浪墙应力稳定计算及配筋计算 1.2.1 防浪墙应力稳定计算 将防浪墙简化为等厚度的“L”型墙，作用在其上荷载有：自重
7、、土压力。 水工建筑物教材提到,欲使填土发生被动破坏，挡土墙位移量需达到墙高的 5%10%。这是工程中所不允许的，所以计算不考虑被动土压力情况。同时,假设 防浪墙发生向上游的位移,此时应考虑主动土压力。下面分别讨论不蓄水状况、 设计洪水状况和校核洪水位三种情况。 考虑防浪墙最危险截面为竖墙与底版交界面，取竖墙部分，一端固定，相当 于悬臂梁结构，受力分析图见图 1-3。 G 力力力Ea 图 1-3 防浪墙应力计算受力分析图 参考土力学,压实填土取 Ko=1.3；由水工建筑物知，主动、静止土 压力的作用分项系数皆为 1.2;又由于 Ka=tg (45-/2)K=1.6 由于抗倾弯矩没包括重力,故所
8、得系数已偏于安全,因而防浪墙不会发生倾覆 破坏。 1.2.2 防浪墙配筋计算 对 L 型挡墙的竖直部分简化为悬臂梁计算 最不利荷载组合为不蓄水情况，静止土压力作用。 安全级别级，o=1.1 短暂状况，=1.0 故最大弯矩设计值 M=58.278KNm 0 1.1 1.0 52.98M 环境级别为：三类 保护层厚度 a=30mm h0=h-a=500-30=470mm 取单位宽度 1m 进行计算，混凝土采用 C25，则轴心抗压强度设计值 。钢筋采用级钢筋， 2 / 5 . 12mmNfc 2 /103mmNfy 截面抵抗矩系数： （2-8） 2 0 hbf M c d s 6 2 1.2 58.
9、278 10 0.0025 12.5 1000 (50030) （2-9） s 211 =0.025 0.9,故设计合理满足要求。 22 趾板配筋 根据混凝土面板堆石坝设计规范DL/T 50161999 要求，趾板宜采用单 层双向钢筋，每向配筋率采用 0.3%0.4%。岩基上趾板钢筋的保护层厚度为 10cm15cm。趾板应用锚筋与基岩连接，锚筋参数可按经验确定。 2.2.1 连接板 连接板混凝土采用 C25，连接板混凝土为双层双向配筋，每向配筋率 0.4%。 参照设计报告，选用 22 钢筋，连接板混凝土含钢量约 106kg/m。 2.2.2 趾板 本工程设计趾板混凝土采用 C25，按配温度筋考
10、虑，采用单层双向钢筋，布 置于趾板表面，每向配筋率 0.4%，顶面保护层 15cm。 As=bh=0.4%1000500=2000mm 式中: As钢筋面积，mm2； 配筋率； b 板单宽，mm； h0 有效宽度，因趾板相对较薄，故计算配筋时 h0仍采用其设 计厚度 h 来考虑。 趾板纵横筋均采用 22 200mm 的钢筋（As=1902mm，小于 5%少筋控制范围,满 足） 。加强筋采用 14200mm，趾板混凝土配筋率约 85.5kg/m。 趾板用砂浆锚杆与基岩连接，锚筋用 25，长 4.0m，间距 1.0m，每排三根。 配筋简图见图 2-3： 课件之家精心整理资料-欢迎你的欣赏 课件之家
11、精心整理资料-欢迎你的欣赏 图 2-3 趾板配筋简图 23 结论 根据计算结果并结合已建工程经验，确定面板厚度为 0.3m，采用单层双向配 筋，钢筋布置在面板中部。连接板宽取 6m，趾板宽左岸取为 4m，右岸 7m，相应 厚度取为 0.6m 和 0.5m。配筋图见图 2-3。 第三章 坝体稳定和变形计算 根据面板堆石坝设计导则DL501693 的规定，结合本工程地质，填料和 地震烈度条件，本工程上水库主副坝可不作稳定分析和应力应变计算。但对于坝 体设计而言，为了更加全面掌握设计程序，这里应对坝体的稳定和变形进行初步 分析。 31 边坡稳定计算 钢筋混凝土面板堆石坝的稳定分析根据工程经验及设计理
12、论,应进行施工期 和蓄水发电期稳定分析,对于本工程而言,由于坝体施工采用的是分层碾压法,整 个坝高在施工期都在不断变化,而且坝体受的汽车荷载、施工机器荷载和人群荷 载，具体考虑起来比较复杂，并且数量级相对较小。要准确分析整个施工期坝体 稳定性比较困难。所以，在这里只计算竣工期的坝体稳定，此时期的坝体稳定分 析应分为上游坡、下游坡两种情况，但对于两种情况进行比较可知，由于上游坡 有钢筋混凝土面板的作用，其整体性要好于干砌块石及次堆石区组成的下游坡， 再加上上游坡库盆需用弃石碴回填至 375m 高程，由于这部分土的作用，上游坡 的稳定性也要好于下游坡。另外，对于钢筋混凝土面板坝其防渗性能很好，其面
13、 板后的浸润线很低，这对上、下游坡的影响都不大，不是控制性工况。而对蓄水 泰安抽水蓄能电站水利枢纽-上水库库盆及导流建筑物设计 - 14 - 期，在施工质量得到保证的前提下，坝体的上游坡由于水的作用，其稳定性会增 大，对于下游坡的稳定性变化会很小。 根据上述分析：上游坡的稳定性要强于下游坡，所以这里只对竣工期下游坡 的稳定性进行分析。 目前对于土石坝边坡稳定性，目前的稳定分析仍基于极限平衡理论，采用假 定滑动面的方法。依据滑弧的不同型式，可分为圆弧滑动法，折线滑动法和复合 滑动法。对于非粘性土坝坡，这里主要采用折线滑动面法，用自编的 FORTRAN 程 序计算抗滑稳定系数。 311 计算公式
14、3111 正常运用 计算示意图见图 3-1： 1 1 2 1: 2 1: 1 1 图 3-1 正常运用非粘性土坝坡稳定计算示意图 DEBC 土块的平衡式为： （3-0cos 1 sin 111111 tgW K WP c 1） 式中意义见图 4-1 21 ， ADE 土块的平衡式为： 课件之家精心整理资料-欢迎你的欣赏 课件之家精心整理资料-欢迎你的欣赏 0)cos(sin)sin( 1 cos 1 211222211222 PWtgP K tgW K cc （3- 2） 式中： 条块间作用力 1 P 各条块重量 21,W W 土体底面的抗剪强度 21 ， 3112 非常运用（地震荷载） 这里
15、主要考虑水平向地震惯性力对坝坡稳定的影响。 水平向地震惯性力代表值可统一用下式表示： （3-3 ） g aGa F iEih i 式中：为作用在质点 的水平向地震惯性力代表值 i Fi 为水平向设计地震加速度代表值，参照水工建筑物表 2- h a 24。对本工程，设计烈度 7 度，gah1 . 0 为重力加速度g 为地震作用的效应折减系数，一般25 . 0 为集中在质点 的重力作用标准值 Ei G i 为质点 的动态分布系数。对于土石坝，按图 3-2 取值。图中 i a i i a 在设计烈度为 7 度，8 度，9 度时，分别为 3.0，2.5，2.0。 m a 泰安抽水蓄能电站水利枢纽-上水
16、库库盆及导流建筑物设计 - 16 - (a)坝高40m (b) 坝高40m 图 3-2 土石坝动态系数分布 i a 具体计算示意图见图 4-3 1 1 2 1: 2 1: 1 1 1 2 图 3-3 非常运用非粘性土坝坡稳定计算示意图 DEBC 土块的平衡式为： （3-0cos)sincos( 1 sin 1111111111 FtgFW K WP c 4） ADE 土块的平衡式为： 课件之家精心整理资料-欢迎你的欣赏 课件之家精心整理资料-欢迎你的欣赏 0cos)cos(sin)sin( 1 )sincos 1 2221122221122222 FPWtgP K tgFW K cc （ （3
17、-5） 式中： ，同正常运用情况 1 P 21,W W 21 ， 作用于各条块上的最不利水平地震荷载（图 4-3 中均指向迎水 21,F F 面） 3.1.2 计算程序 边坡稳定计算的 Fortran 语言程序 折线法计算程序如下: *PROGRAM MAIN *变量说明 *T-坝顶宽度 ， DELT-A 的步长 ， DELB-B 在 Y 方向 的步长 ， DELA-B 在 X 方向的步长 *A（50）存放 A ， B（50，500，2）存放 B， K（50，500）存 放 K， *ZXIA0K（50）存放每个 A 对应下的最小 K， BK（50，2）用来存 放 ZXIAOK 对应下 B 的座
18、标， *（CX，CY）-坝的最低点的座标 H1，H2-坝顶点及坝 下游最低点所对应的高程 * H3-坝底的控制高程 * R1，R2-坝体材料的容重 MOJIAO1，MOJIAO2- -主堆石及次堆石的内摩擦角 * C-堆石的粘聚力 GFK（10）每一分块 的重量 * XFK（10） ，YFK（10）-分块界点 FKJIAO（10）-每分块 的倾角 * FKMO（10）-每分块的内摩擦角 LENGTH（10）-每分 块的底边长度 泰安抽水蓄能电站水利枢纽-上水库库盆及导流建筑物设计 - 18 - * AA（4） ，BB（4）-四条直线的的系数 * Y1，Y2，Y3，Y4-通过 B 点的垂直线与四
19、条直线的交点的纵座 标 * N-分块数 JSK-用来计数 * I，J -循环变量 M-用来表示 A 的个数 * PO-下游坡度 * KH，CZ，AI 是用来表示地震力系数 REAL T,T1,DELT,DELB,DELA REAL A(50),B(50,500,2),K(50,500) REAL ZXIAOK(50),BK(50,2) REAL XFK(0:10),YFK(0:10) REAL LENGTH(10),FKMO(10),FKJIAO(10),GFK(10) REAL CX,CY,H1,H2,H3,PO REAL R1,R2,MOJIAO1,MOJIAO2,C REAL AA(4)
20、,BB(4) REAL Y1,Y2,Y3,Y4 REAL KH,CZ,AI INTEGER I,J,M,JSK,N,KK * READ(*,*)KK *如果 KK 等于 0，则表示不考虑地震力，如果 KK 等于 1，则表示考 虑地震力。 IF(KK.EQ.0) THEN OPEN(0,FILE=W1.DAT) OPEN(1,FILE=W11.OUT) OPEN(2,FILE=W12.OUT) END IF IF(KK.EQ.1) THEN 课件之家精心整理资料-欢迎你的欣赏 课件之家精心整理资料-欢迎你的欣赏 OPEN(0,FILE=W2.DAT) OPEN(1,FILE=W21.OUT) O
21、PEN(2,FILE=W22.OUT) END IF * READ(0,*)T,DELT,DELA,DELB M=INT(T/DELT) T1=M*DELT IF(T.GT.T1) M=M+2 IF(T.EQ.T1) M=M+1 READ(0,*) (A(I),I=1,M) READ(0,*) H1,H2,H3 READ(0,*) R1,R2,MOJIAO1 ,MOJIAO2,C do i=1,2 READ(0,*)AA(I),BB(I) end do READ(0,*)KH,CZ,AI * PO=1.0/AA(1) CY=H1-H3 CX=PO*CY+T AA(3)=(H2-H3)/CX B
22、B(3)=H1-H2 WRITE(*,*)(A(I),I=1,M), CX,CY MOJIAO1=MOJIAO1*3.1415926/180 MOJIAO2=MOJIAO2*3.1415926/180 * DO 10 I=1,M 泰安抽水蓄能电站水利枢纽-上水库库盆及导流建筑物设计 - 20 - AA(4)=CY/(CX-A(I) BB(4)=-A(I)*AA(4) JSK=1 B(I,JSK,1)=A(I) 110 B(I,JSK,2)=-DELB B(I,JSK,1)=B(I,JSK,1)+DELA IF(B(I,JSK,1).GE.CX) THEN JSK=JSK-1 WRITE(*,*
23、) I,JSK GOTO 140 END IF B(I,JSK,2)=B(I,JSK,2)+DELB 130 Y1=AA(1)*B(I,JSK,1)+BB(1) Y2=AA(2)*B(I,JSK,1)+BB(2) Y3=AA(3)*B(I,JSK,1)+BB(3) Y4=AA(4)*B(I,JSK,1)+BB(4) 100 IF(B(I,JSK,2).LT.Y4) THEN B(I,JSK,2)=B(I,JSK,2)+DELB GOTO 100 ELSE GOTO 120 END IF 120 IF(B(I,JSK,2).GT.Y3) THEN GOTO 110 END IF * *在 FEN
24、KUAI 子程序中，利用坝顶的点 A、折线与主次堆石分界线的 交点（XJIAO，YJIAO） 课件之家精心整理资料-欢迎你的欣赏 课件之家精心整理资料-欢迎你的欣赏 *折点 B、坝脚点 C 将滑动体分块，分别求出每分块点的 X、Y 座标、 每分块体底的长度、 * 每分块体底的内摩擦角及与水平线的夹角 CALL FENKUAI(B(I,JSK,1),B(I,JSK,2),A(I),CX,CY,N,AA, $BB,T,Y2,MOJIAO1,MOJIAO2, FKMO,FKJIAO,LENGTH,XFK,YFK) *在 QIUZHOHGLI 子程序中求出了每分块体的重力 CALL QIUZHONGL
25、I(GFK,XFK,YFK,R1,R2,B(I,JSK,1),B(I,JSK,2), $Y2,T,AA,BB,N,A(I) *在 QIUK 子程序中求出了折点为 B 时的安全系数 CALL QIUK(K(I,JSK),GFK,N,FKMO,FKJIAO,LENGTH,C,JSK,I, $B(I,JSK,1),B(I,JSK,2),KH,CZ,AI) * JSK=JSK+1 B(I,JSK,1)=B(I,JSK-1,1) B(I,JSK,2)=B(I,JSK-1,2)+DELB GOTO 130 * *此段用来求出在 A 点不动、B 点循环所的的最小安全系数 140 ZXIAOK(I)=K(I,
26、1) BK(I,1)=B(I,1,1) BK(I,2)=B(I,1,2) DO 20 J=2,JSK IF(K(I,J).LT.ZXIAOK(I) THEN ZXIAOK(I)=K(I,J) BK(I,1)=B(I,J,1) BK(I,2)=B(I,J,2) END IF 20 CONTINUE 泰安抽水蓄能电站水利枢纽-上水库库盆及导流建筑物设计 - 22 - WRITE(1,*) ZXIAOK(I), A(I),BK(I,1),BK(I,2) *88 10CONTINUE CLOSE(1) CLOSE(2) CLOSE(3) END * * *分块子程序 *变量说明（与主程序对应的变量不说
27、明） * （XJIAO，YJIAO）-折线与主次堆石分界线的交点座标， DELX、DELY 分块点的步长 * I-循环变量 SUBROUTINE FENKUAI(X,Y,A,CX,CY,N,AA,BB,T,Y2,MOJIAO1, $MOJIAO2,FKMO,FKJIAO,LENGTH,XFK,YFK) REAL X,Y,A,CX,CY,AA(3),BB(3),T,Y2,MOJIAO1,MOJIAO2 REAL FKMO(10),FKJIAO(10),LENGTH(10),XFK(0:10),YFK(0:10) REAL XJIAO,YJIAO,DELX,DELY INTEGER N,I * *
28、 * * IF(Y.GT.Y2.AND.A.LE.T) THEN CALL QIUZUOBIAO(XJIAO,YJIAO,FKMO,FKJIAO,N,A,X,Y $,CX,CY,MOJIAO1,MOJIAO2,Y2,T,AA,BB) XFK(0)=A 课件之家精心整理资料-欢迎你的欣赏 课件之家精心整理资料-欢迎你的欣赏 YFK(0)=0 XFK(1)=X YFK(1)=Y XFK(2)=XJIAO YFK(2)=YJIAO DELX=(CX-XJIAO)/(N-2) DELY=(CY-YJIAO)/(N-2) DO 210 I=3,N XFK(I)=XFK(I-1)+DELX YFK(I)=
29、YFK(I-1)+DELY 210 CONTINUE END IF * * IF(Y.EQ.Y2.AND.A.LE.T) THEN CALL QIUZUOBIAO(XJIAO,YJIAO,FKMO,FKJIAO,N,A,X,Y $,CX,CY,MOJIAO1,MOJIAO2,Y2,T,AA,BB) XFK(0)=A YFK(0)=0 XFK(1)=X YFK(1)=Y DELX=(CX-XJIAO)/(N-1) DELY=(CY-YJIAO)/(N-1) DO 220 I=2,N XFK(I)=XFK(I-1)+DELX YFK(I)=YFK(I-1)+DELY 220 CONTINUE EN
30、D IF 泰安抽水蓄能电站水利枢纽-上水库库盆及导流建筑物设计 - 24 - * * IF(Y.LT.Y2.AND.A.LT.T) THEN CALL QIUZUOBIAO(XJIAO,YJIAO,FKMO,FKJIAO,N,A,X,Y $,CX,CY,MOJIAO1,MOJIAO2,Y2,T,AA,BB) XFK(0)=A YFK(0)=0 XFK(1)=XJIAO YFK(1)=YJIAO * * IF(X.LE.2*XJIAO/3.0+CX/3.0) THEN XFK(2)=X YFK(2)=Y DELX=(CX-X)/(N-2) DELY=(CY-Y)/(N-2) DO 230 I=3
31、,N XFK(I)=XFK(I-1)+DELX YFK(I)=YFK(I-1)+DELY 230 CONTINUE END IF * * IF(CX/3.0+2*XJIAO/3.0.LT.X.AND.X.LT.2*CX/3.0+XJIAO/3.0) THEN DELX=(X-XJIAO)/3.0 DELY=(Y-YJIAO)/3.0 DO 240 I=2,N-3 XFK(I)=XFK(I-1)+DELX YFK(I)=YFK(I-1)+DELY 240 CONTINUE 课件之家精心整理资料-欢迎你的欣赏 课件之家精心整理资料-欢迎你的欣赏 DELX=(CX-X)/3.0 DELY=(CY-Y
32、)/3.0 DO 250 I=N-2,N XFK(I)=XFK(I-1)+DELX YFK(I)=YFK(I-1)+DELY 250 CONTINUE END IF * * IF(X.GE.2*CX/3.0+XJIAO/3.0) THEN DELX=(X-XJIAO)/4.0 DELY=(Y-YJIAO)/4.0 DO 260 I=2,N-1 XFK(I)=XFK(I-1)+DELX YFK(I)=YFK(I-1)+DELY 260 CONTINUE XFK(N)=CX YFK(N)=CY END IF * * END IF * * IF(Y.LT.Y2.AND.A.EQ.T) THEN CA
33、LL QIUZUOBIAO(XJIAO,YJIAO,FKMO,FKJIAO,N,A,X,Y $,CX,CY,MOJIAO1,MOJIAO2,Y2,T,AA,BB) XFK(0)=A YFK(0)=0 IF(X.LE.2*XJIAO/3.0+CX/3.0) THEN 泰安抽水蓄能电站水利枢纽-上水库库盆及导流建筑物设计 - 26 - XFK(1)=X YFK(1)=Y DELX=(CX-X)/(N-1) DELY=(CY-Y)/(N-1) DO 270 I=2,N XFK(I)=XFK(I-1)+DELX YFK(I)=YFK(I-1)+DELY 270 CONTINUE END IF * *
34、IF(CX/3.0+2*XJIAO/3.0.LT.X.AND.X.LT.2*CX/3.0+XJIAO/3.0) THEN DELX=(X-XJIAO)/3.0 DELY=(Y-YJIAO)/3.0 DO 280 I=1,N-3 XFK(I)=XFK(I-1)+DELX YFK(I)=YFK(I-1)+DELY 280 CONTINUE DELX=(CX-X)/3.0 DELY=(CY-Y)/3.0 DO 290 I=N-2,N XFK(I)=XFK(I-1)+DELX YFK(I)=YFK(I-1)+DELY 290 CONTINUE END IF * * IF(X.GE.2*CX/3.0+X
35、JIAO/3.0) THEN DELX=(X-XJIAO)/4.0 课件之家精心整理资料-欢迎你的欣赏 课件之家精心整理资料-欢迎你的欣赏 DELY=(Y-YJIAO)/4.0 DO 200 I=1,N-1 XFK(I)=XFK(I-1)+DELX YFK(I)=YFK(I-1)+DELY 200 CONTINUE XFK(N)=CX YFK(N)=CY END IF * * END IF DO 201 I=1,N LENGTH(I)=SQRT(XFK(I)-XFK(I-1)*2+(YFK(I)-YFK(I-1)*2) 201 CONTINUE END * * *求交点，求每分块的倾角，求每分
36、块的内摩擦角 *变量说明 与 FENKUAI 子程序所对应的变量不说明 * JIAO1、JIAO2 分别用来代表折线与水平线的交角 * DELX1、DELX2、DELY1、DELY2 分别用来代表 A 点、B 点、C 点之间 的 X、Y 座标的差值 * DAIA、DAIX、DAICX 分别表示用二分法求（XJIAO、YJIAO）时的中 间变量 SUBROUTINE QIUZUOBIAO(XJIAO,YJIAO,FKMO,FKJIAO,N,A,X,Y $,CX,CY,MOJIAO1,MOJIAO2,Y2,T,AA,BB) REAL XJIAO,YJIAO,FKMO(10),FKJIAO(10),
37、A,X,Y REAL CX,CY,MOJIAO1,MOJIAO2,Y2,T,AA(3),BB(3) REAL JIAO1,JIAO2,DELX1,DELX2,DELY1,DELY2,YY1,YY2 泰安抽水蓄能电站水利枢纽-上水库库盆及导流建筑物设计 - 28 - REAL DAIA,DAIX,DAICX INTEGER N * * DELX1=X-A DELX2=CX-X DELY1=Y DELY2=CY-Y JIAO1=ATAN(DELY1/DELX1) JIAO2=ATAN(DELY2/DELX2) DAIA=A DAIX=X DAICX=CX * * IF(Y.GT.Y2.AND.A.
38、LE.T) THEN N=6 FKMO(1)=MOJIAO1 FKMO(2)=MOJIAO1 FKJIAO(1)=JIAO1 FKJIAO(2)=JIAO2 DO 310 I=3,N FKMO(I)=MOJIAO2 FKJIAO(I)=JIAO2 310 CONTINUE 302 XJIAO=(DAICX+DAIX)/2.0 YY1=AA(2)*XJIAO+BB(2) YY2=TAN(JIAO2)*(XJIAO-CX)+CY 课件之家精心整理资料-欢迎你的欣赏 课件之家精心整理资料-欢迎你的欣赏 IF(YY1-YY2).GE.0.001) THEN DAICX=XJIAO GOTO 302 E
39、ND IF IF(YY1-YY2).LE.-0.001) THEN DAIX=XJIAO GOTO 302 END IF YJIAO=YY1 END IF * IF(Y.EQ.Y2.AND.A.LE.T) THEN N=5 XJIAO=X XJIAO=Y IF(A.EQ.T) FKMO(1)=(MOJIAO1+MOJIAO2)/2 IF(A.LT.T) FKMO(1)=MOJIAO1 FKJIAO(1)=JIAO1 DO 320 I=2,N FKMO(I)=MOJIAO2 FKJIAO(I)=JIAO2 320 CONTINUE END IF * IF(Y.LT.Y2.AND.A.LT.T)
40、THEN 301 XJIAO=(DAIA+DAIX)/2.0 YY1=AA(2)*XJIAO+BB(2) YY2=TAN(JIAO1)*(XJIAO-X)+Y 泰安抽水蓄能电站水利枢纽-上水库库盆及导流建筑物设计 - 30 - IF(YY1-YY2).GT.0.01) THEN DAIX=XJIAO GOTO 301 END IF IF(YY1-YY2).LT.-0.01) THEN DAIA=XJIAO GOTO 301 END IF YJIAO=YY1 FKMO(1)=MOJIAO1 FKJIAO(1)=JIAO1 * * IF(X.LE.2*XJIAO/3.0+CX/3.0) THEN
41、N=6 FKMO(2)=MOJIAO2 FKJIAO(2)=JIAO1 DO 330 I=3,N FKMO(I)=MOJIAO2 FKJIAO(I)=JIAO2 330 CONTINUE END IF * * IF(CX/3.0+2*XJIAO/3.0.LT.X.AND.X.LT.2*CX/3.0+XJIAO/3.0) THEN N=7 DO 340 I=2,N-3 FKMO(I)=MOJIAO2 FKJIAO(I)=JIAO1 课件之家精心整理资料-欢迎你的欣赏 课件之家精心整理资料-欢迎你的欣赏 340 CONTINUE DO 350 I=N-2,N FKMO(I)=MOJIAO2 FKJ
42、IAO(I)=JIAO2 350 CONTINUE END IF * * IF(X.GE.2*CX/3.0+XJIAO/3.0) THEN N=6 DO 360 I=2,N-1 FKMO(I)=MOJIAO2 FKJIAO(I)=JIAO1 360 CONTINUE FKMO(N)=MOJIAO2 FKJIAO(N)=JIAO2 END IF * * END IF * * IF(Y.LT.Y2.AND.A.EQ.T) THEN XJIAO=A YJIAO=0 * * IF(X.LE.2*XJIAO/3.0+CX/3.0) THEN N=5 FKMO(1)=MOJIAO2 FKJIAO(1)=JIAO1 泰安抽水蓄能电站水利枢纽-上水库库盆及导流建筑物设计 - 32 - DO 370 I=2,N FKMO(I)=MOJIAO2 FKJIAO(I)=JIAO2 370 CONTINUE END IF * * IF(